Search

In this paper we present a conservative cell-centered Lagrangian finite volume scheme for the solution of the hyper-elasticity equations on unstructured multidimensional grids. The starting point of the new method is the Eucclhyd scheme, which is here combined with the a posteriori Multidimensional Optimal Order Detection (MOOD) limiting strategy to ensure robustness and stability at shock waves with piece-wise linear spatial reconstruction. The ADER (Arbitrary high order schemes using DERivatives) approach is adopted to obtain second-order of accuracy in time as well. This method has been tested in an hydrodynamics context and the present work aims at extending it to the case of hyper-elasticity models. Such models are presented in a fully Lagrangian framework and the dedicated Lagrangian numerical scheme is derived in terms of nodal solver, GCL compliance, subcell forces and compatible discretization. The Lagrangian numerical method is implemented in 3D under MPI parallelization framework allowing to handle genuinely large meshes. A relative large set of numerical test cases is presented to assess the ability of the method to achieve effective second order of accuracy on smooth flows, maintaining an essentially non-oscillatory behavior and general robustness across discontinuities and ensuring at least physical admissibility of the solution where appropriate. Pure elastic neo-Hookean and non-linear materials are considered for our benchmark test problems in 2D and 3D. These test cases feature material bending, impact, compression, non-linear deformation and further bouncing/detaching motions.

We propose an adaptive stencil construction for high order accurate finite volume schemes aposteriori stabilized devoted to solve one-dimensional steady-state hyperbolic equations. High-accuracy (up to the sixth-order presently) is achieved thanks to polynomial reconstructions while stability is provided with an aposteriori MOOD method which controls the cell polynomial degree for eliminating non-physical oscillations in the vicinity of discontinuities. We supplemented this scheme with a stencil construction allowing to reduce even further the numerical dissipation. The stencil is shifted away from troubles (shocks, discontinuities, etc.) leading to less oscillating polynomial reconstructions. Experimented on linear, B\"urgers', and Euler equations, we demonstrate that the adaptive stencil technique manages to retrieve smooth solutions with optimal order of accuracy but also irregular ones without spurious oscillations. Moreover we numerically show that the approach allows to reduce the dissipation still maintaining the essentially non-oscillatory behavior.

We present a new high-accurate, stable and low-dissipative Smooth Particle Hydrodynamics (SPH) method based on Riemann solvers. The method derives from the SPH-ALE formulation first proposed by Vila and Ben Moussa. Moving Least Squares approximations are used for the reconstruction of the variables and the computation of Taylor expansions. The stability of the scheme is achieved by the a posteriori Multi-dimensional Optimal Order Detection (MOOD) paradigm. Such a procedure enables to provide genuine gains in accuracy both for one- and two-dimensional problems involving non-smooth flows when compared to classical SPH methods.

We present a new high-accurate, stable and low-dissipative Smooth Particle Hydrodynamics (SPH) method based on Riemann solvers. The method derives from the SPH-ALE formulation first proposed by Vila and Ben Moussa. Moving Least Squares approximations are used for the reconstruction of the variables and the computation of Taylor expansions. The stability of the scheme is achieved by the a posteriori Multi-dimensional Optimal Order Detection (MOOD) paradigm. Such a procedure enables to provide genuine gains in accuracy both for one- and two-dimensional problems involving non-smooth flows when compared to classical SPH methods.

In this paper, we describe a cell-centered Lagrangian scheme devoted to the numerical simulation of solid dynamics on two-dimensional unstructured grids in planar geometry. This numerical method, utilizes the classical elastic-perfectly plastic material model initially proposed by Wilkins [M.L. Wilkins, Calculation of elastic-plastic flow, Meth. Comput. Phys. (1964)]. In this model, the Cauchy stress tensor is decomposed into the sum of its deviatoric part and the thermodynamic pressure which is defined by means of an equation of state. Regarding the deviatoric stress, its time evolution is governed by a classical constitutive law for isotropic material. The plasticity model employs the von Mises yield criterion and is implemented by means of the radial return algorithm. The numerical scheme relies on a finite volume cell-centered method wherein numerical fluxes are expressed in terms of sub-cell force. The generic form of the sub-cell force is obtained by requiring the scheme to satisfy a semi-discrete dissipation inequality. Sub-cell force and nodal velocity to move the grid are computed consistently with cell volume variation by means of a node-centered solver, which results from total energy conservation. The nominally second-order extension is achieved by developing a two-dimensional extension in the Lagrangian framework of the Generalized Riemann Problem methodology, introduced by Ben-Artzi and Falcovitz [M. Ben-Artzi and J. Falcovitz, Generalized Riemann Problems in Computational Fluid Dynamics, Cambridge Monographs on Applied and Computational Mathematics, 2003]. Finally, the robustness and the accuracy of the numerical scheme are assessed through the computation of several test cases.