Meccanica dei Fluidi
Alessandro Bottaro
Dipartimento di Ingegneria Civile,
Chimica e Ambientale (DICCA)
Secondo Semestre 2014/2015
Textbook
Fluid Mechanics. Fundamentals and
Applications, McGraw-Hill, 2006
Yunus A. Çengal (UNV Reno) and John M. Cimbala
(Penn State)
Includes DVD with movies created at PSU by Prof.
Gary Settles
Available at
Amazon.com, ~ $90.00 (paperback)
Libreria Frasconi, Corso Gastaldi 193r
A version in Italian exists …
Check the following web site for movies of fluid
motion, simulations, exercises, sample exams, etc.:
http://highered.mcgrawhill.com/sites/0072472367/information_center_view0/
(Online Learning Center, Student Edition)
Meccanica dei Fluidi I
2
Information and Introduction
CD-ROM
MultiMedia Fluid Mechanics,
by G.M. Homsy et al., Cambridge U., 2004
Available at
Amazon.com, $21.00
Can be loaned at the DICCA library
(Villa Cambiaso, Mrs. Tina D’Agostino)
Meccanica dei Fluidi I
3
Information and Introduction
Web site
All class material and announcements will be posted on
the course web site: www.dicat.unige.it/bottaro/fmnew.html
Syllabus
Schedule/Calendar
Lecture notes
Message boards
Past mid terms and finals
Exam rules
Grades
Other interesting links can be found at
www.dicat.unige.it/bottaro/teaching.html
Roberto Verzicco
Jean Pierre Petit
Meccanica dei Fluidi I
4
Information and Introduction
IL VOLO, Jean Pierre Petit
Meccanica dei Fluidi I
5
Information and Introduction
Grading
Mid-term exam:
50%
Strongly encouraged!!!!
Final exam:
50%
For those doing “mid-term + final” the oral exam is
optional (to be done in June/July, 2015)
The grade of the oral exam (required for those with a
grade G with 15  G < 18) averages out with the written
tests. Cut-off grade: 15.
For those not doing a “mid-term+final”:
Comprehensive written exam
compulsory oral exam
Meccanica dei Fluidi I
6
+
Information and Introduction
Dates for the “mid-term + final” exams
mid term: April 3, 2015; 14h00 - 17h00, B2
final:
June 3, 2015; 9h00 - 12h00, B2
Optional oral exam for those under this option
is in June/July 2015.
The grade is registered one week after having being
posted on the instructor’s web site. If you do not accept
the grade send an email to: [email protected]
Meccanica dei Fluidi I
7
Information and Introduction
Future indicative exam periods
July 2015
September 2015
EXACT DATES WILL BE POSTED ON THE
INSTRUCTOR’S WEB SITE
Meccanica dei Fluidi I
8
Information and Introduction
Exam policies
Philosophy
One of the best ways to learn something is
through practice and repetition
Therefore, exercises are extremely important
in this class!
If you study and understand the exercises in
the book and elsewhere, you should not have
to struggle with the exams
Meccanica dei Fluidi I
9
Information and Introduction
BONUS!!
The UNIGE-ME class photo competition
Keep an eye on fluid flow phenomena, and
take pictures!
Send me your best original shots, with
indications of date/location/brief description
(about 100 words) of phenomenon you are
displaying.
The three best photographs will gain 3/2/1
points to be added to your final grade.
Only one entry per student. No group entries.
Meccanica dei Fluidi I
10
Information and Introduction
BONUS!!
The UNIGE-ME class photo competition
All photos will be judged by the instructor (with
the possible help of two colleagues) on three
criteria: aesthetic appeal, uniqueness of the
phenomenon, and quality of explanation of the
observed phenomenon.
All photos will be published in a special section
of the instructor’s web site.
Meccanica dei Fluidi I
11
Information and Introduction
EXAMPLES
Dye Droplets at an Oil-Water Interface
This image shows the portion of a glass filled with water (bottom, higher density) and coconut oil
(top, lower density), and the droplets of food dye that rest on the interface between the oil and the
water. It illustrates the effects of surface tension, both at the oil-water interface and at the surface
of the droplets. The droplets are supported by the surface tension at the oil-water interface (they
are on the oil side of the interface and thus they will not mix with the water just yet). Interfacial
tension at the droplet surfaces means that they take on a spherical shape that minimizes their
surface area. When the droplets diffuse through the interface and enter the water (with which they
are miscible), they burst. Just below the oil interface, the different colors have not diffused into
each other yet, but they have on the bottom of the water layer (as indicated by the darker color).
1st prize MIT photo contest 2014
Meccanica dei Fluidi I
12
Information and Introduction
EXAMPLES
Smoke Ring
Smoke rings are possible through the use of toroidal vortices. A toroidal vortex occurs when a
fast-moving parcel of fluid is injected into a stationary fluid. Different parameters, such as
temperature, relative speed, and size of the moving fluid all affect the “crispness” of a smoke ring.
Normally, a vortex is a parcel of fluid spinning around a linear axis, like a tornado or hurricane. In
a toroidal vortex, the axis is still there, but it loops and closes on itself so that the vortex forms a
donut shape. Thus the spinning air traps the smoke inside the vortex, forming a barrier with the
surrounding, stationary fluid. This spinning flow decreases the friction between this parcel of air
and the stationary air around it. Thus the ring can travel for long distances and remain intact, while
other smoke trails blown out with it dissipate.
2nd prize MIT photo contest 2014
Meccanica dei Fluidi I
13
Information and Introduction
EXAMPLES
1st prize UWA photo contest 2013
Paint on a speaker
2nd prize UWA photo contest 2013
Dye flowing into a syphon
3rd prize UWA photo contest 2013
Water balloon to the face
Meccanica dei Fluidi I
14
Information and Introduction
Motivation for Studying Fluid Mechanics
Fluid Mechanics is present almost everywhere
Aerodynamics
Bioengineering and biological systems
Combustion
Energy generation
Geology
Hydraulics and Hydrology
Hydrodynamics
Meteorology
Ocean and Coastal Engineering
Water Resources
…numerous other examples…
Fluid Mechanics is beautiful
Meccanica dei Fluidi I
15
Information and Introduction
Aerodynamics
Meccanica dei Fluidi I
16
Information and Introduction
Bioengineering
Meccanica dei Fluidi I
17
Information and Introduction
Energy generation
Meccanica dei Fluidi I
18
Information and Introduction
Geology
Meccanica dei Fluidi I
19
Information and Introduction
River Hydraulics
Meccanica dei Fluidi I
20
Information and Introduction
Hydraulic Structures
Meccanica dei Fluidi I
21
Information and Introduction
Hydrodynamics
Meccanica dei Fluidi I
22
Information and Introduction
Meteorology
Meccanica dei Fluidi I
23
Information and Introduction
Water Resources
Meccanica dei Fluidi I
24
Information and Introduction
Fluid Mechanics is Beautiful
Meccanica dei Fluidi I
25
Information and Introduction
Tsunamis
Tsunami: Japanese for “Harbour Wave”
Created by earthquakes, land slides, volcanoes,
asteroids/meteors
Pose infrequent but high risk for coastal regions.
Meccanica dei Fluidi I
26
Information and Introduction
Tsunamis: role in religion, evolution, and
apocalyptic events?
Most cultures have deep at
their core a flood myth in which
the great bulk of humanity is
destroyed and a few are left to
repopulate and repurify the
human race. In most of these
stories, God is meting out
retribution, punishing those
who have strayed from his path
Were these “local” floods due
to a tsunami instead of global
events?
Meccanica dei Fluidi I
27
Information and Introduction
Tsunamis: role in religion, evolution, and
apocalyptic events?
Scientists now widely accept that the
worldwide sequence of mass extinctions at
the Cretaceous-Tertiary (K/T) boundary 65.5
million years ago was directly caused by the
collision of an asteroid or comet with Earth.
Evidence for this, all dated to the same
epoch as the extinction event, includes:
the large (200-km diameter) buried
impact structure at Chicxulub in Mexico's
Yucatan Peninsula,
the worldwide iridium-enriched layer at
the K/T boundary,
and the tsunamic deposits well inland in
North America, all dated to the same
epoch as the extinction event.
Meccanica dei Fluidi I
28
Information and Introduction
Tsunamis: role in religion, evolution, and
apocalyptic events?
La Palma Mega-Tsunami = geologic time bomb?
Cumbre Vieja volcano eruption could cause
western half of La Palma (Canary islands) to
collapse into the Atlantic and send a 100 m
tsunami crashing into Eastern coast of U.S.
Meccanica dei Fluidi I
29
Information and Introduction
Methods for Solving Fluid Dynamics
Problems
Analytical Fluid Dynamics (AFD)
Mathematical analysis of governing
equations, including exact and
approximate solutions. This is the primary
focus of this course
Computational Fluid Dynamics (CFD)
Numerical solution of the governing
equations
Experimental Fluid Dynamics (EFD)
Observation and data acquisition.
Meccanica dei Fluidi I
30
Information and Introduction
Analytical Fluid Dynamics
How fast do tsunamis travel in the deep ocean?
Incompressible Navier-Stokes equations
Linearized wave equation for inviscid, irrotational flow
Shallow-water approximation, l/h >> 1 (also kh << 1)
For g = 9.8 m/s2 and h = 3000 m, c = 171 m/s = 617 km/h
Meccanica dei Fluidi I
31
Information and Introduction
Computational Fluid Dynamics
In comparison to analytical
methods, which are good
for providing solutions for
simple geometries or
behavior for limiting
conditions (such as
linearized shallow water
waves), CFD provides a
tool for solving problems
with nonlinear physics and
complex geometry.
Animation by Vasily V. Titov, Tsunami
Inundation Mapping Efforts, NOAA/PMEL
Meccanica dei Fluidi I
32
Information and Introduction
Experimental Fluid Dynamics
Oregon State University
Wave Research
Laboratory
Model-scale experimental
facilities
Tsunami Wave Basin
Large Wave Flume
Dimensional analysis
(Chapter 7 of C&C) is
very important in
designing a model
experiment which
represents physics of
actual problem
Meccanica dei Fluidi I
33
Information and Introduction
Experimental Fluid Dynamics
Experiments are sometimes
conducted in the field or at full
scale
For tsunamis, data acquisition is
used for warning
DART: Deep-ocean Assessment
and Reporting of Tsunamis
(U.S. National Tsunami Hazard
Mitigation Program)
Primary sensor: Bourdon tube for
measuring hydrostatic pressure
Meccanica dei Fluidi I
34
Information and Introduction
THE FLUID DYNAMICS GROUP AT DICAT
Dipartimento di Ingegneria delle Costruzioni, dell’Ambiente e del Territorio,
Università degli Studi di Genova.
Meccanica dei Fluidi I
35
Information and Introduction
BRIEF OVERVIEW OF RECENT ACTIVITIES
1. FUNDAMENTAL STUDIES ON STABILITY AND TRANSITION
- LINEAR STABILITY ANALYSES IN DUCTS, BOUNDARY
LAYERS, FLOWS PAST COMPLIANT WALLS, ETC.
- DIRECT NUMERICAL SIMULATIONS
2. SIMULATION OF 2D AND 3D COMPLEX MOVING GEOMETRIES
- CONVENTIONAL AND CHIMERA GRID APPROACHES
WITH AMR
- FLAPPING FOILS AND WINGS, OPTIMISATION
- BEATING CILIA (ACTIVELY AND PASSIVELY)
3. DETONATION AND SHOCK WAVES
NOTE:
ALL CODES DEVELOPED AND SHOWN HERE ARE HOME-MADE
Meccanica dei Fluidi I
36
Information and Introduction
1.
FUNDAMENTAL STUDIES ON
STABILITY AND TRANSITION
Meccanica dei Fluidi I
37
Information and Introduction
EFFECT OF WALL COMPLIANCE
Study of the effect of compliant walls on the onset of hydrodynamic
instability modes, and on the creation of hydroelastic modes
Wall: shell (or plate) theory
Disturbance field behaves like exp(iax + ibz + st)
Eigensolver: QZ-global or Arnoldi approaches
Initial value solver to treat the initial transient stages
Meccanica dei Fluidi I
38
Information and Introduction
BY-PASS TRANSITION IN A SQUARE DUCT
secondary mean
cross-flow
link with optimal
pertubations?
• Always linearly stable …
• DNS shows a mean
cross-flow motion
• Link with optimal perturbations?
Meccanica dei Fluidi I
39
Information and Introduction
2.
SIMULATION OF 2D AND 3D
COMPLEX MOVING GEOMETRIES
Meccanica dei Fluidi I
40
Information and Introduction
PROPULSION BY BEATING CILIA: THE PLEUROBRACHIA
“Planar” beat patterns of combplates of this small marine invertebrate generate surface waves:
Antiplectic metachronal wave
Propulsion of
the organism
Propulsion
NUMERICAL TECHNIQUE: IMMERSED BOUNDARIES
9 combplates beating at 15Hz
Velocity vectors (arrows) together with velocity magnitude (contours)
3D EFFECTS
The horizontal velocity vectors are magnified 10 times compared to vectors in the vertical plane
3D EFFECTS
Advection of a passive scalar
PASSIVE CILIA: TWO-WAY COUPLING
Meccanica dei Fluidi I
45
Information and Introduction
STRUCTURAL MODELING
Meccanica dei Fluidi I
46
Information and Introduction
Testing the model in the fluid (a vortex pair in a periodic box)
Amplitude of the velocity
Testing the model in the fluid (flow in a channel with one hairy wall)
Amplitude of the velocity
AN ALTERNATIVE TO THE IMMERSED BOUNDARY TECHNIQUE:
THE OVERLAPPING (or CHIMERA) GRIDS METHOD
Meccanica dei Fluidi I
49
Information and Introduction
OVERVIEW OF THE OVERLAPPING GRIDS METHOD
MOVING BOUNDARIES
Meccanica dei Fluidi I
50
Information and Introduction
OVERVIEW OF THE OVERLAPPING GRIDS METHOD
EASILY EXTENDABLE TO MULTIPLE MOVING BODIES
Meccanica dei Fluidi I
51
Information and Introduction
FLAPPING FOILS, i.e. ICARUS DREAM
July 8, 2006
Meccanica dei Fluidi I
52
Information and Introduction
heaving:
pitching:
Meccanica dei Fluidi I
53
Information and Introduction
SIMPLE ENOUGH: CAN USE SINGLE GRID APPROACH
Coarse grid: 64*64 nodes and and domain radius equal to 10 chord lengths
Meccanica dei Fluidi I
54
Information and Introduction
… OR
STILL USE 2 OVERLAPPING GRIDS
Meccanica dei Fluidi I
55
Information and Introduction
HEAVING AIRFOIL, VALIDATION OF START UP
Heathcote, S., Gursul. I., Jet Switching phenomena for a
Plunging Airfoil, 34th Fluid Dynamics Conference and
Exhibit, Portland, Oregon, AIAA-2004-2150, 2004.
Meccanica dei Fluidi I
56
Information and Introduction
HEAVING AIRFOIL, THRUST PRODUCING CASE
Meccanica dei Fluidi I
57
Information and Introduction
HEAVING AIRFOIL, REVERSED VON KARMAN STREET
Heathcote, S., Gursul. I., Jet Switching phenomena for a
Plunging Airfoil, 34th Fluid Dynamics Conference and
Exhibit, Portland, Oregon, AIAA-2004-2150, 2004.
Meccanica dei Fluidi I
58
Information and Introduction
HEAVING AIRFOIL, neutral wake at reduced frequency
Meccanica dei Fluidi I
59
Information and Introduction
THE DRAGONFLY
Vorticity contours on different planes
Instantaneous streamlines
Meccanica dei Fluidi I
60
Information and Introduction
3.
SHOCK AND DETONATION WAVES
Meccanica dei Fluidi I
61
Information and Introduction
SOME NUMERICAL EXAMPLES
COMPRESSIBLE FLOW PAST A CYLINDER WITH AMR
(AUTOMATIC MESH REFINEMENT)
GRID SUMMARY
ITER 1
NG = 5
NP = 11433
IP = 233
ITER 30
NG = 43
NP = 40443
IP = 691
INITIAL CONDITIONS
BEHIND THE SHOCK
r = 2.6
U = 1.250
V=0
P = 3.214
MACH = 0.962
Meccanica dei Fluidi I
INITIAL CONDITIONS
IN FRONT THE SHOCK
r =1
U=0
V=0
P = 0.714
MACH = 0
62
Information and Introduction
SOME NUMERICAL EXAMPLES
COMPRESSIBLE FLOW PAST A CYLINDER WITH AMR
Density
contours
AMR levels
Meccanica dei Fluidi I
63
Information and Introduction
SOME NUMERICAL EXAMPLES
SHOCK/BUBBLE INTERACTION WITH AMR
GRID SUMMARY
ITER 1
NG = 12
NP = 24520
IP = 0
ITER 3000
NG = 12
NP = 27282
IP = 0
INITIAL CONDITIONS
BEHIND THE SHOCK
r = 3.81
U = 2.85
V=0
P = 10
MACH = 0.1485
Meccanica dei Fluidi I
INITIAL CONDITIONS
IN FRONT THE SHOCK
r =1
U=0
V=0
P=1
MACH = 0
64
INITIAL CONDITIONS
WITHIN THE BUBBLE
r = 0.1
U=0
V=0
P=1
MACH = 0
Information and Introduction
SOME NUMERICAL EXAMPLES
SHOCK/BUBBLE INTERACTION WITH AMR
Schlieren countours
Density contours and AMR levels
Meccanica dei Fluidi I
65
Information and Introduction
SOME NUMERICAL EXAMPLES
EXPLOSION IN A CLOSED BOX
Meccanica dei Fluidi I
66
Information and Introduction
Contents of the Fluid Mechanics course
1.
2.
3.
4.
5.
6.
Introductions and basic concepts
Properties of fluids
Pressure and fluid statics
Fluid kinematics
Mass, Bernoulli and energy equation
Momentum analysis of flow systems
7.
8.
9.
10.
11.
Dimensional analysis and p theorem
Internal flows and Moody chart …
(TECH. PHYSICS ?!)
Differential analysis of fluid flows
Approximate solutions of the Navier-Stokes equations
External flows: drag and lift
(IF TIME PERMITS …)
Everybody, including “old” students can take mid-term and final …
Meccanica dei Fluidi I
67
Information and Introduction
Scarica

Introduction - DICAT - Università Degli Studi Di Genova