Solid-fluid dynamics of yield-stress fluids

On the example of two-phase continua experiencing stress induced solid-fluid phase transitions we explore the use of the Euler structure in the formulation of the governing equations. The Euler structure guarantees that solutions of the time evolution equations possessing it are compatible with mechanics and with thermodynamics. The former compatibility means that the equations are local conservation laws of the Godunov type and the latter compatibility means that the entropy does not decrease during the time evolution. In numerical illustrations, in which the one-dimensional Riemann problem is explored, we require that the Euler structure is also preserved in the discretization.Comment: 51 pages, 7 figure

Coupling techniques for nonlinear hyperbolic equations. III. The well-balanced approximation of thick interfaces

We continue our analysis of the coupling between nonlinear hyperbolic problems across possibly resonant interfaces. In the first two parts of this series, we introduced a new framework for coupling problems which is based on the so-called thin interface model and uses an augmented formulation and an additional unknown for the interface location; this framework has the advantage of avoiding any explicit modeling of the interface structure. In the present paper, we pursue our investigation of the augmented formulation and we introduce a new coupling framework which is now based on the so-called thick interface model. For scalar nonlinear hyperbolic equations in one space variable, we observe that the Cauchy problem is well-posed. Then, our main achievement in the present paper is the design of a new well-balanced finite volume scheme which is adapted to the thick interface model, together with a proof of its convergence toward the unique entropy solution (for a broad class of nonlinear hyperbolic equations). Due to the presence of a possibly resonant interface, the standard technique based on a total variation estimate does not apply, and DiPerna's uniqueness theorem must be used. Following a method proposed by Coquel and LeFloch, our proof relies on discrete entropy inequalities for the coupling problem and an estimate of the discrete entropy dissipation in the proposed scheme.Comment: 21 page

ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement

We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.Comment: With updated bibliography informatio

On selection criteria for problems with moving inhomogeneities

We study mechanical problems with multiple solutions and introduce a thermodynamic framework to formulate two different selection criteria in terms of macroscopic energy productions and fluxes. Studying simple examples for lattice motion we then compare the implications for both resting and moving inhomogeneities.Comment: revised version contains new introduction, numerical simulations of Riemann problems, and a more detailed discussion of the causality principle; 18 pages, several figure