Abstract
This paper aims to develop an accurate and efficient numerical model for three-dimensional transient thermal groundwater flow problems. The modified Richards equation and heat transport equation are considered to govern the thermal groundwater flows. For modeling water flows in the subsurface saturated-unsaturated porous media, we combined the method of characteristics (MOC) and a meshless localized radial basis function collocation method (LRBFCM). In order to implement the MOC scheme, the modified Richards equation is reformulated to an advection form and then computed by the particle tracking technique via MOC. We will then solve the heat equation and remaining terms of the temporal Richards equation by the LRBFCM. Seven benchmark subsurface flow problems with and without temperature effects are simulated and discussed to verify the feasibility and efficiency of this novel three-dimensional (3D) numerical model.
Get full access to this article
View all available purchase options and get full access to this article.
Data Availability Statement
Some or all data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request.
Acknowledgments
This project was financially supported by the Ministry of Sciences and Technology (MOST) of Taiwan under the Grant No. MOST 108-2221-E-002-007-MY2 and 107-2221-E-002-023-MY2, it is greatly appreciated.
References
Anibas, C., A. Tolche, G. Ghysels, J. Nossent, U. Schneidewind, M. Huysmans, and O. Batelaan. 2017. “Delineation of spatial-temporal patterns of groundwater/surface-water interaction along a river reach (Aa River, Belgium) with transient thermal modeling.” Hydrogeol. J. 26 (3): 819–835. https://doi.org/10.1007/s10040-017-1695-9.
Bear, J., and A. H.-D. Cheng. 2010. Modeling groundwater flow and contaminant transport. Berlin: Springer.
Bensabat, J., Q. Zhou, and J. Bear. 2000. “An adaptive pathline-based particle tracking algorithm for the Eulerian-Lagrangian method.” Adv. Water Resour. 23 (4): 383–397. https://doi.org/10.1016/S0309-1708(99)00025-1.
Blöcher, M. G., G. Zimmermann, I. Moeck, W. Brandt, A. Hassanzadegan, and F. Magri. 2010. “3D numerical modeling of hydrothermal processes during the lifetime of a deep geothermal reservoir.” Geofluids 10 (3): 406–421. https://doi.org/10.1111/j.1468-8123.2010.00284.x.
Canovas, M., I. Alhama, G. Garcia, E. Trigueros, and F. Alhama. 2017. “Numerical simulation of density-driven flow and heat transport processes in porous media using the network method.” Energies 10 (9): 1359. https://doi.org/10.3390/en10091359.
Canovas, M., I. Alhama, E. Trigueros, and F. Alhama. 2016. “A review of classical dimensionless numbers for the Yusa problem based on discriminated non-dimensionalization of the governing equations.” Hydrol. Processes 30 (22): 4101–4112. https://doi.org/10.1002/hyp.10878.
Chen, W., Z. J. Fu, and C. S. Chen. 2014. Recent advances in radial basis function collocation methods. Berlin: Springer.
Cheng, H. P., J. R. Cheng, and G. T. Yeh. 1996. “A particle tracking technique for the Lagrangian-Eulerian finite element method in multi-dimensions.” Int. J. Numer. Methods Eng. 39 (7): 1115–1136. https://doi.org/10.1002/(SICI)1097-0207(19960415)39:7%3C1115::AID-NME895%3E3.0.CO;2-4.
Hsiang, C. C., C. K. Chou, D. L. Young, J. Sladek, and V. Sladek. 2018. “Applying the method of characteristics and the meshless localized radial basis function collocation method to solve shallow water equations.” J. Eng. Mech. 144 (7): 04018047. https://doi.org/10.1061/(ASCE)EM.1943-7889.0001460.
Huang, K., R. Zhang, and M. T. van Genuchten. 1994. “An Eulerian-Lagrangian approach with an adaptively corrected method of characteristics to simulate variably saturated water flow.” Water Resour. Res. 30 (2): 499–507. https://doi.org/10.1029/93WR02881.
Konikow, L. F. 2011. “The secret to successful solute-transport modeling.” Ground Water 49 (2): 144–159. https://doi.org/10.1111/j.1745-6584.2010.00764.x.
Konikow, L. F., D. J. Goode, and G. Z. Hornberger. 1996. A three-dimensional method-of-characteristics solute-transport model (MOC3D). Washington, DC: US Geology Survey.
Kumar, S., A. Kumar, and D. Baleanu. 2016. “Two analytical methods for time-fractional nonlinear coupled Boussinesq–Burger’s equations arise in propagation of shallow water waves.” Nonlinear Dyn. 85 (2): 699–715. https://doi.org/10.1007/s11071-016-2716-2.
Kumar, S., and M. M. Rashidi. 2014. “New analytical method for gas dynamics equation arising in shock fronts.” Comput. Phys. Commun. 185 (7): 1947–1954. https://doi.org/10.1016/j.cpc.2014.03.025.
Lee, J.-Y., and G.-B. Kim. 2022. “An advanced mixed Lagrangian-Eulerian and finite element method to simulate 3-D subsurface variably saturated flows.” Geosci. J. 2022 (Jan): 21–39. https://doi.org/10.1007/s12303-021-0039-x.
Li, M.-H., H.-P. Cheng, and G.-T. Yeh. 2000. “Solving 3D subsurface flow and transport with adaptive multigrid.” J. Hydrol. Eng. 5 (1): 74–81. https://doi.org/10.1061/(ASCE)1084-0699(2000)5:1(74).
Limberger, J., T. Boxem, M. Pluymaekers, D. Bruhn, A. Manzella, P. Calcagno, F. Beekman, S. Cloetingh, and J.-D. van Weesa. 2018. “Geothermal energy in deep aquifers: A global assessment of the resource base for direct heat utilization.” Renewable Sustainable Energy Rev. 82 (Sep): 961–975. https://doi.org/10.1016/j.rser.2017.09.084.
Lin, H. C., D. R. Richards, G. T. Yeh, J. R. Cheng, H. P. Cheng, and N. L. Jones. 1997. FEMWATER: A three-dimensional finite element computer model for simulating density-dependent flow and transport in variably saturated media. Vicksburg, MS: US Army Corps of Engineers.
Lu, N. 1994. “A semianalytical method of path line computation for transient finite-difference groundwater flow models.” Water Resour. Res. 30 (8): 2449–2459. https://doi.org/10.1029/94WR01219.
Noorizadegan, A., D. L. Young, and C. S. Chen. 2021. “A novel local radial basis function collocation method for multi-dimensional piezoelectric problems.” J. Intell. Mater. Syst. Struct. 33 (12): 1574–1587. https://doi.org/10.1177/1045389X211057207.
Pollock, D. 1988. “Semianalytical computation of path lines for finite-difference models.” Ground Water 26 (6): 743–750. https://doi.org/10.1111/j.1745-6584.1988.tb00425.x.
Saatsaz, M., and S. Eslamian. 2020. “Groundwater modeling and its concepts, classifications, and applications for solute transport simulation in saturated porous media.” In Advances in hydrogeochemistry research. Vicksburg, MS: US Army Corps of Engineers.
Schafer-Perini, A., and J. Wilson. 1991. “Efficient and accurate front tracking for two-dimensional groundwater flow models.” Water Resour. Res. 27 (7): 1471–1485. https://doi.org/10.1029/91WR00720.
Suk H. J., G. T. Yeh, and J. S. Chen. 2021. “Achieving local mass conservation when using continuous Galerkin finite element methods to solve solute transport equations with spatially variable coefficients in a transient state.” J. Hydrol. 595 (1): 126038. https://doi:10.1016/j.jhydrol.2021.126038.
Vienken, T., M. Kreck, and P. Dietrich. 2019. “Monitoring the impact of intensive shallow geothermal energy use on groundwater temperatures in a residential neighborhood.” Geotherm. Energy 7 (8): 1–14. https://doi.org/10.1186/s40517-019-0123-x.
Yeh, G. T. 1999. Computational subsurface hydrology fluid flows. Norwell, MA: Kluwer Academic Publishers.
Yeh, G. T. 2000. Computational subsurface hydrology reactions. Norwell, MA: Kluwer Academic Publishers.
Yeh, G. T., H. P. Cheng, G. B. Huang, F. Zhang, H. C. Lin, E. Edris, and D. Richards. 2004a. A numerical model of flow, thermal transport, and salinity, sediment, and water quality transport in watershed systems of 1-D Stream-river network, 2-D Overland regime, and 3-D subsurface media (WASH123D: Version 2.0). Vicksburg, MS: US Army Corps of Engineers.
Yeh, G. T., Y. Fang, F. Zhang, J. Sun, Y. Li, M. H. Li, and M. D. Siegel. 2010. “Numerical modeling of coupled fluid flow and thermal and reactive biogeochemical transport in porous and fractured media.” Comput. Geosci. 14 (1): 149–170. https://doi.org/10.1007/s10596-009-9140-3.
Yeh, G. T., and H. Shan. 2008. “A mixed Lagrangian-Eulerian and finite element approach to modeling variably saturated flows in three dimensions.” In Proc., World Environmental and Water Resources Congress 2008. New York: Environmental and Water Resources Institute of American Society of Civil Engineers. https://doi.org/10.1061/40976(316)62.
Yeh, G. T., J. Sun, P. M. Jardine, W. D. Burgos, Y. Fang, M. H. Li, and M. D. Siegel. 2004b. HYDROGEOCHEM 5.0: A Three-Dimensional model of coupled fluid flow, thermal transport, and HYDROGEOCHEMICAL transport through variably saturated conditions. Oak Ridge, TN: Oak Ridge National Laboratory.
Yeh, G. T., and C. H. Tsai. 2014. “Mixed Lagrangian-Eulerian and Eulerian approach to discretizing Richards’ equation.” In Proc., American Geophysical Union, Fall Meeting. Washington, DC: American Geophysical Union.
Zanchini, E., S. Lazzari, and A. Priarone. 2012. “Long-term performance of large borehole heat exchanger fields with unbalanced seasonal loads and groundwater flow.” Energy 38 (Dec): 66–77. https://doi.org/10.1016/j.energy.2011.12.038.
Zhou, X., J. Zhao, J. Li, C. Ruan, L. Lin, and Q. Zhang. 2015. “Numerical modeling of 3D transient thermal groundwater flow and heat transport of geothermal fields of low to moderate temperature.” In Proc., World Geothermal Congress Melbourne, 19–25. Melbourne, Australia: International Geothermal Association.
Information & Authors
Information
Published In
Journal of Hydrologic Engineering
Volume 27 • Issue 12 • December 2022
Copyright
© 2022 American Society of Civil Engineers.
History
Received: Aug 6, 2021
Accepted: Jul 18, 2022
Published online: Sep 21, 2022
Published in print: Dec 1, 2022
Discussion open until: Feb 21, 2023
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.