Digital Shaping and Optimization of Fuel Injection Pattern for a Common Rail Automotive Diesel Engine through Numerical Simulation European GT Conference 2017 - Frankfurt am Main Politecnico di Torino: Sapio, F., Piano, A., Millo, F. General Motors GPS: Pesce, C. F. Oct - 2017
Agenda Introduction Injector Modeling Combustion Modeling Injection Strategies Optimization Conclusions Ongoing Work 2
Agenda Introduction Injector Modeling Combustion Modeling Injection Strategies Optimization Conclusions Ongoing Work 3
Aim of the Work Exploiting the fuel injection pattern and shaping potential in reducing BSFC and combustion noise, without exceeding the target emissions level, by means of a virtual test-rig. NVH Engine type Displacement Bore stroke DI Turbocharged Diesel EURO6 1598 cm3 79.7 mm 80.1 mm Compression ratio 16:1 CARBON DIOXIDE CO 2 PARTICULATE MATTER NITROGEN OXIDE NOx Turbocharger Fuel injection system Maximum power and torque Single-stage with VGT Common Rail 100 kw @ 4000rpm 320 Nm @ 2000rpm Source: http://www.autosupermarket.it/magazine/gm-powertrain-europe-realta-italiana/ Johnson, J., Bagley, S., Gratz, L., and Leddy, D., "A Review of Diesel Particulate Control Technology and Emissions Effects - 1992 Horning Memorial Award Lecture," SAE Technical Paper 940233, 1994, doi:10.4271/940233 4
Virtual Test Rig Source: http://www.boschautoparts.com 5
Agenda Introduction Injector Modeling Combustion Modeling Injection Strategies Optimization Conclusions Ongoing Work 6
Experimental Setup In order to build a 1D-CFD model of a CR injector, 2 crucial issues need to be addressed. Internal geometry detection Extensive dataset of experimental injection rate 3D Computed Tomography STS Injection Analyzer Source: http://www.boschautoparts.com 7
Injector Model Results EMI Single Injections EMI Curves Rail Pressure = 400 bar Rail Pressure = 600 bar Rail Pressure = 1000 bar Rail Pressure = 1200 bar 8
Injector Model Results IR Single\Double + Close Pilot Rail Pressure = 400 bar Rail Pressure = 1000 bar 3 Pilots + Main + After Rail Pressure = 400 bar 9
Injection Rate Map (IRM) Engine Model CPU Increase Factor due to adding Detailed Injection* 4 Detailed injectors 17.3 1 Master injector + 3 slaves 6.3 Injection Rate Map 1.1 Advantages Best trade-off between accuracy and CPU time Drawbacks Unable to detect pulse interaction *Source: courtesy of Gamma Technologies Inputs Injected Mass per Pulse or ET Rail Pressure Output Injection Rate (Interpolated) 10
Agenda Introduction Injector Modeling Combustion Modeling Injection Strategies Optimization Conclusions Ongoing Work 11
Bmep [bar] DIPulse Combustion Model Predictive Multi-zone Main Unburned Zone (MUZ) Spray Burned Zone (SBZ) Spray Unburned Zone (SUZ) Calibrated on 28 engine operating points Part Load Points Full Load Points *Source: GT Power Engine Performance User s Manual DIPulse Calibration Parameters* 25 20 Entrainment Rate Multiplier 0.95 2.8 15 10 5 0 1000 1500 2000 2500 3000 3500 4000 Engine Speed [rpm] Ignition Delay Multiplier 0.3 1.7 Premixed Combustion Rate Multiplier 0.05 2.5 Diffusion Combustion Rate Multiplier 0.4 1.4 12
DIPulse Calibration Results Pressure Measured Pressure Predicted Burn Rate Measured Burn Rate Predicted Injection Rate 13
DIPulse Calibration Results Engine map Results in terms of main global indexes are shown, related to the whole engine map (337 operating points). Good agreement with experimental data is obtained. ± 5 % ± 2 deg Avg IMEP error < 5% Avg Pmax error < 5 bar Avg CA of Pmax error < 2 deg Avg CA of MFB50 error < 2 deg ± 5 bar ± 2 deg 14
DIPulse Calibration Results Calibration points Results in terms of emissions are shown, related to the 28 chosen calibration points. EGR sweep points EGR sweep results in terms of NOx emissions are also shown, related to 7 PL engine operating points. 15
Agenda Introduction Injector Modeling Combustion Modeling Injection Strategies Optimization Conclusions Ongoing Work 16
Model Setup Detailed Configuration 4 cylinder based Detailed intake/exhaust geometry Turbocharger EGR circuit Advantages Accurate Drawbacks Computationally expensive 17
Model Setup Simplified Configuration (TPA) Single cylinder based Imposed intake/exhaust pressures and temperatures Imposed intake air residual fraction (EGR rate) Injection Rate Map Injection Timing Controller Injection Number Controller 18
BSFC [fraction] Computational Time [min] Base Model - Validation Simplified configuration achievements Similar results compared with the detailed model *(CPU Intel Xeon E31245 @ 3.30GHz, 4 core) Lower computational time* X 1 6 Simplified Detailed Simplified Detailed 1,1 1,05 1 0,95 0,9 0,85 0,8 0,75 5 4,5 4 3,5 3 2,5 2 1,5 1 0,5 0,7 1000x2 1000x5 1500x2 1500x5 1500x8 2000x2 2000x6 2000x8 2000x12 0 1000x2 1000x5 1500x2 1500x5 1500x8 2000x2 2000x6 2000x8 2000x12 Operating Point [RPMxBMEP] Operating Point [RPMxBMEP] 19
Simulated [db] Combustion Noise in GT-SUITE Evaluation procedure Validation In-cylinder pressure signal is taken during the simulation Engine Map (337 operating points) Cylinder Pressure FFT is performed A-filter + CAVTAB-filter attenuation is introduced Overall octave bands content in db is evaluated ± 0.3 db Experimental [db] 20
Bmep [bar] y-value DoE Approach Operating Key Points 1500 rpm x 2 bar 1500 rpm x 5 bar 2000 rpm x 8 bar 25 20 15 10 5 10 8 6 4 A DoE approach is adopted 1. Latin Hypercube 20 000 experiments per case: exploration of the results space 0 1000 1500 2000 2500 3000 3500 4000 Engine Speed [rpm] Input Variables Number of Injections 2 0 0 2 4 6 8 10 x-value Output Variables 2. Full Factorial 500 000 experiments per case: starting from Latin Hypercube results, search refinement Energizing Time Dwell Time Rail Pressure Injection Timing (SOI Main) EGR Rate Brake Specific Fuel Consumption (BSFC) Brake Specific NOx (BSNOx) Combustion Noise (CN) 21
CNF DoE Post Processing Normalized Output Variables Target emission level is set to the baseline value Normalized BSFC Normalized BSNOx BSFC BSFC t BSNO x BSNO x,t Normalized CN Factor 10 db db t 20 ± 5% 1 1 BSFC 22
Results 1500x2 (RPMxBMEP) 1500 x 2 N-BSFC N-BSNOx N-CNF 1 Pil 0.97 1.00 1.06 2 Pil 0.96 0.99 0.99 3 Pil 0.96 0.98 0.97 4 Pil 0.96 0.95 0.73 5 Pil 0.97 0.99 0.85 1,3 1 Pil 2 Pil 3 Pil 4 Pil 5 Pil 1,2 1,1 1,0 0,9 0,8 0,7 0,6 0,5 N-BSFC N-BSNOx N-CNF 23
Results 1500x5 (RPMxBMEP) 1500 x 5 N-BSFC N-BSNOx N-CNF 1 Pil 0.97 0.95 1.26 2 Pil 0.96 1.01 1.03 3 Pil 0.95 0.99 0.93 4 Pil 0.95 1.00 0.84 5 Pil 0.95 1.02 0.73 1,3 1 Pil 2 Pil 3 Pil 4 Pil 5 Pil 1,2 1,1 1,0 0,9 0,8 0,7 0,6 0,5 N-BSFC N-BSNOx N-CNF 24
Results 2000x8 (RPMxBMEP) 2000 x 8 N-BSFC N-BSNOx N-CNF 1 Pil 0.97 0.99 1.10 2 Pil 0.96 1.05 0.85 3 Pil 0.96 0.96 0.90 4 Pil 0.97 0.95 0.79 5 Pil 0.98 0.98 0.80 1,3 1 Pil 2 Pil 3 Pil 4 Pil 5 Pil 1,2 1,1 1,0 0,9 0,8 0,7 0,6 0,5 N-BSFC N-BSNOx N-CNF 25
Agenda Introduction Injector Modeling Combustion Modeling Injection Strategies Optimization Conclusions Ongoing Work 26
Conclusions Injection Strategies Optimization A methodology for developing a 1D-CFD virtual engine test rig was presented. The virtual test rig includes the injection system model, the engine model with a calibrated predictive combustion model and the innovative control strategy of the ECU for the injection pattern. A large number of fuel injection patterns were tested and the best ones in terms of BSFC and CN reduction, without exceeding the target NOx emissions level, were selected. Some common trends can be highlighted at the end of this work: Compact multi-injection patterns High Rail Pressure Short Dwell Times Progressive Burn Rate Achievements: improvements in BSFC up to 5%, CNF up to 30%, without exceeding in NOx emissions 27
Agenda Introduction Injector Modeling Combustion Modeling Injection Strategies Optimization Conclusions Ongoing Work 28
CN [db] Multi-Objective Results (1500x5) Engine Operating Point 1500 X 5 Complete Dataset Cases tested by the optimization process in GT- SUITE Constrained Dataset BSNOx Constraint added on the complete dataset Complete Dataset Constrained Dataset Pareto Front DoE Baseline Baseline Baseline injection pattern configuration The RED area represents the usable area where: BSFC actual BSFC baseline CN actual CN baseline BSFC [g/kwh] 29
Normalized value [-] Normalized value [-] Normalized value [-] Multi-Objective Resume Engine Operating Point - 1500 X 2 Engine Operating Point - 1500 X 5 Engine Operating Point - 2000 X 8 1,0 0,9 0,8 1,0 0,9 0,8 1,0 0,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0,0 BSFC CN Baseline 1 1 DoE - Optimum 0,962 0,72 MO - Minimum BSFC 0,955 0,66 MO - Minimum CN 0,961 0,43 0,0 BSFC CN Baseline 1 1 DoE - Optimum 0,945 0,69 MO - Minimum BSFC 0,932 0,63 MO - Minimum CN 0,943 0,39 0,0 BSFC CN Baseline 1 1 DoE - Optimum 0,936 0,58 MO - Minimum BSFC 0,928 0,58 MO - Minimum CN 0,936 0,50 By means of GA is possible to reduce computational time of about 2 orders of magnitude with respect to a Full Factorial DoE approach. 30
Acknowledgements Prof. Federico MILLO Eng. Andrea PIANO Dr. Francesco C. PESCE 31
Digital Shaping and Optimization of Fuel Injection Pattern for a Common Rail Automotive Diesel Engine through Numerical Simulation European GT Conference 2017 - Frankfurt am Main Politecnico di Torino: Sapio, F., Piano, A., Millo, F. General Motors GPS: Pesce, C. F. Oct - 2017