Le code CAVIF

Le code CAVIF

Après avoir fait le choix d’une modélisation de type mélange homogène pour aborder le problème du fractionnement primaire d’un jet de carburant dans un moteur Diesel, il convient à présent d’implanter ce modèle dans le code de calcul le plus adapté à nos besoins. Le code CAVIF (Cavitating Internal Flow), développé à l’IFP par Dumont et al. [56, 54], traite déjà de manière eulérienne l’écoulement cavitant à l’intérieur de l’injecteur. C’est un code diphasique tridimensionnel, dérivé du code KMB [73], qui utilise un modèle de mélange pour représenter une phase notée l (carburant liquide) et une phase notée v (vapeur représentant les poches de cavitation). En outre, ce code permet, par la prise en compte de maillages multi-blocs, de calculer des géométries complexes analogues à la réalité, avec également la possibilité de simuler la levée et descente d’aiguille de l’injecteur grâce à un module de mouvement de maillage. Si l’on considère ce qui a été mis en lumière dans la partie 1.1, il est primordial de prendre en compte les phénomènes intervenant à l’intérieur de l’injecteur pour pouvoir représenter de manière satisfaisante la formation du spray. Une approche eulérienne utilisant un modèle de mélange peut être envisagée simplement dans le code CAVIF. De plus, en étudiant l’atomisation des jets directement avec CAVIF, le couplage avec l’écoulement sortant de l’injecteur sera effectué implicitement. La première partie de ce chapitre est consacrée à la description du code CAVIF et de ses schémas numériques. Ensuite, quelques tests numériques, visant à évaluer la diffusion numérique des méthodes utilisées, sont présentés, avant de détailler les modifications apportées au code afin de permettre la prise en compte de la phase gazeuse g. 

Description générale du code

Une présentation succincte du code CAVIF est préalablement nécessaire. Dans la version originale du code, développé par Dumont [54], les équations de Navier-Stokes tridimensionnelles sont résolues avec les termes de compressibilité et de viscosité (équations de conservation de la masse et de la quantité de mouvement), selon une architecture numérique proche de celle du code KMB [73]. J.B. Moreau, Thèse de Doctorat, INP Toulouse, 2005 Chapitre 3 Le code CAVIF 88 Chapitre 3 : Le code CAVIF 3.1.1 Système d’équations résolu On rappelle le système d’équations à résoudre, détaillé à la partie 2.4, pour un mélange de carburant f : ⎧ ⎪⎪⎪⎪⎨ ⎪⎪⎪⎪⎩ ∂ρf ∂t + ∂ρfuj ∂xj = 0 ∂ρfui ∂t + ∂ρfuiuj ∂xj = −∂pf ∂xi + ∂τf,ij ∂xj pf = f(ρf ) (3.1) Le système 3.1 est résolu en exprimant, pour tout point du maillage, la variation temporelle de la masse volumique du mélange f ∂ρf ∂t  , celle de la quantité de mouvement dans les trois directions de l’espace ∂ρfu ∂t  , ∂ρf v ∂t  et ∂ρfw ∂t  . On en déduit alors aisément les grandeurs physiques ρf , u, v, w et ensuite la pression pf grâce à l’équation d’état pf = f(ρf ). 3.1.2 Structure numérique Comme dans le code KMB [73], CAVIF utilise une formulation ALE1 [80, 134] pour laquelle les éléments fluides se déplacent avec une vitesse arbitraire, dite « vitesse de grille ». Outre sa capacité à traiter des domaines déformables, cette méthode présente l’avantage de décomposer l’opérateur différentiel temporel en trois étapes : – Étape A : contribution des termes sources (chimie, évaporation du spray). – Étape B : résolution du gradient de pression et de la diffusion visqueuse. – Étape C : prise en compte des termes convectifs. Les étapes A et B sont purement lagrangiennes, puisque le maillage suit le fluide avec sa vitesse (les termes convectifs n’interviennent pas). En revanche, l’étape C est eulérienne car le maillage se déplace avec sa vitesse propre et les termes convectifs, résultant de la différence de vitesse entre le fluide et le maillage, sont pris en compte dans le calcul de la solution. Le modèle diphasique, présenté au chapitre précédent, traite un seul « pseudofluide » dont les caractéristiques physiques évoluent. Ainsi, quand bien même l’écoulement présenterait des changements de phase physiques (« vaporisation » du carburant liquide, par cavitation, ou « recondensation » de la vapeur de carburant, par implosion), le code ne prendrait en compte qu’un seul et même fluide, du fait que les changements de phase sont simplement déduits des variations de la masse volumique du fluide f. L’étape A se limite donc dans CAVIF au calcul de pf = f(ρf ).

Discrétisation spatiale

La méthode de discrétisation spatiale utilisée est de type volumes finis. On suppose alors que toute grandeur a une valeur uniforme à l’intérieur d’une maille (ou cellule élémentaire de calcul). On résout donc les équations de grandeurs moyennées 1ALE : Arbitrary Lagrangian-Eulerian Modélisation de l’écoulement polyphasique à l’intérieur et en sortie des injecteurs Diesel 3.1 Description générale du code 89 sur un volume de contrôle donné, noté Ω, comme on a pu le voir dans la partie 2.2. Ainsi, pour calculer la divergence d’une grandeur Φ en un point, on utilise le théorème de Gauss-Ostrogradski :  Ω ∂Φ ∂xi dω =  ∂Ω Φ.n dS (3.2) avec n la normale unitaire à la surface délimitant le volume Ω et dirigée vers l’extérieur. CAVIF utilise des maillages structurés, constitués d’hexaèdres (mailles à six faces), de sorte que l’on peut décomposer cette relation comme la somme des flux traversant chaque face, comme on le verra dans la partie 

Repérage local des nœuds et faces

Afin de repérer les huit nœuds (ou sommets) de chaque maille, on associe à chacun un numéro, conformément à la figure 3.1a. Pour les calculs de flux, il est également utile de repérer chacune des six faces de la maille : la face l (left) est la face délimitée par les nœuds 4378, la face r (right) par les nœuds 1265, la face b (bottom) par les nœuds 4123, la face t (top) par les nœuds 8567, la face f (front) par les nœuds 4158 et enfin la face d (derrière) par les nœuds 3267. Étant donné que la face r, par exemple, est aussi la face l de la maille voisine de droite, et que le même raisonnement peut être appliqué aux faces t et d, on n’utilise finalement que les trois faces caractéristiques l, b et f pour calculer les flux. (a) Numérotation locale des nœuds d’une maille et repérage des faces utilisées pour le calcul des flux (b) Volume de contrôle de quantité de mouvement pour le nœud i4, dans une configuration 2D Fig. 3.1 – Discrétisation spatiale dans CAVIF Repérage des mailles et des nœuds Les maillages utilisés étant structurés, on repère une maille par l’indice correspondant à son nœud local n˚4, noté i4. Les indices de maille sont incrémentés d’abord selon l’axe I, puis selon J et enfin K. L’indice i4 d’une maille quelconque est donc déterminé de la manière suivante : i4 = (k-1)*nyp + (j-1)*nxp + i  : Le code CAVIF où les indices i, j et k sont les indices du nœud i4 selon les axes I, J et K, nxp et nyp sont les nombres de nœuds du maillage dans les directions I et J. Flagage des nœuds et mailles A chaque nœud du maillage correspond un certain nombre de flags : il s’agit d’indices prenant des valeurs qui renseignent sur la nature de ce nœud. Ainsi, le flag de maille, noté f, prend une valeur différente selon que le nœud considéré désigne une maille réelle ou fictive. Le flag f renseigne également sur les conditions limites appliquées. Le flag de nœud, noté ifv, spécifie la nature du nœud (fluide, paroi, entrée, sortie). Enfin, le flag de repère, noté ifrez, sert essentiellement à gérer les interfaces entre maillages, en multi-blocs [73]. 3.1.4 Les volumes de contrôle Les vitesses sont localisées aux nœuds, alors que les quantités scalaires (ρf , pf , μf ) sont définies au centre des mailles. On utilise donc deux volumes de contrôle différents : pour les scalaires il s’agit de la maille elle-même (hexaèdre à 6 faces), alors que le volume de contrôle de la quantité de mouvement est en fait constitué des huit parties d’hexaèdres entourant le nœud concerné (cf. figure 3.1b). Ce volume de contrôle possède donc 24 faces en 3D.

Les conditions aux limites

L’étude d’écoulements fortement instationnaires et le modèle implanté dans le code CAVIF impliquent, lors des simulations, la propagation de fortes ondes de pression (l’équation d’état considère une phase liquide compressible). Dans le cas du calcul de l’écoulement à l’intérieur de l’injecteur seulement [54], il était alors indispensable de choisir des conditions aux limites capables de prendre en considération les différentes ondes traversant les frontières du domaine de calcul.