/ *****************************************************************
/ 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="atube.FDNEUT" )
OUTPUT( DELETE )
END
/
TITLE
atube
/
$Re = 100
$vis = 1.0/$Re
FIPREP(  )
DATAPRINT( ADD, CONT )
PROBLEM( ADD, AXI-, INCO, STEA, LAMI, NONL, NEWT, MOME, ENER, FIXE, SING )
SOLUTION( ADD, S.S. = 100, VELC = 0.0001, RESC = 0.0001, ACCF = 0 )
OPTIONS(UPWINDING,NODISSIPATION)
UPWINDING
1 1 0 0 2 0 1 1
PRINTOUT( ADD, ALL, BOUN )
EXECUTION( ADD, NEWJ )
ENTITY( ADD, NAME = "fluid", FLUI )
ENTITY( ADD, NAME = "wall", PLOT )
ENTITY( ADD, NAME = "inlet", PLOT )
ENTITY( ADD, NAME = "outlet", PLOT )
ENTITY( ADD, NAME = "center", PLOT )
BCNODE( ADD, VELO, ENTI = "wall", ZERO )
BCNODE( ADD, UZC, ENTI = "inlet", POLY = 1)
    2,    -8,     0,     2,     0
BCNODE( ADD, URC, ENTI = "inlet", ZERO )
BCNODE( ADD, URC, ENTI = "center", ZERO )
BCNODE( ADD, TEMP, ENTI = "inlet", ZERO )
BCFLUX( ADD, HEAT, ENTI = "wall", CONS = 1 )
/ Set the proper outlet heat flux by energy balance = 4q_s/(D rho c_p U_m)
BCFLUX( ADD, HEAT, ENTI = "outlet", CONS = 4 )
DENSITY( ADD, SET = 1, CONS = 1 )
VISCOSITY( ADD, SET = 1, CONS = $vis )
CONDUCTIVITY( ADD, SET = 1, CONS = 1 )
SPECIFICHEAT( ADD, SET = 1, CONS = 1 )
END(  )