seongwoohan/heartmodel

Computer Models for Optimizing Interventions in Refractory Pulmonary Hypertension

Introduction

We provide computer models, optimizing three proposed surgical treatments of refractory pulmonary hypertension. The models set up the methodology to examine the treatment results without actually performing the surgery. Numerical simulations are performed to compare the change of pressure and oxygen saturation in pulmonary artery and systemic artery to examine the impact of the treatments.

This software is build upon https://www.math.nyu.edu/faculty/peskin/ModSimPrograms/ch1/. See here for original code of computer simulation in heart and circulation.

Experiment

Run the below files to see the numerical simulation in each intervention:

1. circ_asd.m
2. circ_vsd.m
3. circ_d.m

There is a new python version for the model. Run the below files as well to see the numerical simulation in each intervention: (https://github.com/seongwoohan/HeartModel/tree/master/src/python_ver/original).

1. circ_asd.py
2. circ_vsd.py
3. circ.d.py

Refractory Pulmonary Hypertension

Code and computer models for "Computer simulation of surgical interventions for the treatment of refractory pulmonary hypertension" appeard at Mathematical Medicine and Biology, 2022. Paper at: https://arxiv.org/abs/2101.12191

See the READMEs in subdirectories for directions on running the code.

heartmodel's Issues

make max elastance needs to be smaller

HeartModel/src/matlab_ver/in_circ.m

Line 19 in 8c8f681

EmaxLV = 30.00*1000 % (mmHg/L)

I think the max elastance should be more like ~3e3 mmHg / L

see this paper that we cite in ours: https://link.springer.com/article/10.1007/s11517-009-0449-9