Skip to content

Commit

Permalink
Updated JOSS paper: updated the chapter reference to in press, added …
Browse files Browse the repository at this point in the history
…a few more references, and clarified the faster nature and future plans for pymecht
  • Loading branch information
ankushaggarwal committed Dec 19, 2024
1 parent f38be0d commit bf9670d
Show file tree
Hide file tree
Showing 2 changed files with 52 additions and 11 deletions.
57 changes: 49 additions & 8 deletions joss-paper/paper.bib
Original file line number Diff line number Diff line change
Expand Up @@ -20,11 +20,52 @@ @article{AGGARWAL2023105657
author={Aggarwal, Ankush and Hudson, Luke T and Laurence, Devin W and Lee, Chung-Hao and Pant, Sanjay}
}

@bookchapter{AGGARWAL2024,
title = {Heterogeneity and multi-scale modeling in vascular biomechanics},
booktitle = {Integration and Bridging of Multiscale Bioengineering Designs and Tissue Biomechanics},
publisher = {Springer},
year = {under review},
author={Aggarwal, Ankush and Coccarelli, Alberto and McGinty, Sean},
editor={Liao, Jun and Wong, Joyce Y}
}
@bookchapter{Aggarwal2025,
author = {Aggarwal, Ankush and Coccarelli, Alberto and McGinty, Sean},
title = {Heterogeneity and multi-scale modeling in vascular biomechanics},
booktitle = {Integration and Bridging of Multiscale Bioengineering Designs and Tissue Biomechanics},
editor = {Jun Liao and Joyce Y. Wong},
publisher = {Springer Nature Switzerland},
year = {2025},
%pages = {xx--yy},
isbn = {978-3-031-81742-7},
note = {In press},
}

@article{coccarelli2021framework,
title={A framework for incorporating 3D hyperelastic vascular wall models in 1D blood flow simulations},
author={Coccarelli, Alberto and Carson, Jason M and Aggarwal, Ankush and Pant, Sanjay},
journal={Biomechanics and Modeling in Mechanobiology},
volume={20},
number={4},
pages={1231--1249},
year={2021},
publisher={Springer},
doi = {10.1007/s10237-021-01437-5}
}

@software{alberto_coccarelli_2021_4522152,
author = {Alberto Coccarelli and
Jason M Carson and
Ankush Aggarwal and
Sanjay Pant},
title = {1D-Hyperelastic-Haemodynamics: Version1},
month = feb,
year = 2021,
publisher = {Zenodo},
version = {V1},
doi = {10.5281/zenodo.4522152},
}

@article{coccarelli2024new,
title={A new model for evaluating pressure-induced vascular tone in small cerebral arteries},
author={Coccarelli, Alberto and Pant, Sanjay and Polydoros, Ioannis and Harraz, Osama F},
journal={Biomechanics and Modeling in Mechanobiology},
volume={23},
number={1},
pages={271--286},
year={2024},
publisher={Springer},
doi={10.1007/s10237-023-01774-7}
}

6 changes: 3 additions & 3 deletions joss-paper/paper.md
Original file line number Diff line number Diff line change
Expand Up @@ -113,7 +113,7 @@ Detailed documentation is hosted on [`readthedocs`](https://pymecht.readthedocs.
# Advantages over finite element simulation
In principle, the problems that can be solved using `pyMechT` can also be solved using any finite element simulation software. However, `pyMechT` offers the following advantages:

- Geometry and mesh creation would be required for a finite element simulation, which usually takes some time. However, the pre-defined geometrical features in `pyMechT` means that one only needs to choose the right class and parameters. In addition, no meshing is required. This means that setting up the problem is much faster.
- Geometry and mesh creation would be required for a finite element simulation, which usually takes some time. However, the pre-defined geometrical features in `pyMechT` means that one only needs to choose the right class and parameters. In addition, no meshing is required. This means that setting up the problem (i.e., defining the geometry/mesh and loads) is much faster in `pyMechT`. Once the model has been setup, the computational time required to solve it is comparable, depending on the finite element mesh density.

- Enforcing incompressibility in a finite element simulation can be numerically challenging, necessitating approaches such as Lagrange multiplier with a three-field formulation. Instead, in `pyMechT`, the incompressibility is analytically enforced *exactly*, thus making the results more robust.

Expand All @@ -122,9 +122,9 @@ In principle, the problems that can be solved using `pyMechT` can also be solved
- The reference zero-stress state of biological tissues can be unknown or ambiguous. Moreover, the biological tissues are heterogeneous, with multiple layers each of varying properties. These aspects are non-trivial to incorporate in a finite element simulation, due to the need for recreating the geometry and/or incompatability of the initial state. However, it is straightforward to simulate these in `pyMechT`.

# Uses in literature
`pyMechT` has been used for Bayesian model selection based on extensive planar biaxial extension data [@AGGARWAL2023105657]. This work required rapid simulation of varied constitutive models, which was facilitated by `pyMechT`. Similarly, the Bayesian inference via Markov Chain Monte Carlo in `pyMechT` was used to infer the distribution of aortic biomechanical and geometrical properties based on in-vivo measurements (as likelihood) and ex-vivo biaxial extension data (as prior distribution) [@AGGARWAL2024]. Moreover, data-driven model developed in @AGGARWAL2023115812 has been used in `pyMechT` via the `splineI1` and `splineI1I4` material models.
`pyMechT` has been used for Bayesian model selection based on extensive planar biaxial extension data [@AGGARWAL2023105657]. This work required rapid simulation of varied constitutive models, which was facilitated by `pyMechT`. Similarly, the Bayesian inference via Markov Chain Monte Carlo in `pyMechT` was used to infer the distribution of aortic biomechanical and geometrical properties based on in-vivo measurements (as likelihood) and ex-vivo biaxial extension data (as prior distribution) [@Aggarwal2025]. Moreover, data-driven model developed in @AGGARWAL2023115812 has been used in `pyMechT` via the `splineI1` and `splineI1I4` material models.

# Conclusion and future plans
`pyMechT` fills an important gap and allows soft tissue biomechanics researchers to model ex-vivo testing setups in a fast, robust, and flexible manner. The package is numerically efficient and extensively documented. It has facilitated several publications, and we believe that it can benefit the wider community. In the future, we plan to extend the capabilities of the package to include more material models, such as inelastic (viscoelastic, plastic, damage, growth & remodeling), and other ex-vivo setups (such as microindentation). Lastly, the package could be coupled with others to allow multi-physics simulations, such as for hemodynamics and biochemical regulation.
`pyMechT` fills an important gap and allows soft tissue biomechanics researchers to model ex-vivo testing setups in a fast, robust, and flexible manner. The package is numerically efficient and extensively documented. It has facilitated several publications, and we believe that it can benefit the wider community. In the future, we plan to extend the capabilities of the package to include more material models, such as inelastic (viscoelastic, plastic, damage, growth & remodeling) pre-defined formulations, and other ex-vivo setups (such as microindentation using Hertz contact model). Lastly, the package could be coupled with others to allow multi-physics simulations, such as for hemodynamics [@coccarelli2021framework, @alberto_coccarelli_2021_4522152] and biochemical regulation [@coccarelli2024new].

# References

0 comments on commit bf9670d

Please sign in to comment.