Second Two-Day Meeting on Internal Combustion Engine Simulations Using the OpenFOAM technology, Milan 26-27 th November 216. In-cylinder flows and combustion modeling: application and validation to real and enginelike configurations T. Lucchini, G. D'Errico, A. Della Torre, T. Cerri, A. Maghbouli, E. A. Tahmasebi, L. Sforza, D. Paredi Politecnico di Milano, Department of Energy
Topics 2 In-cylinder flows and combustion modeling using OpenFOAM technology OpenFOAM Lib-ICE CFD methodologies Validation Application Next steps
Lib-ICE 3 Internal combustion engine modeling using the OpenFOAM technology OpenFOAM-x.x.x Engine simulation workflow Mesh generation Development/validation Lib-ICE Fuel-air mixing Spray modeling Library: physical models, mesh management 1 2 3 Engine flows Applications: solvers (cold flow, SI, Diesel, aftertreatment), utilities Combustion Diesel combustion SI combustion
Engine simulation workflow 4 Methodology Diesel engines Mesh management Cold flow Fuel-air mixing Combustion SI engines
Engine simulation workflow 5 Methodology Diesel engines Mesh management Cold flow Fuel-air mixing Combustion SI engines Automatic mesh generation Mesh motion Topological changes Discretization Turbulence models Mesh quality Lagrangian spray Sub-models Nozzle flow Simplified or detailed kinetics Ignition Pollutants
Engine simulation workflow 6 Methodology Diesel engines SI engines Fully integrated approaches Mesh management Automatic mesh generation Mesh motion Topological changes Cold flow Discretization Turbulence models Mesh quality Fuel-air mixing Lagrangian spray Sub-models Nozzle flow Open-Source code Combustion Simplified or detailed kinetics Ignition Pollutants
Diesel Engines
Diesel engines 8 Methodology / model development Engine CAD data Nozzle flow Combustion models Automatic mesh generation Mesh handling Cold flow Spray models
Diesel engines 9 Mesh management Automatic mesh generation a) Main engine data b) Spray oriented block mesh Mesh handling Dynamic layering d) Spray-oriented mesh c) Combustion chamber points fit
Diesel engines 1 Mesh management Automatic mesh generation a) Main engine data b) Spray oriented block mesh Mesh handling Dynamic layering Bowl #1 Bowl #2 Bowl #3 Bowl #4 Bowl #5 Bowl #6 d) Spray-oriented mesh c) Combustion chamber points fit From CAD to SIMULATION: 1 minutes
Average Liquid Length [mm] 18 mm Mixture Fraction [-] 2d computational mesh Diesel Engines: spray modeling 11 Spray A from Engine Combustion Network Spray model setup Injection: blob Breakup: KHRT Evaporation: Spalding CFD setup Turbulence model: standard k-e with modified C 1 Test conditions Fuel n-dodecane Nozzle diameter: 9 mm p inj : 5-15 MPa T amb : 7-12 K r amb : 7.6 22.8 kg/m 3 Baseline (p inj = 15 MPa; T amb = 9 K; r amb = 22.8 kg/m 3 ) Dashed: exp Continuous: calc. Liquid Vapor 24 2 16 12 8 4.16.12.8.4 z = 22.5 mm Exp PoliMI z = 45 mm.. 5. 1. Distance from axis [mm] Parametric variations (T amb, r amb ) Experimental Computed 7 8 9 1 11 Ambient Temperature [K]
Diesel Engines: spray modeling 12 FPT C11 Spray (collaboration with TU/e DI N. Maes and Prof. B. Somers, FPT support) TU/e optical vessel 3D mesh (a) CFD Setup (b) Validation Liquid: exp vs computed extinction profiles Liquid penetration: DBI Vapor penetration: Schlieren Operating conditions Nozzle diameter:.25 mm Ambient temperature: 9 K ANR C1 C2 P inj [MPa] 15 8 15 r amb [kg/m 3 ] 22.8 4 4 Consistent with the engine grid in size and structure Spray Huh-Gosman atomization Pilch-Erdman breakup Turbulence modeling k-e with modified C 1
Diesel Engines: spray modeling 13 FPT C11 Spray (collaboration with TU/e DI N. Maes and Prof. B. Somers, FPT support) TU/e optical vessel 3D mesh (a) CFD Setup (b) Validation Liquid: exp vs computed extinction profiles Liquid penetration: DBI Vapor penetration: Schlieren Operating conditions Nozzle diameter:.25 mm Ambient temperature: 9 K ANR C1 C2 P inj [MPa] 15 8 15 r amb [kg/m 3 ] 22.8 4 4 Consistent with the engine grid in size and structure Spray Huh-Gosman atomization Pilch-Erdman breakup Turbulence modeling k-e with modified C 1
Diesel Engines: spray modeling 14 FPT C11 Spray (collaboration with TU/e DI N. Maes and Prof. B. Somers, FPT support) TU/e optical vessel 3D mesh (a) CFD Setup (b) Validation Liquid: exp vs computed extinction profiles Liquid penetration: DBI Vapor penetration: Schlieren Operating conditions Nozzle diameter:.25 mm Ambient temperature: 9 K ANR C1 C2 P inj [MPa] 15 8 15 r amb [kg/m 3 ] 22.8 4 4 Consistent with the engine grid in size and structure Spray Huh-Gosman atomization Pilch-Erdman breakup Turbulence modeling k-e with modified C 1 Vapor: exp (Schlieren), calc. (mixture fraction threshold)
Diesel Engines: cold flow 15 Steady state flow bench simulations (FPT Industrial, Gilles Hardy) Mesh generation CFD Setup Steady-state flow Turbulence model: standard k-e Porous media acts as a flow straigthner in the computational mesh: computation of the swirl torque Verification of angular momentum conservation across the mesh boundaries. Cartesian, body fitted grids with boundary layer generated with cartesianmesh (cfmesh) or snappyhexmesh (OpenFOAM) This work was part of Davide Paredi MSc Thesis
Diesel Engines: cold flow Steady state flow bench simulations (FPT Industrial, Gilles Hardy) Computed flow field.8.6.4.2 Lift = 5 mm Flow coeff. C.8 d [-].6.4.2 Lift = 1 mm Exp. Calc. 16.8.6.4.2 Swirl number N s [-] from Torque Lift = 5 mm Lift = 1 mm Exp. Calc. This work was part of Davide Paredi MSc Thesis
Diesel Engines 17 Combustion modeling Characteristic time-scale model (CTC) Well-mixed model Medium-duty engine operating at low-load with iso-octane fuel. Passenger-car Diesel engine mrif model N-dodecane spray combustion at constant volume conditions ECN Spray-A experiment Variation of ambient temperature, density, oxygen and injection pressure Kinetic mechanism: Luo et al. (111 species, 467 reactions) Ignition delay Density Inj. Press. O 2 conc. Amb. Temp Computed Experimental
Diesel Engines 18 Combustion modeling Characteristic time-scale model (CTC) Well-mixed model Medium-duty engine operating at low-load with iso-octane fuel. Passenger-car Diesel engine mrif model N-dodecane spray combustion at constant volume conditions ECN Spray-A experiment Variation of ambient temperature, density, oxygen and injection pressure Kinetic mechanism: Luo et al. (111 species, 467 reactions) Flame lift-off Density Inj. Press. O 2 conc. Amb. Temp Computed Experimental
Diesel Engines: validation 19 Spray B in Engines (from ECN, collaboration with Dr. Eagle, Dr. Malbec, Dr. Musculus) Optically accessible engine with three hole injector CFD setup and tested points Validation: non reacting Liquid length Vapor distribution CALC. CALC. Nozzle design very similar to spray A: possibility to perform same studies in engine and constantvolume vessel. This work is part of Amin Maghbouli PhD project Same CFD setup used in Spray A simulations, two different combustion models tested: mrif and well-mixed. Tested conditions T@SOI: 8-1 K P inj [bar]: 5 15 bar O 2 [%]: 13-21 % r @ SOI: 15.2 22.8 kg/m 3 Fuel: n-dodecane EXP EXP
Diesel Engines: validation 2 Spray B in Engines (from ECN, collaboration with Dr. Eagle, Dr. Malbec, Dr. Musculus) Optically accessible engine with three hole injector CFD setup and tested points Validation: reacting Nozzle design very similar to spray A: possibility to perform same studies in engine and constantvolume vessel. This work is part of Amin Maghbouli PhD project Automatically generated grid Same CFD setup used in Spray A simulations, two different combustion models tested: mrif and well-mixed. Tested conditions T@SOI: 8-1 K P inj [bar]: 5 15 bar O 2 [%]: 13-21 % r @ SOI: 15.2 22.8 kg/m 3 Fuel: n-dodecane
Diesel Engines: validation 21 EU 6 FPT Engine MD (support from FPT, DI Gilles Hardy) Cylinder pressure validation (mrif and well-mixed) 75% load, 15% EGR, 2 injections 4 14 Exp. m-rif Exp. 4 Grid automatically generated 25 8 2 6 15 4 1 2 5-5 -3-1 1 3 Crank Angle [deg] Initial conditions from 1D simulations 5 Well-mixed 35 7 Pressure [bar] 3 m-rif 8 35 1 [rpm] 4 9 Well-mixed 12 Pressure [bar] 8 25% load, 2% EGR, 3 injections Heat release rate [J/kg] bmep/bmepmax 1 3 6 25 5 2 4 15 3 1 2 5 1-5 -3-1 1 3 Crank Angle [deg] 5 Heat release rate [J/kg] Five operating points: with different EGR levels
Diesel Engines: validation 22 EU 6 FPT Engine MD (support from FPT, DI Gilles Hardy) Cylinder pressure validation (mrif and well-mixed) 1% load, 35% EGR, 3 injections 4 7 Exp. m-rif Exp. 3 4 Grid automatically generated 25 4 2 3 15 2 1 1 5-5 -3-1 1 3 Crank Angle [deg] Initial conditions from 1D simulations 5 Pressure [bar] 35 5 [rpm] 4 16 Well-mixed 6 Pressure [bar] 8 Full load, % EGR, 1 injection Heat release rate [J/kg] bmep/bmepmax 1 m-rif Well-mixed 14 35 12 3 1 25 8 2 6 15 4 1 2 5-5 -3-1 1 3 Crank Angle [deg] 5 Heat release rate [J/kg] Five operating points: with different EGR levels
Cyl. Pressure [bar] Heat Release Rate [J/deg] Diesel Engines: validation 23 PCCI engine (support from FPT, DI Gilles Hardy) PCCI combustion conditions Combustion model: Well-mixed High EGR rate (5%) 3 rpm 14 12 Exp. 35 3 Exp. Computational mesh generated automatically with the Polimi tool 1 8 6 Well-mixed 25 2 15 Well-mixed 4 1 2 5-3 -1 1 3 5-3 -1 1 3 5 Crank Angle [deg] Crank Angle [deg] Diesel fuel: n-dodecane Kinetic mechanism: Faravelli et al. (11 species, 3 reactions)
PCCI combustion conditions High EGR rate (5%) 3 rpm Computational mesh generated automatically with the Polimi tool Diesel Engines: validation 24 PCCI engine (support from FPT, DI Gilles Hardy)
NOx (normalized) soot (normalized) Diesel Engines: validation 25 Heavy Duty engines (support from GE Global Research, Dr. Pasunurthi) Automatic mesh generation Operating conditions SOI variation Soot/NO x trade-off Dual-fuel combustion Diesel: mrif Dual-fuel: PaSR NOx: Zeldovich, soot: Moss P cyl /P max Conventional Diesel Combustion Engine 1: SOI variation SOI #1, Exp. SOI #1, Calc. SOI #2, Exp. SOI #2, Calc. SOI #3, Exp. SOI #3, Calc. -2 2 4 Crank Angle [deg] p/p max Engine 2: pollutant prediction OP.1-2 -1 1 2 3 Crank Angle [deg] 1.2 1.8.6.4.2 Exp. Calc. p/p max OP.3 Calc. Exp. OP.1 OP.2 OP.3-2 -1 1 2 3 Crank Angle [deg] 1.4 1.2 1.8.6.4.2 Exp. Calc.
Diesel Engines validation 26 Heavy Duty engines (support from GE Global Research, Dr. Pasunurthi) Setup Reduced mechanism for Diesel (C 7 H 16 ) and natural gas Combustion model: PaSR Dual-fuel combustion RoHR/RoHR max OP 1 Validation Exp. Calc. RoHR/RoHR max OP 2 Exp. Calc. Pollutants NO x : Zeldovich Soot: Moss Two operating points: same air/fuel ratio, different amount of injected Diesel fuel -2 2 4-2 2 4 Crank Angle [deg] Crank Angle [deg] 1.2 Exp. 1 Calc..8 NO x / NO x,max.6.4.2 OP 1 OP 2
Internal nozzle flow modeling (RANS, cavitation) Spray C & Spray D These two single hole injectors are defined by ECN to compare the effect of geometry on flow behavior. Spray C Spray D Spray C Diesel Engines 27 Next steps 12 11.5 11 1.5 1 9.5 9 ECN contributors 11.72 spray results C mdot (g/s) spray D mdot (g/s) 1.92 1.7 Experiments 1.51 POLIMI Presssure across the hole axis and more mrif with multiple injections Tabulated kinetics: Homogeneous reactors Unsteady RIF blockmesh grid Spray C k factor=, sharp edge Spray D k factor=1.5, rounded edge Spray C cavitated Spray D Not cavitated Transported PDF combustion modeling: Eulerian Stochastic fields with tabulated kinetics! This work is part of Ehsan Tahmasebi PhD project
SI Engines
SI Engines 29 Methodology / model development Engine CAD data Nozzle flow Combustion models Automatic mesh generation Mesh handling Cold flow Spray models
SI Engines: cold flow 3 Darmstadt optical engine (collaboration with Dr. B. Bohm and DI C. P. Ding) 1 2 3 Velocity Magnitude (m/s) Four-valve engine, fully optically accessible PIV measurement techniques: low repetition rate planar PIV high-speed PIV (HS-PIV) stereoscopic PIV (SPIV) tomographic PIV (TPIV).
Lift [mm] SI Engines: cold flow 31 Full cycle SI: Darmstadt optical engine automatic mesh generation STL + data Geometry-oriented block structured mesh snappyhexmesh 1 7.5 Exhaust Intake 5 2.5 18 36 54 72 Crank Angle [deg] Bore 8 mm Stroke 6 mm Conrod 16 mm The user provides combustion chamber geometry and data (bore, stroke, valve lifts) A python program automatically recognizes the direction of ducts, cylinder and valves and generates a geometryoriented background grid snappyhexmesh is then run using the geometry-oriented background mesh
Lift [mm] SI Engines: cold flow 32 Full cycle SI: Darmstadt optical engine automatic mesh generation STL + data Geometry-oriented block structured mesh snappyhexmesh 1 7.5 Exhaust Intake 5 2.5 18 36 54 72 Crank Angle [deg] Bore 8 mm Stroke 6 mm Conrod 16 mm The user provides combustion chamber geometry and data (bore, stroke, valve lifts) A python program automatically recognizes the direction of ducts, cylinder and valves and generates a geometryoriented background grid snappyhexmesh is then run using the geometry-oriented background mesh
SI Engines: cold flow 33 Full cycle SI: Darmstadt optical engine mesh management Full-cycle simulations: - Multiple meshes - Mesh to mesh interpolation strategy. Duration of each mesh: - User defined + quality criteria Initial mesh at Crank angle q q = q curr Generate a new mesh with snappyhexmesh q curr = q Move mesh for Dq Mesh quality and duration satisfied? NO q curr = q curr + Dq YES Move surface geometry to current crank angle q curr NO YES q curr = q end? End of meshing
SI Engines: cold flow 34 Full cycle SI: Darmstadt optical engine case setup Models and boundary conditions Engine geometry data Bore Stroke 86 mm 86 mm Compression ratio 8.5 IVO IVC EVO EVC Speed Combustion 325 CAD 485 CAD 15 CAD 345 CAD 85 rpm no Boundary conditions Unsteady boundary conditions (from exp. data) imposed at inlet and outlet ports. CFD models Second-order discretization (TVD) Turbulence model: standard k-e
SI Engines: cold flow 35 Full cycle SI: Darmstadt optical engine case setup Computational mesh Mesh layouts 4 mm mesh size in the ducts region; 2 mm mesh size in the cylinder and valve region; local refinement up to 1 mm close to cylinder head, piston and liner boundaries; local refinement up to.25 mm close to the valves boundaries. a) Cartesian mesh (automatically generated + snappyhexmesh) b) Flow-oriented mesh (automatically generated + Polimi tool + snappyhexmesh)
SI Engines: cold flow 36 Full cycle SI: Darmstadt optical engine validation
U y [m/s] SI Engines: cold flow 37 Full cycle SI: Darmstadt optical engine validation x and y velocity components along four different measurement lines located at different distances from the cylinder head. 45 CAD mid-intake Similar behavior between the two grids, flow dominated by the incoming air jet Experimental Cartesian Flow-oriented Y = mm Y = -1 mm Y = -2 mm Y = -3 mm -1 y -2-3 x Flow-oriented Cartesian Experimental
U x [m/s] SI Engines: cold flow 38 Full cycle SI: Darmstadt optical engine validation x and y velocity components along four different measurement lines located at different distances from the cylinder head. 45 CAD mid-intake Similar behavior between the two grids, flow dominated by the incoming air jet Experimental Cartesian Flow-oriented Y = mm Y = -1 mm Y = -2 mm Y = -3 mm -1 y -2-3 x Flow-oriented Cartesian Experimental
U y [m/s] SI Engines: cold flow 39 Full cycle SI: Darmstadt optical engine validation x and y velocity components along four different measurement lines located at different distances from the cylinder head. 54 CAD BDC, intake Better prediction of the flow oriented grid: Vortex location, distribution of the two main streams Experimental Cartesian Flow-oriented Y = mm Y = -1 mm Y = -2 mm Y = -3 mm -1 y -2-3 x Flow-oriented Cartesian Experimental
U x [m/s] SI Engines: cold flow 4 Full cycle SI: Darmstadt optical engine validation x and y velocity components along four different measurement lines located at different distances from the cylinder head. 54 CAD BDC, intake Better prediction of the flow oriented grid: Vortex location, distribution of the two main streams Experimental Cartesian Flow-oriented Y = mm Y = -1 mm Y = -2 mm Y = -3 mm -1 y -2-3 x Flow-oriented Cartesian Experimental
U y [m/s] SI Engines: cold flow 41 Full cycle SI: Darmstadt optical engine validation x and y velocity components along four different measurement lines located at different distances from the cylinder head. 66 CAD mid-compression Flow-oriented grid better describe the tumble vortex structure Experimental Cartesian Flow-oriented Y = mm Y = -1 mm Y = -2 mm Y = -3 mm -1 y -2-3 x Flow-oriented Cartesian Experimental
U x [m/s] SI Engines: cold flow 42 Full cycle SI: Darmstadt optical engine validation x and y velocity components along four different measurement lines located at different distances from the cylinder head. 66 CAD mid-compression Flow-oriented grid better describe the tumble vortex structure Experimental Cartesian Flow-oriented Y = mm Y = -1 mm Y = -2 mm Y = -3 mm -1 y -2-3 x Flow-oriented Cartesian Experimental
SI Engines: GDI fuel-air mixing 43 Stratified engine (collaboration with IM-CNR, Ing. Sementa, Ing. Montanaro) Optically accessible GDI engine Injection pressure [bar] Operating points SOI [ BTDC] Charge stratification 1 1 6 High 2 1 11 Low 3 6 6 High 4 6 12 Low bmep = 7.2 bar, SA = 13 BTDC; l = 1.15 (lean)
SI Engines: GDI fuel-air mixing 44 Stratified engine Optically accessible GDI engine: optical/computed data correlation Fuel m.f. P inj = 1 bar; SOI = 11 BTDC Soot chemiluminescence P inj.12 = 1 bar; SOI = 11 BTDC Fuel m.f. P inj = 6 bar; SOI = 6 BTDC Soot chemiluminescence P inj.12 = 6 bar; SOI = 7 BTDC Possible sources of soot: Rich pockets Wall-film Correlated with optical soot chemiluminescence Wall-film l distribution Wall-film l distribution
SI Engines: combustion 45 Modeling Comprehensive combustion model Detailed description of the flame kernel growth process via Lagrangian approach and suitable sub-models (breakdown, electrical circuit, wrinking) Coherent flamelet model for flame propagation in the Eulerian phase (gas) Strict coupling between Eulerian and Lagrangian phases Sub-models Secondary circuit energy transfer R p R s L p L s Spark gap hq el r i Breakdown T i (>1 K) T u (3-6 K) Flame surface density spk N f i 1 i S i V cell Lagrangian particles for the spark-channel Mass and energy equations solved for the flame kernel particles Particles convected by the flow, possibility to predict restrike Possibility to predict the effects of local flow, mixture conditions, turbulence, electrical circuit properties (voltage, current). Prediction of misfire also possible.
SI Engines: combustion 46 Assessment and validation Applications Results Pressurized vessels Multi-ignition systems Fan-generated flow velocity and turbulence fields at the spark-gap Initial combustion stage (Lagrangian Eulerian coupling) Fully turbulent flame (Eulerian model only) This work is part of Lorenzo Sforza PhD project
Enflammed volume [cm 3 ] Enflammed volume [cm 3 ] SI Engines: combustion 47 Experimental validation: Herweg and Maly Engine Flow field, flame propagation and plasma temperature details 2 1.6 1.2 Calc., 1 rpm, l = 1. Calc., 1 rpm, l = 1. Exp., 1 rpm, l = 1. Exp., 1 rpm, l = 1. 2 1.6 1.2 Calc., 1 rpm, l = 1. Calc., 3 rpm, l = 1. Calc., 125 rpm, l = 1. Exp., 1 rpm, l = 1. Exp., 3 rpm, l = 1. Exp., 125 rpm, l = 1..8 Peripheral.8.4.4 Central.2.4.6.8 1 Time after spark [ms] This work is part of Lorenzo Sforza PhD project.2.4.6.8 1 Time after spark [ms]
Conclusions 48 CFD modeling of in-cylinder phenomena at Polimi with OpenFOAM Detailed models continuously validated and improved: Fundamental studies Applied research Consolidated methodologies, currently applied in the context of industrial collaborations but there is still a lot to do!
Thanks for your attention!