Modelling of droplet heating and evaporation in computational fluid dynamics codes

# Modelling of droplet heating and evaporation in computational fluid dynamics codes

## Modelling of droplet heating and evaporation in computational fluid dynamics codes

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -
##### Presentation Transcript

1. Modelling of droplet heating and evaporation in computational fluid dynamics codes Sergei SAZHIN*, Irina SHISHKOVA** , Vladimir LEVASHOV **, Morgan HEIKAL* *Sir Harry Ricardo Laboratory, School of Environment and Technology, Faculty of Science and Engineering, University of Brighton, Brighton BN2 4GJ, UK ** Low Temperature Departments, Moscow Power Engineering Institute, Moscow 111250, Russia

2. Presentation overview • • INTRODUCTION • • LIQUID PHASE MODELS • • GAS PHASE MODELS • •TEST CASES • • MODELLING VERSUS EXPERIMENTS • • KINETIC MODELLING • • RELATED DEVELOPMENTS

3. INTRODUCTION:Processes in Diesel engines • •Formation of a liquid fuel spray • • Fuel spray penetration • • Droplet break-up • • Heating of fuel droplets • • Evaporation of fuel droplets • • Ignition of fuel vapour / air mixture

4. Ignition (experiment) 160 MPa injection into 10 MPa gas 140 MPa injection into 10 MPa gas 100 MPa injection into 10 MPa gas

5. LIQUID PHASE MODELS The Infinite Thermal Conductivity (ITC) models are based on the energy balance equation of the droplet as a whole. The solution to this equation: ,

6. Conduction limit. Effective conductivity

7. where θR is the radiation temperature, Rd is the droplet radius, θRcan be assumed equal to the external temperature Text in the case of an optically thin gas in the whole domain.

8. Analytical solution for h=const

9. Numerical algorithms • numerical algorithm based on the analytical solution (analytical solution at the end of the time step is considered as the initial condition for the next time step) • numerical solution of the discretised heat conduction equation (fully implicit approach) • numerical solution based on the assumption of no temperature gradient inside the droplet (conventional approach currently used in CFD codes)

10. Main results (Liquid phase models) • The numerical algorithm based on the analytical solution is recommended when a compromise between high accuracy and CPU requirements is essential. • The contribution of thermal radiation can be taken into account via the simplified form of the radiation term.

11. GAS PHASE MODELS Model 4 (Abramzon and Sirignano, 1989): , .

12. TEST CASE 1 (Zero-dimensional code) , .

13. Plots for Tg0= 880K, pg0=3 MPa, Td0= 300 K, Rd0=10 μm and vd0=1 m/s. The overall volume of injected liquid fuel was taken equal to 1 mm3, and the volume of air, where the fuel was injected, was taken equal to 883 mm3. The results were obtained based on the ETC model and using seven gas phase models. Text=2000K

14. The same as the previous figure but using three liquid phase numerical algorithms: the algorithm based on the analytical solution of the heat conduction equation inside the droplet (1), the algorithm based on the numerical solution of the heat conduction equation inside the droplet (2), the algorithm based on the assumption that the thermal conductivity inside droplets is infinitely large (3).

15. TEST CASE 2 (KIVA 2 CFD code) , .

16. The total autoignition delay times observed experimentally and computed using the customised version of the KIVA 2 CFD code at three initial in-cylinder pressures. The initial injected liquid fuel temperature was assumed equal to 375 K. The injection pressure was equal to 160 MPa. The modified WAVE model, two liquid phase models (ETC and ITC) and two gas phase models (model 0 and the AS model) were used for computations.

17. Main results(Hydrodynamic Heating and Evaporation Models)The choice of the gas phase model is essential for predicting the evaporation time The choice of the liquid phase model is essential for predicting the initial temperatures and autoignition delaySazhin, S.S., Kristyadi, T., Abdelghaffar, W.A. and Heikal, M.R. (2006) Models for fuel droplet heating and evaporation: comparative analysis, Fuel, 85(12-13), 1613-1630. Sazhin, S.S., Martynov, S.B., Kristyadi, T., Crua, C., Heikal, M.R. (2008) Diesel fuel spray penetration, heating, evaporation and ignition: modelling versus experimentation, International J of Engineering Systems Modelling and Simulation, 1(1) 1-19.

18. ,

19. , Plots of ethanol droplet temperature Td, measured experimentally (solid triangles) and predicted by the model (Tds droplet temperatures at the surface of the droplet, Tdav average droplet temperature, and Tdc droplet temperature at the centre of the droplet) and gas temperature Tg for the initial conditions Rdo= 118.65 mm, Tdo=294 K, C=3.97

20. , Plots of ethanol droplet temperature Td and radius Rd, measured experimentally (solid triangles and squares) and predicted by the model (Tds droplet temperatures at the surface of the droplet, Tdav average droplet temperature, and Tdc droplet temperature at the centre of the droplet) for gas average temperature Tg equal to 1270 K and for the initial conditions Rdo= 52.25 mm, Tdo=309 K, C=10.5

21. Main results(MODELLING VERSUS EXPERIMENTS) For relatively small droplets (initial radii about 65 μm) the experimentally measured droplet temperatures are close to the predicted average droplet temperatures. These temperatures are closer to the temperatures predicted at the centre of the droplets when the droplet diameter is larger than the probe volume size of the two-colour LIF thermometry.C. Maqua, G. Castanet, F. Grish, F. Lemoine, T. Kristyadi , S. S. Sazhin (2008)Monodisperse droplet heating and evaporation: experimental study and modelling, Int J Heat and Mass Transfer,51(15-16), 3932-3945, 2008

22. r T r T s s , Rd Rd , Kinetic Hydrodynamic region region j V q x d Rd 1 2 de

23. APPROXIMATION , .

24. Main results(Kinetic Modelling) • The kinetic effects predicted by the numerical algorithm turned out to be noticeable and cannot be a priori ignored when modelling droplet evaporation. • Shishkova, I.N. and Sazhin, S.S. (2006) A numerical algorithm for kinetic modelling of evaporation processes, J Computational Physics218 (2), 635-653. • Sazhin, S.S., Shishkova, I.N., Kryukov, A.P., Levashov, V.Yu. and Heikal, M.R. (2007) Evaporation of droplets into a background gas: kinetic modelling, Int J Heat Mass Transfer50, 2675-2691. • Sazhin, S.S., Shishkova, I.N. (2008) A kinetic algorithm for modelling the droplet evaporation process in the presence of heat flux and background gas, Atomization and Sprays (in press).

25. Some Related Developments • Development of the new model for transient stationary droplet heating [Int J Thermal Science(2007), V. 46 (4), pp. 444-457]. • Development of the dynamic decomposition method for numerical solution of the system of stiff ODEs [Computers and Fluids (2007), V. 36, pp. 601-610]. • Development of the new model of particle grouping in oscillating flows [European J of Mechanics B/Fluids (2008), V. 27, pp. 131-149; International J Heat and Fluid Flow (2008), V. 29, pp. 415-426]. • Development of droplet break-up models and their application to modelling transient Diesel fuel sprays [International J of Engineering Systems Modelling and Simulation (2008), V. 1(1), pp. 1-19 • Development of vortex ring models and their application to gasoline engines [submitted to J Fluid Mechanics]

26. Unsolved problems • Grid dependence of gas/droplets exchange of heat and mass under the Eulerian/Lagrangian approach in CFD codes. • Effects of turbulence on droplet heating. • Droplet heating during the break-up processes. • Heating of non-spherical droplets. • Heating of small droplets (geometrical optics approximation is not valid; surface tension is important). • Transient heating of moving droplets (boundary layer around the droplet is not fully developed) • Evaporation coefficient and inelastic collisions in kinetic models

27. Acknowledgements • The original results were obtained in collaboration with our colleagues • W. Abdelghaffar, S. Begg, V. Bykov, C. Crua, I. Goldfarb, • V. Gol’dshtein, D.Katoshevski, F. Kaplanskii, T. Kristyadi, A. Kryukov, P. Krutitski, E. Sazhina, T.Shakked, V. Sobolev. • The authors are grateful to the Royal Society • and EPSRC (Project EP/E02243X/1)for financial support.

28. Thank you for your attentionAny comments or suggestionswould be highly appreciated

29. ,

30. Modelling of droplet heating and evaporation in computational fluid dynamics codes Sergei SAZHIN*, Irina SHISHKOVA** , Vladimir LEVASHOV **, Morgan HEIKAL* *Sir Harry Ricardo Laboratory, School of Environment and Technology, Faculty of Science and Engineering, University of Brighton, Brighton BN2 4GJ, UK ** Low Temperature Departments, Moscow Power Engineering Institute, Moscow 111250, Russia

31. whenR Rd COUPLED TRANSIENT SOLUTION where , whenRd<R<Rg

32. when and when

33. h=χh0Fo= t κg /R R=Rg/Rd (no radiation)

34. Main results(Droplet Transient Heating) • The correction for convective heat transfer coefficient needs to be taken into account for all Fo. • The radiative effects of these corrections are negligibly small for small Fo, but it can become significant for large Fo. • Sazhin, S.S., Krutitskii, P.A., Martynov, S.B., Mason, D., Heikal, M.R., Sazhina, E.M. (2006) Transient heating of a semitransparent spherical body, Int J Thermal Science46 (4), 444-457.

35. when when

36. A countable set of positive eigenvalues λn is found from the solution of the equation These are arranged in ascending order 0<λ1 < λ2 < ..... , , .

37. Generalisation of the expression for h

38. C (ratio of the distance between droplets and their diameters), where , and C >3.

39. ,

40. ,

41. Collision of molecules

42. Boltzmann equations

43. (α=a,v; β=a,v)are collision integrals defined as

44. APPROXIMATION , .