A structure-preserving and thermodynamically compatible cell-centered Lagrangian finite volume scheme for continuum mechanics

TL;DR

This paper introduces a novel structure-preserving, thermodynamically compatible cell-centered Lagrangian finite volume scheme for continuum mechanics that exactly maintains key physical laws and properties at the discrete level.

Contribution

It develops a new semi-discrete scheme that preserves energy, entropy, and curl-free conditions simultaneously on moving unstructured meshes, advancing the numerical discretization of the GPR model.

Findings

Achieves exact discrete total energy conservation.

Ensures discrete entropy inequality.

Preserves curl-free properties of fields.

Abstract

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.