6,203 research outputs found
Numerical simulation of conservation laws with moving grid nodes: Application to tsunami wave modelling
In the present article we describe a few simple and efficient finite volume
type schemes on moving grids in one spatial dimension combined with appropriate
predictor-corrector method to achieve higher resolution. The underlying finite
volume scheme is conservative and it is accurate up to the second order in
space. The main novelty consists in the motion of the grid. This new dynamic
aspect can be used to resolve better the areas with large solution gradients or
any other special features. No interpolation procedure is employed, thus
unnecessary solution smearing is avoided, and therefore, our method enjoys
excellent conservation properties. The resulting grid is completely
redistributed according the choice of the so-called monitor function. Several
more or less universal choices of the monitor function are provided. Finally,
the performance of the proposed algorithm is illustrated on several examples
stemming from the simple linear advection to the simulation of complex shallow
water waves. The exact well-balanced property is proven. We believe that the
techniques described in our paper can be beneficially used to model tsunami
wave propagation and run-up.Comment: 46 pages, 7 figures, 7 tables, 94 references. Accepted to
Geosciences. Other author's papers can be downloaded at
http://www.denys-dutykh.com
Arbitrary Lagrangian-Eulerian form of flowfield dependent variation (ALE-FDV) method for moving boundary problems
Flowfield Dependent Variation (FDV) method is a mixed explicit-implicit numerical
scheme that was originally developed to solve complex flow problems through the use
of so-called implicitness parameters. These parameters determine the implicitness of
FDV method by evaluating local gradients of physical flow parameters, hence vary
across the computational domain. The method has been used successfully in solving
wide range of flow problems. However it has only been applied to problems where the
objects or obstacles are static relative to the flow. Since FDV method has been proved
to be able to solve many complex flow problems, there is a need to extend FDV
method into the application of moving boundary problems where an object
experiences motion and deformation in the flow. With the main objective to develop a
robust numerical scheme that is applicable for wide range of flow problems involving
moving boundaries, in this study, FDV method was combined with a body
interpolation technique called Arbitrary Lagrangian-Eulerian (ALE) method. The
ALE method is a technique that combines Lagrangian and Eulerian descriptions of a
continuum in one numerical scheme, which then enables a computational mesh to
follow the moving structures in an arbitrary movement while the fluid is still seen in a
Eulerian manner. The new scheme, which is named as ALE-FDV method, is
formulated using finite volume method in order to give flexibility in dealing with
complicated geometries and freedom of choice of either structured or unstructured
mesh. The method is found to be conditionally stable because its stability is dependent
on the FDV parameters. The formulation yields a sparse matrix that can be solved by
using any iterative algorithm. Several benchmark stationary and moving body
problems in one, two and three-dimensional inviscid and viscous flows have been
selected to validate the method. Good agreement with available experimental and
numerical results from the published literature has been obtained. This shows that the
ALE-FDV has great potential for solving a wide range of complex flow problems
involving moving bodies
Numerical analysis of conservative unstructured discretisations for low Mach flows
This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving. https://authorservices.wiley.com/author-resources/Journal-Authors/licensing-and-open-access/open-access/self-archiving.htmlUnstructured meshes allow easily representing complex geometries and to refine in regions of interest without adding control volumes in unnecessary regions.
However, numerical schemes used on unstructured grids have to be properly defined in order to minimise numerical errors.
An assessment of a low-Mach algorithm for laminar and turbulent flows on unstructured meshes using collocated and staggered formulations is presented. For staggered formulations using cell centred velocity reconstructions the standard first-order method is shown to be inaccurate in low Mach flows on unstructured grids. A recently proposed least squares procedure for incompressible flows is extended to the low Mach regime and shown to significantly improve the behaviour of the algorithm.
Regarding collocated discretisations, the odd-even pressure decoupling is handled through a kinetic energy conserving flux interpolation scheme. This approach is shown to efficiently handle variable-density flows.
Besides, different face interpolations schemes for unstructured meshes are analysed.
A kinetic energy preserving scheme is applied to the momentum equations, namely the Symmetry-Preserving (SP) scheme. Furthermore, a new approach to define the far-neighbouring nodes of the QUICK scheme is presented and analysed. The method is suitable for both structured and unstructured grids, either uniform or not.
The proposed algorithm and the spatial schemes are assessed against a function reconstruction, a differentially heated cavity and a turbulent self-igniting diffusion flame. It is shown that the proposed algorithm accurately represents unsteady variable-density flows. Furthermore, the QUICK schemes shows close to second order behaviour on unstructured meshes and the SP is reliably used in all computations.Peer ReviewedPostprint (author's final draft
- …