Numerical Method of Characteristics for One-Dimensional Blood Flow

Area of Science:

• Computational fluid dynamics
• Biomedical engineering
• Mathematical modeling of physiological systems

Background:

• Accurate cardiovascular system modeling requires solving complex one-dimensional nonlinear hyperbolic systems.
• Existing methods like finite elements and discontinuous Galerkin are computationally intensive.
• Efficient methods are needed for real-time applications, multi-cycle analysis, and uncertainty quantification.

Purpose of the Study:

• To present an efficient and unconditionally stable numerical method for approximating solutions to diagonal nonlinear hyperbolic systems.
• To address computational cost and stability limitations in cardiovascular flow modeling.
• To enable advanced applications such as real-time clinical decision support and uncertainty quantification.

Main Methods:

• Development and theoretical analysis of a novel numerical algorithm for hyperbolic systems.
• Comparison of the proposed method against a discontinuous Galerkin implementation.
• Implementation of the method on physiologically relevant small and large arterial networks.

Main Results:

• The proposed method demonstrates unconditional stability, overcoming time-step restrictions.
• The new algorithm offers improved computational efficiency compared to traditional methods.
• Successful application to arterial networks validates its utility for complex cardiovascular simulations.

Conclusions:

• The developed numerical method provides an efficient and stable solution for cardiovascular flow modeling.
• This advancement facilitates more complex and real-time applications in cardiovascular research and clinical practice.
• The method's robustness is confirmed through theoretical analysis and practical implementation on arterial networks.