3,635 research outputs found
Numerical solution of the Richards equation based catchment runoff model with dd-adaptivity algorithm and Boussinesq equation estimator
This paper presents a pseudo-deterministic catchment runoff model based on the Richards equation model - the governing equation for subsurface flow. The subsurface flow in a catchment is described here by two-dimensional variably saturated flow (unsaturated and saturated). The governing equation is the Richards equation with a slight modification of the time derivative term, as considered e.g. by Neuman. The nonlinear nature of this problem appears in the unsaturated zone only, so it was possible to make use of adaptive domain decomposition algorithm. However delineating of the saturated zone boundary is a nonlinear computationally expensive issue. The simple one-dimensional Boussinesq equation was used here as a rough estimator of the saturated zone boundary. With this estimate the adaptive domain decomposition could always start with an optimal subdomain split, and thus it is now possible to avoid solving huge systems of linear equations in the initial iteration level.
With this measure it is possible to construct an efficient two-dimensional pseudodeterministic catchment runoff model. Finally, the model is tested against real data originating from the Modrava 2 experimental catchment, Czech Republic
A linear domain decomposition method for partially saturated flow in porous media
The Richards equation is a nonlinear parabolic equation that is commonly used
for modelling saturated/unsaturated flow in porous media. We assume that the
medium occupies a bounded Lipschitz domain partitioned into two disjoint
subdomains separated by a fixed interface . This leads to two problems
defined on the subdomains which are coupled through conditions expressing flux
and pressure continuity at . After an Euler implicit discretisation of
the resulting nonlinear subproblems a linear iterative (-type) domain
decomposition scheme is proposed. The convergence of the scheme is proved
rigorously. In the last part we present numerical results that are in line with
the theoretical finding, in particular the unconditional convergence of the
scheme. We further compare the scheme to other approaches not making use of a
domain decomposition. Namely, we compare to a Newton and a Picard scheme. We
show that the proposed scheme is more stable than the Newton scheme while
remaining comparable in computational time, even if no parallelisation is being
adopted. Finally we present a parametric study that can be used to optimize the
proposed scheme.Comment: 34 pages, 13 figures, 7 table
A localized orthogonal decomposition method for semi-linear elliptic problems
In this paper we propose and analyze a new Multiscale Method for solving
semi-linear elliptic problems with heterogeneous and highly variable
coefficient functions. For this purpose we construct a generalized finite
element basis that spans a low dimensional multiscale space. The basis is
assembled by performing localized linear fine-scale computations in small
patches that have a diameter of order H |log H| where H is the coarse mesh
size. Without any assumptions on the type of the oscillations in the
coefficients, we give a rigorous proof for a linear convergence of the H1-error
with respect to the coarse mesh size. To solve the arising equations, we
propose an algorithm that is based on a damped Newton scheme in the multiscale
space
Nonlinear multigrid based on local spectral coarsening for heterogeneous diffusion problems
This work develops a nonlinear multigrid method for diffusion problems
discretized by cell-centered finite volume methods on general unstructured
grids. The multigrid hierarchy is constructed algebraically using aggregation
of degrees of freedom and spectral decomposition of reference linear operators
associated with the aggregates. For rapid convergence, it is important that the
resulting coarse spaces have good approximation properties. In our approach,
the approximation quality can be directly improved by including more spectral
degrees of freedom in the coarsening process. Further, by exploiting local
coarsening and a piecewise-constant approximation when evaluating the nonlinear
component, the coarse level problems are assembled and solved without ever
re-visiting the fine level, an essential element for multigrid algorithms to
achieve optimal scalability. Numerical examples comparing relative performance
of the proposed nonlinear multigrid solvers with standard single-level
approaches -- Picard's and Newton's methods -- are presented. Results show that
the proposed solver consistently outperforms the single-level methods, both in
efficiency and robustness
A reduced basis localized orthogonal decomposition
In this work we combine the framework of the Reduced Basis method (RB) with
the framework of the Localized Orthogonal Decomposition (LOD) in order to solve
parametrized elliptic multiscale problems. The idea of the LOD is to split a
high dimensional Finite Element space into a low dimensional space with
comparably good approximation properties and a remainder space with negligible
information. The low dimensional space is spanned by locally supported basis
functions associated with the node of a coarse mesh obtained by solving
decoupled local problems. However, for parameter dependent multiscale problems,
the local basis has to be computed repeatedly for each choice of the parameter.
To overcome this issue, we propose an RB approach to compute in an "offline"
stage LOD for suitable representative parameters. The online solution of the
multiscale problems can then be obtained in a coarse space (thanks to the LOD
decomposition) and for an arbitrary value of the parameters (thanks to a
suitable "interpolation" of the selected RB). The online RB-LOD has a basis
with local support and leads to sparse systems. Applications of the strategy to
both linear and nonlinear problems are given
NUMERICAL SIMULATION OF A SINGLE RING INFILTRATION EXPERIMENT WITH hp-ADAPTIVE SPACE-TIME DISCONTINUOUS GALERKIN METHOD
We present a novel hp-adaptive space-time discontinuous Galerkin (hp-STDG) method for the numerical solution of the nonstationary Richards equation equipped with Dirichlet, Neumann and seepage face boundary conditions. The hp-STDG method presented in this paper is a generalization of a hp-STDG method which was developed for time dependent non-linear convective-diffusive problems. We describe the method and the single ring experiment, and then we present a numerical experiment which clearly demonstrates the superiority of the hp-STDG method over a discontinuous Galerkin method based on a static fine mesh
A study on iterative methods for solving Richards` equation
This work concerns linearization methods for efficiently solving the
Richards` equation,a degenerate elliptic-parabolic equation which models flow
in saturated/unsaturated porous media.The discretization of Richards` equation
is based on backward Euler in time and Galerkin finite el-ements in space. The
most valuable linearization schemes for Richards` equation, i.e. the
Newtonmethod, the Picard method, the Picard/Newton method and theLscheme are
presented and theirperformance is comparatively studied. The convergence, the
computational time and the conditionnumbers for the underlying linear systems
are recorded. The convergence of theLscheme is theo-retically proved and the
convergence of the other methods is discussed. A new scheme is
proposed,theLscheme/Newton method which is more robust and quadratically
convergent. The linearizationmethods are tested on illustrative numerical
examples
- …