STEADY VISCOUS FLOW THROUGH A VENTURI TUBE
|
|
- Derick Douglas
- 5 years ago
- Views:
Transcription
1 CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 12, Number 2, Summer 2004 STEADY VISCOUS FLOW THROUGH A VENTURI TUBE K. B. RANGER ABSTRACT. Steady viscous flow through an axisymmetric convergent-divergent nozzle is discussed for a wide range of viscosities. It is shown that in general the Stokes flow represents a uniform approximation to the flow for arbitrary viscosity. Introduction The present article discusses the steady axisymmetric viscous flow through a convergent-divergent nozzle. The model for the nozzle is determined by a suitable conformal mapping and is represented by one family of coordinate curves depicted by β = constant. In the first section the Stokes flow is found to be a similarity solution of the coordinate β. This result was previously derived by Sampson [5], Roscoe [4] and is described in the book by Happel and Brenner [3]. The method of derivation given here is different from other authors. It is possible to calculate the pressure drop and wall stress over a section of the nozzle. The main thrust in this discussion is to show that even though the Stokes flow is symmetric it represents a uniform approximation to the motion for a wide range of viscosities. This is achieved by considering coordinate expansions for effectively small and large distances from the nozzle centre. The expansion for large distances is similar to that in [6]. The wall stress and pressure drop in a section of the channel are also valid for a wide range of viscosity and not just for Stokes flow. The applications of the Venturi tube have been discussed in [1]. 1 Axisymmetric Stokes flow through a constricted channel If (z, ρ, φ are cylindrical polar coordinates and ψ = ψ(z, ρ is the stream function for steady axisymmetric motion, the fluid velocity is described AMS subject classification: 34K20, 92D25. Copyright c Applied Mathematics Institute, University of Alberta. 199
2 200 K. B. RANGER by (1 q = curl { ψ } ρ φ = 1 ρ ψ ρ k + 1 ρ ψ z ρ, and for Stokes flow which plays a central role in this analysis the momentum equations can be written in the form (2 p z = µ ρ ρ (L 1(ψ, p ρ = µ ρ z (L 1(ψ, where p = p(z, ρ is the pressure field, and µ the coefficient of viscosity. The Stokes operator is defined by ( 2 (3 L 1 (ψ = z ρ 2 1 ψ. ρ ρ If the pressure p is eliminated from equations (2, then the stream function ψ satisfies the Stokes or repeated operator equation (4 L 2 1(ψ = L 1 (L 1 (ψ = 0. There is a decomposition formula for the solution of iterated operator equations first given by Weinstein [2] and is expressible in the form (5 ψ = v ( 1 + v ( 3, where v ( k, k = 1, 3, satisfy the equation { (6 L k v ( k} = v zz ( k + v ρρ ( k k ρ v( k ρ = 0. v ( 1 and v ( 3 may be interpreted as potential flow stream functions in axisymmetric spaces of three and five dimensions, respectively. Now if z + iρ = F (α + iβ, F (α + iβ 0 is a conformal mapping of the meridional plane (z, ρ into a region of the (α, β plane, then v ( k satisfies the equation (7 ρ k {ρ α k v( k α } + ρ k {ρ β k v( k β } = 0,
3 STEADY VISCOUS FLOW THROUGH A VENTURI TUBE 201 and (8 L 1 (ψ 1 = F (α + iβf (α iβ { ρ α ρ ψ + ρ α β Also the momentum equations (2 transform in the form (9 p α = µ ρ β (L 1(ψ, Oblate ellipsoidal coordinates ρ p β = µ ρ α (L 1(ψ. The conformal mapping (10 z + iρ = F (α + iβ = c sinh(α + iβ produces the real equations (11 z = c sinh α cos β, ρ = c cosh α sin β and the surfaces β = constant, described by ( 2 ( 2 ρ z (12 = c 2 sin β cos β } ψ. β represent a system of confocal hyperboloids of one sheet with common focal ring located at z = 0, ρ = c. The surfaces α = constant, are the orthogonal trajectories and represent a confocal set of planetary or oblate ellipsoids also with common ring at z = 0, ρ = c. The range of the coordinates (α, β are given by 0 β < π/2, < α <. This coordinate system is appropriate to the description of flow through a convergent-divergent nozzle, and for the purposes of application describes a section of stenotic artery. A diagram of the coordinate system in the meridional plane together with a sketch of the flow model are given below in Figure 1 and Figure 2. 2 The boundary value problem for steady Stokes flow The stream function ψ satisfies equation (4 and for the present purposes an appropriate decomposition formula for the solution is described by (5. The boundary conditions of no seepage and no slip on the boundary together with finite flow on the axis are represented by (13 (14 ψ = M, ψ β = 0 at β = β 0, ψ = 0, at β = 0 for < α <.
4 202 K. B. RANGER ρ β=const. α=const. z FIGURE 1: Figure 1. Sketch of typical curves in the meridional plane x z y FIGURE 2: Figure 2. Sketch of the convergent-divergent nozzle
5 STEADY VISCOUS FLOW THROUGH A VENTURI TUBE 203 M is the constant volumetric flux rate through the nozzle. The equations satisfied by v ( 1, v ( 3 are; (15 (16 cosh α cosh α α cosh 3 α α cosh 3 α v ( 1 α v ( 3 α + sin β β + sin 3 β β sin β sin 3 β v ( 1 β v ( 3 β = 0, = 0. Suitable solutions for v ( 1, v ( 3 are forced by the boundary conditions, and in this case v ( 1 = v ( 1 (β, v ( 3 = v ( 3 (β so that (17 d dβ sin β dv ( 1 dβ The solutions are given by = 0, d dβ sin 3 β dv ( 3 dβ = 0. (18 v ( 1 (β = a cos β + b, v ( 3 (β = c cos β c 3 cos3 β + d, where a, b, c and d are constants. The stream function ψ can then be more conveniently written as; (19 ψ = A cos β + C cos 3 β + D where A, C andd are constants to be determined by the boundary conditions. These conditions require that A, C and D satisfy (20 A cos β 0 + C cos 3 β 0 + D = M, A + 3C cos 2 β 0 = 0, A + C + D = 0 and the complete Stokes solution is defined by (21 ψ = M [ cos 3 β 1 3 cos 2 β 0 (cos β 1 ] [3 cos 2 β cos 3. β 0 ] The stream function is a similarity solution of the Stokes equations since the streamlines coincide with the curves β = constant. This flow was previously derived by Sampson, Roscoe and is also contained in the book by Happel and Brenner. The method of derivation in [3] is different from that given here, and the analysis presented in this section appears to be the first application of the decomposition formula (5 to this problem.
6 204 K. B. RANGER The components of fluid velocity u, v along and perpendicular to the symmetry axis respectively are described by (22 (23 u = 1 ρ ψ ρ = 1 ρ β ρψ β = 3M cos β ( cos 2 β 0 cos 2 β c 2 ( sinh 2 α + cos 2 β (3 cos 2 β cos 3 β 0, v = 1 ρ ψ z = 1 ρ β zψ β = 3M sinh α ( cos 2 β cos 2 β 0 sin β c 2 cosh α ( sinh 2 α + cos 2 β [3 cos 2 β cos 3 β 0 ]. The velocity on the axis at the origin α = β = 0 is given by 3M ( 1 cos 2 β 0 (24 u = c 2 [3 cos 2 β cos 3 β 0 ] and decays like e 2 α or 1/r 2 as α. defined by the formula The fluid vorticity ω is (25 ω = L 1(ψ ρ 1 = c 3 cosh α sin β(sinh 2 α + cos 2 β { cosh ψ α cosh α α = sin β β sin β } ψ β 6M cos β sin β c 3 cosh α(sinh 2 α + cos 2 β [3 cos 2 β cos 3 β 0 ]. In particular the vorticity on the boundary β = β 0 is expressed by (26 ω = 6M sin β 0 cos β 0 c 3 cosh α ( sinh 2 α + cos 2 β 0 [3 cos2 β cos 3 β 0 ], and is a maximum at α = 0, and vanishes like O(e 3 α as α. The maximum is given by (27 ω = 6M tan β 0 c 3 [3 cos 2 β cos 3 β 0 ],
7 STEADY VISCOUS FLOW THROUGH A VENTURI TUBE 205 which is an increasing function of β 0 in the range 0 < β 0 < π/2. Since the fluid velocity vanishes at β = β 0, it follows ψ ββ ω = ρc 2 (sinh 2, α + cos 2 β β=β 0 so that the boundary vorticity is related to the shear stress on the boundary in a simple way. The shear stress is negative for all α on β = β 0, so that there is no possibility of streamline separation attached to the boundary. Finally it is noted that the numerical comparison can be substantiated by observing that the pressure force, A p = 1 (sec β The pressure field It has been shown that in (α, β coordinates the momentum equations are equivalent to (28 where (29 L 1 (ψ = p α = µ c cosh α sin β β L 1(ψ, p β = µ c cosh α sin β α L 1(ψ 6M cos β sin 2 β c 2 (sinh 2 α + cos 2 β[3 cos 2 β cos 3 β 0 ]. By elimination of L 1 (ψ from these equations, the pressure field satisfies the Laplace equation and in (α, β coordinates ( 1 (30 cosh α p + 1 ( sin β p = 0. cosh α α α sin β β β Now if s = cos β, then from (28 it is found that p (31 s = 6Mµ ( s 1 c 3 cosh α α sinh 2 α + s 2 (3 cos 2 β cos 3 β 0. Integration shows that (32 p = 3Mµ c 3 cosh α[3 cos 2 β cos 3 β 0 ] { ln(sinh 2 α + s 2 } + A(α α 6Mµ sinh α = c 3 [3 cos 2 β cos 3 β 0 ][sinh 2 α + cos 2 β] + A(α,
8 206 K. B. RANGER where A(α can be taken as a constant since it is readily verified that { cosh α + 1 ( sin β } (33 cosh α α α sin β β β sinh α (sinh 2 α + cos 2 β = 0. The axial pressure drop in a section of nozzle from α = +α 0 to α = α 0, z 0 = c sinh α 0, is given by (34 p = p α=α0 p α= α0 = 12Mµz 0 c 2 (z c2 [2 cos 3 β cos 2 β 0 ]. In particular the pressure drops in the channel along the axis for β 0 = π/6, β 0 = π/4, β 0 = π/3 are given by (35 p π/6 = k 0.049, p π/4 = k 0.207, p π/3 = k 0.5, where the constant k = 12Mz 0 µ/(c 2 (z c2. The ratio of these drops are represented by (36 p π/6 p π/4 = 4.224, p π/6 p π/3 = The minimum area of the cross-sections corresponding to these drops are: (37 A π/6 = 1 4 πc2, A π/4 = 1 2 πc2, A π/3 = 3 4 πc2 and the ratios are given by (38 A π/3 A π/6 = 3, A π/4 A π/6 = 2. These results indicate that the ratio of the pressure drops are slightly greater than the square of the cross-sectional ratios. Again the volume of liquid in the channel bounded by z 0 z z 0, the axis β = 0 and the meridian curve β = β 0, is represented by z0 (39 V = π ρ 2 dz, ρ 2 = sin 2 β 0 (c 2 + z2 z 0 cos 2 β 0 = 2π sin 2 β 0 (z 0 c 2 z cos 2 β 0
9 STEADY VISCOUS FLOW THROUGH A VENTURI TUBE 207 ψ=μ Fig. 1 z=-z 0 ψ=0 z=z 0 β 0 =π/6 Fig. 2 z=-z 0 ψ=0 z=z 0 β 0 =π/4 Fig. 3 z=-z 0 ψ=0 z=z 0 β 0 =π/3
10 208 K. B. RANGER and the pressure drop along the axis from z = z 0 to z = z 0, is given by (40 p = p z=z0 p z= z0 = 12Mz 0 µ c 2 (z c2 [2 cos 3 β cos 2 β 0 ]. This result can be contrasted with the corresponding values determined from steady Poiseulle flow through a straight pipe where the volumetric flux rate is the same and is applied to a section of the pipe with the same volume of liquid. If the radius of the pipe is R, then (41 V = 2πR 2 z 0, R 2 = sin 2 β 0 (c 2 z cos 2 β 0 and from the known results for Poiseulle flow available in the standard texts (42 M = πp R4 9µ, P = dp dz where the pressure gradient is constant. The pressure drop is then 16Mz 0 µ (43 p = 2P z 0 = 2, π sin 4 β 0 (c 2 + z2 0 3 cos 2 β 0 and may be compared with the motion in a constricted channel by the ratio (44 S = 12Mπz 0 µ sin 4 β 0 (c 2 + z 2 0 /3 cos2 β 0 2 c 2 (c 2 + z Mµz 0(2 cos 3 β cos 2 β 0, which is an increasing function of β 0 in the range 0 < β 0 < π/2. In fact as β 0 π/2, S. 4 Steady state Navier-Stokes equations For axisymmetric flow the fluid velocity field is defined by (45 q = curl { ψ } ρ φ = 1 ρ ψ ρ k + 1 ρ ψ z ρ
11 STEADY VISCOUS FLOW THROUGH A VENTURI TUBE 209 and the momentum equations can be written in the form (46 (47 1 ρ 2 ψ zl 1 (ψ = B z ν ρ 1 ρ 2 ψ ρl 1 (ψ = B ρ + ν ρ ρ L 1(ψ, z L 1(ψ where B = p/ρ 0 + (1/2 q 2 is the Bernoulli function or total head of pressure, and ν is the kinematic viscosity. Elimination of B results in the vorticity equation (48 ρ { ψ, L 1 (ψ/ρ 2} (z, ρ = νl 2 1 (ψ. In oblate ellipsoidal coordinates z + iρ = c sinh(α + iβ equations (46 (47 transform into: (49 (50 1 ρ 2 ψ αl 1 (ψ = B α ν ρ 1 ρ 2 ψ βl 1 (ψ = B β + ν ρ β (L 1(ψ, α (L 1(ψ and also there is corresponding mapping for the vorticity equation (51 ρ {ψ, L 1(ψ/ρ} (α, β = ν { cosh α α cosh α α + sin β ( } 1 L 1(ψ β sin β β where (52 { 1 L 1 (ψ = c 2 (sinh 2 α + cos 2 β cosh α α cosh α ψ α + sin β ( } 1 ψ. β sin β β
12 210 K. B. RANGER The coordinate expansion for large α In this situation the stream function ψ is expanded in the asymptotic series (53 ψ F n (βe nα, α, with boundary conditions (54 n=0 F 0 (β 0 = M, F n (β 0 = 0, n 1, F n(β 0 = 0, n 0 and F n (0 = 0 for n 0. For large α the leading order approximations show that (55 ρ { ψ, L 1 (ψ/ρ 2} = O(e 3α, α (α, β and (56 [ cosh α α cosh α + sin β α β sin β ] L 1 (ψ β = O(e 2α, α. Similar results are also determined as α, and in both cases the leading term in the expansion is the Stokes flow which satisfies the equation (57 [ sin β β sin β ] + 6 sin β ( 1 β β sin β F 0 (β = 0. The solution which produces finite flow on the axis is given by (58 F 0 (β = M[cos3 β 1 3 cos 2 β 0 (cos β 1] (3 cos 2 β cos 3. β 0 Higher order approximations satisfy a system of inhomogenous linear partial differential equations of fourth order which in principle can be constructed successively. It is noted from equations (49, (50 that for consistency the components of convective-acceleration show that 1 ρ 2 ψ αl 1 (ψ + 1 [ (ψα 2 + ψ β 2 ] (59 2 α ρ 2 c 2 (sinh 2 = O(e 4α, α + cos 2 β (60 [ 1 ρ 2 ψ βl 1 (ψ β ψ 2 α + ψ2 β ρ 2 c 2 (sinh 2 α + cos 2 β ] = O(e 4α
13 STEADY VISCOUS FLOW THROUGH A VENTURI TUBE 211 and the components of viscous diffusion: (61 ν ρ β (L 1(ψ = O(e 3α, ν ρ α (L 1(ψ = O(e 3α for large α. This demonstrates that the ratio convection/diffusion = O(e α, α and the Stokes approximation is valid at large distances from the nozzle. Similar results apply to higher order approximations since the motion is finite in the flow region. A similar result can be obtained as α. It is not clear whether the leading order term is a good approximation for small α. However it is observed from the exact vorticity equation that on the boundary β = β 0, convectiveacceleration is absent and the governing equation is viscous diffusion, that is L 2 1 (ψ = 0. It then appears that the complete Stokes solution is a leading order approximation to the flow in a neighbourhood of the nozzle wall and is independent of viscosity. A second point is that if ψ = F 0 (β, then the convective-acceleration term from (48 is given by (62 ρ [ψ, L 1(ψ/ρ 2 ] (α, β Also the diffusion term [ (63 cosh α α cosh α = 0 at α = 0. + sin β α β sin β [ F ] β 0 (β cot βf 0(β sinh 2 α + cos 2 β vanishes independent of viscosity. It follows from the vorticity equation that the Stokes solution represents the leading order term in a coordinate expansion for small α described by (64 ψ F 0 (β + where α n G m (β, α 0 m=1 (65 G m (β 0 = G m (β 0 = G m (0 = 0 for m 1. Finally it is reasonable to assume that although the flow is asymmetric about the plane α = 0, that the Stokes solution is a uniformily valid leading order approximation for the flow independent of α and valid for arbitrary viscosity. ]
14 212 K. B. RANGER 5 Conclusions The results for the pressure drop and wall-stress over a section of the flow through a convergent-divergent nozzle appear to be new as the analysis suggest they are valid over a wide range of viscosity. The possibility of application to the flow of blood through a large artery is also included but since the motion is pulsatile the results would have to be interpreted in some sense as mean values over a suitable time period. Since the leading order approximation is symmetric about the plane of symmetry, but the complete flow is asymmetric there is the possibility of reverse flow occuring on a region of the wall downstream of the nozzle. However this effect is likely to be small compared with the main forcing flow. REFERENCES 1. G. Batchelor, An Introduction to Fluid Dynamics, Cambridge University Press, p. 484, J. Happel and H. Brenner, Low Reynolds number hydrodynamics, Martenus Hijoff Publishers, , J. Happel and and H. Brenner, Low Reynolds number hydrodynamics, Martenus Hijhoff Publishers, p. 153, R. Roscoe, Phil. Mag. 40 (1949, R. A. Sampson, Phil Trans. Roy. Soc. A182 (1891, A. Weinstein, Annali Di Mathematics Pure Et Applicata 1959 iv., 39. Department of Mathematics, University of Toronto, Toronto, Ontario, Canada, M5S 3G3 address: ranger@math.utoronto.ca
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 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 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 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 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 information3.5 Vorticity Equation
.0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous
More informationAA210A 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 informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
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 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 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 informationDetailed Outline, M E 320 Fluid Flow, Spring Semester 2015
Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous
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 information6.1 Momentum Equation for Frictionless Flow: Euler s Equation The equations of motion for frictionless flow, called Euler s
Chapter 6 INCOMPRESSIBLE INVISCID FLOW All real fluids possess viscosity. However in many flow cases it is reasonable to neglect the effects of viscosity. It is useful to investigate the dynamics of an
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 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 informationPage 1. Neatly print your name: Signature: (Note that unsigned exams will be given a score of zero.)
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. Vlachos Prof. Ardekani
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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationLecture 3: The Navier-Stokes Equations: Topological aspects
Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without
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 informationV/ t = 0 p/ t = 0 ρ/ t = 0. V/ s = 0 p/ s = 0 ρ/ s = 0
UNIT III FLOW THROUGH PIPES 1. List the types of fluid flow. Steady and unsteady flow Uniform and non-uniform flow Laminar and Turbulent flow Compressible and incompressible flow Rotational and ir-rotational
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationMicroscopic Momentum Balance Equation (Navier-Stokes)
CM3110 Transport I Part I: Fluid Mechanics Microscopic Momentum Balance Equation (Navier-Stokes) Professor Faith Morrison Department of Chemical Engineering Michigan Technological University 1 Microscopic
More informationME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics
ME 431A/538A/538B Homework 22 October 2018 Advanced Fluid Mechanics For Friday, October 26 th Start reading the handout entitled Notes on finite-volume methods. Review Chapter 7 on Dimensional Analysis
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
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 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 informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationChemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017
Chemical and Biomolecular Engineering 150A Transport Processes Spring Semester 2017 Objective: Text: To introduce the basic concepts of fluid mechanics and heat transfer necessary for solution of engineering
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
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 informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationQuick 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 informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationVorticity Equation Marine Hydrodynamics Lecture 9. Return to viscous incompressible flow. N-S equation: v. Now: v = v + = 0 incompressible
13.01 Marine Hydrodynamics, Fall 004 Lecture 9 Copyright c 004 MIT - Department of Ocean Engineering, All rights reserved. Vorticity Equation 13.01 - Marine Hydrodynamics Lecture 9 Return to viscous incompressible
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 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 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 informationCHAPTER 3 BASIC EQUATIONS IN FLUID MECHANICS NOOR ALIZA AHMAD
CHAPTER 3 BASIC EQUATIONS IN FLUID MECHANICS 1 INTRODUCTION Flow often referred as an ideal fluid. We presume that such a fluid has no viscosity. However, this is an idealized situation that does not exist.
More informationvector H. If O is the point about which moments are desired, the angular moment about O is given:
The angular momentum A control volume analysis can be applied to the angular momentum, by letting B equal to angularmomentum vector H. If O is the point about which moments are desired, the angular moment
More informationUNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence
UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: Day of examination: Friday 15. June 212 Examination hours: 9. 13. This problem set consists of 5 pages. Appendices: Permitted
More information0 = p. 2 x + 2 w. z +ν w
Solution (Elliptical pipe flow (a Using the Navier Stokes equations in three dimensional cartesian coordinates, given that u =, v = and w = w(x,y only, and assuming no body force, we are left with = p
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 informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationLecture 30 Review of Fluid Flow and Heat Transfer
Objectives In this lecture you will learn the following We shall summarise the principles used in fluid mechanics and heat transfer. It is assumed that the student has already been exposed to courses in
More informationCE Final Exam. December 12, Name. Student I.D.
CE 100 - December 12, 2009 Name Student I.D. This exam is closed book. You are allowed three sheets of paper (8.5 x 11, both sides) of your own notes. You will be given three hours to complete four problems.
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 informationVisualization of flow pattern over or around immersed objects in open channel flow.
EXPERIMENT SEVEN: FLOW VISUALIZATION AND ANALYSIS I OBJECTIVE OF THE EXPERIMENT: Visualization of flow pattern over or around immersed objects in open channel flow. II THEORY AND EQUATION: Open channel:
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 informationChapter 4 DYNAMICS OF FLUID FLOW
Faculty Of Engineering at Shobra nd Year Civil - 016 Chapter 4 DYNAMICS OF FLUID FLOW 4-1 Types of Energy 4- Euler s Equation 4-3 Bernoulli s Equation 4-4 Total Energy Line (TEL) and Hydraulic Grade Line
More informationMeasurements using Bernoulli s equation
An Internet Book on Fluid Dynamics Measurements using Bernoulli s equation Many fluid measurement devices and techniques are based on Bernoulli s equation and we list them here with analysis and discussion.
More informationMAE 3130: Fluid Mechanics Lecture 7: Differential Analysis/Part 1 Spring Dr. Jason Roney Mechanical and Aerospace Engineering
MAE 3130: Fluid Mechanics Lecture 7: Differential Analysis/Part 1 Spring 2003 Dr. Jason Roney Mechanical and Aerospace Engineering Outline Introduction Kinematics Review Conservation of Mass Stream Function
More informationChapter 3 Bernoulli Equation
1 Bernoulli Equation 3.1 Flow Patterns: Streamlines, Pathlines, Streaklines 1) A streamline, is a line that is everywhere tangent to the velocity vector at a given instant. Examples of streamlines around
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 informationFluid Mechanics Qualifying Examination Sample Exam 2
Fluid Mechanics Qualifying Examination Sample Exam 2 Allotted Time: 3 Hours The exam is closed book and closed notes. Students are allowed one (double-sided) formula sheet. There are five questions on
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 informationVorticity Stream function Solver in Cylindrical Coordinates
Vticity Stream function Solver in Cylindrical Codinates L. Brieda May 27, 2016 This document summarizes equations used to solve flow in a cylindrical pipe using the stream function approach as seen in
More informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
More informationCOURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics
COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour Basic Equations in fluid Dynamics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Description of Fluid
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 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 informationVortex motion. Wasilij Barsukow, July 1, 2016
The concept of vorticity We call Vortex motion Wasilij Barsukow, mail@sturzhang.de July, 206 ω = v vorticity. It is a measure of the swirlyness of the flow, but is also present in shear flows where the
More informationCBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates
CBE 6333, R. Levicky 1 Orthogonal Curvilinear Coordinates Introduction. Rectangular Cartesian coordinates are convenient when solving problems in which the geometry of a problem is well described by the
More information2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;
Accepted for publication in J. Fluid Mech. 1 Viscous Potential Flow By D.D. Joseph Department of Aerospace Engineering and Mechanics, University of Minnesota, MN 55455 USA Email: joseph@aem.umn.edu (Received
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationChapter 2 CONTINUUM MECHANICS PROBLEMS
Chapter 2 CONTINUUM MECHANICS PROBLEMS The concept of treating solids and fluids as though they are continuous media, rather thancomposedofdiscretemolecules, is one that is widely used in most branches
More informationThe effective slip length and vortex formation in laminar flow over a rough surface
The effective slip length and vortex formation in laminar flow over a rough surface Anoosheh Niavarani and Nikolai V. Priezjev Movies and preprints @ http://www.egr.msu.edu/~niavaran A. Niavarani and N.V.
More informationHomework #4 Solution. μ 1. μ 2
Homework #4 Solution 4.20 in Middleman We have two viscous liquids that are immiscible (e.g. water and oil), layered between two solid surfaces, where the top boundary is translating: y = B y = kb y =
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 informationSignature: (Note that unsigned exams will be given a score of zero.)
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. Dabiri Prof. Wassgren Prof.
More informationReynolds, an engineering professor in early 1880 demonstrated two different types of flow through an experiment:
7 STEADY FLOW IN PIPES 7.1 Reynolds Number Reynolds, an engineering professor in early 1880 demonstrated two different types of flow through an experiment: Laminar flow Turbulent flow Reynolds apparatus
More informationProblem 4.3. Problem 4.4
Problem 4.3 Problem 4.4 Problem 4.5 Problem 4.6 Problem 4.7 This is forced convection flow over a streamlined body. Viscous (velocity) boundary layer approximations can be made if the Reynolds number Re
More informationBERNOULLI EQUATION. The motion of a fluid is usually extremely complex.
BERNOULLI EQUATION The motion of a fluid is usually extremely complex. The study of a fluid at rest, or in relative equilibrium, was simplified by the absence of shear stress, but when a fluid flows over
More informationPROBLEM SET 6. SOLUTIONS April 1, 2004
Harvard-MIT Division of Health Sciences and Technology HST.54J: Quantitative Physiology: Organ Transport Systems Instructors: Roger Mark and Jose Venegas MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departments
More informationContents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1
Contents 1. Introduction 2. Fluids 3. Physics of Microfluidic Systems 4. Microfabrication Technologies 5. Flow Control 6. Micropumps 7. Sensors 8. Ink-Jet Technology 9. Liquid Handling 10.Microarrays 11.Microreactors
More informationExam in Fluid Mechanics 5C1214
Eam in Fluid Mechanics 5C1214 Final eam in course 5C1214 13/01 2004 09-13 in Q24 Eaminer: Prof. Dan Henningson The point value of each question is given in parenthesis and you need more than 20 points
More informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationFORMULA 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 informationFluid dynamics is the study of flows of fluid. However, before we begin discussing the dynamics of
Hydrodynamics of Newtonian Fluids Cameron Stewart Universität Stuttgart cameron.n.stewart@gmail.com April 5, 2017 Abstract In this review we provide an overview of the hydrodynamics of viscous, Newtonian
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 informationBoundary 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 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 informationFLUID MECHANICS AND HEAT TRANSFER
AN INTRODUCTION TO FLUID MECHANICS AND HEAT TRANSFER AN INTRODUCTION TO FLUID MECHANICS AND HEAT TRANSFER WITH APPLICATIONS IN CHEMICAL & MECHANICAL PROCESS ENGINEERING BY J. M. KAY AND R. M. NEDDERMAN
More informationThe role of micro-scale inertia in transport processes. - Ganesh Subramanian (JNCASR, Bangalore)
The role of micro-scale inertia in transport processes - Ganesh Subramanian (JNCASR, Bangalore) Micro-scale inertia (suspensions and emulsions) Navier-Stokes equations : Re u t + Re u u = Inertial acceleration
More informationFor more info
Characteristic of Ideal fluid:- (a) It is incompressible (b) It is non-viscous (c) Flow of ideal fluid is irrational (d) It is capable of exhibiting steady flow Stream line flow:- Flow of a liquid fluid
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar
More informationConvective 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 informationFree 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 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 information2 Law of conservation of energy
1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion
More informationFluid mechanics and living organisms
Physics of the Human Body 37 Chapter 4: In this chapter we discuss the basic laws of fluid flow as they apply to life processes at various size scales For example, fluid dynamics at low Reynolds number
More informationMAE 101A. Homework 7 - Solutions 3/12/2018
MAE 101A Homework 7 - Solutions 3/12/2018 Munson 6.31: The stream function for a two-dimensional, nonviscous, incompressible flow field is given by the expression ψ = 2(x y) where the stream function has
More informationRate of Flow Quantity of fluid passing through any section (area) per unit time
Kinematics of Fluid Flow Kinematics is the science which deals with study of motion of liquids without considering the forces causing the motion. Rate of Flow Quantity of fluid passing through any section
More informationBACHELOR OF TECHNOLOGY IN MECHANICAL ENGINEERING (COMPUTER INTEGRATED MANUFACTURING)
No. of Printed Pages : 6 BME-028 BACHELOR OF TECHNOLOGY IN MECHANICAL ENGINEERING (COMPUTER INTEGRATED MANUFACTURING) Term-End Examination December, 2011 00792 BME-028 : FLUID MECHANICS Time : 3 hours
More informationThe Kelvin transformation in potential theory and Stokes flow
IMA Journal of Applied Mathematics (2009 74, 427 438 doi:10.1093/imamat/hxn027 Advance Access publication on September 9, 2008 The Kelvin transformation in potential theory and Stokes flow GEORGE DASSIOS
More information