Issue ESAIM: Proc. Volume 28, 2009 CEMRACS 2008 - Modelling and Numerical Simulation of Complex Fluids Page(s) 80 - 99 DOI http://dx.doi.org/10.1051/proc/2009040 Published online 23 November 2009

ESAIM: Proc., 2009, Vol. 28, pp. 80-99
DOI: 10.1051/proc/2009040

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

Julien Dambrine¹, Philippe Hoch², Raphaël Kuate³, Jérôme Lohéac⁴, Jérôme Métral², Bernard Rebourcet² and Lisl Weynans<sup>1

1</sup>  Université Bordeaux 1 351, cours de la Libération F-33405 Talence cedex, France
²  CEA-DAM, DIF, F 91297, Arpajon, France
³  Université de Valenciennes, LAMAV, LME. Mont Houy 59313 Valenciennes cedex 9, France
⁴  LJK (Laboratoire Jean Kuntzmann) et LSP (Laboratoire de Spectrométrie Physique) Université Joseph Fourier - BP 53 - 38041 Grenoble Cedex 9, France

Published online: 23 November 2009

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.

Résumé
De nombreux systèmes de lois de conservation sont intégrés à l'aide du formalisme lagrangien où les sommets des mailles voient leurs position varier au cours du temps. La forme des mailles ne permet pas toujours d'assurer une bonne précision du calcul et les techniques ALE sont nécessaires. Nous nous intéressons ici à des phénomènes de conduction dépendant du gradient de vitesse couplés à la dynamique des gaz. Pour cela, nous proposons l'étude d'un simulacre de modèle de mélange turbulent du second ordre construit pour combiner terme elliptique et noyau hyperbolique. La partie hyperbolique est résolue par des schéma centrés volumes-finis et le remaillage de la phase ALE est celui décrit dans [7]. Une partie importante de ce papier est la discrétisation de l'équation parabolique anisotrope sur maillage polygonal non régulier. Celle-ci s'inspire du schéma dans [9] qui assure la positivité de la solution numérique. Nous proposons une variante basée sur le découpage en triangles des mailles polygonales. Nous montrons quelques résultats préliminaires sur un couplage faible hydrodynamique et équation parabolique dont le tenseur de diffusion dépend des tensions de Reynolds.

© EDP Sciences, ESAIM 2009

What is OpenURL?