390 likes | 508 Views
This document details the validation process for Eulerian spray modeling techniques applied to the characterization of diesel fuel injection. The study, conducted by Milan Vujanovic at the University of Zagreb's Faculty of Mechanical Engineering, focuses on Nozzle D (205 µm diameter) under various conditions, including rail pressures of 500 and 1200 bar, gas chamber pressures of 54-800 bar, and gas temperatures of 900 K. Key findings highlight comparisons between modeled penetration rates for liquid and vapor phases against experimental data, accounting for turbulence and dissipation coefficients in the context of power engineering and energy management.
E N D
University of Zagreb Faculty of Mechanical Engineering and Naval Architecture Department of Energy, Power Engineering and Environment Chair of Power Engineering and Energy Management Status – Validation of Eulerian Spray Modelling Milan Vujanovic May, 2006
Validation: I-Level project Version v8.5006 vs. Version v8.5014 Nozzle D – 205 micron diameter Rail pressure – 500 bar Gas chamber pressure – 72 bar Gas temperature in chamber - 900 K
Test Case - Nozzle D – 205 micron diameter Experimental data – injection rate:
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00261.0e-06 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results 8.5006 8.5014
Validation: I-Level project Impact of initial k and epsilon values Nozzle D – 205 micron diameter Rail pressure – 500 bar Gas chamber pressure – 72 bar Gas temperature in chamber - 900 K Case 1_1 Turb. kin. energy – 10 m2/s2 Turb. length scale – 2e-05 m Turb. diss. rate – 259 808 m2/s3 Case 6_1 Turb. kin. energy – 250 m2/s2 Turb. length scale – 2e-05 m Turb. diss. rate – 3.247e+07 m2/s3
Penetration for liquid phase and vapour phase compared with experimental results Case 1_1 Turb. kin. energy – 10 m2/s2 Turb. length scale – 2e-05 m Turb. diss. rate – 259 808 m2/s3 Case 6_1 Turb. kin. energy – 250 m2/s2 Turb. length scale – 2e-05 m Turb. diss. rate – 3.247e+07 m2/s3
Validation: I-Level project Impact of constant cε2 Nozzle D – 205 micron diameter Rail pressure – 500 bar Gas chamber pressure – 72 bar Gas temperature in chamber - 900 K The constant cε2 in the transport equation for the dissipation rate of the turbulent kinetic energy was set tocε2= 1.8 instead cε2=1.92
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00261.0e-06 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00261.0e-06 / 5.0e-07 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8
Validation: I-Level project Impact of constant cε2 Nozzle D – 205 micron diameter Rail pressure – 1200 bar Gas chamber pressure – 72 bar Gas temperature in chamber - 900 K The constant cε2 in the transport equation for the dissipation rate of the turbulent kinetic energy was set tocε2= 1.8 instead cε2=1.92
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00265.0e-07 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8
Validation: I-Level project Impact of constant cε2 Nozzle D – 205 micron diameter Rail pressure – 500 bar Gas chamber pressure – 54 bar Gas temperature in chamber - 900 K The constant cε2 in the transport equation for the dissipation rate of the turbulent kinetic energy was set tocε2= 1.8 instead cε2=1.92
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00261.0e-06 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8
Validation: I-Level project Impact of constant cε2 Nozzle D – 205 micron diameter Rail pressure – 800 bar Gas chamber pressure – 54 bar Gas temperature in chamber - 900 K The constant cε2 in the transport equation for the dissipation rate of the turbulent kinetic energy was set tocε2= 1.8 instead cε2=1.92
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00265.0e-07 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 4.5
Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8
Validation: I-Level project Impact of constant cε2 Nozzle D – 205 micron diameter Rail pressure – 1200 bar Gas chamber pressure – 54 bar Gas temperature in chamber - 900 K The constant cε2 in the transport equation for the dissipation rate of the turbulent kinetic energy was set tocε2= 1.8 instead cε2=1.92
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00265.0e-07 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8
Validation: I-Level project k – zeta – f turbulence model Nozzle D – 205 micron diameter Rail pressure – 500 bar Gas chamber pressure – 72 bar Gas temperature in chamber - 900 K
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00261.0e-06 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results k – epsilon k –zeta - f
Validation: I-Level project k – zeta – f turbulence model Nozzle D – 205 micron diameter Rail pressure – 1200 bar Gas chamber pressure – 54 bar Gas temperature in chamber - 900 K
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.00265.0e-07 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results k – epsilon k –zeta - f
Validation: I-Level project Calculation with nozzle interface Coupling internal nozzle flow simulation and initialisation of spray calculation Nozzle D – 205 micron diameter Rail pressure – 500 bar Gas chamber pressure – 72 bar Gas temperature in chamber - 900 K Using the data of the two phase flow calculation inside the nozzle as a start and boundary condition for Eulerian spray calculation
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.0026 1.0e-06 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results without nozzle interface with nozzle interface
Validation: I-Level project Calculation with nozzle interface Coupling internal nozzle flow simulation and initialisation of spray calculation Nozzle D – 205 micron diameter Rail pressure – 1200 bar Gas chamber pressure – 72 bar Gas temperature in chamber - 900 K Using the data of the two phase flow calculation inside the nozzle as a start and boundary condition for Eulerian spray calculation
Calculation settings Time discretisation: Upto Time [s] Δt upto 1.0e-6 2.5e-08 upto 1.0e-4 2.5e-07 upto 2.0e-4 5.0e-07 upto 0.0026 1.0e-06 The liquid → Diesel →T=373 K Eulerian spray with 6 phases Primary brake-up model: Dies.Core Injection Secondary brake-up model: Wave model Evaporation model: Abramzon-Sirignano model Turbulent dispersion coefficient = 6
Penetration for liquid phase and vapour phase compared with experimental results without nozzle interface with nozzle interface
University of Zagreb Faculty of Mechanical Engineering and Naval Architecture Department of Energy, Power Engineering and Environment Chair of Power Engineering and Energy Management The end
2nd phase of validation: I-Level project Nozzle D – 205 micron diameter Experimental data – injection rate:
Test Case: I-Level project Nozzle D – 205 micron diameter Experimental data – injection rate: