Robust numerical schemes for anisotropic diffusion problems, a first step for turbulence modeling in Lagrangian hydrodynamics

Abstract

Numerous systems of conservation laws are discretized on Lagrangian meshes where cells nodes move with matter. For complex applications, cells shape or aspect ratio often do not insure sufficient accuracy to provide an acceptable numerical solution and use of ALE technics is necessary. Here we are interested with conduction phenomena depending on velocity derivatives coming from the resolution of gas dynamics equations. For that, we propose the study of a mock of second order turbulent mixing model combining an elliptical part and an hyperbolic kernel. The hyperbolic part is approximated by finite-volume centered scheme completed by a remapping step see [7]. A major part of this paper is the discretization of the anisotropic parabolic equation on polygonal distorted mesh. It is based on the scheme described in [9] ensuring the positivity of the numerical solution. We propose an alternative based on the partitioning of polygons in triangles. We show some preliminary results on a weak coupling of hydrodynamics and parabolic equation whose tensor diffusion coefficient depends on Reynolds stresses.