6,936 research outputs found
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
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
Three-dimensional CFD simulations with large displacement of the geometries using a connectivity-change moving mesh approach
This paper deals with three-dimensional (3D) numerical simulations involving 3D moving geometries with large displacements on unstructured meshes. Such simulations are of great value to industry, but remain very time-consuming. A robust moving mesh algorithm coupling an elasticity-like mesh deformation solution and mesh optimizations was proposed in previous works, which removes the need for global remeshing when performing large displacements. The optimizations, and in particular generalized edge/face swapping, preserve the initial quality of the mesh throughout the simulation. We propose to integrate an Arbitrary Lagrangian Eulerian compressible flow solver into this process to demonstrate its capabilities in a full CFD computation context. This solver relies on a local enforcement of the discrete geometric conservation law to preserve the order of accuracy of the time integration. The displacement of the geometries is either imposed, or driven by fluid–structure interaction (FSI). In the latter case, the six degrees of freedom approach for rigid bodies is considered. Finally, several 3D imposed-motion and FSI examples are given to validate the proposed approach, both in academic and industrial configurations
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
An efficient parallel immersed boundary algorithm using a pseudo-compressible fluid solver
We propose an efficient algorithm for the immersed boundary method on
distributed-memory architectures, with the computational complexity of a
completely explicit method and excellent parallel scaling. The algorithm
utilizes the pseudo-compressibility method recently proposed by Guermond and
Minev [Comptes Rendus Mathematique, 348:581-585, 2010] that uses a directional
splitting strategy to discretize the incompressible Navier-Stokes equations,
thereby reducing the linear systems to a series of one-dimensional tridiagonal
systems. We perform numerical simulations of several fluid-structure
interaction problems in two and three dimensions and study the accuracy and
convergence rates of the proposed algorithm. For these problems, we compare the
proposed algorithm against other second-order projection-based fluid solvers.
Lastly, the strong and weak scaling properties of the proposed algorithm are
investigated
- …