An official website of the United States government
In biomechanics, local phenomena, such as tissue perfusion, are strictly related to the global features of the surrounding blood circulation. In this paper, we propose a heterogeneous model where a local, accurate, 3D description of tissue perfusion by means of fluid flows through deformable porous media equations is coupled with a systemic, 0D, lumped model of the remainder of the circulation, where the fluid flow through a vascular network is described via its analog with a current flowing through an electric circuit. This represents a multiscale strategy, which couples an initial boundary value problem to be used in a specific tissue region with an initial value problem in the surrounding circulatory system. This PDE/ODE coupling leads to interface conditions enforcing the continuity of mass and the balance of stresses across models at different scales, and careful consideration is taken to address this interface mismatch. The resulting system involves PDEs of mixed type with interface conditions depending on nonlinear ODEs. A new result on local existence of solutions for this multiscale interface coupling is provided in this article.
more »
« less
This content will become publicly available on February 1, 2026
Comparing Interface Conditions for a 3D–0D Multiscale Interface Coupling With Applications in Tissue Perfusion
Many pathologies are related to hemodynamic changes occurring at the microvascular level, where small vessels pierce the tissue, perfusing it with blood. Since there is a large number of vessels of small caliber, it is impractical to model the fluid flow through each one of them separately, as it is done in the case of large arteries using, for example, the Navier–Stokes equations. As an alternative, tissue perfusion is modeled here via three‐dimensional (3D) partial differential equations (PDEs) for fluid flow through deformable porous media, where blood vessels are modeled as pores within a deformable solid representing the tissue. Since it is known that the local perfusion is related to the systemic features of surrounding blood circulation, we couple the PDE system with a zero‐dimensional (0D) lumped circuit model, obtained by the analogy between fluid flows in hydraulic networks and current flowing in electrical circuits. An important feature in this multiscale 3D–0D coupling is the specification of interface conditions between the 3D and the 0D parts of the system. In this article, we focus on two types of interface conditions driven by physical considerations, and compare the behavior of the solutions for the two different scenarios.
more »
« less
- PAR ID:
- 10639179
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Biomedical Engineering
- Volume:
- 41
- Issue:
- 2
- ISSN:
- 2040-7939
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Numerous disease conditions involve the sudden or progressive loss of blood flow. Perfusion restoration is vital for returning affected organs to full health. While a range of clinical interventions can successfully restore flow to downstream tissues, the microvascular responses after a loss-of-flow event can vary over time and may involve substantial microvessel instability. Increased insight into perfusion-mediated capillary stability and access-to-flow is therefore essential for advancing therapeutic reperfusion strategies and improving patient outcomes. To that end, we developed a tissue-based microvascular fluidics model to better understand (i) microvascular stability and access-to-flow over an acute time course post-ischemia, and (ii) collateral flow in vessels neighboring an occlusion site. We utilized murine intestinal tissue regions by catheterizing a feeder artery and introducing perfusate at physiologically comparable flow-rates. The cannulated vessel as well as a portion of the downstream vessels and associated intestinal tissue were cultured while constant perfusion conditions were maintained. An occlusion was introduced in a selected arterial segment, and changes in perfusion within areas receiving varying degrees of collateral flow were observed over time. To observe the microvascular response to perfusion changes, we incorporated (i) tissues harboring cell-reporter constructs, specifically Ng2-DsRed labeling of intestinal pericytes, and (ii) different types of fluorescent perfusates to quantify capillary access-to-flow at discrete time points. In our model, we found that perfusion tracers could enter capillaries within regions downstream of an occlusion upon the initial introduction of perfusion, but at 24 h tissue perfusion was severely decreased. However, live/dead cell discrimination revealed that the tissue overall did not experience significant cell death, including that of microvascular pericytes, even after 48 h. Our findings suggest that altered flow conditions may rapidly initiate cellular responses that reduce capillary access-to-flow, even in the absence of cellular deterioration or hypoxia. Overall, this ex vivo tissue-based microfluidics model may serve as a platform upon which a variety of follow-on studies may be conducted. It will thus enhance our understanding of microvessel stability and access-to-flow during an occlusive event and the role of collateral flow during normal and disrupted perfusion.more » « less
-
Microfluidic devices offer well-defined physical environments that are suitable for effective cell seeding and in vitro three-dimensional (3D) cell culture experiments. These platforms have been employed to model in vivo conditions for studying mechanical forces, cell–extracellular matrix (ECM) interactions, and to elucidate transport mechanisms in 3D tissue-like structures, such as tumor and lymph node organoids. Studies have shown that fluid flow behavior in microfluidic slides (µ-slides) directly influences shear stress, which has emerged as a key factor affecting cell proliferation and differentiation. This study investigates fluid flow in the porous channel of a µ-slide using computational fluid dynamics (CFD) techniques to analyze the impact of perfusion flow rate and porous properties on resulting shear stresses. The model of the µ-slide filled with a permeable biomaterial is considered. Porous media fluid flow in the channel is characterized by adding a momentum loss term to the standard Navier–Stokes equations, with a physiological range of permeability values. Numerical simulations are conducted to obtain data and contour plots of the filtration velocity and flow-induced shear stress distributions within the device channel. The filtration flow is subsequently measured by performing protein perfusions into the slide embedded with native human-derived ECM, while the flow rate is controlled using a syringe pump. The relationships between inlet flow rate and shear stress, as well as filtration flow and ECM permeability, are analyzed. The findings provide insights into the impact of shear stress, informing the optimization of perfusion conditions for studying tissues and cells under fluid flow.more » « less
-
Abstract We consider the finite element approximation of a coupled fluid‐structure interaction (FSI) system, which comprises a three‐dimensional (3D) Stokes flow and a two‐dimensional (2D) fourth‐order Euler–Bernoulli or Kirchhoff plate. The interaction of these parabolic and hyperbolic partial differential equations (PDE) occurs at the boundary interface which is assumed to be fixed. The vertical displacement of the plate dynamics evolves on the flat portion of the boundary where the coupling conditions are implemented via the matching velocities of the plate and fluid flow, as well as the Dirichlet boundary trace of the pressure. This pressure term also acts as a coupling agent, since it appears as a forcing term on the flat, elastic plate domain. Our main focus in this work is to generate some numerical results concerning the approximate solutions to the FSI model. For this, we propose a numerical algorithm that sequentially solves the fluid and plate subsystems through an effective decoupling approach. Numerical results of test problems are presented to illustrate the performance of the proposed method.more » « less
-
Abstract. We consider a nonlinear, moving boundary, fluid-structure interaction problem between a time dependent incompressible, viscous fluid flow, and an elastic structure composed of a cylindrical shell supported by a mesh of elastic rods. The fluid flow is modeled by the time-dependent Navier- Stokes equations in a three-dimensional cylindrical domain, while the lateral wall of the cylinder is modeled by the two-dimensional linearly elastic Koiter shell equations coupled to a one-dimensional system of conservation laws defined on a graph domain, describing a mesh of curved rods. The mesh supported shell allows displacements in all three spatial directions. Two-way coupling based on kinematic and dynamic coupling conditions is assumed between the fluid and composite structure, and between the mesh of curved rods and Koiter shell. Problems of this type arise in many ap- plications, including blood flow through arteries treated with vascular prostheses called stents. We prove the existence of a weak solution to this nonlinear, moving boundary problem by using the time discretization via Lie operator splitting method combined with an Arbitrary Lagrangian-Eulerian approach, and a non-trivial extension of the Aubin-Lions-Simon compactness result to problems on moving domains.more » « less
