Title: Resolving hydrodynamic singularities through molecular interactions
1Resolving hydrodynamic singularities through
molecular interactions
- Len Pismen
- Technion, Haifa, Israel
-
- Supported by Israel Science Foundation
- Minerva Center for Nonlinear Physics of
Complex Systems
2Cusp on a free interface
3Hydrodynamic problem
kinematic condition
boundary conditions
?
tangential stress normal stress
?????? ?
complex flow potential
Stokes solution
boundary conditions
4Conformal transformation
conformal transformation
???
????? ???
? ?????
Ca???at ? ????? ?
? ?????
5Hydrodynamic solution
must vanish at q-p/2
defines parameter ? as a function of Cau0????
Ca
?
6Hydrodynamic solution flow pattern
Ca 0.1?
Ca ????
7Flow near the tip inflow or stagnation point?
Ca
log10(??L)
8Anatomy of a cusp
Hydrodynamic picture the interface of zero
thickness terminates in a cusp singularity Shikhmu
rzaevs picture there is an interfacial layer
(1) of unspecified nature its properties change
starting from equilibrium values the cusp widens
into a transition zone (2), and a
surface-tension-relaxation tail (3) becomes a
gradually disappearing internal
interface. Diffuse interface picture A narrow
interfacial zone (1) is close to local
equilibrium vapor condenses near the cusp (2)
bulk fluid relaxes to equilibrium in the
diffusion zone (3)
9Equilibrium density functional theory
Free energy of inhomogeneous fluid
e.g. hard-core potential
Euler Lagrange equation
chemical potential
static equation of state
10Dynamic diffuse interface theory
Cahn Hilliard equation (in Galilean frame)
dynamic equation of state
G mobility
coupling to hydrodynamics
Stokes equation
continuity equation
compressible flow in boundary region
11Multiscale perturbation scheme
- Find an equilibrium density profile
- liquid at z ???? ????? at z ???
- Perturbations
- Weak gradient of chemical potential
- Propagation of interphase boundary
- Disjoining/conjoining potential
- Interfacial curvature
- Expand in scale ratio
- Use solvability condition
- Match to outer solution
12Equilibrium density profile
1D static density functional equation
Q(z)
Density profile for a van der Waals fluid
h3 tail
h3 tail
Gibbs surface (defines the nominal thickness)
13Perturbation scheme inner solution
- Inner equation for chemical potential
Inner chemical potential
- Material balance across the layer
- Dynamic shift of chemical potential
14Interaction of interfaces
liquid
Change of surface tension
h
vapor
liquid
Shift of chemical potential
equilibrium
valid at hgtgtd significant at hd
shifted equilibrium
151d solutions of the nonlocal equations
h/d
g
h/d
16Variable surface tension
surface tension vs. gap width
g
h lt d
h gt d
h/d
Ca
Cusp appears at a finite Ca it is lower than
Ca in the hydrodynamic solution with the same
tip curvature
Ca(?)
Ca(d)
log10(?L or L/d)
17Effect of variable surface tension velocity
?
?d
velocity
d 108
velocity
Ca??? ???????
?
blow-up near the tip
q
d 108
velocity
velocity
d 108
d 106
Ca??? ???????
q
18Larger Ca? Higher-order singularity
1/3 gt a gt 8/19 0 gt b gt 1/95 close
together
tail close together
a1/3 b0
a8/19 b1/95
tip diverge
19Outer solution gas phase
longitudinal vapor flux
liquid
condensation zone
longitudinal air flux
gas flow
longitudinal condensation flux
liquid
local condensation flux
condensation starts when the gap narrows to
provide necessary driving force
gap h at the entrance of condensation zone
107 m
20Outer solution liquid phase
density depletion
convective diffusion equation ( u velocity
along the cusp)
equilibrium
shifted equilibrium
potential difference driving the condensation flux
Weak but slowly decaying density depletion
downstream
h/d
21Conclusions
- Cusp singularity is resolved trough molecular
interactions on nanoscale - Effect of variable surface tension
- Cusp is formed at finite capillary number
- Weak Marangoni flow near the cusp
- Effect of condensation
- Air is removed from the cusp by diffusion
- Weak depletion of liquid density near the cusp
- Relaxation effect
- Weak depletion tail downstream from the cusp