11,404 research outputs found
Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes
We present a new family of high order accurate fully discrete one-step
Discontinuous Galerkin (DG) finite element schemes on moving unstructured
meshes for the solution of nonlinear hyperbolic PDE in multiple space
dimensions, which may also include parabolic terms in order to model
dissipative transport processes. High order piecewise polynomials are adopted
to represent the discrete solution at each time level and within each spatial
control volume of the computational grid, while high order of accuracy in time
is achieved by the ADER approach. In our algorithm the spatial mesh
configuration can be defined in two different ways: either by an isoparametric
approach that generates curved control volumes, or by a piecewise linear
decomposition of each spatial control volume into simplex sub-elements. Our
numerical method belongs to the category of direct
Arbitrary-Lagrangian-Eulerian (ALE) schemes, where a space-time conservation
formulation of the governing PDE system is considered and which already takes
into account the new grid geometry directly during the computation of the
numerical fluxes. Our new Lagrangian-type DG scheme adopts the novel a
posteriori sub-cell finite volume limiter method, in which the validity of the
candidate solution produced in each cell by an unlimited ADER-DG scheme is
verified against a set of physical and numerical detection criteria. Those
cells which do not satisfy all of the above criteria are flagged as troubled
cells and are recomputed with a second order TVD finite volume scheme. The
numerical convergence rates of the new ALE ADER-DG schemes are studied up to
fourth order in space and time and several test problems are simulated.
Finally, an application inspired by Inertial Confinement Fusion (ICF) type
flows is considered by solving the Euler equations and the PDE of viscous and
resistive magnetohydrodynamics (VRMHD).Comment: 39 pages, 21 figure
Nonlinear Methods for Model Reduction
The usual approach to model reduction for parametric partial differential
equations (PDEs) is to construct a linear space which approximates well
the solution manifold consisting of all solutions with
the vector of parameters. This linear reduced model is then used for
various tasks such as building an online forward solver for the PDE or
estimating parameters from data observations. It is well understood in other
problems of numerical computation that nonlinear methods such as adaptive
approximation, -term approximation, and certain tree-based methods may
provide improved numerical efficiency. For model reduction, a nonlinear method
would replace the linear space by a nonlinear space . This idea
has already been suggested in recent papers on model reduction where the
parameter domain is decomposed into a finite number of cells and a linear space
of low dimension is assigned to each cell.
Up to this point, little is known in terms of performance guarantees for such
a nonlinear strategy. Moreover, most numerical experiments for nonlinear model
reduction use a parameter dimension of only one or two. In this work, a step is
made towards a more cohesive theory for nonlinear model reduction. Framing
these methods in the general setting of library approximation allows us to give
a first comparison of their performance with those of standard linear
approximation for any general compact set. We then turn to the study these
methods for solution manifolds of parametrized elliptic PDEs. We study a very
specific example of library approximation where the parameter domain is split
into a finite number of rectangular cells and where different reduced
affine spaces of dimension are assigned to each cell. The performance of
this nonlinear procedure is analyzed from the viewpoint of accuracy of
approximation versus and
Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion
We present a class of spline finite element methods for time-domain wave
propagation which are particularly amenable to explicit time-stepping. The
proposed methods utilize a discontinuous Galerkin discretization to enforce
continuity of the solution field across geometric patches in a multi-patch
setting, which yields a mass matrix with convenient block diagonal structure.
Over each patch, we show how to accurately and efficiently invert mass matrices
in the presence of curved geometries by using a weight-adjusted approximation
of the mass matrix inverse. This approximation restores a tensor product
structure while retaining provable high order accuracy and semi-discrete energy
stability. We also estimate the maximum stable timestep for spline-based finite
elements and show that the use of spline spaces result in less stringent CFL
restrictions than equivalent piecewise continuous or discontinuous finite
element spaces. Finally, we explore the use of optimal knot vectors based on L2
n-widths. We show how the use of optimal knot vectors can improve both
approximation properties and the maximum stable timestep, and present a simple
heuristic method for approximating optimal knot positions. Numerical
experiments confirm the accuracy and stability of the proposed methods
- …