What Is CLSVOF
What Is CLSVOF
What Is CLSVOF
Yan Zhan
8-18-2011
Outline
• What is CLSVOF?
• Why CLSVOF?
• How to implement CLSVOF?
What is CLSVOF?
Volume of Fluid (VOF) Level Set (LS)
• Volumetric phase fraction F • Level-Set function Ф
Phase 1 F=1 Phase 1 Ф>0
Phase 2 F=0 Phase 2 Ф<0
Interface 0<F<1 Interface Ф=0
• Transport of F: • Transport of F:
~ ~
( ρF ) + ∇ ⋅ ( ρUF ) = 0
t
( ρφ ) + ∇ ⋅ ( ρUφ ) = 0
t
CLSVOF
Why CLSVOF?
Development of the liquid jet (time step is 2.5 μm) (Menard, 2007)
Why CLSVOF?
Liquid parcels
How to implement CLSVOF?
• Couple LS with VOF within the CFD code FLUENT by
implementing user defined functions (UDF)
• UDF
– User written program that can be linked with FLUENT at
run-time
– Programmed in C and FLUENT defined macros
– User-defined scalar (UDS) transport modeling customize
FLUENT for level set equation
How to implement CLSVOF?
How to implement CLSVOF?
(t = 0.2s)
(t = 0.2s)
How to implement CLSVOF?
• Setup UDS for LS in FLUENT
Scalar ø Transport Equation
∂ρφ ~
+ ∇ ⋅ ( ρUφ ) = 0
∂t
∂ρφ
– Unsteady term
∂t
– Convection term ~
∇ ⋅ ( ρUφ )
– Diffusive term
0
– Source term
0
∂φ
Additional term appear for turbulent flow such as − ρ u ′φ ′ = Γ
∂x
j t
j
How to implement CLSVOF?
• Setup UDS for LS in FLUENT
– Set number of UDS
– Set UDS terms (Appendix A)
• DEFINE_UDS_UNSTEADY
– Get unsteady term for scalar equation
• DEFINE_UDS_FLUX
– Returns user specified flux
• DEFINE_DIFFUSIVITY
– Returns user diffusion coefficient (Γ)
• DEFINE_SOURCE
– Set UDS boundary conditions
• Constant
• UDF: DEFINE_PROFILE
Appendix Equations
• Incompressible two-phase flow
∇ ⋅U = 0
∇p 1 1
U + U ⋅ ∇U = − + ∇ ⋅ (2µ (φ ) D ) − γκ (φ )∇H (φ ) + F
t
ρ (φ ) ρ (φ ) ρ (φ )
φ + U ⋅ ∇φ = 0
t
F + ∇ ⋅ (UF ) = 0
t
Density ρ(Ф), viscosity μ (Ф), and curvature κ (Ф) are written as,
ρ (φ ) = ρ (1 − H (φ )) + ρ H (φ )
g l
µ (φ ) = µ (1 − H (φ )) + µ H (φ )
g l
∇φ
κ (φ ) = ∇ ⋅
∇φ
D is defined as the rate of deformation tensor
D = (∇U ) + (∇U ) T
Appendix Equations
• Incompressible two-phase flow
The surface tension force is
1
γκ (φ )∇H (φ )
ρ (φ )
where H is the Heaviside function,
1 if φ > 0
H (φ ) =
0 otherwise.
F will be initialized in each computational cell Ωij
1
F = ∫ H (φ ( r , z ,0))rdrdz
∆r∆z
ij Ω ij
where Ωij is
Ω = (r , z ) r ≤ r ≤ r
ij i i +1
and z ≤ z ≤ z
j j +1
Appendix Equations
• Re-Initialization
– Reinitialize ф
where w is the characteristic velocity pointing outward from the free surface