Several user-controlled options have been provided as an aid
to modifying Wind-US.
These options all may be selected by using the keyword TEST
in the input data file.
The test options typically control program features which are under test,
and have not been accepted for production use.
Each test option is of the form
As test options are accepted, they are "hard wired" into the code
and the test option described here becomes meaningless.
Test options may also be rejected based on trial runs.
In that case, the option code described here will also become meaningless.
For this reason, the list of valid options is not consecutive.
Selection of an invalid option will be accepted by Wind-US,
but will have no effect.
The various test options and modes are described below.
For each option, the subroutines referencing that option are listed in
parentheses.
Non-Production Test Options
|
| number |
| Description
|
|---|
|
| 2 | | Designed for parallel processing data
transfer debugging
| |
mode | |
Result
|
|---|
|
1 | |
Don't read zonal boundary data
(evrwbd)
| |
2 | |
Not used
| |
4 | |
Don't solve zone
(evsolv)
| |
8 | |
Don't write zonal data
(evwzon)
| |
16 | |
Don't write zonal boundary data
(evrwbd)
| |
32 | |
Not used
| |
64 | |
Not used
| |
128 | |
Not used
| |
256 | |
Don't update boundary conditions
(lpschm)
|
Set mode equal to the sum of the desired actions.
|
| 3 | | Parallel processing
task tracing
| |
mode | |
Action
|
|---|
|
1 | | Trace event reads/writes
(rwev)
| |
2 | | Trace file I/O
(rwnh, rwd, rwi, rwr)
| |
4 | | Trace network traffic
(psexit, psrwev, psrwgv, psrwnh, psrwc, psrwd, psrwi,
psrwr, rwbc)
| |
8 | | Trace task begin/end
(psspwn, tskbeg, tskidl)
| |
16 | | Print task queue for debugging
(psqprt)
|
Set mode equal to the sum of the desired actions.
I.e., setting mode = 5 will trace both event reads/writes
and network traffic.
|
| 5 | | Flag for procedure used in gas1
to compute effective specific heat ratio β and sound
speed a for frozen and finite-rate chemistry.
| |
mode | |
Procedure
|
|---|
|
|
> 0 | | Iterate on pressure,
with a maximum of mode iterations
| |
< 0 | | Iterate on temperature,
with a maximum of 20 iterations
|
The default is to iterate on pressure, with a maximum of 20
iterations.
(gas1)
|
| 6 | | Write .cfl file compatible with the
"old" code
(asnsx, asvisc)
|
| 7 | | Do not use high
performance C I/O interface
(openf)
|
| 8 | | Use Version 2
common files
(cftopn, mpinit, openf, zhinit)
|
| 9 | | For mode = 1,
run a verification case using the method of manufactured
solutions.
If mode = 2, in addition to running the case, Plot3d
function files named mms.exact.fxn, mms.cfd.fxn,
and mms.error.fxn are written containing the exact
solution, the Wind-US solution, and the error.
The functional forms for the primitive variables are combinations
of sine and cosine functions, and are coded in subroutines
mms_3d_euler_bc and mms_3d_ns_bc for the Euler
and Navier-Stokes equations, respectively.
The corresponding source terms for the governing equations
are coded in subroutines mms_3d_euler_src and
mms_3d_ns_src.
The choice of an Euler or Navier-Stokes solution is determined by
the TURBULENCE keyword.
If the flow is inviscid in all zones, an Euler solution is
computed; otherwise, a constant-viscosity laminar Navier-Stokes
solution is computed.
The flow may be 2-D or 3-D, as determined by the grid file.
Multi-zone grids may be used, but the use of
COUPLING MODE
CHARACTERISTIC is recommended for Navier-Stokes solutions.
(Viscous terms are neglected in Roe boundary coupling, and they
are large in the method of manufactured solutions.)
Frozen boundary conditions must
be specified in the grid file at all non-coupled boundaries.
The appropriate boundary conditions for the functional forms being
used will then be set automatically.
See Roach [Roach, P. J. (1998) Verification and Validation in
Computational Science and Engineering, Hermosa Publishers,
New Mexico], and Salari and Knupp [Salari, K., and Knupp,
P. (2000) "Code Verification by the Method of Manufactured
Solutions," SAND 2000-1444, Sandia National Laboratories,
Albuquerque, New Mexico], for more detailed information on the
method of manufactured solutions.
(l2norm1d, lpschm, mms_test_soln, nsrhss, pstinp, sins, tdbcni)
|
| 10 | | Print the time
step information (i.e., minimum/maximum CFL and Δt)
into the .lis file every mode cycles, instead of
just on the first cycle.
(NSzsolv)
|
| 11 | | Normally, the flowfield is rotated to be
consistent with changes in the global angles of attack and yaw.
This TEST option prevents that rotation.
(Changes in the global angles may occur at the start of a restart
run, if the user specifies angles that are different from the
values in the .cfl file, or during a run when the
FIXED_CL keyword is used.)
(lpgrp)
|
| 12 | | In the
HLLC scheme, use a more accurate
method for computing the wave speeds.
(hllc, US_HLLC)
|
| 14 | | Add pressure to
the .cfl file for structured grids.
(It's always there for unstructured grids.)
(axflow)
|
| 15 | | For the SST model
with an ε limiter (the
LESB
keyword), modify the grid filter width Δ used in the limiter.
This test option only applies to structured grids.
(sst1, sst2)
| |
mode | |
Grid Filter Width
|
|---|
|
0 | |
max(ds1, ds2,
ds3, V dt,
k1/2 dt)
| |
1 | |
max(ds1, ds2,
ds3,
k1/2 dt)
| |
2 | |
max(ds1, ds2,
ds3, V dt)
| |
3 | |
max(ds1, ds2,
ds3)
|
where ds1, ds2, and
ds3 are the distances along the grid lines in
the three directions, V is the velocity, k is the
turbulent kinetic energy, and dt is the time step size.
|
| 17 | | Use "new" Baldwin-Barth turbulence model
(bbarth)
|
| 20 | | Non-dimensionalize k and ω
in the SST model the "old" way
(aijkrg, aikeps, sstprtinp, sstpstprt)
|
| 21 | | Spalart-Allmaras
turbulence model.
(See the
TURBULENCE keyword.)
(goldbergprtinp, goldbergpstprt, redimsa, sabound, saprtinp,
sapstprt, sinut, spalart, TURB_goldbergboundUS,
TURB_saboundUS)
| |
mode | |
Result
|
|---|
|
0 | | Use original 1992 model,
with an ft2 term for laminar
stabilization and a default freestream value for
νt of 5.0.
| |
1 | | Like mode 0, except with
"corrections" to the production and destruction terms.
This is equivalent to the default model for WIND beta versions
4.15 to 4.92.
| |
2 | | Like mode 1, except without
the ft2 term, and with a default
freestream value for νt of 0.1.
In addition, the initial value of the dependent variable
is set to the freestream νt.
This is equivalent to the default model for WIND versions
prior to WIND beta 4.15, and includes a slight error that
makes the model overly dissipative.
|
|
| 25 | | In the Baldwin-Lomax model, use
y+ based on wall vorticity
(blomax)
|
| 26 | | Use local values in y+
damping for the Baldwin-Lomax, Cebeci-Smith, Baldwin-Barth, and
k-ε models
(bbdamp, blomax, cebeci, kepy2)
|
| 29 | | For the Cebeci-Smith
model, use y+ based on wall vorticity.
(algtur, cebeci)
For the k-ε models, TEST 29 is
a production limiter
(kelhssch, kelhssrg, keprod, kerhssch, kerhssrg,
redimkeps3d)
| |
mode | |
Result
|
|---|
|
0 | | Production limited to
20 × dissipation
| |
1 | | Production computed from
vorticity, but not limited
| |
2 | | Production not limited
| |
3 | | Production computed from
vorticity, then limited
| |
4 | | Production computed using
2μtSijSij, but not
limited
| |
5 | | Production computed using
2μtSijSij, then
limited
|
Modes 4 and 5 only apply to the Rumsey-Gatski ASM
k-ε model.
|
| 30 | | Irrotational boundary condition at freestream
inflow boundaries
(BC_bcfree)
|
| 31 | | Use quadratic fit
for normal derivative boundary conditions
(BC_bcwall)
|
| 35 | | Flag for LOADS
summary file
(frcsum)
| |
mode | |
Result
|
|---|
|
1 | | Write file in "old format"
| |
2 | | Write loads on unstructured
surfaces
|
|
| 37 | | In wbnd2, deallocate memory used
for boundary coupling data before returning
(wbnd2)
|
| 40 | | Old switch for compressibility correction
in k-ε model; use
K-E
COMPRESSIBILITY CORRECTION instead.
(keppstin)
|
| 46 | | In the
SST turbulence model,
in blowing regions and
bleed regions with a specified
negative bleed flow rate, set
μturb = 10μlam
along the wall.
(sstbound, TURB_sstboundUS)
|
| 47 | | For the algebraic turbulence models, smooth
turbulent viscosity in each i-plane using simple averaging.
mode = number of smoothing passes.
(smtvis)
|
| 48 | | For the algebraic turbulence models, smooth
turbulent viscosity in three dimensions using simple averaging.
mode = number of smoothing passes.
(smtvis)
|
| 49 | | Modified Runge-Kutta
smoothing (see the
SMOOTHING keyword)
(bdload, dampi, dampj, dampk)
| |
mode | |
Result
|
|---|
|
2 | | No pressure switch on
second-order dissipation
| |
3 | | No pressure switch, and an
LES type filtering of nonlinear terms
|
|
| 51 | | Limit the
turbulent viscosity μT, such that the
maximum value of
μT / (μL)∞ =
mode × 1000.
Suggested range is 50 < mode < 100.
Do not use this option with the k-ε models;
use K-E
MAXIMUM TURBULENT VISCOSITY instead.
(keppstin, mutlim, US_UpdateTurb)
|
| 52 | | When using
BLOW PLENUM, print a
warning when the plenum total pressure is automatically raised
because it was less than the local static pressure
(BC_bcbled)
|
| 54 | | Reserved for use at Boeing
|
| 55 | | Reserved for use at Boeing
|
| 56 | | No energy addition to fluid due to MFD equations
(emdef)
|
| 57 | | Implicit terms on for the Spalart-Allmaras
and SST turbulence models
(spalart, sst1, sst2)
|
| 58 | | Store the Lorentz force in the .cfl
file instead of the electric field
(emdef, emzwrt)
|
| 59 | | Apply
SMOOTHING keyword
values to the Spalart-Allmaras model as well as the mean flow solver.
(Currently deactivated.)
(spalart)
|
| 60 | | Print warning when
local and coupled normal velocities have inconsistent directions at
a coupled boundary
(postrbs)
|
| 61 | | When mode = 2,
all boundary conditions are applied, whether or not they're
consistent with the IBLANK values.
This only affects corners, where there are usually multiple
boundary conditions.
So, if a wall boundary at j = 1 meets an
outflow boundary at i = imax,
if TEST 61 2 is specified, both boundary
conditions are applied.
(tdbcgs)
|
| 62 | | Compute cell
areas and volumes using procedure from WIND 4
(mphzmet, nsrhsv, tdarea, vismet)
|
| 63 | | Eliminate the
"fat" boundary cells in any coordinate direction.
mode = 1, 2, or 4 indicates the i, j,
and k direction, respectively.
Set mode equal to the sum of the desired directions.
I.e., setting mode = 5 will eliminate the "fat"
boundary cells in the i and k directions.
(mphzmet, nsrhsv, tdarea, tdarea2, vismet)
|
| 64 | | Remove dt from dq when
computing residuals
(l2norm1d)
|
| 65 | | In marching solutions, lower the CFL
number for the last marching step
(NSzsolv)
|
| 66 | | Don't update β in gas3, for
ireal = 2 and ispec = 2.
This test option is not recommended but will decrease run time.
(gas3)
|
| 67 | | Bleed/blowing
boundaries in turbulence model solution
(kebc, sabound, sstbound, TURB_goldbergboundUS,
TURB_saboundUS, TURB_sstboundUS, vprtbound)
| |
mode | |
Result
|
|---|
|
0 | | Treat bleed/blowing
boundaries as no-slip walls.
| |
1 | | Treat bleed/blowing
boundaries as slip walls in the Spalart-Allmaras, Chien
k-ε, and SST models.
This was the default behavior prior to WIND 5.101.
| |
2 | | Same as mode 0, except
do not use the procedure of Wilcox [Wilcox, D. C.,
Turbulence Modeling for CFD, DCW Industries, 2000]
to compute the boundary condition for omega in the SST
model.
|
|
| 68 | | If the density is zero at a coupled boundary, issue a warning, ignore
the coupling data, and continue.
The default is issue an error message and abort.
(postrbs)
|
| 70 | | Tolerance level for converging gas properties
P, ρ, or T in gas2.
Tolerance level = 0.1mode.
(gas2)
|
| 71 | | Calculation
of thermodynamic properties from curve fit equations in
.chm file.
(aichem, aijkrg, aikeps, airgun, aixyzrg, BC_bcfreebc,
BC_bcfreechar, BC_ijkfreechar, BC_usrffreechar, chpstin1,
chrhsv, cpfun, genturb, gibfn, hfun, kwinflow, phinit,
pstinp, tdutv1, therm1, transp, turbupd, uspeci)
| |
mode | |
Result
|
|---|
|
0 | |
Abort if temperature exceeds the maximum for the curve
fits
| |
1 | |
Extrapolate if temperature exceeds the maximum for the curve
fits, and write a warning message to the .lis file
| |
2 | |
Same as mode 1, except don't write a warning message
| |
3 | |
If temperature is outside the range of the curve fits,
use the values at the nearest limit
| |
4 | |
Use constant values consistent with a gas constant of
287 m2/sec2-K, and a specific heat
at constant pressure of 1004.5 m2/sec2-K
| |
5 | |
If temperature is outside the range of the curve fits,
set specific heat to value at the nearest limit and
extrapolate for remaining properties
|
|
| 72, 73 | | For
structured grids, periodic boundaries may be established
by setting up double (or larger) fringes at the boundaries
using GMAN, with the boundary condition type specified as
frozen.
TEST options 72 and 73 provide the additional information
needed to apply the boundary condition.
TEST option 72 specifies the direction(s) of
periodicity, where values of 1, 2, and 4 indicate the i,
j, and k direction, respectively.
Set mode equal to the sum of the desired directions.
E.g., setting mode = 5 means periodicity in both
the i and k directions.
TEST option 73 specifies the depth of the fringes.
As an example, with TEST 72 1 (periodicity in the
i direction) and TEST 73 2 (a double fringe), we
have
| |
Conditions at | |
Come from
|
|---|
|
I = 1 | | I = IDIM - 3
| |
I = 2 | | I = IDIM - 2
| |
I = IDIM - 1 | | I = 3
| |
I = IDIM | | I = 4
|
This procedure may only be used for periodic boundaries within a
single structured zone, with point-matched grids in the source
and target regions.
It's primarily useful for cases with higher-order differencing
schemes, since the order of the scheme is preserved across the
boundary.
(kebc, sabound, sstbound, tdbcni, vprtbound, vsctbound)
|
| 74 | | Include
k-direction spacing when computing minimum time step in
2-D flows
(tdtmst)
|
| 84 | | Use "old" viscous metric calculation
(dsolv, vismet)
|
| 85 | | Check for zero volumes when computing
viscous metrics
(dsolv, vismet)
| |
mode | |
Result
|
|---|
|
1 | | Check; if ≤ 0 print
message and continue
| |
2 | | Check; if ≤ 0 print
message and stop
|
|
| 87 | | Freezes supersonic inflow at initial
conditions
(BC_bcfreechar)
|
| 88 | | Bypass negative T check in
tdgas.
This is needed for chemistry if SHF (heat of formation) varies widely
since we only have an old SHF to use to estimate T.
(lpschm, tdgas)
|
| 89 | | Use "old" species flux correction method
(gas1)
|
| 90 | | Chemistry stuff
(chinv)
| |
mode | |
Result
|
|---|
|
0 | | Analytic chemistry Jacobian
(ns = 5 only)
| |
1 | | Householder chemistry
Jacobian (ns > 5)
| |
2 | | Solves chemistry source term
explicitly
|
|
| 91 | | Gas constant
(BC_bcfreebc, chpstin1, therm1)
| |
mode | |
Result
|
|---|
|
1 | | β = γ =
β∞
| |
2 | | β = γ = 1.4
|
|
| 92 | | Utilize operator
splitting for the reacting chemistry source terms to increase the
stability of the integration, allowing more efficient solution of
the coupled system.
A 4th-order Pade approximation is used to integrate the reaction
source terms, with mode setting the number of subiterations.
Setting mode = 0 indicates no operator splitting.
(chimplicit, chrhss, US_chemsrc, US_chimplicit,
US_GaussSeidel)
|
| 93 | | In chemistry
species diffusion terms, omit the diffusion gradient in the
conduction term.
This test option only applies to structured grids.
(chrhsv, prtinp, rhsvfl, tdutv1)
|
| 94 | | Turn off implicit chemistry terms.
This test option only applies to structured grids.
(tdimafk, tdimfu, tdimja, tdutaa)
|
| 95 | | Turn off chemistry source term
(rates, rates1, rates1a, rates2, rates2a, rates3, rates3a,
rates4, rates4a, ratesa, ratesaa, ratesadl, ratesadla,
ratesarr, ratesb, ratesba, ratesbe, ratesbea, ratesf,
ratesfa, ratesg, ratesga)
|
| 96 | | Apply chemistry source term over mode
iterations for finite-rate non-equilibrium chemistry
(US_chemsrc, chrhss)
|
| 97 | | P. D. Thomas turbulence model scanning direction.
By default, Wind-US starts at viscous walls and moves into the field.
This test option forces the code to calculate turbulent parameters
from any boundary, in addition to walls.
(algtur)
| |
mode | |
Result
|
|---|
|
0 | | use j lines
| |
1 | | use k lines
| |
2 | | use j and k lines
|
|
| 99 | | Initialize finite rate chemistry with Liu
and Vinokur curve fits
(gas2)
| |
mode | |
Result
|
|---|
|
0 | | Do not track the species
(valid to 50K?)
| |
1 | | Track the species
(valid to 10K?)
|
|
| 100 | | Characteristic time-stepping boundary condition
(BC_bcfree, BC_test100)
| |
mode | |
Result
|
|---|
|
0 | | Second order, with
limit of ΔQ ≤ Q/2
| |
1 | | 1st-order, original
characteristic treatment
| |
2 | | 2nd-order, original
characteristic treatment (only available for structured
grids)
| |
3 | | 1st-order, Roe's average
characteristic treatment
| |
4 | | 2nd-order, Roe's average
characteristic treatment (only available for structured
grids)
|
|
| 102 | | Use time-averaged back pressure for mass flow
boundary condition
(BC_pdsmfr)
|
| 104 | | Treatment of implicit viscous terms
(rcutv1, tdutv1)
| |
mode | |
Result
|
|---|
|
0 | | Use simple implicit
viscous terms
| |
1 | | Turn off implicit
viscous terms
| |
2 | | Use fully-implicit
viscous terms
|
|
| 105 | | Time step type
(iterprtinp, tdtmst)
| |
mode | |
Time Step Type
|
|---|
|
0 | | Δt =
CFL / max (λξ, λη,
λζ)
| |
1 | | Flow angle scaling,
Δt = CFL × (fξ Δξ +
fη Δη +
fζ Δζ), where
fξ =
[1 + tan θ + tan ψ]1/2
fη = fξ tan θ
fζ = fξ tan ψ
| |
2 | | Velocity scaling,
Δt = CFL × min (fξ Δξ,
fη Δη,
fζ Δζ), where
fξ = u / uξ /
|uξ + c|
fη = u / uη /
|uη + c|
fζ = u / uζ /
|uζ + c|
| |
3 | | Δt =
CFL × min (Δξ, Δη,
Δζ) / (|u| + c)
| |
4 | | Δt =
CFL / (λξ + λη +
λζ)
|
|
| 106 | | Compute
the time step at the start of every cycle (even when Newton
time stepping is being used), instead of at the start of every
iteration.
(axflow, lpgrp, lpmg, redim)
|
| 108 | | Extrapolate freestream outflow
(BC_bcfree, BC_bcfreebc)
| |
mode | |
Mach | |
Outflow Conditions
|
|---|
|
0 | | < 1 | | Hold upstream
running characteristic at freestream
| |
| | > 1 | | Extrapolate all, even in
boundary layer
| |
1 | | All | | Extrapolate all, even in
boundary layer
|
|
| 109 | | Boundary flux treatment
(roewal, tdup1)
For tdup1:
| |
mode | |
Result
|
|---|
|
0 | | Conservative
| |
1 | | Upwind extrapolation from
interior
|
For roewal:
| |
mode | |
Result
|
|---|
|
0 | | Characteristic inflow,
conservation if flow parallel to wall
| |
1 | | Characteristic regardless
| |
1000 | | Use conservative wall
treatment at all boundaries
|
|
| 110 | | Grid area variation
limiting.
Not allowed for iorder > 24, i.e., for the following
Roe and Van Leer explicit operators: third-order fully upwind,
fourth-order upwind-biased, fourth-order central, and fifth-order
upwind-biased.
(roecof)
| |
mode | |
A2 / A1
|
|---|
|
0 | | ∞
| |
1 | | 2.0
| |
2 | | 1.5
| |
3 | | 1.33
| |
4 | | 1.1
|
|
| 111 | | Singular matrix check
(jacpr4, jacpr5, jacprg, jacprg4, tdsol4, tdsol9, tdsol11,
tdsolg, tdsolg4, tdsolv)
| |
mode | |
Result
|
|---|
|
1 | | Check, but don't print
results
| |
2 | | Don't check
|
|
| 112 | | Corrected upwind scheme at boundaries.
Defaults to corrected scheme, mode > 0 uses second order
smoothing with mode / 1000 as the smoothing level.
Users should not use this option.
(rhsupw)
|
| 113 | | Check for reverse flow at inflow
and outflow boundaries
(BC_bcconf, BC_bcfree)
| |
mode | |
Result
|
|---|
|
0, 1 | | Print a warning message
and continue
| |
2 | | Print an error message
and stop
| |
3 | | Set velocity to zero
and continue, with no warning message.
This mode only applies at
outflow boundaries.
|
|
| 114 | | Central difference ζ operator.
(tdup1)
| |
mode | |
Result
|
|---|
|
0 | | Upwind
| |
n | | Central smoothing
coefficient = n / 1000
|
|
| 115 | | Do not rescale inviscid wall total velocity
to equal adjacent value, just subtract the normal component from the
adjacent value
(BC_bcvel)
|
| 116 | | Set inward pointing normal to zero in
tdbcm1 at unknown grid topology points.
(tdbcm1)
|
| 117 | |
Freeze inflow boundaries, even in subsonic flow
(BC_bcfreebc, BC_tst100bc)
| |
mode | |
Result
|
|---|
|
1 | | Freeze all inflow
| |
2 | | Freeze only arbitrary inflow
points
| |
3 | | Freeze characteristics on all
i = 1 boundaries
|
|
| 118 | | Singular axis fix
(radavg)
| |
mode | |
Result
|
|---|
|
0 | | Average density, momentum
components, and pressure
| |
1 | | Average density, velocity
components, and pressure
|
|
| 121 | | Under-relaxation of points adjacent to singular
axis
(bcpinw, bcsing, kebc, relsng, sabound, sngthrm,
sstbound)
| |
mode | |
Result
|
|---|
|
0 | | Value on axis is a
radius-weighted average of the values at the adjacent
points; values at the adjacent points are unchanged
| |
n | | Value on axis
computed as for mode 1; values at the adjacent points are
computed from
Fadj = (1 − r) Fadj +
r Faxis
where Fadj is the value at the adjacent point,
Faxis is the axis value, and
r = n / 1000.
|
|
| 122 | | Allow left-handed coordinates
(tdarea1, tdarea2)
|
| 123 | | Track back
pressure, mass flow, and integrated total pressure for outflow
boundary conditions, writing the values into the .lis
file.
The back pressure will be in psi, and the integrated total
pressure will be the ratio of the boundary value to the
freestream value.
When the MASS FLOW
keyword is used, the mass flow value will actually be the ratio
of the computed value to the desired value.
For the other outflow boundary conditions the dimensional mass
flow is written, in lbm/sec.
Note that when the
MASS FLOW keyword
is used, the back pressure and the mass flow are automatically
written into the .lis file, whether this test option is
used or not.
Specifying TEST 123 will add the integrated total
pressure.
Also note that when the
MASS FLOW keyword
is used, all three values may be extracted from the .lis
file using the resplt
(or resplt.pl)
utility.
For the other outflow boundary conditions, only the integrated
total pressure can be extracted.
(BC_IntgrtBnd, BC_pdsmfr)
|
| 124 | | Write convergence
information to list output (.lis) file every iteration
instead of every cycle
(print_res)
|
| 126 | | Compressibility correction to Baldwin-Lomax
turbulence model
(blomax)
| |
mode | |
Result
|
|---|
|
0 | | No compressibility correction
| |
1 | | κ = 0.0180 for
Baldwin-Lomax model (CFL3D uses this)
|
|
| 127 | | Scale printed residual by maximum residual
over all time steps
(print_res)
|
| 128 | | Check the L2 norm
of the residual for convergence instead of the maximum residual
(l2norm1d)
|
| 131 | | For boundary layers on
j = 1 walls, set the time step in the boundary
layer to a (larger) "outer" value, defined as the value at
j = mode.
I.e., (Δt)j =
(Δt)mode for
j < mode.
(tdtmst)
|
| 132 | | Renormalize, changing from total to
static values.
Normalizing values in the .cfl file are unchanged.
(redim)
|
| 134 | |
Second order characteristic extrapolation for adjacent conditions
(BC_bcptcond)
| |
mode | |
Result
|
|---|
|
2 | | First order
| |
0, 1 | | Second order, using
a minmod limiter
|
|
| 135 | | Resets the time step using a weighting
function between the ordinary Euler CFL number and a new "viscous
CFL number", for convergence acceleration in viscous layers.
The viscous CFL number is set to mode/1000.
Limited testing indicates that a value of mode = 50 is
stable and increases the time step near the wall by at least an
order of magnitude.
(tdtmst)
|
| 136 | | Divergence checker, mode =
n1 + 10 n2, where
(lpgrp)
| |
n1 | |
Divergence Definition
|
|---|
|
1 | | Max residual > 1.0,
L2 norm increasing
| |
2 | | Max residual > 5.0,
L2 norm increasing
| |
3 | | Max residual > 10.0,
L2 norm increasing
|
and
| |
n2 | |
Action Taken When Diverging
|
|---|
|
1 | | Terminate iteration for
current cycle
| |
2 | | Abort run
| |
3 | | Reduce CFL number by 1/2 and
terminate iteration for current cycle
|
|
| 137 | | Butt line interpolation region for
USERSPEC;
smear USERSPEC conditions over
0.001 × butt line at minimum and maximum butt line
(uspeci)
| |
mode | |
Result
|
|---|
|
0 | | No interpolation
| |
n | | n = 0.001 ×
butt line for interpolation
|
|
| 138 | | Use large cell
Jacobians at boundaries
(BC_bcwall, chrhsv, mphzmet, nsrhsv, tdarea, tdarea1,
tdarea2, PreWallBC, US_tdbcg, vismet)
| |
mode | |
Result
|
|---|
|
≤ 1 | | Use large cells
| |
2 | | Use large cells,
central difference Jacobian
| |
3 | | Throw out half cell at
boundaries
| |
5 | | Solve
∂P/∂n equation at walls
|
|
| 139 | | Turn on
grid-based flux limiter
(tdup1)
|
| 140 | | Use
first-order differencing for computing
∂(u,v,w) / ∂ξ
term in vorticity used in turbulence models
(vortcy)
|
| 141 | | 2nd-order wall
boundary conditions (explicit)
(BC_bcvel, BC_bcwall, US_tdbcg)
| |
mode | |
Result
|
|---|
|
1 | | Second order
∂p/∂n, ∂T/∂n, and
∂utan/∂n
| |
2 | | Second order
∂p/∂n and ∂T/∂n, but not
∂utan/∂n
|
|
| 150 | | Singular axis
on symmetry planes.
When symmetry plane test fails, zero this component of velocity.
(bcsing)
| |
mode | |
Result
|
|---|
|
1 | | u = 0
| |
2 | | v = 0
| |
3 | | w = 0
| |
4 | | do not zero any component
(use average)
| |
5 | | v = w = 0
| |
6 | | u = w = 0
| |
7 | | u = v = 0
|
|
| 151 | | For a singular axis, the value on the
axis is a radius-weighted average of the values two points away
from the axis, instead of the values at the adjacent points;
values at the adjacent points are set to the average of the axis
value and the value two points away from the axis.
E.g., for a singular axis at j = 1, with the
k direction "circumferential", the value on the axis is a
radius-weighted average of the values at j = 3,
instead of at j = 2.
Then for each k, the value at j = 2 is
set to the average of the values at j = 1 and
j = 3.
This test option overrides TEST 121.
(bcpinw, bcsing, kebc, kerot, radavg, relsng, sabound,
sngthrm, sstbound, vprtbound, vsctbound)
|
| 153 | | Iteration
frequency for updating pressure at outflow boundaries.
The default is 5.
(BC_IntgrtBnd)
|
| 154 | | When computing values at "undefined"
boundary points, and no neighboring non-hole, non-fringe, points
are found, average over neighboring fringe points
(bcundef, kebc, sabound, sstbound, vprtbound, vsctbound)
|
| 155 | | For unstructured
grids, extrapolate at freestream characteristic boundaries.
(BC_bcfree)
|
| 157 | |
USERSPEC inflow
(uspeci)
| |
mode | |
Result
|
|---|
|
1 | | Put USERSPEC inflow
at all points in the buttline range.
Do not check for above/below vehicle.
| |
2 | | Same as mode 1, but also
ignore fuselage station check.
|
|
| 158 | | Write various unstructured grid info to
.lis file
(opngrd)
|
| 160 | | Pressure correction factor = mode / 1000,
for specified mass flow boundary
(BC_pdsmfr)
|
| 162 | | Cebeci-Smith boundary layer edge definition
(cebeci)
| |
mode | |
Result
|
|---|
|
0 | | 1.0% change in
Utotal
| |
1 | | 0.995 Ht
| |
2 | | 0.99 Utotal
| |
3 | | 0.9999 Utotal
|
|
| 163 | | Criteria for defining Fmax
in Baldwin-Lomax model
| |
mode | |
Result
|
|---|
|
> 0 | | Search outward from
wall; use first peak in F that is followed by a
fractional decrease in F of mode / 1000.
| |
< 0 | | Use the max F
value between the wall and the | mode |'th grid
point from the wall
|
The default value is +700.
(blomax)
|
| 164 | | Iteration interval for updating gas
properties and species for ireal = 2 (i.e.,
equilibrium air).
The default value is 1.
(tdgas)
|
| 165 | | Sets the
tolerance for the distance between grid points used to determine
a singular direction to 10mode/1000.
The default is a tolerance of 10−8 (i.e.,
mode = −8000).
(bcsing, direct)
|
| 168 | | In a marching solution with the algebraic
turbulence models, begin turbulent flow at streamwise station
mode.
(tdvis1)
|
| 170 | | NASA Ames time step formula.
(CFL increases as 1/(Δy)1/2 near the wall.
Thus, Δt decreases as (Δy)1/2,
not Δy as the default.)
CA is scalar coefficient on CFL; i.e.,
CFLwall is proportional to CA.
This test option has no effect if
TEST 105 mode 1, 2, or 3 is used.
(tdtmst)
|
| 172 | | Turn off base energy for ireal = 3
(i.e., frozen and finite-rate chemistry).
(aichem, aijkrg, aikeps, aixyzrg, BC_bcfreechar,
BC_ijkfreecha, BC_usrffreechar, dqlim1, gas1, gas2, gas3, gas4,
gasint, genturb, redimchem, stomp, tdimfp, tdroe4, tdroe5,
turbupd, US_DqLimit, US_EnsightGama, US_EnsightMach,
uspeci)
|
| 174 | | For the algebraic turbulence models, the
iteration interval for updating the turbulent viscosity.
The default is 1.
(tdvis1)
|
| 175 | | Boundary conditions at freestream radial
outer boundaries (nzn = −6)
(BC_bcfree, BC_bcfreechar, tdbcgs)
| |
mode | |
Result
|
|---|
|
0 | | Compute characteristics
from freestream conditions
| |
1 | | Compute characteristics
from conditions at i = 1
along boundary
| |
2 | | Extrapolate without testing
at k boundaries; treat i and j
boundaries as in mode 1
|
|
| 177 | | Freeze maximum residual
(lpgrp)
|
| 178 | | Insert
freestream species buffer during
BLOW VALVE relaxation
(BC_bcbled)
|
| 179 | | Modify solidbody rotation radius to get
linear swirl profile in r, but with zero velocity not at
center of rotation.
mode = 1000 r0, where r0
is the radius (from the point xc,
yc, zc specified using the
SOLIDBODY keyword in the
ARBITRARY INFLOW
keyword block) for zero velocity.
(rotat)
|
| 180 | | Defines the
radius of the solidbody rotation core for actuator disk free
vortex modeling.
mode = 1000 rcore,
where rcore is the solidbody core radius.
(rotat1)
|
| 182 | | For Roe coupling, don't modify boundary
states at interior face points for consistency with boundary values
(tdbcgs)
|
| 185 | | ??
(gasint)
|
| 187 | | mode / 1000 = factor for suppression
of streamwise pressure gradient in a marching solution.
The default is 950, and values below 800 are not recommended.
When separation or strong adverse pressure gradients are
causing problems, values between 800 and 900 will really help.
(rhsmar, tdimfp)
|
| 189 | | If
a first-order upwind explicit operator modified for stretched
grids is used (e.g.,
RHS ROE FIRST PHYSICAL),
then TEST 189 1 must also be specified.
(numprtinp)
|
| 190 | | Outgoing wave Roe boundary treatment
(pstinp, roecof, roeht, tdbcgs)
| |
mode | |
Result
|
|---|
|
0, 1 | | Use normal Roe
boundary treatment (uses boundary point in formulation)
| |
2 | | Lower order by one (does
not use boundary point in formulation).
This option cannot be used with
TVD in the same zone.
| |
3 | | Use zero-order extrapolation
|
|
| 192 | | Save metrics in a temporary file.
After the first cycle, metrics will be read rather than computed.
This eliminates the CPU resources required to re-compute the
metrics each cycle, but adds significant I/O to the computation.
In the past, on at least some Cray systems, this reduced the CPU
time by approximately 40.8 micro-seconds / (node-cycle).
On the more common platforms, however, it is generally faster to
re-compute the metrics rather than store them.
(lpcycl)
|
| 193 | | Do not stop if a singular line
is encountered normal to a wall
(bbdamp, blinit, kepy2)
|
| 194 | | Bypass singular viscous metric check
(dsolv, mphzmet, vismet)
|
| 195 | | When using
BLOW SURFACE, print a
warning when the flowfield static pressure becomes larger than
the plenum total pressure, causing blowing to be turned off
at that point.
Note that this is a five-line message written for each iteration
and each "closed" node, and could cause the .lis file to
become very large very quickly.
(BC_bcbled)
|
| 196 | | Overlapping grid: print an error message
if there are no points adjacent to a fringe point.
(BC_norot, srfpar, tdbcgs)
|
| 197 | | Roe self-coupling mode
(pstinp)
| |
mode | |
Result
|
|---|
|
0 | | Once per iteration, using
bcself
| |
1 | | Once per cycle, using
standard zone coupling
|
|
| 199 | | Singular axis
averaging; average from 1 to (max − 1), not
1 to max.
(bcsing, linzero, radavg, relsng, sngthrm)
|
| 200 | | Don't bomb for negative speed of sound in
tdroe3A
(tdroe3A)
|
