A stable high-order FC-based methodology for hemodynamic wave propagation

JOURNAL OF COMPUTATIONAL PHYSICS

卷:405

Article

Amlani, Faisal¹  Pahlevan, Niema M.^1,2 

[1] Univ Southern Calif, Dept Aerosp & Mech Engn, Los Angeles, CA 90007 USA

[2] Univ Southern Calif, Sch Med, Los Angeles, CA 90007 USA

关键词: Fourier continuation;    Hemodynamic wave simulation;    Mathematical physiology;    High-order methods;   

DOI  :  10.1016/j.jcp.2019.109130

来源: Elsevier

【 摘 要 】

A high-order numerical algorithm is proposed for the solution of one-dimensional arterial pulse wave propagation problems based on use of an accelerated Fourier continuation (FC) methodology for accurate Fourier expansion of non-periodic functions. The solver provides high-order accuracy, mild Courant-Friedrichs-Lewy (CFL) constraints on the time discretization and, importantly, results that are essentially free of spatial dispersion errors-enabling fast and accurate resolution of clinically-relevant problems requiring simulation of many cardiac cycles or vascular segments. The left ventricle-arterial model that is employed presents a particularly challenging case of ordinary differential equation (ODE)-governed boundary conditions that include a hybrid ODE-Dirichlet model for the left ventricle and an ODE-based Windkessel model for truncated vasculature. Results from FC-based simulations are shown to capture the important physiological features of pressure and flow waveforms in the systemic circulation. The robustness of the proposed solver is demonstrated through a number of numerical examples that include performance studies and a physiologically-accurate case study of the coupled left ventricle-arterial system. The results of this paper imply that the FC-based methodology is straightforwardly applicable to other biological and physical phenomena that are governed by similar hyperbolic partial differential equations (PDEs) and ODE-based time-dependent boundaries. (C) 2019 Elsevier Inc. All rights reserved.

【 授权许可】

Free