TALK=T;RUN(1,1)
DISPLAY
A fluid (liquid steel) flows steadily in a direction normal to$
a
thick plate in channel. The Reynolds Number (i.e. velocity times
channel width divided by kinematic viscosity) equals 10.
For calculation deformation of plate is used CDM or SDM models.
ENDDIS
PHOTON USE
p;phi
0.4 1 1
set prop off
msg
msg PRPS contours
msg
cont PRPS z 1 fil;.0001
gr OU Z 1
pause
msg
msg P1 contours
msg
GR OU Z 1;
cont P1 z 1 fil;.0001
pause
cl
msg
msg
msg Velocity vector
GR OU Z 1;
vec z 1 SH DASH 0
pause
cl
msg
msg U1 contours
msg
GR OU Z 1;
cont U1 z 1 fil;.0001
pause
cl
msg
msg V1 contours
msg
GR OU Z 1;
cont V1 z 1 fil;.0001
pause
cl
msg
msg EPSX contours
msg
GR OU Z 1;
cont EPSX z 1 x 11 20 y 1 20 fil;.0001
pause
cl
msg
msg EPSY contours
msg
GR OU Z 1;
cont EPSY z 1 x 11 20 y 1 20 fil;.0001
pause
cl
msg
msg STRX contours
msg
GR OU Z 1;
cont STRX z 1 x 11 20 y 1 20 fil;.0001
pause
cl
msg
msg STRY contours
msg
GR OU Z 1;
cont STRY z 1 x 11 20 y 1 20 fil;.0001
pause
cl
set prop off
msg
msg X-displacement contours
msg
GR OU Z 1;
cont D_X z 1 x 11 20 y 1 20 fil;.0001
pause
cl
set prop off
msg
msg Y-displacement contours
msg
GR OU Z 1;
cont D_Y z 1 x 11 20 y 1 20 fil;.0001
pause
cl
set prop off
msg
msg Displacements vector
msg
GR OU Z 1;
SET VEC COMP - - -
SET VEC COMP D_X D_Y -
VECT z 1 x 11 20 y 1 20 SH DASH 0
pause
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(2D Fluid flow around an beam; S801 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 14
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 50
XULAST =4.2
XFRAC(1)=0.047619 ;XFRAC(2)=0.095238
XFRAC(3)=0.142857 ;XFRAC(4)=0.190476
XFRAC(5)=0.238095 ;XFRAC(6)=0.285714
XFRAC(7)=0.333333 ;XFRAC(8)=0.380952
XFRAC(9)=0.428571 ;XFRAC(10)=0.476191
XFRAC(11)=0.480952 ;XFRAC(12)=0.485714
XFRAC(13)=0.490476 ;XFRAC(14)=0.495238
XFRAC(15)=0.5 ;XFRAC(16)=0.504762
XFRAC(17)=0.509524 ;XFRAC(18)=0.514286
XFRAC(19)=0.519048 ;XFRAC(20)=0.52381
XFRAC(21)=0.539683 ;XFRAC(22)=0.555556
XFRAC(23)=0.571429 ;XFRAC(24)=0.587302
XFRAC(25)=0.603175 ;XFRAC(26)=0.619048
XFRAC(27)=0.634921 ;XFRAC(28)=0.650794
XFRAC(29)=0.666667 ;XFRAC(30)=0.68254
XFRAC(31)=0.698413 ;XFRAC(32)=0.714286
XFRAC(33)=0.730159 ;XFRAC(34)=0.746032
XFRAC(35)=0.761905 ;XFRAC(36)=0.777778
XFRAC(37)=0.793651 ;XFRAC(38)=0.809524
XFRAC(39)=0.825397 ;XFRAC(40)=0.84127
XFRAC(41)=0.857143 ;XFRAC(42)=0.873016
XFRAC(43)=0.888889 ;XFRAC(44)=0.904762
XFRAC(45)=0.920635 ;XFRAC(46)=0.936508
XFRAC(47)=0.952381 ;XFRAC(48)=0.968254
XFRAC(49)=0.984127 ;XFRAC(50)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 40
YVLAST =1.
YFRAC(1)=0.025 ;YFRAC(2)=0.05
YFRAC(3)=0.075 ;YFRAC(4)=0.1
YFRAC(5)=0.125 ;YFRAC(6)=0.15
YFRAC(7)=0.175 ;YFRAC(8)=0.2
YFRAC(9)=0.225 ;YFRAC(10)=0.25
YFRAC(11)=0.275 ;YFRAC(12)=0.3
YFRAC(13)=0.325 ;YFRAC(14)=0.35
YFRAC(15)=0.375 ;YFRAC(16)=0.4
YFRAC(17)=0.425 ;YFRAC(18)=0.45
YFRAC(19)=0.475 ;YFRAC(20)=0.5
YFRAC(21)=0.525 ;YFRAC(22)=0.55
YFRAC(23)=0.575 ;YFRAC(24)=0.6
YFRAC(25)=0.625 ;YFRAC(26)=0.65
YFRAC(27)=0.675 ;YFRAC(28)=0.7
YFRAC(29)=0.725 ;YFRAC(30)=0.75
YFRAC(31)=0.775 ;YFRAC(32)=0.8
YFRAC(33)=0.825 ;YFRAC(34)=0.85
YFRAC(35)=0.875 ;YFRAC(36)=0.9
YFRAC(37)=0.925 ;YFRAC(38)=0.95
YFRAC(39)=0.975 ;YFRAC(40)=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(138)=D_Y
NAME(139)=D_X ;NAME(140)=DRH1
NAME(141)=DVO1 ;NAME(142)=ENUL
NAME(143)=DEN1 ;NAME(144)=EPSZ
NAME(145)=EPSX ;NAME(146)=EPSY
NAME(147)=STXY ;NAME(148)=STRY
NAME(149)=STRX ;NAME(150)=PRPS
* 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(D_Y,Y,N,N,N,N,Y)
SOLUTN(D_X,Y,N,N,N,N,Y)
SOLUTN(DRH1,Y,N,N,N,N,Y)
SOLUTN(DVO1,Y,N,N,N,N,Y)
SOLUTN(ENUL,Y,N,N,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(EPSZ,Y,N,N,N,N,Y)
SOLUTN(EPSX,Y,N,N,N,N,Y)
SOLUTN(EPSY,Y,N,N,N,N,Y)
SOLUTN(STXY,Y,N,N,N,N,Y)
SOLUTN(STRY,Y,N,N,N,N,Y)
SOLUTN(STRX,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
DEN1 = 143
VISL = 142
PRPS = 150
************************************************************
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.
RADIA =0. ;RADIB =0.
ENUL =1.0E-05 ;ENUT =0.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRT(U1)=1. ;PRT(V1)=1.
CP1 =1. ;CP2 =1.
* List of user-defined materials to be read by EARTH
MATFLG=T;IMAT=2
* Name
*Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr.
*
99 7800.0 1.0 473.0 43.0 0.0 0.0
*
100 2700.0 0.35 896.0 204.0 2.35E-5 1.47E-11
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
FIINIT(V1)=1.0E-10 ;FIINIT(D_Y)=0.
FIINIT(D_X)=0. ;FIINIT(DRH1)=1.0E-10
FIINIT(DVO1)=1.0E-10 ;FIINIT(ENUL)=1.0E-10
FIINIT(DEN1)=1.0E-10 ;FIINIT(EPSZ)=1.0E-10
FIINIT(EPSX)=1.0E-10 ;FIINIT(EPSY)=1.0E-10
FIINIT(STXY)=1.0E-10 ;FIINIT(STRY)=1.0E-10
FIINIT(STRX)=1.0E-10 ;FIINIT(PRPS)=99.
PATCH(BODY ,INIVAL, 11, 20, 1, 20, 1, 1, 1, 1)
INIT(BODY ,PRPS,1. ,100. )
INIADD = F
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(IN1 ,WEST , 1, 1, 1, 40, 1, 1, 1, 1)
COVAL(IN1 ,P1 , FIXFLU ,7.8E+04 )
COVAL(IN1 ,U1 ,0. ,10. )
COVAL(IN1 ,V1 ,0. ,0. )
PATCH(OUT ,EAST , 50, 50, 1, 40, 1, 1, 1, 1)
COVAL(OUT ,P1 ,1.0E+05 ,0. )
COVAL(OUT ,U1 ,0. ,10. )
PATCH(TOP ,NWALL , 1, 50, 40, 40, 1, 1, 1, 1)
COVAL(TOP ,U1 ,1. ,0. )
COVAL(TOP ,V1 ,1. ,0. )
PATCH(BOTTOM1 ,SWALL , 1, 10, 1, 1, 1, 1, 1, 1)
COVAL(BOTTOM1 ,U1 ,1. ,0. )
COVAL(BOTTOM1 ,V1 ,1. ,0. )
PATCH(BOTTOM2 ,SWALL , 21, 50, 1, 1, 1, 1, 1, 1)
COVAL(BOTTOM2 ,U1 ,1. ,0. )
COVAL(BOTTOM2 ,V1 ,1. ,0. )
PATCH(BODYFIX ,CELL , 11, 20, 1, 1, 1, 1, 1, 1)
COVAL(BODYFIX ,U1 ,1.0E+05 ,0. )
COVAL(BODYFIX ,V1 ,1.0E+05 ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 4000 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-07
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(U1)=100
LITER(V1)=100
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,100.)
RELAX(V1,FALSDT,100.)
RELAX(D_Y,LINRLX,1.)
RELAX(D_X,LINRLX,1.)
RELAX(DRH1,LINRLX,1.)
RELAX(DVO1,LINRLX,1.)
RELAX(ENUL,LINRLX,1.)
RELAX(DEN1,LINRLX,1.)
RELAX(EPSZ,LINRLX,1.)
RELAX(EPSX,LINRLX,1.)
RELAX(EPSY,LINRLX,1.)
RELAX(STXY,LINRLX,1.)
RELAX(STRY,LINRLX,1.)
RELAX(STRX,LINRLX,1.)
RELAX(PRPS,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(D_Y)=1.0E+10 ;VARMIN(D_Y)=-1.0E+10
VARMAX(D_X)=1.0E+10 ;VARMIN(D_X)=-1.0E+10
VARMAX(DRH1)=1.0E+10 ;VARMIN(DRH1)=-1.0E+10
VARMAX(DVO1)=1.0E+10 ;VARMIN(DVO1)=-1.0E+10
VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
VARMAX(EPSZ)=1.0E+10 ;VARMIN(EPSZ)=-1.0E+10
VARMAX(EPSX)=1.0E+10 ;VARMIN(EPSX)=-1.0E+10
VARMAX(EPSY)=1.0E+10 ;VARMIN(EPSY)=-1.0E+10
VARMAX(STXY)=1.0E+10 ;VARMIN(STXY)=-1.0E+10
VARMAX(STRY)=1.0E+10 ;VARMIN(STRY)=-1.0E+10
VARMAX(STRX)=1.0E+10 ;VARMIN(STRX)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
STRA = T
PARSOL = F
ISG52 = 3
ISG62 = 1
CSG10 ='Q1'
SPEDAT(SET,RLXFAC,RLXU1D,R,0.5)
SPEDAT(SET,RLXFAC,RLXV1D,R,0.5)
SPEDAT(SET,STORED,D_X,C,=U1!IMAT>99!ZSLFIN)
SPEDAT(SET,STORED,D_Y,C,=V1!IMAT>99!ZSLFIN)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,MATERIAL,99,L,T)
SPEDAT(SET,MATERIAL,100,L,T)
************************************************************
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(D_Y,Y,N,Y,N,N,N)
OUTPUT(D_X,Y,N,Y,N,N,N)
OUTPUT(DRH1,Y,N,Y,N,N,N)
OUTPUT(DVO1,Y,N,Y,N,N,N)
OUTPUT(ENUL,Y,N,Y,N,N,N)
OUTPUT(DEN1,Y,N,Y,N,N,N)
OUTPUT(EPSZ,Y,N,Y,N,N,N)
OUTPUT(EPSX,Y,N,Y,N,N,N)
OUTPUT(EPSY,Y,N,Y,N,N,N)
OUTPUT(STXY,Y,N,Y,N,N,N)
OUTPUT(STRY,Y,N,Y,N,N,N)
OUTPUT(STRX,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 12 ;IYMON = 18 ;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 = 8 ;IXPRL = 22
NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 22
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