Articles | Volume 16, issue 10 
            
                
                    
            
            
            https://doi.org/10.5194/se-16-1227-2025
                    © Author(s) 2025. This work is distributed under 
the Creative Commons Attribution 4.0 License.
                the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/se-16-1227-2025
                    © Author(s) 2025. This work is distributed under 
the Creative Commons Attribution 4.0 License.
                the Creative Commons Attribution 4.0 License.
Revisiting Gassmann-type relationships within Biot poroelastic theory
Yury Alkhimenkov
CORRESPONDING AUTHOR
                                            
                                    
                                            Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
                                        
                                    Yury Y. Podladchikov
                                            Institute of Earth Sciences, University of Lausanne, Lausanne, Switzerland
                                        
                                    Related authors
Yury Alkhimenkov and Yury Y. Podladchikov
                                    Geosci. Model Dev., 18, 563–583, https://doi.org/10.5194/gmd-18-563-2025, https://doi.org/10.5194/gmd-18-563-2025, 2025
                                    Short summary
                                    Short summary
                                            
                                                The accelerated pseudo-transient (APT) method is an efficient way to solve partial differential equations, particularly well-suited for parallel computing. This paper explores the APT method's effectiveness in solving elastic, viscoelastic, and hydromechanical problems, focusing on quasi-static conditions in 1D, 2D, and 3D. The study examines the best numerical settings for fast and accurate solutions. The paper shows how the APT method can handle complex problems in high-resolution models.
                                            
                                            
                                        Yury Alkhimenkov, Lyudmila Khakimova, and Yury Podladchikov
                                        EGUsphere, https://doi.org/10.5194/egusphere-2024-3237, https://doi.org/10.5194/egusphere-2024-3237, 2024
                                    Short summary
                                    Short summary
                                            
                                                This study examines stress drops and earthquake nucleation in elasto-plastic media using 2D simulations, highlighting the importance of high temporal and spatial resolutions in capturing stress evolution and strain fields. Stress drops reflect fault rupture mechanics and emulate earthquake behavior. The non-Gaussian distribution of stress drop amplitudes resembles "solid turbulence." Elasto-plastic models simulate key earthquake processes and could improve seismic hazard assessment.
                                            
                                            
                                        Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
                                    Geosci. Model Dev., 14, 7749–7774, https://doi.org/10.5194/gmd-14-7749-2021, https://doi.org/10.5194/gmd-14-7749-2021, 2021
                                    Short summary
                                    Short summary
                                            
                                                We propose an implementation of the material point method using graphical processing units (GPUs) to solve elastoplastic problems in three-dimensional configurations, such as the granular collapse or the slumping mechanics, i.e., landslide. The computational power of GPUs promotes fast code executions, compared to a traditional implementation using central processing units (CPUs). This allows us to study complex three-dimensional problems tackling high spatial resolution.
                                            
                                            
                                        Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
                                    Geosci. Model Dev., 13, 6265–6284, https://doi.org/10.5194/gmd-13-6265-2020, https://doi.org/10.5194/gmd-13-6265-2020, 2020
                                    Short summary
                                    Short summary
                                            
                                                In this work, we present an efficient and fast material point method (MPM) implementation in MATLAB. We first discuss the vectorization strategies to adapt this numerical method to a MATLAB implementation. We report excellent agreement of the solver compared with classical analysis among the MPM community, such as the cantilever beam problem. The solver achieves a performance gain of 28 compared with a classical iterative implementation.
                                            
                                            
                                        Yury Alkhimenkov and Yury Y. Podladchikov
                                    Geosci. Model Dev., 18, 563–583, https://doi.org/10.5194/gmd-18-563-2025, https://doi.org/10.5194/gmd-18-563-2025, 2025
                                    Short summary
                                    Short summary
                                            
                                                The accelerated pseudo-transient (APT) method is an efficient way to solve partial differential equations, particularly well-suited for parallel computing. This paper explores the APT method's effectiveness in solving elastic, viscoelastic, and hydromechanical problems, focusing on quasi-static conditions in 1D, 2D, and 3D. The study examines the best numerical settings for fast and accurate solutions. The paper shows how the APT method can handle complex problems in high-resolution models.
                                            
                                            
                                        Yury Alkhimenkov, Lyudmila Khakimova, and Yury Podladchikov
                                        EGUsphere, https://doi.org/10.5194/egusphere-2024-3237, https://doi.org/10.5194/egusphere-2024-3237, 2024
                                    Short summary
                                    Short summary
                                            
                                                This study examines stress drops and earthquake nucleation in elasto-plastic media using 2D simulations, highlighting the importance of high temporal and spatial resolutions in capturing stress evolution and strain fields. Stress drops reflect fault rupture mechanics and emulate earthquake behavior. The non-Gaussian distribution of stress drop amplitudes resembles "solid turbulence." Elasto-plastic models simulate key earthquake processes and could improve seismic hazard assessment.
                                            
                                            
                                        Ludovic Räss, Ivan Utkin, Thibault Duretz, Samuel Omlin, and Yuri Y. Podladchikov
                                    Geosci. Model Dev., 15, 5757–5786, https://doi.org/10.5194/gmd-15-5757-2022, https://doi.org/10.5194/gmd-15-5757-2022, 2022
                                    Short summary
                                    Short summary
                                            
                                                Continuum mechanics-based modelling of physical processes at large scale requires huge computational resources provided by massively parallel hardware such as graphical processing units. We present a suite of numerical algorithms, implemented using the Julia language, that efficiently leverages the parallelism. We demonstrate that our implementation is efficient, scalable and robust and showcase applications to various geophysical problems.
                                            
                                            
                                        Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
                                    Geosci. Model Dev., 14, 7749–7774, https://doi.org/10.5194/gmd-14-7749-2021, https://doi.org/10.5194/gmd-14-7749-2021, 2021
                                    Short summary
                                    Short summary
                                            
                                                We propose an implementation of the material point method using graphical processing units (GPUs) to solve elastoplastic problems in three-dimensional configurations, such as the granular collapse or the slumping mechanics, i.e., landslide. The computational power of GPUs promotes fast code executions, compared to a traditional implementation using central processing units (CPUs). This allows us to study complex three-dimensional problems tackling high spatial resolution.
                                            
                                            
                                        Martin Franz, Michel Jaboyedoff, Ryan P. Mulligan, Yury Podladchikov, and W. Andy Take
                                    Nat. Hazards Earth Syst. Sci., 21, 1229–1245, https://doi.org/10.5194/nhess-21-1229-2021, https://doi.org/10.5194/nhess-21-1229-2021, 2021
                                    Short summary
                                    Short summary
                                            
                                                A landslide-generated tsunami is a complex phenomenon that involves landslide dynamics, wave dynamics and their interaction. This phenomenon threatens numerous lives and infrastructures around the world. To assess this natural hazard, we developed an efficient numerical model able to simulate the landslide, the momentum transfer and the wave all at once. The good agreement between the numerical simulations and physical experiments validates our model and its novel momentum transfer approach.
                                            
                                            
                                        Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
                                    Geosci. Model Dev., 13, 6265–6284, https://doi.org/10.5194/gmd-13-6265-2020, https://doi.org/10.5194/gmd-13-6265-2020, 2020
                                    Short summary
                                    Short summary
                                            
                                                In this work, we present an efficient and fast material point method (MPM) implementation in MATLAB. We first discuss the vectorization strategies to adapt this numerical method to a MATLAB implementation. We report excellent agreement of the solver compared with classical analysis among the MPM community, such as the cantilever beam problem. The solver achieves a performance gain of 28 compared with a classical iterative implementation.
                                            
                                            
                                        Cited articles
                        
                        Alkhimenkov, Y. and Podladchikov, Y.: Biot_Gassmann, Zenodo [code], https://doi.org/10.5281/zenodo.13942952, 2024. a, b
                    
                
                        
                        Alkhimenkov, Y. and Podladchikov, Y.: Revisiting Gassmann-Type Relationships within Biot Poroelastic Theory – Symbolic Derivations and Thermodynamic Validation, Zenodo [code], https://doi.org/10.5281/zenodo.15777522 (last access: 30 June 2025), 2025. a, b, c
                    
                
                        
                        Alkhimenkov, Y., Caspari, E., Lissa, S., and Quintal, B.: Azimuth-, angle- and frequency-dependent seismic velocities of cracked rocks due to squirt flow, Solid Earth, 11, 855–871, https://doi.org/10.5194/se-11-855-2020, 2020b. a, b
                    
                
                        
                        Berryman, J. G.: Origin of Gassmann’s equations, Geophysics, 64, 1627–1629, 1999. a
                    
                
                        
                        Berryman, J. G. and Milton, G. W.: Exact results for generalized Gassmann’s equations in composite porous media with two constituents, Geophysics, 56, 1950–1960, 1991. a
                    
                
                        
                        Biot and Willis, D.: The elastic coeff cients of the theory of consolidation, J. Appl. Mech., 15, 594–601, 1957. a
                    
                
                        
                        Biot, M.: Theory of Propagation of Elastic Waves in a Fluid-Saturated Porous Solid. I. Low-Frequency Range, The Journal of the Acoustical Society of America, 28, 168–178, 1956. a
                    
                
                        
                        Burridge, R. and Keller, J. B.: Poroelasticity equations derived from microstructure, The Journal of the Acoustical Society of America, 70, 1140–1146, 1981. a
                    
                
                        
                        Cheng, A. H.-D.: Poroelasticity, Springer, 2016. a
                    
                
                        
                        Coussy, O.: Poromechanics, John Wiley & Sons, ISBN 9780470849200, https://doi.org/10.1002/0470092718, 2004. a
                    
                
                        
                        Coussy, O.: Mechanics and physics of porous solids, John Wiley & Sons, 296 pp., ISBN 978-0-470-72135-3, 2011. a
                    
                
                        
                        Coussy, O., Dormieux, L., and Detournay, E.: From mixture theory to Biot’s approach for porous media, International Journal of Solids and Structures, 35, 4619–4635, 1998. a
                    
                
                        
                        Dormieux, L., Kondo, D., and Ulm, F.-J.: Microporomechanics, John Wiley & Sons, ISBN 9780470031889, 2006. a
                    
                
                        
                        Fortin, J. and Guéguen, Y.: Porous and cracked rocks elasticity: Macroscopic poroelasticity and effective media theory, Mathematics and Mechanics of Solids, 26, 1158–1172, 2021. a
                    
                
                        
                        Frenkel, J.: On the theory of seismic and seismoelectric phenomena in a moist soil, Journal of Physics, III, 230–241, 1944. a
                    
                
                        
                        Guéguen, Y. and Boutéca, M.: Mechanics of fluid-saturated rocks, Elsevier, ISBN 9780123053558, 2004. a
                    
                
                        
                        Gurevich, B.: Comparison of the low-frequency predicitons of Biot's and de Boer's poroelasticity theories with Gassmann’s equation, Applied physics letters, 91, 091919, https://doi.org/10.1063/1.2778763, 2007. a
                    
                
                        
                        Gyarmati, I.: Non-equilibrium thermodynamics, Springer, 1970. a
                    
                
                        
                        Jou, D., Casas-Vázquez, J., Lebon, G., Jou, D., Casas-Vázquez, J., and Lebon, G.: Extended irreversible thermodynamics, Springer, https://doi.org/10.1007/978-90-481-3074-0, 1996. a
                    
                
                        
                        Lopatnikov, S. L. and Cheng, A.-D.: Macroscopic Lagrangian formulation of poroelasticity with porosity dynamics, Journal of the Mechanics and Physics of Solids, 52, 2801–2839, 2004. a
                    
                
                        
                        Mavko, G., Mukerji, T., and Dvorkin, J.: The rock physics handbook, Cambridge University Press, https://doi.org/10.1017/9781108333016, 2020. a
                    
                
                        
                        Müller, T. M. and Sahay, P. N.: Elastic potential energy in linear poroelasticity, Geophysics, 84, W1–W20, 2019. a
                    
                
                        
                        Nur, A. and Byerlee, J.: An exact effective stress law for elastic deformation of rock with fluids, J. Geophys. Res., 76, 6414–6419, 1971. a
                    
                
                        
                        Pelissier, M. A., Hoeber, H., van de Coevering, N., and Jones, I. F.: Classics of elastic wave theory, Society of Exploration Geophysicists, 24, ISBN 9781560801429, https://doi.org/10.1190/1.9781560801931, 2007. a
                    
                
                        
                        Pride, S. R. and Garambois, S.: Electroseismic wave theory of Frenkel and more recent developments, Journal of Engineering Mechanics, 131, 898–907, 2005. a
                    
                
                        
                        Sahay, P. N.: Biot constitutive relation and porosity perturbation equation, Geophysics, 78, L57–L67, 2013. a
                    
                
                        
                        Schenk, O. and Gärtner, K.: Solving unsymmetric sparse systems of linear equations with PARDISO, Future Generation Computer Systems, 20, 475–487, 2004. a
                    
                
                        
                        Sevostianov, I.: Gassmann equation and replacement relations in micromechanics: A review, International Journal of Engineering Science, 154, 103344, https://doi.org/10.1016/j.ijengsci.2020.103344, 2020.  a
                    
                
                        
                        Smith, T. M., Sondergeld, C. H., and Rai, C. S.: Gassmann fluid substitutions: A tutorial, Geophysics, 68, 430–440, 2003. a
                    
                
                        
                        Tarokh, A. and Makhnenko, R. Y.: Remarks on the solid and bulk responses of fluid-filled porous rock, Geophysics, 84, WA83–WA95, 2019. a
                    
                
                        
                        Thomsen, L.: Discussion on the logical error in Gassmann poroelasticity: Numerical computations and generalizations, https://doi.org/10.1111/1365-2478.13290, 2023a. a, b, c
                    
                
                        
                        Ulm, F. and Coussy, O.: Mechanics and Durability of Solids, no. v. 1 in MIT/Prentice Hall series on civil, environmental, and systems engineering, Prentice Hall, ISBN 9780130479570, 2003. a
                    
                
                        
                        Wang, H.: Theory of linear poroelasticity with applications to geomechanics and hydrogeology, Princeton university press, ISBN 9780691037462, 2000. a
                    
                Short summary
            We present a thermodynamically consistent derivation of extended Biot poroelasticity, showing that Gassmann, Brown–Korringa, and related models emerge as special cases. Our formulation clarifies the conditions under which Gassmann’s relation holds and extends it by incorporating off-diagonal Hessian terms. Symbolic Maple code with consistency checks ensures full transparency, reproducibility, and accessibility for further research.
            We present a thermodynamically consistent derivation of extended Biot poroelasticity, showing...
            
         
 
                        
                                         
                        
                                         
                        
                                         
                        
                                         
                        
                                         
             
             
            