70,752 research outputs found
Krylov implicit integration factor discontinuous Galerkin methods on sparse grids for high dimensional reaction-diffusion equations
Computational costs of numerically solving multidimensional partial
differential equations (PDEs) increase significantly when the spatial
dimensions of the PDEs are high, due to large number of spatial grid points.
For multidimensional reaction-diffusion equations, stiffness of the system
provides additional challenges for achieving efficient numerical simulations.
In this paper, we propose a class of Krylov implicit integration factor (IIF)
discontinuous Galerkin (DG) methods on sparse grids to solve reaction-diffusion
equations on high spatial dimensions. The key ingredient of spatial DG
discretization is the multiwavelet bases on nested sparse grids, which can
significantly reduce the numbers of degrees of freedom. To deal with the
stiffness of the DG spatial operator in discretizing reaction-diffusion
equations, we apply the efficient IIF time discretization methods, which are a
class of exponential integrators. Krylov subspace approximations are used to
evaluate the large size matrix exponentials resulting from IIF schemes for
solving PDEs on high spatial dimensions. Stability and error analysis for the
semi-discrete scheme are performed. Numerical examples of both scalar equations
and systems in two and three spatial dimensions are provided to demonstrate the
accuracy and efficiency of the methods. The stiffness of the reaction-diffusion
equations is resolved well and large time step size computations are obtained
Optimal Local Multi-scale Basis Functions for Linear Elliptic Equations with Rough Coefficient
This paper addresses a multi-scale finite element method for second order
linear elliptic equations with arbitrarily rough coefficient. We propose a
local oversampling method to construct basis functions that have optimal local
approximation property. Our methodology is based on the compactness of the
solution operator restricted on local regions of the spatial domain, and does
not depend on any scale-separation or periodicity assumption of the
coefficient. We focus on a special type of basis functions that are harmonic on
each element and have optimal approximation property. We first reduce our
problem to approximating the trace of the solution space on each edge of the
underlying mesh, and then achieve this goal through the singular value
decomposition of an oversampling operator. Rigorous error estimates can be
obtained through thresholding in constructing the basis functions. Numerical
results for several problems with multiple spatial scales and high contrast
inclusions are presented to demonstrate the compactness of the local solution
space and the capacity of our method in identifying and exploiting this compact
structure to achieve computational savings
Comparison of data-driven uncertainty quantification methods for a carbon dioxide storage benchmark scenario
A variety of methods is available to quantify uncertainties arising with\-in
the modeling of flow and transport in carbon dioxide storage, but there is a
lack of thorough comparisons. Usually, raw data from such storage sites can
hardly be described by theoretical statistical distributions since only very
limited data is available. Hence, exact information on distribution shapes for
all uncertain parameters is very rare in realistic applications. We discuss and
compare four different methods tested for data-driven uncertainty
quantification based on a benchmark scenario of carbon dioxide storage. In the
benchmark, for which we provide data and code, carbon dioxide is injected into
a saline aquifer modeled by the nonlinear capillarity-free fractional flow
formulation for two incompressible fluid phases, namely carbon dioxide and
brine. To cover different aspects of uncertainty quantification, we incorporate
various sources of uncertainty such as uncertainty of boundary conditions, of
conceptual model definitions and of material properties. We consider recent
versions of the following non-intrusive and intrusive uncertainty
quantification methods: arbitary polynomial chaos, spatially adaptive sparse
grids, kernel-based greedy interpolation and hybrid stochastic Galerkin. The
performance of each approach is demonstrated assessing expectation value and
standard deviation of the carbon dioxide saturation against a reference
statistic based on Monte Carlo sampling. We compare the convergence of all
methods reporting on accuracy with respect to the number of model runs and
resolution. Finally we offer suggestions about the methods' advantages and
disadvantages that can guide the modeler for uncertainty quantification in
carbon dioxide storage and beyond
- …