Syllabus

---

Homework/Project

Homework #1 , due 1:30 PM, Thursday, September 13
Solution

Project #1 , due 1:30 PM, Tuesday, October 2

Update (9/25): A clarification is made in Task 4. See the strikethrough text.
Please use Fluent's built-in function to compute the rate of heat transfer 
at bottom plate. The "alternative" approach similar to that used in HW1 
will actually produce a nontrivial discrepancy in energy balance in this case.
(Roughly speaking, this is due to the effect of heat conduction at the inlet,
which is unaccounted for if using the alternative approach. This was noted
in the additional remark in HW1. This effect is small in HW1 but becomes more
significant in the current problem. We will explain the detail in future lectures.)

Further background:

In case it might be useful, the following is a link to a paper that provides
the empirically determined relation of ρ as a function of T for water.
(See Eq. (1) or Eq. (3) in the paper.) 

NIST 1992 formula for ρ(T) 

By taking the derivative of ρ(T) one can readily obtain the thermal expansion
coefficient. Note that the thermal expansion coefficient (as used in Fluent) 
is defined by β ≡ α−1(∂α/∂T) where α ≡ 1/ρ is specific volume.

Alternatively, the operating temperature, operating density, and thermal
expansion coefficient can also be obtained from look-up tables in 
standard engineering handbooks (some are available online).

Discussion on solution for Project 1

Reference solution #1  (thanks to Doni Tapederi)

Reference solution #2  (thanks to Dhivakar Murugesan)

Reference solution #3  (thanks to Hunter Halversen)

Project #2 , due 1:30 PM, Thursday, October 18

Discussion on solution for Project 2

Reference solution #1  (thanks to Conner Cameron)

Reference solution #2  (thanks to Saurabh Prabhu)

Project #3 , due 1:30 PM, Thursday, November 8

Data for geometry:
faucet_data.txt 
For reference only, Matlab code for generating the data file:
waterfaucet3.m 

Discussion on solution for Project 3

Reference solution #1  (thanks to Richard Nile)

Project #4 , due 1:30 PM, Tuesday, November 27

Important! Correction made on Task 3: "x-velocity" should be "y-velocity".
See corrections highlighted in yellow. An additional minor correction of
typo was made: "steady pressure" changed to "static pressure" in deliverable (i)
of Task 3a.

Additional note: If any of the cases in Task 2 and 3 is in the oscillatory or chaotic
regime, the lift and drag might fluctuate with the number of iteration. In that
case, you can use the averaged value for the report. Alternatively, you may
switch to transient solution and use the time-averaged values of lift and
drag in the report. (This note does not guarantee or imply that any of the
cases in those tasks is oscillatory or chaotic. It is provided here 
"just in case".) This remark (about using averaged values) applies to only
lift and drag. For the contour plots, it is acceptable to show just one plot 
from any phase of the oscillation (if there is an oscillation) or the last point
of iteration. For the deliverables in Task 3 about the pressure and viscous
components of the drag, it is also acceptable to just report the values from 
the last point of iteration.

profile data file, flyingsaucer2DH.txt 

(For reference only. Matlab code that generates flyingsaucer2DH.txt) 

Discussion on solution for Project 4

Reference solution #1  (thanks to Brent Skabelund)

Reference solution #2  (thanks to Marielle Debeurre)

Reference solution #3  (thanks to Ramteja Kondakindi)

Project #5 , due 12:10 PM, Thursday, December 6

Slides (typed)

Slides from the first lecture (8/16)

---

Matlab

Octave (free and open alternative to Matlab)

Basic Matlab programming

Links to individual matlab example codes

---