/ *****************************************************************
/ Disclaimer: This file was written by GAMBIT and contains
/ all the continuum and boundary entities and coordinate systems
/ defined in GAMBIT. Additionally, some frequently used FIPREP
/ commands are added. Modify/Add/Uncommment any necessary commands.
/ Refer to FIPREP documentation for complete listing of commands.
/ *****************************************************************
/
/          CONVERSION OF NEUTRAL FILE TO FIDAP Database
/
FICONV( NEUTRAL )
INPUT( FILE="step.FDNEUT" )
OUTPUT( DELETE )
END
/
TITLE
step
/
FIPREP
/
/          PROBLEM SETUP
/
PROBLEM (2-D, TURB, NONLINEAR, ISOTHERMAL)
EXECUTION( NEWJOB )
PRINTOUT( NONE )
DATAPRINT( NODES=3, ELEM=1 )
/ Control the echoing of the input data
/ No printout of node data
/ Print element group control data only
/
/          CONTINUUM ENTITIES
/
ENTITY ( NAME = "fluid", FLUID, PROPERTY = "fluid" )
/
/          BOUNDARY ENTITIES
/
ENTITY ( NAME = "inlet", PLOT )
ENTITY ( NAME = "wall1", WALL )
ENTITY ( NAME = "wall2", WALL )
ENTITY ( NAME = "wall3", WALL )
ENTITY ( NAME = "wall4", WALL )
ENTITY ( NAME = "outlet", PLOT )
/
/          SOLUTION PARAMETERS
/
SOLUTION( S.S. = 100, VELCONV = .001, RESCONV = .001)
RELAXATION
0.5 0.5 0 0 0 0 0.5 0.5
/ RElaxation factors for
/ u1,u2,u3,p,T,s,K,Epsilon
PRESSURE( MIXE, DISCONTINUOUS )
OPTIONS( UPWINDING)
/SCALE( VALUE = 1 )
/TIMEINTEGRATION( BACKWARD,NSTEPS = ,DT = 0,,, )
/POSTPROCESS( NBLOCKS =  )
/
/          MATERIAL PROPERTIES
/
/ Partial list of Material Properties data
/
DENSITY( SET = "fluid", CONSTANT = 1 )
/
VISCOSITY( SET = "fluid", CONSTANT = 2.2222e-5 , TWO-)
/
EDDYVISCOSITY(BOUSSINESQ)
/ Standard
TURBOPTION(ADD,STANDARD)
    1.44,     1.92,     0.09,     1.00,     1.30,     0.90,     0.41,     26,     0.80,     0.9,     0.9,
    0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,
    0.9,     0.9,     0.9,     0
/ Extended, Original
/TURBOPTION(ADD,EXTE)
/   1.15,     1.9,     0.09,     0.75,     1.15,     0.90,     0.41,     26,     0.80,     0.9,     0.9,
/   0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,
/   0.9,     0.9,     0.9,     0.25
/ Extended, revised
/TURBOPTION(ADD,EXTE)
/   1.35,     1.9,     0.09,     0.75,     1.15,     0.90,     0.41,     26,     0.80,     0.9,     0.9,
/   0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,
/   0.9,     0.9,     0.9,     0.05
/ RNG Original
/TURBOPTION(ADD,RNG)
/   1.42,     1.68,     0.085,     0.7179,     0.7179,     0.90,     0.3875,     26,     0.80,     0.9,     0.9,
/   0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,
/   0.9,     0.9,     0.9,     0
/ RNG revised
/TURBOPTION(ADD,RNG)
/   1.45,     1.83,     0.0865,     0.80,     1.15,     0.90,     0.41,     26,     0.80,     0.9,     0.9,
/   0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,     0.9,
/   0.9,     0.9,     0.9,     0
/ RE number based on step height = 45,000
/
/          INITIAL AND BOUNDARY CONDITIONS
/
BCNODE( UX, CONSTANT = 1.0, ENTITY = "inlet" )
BCNODE( UY, CONSTANT = 0.0, ENTITY = "inlet" )
$energy = 0.003*1.0^2
BCNODE( KINE, CONSTANT = $energy, ENTITY = "inlet" )
$dissp = 0.007*1.0^3/2.0
BCNODE( DISS, CONSTANT = $dissp, ENTITY = "inlet" )
BCNODE( VELO, CONSTANT = 0, ENTITY = "wall1" )
BCNODE( VELO, CONSTANT = 0, ENTITY = "wall2" )
BCNODE( VELO, CONSTANT = 0, ENTITY = "wall3" )
BCNODE( VELO, CONSTANT = 0, ENTITY = "wall4" )
/
ICNODE( VELO, CONSTANT = 0, ALL )
ICNODE( KINE, CONSTANT = $energy, ALL )
ICNODE( DISS, CONSTANT = $dissp, ALL )
/
RENUMBER
/ Turn on  node number optimal renumbering option
END
/