Related papers: Generalized Green-Ampt Approach to 1D Oscillatory …
Richards equation is often used to represent two-phase fluid flow in an unsaturated porous medium when one phase is much heavier and more viscous than the other. However, it cannot describe the fully saturated flow for some capillary…
In this study, a novel semi-implicit second-order temporal scheme combined with the finite element method for space discretization is proposed to solve the coupled system of infiltration and solute transport in unsaturated porous media. The…
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…
Analysing two-dimensional shallow water equations with idealised bottom topographies have many applications in the atmospheric and oceanic sciences; however, restrictive flow pattern assumptions have been made to achieve explicit solutions.…
In this paper, we present a multiscale method for simulations of the multicontinua unsaturated flow problems in heterogeneous fractured porous media. The mathematical model is described by the system of Richards equations for each continuum…
The derivation of a quasi-geostrophic (QG) system from the rotating shallow water equations on a midlatitude beta-plane coupled with moisture is presented. Condensation is prescribed to occur whenever the moisture at a point exceeds a…
In this work, we study the behaviour of saturation fronts for two phase flow through a long homogeneous porous column. In particular, the model includes hysteresis and dynamic effects in the capillary pressure and hysteresis in the…
This paper explores static non-uniform grid solvers that adapt three raster-based flood models on an optimised non-uniform grid: the second-order discontinuous Galerkin (DG2) model representing the modelled data as piecewise-planar fields,…
In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced-order methods are demanded. Combining these…
Multi-component fluid flow simulations in multi-scale porous structures often involve regions that are under-resolved at practical computational resolutions. Accurately capturing the contributions from these unresolved regions is critical.…
Modeling unsaturated flow through soils with water uptake by plan root has many applications in agriculture and water resources management. In this study, our aim is to develop efficient numerical techniques for solving the Richards…
The Richards equation, a nonlinear elliptic parabolic equation, is widely used to model infiltration in porous media. We develop a finite element method for solving the Richards equation by introducing a new bounded auxiliary variable to…
We propose an analytical framework to model the effect of single and multiple mechanical surface oscillators on the dynamics of vertically polarized elastic waves propagating in a semi-infinite medium. The formulation extends the canonical…
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…
Simulations of immiscible flows involving surface tension (ST) require a robust high-fidelity framework. State-of-the-art multi-phase models, such as the Conservative Level-Set (CLS) approach, rely on Eulerian representations of the fluids…
We propose a new model to approximate the wave response of waveguides containing an arbitrary number of small inclusions. The theory is developed to consider any one-dimensional waveguide (longitudinal, flexural, shear, torsional waves or a…
Analysing the impact of bottom friction on shallow water waves over bottom terrains is important in areas including environmental and coastal engineering as well as the oceanic and atmospheric sciences. However, current theoretical…
The unsteady response of a water free surface to a localized pressure source moving at constant speed $U$ in the range $0.95c_\mathrm{min} \lesssim U \leq 1.02 c_\mathrm{min}$, where $c_\mathrm{min}$ is the minimum phase speed of linear…
In this paper we propose a method for the construction of locally conservative flux fields from Generalized Multiscale Finite Element Method (GMsFEM) pressure solutions. The flux values are obtained from an element-based postprocessing…
In this paper, we combine discrete empirical interpolation techniques, global mode decomposition methods, and local multiscale methods, such as the Generalized Multiscale Finite Element Method (GMsFEM), to reduce the computational…