shockwave.afl

$GLOBAL,
	TITLE    = 'shockwave',
	SIMTYPE  = 2,
	FOLDER   = 'shockwave',
	MACH     = 1.0,
	REYNOLDS = 300000.0,
	ORDER    = 1,
	$END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE TURBULENCE MODEL
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

COMMENT:~~~~~~~~~~~~SPALART-ALLMARAS~~~~~~~~~~~~~~~
$TURB,
    TURBTYPE = 0,
    $END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE SPATIAL SCHEME
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

$SPATIAL,
	SCHEME = 5,	-999,	-999,
	METRIC = 0,
	VISCOUS = 0,  	0,	0,
	VISCHEME = 0,
	VISCMON = 0,
	CUTOFFI = 0.00001,
	ISOTROPI = 1,
	ICUT = 0, 0, 0,
	SUTHERLN = 0.38,
	$END

SPATIAL,
    SCHEME   = 33, -999, -999,
    METRIC   = 13,
    VISCOUS  = 0,  0,  0,
    VISCHEME = 0,
    VISCMON  = 0,
    SUTHERLN = 0.38,
    $END	
	
SPATIAL,
    SCHEME   = 34, -999, -999,
    METRIC   = 13,
    VISCOUS  = 0,  0,  0,
    VISCHEME = 0,
    VISCMON  = 0,
    SUTHERLN = 0.38,
    $END	
	
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE TIME INTEGRATION SCHEME
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
$TIMESTEP,
	MAXITER  = 178,
	SCHEME   = 'RK4',
	KTVDRK   = 0,
	GLOBALDT = 1.E-3,
	SUBDT    = 1.E-3,
	SUBON    = 1,
	MAXSUB   = 2,
	$END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE DAMPING/FILTERING
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

 $DAMPING,
        TYPE =            3 ,
        ES2 =   0.2000000     ,
        ES4 =   0.02 ,
        FES2 =    10.000000     ,
        FES4 =    2.000000     ,
        MAXRED =            1 ,
        OMGAV =   0.5 ,
        SRCONST =   0.0000000E+00 ,
        FILTER =            0 ,           0 ,           0 ,
        NFILTER =            1 ,           0 ,           0 ,
        ALPHAN =   0.45000     ,  0.45000     ,  0.45000     ,
        ALPHA1 =   0.0E+00     ,  0.0E+00     ,  0.0E+00 ,
        ALPHA2 =   0.49000     ,  0.49000     ,  0.49000     ,
        ALPHA3 =   0.45000     ,  0.45000     ,  0.45000     ,
        ALPHA4 =   0.45000     ,  0.45000     ,  0.45000     ,
        ALPHA5 =   0.45000     ,  0.45000     ,  0.45000     ,
        ORDERN =            8 ,           8 ,           8 ,
        ORDER1 =            0 ,           0 ,           0 ,
        ORDER2 =            2 ,           2 ,           2 ,
        ORDER3 =            4 ,           4 ,           4 ,
        ORDER4 =            6 ,           6 ,           6 ,
        ORDER5 =            8 ,           8 ,           8 ,
      $END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    OUTPUT AND DEBUG SETUP
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

$OUTPUT,
	FORMAT = 'TECP',
	PRINTFRQ = 50,
    TECANIM  = 1,
    TRANGE   = 0,1000000,2,
	$END

$DEBUG,
	POINT = 4, 4, 2,
	$END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:    DEFINE GRID
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
$BLOCK,
	NAME     = 'block01',
	MESHFILE = 'shock.in',
	FILETYPE = 'AEROFLO',
	ORDER    = 1,
	ISIZE    = 192,
	JSIZE    = 3,
	KSIZE    = 3,
	PERIODIC = 0, 1, 1,
	$END

COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
COMMENT:      DEFINE BOUNDARY CONDITIONS
COMMENT:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

COMMENT:~~~~~~~~~~~~I-DIRECTION~~~~~~~~~~~~~~~~~~~~

comment:	Outflow condition at i= IE
$FACEBC,
	BLOCK = 'block01',
	TYPE = 31,
	SIDE = 2,
	$END