mjm@[Link] [Link] [Link]

edu/~chegdept
Use of Mathematica as a Teaching Tool
for (Computational) Fluid Dynamics and
Transport Phenomena
Mark J. McCready
University of Notre Dame
mjm@[Link] [Link] [Link]
Outline
Mixing
Importance of dimensionless groups
Experimental study of mixing of viscous materials
Numerical solution to flow in a rectangular duct
Mathematica used to show finite difference vs finite element
Mathematica notebooks on computational fluid
flow and heat transfer problems
Boundary-layer flow, Falkner-Skan problem
Natural convection thermal boundary layer
Mathematica notebooks for other fluid flow
problems
Creeping flow past a sphere
Introduction to multiphase flows
mjm@[Link] [Link] [Link]
Dimensionless groups
A big theme throughout the Junior-Level
Fluid Dynamics course is the importance of
comparing competing or cooperating
effects and how dimensionless groups
inherently do this.
To make the point we did a laboratory
exercise on mixing and combined with
dimensional analysis.
mjm@[Link] [Link] [Link]
Mixing Experiment
mjm@[Link] [Link] [Link]
Coloring
Liquid
Soap
mjm@[Link] [Link] [Link]
Toothpaste and Karo Syrup
mjm@[Link] [Link] [Link]
They can be mixed
mjm@[Link] [Link] [Link]
Mixing questions using
Dimensional analysis
mjm@[Link] [Link] [Link]
Mixing Answers
Viscous mixing

1
2
3 2
, ,
L
R
P
L
mjm@[Link] [Link] [Link]
Dimensionless groups do not need
to be on technical subjects
Cr
How Smart You Are
How Smart You Think You Are

mjm@[Link] [Link] [Link]
Dimensionless
Confucius Proverb
He who knows not and knows he knows not
is a child, teach him, Cr~1
He who knows not and knows not he knows
not is a fool, shun him, Cr<<1
He who knows and knows not he knows is
asleep, awaken him, Cr>>1
He who knows and knows he knows is wise,
follow him Cr~1
mjm@[Link] [Link] [Link]
Dimensionless Proverb
Child ~ Wise person
mjm@[Link] [Link] [Link]
Dimensional Analysis NoteBook
mjm@[Link] [Link] [Link]
Mathematica Note Books
Many of these (and other useful materials) are also available from MathSource, at the Wolfram Research website.
This and other courses that use Mathematica materials can be found at the Mathematica Courseware web site.
A simple Mathematica primer,
Mathematica_primer.[Link]. (Notebook)
Mathematica Primer(html)
A basic introduction to dimensional analysis including physical motivation and how to solve pipe flow.
[Link]. (Notebook)
[Link] (html)
A simple primer on why we use log-log plots and what they mean,
Primer on log-log and semilog plots. (Notebook)
Primer on log-log and semilog plots(html)
An exhaustive solution of the lubricated flow example ("core-annular flow") from Middleman 3.2.3, pp79-82).
It demonstrates a number of Mathematica features and several important basic ideas from this course,
[Link] (notebook format)
[Link] (html, this is not as good as the Mathematica version, but you don't need MathReader.)
This one shows how to use the chain rule to nondimensionalize differential equations. It also makes a point that the
Resulting dimensionless terms are of order 1.
Making a differential equation dimensionless (Notebook format)
Making a differential equation dimensionless (html)
mjm@[Link] [Link] [Link]
Finite Difference and Finite Element
solution to flow in a Rectangular Duct
[Link]
mjm@[Link] [Link] [Link]
Objectives for CFD
BS engineers, whether we like it or not, will
increasingly be using computational
packages and analytical instruments that
are turn key (they dont understand how
they work). We need to instill in them
both a healthy skepticism that they need
to verify the answers, and enough
fundamental understanding of the
different subject so that they can.
mjm@[Link] [Link] [Link]
As the aspect ratio is varied from one, the
solution become increasingly inaccurate. The
finite different method works better
mjm@[Link] [Link] [Link]
Comparison of solutions
Analytical
Finite diff.
Finite element
mjm@[Link] [Link] [Link]
Boundary-Layer flows
[Link]
mjm@[Link] [Link] [Link]
Boundary-Layer Flow (cont)
mjm@[Link] [Link] [Link]
Boundary Layer Flow
mjm@[Link] [Link] [Link]
Problem set-up
mjm@[Link] [Link] [Link]
Problem set-up continued
mjm@[Link] [Link] [Link]
Flow over
a flat
plate
solution
mjm@[Link] [Link] [Link]
Flow past a wedge
mjm@[Link] [Link] [Link]
Flow
past a
wedge
mjm@[Link] [Link] [Link]
Solution
code
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
Plots at different angles
mjm@[Link] [Link] [Link]
Solution of natural convection boundary - layer
flow near a heated flat plate
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
Demonstration of the effect of flow regime on
pressure drop in multifluid flows
[Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
Surprising result
mjm@[Link] [Link] [Link]
Creeping flow past a stationary sphere
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
Flow past a sphere problem
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
mjm@[Link] [Link] [Link]
Velocity field magnitude
mjm@[Link] [Link] [Link]
Conclusions
Some elements of mixing are incorporated
into the fluid dynamics course
Laboratory experiment
Dimensional analysis
The main idea we attempt to convey about
computational fluid dynamics is that it is
wonderful if it works, but make sure your
solution is correct.
Strategy is like using different excess Gibbs Free
Energy models to design distillation columns with
a process simulator
mjm@[Link] [Link] [Link]
Conclusions (cont.)
Mathematica notebooks can be used to
show students
Computations
To do algebra that is too tedious for them to do
To allow them to explore the solution
To incorporate other media
Questions remain as to if our approach
gives a significant or incremental benefit
to the students.