Research Article  Open Access
Roberto Palma, Antonio M. Callejas, "Extended Poroelasticity: An Analytical Solution and Its Application to pWave Propagation in Cervical Tissues", Mathematical Problems in Engineering, vol. 2018, Article ID 5280141, 6 pages, 2018. https://doi.org/10.1155/2018/5280141
Extended Poroelasticity: An Analytical Solution and Its Application to pWave Propagation in Cervical Tissues
Abstract
This work presents a semianalytical solution based on Laplace transform to study the behaviour of poroelastic materials in the context of the Extended Nonequilibrium Thermodynamics. In this framework, the fluid phase incorporates a relaxation time and, consequently, a frequencydependence appears. This rheological behaviour could explain the frequencydependence experimentally observed in biological tissues, which has traditionally attributed to the solid phase of tissues. In particular, the analytical solution is applied to two cases, heaviside and sinusoidal inputs, of a semiinfinite domain, which is filled with a material such as the human cervix. From the results, it is observed that the frequencydependence of the fluid phase could be relevant to high relaxation times while for null relaxation times the classical poroelastic theory is recovered. Finally, the present analytical solution could be used to validate future computational codes and experimental settings.
1. Introduction
Nowadays, understanding mechanical behaviour of the human cervix is a challenge for theoretical, computational, and experimental communities since it could be used to develop functional anticipation diagnostic tools, which will be applied to reduce the main cause of infant mortality: preterm birth, according to [1]. Notice that the understanding of cervical tissues is considered one of the most pressing problems in obstetrics [2].
According to [3], the human cervix is composed of a distribution of cells embedded in an extracellular matrix of fibrillar collagen, which can be considered as the solid phase and represent 2030% of the tissue and water with dissolved salts (fluid phase). The latter phase is the responsible for the cervix incompressibility [4].
There exist many models to study dynamic responses of tissues in specialised literature. However, most of them are contradictory since they are based on phenomenological equations, which lack robustness. In particular, measured values in cervix tissues differ several orders of magnitude for medium and high frequencies [5]. In this connection, there are several viscoelastic approaches to properly characterise the human cervix: a linear viscoelastic model was proposed in [6], an inhomogenous model in [7], and several approaches by the Rus’ group; see [8–10].
Despite the fact that the cervix is composed of solid and fluid, all the previous works are only focused on the solid contribution. On the contrary, the Extended Nonequilibrium Thermodynamics (ENET) [11] incorporates a viscouslike term to the fluid phase, which could explain the frequencydependence of the cervix. For instance, R. Palma et al. have developed analytical and numerical solutions to study the second sound in thermoelasticity (see [12, 13]), which can be considered as a viscothermal effect. Also, ENET is applied in [14] to study the Debye relaxation: a viscoelectric effect.
On this ground, this paper presents a poroelastic formulation based on the ENET to study the frequencydependence of cervix tissues due to the fluid phase, and it could be considered the main novelty of the present work. Then, the formulation is expressed in onedimension in order to obtain a dynamical solution by a semianalytical approach based on Laplace transform. Finally, two numerical examples are developed to highlight the importance of the fluidrelaxation in the response of the cervix.
The present work assumes small strains and linear elasticity for the solid phase. For the fluid, constant properties are considered and the effects of temperature are neglected.
2. ThreeDimensional Governing Equations
Consider an arbitrary domain and boundary for which the governing equations are composed of equilibrium and constitutive equations and of the boundary conditions.
With regard to the equilibrium equation and since the domain contains solid and fluid constituents, two equations must be enforced: linear momentum and mass conservation. Mathematically, they are expressed in local form as follows:where , , and denote mass density (including solid and fluid), acceleration, and body forces, respectively; is the Cauchy stress tensor; and denote the rate of change in fluid mass through the boundary and the production of fluid from an external source, respectively. Finally, , , , and denote volumetric strain, fluid pressure, and two constants closely related to the porosity and fluid bulk modulus, respectively.
Obviously, two constitutive equations, solid and fluid phases, are required to model poroelasticity; these equations readAt this point, it is necessary to define all terms in (2): (i) denotes the elastic fourth order tensor, which is composed of matrix and of fibre (both solid phases) by the rule of mixture: where is the percentage of matrix, denotes the fourthorder identity tensor, and , are the Lamé parameters for matrix and fibre constituents, respectively.(ii) is the small strain second order tensor and denotes the symmetric part of the gradient of displacements.(iii) denotes solid permeability and it is closely related to the fluid viscosity .(iv) is the relaxation times, which is introduced by the assumption of a mixed entropy; see [15]. This empirical parameter is responsible for viscosity in the fluid phase and, consequently, for the frequencydependence. Notice that the classical poroelasticity theory is recovered by imposing .
Finally, the Dirichlet and Neumann boundary conditions for the extended poroelastic problem readwhere and denote prescribed displacements and pressure, respectively, and and prescribed traction and fluid flux on the boundary with outward normal .
3. Analytical Solution
This section presents a onedimensional, semianalytical, and dynamical solution for a halfspace filled with a poroelastic material. For this purpose, the threedimension extended poroelastic equations reported in Section 2 are rewritten along the axis in order to apply the state space technique; see [16]. The threedimensional Euclidean coordinates becomeand the components of strain tensor are reduced toand, finally, the constitutive equation of (2) (upper) becomes
In absence of body forces and mass sources , the equilibrium equations of (1) are reduced to
Now, equations (6), (7), and (8) are rewritten in the Laplace domain by applying the transformation :These equations can be expressed in compact form by introducing the coefficients , , , and :to giveand the closed solution of this system of two couple equations can be expressed as follows:
The solution of this system is obtained by applying the CayleyHamilton theorem [16] to givewhere and are the boundary conditions and the coefficient and are explicitly given bywhere and are the solutions of the following characteristic equations:
Finally, the semianalytical solution is attained by imposing boundary conditions and by inverting the Laplace transform using Riemannsum approximations, as in [12].
4. Results
This section presents two analytical solutions, called cases, in order to highlight the main features of the present formulation. For this purpose, the material properties are obtained from the literature (see Table 1), and they are real measured variables of the human cervix. In particular, Lamé parameters for fibre and matrix phases are obtained from [8, 9], bulk modulus of fluid from [17], solid permeability from [18], and fluid viscosity from [19]. Finally, the coefficient obeys the composition of the cervix, namely, 8070% of fluid phase; see [20].

4.1. Case I
For case I, the boundary conditions are assumed to be a heaviside unit step function :
Introducing (16) in (12) and taking into account (14), the solutions for and in the Laplace domain readFurthermore, the mechanical displacement and the flux can be obtained taking into accountto give
Consider a onedimensional and semiinfinite domain, which is filled with a material as that of cervix. At , a fluid pressure of heaviside type is applied and, then, the evolution of all variables of (17) and (18) is numerically studied. In particular, Figure 1 shows the time evolution of (left column), (middle column), and (right column) for three different relaxation times: 0.5, 0,1, 0 [s] (each one in a row). In addition, each figure contains three different curves at three time instants 0.3, 0.5, 1 [s].
For [s] (extended poroelasticity) a hyperbolic behaviour is attained in the evolution of and, consequently, in . It is observed in the wave front, which is moving at each time instant. On the contrary, for the case a parabolic behaviour is observed, absence of wave front, and the classical poroelastic theory is recovered. In this sense, the extended poroelasticity removes the paradox of infinite velocity, which lacks physical meaning: in nature, there are no infinity velocities.
Since the problem is coupled, the stresses for [s] show a viscouslike behaviour that could explain the frequencydependence experimentally observed in cervical tissues. On the contrary, for [s], a linear behaviour is observed and, therefore, it is concluded that the classical poroelasticity can not be used to model the human cervix.
4.2. Case II
For case II, a sinusoidal function of frequency is prescribed:Following a procedure similar to that of case I, the solutions for and in the Laplace domain areand the mechanical displacement and the flux are given by
Consider the same semiinfinite domain as that in case I. Now, a fluid pressure of sinusoidal type is applied at and all variables of (21) and (22) are shown in Figure 2. Again, the same conclusion as those in case I can be observed for this sinusoidal input. Nevertheless, the curves are smoother due to the nature of the sinusoidal signal. Therefore, these solutions are more amenable for future computational validations since it is not necessary to use regularisation schemes.
5. Conclusions
This work has presented a theoretical approach based on Nonequilibrium Thermodynamics to study the behaviour of poroelastic materials taking into account the frequencydependence of the fluid phase. In this connection, the main novelty of the present work is the incorporation of relaxation times for the fluid phase to perform a material constitution applied to biological tissues. Then, the threedimensional governing equations are reduced to one dimension in order to obtain a semianalytical and dynamical solution based on Laplace transform. In particular, the solution is applied to simulate a semiinfinite domain, which is filled with a material such as the human cervix, and it is observed that the frequencydependence also could be due to the fluid phase.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the Ministry of Education DPI201783859 R, DPI201451870 R, DPI201017065, and UNGR15CE 3664, Ministry of Health DTS15/00093, and PI16/00339, PI01072017, PIN00302017, and Junta de Andalucía P11CTS 8089 projects.
References
 M. House and S. Socrate, “The cervix as a biomechanical structure,” Ultrasound in Obstetrics & Gynecology, vol. 28, no. 6, pp. 745–749, 2006. View at: Publisher Site  Google Scholar
 N. Uldbjerg, G. Ekman, A. Malmström, B. Sporrong, U. Ulmsten, and L. Wingerup, “Biochemical and morphological changes of human cervix after local application of prostaglandin e2 in pregnancy,” The Lancet, vol. 31, pp. 267268, 1981. View at: Publisher Site  Google Scholar
 M. Winkler and W. Rath, “Changes in the cervical extracellular matrix during pregnancy and parturition,” Journal of Perinatal Medicine, vol. 27, no. 1, pp. 45–60, 1999. View at: Google Scholar
 S. Weiss, T. Jaermann, P. Schmid et al., “Threedimensional fiber architecture of the nonpregnant human uterus determined ex vivo using magnetic resonance diffusion tensor imaging,” Anatomical Record  Part A Discoveries in Molecular, Cellular, and Evolutionary Biology, vol. 288, no. 1, pp. 84–90, 2006. View at: Publisher Site  Google Scholar
 E. Turgay, S. Salcudean, and R. Rohling, “Identifying the mechanical properties of tissue by ultrasound strain imaging,” Ultrasound in Medicine & Biology, vol. 32, no. 2, pp. 221–235, 2006. View at: Publisher Site  Google Scholar
 J. M. Pereira, J. M. Mansour, and B. R. Davis, “Dynamic measurement of the viscoelastic properties of skin,” Journal of Biomechanics, vol. 24, no. 2, pp. 157–162, 1991. View at: Publisher Site  Google Scholar
 E. A. Barannik, A. Girnyk, V. Tovstiak, A. I. Marusenko, S. Y. Emelianov, and A. P. Sarvazyan, “Doppler ultrasound detection of shear waves remotely induced in tissue phantoms and tissue in vitro,” Ultrasonics, vol. 40, no. 18, pp. 849–852, 2002. View at: Publisher Site  Google Scholar
 L. Peralta, G. Rus, N. Bochud, and F. S. Molina, “Assessing viscoelasticity of shear wave propagation in cervical tissue by multiscale computational simulation,” Journal of Biomechanics, vol. 48, no. 9, pp. 1549–1556, 2015. View at: Publisher Site  Google Scholar
 L. Peralta, G. Rus, N. Bochud, and F. S. Molina, “Mechanical assessment of cervical remodelling in pregnancy: Insight from a synthetic model,” Journal of Biomechanics, vol. 48, no. 9, pp. 1557–1565, 2015. View at: Publisher Site  Google Scholar
 A. Callejas, A. Gomez, J. Melchor et al., “Performance Study of a Torsional Wave Sensor and Cervical Tissue Characterization,” Sensors, vol. 17, no. 9, p. 2078, 2017. View at: Publisher Site  Google Scholar
 D. Jou, J. CasasVázquez, and G. Lebon, Extended Irreversible Thermodynamics, SpringerVerlag, Berlin, Heidelberg, Germany, 1996. View at: Publisher Site  MathSciNet
 R. Palma, J. L. PérezAparicio, and R. L. Taylor, “Nonlinear finite element formulation applied to thermoelectric materials under hyperbolic heat conduction model,” Computer Methods Applied Mechanics and Engineering, vol. 213–216, pp. 93–103, 2012. View at: Publisher Site  Google Scholar
 R. Palma, E. Moliner, and J. L. PérezAparicio, “Elastothermoelectric beam formulation for modeling thermoelectric devices,” Finite Elements in Analysis and Design, vol. 129, pp. 32–41, 2017. View at: Publisher Site  Google Scholar
 R. Palma, J. L. PerezAparicio, and R. L. Taylor, “Dissipative finite element formulation applied to piezoelectric materials with Debye memory,” IEEE/ASME Transactions on Mechatronics, 2018. View at: Google Scholar
 G. Lebon, T. Desaive, and P. Dauby, “A unified extended thermodynamic description of diffusion, thermodiffusion, suspensions, and porous media,” Transaction of the ASME, vol. 73, no. 1, pp. 16–20, 2006. View at: Publisher Site  Google Scholar
 H. M. Youssef and E. A. AlLehaibi, “Statespace approach of twotemperature generalized thermoelasticity of onedimensional problem,” International Journal of Solids and Structures, vol. 44, no. 5, pp. 1550–1562, 2007. View at: Publisher Site  Google Scholar
 J. B. Pendry and J. Li, “An acoustic metafluid: Realizing a broadband acoustic cloak,” New Journal of Physics , vol. 10, 2008. View at: Google Scholar
 F. J. O'Brien, B. A. Harley, M. A. Waller, I. V. Yannas, L. J. Gibson, and P. J. Prendergast, “The effect of pore size on permeability and cell attachment in collagen scaffolds for tissue engineering,” Technology and Health Care, vol. 15, no. 1, pp. 3–17, 2007. View at: Google Scholar
 L. Colla, L. Fedele, M. Scattolini, and S. Bobbo, “Waterbased Fe 2O 3 nanofluid characterization: Thermal conductivity and viscosity measurements and correlation,” Advances in Mechanical Engineering, vol. 2012, 2012. View at: Google Scholar
 P. Fratzl, “Collagen: structure and mechanics, an introduction,” Collagen: Structure and Mechanics, pp. 1–13, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Roberto Palma and Antonio M. Callejas. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.