Reply to "fixation and perception".

P. Boyce,D. C. West

DOI: https://doi.org/10.1080/713818086

1968-05-01

Abstract:

What problem does this paper attempt to address?

Physics-constrained coupled neural differential equations for one dimensional blood flow modeling

Hunor Csala,Arvind Mohan,Daniel Livescu,Amirhossein Arzani

2024-11-08

Abstract:Computational cardiovascular flow modeling plays a crucial role in understanding blood flow dynamics. While 3D models provide acute details, they are computationally expensive, especially with fluid-structure interaction (FSI) simulations. 1D models offer a computationally efficient alternative, by simplifying the 3D Navier-Stokes equations through axisymmetric flow assumption and cross-sectional averaging. However, traditional 1D models based on finite element methods (FEM) often lack accuracy compared to 3D averaged solutions. This study introduces a novel physics-constrained machine learning technique that enhances the accuracy of 1D blood flow models while maintaining computational efficiency. Our approach, utilizing a physics-constrained coupled neural differential equation (PCNDE) framework, demonstrates superior performance compared to conventional FEM-based 1D models across a wide range of inlet boundary condition waveforms and stenosis blockage ratios. A key innovation lies in the spatial formulation of the momentum conservation equation, departing from the traditional temporal approach and capitalizing on the inherent temporal periodicity of blood flow. This spatial neural differential equation formulation switches space and time and overcomes issues related to coupling stability and smoothness, while simplifying boundary condition implementation. The model accurately captures flow rate, area, and pressure variations for unseen waveforms and geometries. We evaluate the model's robustness to input noise and explore the loss landscapes associated with the inclusion of different physics terms. This advanced 1D modeling technique offers promising potential for rapid cardiovascular simulations, achieving computational efficiency and accuracy. By combining the strengths of physics-based and data-driven modeling, this approach enables fast and accurate cardiovascular simulations.

Fluid Dynamics,Machine Learning
A one-dimensional computational model for blood flow in an elastic blood vessel with a rigid catheter

Aseem Milind Pradhan,Fernando Mut,Juan Raul Cebral

DOI: https://doi.org/10.1002/cnm.3834

Abstract:Strokes are one of the leading causes of death in the United States. Stroke treatment involves removal or dissolution of the obstruction (usually a clot) in the blocked artery by catheter insertion. A computer simulation to systematically plan such patient-specific treatments needs a network of about 105 blood vessels including collaterals. The existing computational fluid dynamic (CFD) solvers are not employed for stroke treatment planning as they are incapable of providing solutions for such big arterial trees in a reasonable amount of time. This work presents a novel one-dimensional mathematical formulation for blood flow modeling in an elastic blood vessel with a centrally placed rigid catheter. The governing equations are first-order hyperbolic partial differential equations, and the hypergeometric function needs to be computed to obtain the characteristic system of these hyperbolic equations. We employed the Discontinuous Galerkin method to solve the hyperbolic system and validated the implementation by comparing it against a well-established 3D CFD solver using idealized vessels and a realistic truncated arterial network. The results showed clinically insignificant differences in steady flow cases, with overall variations between 1D and 3D models remaining below 10%. Additionally, the solver accurately captured wave reflection phenomena at domain discontinuities in unsteady cases. A primary advantage of this model over 3D solvers is its ease in obtaining a discretized geometry of complex vasculatures with multiple arterial branches. Thus, the 1D computational model offers good accuracy and applicability in simulating complex vasculatures, demonstrating promising potential for investigating patient-specific endovascular interventions in strokes.
A modular framework for implicit 3D–0D coupling in cardiac mechanics

Aaron L. Brown,Matteo Salvador,Lei Shi,Martin R. Pfaller,Zinan Hu,Kaitlin E. Harold,Tzung Hsiai,Vijay Vedula,Alison L. Marsden

DOI: https://doi.org/10.1016/j.cma.2024.116764

IF: 6.588

2024-01-20

Computer Methods in Applied Mechanics and Engineering

Abstract:In numerical simulations of cardiac mechanics , coupling the heart to a model of the circulatory system is essential for capturing physiological cardiac behavior. A popular and efficient technique is to use an electrical circuit analogy, known as a lumped parameter network or zero-dimensional (0D) fluid model, to represent blood flow throughout the cardiovascular system. Due to the strong physical interaction between the heart and the blood circulation, developing accurate and efficient numerical coupling methods remains an active area of research. In this work, we present a modular framework for implicitly coupling three-dimensional (3D) finite element simulations of cardiac mechanics to 0D models of blood circulation. The framework is modular in that the circulation model can be modified independently of the 3D finite element solver, and vice versa. The numerical scheme builds upon a previous work that combines 3D blood flow models with 0D circulation models (3D fluid–0D fluid). Here, we extend it to couple 3D cardiac tissue mechanics models with 0D circulation models (3D structure–0D fluid), showing that both mathematical problems can be solved within a unified coupling scheme. The effectiveness, temporal convergence, and computational cost of the algorithm are assessed through multiple examples relevant to the cardiovascular modeling community. Importantly, in an idealized left ventricle example, we show that the coupled model yields physiological pressure–volume loops and naturally recapitulates the isovolumic contraction and relaxation phases of the cardiac cycle without any additional numerical techniques. Furthermore, we provide a new derivation of the scheme inspired by the Approximate Newton Method of Chan (1985), explaining how the proposed numerical scheme combines the stability of monolithic approaches with the modularity and flexibility of partitioned approaches.

mechanics,engineering, multidisciplinary,mathematics, interdisciplinary applications
A Modular Framework for Implicit 3D-0D Coupling in Cardiac Mechanics

Aaron L. Brown,Matteo Salvador,Lei Shi,Martin R. Pfaller,Zinan Hu,Kaitlin E. Harold,Tzung Hsiai,Vijay Vedula,Alison L. Marsden

2023-10-21

Abstract:In numerical simulations of cardiac mechanics, coupling the heart to a model of the circulatory system is essential for capturing physiological cardiac behavior. A popular and efficient technique is to use an electrical circuit analogy, known as a lumped parameter network or zero-dimensional (0D) fluid model, to represent blood flow throughout the cardiovascular system. Due to the strong physical interaction between the heart and the blood circulation, developing accurate and efficient numerical coupling methods remains an active area of research. In this work, we present a modular framework for implicitly coupling three-dimensional (3D) finite element simulations of cardiac mechanics to 0D models of blood circulation. The framework is modular in that the circulation model can be modified independently of the 3D finite element solver, and vice versa. The numerical scheme builds upon a previous work that combines 3D blood flow models with 0D circulation models (3D fluid - 0D fluid). Here, we extend it to couple 3D cardiac tissue mechanics models with 0D circulation models (3D structure - 0D fluid), showing that both mathematical problems can be solved within a unified coupling scheme. The effectiveness, temporal convergence, and computational cost of the algorithm are assessed through multiple examples relevant to the cardiovascular modeling community. Importantly, in an idealized left ventricle example, we show that the coupled model yields physiological pressure-volume loops and naturally recapitulates the isovolumic contraction and relaxation phases of the cardiac cycle without any additional numerical techniques. Furthermore, we provide a new derivation of the scheme inspired by the Approximate Newton Method of Chan (1985), explaining how the proposed numerical scheme combines the stability of monolithic approaches with the modularity and flexibility of partitioned approaches.

Numerical Analysis,Biological Physics
A multicore parallel algorithm for multiscale modelling of an entire human blood circulation network

Jiawei Liu,Hiroshi Suito

2024-04-02

Abstract:The presented multi-scale, closed-loop blood circulation model includes arterial, venous, and portal venous systems, heart-pulmonary circulation, and micro-circulation in capillaries. One-dimensional models simulate large blood vessel flow, whereas zerodimensional models are used for simulating blood flow in vascular subsystems corresponding to peripheral arteries and organs. Transmission conditions at bifurcation and confluence are solved using Riemann invariants. Blood circulation simulation in the portal venous system and related organs (liver, stomach, spleen, pancreas, intestine) is particularly targeted. Those organs play important roles in metabolic system dynamics. The proposed efficient parallel algorithms for multicore environments solve these equations much faster than serial computations.

Medical Physics,Numerical Analysis
A simplified particulate model for coarse-grained hemodynamics simulations

Florian Janoschek,Federico Toschi,Jens Harting

DOI: https://doi.org/10.1103/PhysRevE.82.056710

2010-09-28

Abstract:Human blood flow is a multi-scale problem: in first approximation, blood is a dense suspension of plasma and deformable red cells. Physiological vessel diameters range from about one to thousands of cell radii. Current computational models either involve a homogeneous fluid and cannot track particulate effects or describe a relatively small number of cells with high resolution, but are incapable to reach relevant time and length scales. Our approach is to simplify much further than existing particulate models. We combine well established methods from other areas of physics in order to find the essential ingredients for a minimalist description that still recovers hemorheology. These ingredients are a lattice Boltzmann method describing rigid particle suspensions to account for hydrodynamic long range interactions and---in order to describe the more complex short-range behavior of cells---anisotropic model potentials known from molecular dynamics simulations. Paying detailedness, we achieve an efficient and scalable implementation which is crucial for our ultimate goal: establishing a link between the collective behavior of millions of cells and the macroscopic properties of blood in realistic flow situations. In this paper we present our model and demonstrate its applicability to conditions typical for the microvasculature.

Soft Condensed Matter,Biological Physics,Computational Physics
Simulation of unsteady blood flows in a patient-specific compliant pulmonary artery with a highly parallel monolithically coupled fluid-structure interaction algorithm

Fande Kong,Vitaly Kheyfets,Ender Finol,Xiao-Chuan Cai

DOI: https://doi.org/10.48550/arXiv.1810.04289

2018-10-10

Abstract:Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient-specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems, but need improvements when the problems become harder in the case when (1) the geometry becomes complex (from a few branches to a full pulmonary artery), (2) the model becomes more complex (from fluid-only calculation to coupled fluid-structure interaction calculation), (3) both the fluid and wall models become highly nonlinear, and (4) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid-structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient-specific, full size pulmonary artery obtained from CT (Computed Tomography) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier-Stokes equations and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. The proposed numerical algorithm and software scale well beyond 10,000 processor cores on a supercomputer for solving the fluid-structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.

Computational Physics,Computational Engineering, Finance, and Science,Distributed, Parallel, and Cluster Computing,Numerical Analysis
Towards Blood Flow in the Virtual Human: Efficient Self-Coupling of HemeLB

J. W. S. McCullough,R. A. Richardson,A. Patronis,R. Halver,R. Marshall,M. Ruefenacht,B. J. M. Wylie,T. Odaker,M. Wiedemann,B. Lloyd,E. Neufeld,G. Sutmann,A. Skjellum,D. Kranzlmüller,P. V. Coveney

DOI: https://doi.org/10.48550/arXiv.2010.04144

2020-10-09

Abstract:Many scientific and medical researchers are working towards the creation of a virtual human - a personalised digital copy of an individual - that will assist in a patient's diagnosis, treatment and recovery. The complex nature of living systems means that the development of this remains a major challenge. We describe progress in enabling the HemeLB lattice Boltzmann code to simulate 3D macroscopic blood flow on a full human scale. Significant developments in memory management and load balancing allow near linear scaling performance of the code on hundreds of thousands of computer cores. Integral to the construction of a virtual human, we also outline the implementation of a self-coupling strategy for HemeLB. This allows simultaneous simulation of arterial and venous vascular trees based on human-specific geometries.

Medical Physics
A cardiac electromechanics model coupled with a lumped parameters model for closed-loop blood circulation. Part II: numerical approximation

Francesco Regazzoni,Matteo Salvador,Pasquale Claudio Africa,Marco Fedele,Luca Dede',Alfio Quarteroni

DOI: https://doi.org/10.48550/arXiv.2011.15051

2020-11-30

Numerical Analysis

Abstract:In the framework of accurate and efficient segregated schemes for 3D cardiac electromechanics and 0D cardiovascular models, we propose here a novel numerical approach to address the coupled 3D-0D problem introduced in Part I of this two-part series of papers. We combine implicit-explicit schemes to solve the different cardiac models in a multiphysics setting. We properly separate and manage the different time and space scales related to cardiac electromechanics and blood circulation. We employ a flexible and scalable intergrid transfer operator that enables to interpolate Finite Element functions among different meshes and, possibly, among different Finite Element spaces. We propose a numerical method to couple the 3D electromechanical model and the 0D circulation model in a numerically stable manner within a fully segregated fashion. No adaptations are required through the different phases of the heartbeat. We also propose a robust algorithm to reconstruct the stress-free reference configuration. Due to the computational cost associated with the numerical solution of this inverse problem, the reference configuration recovery algorithm comes along with a novel projection technique to precisely recover the unloaded geometry from a coarser representation of the computational domain. We show the convergence property of our numerical schemes by performing an accuracy study through grid refinement. To prove the biophysical accuracy of our computational model, we also address different scenarios of clinical interest in our numerical simulations by varying preload, afterload and contractility. Indeed, we simulate physiologically relevant behaviors and we reproduce meaningful results in the context of cardiac function.
Simulating Cardiovascular Haemodynamics By Bond Graph Technique

kaiyong dai,gangmin ning,shuzhen chen,juan du,xiaoxiang zheng

DOI: https://doi.org/10.1007/978-3-540-36841-0_46

2007-01-01

Abstract:A bond graph model for simulating cardiovascular haemodynamics is proposed. The model is given by an elastance variable description of heart, head and neck circulation, pulmonary circulation, coronary circulation, abdomen circulation and extremity circulation. Also, the carotid baroreflex control mechanism for circulation is embeded. The model responds to sympathetic and vagal activity by modifying systemic peripheral resistance, heart rate and ventricular end-systolic elastance. It is able to simulate the blood pressure and blood flow wave on arbitrary segment of circulation system, subsequently, the cardiovascular parameters such as stroke volume, cardiac output and eject fraction of heart can be computed. The simulation results under physiological condition are in agreement with human physiological experimental data. By modifying the arterial compliance and balance point of barorecpter, the hypertension is simulated and the results are consistent with the pathological features. In conclusion, the model could simulate the performance the cardiovascular system validly.
Cellular blood flow modelling with HemoCell

Gabor Zavodszky,Christian Spieker,Benjamin Czaja,Britt van Rooij

2023-05-04

Abstract:Many of the intriguing properties of blood originate from its cellular nature. Bulk effects, such as viscosity, depend on the local shear rates and on the size of the vessels. While empirical descriptions of bulk rheology are available for decades, their validity is limited to the experimental conditions they were observed under. These are typically artificial scenarios (e.g., perfectly straight glass tube or in pure shear with no gradients). Such conditions make experimental measurements simpler, however, they do not exist in real systems (i.e., in a real human circulatory system). Therefore, as we strive to increase our understanding on the cardiovascular system and improve the accuracy of our computational predictions, we need to incorporate a more comprehensive description of the cellular nature of blood. This, however, presents several computational challenges that can only be addressed by high performance computing. In this chapter we describe HemoCell , an open-source high performance cellular blood flow simulation, which implements validated mechanical models for red blood cells and is capable of reproducing the emergent transport characteristics of such a complex cellular system. We discuss the accuracy, the range of validity, and demonstrate applications on a series of human diseases.

Computational Engineering, Finance, and Science,Biological Physics
Hybrid models for simulating blood flow in microvascular networks

Ettore Vidotto,Timo Koch,Tobias Köppl,Rainer Helmig,Barbara Wohlmuth

DOI: https://doi.org/10.1137/18M1228712

2019-06-27

Abstract:In this paper, we are concerned with the simulation of blood flow in microvascular networks and the surrounding tissue. To reduce the computational complexity of this issue, the network structures are modeled by a one-dimensional graph, whose location in space is determined by the centerlines of the three-dimensional vessels. The surrounding tissue is considered as a homogeneous porous medium. Darcy's equation is used to simulate flow in the extra-vascular space, where the mass exchange with the blood vessels is accounted for by means of line source terms. However, this model reduction approach still causes high computational costs, in particular, when larger parts of an organ have to be simulated. This observation motivates the consideration of a further model reduction step. Thereby, we homogenize the fine scale structures of the microvascular networks resulting in a new hybrid approach modeling the fine scale structures as a heterogeneous porous medium and the flow in the larger vessels by one-dimensional flow equations. Both modeling approaches are compared with respect to mass fluxes and averaged pressures. The simulations have been performed on a microvascular network that has been extracted from the cortex of a rat brain.

Numerical Analysis,Fluid Dynamics
Development of a 3D vascular network visualization platform for one-dimensional hemodynamic simulation

Yan Chen,Masaharu Kobayashi,Changyoung Yuhn,Marie Oshima

DOI: https://doi.org/10.1101/2024.01.08.574602

2024-01-08

Abstract:Recent advancements in computational performance and medical simulation technology, have made significant strides, particularly in predictive diagnosis. This study focuses on the blood flow simulation reduced-order models, which provide swift and cost-effective solutions for complex vascular systems, positioning them as practical alternatives to 3D simulations in resource-limited medical settings. The paper introduces a visualization platform for patient-specific and image-based 1D-0D simulations. This platform covers the entire workflow, from modeling to dynamic 3D visualization of simulation results. Two case studies on respectively carotid stenosis and arterial remodeling demonstrate its utility in surgical planning and preventing complications like cerebral hyperperfusion syndrome.

Bioengineering
A Meshless Solver for Blood Flow Simulations in Elastic Vessels Using Physics-Informed Neural Network

Han Zhang,Raymond Chan,Xue-Cheng Tai

2024-06-06

Abstract:Investigating blood flow in the cardiovascular system is crucial for assessing cardiovascular health. Computational approaches offer some non-invasive alternatives to measure blood flow dynamics. Numerical simulations based on traditional methods such as finite-element and other numerical discretizations have been extensively studied and have yielded excellent results. However, adapting these methods to real-life simulations remains a complex task. In this paper, we propose a method that offers flexibility and can efficiently handle real-life simulations. We suggest utilizing the physics-informed neural network (PINN) to solve the Navier-Stokes equation in a deformable domain, specifically addressing the simulation of blood flow in elastic vessels. Our approach models blood flow using an incompressible, viscous Navier-Stokes equation in an Arbitrary Lagrangian-Eulerian form. The mechanical model for the vessel wall structure is formulated by an equation of Newton's second law of momentum and linear elasticity to the force exerted by the fluid flow. Our method is a mesh-free approach that eliminates the need for discretization and meshing of the computational domain. This makes it highly efficient in solving simulations involving complex geometries. Additionally, with the availability of well-developed open-source machine learning framework packages and parallel modules, our method can easily be accelerated through GPU computing and parallel computing. To evaluate our approach, we conducted experiments on regular cylinder vessels as well as vessels with plaque on their walls. We compared our results to a solution calculated by Finite Element Methods using a dense grid and small time steps, which we considered as the ground truth solution. We report the relative error and the time consumed to solve the problem, highlighting the advantages of our method.

Numerical Analysis,Fluid Dynamics
An Arbitrarily High-Order Fully Well-balanced Hybrid Finite Element-Finite Volume Method for a One-dimensional Blood Flow Model

Yongle Liu,Wasilij Barsukow

2024-04-28

Abstract:In this paper, we propose an arbitrarily high-order accurate fully well-balanced numerical method for the one-dimensional blood flow model. The developed method is based on a continuous representation of the solution and a natural combination of the conservative and primitive formulations of the studied PDEs. The degrees of freedom are defined as point values at cell interfaces and moments of the conservative variables inside the cell, drawing inspiration from the discontinuous Galerkin method. The well-balanced property, in the sense of an exact preservation of both the zero and non-zero velocity equilibria, is achieved by a well-balanced approximation of the source term in the conservative formulation and a well-balanced residual computation in the primitive formulation. To lowest (3rd) order this method reduces to the method developed in [Abgrall and Liu, A New Approach for Designing Well-Balanced Schemes for the Shallow Water Equations: A Combination of Conservative and Primitive Formulations, arXiv preprint,

Numerical Analysis
A monolithic 3D-0D coupled closed-loop model of the heart and the vascular system: Experiment-based parameter estimation for patient-specific cardiac mechanics

Marc Hirschvogel,Marina Bassilious,Lasse Jagschies,Stephen M Wildhirt,Michael W Gee

DOI: https://doi.org/10.1002/cnm.2842

Abstract:A model for patient-specific cardiac mechanics simulation is introduced, incorporating a 3-dimensional finite element model of the ventricular part of the heart, which is coupled to a reduced-order 0-dimensional closed-loop vascular system, heart valve, and atrial chamber model. The ventricles are modeled by a nonlinear orthotropic passive material law. The electrical activation is mimicked by a prescribed parameterized active stress acting along a generic muscle fiber orientation. Our activation function is constructed such that the start of ventricular contraction and relaxation as well as the active stress curve's slope are parameterized. The imaging-based patient-specific ventricular model is prestressed to low end-diastolic pressure to account for the imaged, stressed configuration. Visco-elastic Robin boundary conditions are applied to the heart base and the epicardium to account for the embedding surrounding. We treat the 3D solid-0D fluid interaction as a strongly coupled monolithic problem, which is consistently linearized with respect to 3D solid and 0D fluid model variables to allow for a Newton-type solution procedure. The resulting coupled linear system of equations is solved iteratively in every Newton step using 2 × 2 physics-based block preconditioning. Furthermore, we present novel efficient strategies for calibrating active contractile and vascular resistance parameters to experimental left ventricular pressure and stroke volume data gained in porcine experiments. Two exemplary states of cardiovascular condition are considered, namely, after application of vasodilatory beta blockers (BETA) and after injection of vasoconstrictive phenylephrine (PHEN). The parameter calibration to the specific individual and cardiovascular state at hand is performed using a 2-stage nonlinear multilevel method that uses a low-fidelity heart model to compute a parameter correction for the high-fidelity model optimization problem. We discuss 2 different low-fidelity model choices with respect to their ability to augment the parameter optimization. Because the periodic state conditions on the model (active stress, vascular pressures, and fluxes) are a priori unknown and also dependent on the parameters to be calibrated (and vice versa), we perform parameter calibration and periodic state condition estimation simultaneously. After a couple of heart beats, the calibration algorithm converges to a settled, periodic state because of conservation of blood volume within the closed-loop circulatory system. The proposed model and multilevel calibration method are cost-efficient and allow for an efficient determination of a patient-specific in silico heart model that reproduces physiological observations very well. Such an individual and state accurate model is an important predictive tool in intervention planning, assist device engineering and other medical applications.
A geometric multiscale model for the numerical simulation of blood flow in the human left heart

Alberto Zingaro,Ivan Fumagalli,Luca Dede,Marco Fedele,Pasquale C. Africa,Antonio F. Corno,Alfio Quarteroni

DOI: https://doi.org/10.3934/dcdss.2022052

2022-01-01

Abstract:We present a new computational model for the numerical simulation of blood flow in the human left heart. To this aim, we use the Navier-Stokes equations in an Arbitrary Lagrangian Eulerian formulation to account for the endocardium motion and we model the cardiac valves by means of the Resistive Immersed Implicit Surface method. To impose a physiological displacement of the domain boundary, we use a 3D cardiac electromechanical model of the left ventricle coupled to a lumped-parameter (0D) closed-loop model of the remaining circulation. We thus obtain a one-way coupled electromechanics-fluid dynamics model in the left ventricle. To extend the left ventricle motion to the endocardium of the left atrium and to that of the ascending aorta, we introduce a preprocessing procedure according to which an harmonic extension of the left ventricle displacement is combined with the motion of the left atrium based on the 0D model. To better match the 3D cardiac fluid flow with the external blood circulation, we couple the 3D Navier-Stokes equations to the 0D circulation model, obtaining a multiscale coupled 3D-0D fluid dynamics model that we solve via a segregated numerical scheme. We carry out numerical simulations for a healthy left heart and we validate our model by showing that meaningful hemodynamic indicators are correctly reproduced.

English Else
Blood flow distribution in an anatomically detailed arterial network model: criteria and algorithms

Pablo J. Blanco,Sansuke M. Watanabe,Enzo A. Dari,Marco Aurélio R. F. Passos,Raúl A. Feijóo

DOI: https://doi.org/10.1007/s10237-014-0574-8

IF: 3.5

2014-03-30

Biomechanics and Modeling in Mechanobiology

Abstract:Development of blood flow distribution criteria is a mandatory step toward developing computational models and numerical simulations of the systemic circulation. In the present work, we (i) present a systematic approach based on anatomical and physiological considerations to distribute the blood flow in a 1D anatomically detailed model of the arterial network and (ii) develop a numerical procedure to calibrate resistive parameters in terminal models in order to effectively satisfy such flow distribution. For the first goal, we merge data collected from the specialized medical literature with anatomical concepts such as vascular territories to determine blood flow supply to specific (encephalon, kidneys, etc.) and distributed (muscles, skin, etc.) organs. Overall, 28 entities representing the main specific organs are accounted for in the detailed description of the arterial topology that we use as model substrate. In turn, 116 vascular territories are considered as the basic blocks that compose the distributed organs throughout the whole body. For the second goal, Windkessel models are used to represent the peripheral beds, and the values of the resistive parameters are computed applying a Newton method to a parameter identification problem to guarantee the supply of the correct flow fraction to each terminal location according to the given criteria. Finally, it is shown that, by means of the criteria developed, and for a rather standard set of model parameters, the model predicts physiologically realistic pressure and flow waveforms.

engineering, biomedical,biophysics
Numerical simulation and fast method for the 0D‐1D multi‐scale coupled model and its application in ischemic brain tissue blood flow problems

Yi Liu,Junqing Jia,Fanhai Zeng,Xiaoyun Jiang

DOI: https://doi.org/10.1002/cnm.3828

IF: 2.1

2024-04-24

International Journal for Numerical Methods in Biomedical Engineering

Abstract:Abstract As living standards rise, more and more people are paying attention to their own health, especially issues such as cerebral thrombosis, cerebral infarction, and other cerebral blood flow problems. An accurate simulation of blood flow within cerebral vessels has emerged as a crucial area of research. In this study, we focus on microcirculatory blood flow in ischemic brain tissue and employ a 0D‐1D geometric multi‐scale coupled model to characterize this process. Given the intricate nature of human cerebral vessels, we apply a numerical method combining the finite element method and the third‐order Runge–Kutta method to resolve the coupled model. To enhance computational efficiency, we introduce a fast method based on the reduced‐order extrapolation algorithm. Our numerical example underscores the stability of the method and convergence accuracy to , while significantly improving the accuracy and efficiency of blood flow simulation, making the mechanism analysis more accurate. Additionally, we present examples detailing variations and distribution of intracranial pressure and blood flow in ischemic brain tissue throughout a cardiac cycle. Both reduced vascular compliance and vascular stenosis can have adverse effects on intracranial cerebral pressure and blood flow, leading to insufficient local oxygen supply and negative effects on brain function. Meanwhile, there will also be corresponding changes in volume flow and pulsatile blood pressure.

engineering, biomedical,mathematics, interdisciplinary applications,mathematical & computational biology
Object‐Oriented Lumped‐Parameter Modeling of the Cardiovascular System for Physiological and Pathophysiological Conditions

Luca Rosalia,Caglar Ozturk,David Van Story,Markus A. Horvath,Ellen T. Roche

DOI: https://doi.org/10.1002/adts.202000216

2021-02-11

Advanced Theory and Simulations

Abstract:In this work, a lumped‐parameter Windkessel model of the cardiovascular system that simulates biomechanical parameters of the human physiology is presented. The object‐oriented platform provided by the MATLAB‐based modeling environment SIMSCAPE is employed to compute blood pressures and flows in each heart chamber and at various sites of the vascular tree. The hydraulic domain allows the determination of cardiovascular hemodynamics intuitively from geometrical and mechanical properties of the system, while custom elements model the pumping action of the heart and the effects of respiration on blood flow. The model is validated by comparing predicted hemodynamics with normal physiology during both systole and diastole, demonstrating that changes in arterial pressures with breathing are consistent with reported physiological effects of cardiorespiratory coupling. The capabilities of this platform are explored through two exemplary case studies: i) pressure‐overload heart failure due to aortic constriction, validated in vitro and via finite element analysis, and ii) single‐ventricle Fontan physiology, validated in vitro and compared with the clinical literature. This platform provides a practical tool for the calculation of cardiovascular hemodynamics from hydraulic parameters, enabling the intuitive creation of in silico representations of complex circulatory loops, the planning and optimization of medical interventions, and the prediction of clinically relevant patient‐specific hemodynamics. A lumped‐parameter model of the cardiovascular system is developed to determine cardiovascular hemodynamics from hydraulic parameters in physiological and pathophysiological conditions, including left ventricular pressure overload and single‐ventricle physiology. This provides a practical tool for the simulation of complex circulatory loops, the optimization of medical interventions, and the prediction of patient‐specific hemodynamics of clinical relevance.

multidisciplinary sciences

Reply to "fixation and perception".

Physics-constrained coupled neural differential equations for one dimensional blood flow modeling

A one-dimensional computational model for blood flow in an elastic blood vessel with a rigid catheter

A modular framework for implicit 3D–0D coupling in cardiac mechanics

A Modular Framework for Implicit 3D-0D Coupling in Cardiac Mechanics

A multicore parallel algorithm for multiscale modelling of an entire human blood circulation network

A simplified particulate model for coarse-grained hemodynamics simulations

Simulation of unsteady blood flows in a patient-specific compliant pulmonary artery with a highly parallel monolithically coupled fluid-structure interaction algorithm

Towards Blood Flow in the Virtual Human: Efficient Self-Coupling of HemeLB

A cardiac electromechanics model coupled with a lumped parameters model for closed-loop blood circulation. Part II: numerical approximation

Simulating Cardiovascular Haemodynamics By Bond Graph Technique

Cellular blood flow modelling with HemoCell

Hybrid models for simulating blood flow in microvascular networks

Development of a 3D vascular network visualization platform for one-dimensional hemodynamic simulation

A Meshless Solver for Blood Flow Simulations in Elastic Vessels Using Physics-Informed Neural Network

An Arbitrarily High-Order Fully Well-balanced Hybrid Finite Element-Finite Volume Method for a One-dimensional Blood Flow Model

A monolithic 3D-0D coupled closed-loop model of the heart and the vascular system: Experiment-based parameter estimation for patient-specific cardiac mechanics

A geometric multiscale model for the numerical simulation of blood flow in the human left heart

Blood flow distribution in an anatomically detailed arterial network model: criteria and algorithms

Numerical simulation and fast method for the 0D‐1D multi‐scale coupled model and its application in ischemic brain tissue blood flow problems

Object‐Oriented Lumped‐Parameter Modeling of the Cardiovascular System for Physiological and Pathophysiological Conditions