Module 3: "Thin Film Hydrodynamics" Lecture 11: "" The Lecture Contains: Micro and Nano Scale Hydrodynamics with and without Free Surfaces

Similar documents
REE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics

Differential relations for fluid flow

Module 3: "Thin Film Hydrodynamics" Lecture 12: "" The Lecture Contains: Linear Stability Analysis. Some well known instabilities. Objectives_template

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

Lecture 2: Hydrodynamics at milli micrometer scale

contact line dynamics

Chapter 1: Basic Concepts

PHYSICAL MECHANISM OF CONVECTION

ISCST shall not be responsible for statements or opinions contained in papers or printed in its publications.

Boundary Conditions in Fluid Mechanics

Chapter 9: Differential Analysis of Fluid Flow

University of Hail Faculty of Engineering DEPARTMENT OF MECHANICAL ENGINEERING. ME Fluid Mechanics Lecture notes. Chapter 1

7 The Navier-Stokes Equations

THE EFFECT OF SLIP CONDITION ON UNSTEADY MHD OSCILLATORY FLOW OF A VISCOUS FLUID IN A PLANER CHANNEL

CLASS SCHEDULE 2013 FALL

Shell Balances in Fluid Mechanics

Chapter 9: Differential Analysis

An-Najah National University Civil Engineering Department. Fluid Mechanics. Chapter 1. General Introduction


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

Module 9: Packed beds Lecture 29: Drag, particles settling. Flow through a packed bed of solids. Drag. Criteria of settling.

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

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fluid Mechanics Theory I

Introduction and Fundamental Concepts (Lectures 1-7)

MYcsvtu Notes HEAT TRANSFER BY CONVECTION

dynamics of f luids in porous media

4.2 Concepts of the Boundary Layer Theory

Figure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m

HEAT TRANSFER BY CONVECTION. Dr. Şaziye Balku 1

Fluid Mechanics II Viscosity and shear stresses

Module17: Intermolecular Force between Surfaces and Particles. Lecture 23: Intermolecular Force between Surfaces and Particles

1. Introduction, tensors, kinematics

Nanoscale Energy Transport and Conversion A Parallel Treatment of Electrons, Molecules, Phonons, and Photons

Contents. Preface XI Symbols and Abbreviations XIII. 1 Introduction 1

Contents. Preface XIII. 1 General Introduction 1 References 6

Contents. I Introduction 1. Preface. xiii

arxiv: v1 [physics.flu-dyn] 14 Jul 2015

Studies on flow through and around a porous permeable sphere: II. Heat Transfer

Introduction to Marine Hydrodynamics

Principles of Convection

A new approach for local similarity solutions of an unsteady hydromagnetic free convective heat transfer flow along a permeable flat surface

Sliding Contact Bearings

ME 262 BASIC FLUID MECHANICS Assistant Professor Neslihan Semerci Lecture 4. (Buoyancy and Viscosity of water)

FLOWS IN LIQUID FOAMS

Slow viscous flow in a microchannel with similar and different superhydrophobic walls

Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017


The effective slip length and vortex formation in laminar flow over a rough surface

Module 6: Free Convections Lecture 26: Evaluation of Nusselt Number. The Lecture Contains: Heat transfer coefficient. Objectives_template

Homework #4 Solution. μ 1. μ 2

UNIT II CONVECTION HEAT TRANSFER

Further Applications of Newton s Laws - Friction Static and Kinetic Friction

centrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration

Solvability condition for the moving contact line

Parash Moni Thakur. Gopal Ch. Hazarika

Modeling the combined effect of surface roughness and shear rate on slip flow of simple fluids

Driven large contact angle droplets on chemically heterogeneous substrates

Fig. 1. The coordinate system and the directions of the velocities.

ME3560 Tentative Schedule Spring 2019

FLUID MECHANICS. Chapter 9 Flow over Immersed Bodies

Lecture 2 Fluid dynamics in microfluidic systems

MHD Non-Newtonian Power Law Fluid Flow and Heat Transfer Past a Non-Linear Stretching Surface with Thermal Radiation and Viscous Dissipation

Instabilities in the Flow of Thin Liquid Films

ME3560 Tentative Schedule Fall 2018

Microfluidics 1 Basics, Laminar flow, shear and flow profiles

In the name of Allah the most beneficent the most merciful

Application of Reconstruction of Variational Iteration Method on the Laminar Flow in a Porous Cylinder with Regressing Walls

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

Chapter 2 CONTINUUM MECHANICS PROBLEMS

Lecture 8: Tissue Mechanics

Turbulent boundary layer

Table of Contents. Foreword... xiii. Preface... xv

Module 6: Free Convections Lecture 25: Buoyancy Driven Flows. The Lecture Contains: Free Convection. Objectives_template

Candidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used.

Boundary-Layer Theory

Convection Heat Transfer. Introduction

Mathematical modeling of three-layer beam hydroelastic oscillations

Exercise 5: Exact Solutions to the Navier-Stokes Equations I

Fluid Mechanics. Spring 2009

Instability & Patterning of Thin Polymer Films Prof. R. Mukherjee Department of Chemical Engineering Indian Institute Of Technology, Kharagpur

Detailed Outline, M E 521: Foundations of Fluid Mechanics I

On the breakup of fluid films of finite and infinite extent

MULTIPLE-CHOICE PROBLEMS :(Two marks per answer) (Circle the Letter Beside the Most Correct Answer in the Questions Below.)

Prediction of Coating Thickness

DIMENSIONAL ANALYSIS IN MOMENTUM TRANSFER

What s important: viscosity Poiseuille's law Stokes' law Demo: dissipation in flow through a tube

ME 309 Fluid Mechanics Fall 2010 Exam 2 1A. 1B.

Lecture-4. Flow Past Immersed Bodies

This chapter is a study of the shear stress as a function of the shear rate for Newtonian and non-newtonian biological materials.

Hydrodynamics. l. mahadevan harvard university

arxiv: v1 [physics.class-ph] 13 Sep 2008

Laminar Boundary Layers. Answers to problem sheet 1: Navier-Stokes equations

Foundations of. Colloid Science SECOND EDITION. Robert J. Hunter. School of Chemistry University of Sydney OXPORD UNIVERSITY PRESS

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

ch-01.qxd 8/4/04 2:33 PM Page 1 Part 1 Basic Principles of Open Channel Flows

Iran University of Science & Technology School of Mechanical Engineering Advance Fluid Mechanics

Heat Transfer Convection

Transfer Rate in a Reacting Flow for the Refined Sodium Bicarbonate Production Process

Problem 4.3. Problem 4.4

Transcription:

The Lecture Contains: Micro and Nano Scale Hydrodynamics with and without Free Surfaces Order of Magnitude Analysis file:///e /courses/colloid_interface_science/lecture11/11_1.htm[6/16/2012 1:39:56 PM]

Micro and Nano Scale Hydrodynamics with and without Free Surfaces When studying free surface flows, the dynamics of interfaces comes into play. So the equations of motion are important here. These equations are mentioned below. Equation of continuity Navier-Stokes equation (3.1) (3.2) There are four unknowns pressure field and three components of the velocity vector, and four equations. Therefore, theoretically, this system of equations can be solved. But due to their complexity, these equations have not been analytically solved yet. Figure 3.1. Free surface flow file:///e /courses/colloid_interface_science/lecture11/11_2.htm[6/16/2012 1:39:56 PM]

Consider a free fluid surface of variable height (Fig. 3.1) is the length scale along the -direction. Since is the rate of change in the film thickness recorded at a point. (3.3) where, and are the velocity components along and, respectively. The above-mentioned relation is called the kinematic condition file:///e /courses/colloid_interface_science/lecture11/11_3.htm[6/16/2012 1:39:56 PM]

Alternatively, the kinematic condition can also be derived as follows. Consider a pill box (Fig. 3.1). Mass balance on the pill-box leads to which can be simplified using the Leibnitz rule as the last term on the LHS vanishes because of the no slip condition at the solid boundary. file:///e /courses/colloid_interface_science/lecture11/11_4.htm[6/16/2012 1:39:57 PM]

From equation of continuity, we have therefore, In case of an impermeable solid boundary, as in the present case, which leads to. which is same as Eq (39). file:///e /courses/colloid_interface_science/lecture11/11_5.htm[6/16/2012 1:39:57 PM]

Order of magnitude analysis In this case, where, is the length scale and is the mean film thickness. Equation of continuity If we assume, then where Now if then file:///e /courses/colloid_interface_science/lecture11/11_6.htm[6/16/2012 1:39:57 PM]

Kinematic condition On LHS, the first term is, and the second is Therefore, for consistency, t should be of Navier-Stokes equation (x-component) or This means that the unsteady term is as important as the inertial terms. And, file:///e /courses/colloid_interface_science/lecture11/11_7.htm[6/16/2012 1:39:58 PM]

Navier-Stokes equation (z-component) After considering the orders of all the terms and taking the limit we get This equation shows that the pressure distribution is a hydrostatic pressure distribution with respect to z-direction. It further means that the flow is largely in the x-direction. file:///e /courses/colloid_interface_science/lecture11/11_8.htm[6/16/2012 1:39:58 PM]

Example 1. Squeezing of fluid between two flat surfaces (Fig. 3.2). Figure 3.2: Flow between two flat plates We want to obtain expressions for the velocity with which the two surfaces approach each other and the time it takes to do so. Here the distance 2h between the two plates is not a function of x and is only a function of time. file:///e /courses/colloid_interface_science/lecture11/11_9.htm[6/16/2012 1:39:58 PM]

Using we get, is to be obtained vanishes due to the no slip condition and vanishes because is independent of. file:///e /courses/colloid_interface_science/lecture11/11_10.htm[6/16/2012 1:39:58 PM]

Therefore, we get, is still an unknown. Note that is a positive quantity. file:///e /courses/colloid_interface_science/lecture11/11_11.htm[6/16/2012 1:39:59 PM]

Integrate once to get Symmetry of about mid-plane i.e. x = 0 implies that At, surrounding fluid pressure. This is what causes the symmertic drainage. file:///e /courses/colloid_interface_science/lecture11/11_12.htm[6/16/2012 1:39:59 PM]

Obtaining P in Terms of Velocity of Thinning A quasi-steady analysis leads to is the force applied by the fluid on the plate. file:///e /courses/colloid_interface_science/lecture11/11_13.htm[6/16/2012 1:39:59 PM]

Equating the two pressures gives Integrate, to get file:///e /courses/colloid_interface_science/lecture11/11_14.htm[6/16/2012 1:40:00 PM]

Analysis of the solution 1. This means that if the length is large, drainage would slow down. This is because of the viscous resistance. 2. As This may lead us to believe that the two surfaces will not come together. But actually this isn t true because when the surfaces are close together, Van-Der Waals forces will come into play and they will affect the system. When the film thickness reduces to molecular scale, the intermolecular interactions become significant. (3.6) Substituting this in the differential equation Now we can observe that as h becomes very small, velocity becomes constant. This implies that coalescence would occur in a finite time. If the exponent of h is smaller it a leads to deceleration as the surfaces approach each other and contact occurs in infinite time. On the other hand, if the exponent of h is larger it results in acceleration and the two surfaces crash into each other. Thus we can see that Lipshitz-van der Waals interactions lead to macroscopically verifiable interactions. file:///e /courses/colloid_interface_science/lecture11/11_15.htm[6/16/2012 1:40:00 PM]

Example. Liquid coating on a solid substrate (Fig. 3.3). Fig. 3.4 :Liquid film on a solid surface Navier-Stokes equation reduces to- From the no-slip condition, at We take another assumption that the bounding fluid is less viscous than the film. This would mean that the bulk fluid does not offer any drag or shear force on the film. from the no-slip condition. file:///e /courses/colloid_interface_science/lecture11/11_16.htm[6/16/2012 1:40:00 PM]

Using the continuity equation, We will use this equation to find the normal component of the velocity. Substrate (a solid) is impermeable, therefore (3.7) Conservative form of the equation (3.8) This is also the governing equation for the film thickness. Conservative form is more suitable for solving this problem numerically. file:///e /courses/colloid_interface_science/lecture11/11_17.htm[6/16/2012 1:40:00 PM]

s Why should the surface deform? 1. Because a heavier fluid is on top of a lighter fluid, 2. Because of spreading of a fluid on a non-wettable surface. A point to be noted is that these are not equilibrium situations. Figure 3.5. Thin film being coated on a surface Weak flow in the bounding fluid means that pressure distribution is hydrostatic. Also note that the mean film thickness is constant because of conservation of volume. is a constant and file:///e /courses/colloid_interface_science/lecture11/11_18.htm[6/16/2012 1:40:01 PM]

If the film thickness is very small, the potential energy due to interaction between the substrate and film needs to be taken into account. Electrostatic interactions can also be important. where accounts for the other interaction potentials. It is a function of and hence, of (3.9) This is a non-linear PDE of 4th order, which can generally be solved only using numerical techniques. file:///e /courses/colloid_interface_science/lecture11/11_19.htm[6/16/2012 1:40:01 PM]