Computational Fluid Dynamics Study of the Effect of Gravity on the Impact of a Drop onto Dry and Wet Surfaces

Similar documents
Simulation of the Spray Characteristics of Black Liquor. Colloquium on BL Combustion and Gasification Denise Levesque

Splat formation in plasma-spray coating process*

Study of water falling film over horizontal drop shaped and inverted drop shaped tubes

CFD ANALYSIS OF CONVECTIVE FLOW IN A SOLAR DOMESTIC HOT WATER STORAGE TANK

Formation of Porosities during Spot Laser Welding: Case of Tantalum Joining

Computational Fluid Dynamics (CFD) Simulations of Dispersed Flow of Viscous Oil in Water in a Horizontal Pipe

PERFORMANCE ANALYSIS OF NATURAL DRAFT WET COOLING TOWER AT OPTIMIZED INJECTION HEIGHT

Deformation behavior of a liquid droplet impacting a solid surface

SIMULATION OF ZINC FILM FORMATION DURING CONTINUOUS WITHDRAWAL OF STEEL STRIPS FROM GALVANIZING BATHS

Effect of Rotation Speed of a Rotary Thermal Wheel on Ventilation Supply Rates of Wind Tower System

Numerical Simulation of the Collision of a Droplet with a Heated Solid Surface

COMSOL Multiphysics Simulation of Flow in a Radial Flow Fixed Bed Reactor (RFBR)

Numerical Simulation on Effects of Electromagnetic Force on the Centrifugal Casting Process of High Speed Steel Roll

Modelling of Material Removal in Abrasive Flow Machining Process Using CFD Simulation

CFD SIMULATIONS OF A CAPILLARY FORCE DRIVEN TWO-PHASE FLOW IN THE ANODE FLOW-FIELD OF PASSIVE-FEED µdmfc.

Multiphase Flow and Heat Transfer

Numerical Modeling of Buoyancy-driven Natural Ventilation in a Simple Three Storey Atrium Building

Numerical modelling of the filling of formworks with self-compacting concrete

Topography of laminar and oscillatory thermal and thermocapillary convection in a Czochralski process

VI.3 Application of Computational Fluid Dynamics to Simulate the Behaviour of Fluids inside Vertical Stack of Building Drainage System

Numerical and Experimental Modeling of Producer Gas Carburettor

Numerical modelling of the filling of formworks with self-compacting concrete

Liquid-Solid Phase Change Modeling Using Fluent. Anirudh and Joseph Lam 2002 Fluent Users Group Meeting

Task 1 For Task 1, the outlet was set as a zero-gauge pressure outlet, which is the same outlet condition as the five smaller pipes.

By Kimbal A. Hall, PE. Submitted to: WESTFALL MANUFACTURING COMPANY. May ALDEN RESEARCH LABORATORY, INC. 30 Shrewsbury Street Holden, MA 01520

ANSYS, Inc. March 31, Phase Change Models in FLUENT

EFFECT OF COMBINED HOLE CONFIGURATION ON FILM COOLING WITH AND WITHOUT MIST INJECTION

CFD Analysis of Pelton Runner

Modeling of micrometre-sized molten metallic droplet impact on a solid wall

Analysis of Side Sluice in Trapezoidal Channel

INVESTIGATIONS ON PERFORMANCE OF A SAVONIUS HYDROKINETIC TURBINE

Numerical Simulation of the Flue Gas Desulfurization in the Sieve Tray Tower

Numerical Investigation of Swirl's Effects in the Outer Annulus of a Reverse-flow Gas Turbine Combustor

MODELLING & SIMULATION OF LASER MATERIAL PROCESSING: PREDICTING MELT POOL GEOMETRY AND TEMPERATURE DISTRIBUTION

Optimization of Wind Loads Affecting Skyscrapers using ANSYS

DRUG DELIVERY PARTICLES FABRICATED USING THE ELECTROHYDRODYNAMIC ATOMIZATION (EHDA) METHOD

Splashing of a Drop Impacting on a Thin Liquid Film

NUMERICAL SIMULATION OF AIR NATURAL CIRCULATION AND THERMAL RADIATION IN PASSIVE CONTAINMENT COOLING SYSTEM

A Numerical Investigation on bubble formation in Coolant channel with various Ethylene Glycol and water compositions and change in contact angle

Bubble Penetration Through a Single Layer Sphere Bed

INVESTIGATION ON DRAG COEFFICIENT OF SUPERCIRTICAL WATER CROSS-FLOW PAST CYLINDER BIOMASS PARTICEL AT LOW REYNOLDS NUMBERS

On the Vortex Formation Effect During the Application of a Nitrogen-Gas Assisted Laser-Fusion Cutting Technique to Stainless Steel

Analysis on the influence of rotational speed to aerodynamic performance of vertical axis wind turbine

heat exchanger modelling for generator ventilation systems

INVESTIGATION OF FLOW CHARACTERISTICS OF HELLER-TYPE COOLING TOWERS WITH DIFFERENT COOLING DELTA ANGLES

Numerical investigations on the effect of slenderness ratio of matrix elements in cryogenic chill down process

Droplet formation mechanisms in metallurgical processes

Increase the Cooling Efficiency of Casting Furnace Stave Cooler with Use of Experimental and Analytical Analysis

CFD Modeling of Crankcase Inertial Oil Separators. Arun P. Janakiraman Anna Balazy Saru Dawar

SIMULATION OF SOLIDIFICATION PROCESS USED IN TWIN ROLL CASTING

COMPUTATIONAL STUDY OF TRANSITION OF OIL-WATER FLOW MORPHOLOGY DUE TO SUDDEN CONTRACTION IN MICROFLUIDIC CHANNEL

CFD and Wind Tunnel Study of the Performance of a Multi- Directional Wind Tower with Heat Transfer Devices

Computational Fluid Dynamic Analysis in De-staging of Centrifugal Pumps

DESIGN AND ANALYSIS OF GAS TURBINE INTERNAL COOLING PASSAGE P.Sethunathan, S.Prathap, M. Prabakaran, S. Pawanraj, R. Siddharth

Heat transfer enhancement in fire tube boiler using hellically ribbed tubes

CFD analysis of double-chambered crematories using biomass producer gas as a fuel source

International Conference on Materials Engineering and Information Technology Applications (MEITA 2015)

JJMIE Jordan Journal of Mechanical and Industrial Engineering

NUMERICAL STUDY ON FILM COOLING AND CONVECTIVE HEAT TRANSFER CHARACTERISTICS IN THE CUTBACK REGION OF TURBINE BLADE TRAILING EDGE

EVALUATE THE RISK OF HYDRATE BLOCKAGE INBEND FLOWLINES

Project Background (long term goals)

International Journal Of Scientific Research And Education Volume 3 Issue 3 Pages March-2015 ISSN (e): Website:

Numerical Analysis to Predict the Aerodynamic Performance of Tilt Wing of a Solar Powered UAV

ScienceDirect. CFD Analysis of a Kerosene Fuel Tank to Reduce Liquid Sloshing

Modelling Turbulent Non-Premixed Combustion in Industrial Furnaces

PLANNED COURSE ROTATION FOR CAE COURSES

Numerical Simulation of Airflow Organization and Temperature Distribution in Bulk Curing Barn Yong-qiao WANG, Yan-qiong LI, Yan ZHAO

ME 239: Rocket Propulsion. Real Nozzles. J. M. Meyers, PhD

Comparative Analysis of Different Orifice Geometries for Pressure Drop

Study on the Engineering Simulation of the Electric Hot Air De-icing. Device of the Substation Equipment. Jian Li1,a

INVESTIGATION OF INTERFACIAL INSTABILITIES WITH A TWO- LAYER SLIDE-COATING PROCESS

Numerical Investigation of a Single Phase Air Flow inside Porous Media of a Cross Flow Cooling Tower

ISSN: ISO 9001:2008 Certified International Journal of Engineering Science and Innovative Technology (IJESIT) Volume 3, Issue 2, March 2014

Simulation of Binary Mixture Droplet Evaporation using VOF Methods

FLUID STRUCTURE INTERACTION MODELLING OF WIND TURBINE BLADES BASED ON COMPUTATIONAL FLUID DYNAMICS AND FINITE ELEMENT METHOD

SOLIDIFICATION OF IMPINGING MOLTEN METAL DROPLET ON A COLD SUBSTRATE

Research Article Influence of the Slag Density on the Splashing Process in a Steelmaking Converter

Abstract. Introduction. PS Modeling non-uniform airflow distribution in large grain silos using Fluent. D. Garg 1, D.E.

A Study of Twin Co- and Counter-Rotating Vertical Axis Wind Turbines with Computational Fluid Dynamics

ESTIMATION OF RESIN FLOW FOR FRP BASE ON MPS METHOD

Lateral Outflow from Supercritical Channels

Prediction of Pollutant Emissions from Industrial Furnaces Using Large Eddy Simulation

Industrial Coatings: HIGHLY DURABLE, SUPERHYDROPHOBIC AND OLEOPHOBIC. (White Paper) NeverWet LLC

Numerical Analysis of Solar Updraft Tower with Solar Load applied to Discrete Ordinate Model

Evacuation Safety Evaluation of Inundated Stairs Using 3D Numerical Simulation

ANALYSIS OF PERFORMANCE OF NATURAL DRAFT COOLING TOWER AT OPTIMIZED INJECTION HEIGHT

CFD/FEM Based Analysis Framework for Wind Effects on Tall Buildings in Urban Areas

Texture and Wettability of Metallic Lotus Leaves

Analysis of Shear Wall Transfer Beam Structure LEI KA HOU

A New Analytical Approach to Improving the Aerodynamic Performance of the Gyromill Wind Turbine

Comparison between micrometer and millimetre sized metal and ceramic lamellas, for a better understanding of the spray process.

PS Modeling airflow in outdoor grain pile aeration systems. Introduction

Computational Analyses of Combustive Vortex Flows in Liquid Rocket Engines

Wind Loads Acting on Solar Panels in a Row by CFD Analysis

CFD Modelling of an Aerosol Exposure Chamber for Medical Studies G. Manenti, M. Derudi, G. Nano, R. Rota

Numerical method for modelling spray quenching of cylindrical forgings

Rapid DNA amplification in a capillary tube by natural convection with a single isothermal heater

Performance Analysis for Natural Draught Cooling Tower & Chimney through Numerical Simulation

Investigating Two Configurations of a Heat Exchanger in an Indirect Heating Integrated Collector Storage Solar Water Heating System

Modeling the after-treatment system of Diesel and S.I. internal combustion engines by means of OpenFOAM

Transcription:

Computational Fluid Dynamics Study of the Effect of Gravity on the Impact of a Drop onto Dry and Wet Surfaces Presenter : Murat Dinc MURAT DINC and DONALD D. GRAY, Dept. of Civil and Environmental Engineering, NICHOLAS HILLEN, J. STEPHEN TAYLOR, and JOHN KUHLMAN, Dept. of Mechanical and Aerospace Engineering, West Virginia University, Morgantown, WV 26506 3/20/2014 WVAS 2012 1

Background The impact of liquid drops onto both dry surfaces and wet surfaces has been studied by many researchers for more than a century The interaction between a drop and a surface may involve inertial, surface tension, viscous, and gravitational forces Experimental Sprays include millions of drop impingement on a surface Spray cooling is the key part in thermal management of the design of next generation electronic, aircraft and spacebased systems Computational 3/20/2014 2

Goal of the Spray Cooling Project The Project involves both experimental and computational studies aiming to develop correlations for the previously developed computationally effective Monte Carlo Spray Simulation Model In this presentation; contact angle, gravity effects will be investigated quantitatively for: A single drop impact on a dry surface A single drop impact on a wet surface 3/20/2014 WVAS 2012 3

Contact Angle/Wettability Lotus Leaf http://www.ramehart.com/contactangle.htm Surfaces with a contact angle < 90 are referred as hydrophilic, with an angle > 90 are referred as hydrophobic The higher the contact angle the higher the hydrophobicity of a surface, and results in less wettability Some plants show contact angles up to 160 and are called super-hydrophobic Plants like the lotus leaf can reach a contact angle of 170 Contact Angle Values Surface Glass Polymer Steel Teflon Liquid Water Water Water Water Contact angle, θ ( 0 ) 0 77 90 107 3/20/2014 WVAS 2012 4

Gravity Gravity (g) is a natural phenomenon by which physical bodies attract with a force proportional to their masses. Newton s Theory of Gravity: F g = m g Gravity Values Earth Gravity (m/s 2 ) Lunar Gravity (m/s 2 ) Zero Gravity Negative Gravity (m/s 2 ) 9.81 1.618 0-9.81 3/20/2014 WVAS 2012 5

Computational Approach Computational Model: CFD code: ANSYS 14 FLUENT The Navier-Stokes and continuity equations are solved in axisymmetric coordinates using the Finite Volume Method (FVM) for unsteady, incompressible and laminar flow Volume of Fluid (VOF) multiphase model is implemented with three phases (liquid drop, liquid layer and gas) to track the interface between the liquid phases and the surrounding gas 3/20/2014 WVAS 2012 6

Initial and Boundary Conditions Pressure Outlet 0.05 x 0.05 [m 2 ] axisymmetric domain Pressure Outlet Boundary Conditions Axis Pressure Outlet Wall Initial Conditions Smallest cell size= Do/168 Drop Liquid Layer 3/20/2014 WVAS 2012 7

Definition of Drop Impact Dry Surface Initial Stage Advancing Stage Receding Stage First Impact at t 0 Maximum Spread Drop receding D o V 0 V r V r Wet Surface D First Impact at t 0 Maximum crown height generation Crown collapse D o V 0 h H max V r V r 3/20/2014 WVAS 2012 8

Validation of the Computational Model t=1 ms t=3 ms t=7.5 ms t=10 ms Comparison of results for D = 4.2 mm droplet impacting on a thin liquid film with h = 2.1 mm: We = 2009, and h * = 0.5 (Experimental results of Wang and Chen (2000) in the left column of images, numerical results of Asadi and Passandideh-Fard (2009) in the middle column, numerical results of the current study with uniform-fine mesh elements in the right column) 3/20/2014 WVAS 2012 9

Non-Dimensional Parameters Dimensionless numbers in engineering and physics are used for data reduction of similar problems. Important parameters used in drops are: Reynolds Number: Froude Number: Re = U D o υ Fr = V2 gd Bond Number: Bo = gρd 2 σ Weber Number: Non-dimensional Film thickness: Non-Dimensional Numbers We = ρu2 D o σ h = h D o Surface Liquid We Bo Re Fr (Earth gravity) h/do Dry Water 28 0.55 2000 51 0 Dry Water 39 2.15 2800 36.4 0 Wet Water 137 2.7 6690 51.5 0.837 3/20/2014 WVAS 2012 10

Results/Comparisons Results: Gravity and Contact Angle effect on dry surfaces 3/20/2014 WVAS 2012 11

Normal Gravity (m/s2) Lunar Gravity (m/s2) Positive Earth Gravity 0 ms 5 ms Dry Surface-Gravity Effect Zero Gravity Negative Gravity (m/s2) 9.81 1.618 0-9.81 Negative Earth Gravity 0 ms 5 ms Contours of Volume Fraction, 90 Contact Angle (Blue shows liquid drop while red shows air) Liquid We Re Fr h/do Do (mm) Vo (m/s) Water 28 2000 51 0 2 1 10 ms 15 ms 10 ms 15 ms 20 ms 25 ms 20 ms 25 ms 30 ms 30 ms 3/20/2014 WVAS 2012 12

Dry Surface-Gravity Effect Contours of Volume Fraction (Blue shows liquid drop while red shows air) Zero Gravity Lunar Gravity 0 ms 5 ms 0 ms 5 ms 10 ms 15 ms 10 ms 15 ms 20 ms 25 ms 20 ms 25 ms 30 ms 30 ms 3/20/2014 WVAS 2012 13

Dry Surface-Contact Angle Effect Contours of Volume Fraction Earth Gravity 77 0 ms 5 ms Liquid We Re Fr h/do Do (mm) 90 107 0 ms 5 ms 0 ms 5 ms Vo (m/s) Water 39 2800 36.5 0 2.8 1 10 ms 15 ms 10 ms 15 ms 10 ms 15 ms 20 ms 25 ms 20 ms 25 ms 20 ms 25 ms 30 ms 35 ms 30 ms 35 ms 30 ms 35 ms 40 ms 40 ms 40 ms 3/20/2014 WVAS 2012 14

Results/Comparisons Results: Gravity and Contact Angle effect on wet surfaces 3/20/2014 WVAS 2012 15

Normal Earth Gravity (0 Contact Angle) Wet Surface-Contact Angle Effect Contours of Volume Fraction (Red shows liquid drop, green shows liquid layer while blue shows air) Liquid We Re Fr h/do Do (mm) Vo (m/s) Water 137 6690 51. 0.837 4.48 1.5 Normal Earth Gravity (180 Contact Angle) 0 ms 5 ms 0 ms 5 ms 8 ms 16 ms 8 ms 16 ms 25 ms 35 ms 25 ms 35 ms 45 ms 50 ms 45 ms 60 ms 3/20/2014 WVAS 2012 16

Wet Surface-Gravity Effect Contours of Volume Fraction (Red shows liquid drop, green shows liquid layer while blue shows air) Negative Earth Gravity 0 ms 5 ms Zero Gravity 0 ms 5 ms 8ms 16ms 8ms 16ms 25ms 35ms 25ms 35ms 45ms 60ms 45ms 60ms 82ms 80ms 3/20/2014 WVAS 2012 17

Conclusions-1 At advancing stage flow behavior is obtained very similar since Fr >> 1 (intertial forces are dominant over gravity) However, small differences start to appear at receding stage Generally, differences are not so much since Bond (Bo) number ~ 1 that makes gravitational forces comparable to surface tension forces Chandra and Qiao (1996) also showed same phenomenon where they experimentally tested a single water drop impact on a dry steel surface in normal and low gravity 3/20/2014 WVAS 2012 18

Conclusions-2 Results show that higher surface contact angles (hydrophobic), and negative gravity environment cannot be effective for spray cooling due to none or less drop liquid remaining on the surface However, analysis with heated surface is needed in order to study the heat transfer effectiveness of drops on dry and wet surfaces on low and normal earth gravity Computational approaches can be useful tool to study multiphase flow simulations such as drops and sprays, and can be correlated into the previously developed Monte Carlo Simulation Model 3/20/2014 WVAS 2012 19

Acknowledgments We gratefully acknowledge the financial support of this project under NASA Cooperative Agreement NNX10AN0YA. 3/20/2014 WVAS 2012 20

QUESTIONS/COMMENTS 3/20/2014 WVAS 2012 21

THANK YOU 3/20/2014 WVAS 2012 22