Kinetic Modeling and Monte-Carlo Simulations of Droplet Coalescence in a Turbulent Gas Flow PUBLIC ACCESS

[+] Author Affiliations
Philippe Villedieu

ONERA, Toulouse Cedex, France

Olivier Simonin

IMFT, Toulouse, France

Paper No. FEDSM2002-31385, pp. 673-674; 2 pages
  • ASME 2002 Joint U.S.-European Fluids Engineering Division Conference
  • Volume 1: Fora, Parts A and B
  • Montreal, Quebec, Canada, July 14–18, 2002
  • Conference Sponsors: Fluids Engineering Division
  • ISBN: 0-7918-3615-0 | eISBN: 0-7918-3600-2
  • Copyright © 2002 by ASME


Two-phase gas-droplet flows are involved in a lot of industrial applications, especially in the combustion field (Diesel engine, turbomachinery, rocket engine,[[ellipsis]]). Among all the characteristics of the spray, the droplet size distribution generally has a major influence on the global performances of the system and must be accurately taken into account in a numerical simulation code. This is a difficult task because the carrier gas flow is very often turbulent. Hence, droplets located in the vicinity of the same point may have different velocities and coalesce, leading at the end to a strong modification of the initial droplet size distribution. The first part of our contribution will be devoted to the presentation of a new kinetic model for droplet coalescence in turbulent gas flows. This model is an extension, to the case of sprays, of the ideas introduced by Simonin, Deutsch and Laviéville in [1]. The key ingredient is the use of the “joint density function”, fgp (t, x , r, v , u ), representing the density of droplets at time t, located at point x , with radius r and velocity v and “viewing” an instantaneous turbulent gas velocity u . The great advantage of using fgp (t, x , r, v , u ) instead of the usual density function fp (t, x , r, v ) is the possibility to close the collision operator, in the governing kinetic equation, with less restrictive assumptions on the velocity correlations of two colliding droplets. The link between this model and the usual one (relying on the so-called “chaos assumption”) will be discussed. In the second part of our contribution, we shall present a new Monte-Carlo algorithm derived from our kinetic model. Numerical simulation results, for some academic test cases (homogeneous isotropic turbulence), will be shown and compared to the results obtained with a classical algorithm for droplet collision, based on the chaos assumption (see for example [2] or [3]). The figure 1 below shows a comparison between the temporal evolution of the mass mean radius computed by a classical collision model (neglecting the influence of gas and droplet velocity correlation) and by the “joint-pdf” based model. In the first case, the growth rate of the droplet, due to coalescence phenomena, is overestimated. Moreover, figure 2 shows that the droplet kinetic energy, induced by the turbulent gas motion, decays rapidly with the chaos assumption based model, as already noticed by Laviéville et al [1] in the case of solid particle collisions.

Copyright © 2002 by ASME
This article is only available in the PDF format.



Interactive Graphics


Country-Specific Mortality and Growth Failure in Infancy and Yound Children and Association With Material Stature

Use interactive graphics and maps to view and sort country-specific infant and early dhildhood mortality and growth failure data and their association with maternal

Citing articles are presented as examples only. In non-demo SCM6 implementation, integration with CrossRef’s "Cited By" API will populate this tab (http://www.crossref.org/citedby.html).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In