A structure-preserving and thermodynamically compatible cell-centered Lagrangian finite volume scheme for continuum mechanics
By: Walter Boscheri , Michael Dumbser , Raphael Loubère and more
Potential Business Impact:
Makes computer models of materials more accurate.
In this work we present a novel structure-preserving scheme for the discretization of the Godunov-Peshkov-Romenski (GPR) model of continuum mechanics written in Lagrangian form. This model admits an extra conservation law for the total energy (first principle of thermodynamics) and satisfies the entropy inequality (second principle of thermodynamics). Furthermore, in the absence of algebraic source terms, the distortion field of the continuum and the specific thermal impulse satisfy a curl-free condition, provided the initial data are curl-free. Last but not least, the determinant of the distortion field is related to the density of the medium, i.e. the system is also endowed with a nonlinear algebraic constraint. The objective of this work is to construct and analyze a new semi-discrete thermodynamically compatible cell-centered Lagrangian finite volume scheme on moving unstructured meshes that satisfies the following structural properties of the governing PDE exactly at the discrete level: i) compatibility with the first law of thermodynamics, i.e. discrete total energy conservation; ii) compatibility with the second law of thermodynamics, i.e. discrete entropy inequality; iii) exact discrete compatibility between the density and the determinant of the distortion field; iv) exact preservation of the curl-free property of the distortion field and of the specific thermal impulse in the absence of algebraic source terms. We show that it is possible to achieve all above properties simultaneously. Unlike in existing schemes, we choose to directly discretize the entropy inequality, hence obtaining total energy conservation as a consequence of an appropriate and thermodynamically compatible discretization of all the other equations.
Similar Papers
High order well-balanced and total-energy-conserving local discontinuous Galerkin methods for compressible self-gravitating Euler equations
Numerical Analysis
Simulates gravity and gas flow accurately.
Variational discretizations of viscous and resistive magnetohydrodynamics using structure-preserving finite elements
Numerical Analysis
Simulates space storms more accurately.
A decoupled and structure-preserving direct discontinuous Galerkin method for the Keller-Segel Model
Numerical Analysis
Models how cells move and group together.