Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEfficient Solvers for Nonstandard Models for Flow and Transport in Unsaturated Porous Media
380
0
0.0
(
0
)
Added by
Davide Illiano
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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 to incorporate nonstandard effects like dynamic capillarity and hysteresis, and a reactive transport equation for the solute. The two model components are strongly coupled. On one hand, the flow affects the concentration of the solute; on the other hand, the surface tension is a function of the solute, which impacts the capillary pressure and, consequently, the flow. After applying an Euler implicit scheme, we consider a set of iterative linearization schemes to solve the resulting nonlinear equations, including both monolithic and two splitting strategies. The latter include a canonical nonlinear splitting and an alternate linearized splitting, which appears to be overall faster in terms of numbers of iterations, based on our numerical studies. The (time discrete) system being nonlinear, we investigate different linearization methods. We consider the linearly convergent L-scheme, which converges unconditionally, and the Newton method, converging quadratically but subject to restrictions on the initial guess. Whenever hysteresis effects are included, the Newton method fails to converge. The L-scheme converges; nevertheless, it may require many iterations. This aspect is improved by using the Anderson acceleration. A thorough comparison of the different solving strategies is presented in five numerical examples, implemented in MRST, a toolbox based on MATLAB.
rate research
Read More
In this article, we present new random walk methods to solve flow and transport problems in unsaturated/saturated porous media, including coupled flow and transport processes in soils, heterogeneous systems modeled through random hydraulic conductivity and recharge fields, processes at the field and regional scales. The numerical schemes are based on global random walk algorithms (GRW) which approximate the solution by moving large numbers of computational particles on regular lattices according to specific random walk rules. To cope with the nonlinearity and the degeneracy of the Richards equation and of the coupled system, we implemented the GRW algorithms by employing linearization techniques similar to the $L$-scheme developed in finite element/volume approaches. The resulting GRW $L$-schemes converge with the number of iterations and provide numerical solutions that are first-order accurate in time and second-order in space. A remarkable property of the flow and transport GRW solutions is that they are practically free of numerical diffusion. The GRW solutions are validated by comparisons with mixed finite element and finite volume solutions in one- and two-dimensional benchmark problems. They include Richards equation fully coupled with the advection-diffusion-reaction equation and capture the transition from unsaturated to saturated flow regimes. For completeness, we also consider decoupled flow and transport model problems for saturated aquifers.
Flow and multicomponent reactive transport in saturated/unsaturated porous media are modeled by ensembles of computational particles moving on regular lattices according to specific random walk rules. The occupation number of the lattice sites is updated with a global random walk (GRW) procedure which simulates the evolution of the ensemble with computational costs comparable to those for a single random walk simulation in sequential procedures. To cope with the nonlinearity and the degeneracy of the Richards equation the GRW flow solver uses linearization techniques similar to the $L$-scheme developed in finite element/volume approaches. Numerical schemes for reactive transport, coupled with the flow solver via numerical solutions for saturation and water flux, are implemented in splitting procedures. Diffusion-advection steps are solved by GRW algorithms using either biased or unbiased random walk probabilities. Since the number of particles in GRW simulations can be as large as the number of molecules involved in chemical reactions, one avoids the cumbersome problem of rescaling particle densities to approximate concentrations. Reaction steps are therefore formulated in terms of concentrations, as in deterministic approaches. The numerical convergence of the new schemes is demonstrated by comparisons with manufactured analytical solutions. Coupled flow and reactive transport problems of contaminant biodegradation described by the Monod model are further solved and the influence of flow nonlinearity/degeneracy and of the spatial heterogeneity of the medium is investigated numerically.
In this work, we consider a mathematical model for flow in a unsaturated porous medium containing a fracture. In all subdomains (the fracture and the adjacent matrix blocks) the flow is governed by Richards equation. The submodels are coupled by physical transmission conditions expressing the continuity of the normal fluxes and of the pressures. We start by analyzing the case of a fracture having a fixed width-length ratio, called $varepsilon > 0$. Then we take the limit $varepsilon to 0$ and give a rigorous proof for the convergence towards effective models. This is done in different regimes, depending on how the ratio of porosities and permeabilities in the fracture, respectively matrix scale with respect to $varepsilon$, and leads to a variety of effective models. Numerical simulations confirm the theoretical upscaling results.
In this work we consider the transport of a surfactant in a variably saturated porous media. The water flow is modelled by the Richards equations and it is fully coupled with the transport equation for the surfactant. Three linearization techniques are discussed: the Newton method, the modified Picard and the L-scheme. Based on these, monolithic and splitting schemes are proposed and their convergence is analyzed. The performance of these schemes is illustrated on four numerical examples. For these examples, the number of iterations and the condition numbers of the linear systems emerging in each iteration are presented.
In this paper, we study a model for the transport of an external component, e.g., a surfactant, in variably saturated porous media. We discretize the model in time and space by combining a backward Euler method with the linear Galerkin finite elements. The Newton method and the L-Scheme are employed for the linearization and the performance of these schemes is studied numerically. A special focus is set on the effects of dynamic capillarity on the transport equation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
