Dynamics of Droplet-Droplet and Droplet-Film Collision. C. K. Law Princeton University

Similar documents
Report Documentation Page

Attribution Concepts for Sub-meter Resolution Ground Physics Models

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)

Thermo-Kinetic Model of Burning for Polymeric Materials

System Reliability Simulation and Optimization by Component Reliability Allocation

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.

A report (dated September 20, 2011) on. scientific research carried out under Grant: FA

Report Documentation Page

Diagonal Representation of Certain Matrices

Broadband matched-field source localization in the East China Sea*

SW06 Shallow Water Acoustics Experiment Data Analysis

Report Documentation Page

Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures

Crowd Behavior Modeling in COMBAT XXI

On Applying Point-Interval Logic to Criminal Forensics

CRS Report for Congress

Marginal Sea - Open Ocean Exchange

Improvements in Modeling Radiant Emission from the Interaction Between Spacecraft Emanations and the Residual Atmosphere in LEO

Predictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER

Direct Numerical Simulation of Aeolian Tones

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects

Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST)

Molecular Characterization and Proposed Taxonomic Placement of the Biosimulant 'BG'

Aviation Weather Routing Tool: A Decision Aid for Manned/Unmanned Aircraft Routing

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION

Scattering of Internal Gravity Waves at Finite Topography

Experimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin

Periodic Magnetoresistance Oscillations in Side-Gated Quantum Dots

DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited.

STUDY OF DETECTION LIMITS AND QUANTITATION ACCURACY USING 300 MHZ NMR

A GENERAL PROCEDURE TO SET UP THE DYADIC GREEN S FUNCTION OF MULTILAYER CONFORMAL STRUCTURES AND ITS APPLICATION TO MICROSTRIP ANTENNAS

Wavelets and Affine Distributions A Time-Frequency Perspective

PIPS 3.0. Pamela G. Posey NRL Code 7322 Stennis Space Center, MS Phone: Fax:

Real-Time Environmental Information Network and Analysis System (REINAS)

REPORT DOCUMENTATION PAGE

Predicting Tropical Cyclone Formation and Structure Change

REPORT DOCUMENTATION PAGE. Theoretical Study on Nano-Catalyst Burn Rate. Yoshiyuki Kawazoe (Tohoku Univ) N/A AOARD UNIT APO AP

Volume 6 Water Surface Profiles

Catalytic Oxidation of CW Agents Using H O in 2 2 Ionic Liquids

Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts

INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS

Ocean Acoustics Turbulence Study

Extension of the BLT Equation to Incorporate Electromagnetic Field Propagation

Comparative Analysis of Flood Routing Methods

An Examination of 3D Environmental Variability on Broadband Acoustic Propagation Near the Mid-Atlantic Bight

Determining the Stratification of Exchange Flows in Sea Straits

Experimental and Theoretical Studies of Ice-Albedo Feedback Processes in the Arctic Basin

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Modulation Instability of Spatially-Incoherent Light Beams and Pattern Formation in Incoherent Wave Systems

Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector

VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ

HIGH-POWER SOLID-STATE LASER: LETHALITY TESTING AND MODELING

LAGRANGIAN MEASUREMENTS OF EDDY CHARACTERISTICS IN THE CALIFORNIA CURRENT

High Resolution Surface Characterization from Marine Radar Measurements

Effects of Constrictions on Blast-Hazard Areas

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data.

Playing Abstract games with Hidden States (Spatial and Non-Spatial).

Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method

Coupled Ocean-Atmosphere Modeling of the Coastal Zone

DSMC Calculation of Vortex Shedding behind a Flat Plate with a New Intermolecular Collision Scheme

Generation and Propagation of Internal Solitary Waves on the Continental Shelf and Slope

FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES

The Bernoulli equation in a moving reference frame

An Observational and Modeling Study of Air-Sea Fluxes at Very High Wind Speeds

Regional Stratification and Shear of the Various Streams Feeding the Philippine Straits ESR Component

Optimizing Robotic Team Performance with Probabilistic Model Checking

Coastal Engineering Technical Note

MTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS

Characterization of Caribbean Meso-Scale Eddies

SMA Bending. Cellular Shape Memory Structures: Experiments & Modeling N. Triantafyllidis (UM), J. Shaw (UM), D. Grummon (MSU)

HYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling. Naval Research Laboratory, Stennis Space Center, USA

Mass Transport by Second Mode Internal Solitary Waves

Coastal Mixing and Optics

USMC Enlisted Endstrength Model

Rogue Wave Statistics and Dynamics Using Large-Scale Direct Simulations

Nonlinear Internal Waves: Test of the Inverted Echo Sounder

Modeling the Impact of Extreme Events on Margin Sedimentation

EFFECTS OF LOCAL METEOROLOGICAL VARIABILITY ON SURFACE AND SUBSURFACE SEISMIC-ACOUSTIC SIGNALS

Satellite Observations of Surface Fronts, Currents and Winds in the Northeast South China Sea

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Improved Parameterizations Of Nonlinear Four Wave Interactions For Application In Operational Wave Prediction Models

Computer Simulation of Sand Ripple Growth and Migration.

Design for Lifecycle Cost using Time-Dependent Reliability

Reverse Ion Acceleration by Laser-Matter Interaction

Understanding Near-Surface and In-cloud Turbulent Fluxes in the Coastal Stratocumulus-topped Boundary Layers

Surface Fluxes and Wind-Wave Interactions in Weak Wind Conditions

The Extratropical Transition of Tropical Cyclones

INTERACTION AND REMOTE SENSING OF SURFACE WAVES AND TURBULENCE

Internal Waves and Mixing in the Aegean Sea

Maximizing the Bandwidth from Supercontinuum Generation in Photonic Crystal Chalcogenide Fibers

Vortex Rossby Waves and Hurricane Evolution in the Presence of Convection and Potential Vorticity and Hurricane Motion

uniform distribution theory

Abyssal Current Steering of Upper Ocean Current Pathways in an Ocean Model with High Vertical Resolution

Contract No. N C0123

Transcription:

Dynamics of Droplet-Droplet and Droplet-Film Collision C. K. Law Princeton University The physical phenomena of droplet-droplet and droplet-film collision in the head-on orientation were studied experimentally and computationally, with emphasis on the transition between bouncing and merging of the liquid surfaces. Experimentally, the droplets (~300 µm diameter) were generated using the ink-jet printing technique, and imaged using stroboscopy and high-speed cine-photography for the droplet-droplet and droplet-film collision events respectively. Computationally, the collision event was simulated using the front-tracking technique. For the study of droplet-droplet collision, the instant of merging was experimentally determined and then used as an input in the computational simulation of the entire collision event. The simulation identified the differences between collision and merging at small and large Weber numbers, and satisfactorily described the dynamics of the interdroplet gap including the role of the van der Waals force in effecting surface rupture. For the study of droplet-film collision, extensive experimental mapping showed that the collision dynamics is primarily affected by the droplet Weber number (We) and the film thickness scaled by the droplet radius (H), that while droplet absorption by the film is facilitated with increasing droplet Weber number, the boundary of transition is punctuated by an absorption peninsula, in the We-H space, within which absorption is facilitated for smaller Weber numbers. Results from computation simulation revealed the essential dependence of the collision dynamics on the restraining nature of the solid surface, the energy exchange between the droplet and the film, and the coherent motion of the gas-liquid interfaces. Partial absorption with the emission of a secondary droplet of smaller size was also observed and explained. 1

Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 1. REPORT DATE 22 JUN 2004 2. REPORT TYPE N/A 3. DATES COVERED - 4. TITLE AND SUBTITLE Dynamics of Droplet-Droplet and Droplet-Film Collision 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) Princeton University 8. PERFORMING ORGANIZATION REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR S ACRONYM(S) 12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release, distribution unlimited 11. SPONSOR/MONITOR S REPORT NUMBER(S) 13. SUPPLEMENTARY NOTES See also ADM001793, International Symposium on Energy Conversion Fundamentals Held in Istanbul, Turkey on 21-25 June 2005., The original document contains color images. 14. ABSTRACT 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT UU a. REPORT unclassified b. ABSTRACT unclassified c. THIS PAGE unclassified 18. NUMBER OF PAGES 18 19a. NAME OF RESPONSIBLE PERSON Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18

Dynamics of Droplet-Droplet and Click to edit Master title style Second level Droplet-Film Collision C. K. Law, Princeton University Fourth International level Symposium on Fifth Energy levelconversion Fundamentals June 21-25, 2004, Istanbul, Turkey 1 111 1

Practical Relevance Click to edit Master title style Droplet-droplet collision: - Spray Click processes to edit Master text styles - Rain drop formation - Model Second for nuclear level fusion Fourth chamber level - Raindrop Fifth level interaction with ocean and soil Droplet-surface collision: - Spray impingement on surface of combustion - Meteorite impaction - Ink-jet printing - Spray-painting/coating 2 222 2

Fundamental Interests Click to edit Master title style 3 Dynamics of deformable bodies Large deformation implies nonlinear dynamics Many decades of length scales: Second level - Droplet size: 10 ~ 10 3 µm - Molecular Third level force distance: 10-2 µm - Continuum flows Wide range of physical regimes: - Rarefied flows - van der Waals force Computationally & experimentally challenging Description of state of merging 333 3

Click to edit Master title style Second level 4 444 4

Computational Simulation Click to edit Master title style Front tracking method of Tryggvason & immersed Click to edit boundary Master method text styles of Peskin Second level Continuum Third levelformulation ( ρfifth V) ( level 1 T 8 + ρvv) = p + g + ( V + V ) t Re We s κnδ ( r r f ) da 5 555 5

Dynamics of Colliding Interfaces Click to edit Master title style 0.5 0-0.5 0.02 Click to edit Master t <= 0.45 text styles 0.76 t > 0.45 Second level 0.92 0.015 1.07 h 1.11 0.01 0.005 t = 0.13 0.45 0.27 0 0.002 0.004 0.006 0.008 0.01 0 0 0.2 0.4 0.6 0.8 1 1.2 t h center h rim 6 666 6

Matched Instant of Merging Click to edit Master title style Identify experimental instant of curvature reversal of interface 0.298 (0.124 ms) cusp Second as state level of surface Third level rupture 0.622 (0.258 ms) simulation Fourth level 0.916 (0.380 ms) Use this instant in Close agreement with experimental images t = 0 0.953 (0.395 ms) 0.984 (0.408 ms) 1.076 (0.446 ms) 1.362 (0.565 ms) 1.637 (0.679 ms) 1.944 (0.806 ms) 2.390 (0.990 ms) 2.923 (1.212 ms) 7 1.013 (0.420 ms) 3.280 (1.360 ms) 777 7

Mis-Matched Instants of Merging Click to edit Master title style Second level t = 0 1.076 0.298 0.622 0.916 0.953 1.362 1.637 1.944 2.390 0.984 2.923 8 1.013 3.280 Advanced rupture Delayed rupture 888 8

Surface Rupture Induced by Click Molecular to edit Master Force title (1/3) style To effect surface rupture, add van der 0.01 Waals force in formulation: Second level A/(6πδ 3 0.008 ) h 0.006 divergence of 0.004 interface Fourth for level certain A 0.002 A* = 1.6 10-5 Simulation shows Adjust A to agree with empirical instant of rupture 0.012 h rim (A* 0) h center (A* = 0) h rim (A* = 0) A* = 2.4 10 - Rupture 0 5 0.4 0.6 0.8 1 1.2 1.4 1.6 t A* = 0.8 10-5 9 999 9

Surface Rupture Induced by Click Molecular to edit Master Force title (2/3) style 0.012 0.012 h 0.01 0.008 0.006 0.004 0.002 0.01 Click to hedit rim (A* 0) Master text styles 0.008 Second h center (A* = 0) level 0.006 h rim (A* = 0) 0.004 A* = 0.8 10-5 A* = 1.6 10-5 0.002 A* = 2.4 10-5 Rupture 0 0 0.4 0.6 0.8 1 1.2 1.4 1.6 Soft merging occurs late in collision, at rim t h center (A* = 0) h rim (A* 0) Rupture A* = 0.8 10-6 A* = 2.4 10-6 A* = 2.0 10-6 0.2 0.4 0.6 0.8 1 1.2 Hard merging occurs early in collision, at rim t h rim (A* = 0) 10 10 10

Surface Rupture Induced by Click Molecular to edit Master Force title (3/3) style Augmented Hamaker constant is 10 3 to 10Click 6 larger to edit than Master the real text one styles Implying Second the level interfacial gap width from simulation Third level is one to two orders too large Consistent Fourth level with: Simulated Fifth levelgap width ~ 10-1 µm van der Waals force range ~ 10-3 -10-2 µm 11 11 11

Missing Physics Click to edit Master title style Additional physics needed so that the interface Click to can edit reach Master the text range styles of the molecular Second level force Need Third rarefied level flow treatment Need Fourth variable level density treatment Computationally Fifth level challenging Possibly described through analysis 12 12 12

Dynamics of Droplet-Film Collision Click to edit (Thin Master film) title style Second level (a) (a) Bouncing (b) Partial absorption (c) Total absorption (b) 13 (c) 13 13

Dynamics of Droplet-Film Collision Click to edit (Thick Master film) title style 14 Second level (a) (b) (c) (a) Bouncing (b) Partial absorption (c) Total absorption 14 14

Regime Diagram of Click Droplet-Film to edit Master Collision title style Controlling parameters: Weber 2 number & film thickness Second (H=h/R) level 1.5 Bouncing 1 peninsula around H 1 0.5 Unquantified boundary Existence of merging Existence of triple reversal in absorption & bouncing with increasing H H = h /R 2.5 0 0 5 10 15 20 25 30 We = 2R ρ V 2 /σ Merging R = 0.165 mm R = 0.250 mm R = 0.305 mm 15 15 15

Factors Affecting Click Bouncing to edit Master vs Absorption title style Solid surface restrains droplet motion Impact inertia narrows interfacial gap Droplet Second flattening level increases resistance Droplet deformation increases viscous loss Fourth level Energy Fifth level transfer to film Oscillatory resonance around H 1 16 16 16

Summary Click to edit Master title style Instant of droplet merging can be empirically determined and used to simulate complete droplet collision sequence Second level augmented Third level van der Waals force Surface rupture can be induced through an Proper description of interfacial dynamics, allowing for rarefied flow and variable density, is essential For droplet-film collision, existence of merging peninsula around H 1 17 17 17