```
TALK=T;RUN(1,1)
DISPLAY
A fluid flows steadily in a direction normal to a
flat plate in channel. The Reynolds Number (i.e. velocity times
channel width divided by kinematic viscosity) equals 100.0.
For calculation deformation of plate is use "Bent of Solid Plate"\$
model.
ENDDIS
PHOTON USE
p;;;;

msg
msg P1 contours
msg
GR OU Z 1;
cont P1 z 1 fil;.0001
GR OU X 21 Y 1 10 Z 1 1 COL 1
pause

cl
msg
msg
msg Velocity vector
GR OU Z 1;
GR OU X 21 Y 1 10 Z 1 1 COL 1
vec z 1 SH DASH 0
pause

cl
msg
msg
msg Deflection plate vector
GR OU Z 1;
GR OU X 21 Y 1 10 Z 1 1 COL 1
SET VEC COMP - - -
SET VEC COMP STDP - -
SET VEC RE 0.02
VECT IZ 1 SH DASH 0
pause

ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************

TEXT(2D Flow Past A Bent Solid Plate; S401   )

************************************************************
************************************************************

IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 70
XULAST =14.
XFRAC(1)=0.014286 ;XFRAC(2)=0.028571
XFRAC(3)=0.042857 ;XFRAC(4)=0.057143
XFRAC(5)=0.071429 ;XFRAC(6)=0.085714
XFRAC(7)=0.1 ;XFRAC(8)=0.114286
XFRAC(9)=0.128571 ;XFRAC(10)=0.142857
XFRAC(11)=0.157143 ;XFRAC(12)=0.171429
XFRAC(13)=0.185714 ;XFRAC(14)=0.2
XFRAC(15)=0.214286 ;XFRAC(16)=0.228571
XFRAC(17)=0.242857 ;XFRAC(18)=0.257143
XFRAC(19)=0.271429 ;XFRAC(20)=0.285714
XFRAC(21)=0.3 ;XFRAC(22)=0.314286
XFRAC(23)=0.328571 ;XFRAC(24)=0.342857
XFRAC(25)=0.357143 ;XFRAC(26)=0.371429
XFRAC(27)=0.385714 ;XFRAC(28)=0.4
XFRAC(29)=0.414286 ;XFRAC(30)=0.428571
XFRAC(31)=0.442857 ;XFRAC(32)=0.457143
XFRAC(33)=0.471429 ;XFRAC(34)=0.485714
XFRAC(35)=0.5 ;XFRAC(36)=0.514286
XFRAC(37)=0.528571 ;XFRAC(38)=0.542857
XFRAC(39)=0.557143 ;XFRAC(40)=0.571429
XFRAC(41)=0.585714 ;XFRAC(42)=0.6
XFRAC(43)=0.614286 ;XFRAC(44)=0.628571
XFRAC(45)=0.642857 ;XFRAC(46)=0.657143
XFRAC(47)=0.671429 ;XFRAC(48)=0.685714
XFRAC(49)=0.7 ;XFRAC(50)=0.714286
XFRAC(51)=0.728571 ;XFRAC(52)=0.742857
XFRAC(53)=0.757143 ;XFRAC(54)=0.771429
XFRAC(55)=0.785714 ;XFRAC(56)=0.8
XFRAC(57)=0.814286 ;XFRAC(58)=0.828571
XFRAC(59)=0.842857 ;XFRAC(60)=0.857143
XFRAC(61)=0.871429 ;XFRAC(62)=0.885714
XFRAC(63)=0.9 ;XFRAC(64)=0.914286
XFRAC(65)=0.928571 ;XFRAC(66)=0.942857
XFRAC(67)=0.957143 ;XFRAC(68)=0.971429
XFRAC(69)=0.985714 ;XFRAC(70)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 20
YVLAST =4.
YFRAC(1)=0.05 ;YFRAC(2)=0.1
YFRAC(3)=0.15 ;YFRAC(4)=0.2
YFRAC(5)=0.25 ;YFRAC(6)=0.3
YFRAC(7)=0.35 ;YFRAC(8)=0.4
YFRAC(9)=0.45 ;YFRAC(10)=0.5
YFRAC(11)=0.55 ;YFRAC(12)=0.6
YFRAC(13)=0.65 ;YFRAC(14)=0.7
YFRAC(15)=0.75 ;YFRAC(16)=0.8
YFRAC(17)=0.85 ;YFRAC(18)=0.9
YFRAC(19)=0.95 ;YFRAC(20)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(142)=MO/T
NAME(143)=DP/T ;NAME(144)=MOTH
NAME(145)=DPTH ;NAME(146)=STMO
NAME(147)=STDP ;NAME(149)=NPOR
NAME(150)=EPOR
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(MO/T,Y,N,N,N,N,Y)
SOLUTN(DP/T,Y,N,N,N,N,Y)
SOLUTN(MOTH,Y,N,N,N,N,Y)
SOLUTN(DPTH,Y,N,N,N,N,Y)
SOLUTN(STMO,Y,N,N,N,N,Y)
SOLUTN(STDP,Y,N,N,N,N,Y)
SOLUTN(NPOR,Y,N,N,N,N,Y)
SOLUTN(EPOR,Y,N,N,N,N,Y)
EPOR = 150 ;HPOR = 0 ;NPOR = 149 ;VPOR = 0
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1. ;TMP1 =0.
EL1 =0.
TSURR =0. ;TEMP0 =0.
PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
ENUL =0.04 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRT(U1)=1. ;PRT(V1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=0. ;FIINIT(U1)=1.
FIINIT(V1)=1.0E-10 ;FIINIT(MO/T)=0.
FIINIT(DP/T)=0. ;FIINIT(MOTH)=0.
FIINIT(DPTH)=0. ;FIINIT(STMO)=0.
FIINIT(STDP)=0. ;FIINIT(NPOR)=1.
FIINIT(EPOR)=1.

PATCH(BEFORE1 ,INIVAL, 1, 20, 1, 20, 1, 1, 1, 1)
INIT(BEFORE1 ,P1 , FIXVAL ,3. )
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources

PATCH(BEFORE ,CELL , 1, 20, 1, 20, 1, 1, 1, 1)
COVAL(BEFORE ,P1 , FIXVAL ,3. )

PATCH(AFTER ,CELL , 21, 70, 1, 20, 1, 1, 1, 1)
COVAL(AFTER ,P1 , FIXVAL ,0. )

PATCH(IN ,WEST , 1, 1, 1, 20, 1, 1, 1, 1)
COVAL(IN ,P1 , FIXFLU ,1. )
COVAL(IN ,U1 ,0. ,1. )
COVAL(IN ,V1 ,0. ,0. )

PATCH(OUT ,EAST , 70, 70, 1, 20, 1, 1, 1, 1)
COVAL(OUT ,P1 ,1.0E+05 ,0. )
COVAL(OUT ,U1 ,0. ,1. )

PATCH(TOP ,NWALL , 1, 70, 20, 20, 1, 1, 1, 1)
COVAL(TOP ,U1 ,1. ,0. )
COVAL(TOP ,V1 ,1. ,0. )

PATCH(BOTTOM ,SWALL , 1, 70, 1, 1, 1, 1, 1, 1)
COVAL(BOTTOM ,U1 ,1. ,0. )
COVAL(BOTTOM ,V1 ,1. ,0. )

PATCH(PLT*1 ,EWALL , 20, 20, 1, 10, 1, 1, 1, 1)
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 2 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=-50 ;LITER(U1)=10
LITER(V1)=10
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,1.)
RELAX(V1,FALSDT,1.)
RELAX(MO/T,LINRLX,1.)
RELAX(DP/T,LINRLX,1.)
RELAX(MOTH,LINRLX,1.)
RELAX(DPTH,LINRLX,1.)
RELAX(STMO,LINRLX,1.)
RELAX(STDP,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(MO/T)=1.0E+10 ;VARMIN(MO/T)=-1.0E+10
VARMAX(DP/T)=1.0E+10 ;VARMIN(DP/T)=-1.0E+10
VARMAX(MOTH)=1.0E+10 ;VARMIN(MOTH)=-1.0E+10
VARMAX(DPTH)=1.0E+10 ;VARMIN(DPTH)=-1.0E+10
VARMAX(STMO)=1.0E+10 ;VARMIN(STMO)=-1.0E+10
VARMAX(STDP)=1.0E+10 ;VARMIN(STDP)=-1.0E+10
VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10
VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG52 = 3
ISG62 = 1
SPEDAT(SET,HPLATE,PLT*1,R,2.0E-03)
SPEDAT(SET,YOUNG,PLT*1,R,2.0E+11)
SPEDAT(SET,POISSON,PLT*1,R,0.3)
SPEDAT(SET,FIXEDEDGE,PLT*1,C,BEG)
SPEDAT(SET,POROSITY,PLT*1,R,0.)
SPEDAT(SET,STORED,DPTH!PLT*1,C,=8.531248E-04*(YG^4-4*YG^3*2.+6*YG\$)
SPEDAT(SET,STORED,DPTH!PLT*1,C,^2*4.)!ZSLFIN)
SPEDAT(SET,STORED,MOTH!PLT*1,C,=-1.5*(2.-YG)^2!ZSLFIN)
SPEDAT(SET,STORED,DP/T!PLT*1,C,=STDP/(DPTH+1.E-20)!ZSLFIN)
SPEDAT(SET,STORED,MO/T!PLT*1,C,=STMO/(MOTH+1.E-20)!ZSLFIN)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(MO/T,Y,N,Y,N,N,N)
OUTPUT(DP/T,Y,N,Y,N,N,N)
OUTPUT(MOTH,Y,N,Y,N,N,N)
OUTPUT(DPTH,Y,N,Y,N,N,N)
OUTPUT(STMO,Y,N,Y,N,N,N)
OUTPUT(STDP,Y,N,Y,N,N,N)
OUTPUT(NPOR,Y,N,Y,N,N,N)
OUTPUT(EPOR,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 25 ;IYMON = 10 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = 1 ;IXPRF = 18 ;IXPRL = 22
NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP
```