Présentation 4

Bitte decken Sie die schraffierte Fläche mit einem Bild ab.
Please cover the shaded area with a picture.
(24,4 x 11,0 cm)
ETUDE D’UNE ÉTINCELLE DANS UNE CONFIGURATION
SIMPLIFIÉE: SIMULATION NUMÉRIQUE
Par:
Malyk Benmouffok, Amath Lo, Pierre Freton
Philippe Teulet et Jean-Jacques Gonzalez
www.continental-corporation.com
CAE XII, Clermont-Ferrand
16 et 17 mars 2015
ETUDE D’UNE ÉTINCELLE DANS UNE CONFIGURATION
SIMPLIFIÉE: SIMULATION NUMÉRIQUE
Contexte de l’étude:
 Diminution des ressources fossiles
 Augmentation de la pollution de l’air
Conséquences:
 Réglementation (normes Euro5, Euro6…)
 Obligation d’amélioration des moteurs pour les constructeurs
 Moteur à allumage commandé.
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
2
Principe de fonctionnement d’un moteur à allumage
commandé
Système d’allumage:
Tension alimentation,
bobine et câbles
Phase d’action de la
bougie d’allumage,
déclenchement de
l’étincelle
A. LALLEMAND, Convertisseurs Thermomecaniques. Ed. Techniques Ingénieur
Injection indirecte*
Injection directe*
*C. Ternel,
Contribution au développement de l’allumage par laser pour les moteurs à combustion interne
Cycle moteur 4 temps
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
3
Motivations
* Véhicule léger et petit utilitaire (-6T), estimation mondiale
› La part de marché des
véhicules essence
(Gazoline) est
considérable.
Fort potentiel d’amélioration avec notamment :
› La diminution de la taille des moteurs
(downsizing)
Performances équivalentes
› La recirculation des gaz d’échappement (EGR)
+ EGR
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
4
FAMAC
Contraintes des moteurs essence « downsizés »:
› Cliquetis
› Difficultés d’initier la combustion dans des mélanges
fortement dilués (air+carburant+EGR)
Projet ANR : Fondamentaux d’Allumage pour les Moteurs à
Allumage Commandé (FAMAC)
› Problématique : Comprendre la physique fondamentale de
l’étincelle dans le but d’optimiser le transfert d’énergie du
plasma généré vers le mélange réactif du milieu
environnant
Simulations
Expérimentations
Configuration simplifiée et configuration réelle
dans l’air et dans un mélange air/propane
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
5
 Données d’entrée nécessaires pour la simulation
 Fonctionnement du modèle
 Données utilisées
 Géométrie
 Conditions initiales
 Résultats et comparaison avec l’expérience
 Conclusions
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
6
Données d’entrée nécessaires pour la simulation
R. Maly (1984)*
L’étincelle est définie par la succession de 4 phases:

Pre-breakdown

Breakdown

Arc

Glow
Conditions nécessaires pour la simulation:
-ETL
-Connaitre le plus tôt possible après le claquage:
 le rayon
 La température
du canal conducteur
 La pression
 Le courant
* Dans l’ouvrage Fuel Economy: In Road Vehicles Powered by Spark Ignition Engines. Springer, chap3, 1984.
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
7
Fonctionnement du modèle
Simulation:
2D axisymétrique, transitoire avec ANSYS Fluent:
Air
Propriétés thermodynamiques:
(f(T,P); 300-60 000K; 1-400bar)
Hypothèses:
E.T.L.
Laminaire
Parois d’électrodes adiabatiques
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
8
Données utilisées
Géométrie
Géométrie basée sur configuration expérimentale du projet FAMAC(durée de décharge <700ns,
meilleure stabilité de la décharge):
 Distance inter-électrodes : 3mm
 Rayon de courbure des électrodes : 150µm
 Dimensions du domaine de calcul : 10 x10 mm
Nombre de mailles: 204676 cellules
Taille minimale des cellules: 2µm
Pas de temps: 10ns
Durée de simulation: 3.5 jours
Densité de courant
Jz
• Domaine ouvert
• Pext : 1bar
V=0V
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
9
Conditions initiales
Détermination expérimentale
Conditions initiales à t=15ns
Profil de température N 2
Profil de Température électronique
40000
1200
1100
36000
1000
Température (K)
34000
32000
30000
28000
26000
900
800
700
600
24000
500
22000
Arc
Streamer
400
20000
300
18000
0
20
40
60
80
100
120
0
100
200
300
400
Rayon (µm)
Rayon (µm)
 La température initiale dans le canal est déterminée à partir des
mesures expérimentales de température électronique dans la
phase d’arc et des molécules de N2 dans la phase de streamer
 Le canal initial de la simulation présentera donc une partie très
chaude et aura un préchauffage périphérique dû au streamer
CAE XII
Clermont-Ferrand
16 et 17 mars
500
600
700
-
Température électronique (K)
38000
+
3/19/2015
Benmouffok Malyk, © Continental AG
10
Conditions initiales
Détermination expérimentale
Conditions initiales à t=15ns
-
100
Profil de densité électronique
17
-3
Densité électronique (10 cm )
90
+
80
70
Plus de signal au-delà de
140µm.
60
Le rayon choisi est
50
donc de 140µm
40
0
20
40
60
80
100
120
140
160
Rayon (µm)
 Rayon déterminé à l’aide des mesures de profil
de densité électronique
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
11
Conditions initiales
Détermination expérimentale
Conditions initiales à t=15ns
• On relève la position de l’onde
(Schlieren)
• On ajoute le point obtenu à l’aide
des mesure de Ne
• On obtient la vitesse => la pression
• Braginski
• Rankine-Hugoniot
 La vitesse de déplacement de l’onde de pression permet de
remonter à la pression initiale dans les premiers instants de l’arc,
juste après le claquage.
CAE XII
Clermont-Ferrand
16 et 17 mars
P100bar
3/19/2015
Benmouffok Malyk, © Continental AG
12
Conditions initiales
Application au modèle
Conditions initiales à t=15ns
-
-
Température = 36kK->1200K
Préchauffage = 1200K->300K
Pression = 100bar
Rayon = 140µm
Mesures réalisées au
centre :
-Profil de température
-Densité électronique
Reste du domaine:
1bar/300K
Energie déposée en
condition initiale: 15mJ
180
+
+
160
Current measurement
Interpolation
140
Current (A)
120
100
Ces conditions initiales sont propres à
la configuration utilisée
80
60
40
20
0
-20
0
100
200
300
400
500
600
700
800
Time (ns)
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
13
Résultats
Champs de pression et de température
Température
Pression
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
14
Résultats
Comparaison Expérience/modèle
CORIA mesures expérimentales - r=0mm; z=1.5mm
LAPLACE simulation - r=0mm; z=1.5mm
50000
-
Temperature (K)
45000
40000
35000
+
30000
25000
0
50
100
150
200
250
300
350
400
Temps (ns)
Evolution de la température sur l’axe de la décharge, à mi-distance des électrodes.
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
15
Résultats
Comparaison Expérience/modèle
t=30µs
La figure de droite
correspond à un gradient
de densité reconstruit
numériquement
Simulation
Image Schlieren
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
16
Résultats
Comparaison Expérience/modèle
t=30µs
Simulation
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
17
Résultats
Comparaison Expérience/modèle
t=30µs
Simulation
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
18
Résultats
Comparaison Expérience/modèle
t=30µs
Simulation
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
19
Résultats
Comparaison Expérience/modèle
t=30µs
Simulation
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
20
Résultats
Comparaison Expérience/modèle
t=30µs
Simulation
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
21
Résultats
Comparaison Expérience/modèle
t=30µs
Simulation
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
22
Résultats
Comparaison Expérience/modèle
t=30µs
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
23
Résultats
Comparaison Expérience/modèle
3,0
EM2C mesures Schlieren
LAPLACE simulation
Distance radiale (mm)
2,5
2,0
Comparaison de la
position de l’onde de
pression schlieren avec
le modèle.
1,5
1,0
0,5
0,0
0
500
1000
1500
2000
2500
3000
3500
Temps (ns)
CAE XII
Clermont-Ferrand
4000
Donnée Schlieren obtenues par:
Da Xu, Sara LOVASCIO, Christophe Laux
Ecole Centrale Paris
Bâtiment PECLET
Grande Voie des Vignes 92295
Châtenay-Malabry Cedex
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
24
Conclusions
 Comparaison et validation du modèle dans le cas d’un arc
nanoseconde:
•
Très bon accord de l’évolution de la température et de la
position de l’onde.
•
Les volumes des tores de gaz chauds sont également
comparables
Perspectives:
Détermination de nouvelles conditions initiales pour une décharge
d’allumage conventionnelle
Observer le comportement de l’arc dans une géométrie réelle de
bougie, en particulier lorsque ce dernier est soumis à un soufflage
latéral
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
25
Merci
de votre attention!
CAE XII
Clermont-Ferrand
16 et 17 mars
3/19/2015
Benmouffok Malyk, © Continental AG
26