Discontinuous Galerkin method for a three-dimensional coupled fluid-poroelastic model with applications to brain fluid mechanics
The modeling of the interaction between a poroelastic medium and a fluid in a hollow cavity is crucial for understanding, e.g., the multiphysics flow of blood and Cerebrospinal Fluid (CSF) in the brain, the supply of blood by the coronary arteries in heart perfusion, or the interaction between groun...
Saved in:
| Main Author: | |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
AIMS Press
2025-04-01
|
| Series: | Mathematics in Engineering |
| Subjects: | |
| Online Access: | https://www.aimspress.com/article/doi/10.3934/mine.2025006 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1850125182705860608 |
|---|---|
| author | Ivan Fumagalli |
| author_facet | Ivan Fumagalli |
| author_sort | Ivan Fumagalli |
| collection | DOAJ |
| description | The modeling of the interaction between a poroelastic medium and a fluid in a hollow cavity is crucial for understanding, e.g., the multiphysics flow of blood and Cerebrospinal Fluid (CSF) in the brain, the supply of blood by the coronary arteries in heart perfusion, or the interaction between groundwater and rivers or lakes. In particular, the cerebral tissue's elasticity and its perfusion by blood and interstitial CSF can be described by Multi-compartment Poroelasticity (MPE), while CSF flow in the brain ventricles can be modeled by the (Navier-)Stokes equations, the overall system resulting in a coupled MPE-(Navier-)Stokes system. The aim of this paper is three-fold. First, we aim to extend a recently presented discontinuous Galerkin method on polytopal grids (PolyDG) to incorporate three-dimensional geometries and physiological interface conditions. Regarding the latter, we consider here the Beavers-Joseph-Saffman (BJS) conditions at the interface: These conditions are essential to model the friction between the fluid and the porous medium. Second, we quantitatively analyze the computational efficiency of the proposed method on a domain with small geometrical features, thus demonstrating the advantages of employing polyhedral meshes. Finally, by a comparative numerical investigation, we assess the fluid-dynamics effects of the BJS conditions and of employing either Stokes or Navier-Stokes equations to model the CSF flow. The semidiscrete numerical scheme for the coupled problem is proved to be stable and optimally convergent. Temporal discretization is obtained using Newmark's $ \beta $-method for the elastodynamics equation and the $ \theta $-method for the remaining equations of the model. The theoretical error estimates are verified by numerical simulations on a test case with a manufactured solution, and a numerical investigation is carried out on a three-dimensional geometry to assess the effects of interface conditions and fluid inertia on the system. |
| format | Article |
| id | doaj-art-28c7bd7a49c8421699e5a96bb3bae1cc |
| institution | OA Journals |
| issn | 2640-3501 |
| language | English |
| publishDate | 2025-04-01 |
| publisher | AIMS Press |
| record_format | Article |
| series | Mathematics in Engineering |
| spelling | doaj-art-28c7bd7a49c8421699e5a96bb3bae1cc2025-08-20T02:34:10ZengAIMS PressMathematics in Engineering2640-35012025-04-017213016110.3934/mine.2025006Discontinuous Galerkin method for a three-dimensional coupled fluid-poroelastic model with applications to brain fluid mechanicsIvan Fumagalli0MOX, Dipartimento di Matematica, Politecnico di Milano, piazza Leonardo da Vinci 32, Milan 20133, ItalyThe modeling of the interaction between a poroelastic medium and a fluid in a hollow cavity is crucial for understanding, e.g., the multiphysics flow of blood and Cerebrospinal Fluid (CSF) in the brain, the supply of blood by the coronary arteries in heart perfusion, or the interaction between groundwater and rivers or lakes. In particular, the cerebral tissue's elasticity and its perfusion by blood and interstitial CSF can be described by Multi-compartment Poroelasticity (MPE), while CSF flow in the brain ventricles can be modeled by the (Navier-)Stokes equations, the overall system resulting in a coupled MPE-(Navier-)Stokes system. The aim of this paper is three-fold. First, we aim to extend a recently presented discontinuous Galerkin method on polytopal grids (PolyDG) to incorporate three-dimensional geometries and physiological interface conditions. Regarding the latter, we consider here the Beavers-Joseph-Saffman (BJS) conditions at the interface: These conditions are essential to model the friction between the fluid and the porous medium. Second, we quantitatively analyze the computational efficiency of the proposed method on a domain with small geometrical features, thus demonstrating the advantages of employing polyhedral meshes. Finally, by a comparative numerical investigation, we assess the fluid-dynamics effects of the BJS conditions and of employing either Stokes or Navier-Stokes equations to model the CSF flow. The semidiscrete numerical scheme for the coupled problem is proved to be stable and optimally convergent. Temporal discretization is obtained using Newmark's $ \beta $-method for the elastodynamics equation and the $ \theta $-method for the remaining equations of the model. The theoretical error estimates are verified by numerical simulations on a test case with a manufactured solution, and a numerical investigation is carried out on a three-dimensional geometry to assess the effects of interface conditions and fluid inertia on the system.https://www.aimspress.com/article/doi/10.3934/mine.2025006navier-stokes equationsmultiple-network poroelasticity theorybeavers-joseph-saffman interface conditionspolyhedral meshcerebrospinal fluid |
| spellingShingle | Ivan Fumagalli Discontinuous Galerkin method for a three-dimensional coupled fluid-poroelastic model with applications to brain fluid mechanics Mathematics in Engineering navier-stokes equations multiple-network poroelasticity theory beavers-joseph-saffman interface conditions polyhedral mesh cerebrospinal fluid |
| title | Discontinuous Galerkin method for a three-dimensional coupled fluid-poroelastic model with applications to brain fluid mechanics |
| title_full | Discontinuous Galerkin method for a three-dimensional coupled fluid-poroelastic model with applications to brain fluid mechanics |
| title_fullStr | Discontinuous Galerkin method for a three-dimensional coupled fluid-poroelastic model with applications to brain fluid mechanics |
| title_full_unstemmed | Discontinuous Galerkin method for a three-dimensional coupled fluid-poroelastic model with applications to brain fluid mechanics |
| title_short | Discontinuous Galerkin method for a three-dimensional coupled fluid-poroelastic model with applications to brain fluid mechanics |
| title_sort | discontinuous galerkin method for a three dimensional coupled fluid poroelastic model with applications to brain fluid mechanics |
| topic | navier-stokes equations multiple-network poroelasticity theory beavers-joseph-saffman interface conditions polyhedral mesh cerebrospinal fluid |
| url | https://www.aimspress.com/article/doi/10.3934/mine.2025006 |
| work_keys_str_mv | AT ivanfumagalli discontinuousgalerkinmethodforathreedimensionalcoupledfluidporoelasticmodelwithapplicationstobrainfluidmechanics |