The effect of turbulence, gravity, and non-continuum hydrodynamic interactions on the drop size distribution in clouds

Basic Information

• Paper ID: 2501.01086
• Title: The effect of turbulence, gravity, and non-continuum hydrodynamic interactions on the drop size distribution in clouds
• Authors: Johnson Dhanasekaran, Donald L. Koch (Cornell University)
• Classification: physics.flu-dyn
• Publication Date: January 2, 2025
• Paper Link: https://arxiv.org/abs/2501.01086

Abstract

This study focuses on the evolution of micrometer-scale water droplets in clouds, particularly examining the "size gap" region at 15-40 μm radius, where condensation and differential sedimentation are least effective in promoting growth. This bottleneck leads to inaccurate growth models, while turbulence may correct discrepancies with in-situ cloud measurements. The research provides detailed analysis of the roles of turbulent collisions, droplet mixing, and water vapor fluctuations in bridging the "size gap." Collisions driven by the coupled effects of turbulent shear and differential sedimentation can grow drizzle-sized droplets. Turbulence-induced water vapor fluctuations also promote growth, maintaining polydispersity during the initial condensation-driven growth phase and facilitating subsequent differential sedimentation-driven coalescence growth.

Research Background and Motivation

Core Problem

The evolution of cloud droplet size distribution exhibits a critical "size gap" problem:

1. Small droplets (<15 μm): Grow primarily through condensation
2. Large droplets (>40 μm): Grow through differential sedimentation collision and coalescence
3. Medium droplets (15-40 μm): Neither mechanism is sufficiently effective, forming a growth bottleneck

Problem Significance

• Droplet size distribution affects radiative heat exchange, and its evolution controls precipitation
• Existing models predict rain formation times much longer than the tens of minutes observed in warm clouds
• This discrepancy has important implications for atmospheric research and climate modeling

Limitations of Existing Approaches

1. DNS Study Limitations: Direct numerical simulations by Li et al. did not include hydrodynamic interactions
2. Insufficient Hydrodynamic Treatment: Studies by Xue et al. and Grabowski et al. included hydrodynamic interactions but did not properly couple turbulence and gravity
3. Neglect of Non-continuum Effects: Most studies did not consider the breakdown of the continuum assumption as droplets approach
4. Reynolds Number Constraints: Taylor microscale Reynolds numbers in DNS studies are far below typical cloud values (O(10⁴))

Core Contributions

1. Comprehensive Physical Model: Developed a complete model simultaneously considering turbulent collisions, gravitational sedimentation, non-continuum hydrodynamic interactions, and inertial clustering
2. Monte Carlo Simulation Scheme: Established a Monte Carlo evolution scheme capturing turbulent intermittency and mixing
3. Quantification of Key Mechanisms: Quantitatively assessed the relative importance of various physical mechanisms in bridging the "size gap"
4. Parameter Sensitivity Analysis: Systematically studied the effects of mean free path, turbulence intensity, temperature, and pressure
5. Industrial and Planetary Applications: Provided theoretical foundation for industrial aggregators and extraterrestrial precipitation research

Methodology Details

Task Definition

Study the complete evolution of cloud water droplets from submicron condensation nuclei through condensation, turbulent collisions, and gravitational coalescence to drizzle size (>40 μm).

Model Architecture

1. Condensation Growth Model

Water droplet radius growth rate controlled by diffusion:

da/dt = D/(ρa)[ρ∞ - ρa(Ta)]

where differences between droplet temperature and surrounding air temperature are considered, determined through energy balance.

2. Turbulent Water Vapor Fluctuation Model

• Discretize air parcels into Np equal-volume parcels
• Water vapor content in each parcel evolves as an Ornstein-Uhlenbeck process
• Characteristic time scale is the Eulerian integral time scale τE
• Standard deviation is a fixed fraction f (typically 1%) of mean water vapor content

3. Collision Rate Model

Collision rate constant expressed as:

Gij = gij(r)|r=ai+aj × Kij

where:

• Kij is the collision kernel at zero droplet inertia
• gij(r) is the radial distribution function, capturing inertial clustering effects

4. Monte Carlo Evolution Scheme

• Use Monte Carlo droplets to represent multiple real droplets
• Collision probability: P = Gij·Nd,j·dt/V
• Modified probability formula accounting for computational constraints

Technical Innovations

1. Coupled Physical Processes: First complete coupling of turbulent shear, differential sedimentation, non-continuum hydrodynamics, and inertial clustering
2. Non-continuum Effects: Accurate treatment of the breakdown of continuum assumption as droplets approach
3. Turbulent Intermittency: Use log-normal distribution to simulate intermittency of Kolmogorov shear rate
4. Multi-parcel System: Capture different turbulent histories and water vapor fluctuations through parcel system

Experimental Setup

Default Parameters

• Turbulent dissipation rate: ε = 0.01 m²/s³
• Taylor microscale Reynolds number: Reλ = 2500
• Water vapor fluctuation intensity: f = 1%
• Surface temperature: T₀ = 293 K
• Surface pressure: P₀ = 1 bar
• Updraft velocity: 1 m/s

Simulation Parameters

• Volume: V = 1 m³
• Number of parcels: Np = 100
• Initial Monte Carlo droplets per parcel: 100
• Initial droplet radius: uniformly distributed 0.1-1 μm

Evaluation Metrics

1. Volume-averaged radius: av = Σa⁴ᵢ/Σa³ᵢ
2. Dispersion: D = σa/ (ratio of standard deviation to mean radius)

Experimental Results

Main Results

1. Complete Model Evolution

• Early phase: Condensation dominates, forming near-Gaussian distribution
• Middle phase: Entering "size gap," growth slows
• Late phase: Collision dominates, forming exponential tail distribution

2. Turbulence Intensity Effects

Under different turbulent dissipation rates (ε = 0.001-0.1 m²/s³):

• High turbulence significantly accelerates bridging the "size gap"
• Water vapor fluctuations more critical at low turbulence

3. Relative Importance of Mechanisms

By progressively disabling different mechanisms:

• Non-continuum Effects: Most critical; ignoring severely overestimates growth rate
• Inertial Clustering: Moderate enhancement (approximately 15%)
• Turbulent Intermittency: Weak influence
• Water Vapor Fluctuations: Critical at low turbulence; may slightly hinder growth at high turbulence

Ablation Studies

1. Water Vapor Fluctuation Effects

• f = 0.1% to 1%: Enhanced growth rate
• May slightly slow growth at high turbulence by increasing polydispersity

2. Mean Free Path Sensitivity

• Tested in range λg = 50-500 nm
• Larger mean free path significantly accelerates growth
• Earth conditions: λg ≈ 70-100 nm, collision efficiency β = 0.1-0.3

3. Initial Condition Effects

• Pressure: Primarily affects through mean free path
• Temperature: Primarily affects through water vapor carrying capacity

Important Findings

1. "Lucky Droplet" Hypothesis Invalid: Few large droplets produced by turbulent intermittency insufficient to affect overall distribution evolution
2. Limited Inertial Clustering: While enhancing, far from order-of-magnitude improvement
3. Non-continuum Effects Dominant: Collision efficiency can be as low as 0.01-0.3 under typical cloud conditions
4. Dual Role of Polydispersity: Promotes differential sedimentation but hinders turbulent collisions

Related Work

Main Research Directions

1. DNS Studies: Direct numerical simulations by Li et al., Ayala et al.
2. Collision Kernel Theory: Saffman-Turner turbulent collision theory
3. Inertial Clustering: Caustic effect research by Wilkinson et al.
4. Non-continuum Fluid Mechanics: Lubrication theory by Sundararajakumar and Koch

Advantages of This Work

1. First comprehensive consideration of all key physical processes
2. Accurate treatment of high Reynolds number turbulent conditions
3. Complete non-continuum hydrodynamics included
4. Systematic parameter sensitivity analysis

Conclusions and Discussion

Main Conclusions

1. Turbulent collisions sufficient to bridge the "size gap" without requiring water vapor fluctuations
2. Non-continuum hydrodynamics is the critical limiting factor
3. Water vapor fluctuations important under low turbulence conditions
4. Inertial clustering provides moderate enhancement
5. Turbulent intermittency has limited influence

Limitations

1. Does not consider "slingshot effect" in strong inertia regime
2. Assumes viscous forces dominate hydrodynamic interactions
3. Neglects droplet deformation and van der Waals forces
4. Simplified turbulent statistics model

Future Directions

1. Extension to industrial aggregator applications
2. Extraterrestrial planetary precipitation research
3. More accurate turbulence-gravity coupling models
4. Experimental validation and comparison with field observations

In-Depth Evaluation

Strengths

1. Physical Completeness: First comprehensive consideration of all key processes in cloud microphysics
2. Advanced Numerical Methods: Monte Carlo scheme effectively handles discrete collisions and continuous condensation
3. Comprehensive Parameter Study: Systematic exploration of broad parameter space
4. High Practical Value: Provides practical tools for cloud physics and industrial applications
5. Significant Theoretical Contribution: Clarifies relative importance of multiple disputed mechanisms

Weaknesses

1. Computational Complexity: Model is relatively complex with potentially high computational cost
2. Lack of Experimental Validation: Primarily based on theoretical and numerical studies, lacking direct experimental comparison
3. Simplified Parameterization: Some turbulent statistics still use simplified models
4. Limited Applicability: Primarily focused on warm clouds, does not consider ice phase processes

Impact

1. Academic Contribution: Provides new theoretical framework for cloud physics
2. Practical Value: Can improve weather forecasting and climate models
3. Interdisciplinary Applications: Provides tools for industrial processes and planetary science
4. Methodological Significance: Monte Carlo method can be generalized to other multiphase flow systems

Applicable Scenarios

1. Cloud Physics Research: Modeling water droplet growth in warm clouds
2. Weather Forecasting: Improving precipitation forecasting microphysics schemes
3. Industrial Applications: Optimizing aerosol aggregator design
4. Planetary Science: Precipitation processes in extraterrestrial atmospheres

References

The paper cites 28 important references covering multiple fields including cloud physics, turbulence theory, and numerical methods, providing a solid theoretical foundation for this research.