Numerical simulation of Richards equation in partially saturated porous media: under-relaxation and mass balance

Abstract

Numerical simulation of Richards equation in unsaturated soil is known to be difficult because of the highly non-linear material properties involved. An extensive study was conducted to clarify the role of mass balance and under-relaxation on the rate in which the correct solution is approached. Both finite element and finite difference techniques were considered. For the former, the hbased formulation was compared with a massconservative mixed form. The conductivity function (K) was under-relaxed in two ways while the capacity function (C) is computed following the standard mass or non-mass conservative schemes recommended in literature. For fairly coarse discretisation, it was found that large errors were produced when K was under-relaxed by evaluating it at the average of heads from current and previous time step (UR1), regardless of the numerical scheme used. Maintaining global mass balance is found to have little impact on the accuracy. All numerical schemes that underrelaxed K by computing it at the average of two most recent iterations in the current time step (UR2) converged quicker to the correct solution with increasing discretisation, although more iterations per time step than UR1 is needed to achieve a stable solution. An important practical ramification is that it appears to be possible to achieve reasonably accurate and oscillation-free results using fairly coarse discretisation by making only minor modifications (namely, using UR2 for conductivity function) to the h-based finite element formulation and applying some minimum time step criteria.

Numerical modelling on fate and transport of petroleum hydrocarbons in an unsaturated subsurface system for varying source scenario

“…The porous media hydraulic properties, initial and boundary conditions for this validation test are shown in table 4. Figure 3 provides the comparison of the results obtained from present flow model with that of Phoon et al (2007). It is observed that the model predictions are in good agreement with the reported values of pressure head.…”

“…To examine the performance of the flow model (equation 1), numerical simulations were carried out using data from the one-dimensional numerical experiment of Phoon et al (2007), where water infiltrates into a homogeneous soil column which is initially dry. The porous media hydraulic properties, initial and boundary conditions for this validation test are shown in table 4.…”

Numerical modeling of biological clogging on transport of nitrate in an unsaturated porous media

“…Since the analytical solutions are not available for the present study, the numerical model has been validated separately for water flow in unsaturated zone and bioclogging in saturated zone with their appropriate published analytical/numerical solution. The numerical solution for water flow in an unsaturated zone has been provided by Mitchell and Mayer (1998), Lehmann and Ackerer (1998) and Phoon et al (2007). It is observed from Fig.…”

Avaliação dos algoritmos de Picard-Krylov na solução da equação de Richards

“…Para validar o código computacional foram resolvidos alguns problemas clássicos da literatura, a saber: Haverkamp et al [11], Phoon et al [7], Vasconcellos [2], Celia et al [8] e Juncu et al [5]. Igualmente foram utilizadas as soluções semi-analíticas de Philip [6], Warrick's et al…”