Computational & Technology Resources
an online resource for computational,
engineering & technology publications
PROCEEDINGS OF THE ELEVENTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY
Edited by: B.H.V. Topping
Hierarchical Numerical Modelling of Nested Poroelastic Media
E. Rohan1 and R. Cimrman2
1Department of Mechanics, Faculty of Applied Sciences, 2New Technologies Research Centre,
E. Rohan, R. Cimrman, "Hierarchical Numerical Modelling of Nested Poroelastic Media", in B.H.V. Topping, (Editor), "Proceedings of the Eleventh International Conference on Computational Structures Technology", Civil-Comp Press, Stirlingshire, UK, Paper 250, 2012. doi:10.4203/ccp.99.250
Keywords: poroelasticity, homogenization, double porosity, micromechanics, multiscale modelling.
In nature as well as in technical practice one can often find materials with different levels of porosity at different scales. Considering poroelastic  description at the smallest scale, the theory of homogenization [3,8] provides a natural way of upscaling this description to higher levels of porosities in a sense that effective poroelastic material coefficients (consistent with the Biot model) at a higher level are obtained by applying homogenization to the lower level. This leads to a suitable hierarchical description of the porous medium, where different porosities associated with different scale levels are taken into account.
In this paper homogenization-based formulae are provided which enable the computation of the poroelasticity coefficients for a given geometry and topology at the micro- and meso-scopic levels. The homogenization at each scale level proceeds in two steps: 1) Find effective (homogenized) coefficients by solving auxiliary problems for several characteristic (or corrector) functions, cf. [7,5]. 2) Compute the homogenized coefficients that can be used for the higher level and, or "global" (homogenized) model of the current level. 3) Solve the homogenized problem with given loading and boundary conditions at the macro-scale (the highest level). As a result of the linearity of the problems, those steps are decoupled in a sense that the computation of the homogenized coefficients for the global level is valid for any point having the corresponding "microstructure".
The two-level upscaling "micro-meso-macro" is presented and the influence of the pore geometry and topology is illustrated. In  systems of interconnected pores at different scales were considered; in a steady state there is only one pore fluid pressure. In this paper, a different arrangement of porosities is described which can be mutually separated by a semipermeable interface, however, each porosity can form a separate connected system. In this situation, the homogenized problem results in two different pressures. At the mesoscopic scale we take into account the Darcy flow in the poroelastic matrix, although in the mesoscopic fractures (called channels in our terminology) the fluid is assumed to be static with no pressure gradients.
In the literature, homogenization in poroelasticity is a frequently discussed issue , but to our knowledge, the results reported in this paper are novel, namely due to the numerical feedback and computer implementation .
purchase the full-text of this paper (price £20)