Thin flow over a sphere
|
|
- Winfred Chandler
- 6 years ago
- Views:
Transcription
1 1 University of Waterloo, Faculty of Mathematics, Waterloo, Canada 2 Ryerson University, Department of Mathematics, Toronto, Canada IMA8 2016, June 12-16, Bad Honnef, Germany
2 Problem description Related problem Previous work Governing equations Scaling Dimensionless equations Boundary conditions Small δ expansion Leading-order problem O(δ) problem
3 Problem description Related problem Previous work The steady flow of a thin fluid layer over a sphere occurs naturally in a globe fountain as shown. A globe fountain consists of a sphere with a hole at the top whereby water is pumped out at a constant rate. An analytical study of the flow has been carried out in this study.
4 A striking example of levitation is encountered in the kugel fountain where a granite sphere, Outline sometimes weighing over a tonne, is kept aloft by a thin film of flowing water. In this paper we Problem description explain the working principle Mathematical behind this formulation levitation. We show that the fountain can be viewed as Related problem a giant roller bearing and thus Analytical forms a solution prime example Previous of lubrication work theory. It is demonstrated how the viscosity and flow rate of the fluid determine (i) the remarkably small thickness of the film supporting the sphere and (ii) the surprisingly long time it takes for rotations to damp out. The theoretical results compare well with measurements on a fountain holding a granite sphere of one meter in diameter. We close by discussing several related cases of levitation by lubrication. I. INTRODUCTION ite sphere or kugel fountains (see Fig. 1) are The kugel shown here iar sight in town squares and science museums, aller ones involves often with a massive a marble granite sphere decorate rivate homes and gardens. These fountains conperfectly polished ball floating in a socket that sphere which floats and spins cisely around on a thin it. The filmfluid ofthat flowing wells up water around of the socket is pumped into the fountain via a which is out of a hole the base. In spite of its considerable weight, the is easily atbrought the base into of a spinning the fountain. motion, which tractive sight especially when the surface of the is engraved with the Earth s map, a soccer ball, ht sky, or something of the kind. The fluid layer the socket and sphere is very thin (thinner than card 1 ), which is important for any kugel on disa public place, since it means there is no risk of s fingers being caught under the spinning sphere. ite its popularity, the granite sphere fountain is understood by most people. When we asked vis- FIG. 1: One of the largest granite sphere fountains i world, the Grand Kugel at the Science Museum of Vir Richmond. The sphere has a diameter of 2.65 m and a of about 27 tonnes.
5 Problem description Related problem Previous work Takagi & Huppert (JFM, 2010) studied a constant volume of fluid released at the top of a sphere using lubrication theory. Their analytical results agreed well with their experiments. They also investigated the onset of instability of the advancing front as it split into a series of rivulets. The dynamics of the kugel fountain was recently analysed by Snoeijer & van der Weele (Am. J. Phys., 2014), again using lubrication theory.
6 Governing equations Scaling Dimensionless equations Boundary conditions Formulated in spherical coordinates (r, θ, φ) with the hole oriented about the polar axis θ = 0 and assuming azimuthal symmetry, the governing steady-state Navier-Stokes equations become: + ν r 2 [ r + ν r 2 [ r ( r 2 v r ( r 2 v sinθ ) + (ru sinθ) = 0, r θ v v r + u r ) + 1 sinθ v u r + u r ( r 2 u r v θ u2 r θ u θ + uv r ) + 1 sinθ = 1 P ρ ( sinθ v θ θ r g cosθ ) 2 sinθ = 1 P ρr θ + g sinθ ( sinθ u ) + 2 v θ θ (u sinθ) 2v θ u ] sin 2 θ. ],
7 Governing equations Scaling Dimensionless equations Boundary conditions Here, u, v denote the velocity components in the θ and radial directions, respectively, P refers to the pressure, ν = µ/ρ represents the kinematic viscosity while µ is the dynamic viscosity and ρ is the fluid density whereas g is the acceleration due to gravity. The thickness of the fluid layer is scaled by the Nusselt thickness, H, which for a vertical incline is given by H 3 = 3νQ g, where Q is the constant flow rate exiting the small hole. For a fluid layer having a width of unity Q = UH where U is the velocity scale.
8 Governing equations Scaling Dimensionless equations Boundary conditions The coordinate y is introduced which is related to r through the relation r = R + y with R referring to the radius of the sphere. This can be scaled as r ( y ) R = 1 + δ, H where the dimensionless parameter δ = H/R 1 denotes the shallowness parameter. With this scaling the dimensionless flow variables and coordinate y are given by (u, v, P, y) (Uũ, δuṽ, ρu 2 P, Hỹ), where the tilde denotes a dimensionless quantity.
9 Governing equations Scaling Dimensionless equations Boundary conditions In dimensionless form and suppressing the tildes the governing equations become: [ (1 + δy) 2 v sinθ ] + [(1 + δy)u sinθ] = 0, y θ δ 2 v v y + δ2 u v (1 + δy) θ δu2 (1 + δy) = P y 3 cosθ δ + Re Re(1 + δy) 2 [ ( (1 + δy) 2 v ) ( + δ2 sinθ v ) 2δ ] y y sinθ θ θ sinθ θ (u sinθ) 2δ2 v δv u y + δu u (1 + δy) θ + δ2 uv (1 + δy) = δ (1 + δy) [ ( (1 + δy) 2 u ) ( + δ2 sinθ u y y sinθ θ θ where Re = Q/ν denotes the Reynolds number. P sinθ +3 θ Re + 1 Re(1 + δy) 2 ], ) + 2δ 3 v θ δ2 u sin 2 θ,
10 Governing equations Scaling Dimensionless equations Boundary conditions Along the air-fluid interface dynamic conditions along the steady free surface y = η(θ) are applied: [ 3δ 2 (η ) 2 + 2(1 + δη) 2 δ(1 + δη)η P P a = δwe [(1 + δη) 2 + δ 2 (η ) 2 ] 3/2 ] δη [ cotθ (1 + δη) 2δ + (1 + δη) 2 + δ 2 (η ) 2 Re[(1 + δη) 2 + δ 2 (η ) 2 (1 + δη) 2 v ] y + δ2 (η ) 2 ( δv + u ) ( ( u (1 + δη)η (1 + δη) θ y + δ δ v ))] (1 + δη) θ u, ( ( u [(1 + δη) 2 δ 2 (η ) 2 ] y + δ δ v )) (1 + δη) θ u ( ( )) v +2δ(1 + δη)η y δ u (1 + δη) θ + δv = 0.
11 Governing equations Scaling Dimensionless equations Boundary conditions Here, We = σh/(ρq 2 ) is the Weber number with σ denoting surface tension, P a the constant ambient air pressure, η = h/h with h being the thickness of the fluid layer, and the prime denotes differentiation with respect to θ. The kinematic condition along the free surface y = η(θ) is given by v = u (1 + δη) η, and the no-slip and impermeability conditions on the surface of the sphere are u = v = 0 at y = 0.
12 Small δ expansion Leading-order problem O(δ) problem For small δ an approximate analytical solution can be constructed by expanding the flow variables in the following series u(y, θ) = u 0 (y, θ) + δu 1 (y, θ) +, v(y, θ) = v 0 (y, θ) + δv 1 (y, θ) +, P(y, θ) = P 0 (y, θ) + δp 1 (y, θ) +, η(θ) = η 0 (θ) + δη 1 (θ) +. Substituting these expansions into the equations of motion and expanding the dynamic conditions in powers of δ leads to a hierarchy of problems.
13 Small δ expansion Leading-order problem O(δ) problem The leading-order problem is governed by the system P 0 y = 3 cosθ Re, subject to 2 u 0 y ( v0 sinθ y + u 0 θ 2 = 3 sinθ, ) + u 0 cosθ = 0, P 0 = P a, u 0 y = 0 at y = η 0, u 0 = v 0 = 0 at y = 0.
14 Small δ expansion Leading-order problem O(δ) problem The solutions are easily found to be P 0 (y, θ) = P a + 3 cosθ Re (η 0 y), u 0 (y, θ) = 3 2 y sinθ(2η 0 y), v 0 (y, θ) = y 2 2 (6η 0 cosθ + 3η 0 sinθ 2y cosθ). Thus, to leading order the pressure is hydrostatic and the velocity in the θ direction, u 0, has a parabolic profile in y which is consistent with flow down an inclined surface at an angle of θ with the horizontal. Further, the velocity component u 0 is symmetric about the plane θ = π/2.
15 Small δ expansion Leading-order problem O(δ) problem The leading-order term for the unknown free surface can be determined by applying the kinematic condition which when transferred from y = η to y = η 0 takes the form v 0 = u 0 η 0 at y = η 0. This leads to the differential equation and solution given by η 0 = 2η 0 cosθ 3 sinθ η 0 (θ) = C sin 2/3 θ. The constant C is found by applying η 0 = h 0 at θ = θ 0 where h 0 and θ 0 are free dimensionless parameters. Thus, C = h 0 sin 2/3 θ 0 and the singularity at θ = 0 is removed. Since we expect the flow to separate from the surface before reaching the bottom of the sphere, the singularity at θ = π is also resolved.
16 Small δ expansion Leading-order problem O(δ) problem The O(δ) problem satisfies the system subject to 2 u 1 y 2 sinθ = Re θ sinθ ( ) P 0 + u2 0 + (Rev 0 2) u 0 2 y, ( v1 y + u 1 θ ( 2y v 0 y + y u 0 θ + 2v 0 ) + u 1 cosθ = ) yu 0 cosθ, u 1 y + η 2 u 0 1 y 2 = u u 0 0 2η 0 y v 0 2η 0 y at y = η 0, u 1 = v 1 = 0 at y = 0.
17 Small δ expansion Leading-order problem O(δ) problem The kinematic condition applied at y = η 0 yields the differential equation for the free surface correction, η 1, given by ( u 0 η 1 + η 0 u 0 y v ) 0 η 1 = v 1 + η 0 η y 0u 0 η 0u 1, subject to η 1 = 0 at θ = θ 0. The solution to the O(δ) problem is significantly more complicated and was obtained using the Maple computer algebra system. For example, the solution for u 1 is:
18 Small δ expansion Leading-order problem O(δ) problem u 1 (y, θ) = y ( Rey 5 sinθ cosθ + 6CRey 4 sin 1/3 θ cosθ 40 ) +15CReη 0y 3 sin 4/3 θ + 60y 2 sinθ + 60y[η 0 cosθ 3C sin 1/3 θ] + F (θ) where F (θ) = 120η 1 sinθ+240c 2 sin 4/3 θ[η 0 cosθ+sinθ] 60C 4 Reη 0 sin 2/3 θ +240C(η 0) 2 sin 1/3 θ 120Cη 0 sin 2/3 θ cosθ 24C 5 Re sin 7/3 θ cosθ. The solution for η 1 can be expressed in the form η 1 (θ) = 512C 2 θ 945 sin 2/3 f (α)dα, θ θ 0 and can be integrated numerically.
19 Small δ expansion Leading-order problem O(δ) problem The function f (α) is defined as: ( f (α) = sin 23/3 α 180C 3 Re sin 2/3 α cos2α 396C 3 Re sin 2/3 α ) +280C cos3α C cosα sin 2/3 α cos3α 2625 sin 2/3 α cosα / (462+cos12α 12 cos10α+66 cos8α 220 cos6α+495 cos4α 792 cos2α).
20 0.55 The leading-order, η 0, and first-order correction, η 1, to the free surface with Re = 1, δ = 0.1, θ 0 = 0.2 and h 0 = η η 0 η 0 + δ η θ
21 1.5 The free surface η = η 0 + δη 1 shown in Cartesian coordinates with Re = 1, δ = 0.1, θ 0 = 0.2 and h 0 = 0.5. y sphere free surface x
22 The dimensionless flow rate per unit width, ˆQ, and dimensionless average streamwise velocity, Û, can be computed using the leading-order solution as follows ˆQ = η0 0 u 0 (y, θ)dy = η 3 0 sinθ, Û = ˆQ η 0 = η 2 0 sinθ. Substituting η 0 = C/ sin 2/3 θ yields ˆQ = C 3 / sinθ and Û = C 2 / sin 1/3 θ. Both the flow rate and average speed decrease as θ increases from θ 0 to π/2 which leads to a decrease in fluid thickness. Note that ˆQ sinθ is constant!
23 The leading-order, u 0, and first-order correction, u 1, to the velocity profile with Re = 1, δ = 0.1, θ 0 = 0.2 and h 0 = 1 at θ = π/4. u u 0 u 0 + δ u y
24 The velocity profile u = u 0 + δu 1 using Re = 1, δ = 0.1, θ 0 = 0.2 and h 0 = 1 at θ = π/4 and θ = π/2. u θ = π/2 θ = π/ y
25 The point of separation can be estimated by using the zero-stress condition: u y = 0 at y = 0. Using u = u 0 leads to 3η 0 (θ) sinθ = 0 and yields the separation angle θ s = π. Including the first-order correction u = u 0 + δu 1 requires solving the equation 3η 0 (θ) sinθ + δ 40 F (θ) = 0, to determine θ s. Here, F (θ) was previously defined. As expected, θ s occurs near the bottom of the sphere with little dependence on Re and δ. When h 0 = 1 and θ 0 = 0.2 separation occurs at θ s while when h 0 = 0.5 and θ 0 = 0.3 separation occurs at θ s
26 The calculation was repeated for the problem of thin flow over a cylinder and very similar streamwise velocity profiles were obtained. However, the free surface was found to vary more rapidly with θ for the sphere than it does for the cylinder. Expressions for η 0 (θ) and η 1 (θ) for the cylinder are given by η 0 (θ) = η 1 (θ) = 7560 sin 1/3 θ D sin 1/3 θ where D = h 0 sin 1/3 θ 0, D 2 θ θ 0 1 sin 11/3 α ( 245 sin 1/3 α cos3α 2205 sin 1/3 α cosα + 70D cos3α D cosα ) +72D 3 Re sin 4/3 α cos2α + 504D 3 Re sin 4/3 α dα.
27 For the cylinder the dimensionless flow rate per unit width, ˆQ, and dimensionless average streamwise velocity, Û, are given by ˆQ = D 3 and Û = ˆQ/η 0 = D 2 sin 1/3 θ. Here, the average streamwise velocity increases as θ increases from θ 0 to π/2 and since the flow rate remains constant the fluid thickness must decrease accordingly.
28 The free surface η = η 0 + δη 1 for the cylinder and sphere with Re = 1, δ = 0.1, θ 0 = 0.2 and h 0 = η cylinder sphere θ
29 Discussed in this work was an analytical investigation of the steady flow of a thin fluid layer over a sphere resulting from a constant discharge from a small hole at the top of the sphere. As an extension the cylindrical case was also solved. The variation in fluid layer thickness is more pronounced for the sphere than it is for the cylinder. This is because the average streamwise velocity increases as the fluid flows over the cylinder and since the flow rate is constant the thickness must therefore decrease. For the sphere both the flow rate and the average streamwise velocity decrease as the fluid flows over the surface which leads to a more rapid decrease in fluid layer thickness. The technique and approach adopted here can be used to model other thin flows that occur in similar settings.
Exercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationCandidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used.
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2011 2012 FLUID DYNAMICS MTH-3D41 Time allowed: 3 hours Attempt FIVE questions. Candidates must show on each answer book the type
More informationFluid 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ρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationFluid 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 informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationMA3D1 Fluid Dynamics Support Class 5 - Shear Flows and Blunt Bodies
MA3D1 Fluid Dynamics Support Class 5 - Shear Flows and Blunt Bodies 13th February 2015 Jorge Lindley email: J.V.M.Lindley@warwick.ac.uk 1 2D Flows - Shear flows Example 1. Flow over an inclined plane A
More informationFigure 11.1: A fluid jet extruded where we define the dimensionless groups
11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates
More informationFigure 3: Problem 7. (a) 0.9 m (b) 1.8 m (c) 2.7 m (d) 3.6 m
1. For the manometer shown in figure 1, if the absolute pressure at point A is 1.013 10 5 Pa, the absolute pressure at point B is (ρ water =10 3 kg/m 3, ρ Hg =13.56 10 3 kg/m 3, ρ oil = 800kg/m 3 ): (a)
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationDifferential 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 informationPF BC. When is it Energetically Favourable for a Rivulet of Perfectly Wetting Fluid to Split?
PF 05-0109-BC When is it Energetically Favourable for a Rivulet of Perfectly Wetting Fluid to Split? S. K. Wilson 1 and B. R. Duffy 2 Department of Mathematics, University of Strathclyde, Livingstone Tower,
More informationAny correspondence concerning this service should be sent to The Strathprints Administrator:
Sullivan, J.M. and Wilson, S.K. and Duffy, B.R. 2008 Air-blown rivulet flow of a perfectly wetting fluid on an inclined substrate. In: Progress in Industrial Mathematics at ECMI 2006. Springer, pp. 774-778.
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationFluid Flow in Open Helical Channels
Hayden Tronnolone Supervisor: Yvonne Stokes School of Mathematical Sciences September 24, 2010 Much research has been conducted into flows in closed, helical pipes. Such flows have application to modelling
More informationcontact line dynamics
contact line dynamics part 2: hydrodynamics dynamic contact angle? lubrication: Cox-Voinov theory maximum speed for instability corner shape? dimensional analysis: speed U position r viscosity η pressure
More informationDouble-diffusive lock-exchange gravity currents
Abstract Double-diffusive lock-exchange gravity currents Nathan Konopliv, Presenting Author and Eckart Meiburg Department of Mechanical Engineering, University of California Santa Barbara meiburg@engineering.ucsb.edu
More informationFlow past a slippery cylinder
Faculty of Mathematics, University of Waterloo, Canada EFMC12, September 9-13, 2018, Vienna, Austria Problem description and background Conformal mapping Boundary conditions Rescaled equations Asymptotic
More informationHomework of chapter (1) (Solution)
بسم اهلل الرمحن الرحيم The Islamic University of Gaza, Civil Engineering Department, Fluid mechanics-discussion, Instructor: Dr. Khalil M. Al Astal T.A: Eng. Mohammed H El Nazli Eng. Sarah R Rostom First
More informationSimulating 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 informationWhat s important: viscosity Poiseuille's law Stokes' law Demo: dissipation in flow through a tube
PHYS 101 Lecture 29x - Viscosity 29x - 1 Lecture 29x Viscosity (extended version) What s important: viscosity Poiseuille's law Stokes' law Demo: dissipation in flow through a tube Viscosity We introduced
More informationChapter 10 - Mechanical Properties of Fluids. The blood pressure in humans is greater at the feet than at the brain
Question 10.1: Explain why The blood pressure in humans is greater at the feet than at the brain Atmospheric pressure at a height of about 6 km decreases to nearly half of its value at the sea level, though
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationContents. 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 informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More informationAnalysis of Turbulent Free Convection in a Rectangular Rayleigh-Bénard Cell
Proceedings of the 8 th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows Lyon, July 2007 Paper reference : ISAIF8-00130 Analysis of Turbulent Free Convection
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More informationCandidates 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 informationUNIVERSITY 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 informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationIntroduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as:
Symbolic Calculation of Free Convection for Porous Material of Quadratic Heat Generation in a Circular Cavity Kamyar Mansour Amirkabir University of technology, Tehran, Iran, 15875-4413 mansour@aut.ac.ir
More information2. 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 informationFluid 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 informationEinstein on the beach
Einstein on the beach fluid flow on a curved space Jean-Luc Thiffeault Department of Mathematics University of Wisconsin Madison with Andrew Belmonte (Penn State), Khalid Kamhawi, and Jay Johnson (UW)
More information2 Navier-Stokes Equations
1 Integral analysis 1. Water enters a pipe bend horizontally with a uniform velocity, u 1 = 5 m/s. The pipe is bended at 90 so that the water leaves it vertically downwards. The input diameter d 1 = 0.1
More information3.1 Definition Physical meaning Streamfunction and vorticity The Rankine vortex Circulation...
Chapter 3 Vorticity Contents 3.1 Definition.................................. 19 3.2 Physical meaning............................. 19 3.3 Streamfunction and vorticity...................... 21 3.4 The Rankine
More informationBoundary Layers. Lecture 2 by Basile Gallet. 2i(1+i) The solution to this equation with the boundary conditions A(0) = U and B(0) = 0 is
Boundary Layers Lecture 2 by Basile Gallet continued from lecture 1 This leads to a differential equation in Z Z (A ib) + (A ib)[ C + Uy Uy 2i(1+i) ] = 0 with C = 2. The solution to this equation with
More information10.52 Mechanics of Fluids Spring 2006 Problem Set 3
10.52 Mechanics of Fluids Spring 2006 Problem Set 3 Problem 1 Mass transfer studies involving the transport of a solute from a gas to a liquid often involve the use of a laminar jet of liquid. The situation
More informationLaminar Boundary Layers. Answers to problem sheet 1: Navier-Stokes equations
Laminar Boundary Layers Answers to problem sheet 1: Navier-Stokes equations The Navier Stokes equations for d, incompressible flow are + v ρ t + u + v v ρ t + u v + v v = 1 = p + µ u + u = p ρg + µ v +
More informationcontact line dynamics
contact line dynamics Jacco Snoeijer Physics of Fluids - University of Twente sliding drops flow near contact line static contact line Ingbrigtsen & Toxvaerd (2007) γ γ sv θ e γ sl molecular scales macroscopic
More informationPHYS 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 informationProf. Scalo Prof. Vlachos Prof. Ardekani Prof. Dabiri 08:30 09:20 A.M 10:30 11:20 A.M. 1:30 2:20 P.M. 3:30 4:20 P.M.
Page 1 Neatly print your name: Signature: (Note that unsigned exams will be given a score of zero.) Circle your lecture section (-1 point if not circled, or circled incorrectly): Prof. Scalo Prof. Vlachos
More informationHydrodynamic Lubrication
Hydrodynamic Lubrication Experiment with 'Floating' Drops Jaywant H Arakeri and K R Sreenivas This article gives the principle of hydrodynamic lubrication and also presents the new phenomenon of levitating
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More information8.6 Drag Forces in Fluids
86 Drag Forces in Fluids When a solid object moves through a fluid it will experience a resistive force, called the drag force, opposing its motion The fluid may be a liquid or a gas This force is a very
More informationMicrofluidics 1 Basics, Laminar flow, shear and flow profiles
MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,
More informationThe 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 informationDesign and Modeling of Fluid Power Systems ME 597/ABE Lecture 7
Systems ME 597/ABE 591 - Lecture 7 Dr. Monika Ivantysynova MAHA Professor Fluid Power Systems MAHA Fluid Power Research Center Purdue University Content of 6th lecture The lubricating gap as a basic design
More informationPart II. Fluid Dynamics II. Year
Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 Paper 2, Section II 38C 41 An initially unperturbed two-dimensional inviscid jet in h < y < h has uniform speed U
More informationINTEGRAL ANALYSIS OF LAMINAR INDIRECT FREE CONVECTION BOUNDARY LAYERS WITH WEAK BLOWING FOR SCHMIDT NO. 1
INTEGRA ANAYSIS OF AMINAR INDIRECT FREE CONVECTION BOUNDARY AYERS WITH WEAK BOWING FOR SCHMIDT NO. Baburaj A.Puthenveettil and Jaywant H.Arakeri Department of Mechanical Engineering, Indian Institute of
More informationTALLINN UNIVERSITY OF TECHNOLOGY, DIVISION OF PHYSICS 13. STOKES METHOD
13. STOKES METHOD 1. Objective To determine the coefficient of viscosity of a known fluid using Stokes method.. Equipment needed A glass vessel with glycerine, micrometer calliper, stopwatch, ruler. 3.
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationAdvanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell
Laminar external natural convection on vertical and horizontal flat plates, over horizontal and vertical cylinders and sphere, as well as plumes, wakes and other types of free flow will be discussed in
More informationStudent name: This is a closed book examination. You are allowed 1 sheet of 8.5 x 11 paper with notes.
13.012 Marine Hydrodynamics for Ocean Engineers Fall 2004 Quiz #2 Student name: This is a closed book examination. You are allowed 1 sheet of 8.5 x 11 paper with notes. For the problems in Section A, fill
More informationCHME 302 CHEMICAL ENGINEERING LABOATORY-I EXPERIMENT 302-V FREE AND FORCED CONVECTION
CHME 302 CHEMICAL ENGINEERING LABOATORY-I EXPERIMENT 302-V FREE AND FORCED CONVECTION OBJECTIVE The objective of the experiment is to compare the heat transfer characteristics of free and forced convection.
More information1. Introduction - Tutorials
1. Introduction - Tutorials 1.1 Physical properties of fluids Give the following fluid and physical properties(at 20 Celsius and standard pressure) with a 4-digit accuracy. Air density : Water density
More informationNonlinear Analysis: Modelling and Control, 2008, Vol. 13, No. 4,
Nonlinear Analysis: Modelling and Control, 2008, Vol. 13, No. 4, 513 524 Effects of Temperature Dependent Thermal Conductivity on Magnetohydrodynamic (MHD) Free Convection Flow along a Vertical Flat Plate
More information- 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 informationThermocapillary 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 informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More information1/3/2011. This course discusses the physical laws that govern atmosphere/ocean motions.
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationPhys 1401: General Physics I
1. (0 Points) What course is this? a. PHYS 1401 b. PHYS 1402 c. PHYS 2425 d. PHYS 2426 2. (0 Points) Which exam is this? a. Exam 1 b. Exam 2 c. Final Exam 3. (0 Points) What version of the exam is this?
More informationChapter 10. Rotation
Chapter 10 Rotation Rotation Rotational Kinematics: Angular velocity and Angular Acceleration Rotational Kinetic Energy Moment of Inertia Newton s nd Law for Rotation Applications MFMcGraw-PHY 45 Chap_10Ha-Rotation-Revised
More informationHeat and Mass Transfer
1 Comments on six papers published by S.P. Anjali Devi and R. Kandasamy in Heat and Mass Transfer, ZAMM, Mechanics Research Communications, International Communications in Heat and Mass Transfer, Communications
More informationMM303 FLUID MECHANICS I PROBLEM SET 1 (CHAPTER 2) FALL v=by 2 =-6 (1/2) 2 = -3/2 m/s
MM303 FLUID MECHANICS I PROBLEM SET 1 (CHAPTER ) FALL 018 1) For the velocity fields given below, determine: i) Whether the flow field is one-, two-, or three-dimensional, and why. ii) Whether the flow
More informationthe Cartesian coordinate system (which we normally use), in which we characterize points by two coordinates (x, y) and
2.5.2 Standard coordinate systems in R 2 and R Similarly as for functions of one variable, integrals of functions of two or three variables may become simpler when changing coordinates in an appropriate
More information15 The coffee cup J Nonlinear Dynamics II: Continuum Systems Lecture 15 Spring 2015
18.354J Nonlinear Dynamics II: Continuum Systems Lecture 15 Spring 2015 15 The coffee cup Let s try and apply our knowledge of fluid dynamics to a real observation, to test whether the theory actually
More informationBasic 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 informationThe Navier-Stokes equations for an incompressible fluid in indicial notation are. i u i = 0 (9.1)
Chapter 9 Viscous flow The Navier-Stokes equations for an incompressible fluid in indicial notation are i u i = 0 9.1 ρ t + u j j u i = i p + µ 2 j u i 9.2 When these are applied to practical problems,
More informationShallow free-surface flows.
Shallow free-surface flows. Salil Gogte Department of Mechanical Engineering, University of New Mexico Albuquerque, NM 87131, USA March 4, 2005 Abstract The analysis of thin film flows with large distortions
More informationAP Physics First Nine Weeks Review
AP Physics First Nine Weeks Review 1. If F1 is the magnitude of the force exerted by the Earth on a satellite in orbit about the Earth and F2 is the magnitude of the force exerted by the satellite on the
More informationChapter 5. Shallow Water Equations. 5.1 Derivation of shallow water equations
Chapter 5 Shallow Water Equations So far we have concentrated on the dynamics of small-scale disturbances in the atmosphere and ocean with relatively simple background flows. In these analyses we have
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationTOPIC D: ROTATION EXAMPLES SPRING 2018
TOPIC D: ROTATION EXAMPLES SPRING 018 Q1. A car accelerates uniformly from rest to 80 km hr 1 in 6 s. The wheels have a radius of 30 cm. What is the angular acceleration of the wheels? Q. The University
More informationMoments of Inertia (7 pages; 23/3/18)
Moments of Inertia (7 pages; 3/3/8) () Suppose that an object rotates about a fixed axis AB with angular velocity θ. Considering the object to be made up of particles, suppose that particle i (with mass
More informationSimulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions
Chapter 1 Earth Science Simulation Study on the Generation and Distortion Process of the Geomagnetic Field in Earth-like Conditions Project Representative Yozo Hamano Authors Ataru Sakuraba Yusuke Oishi
More informationREE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics
REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationComputation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
More informationSpinning of a molten threadline Steady-state isothermal viscous flows
Spinning of a molten threadline Steady-state isothermal viscous flows Jet equations and shape M.A. Matovich and J.R.A. Pearson (1969) Matthieu Varagnat July 13, 2007 Summer Reading Group The authors J.R.A.
More informationCENG 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 informationFlow and instability of thin films on a cylinder and sphere
J. Fluid Mech. (2010), vol. 647, pp. 221 238. c Cambridge University Press 2010 doi:10.1017/s0022112009993818 221 Flow and instability of thin films on a cylinder and sphere DAISUKE TAKAGI AND HERBERT
More informationExercise solutions: concepts from chapter 10
Reading: Fundamentals of Structural Geology, Ch 0 ) The flow of magma with a viscosity as great as 0 0 Pa s, let alone that of rock with a viscosity of 0 0 Pa s, is difficult to comprehend because our
More informationSuppression of 3D flow instabilities in tightly packed tube bundles
Suppression of 3D flow instabilities in tightly packed tube bundles Nicholas Kevlahan kevlahan@mcmaster.ca Department of Mathematics & Statistics CSFD, June 13 15 2004 p.1/33 Collaborators CSFD, June 13
More informationINTRODUCTION TO FLUID MECHANICS June 27, 2013
INTRODUCTION TO FLUID MECHANICS June 27, 2013 PROBLEM 3 (1 hour) A perfect liquid of constant density ρ and constant viscosity µ fills the space between two infinite parallel walls separated by a distance
More informationTHERMOCAPILLARY CONVECTION IN A LIQUID BRIDGE SUBJECTED TO INTERFACIAL COOLING
THERMOCAPILLARY CONVECTION IN A LIQUID BRIDGE SUBJECTED TO INTERFACIAL COOLING Melnikov D. E. and Shevtsova V. M. Abstract Influence of heat loss through interface on a supercritical three-dimensional
More informationModule 3: "Thin Film Hydrodynamics" Lecture 11: "" The Lecture Contains: Micro and Nano Scale Hydrodynamics with and without Free Surfaces
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]
More informationMECHANICAL PROPERTIES OF FLUIDS:
Important Definitions: MECHANICAL PROPERTIES OF FLUIDS: Fluid: A substance that can flow is called Fluid Both liquids and gases are fluids Pressure: The normal force acting per unit area of a surface is
More informationQuiz #8. Vector. 2) Given A( 1, 4, 3), and B( 3, 4, 1), calculate A B
Quiz #8 Vector 1) Given A(1, 2), and B( 3, 4), calculate A B 2) Given A( 1, 4, 3), and B( 3, 4, 1), calculate A B 3) Given the following magnitude of forces in Figure 1: α = 20, θ = 60, β = 30, F 1 = 1N,
More informationConvection. 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 informationInfluence of Lock Aspect Ratio upon the Evolution of an Axisymmetric Intrusion
Under consideration for publication in J. Fluid Mech. 1 Influence of Lock Aspect Ratio upon the Evolution of an Axisymmetric Intrusion By AMBER M. HOLDSWORTH and BRUCE R. SUTHERLAND Departments of Physics
More informationDepartment of Civil Engineering Hydraulics and Water Resources Division Application and Solution I
Question 1: The Specific weight of water is 1000 /. Using this given value, find the specific mass of water in SI units (g= m/s ). Solution 1: The specific mass of water in SI units: 1 N 1000 m 9810 Nm
More informationHydrodynamics of wetting phenomena. Jacco Snoeijer PHYSICS OF FLUIDS
Hydrodynamics of wetting phenomena Jacco Snoeijer PHYSICS OF FLUIDS Outline 1. Creeping flow: hydrodynamics at low Reynolds numbers (2 hrs) 2. Thin films and lubrication flows (3 hrs + problem session
More informationChapter 6 Laminar External Flow
Chapter 6 aminar Eternal Flow Contents 1 Thermal Boundary ayer 1 2 Scale analysis 2 2.1 Case 1: δ t > δ (Thermal B.. is larger than the velocity B..) 3 2.2 Case 2: δ t < δ (Thermal B.. is smaller than
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More information