High-Resolution Simulation of Turbulent Collision of Cloud Droplets

A novel parallel implementation of hybrid DNS (Direct Numerical Simulation) code for simulating collision-coalescence of aerodynamically interacting particles in a turbulent flow has been developed. An important application of this code is to quantify turbulent collision-coalescence rate of cloud droplets, relevant to warm rain formation, under physically realistic conditions. The code enables performing high-resolution DNS of turbulent collisions so the simulation results can be used to begin addressing the question of Reynolds number dependence of pair and collision statistics. The new implementation is based on MPI (Message Passing Interface) library, and thus the code can run on computers with distributed memory. This development enables to conduct hybrid DNS with flow field solved at grid resolutions up to 5123 while simultaneously track up to several million aerodynamically-interacting droplets. In this paper we discuss key elements of the MPI implementation and present preliminary results from the high resolution simulations. The key conclusion is that, for small cloud droplets, the results on pair statistics and collision kernel appear to reach their saturation values as the flow Reynolds number is increased.

[1]  Wojciech W. Grabowski,et al.  Theoretical Formulation of Collision Rate and Collision Efficiency of Hydrodynamically Interacting Cloud Droplets in Turbulent Atmosphere , 2005 .

[2]  Diego Donzis,et al.  Dissipation and enstrophy in isotropic turbulence: Resolution effects and scaling in direct numerical simulations , 2008 .

[3]  L. Biferale,et al.  Acceleration statistics of heavy particles in turbulence , 2006, Journal of Fluid Mechanics.

[4]  Bogdan Rosa,et al.  Effects of turbulence on the geometric collision rate of sedimenting droplets. Part 2. Theory and parameterization , 2008 .

[5]  Y. Kaneda,et al.  Study of High-Reynolds Number Isotropic Turbulence by Direct Numerical Simulation , 2009 .

[6]  Lian-Ping Wang,et al.  Parallel Implementation of Particle Tracking and Collision in a Turbulent Flow , 2009, PPAM.

[7]  Lance R. Collins,et al.  Reynolds number scaling of particle clustering in turbulent aerosols , 2004 .

[8]  Lian-Ping Wang,et al.  Scalable parallel FFT for spectral simulations on a Beowulf cluster , 2001, Parallel Comput..

[9]  A. Wexler,et al.  Statistical mechanical description and modelling of turbulent collision of inertial particles , 1998, Journal of Fluid Mechanics.

[10]  Paul A. Vaillancourt,et al.  Statistics and Parameterizations of the Effect of Turbulence on the Geometric Collision Kernel of Cloud Droplets , 2007 .

[11]  Wojciech W. Grabowski,et al.  Effects of stochastic coalescence and air turbulence on the size distribution of cloud droplets , 2006 .

[12]  Y. Kaneda,et al.  Small-scale statistics in high-resolution direct numerical simulation of turbulence: Reynolds number dependence of one-point velocity gradient statistics , 2007, Journal of Fluid Mechanics.

[13]  Bogdan Rosa,et al.  A spurious evolution of turbulence originated from round-off error in pseudo-spectral simulation , 2009 .

[14]  Lian-Ping Wang,et al.  A hybrid approach for simulating turbulent collisions of hydrodynamically-interacting particles , 2007, J. Comput. Phys..

[15]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[16]  M. Maxey,et al.  Settling velocity and concentration distribution of heavy particles in homogeneous isotropic turbulence , 1993, Journal of Fluid Mechanics.

[17]  Lian-Ping Wang On the dispersion of heavy particles by turbulent motion , 1993 .

[18]  Bogdan Rosa,et al.  Effects of turbulence on the geometric collision rate of sedimenting droplets. Part 1. Results from direct numerical simulation , 2008 .

[19]  K. Sreenivasan On the universality of the Kolmogorov constant , 1995 .

[20]  Bogdan Rosa,et al.  Turbulent collision of inertial particles: point-particle based, hybrid simulations and beyond , 2009 .