21 research outputs found
Arbitrary-Lagrangian-Eulerian One-Step WENO Finite Volume Schemes on Unstructured Triangular Meshes
In this article we present a new class of high order accurate
Arbitrary-Eulerian-Lagrangian (ALE) one-step WENO finite volume schemes for
solving nonlinear hyperbolic systems of conservation laws on moving two
dimensional unstructured triangular meshes. A WENO reconstruction algorithm is
used to achieve high order accuracy in space and a high order one-step time
discretization is achieved by using the local space-time Galerkin predictor.
For that purpose, a new element--local weak formulation of the governing PDE is
adopted on moving space--time elements. The space-time basis and test functions
are obtained considering Lagrange interpolation polynomials passing through a
predefined set of nodes. Moreover, a polynomial mapping defined by the same
local space-time basis functions as the weak solution of the PDE is used to map
the moving physical space-time element onto a space-time reference element. To
maintain algorithmic simplicity, the final ALE one-step finite volume scheme
uses moving triangular meshes with straight edges. This is possible in the ALE
framework, which allows a local mesh velocity that is different from the local
fluid velocity. We present numerical convergence rates for the schemes
presented in this paper up to sixth order of accuracy in space and time and
show some classical numerical test problems for the two-dimensional Euler
equations of compressible gas dynamics.Comment: Accepted by "Communications in Computational Physics
WENO schemes on unstructured meshes using a relaxed a posteriori MOOD limiting approach
In this paper a relaxed formulation of the a posteriori Multi-dimensional Optimal Order Detection (MOOD) limiting approach is introduced for weighted essentially non-oscillatory (WENO) finite volume schemes on unstructured meshes. The main goal is to minimise the computational footprint of the MOOD limiting approach by employing WENO schemes—by virtue of requiring a smaller number of cells to reduce their order of accuracy compared to an unlimited scheme. The key characteristic of the present relaxed MOOD formulation is that the Numerical Admissible Detector (NAD) is not uniquely defined for all orders of spatial accuracy, and it is relaxed when reaching a 2nd-order of accuracy. The augmented numerical schemes are applied to the 2D unsteady Euler equations for a multitude of test problems including the 2D vortex evolution, cylindrical explosion, double-Mach reflection, and an implosion. It is observed that in many events, the implemented MOOD paradigm manages to preserve symmetry of the forming structures in simulations, an improvement comparing to the non-MOOD limited counterparts which cannot be easily obtained due to the multi-dimensional reconstruction nature of the schemes
High Order Cell-Centered Lagrangian-Type Finite Volume Schemes with Time-Accurate Local Time Stepping on Unstructured Triangular Meshes
We present a novel cell-centered direct Arbitrary-Lagrangian-Eulerian (ALE)
finite volume scheme on unstructured triangular meshes that is high order
accurate in space and time and that also allows for time-accurate local time
stepping (LTS). The new scheme uses the following basic ingredients: a high
order WENO reconstruction in space on unstructured meshes, an element-local
high-order accurate space-time Galerkin predictor that performs the time
evolution of the reconstructed polynomials within each element, the computation
of numerical ALE fluxes at the moving element interfaces through approximate
Riemann solvers, and a one-step finite volume scheme for the time update which
is directly based on the integral form of the conservation equations in
space-time. The inclusion of the LTS algorithm requires a number of crucial
extensions, such as a proper scheduling criterion for the time update of each
element and for each node; a virtual projection of the elements contained in
the reconstruction stencils of the element that has to perform the WENO
reconstruction; and the proper computation of the fluxes through the space-time
boundary surfaces that will inevitably contain hanging nodes in time due to the
LTS algorithm. We have validated our new unstructured Lagrangian LTS approach
over a wide sample of test cases solving the Euler equations of compressible
gasdynamics in two space dimensions, including shock tube problems, cylindrical
explosion problems, as well as specific tests typically adopted in Lagrangian
calculations, such as the Kidder and the Saltzman problem. When compared to the
traditional global time stepping (GTS) method, the newly proposed LTS algorithm
allows to reduce the number of element updates in a given simulation by a
factor that may depend on the complexity of the dynamics, but which can be as
large as 4.7.Comment: 31 pages, 13 figure
Lagrangian ADER-WENO Finite Volume Schemes on Unstructured Triangular Meshes Based On Genuinely Multidimensional HLL Riemann Solvers
In this paper we use the genuinely multidimensional HLL Riemann solvers
recently developed by Balsara et al. to construct a new class of
computationally efficient high order Lagrangian ADER-WENO one-step ALE finite
volume schemes on unstructured triangular meshes. A nonlinear WENO
reconstruction operator allows the algorithm to achieve high order of accuracy
in space, while high order of accuracy in time is obtained by the use of an
ADER time-stepping technique based on a local space-time Galerkin predictor.
The multidimensional HLL and HLLC Riemann solvers operate at each vertex of the
grid, considering the entire Voronoi neighborhood of each node and allows for
larger time steps than conventional one-dimensional Riemann solvers. The
results produced by the multidimensional Riemann solver are then used twice in
our one-step ALE algorithm: first, as a node solver that assigns a unique
velocity vector to each vertex, in order to preserve the continuity of the
computational mesh; second, as a building block for genuinely multidimensional
numerical flux evaluation that allows the scheme to run with larger time steps
compared to conventional finite volume schemes that use classical
one-dimensional Riemann solvers in normal direction. A rezoning step may be
necessary in order to overcome element overlapping or crossing-over. We apply
the method presented in this article to two systems of hyperbolic conservation
laws, namely the Euler equations of compressible gas dynamics and the equations
of ideal classical magneto-hydrodynamics (MHD). Convergence studies up to
fourth order of accuracy in space and time have been carried out. Several
numerical test problems have been solved to validate the new approach
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