170 likes | 179 Views
This study investigates the dynamic motion of a levitated droplet in an oscillating flow field using numerical simulations. The droplet's motion, flow field, and pressure distribution are analyzed and compared with experimental results.
E N D
Multiphysics 2009, Dec. 12, 2009 Numerical simulation of droplet motion and two-phase flow field in an oscillating container Tadashi Watanabe Center for Computational Science and e-Systems Japan Atomic Energy Agency o Background and Objectives o Numerical Method o Flow Field o Comparison o Summary
Background and Objectives Levitated Droplet : Free from effects of container wall Measurement of material properties of high-temperature molten metal,,, Levitation : electromagnetic, ultrasonic,,, Rotation : acoustic,,, Surface tension--- Oscillation frequency, Rotating shape,,, Viscosity --- Damping, Shape deformation,,, Oscillation Rotation Numerical simulations are performed to study the dynamic motion of the droplet in the oscillating flow fields.
Numerical Method (1) Oscillating Boundary Arbitrary Lagrangian-Eulerian mesh with oscillation speed of boundary Liquid Droplet Slip Boundary Incompressible + pseudo compressible Gas Oscillating Boundary
Numerical Method (2) Governing Equations for Fluid Motion Continuity Navier-Stokes Surface Tension Force Curvature Interpolation Pseudo Compressibility
Numerical Method (3) Governing Equations for Level Set Function interface Transport Reinitialization Mass Conservation
Numerical Method (4) FDM: 2nd order Adams-Bashforth method 2nd order upwind difference SMAC method for pressure and velocity Bi-CGSTAB method for Poisson equation Parallelization Simulation region : 10 mm x 17 mm (100x170) Droplet radius : 2 mm Time step : 1.0e-6 s Oscillation frequency : 20 kHz Sound pressure : 0.25~0.5 kPa Droplet : density = 998.2 kg/m3 viscosity = 0.998e-3 Ns/m2 surface tension = 0.0145 N/m Gas : density =1.166 kg/m3 viscosity = 1.819e-5 Ns/m2 sound speed = 340 m/s Liquid droplet Gas Oscillation
Numerical Method (5) Validation : Sloshing Experiment Liu and Lin, J. Comp. Phys. 227(2008)p3921 x=-5.0sinwt : w=6.0578s-1 190x100 probe3 probe2 probe1
Flow Field (1) Example of Pressure Distribution/Variation Vertical Position 0.0085 0.0 -0.0085 pressure node : 0.0 pressure node : -0.0085
Flow Field (2) Velocity Field and Droplet Motion Pressure node
Flow Field (3) Velocity Field and Droplet Motion t=0.00 s 0.05s 0.1 s 0.15 s 0.2 s 0.25 s Pressure node : middle
Flow Field (4) Velocity Field and Droplet Motion t=0.00 s 0.05s 0.1 s 0.15 s 0.2 s 0.25 s Pressure node : bottom
Flow Field (5) Velocity Field and Droplet Motion t=0.00 s 0.05s 0.1 s 0.15 s 0.2 s 0.25 s Pressure node : top
Flow Field (6) Vertical Position Droplet Position 0.0085 Pressure node Top 0.0 Middle Bottom -0.0085
Comparison (1) Incompressible Case Pressure node : bottom Pressure node : top t=0.00 s 0.05s 0.1 s 0.15 s 0.2 s 0.25 s
Comparison (2) Stationary/Oscillating Droplet Stationary Droplet (oscillating container) Oscillating Droplet (stationary container) scale x4 t=0.00 s 0.05s 0.1 s 0.15 s 0.2 s 0.25 s
Comparison (3) with Experiment Oscillating Circular Cylinder Tatsuno, Bull. Kyushu Univ. Appl. Mech., 128(2005)p23
Summary Motions of the droplet and the flow field in an oscillating container have been simulated numerically using the coupled level set and ALE method. ・Upward and downward flows from the droplet surface to the container wall appeared in the oscillating direction. ・ The droplet moved toward the pressure node, but this is not the case for incompressible case. ・ Induced flow field was similar to the flow field around an oscillating droplet/cylinder.