Abstract
A new mathematical model for the macroscopic behavior of a material composed of a poroelastic solid embedding a Newtonian fluid network phase (also referred to as vascularized poroelastic material), with fluid transport between them, is derived via asymptotic homogenization. The typical distance between the vessels/channels (microscale) is much smaller than the average size of a whole domain (macroscale). The homogeneous and isotropic Biot’s equation (in the quasi-static case and in absence of volume forces) for the poroelastic phase and the Stokes’ problem for the fluid network are coupled through a fluid-structure interaction problem which accounts for fluid transport between the two phases; the latter is driven by the pressure difference between the two compartments. The averaging process results in a new system of partial differential equations that formally reads as a double poroelastic, globally mass conserving, model, together with a new constitutive relationship for the whole material which encodes the role of both pore and fluid network pressures. The mathematical model describes the mutual interplay among fluid filling the pores, flow in the network, transport between compartments, and linear elastic deformation of the (potentially compressible) elastic matrix comprising the poroelastic phase. Assuming periodicity at the microscale level, the model is computationally feasible, as it holds on the macroscale only (where the microstructure is smoothed out), and encodes geometrical information on the microvessels in its coefficients, which are to be computed solving classical periodic cell problems. Recently developed double porosity models are recovered when deformations of the elastic matrix are neglected. The new model is relevant to a wide range of applications, such as fluid in porous, fractured rocks, blood transport in vascularized, deformable tumors, and interactions across different hierarchical levels of porosity in the bone.
![](https://cdn.statically.io/img/media.springernature.com/m312/springer-static/image/art%3A10.1007%2Fs11012-017-0625-1/MediaObjects/11012_2017_625_Fig1_HTML.gif)
Similar content being viewed by others
Notes
These two paper differ in some scaling assumptions concerning the drug transport analysis, but the double porosity fluid transport models, for a purely Newtonian blood vessels’ rheology and for macroscopically uniform media, coincide, as actually demonstrated in [35].
The vector field \(\varvec{w}_p\) that satisfies Darcy’s law (8) actually represents the average relative fluid velocity in the porous medium and is also referred to as the discharge, flux, or filtration velocity. It is indeed related to the specific relative velocity \(\varvec{v}_p-\dot{\varvec{u}}\) via the interstitial porosity \(\phi\). However, since we conveniently carry out our analysis in terms of \(\varvec{w}_p\) hereinafter, we simply refer to \(\varvec{w}_p\) as the relative interstitial velocity.
Our admissibility constraint (20) is equivalent to equation (2.a), page 7, [43], setting their \(\beta =0\) and considering that \(c_0\), \(\varvec{v}^1\), \(\varvec{v}^2\), and \(\varvec{q}\) are denoted by \(\tilde{\alpha }\), \(\dot{\varvec{u}}\), \(\varvec{u}_n\), \(\varvec{w}_p\) in our manuscript.
Note that in the appendix reported in [13] the third rank tensor \(v^{ij}_p\) is equivalent to our \(\mathbb {{\mathcal {A}}}\) and its gradient to our \({\mathbb {M}}\). However, the latter is always identified to \({\mathbb {LC}}\), although relationship (138) actually rigorously holds only when \({\mathbb {C}}\) is locally constant.
See Section 4.3.1, pages 24–25, [35], and the model (4.100-4.103) therein, which is in turn proved to be equivalent to Eqs. (54), (56), (70), (75) in [42]. In [35], the quantities \({\tilde{{\mathsf {K}}}}\), \({\tilde{{\mathsf {G}}}}\), \(p_p\), \(\varvec{v}_n\), \(\phi \varvec{v}_p\), \(|\varOmega _p|\), \(|\varGamma |\) are denoted by \({\mathsf {K}}\), \({\mathsf {E}}\), \(p_t\), \(\varvec{u}_n\), \(\varvec{u}_t\), \(|\varOmega _t|\), S, respectively. Both in [42] and [35], only specific averages are used, and they are denoted as we denote non-specific averages, i.e. \(\displaystyle \left\langle \,\bullet \,\right\rangle _k\).
References
Auriault JL, Boutin C, Geindreau C (2010) Homogenization of coupled phenomena in heterogenous media, vol 149. Wiley, Hoboken
Bai M, Meng F, Elsworth D, Abousleiman Y, Roegiers JC (1999) Numerical modelling of coupled flow and deformation in fractured rock specimens. Int J Numer Anal Methods Geomech 23(2):141–160
Bakhvalov N, Panasenko G (1989) Homogenisation averaging processes in periodic media. Springer, Berlin
Beavers GS, Joseph DD (1967) Boundary conditions at a naturally permeable wall. J Fluid Mech 30:197–207
Berryman JG (2005) Comparison of upscaling methods in poroelasticity and its generalizations. J Eng Mech 131(9):928–936
Biot MA (1941) General theory of three-dimensional consolidation. J Appl Phys 12(2):155–164
Biot MA (1955) Theory of elasticity and consolidation for a porous anisotropic solid. J Appl Phys 26(2):182–185
Biot M (1956) General solutions of the equations of elasticity and consolidation for a porous material. J Appl Mech 23(1):91–96
Biot MA (1956) Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low-frequency range. J Acoust Soc Am 28(2):168–178
Biot MA (1962) Mechanics of deformation and acoustic propagation in porous media. J Appl Phys 33(4):1482–1498
Bottaro A, Ansaldi T (2012) On the infusion of a therapeutic agent into a solid tumor modeled as a poroelastic medium. J Biomech Eng 134:1–6
Bruna M, Chapman SJ (2015) Diffusion in spatial varying porous media. SIAM J Appl Math 75(4):1648–1674
Burridge R, Keller J (1981) Poroelasticity equations derived from microstructure. J Acoust Soc Am 70:1140–1146
Coussy O (2011) Mechanics and physics of porous solids. Wiley, Hoboken
Cowin SC (1999) Bone poroelasticity. J Biomech 32(3):217–238
Dalwadi MP, Griffiths IM, Bruna M (2015) Understanding how porosity gradients can make a better filter using homogenization theory. In: Proceedings of Royal Society A, vol 471, p 20150464
Detournay E, Cheng AD (2004) Fundamentals of poroelasticity. In: Fairhurst C (ed) Comprehensive rock engineering: principles, practice and projects, analysis and design method, vol II. Pergamon Press, Oxford, pp 113–171
Discacciati M, Quarteroni A (2009) Navier-stokes/darcy coupling: modeling, analysis, and numerical approximation. Rev Matem Complut 22(2):315–426
Eshelby J (1957) The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proc R Soc Lond Ser A Math Phys Sci 241:376–396
Gurtin M, Fried E, Anand L (2010) The mechanics and thermodynamics of continua. Cambridge University Press, Cambridge
Holmes M (1995) Introduction to perturbation method. Springer, Berlin
Hornung U (1997) Homogenization and porous media. Springer, Berlin
Jain RK (1987) Transport of molecules across tumor vasculature. Cancer Metastasis Rev 6:559–594
Jain RK (1990) Physiological barriers to delivery of monoclonal antibodies and other macromolecules in tumors. Cancer Res 50:814–819
Jain RK, Baxter LT (1988) Mechanisms of heterogeneous distribution of monoclonal antibodies and other macromolecules in tumors: significance of elevated interstitial pressure. Cancer Res 48:7022–7032
Jones IP (1973) Low reynolds number flow past a porous spherical shell. Math Proc Camb Philos Soc 73:231–238
Kedem O, Katchalsky A (1958) Thermodynamic of the permeability of biological membranes to non-electrolytes. Biochim Biophys Acta 27:229–246
Mascheroni P, Penta R (2016) The role of the microvascular network structure on diffusion and consumption of anti-cancer drugs. Int J Numer Methods Biomed Eng. doi:10.1002/cnm.2857
Mei CC, Vernescu B (2010) Homogenization methods for multiscale mechanics. World Scientific, Singapore
Meirmanov A (2010) Double porosity models for liquid filtration in incompressible poroelastic media. Math Models Methods Appl Sci 20(04):635–659
Mikelić A, Wheeler MF (2013) Convergence of iterative coupling for coupled flow and geomechanics. Comput Geosci 17(3):455–461
Netti PA, Baxter LT, Boucher Y, Skalak R, Jain RK (1995) Time-dependent behavior of interstitial fluid pressure in solid tumors: implications for drug delivery. Cancer Res 55(22):5451–5458
Papanicolau G, Bensoussan A, Lions JL (1978) Asymptotic analysis for periodic structures. Elsevier, Amsterdam
Penta R, Ambrosi D (2015) The role of the microvascular tortuosity in tumor transport phenomena. J Theor Biol 364:80–97
Penta R, Ambrosi D, Quarteroni A (2015) Multiscale homogenization for fluid and drug transport in vascularized malignant tissues. Math Models Methods Appl Sci 25(1):79–108
Penta R, Ambrosi D, Shipley RJ (2014) Effective governing equations for poroelastic growing media. Q J Mech Appl Math 67(1):69–91
Penta R, Gerisch A (2016) Investigation of the potential of asymptotic homogenization for elastic composites via a three-dimensional computational study. Comput Vis Sci 17(4):185–201
Penta R, Raum K, Grimal Q, Schrof S, Gerisch A (2016) Can a continuous mineral foam explain the stiffening of aged bone tissue? A micromechanical approach to mineral fusion in musculoskeletal tissues. Bioinspir Biomim 11(3):035004
Saffman PG (1971) On the boundary condition at the surface of a porous medium. Stud Appl Math 50(2):93–101
Sanchez-Palencia E (1980) Non-homogeneous media and vibration theory. Lecture Notes in Physics. Springer, Berlin
Sanchez-Palencia E (1983) Homogenization method for the study of composite media. Asymptotic analysis II. Lecture notes in mathematics. vol 985, pp 192–214
Shipley RJ, Chapman J (2010) Multiscale modelling of fluid and drug transport in vascular tumors. Bull Math Biol 72:1464–1491
Showalter RE (2005) Poroelastic filtration coupled to stokes flow. Lecture Notes in Pure and Applied Mathematics .vol 242, pp 229–237
Taffetani M, de Falco C, Penta R, Ambrosi D, Ciarletta P (2014) Biomechanical modelling in nanomedicine: multiscale approaches and future challenges. Arch Appl Mech 84(9–11):1627–1645
Wang H (2000) Theory of linear poroelasticity with applications to geomechanics and hydrogeology. Princeton University Press, Princeton
Weiner S, Wagner HD (1998) The material bone: structure-mechanical function relations. Ann Rev Mater Sci 28:271–298
Acknowledgements
The authors are indebted to Davide Ambrosi, Alf Gerisch, and Rebecca Shipley for useful hints and discussions concerning the content of this manuscript.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Funding
This study was funded by the Ministry of Economy in Spain, under the project reference DPI2014-58885-R.
Rights and permissions
About this article
Cite this article
Penta, R., Merodio, J. Homogenized modeling for vascularized poroelastic materials. Meccanica 52, 3321–3343 (2017). https://doi.org/10.1007/s11012-017-0625-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11012-017-0625-1