This work presents a theoretical and numerical study of the flow of the interstitial fluid saturating a porous medium, principally aimed at modeling a bio-mimetic material and assumed to experience a dynamic regime different from the Darcian one, as is typically hypothesized in biomechanical scenarios. The main aspect of our research is the conjecture according to which, for a particular mechanical state of the porous medium, the fluid exhibits two types of deviation from Darcy’s law. One is due to the inertial forces characterizing the pore scale dynamics of the fluid. This aspect can be resolved by turning to the Forchheimer correction to Darcy’s law, which introduces non-linearities in the relationship between the fluid filtration velocity and the dissipative forces describing the interactions between the fluid and the solid matrix. The second source of discrepancies from classical Darcy’s law emerges, for example, when pore scale disturbances to the flow, such as obstructions of the fluid path or clogging of the pores, result in a time delay between drag force and filtration velocity. Recently, models have been proposed in which such delay is described through constitutive laws featuring fractional operators. Whereas, to the best of our knowledge, the aforementioned behaviors have been studied separately or in small deformations, we present a model of fluid flow in a deformable porous medium undergoing large deformations in which the fluid motion obeys a fractionalized Forchheimer’s correction. After reviewing Forchheimer’s formulation, we present a fractionalization of the Darcy–Forchheimer law, and we explain the numerical procedure adopted to solve the highly non-linear boundary value problem resulting from the presence of the two considered deviations from the Darcian regime. We complete our study by highlighting the way in which the fractional order of the model tunes the magnitude of the pore pressure and fluid filtration velocity.
AteshianG.On the theory of reactive mixtures for modeling biological growth. Biomech Model Mechan2007; 6(6): 423–445.
2.
HolmesMMowV.The nonlinear characteristics of soft gels and hydrated connective tissues in ultrafiltration. J Biomech1990; 23: 1145–1156.
3.
HuygheJMLoonRVBaaijensF.Fluid-solid mixtures and electrochemomechanics: the simplicity of lagrangian mixture theory. Comput Appl Math2004; 23(2–3): 235–258.
4.
LoretBSimõesF.A framework for deformation, generalized diffusion, mass transfer and growth in multi-species multi-phase biological tissues. Eur J Mech A2005; 24: 757–781.
5.
GrilloAFedericoSWittumG.Growth, mass transfer, and remodeling in fiber-reinforced, multi-constituent materials. Int J Nonlinear Mech2012; 47: 388–401.
6.
MaroudasABulloughP.Permeability of articular cartilage. Nature1968; 219(5160): 1260–1261.
7.
GuilakF.Compression-induced changes in the shape and volume of the chondrocyte nucleus. J Biomech1995; 28(12): 1529–1541.
8.
HolzapfelGGasserTOgdenR.A new constitutive framework for arterial wall mechanics and a comparative study of material models. J Elast2000; 61(1–3): 1–48.
9.
MollenhauerJAurichMMuehlemanC, et al. X-ray diffraction of the molecular substructure of human articular cartilage. Connect Tissue Res2003; 44(5): 201–207.
10.
MerodioJOgdenR.Mechanical response of fiber-reinforced incompressible non-linearly elastic solids. Int J Nonlinear Mech2005; 40(2–3): 213–227.
LanirY.Constitutive equations for fibrous connective tissues. J Biomech1983; 16: 1–12.
13.
GasserTOgdenRHolzapfelG.Hyperelastic modelling of arterial layers with distributed collagen fibre orientations. J R Soc Interface2006; 3: 15–35.
14.
FedericoSHerzogW.Towards an analytical model of soft biological tissues. J Biomech2008; 41: 3309–3313.
15.
FedericoSHerzogW.On the permeability of fibre-reinforced porous materials. Int J Solids Struct2008; 45: 2160–2172.
16.
FedericoSGasserT.Non-linear elasticity of biological tissues with statistical fibre orientation. J R Soc Interface2010; 7: 955–966.
17.
CarfagnaMGrilloA.The spherical design algorithm in the numerical simulation of biological tissues with statistical fibre-reinforcement. Comput Vis Sci2017; 18: 157–184.
18.
HolzapfelGAOgdenRWSherifovaS.On fibre dispersion modelling of soft biological tissues: a review. Proc R Soc A2019; 475(2224): 1–22.
19.
FedericoSGrilloA.Elasticity and permeability of porous fibre-reinforced materials under large deformations. Mech Mater2012; 44: 58–71.
20.
TomicAGrilloAFedericoS.Poroelastic materials reinforced by statistically oriented fibres — numerical implementation and application to articular cartilage. IMA J Appl Math2014; 79: 1027–1059.
21.
QuinnTMDierickxPGrodzinskyAJ.Glycosaminoglycan network geometry may contribute to anisotropic hydraulic permeability in cartilage under compression. J Biomech2001; 34(11): 1483–1490.
HolzapfelGAOgdenRW.Constitutive modelling of arteries. Proc R Soc A Math Phys Eng Sci2010; 466(2118): 1551–1597.
24.
PierceDMRickenTHolzapfelGA.A hyperelastic biphasic fibre-reinforced model of articular cartilage considering distributed collagen fibre orientations: continuum basis, computational aspects and applications. Comput Methods Biomech Biomed Engin2013; 16(12): 1344–1361.
25.
HashlamounKGrilloAFedericoS.Efficient evaluation of the material response of tissues reinforced by statistically oriented fibres. Z Angew Math Phys2016; 67: 113–145.
26.
LiKOgdenRHolzapfelG.Computational method for excluding fibers under compression in modeling soft fibrous solids. Eur J Mech A Solids2016; 57: 178–193.
27.
GizziAPandolfiAVastaM.A generalized statistical approach for modeling fiber-reinforced materials. J Eng Math2017; 109(1): 211–226.
28.
MowVCRatcliffeAPooleAR.Cartilage and diarthrodial joints as paradigms for hierarchical materials and structures. Biomaterials1992; 13(2): 67–97.
29.
HanSKFedericoSGrilloA, et al. The mechanical behaviour of chondrocytes predicted with a micro-structural model of articular cartilage. Biomech Model Mechan2006; 6(3): 139–150.
30.
GrilloACarfagnaMFedericoS.Non-Darcian flow in fibre-reinforced biological tissues. Meccanica2017; 52: 3299–3320.
31.
AlaimoGPiccoloVCutoloA, et al. A fractional order theory of poroelasticity. Mech Res Commun2019; 100: 103395.
32.
HassanizadehS. Derivation of basic equations of mass transp. porous med., Part 2. Generalized Darcy’s and Fick’s laws. Adv Water Resour1986; 9: 207–222.
33.
BearJBachmatY.Introduction to modeling of transport phenomena in porous media. Dordrecht: Kluwer, 1990.
34.
BennethumLMuradMCushmanJ.Macroscale thermodynamics and the chemical potential for swelling porous media. Transp Porous Media2000; 39(2): 187–225.
35.
KhaledARVafaiK.The role of porous media in modeling flow and heat transfer in biological tissues. Int J Heat Mass Transf2003; 46(26): 4989–5003.
36.
GrilloACarfagnaMFedericoS.The Darcy-Forchheimer law for modelling fluid flow in biological tissues. Theor Appl Mech2014; 41(4): 283–322.
37.
ThauvinFMohantyKK.Network modeling of non-darcy flow through porous media. Transp Porous Media1998; 31(1): 19–37.
38.
DeseriLZingalesM.A mechanical picture of fractional-order Darcy equation. Commun Nonlinear Sci Numer Simul2015; 20(3): 940–949.
39.
JangJChenJ.Variable porosity and thermal dispersion effects on vortex instability of a horizontal natural convection flow in a saturated porous medium. Wärme Stoffübertragung1994; 29: 153–160.
40.
BulleRAlottaGMarchioriG, et al. The human meniscus behaves as a functionally graded fractional porous medium under confined compression conditions. Appl Sci2021; 11(20).
41.
IaffaldanoGCaputoMMartinoS.Experimental and theoretical memory diffusion of water in sand. Hydrol Earth Syst Sc2005; 10: 93–100.
42.
CaputoM.Models of flux in porous media with memory. Water Resour Res2000; 36(3): 693–705.
43.
SaporaACornettiPChiaiaB, et al. Nonlocal diffusion in porous media: a spatial fractional approach. J Eng Mech2016; 143: D4016007.
44.
GundaSNatarajanSBarreraO.On the fractional transversely isotropic functionally graded nature of soft biological tissues: application to the meniscal tissue. J Mech Behav Biomed Mater2023; 143: 105855.
45.
BarreraO.A unified modelling and simulation for coupled anomalous transp porous media and its finite element implementation. Comput Mech2021; 68: 1267–1282.
46.
CaputoMPlastinoW.Diffusion in porous layers with memory. Geophys J Int2004; 158(1): 385–396.
47.
UzunerSLiLKucukS, et al. Changes in knee joint mechanics after medial meniscectomy determined with a poromechanical model. J Biomech Eng2020; 142(10).
48.
AmiriFBolognaENuzzoG, et al. Fractional-order poromechanics for a fully saturated biological tissue: biomechanics of meniscus. Int J Numer Meth Bio2023; .
49.
CarpinteriAMainardiF (eds) Fractals and fractional calculus in continuum mechanics. Vienna: Springer, 1997.
50.
PaolaMDPirrottaAValenzaA.Visco-elastic behavior through fractional calculus: an easier method for best fitting experimental results. Mech Mater2011; 43(12): 799–806.
51.
IonescuCLopesACopotD, et al. The role of fractional calculus in modeling biological phenomena: a review. Commun Nonlinear Sci Numer Simul2017; 51: 141–159.
52.
MaesJADonahueTH.Time dependent properties of bovine meniscal attachments: stress relaxation and creep. J Biomech2006; 39(16): 3055–3061.
53.
CamardaLBolognaEPavanD, et al. Posterior meniscal root repair: a biomechanical comparison between human and porcine menisci. Muscles Ligaments Tendons J2019; 09(01): 76.
54.
BolognaEPaolaMDDayalK, et al. Fractional-order nonlinear hereditariness of tendons and ligaments of the human knee. Philos Trans R Soc Math Phys Eng Sci2020; 378(2172): 20190294.
55.
NorbergCFilipponeGAndreopoulosF, et al. Viscoelastic and equilibrium shear properties of human meniscus: relationships with tissue structure and composition. J Biomech2021; 120: 110343.
56.
PajerskiJ.Nonlinear biphasic microstructural numerical analysis of articular cartilage and chondrocytes. MSc Thesis, The University of Calgary, Calgary, AB, Canada, 2010.
57.
MaginRLRoystonTJ.Fractional-order elastic models of cartilage: a multi-scale approach. Commun Nonlinear Sci Numer Simul2010; 15(3): 657–664.
58.
SmythPAGreenI.Fractional calculus model of articular cartilage based on experimental stress-relaxation. Mech Time Depend Mat2015; 19(2): 209–228.
59.
QuiligottiS.On bulk growth mechanics of solid-fluid mixtures: kinematics and invariance requirements. Theoret Appl Mech2002; 28–29: 277–288.
60.
QuiligottiSMauginGdell’IsolaF.An eshelbian approach to the nonlinear mechanics of constrained solid-fluid mixtures. Acta Mech2003; 160: 45–60.
61.
FedericoSHerzogW.On the anisotropy and inhomogeneity of permeability in articular cartilage. Biomech Model Mechanobiol2008; 7: 367–378.
62.
HassanizadehMGrayWG.General conservation equations for multi-phase systems: 1. averaging procedure. Adv Water Resour1979; 2: 131–144.
63.
MarsdenJHughesT.Mathematical foundations of elasticity. Mineola, NY: Dover Publications, Inc., 1983.
64.
FedericoSAlhasadiMFGrilloA.Eshelby’s inclusion theory in light of Noether’s theorem. Math Mech Complex Syst2019; 7(3): 247–285.
65.
FedericoSGrilloASegevR.Material description of fluxes in terms of differential forms. Continuum Mech Therm2015; 28(1–2): 379–390.
66.
GrilloAGiversoCFavinoM, et al. Mass transp porous media with variable mass. In: DelgadoJMPQBarbosa de LimaAGVázquez da SilvaM (eds) Advanced structured materials. Berlin: Springer, 2012, pp. 27–61. DOI:10.1007/978-3-642-30532-0)_2.
67.
BonetJWoodR.Nonlinear continuum mechanics for finite element analysis. New York: Cambridge University Press, 2008.
68.
HassanizadehSMLeijnseA.A non-linear theory of high-concentration-gradient dispersion in porous media. Adv Water Resour1995; 18(4): 203–215.
69.
HassanizadehSM.Derivation of basic equations of mass Transp Porous Media, Part 1. Macroscopic balance laws. Adv Water Resour1986; 9(4): 196–206.
70.
BennethumLGiorgiT.Generalized Forchheimer equation for two-phase flow based on hybrid mixture theory. Transp Porous Media1997; 26: 261–275.
71.
KöpfMCorinthCHaferkampO, et al. Anomalous diffusion of water in biological tissues. Biophys J1996; 70(6): 2950–2958.
72.
ÖzarslanEBasserPJShepherdTM, et al. Observation of anomalous diffusion in excised tissue by characterizing the diffusion-time dependence of the MR signal. J Magn Reson2006; 183(2): 315–323.
73.
MaginRLAbdullahOBaleanuD, et al. Anomalous diffusion expressed through fractional order differential operators in the bloch–torrey equation. J Magn Reson2008; 190(2): 255–270.
74.
MaginRLIngoCColon-PerezL, et al. Characterization of anomalous diffusion in porous biological tissues using fractional order derivatives and entropy. Micropor Mesopor Mat2013; 178: 39–43.
75.
LiuI.Method of lagrange multipliers for exploitation of the entropy principle. Archive Rational Mech Anal1972; 46: 131–148.
76.
AteshianGHumphreyJ.Continuum mixture models of biological growth and remodeling: past successes and future opportunities. Annu Rev Biomed Eng2012; 14(1): 97–111.
77.
AtanackovićTKonjikSOparnicaL, et al. The Cattaneo type space-time fractional heat conduction equation. Continuum Mech Therm2011; 24(4–6): 293–311.
78.
ŽeliVZoricaD.Analytical and numerical treatment of the heat conduction equation obtained via time-fractional distributed-order heat conduction law. Physica A2018; 492: 2316–2335.
79.
FedericoS.The Truesdell rate in continuum mechanics. Z Angew Math Phys2022; 73(3): 109.
80.
PodlubnyI.Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. New York: Academic Press, 1998.
GrilloALogashenkoDStichelS, et al. Forchheimer’s correction in modelling flow and transport in fractured porous media. Comput Vis Sci2012; 15: 169–190.
83.
TekMCoatsKKatzD.The effect of turbulence on flow of natural gas through porous reservoirs. J Petrol Technol1962; 14(7): 799–806.
84.
GeertsmaJ.Estimating the coefficient of inertial resistance in fluid flow through porous media. Soc Petrol Eng J1974; 14(5): 445–450.
85.
WangXThauvinFMohantyK.Non-darcy flow through anisotropic porous media. Chem Eng Sci1999; 54(12): 1859–1869.
86.
HuangHAyoubJ.Applicability of the Forchheimer equation for non-Darcy flow in porous media. SPE J2008; 13(01): 112–122.
87.
InnocentiniMDMFaleirosRKPisaniR, et al. Permeability of porous gelcast scaffolds for bone tissue engineering. J Porous Mat2009; 17(5): 615–627.
88.
GrilloAProhlRWittumG.A poroplastic model of structural reorganisation in porous media of biomechanical interest. Continuum Mech Therm2016; 28: 579–601.
89.
KnodelMDi StefanoSNägelA, et al. An efficient algorithm for biomechanical problems based on a fully implicit nested newton solver. Theor Appl2022; 49(2): 183–221.
90.
GrilloAProhlRWittumG.A generalised algorithm for anelastic processes in elastoplasticity and biomechanics. Math Mech Solids2017; 22(3): 502–527.
91.
SimoJHughesT.Computational inelasticity. New York: Springer, 1998.
92.
ABAQUS Theory Manual. Abaqus 6.11 Theory Manual. Johnston, RI: Dassault Systèmes Simulia Corp, 2011.
93.
GiversoCStefanoSDGrilloA, et al. A three dimensional model of multicellular aggregate compression. Soft Matter2019; 15(48): 10005–10019.
94.
HolmesMH.Finite deformation of soft tissue: analysis of a mixture model in uni-axial compression. J Biomech Eng1986; 108(4): 372–381.
95.
Ramírez-TorresAPentaRGrilloA.Effective properties of fractional viscoelastic composites via two⣳scale asymptotic homogenization. Math Methods Appl Sci2023; 46(16): 16500–16520.
96.
WaghorneJBonomoFRabbaniA, et al. On the characteristics of natural hydraulic dampers: an image-based approach to study the fluid flow behaviour inside the human meniscal tissue. Acta Biomater2023; 175: 157–169.
97.
AgustoniGMaritzJKennedyJ, et al. High resolution micro-computed tomography reveals a network of collagen channels in the body region of the knee meniscus. Ann Biomed Eng2021; 49: 2273–2281.
98.
VetriVDragnevskiKTkaczykM, et al. Advanced microscopy analysis of the micro-nanoscale architecture of human menisci. Sci Rep2019; 9(1): 18732.
99.
CompteAMetzlerR.The generalized Cattaneo equation for the description of anomalous transport processes. J Phys A Math Gen1997; 30(21): 7277–7289.
100.
Ramrez-TorresADi StefanoSGrilloA.Influence of non-local diffusion in avascular tumour growth. Math Mech Solids2021; 26(9): 1264–1293.
101.
BarreraOTarletonETangH, et al. Modelling the coupling between hydrogen diffusion and the mechanical behaviour of metals. Comput Mater Sci2016; 122: 219–228.
102.
SreejithPSrikanthKKannanK, et al. A thermodynamic framework for additive manufacturing of crystallizing polymers, Part II: simulation of the printing of a stent. Int J Eng Sci2023; 184: 103790.
