TY - JOUR

T1 - A semianalytical approach for solving first-order perturbation equations of dissolution-timescale reactive infiltration instability problems in fluid-saturated rocks

AU - Zhao, Chongbin

AU - Hobbs, Bruce E.

AU - Ord, Alison

PY - 2020/10/25

Y1 - 2020/10/25

N2 - This paper presents a semianalytical approach for solving first-order perturbation (FOP) equations, which are used to describe dissolution-timescale reactive infiltration instability (RII) problems in fluid-saturated rocks. The proposed approach contains two parts because the chemical dissolution reaction divides the whole problem domain into two subdomains. In the first part, the interface-condition substitution strategy is used to derive the analytical expressions of purely mathematical solutions for the FOP equations in the upstream subdomain, where the dissolution chemical reaction is ceased and the FOP equations are weakly coupled. In the second part, the finite element method (FEM) is used to derive the analytical expressions of numerical solutions for the FOP equations in the downstream subdomain, where the dissolution chemical reaction needs to be considered and the FOP equations are strongly coupled so that it is impossible to derive purely mathematical solutions for them. Particular attention is paid to the development of the element-by-element forward marching strategy, which is associated with the use of the FEM for solving this new kind of scientific problem. The related analytical results demonstrated that (1) both the dynamic characteristic of a reactive infiltration system and the dimensionless wavenumber can have pronounced influences on the distribution of the FOP dimensionless acid concentration within the entire domain of the dissolution-timescale RII problems in fluid-saturated rocks and (2) the FOP dimensionless acid concentration distribution exhibits two significantly different patterns in the upstream and downstream subdomains of the dissolution-timescale RII system.

AB - This paper presents a semianalytical approach for solving first-order perturbation (FOP) equations, which are used to describe dissolution-timescale reactive infiltration instability (RII) problems in fluid-saturated rocks. The proposed approach contains two parts because the chemical dissolution reaction divides the whole problem domain into two subdomains. In the first part, the interface-condition substitution strategy is used to derive the analytical expressions of purely mathematical solutions for the FOP equations in the upstream subdomain, where the dissolution chemical reaction is ceased and the FOP equations are weakly coupled. In the second part, the finite element method (FEM) is used to derive the analytical expressions of numerical solutions for the FOP equations in the downstream subdomain, where the dissolution chemical reaction needs to be considered and the FOP equations are strongly coupled so that it is impossible to derive purely mathematical solutions for them. Particular attention is paid to the development of the element-by-element forward marching strategy, which is associated with the use of the FEM for solving this new kind of scientific problem. The related analytical results demonstrated that (1) both the dynamic characteristic of a reactive infiltration system and the dimensionless wavenumber can have pronounced influences on the distribution of the FOP dimensionless acid concentration within the entire domain of the dissolution-timescale RII problems in fluid-saturated rocks and (2) the FOP dimensionless acid concentration distribution exhibits two significantly different patterns in the upstream and downstream subdomains of the dissolution-timescale RII system.

KW - chemical dissolution

KW - dissolution timescale

KW - forward marching strategy

KW - perturbation equations

KW - porous rocks

KW - reactive infiltration instability

UR - http://www.scopus.com/inward/record.url?scp=85088580559&partnerID=8YFLogxK

U2 - 10.1002/nag.3119

DO - 10.1002/nag.3119

M3 - Article

AN - SCOPUS:85088580559

VL - 44

SP - 2070

EP - 2092

JO - International Journal of Numerical and Analytical Methods in Geomechanics

JF - International Journal of Numerical and Analytical Methods in Geomechanics

SN - 0363-9061

IS - 15

ER -