Mathematical analysis and numerical approximation of a general linearized poro-hyperelastic model

Published on Jun 1, 2021in Computers & Mathematics With Applications3.37
· DOI :10.1016/J.CAMWA.2020.07.025
Nicolás Barnafi1
Estimated H-index: 1
(Polytechnic University of Milan),
Nicolas A. Barnafi2
Estimated H-index: 2
(Polytechnic University of Milan)
+ 1 AuthorsAlfio Quarteroni83
Estimated H-index: 83
(Polytechnic University of Milan)
Sources
Abstract
Abstract We describe the behavior of a deformable porous material by means of a poro-hyperelastic model that has been previously proposed in Chapelle and Moireau (2014) under general assumptions for mass and momentum balance and isothermal conditions for a two-component mixture of fluid and solid phases. In particular, we address here a linearized version of the model, based on the assumption of small displacements. We consider the mathematical analysis and the numerical approximation of the problem. More precisely, we carry out firstly the well-posedness analysis of the model. Then, we propose a numerical discretization scheme based on finite differences in time and finite elements for the spatial approximation; stability and numerical error estimates are proved. Particular attention is dedicated to the study of the saddle-point structure of the problem, that turns out to be interesting because velocities of the fluid phase and of the solid phase are combined into a single quasi-incompressibility constraint. Our analysis provides guidelines to select the componentwise polynomial degree of approximation of fluid velocity, solid displacement and pressure, to obtain a stable and robust discretization based on Taylor–Hood type finite element spaces. Interestingly, we show how this choice depends on the porosity of the mixture, i.e. the volume fraction of the fluid phase.
References38
Newest
#1Rana Zakerzadeh (University of Texas at Austin)H-Index: 6
#2Paolo Zunino (Polytechnic University of Milan)H-Index: 28
We study the effect of poroelasticity on fluid–structure interaction. More precisely, we analyze the role of fluid flow through a deformable porous matrix in the energy dissipation behavior of a poroelastic structure. For this purpose, we develop and use a nonlinear poroelastic computational model and apply it to the fluid–structure interaction simulations. We discretize the problem by means of the finite element method for the spatial approximation and using finite differences in time. The nume...
5 CitationsSource
#1Bruno Burtschell (Université Paris-Saclay)H-Index: 2
#2Philippe Moireau (Université Paris-Saclay)H-Index: 21
Last. Dominique Chapelle (Université Paris-Saclay)H-Index: 5
view all 3 authors...
We consider a previously proposed general nonlinear poromechanical formulation, and we derive a linearized version of this model. For this linearized model, we obtain an existence result and we propose a complete discretization strategy - in time and space - with a special concern for issues associated with incompressible or nearly-incompressible behavior. We provide a detailed mathematical analysis of this strategy, the main result being an error estimate uniform with respect to the compressibi...
2 CitationsSource
#1Ilona Ambartsumyan (University of Pittsburgh)H-Index: 6
#2Eldar Khattatov (University of Pittsburgh)H-Index: 5
Last. Paolo Zunino (Polytechnic University of Milan)H-Index: 28
view all 4 authors...
We study a finite element computational model for solving the coupled problem arising in the interaction between a free fluid and a fluid in a poroelastic medium. The free fluid is governed by the Stokes equations, while the flow in the poroelastic medium is modeled using the Biot poroelasticity system. Equilibrium and kinematic conditions are imposed on the interface. A mixed Darcy formulation is employed, resulting in continuity of flux condition of essential type. A Lagrange multiplier method...
32 CitationsSource
#1Dietmar Gallistl (UT: University of Twente)H-Index: 16
4 CitationsSource
#1Bruno Burtschell (University of Paris)H-Index: 2
#2Dominique Chapelle (University of Paris)H-Index: 5
Last. Philippe Moireau (University of Paris)H-Index: 21
view all 3 authors...
Proposed time discretization scheme for general two-phase poromechanical model.Discrete energy estimate with the total free energy of the mixture, in a general nonlinear framework.Numerical examples with representative test problems. We consider a general nonlinear poromechanical model, formulated based on fundamental thermodynamics principle, suitable for representing the coupling of rapid internal fluid flows with large deformations of the solid, and compatible with a wide class of constitutiv...
5 CitationsSource
#1Michelle Hine Armstrong (UA: University of Arizona)H-Index: 2
#2Adrian Buganza Tepole (Stanford University)H-Index: 16
Last. Jonathan P. Vande GeestH-Index: 22
view all 5 authors...
The purpose of this manuscript is to establish a unified theory of porohyperelasticity with transport and growth and to demonstrate the capability of this theory using a finite element model developed in MATLAB. We combine the theories of volumetric growth and mixed porohyperelasticity with transport and swelling (MPHETS) to derive a new method that models growth of biological soft tissues. The conservation equations and constitutive equations are developed for both solid-only growth and solid/f...
12 CitationsSource
#1Radomir Chabiniok ('KCL': King's College London)H-Index: 16
#2Vicky Y. Wang (University of Auckland)H-Index: 11
Last. David Nordsletten ('KCL': King's College London)H-Index: 25
view all 12 authors...
With heart and cardiovascular diseases continually challenging healthcare systems worldwide, translating basic research on cardiac (patho)physiology into clinical care is essential. Exacerbating this already extensive challenge is the complexity of the heart, relying on its hierarchical structure and function to maintain cardiovascular flow. Computational modelling has been proposed and actively pursued as a tool for accelerating research and translation. Allowing exploration of the relationship...
106 CitationsSource
#1Martin Costabel (University of Rennes)H-Index: 42
#2Michel Crouzeix (University of Rennes)H-Index: 20
Last. Yvon Lafranche (University of Rennes)H-Index: 5
view all 4 authors...
The inf-sup constant for the divergence, or LBB constant, is related to the Cosserat spectrum. It has been known for a long time that on non-smooth domains the Cosserat operator has a non-trivial essential spectrum, which can be used to bound the LBB constant from above. We prove that the essential spectrum on a plane polygon consists of an interval related to the corner angles and that on three-dimensional domains with edges, the essential spectrum contains such an interval. We obtain some nume...
10 CitationsSource
#1Alain Goriely (University of Oxford)H-Index: 58
#2Mgd Marc Geers (TU/e: Eindhoven University of Technology)H-Index: 65
Last. Ellen Kuhl (Stanford University)H-Index: 71
view all 10 authors...
The human brain is the continuous subject of extensive investigation aimed at understanding its behavior and function. Despite a clear evidence that mechanical factors play an important role in regulating brain activity, current research efforts focus mainly on the biochemical or electrophysiological activity of the brain. Here, we show that classical mechanical concepts including deformations, stretch, strain, strain rate, pressure, and stress play a crucial role in modulating both brain form a...
190 CitationsSource
#1Anh-Tu Vuong (TUM: Technische Universität München)H-Index: 5
#2Lena Yoshihara (TUM: Technische Universität München)H-Index: 7
Last. Wolfgang A. Wall (TUM: Technische Universität München)H-Index: 62
view all 3 authors...
Abstract In the last few decades modeling deformation and flow in porous media has been of great interest due to its possible application areas in various fields of engineering such as biomechanics, soil mechanics, geophysics, physical chemistry and material sciences. Due to the high complexity and in most cases also unknown geometry of porous media on the microscale, a fully resolved model is nearly impossible to obtain, but most of the times also not necessary to answer important questions. As...
35 CitationsSource
Cited By2
Newest
Source
#1Jakub Wiktor BothH-Index: 5
#2Nicolas A. BarnafiH-Index: 2
Last. Alfio QuarteroniH-Index: 83
view all 5 authors...
We address numerical solvers for a poromechanics model particularly adapted for soft materials, as it generally respects thermodynamics principles and energy balance. Considering the multi-physics nature of the problem, which involves solid and fluid species, interacting on the basis of mass balance and momentum conservation, we decide to adopt a solution strategy of the discrete problem based on iterative splitting schemes. As the model is similar (but not equivalent to) the Biot poromechanics ...