svZeroDSolver is a Python code that simulates the hemodynamics in zero-dimensional (0D) lumped parameter models of vascular networks. These 0D models are governed by differential algebraic equations (DAEs).
The solver uses a highly modular framework to model the vascular anatomy, using individual 0D elements to represent different parts of
the vascular anatomy (and boundary conditions). The individual 0D elements and their associated governing equations defined in blocks.py.
In solver.py, the blocks are assembled and simulated using the generalized-alpha time-stepping method defined in time_integration.py.
The svZeroDSolver Python files are in the svzerodsolver Python package directory.
svZeroDSolver currently supports the following vascular 0D modeling options and boundary conditions:
- Resistor
- Resistor-capacitor
- Resistor-inductor
- Resistor-capacitor-inductor
- Pressure
- Resistor
- RCR
- Coronary
- Flow
The following software is required:
- Python 3
The following Python packages are required:
- os
- re
- sys
- pdb
- copy
- numpy
- argparse
- tqdm.tqdm
- importlib
- matplotlib.pyplot
- collections.defaultdict
- scipy
- scipy.interpolate
- scipy.sparse.linalg
- scipy.sparse.csr_matrix
The solver can be executed three ways.
- Execute from the top level Git repository using
python -m svzerodsolver.solver SOLVER_INPUT_FILE
- Setting the
PYTHONPATHenvironent variable
The solver can be executed from any directory by setting the PYTHONPATH environent variable
to the top level Git repository
export PYTHONPATH=$PYTHONPATH:/$HOME/svZeroDSolver/
- Installing the
svzerodsolverPython package
The svzerodsolver Python package is installed using the setup.py script
python setup.py install
The options supported by the solver are listed using
python -m svzerodsolver.solver --help