Absorption of gas by a falling liquid film

Size: px
Start display at page:

Download "Absorption of gas by a falling liquid film"

Transcription

1 Absorption of gas by a falling liquid film Christoph Albert Dieter Bothe Mathematical Modeling and Analysis Center of Smart Interfaces/ IRTG 1529 Darmstadt University of Technology 4th Japanese-German International Workshop on Mathematical Fluid Dynamics, Tokyo Mathematical Modeling and Analysis Christoph Albert

2 Wavy Falling Films Advantages Good heat transfer due to small thickness Large Interface Applications Evaporation Cooling Absorption Tokyo Mathematical Modeling and Analysis Christoph Albert 2

3 Hydrodynamic Model Assumptions Incompressible, Newtonian two-phase flow No phase transition Constant surface tension t (ρu)+ (ρu u)+ p= S+ρ g, Ω Σ u=0, Ω Σ p I S n Σ =σ κ n Σ, Σ u =0, Σ S=η( u+( u) T ) Tokyo Mathematical Modeling and Analysis Christoph Albert 3

4 The Volume of Fluid Method One Fluid Formulation: t (ρu)+ (ρu u)+ p= S+ρ g+δ f Σ u=0 t f +u f = Density and viscosity depend on volume fraction Volume fraction has to be transported Additional surface tension force term f Σ Tokyo Mathematical Modeling and Analysis Christoph Albert 4

5 Parasitic Currents VOF simulations suffer from unphysical oscillation of velocity, so-called parasitic currents Stem from numerical treatment of interfacial jump condition for stress: p I S n Σ =σ κ n Σ Especially serious in stagnant flow situations Parasitic Currents in a VOF simulation of a static droplet Also problematic in simulations of falling films, which are convection dominated Tokyo Mathematical Modeling and Analysis Christoph Albert 5

6 Numerical Setup u(t, y) x=0 =(1+ϵsin(2 π ω t )) ( ρ l μ L g δ 0 2 [ y δ ( y δ 0 ) 2 ], 0 ) Tokyo Mathematical Modeling and Analysis Christoph Albert 6

7 Continuum Surface Stress (CSS) ( f σ(i n Σ n Σ )) Body Force, Approximate n Σ by differentiating a smoothed f-field Momentum conservative Standard Surface Tension model in FS3D; delivers good results in many twophase flow situations Tokyo Mathematical Modeling and Analysis Christoph Albert 7

8 Falling Films with CSS Water/Air, Re cells per mean film thickness (0.265mm) A = 30%, f = 20 Hz Tokyo Mathematical Modeling and Analysis Christoph Albert 8

9 Continuum Surface Force (CSF) Body Force σ κ f In the case of a sphere and constant curvature, an exact balance between surface tension and pressure can be achieved: Balanced Force Renardy, J. Comput. Phys. 183 (2002) h i Greatly reduces parasitic currents Relies on good curvature information Easy for a falling film: κ= h xx 3 2 (1+h x 2 ) i 1 i i Tokyo Mathematical Modeling and Analysis Christoph Albert 9

10 Falling Films with CSF Water/Air, Re cells per mean film thickness (0.265mm) A = 30%, f = 20 Hz Tokyo Mathematical Modeling and Analysis Christoph Albert 10

11 Comparison to experiment: Dietze Experimental data from DMSO/Air = m 2 Re = 8.6 f = 16 Hz s G.F. Dietze, Flow Separation in Falling Liquid Films, 2010, = kg m 3, = N m A = 40% Resolution Tokyo Mathematical Modeling and Analysis Christoph Albert 11

12 Comparison to experiment: Film thickness Film thickness at x = 56 mm over a time span of 0.3 s CSS overshoots Tokyo Mathematical Modeling and Analysis Christoph Albert 12

13 Species Transport in Falling Films Data from Yoshimura et.al., 1996 Transport of Oxygen into a water film at 18 C Sc= ν L D =570 Sh= k Lδ 0 D k L = Γ L ln ( C S C in C S C out ) Tokyo Mathematical Modeling and Analysis Christoph Albert 13

14 Species Transport in Falling Films: Model Assumptions: Dilute system => species bears no mass or momentum, and does not affect viscosity No adsorption at the interface => surface tension stays constant No chemical reaction Local thermodynamic equilibrium at the interface => Henry's law holds Constant Henry coefficient t c+u c=d Δ c, D c n Σ =0, c L =H c G, Ω G (t ) Ω L (t) Σ Σ Tokyo Mathematical Modeling and Analysis Christoph Albert 14

15 Numerical Approach Two scalar approach with one-sided concentration gradient transfer flux Concentration is advected in the same way as volume fraction Two concentrations stored in each interfacial cell Both concentrations and Henry's law yield interfacial concentration Concentration gradient in liquid phase computed between interface and some value in bulk by subgrid model Diffuse flux from gas to liquid according to Fick's Law ϕ G =c χ ΩG ϕ L =c χ ΩL Tokyo Mathematical Modeling and Analysis Christoph Albert 15

16 Calculations Water film at 18 C in Oxygen atmosphere δ 0 Resolution 16 Re=31 Sc=50 k H,cc = c G =4.19e-5 c L,0 =2.37e-7 D=2.14e Tokyo Mathematical Modeling and Analysis Christoph Albert 16

17 Validation: Method Solve stationary advection diffusion problem on domain with boundary conditions. u( y) x c=d 2 y c [0, x max ] [0,δ 0 ] c x=0 =c L,0 y c y=0 =0 x c x=10cm =0 Solved with Matlab ODE Solver, by defining x as pseudotime. y c y=δ0 =k H,cc c G Tokyo Mathematical Modeling and Analysis Christoph Albert 17

18 Validation: Result Tokyo Mathematical Modeling and Analysis Christoph Albert 18

19 0 Hz 0.46mm 0mm 5.90cm 8.53cm Planar Interface Pure Diffusion Tokyo Mathematical Modeling and Analysis Christoph Albert 19

20 20 Hz 0.53mm mm 5.25cm 9.19cm Time-periodic wave structures appear Large Wave humps, preceeded by several smaller capillary waves Tokyo Mathematical Modeling and Analysis Christoph Albert 20

21 20 Hz 0.53mm mm 6.33cm 6.70cm Filaments of high concentration in the large wave humps Develop along the streamlines of the large vortex Touch the interface at a hyperbolic point Tokyo Mathematical Modeling and Analysis Christoph Albert 21

22 20Hz Highly nonmonotonous concentration profiles Tokyo Mathematical Modeling and Analysis Christoph Albert 22

23 20 Hz Strong contribution from the capillary wave region Tokyo Mathematical Modeling and Analysis Christoph Albert 23

24 20 Hz 0.53mm 0mm 6.59cm 6.77cm Flow up the wall in the reference frame of the wall Tokyo Mathematical Modeling and Analysis Christoph Albert 24

25 20 Hz 0.33mm 0mm 6.55cm Pressure in film according to Young- Laplace: Δ p σ κ 6.79cm Increase in pressure large enough to drive water up the wall Compare Dietze et al., J. Fluid Mech., 637, Tokyo Mathematical Modeling and Analysis Christoph Albert 25

26 30 Hz 0.53mm 0mm 7.22cm 9.84cm Less capillary waves Wave length, peak height, and wave velocity decrease Tokyo Mathematical Modeling and Analysis Christoph Albert 26

27 40 Hz 0.53mm 0mm 7.39cm 9.65cm No capillary waves; waves become sinusoidal Wave length, peak height, and wave velocity decrease further Tokyo Mathematical Modeling and Analysis Christoph Albert 27

28 40 Hz 0.39mm 0mm 7.91cm 8.03cm At this Reynolds Number, there exists backflow even when Capillary Waves are absent Tokyo Mathematical Modeling and Analysis Christoph Albert 28

29 Re 15 / 15Hz 0.43mm 0mm 5.84cm 7.87cm At this Reynolds Number, concentration profiles are monotonous Tokyo Mathematical Modeling and Analysis Christoph Albert 29

30 Re 15 / 15Hz Tokyo Mathematical Modeling and Analysis Christoph Albert 30

31 Re 15 / 15Hz 0.43mm mm 6.98cm 7.12cm No Vortex relative to wave velocity Tokyo Mathematical Modeling and Analysis Christoph Albert 31

32 Thank you for your attention! Tokyo Mathematical Modeling and Analysis Christoph Albert 32

Reduction of parasitic currents in the DNS VOF code FS3D

Reduction of parasitic currents in the DNS VOF code FS3D M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring

More information

Mass Transfer Operations I Prof. Bishnupada Mandal Department of Chemical Engineering Indian Institute of Technology, Guwahati

Mass Transfer Operations I Prof. Bishnupada Mandal Department of Chemical Engineering Indian Institute of Technology, Guwahati Mass Transfer Operations I Prof. Bishnupada Mandal Department of Chemical Engineering Indian Institute of Technology, Guwahati Module - 2 Mass Transfer Coefficient Lecture - 4 Boundary Layer Theory and

More information

Introduction to Heat and Mass Transfer. Week 10

Introduction to Heat and Mass Transfer. Week 10 Introduction to Heat and Mass Transfer Week 10 Concentration Boundary Layer No concentration jump condition requires species adjacent to surface to have same concentration as at the surface Owing to concentration

More information

Convective Mass Transfer

Convective Mass Transfer Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface

More information

CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM EVAPORATOR

CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM EVAPORATOR Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM

More information

Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport

Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport Christoph Lehrenfeld and Arnold Reusken Preprint No. 413 December 2014 Key words: Two-phase flow, mass transport,

More information

Contents. I Introduction 1. Preface. xiii

Contents. I Introduction 1. Preface. xiii Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................

More information

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum) 2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities

More information

Convection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.

Convection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds. Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,

More information

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit

More information

Simulation of Liquid Jet Breakup Process by Three-Dimensional Incompressible SPH Method

Simulation of Liquid Jet Breakup Process by Three-Dimensional Incompressible SPH Method Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 212 ICCFD7-291 Simulation of Liquid Jet Breakup Process by Three-Dimensional Incompressible SPH

More information

Direct Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid

Direct Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid Direct Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid Nima. Samkhaniani, Azar. Ajami, Mohammad Hassan. Kayhani, Ali. Sarreshteh Dari Abstract In this paper, direct numerical simulation

More information

Outline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer

Outline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer Diffusion 051333 Unit operation in gro-industry III Department of Biotechnology, Faculty of gro-industry Kasetsart University Lecturer: Kittipong Rattanaporn 1 Outline Definition and mechanism Theory of

More information

Mass Transfer Fundamentals. Chapter#3

Mass Transfer Fundamentals. Chapter#3 Mass Transfer Fundamentals Chapter#3 Mass Transfer Co-efficient Types of Mass Transfer Co-efficient Convective mass transfer can occur in a gas or liquid medium. Different types of mass transfer coefficients

More information

A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows

A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows M. Raessi, M. Bussmann*, and J. Mostaghimi Department of Mechanical and Industrial Engineering,

More information

NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE

NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00786 NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A

More information

Basic concepts in viscous flow

Basic concepts in viscous flow Élisabeth Guazzelli and Jeffrey F. Morris with illustrations by Sylvie Pic Adapted from Chapter 1 of Cambridge Texts in Applied Mathematics 1 The fluid dynamic equations Navier-Stokes equations Dimensionless

More information

Height function interface reconstruction algorithm for the simulation of boiling flows

Height function interface reconstruction algorithm for the simulation of boiling flows Computational Methods in Multiphase Flow VI 69 Height function interface reconstruction algorithm for the simulation of boiling flows M. Magnini & B. Pulvirenti Dipartimento di Ingegneria Energetica, Nucleare

More information

n v molecules will pass per unit time through the area from left to

n v molecules will pass per unit time through the area from left to 3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular

More information

Derivation of continuum models for the moving contact line problem based on thermodynamic principles. Abstract

Derivation of continuum models for the moving contact line problem based on thermodynamic principles. Abstract Derivation of continuum models for the moving contact line problem based on thermodynamic principles Weiqing Ren Courant Institute of Mathematical Sciences, New York University, New York, NY 002, USA Weinan

More information

Introduction to Heat and Mass Transfer. Week 12

Introduction to Heat and Mass Transfer. Week 12 Introduction to Heat and Mass Transfer Week 12 Next Topic Convective Heat Transfer» Heat and Mass Transfer Analogy» Evaporative Cooling» Types of Flows Heat and Mass Transfer Analogy Equations governing

More information

Part I.

Part I. Part I bblee@unimp . Introduction to Mass Transfer and Diffusion 2. Molecular Diffusion in Gasses 3. Molecular Diffusion in Liquids Part I 4. Molecular Diffusion in Biological Solutions and Gels 5. Molecular

More information

Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,

Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain

More information

A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations

A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm

More information

Numerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand

Numerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand Institute of Aerospace Thermodynamics Numerical modelling of phase change processes in clouds Challenges and Approaches Martin Reitzle Bernard Weigand Introduction Institute of Aerospace Thermodynamics

More information

Fluid Mechanics Theory I

Fluid Mechanics Theory I Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to

More information

Spri ringer. INTERFACIAL TRANSPORT PHENOMENA 2 nd Edition. John C. Slattery Department ofaerospace Engineering Texas A&M University

Spri ringer. INTERFACIAL TRANSPORT PHENOMENA 2 nd Edition. John C. Slattery Department ofaerospace Engineering Texas A&M University INTERFACIAL TRANSPORT PHENOMENA 2 nd Edition John C. Slattery Department ofaerospace Engineering Texas A&M University Leonard Sagis Department of Agrotechnology & Food Science Wageningen University Eun-Suok

More information

The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011

The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 Derivation of Navier-Stokes Equation 1 The total stress tensor

More information

Simulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework

Simulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,

More information

A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant

A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant Ashley J. James Department of Aerospace Engineering and Mechanics, University of Minnesota John Lowengrub

More information

CHEN 7100 FA16 Final Exam

CHEN 7100 FA16 Final Exam CHEN 7100 FA16 Final Exam Show all work and state all assumptions for full credit. The exam is closed book, notes, and homework. Only the course reader and your mind should be open. No electronic devices

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

More information

AC & DC Magnetic Levitation and Semi-Levitation Modelling

AC & DC Magnetic Levitation and Semi-Levitation Modelling International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, March 24-26, 2003 AC & DC Magnetic Levitation and Semi-Levitation Modelling V. Bojarevics, K. Pericleous Abstract

More information

Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer

Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer 1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic

More information

Numerical simulation of polyurethane foaming processes on bubble scale

Numerical simulation of polyurethane foaming processes on bubble scale Numerical simulation of polyurethane foaming processes on bubble scale 7th OpenFOAM Workshop Darmstadt, 25-28 June 2012 Stephanie Geier and Manfred Piesche Institute of Mechanical Process Engineering University

More information

Mechanical Engineering. Postal Correspondence Course HEAT TRANSFER. GATE, IES & PSUs

Mechanical Engineering. Postal Correspondence Course HEAT TRANSFER. GATE, IES & PSUs Heat Transfer-ME GATE, IES, PSU 1 SAMPLE STUDY MATERIAL Mechanical Engineering ME Postal Correspondence Course HEAT TRANSFER GATE, IES & PSUs Heat Transfer-ME GATE, IES, PSU 2 C O N T E N T 1. INTRODUCTION

More information

meters, we can re-arrange this expression to give

meters, we can re-arrange this expression to give Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.

More information

Chapter 5. The Differential Forms of the Fundamental Laws

Chapter 5. The Differential Forms of the Fundamental Laws Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations

More information

Free energy concept Free energy approach LBM implementation Parameters

Free energy concept Free energy approach LBM implementation Parameters BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH

More information

Fluid Dynamics Exercises and questions for the course

Fluid Dynamics Exercises and questions for the course Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r

More information

Mush liquid interfaces with cross flow

Mush liquid interfaces with cross flow Mush liquid interfaces with cross flow Devin Conroy March 15, 27 1 Introduction The solidification of a binary melt growing into a supercooled region may lead to the formation of a mushy layer as a result

More information

6.2 Governing Equations for Natural Convection

6.2 Governing Equations for Natural Convection 6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed

More information

Transport processes. 7. Semester Chemical Engineering Civil Engineering

Transport processes. 7. Semester Chemical Engineering Civil Engineering Transport processes 7. Semester Chemical Engineering Civil Engineering 1 Course plan 1. Elementary Fluid Dynamics 2. Fluid Kinematics 3. Finite Control Volume nalysis 4. Differential nalysis of Fluid Flow

More information

Principles of Convective Heat Transfer

Principles of Convective Heat Transfer Massoud Kaviany Principles of Convective Heat Transfer Second Edition With 378 Figures Springer Contents Series Preface Preface to the Second Edition Preface to the First Edition Acknowledgments vii ix

More information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed

More information

Differential equations of mass transfer

Differential equations of mass transfer Differential equations of mass transfer Definition: The differential equations of mass transfer are general equations describing mass transfer in all directions and at all conditions. How is the differential

More information

Boundary Conditions in Fluid Mechanics

Boundary Conditions in Fluid Mechanics Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial

More information

MOLECULAR SIMULATION OF THE MICROREGION

MOLECULAR SIMULATION OF THE MICROREGION GASMEMS2010-HT01 MOLECULAR SIMULATION OF THE MICROREGION E.A.T. van den Akker 1, A.J.H. Frijns 1, P.A.J. Hilbers 1, P. Stephan 2 and A.A. van Steenhoven 1 1 Eindhoven University of Technology, Eindhoven,

More information

2. FLUID-FLOW EQUATIONS SPRING 2019

2. FLUID-FLOW EQUATIONS SPRING 2019 2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid

More information

Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition

Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow

More information

Silica deposition in superheated geothermal systems

Silica deposition in superheated geothermal systems PROCEEDINGS, 43rd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 1-14, 018 SGP-TR-13 Silica deposition in superheated geothermal systems Vijay Chauhan,

More information

Chapter 2 Mass Transfer Coefficient

Chapter 2 Mass Transfer Coefficient Chapter 2 Mass Transfer Coefficient 2.1 Introduction The analysis reported in the previous chapter allows to describe the concentration profile and the mass fluxes of components in a mixture by solving

More information

2.3 The Turbulent Flat Plate Boundary Layer

2.3 The Turbulent Flat Plate Boundary Layer Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

6. Laminar and turbulent boundary layers

6. Laminar and turbulent boundary layers 6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary

More information

Math 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations

Math 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The

More information

Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk

Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,

More information

x j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2

x j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2 Merging of drops to form bamboo waves Yuriko Y. Renardy and Jie Li Department of Mathematics and ICAM Virginia Polytechnic Institute and State University Blacksburg, VA -, U.S.A. May, Abstract Topological

More information

Chapter 6 Fundamental Concepts of Convection

Chapter 6 Fundamental Concepts of Convection Chapter 6 Fundamental Concepts of Convection 6.1 The Convection Boundary Layers Velocity boundary layer: τ surface shear stress: s = μ u local friction coeff.: C f y y=0 τ s ρu / (6.) (6.1) Thermal boundary

More information

Heat and Mass Transfer Unit-1 Conduction

Heat and Mass Transfer Unit-1 Conduction 1. State Fourier s Law of conduction. Heat and Mass Transfer Unit-1 Conduction Part-A The rate of heat conduction is proportional to the area measured normal to the direction of heat flow and to the temperature

More information

Unit operations of chemical engineering

Unit operations of chemical engineering 1 Unit operations of chemical engineering Fourth year Chemical Engineering Department College of Engineering AL-Qadesyia University Lecturer: 2 3 Syllabus 1) Boundary layer theory 2) Transfer of heat,

More information

PHYS 432 Physics of Fluids: Instabilities

PHYS 432 Physics of Fluids: Instabilities PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible

More information

Sound Pressure Generated by a Bubble

Sound Pressure Generated by a Bubble Sound Pressure Generated by a Bubble Adrian Secord Dept. of Computer Science University of British Columbia ajsecord@cs.ubc.ca October 22, 2001 This report summarises the analytical expression for the

More information

Introduction to Heat and Mass Transfer. Week 9

Introduction to Heat and Mass Transfer. Week 9 Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional

More information

Measurement of Liquid Film Thickness in Micro Square Channel

Measurement of Liquid Film Thickness in Micro Square Channel Measurement of Liquid Film Thickness in Micro Square Channel Youngbae Han and Naoki Shikazono Department of Mechanical Engineering, The University of Tokyo Hongo 7-3-1, Bunkyo-ku, Tokyo 113-8656, Japan

More information

Thermocapillary Migration of a Drop

Thermocapillary Migration of a Drop Thermocapillary Migration of a Drop An Exact Solution with Newtonian Interfacial Rheology and Stretching/Shrinkage of Interfacial Area Elements for Small Marangoni Numbers R. BALASUBRAMANIAM a AND R. SHANKAR

More information

Development and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations

Development and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Development and analysis of a Lagrange-Remap sharp interface solver

More information

Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow

Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow OCEN 678-600 Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow Date distributed : 9.18.2005 Date due : 9.29.2005 at 5:00 pm Return your solution either in class or in my mail

More information

1 Introduction to Governing Equations 2 1a Methodology... 2

1 Introduction to Governing Equations 2 1a Methodology... 2 Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................

More information

CFD in COMSOL Multiphysics

CFD in COMSOL Multiphysics CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective

More information

Boundary Conditions for the Moving Contact Line Problem. Abstract

Boundary Conditions for the Moving Contact Line Problem. Abstract Boundary Conditions for the Moving Contact Line Problem Weiqing Ren Courant Institute of Mathematical Sciences, New York University, New York, NY 10012, USA Weinan E Department of Mathematics and PACM,

More information

Solving the Euler Equations!

Solving the Euler Equations! http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p

More information

Numerical Studies of Droplet Deformation and Break-up

Numerical Studies of Droplet Deformation and Break-up ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Break-up B. T. Helenbrook Department of Mechanical and

More information

UNIVERSITY of LIMERICK

UNIVERSITY of LIMERICK UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF

More information

IHTC DRAFT MEASUREMENT OF LIQUID FILM THICKNESS IN MICRO TUBE ANNULAR FLOW

IHTC DRAFT MEASUREMENT OF LIQUID FILM THICKNESS IN MICRO TUBE ANNULAR FLOW DRAFT Proceedings of the 14 th International Heat Transfer Conference IHTC14 August 8-13, 2010, Washington D.C., USA IHTC14-23176 MEASUREMENT OF LIQUID FILM THICKNESS IN MICRO TUBE ANNULAR FLOW Hiroshi

More information

arxiv: v1 [physics.flu-dyn] 1 Mar 2017

arxiv: v1 [physics.flu-dyn] 1 Mar 2017 Direct numerical simulation of variable surface tension flows using a Volume-of-Fluid method arxiv:1703.00327v1 [physics.flu-dyn] 1 Mar 2017 Ivana Seric, Shahriar Afkhami, Lou Kondic Department of Mathematical

More information

LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE

LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE Proceedings of the ASME/JSME 2011 8th Thermal Engineering Joint Conference AJTEC2011 March 13-17, 2011, Honolulu, Hawaii, USA AJTEC2011-44190 LIQUID FILM THICKNESS OF OSCILLATING FLOW IN A MICRO TUBE Youngbae

More information

P = 1 3 (σ xx + σ yy + σ zz ) = F A. It is created by the bombardment of the surface by molecules of fluid.

P = 1 3 (σ xx + σ yy + σ zz ) = F A. It is created by the bombardment of the surface by molecules of fluid. CEE 3310 Thermodynamic Properties, Aug. 27, 2010 11 1.4 Review A fluid is a substance that can not support a shear stress. Liquids differ from gasses in that liquids that do not completely fill a container

More information

Candidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.

Candidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request. UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book

More information

Mass Transfer Operations

Mass Transfer Operations College of Engineering Tutorial # 1 Chemical Engineering Dept. 14/9/1428 1. Methane and helium gas mixture is contained in a tube at 101.32 k Pa pressure and 298 K. At one point the partial pressure methane

More information

FLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes

FLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam

More information

Diffusion and Adsorption in porous media. Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad

Diffusion and Adsorption in porous media. Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad Diffusion and Adsorption in porous media Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad Contents Introduction Devices used to Measure Diffusion in Porous Solids Modes of transport in

More information

Level Set Tumor Growth Model

Level Set Tumor Growth Model Level Set Tumor Growth Model Andrew Nordquist and Rakesh Ranjan, PhD University of Texas, San Antonio July 29, 2013 Andrew Nordquist and Rakesh Ranjan, PhD (University Level Set of Texas, TumorSan Growth

More information

CFD Simulation Of Hydrodynamics And Heat Transfer In Flow Of Liquids On Inclined Plates

CFD Simulation Of Hydrodynamics And Heat Transfer In Flow Of Liquids On Inclined Plates CFD Simulation Of Hydrodynamics And Heat Transfer In Flow Of Liquids On Inclined Plates Veeranna Modi 1*, Omprakash Hebbal 2 1* PG Student, Thermal Power Engineering, PDA College of Engineering, Gulbarga-585102,

More information

LINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION

LINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION LINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION É. FONTANA 1, E. MANCUSI 1,2, A. A. ULSON DE SOUZA 1, S. M. A. GUELLI U. SOUZA 1 1 Universidade Federal

More information

5.2 Surface Tension Capillary Pressure: The Young-Laplace Equation. Figure 5.1 Origin of surface tension at liquid-vapor interface.

5.2 Surface Tension Capillary Pressure: The Young-Laplace Equation. Figure 5.1 Origin of surface tension at liquid-vapor interface. 5.2.1 Capillary Pressure: The Young-Laplace Equation Vapor Fo Fs Fs Fi Figure 5.1 Origin of surface tension at liquid-vapor interface. Liquid 1 5.2.1 Capillary Pressure: The Young-Laplace Equation Figure

More information

Numerical Simulation of Film Flow over an Inclined Plate: Effects of Solvent Properties and Contact Angle

Numerical Simulation of Film Flow over an Inclined Plate: Effects of Solvent Properties and Contact Angle Numerical Simulation of Film Flow over an Inclined Plate: Effects of Solvent Properties and Contact Angle Janine Carney and Rajesh Singh Multiphase Flow Science Workshop August 5-6, 214 Lakeview Golf Resort

More information

Surface and Interfacial Tensions. Lecture 1

Surface and Interfacial Tensions. Lecture 1 Surface and Interfacial Tensions Lecture 1 Surface tension is a pull Surfaces and Interfaces 1 Thermodynamics for Interfacial Systems Work must be done to increase surface area just as work must be done

More information

FORMULA SHEET. General formulas:

FORMULA SHEET. General formulas: FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to

More information

Understanding Magnetic Field Gradient Effect From a Liquid Metal Droplet Movement

Understanding Magnetic Field Gradient Effect From a Liquid Metal Droplet Movement Donghong Gao* Neil B. Morley Vijay Dhir Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA 90095 Understanding Magnetic Field Gradient Effect From a Liquid Metal

More information

BAE 820 Physical Principles of Environmental Systems

BAE 820 Physical Principles of Environmental Systems BAE 820 Physical Principles of Environmental Systems Inter phase mass transfer Dr. Zifei Liu Mass transfer between two phases For a two phase system not at equilibrium (concentrations in the two phases

More information

Computational Fluid Dynamics 2

Computational Fluid Dynamics 2 Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2

More information

Quick Recapitulation of Fluid Mechanics

Quick Recapitulation of Fluid Mechanics Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change

More information

CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer

CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic

More information

PORE-SCALE SIMULATION OF MASS TRANSFER ACROSS IMMISCIBLE INTERFACES

PORE-SCALE SIMULATION OF MASS TRANSFER ACROSS IMMISCIBLE INTERFACES PORE-SCALE SIMULATION OF MASS TRANSFER ACROSS IMMISCIBLE INTERFACES A report submitted to the department of Energy Resources Engineering of Stanford University in partial fulfillment of the requirements

More information

Differential relations for fluid flow

Differential relations for fluid flow Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow

More information

TABLE OF CONTENTS CHAPTER TITLE PAGE

TABLE OF CONTENTS CHAPTER TITLE PAGE v TABLE OF CONTENTS CHAPTER TITLE PAGE TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES LIST OF SYMBOLS LIST OF APPENDICES v viii ix xii xiv CHAPTER 1 INTRODUCTION 1.1 Introduction 1 1.2 Literature Review

More information

Finite Element Model of a complex Glass Forming Process as a Tool for Control Optimization

Finite Element Model of a complex Glass Forming Process as a Tool for Control Optimization Excerpt from the Proceedings of the COMSOL Conference 29 Milan Finite Element Model of a complex Glass Forming Process as a Tool for Control Optimization Felix Sawo and Thomas Bernard Fraunhofer Institute

More information

ACOUSTICS SIMULATION IN THE PRESENCE OF MOVING INTERFACES IN MULTIPHASE FLOWS

ACOUSTICS SIMULATION IN THE PRESENCE OF MOVING INTERFACES IN MULTIPHASE FLOWS 6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 18, Glasgow, UK ACOUSTICS SIMULATION IN THE PRESENCE OF MOVING INTERFACES

More information

I. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010

I. Borsi. EMS SCHOOL ON INDUSTRIAL MATHEMATICS Bedlewo, October 11 18, 2010 : an : an (Joint work with A. Fasano) Dipartimento di Matematica U. Dini, Università di Firenze (Italy) borsi@math.unifi.it http://web.math.unifi.it/users/borsi porous EMS SCHOOL ON INDUSTRIAL MATHEMATICS

More information