Skip to content

Correction viscosity term in ChamberSphere block#217

Merged
KatrinKoesler merged 6 commits intoSimVascular:masterfrom
KatrinKoesler:chamber_sphere
Feb 27, 2026
Merged

Correction viscosity term in ChamberSphere block#217
KatrinKoesler merged 6 commits intoSimVascular:masterfrom
KatrinKoesler:chamber_sphere

Conversation

@KatrinKoesler
Copy link
Collaborator

Correction viscosity term in ChamberSphere block

Current situation

This addresses the problem highlighted in #216. Calculations indicated that the implementation of the viscosity term might previously have been incorrect. Further calculations have shown this is the case.

Release Notes

The viscosity term has been updated from 2 * eta * dCG * (1 - 2 * C**(-6)) to eta * dCG * (1 + 2 * C**(-6)). The test file for the ChamberSphere test case has been adjusted according to the new outputs. The overall behavior has not changed significantly due to this change, which is expected since the term has a much lower contribution than the elastic strain energy term.

Documentation

Has been changed accordingly with the new formulation of the viscosity term.

Testing

Test case has been adjusted accordingly, all tests run.

Code of Conduct & Contributing Guidelines

@KatrinKoesler KatrinKoesler marked this pull request as ready for review February 26, 2026 20:40
@mrp089 mrp089 self-requested a review February 26, 2026 22:34
@mrp089
Copy link
Member

mrp089 commented Feb 26, 2026

@KatrinKoesler, nicely done! Can you add a screenshot of your before/after comparison for reference? Thanks!

@KatrinKoesler
Copy link
Collaborator Author

Comparison_new_viscosity pv_loop_new_viscosity

Plotting a pv-loop and pressures/volume/in- and outflow over time in a cardiac cycle

@KatrinKoesler KatrinKoesler merged commit 991fa17 into SimVascular:master Feb 27, 2026
10 checks passed
@KatrinKoesler KatrinKoesler deleted the chamber_sphere branch February 27, 2026 02:20
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants