Test Options

Several user-controlled options have been provided as an aid to modifying the WIND code. 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

   TEST number mode

If mode is not described for a given test option, the user should use mode = 1 to activate that option.

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 the WIND code, 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
1		Memory size check. Stop the run after displaying the memory requirements but before memory is allocated. (mpinit, mpterm, opngrd, wind, zinit)     mode    Result n < 0 Retrieve dimensions from the actual grid (.cgd) file n > 0 Retrieve dimensions from a text file with the following format: #zones idim1 jdim1 kdim1 nfrpts1 ... idimn jdimn kdimn nfrptsn Each of the nfrpts values is 0 except for overlapping grids. This mode allows one to determine how much memory WIND will require before generating the grid. Note that the name of the text file should end in .cgd, like the actual grid file.
2		Designed for parallel processing data transfer debugging     mode    Result 1 Not used 2 Don't read zonal boundary data (evrwbd) 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, rwstat, rwd, rwi, rwr) 4 Trace network traffic (psexit, psrwev, psrwgv, psrwnh, psrwd, psrwi, psrwr) 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.
4		Reserved for internal Boeing use in routines for hybrid unstructured grid capability (opngrd)
5		Compute the L2 norm of the residual for the Navier-Stokes equations in the same way as the "old" code (l2norm)
6		Write .cfl file compatible with the "old" code (asnsx, asvisc, mpini2)
7		Do not use high performance C I/O interface (openf, rwstat)
8		Use Version 2 common files (mpinit, openf, rwstat)
10		Print the time step information (i.e., minimum/maximum CFL and Delta t) into the .lis file every mode cycles, instead of just on the first cycle. (zsolv)
17		Use "new" Baldwin-Barth turbulence model (bbarth)
20		Non-dimensionalize k and omega in the SST model the "old" way (aijkrg, aikeps, sst, sstprtinp, sstpstprt)
21		Spalart-Allmaras turbulence model. (See the TURBULENCE keyword.) (redimsa, sabound, sapreinp, sinut, spalart)     mode    Result 0 Use original 1992 model, with an ft2 term for laminar stabilization and a default freestream value for nut 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 nut of 0.1. In addition, the initial value of the dependent variable is set to the freestream nut. 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 Chien k-epsilon models (bbdamp, blomax, cebeci, kepy2)
29		For the Cebeci-Smith model, use y+ based on wall vorticity. (algtur, cebeci) For the Chien k-epsilon model, TEST 29 is a production limiter (kelhssch, keprod, kerhssch)     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
30		Irrotational boundary condition at freestream inflow boundaries (bcfree)
40		Old switch for compressibility correction in k-epsilon 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 muturb = 10 mulam along the wall. (sstbound)
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) (dampi, dampj, dampk, bdload)     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 muT, such that the maximum value of muT / (muL)inf = mode × 1000. Suggested range is 50 < mode < 100. Do not use this option with the Chien k-epsilon model; use K-E MAXIMUM TURBULENT VISCOSITY instead. (keppstin, mutlim)
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 (bcbled)
53		Allows no-slip in inviscid flow at "viscous wall" boundaries (bcwall)
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, sst)
58		Store the Lorentz force in the .cfl file instead of the electric field (emdef)
59		Apply SMOOTHING keyword values to the Spalart-Allmaras model as well as the mean flow solver. (Currently deactivated.) (spalart)
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)
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. (tdarea)
64		Remove dt from dq when computing residuals (l2norm)
65		In marching solutions, lower the CFL number for the last marching step (zsolv)
66		Don't update beta in gas3, for ireal = 2 and ispec = 2. This test option is not recommended but will decrease run time. (gas3)
67		When solving turbulence model equations, treat bleed boundaries as slip walls instead of no-slip walls. This was the default behavior prior to WIND 5.101. (kebc, sabound, sstbound)
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, rho, or T in the gas routines. Tolerance level = (0.1)mode. (gas1, gas2, gas3, gas4)
71		Apply curve fit equations for thermodynamic properties even when the temperature is outside the applicable range. If mode = 1, an error message will be written whenever this occurs. For other non-zero values of mode, no error message is written. If this test option is not set, the run will be terminated if the temperature is outside the range of the curve fits. (cpfun, gibfn, hfun, proper)
75		?? (prpold, tdimafk, tdimfu)
77		Set bleed limits for a specific F/A-18 15% scale sting and distortion model, for use with the BLEED AEDC keyword (bcaedc)     mode    Result 0 F/A-18E 1 F/A-18A
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 (bcfree)
88		Bypass negative T check in tdgas1. 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, tdgas1)
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 (chprtinp, therm1, thermp)     mode    Result 1 beta = gamma = betainf 2 beta = gamma = 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, prpold, rhssrc)
93		Turn on/off chemistry species diffusion (chrhsv, rhsvfl, tdutv1)     mode    Result 0 Include all chemistry species diffusion terms 1 Include chemistry species diffusion terms, except for the diffusion gradient in the conduction term 2 Same as 0 3 Ignore chemistry species diffusion
94		Turn off implicit chemistry terms (tdimpl, tdutaa)
95		Turn off chemistry source term (rates, rates1, rates2, rates3, rates4, ratesa, ratesadl, ratesb, ratesf)
96		Apply chemistry source term over mode iterations for finite-rate non-equilibrium chemistry (chrhss)
97		P. D. Thomas turbulence model scanning direction. By default, WIND 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
98		mode is the relaxation factor for the species mass fraction from the Liu and Vinokur (real gas) model. By increasing mode, the mass fraction relaxes faster. Default is mode = 5. mode = 1 corresponds to instantaneous. (gas1, gas2, gas3, gas4)
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 (bcfree)     mode    Result 0 Second order, with limit of Delta Q <= Q/2 1 1st-order, original characteristic treatment 2 2nd-order, original characteristic treatment 3 1st-order, Roe's average characteristic treatment 4 2nd-order, Roe's average characteristic treatment
102		Use time-averaged back pressure for mass flow boundary condition (pdsmfr)
104		Turn off implicit viscous terms (tdutv1)
105		Time step type (prtinp, tdtmst)     mode    Time Step Type 0 Delta t = CFL / max (lambdaxi, lambdaeta, lambdazeta) 1 Flow angle scaling, Delta t = CFL × (fxi Delta xi + feta Delta eta + fzeta Delta zeta), where fxi = [1 + tan (theta) + tan (psi)]1/2 feta = fxi tan (theta) fzeta = fxi tan (psi) 2 Velocity scaling, Delta t = CFL × min (fxi Delta xi, feta Delta eta, fzeta Delta zeta), where fxi = u / uxi / | uxi + c | feta = u / ueta / | ueta + c | fzeta = u / uzeta / | uzeta + c | 3 Delta t = CFL × min (Delta xi, Delta eta, Delta zeta) / (| u | + c) 4 Delta t = CFL / (lambdaxi + lambdaeta + lambdazeta)
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 (bcfree)     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 infinity 1 2.0 2 1.5 3 1.33 4 1.1
111		Singular matrix check (jacpr4, jacpr5, jacprg, tdsol4, tdsol9, tdsol11, tdsolg, 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 (bcconf, bcfree)     mode    Result 0, 1 Print a warning message and continue 2 Print an error message and stop
114		Central difference zeta 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 (bcwall)
116		Set inward pointing normal to zero in tdbcm1 at unknown grid topology points. (tdbcm1)
117		Freeze inflow boundaries, even in subsonic flow (bcfree)     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 (tdarea)
123		Track back pressure, mass flow, and integrated total pressure with outflow boundary conditions (bcconf, pdsmfr)
124		Write convergence information to list output (.lis) file every iteration instead of every cycle (lpschm)
126		Compressibility correction to Baldwin-Lomax turbulence model (blomax)     mode    Result 0 No compressibility correction 1 kappa = 0.0180 for Baldwin-Lomax model (CFL3D uses this)
127		Scale printed residual by maximum residual over all time steps (lpschm)
128		Check the L2 norm of the residual for convergence instead of the maximum residual (l2norm)
130		Roe scheme physical space extrapolation (bcfree)
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., (Delta t)j = (Delta t)mode for j < mode. (tdtmst)
132		Renormalize, changing from total to static values. Normalizing values in the .cfl file are unchanged. (redim)
133		Print multi-grid sub-iteration convergence data (currently deactivated) (lpschm)
134		Second order characteristic extrapolation for adjacent conditions (bcfree)     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 4 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 (bcwall, chrhsv, nsrhsv, tdarea, vismet)     mode    Result <= 1 Use large cells 2 Use large cells, central difference Jacobian 3 Throw out half cell at boundaries 5 Solve dP / dn equation at walls
139		Turn on grid-based flux limiter (tdup1)
140		Use first-order differencing for computing d(u,v,w)/dxi term in vorticity used in turbulence models (vortcy)
141		2nd-order wall boundary conditions (explicit) (bcwall)     mode    Result 1 Second order dp/dn, dT/dn, and dutan/dn 2 Second order dp/dn and dT/dn, but not dutan/dn
147		?? (currently deactivated) (turpdt)
148		?? (currently deactivated) (turpdt)
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)
153		Iteration frequency for updating pressure at outflow boundaries. The default is 5. (bcconf)
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)
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		Reserved for internal Boeing use in routines for hybrid unstructured grid capability (opngrd, pstinp)
160		Pressure correction factor = mode / 1000, for specified mass flow boundary (pdsmfr)
162		Cebeci-Smith boundary layer edge definition (cebeci)     mode    Result 0 1.0% change in Utotal 1 0.995Ht 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 = 1 and 2. The default value is 1. (tdgas1)
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		For the algebraic turbulence models, begin turbulent flow at i = mode (tdvis1)
170		NASA Ames time step formula. (CFL increases as 1 / (Delta y)1/2 near the wall. Thus, Delta t decreases as (Delta y)1/2, not Delta 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)     mode    CA 1 0.01 n 0.001n
172		Turn off base energy for ireal = 1. (aichem, aijkrg, gas1, gas2, gas3, gas4, gasint, pstinp, stomp, tdimfp, tdroe4, 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) (bcfree)     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 (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)
181		?? (evrwbd, wbnd2)
182		?? (tdbcgs)
185		?? (gasint, mpinit)
187		mode / 1000 = factor for suppression of streamwise pressure gradient when marching. 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. (prtinp)
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
191		Options used with the BLOW and BLOW SURFACE keywords.     mode    Result 0 Angles are defined with respect to the surface normal. 1 Angles are defined with respect to the projection of the surface normal onto a constant z plane. 2 Angles are defined as in mode 0, and equilibrium air chemistry is allowed. The surface values of the effective gamma and compressibility factor are approximated by using the values at the adjacent flowfield point. 3 Angles are defined as in mode 1, and equilibrium air chemistry is allowed as in mode 2. 4 Backward compatibility mode for the BLOW keyword. With this mode, the syntax is "BLOW ibreg mdot Tinj angle", where angle is the blowing angle relative to the x-y plane, in degrees. Note, though, that this mode is not recommended, as it does not provide the requested mass flux for surfaces whose normal vector contains components in the z direction (i.e., surfaces with transverse curvature). Modes 0 and 1 apply to both the BLOW and BLOW SURFACE forms; modes 2 and 3 only apply to the BLOW SURFACE form; and mode 4 only applies to the basic BLOW form. (bcbled)
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, 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. (bcbled)
196		Overlapping grid: print an error message if there are no points adjacent to a fringe point. (inttur, tdbcgs)
197		Roe self-coupling mode (pstinp)     mode    Result 0 Once per iteration, using bcself 1 Once per cycle, using standard zone coupling
198		Sets the tolerance level for fuselage station USERSPEC inflow to mode/1000. The default tolerance is 0.01 (i.e., mode = 10). (Currently deactivated) (uspeci)
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 tdroe3 (tdroe3)