Introduction to Fluid Mechanics*

Fred Stern, Tao Xing, Jun Shao, Surajeet Ghosh

8/26/2005

AFD EFD CFD

ѷ U=0
DU 1
=-�p+�+ѷ 2
U ui u j
Dt Re

57:020 Fluid Mecha 1

 Fluid Mechanics
• Fluids essential to life
• Human body 65% water
• Earth’s surface is 2/3 water
• Atmosphere extends 17km above the earth’s surface
• History shaped by fluid mechanics
• Geomorphology
• Human migration and civilization
• Modern scientific and mathematical theories and methods
• Warfare
• Affects every part of our lives

57:020 Fluid Mecha 2

 History
Faces of Fluid Mechanics

Archimedes Newton Leibniz Bernoulli Euler

(C. 287-212 BC) (1642-1727) (1646-1716) (1667-1748) (1707-1783)

Navier Stokes Reynolds Prandtl Taylor

(1785-1836) (1819-1903) (1842-1912) (1875-1953) (1886-1975)
57:020 Fluid Mecha 3
 Significance
• Fluids omnipresent
• Weather & climate
• Vehicles: automobiles, trains, ships, and
planes, etc.
• Environment
• Physiology and medicine
• Sports & recreation
• Many other examples!

57:020 Fluid Mecha 4

 Weather & Climate
Tornadoes Thunderstorm

Global Climate Hurricanes

57:020 Fluid Mecha 5

 Vehicles
Aircraft Surface ships

High-speed rail Submarines

57:020 Fluid Mecha 6

 Environment

Air pollution River hydraulics

57:020 Fluid Mecha 7

 Physiology and Medicine

Blood pump Ventricular assist device

57:020 Fluid Mecha 8

 Sports & Recreation

Water sports Cycling Offshore racing

Auto racing Surfing

57:020 Fluid Mecha 9

 Fluids Engineering

Reality

Fluids Engineering System Components Idealized

EFD, U D = B2 + P2 Mathematical Physics Problem Formulation

AFD, U m CFD, U s = U SM + U SN
2 2

57:020 Fluid Mecha 10

 Analytical Fluid Dynamics
• The theory of mathematical physics
problem formulation
• Control volume & differential analysis
• Exact solutions only exist for simple
geometry and conditions
• Approximate solutions for practical
applications
• Linear
• Empirical relations using EFD data

57:020 Fluid Mecha 11

 Analytical Fluid Dynamics
• Lecture Part of Fluid Class
• Definition and fluids properties
• Fluid statics
• Fluids in motion
• Continuity, momentum, and energy principles
• Dimensional analysis and similitude
• Surface resistance
• Flow in conduits
• Drag and lift

57:020 Fluid Mecha 12

 Analytical Fluid Dynamics
• Example: laminar pipe flow
rUD
Assumptions: Fully developed, Low Re = < 2000
m
Approach: Simplify momentum equation,
integrate, apply boundary conditions to Schematic
determine integration constants and use
energy equation to calculate head loss
0
Du 0 p   u  2u 
2 0
=- + m  2 + 2  + gx
Dt x  x y 
Exact solution :
u(r) = 1 (- �
p)(R2 - r 2)
4m � x
8m du
f = 8t w = dy w = 64
Friction factor:
rV 2 rV 2 Re
p1 p2 L V 2 32 m LV
Head loss: g + z1 = g + z2 + h f hf = f =
D 2g g D2
57:020 Fluid Mecha 13
 Analytical Fluid Dynamics
• Example: turbulent flow in smooth pipe( Re > 3000)
Three layer concept (using dimensional analysis)
+
u =u u *
y = yu n
+ * u* = t w r
1. Laminar sub-layer (viscous shear dominates)
u+ = y+ 0 < y+ < 5
2. Overlap layer (viscous and turbulent shear important)
1
u+ = ln y + + B 20 < y + < 105 (k=0.41, B=5.5)
k
U -u � r� +
3. Outer layer (turbulent shear dominates) *
= f 1 - �y > 10
�
5

u � r0 �
u( r) 1 ( r - r) u
*
Assume log-law is valid across entire pipe: = ln 0 +B
u* k n

Integration for average velocity and using EFD data to adjust constants:
1
= 2log ( Re f 1 2 ) - .8
f
57:020 Fluid Mecha 14
 Analytical Fluid Dynamics
• Example: turbulent flow in rough pipe
Both laminar sublayer and overlap
layer
Inner
are affected +
= u+ ( y k )
layer:by uroughness

Outer layer: unaffected

1 y
Overlap layer: u+ = ln + constant
k k
Three regimes of flow depending on k+
1. K+<5, hydraulically smooth (no effect of roughness)
2. 5 < K+< 70, transitional roughness (Re dependent)
3. K+> 70, fully rough (independent Re)
For 3, using EFD data to adjust constants:
1 y 1 k D
u+ = ln + 8.5 �f ( Re ) Friction factor: = -2log
k k f 3.7

57:020 Fluid Mecha 15

 Analytical Fluid Dynamics
• Example: Moody diagram for turbulent pipe flow
Composite Log-Law for smooth and rough pipes is given by the Moody diagram:

1 �k D 2.51 �
1
= -2log � + 12�
f 2 �3.7 Re f �

57:020 Fluid Mecha 16

Experimental Fluid Dynamics (EFD)

Definition:
Use of experimental methodology and procedures for solving fluids engineering
systems, including full and model scales, large and table top facilities,
measurement systems (instrumentation, data acquisition and data reduction),
uncertainty analysis, and dimensional analysis and similarity.

EFD philosophy:
• Decisions on conducting experiments are governed by the ability of the expected
test outcome, to achieve the test objectives within allowable uncertainties.
• Integration of UA into all test phases should be a key part of entire experimental
program
• test design
• determination of error sources
• estimation of uncertainty
• documentation of the results

57:020 Fluid Mecha 17

 Purpose

• Science & Technology: understand and investigate a

phenomenon/process, substantiate and validate a theory
(hypothesis)
• Research & Development: document a process/system,
provide benchmark data (standard procedures,
validations), calibrate instruments, equipment, and
facilities
• Industry: design optimization and analysis, provide data
for direct use, product liability, and acceptance
• Teaching: instruction/demonstration

57:020 Fluid Mecha 18

 Applications of EFD

Application in science & technology Application in research & development

Picture of Karman vortex shedding Tropic Wind Tunnel has the ability to create
temperatures ranging from 0 to 165 degrees
Fahrenheit and simulate rain

57:020 Fluid Mecha 19

 Applications of EFD (cont’d)

Example of industrial application

NASA's cryogenic wind tunnel simulates flight

conditions for scale models--a critical tool in
designing airplanes.

Application in teaching

Fluid dynamics laboratory

57:020 Fluid Mecha 20
 Full and model scale

• Scales: model, and full-scale

• Selection of the model scale: governed by dimensional analysis and similarity

57:020 Fluid Mecha 21

 Measurement systems
• Instrumentation
• Load cell to measure forces and moments
• Pressure transducers
• Pitot tubes
• Hotwire anemometry
• PIV, LDV
• Data acquisition
• Serial port devices
• Desktop PC’s
• Plug-in data acquisition boards
• Data Acquisition software - Labview
• Data analysis and data reduction
• Data reduction equations
• Spectral analysis

57:020 Fluid Mecha 22

 Instrumentation

Pitot tube

Load cell

Hotwire 3D - PIV
57:020 Fluid Mecha 23
 Data acquisition system

Hardware

Software - Labview

57:020 Fluid Mecha 24

 Data reduction methods
TEM PERATURE
W ATER
TEM PERATURE
A IR
P IP E
PRESSURE
VENTURI
PRESSURE
IN D IV ID U
M EASUREM
SYSTEM

• Data reduction equations

T T z z M EASUREM
w a SM DM
B , Pz B , Pz O F IN D IV ID
B , PT B , PT
• Spectral analysisTw
w
T a
a
zSM SM zDM DM V A R IA B L

r = F (T )
w w
ra = F (T a )
D ATA R E D U C
Q = F ( z DM
) E Q U A T IO
2
g   D
5 rw
f = F (r , r , z
w a SM
,Q ) = 2
(z - z
ra SM i SM j
)
8LQ

 2  gr w
u (r ) =   
 z SM Stag ( r ) - z SM Stat  
 ra 
f E X P E R IM E N
B f , P fequations
Example of data reduction R E S U LT

57:020 Fluid Mecha 25

 Spectral analysis
Aim: To analyze the natural FFT: Converts a function from amplitude as function
unsteadiness of the separated flow, of time to amplitude as function of frequency
around a surface piercing
strut, using FFT.

Fast Fourier Transform

Free-surface wave Time history of wave

elevation contours elevation

Surface piercing strut FFT of wave elevation Power spectral density

57:020 Fluid Mecha 26 of wave elevation
 Uncertainty analysis

Rigorous methodology for uncertainty assessment

using statistical and engineering concepts

E LE M EN TA L
ERRO R SO UR CES

IN D IV ID U A L
1 2 J M EASUREM ENT
SYSTEM S

X X X M EASUREM ENT
1 2 J O F IN D IV ID U A L
B , P B , P B , P V A R IA B L E S
1 1 2 2 J J

D A T A R E D U C T IO N
r = r ( X , X ,......, X ) E Q U A T IO N
1 2 J

r E X P E R IM E N T A L
B , P R E S U LT
r r

57:020 Fluid Mecha 27

 Dimensional analysis
• Definition : Dimensional analysis is a process of formulating fluid mechanics problems in
in terms of non-dimensional variables and parameters.
• Why is it used :
• Reduction in variables ( If F(A1, A2, … , An) = 0, then f(1, 2, … r < n) = 0,
where, F = functional form, Ai = dimensional variables, j = non-dimensional
parameters, m = number of important dimensions, n = number of dimensional variables, r
= n – m ). Thereby the number of experiments required to determine f vs. F is reduced.
• Helps in understanding physics
• Useful in data analysis and modeling
• Enables scaling of different physical dimensions and fluid properties
Example Drag = f(V, L, r, m, c, t, e, T, etc.)
From dimensional analysis,

Examples of dimensionless quantities : Reynolds number, Froude

Vortex shedding behind cylinder
Number, Strouhal number, Euler number, etc.
57:020 Fluid Mecha 28
 Similarity and model testing
• Definition : Flow conditions for a model test are completely similar if all relevant
dimensionless parameters have the same corresponding values for model and prototype.
• i model = i prototype i = 1
• Enables extrapolation from model to full scale
• However, complete similarity usually not possible. Therefore, often it is necessary to
use Re, or Fr, or Ma scaling, i.e., select most important and accommodate others
as best possible.
• Types of similarity:
• Geometric Similarity : all body dimensions in all three coordinates have the same
linear-scale ratios.
• Kinematic Similarity : homologous (same relative position) particles lie at homologous
points at homologous times.
• Dynamic Similarity : in addition to the requirements for kinematic similarity the model
and prototype forces must be in a constant ratio.

57:020 Fluid Mecha 29

 EFD process
• “EFD process” is the steps to set up an experiment and
take data

57:020 Fluid Mecha 30

 EFD – “hands on” experience

Lab1: Measurement of density and

kinematic viscosity of a fluid
C h o rd - w is e
P re s s u re

Tygon
Tap s
L L o a d C e ll
Lab2: Measurement of
flow rate, friction factor and
T u b in g D

To

velocity profiles in smooth and

L o a d C e ll
S c a n iv a lv e

rough pipes.

Lab3: Measurement of surface pressure

Distribution, lift and drag coefficient for an airfoil
57:020 Fluid Mecha 31
 Computational Fluid Dynamics
• CFD is use of computational methods for
solving fluid engineering systems, including
modeling (mathematical & Physics) and
numerical methods (solvers, finite differences,
and grid generations, etc.).
• Rapid growth in CFD technology since advent
of computer

ENIAC 1, 1946 IBM WorkStation

57:020 Fluid Mecha 32

 Purpose
• The objective of CFD is to model the continuous fluids
with Partial Differential Equations (PDEs) and
discretize PDEs into an algebra problem, solve it,
validate it and achieve simulation based design
instead of “build & test”

• Simulation of physical fluid phenomena that are

difficult to be measured by experiments: scale
simulations (full-scale ships, airplanes), hazards
(explosions,radiations,pollution), physics (weather
prediction, planetary boundary layer, stellar
evolution).

57:020 Fluid Mecha 33

 Modeling
• Mathematical physics problem formulation of fluid
engineering system
• Governing equations: Navier-Stokes equations (momentum),
continuity equation, pressure Poisson equation, energy
equation, ideal gas law, combustions (chemical reaction
equation), multi-phase flows(e.g. Rayleigh equation), and
turbulent models (RANS, LES, DES).
• Coordinates: Cartesian, cylindrical and spherical coordinates
result in different form of governing equations
• Initial conditions(initial guess of the solution) and Boundary
Conditions (no-slip wall, free-surface, zero-gradient,
symmetry, velocity/pressure inlet/outlet)
• Flow conditions: Geometry approximation, domain, Reynolds
Number, and Mach Number, etc.

57:020 Fluid Mecha 34

 Modeling (examples)
Developing flame surface (Bell et al., 2001)
Free surface animation for ship in
regular waves

Evolution of a 2D mixing layer laden with particles of Stokes

Number 0.3 with respect to the vortex time scale (C.Narayanan)

57:020 Fluid Mecha 35

 Modeling (examples, cont’d)

3D vortex shedding behind a circular cylinder (Re=100,DNS,J.Dijkstra)

DES,
Re=105, Iso-
surface of Q
criterion (0.4)
for turbulent
flow around
NACA12 with
angle of attack
60 degrees

LES of a turbulent jet. Back wall shows a slice of the dissipation rate and the
bottom wall shows a carpet plot of the mixture fraction in a slice through the jet
centerline, Re=21,000 (D. Glaze).
57:020 Fluid Mecha 36
 Numerical methods
• Finite difference methods: using y
x
numerical scheme to approximate jmax
the exact derivatives in the PDEs
j+1
�2
P Pi +1 - 2 Pi + Pi -1
= j
�x2 x 2 y
j-1
�P2
Pj +1 - 2 Pj + Pj -1
=
�y2 y 2

• Finite volume methods

o i-1 i i+1 imax x
• Grid generation: conformal
mapping, algebraic methods and
differential equation methods
• Grid types: structured,
unstructured
• Solvers: direct methods (Cramer’s
rule, Gauss elimination, LU
decomposition) and iterative
methods (Jacobi, Gauss-Seidel,
SOR)
Slice of 3D mesh of a fighter aircraft
57:020 Fluid Mecha 37
 CFD process
Geometry Physics Mesh Solve Reports Post-
Processing

Select Heat Transfer Unstructured Steady/ Forces Report Contours

Geometry ON/OFF (automatic/ Unsteady (lift/drag, shear
manual) stress, etc)

Compressible Structured Iterations/ XY Plot Vectors

Geometry ON/OFF (automatic/ Steps
Parameters manual)

Domain Shape Flow Convergent Verification Streamlines

and Size properties Limit

Viscous Model Precisions Validation

(single/
double)

Boundary Numerical
Conditions Scheme

Initial
Conditions

57:020 Fluid Mecha 38

 Commercial software
• CFD software
1. FLUENT: http://www.fluent.com
2. FLOWLAB: http://www.flowlab.fluent.com
3. CFDRC: http://www.cfdrc.com
4. STAR-CD: http://www.cd-adapco.com
5. CFX/AEA: http://www.software.aeat.com/cfx
• Grid Generation software
1. Gridgen: http://www.pointwise.com
2. GridPro: http://www.gridpro.com
• Visualization software
1. Tecplot: http://www.amtec.com
2. Fieldview: http://www.ilight.com

57:020 Fluid Mecha 39

 “Hands-on” experience using CFD
Educational Interface (pipe template)

57:020 Fluid Mecha 40

 “Hands-on” experience using CFD
Educational Interface (airfoil template)

57:020 Fluid Mecha 41

 57:020 Fluid Mechanics
• Lectures cover basic concepts in fluid statics,
kinematics, and dynamics, control-volume, and
differential-equation analysis methods. Homework
assignments, tests, and complementary EFD/CFD
labs
• This class provides an introduction to all three tools:
AFD through lecture and CFD and EFD through labs
• ISTUE Teaching Modules
(http://www.iihr.uiowa.edu/~istue) (next two slides)

57:020 Fluid Mecha 42

 TM Descriptions
Table 1: ISTUE Teaching Modules for Introductory Level Fluid Mechanics at Iowa
Teaching Modules TM for Fluid Property TM for Pipe Flow TM for Airfoil Flow

Overall Purpose Hands-on student Hands-on student experience Hands-on student experience with
experience with table-top with complementary EFD, CFD, complementary EFD, CFD, and UA
facility and simple MS for and UA for Introductory Pipe for Introductory Airfoil Flow,
fluid property Flow, including friction factor including lift and drag, surface
measurement, including and mean velocity measurements pressure, and mean and turbulent
comparison manufacturer and comparisons benchmark wake velocity profile measurements
values and rigorous data, laminar and turbulent flow and comparisons benchmark data,
implementation standard CFD simulations, modeling and inviscid and turbulent flow
EFD UA verification studies, and simulations, modeling and verification
validation using AFD and EFD. studies, and validation using AFD and
EFD.

Educational Materials FM and EFD lecture; lab FM, EFD and CFD lectures; lab FM, EFD and CFD lectures; lab
report instructions; pre lab report instructions; pre lab report instructions; pre lab questions,
questions, and EFD questions, and EFD and CFD and EFD and CFD exercise notes.
exercise notes. exercise notes.

ISTUE ASEE papers Paper 1 Paper2

FM Lecture Introduction to Fluid Mechanics

Lab Report Instructions EFD lab report Instructions CFD lab report Instructions

Continued in next slide…

http://css.engineering.uiowa.edu/~fluids
57:020 Fluid Mecha 43
 TM Descriptions, cont’d
Teaching Modules TM for Fluid TM for Pipe Flow TM for Airfoil Flow
Property
CFD CFD Lecture Introduction to CFD

Exercise Notes None CFD Prelab1 CFD Prelab2

PreLab1 Questions PreLab 2 Questions
CFD Lab 1 CFD Lab2
Lab Concepts Lab Concepts

EFD Data EFD Data

EFD EFD EFD and UA

Lecture
Exercise Notes PreLab1 Questions PreLab2 Questions PreLab3 Questions
Lab1 Lecture Lab2 Lecture Lab3 Lecture
Lab 1 exercise notes Lab 2 exercise notes Lab 3 exercise notes
Lab 1 data Lab2 Data reduction Data Reduction Sheet
reduction sheet sheet (smooth & rough) Lab concepts
Lab concepts Lab concepts

UA(EFD) Lecture EFD UA Report

Exercise Notes Instructions UA

UA(CFD) Lecture

Exercise Notes None

57:020 Fluid Mecha 44