Nicolas Favrie

Solid-fluid diffuse interface model in cases of extreme deformations

Diffuse interface methods have been recently proposed and successfully used for accurate compressible multi-fluid computations Abgrall [1]; Kapila et al. [20]; Saurel et al. [30]. These methods deal with extended systems of hyperbolic equations involving a non-conservative volume fraction equation and relaxation terms. Following the same theoretical frame, we derive here an Eulerian diffuse interface model for elastic solid-compressible fluid interactions in situations involving extreme deformations. Elastic effects are included following the Eulerian conservative formulation proposed in Godunov [16], Miller and Colella [23], Godunov and Romenskii [17], Plohr and Plohr [27] and Gavrilyuk et al. [14]. We apply first the Hamilton principle of stationary action to derive the conservative part of the model. The relaxation terms are then added which are compatible with the entropy inequality. In the limit of vanishing volume fractions the Euler equations of compressible fluids and a conservative hyperelastic model are recovered. It is solved by a unique hyperbolic solver valid at each mesh point (pure fluid, pure solid and mixture cell). Capabilities of the model and methods are illustrated on various tests of impacts of solids moving in an ambient compressible fluid.

Multi-solid and multi-fluid diffuse interface model

We extend the model of diffuse solid-fluid interfaces developed earlier by authors of this paper to the case of arbitrary number of interacting hyperelastic solids. Plastic transformations of solids are taken into account through a Maxwell type model. The specific energy of each solid is given in separable form: it is the sum of a hydrodynamic part of the energy depending only on the density and the entropy, and an elastic part of the energy which is unaffected by the volume change. It allows us to naturally pass to the fluid description in the limit of vanishing shear modulus. In spite of a large number of governing equations, the model has a quite simple mathematical structure: it is a duplication of a single visco-elastic model. The model is well posed both mathematically and thermodynamically: it is hyperbolic and compatible with the second law of thermodynamics. The resulting model can be applied in the situations involving an arbitrary number of fluids and solids. In particular, we show the ability of the model to describe spallation and penetration phenomena occurring during high velocity impacts.

Modelling dynamic and irreversible powder compaction

A multiphase hyperbolic model for dynamic and irreversible powder compactionis built. Four important points have to be addressed in this case. The first one isrelated to the irreversible character of powder compaction. When a granular media issubjected to a loading–unloading cycle, the final volume is lower than the initial one.To deal with this hysteresis phenomenon, a multiphase model with relaxation is built.During loading, mechanical equilibrium is assumed corresponding to stiff mechanicalrelaxation, while during unloading non-equilibrium mechanical transformation isassumed. Consequently, the sound speed of the limit models are very different duringloading and unloading. These differences in acoustic properties are responsible forirreversibility in the compaction process. The second point is related to dynamiceffects, where pressure and shock waves play an important role. Wave dynamics isguaranteed by the hyperbolic character of the equations. Phase compressibility aswell as configuration energy are taken into account. The third point is related tomulti-dimensional situations that involve material interfaces. Indeed, most processeswith powder compaction entail free surfaces. Consequently, the model should be ableto solve interfaces separating pure fluids and granular mixtures. Finally, the fourthpoint is related to gas permeation that may play an important role in some specificpowder compaction situations. This poses the difficult question of multiple-velocitydescription. These four points are considered in a unique model fitting the frameof multiphase theory of diffuse interfaces (Saurel & Abgrall, J. Comput. Phys.,vol. 150, 1999, p. 425; Kapila et al., Phys. Fluids, vol. 13, 2001, p. 3002; Saurel et al.,J. Comput. Phys., vol. 228, 2009, p. 1678). The ability of the model to deal with thesevarious effects is validated on basic situations, where each phenomenon is consideredseparately. Except for the material EOS (hydrodynamic and granular pressures andenergies), which are determined on the basis of separate experiments found in theliterature, the model is free of adjustable parameter.Key words: granular media, particle/fluid flows, shock waves

Mathematical and numerical model for nonlinear viscoplasticity.

A macroscopic model describing elastic–plastic solids is derived in a special case of the internal specific energy taken in separable form: it is the sum of a hydrodynamic part depending only on the density and entropy, and a shear part depending on other invariants of the Finger tensor. In particular, the relaxation terms are constructed compatible with the von Mises yield criteria. In addition, Maxwell-type material behaviour is shown up: the deviatoric part of the stress tensor decays during plastic deformations. Numerical examples show the ability of this model to deal with real physical phenomena.

A thermodynamically compatible splitting procedure in hyperelasticity

Abstract A material is hyperelastic if the stress tensor is obtained by variation of the stored energy function. The corresponding 3D mathematical model of hyperelasticity written in the Eulerian coordinates represents a system of 14 conservative partial differential equations submitted to stationary differential constraints. A classical approach for numerical solving of such a 3D system is a geometrical splitting: the 3D system is split into three 1D systems along each spatial direction and solved then by using a Godunov type scheme. Each 1D system has 7 independent eigenfields corresponding to contact discontinuity, longitudinal waves and shear waves. The construction of the corresponding Riemann solvers is not an easy task even in the case of isotropic solids. Indeed, for a given specific energy it is extremely difficult, if not impossible, to check its rank-one convexity which is a necessary and sufficient condition for hyperbolicity of the governing equations. In this paper, we consider a particular case where the specific energy is a sum of two terms. The first term is the hydrodynamic energy depending only on the density and the entropy, and the second term is the shear energy which is unaffected by the volume change. In this case a very simple criterion of hyperbolicity can be formulated. We propose then a new splitting procedure which allows us to find a numerical solution of each 1D system by solving successively three 1D sub-systems. Each sub-system is hyperbolic, if the full system is hyperbolic. Moreover, each sub-system has only three waves instead of seven, and the velocities of these waves are given in explicit form. The last property allows us to construct reliable Riemann solvers. Numerical 1D tests confirm the advantage of the new approach. A multi-dimensional extension of the splitting procedure is also proposed.

A model and numerical method for compressible flows with capillary effects

A new model for interface problems with capillary effects in compressible fluids is presented together with a specific numerical method to treat capillary flows and pressure waves propagation. This new multiphase model is in agreement with physical principles of conservation and respects the second law of thermodynamics. A new numerical method is also proposed where the global system of equations is split into several submodels. Each submodel is hyperbolic or weakly hyperbolic and can be solved with an adequate numerical method. This method is tested and validated thanks to comparisons with analytical solutions (Laplace law) and with experimental results on droplet breakup induced by a shock wave.

The piston problem in hyperelasticity with the stored energy in separable form

The piston problem for a hyperelastic hyperbolic conservative model where the stored energy is given in separable form is studied. The eigenfields corresponding to the hyperbolic system are of three types: linearly degenerate fields (corresponding to the contact characteristics), the fields which are genuinely nonlinear in the sense of Lax (corresponding to longitudinal waves), and, finally, nonlinear fields which are not genuinely nonlinear (corresponding to transverse waves). Taking the initial state free of stresses, we presented possible auto-similar solutions to the piston problem. In particular, we have shown that the equations admit transverse shock waves having a remarkable property: the solid density is decreasing through such a shock, it is thus a ‘rarefaction’ shock.

Modeling hyperelasticity in non-equilibrium multiphase flows

The aim of this article is the construction of a multiphase hyperelastic model. The Eulerian formulation of the hyperelasticity represents a system of 14 conservative partial differential equations submitted to stationary differential constraints. This model is constructed with an elegant approach where the specific energy is given in separable form. The system admits 14 eigenvalues with 7 characteristic eigenfields. The associated Riemann problem is not easy to solve because of the presence of 7 waves. The shear waves are very diffusive when dealing with the full system. In this paper, we use a splitting approach to solve the whole system using 3 sub-systems. This method reduces the diffusion of the shear waves while allowing to use a classical approximate Riemann solver. The multiphase model is obtained by adapting the discrete equations method. This approach involves an additional equation governing the evolution of a phase function relative to the presence of a phase in a cell. The system is integrated over a multiphase volume control. Finally, each phase admits its own equations system composed of three sub-systems. One and three dimensional test cases are presented.

Dynamic compaction of granular materials

An Eulerian hyperbolic multiphase flow model for dynamic and irreversible compaction of granular materials is constructed. The reversible model is first constructed on the basis of the classical Hertz theory. The irreversible model is then derived in accordance with the following two basic principles. First, the entropy inequality is satisfied by the model. Second, the corresponding ‘intergranular stress’ coming from elastic energy owing to contact between grains decreases in time (the granular media behave as Maxwell-type materials). The irreversible model admits an equilibrium state corresponding to von Mises-type yield limit. The yield limit depends on the volume fraction of the solid. The sound velocity at the yield surface is smaller than that in the reversible model. The last one is smaller than the sound velocity in the irreversible model. Such an embedded model structure assures a thermodynamically correct formulation of the model of granular materials. The model is validated on quasi-static experiments on loading–unloading cycles. The experimentally observed hysteresis phenomena were numerically confirmed with a good accuracy by the proposed model.

Internal-variable modeling of solids with slow dynamics: Wave propagation and resonance simulations

Rocks and concrete are known to soften under a dynamic loading, i.e. the speed of sound diminishes with forcing amplitudes. To reproduce this behavior, an internal-variable model of continuum is proposed. It is composed of a constitutive law for the stress and an evolution equation for the internal variable. Nonlinear viscoelasticity of Zener type is accounted for by using additional internal variables. The proposed system of partial differential equations is solved numerically using finite-volume methods. The numerical tool is used to reproduce qualitatively Nonlinear Resonance Ultrasound Spectroscopy (NRUS) and Dynamic Acoustoelastic Testing (DAET) experiments. A frequency-domain approach based on finite elements, harmonic balance and numerical continuation is compared to the time-domain method. This approach is promising for upcoming experimental validations with respect to resonance experiments.