Computational Transport Phenomena
Diskette files (file \library\readme.1st)
(see also file \library\files.txt)
These files can also be downloaded from:
http://www.lehigh.edu/~wes1/wes1.html
The files are organized on the diskette, and in this Web page,
by chapters, e.g., the files for chapter m1 are in subdirectory
m1, etc. All of the code is written in transportable Fortran 77,
and should therefore run on any computer with a Fortran 77 compiler.
Conversion to Fortran 90 is easily accomplished with the following
steps:
(1) Continuation characters in column 6 should be
removed and replaced with & at the end of the preceding
line.
(2) Comment characters C or * in column 1 should be
replaced with !.
Since, with the exception of (1) and (2), Fortran 90 is fully
compatible with Fortran 77, no other changes should be necessary
to run the programs under Fortran 90.
All of the files are thoroughly commented, so only a brief descrip-
tion of their content is given below. The individual files can be
identifed by type according to the extensions
.for subordinate Fortran subroutines called by the main
program for the specifc problem (in a few cases, e.g.,
h1, the main program may be contained in this Fortran
.for file, but in most cases, the main program is
main.for in subdirectory /library)
.dat data for the problem (in a few cases for which
multiple data files are provided, they are given
individual names, e.g., m3.da1, m3.da2, m3.da3)
.m Matlab file for plotting the solutions (in a few
cases for which multiple Matlab files are provided,
they are given individual names, e.g., ma1.m, ma2.m
ma3.m)
The main programs (named main.for in subdirectory /library or with
the .for extension for a specific chapter) contain OPEN statements
that define the input/output file names. Preceding the OPEN
statements are the assignment statements which define the input/
output unit numbers. These may have to be changed to meet the
local computing requirements. The OPEN statements can easily
be found with an editor by searching for the characters "OPEN("
(with no quotes).
In general, the input/output files named in the OPEN statements
are "data" and "output". Thus, the data file for a particular
problem should be renamed "data", e.g., m1.dat is renamed "data"
before executing the program consisting of main.for + m1.for +
m1.dat.
Spatial derivatives in PDEs are, in some cases, computed with DSS
routines (Differentiation in Space Subroutines) contained in a
single file dss.for in subdirectory /library. These routines are
numbered DSS002 to DSS044, and have a *DECK line at the beginning
of each routine, e.g., subroutine DSS002 starts with the line
*DECK DSS002
Therefore, specific routines can easily be located within file dss.for
by using an editor to search for the characters "*DECK DSSXXX" where
XXX is the number of a specifc routine.
Questions concerning the programs on the diskette or in the Web page
can be directed to:
W. E. Schiesser
Iacocca D307
Lehigh University
111 Research Drive
Bethlehem, PA 18015 USA
(610) 758-4264 (office)
(610) 758-5057 (fax)
wes1@lehigh.edu (Internet)
http://www.lehigh.edu/~wes1/wes1.html
/m1:
Laminar boundary layer flow
main.for, m1.for, m1.dat
m1ser.for - series solution for Blasius equation
/m2:
Unsteady laminar flow in a tube
main.for, m2.for, m2.dat, m2a.m, mab.m
m2cl.for - infinite series solution for the centerline
m2ex.for - explicit finite difference solution (equi-
valent to Euler time integration)
/m3:
Nonlinear, front-sharpening convective flow systems
main.for, m3.for, m3.da1, m3.da2, m3.da3, m3.m
(three data sets for NCASE = 1, 2, 3 selected
in m3.for; these data file names must be used
in the OPEN statement in main.for)
/h1:
Heat conduction in a semi infinite system
h1.for, h1.m
/h2:
Heat conduction in a semi infinite system
main.for, h2.for, h2cart.dat, h2cart.m
h2cylr.dat, h2cylr.m
h2sphr.dat, h2sphr.m
"cart", "cylr" and "sphr" denote Cartesian,
cylindrical and spherical coordinates, respectively
The extensions to two and three dimensions are included
in the following files (main.for can be used as the main
program):
Two dimensions - h22d.for, h22d.dat
h2series.for - series solution for Cartesian coordinates
/h3:
Heat transfer in a circular fin
main.for, h3.for, h3.dat
h3ss1.for - steady state solution, version 1 (based
on Bessel functions)
h3ss2.for - steady state solution, version 2 (based
on hyperbolic functions)
/h4:
The Graetz problem with constant wall heat flux
main.for, h4.for, h4.das, h4dal, h4.m
h4.das - data file for short times
h4.dal - data file for long times
/h5:
The Graetz problem with constant wall temperature
main.for, h5.for, h5.dat, h5.m
main.for, h5ivz.for, h5ivz.dat, h5ivz.m - steady
state solution with axial position z as the initial
value variable
h5cl.for - series solution for the centerline temperature
/h6:
Heat exchanger dynamics
main.for, h6.for, h6.dat, h6.m
/ma1:
A dynamic mass transfer model
main.for, ma1.for, ma1.dat, ma1.m
/ma2:
Mass transfer with simultaneous convection and diffusion
main.for, ma2.for, ma2.dat, ma2.m
/ma3:
Transient multicomponent diffusion
main.for, ma3.for, ma3.dat, ma3.m
WES
01AUG97