Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

Zhao, CB, Hobbs, BE, Walshe, JL, Muhlhaus, HB and Ord, A (2001) Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. Computer Methods In Applied Mechanics And Engineering, 190 18-19: 2277-2293. doi:10.1016/S0045-7825(00)00304-2


Author Zhao, CB
Hobbs, BE
Walshe, JL
Muhlhaus, HB
Ord, A
Title Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins
Journal name Computer Methods In Applied Mechanics And Engineering   Check publisher's open access policy
ISSN 0045-7825
Publication date 2001-01-01
Sub-type Article (original research)
DOI 10.1016/S0045-7825(00)00304-2
Open Access Status DOI
Volume 190
Issue 18-19
Start page 2277
End page 2293
Total pages 17
Place of publication Lausanne
Publisher Elsevier Science Sa
Language eng
Abstract In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.
Keyword Mathematics, Interdisciplinary Applications
Engineering, Multidisciplinary
Mechanics
Fluid-rock Interaction
Generalized Concentration Concept
Term Splitting Algorithm
Equivalent Source/sink Term
Finite Element Modeling
Mass-transport Problems
Porous-media
Heat-transfer
Systems
Mineralization
Flow
Throughflow
Gradient
Q-Index Code C1
Q-Index Status Provisional Code
Institutional Status Unknown

Document type: Journal Article
Sub-type: Article (original research)
Collection: Earth Systems Science Computational Centre Publications
 
Versions
Version Filter Type
Citation counts: TR Web of Science Citation Count  Cited 31 times in Thomson Reuters Web of Science Article | Citations
Scopus Citation Count Cited 39 times in Scopus Article | Citations
Google Scholar Search Google Scholar
Created: Mon, 13 Aug 2007, 22:13:53 EST