Verification of a sequential thermo-poroelasticity formulation in PFLOTRAN
Pith reviewed 2026-07-02 01:26 UTC · model grok-4.3
The pith
A sequential non-iterative fixed-stress split for thermo-poroelastic equations matches analytical benchmark solutions for pressure, temperature, and deformation.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The thermo-hydrologic equations are solved on control-volume blocks while the quasi-static momentum balance is solved on an element-based dual mesh, with coupling via a strictly sequential non-iterative fixed-stress split. This produces solutions that agree with analytical or semi-analytical responses for pressure diffusion, temperature field, and mechanical deformation in poroelastic and thermo-poroelastic benchmarks. A mapping treatment for discontinuities is also validated against an analytical solution.
What carries the argument
The strictly sequential, non-iterative fixed-stress split that solves the TH system implicitly then updates mechanics.
If this is right
- Accurate modeling of coupled processes in geologic porous media becomes feasible without iterative coupling.
- The discontinuity treatment enables simulation of fractures in THM contexts.
- The approach supports applications involving enhanced geothermal systems and subsurface energy storage.
- Verification establishes reliability for pressure, temperature, and deformation fields.
Where Pith is reading between the lines
- For problems with moderate coupling strength, avoiding iteration may reduce computational cost without loss of accuracy.
- The dual mesh approach for mechanics could be adapted to other simulation frameworks.
- Testing the method on problems with stronger nonlinearities or longer time scales would further confirm its robustness.
Load-bearing premise
The non-iterative fixed-stress split produces accurate enough coupled solutions for the selected benchmark problems and time scales without needing iteration or extra stabilization.
What would settle it
Observation of significant deviation between the computed pressure diffusion or displacement fields and the corresponding analytical solutions in any of the verification benchmarks.
read the original abstract
We present the verification of a thermo--hydrologic--mechanical capability implemented within the PFLOTRAN framework, with emphasis on benchmark-based assessment of the THM implementation. The thermal--hydrologic (TH) equations for mass and energy balance are solved on control-volume blocks or Voronoi cells, while the quasi-static momentum balance is solved on an element-based dual mesh. The coupling is achieved using a strictly sequential, non-iterative fixed-stress split strategy in which the TH system is solved implicitly for pressure and temperature, followed by a mechanics update for the displacement unknowns. Several verification problems are set up against poroelastic and thermo-poroelastic benchmarks, demonstrating agreement with analytical or semi-analytical benchmark responses for pressure diffusion, the temperature field, and mechanical deformation. In addition, we propose a treatment for discontinuities (e.g., fractures) based on mapping between mechanical and flow degrees of freedom, and validate the approach by comparison to an analytical solution. This work establishes the basis for thermo-poroelastic coupling in PFLOTRAN and provides a solid modeling foundation for a range of applications (e.g., enhanced geothermal systems and other subsurface energy storage) involving coupled thermal--hydrologic--mechanical (THM) processes in geologic porous media.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents the verification of a thermo-hydrologic-mechanical (THM) capability in PFLOTRAN. The TH equations are solved on control-volume or Voronoi cells and mechanics on a dual mesh, coupled via a strictly sequential non-iterative fixed-stress split (TH solved implicitly, followed by mechanics update). Verification is shown via agreement with analytical/semi-analytical benchmarks for pressure diffusion, temperature, deformation, and a discontinuity treatment.
Significance. If the non-iterative split accuracy holds beyond the tested cases, the work supplies a verified THM foundation in PFLOTRAN for applications such as enhanced geothermal systems. Credit is due for performing verification against external analytical benchmarks rather than self-referential fitting.
major comments (2)
- [Verification results / Abstract] Verification sections: claims of agreement with benchmarks are stated but no quantitative error metrics (L2 norms, relative errors, maximum deviations), mesh convergence studies, or time-step refinement data are reported. This leaves the central verification claim unquantified.
- [Coupling strategy / Verification problems] Coupling description and results: the strictly sequential non-iterative fixed-stress split is adopted without analysis of splitting error (O(Δt) truncation error known to depend on Biot coefficient, hydraulic diffusivity, and consolidation time). No comparison to an iterative or monolithic solver is provided, which is load-bearing for asserting that the formulation is verified for the stated applications.
minor comments (1)
- [Abstract] Abstract: states 'several verification problems' without enumerating them or indicating which benchmark addresses which field (pressure, temperature, mechanics).
Simulated Author's Rebuttal
We thank the referee for their constructive feedback on our manuscript. We address each major comment below.
read point-by-point responses
-
Referee: [Verification results / Abstract] Verification sections: claims of agreement with benchmarks are stated but no quantitative error metrics (L2 norms, relative errors, maximum deviations), mesh convergence studies, or time-step refinement data are reported. This leaves the central verification claim unquantified.
Authors: We agree that quantitative error metrics would strengthen the presentation. In the revised manuscript we will report L2 norms, relative errors, and maximum deviations for pressure, temperature, and displacement where analytical solutions exist. We will also add a mesh-convergence study for one benchmark and time-step refinement results for the diffusion problems. revision: yes
-
Referee: [Coupling strategy / Verification problems] Coupling description and results: the strictly sequential non-iterative fixed-stress split is adopted without analysis of splitting error (O(Δt) truncation error known to depend on Biot coefficient, hydraulic diffusivity, and consolidation time). No comparison to an iterative or monolithic solver is provided, which is load-bearing for asserting that the formulation is verified for the stated applications.
Authors: The verification is performed directly against analytical solutions of the coupled THM system; the observed agreement therefore already bounds the practical splitting error for the tested parameters and time scales. We will add a short discussion of the known O(Δt) truncation error and its dependence on the Biot coefficient. A comparison to monolithic or iterative solvers is not possible within the current PFLOTRAN implementation and lies outside the scope of this verification study. revision: partial
Circularity Check
Verification against external analytical benchmarks shows no circularity
full rationale
The paper implements a sequential non-iterative fixed-stress THM solver in PFLOTRAN and verifies the implementation by direct comparison to independent analytical and semi-analytical benchmark solutions for pressure diffusion, temperature, and deformation (including a discontinuity treatment). No parameters are fitted to the verification data and then re-used as 'predictions'; the benchmarks are external references. Any self-citations are incidental and not load-bearing for the verification claims. The central result is numerical agreement with outside solutions, rendering the reported derivation chain self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Heinemann, N., Alcalde, J., Miocic, J.M., Hangx, S.J.T., Kallmeyer, J., Ostertag-Henning, C., Hassanpouryouzband, A., Thaysen, E.M., Strobel, G.J., Schmidt-Hattenberger, C., Edlmann, K., Wilkinson, M., Bentham, M., Stuart Haszeldine, R., Carbonell, R., Rudloff, A.: Enabling large- scale hydrogen storage in porous media – the scientific challenges. Energy ...
-
[2]
Geotechnical and Geological Engineering30, 525–551 (2012) https://doi.org/10.1007/s10706-011-9491-0
Rutqvist, J.: The geomechanics of co 2 storage in deep sedimentary formations. Geotechnical and Geological Engineering30, 525–551 (2012) https://doi.org/10.1007/s10706-011-9491-0
-
[3]
McClure, M.W., Horne, R.N.: An investigation of stimulation mechanisms in enhanced geothermal systems. International Journal of Rock Mechanics and Mining Sciences72, 242–260 (2014) https: //doi.org/10.1016/j.ijrmms.2014.07.011
-
[4]
General theory of three-dimensional consolidation
Biot, M.A.: General theory of three-dimensional consolidation. Journal of Applied Physics12(2), 155–164 (1941) https://doi.org/10.1063/1.1712886
-
[5]
Journal of Applied Physics26(2), 182–185 (1955) https://doi.org/10.1063/1.1721956
Biot, M.A.: Theory of elasticity and consolidation for a porous anisotropic solid. Journal of Applied Physics26(2), 182–185 (1955) https://doi.org/10.1063/1.1721956
-
[6]
Princeton University Press, Princeton, NJ (2000).http://www.jstor.org/stable/j.ctt1jktrr4
Wang, H.F.: Theory of Linear Poroelasticity with Applications to Geomechanics and Hydrogeology. Princeton University Press, Princeton, NJ (2000).http://www.jstor.org/stable/j.ctt1jktrr4
2000
-
[7]
Coussy, O.: Thermoporoelasticity, pp. 71–112. John Wiley & Sons, Ltd, Chichester, UK (2003). Chap. 4. https://doi.org/10.1002/0470092718.ch4 . https://onlinelibrary.wiley.com/doi/abs/10.1002/0470092718.ch4
-
[8]
Lewis, R.W., Schrefler, B.A.: The Finite Element Method in the Static and Dynamic Deformation 15 Table A1: RelativeL 2 errors for the numerical experiments. Experiment Evaluation Rel.L 2 error V ariable Terzaghi’s problem (4.1) t= 0.005 s 1.2×10 −3 pressure t= 0.05 s 3.7×10 −4 pressure t= 0.2 s 4.7×10 −4 pressure t= 0.5 s 1.4×10 −3 pressure t= 0.005 s 9.0...
1998
-
[9]
Earth and Planetary Science Letters240(3), 539–558 (2005) https://doi.org/10.1016/j.epsl.2005.09.017
Steefel, C.I., DePaolo, D.J., Lichtner, P.C.: Reactive transport modeling: An essential tool and a new research approach for the earth sciences. Earth and Planetary Science Letters240(3), 539–558 (2005) https://doi.org/10.1016/j.epsl.2005.09.017
-
[10]
Hydrogeology Journal11(1), 7–40 (2003)
Rutqvist, J., Stephansson, O.: The role of hydromechanical coupling in fractured rock engineering. Hydrogeology Journal11(1), 7–40 (2003)
2003
-
[12]
SPE Journal 3(3), 219–226 (1998) https://doi.org/10.2118/50939-PA
Settari, A., Mourits, F.M.: A coupled reservoir and geomechanical simulation system. SPE Journal 3(3), 219–226 (1998) https://doi.org/10.2118/50939-PA
-
[13]
Convergence of iterative coupling for coupled flow and geome- chanics
Mikeli´ c, A., Wheeler, M.F.: Convergence of iterative coupling for coupled flow and geomechanics. Computational Geosciences17(3), 479–496 (2013) https://doi.org/10.1007/s10596-012-9318-y
-
[14]
Journal of Open Source Software 9(102), 6973 (2024) https://doi.org/10.21105/joss.06973
Settgast, R.R., Aronson, R.M., Besset, J.R., Borio, A., Bui, Q.M., Byer, T.J., Castelletto, N., Cit- rain, A., Corbett, B.C., Corbett, J., Cordier, P., Cremon, M.A., Crook, C.M., Cusini, M., Fei, F., Frambati, S., Franc, J., Franceschini, A., Frigo, M., Fu, P., Gazzola, T., Gross, H., Hamon, F., Han, B.M., Hao, Y., Hasanzade, R., Homel, M., Huang, J., Jin...
-
[15]
14: Kirsch benchmark results for experiment 4.5
Bilke, L., Naumov, D., Fischer, T., Wang, W., Mollaali, M., Buchwald, J., Max, J., Lehmann, C., 18 10−1 100 101 Radius (m) −1.6 −1.4 −1.2 −1.0 −0.8 −0.6 −0.4 −0.2 0.0 Stress (Pa) ×107 Hoop stress Radial stress Shear stress (a) 10−1 100 101 Radius (m) −2.0 −1.5 −1.0 −0.5 0.0 Stress (Pa) ×107 Hoop stress Radial stress Shear stress (b) 10−1 100 101 Radius (m...
2025
-
[16]
SoftwareX31, 102264 (2025) https://doi.org/10.1016/j.softx.2025.102264
Harbour, L., Giudicelli, G., Lindsay, A.D., Hansel, P.G., Icenhour, C., Li, M., Miller, J.M., Stogner, R.H., Behne, P., Yankura, D., Prince, Z.M., DeChant, C., Schwen, D., Spencer, B.W., Tano, M., Choi, N., Wang, Y., Nezdyur, M., Miao, Y., Hu, T., Kumar, S., Matthews, C., Langley, B., Nobre, N., Blair, A., MacMackin, C., Rocha, H.B., Palmer, E., Carter, J...
-
[17]
Hammond, G.E., Lichtner, P.C., Mills, R.T.: Evaluating the performance of parallel subsurface simulators: An illustrative example with pflotran. Water Resources Research50(1), 208–228 (2014) https://doi.org/10.1002/2012WR013483 https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1002/2012WR013483 19
-
[18]
Technical report, Los Alamos National Lab
Lichtner, P.C., Hammond, G.E., Lu, C., Karra, S., Bisht, G., Andre, B., Mills, R., Kumar, J.: Pflotran user manual: A massively parallel reactive flow and transport model for describing surface and subsurface processes. Technical report, Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Sandia National Lab. (SNL-NM), Albuquerque, NM (United...
-
[19]
Los Alamos, NM, Los Alamos National Laboratory (2001)
Lichtner, P.C., Hammond, G., Lu, C., Karra, S., Bisht, G., Andre, B., Mills, R., Kumar, J., Frederick, J.: Pflotran user manual. Los Alamos, NM, Los Alamos National Laboratory (2001)
2001
-
[20]
John Wiley & Sons, Chichester (2000)
Okabe, A., Boots, B., Sugihara, K., Chiu, S.N.: Spatial Tessellations: Concepts and Applications of Voronoi Diagrams, 2nd edn. John Wiley & Sons, Chichester (2000). https://doi.org/10.1002/ 9780470317013
2000
-
[21]
deuxi` eme m´ emoire: Recherches sur les parall´ elo` edres primitifs
Voronoi, G.: Nouvelles applications des param` etres continus ` a la th´ eorie des formes quadratiques. deuxi` eme m´ emoire: Recherches sur les parall´ elo` edres primitifs. Journal f¨ ur die reine und angewandte Mathematik134, 198–287 (1908)
1908
-
[22]
Bulletin de l’Acad´ emie des Sciences de l’URSS, Classe des sciences math´ ematiques et naturelles (6), 793–800 (1934)
Delaunay, B.: Sur la sph` ere vide. Bulletin de l’Acad´ emie des Sciences de l’URSS, Classe des sciences math´ ematiques et naturelles (6), 793–800 (1934)
1934
-
[23]
Cambridge University Press, ??? (2019)
Lie, K.-A.: An Introduction to Reservoir Simulation Using MATLAB/GNU Octave: User Guide for the MATLAB Reservoir Simulation Toolbox (MRST). Cambridge University Press, ??? (2019)
2019
-
[24]
Compu- tational Geosciences6(3-4), 405–432 (2002) https://doi.org/10.1023/A:1021291114475
Aavatsmark, I.: An introduction to multipoint flux approximations for quadrilateral grids. Compu- tational Geosciences6(3-4), 405–432 (2002) https://doi.org/10.1023/A:1021291114475
-
[25]
Palagi, C.L., Aziz, K.: Use of voronoi grid in reservoir simulation. SPE Advanced Technology Series2(02), 69–77 (1994) https://doi.org/10.2118/22889-PA https://onepetro.org/ATS/article- pdf/2/02/69/2606005/spe-22889-pa.pdf
-
[26]
Computational Geosciences2(4), 259–290 (1998) https://doi.org/ 10.1023/A:1011549407709
Edwards, M.G., Rogers, C.F.: Finite volume discretization with imposed flux continuity for the general tensor pressure equation. Computational Geosciences2(4), 259–290 (1998) https://doi.org/ 10.1023/A:1011549407709
-
[27]
In: Ciarlet, P.G., Lions, J.-L
Eymard, R., Gallou¨ et, T., Herbin, R.: Finite volume methods. In: Ciarlet, P.G., Lions, J.-L. (eds.) Handbook of Numerical Analysis vol. 7, pp. 713–1020. Elsevier, Amsterdam (2000)
2000
-
[28]
Kim, J., Tchelepi, H.A., Juanes, R.: Stability and convergence of sequential methods for coupled flow and geomechanics: Fixed-stress and fixed-strain splits. Computer Methods in Applied Mechanics and Engineering200(13-16), 1591–1606 (2011) https://doi.org/10.1016/j.cma.2010.12.022
-
[29]
http://www.mcs.anl.gov/petsc (2015)
Balay, S., Abhyankar, S., Adams, M.F., Brown, J., Brune, P., Buschelman, K., Dalcin, L., Eijkhout, V., Gropp, W.D., Kaushik, D., Knepley, M.G., McInnes, L.C., Rupp, K., Smith, B.F., Zampini, S., Zhang, H.: PETSc Web page. http://www.mcs.anl.gov/petsc (2015). http://www.mcs.anl.gov/ petsc
2015
-
[30]
Technical Report ANL-95/11 - Revision 3.6, Argonne National Laboratory (2015)
Balay, S., Abhyankar, S., Adams, M.F., Brown, J., Brune, P., Buschelman, K., Dalcin, L., Eijkhout, V., Gropp, W.D., Kaushik, D., Knepley, M.G., McInnes, L.C., Rupp, K., Smith, B.F., Zampini, S., Zhang, H.: PETSc users manual. Technical Report ANL-95/11 - Revision 3.6, Argonne National Laboratory (2015). http://www.mcs.anl.gov/petsc
2015
-
[31]
(eds.) Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries, pp
Balay, S., Gropp, W.D., McInnes, L.C., Smith, B.F.: In: Arge, E., Bruaset, A.M., Langtangen, H.P. (eds.) Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries, pp. 163–202. Birkh¨ auser Boston, Boston, MA (1997). https://doi.org/10.1007/978-1-4612-1986-68
-
[32]
John Wiley & Sons, New York (1943)
Terzaghi, K.: Theoretical Soil Mechanics. John Wiley & Sons, New York (1943)
1943
-
[33]
20 Highway Research Board Proceedings37, 584–617 (1958)
Schiffman, R.L.: Consolidation of soil under time-dependent loading and varying permeability. 20 Highway Research Board Proceedings37, 584–617 (1958)
1958
-
[34]
G´ eotechnique3(7), 287–299 (1953)
Mandel, J.: Consolidation des sols (´ etude math´ ematique). G´ eotechnique3(7), 287–299 (1953)
1953
-
[35]
Gongcheng Lixue/Engineering Mechanics22, 186–191 (2005)
Bai, B.: One-dimensional thermal consolidation characteristics of geotechnical media under non- isothermal condition. Gongcheng Lixue/Engineering Mechanics22, 186–191 (2005)
2005
-
[36]
Technical Report SAND2025-13033, Sandia National Laboratories (2025)
Sandia National Laboratories: Cubit mesh generation toolkit user documentation. Technical Report SAND2025-13033, Sandia National Laboratories (2025). https://cubit.sandia.gov
2025
-
[37]
https://doi.org/10.11578/dc.20221017.1
Stimpson, C., Sjaardema, G., Dudley, C., Yorgason, R., Hanks, B., White, D., Kerr, R., Mitchell, S., Shepherd, J., Tautges, T., et al.: Cubit v.16.X (2022). https://doi.org/10.11578/dc.20221017.1 . https://www.osti.gov/biblio/code-95109
-
[38]
https://lagrit.lanl.gov
Los Alamos National Laboratory: LaGriT: Los Alamos Grid Toolbox. https://lagrit.lanl.gov. Mesh generation and optimization software (2024)
2024
-
[39]
Karra, S.: satkarra/voronoi. GitHub. [Online; accessed 23-3-2026] (2024). https://github.com/ satkarra/voronoi
2026
-
[40]
https://github.com/nschloe/meshio
Schl¨ omer, N.: meshio: Python tools for I/O of mesh files. https://github.com/nschloe/meshio. Accessed: 2026-03-23 (2023)
2026
-
[41]
Computational Geosciences29, 25 (2025) https://doi.org/10.1007/s10596-025-10364-8
Huang, Y., Huang, X., Hyman, J.D.,et al.: Analytic solutions and field-scale application for verifi- cation of coupled thermo-hydro-mechanical processes in subsurface fractured media. Computational Geosciences29, 25 (2025) https://doi.org/10.1007/s10596-025-10364-8
-
[42]
Castelletto, N., White, J.A., Tchelepi, H.A.: Accuracy and convergence properties of the fixed- stress iterative solution of two-way coupled poromechanics. International Journal for Numerical and Analytical Methods in Geomechanics39(14), 1593–1618 (2015) https://doi.org/10.1002/nag.2400 https://onlinelibrary.wiley.com/doi/pdf/10.1002/nag.2400
-
[43]
G´ eotechnique46(2), 187–195 (1996) https://doi.org/10.1680/geot.1996.46.2.187
Abousleiman, Y., Cheng, A.H.-D., Cui, L., Detournay, E., Roegiers, J.-C.: Mandel’s problem revisited. G´ eotechnique46(2), 187–195 (1996) https://doi.org/10.1680/geot.1996.46.2.187
-
[44]
Rock Mechanics and Rock Engineering53(4) (2019) https://doi.org/10.1007/s00603-019-01989-0
Gao, Q., Ghassemi, A.: Three-dimensional thermo-poroelastic modeling and analysis of flow, heat transport and deformation in fractured rock with applications to a lab-scale geothermal system. Rock Mechanics and Rock Engineering53(4) (2019) https://doi.org/10.1007/s00603-019-01989-0
-
[45]
Pro- ceedings of the Royal Society of London
Sneddon, I.N.: The distribution of stress in the neighbourhood of a crack in an elastic solid. Pro- ceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences187(1009), 229–260 (1946) https://doi.org/10.1098/rspa.1946.0077
-
[46]
Quarterly of Applied Mathematics4, 262–267 (1946)
Sneddon, I.N., Elliot, H.A.: The opening of a griffith crack under internal pressure. Quarterly of Applied Mathematics4, 262–267 (1946)
1946
-
[47]
Inter- national Journal of Rock Mechanics and Mining Sciences44(5), 739–757 (2007) https://doi.org/10
Adachi, J., Siebrits, E., Peirce, A., Desroches, J.: Computer simulation of hydraulic fractures. Inter- national Journal of Rock Mechanics and Mining Sciences44(5), 739–757 (2007) https://doi.org/10. 1016/j.ijrmms.2006.11.006
2007
-
[48]
CRC Press, Boca Raton, FL (2017)
Anderson, T.L.: Fracture Mechanics: Fundamentals and Applications, 4th edn. CRC Press, Boca Raton, FL (2017). https://doi.org/10.1201/9781315370293
-
[49]
Zeitshrift des Vereines deutscher Ingenieure42, 797–807 (1898)
Kirsch, E.G.: Die theorie der elastizit¨ at und die bed¨ urfnisse der festigkeitslehre. Zeitshrift des Vereines deutscher Ingenieure42, 797–807 (1898)
-
[50]
Wiley- Blackwell, Maden, MA (2007)
Jaeger, J.C., Cook, N.G.W., Zimmerman, R.: Fundamentals of Rock Mechanics, 4th edn. Wiley- Blackwell, Maden, MA (2007). https://doi.org/10.1017/CBO9780511735349
-
[51]
Cambridge University Press, Cambridge (2007)
Zoback, M.D.: Reservoir Geomechanics. Cambridge University Press, Cambridge (2007). https:// doi.org/10.1017/CBO9780511586477 . Publication date: 10 December 2009 21
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.