Search

Summary: The purpose of this work is to build a general framework to reconstruct the underlying fields within a finite volume (FV) scheme solving a hyperbolic system of PDEs (Partial Differential Equations). In an FV context, the data are piecewise constants per computational cell and the physical fields are reconstructed taking into account neighbor cell values. These reconstructions are further used to evaluate the physical states usually used to feed a Riemann solver which computes the associated numerical fluxes. The physical field reconstructions must obey some properties linked to the system of PDEs such as the positivity, but also some numerically based ones, like an essentially nonoscillatory behavior. Moreover, the reconstructions should be highly accurate for smooth flows and robust/stable for discontinuous solutions. To ensure a solution property preserving reconstruction, we introduce a methodology to blend high‐/low‐order polynomials and hyperbolic tangent reconstructions. A boundary variation diminishing algorithm is employed to select the best reconstruction in each cell. An a posteriori MOOD detection procedure is employed to ensure the positivity by recomputing the rare problematic cells using a robust first‐order FV scheme. We illustrate the performance of the proposed scheme via the numerical simulations for some benchmark tests which involve vacuum or near vacuum states, strong discontinuities, and also smooth flows. The proposed scheme maintains high accuracy on smooth profile, preserves the positivity and can eliminate the oscillations in the vicinity of discontinuities while maintaining sharper discontinuity with superior solution quality compared to classical highly accurate FV schemes. [ABSTRACT FROM AUTHOR]

SUMMARY The aim of the present work is the 3D extension of a general formalism to derive a staggered discretization for Lagrangian hydrodynamics on unstructured grids. The classical compatible discretization is used; namely, momentum equation is discretized using the fundamental concept of subcell forces. Specific internal energy equation is obtained using total energy conservation. The subcell force is derived by invoking the Galilean invariance and thermodynamic consistency. A general form of the subcell force is provided so that a cell entropy inequality is satisfied. The subcell force consists of a classical pressure term plus a tensorial viscous contribution proportional to the difference between the node velocity and the cell-centered velocity. This cell-centered velocity is an extra degree of freedom solved with a cell-centered approximate Riemann solver. The second law of thermodynamics is satisfied by construction of the local positive definite subcell tensor involved in the viscous term. A particular expression of this tensor is proposed. A more accurate extension of this discretization both in time and space is also provided using a piecewise linear reconstruction of the velocity field and a predictor-corrector time discretization. Numerical tests are presented in order to assess the efficiency of this approach in 3D. Sanity checks show that the 3D extension of the 2D approach reproduces 1D and 2D results. Finally, 3D problems such as Sedov, Noh, and Saltzman are simulated. Copyright © 2012 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]