526 research outputs found
Addressing the challenges of implementation of high-order finite volume schemes for atmospheric dynamics of unstructured meshes
The solution of the non-hydrostatic compressible Euler equations using Weighted Essentially Non-Oscillatory (WENO) schemes in two and three-dimensional unstructured meshes, is presented. Their key characteristics are their simplicity; accuracy; robustness; non-oscillatory properties; versatility in handling any type of grid topology; computational and parallel efficiency. Their defining characteristic is a non-linear combination of a series of high-order reconstruction polynomials arising from a series of reconstruction stencils. In the present study an explicit TVD Runge-Kutta 3rd -order method is employed due to its lower computational resources requirement compared to implicit type time advancement methods. The WENO schemes (up to 5th -order) are applied to the two dimensional and three dimensional test cases: a 2D rising
Adaptive mesh refinement techniques for high-order finite-volume WENO schemes
This paper demonstrates the capabilities of Adaptive Mesh Refinement Techniques (AMR) on 2D hybrid unstructured meshes, for high order finite volume WENO methods. The AMR technique developed is a conformal adapting unstructured hybrid quadrilaterals and triangles (quads & tris) technique for resolving sharp flow features in accurate manner for steady-state and time dependent flow problems. In this method, the mesh can be refined or coarsened which depends on an error estimator, making decision at the parent level whilst maintaining a conformal mesh, the unstructured hybrid mesh refinement is done hierarchically.When a numerical method can work on a fixed conformal mesh this can be applied to do dynamic mesh adaptation. Two Refinement strategies have been devised both following a H-P refinement technique, which can be applied for providing better resolution to strong gradient dominated problems. The AMR algorithm has been tested on cylindrical explosion test and forward facing step problems
Implementation of a low-mach number modification for high-order finite-volume schemes for arbitrary hybrid unstructured meshes
An implementation of a novel low-mach number treatment for high-order finite-volume schemes using arbitrary hybrid unstructured meshes is presented in this paper. Low-Mach order modifications for Godunov type finite-volume schemes have been implemented successfully for structured and unstructured meshes, however the methods break down for hybrid mesh topologies containing multiple element types. The modification is applied to the UCNS3D finite-volume framework for compressible flow configurations, which have been shown as very capable of handling any type of grid topology. The numerical methods under consideration are the Monotonic Upstream-Centered Scheme for Conservation Laws (MUSCL) and the Weighted Essentially Non-Oscillatory (WENO) schemes for two-dimensional mixed-element type unstructured meshes. In the present study the HLLC Approximate Riemann Solver is used with an explicit TVD Runge-Kutta 3rd-order method due to its excellent scalability. These schemes (up to 5th-order) are applied to well established two-dimensional and three-dimensional test cases. The challenges that occur when applying these methods to low-mach flow configurations is thoroughly analysed and possible improvements and further test cases are suggested
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 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
High-Order Unstructured Lagrangian One-Step WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems: Applications to Compressible Multi-Phase Flows
In this article we present the first better than second order accurate
unstructured Lagrangian-type one-step WENO finite volume scheme for the
solution of hyperbolic partial differential equations with non-conservative
products. The method achieves high order of accuracy in space together with
essentially non-oscillatory behavior using a nonlinear WENO reconstruction
operator on unstructured triangular meshes. High order accuracy in time is
obtained via a local Lagrangian space-time Galerkin predictor method that
evolves the spatial reconstruction polynomials in time within each element. The
final one-step finite volume scheme is derived by integration over a moving
space-time control volume, where the non-conservative products are treated by a
path-conservative approach that defines the jump terms on the element
boundaries. The entire method is formulated as an Arbitrary-Lagrangian-Eulerian
(ALE) method, where the mesh velocity can be chosen independently of the fluid
velocity.
The new scheme is applied to the full seven-equation Baer-Nunziato model of
compressible multi-phase flows in two space dimensions. The use of a Lagrangian
approach allows an excellent resolution of the solid contact and the resolution
of jumps in the volume fraction. The high order of accuracy of the scheme in
space and time is confirmed via a numerical convergence study. Finally, the
proposed method is also applied to a reduced version of the compressible
Baer-Nunziato model for the simulation of free surface water waves in moving
domains. In particular, the phenomenon of sloshing is studied in a moving water
tank and comparisons with experimental data are provided
High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes
We present a new family of very high order accurate direct
Arbitrary-Lagrangian-Eulerian (ALE) Finite Volume (FV) and Discontinuous
Galerkin (DG) schemes for the solution of nonlinear hyperbolic PDE systems on
moving 2D Voronoi meshes that are regenerated at each time step and which
explicitly allow topology changes in time.
The Voronoi tessellations are obtained from a set of generator points that
move with the local fluid velocity. We employ an AREPO-type approach, which
rapidly rebuilds a new high quality mesh rearranging the element shapes and
neighbors in order to guarantee a robust mesh evolution even for vortex flows
and very long simulation times. The old and new Voronoi elements associated to
the same generator are connected to construct closed space--time control
volumes, whose bottom and top faces may be polygons with a different number of
sides. We also incorporate degenerate space--time sliver elements, needed to
fill the space--time holes that arise because of topology changes. The final
ALE FV-DG scheme is obtained by a redesign of the fully discrete direct ALE
schemes of Boscheri and Dumbser, extended here to moving Voronoi meshes and
space--time sliver elements. Our new numerical scheme is based on the
integration over arbitrary shaped closed space--time control volumes combined
with a fully-discrete space--time conservation formulation of the governing PDE
system. In this way the discrete solution is conservative and satisfies the GCL
by construction.
Numerical convergence studies as well as a large set of benchmarks for
hydrodynamics and magnetohydrodynamics (MHD) demonstrate the accuracy and
robustness of the proposed method. Our numerical results clearly show that the
new combination of very high order schemes with regenerated meshes with
topology changes lead to substantial improvements compared to direct ALE
methods on conforming meshes
- …