390 likes | 489 Views
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
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: