Related papers: A finite element method using a bounded auxiliary …
We consider here a cell-centered finite difference approximation of the Richards equation in three dimensions, averaging for interface values the hydraulic conductivity $K=K(p)$, a highly nonlinear function, by arithmetic, upstream, and…
This paper proposes two algorithms to impose seepage boundary conditions in the context of Richards' equation for groundwater flows in unsaturated media. Seepage conditions are non-linear boundary conditions, that can be formulated as a set…
We consider implementational aspects of the mixed finite element method for a special class of nonlinear problems. We establish the equivalence of the hybridized formulation of the mixed finite element method to a nonconforming finite…
The Richards' equation is a model for flow of water in unsaturated soils. The coefficients of this (nonlinear) partial differential equation describe the permeability of the medium. Insufficient or uncertain measurements are commonly…
Traditional finite element approaches are well-known to introduce spurious oscillations when applied to advection-dominated problems. We explore alleviation of this issue from the perspective of a generalized finite element formulation,…
This paper presents a mathematical analysis of a doubly degenerate parabolic equation and its application to the Richards equation using a bounded auxiliary variable. We establish the existence of weak solutions using semi-implicit time…
We present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that are modeled using the Richards' equation. We propose a stochastic extension for the empirical models that are…
Simulating infiltration in porous media using Richards' equation remains computationally challenging due to its parabolic structure and nonlinear coefficients. While a wide range of numerical methods for differential equations have been…
We study several iterative methods for fully coupled flow and reactive transport in porous media. The resulting mathematical model is a coupled, nonlinear evolution system. The flow model component builds on the Richards equation, modified…
A new modified Galerkin / Finite Element Method is proposed for the numerical solution of the fully nonlinear shallow water wave equations. The new numerical method allows the use of low-order Lagrange finite element spaces, despite the…
Recently, the nonlinearity continuation method has been used to numerically solve boundary value problems for steady-state Richards equation. The method can be considered as a predictor-corrector procedure with the simplest form which has…
In this note, the existing porousMultiphaseFoam toolbox, developed initially for any two-phase flow in porous media is extended to the specific case of the Richards' equation which neglect the pressure gradient of the non-wetting phase.…
This work proposes a nonlinear finite element method whose nodal values preserve bounds known for the exact solution. The discrete problem involves a nonlinear projection operator mapping arbitrary nodal values into bound-preserving ones…
We propose a new numerical method to solve the Cahn-Hilliard equation coupled with non-linear wetting boundary conditions. We show that the method is mass-conservative and that the discrete solution satisfies a discrete energy law similar…
In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation…
We present a continuous finite element method for some examples of fully nonlinear elliptic equation. A key tool is the discretisation proposed in Lakkis & Pryer (2011, SISC) allowing us to work directly on the strong form of a linear PDE.…
In this paper, we focus on one-dimensional vertical infiltration, assuming constant diffusivity and a quadratic relationship between hydraulic conductivity and water content. Under these assumptions, Richards' equation reduces to Burgers'…
In this paper, we first propose a coupled numerical model of unsaturated flow in soils and plant root water uptake. The Richards equation and different formulations are used in the developed numerical model to describe infiltration in root…
Modeling flow in geosystems with natural fault is a challenging problem due to low permeability of fault compared to its surrounding porous media. One way to predict the behavior of the flow while taking the effects of fault into account is…
We construct a finite element method for the numerical solution of a fractional porous medium equation on a bounded open Lipschitz polytopal domain $\Omega \subset \mathbb{R}^{d}$, where $d = 2$ or $3$. The pressure in the model is defined…