<< Chapter < Page Chapter >> Page >
Velocity vectors for square with top side sliding

Assignment 12.1

Use the code in the /numer/ directory of the CENG501 web site as a starting point. Add the capability to include arbitrarily specified normal component of velocity at the boundaries. The stream function at the boundary should be numerically calculated from the specified velocity at the boundary. Test the code with Couette and Plane-Poiseuille flows first in the x and then in the y directions. Display vector plot of velocity and contour plots of stream function and vorticity. Show your analysis and code.

Assignment 12.2

Form teams of two or three and do one of the following additions to the code. Find suitable problems to test code.

  1. Calculate pressure and stress. Test problems: Couette and Plane-Poiseuille flows.
  2. Calculate transient and finite Reynolds number flows. Test problems: Plate suddenly set in motion and flow in cavity.
  3. Add option for cylindrical polar coordinates. Use the coordinate transformation, z = ln r . Test problems: line source and annular rotational Couette flow.

Form new teams of two or three with one member that developed each part of the code. Include all features into the code. Test the combined features. Work together as a team to do the following simulation assignments.

Assignment 12.3

Boundary layer on flat plate. Compute boundary layer velocity profiles and drag coefficient and compare with the Blasius solution. What assumption is made about the value of the Reynolds in the Blasius solution?

Assignment 12.4

Flow around cylinder. Assume potential flow far from the cylinder. Calculate drag coefficient as a function of Reynolds number and compare with literature.

Assignment 12.5

Flow around cylinder. Find the Reynolds number at which boundary layer separation occurs.

Fortran on owlnet

Modern FORTRAN compilers in the MSWindows are very user friendly. They have built in project workspace routines, documentation, and integrated editor and debugger. If you use the FORTRAN compiler on Owlnet, it is useful to know how to use make files to compile and link FORTRAN code. The following is an example make file. You may give it the name, makefile , with no extension.

f77 -c streamf1.f; f77 -c bc1.f;f77 -c coef1.f; f77 -c asolve.f;f77 -c atimes.f; f77 -c dsprsax.f;f77 -c dsprstx.f; f77 -c linbcg.f;f77 -c snrm.f; f77 -o exe streamf1.o bc1.o coef1.o asolve.o atimes.o dsprsax.o dsprstx.o linbcg.o snrm.o

You have to give yourself permission to execute the make file with the one-time command,

chmod +x makefile

The line with the -c option compiles the source code and makes an object code file with the extension of .o. The line with the -o option links the object files into the executable file called exe . This executable file can be executed by typing its name, exe , on the command line.

This example makefile is not very efficient because it will compile every source code listed. There are instructions to recompile only the source code that has been modified or "touched" since the last time the makefile was invoked. However, I do not recall the instructions.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Transport phenomena. OpenStax CNX. May 24, 2010 Download for free at http://cnx.org/content/col11205/1.1
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Transport phenomena' conversation and receive update notifications?

Ask