sajben_s2u transonic diffuser
Unstructured hexahedral grid
Run 1 - Manual smoothing and convergence monitoring


freestream total 0.90 19.58  500.    0. 0.

cycles 2000
iter_cycle 5 print frequency 5

CFL 10.


/  For the three smoothing parameters below:
/
/     SECOND really means "upwind flux dissipation".
/         0.5 is most stable, 0.1 is most accurate.
/   
/     FOURTH really means "Jacobian dissipation."
/         Always leave set to 1.5.
/
/     SMLIMT means TVD slope limiter.  1.0 is most
/         stable, 4.0 is most accurate.

SMOOTHING SECOND 0.10 FOURTH 1.5 SMLIMT 4.0


/  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

END