Positivity-preserving cell-centered Lagrangian schemes for multi-material compressible flows: From first-order to high-orders

One of the main issues in the field of numerical schemes is to ally robustness with accuracy. Considering gas dynamics, numerical approximations may generate negative density or pressure, which may lead to nonlinear instability and crash of the code. This phenomenon is even more critical using a Lagrangian formalism, the grid moving and being deformed during the calculation. Furthermore, most of the problems studied in this framework contain very intense rarefaction and shock waves. In this paper, the admissibility of numerical solutions obtained by high-order finite-volume-scheme-based methods, such as the discontinuous Galerkin (DG) method, the essentially non-oscillatory (ENO) and the weighted ENO (WENO) finite volume schemes, is addressed in this Lagrangian gas dynamics framework. To this end, we first focus on the one-dimensional case. After briefly recalling how to derive Lagrangian forms of the gas dynamics system of equations, a discussion on positivity-preserving approximate Riemann solvers, ensuring first-order finite volume schemes to be positive, is then given. This study is conducted for both ideal gas and non ideal gas equations of state (EOS), such as the Jones-Wilkins-Lee (JWL) EOS or the Mie-Grüneisen (MG) EOS. It enables us to derive time step conditions ensuring the desired positivity property, as well as L 1 stability of the specific volume and total energy over the domain. Then, making use of the work presented in [74, 75, 15], this positivity study is extended to high-orders of accuracy, where new time step constraints are obtained, and proper limitation is required. This whole analysis is finally applied to the two-dimensional case, and shown to fit a wide range of numerical schemes in the literature, such as the GLACE scheme [12], the EUCCLHYD scheme [55], the GLACE scheme on conical meshes [8], and the LCCDG method [72]. Through this new procedure, scheme robustness is highly improved and hence new problems can be tackled. Numerical results are provided to demonstrate the effectiveness of these methods. Finally, let us emphasize that even if this paper is concerned with purely Lagrangian schemes, the theory developed is of fundamental importance for any methods relying on a purely Lagrangian step, as ALE methods or non-direct Euler schemes

Numerical methods for all-speed flows in fluid-dynamics and non-linear elasticity

In this thesis we are concerned with the numerical simulation of compressible materials flows, including gases, liquids and elastic solids. These materials are described by a monolithic Eulerian model of conservation laws, closed by an hyperelastic state law that includes the different behaviours of the considered materials. A novel implicit relaxation scheme to solve compressible flows at all speeds is proposed, with Mach numbers ranging from very small to the order of unity. The scheme is general and has the same formulation for all the considered materials, since a direct dependence on the state law is avoided via the relaxation. It is based on a fully implicit time discretization, easily implemented thanks to the linearity of the transport operator in the relaxation system. The spatial discretization is obtained by a combination of upwind and centered schemes in order to recover the correct numerical viscosity in different Mach regimes. The scheme is validated with one and two dimensional simulations of fluid flows and of deformations of compressible solids. We exploit the domain discretization through Cartesian grids, allowing for massively parallel computations (HPC) that drastically reduce the computational times on 2D test cases. Moreover, the scheme is adapted to the resolution on adaptive grids based on quadtrees, implementing adaptive mesh refinement techinques. The last part of the thesis is devoted to the numerical simulation of heterogeneous multi-material flows. A novel sharp interface method is proposed, with the derivation of implicit equilibrium conditions. The aim of the implicit framework is the solution of weakly compressible and low Mach flows, thus the proposed multi-material conditions are coupled with the implicit relaxation scheme that is solved in the bulk of the flow. Dans cette thèse on s’intéresse à la simulation numérique d’écoulements des matériaux compressibles, voir fluides et solides élastiques. Les matériaux considérés sont décrits avec un modèle monolithique eulérian, fermé avec une loi d’état hyperélastique qui considère les différents comportéments des matériaux. On propose un nouveau schéma de relaxation qui résout les écoulements compressibles dans des différents régimes, avec des nombres de Mach très petits jusqu’à l’ordre 1. Le schéma a une formulation générale qui est la même pour tous le matériaux considérés, parce que il ne dépend pas directement de la loi d’état. Il se base sur une discrétization complétement implicite, facile à implémenter grâce à la linearité de l’opérateur de transport du système de relaxation. La discrétization en éspace est donnée par la combinaison de flux upwind et centrés, pour retrouver la correcte viscosité numérique dans les différents régimes. L’utilisation de mailles cartésiennes pour les cas 2D s’adapte bien à une parallélisation massive, qui permet de réduire drastiquement le temps de calcul. De plus, le schéma a été adapté pour la résolution sur des mailles quadtree, pour implémenter l’adaptivité de la maille avec des critères entropiques. La dernière partie de la thèse concerne la simulation numérique d’écoulements multi-matériaux. On a proposé une nouvelle méthode d’interface “sharp”, en dérivant les conditions d’équilibre en implicite. L’objectif est la résolution d’interfaces physiques dans des régimes faiblement compressibles et avec un nombre de Mach faible, donc les conditions multi-matériaux sont couplées au schéma implicite de relaxation

A multimaterial Eulerian approach for fluid-solid interaction

This thesis is devoted to understanding and modeling multimaterial interactions, and to develop accordingly a robust scheme taking into account the largest variety of those, with a particular interest in resolving solid/fluid configurations. This very general frame of studies can be tackled with numerous different approaches as several issues arise and need to be addressed before attempting any modelisation of these problems. A first questioning should be the frame of reference to be used for the materials considered. Eulerian shock-capturing schemes have advantages for modeling problems involving complex non-linear wave structures and large deformations. If originally reserved mostly to fluids components, recent work has focused on extending Eulerian schemes to other media such as solid dynamics, as long as the set of equations employed is written under a hyperbolic system of conservation laws. Another matter of interest when dealing with multiple immiscible materials it the necessity to include some means of tracking material boundaries within a numerical scheme. Interface tracking methods based on the use of level set functions are an attractive alternative for problems with sliding interfaces since it allows discontinuous velocity profiles at the material boundaries whilst employing fixed grids. However, its intrinsic lack of variables conservation needs to be circumvented by applying an appropriate fix near the interface, where cells might comprise multiple components. Another requirement is the ability to correctly predict the physical interaction at the interface between the materials. For that purpose, the Riemann problem corresponding to the interfacial conditions needs to be formulated and solved. This implies in turn the need of appropriate Riemann solvers; if they are largely available when the materials are identical (i.e. governed by the same set of equations), a specific Riemann solver will be developed to account for fluid/solid interaction. Eventually, these newly developed methods will be tested on a wide range of different multimaterial problems, involving several materials undergoing large deformations. The materials used, whether modelling fluid/fluid or solid/fluid interactions, will be tested using various initial conditions from both sides of the interface, to demonstrate the robustness of the solver and its flexibility. These testcases will be carried out in 1D, 2D and 3D frames, and compared to exact solutions or other numerical experiments conducted in previous studies.EThOS - Electronic Theses Online ServiceGBUnited Kingdo