1 / 39

Status – Validation of Eulerian Spray Modelling

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

magee
Download Presentation

Status – Validation of Eulerian Spray Modelling

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 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

  2. 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

  3. Test Case - Nozzle D – 205 micron diameter Experimental data – injection rate:

  4. 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

  5. Penetration for liquid phase and vapour phase compared with experimental results 8.5006 8.5014

  6. 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

  7. 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

  8. 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

  9. 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

  10. Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8

  11. 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

  12. Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8

  13. 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

  14. 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

  15. Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8

  16. 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

  17. 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

  18. Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8

  19. 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

  20. 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

  21. Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8

  22. 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

  23. 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

  24. Penetration for liquid phase and vapour phase compared with experimental results cε2=1.92 cε2=1.8

  25. 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

  26. 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

  27. Penetration for liquid phase and vapour phase compared with experimental results k – epsilon k –zeta - f

  28. 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

  29. 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

  30. Penetration for liquid phase and vapour phase compared with experimental results k – epsilon k –zeta - f

  31. 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

  32. 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

  33. Penetration for liquid phase and vapour phase compared with experimental results without nozzle interface with nozzle interface

  34. 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

  35. 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

  36. Penetration for liquid phase and vapour phase compared with experimental results without nozzle interface with nozzle interface

  37. 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

  38. 2nd phase of validation: I-Level project Nozzle D – 205 micron diameter Experimental data – injection rate:

  39. Test Case: I-Level project Nozzle D – 205 micron diameter Experimental data – injection rate:

More Related