Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA direct Eulerian GRP scheme for a blood flow model in arteries
172
0
0.0
(
0
)
Added by
Qinglong Zhang
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this paper, we propose a direct Eulerian generalized Riemann problem (GRP) scheme for a blood flow model in arteries. It is an extension of the Eulerian GRP scheme, which is developed by Ben-Artzi, et. al. in J. Comput. Phys., 218(2006). By using the Riemann invariants, we diagonalize the blood flow system into a weakly coupled system, which is used to resolve rarefaction wave. We also use Rankine-Hugoniot condition to resolve the local GRP formulation. We pay special attention to the acoustic case as well as the sonic case. The extension to the two dimensional case is carefully obtained by using the dimensional splitting technique. We test that the derived GRP scheme is second order accuracy.
rate research
Read More
The paper proposes a second-order accurate direct Eulerian generalized Riemann problem (GRP) scheme for the radiation hydrodynamical equations (RHE) in the zero diffusion limit. The difficulty comes from no explicit expression of the flux in terms of the conservative vector. The characteristic fields and the relations between the left and right states across the elementary-waves are first studied, and then the solution of the one-dimensional Riemann problem is analyzed and given. Based on those, the direct Eulerian GRP scheme is derived by directly using the generalized Riemann invariants and the Runkine-Hugoniot jump conditions to analytically resolve the left and right nonlinear waves of the local GRP in the Eulerian formulation. Several numerical examples show that the GRP scheme can achieve second-order accuracy and high resolution of strong discontinuity.
A new two-dimensional model for blood flows in arteries with arbitrary cross sections is derived. The model consists of a system of balance laws for conservation of mass and balance of momentum in the axial and angular directions. The equations are derived by applying asymptotic analysis to the incompressible Navier-Stokes equations in narrow, large vessels and integrating in the radial direction in each cross section. The main properties of the system are discussed and a positivity-preserving well-balanced central-upwind scheme is presented. The merits of the scheme will be tested in a variety of scenarios. In particular, numerical results of simulations using an idealized aorta model are shown. We analyze the time evolution of the blood flow under different initial conditions such as perturbations to steady states consisting of a bulging in the vessels wall. We consider different situations given by distinct variations in the vessels elasticity.
We are interested in simulating blood flow in arteries with a one dimensional model. Thanks to recent developments in the analysis of hyperbolic system of conservation laws (in the Saint-Venant/ shallow water equations context) we will perform a simple finite volume scheme. We focus on conservation properties of this scheme which were not previously considered. To emphasize the necessity of this scheme, we present how a too simple numerical scheme may induce spurious flows when the basic static shape of the radius changes. On contrary, the proposed scheme is well-balanced: it preserves equilibria of Q = 0. Then examples of analytical or linearized solutions with and without viscous damping are presented to validate the calculations. The influence of abrupt change of basic radius is emphasized in the case of an aneurism.
We performed numerical simulations of blood flow in arteries with a variable stiffness and cross-section at rest using a finite volume method coupled with a hydrostatic reconstruction of the variables at the interface of each mesh cell. The method was then validated on examples taken from the literature. Asymptotic solutions were computed to highlight the effect of the viscous and viscoelastic source terms. Finally, the blood flow was computed in an artery where the cross-section at rest and the stiffness were varying. In each test case, the hydrostatic reconstruction showed good results where other simpler schemes did not, generating spurious oscillations andnonphysical velocities.
We propose an integrated electromechanical model of the human heart, with focus on the left ventricle, wherein biophysically detailed models describe the different physical phenomena concurring to the cardiac function. We model the subcellular generation of active force by means of an Artificial Neural Network, which is trained by a suitable Machine Learning algorithm from a collection of pre-computed numerical simulations of a biophysically detailed, yet computational demanding, high-fidelity model. To provide physiologically meaningful results, we couple the 3D electromechanical model with a closed-loop 0D (lumped parameters) model describing the blood circulation in the whole cardiovascular network. We prove that the 3D-0D coupling of the two models is compliant with the principle of energy conservation, which is achieved in virtue of energy-consistent boundary conditions that account for the interaction among cardiac chambers within the computational domain, pericardium and surrounding tissue. We thus derive an overall balance of mechanical energy for the 3D-0D model. This provides a quantitative insight into the energy utilization, dissipation and transfer among the different compartments of the cardiovascular network and during different stages of the heartbeat. In virtue of this new model and the energy balance, we propose a new validation tool of heart energy usage against relationships used in the daily clinical practice. Finally, we provide a mathematical formulation of an inverse problem aimed at recovering the reference configuration of one or multiple cardiac chambers, starting from the stressed configuration acquired from medical imaging. This is fundamental to correctly initialize electromechanical simulations. Numerical methods and simulations of the 3D-0D model will be detailed in Part II.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
