TALK=T;RUN(1,1)
PHOTON USE
p
n4
up z
msg contours of S at time = 0.1s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n8
up z
msg contours of S at time = 0.20s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n12
up z
msg contours of S at time = 0.35s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n16
up z
msg contours of S at time = 0.40s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n20
up z
msg contours of S at time = 0.50s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n24
up z
msg contours of S at time = 0.60s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n28
up z
msg contours of S at time = 0.70s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n32
up z
msg contours of S at time = 0.80s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n36
up z
msg contours of S at time = 0.90s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
MSG Press Enter to continue
pause
p
n40
up z
msg contours of S at time = 1.0s
con s1 x 1 sh;.49 .52 40
msg velocity vectors
vec x 1
msg -
msg Press e to END
enduse
DISPLAY
This case simulates the slumping of a liquid column that may
arise, for example, when container walls break down. In
addition an obstacle can be errected in the path of the liquid.
liquid
| column
|------ z ^
| - - | |
| - | +----->
| -- | | wall y
| -- | |
|_____|______|___
/////////////////
The run shows that use of PLANT results in compact Q1 coding to
introduce the basic features of SEM method. It should be also
noted that upwind convection scheme may not perform well enough
suffering exessive false diffusion. The low dispersion schemes
give significantly higher accuracy in free surface location
calculation. For more details and comparison with experiment see
the entry SEM of PHOENICS 2.0 and older and references cited
therein. One of them, namely, of C.W.Hirt and B.D. Nichols in J.
Comput. Phys., v. 39, 1981, pp.201- , is openly available.
ENDDIS
PLANTBEGIN
** Efective pressure gradient sources
PATCH(DPDZ,CELL,1,NX,1,NY,1,NZ,1,lstep)
VAL=(1.-1./(1000.*S1+(1.-S1)))*AHIGH*(HIGH(P1)-P1)
COVAL(DPDZ,W1,FIXFLU,GRND)
PATCH(DPDY,CELL,1,NX,1,NY,1,NZ,1,lstep)
VAL=(1.-1./(1000.*S1+(1.-S1)))*ANORTH*(NORTH(P1)-P1)
COVAL(DPDY,V1,FIXFLU,GRND)
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT( Slumping of a liquid column by SEM meth)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 146
************************************************************
Group 2. Time dependence
STEADY = F
* Set overall time and no. of steps
TFIRST =0. ;TLAST =1.
FSTEP = 1 ;LSTEP = 40
Method of pairs used for grid setting.
TFRAC(2)=0.025
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 1
XULAST =1.
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 10
YVLAST =2.
YFRAC(1)=0.1 ;YFRAC(2)=0.2
YFRAC(3)=0.3 ;YFRAC(4)=0.4
YFRAC(5)=0.5 ;YFRAC(6)=0.6
YFRAC(7)=0.7 ;YFRAC(8)=0.8
YFRAC(9)=0.9 ;YFRAC(10)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 10
ZWLAST =0.5
ZFRAC(1)=0.1 ;ZFRAC(2)=0.2
ZFRAC(3)=0.3 ;ZFRAC(4)=0.4
ZFRAC(5)=0.5 ;ZFRAC(6)=0.6
ZFRAC(7)=0.7 ;ZFRAC(8)=0.8
ZFRAC(9)=0.9 ;ZFRAC(10)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(150)=S1
* 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,N,N,N,N)
SOLUTN(V1,Y,Y,N,N,N,N)
SOLUTN(W1,Y,Y,N,N,N,N)
SOLUTN(S1,Y,Y,N,N,N,Y)
************************************************************
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(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(S1,N,Y,N,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = T ;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(V1)=1. ;PRNDTL(W1)=1.
PRT(V1)=1. ;PRT(W1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=0. ;FIINIT(V1)=0.
FIINIT(W1)=0. ;FIINIT(S1)=0.
PATCH(LIQUID ,INIVAL, 1, 1, 1, 5, 1, 5, 1, 1)
INIT(LIQUID ,S1 ,0. ,1. )
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(REFP ,CELL , 1, 1, 1, 10, 10, 10, 1, 40)
COVAL(REFP ,P1 ,1. ,0. )
PATCH(GRAV ,PHASEM, 1, 1, 1, 10, 1, 10, 1, 40)
COVAL(GRAV ,W1 , FIXFLU ,-9.81 )
PATCH(DPDZ ,CELL , 1, 1, 1, 10, 1, 10, 1, 40)
COVAL(DPDZ ,W1 , FIXFLU , GRND )
PATCH(DPDY ,CELL , 1, 1, 1, 10, 1, 10, 1, 40)
COVAL(DPDY ,V1 , FIXFLU , GRND )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 10 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(V1)=10
LITER(W1)=10 ;LITER(S1)=1
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(S1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(V1,FALSDT,0.1)
RELAX(W1,FALSDT,0.1)
RELAX(S1,FALSDT,1.0E+09)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(S1)=1. ;VARMIN(S1)=0.
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=70.82 ;EX(V1)=0.4885
EX(W1)=0.1421 ;EX(S1)=0.134
************************************************************
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,Y,Y,Y,Y,Y)
OUTPUT(V1,Y,Y,Y,Y,Y,Y)
OUTPUT(W1,Y,Y,Y,Y,Y,Y)
OUTPUT(S1,Y,Y,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 5 ;IZMON = 5
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NTPRIN = 20 ;ISTPRF = 1 ;ISTPRL = 100000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
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
IDISPA = 4 ;IDISPB = 0 ;IDISPC = 0
CSG1 ='N'
STOP