/ RE = 200
/
/          CONVERSION OF NEUTRAL FILE TO FIDAP Database
/
FICONV( NEUTRAL )
INPUT( FILE="cylind.FDNEUT" )
OUTPUT( DELETE )
END
/
/
TITLE
cylind
/
FIPREP(  )
/
EXECUTION( ADD, NEWJ )
/
PROBLEM( ADD, 2-D, INCO, TRAN, LAMI, NONL, NEWT, MOME, ISOT, FIXE, SING )
/
TIMEINTEGRATION( ADD, BACK, NSTE = 801, TSTA = 0, DT = 0.20, FIXE )
SOLUTION( S.S. = 50, VELCONV = .001, RESCONV = .001)
/
POSTPROCESS(NBLOCKS=1)
1 801 4
/
DATAPRINT( ADD, CONT )
PRINTOUT( ADD, NONE, BOUN )
/
ENTITY( ADD, NAME = "fluid", FLUI )
ENTITY( ADD, NAME = "bottom", PLOT )
ENTITY( ADD, NAME = "top", PLOT )
ENTITY( ADD, NAME = "inlet", PLOT )
ENTITY( ADD, NAME = "outlet", PLOT )
ENTITY( ADD, NAME = "cylinder", PLOT )
/
DENSITY( ADD, SET = 1, CONS = 1.0)
VISCOSITY( ADD, SET = 1, CONS = 0.005 )
/
/ ICNODE( ADD, VELO, ZERO, ENTI = "fluid", X, Y, Z )
/ To read steady-state flow field as the initial condition:
/ change *.FDPOST to *.FDREST
/ and then
/ For the 1st run
ICNO (VELO, ZERO)
/ For 2nd and third runs
/ ICNO (VELO, READ)
/
BCNODE( ADD, VELO, ENTI = "cylinder", ZERO, X, Y, Z )
BCNODE( ADD, UY, ENTI = "inlet", ZERO )
BCNODE( ADD, UX, ENTI = "inlet", CONS = 1.00 )
BCNODE( ADD, UY, ENTI = "bottom", ZERO )
BCNODE( ADD, UY, ENTI = "top", ZERO )
/
END(  )