Multiscale numerical modelling of perfusion in deformable double porous media described by the Biot-Darcy-Brinkman model¶
We consider a double porous medium which consists of an elastic solid matrix perforated by a system of channels filled with an incompressible fluid with interface . These components are arranged in a periodic lattice at both the micro- and mesoscopic level. Two small scale parameters and chracterize micro- and meso-porosities.
At the mesoscopic scale, the periodic structure is formed by two systems of fluid filled channels occupying domains , whereby and by domain which is constituted by a microporous material. In particular, domain represents micro pores saturated by fluid, whereas is the skeleton, see Fig. 1.
To summarize the decompositions,
The mesoscopic channels are mutually separated by the micro-porous matrix , so that also two disconnected interfaces are defined .
The microporous material is generated as a periodic lattice by repeating representative volume element (RVE) occupying domain . It splits in the solid and fluid part with interface , as seen on Fig. 1.
Similarly, at the mesoscopic level the porous structure is generated by RVE , with decomposition into the part filled by microporous matrix , and fluid channels . The channels further split into two sub-parts , each with interface , see Fig. 1.
The superscripts and denote the quantities oscillating within the heterogeneous structure with the period equal to the size of the micro- and mesoscopic periodic unit. However, we drop the superscript in following text to simplify the notation.
The mechanical behavior of such a structure can be described using the two-level asymptotic homogenization method, (for more detailed explanaition we refer to [RohanTurjanicovaLukes2020]).
The mechanical properties of the deformable matrix are given by elasticity tensor which satisfies the usual symmetries.
The fluid is characterized by viscosity which is given by a piece-wise constant function,
The scaling of the viscosity in micropores is the standart consequence of the non-slip boundary condition on he pore wall.
The problem of the fluid flow in deformable media at microscopic level is given by the following equilibrium equations and boundary conditions governing displacement of the solid and both the fluid pressure and velocity fields :
where is the fluid stress, is the strain in the solid with components , denotes the volume forces in the solid or in the fluid, and is the surfacetraction stresses acting on the solid part. The relative fluid velocity in the fluid-filled pores is defined whit use of a smooth extention of the dislacement field from solid to whole domain .
Due to the double-porous nature of the medium, we performe two levels of homogenization. The 1st-level of homogenization concerns the asymptotic analysis related to the fluid-structure interaction in microporous structure situated in . We apply the standard homogenization techniques to the above problem. It results in the limit model for , where is the scale parameter relating the microscopic and macroscopic length scales. The homogenization process leads to local microscopic problems, defined within a reference periodic cell , and to the mesoscopic problem describing the behavior of the homogenized matrix at the mesoscopic level. The mesoscopic problem involves the homogenized material coefficients which are evaluated using the solutions of the local problems. The 2nd-level of homogenization deals with upsacling from meso- to macroscopic scale. It results in the limit model for and subsequently in the local mesoscopic problems on a reference periodic cell , and in the derivation of the homogenized problem at macroscopic level. Due to linearity of the problem, the microscopic, mesoscopic and macroscopic problems are decoupled.
The local microscopic responses are given by the following sub-problems which are solved within the periodic reference cell , see Fig. 1, that is decomposed similarly to the decomposition of domain :
- Find , such that for all for any
- Find , such that for all , satisfying
The microscopic sub-problems are solved with the periodic boundary conditions and is the interface between solid and fluid part of the cell .
With the characteristic responses obtained by solving local sub-problems, the homogenized material coefficients , , and can be evaluated using the following expressions:
Homogenization - 2nd level¶
At the 2st-level of homogenization, the asymptotic analysis is related to the interaction between the homogenized microporous matrix in and fluid in channels at mesoscopic level. By similar upscaling procedure as in 1st level of homogenization, we obtain local mesoscopic problems, defined within a reference periodic cell , where we enter the homogenized coefficients obtained by 1st level homogenization. However, when compared to the model of a mesoscopic structure containing only one system of channels, the upscaled mesoscopic structures involving two mesoscopic channels yields two macroscopic velocity fields which describe the two parallel flows.
We also arrive to the global problem describing the behavior of the homogenized matrix at the macroscopic level. The homogenized material coefficients describing whichdescribe behavior at macroscopic level are evaluated using the solutions of the local mesoscopic problems.Due to linearity of the problem, the microscopic, mesoscopic and macroscopic problems are decoupled.
The local mesoscopic responses are given by the following sub-problems which are solved within the periodic reference cell , see Fig. 1, that is decomposed similarly to the decomposition of domain :
- Find , such that for all for any
- Find , such that for all satisfying
- Find , such that for all , satisfying
The mesoscopic sub-problems are solved with the periodic boundary conditions and is the interface between the matrix part and system of channels .
With the characteristic responses obtained by solving local sub-problems, the homogenized material coefficients , , , , , and can be evaluated using the following expressions:
All coefficients are symmetric with respect to indices related to strain and strain rate tensors, i.e. .
The global macroscopic problem is defined in terms of the homogenized coefficients as: Find the macroscopic displacements , pressure and velocity fields such that for all , and
The Dirichlet boundary conditions prescribed for and can be imposed,
The pressure fulfils zero-means conditions in the whole domain .
The complementary Neumann-type boundary condition are specfied as follows
For the purpose of this example, we simplify the problem (13). We omitt all volume froces and and surface tractions . Also we will consider closed micropores on the whole boundary , i.e. . The problem (10) becomes: Find the macroscopic displacements , velocity and pressure fields such that for all , and
The initial conditions of fields and are necessary for computation of time dependent problem and are computed using the steady state form of the boundary problem defined above.
To run the numerical simulation, download the archive, unpack it in the main SfePy directory and type:
This invoke the
simply.py script which calculates the macroscopic
problem (13) and calls the homogenization engine that solves the
local subproblems for given parameter , viscosity and elastic tensor . First, it solves subproblems (3) and (4)
on microscopic cell (see Fig. 2 left) and
evaluates the homogenized
coefficients (5). Then, using solution from previous step, it solves subproblems
(6)–(8) on mesoscopic cell (see Fig. 2 right)
and evaluates the homogenized coefficients (9). See [CimrmanLukesRohan2019] for more details related
to the SfePy homogenization engine.
Then, the script computes both the steady state and time evolution of the problem (13). For steady state, the macroscopic sample is fixed on both top and bottom side, so that no displacements are allowed, see Fig. 3. The defromation is induced due to the flow through porous matrix, as the responce to the prescribed velocities , see (13).
For simulation of the time evolution of the macroscopic problem we take the steady state as initial value of displacement and preassure fields. In this case, the gradually increasing displacement is prescribed on the top side of the macroscopic sample.
The resulting macroscopic pressure field , displacement and the velocity fields are depicted in Fig. 4, Fig. 5 and Fig. 4. These figures depict distribution of macroscopic quantities on the deformingin at macroscopic specimen at different computational times (, and ). The nondeformed shape of macroscopic specimen is visualised by its outline.
|[RohanTurjanicovaLukes2020]||Rohan E., Turjanicová J., Lukeš V. Multiscale modelling and simulations of tissue perfusion using the Biot-Darcy-Brinkman model. Submitted to Computers & Structures, 2020|
|[CimrmanLukesRohan2019]||Cimrman R., Lukes V., Rohan E. Multiscale finite element calculations in Python using SfePy. Advances in Computational Mathematics, 45(4):1897-1921, 2019, DOI:10.1007/s10444-019-09666-0|