Volume 24, 2008CEMRACS 2007
|Page(s)||111 - 129|
|Published online||10 October 2008|
Non conformal adaptation and mesh smoothing for compressible Lagrangian fluid dynamics
CEA-DIF, Bruyères-Le-Châtel 91297, Arpajon Cedex, France,
2 Université Henri Poincaré Nancy 1., Institut Elie Cartan Nancy (Mathématiques), B.P. 239, F-54506 Vandoeuvre-lès-Nancy Cedex, France,
3 University of West Bohemia, Faculty of Applied Sciences, Department of Mathematics, Univerzitni 22, Plzen 306 14, Czech Republic,
4 Université Paris-Est Marne-la-Vallée, LETEM, Champs sur Marne, 77454 Marne-la-Vallée, Cedex 2, France,
The context of the present work is the numerical approximation of two-dimensional compressible Lagrangian fluid dynamics. The equation on ρ , U and E respectively the density, the velocity and the total energy describes the conservation of mass, momentum, and total energy of the fluid. For a numerical approximation, the starting point is a first order non direct ALE (Arbitrary Lagrangian Eulerian method) that consists of two stages. The first step is a pure Lagrangian step with a centered scheme  or  (for staggered schemes see , , , but we did not use them in this work). The second step is the union of a rezoning step (mesh is smoothed to cure pathological cells and/or adapted to concentrate nodes near singularities, but the connectivity is fixed here) , , , , and a remapping step, defining the old values on this new mesh in a conservative way (see , , ). In some cases, we need to refine/derefine locally the moving mesh, so that this adaptation step takes place between existing Lagrangian/ALE formulations, making the entire process a three steps scheme. We propose an adaptation strategy based on a variant of local non conformal cells. More precisely, as we will see the notion of “non conformal” is only partial, meaning that the new nodes created by adaptation can be a degree of freedom exactly as the master nodes, we call them semi-slaves. With a different approach for local non conformal quadrangular based refinement (see ), in our case unlike their Lagrangian step, we consider all nodes as degree of freedom meaning that we do not consider hanging nodes. The refinement step is generic with respect to the polygonal mesh, the number of nodes in a cell is arbitrary and these are not restricted to be a constant (only triangles and/or quadrilaterals), the conservation is straightforward. The derefinement step is done by local non conformity and the conservation can be obtain either by self-intersection volumes fluxing (line clipping) or by polygonal clipping (more precise than coarse volume averaging). Our approach naturally leads to define this notion of non conformity as a special object inside the description of a generic unstructured polygonal mesh. The development and the simulations have been done with the GO++ package. keywords: Polygonal Unstructured ALE, Hybrid Mesh Representation, Local non Conforming Cell, Dynamic Adaptation, Centered Lagrangian Hydrodynamic, C++.
© EDP Sciences, ESAIM, 2008
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.