Lecture 1.1 Introduction to Fluid Dynamics

Similar documents
FLUID MECHANICS TESTS

FLOWING FLUIDS AND PRESSURE VARIATION

Chapter 1 - Basic Equations

1 Mathematical Concepts

MAE 3130: Fluid Mechanics Lecture 5: Fluid Kinematics Spring Dr. Jason Roney Mechanical and Aerospace Engineering

Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion

Inviscid Flows. Introduction. T. J. Craft George Begg Building, C41. The Euler Equations. 3rd Year Fluid Mechanics

Chapter 5: Introduction to Differential Analysis of Fluid Motion

CHAPTER 3. Elementary Fluid Dynamics

Module 4: Fluid Dynamics Lecture 9: Lagrangian and Eulerian approaches; Euler's acceleration formula. Fluid Dynamics: description of fluid-motion

Chapter 1 - Basic Equations

2.29 Numerical Marine Hydrodynamics Spring 2007

Essay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions

Grad operator, triple and line integrals. Notice: this material must not be used as a substitute for attending the lectures

Table of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Contents.

Realistic Animation of Fluids

Driven Cavity Example

1. Suppose that the equation F (x, y, z) = 0 implicitly defines each of the three variables x, y, and z as functions of the other two:

Chapter 4 FLUID KINEMATICS

University of California, Berkeley

University of Saskatchewan Department of Mathematics & Statistics MATH Final Instructors: (01) P. J. Browne (03) B. Friberg (05) H.

Name: Class: Date: 1. Use Lagrange multipliers to find the maximum and minimum values of the function subject to the given constraint.

Data Visualization. Fall 2017

MAC2313 Final A. a. The vector r u r v lies in the tangent plane of S at a given point. b. S f(x, y, z) ds = R f(r(u, v)) r u r v du dv.

Vector Visualization. CSC 7443: Scientific Information Visualization

1 Vector Functions and Space Curves

BASICS OF FLUID MECHANICS AND INTRODUCTION TO COMPUTATIONAL FLUID DYNAMICS

Solution for Euler Equations Lagrangian and Eulerian Descriptions

Solution for Euler Equations Lagrangian and Eulerian Descriptions

Solution for Euler Equations Lagrangian and Eulerian Descriptions

Example 13 - Shock Tube

MA FINAL EXAM Green April 30, 2018 EXAM POLICIES

The Divergence Theorem

Solution of final examination

Vector Visualisation 1. global view

Math 5BI: Problem Set 2 The Chain Rule

Chapter 5 Partial Differentiation

18.02 Final Exam. y = 0

This exam will be cumulative. Consult the review sheets for the midterms for reviews of Chapters

Numerical Methods in Aerodynamics. Fluid Structure Interaction. Lecture 4: Fluid Structure Interaction

MA 174: Multivariable Calculus Final EXAM (practice) NO CALCULATORS, BOOKS, OR PAPERS ARE ALLOWED. Use the back of the test pages for scrap paper.

A simple method for potential flow simulation of cascades

PURE MATHEMATICS 212 Multivariable Calculus CONTENTS. Page. 1. Assignment Summary... i 2. Summary Assignments...2

Debojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering

Directional Derivatives and the Gradient Vector Part 2

Particle Velocimetry Data from COMSOL Model of Micro-channels

MATH 2023 Multivariable Calculus

Wall-bounded laminar sink flows

minutes/question 26 minutes

Topic 5.1: Line Elements and Scalar Line Integrals. Textbook: Section 16.2

Module 1 Lecture Notes 2. Optimization Problem and Model Formulation

MA 243 Calculus III Fall Assignment 1. Reading assignments are found in James Stewart s Calculus (Early Transcendentals)

Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods

Development of the Compliant Mooring Line Model for FLOW-3D

SPC 307 Aerodynamics. Lecture 1. February 10, 2018

Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data

Overview. Applications of DEC: Fluid Mechanics and Meshing. Fluid Models (I) Part I. Computational Fluids with DEC. Fluid Models (II) Fluid Models (I)

Flow Visualisation - Background. CITS4241 Visualisation Lectures 20 and 21

MATH 2400, Analytic Geometry and Calculus 3

CGT 581 G Fluids. Overview. Some terms. Some terms

MAT175 Overview and Sample Problems

Vector Field Visualisation

Let and be a differentiable function. Let Then be the level surface given by

Vector Visualization

Section 17.7: Surface Integrals. 1 Objectives. 2 Assignments. 3 Maple Commands. 4 Lecture. 4.1 Riemann definition

Chapter 15 Vector Calculus

2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123

ANALYSIS OF THE MODIFIED OSEEN EQUATIONS USING ORTHOGONAL FUNCTION EXPANSIONS

Axisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows

Faculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia

21. Efficient and fast numerical methods to compute fluid flows in the geophysical β plane

The directional derivative of f x, y in the direction of at x, y u. f x sa y sb f x y (, ) (, ) 0 0 y 0 0

Verification of Moving Mesh Discretizations

Chapter 3 Numerical Methods

Hardware Accelerated Real-Time Fluid Surface Simulation

Surfaces and Integral Curves

CS-184: Computer Graphics Lecture #21: Fluid Simulation II

Exam 2 Preparation Math 2080 (Spring 2011) Exam 2: Thursday, May 12.

Volumes of Solids of Revolution Lecture #6 a

Background for Surface Integration

MATH 234. Excercises on Integration in Several Variables. I. Double Integrals

F dr = f dx + g dy + h dz. Using that dz = q x dx + q y dy we get. (g + hq y ) x (f + hq x ) y da.

PURE MATHEMATICS 212 Multivariable Calculus CONTENTS. Page. 1. Assignment Summary... i 2. Introduction Timetable Assignments...

COPYRIGHTED MATERIAL INTRODUCTION CHAPTER 1

FLUID MECHANICAL DESIGN OF THE POTENTIAL FLOW TURBINE. This paper presents the software package developed at Fluid and Power Research

MESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP

Modeling icing using cartesian grids, penalization & level sets

Math 2130 Practice Problems Sec Name. Change the Cartesian integral to an equivalent polar integral, and then evaluate.

the lines of the solution obtained in for the twodimensional for an incompressible secondorder

Dr. Allen Back. Nov. 21, 2014

DIGITAL THESES PROJECT. I authorise the University to make this thesis available on the Internet under the International digital theses project

Directional Derivatives and the Gradient Vector Part 2

Level Set Methods and Fast Marching Methods

Lecture 2.2 Cubic Splines

A brief introduction to fluidstructure. O. Souček

MAC2313 Test 3 A E g(x, y, z) dy dx dz

Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent

Free-surface flow under a sluice gate from deep water

Characteristic Aspects of SPH Solutions

Vector Visualization Chap. 6 March 7, 2013 March 26, Jie Zhang Copyright

Transcription:

Lecture 1.1 Introduction to Fluid Dynamics 1

Introduction A thorough study of the laws of fluid mechanics is necessary to understand the fluid motion within the turbomachinery components. In this introductory lecture, some relevant concepts of fluid mechanics are revised, before entering the realms of computational fluid dynamics. Fluid mechanics has two parts: Kinematics Dynamics 2

The kinematics is about describing the fluid motion without taking into account the forces that cause the motion. In fluid dynamics, the forces are considered for the analysis of fluid motion. Governing equations are derived by considering the balance of these forces. Computational fluid dynamics is a subject where the fluid motion is studied by solving the governing differential equations using numerical methods. 3

Microscopic versus Macroscopic Approaches Microscopic and macroscopic are the two approaches using which the fluid flows are analyzed and the flow properties like density, velocity, etc. are determined at each point of the flow domain. In the microscopic approach, the molecular mean free path is far larger than the characteristic dimension. In the macroscopic approach, the fluid is considered as a continuum, that is, the fluid in a region is assumed to occupy every geometric point of that region. The point here means a very small volume in the fluid region whose size is zero in the limit. 4

Flow studies in machines working at very low pressures such as a turbo molecular pump require microscopic view. The flows in turbomachines, which are of interest in this course, are all in the continuum and hence the macroscopic approach is appropriate. The course objectives are, therefore, to know The formulation of equations governing the fluid flow in turbomachines Physical modeling of the flow domains in a turbomachine The methodology to solve the governing equations and The analysis of the numerical results. The results obtained for a turbomachine model can be extended to another geometrically similar prototype machine when the kinematic and dynamic similarity rules are respected between the model and the prototype. 5

Eulerian and Lagrangian Approaches In the Eulerian approach, the flow properties are described as functions of space and time, like in cinematography. The complete state of motion is described by a succession of instantaneous states of flow. Thus if is a flow property, then (x,y,z,t)or ( rt, ) is the value of at r, the position occupied by the particle and at the instant of time t. At a later time the fluid particle occupying the position (x, y, z) will be different. If the flow is steady then is independent of t and is a function of x, y and z alone. 6

In the Lagrange method, fluid particles are identified by their position r (or x o, y o, z o,) at any time t o and describe their o trajectories. The trajectory of the line along which a particle moves is called a path line. The instantaneous velocity and acceleration v and can be expressed as v v r t a a ro t 2 r v r o 2 t r const. t r const t r const, and, o o o Except for multi-phase flows such as gas-solid flows, where one or both of the phases may be considered for the Lagrangian description, the Eulerian approach is more appropriate for describing the fluid motion. The governing equations in this course are therefore written using only the Eulerian approach. a 7

Velocity of a Fluid Particle In the Eulerian description, the particles which occupy a chosen point are going to change from time t to t + δt. Let P be the point whose position vector is r at time t with respect to some fixed point O. That is, OP r After a time δt the same particle is moved to P, as shown in the Fig. 1.1.1, having a distance δr from P. Therefore, the particle velocity v at P is given by r dr v lim (1.1.1) t 0 t dt Figure 1.1.1 Velocity 8

That is, the velocity is dependent on r and t. The velocity in Cartesian coordinates is given by v uiˆvj ˆwkˆ (1.1.2) where dx dy dz u, v, w dt dt dt iˆ, ˆjandkˆ are the unit vectors in the coordinate directions x, y, z, respectively. Note that r xiˆ yj ˆzkˆ (1.1.3) 9

Substantial Derivative In order to compute the time rate of change of any property at the chosen point, therefore, the local rate of change of the property and also the change of its position have to be used. That is D u v w v Dt t x y z t (1.1.4) The derivative D/Dt is called by several names such as substantial derivative, material derivative, particle derivative or total derivative. 10

Acceleration of a Fluid Particle The first term on R.H.S. of Eq. (1.1.4), is called the local derivative which indicates the unsteady time variation of fluid property at a point. The sum of the last three terms is called convective derivative of the particle. Using the substantial derivative, the acceleration of a fluid particle is written as Dv v v v v v (1.1.5) a v v u v w Dt t t x y z Note that the acceleration Dv Dt is called total acceleration, comprising of local acceleration (partial derivative w.r.t. time) and convective acceleration (caused due to the fluid motion itself). (x, y, z) is the inertial frame of reference. 11

Streamline Streamline is an imaginary curve in the flow such that the tangent at each point on the curve coincides with the direction of the velocity at that point. When the flow is steady, the streamline coincides with the path line of any fluid particle. The equation of a streamline is dx dy dz u v w (1.1.6) However, in unsteady motion, the flow pattern varies with time. Therefore, the streamlines and path lines differ with each other. (refer Fig. 1.1.2). The turbomachinery flows, in general, are unsteady. 12

Fig. 1.1.2 Streamlines and pathlines 13

Stream Surface and Stream Tube A stream sheet or stream surface is generated by a group of streamlines passing through some given curve, C. If the curve C is closed, the streamlines form a stream tube. (Since a stream tube is like a virtual solid surface, a particle can have no motion perpendicular to the streamline.) The space that the fluid occupies in a turbomachine can be considered as a stream tube. Even the flow between the blades of the turbomachines may be treated approximately in this manner. This approach gives an average behavior of fluid flow and is known as a quasi-one dimensional treatment. 14

Stream Function As there can be no flow across a streamline, stream function is a constant over a given streamline. From the equation of streamline Eq. (1.1.6), the velocity components can be readily derived as derivatives of stream function with respect to the spatial coordinates. For incompressible flow, the change in stream function signifies the change in mass flow rate. For inviscid and incompressible flow, (also irrotational flows) the stream function satisfies the Laplace s equation. A flow net, with the mutually orthogonal stream function and velocity potential functions, can be constructed for such flows. 15

Vorticity and Irrotational Flow The curl of the velocity vector is called vorticity, v If i j k, then, w v u w v u ; ; and y z z x x y (1.1.7) Vorticity is a vector and its magnitude is twice the angular velocity of the fluid particle at a given position in the flow field. 16

A flow field is called irrotational when the vorticity is zero i.e. = = =0. Velocity potential exists in such a flow and hence is called potential flow. The velocity components in this case are given by u ; v ; w x y z (1.1.8) 17

Vortex Line, Vortex Sheet and Vortex Tube Vortex line is an imaginary curve lying in the flow field such that its tangent at any point coincides with the direction of vorticity at that point. The equation of a vortex line is dx dy dz Vortex sheet is a surface composed of vortex lines. (1.1.9) Vortex tube is generated by vortex lines, drawn through each point of a closed curve C. A vortex tube is also referred simply as a vortex. A vortex tube of a very small cross-section is called a vortex filament. 18

Circulation The instantaneous line integral of the tangential velocity around any closed curve C in a flow field is called circulation, denoted by. Fig.1.1.4 Circulation Fig. 1.1.3 Circulation It is taken as positive when C, as shown in Fig. 1.1.3, is traversed such that the area enclosed by C lies to the left. v ds udxvdy wdz C A C v nda (by Stokes theorem) (1.1.10) Stokes theorem links circulation with vorticity. 19

The circulation around any curve C, bounding an area A (simply or multiply connected), is the sum of the circulations around all the lesser areas in to which A may be arbitrarily divided. Strength of a vortex tube is the circulation along a circuit lying on the surface of the vortex tube and passing around it only once. Strength of a vortex is same throughout its length. A vortex cannot have an end within the fluid i.e. vortex filament either forms a closed curve or extends to the fluid boundaries. The above points can be applied to develop a simple theory called isolated blade element theory for turbomachines, where the number of blades is not very large. 20

Summary of Lecture 1.1 In this lecture introductory aspects of fluid dynamics such as Microscopic versus Macroscopic approaches, Eulerian versus Lagrangian approaches and basic definitions of velocity stream function, vorticity and circulation are presented. The understanding of the concept of circulation is required to develop blade element theory used for turbomachines. END OF LECTURE 1.1 21