A Fully Discrete, Mass-Conservative Adaptive Discontinuous Spectral-Element Method for the One-Dimensional Richards Equation in Mixed Form
Abstract
We present and analyze a spatially and temporally adaptive high-order method for the one-dimensional Richards equation describing vertical water transport in unsaturated soil, with an emphasis on resolving sharp wetting fronts while satisfying a fully discrete water mass balance. The spatial discretization uses an arbitrary-order local discontinuous Galerkin spectral-element formulation with collocated Legendre-Gauss-Lobatto quadrature, and the time integration is based on a high-order Rosenbrock-Wanner method. We focus primarily on the mixed formulation of the Richards equation, which advances water content as the conserved variable while evaluating the nonlinear diffusion operator in terms of pressure head; the resulting fully discrete scheme satisfies the water mass balance in exact arithmetic, independent of spatial resolution and time step size, including under spatial and temporal adaptivity. In contrast, the pressure-head formulation advances pressure head through a nonlinear capacity function and satisfies the corresponding conservation identity only at the semi-discrete level. Numerical experiments examine convergence through spatial and temporal refinement studies, with mass conservation assessed under spatial and temporal adaptivity using the solver flux output method (SFOM). The SFOM results confirm machine-precision mass conservation for the mixed formulation, whereas the pressure-head formulation exhibits mass balance errors above roundoff. Moreover, in the temporal refinement studies, the mixed formulation yields more accurate solutions than the pressure-head formulation for a prescribed time step size. The combination of high-order accuracy, spatial and temporal adaptivity, and fully discrete mass conservation makes the proposed mixed-form discontinuous spectral-element method a promising approach for simulating unsaturated flow in next-generation land models.
Related
- A well-balanced path-conservative discontinuous Galerkin method with generalized hydrostatic reconstruction for non-conservative balance laws
- High-order in-cell discontinuous reconstruction path-conservative methods for conservative systems in primitive variables
- Resolvent Expansions at All Spectral Thresholds for Three-Dimensional Discrete Schrödinger Operators
- GPU Performance of an Entropy-Stable Discontinuous Galerkin Euler Solver with Non-Conservative Terms
- Discrete hypocoercive estimates for discontinuous Galerkin methods: application to the Vlasov-Poisson-Fokker-Planck system
BibTeX
@misc{montoya2026fully,
title = {{A Fully Discrete, Mass-Conservative Adaptive Discontinuous Spectral-Element Method for the One-Dimensional Richards Equation in Mixed Form}},
author = {Tristan Montoya and Raymond J. Spiteri},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.17031},
eprint = {2610.17031},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.17031},
}
RIS
TY - UNPB
TI - A Fully Discrete, Mass-Conservative Adaptive Discontinuous Spectral-Element Method for the One-Dimensional Richards Equation in Mixed Form
AU - Tristan Montoya
AU - Raymond J. Spiteri
PY - 2026
DA - 2026/10/12
PB - arXiv
DO - 10.48550/arxiv.2610.17031
UR - https://arxiv.org/abs/2610.17031
AB - We present and analyze a spatially and temporally adaptive high-order method for the one-dimensional Richards equation describing vertical water transport in unsaturated soil, with an emphasis on resolving sharp wetting fronts while satisfying a fully discrete water mass balance. The spatial discretization uses an arbitrary-order local discontinuous Galerkin spectral-element formulation with collocated Legendre-Gauss-Lobatto quadrature, and the time integration is based on a high-order Rosenbrock-Wanner method. We focus primarily on the mixed formulation of the Richards equation, which advances water content as the conserved variable while evaluating the nonlinear diffusion operator in terms of pressure head; the resulting fully discrete scheme satisfies the water mass balance in exact arithmetic, independent of spatial resolution and time step size, including under spatial and temporal adaptivity. In contrast, the pressure-head formulation advances pressure head through a nonlinear capacity function and satisfies the corresponding conservation identity only at the semi-discrete level. Numerical experiments examine convergence through spatial and temporal refinement studies, with mass conservation assessed under spatial and temporal adaptivity using the solver flux output method (SFOM). The SFOM results confirm machine-precision mass conservation for the mixed formulation, whereas the pressure-head formulation exhibits mass balance errors above roundoff. Moreover, in the temporal refinement studies, the mixed formulation yields more accurate solutions than the pressure-head formulation for a prescribed time step size. The combination of high-order accuracy, spatial and temporal adaptivity, and fully discrete mass conservation makes the proposed mixed-form discontinuous spectral-element method a promising approach for simulating unsaturated flow in next-generation land models.
ER -