CFD Combustion Models for IC Engines Rolf D. Reitz Engine Research Center University of Wisconsin-Madison ERC Symposium, June 7, 27 http://www.cae.wisc.edu/~reitz
Combustion and Emission Models at the ERC Turbulence vs. chemistry time scales 1. Characteristic Time Combustion Model (KIVA-CTC) Auto-ignition model (SHELL) ignition with 8 reaction steps Combustion model (CTC) turbulence mixing timescales 2. Representative Interactive Flamelet Model (KIVA-RIF) Flamelet equations/detailed chemistry coupled with CFD (Peters et al.) Effect of turbulence is modeled with average scalar dissipation rate 3. Direct Integration of CFD with Chemistry (KIVA-CHEMKIN) Each computational cell treated as well stirred reactor ERC reduced mechanism for n-heptane chemistry timescales 4. KIVA-CHEMKIN-G Model Explicit modeling of turbulent flame propagation combined with detailed chemistry Eddy breakup timescales Chemistry timescales
Auto-ignition and Detailed Chemistry Models Skeletal Chemistry Mechanism for Multi-Surrogate Fuels Gasoline Iso-octane 857 species 3586 reactions 25 species 51 reactions Diesel Bio-diesel N-heptane Methyl butanoate 57 species 258 reactions 264 species 1219 reactions 3 species 65 reactions 41 species 15 reactions 47 41 Multi-fuel Mechanism (ERC-PRF-Bio) Aromatic Alcohol Toluene Ehthanol surrogate fuel Multi-component fuel Comprehensive mechanism Skeletal mechanism Dr. Youngchul Ra Jessica Brakora SAE 24-1-558, SAE 27-1-19, SAE 27-1-165
Validation - Auto-ignition Model n-heptane-iso-octane/air mixtures, stoichiometric, P=4 Bar) Calculation using ERC-PRF mechanism Ignition Delay Time (ms) 1 1 1.1 1% i-octane 9% i-octane / 1% n-heptane 8% i-octane / 2% n-heptane 6% i-octane / 4% n-heptane 1% n-heptane Shock tube measurements (Fieweger et al. 1997) SAE 27-1-165.1.7.8.9 1. 1.1 1.2 1.3 1.4 1.5 1(K)/T
Validation Cat HD engine 6 5 4 3 2 1 6 5 4 3 2 1 rpm = 1737 eq_ratio =.23 c_ratio = 16.1 inj_timing = -113 ca -4-2 2 4 7 rpm CAD ATDC Experiment eq_ratio =.27 Shell-CTC c_ratio = 16.1 CHEMKIN inj_timing = - 274.5 ca -4-2 2 4 CAD ATDC 6 5 4 3 2 1 CO (g/kgf) 18 16 Experiment 9 14 Shell-CTC Gasoline HCCI CHEMKIN 8 12 7 1 AHRR (J/deg) 9 8 7 6 AHRR (J/deg) 5 4 3 2 1 8 6 4 2 3.5 3. 2.5 NOx (g/kgf) 2. 1.5 1..5. Experiment Shell-CTC Chemkin.135 Experiment.2.23.27.29.39 Shell-CTC Equivalence ratio Chemkin.135.2.23.27.29.39 Equivalence ratio SAE 22-1-418, SAE 27-1-19
Validation - GM Light-duty Diesel Engine KIVA-CHEMKIN SOIC = 39, -32 and 23 SOIC=-32 GM 1.9L 4-cylinder Bore Stroke 82 9.4 mm Compression Ratio 16.6: 1 Engine Speed 2 rev/min IMEP 5.5 Bar EGR 67% PCCI SAE 27-1-193
Advances in SI Engine Combustion Modeling G-equation Combustion Model with Detailed Chemical Kinetics SAE 2-1-289 Discrete particle ignition model Burnt Gas G-Equation Flame propagation End Gas ~ G t + ( v f v vertex ~ ) G = u S T Auto-ignition (detailed kinetics) ~ G D T ~ ~ k G SAE 23-1-722 SAE 27-1-165
2.5 2. 1.5 1. Validation Ford SI Engine EXPT SIMU PFI mode -32 O ATDC 2.5 KIVA-CHEMKIN-G 2. 1.5 1. EXPT SIMU PFI mode -44 O ATDC Ford DISI Bore Stroke 89 79.5 mm Compression Ratio 12 : 1 Engine Speed 15 rev/min PFI Spark sweep.5 3.5 25 EXPT Spk = -2 o ATDC. 3. SIMU -1-5 5 12 2.5 2. 1.5 1. Crank Angle ( o ATDC).5 5 4. 35. 3.5 EXPT MAP = 75 kpa -.5 SIMU 3 3. -8-6 -4-2 2 4 6 8 1 2.5 Crank Angle ( o 25 ATDC) 2. 2 1.5 1. 1.5. 5 -.5-8 -6-4 -2 2 4 6 8 1 Crank Angle ( o ATDC) 15 15 1.5 3.5 25 EXPT Spk = -32 o ATDC. 3. SIMU -1-5 5 1 2 2.5 Heat Release Rate (J/Deg Heat Release Pressure Rate (MPa) (J/Deg 2. 1.5 1. Crank Angle ( o ATDC).5 5. 4. 35 -.5 3.5 EXPT MAP = 1 kpa SIMU 3 3. -8-6 -4-2 2 4 6 8 1 Crank Angle ( o ATDC) 25 2.5 2. 2 1.5 1. 1.5. 5 -.5-8 -6-4 -2 2 4 6 8 1 Crank Angle ( o ATDC) 15 1 15 Heat Release Rate (J/De Heat Release Rate (J/Deg DI Spark sweep DI MAP sweep SAE 27-1-165
Validation - Sandia Cummins Diesel High-T, long-id Condition (Premixed Combustion) O 2 = 21% SOI = -5 AHRR (J/deg) 12 1 8 6 4 2 TEST CHEMKIN G_Model Injection Profile -1-15 -1-5 5 1 15 2 25 CAD ATDC High-T, Long-ID Experiment KIVA-CHEMKIN KIVA-CHEMKIN-G 9 8 7 6 5 4 3 2 1 SAE 26-1-55 Green: OH Red: Soot
Validation - Cummins Dual-Fuel Engine Experimental SAE 27-1-171 KIVA-CHEMKIN-G NOx Emissions (grams) Experimental Measurements 18 Case 1 Case 2 Case 3 15 Case 4 Case 5 12 AHRR (J/deg) 9 6 3.3.25.2.15.1.5. -2-2 -15-1 -5 5 1 15 2 25 3 35 4 Measured Predicted CAD ATDC NOx Case1 Case2 Case3 Case4 Case5 C 16 14 12 1 8 6 4 2 AHRR (J/deg) 18 15 12 9 6 3 Model Predictions BMEP 1 bar SOI -87 BTDC PIVC 2.1 bar NG.44 CR 14.5 Case 1 Case 2 Case 3 Case 4 Case 5-2 -2-15 -1-5 5 1 15 2 25 3 35 4 CAD ATDC Case 1 2 3 4 5 TIVC (K) 392 397.5 43 45 48 16 14 12 1 8 6 4 2 Primary Fuel: Natural Gas Ignition Source: Diesel Pilot (1%)
Validation - Cummins Dual-Fuel Engine KIVA-CHEMKIN vs KIVA-CHEMKIN-G AHRR(J/deg) 18 15 12 9 6 Test KIVA-CHEMKIN Case 1 16 14 12 1 8 6 4 AHRR(J/deg) 18 15 12 9 6 Test KIVA-CHEMKIN-G Case 1 16 14 12 1 8 6 4 3-2 -5-4 -3-2 -1 1 2 3 4 5 CAD ATDC 2 3-2 -5-4 -3-2 -1 1 2 3 4 5 CAD ATDC 2 Without flamepropagation combustion is too slow
Summary and Conclusions Current CFD models describe important physical and chemical processes of engine flows, but validation experiments are needed for model development and improvement CFD combustion models capture engine performance and emissions trends and can be combined with optimization tools for design and evaluation of engine design concepts Methodologies are in place for modeling a wide range of combustion regimes: HCCI PCCI Diesel SI flame propagation dual-fuel Multi-component fuel vaporization models can be coupled with advanced chemical kinetics models Surrogate fuel models are under development to represent representative practical fuels (MB, toluene, ethanol..)
Multi-Component fuel modeling Dr. Youngchul Ra Gasoline* * Adapted from J. Farrel, Exxon Mobil Corp. Diesel Aromatics % 16 Sulfur ppm 7.3 Paraffins % 42 Napthenes % 42 Olefins %.3 e l o m C/H ratio - 6.393 r o n o i t u b i r t s i d 1.6 5 gasoline composition 1.4 Discrete g p (mw i ) iso-octane approximation 4 1.2 Single comp approx 1 3.8.6 Continuous f p (I) 2.4 1.2 5 1 15 2 25 3 molecular weight d i s t r i b u t i o n o r m o l
PRF5: Emissions HC at EVO [g.kg-f] 14 12 1 8 6 4 PRF5 Diesel-sgl CO at EVO [g/kg-f] 8 7 6 5 4 3 2 PRF5 Diesel-sgl 2 HC -5-45 -4-35 -3-25 -2 1 CO -5-45 -4-35 -3-25 -2 start of injection command [deg ATDC] start of injection command [deg ATDC].25.2 PRF5 Diesel-sgl NOx [g/kg-f].15.1 LTC regime (EGR=65%) Injection sweep.5. NOx -5-45 -4-35 -3-25 -2 start of injection command [deg ATDC]
Reduced Biodiesel Mechanism Jessica Brakora LLNL methyl butanoate (MB) * Surrogate to represent the large methyl esters in biodiesel MB chemical formula: C 5 H 1 O 2 264 species, 1219 reactions Mechanism Reduction Method SENKIN constant volume analysis Identify key species and reaction paths based on ignition timing Peak concentration: remove species that never exceed specified mole fraction thresholds (-112 species) XSENKPLOT reaction flux analysis: identify key reaction pathways, remove less important paths (-111 species) Adjust rate constants for key reactions to optimize ignition timing *Proc. Comb. Inst., Vol 28, 2/pp.1579-1586 1.E+ 1.E-1 1.E-2 1.E-3 1.E-4 Large flux to mb2j (C 5 H 9 O 2 ); test removal of other paths Ignition Delay [s] LLNL detailed (264 species) rm spec < 1^-1 (199 spec) rm spec < 1^-9 (179 spec) rm spec < 1^-8 (162 spec) No noticeable change in ignition; these species can be removed 1.E-5 6 7 8 9 1 11 12 13 14 Initial Temperature [K] XSENKPLOT reaction paths: MB fuel breakup
Reduced MB Kinetic Mechanism mb ch3o2h (ch3o, ch2o, ch3o2) mb2ooh mb2oo mb2j These species further break down into lower-level carbon species (ch3oco) mb2o mb2ooh4j mp2d MB reduced coupled with ERC n-heptane: 41 species, 15 reactions me2*o c2h5cho mb2ooh4oo ch2cho 1.E+2 P = 4bar, Phi = 1. LLNL detailed me2j*o (ch2co,ch3oco) (c2h5) mb4ooh2*o mp3j2*o ch2co (ch2o) Ignition Delay [ms] 1.E+1 1.E+ 1.E-1 MB reduced (ch3oco, ch2co) 1.E-2 6 7 8 9 1 11 12 13 14 Initial Temperature [K]
Auto-ignition and Detailed Chemistry Models Diesel Gasoline RH + O 2 ROOH First Stage Ignition OH ORO Olefin RH + O 2 HO 2 H 2 O 2 ERC n-heptane mechanism SAE 24-1-558: 34 species H 2 O 2 + (M) OH + OH + (M) Second Stage Ignition ERC PRF Reduced Chemistry SAE 27-1-19, SAE 27-1-165: 47 species
Multi-component Fuel DMC/PRF Models PRF5 fuel pressure [MPa] 14 12 1 8 6 4 2 n-heptane iso-octane 28 24 2 16 12 8 4 HRR [J/deg] P, inj= -45 P, inj= -42 P, inj= -39 P, inj= -33 P, inj= -27 P, inj= -21 P, ic8 P, nc7 HRR, inj= -45 HRR, inj= -42 HRR, inj= -39 HRR, inj= -33-3 -2-1 1 2 3 4 5 CA [deg ATDC] HRR, inj= -27 HRR, inj= -21 HRR, ic8 HRR, nc7