TY - JOUR
T1 - Numerical issues arising in determination of interlayer conductivities in layered unsaturated soils
AU - Liu, Ruowen
AU - Welfert, Bruno
AU - Houston, Sandra
N1 - Funding Information:
This study was supported by the National Science Foundation (NSF) under Grant CMMI-0825089. The opinions, conclusions, and interpretations expressed in this paper are those of the authors, and not necessarily of the NSF. There are no datasharing issues, because all of the numerical simulations are provided in the figures produced by solving the 1D Richards equation. The initial setups of the numerical experiments are all stated in the paper and are also available from the authors on request.
Publisher Copyright:
© 2016 American Society of Civil Engineers.
Copyright:
Copyright 2020 Elsevier B.V., All rights reserved.
PY - 2017/3/1
Y1 - 2017/3/1
N2 - The one-dimensional (1D) Richards equation is used to model infiltration flow in an unsaturated two-layer soil. Imposing continuities of both flux and pressure head at the interface yields a nonlinear equation for determining interface conductivities. The authors show that multiple solutions of this nonlinear interface equation may exist if the spatial discretization is not fine enough around the interface, in particular as sharp wetting fronts pass through the interface, or for flow across highly dissimilar materials. Three hydraulic models, the Gardner model (G), the Mualem-van Genuchten model (MvG), and the Fredlund-Xing-Leong-Rahardjo model (FXLR), are investigated to demonstrate the nonuniqueness of solutions to the interface problem. For the simplest G model, a full mathematical analysis of the interface problem in terms of two parameters depending on local hydraulic conditions and mesh size is possible. For more advanced models, the interface equation can only be analyzed numerically. In all cases, the interface equation exhibits multiple solutions under suitable conditions, and the impact on the overall numerical solutions is studied. In particular, both a fixed-point iteration and a Newton iteration are used to solve the interface equation numerically.
AB - The one-dimensional (1D) Richards equation is used to model infiltration flow in an unsaturated two-layer soil. Imposing continuities of both flux and pressure head at the interface yields a nonlinear equation for determining interface conductivities. The authors show that multiple solutions of this nonlinear interface equation may exist if the spatial discretization is not fine enough around the interface, in particular as sharp wetting fronts pass through the interface, or for flow across highly dissimilar materials. Three hydraulic models, the Gardner model (G), the Mualem-van Genuchten model (MvG), and the Fredlund-Xing-Leong-Rahardjo model (FXLR), are investigated to demonstrate the nonuniqueness of solutions to the interface problem. For the simplest G model, a full mathematical analysis of the interface problem in terms of two parameters depending on local hydraulic conditions and mesh size is possible. For more advanced models, the interface equation can only be analyzed numerically. In all cases, the interface equation exhibits multiple solutions under suitable conditions, and the impact on the overall numerical solutions is studied. In particular, both a fixed-point iteration and a Newton iteration are used to solve the interface equation numerically.
KW - Interlayer conductivity
KW - Layered soil
KW - Richards equation
KW - Unsaturated flow
UR - http://www.scopus.com/inward/record.url?scp=85012961114&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85012961114&partnerID=8YFLogxK
U2 - 10.1061/(ASCE)GM.1943-5622.0000749
DO - 10.1061/(ASCE)GM.1943-5622.0000749
M3 - Article
AN - SCOPUS:85012961114
SN - 1532-3641
VL - 17
JO - International Journal of Geomechanics
JF - International Journal of Geomechanics
IS - 3
M1 - 04016078
ER -