[email protected] http://www.nd.edu/~mjm http://www.nd.

edu/~chegdept
Use of Mathematica as a Teaching Tool
for (Computational) Fluid Dynamics and
Transport Phenomena
Mark J. McCready
University of Notre Dame
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
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
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
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.
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Mixing Experiment
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Coloring
Liquid
Soap
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Toothpaste and Karo Syrup
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
They can be mixed
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Mixing questions using
Dimensional analysis
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Mixing Answers
Viscous mixing

1
2
3 2
, ,
L
R
P
L
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Dimensionless groups do not need
to be on technical subjects
Cr
How Smart You Are
How Smart You Think You Are

[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
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
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Dimensionless Proverb
Child ~ Wise person
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Dimensional Analysis NoteBook
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
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.1.nb. (Notebook)
Mathematica Primer(html)
A basic introduction to dimensional analysis including physical motivation and how to solve pipe flow.
dimensional.analysis.nb. (Notebook)
dimensional.analysis.html (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,
lubricatedflow.nb (notebook format)
lubricatedflow.html (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)
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Finite Difference and Finite Element
solution to flow in a Rectangular Duct
http://www.nd.edu/~mjm/RectangularDuct.nb
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
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.
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
As the aspect ratio is varied from one, the
solution become increasingly inaccurate. The
finite different method works better
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Comparison of solutions
Analytical
Finite diff.
Finite element
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Boundary-Layer flows
http://www.nd.edu/~mjm/boundarylayer.nb
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Boundary-Layer Flow (cont)
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Boundary Layer Flow
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Problem set-up
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Problem set-up continued
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Flow over
a flat
plate
solution
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Flow past a wedge
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Flow
past a
wedge
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Solution
code
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Plots at different angles
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Solution of natural convection boundary - layer
flow near a heated flat plate
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Demonstration of the effect of flow regime on
pressure drop in multifluid flows
http://www.nd.edu/~mjm/Effect.of.Flow.Regime.pdf
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Surprising result
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Creeping flow past a stationary sphere
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Flow past a sphere problem
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
Velocity field magnitude
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
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
[email protected] http://www.nd.edu/~mjm http://www.nd.edu/~chegdept
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.