sajben_s2u transonic diffuser
Unstructured hexahedral grid
Run 2 - "CONVERGE LOADS" strategy


freestream total 0.90 19.58  500.    0. 0.

cycles 20000
iter_cycle 5 print frequency 5

CFL 10.


/ Use the CONVERGE LOADS strategy: Wind-US monitors convergence of loads,
/ and automatically adjusts SECOND and SMLIMT when changes in forces and
/ moments are smaller than VAL (set to 0.2 here)

CONVERGE LOAD FREQUENCY 10 0.2


/  Spatial operator info.  2nd-order Rusanov is recommended
/  for unstructured grids.

RHS Rusanov Second CELLCENTERED


/  Temporal scheme.  Use Gauss-Seidel implicit operator.
/  Save the Jacobians in memory rather than recompute at each iteration
/  (which would reduce memory usage by about 20%, but adds calculations).
/  Solve as a diagonal matrix system (not full block, which is more stable 
/  but is also about 10% slower).


IMPLICIT UGAUSS SAVE_JACOBIAN DIAGONAL SUBITERATIONS 4


/  Type of time step: local (varies in space), recommended for steady
/  flows, as opposed to global (constant) for time accurate flows.
/  But limit the maximum to minimum step size to DTRATIO.  Use CFL number
/  given above.

TIMEMARCHING VARIABLE-DT DTRATIO 1000.0


/  Whenever the skewness angle where 2 adjacent cells meet is greater
/  than the specified number of degrees, reduce order of RHS to 1st.

GRID LIMITER ANGLE 150


/  Inflow BC: impose the total pressure (psi) given below.
/  Since HOLD_TOTALS is also there, then, for subsonic inflow,
/  hold the total pressure (psi) and temperature (R) given below
/  and whatever flow angles (degrees) are initially on the inflow
/  planes, and extrapolate Mach number.  The Mach no. below will
/  only be used for initializing the flow if there's no CFL file.
/  Each USURFACE, also called a "USURF number", also called a
/  "surface segment", gets its own specifications.

ARBITRARY INFLOW
TOTAL
HOLD_TOTALS
ZONE 1
/             usurf#  Mach   Ptot   Ttot  alpha beta
   USURFACE       1   0.90  19.580  500.0    0. 0.
ENDINFLOW


/  Outflow BC: impose the given static pressure (psi).

DOWNSTREAM PRESSURE 14.1 ZONE 1


/  Turbulence model.  Spalart's default free-stream eddy viscosity nu_t
/  is 5.0, but could be changed via the FREE_ANUT keyword.  The Spalart
/  model is only used within a certain distance of a wall that defaults
/  to 27 times what the Prandtl turbulent boundary layer thickness would be
/  at x=Lref.  You can change this distance by using the MAX_WALL_DISTANCE
/  keyword and specifying the parameter DistMax in whatever length units
/  the grid coordinates are in.  Wall functions are available for both
/  structured and unstructured solvers.

TURBULENCE MODEL SPALART
/FREE_ANUT 5.0
/MAX_WALL_DISTANCE 30
/WALL FUNCTION ON


/  Debug inputs are specified as "DEBUG NUMBER MODE".  Must have some
/  MODE, say 1, even if not used.

/  Save wall distances in the CFL file

Debug  2 1

/  Smooth the initial and restart conditions,
/  using MODE-many smoothing passes

Debug 13 1

/  Include the diagonal viscous flux Jacobian terms.
/  Important for y+ < 0.5.

Debug 18 1

/  Convergence criteron: assume convergence when forces
/  and moments change less than 1/MODE.  Needed for the
/  CONVERGE LOADS strategy.

Debug 24 50


/  Print loads - needed for the CONVERGE LOADS strategy

loads
   print totals drag frequency 10
   reference area 144.0
   reference length 12.0
   reference moment center 0.0 0.0 0.0
   zone 1
   surface u 3 force viscous
endloads


END