/ccoat_p.FDREAD
/This is a parameterized file for a curtain coating problem.
/This file both runs the simulation and postprocesses the results.
/To run this file type the command:
/fidap -id ccoat_p -in ccoat_p.FDREAD -new
/Note: The file ccoat.FDNEUT must also be in your working directory.
/

FICONV(NEUTRAL,NORESULTS,INPUT)
INPUT(FILE= "ccoat.FDNEUT")
END
/
/units are CGS
/
/scale geometry in GAMBIT so that characteristic length at outlet 
/  is .06 cm
/length units are cm
$length = .06
 
/scale factor
$scfac = $length
 
/density units are g/m3
$density = 1.25
 
/viscosity units are g/(cm-s)
$viscosity = 5.27
 
/surface tension units are g/s2
$surftens = 30.0

/velocity units of falling curtain and moving substrate are cm/s
$vel_c = 200
$vel_s = 200

/slip coefficient
$slipcoef = .01

/contact angle units are degrees
$contangle = 20
 
/Reynolds number and Capillary number at "outlet"
$Re = $density * $vel_s * $length / $viscosity
$Ca = $viscosity * $vel_s / $surftens
 
/Relaxation values for velocity, pressure and free surface d.o.f.
$vel_rel = .3
$pre_rel = .5
$sur_rel = .975

/title information
$fr = 100
$fc = 100
$Re1 = 100 * $Re
$Re2 = INT ($Re1)
$Re3 = ($Re2/$fr)

$Ca1 = 100 * $Ca
$Ca2 = INT ($Ca1)
$Ca3 = ($Ca2/$fc)

$head = "Curtain Coating:"
$rtitle = " Re = " + NTOS ($Re3)
$ctitle = " Ca = " + NTOS ($Ca3)

TITLE (String)
($head + $rtitle + $ctitle)

FIPREP
 PROB (2-D, STEA, INCO, LAMI, NONL, NEWT, MOME, ISOT, FREE, SING)
 PRES (MIXE, DISC)
 EXEC (NEWJ)
 SOLU (SEGR = 100, VELC = 0.1E-03, SURF = 0.1E-02, PREC = 21, KINE = 500)
 DATA (NONE)
 RELA (HYBR)
/ U        V       W  P        T S
 $vel_rel $vel_rel 0 $pre_rel  0 $sur_rel
 PRIN (NONE)
 SCAL (VALU = $scfac)                       

 ENTI (NAME = "fluid", FLUI, PROP = "fluid")
 ENTI (NAME = "inlet", PLOT)
 ENTI (NAME = "outlet", PLOT)
 ENTI (NAME = "slip", SLIP, PROP = "slip" )
 ENTI (NAME = "lfree",  SURF, DEPT = 0,  MAPP, X = 1.0)
 ENTI (NAME = "rfree1", SURF, DEPT = -1, MAPP, NORM)
 ENTI (NAME = "rfree2", SURF, DEPT = 0,  MAPP, NORM, ANG1 = 0.0)

 DENS (SET = "fluid", CONS = $density)
 VISC (SET = "fluid", CONS = $viscosity)
 SURF (SET = "fluid", CONS = $surftens)
 SLIP (SET = "slip",  CONS = $slipcoef, UT = $vel_s)

 BCNO (UT, ENTI = "slip",   CONS = $vel_s)
 BCNO (UN, ENTI = "slip",   ZERO)
 BCNO (UX, ENTI = "inlet",  ZERO)
 BCNO (UY, ENTI = "inlet",  CONS = -$vel_c)
 BCNO (UY, ENTI = "outlet", ZERO)
/ Specify the contact angle at node 14
 BCNO (CONT, CONS = $contangle, NODE = 14, X = 1.0, Y = 0.0)
/ The following lines release the constraints of UT on 4 selected nodes
/ so free slip condition will be applied instead.
 BCNO (UX, NODE = 14, FREE)
 BCNO (UX, NODE = 156, FREE)
 BCNO (UX, NODE = 136, FREE)
 BCNO (UX, NODE = 155, FREE)
/ The following 5 lines correct the local normals
 BCNO (COOR, NODE = 1)
 BCNO (COOR, NODE = 2)
 BCNO (COOR, NODE = 14)
 BCNO (COOR, NODE = 27)
 BCNO (COOR, NODE = 28)
/ Fix nodes 1 and 2 on the free surfaces
 BCNO (SURF, NODE = 1, ZERO)
 BCNO (SURF, NODE = 2, ZERO)

 ICNO (UY, CONS = -$vel_c, ENTI = "fluid")
END
CREATE(FIPREP,DELE)
CREATE(FISOLV)
PARAMETER(LIST)
RUN(FISOLV,FORE)