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