Strain-dependent active stress model - #245
Open
KatrinKoesler wants to merge 44 commits into
Open
KatrinKoesler wants to merge 44 commits into
KatrinKoesler wants to merge 44 commits into
Conversation
…lver into ventricle_gr_144
…ype and CMakeLists File. Removed comments from ChamberSphere_StrainDepActStress files
…based on Caruel 2013
…lver into ventricle_gr_144
…d changed runtime error to warning in case of non-convergence of the non-linear solver.
…tress, changed parameters to Caruel paper parameters in test file for strain dependent active stress
modified: scripts/ChamberSphere_StrainDepActStress.yaml
modified: scripts/jacobian_strainDependence.py
modified: src/model/ChamberSphere_StrainDepActStress.cpp
modified: tests/cases/chamber_sphere_closed_loop_straindep.json
ChamberSphere_StrainDepActStress.cpp: now implemented so strain depenedency works and simulation running
Outputs are pv loop like with decreasing cardiac volume, crashes in 10th cycle
…es for closed loop and Caruel circulations
…ss spherical model (results don't change much but timing wise now corresponds better to Caruel et al. 2013); added test cases with piecewise valves and split vessels to change the order of the resistance and capacitance in the circulation (better end diastolic behavior when compliance before resistance)
…lver into ventricle_gr_144
…unction and not inside ChamberSphere
…active stress models now can be chosen for ChamberSphere via input file
… structure. Added result files for strain-dependent test case
…n case activation_function is missing
![claude[bot]](https://avatars.githubusercontent.com/in/1236702?s=60&v=4){kind=small}
mrp089
requested changes
Sep 17, 2026
| // strain-dependent active-stress model (Caruel 2013) | ||
| // ============================================================ | ||
|
|
||
| double PiecewiseRateActivation::compute(double time) { |
Member
There was a problem hiding this comment.
This differs from other activation functions in that it returns values outside [0,1] and has a physical unit (1/s). Is there a way to consolidate these? Or should the rate activation have its own class?
Collaborator
Author
There was a problem hiding this comment.
The activation function is a function summarizing the chemical activation of the cardiomyocytes and needs to be defined in this way so that the model works. I could make its own class because it doesn't really make sense to use it with the strain-independent active stress models, we currently have implemented
Collaborator
Author
There was a problem hiding this comment.
But if this activation function is parametrized, it could also be used to create activation functions that are piecewise linear and between [0,1]
…rain dependent values via break point vectors in the input file
…ve stress case. Change 6 valve_status values in the results file that now differed from the prior results file by around 1e-6 (these might have be affected due to the numerical instabilities in this test case). all other values remain unchanged to prior test case!
This branch has not been deployed
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
2 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Added strain-dependent active stress model for ChamberSphere
Current situation
Active stress generation in physiology is strain-dependent. There is no strain-dependent active stress model implemented in svZeroDSolver yet. This addresses #244.
Release Notes
Documentation
All additions are documented, documentation check passes in Github.
Testing
All tests pass, including the newly instated strain-dependent active stress test case
Code of Conduct & Contributing Guidelines