Skip navigation links
(Wind-US Documentation Home Page) (Wind-US User's Guide) (GMAN User's Guide) (MADCAP User's Guide) (CFPOST User's Guide) (Wind-US Utilities) (Common File User's Guide) (Wind-US Installation Guide) (Wind-US Developer's Reference) (Guidelines Documents)

(Introduction) (Tutorial) (Geometry and Flow Physics Modeling) (Numerical Modeling) (Boundary Conditions) (Convergence Monitoring) (Files) (Scripts) (Parallel Processing) (Keyword Reference) (Test Options)

Files

A number of files are used by Wind-US in the course of a solution. The script file you run to submit Wind-US jobs will assign all the necessary files to their appropriate Fortran unit numbers, so you should not need to do any of that yourself. Each of the support files is described briefly below.

Input Data File (.dat)

Fortran unit number 5

The input data file is the primary control file for Wind-US and must be created for each case you want to run with the code. Input data and code options are entered in this file through the use of descriptive keywords. You should observe the following formatting rules in creating the input data file:

  1. The first three lines of the file are reserved for geometry, flow condition, and arbitrary titles, respectively. Each of these titles may be up to 64 characters long. None of these first three lines may start with the (case-insensitive) word "Include".

  2. Blank lines and comments, beginning with a "/" or "!" character, may be placed anywhere in the file after the first three lines. Both dedicated comment lines and trailing comments are supported. Data file readability may be improved dramatically through the liberal use of comments - for example, separating logical sections of the data file: numerical algorithm, force integrations, test options, etc.

  3. Block data, such as that specified in arbitrary inflow or chemistry cases, must be contiguous. Only keywords corresponding to the block may reside between the beginning and ending block indicators (e.g., CHEMISTRY and ENDCHEMISTRY).

  4. Keywords may be entered in upper or lower case.

  5. Abbreviations for keywords may be used, as long as they are unique. If they are not, you may not get the results you expect. For example, it's not a good idea to use single-letter abbreviations for keywords.

  6. The input data file may include content from other files through the use of the INCLUDE keyword.

The following is an example of a simple input data file:

   Geometric Title
   Flowfield Condition Title
   Optional Title

   !-- Freestream "static" conditions:  Mach, p(psi), T(R), alpha(deg), beta(deg)
   Freestream static 0.9 14.7 530. 4. 0.

   !-- Specify turbulence model
   Turbulence Model Inviscid

   !-- Specify how long to run
   Cycles 15

   End

A copy of the input data file will be written as character data in the root node of the flow file (.cfl). This data can be read with the cflinfo utility.

Grid File (.cgd)

Fortran unit number 11

The computational grid used by Wind-US for a particular case is stored in the grid file. In this file are stored the (x, y, z) coordinates of all computational grid points, zone coupling interpolation factors, and grid reference and scaling data.

This file was originally referred to as the common grid file, so named because the file was formatted according to Boeing's Common File Format (CFF). Wind-US also supports grid files in CGNS (CFD General Notation System)format. [Detailed information on the CGNS standard may be found at the CGNS web site external link.] Both common files and CGNS files are binary, and portable to virtually every hardware platform, except Cray, with no need for explicit data conversions.

Grid files for Wind-US may be created by several mesh-generation codes in either common file or CGNS format. If necessary, the cfcnvt utility may be used to convert a variety of other formats, including PLOT3D xyz format, to common file format.

Zone coupling, reference, and scaling data are added to common grid files using the GMAN program. Common grid files are also used in the CFPOST post-processing package. Neither GMAN nor CFPOST currently support CGNS files, however.

Even when third-party mesh-generation codes are used to create the grid file, GMAN or MADCAP should still be used to perform grid quality checks and confirm that no boundary surface condition remains undefined. Often zone coupling information, particularly for non-point-matched boundaries, is missing from the file.

Keywords: CGNSBASE

Flow File (.cfl)

Fortran unit number 20

The flow file contains the computed flow field. For Navier-Stokes and Euler solutions, the file contains density, momentum, and energy data, and, for viscous solutions, turbulence data. The flow file also records the current solution cycle number to allow the file to be used for solution restarts.

Like the grid file, the flow file was originally referred to as the common flow file, but may now be written in either common file or CGNS format. Several graphical post-processing programs are able to read both common files and CGNS files.

The CFPOST post-processing package may be used with common flow files to produce other files for post-processing and/or to create flowfield plots directly. CFPOST does not currently support CGNS files, however.

Keywords: CGNSBASE

Boundary Data File (.tda)

Fortran unit number 14

Wind-US's boundary data file is used during solution restarts, and results in smoother restarts, especially with higher-order boundary coupling. The file acts as a buffer for the transfer of zone coupling information and a holding bin for data needed throughout a Wind-US run but not stored in the grid or solution file.

Time History File (.cth)

Fortran unit number 19

The time history file is a common file which stores data resulting from the use of Wind-US's history tracking capability. The file contains a lookup table corresponding to the range of computational indices tracked during the run, time stamp data, and flowfield data for each time stamp. Upon completion of a time history run, the auxiliary program thplt may be used to view the contents of the time history file.

Keywords: HISTORY

List Output File (.lis)

Fortran unit number 6

The purpose of Wind-US's list output file is to echo the input from the input data file, track convergence and integration results, record CPU/job statistics, and log code error messages. The auxiliary program resplt may be used to extract convergence and integration data from the list output file and create a GENPLOT-style data file, which may be plotted with the plot command in the CFPOST post-processing package. For parallel-processing runs, this file also contains messages relating to the PVM system and to the allocation and operation of slave processors. Convergence data is written to this file for each of the equation sets solved during the run, including the Euler, Navier-Stokes, and non-algebraic turbulence model sets.

Time Data File (.cft)

Fortran unit number 22

The time data file may be used for storing the computed flow field at the extra time levels required for second-order time differencing, Newton iteration, and/or dual time stepping. Currently by default, the extra time levels for second-order time differencing are stored in the .cfl file, and those for Newton iteration and dual time stepping are stored in the .cft file and linked to the .cfl file. Keywords in the TEMPORAL keyword block may be used to specify where the extra time levels should be stored.

Keywords: IMPLICIT ORDER, NEWTON, TEMPORAL

Edge Data File (.cge)

Fortran unit number 21

The edge data file is used with unstructured grids to store grid-related information about each cell that's required during the solution procedure. The file is created during the initial run, and by default is saved for use in subsequent restart runs, but will be recreated if it's missing. When the file is opened, if necessary it is automatically split into multiple files to keep the size of each file below two gigabytes, with the main file transparently linked to the separate files.

Keywords: DEBUG 8, DEBUG 9

Wind-US Control File (WINDCTRL)

Fortran unit number 1

The control file can be used to modify certain program inputs while the code is running. This is particularly useful on large cluster machines, since small changes can be made without having to stop the job, edit the input (.dat) file, resubmit it, and wait in the queue to run again.

At the beginning of every cycle, the program checks for the WINDCTRL control file. If it exists, the file is opened, checked for recognized commands, and then deleted. Note that the control file must be placed in the directory where the code is running, and that this directory may be different from the one where the wind command was issued. If the -runinplace option to the wind script is not used, the job will be run in a remote directory. The root name of the remote run directory may be specified using the -runroot option, or in response to a command line prompt. The default for the root name of the remote directory is /tmp. The full name of the remote directory will be rootname/userid/basename.scr, where rootname is the root name described above, userid is your userid, and basename is the base name of your .dat file.

Below are the keywords that the control file supports:

    CHECKPOINT [CFLONLY | FULL]   Forces a checkpoint of the file(s). CFLONLY will only update the flow (.cfl) file. FULL is a complete checkpoint, which includes updating the boundary data (.tda) file. If neither option is specified after CHECKPOINT, then CFLONLY will be assumed.
 
CFLMAX cflmax Changes the value of CFL AUTOMATIC CFLMAX for the current run stage.
 
CYCLES end_cycle Changes the number of requested CYCLES for the current run stage. The value end_cycle is the last cycle number to complete before terminating the current SOLVER-STAGE. If the value is less than the current cycle, the stage will terminate with the current cycle.

Example: Starting from a flow solution (.cfl) file at 10000 cycles and an input (.dat) file that requests 5000 cycles, the code would normally complete at 15000 (total) cycles. In the control file, specifying CYCLES 12000 will stop the solution at 12000 (total) cycles, or at the completion of the next cycle if already past 12000. CYCLES 20000 will extend the run to 20000 (total) cycles.
 
TIMELEVELS end_level Changes the number of requested TEMPORAL TIME LEVELS for the current run stage. The value end_level is the last time level to complete before terminating the current SOLVER-STAGE. If the value is less than the current time level, the stage will terminate with the current time level.
 
GLOBALRESIDUAL Activates RESIDUAL_OUTPUT TYPE GLOBAL.
 
ZONALRESIDUAL Activates RESIDUAL_OUTPUT TYPE ZONAL.
 
RELAX value Changes the value of the relaxation factor used with the DQ LIMITER. This value typically does not need to be modified.
 
QUIT-STAGE Terminates the current SOLVER-STAGE and proceeds to the next.

Wind-US Stop File (NDSTOP)

Fortran unit number 1

The stop file is used to stop Wind-US execution cleanly in the middle of a CFD solution. Although a solution may be stopped simply by killing the Wind-US process from the system queue, doing so will not ensure a clean update of all zonal flowfield data to the flowfield file. The stop file - named NDSTOP - provides a means to cleanly shut down a running solution, completing the current cycle or zone, performing zone coupling, updating the flowfield file, and removing all symbolic links to Fortran file unit numbers.

To stop the code, you must place one of two words in the stop file in the directory in which the Wind-US job is running. The word STOP in the stop file signals Wind-US to complete the current cycle (at the end of the last zone) and exit. In single-processing mode, you may also use the word STOPZONE in the stop file to stop Wind-US after completing the zone currently in memory. [In multi-processing mode, STOPZONE does the same thing as STOP.] You may then restart the code in your next run, starting at the next zone or back at zone one. Regardless of the existence of this file, Wind-US will stop if the requested number of cycles has been computed, or if your solution has converged.

During the clean-up procedure at the end of the job, the NDSTOP file is automatically removed.

Note that the directory in which the Wind-US job is running, where the NDSTOP file must be, is not necessarily the one you were in when the wind command was issued. If you don't use the -runinplace option to the wind script, the job will be run in a remote directory. The root name of the remote run directory may be specified using the -runroot option, or in response to a command line prompt. The default for the root name of the remote directory is /tmp. The full name of the remote directory will be rootname/userid/basename.scr, where rootname is the root name described above, userid is your userid, and basename is the base name of your .dat file.

On a Unix system, you might submit a simple "at" job for a later time as follows:

   at 0530 monday
   echo STOP > NDSTOP
   ^D

At 5:30 AM on the next Monday, the system would create the NDSTOP file with the word "STOP" in it. Note that this will not stop the run exactly at 5:30 AM; Wind-US must still complete the current cycle, which may take an hour or more for large cases.

Keywords: RESTART

Temperature and Transition Specification Files

Fortran unit number 45

Wind-US can read temperature and transition files, specifying the temperature or transition to turbulence on any boundary surface, that were created with the older tmptrn utility. This option is included for backward compatibility with WIND 2.0, and is only needed for an initial (i.e., non-restart) run. The wind script copies the files to the run directory, and Wind-US opens the files directly (unit 45) using the file names specified with the TTSPEC keyword. New applications should use the latest tmptrn utility to write the temperature/transition data into the flow (.cfl) file directly.

Keywords: TTSPEC

Chemistry Files (.chm)

Fortran unit number 26

The generalized chemistry files contain all the information Wind-US requires to compute a general chemistry mixture. A chemistry file has a header line, followed by three sections containing data defining the thermodynamic properties, reaction rates, and transport properties.

Header

The first line in the file must contain a single value, ispec, specifying the type of reaction and the format for the finite-rate data. The format for the line is:

   ispec     ISPEC
where

    ispec   Reaction and Format
100, 110, 115 Forward and backward elementary reactions, with two reactants and two products. Backward rate computed using an equilibrium constant. (Format 1)
130, 135, 136, 137 Forward and backward elementary reactions, with two reactants and two products. Forward and backward rates may be specified separately, or one may be computed using an equilibrium constant. For ispec = 137, in addition to elementary reactions with two reactants and two products, "general" and "global" exchange reactions with up to three reactants and three products may also be used. (Format 2)
120, 121 Forward reactions only. Intended for detonation or rapid combustion problems. (Format 3)

(The value ispec is read using a list-directed read statement. The ISPEC label following the value is thus optional, but present in the standard chemistry files supplied with Wind-US. The labels following the values ns, nreq, ndeq, and tfrmin, described in the next two sections, are also optional.)

It should be noted that ispec values of 1, 3, and 4 were used in earlier versions of Wind-US. Values of 1 and 4 are now equivalent to ispec = 100, and 3 is now equivalent to ispec = 130.

Thermodynamic Properties

This section of the chemistry data file contains the information necessary to compute the thermodynamic properties for each species. The general format for this section is:

   THERMODYNAMIC COEFFICIENTS
   Title, line 1
   Title, line 2
   ns     NS
   Curve type
   Information defining species 1
   ...
   Information defining species ns
where ns is the number of species, and Curve type specifies the type of curve fits used to define the thermodynamic properties. The section title THERMODYNAMIC COEFFICIENTS must be the second line in the file, immediately after the line defining ispec. The Curve type line may be omitted (not left blank), or specified as one of SPARKCRV, WINDNASA, and NASA3287. If the line is omitted, the SPARKCRV curve type is assumed.

The information defining each species, and the specific format used, depends on the curve type, as described below.

SPARKCRV

[Note - Because the SPARKCRV format lacks information needed to correctly calculate entropy, it is now considered obsolete and may not be supported in the future. Test 19 must be active to use this format. Users are encouraged to use the NASA3287 format instead.]

This is the original file format used for Wind-US chemistry files. The specific heat at constant pressure Cp for each species is defined by a series of fourth-order polynomials, each valid within a defined temperature range. I.e.,

Cp / R = a1 + a2T + a3T2 + a4T3 + a5T4

where R is the universal gas constant, and T is in K.

Given the above, polynomials may be derived for the enthalpy, entropy, and Gibbs free energy.

H / RT = a1 + a2T/2 + a3T2/3 + a4T3/4 + a5T4/5 + b1/T
S / R = a1 ln T + a2T + a3T2/2 + a4T3/3 + a5T4/4 + b2
G / RT = a1 (1 − ln T) − a2T/2 − a3T2/6 − a4T3/12 − a5T4/20 + b1/Tb2

For the SPARKCRV curve type, the enthalpy reference state is at T = 0 K. I.e., the curve fit for H/RT actually computes

H / RT = (ΔfH(0) + ΔH(0→T)) / RT

where ΔfH(0) is the heat of formation at 0 K, and ΔH(0→T) is the change in enthalpy between 0 K and T.

For each species, the information defining the species and for computing Cp / R, etc., is stored in the .chm file as described in the following table. Note that records 2-3 are repeated for each temperature range.

   
   
Contents of Record   Columns   Format

Record 1
    • Name of species 1-8 a8
• Number of curves (i.e., temperature ranges) defining thermodynamic properties 16-20 i5
• Name, and number per molecule, of constituent elements 25-44 4(a2,f3.0)
• Low-temperature (i.e., constant) value for Cp / R. The default is the value computed from the polynomial curve fit, evaluated at the beginning temperature for the first curve. 51-60 f10.1
• Molecular weight 66-75 f10.3
 
Record 2
• Minimum and maximum temperature for curve 1-30 2e15.5
• Coefficients a1-a3 in equation for thermodynamic properties 31-75 3e15.5
 
Record 3
• Coefficients a4, a5, b1, and b2 in equations for thermodynamic properties 1-60 4e15.5

Example

THERMODYNAMIC COEFFICIENTS
 FROM NASA-RFL-TR-70-3, NASA-CR-111989, MAC LIB FISH NO. N71-38747

    5      NS
SPARKCRV
O2                 3    O  2.   0.   0.   0.             3.5         32.000
        300.000       1000.000    0.37190E+01   -0.25170E-02    0.85840E-05
   -0.83000E-08    0.27080E-11  -0.104419E+04    0.00000E+00
       1000.000       6000.000    0.33160E+01    0.11510E-02   -0.37260E-06
    0.61860E-10   -0.36660E-14  -0.104419E+04    0.00000E+00
       6000.000      15000.000    0.37210E+01    0.42540E-03   -0.28350E-07
    0.60500E-12   -0.51860E-17  -0.104419E+04    0.00000E+00
NO                 3    O  1.N  1.   0.   0.             3.5         30.008
        300.000       1000.000    0.41470E+01   -0.41200E-02    0.96920E-05
   -0.78630E-08    0.22310E-11   0.979001E+04    0.00000E+00
       1000.000       6000.000    0.32210E+01    0.12210E-02   -0.42970E-06
    0.65590E-10   -0.34510E-14   0.979001E+04    0.00000E+00
       6000.000      15000.000    0.38450E+01    0.25210E-03   -0.26580E-07
    0.21620E-11   -0.63810E-16   0.979001E+04    0.00000E+00
... etc.
WINDNASA

[Note - Because the WINDNASA format lacks information needed to correctly calculate entropy, it is now considered obsolete and may not be supported in the future. Test 19 must be active to use this format. Users are encouraged to use the NASA3287 format instead.]

This is the same as the SPARKCRV curve fits described above, except for the addition of a single final record for each species that specifies the heat of formation at 0 K. The same equations are used for Cp/R, etc. However, with the WINDNASA curve fits, the enthalpy (and Gibbs free energy, since G/RT = H/RTS/R) use an enthalpy reference state of 298.15 K instead of 0 K. I.e., the curve fit for enthalpy actually computes

H' / RT = (ΔfH(298) + ΔH(298→T)) / RT

where ΔfH(298) is the heat of formation at 298.15 K, and ΔH(298→T) is the change in enthalpy between 298.15 K and T.

In Wind-US, when ADJUST is specified in the CHEMISTRY keyword block, the values of enthalpy and Gibbs free energy returned by the WINDNASA curve fits are shifted to change the reference state to 0 K. I.e.,

H / RT = H' / RT + ΔH(shift) / RT
= (ΔfH(298) + ΔH(298→T)) / RT + (ΔfH(0) − ΔfH(298) + ΔH(0→298)) / RT
= (ΔfH(0) + ΔH(0→T)) / RT

When ADJUST is specified, ΔfH(0), the heat of formation at 0 K in the above equation, is read from the .chm file. If ADJUST is not specified, the enthalpy reference state is left as 298.15 K, and the line in the .chm file containing the heat of formation at 0 K, if present, is ignored.

For each species, the information defining the species and for computing Cp / R, etc., is stored in the .chm file as described in the following table. Note that records 2-3 are repeated for each temperature range. Except for the last record specifying the heat of formation at 0 K, this is the same as for the SPARKCRV format.

   
   
Contents of Record   Columns   Format

Record 1
    • Name of species 1-8 a8
• Number of curves (i.e., temperature ranges) defining thermodynamic properties 16-20 i5
• Name, and number per molecule, of constituent elements 25-44 4(a2,f3.0)
• Low-temperature (i.e., constant) value for Cp / R. The default is the value computed from the polynomial curve fit, evaluated at the beginning temperature for the first curve. 51-60 f10.1
• Molecular weight 66-75 f10.3
 
Record 2
• Minimum and maximum temperature for curve 1-30 2e15.5
• Coefficients a1-a3 in equations for thermodynamic properties 31-75 3e15.5
 
Record 3
• Coefficients a4, a5, b1, and b2 in equations for thermodynamic properties 1-60 4e15.5
 
Record 4 (only required if ADJUST is used)
• Identifier string starting with the word "Heat" 1-31 a31
• Heat of formation at 0 K, in J/mole 32-46 f15.3

Example

THERMODYNAMIC COEFFICIENTS
 CURVE FITS FROM NASA LEWIS CET86 HIGH TEMPERATURE THERMO DATA BASE

    7      NS
WINDNASA
O2                 2    O  2.   0.   0.   0.             3.5         31.998
        300.000       5000.000  3.1162949E+00  1.5886094E-03 -6.7904360E-07
  1.4714899E-10 -1.1729212E-14 -9.9401794E+02  6.4600671E+00
       5000.000      15000.000  2.5782323E+00  8.5796324E-04 -7.6397647E-08
  1.3412616E-12  3.2564804E-17  1.1504711E+03  1.1400551E+01
Heat of Formation at 0 deg K    0.0
H                  2    H  1.   0.   0.   0.             2.5          1.008
        300.000       5000.000  2.5000000E+00  0.0000000E+00  0.0000000E+00
  0.0000000E+00  0.0000000E+00  2.5474038E+04 -4.5991986E-01
       5000.000      15000.000  2.5000000E+00  0.0000000E+00  0.0000000E+00
  0.0000000E+00  0.0000000E+00  2.5474038E+04 -4.5991986E-01
Heat of Formation at 0 deg K    216024.1
... etc.
NASA3287

This format for the curve fits defining the thermodynamic properties is derived from the one defined in NASA TP-3287. [McBride, B. J., Gordon, S., and Reno, M. A. (2001) "Thermodynamic Data for Fifty Reference Elements," NASA TP-3287/REV1.] A series of curves is again used to define the specific heat at constant pressure, enthalpy, entropy, and Gibbs free energy, with each curve valid within a defined temperature range. I.e.,

Cp / R = a1T−2 + a2T−1 + a3 + a4T + a5T2 + a6T3 + a7T4 + a8T5
H / RT = −a1T−2 + a2T−1 ln T + a3 + a4T/2 + a5T2/3 + a6T3/4 + a7T4/5 + a8T5/6 + b1/T
S / R = −a1T−2/2 − a2T−1 + a3 ln T + a4T + a5T2/2 + a6T3/3 + a7T4/4 + a8T5/5 + b2
G / RT = −a1T−2/2 + a2T−1 (1 + ln T) + a3 (1 − ln T) − a4T/2 − a5T2/6 − a6T3/12 − a7T4/20 − a8T5/30 + b1/Tb2

where R is the universal gas constant, and T is in K.

Like the WINDNASA curve fits, the enthalpy reference state for the NASA3287 curve fits is 298.15 K. When ADJUST is specified in the CHEMISTRY keyword block, the values of enthalpy and Gibbs free energy are shifted to change the reference state to 0 K, as described above for the WINDNASA curve fits.

For each species, the information defining the species and for computing Cp / R, etc., is stored in the .chm file as described in the following table. Note that records 3-5 are repeated for each temperature range.

   
   
Contents of Record   Columns   Format

Record 1
    • Name of species 1-18 a18
• Comments (not used in Wind-US) 19-76 a58
 
Record 2
• Number of curves (i.e., temperature ranges) defining thermodynamic properties 1-2 i2
• Identification code (not used in Wind-US) 4-9 a6
• Name, and number per molecule, of constituent elements 11-50 5(a2,f6.2)
• Flag indicating standard state (not used in Wind-US) 52 i1
• Molecular weight 53-65 f13.5
• Heat of formation at 298.15 K, in J/mole 66-80 f15.3
 
Record 3
• Minimum and maximum temperature for curve 2-21 2f10.3
• Number of coefficients an in equation for Cp / R 23 i1
• Exponents of T in equation for Cp / R 24-63 8f5.1
• Enthalpy difference from 298.15 K to 0 K, in J/mole 66-80 f15.3
 
Record 4
• Coefficients a1-a5 in equations for thermodynamic properties 1-80 5d16.8
 
Record 5
• Coefficients a6-a8 in equations for thermodynamic properties 1-48 3d16.8
• Coefficients b1 and b2 in equations for thermodynamic properties 49-80 2d16.8
 
Record 6 (only required if ADJUST is used)
• Identifier string starting with the word "Heat" 1-31 a31
• Heat of formation at 0 K, in J/mole 32-46 f15.3

Example

THERMODYNAMIC COEFFICIENTS
 CURVE FIT JANAF 6000K DATA EXTRAP. TO 15000K & NASA LEWIS DATA RE-FITTED

    5      NS
NASA3287
CO2               Props & Hf298: TPIS v2,pt1,1991,p27.
 3 l 7/88 C   1.00O   2.00    0.00    0.00    0.00 0     44.00980    -393510.000
    200.000  1000.000 7 -2.0 -1.0  0.0  1.0  2.0  3.0  4.0  0.0         9365.469
  4.94378364D+04 -6.26429208D+02  5.30181336D+00  2.50360057D-03 -2.12470010D-07
 -7.69148680D-10  2.84997991D-13  0.00000000D+00 -4.52818986D+04 -7.04876965D+00
   1000.000  6000.000 7 -2.0 -1.0  0.0  1.0  2.0  3.0  4.0  0.0         9365.469
  1.17696943D+05 -1.78880147D+03  8.29154353D+00 -9.22477831D-05  4.86963541D-09
 -1.89206384D-12  6.33067509D-16  0.00000000D+00 -3.90834501D+04 -2.65268192D+01
   6000.000 20000.000 7 -2.0 -1.0  0.0  1.0  2.0  3.0  4.0  0.0         9365.469
 -1.54440594D+09  1.01683595D+06 -2.56137666D+02  3.36936340D-02 -2.18115756D-06
  6.99132366D-11 -8.84221052D-16  0.00000000D+00 -8.04312703D+06  2.25415288D+03
Heat of Formation at 0 deg K    -393149.56
H2O               CODATA,1989. JRNBS v92,1987,p35. TRC tuv-25,10/88.
 2 l 8/89 H   2.00O   1.00    0.00    0.00    0.00 0     18.01528    -241826.000
    200.000  1000.000 7 -2.0 -1.0  0.0  1.0  2.0  3.0  4.0  0.0         9904.092
 -3.94795999D+04  5.75572977D+02  9.31783351D-01  7.22271091D-03 -7.34255448D-06
  4.95504134D-09 -1.33693261D-12  0.00000000D+00 -3.30397425D+04  1.72420539D+01
   1000.000  6000.000 7 -2.0 -1.0  0.0  1.0  2.0  3.0  4.0  0.0         9904.092
  1.03497224D+06 -2.41269895D+03  4.64611114D+00  2.29199814D-03 -6.83683007D-07
  9.42646842D-11 -4.82238028D-15  0.00000000D+00 -1.38428625D+04 -7.97815119D+00
Heat of Formation at 0 deg K    -238918.95
... etc.

Reaction Rates

This section of the .chm file contains the reaction rate information. The first five lines are:

   FINITE RATE COEFFICIENTS
   Title, line 1
   Title, line 2
   nreq ndeq     NREQ, NDEQ
   tfrmin        TFRMIN
where nreq is the total number of reactions (dissociation + exchange + ionization), ndeq is the number of dissociation reactions (i.e., have a third body), and tfrmin is the temperature in K below which no reactions occur. The section title FINITE RATE COEFFICIENTS must immediately follow the last data line in the Thermodynamic Properties section, with no blank lines in between.

There are three possible formats for the remainder of this section, depending on the value of ispec.

Format 1; ispec = 100, 110, 115

This option may be used for any chemically reacting flow with dissociation and exchange reactions of the form

S1 ↔ S3 + S4
S1 + S2 ↔ S3 + S4

where S1, etc., represent chemical species.

Rate Equations

For dissociation reactions, the forward and backward reaction rate equations are of the form

rf = kf [S1]
rb = kb [S3] [S4]

where [S1] represents the concentration of species S1, etc.

For exchange reactions, the forward and backward reaction rate equations are of the form

rf = kf [S1] [S2]
rb = kb [S3] [S4]
Rate Coefficients

The forward and backward reaction rate coefficients kf and kb are computed using the equations

kf = CTS eD / (KBT)
kb = kf / K

where T is the temperature in K, the ratio D/KB is in K, and KB is the Boltzmann constant. The coefficient C has units (cm3/g-mol) (sec-KS)−1.

For ispec = 100 or 110, the equilibrium constant K is computed using

K = exp(a1 + a2Z + a3Z2 + a4Z3 + a5Z4)
Z = 10,000/T

For ispec = 115 the equilibrium constant K is computed from the change in Gibbs free energy for the reaction. I.e.,

K = e−ΔGo/RT

where ΔGo is the change in Gibbs free energy, R is the gas constant, and T is the temperature. The Gibbs free energy for each species is determined from the information specified in the Thermodynamic Properties section.

Third Body Reaction Rates

For the dissociation reactions the specification of the third body reaction rates depends on the ispec value used.

    ispec   Third Body Reaction Rate
100 Variable third body reaction rates are used. The .chm file contains the names of the third body reactants and the corresponding coefficients C used in the formula for kf. Note that values must be specified for each of the ns species.
110, 115 Variable or average third body reaction rates are used. The .chm file contains the average third body coefficient C, and, for variable reaction rates, the names of the third body reactants and the corresponding third body efficiencies.
File Format

The information needed for each reaction is stored in the chemistry data (.chm) file, with the dissociation reactions listed first. I.e., reaction data is specified in the following order:

   Information defining dissociation reaction 1
   ...
   Information defining dissociation reaction ndeq
   Information defining exchange reaction 1
   ...
   Information defining exchange reaction nreq - ndeq
where nreq is the total number of reactions, and ndeq is the number of dissociation reactions.

The data to be specified in the .chm file is described in detail in the following table.

   
   
Contents of Record   Columns   Format

Record 1 (for all ispec values)
    • Name of first reactant 1-5 a5
• Name of second reactant (blank for dissociation) 9-13 a5
• Name of first product 17-21 a5
• Name of second product 25-29 a5
• Temperature exponent S in equation for kf 33-44 e12.4
D/KB in equation for kf 45-56 e12.4
• Coefficient C in equation for kf (only for ispec = 115). For dissociation reactions, this is the average third body coefficient. 57-68 e12.4
 
Record 2 (for all ispec values except 115)
• Coefficients a1-a5 in equation for K 1-60 5e12.4
 
Records 2+ (for ispec = 110)
For dissociation reactions with variable third body reaction rates, these records specify the third body reactants and the corresponding third body efficiencies. To use an average third body reaction rate with ispec = 110, simply omit these lines.

Only the non-unity efficiencies have to be specified. The first record in this group is

   THIRD [BODY] [EFFICIENCY] nnot1
where nnot1 is the number of third body reactants with a non-unity efficiency. (The brackets are used to indicate that the words BODY and EFFICIENCY are optional.) The names of the third body reactants and the corresponding efficiencies are specified in the next nnot1 records, in free format, as follows:
   species(1)  efficiency(1)
   ...
   species(nnot1)  efficiency(nnot1)
Records 2+ (for ispec = 115)
For dissociation reactions, these records specify the third body reactants and the corresponding third body efficiencies. The format is the same as described above for ispec = 110, with one exception. The first record
   THIRD [BODY] [EFFICIENCY] nnot1
must be specified even when an average third body reaction rate is to be used. To use an average third body reaction rate with ispec = 115, specify a zero value (not blank) for nnot1, and omit the lines listing the species and the corresponding efficiencies.
 
Records 3+ (for ispec = 100)
• Name of third body reactant (only for dissociation reactions) 17-24 a8
• Coefficient C in equation for kf. For the dissociation reactions, this is the coefficient for the specified third body reactant. 25-36 e12.4
 
Record 3 (for ispec = 110)
• Coefficient C in equation for kf. For the dissociation reactions, this is the average third body coefficient. 25-36 e12.4

With Format 1, the data is read in subroutine frtin1, and the rates are computed in rates (for ispec = 100), ratesa (ispec = 110), or ratesadl (ispec = 115).

Example

FINITE RATE COEFFICIENTS
 FROM AIAA 88-0513

    5,3      NREQ,NDEQ
  2000.      TFRMIN
O2              O       O         -1.0       5.9500E+04
   1.335      -4.127      -0.616       0.093      -0.005       5
                O2       2.75E+19
                NO       2.75E+19
                O        8.25E+19
                N        8.25E+19
                N2       2.75E+19
N2              N       N         -1.6       1.1320E+05
   3.898     -12.611       0.683      -0.118       0.006       5
                O2       3.70E+21
                NO       3.70E+21
                O        1.11E+22
                N        1.11E+22
                N2       3.70E+21
NO              N       O         -0.5       7.5500E+04
   1.549      -7.784       0.228      -0.043       0.002       5
                O2       2.30E+17
                NO       2.30E+17
                O        4.60E+17
                N        4.60E+17
                N2       2.30E+17
NO      O       O2      N          1.29      1.9220E+04
   0.215      -3.652       0.843      -0.136       0.007       0
                         2.16E+08
O       N2      NO      N          0.1       3.7700E+04
   2.349      -4.828       0.455      -0.075       0.004       0
                         3.18E+13
Format 2; ispec = 130, 135, 136, 137

This format is the most flexible one in Wind-US, with a variety of options for specifying the necessary information. "Elementary" dissociation and exchange reactions may be specified using ispec = 130, 135, and 136, the same as with Format 1. In addition, for ispec = 137, "general" and "global" exchange reactions may be specified with up to three reactants and products.

Rate Equations

General exchange reactions are of the form

c1S1 + c2S2 + c3S3c4S4 + c5S5 + c6S6

with forward and backward reaction rate equations of the form

rf = kf [S1]a [S2]b [S3]c
rb = kb [S4]d [S5]e [S6]f

Global exchange reactions are forward only, of the form

c1S1 + c2S2 + c3S3c4S4 + c5S5 + c6S6

with corresponding reaction rate equations

rf = kf [S1]a [S2]b [S3]c
rb = 0
Rate Coefficients

For all the ispec values, the forward and backward reaction rate coefficients kf and kb are computed using the equations

kf = Cf TSf eDf / (KBT)
kb = Cb TSb eDb / (KBT)

where T is the temperature in K, the ratios Df/KB and Db/KB are in K, and KB is the Boltzmann constant. The reaction rate coefficients for the j'th reaction, (Cf)j and (Cb)j, have units

(Cf)j :   (cm3/g-mol)Oj−1 (sec-KSf)−1     (Cb)j :   (cm3/g-mol)Oj−1+vj (sec-KSb)−1

where Oj is the order of the reaction (i.e., the total number of moles of reactants), and vj is the number of moles of products minus the number of moles of reactants.

Note that unlike Format 1, with Format 2 the backward reaction rate coefficient kb may be specified independently of the forward reaction rate coefficient. However, if the specified input value of either the forward or backward coefficient Cf or Cb is zero, the corresponding reaction rate coefficient kf or kb is computed from the non-zero reaction rate coefficient and an equilibrium constant that's computed using Gibbs function, just as in Format 1 with ispec = 115. E.g., if the specified Cb = 0, the backward reaction rate coefficient is computed as

kb = kf / K

where the equilibrium constant K is

K = e−ΔGo/RT
Eddy Dissipation Concept

For ispec = 137, the eddy dissipation concept may be used, specified using the EDC keyword in the CHEMISTRY keyword block. Note that although three reactants and products are allowed with ispec = 137, the eddy dissipation concept should only be used with two reactants and forward reactions only, and is only applicable to non-premixed combustion. It's intended for rapid combustion reactions with the fuel as species S1 and the oxidizer as species S2.

With the eddy dissipation concept, a second forward reaction rate is computed in addition to the one described above. The two reaction rates are compared, and the lower rate is used.

The reaction rate equation using the eddy dissipation concept is of the form

r = Cedc (ρ/M1) Cμω min {[S1], [S2] (c1M1) / (c2M2)}

where Cedc is the eddy dissipation coefficient, specified using the EDC keyword in the CHEMISTRY keyword block; ρ is the fluid density; M1 and M2 are the molecular weights of species 1 (fuel) and species 2 (oxidizer); Cμ is a constant equal to 0.09; ω is the specific dissipation rate from the SST turbulence model; [S1] and [S2] are the concentrations of species 1 (fuel) and species 2 (oxidizer); and c1 and c2 are the number of molecules in the reaction for species 1 (fuel) and species 2 (oxidizer).

Note that the eddy dissipation concept may currently be used only with the SST turbulence model.

Damkohler Rate Limiter

Also for ispec = 137, the Damkohler rate limiter may be used, specified using the DAMKOHLER keyword in the CHEMISTRY keyword block.

With the Damkohler limiter, both forward and backward reaction rates are limited such that the ratio Da = τFD / τCH is less than some user-specified maximum Damax. The fluid dynamic time scale τFD is taken as the local time step Δt. The chemical reaction time scale τCH is the value that would allow the concentration of the specie with the smallest negative defect (i.e., the one that will first reach zero concentration) to drop by a ratio of 1/e.

Third Body Reaction Rates

Like in Format 1, for the dissociation reactions the specification of the third body reaction rates depends on the ispec value used.

    ispec   Third Body Reaction Rate
130, 136 Average third body reaction rates are used. The .chm file contains the average third body coefficients Cf and Cb.
135, 137 Variable third body reaction rates are used. The .chm file contains the average third body coefficients Cf and Cb, plus the names of the third body reactants and the corresponding third body efficiencies. Note: currently these ispec values only apply to structured grids.

Note: currently dissociation reactions with variable third body reaction rates (ispec values of 135 and 137) are only available for structured grids.
File Format

The information needed for each reaction is stored in the chemistry data (.chm) file, with the dissociation reactions listed first. I.e., reaction data is specified in the following order:

   Information defining dissociation reaction 1
   ...
   Information defining dissociation reaction ndeq
   Information defining exchange reaction 1
   ...
   Information defining exchange reaction nreq - ndeq
where nreq is the total number of reactions, and ndeq is the number of dissociation reactions.

First, for each of the dissociation reactions, the .chm file contains the following:

   
   
Contents of Record   Columns   Format

Record 1 (for all ispec values)
    • Name of first reactant 1-5 a5
• Name of second reactant (blank for dissociation) 9-13 a5
• Name of first product 17-21 a5
• Name of second product 25-29 a5
• Temperature exponent Sf in equation for kf 33-44 e12.4
Df/KB in equation for kf 45-56 e12.4
• Average third body coefficient Cf in equation for kf 57-68 e12.4
 
Record 2 (for all ispec values)
• Temperature exponent Sb in equation for kb 33-44 e12.4
Db/KB in equation for kb 45-56 e12.4
• Average third body coefficient Cb in equation for kb 57-68 e12.4
 
Records 3+ (for ispec = 135, 137)
These records specify the third body reactants and the corresponding third body efficiencies. Only the non-unity efficiencies have to be specified. The first record in this group must be
   THIRD [BODY] [EFFICIENCY] nnot1
where nnot1 is the number of third body reactants with a non-unity efficiency. (The brackets are used to indicate that the words BODY and EFFICIENCY are optional.) The names of the third body reactants and the corresponding efficiencies are specified in the next nnot1 records, in free format, as follows:
   species(1)  efficiency(1)
   ...
   species(nnot1)  efficiency(nnot1)

Following the data for the dissociation reactions, the .chm file contains the data for each of the exchange reactions. If ispec = 137, this section begins with an optional line containing either the word GENERAL or GLOBAL. I.e.,

   [GENERAL | GLOBAL]
If this line is not present, or if ispec is not 137, only elementary exchange reactions with two reactants and two products may be used, and the data for each reaction is specified as follows:

   
   
Contents of Record   Columns   Format

Record 1
    • Name of first reactant 1-5 a5
• Name of second reactant 9-13 a5
• Name of first product 17-21 a5
• Name of second product 25-29 a5
• Temperature exponent Sf in equation for kf 33-44 e12.4
Df/KB in equation for kf 45-56 e12.4
• Coefficient Cf in equation for kf 57-68 e12.4
 
Record 2
• Temperature exponent Sb in equation for kb 33-44 e12.4
Db/KB in equation for kb 45-56 e12.4
• Coefficient Cb in equation for kb 57-68 e12.4

If ispec = 137, and GENERAL or GLOBAL is specified, general or global exchange reactions with up to three reactants and three products are used, and the data for each reaction is specified as follows:

   
   
Contents of Record   Columns   Format

Record 1
    • Number of molecules and name of first reactant 1-15 f5.1,a10
• Number of molecules and name of second reactant 16-30 f5.1,a10
• Number of molecules and name of third reactant 31-45 f5.1,a10
• Number of molecules and name of first product 46-60 f5.1,a10
• Number of molecules and name of second product 61-75 f5.1,a10
• Number of molecules and name of third product 76-90 f5.1,a10
 
Record 2
• Temperature exponent Sf in equation for kf 3-14 g12.4
Df/KB in equation for kf 17-28 g12.4
• Coefficient Cf in equation for kf 31-42 g12.4
 
Record 2+ (if GENERAL was specified)
• Temperature exponent Sb in equation for kb 3-14 g12.4
Db/KB in equation for kb 17-28 g12.4
• Coefficient Cb in equation for kb 31-42 g12.4
 
Record 3
This section specifies the concentration exponents a-f in the forward and backward reaction rate equations. It starts with a line that must contain either the word STOICHIOMETRIC or EXPONENTS. I.e.,
   {STOICHIOMETRIC | EXPONENTS}
If STOICHIOMETRIC is specified, the exponents a through f are set equal to the stoichiometric coefficients c1 through c6, respectively.

If EXPONENTS is specified, the exponents a through f are specified by the following line, in free format, as follows:

   a  b  c  d  e  f
All six values must be specified. Note that since global exchange reactions are forward only, when GLOBAL is specified the exponents d through f aren't actually used, and may be specified as 0.0.

With Format 2, the data is read in subroutine frtin3, and the rates are computed in ratesbe.

Example

FINITE RATE COEFFICIENTS
 FROM EVAN & SCHEXNAYDER - CONVERTED TO EQUILIBRIUM CONSTANT FORM

    8,4      NREQ,NDEQ
   300.      TFRMIN
O2   1.         O    1. O    1.   -1.0       5.9340E+04   7.20E+18
                                  -1.0       0.0          4.00E+17
H2   1.         H    1. H    1.   -1.0       5.1987E+04   5.50E+18
                                  -1.0       0.0          1.80E+18
H2O  1.         OH   1. H    1.   -1.5       5.9386E+04   5.20E+21
                                  -1.5       0.0          4.40E+20
OH   1.         O    1. H    1.   -1.0       5.0830E+04   8.50E+18
                                  -1.0       0.0          7.10E+18
O2   1. H    1. OH   1. O    1.    0.0       8.4550E+03   2.20E+14
                                   0.0       0.0          1.50E+13
H2   1. O    1. OH   1. H    1.    0.0       5.5860E+03   7.50E+13
                                   0.0       4.4290E+03   3.00E+13
H2O  1. O    1. OH   1. OH   1.    0.0       9.0590E+03   5.80E+13
                                   0.0       5.0300E+02   5.30E+12
H2O  1. H    1. OH   1. H2   1.    0.0       1.0116E+04   8.40E+13
                                   0.0       2.6000E+03   2.00E+13
Format 3; ispec = 120, 121
[Note - This format is considered obsolete, and may be removed in the future. The same capability is available using Format 2 with ispec = 137.]

This option is only for global 1- or 2-reaction chemistry, and only for forward reactions. It is a quick method for simulating detonation problems, for example, in which reactions proceed only forward to completion. It allows dissociation and exchange reactions of the form

c1S1c3S3 + c4S4
c1S1 + c2S2c3S3 + c4S4
Rate Equations

For dissociation reactions, the reaction rate equation is of the form

r = k [S1]a

For exchange reactions, the reaction rate equation is of the form

r = k [S1]a [S2]b
Rate Coefficient

The reaction rate coefficient k is computed using the equation

k = CTS eD / (KBT)

where T is the temperature in K, the ratio D/KB is in K, and KB is the Boltzmann constant. The reaction rate coefficient for the j'th reaction, Cj, has units (cm3/g-mol)Oj−1 (sec-KS)−1, where Oj is the order of the reaction (i.e., the total number of moles of reactants).

Eddy Dissipation Concept

For ispec = 121, the eddy dissipation concept may be used, as described above for Format 2 with ispec = 137.

File Format

The information needed for each reaction is stored in the chemistry data (.chm) file in the following order:

   Information defining reaction 1
   ...
   Information defining reaction nreq
where nreq is the total number of reactions.

The data to be specified in the .chm file is described in detail in the following table.

   
   
Contents of Record   Columns   Format

Record 1
    • Name and number of molecules of first reactant 1-8 a5,f3.1
• Name and number of molecules of second reactant (blank for dissociation) 9-16 a5,f3.1
• Name and number of molecules of first product 17-24 a5,f3.1
• Name and number of molecules of second product 25-32 a5,f3.1
• Temperature exponent S in equation for k 33-44 e12.4
D/KB in equation for k 45-56 e12.4
 
Record 2
• Place holders for future use. Leave blank.
 
Record 3
• Concentration exponents a and b in the reaction rate equation 1-16 2f8.3
• Coefficient C in equation for k 25-36 e12.4

With Format 3, the data is read in subroutine frtin2, and the rates are computed in ratesf (for ispec = 120) or ratesg (for ispec = 121).

Example

FINITE RATE COEFFICIENTS
WESTBROOK-DREYER GLOBAL REACTION MODEL 

    1,0      NREQ,NDEQ
   250.      TFRMIN
C2H4 1. O2   3. CO2  2. H2O  2.    0.0       1.5098E+04
   0.0         0.0         0.0         0.0         0.0         0
0.10    1.65             2.00E+12

Transport Properties

This section of the chemistry data file contains the information used to compute the laminar viscosity and thermal conductivity. The formulation is based on Wilke's law, with either Sutherland's law or a NASA four-coefficient formula used individually for each species. Different constants or coefficients may be specified for different temperature ranges.

For Sutherland's law, the equations are of the following form:

μ/μ0 = (T/(T0)μ) 3/2 ((T0)μ + Sμ) / (T + Sμ)
k/k0 = (T/(T0)k) 3/2 ((T0)k + Sk) / (T + Sk)

For the NASA formula, the equations are [Svehla, R. A. (1995) "Transport Coefficients for the NASA Lewis Chemical Equilibrium Program," NASA TM-4647]:

ln μn = Aμ ln T + Bμ/T + Cμ/T2 + Dμ
ln kn = Ak ln T + Bk/T + Ck/T2 + Dk

The data is read from the chemistry data file, in the following form:

   TRANSPORT COEFFICIENTS [NASA]
   Title, line 1
   Title, line 2
   Information for species 1
   ...
   Information for species ns
The section title TRANSPORT COEFFICIENTS [NASA] must immediately follow the last data line in the Reaction Rates section, with no blank lines in between. If the optional word NASA is in the section title, the NASA formulas are used to compute the viscosity and thermal conductivity; otherwise Sutherland's law is used.

When Sutherland's law is used, the information for each species is stored in the .chm file as described in the following table. The reference viscosity μ0 is in millipoise; the reference conductivity k0 is in BTU/(hour-ft-°R); the reference temperatures and temperature offsets are in °R; and the beginning and ending temperatures for each curve are in K. Note that records 3 and 4 are repeated, in pairs, for each temperature range after the first. I.e., the Sutherland's law constants for both μ and k are listed for the first curve, followed by the constants for both μ and k for the second curve, etc.

   
   
Contents of Record   Columns   Format

Record 1
    • Name of species 1-8 a8
• Number of curves (i.e., temperature ranges) defining μ and k 11-15 i5
• Minimum and maximum temperature for first curve 16-35 2f10.3
• Reference viscosity μ0 for first curve 37-48 e12.4
• Reference temperature (T0)μ for first curve 49-60 e12.4
• Reference temperature offset Sμ for first curve 61-72 e12.4
 
Record 2
• Reference conductivity k0 for first curve 37-48 e12.4
• Reference temperature (T0)k for first curve 49-60 e12.4
• Reference temperature offset Sk for first curve 61-72 e12.4
 
Record 3
• Minimum and maximum temperature for curve 16-35 2f10.3
• Reference viscosity μ0 for curve 37-48 e12.4
• Reference temperature (T0)μ for curve 49-60 e12.4
• Reference temperature offset Sμ for curve 61-72 e12.4
 
Record 4
• Reference conductivity k0 for curve 37-48 e12.4
• Reference temperature (T0)k for curve 49-60 e12.4
• Reference temperature offset Sk for curve 61-72 e12.4

When the NASA formula is used, the format for the information in the .chm file is as described in the following table. With the NASA formula, the viscosity μ is in micropoise, the thermal conductivity k is in microwatts/(cm-K), and the beginning and ending temperatures for each curve are in K. Note that record 2 is repeated for each viscosity temperature range, followed by record 3 repeated for each thermal conductivity temperature range. I.e., the coefficients for all the μ curves are listed first, followed by the coefficients for all the k curves.

   
   
Contents of Record   Columns   Format

Record 1
    • Name of species 1-16 a16
• The letter "V" 35 a1
• Number of curves (i.e., temperature ranges) defining μ 36 i1
• The letter "C" 37 a1
• Number of curves (i.e., temperature ranges) defining k 38 i1
 
Records 2+
• The letter "V" 2 a1
• Minimum and maximum temperature for curve 3-20 2f9.2
• Coefficients Aμ, Bμ, Cμ, and Dμ for curve 21-80 4e15.8
 
Records 3+
• The letter "C" 2 a1
• Minimum and maximum temperature for curve 3-20 2f9.2
• Coefficients Ak, Bk, Ck, and Dk for curve 21-80 4e15.8

If the effective binary diffusion model is being used (specified by the the DIFFUSION EFFECTIVE-BINARY keyword in the CHEMISTRY keyword block), the Lennard-Jones parameters [Lennard-Jones, J. E. (1924) "On the Determination of Molecular Fields. I. From the Variation of the Viscosity of a Gas with Temperature," Proc. of the Royal Society, A106, pp. 441-462; Lennard-Jones, J. E. (1924) "On the Determination of Molecular Fields. II. From the Equation of State of a Gas,'' Proc. of the Royal Society, A106, pp. 463-477] used to compute the effective binary diffusivity for each species must also be included in the .chm file, following the TRANSPORT COEFFICIENTS data described above. Values for these parameters for a variety of species may be found in the literature (e.g., Bird, R. B., Stewart, W. E., and Lightfoot, E. N. (1960) Transport Phenomena, John Wiley & Sons, New York.)

The data is read in the following form:

   LENNARD-JONES PARAMETERS
   Title, line 1
   Title, line 2
   Parameters for species 1
   ...
   Parameters for species ns
The title LENNARD-JONES PARAMETERS must immediately follow the preceding data line, with no blank lines in between.

For each species, the data is read using the format described below. The parameter σ is in angstroms, and ε/KB is in degrees Kelvin.

   
   
Contents of Record   Columns   Format

    • Name of species 1-8 a8
• Lennard-Jones parameter σ 37-48 f12.4
• Lennard-Jones parameter ε/KB 49-60 f12.4

Example

TRANSPORT COEFFICIENTS


O2            1   300.000 15000.000   1.9190E-01  4.9160E+02  2.5000E+02
                                      1.4190E-02  4.9160E+02  4.0000E+02
NO            1   300.000 15000.000   1.3700E-01  4.9160E+02  4.0000E+02
                                      8.4070E-03  4.9160E+02  4.0000E+03
O             1   300.000 15000.000   1.7030E-01  7.5000E+02  1.5500E+03
                                      1.0360E-02  4.9160E+02  2.3000E+03
N             1   300.000 15000.000   1.7890E-01  4.9160E+02  2.3000E+02
                                      6.8900E-02  4.9160E+02  2.3030E+02
N2            1   300.000 15000.000   1.6630E-01  4.9160E+02  1.9200E+02
                                      1.4000E-02  4.9160E+02  3.0000E+02

Keywords: CHEMISTRY

Turbomachinery Data Files

Fortran unit number 77

The effects of turbomachinery in a duct may be simulated using an actuator duct model, originally developed at MIT. [Gong, Y., "A Computational Model for Rotating Stall and Inlet Distortion in Multistage Compressors," Ph. D. Dissertation, Dept. of Aeronautics and Astronautics, Massachusetts Institute of Technology, Cambridge, Massachusetts, 1999.] The user input needed to define the characteristics of the turbomachinery being modeled is read from a set of turbomachinery data files. These are formatted files, with a separate file required for each blade row. Currently each blade row must be in a separate zone. The names of the files, and the zone each file corresponds to, is specified using the TURBOSPEC keyword block. Unless specified otherwise, all input data is non-dimensional.

The first three lines in the file are as follows:

   IROW
   RPM
   IBLADE
where

    IROW   Blade row designation. This value must be the same as the zone number for the zone containing the blade row.
 
RPM Rotational speed of the blade in revolutions per minute.
 
IBLADE Flag defining what data is being specified in the rest of the file, and its format. The possible values are:

0   Input leading and trailing edge blade angles, and force coefficients, both as functions of radial location. Compute the local blade angles using linear interpolation in the axial direction between the leading and trailing edge angles.
2 Input a uniform body force in the axial (i.e., x) direction, plus relaxation factors for the body force source terms in both the x-momentum and energy equations. With this option, no body force source terms are added to the y- and z-momentum equations.
3 Input a uniform body force for a combustor, plus a relaxation factor for the energy equation. With this option, no body force source terms are added to the momentum equations.

The format for the rest of the file depends on the value of the flag IBLADE.

IBLADE = 0

For IBLADE = 0, the format is:
   NBLADE  NANG   NBODY
   RBLE    BLE    RBTE   BTE   (This line occurs NANG times.)
   ...
   RBODY   CKVIN  CKNIN        (This line occurs NBODY times.)
   ...
where

    NBLADE   Number of blades.
 
NANG Number of spanwise (i.e., radial) locations at which the blade leading and trailing edge angles are being specified.
 
NBODY Number of spanwise locations at which the force coefficients are being specified.
 
RBLE, RBTE Spanwise locations (i.e., radii) at which the blade leading and trailing edge angles are being specified.
 
BLE, BTE Blade leading and trailing edge angles, in degrees, at the spanwise locations specified by RBLE and RBTE. The direction of rotation is chosen as the positive direction. Thus, an angle is positive if the flow in that angle's direction has a tangential component in the same direction as the blade rotation. This implies that leading and trailing edge angles will be negative for compressor rotor blades. For stator blades, the leading edge angle will be positive, but the trailing edge angle may be positive or negative.
 
RBODY Spanwise location at which the force coefficients are being specified.
 
CKVIN, CKNIN Parallel (i.e., viscous) and normal force coefficients at the spanwise location specified by RBODY.

IBLADE = 2

For IBLADE = 2, only one line is needed after IBLADE:

   BODYF  RLXFM  RLXFE
where

    BODYF   Uniform body force in the axial (i.e., x) direction
 
RLXFM, RLXFE Relaxation factors for the body force source terms in the x-momentum and energy equations.

IBLADE = 3

For IBLADE = 3, again only one line is needed after IBLADE:

   BODYF  RLXFE
where

    BODYF   Uniform body force in the axial (i.e., x) direction
 
RLXFE Relaxation factor for the body force source term in the energy equation.

Keywords: TURBOSPEC

Memory Log File (memdebug.lis)

Fortran unit number 97

A memory log file, always named memdebug.lis, may be created using DEBUG option 65 to track memory allocation/deallocation requests.

Keywords: DEBUG 65

Reserved Files

Unit numbers 15, 32, and 55 are reserved for proprietary features, and should not be used by Wind-US developers.


Last updated 24 Jul 2018