TALK=T;RUN(1,1)
PHOTON USE
p;;;;;
gr ou z 1 1
gr ou y m z 1 1;gr ou z 2;gr ou x 1 z 1 1 ;gr ou x m z 1 1
gr ou y m z 2 2;gr ou z 3;gr ou x 1 z 2 2 ;gr ou x m z 2 2
gr ou y m z 3 3;gr ou z 4;gr ou x 1 z 3 3 ;gr ou x m z 3 3
gr ou y m z 4 4;gr ou z 5;gr ou x 1 z 4 4 ;gr ou x m z 4 4
vec z 1 sh
msg
msg In slab of IZ=1 the flow is 45 degree
msg from right bottom corner.
msg
pause
vec off;con off;red
vec z 2 sh
msg
msg In slab IZ=2 the flow is a solid body
msg counter clockwise rotation.
msg
pause
vec off;con off;red
con h1 z 3 fil;.001
vec z 3 sh
msg
msg In slab IZ=3 the flow is the superposition
msg of those in slabs 1 and 2.
msg
msg The source dispersion plume is as shown.
msg
msg Hit Enter to get the result of convection
msg fluxes transformation.
pause
vec off;con off;red
con h1 z 4 fil;.001
vec z 4 sh
msg
MSG
msg In slab IZ=4 the dispersion is exactly the same as
msg before in spite of background flow being in opposite
msg direction.
ENDUSE
DISPLAY
The four problems are arranged in four slab cartesian
box. Each slab is devoted to the particular problem.
The PLANTing of initial velocity fields, manipulation
with them and the alteration of convection fluxes by
add-extra-velocity option are demonstrated.
ENDDIS
PLANTBEGIN
PATCH(INITZ2,INIVAL,1,NX,1,NY,2,2,1,1)
VAL= -(YG2D-15.)
COVAL(INITZ2,U1,zero,GRND)
VAL= XG2D-15.
COVAL(INITZ2,V1,zero,GRND)
The velocity field in the second slab is initialised as
solid body rotation.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PATCH(INITZ3,INIVAL,1,NX,1,NY,3,3,1,1)
VAL= U1[,,1]+U1[,,2]+YG2D
COVAL(INITZ3,U1,zero,GRND)
VAL= V1[,,1]+V1[,,2]-XG2D
COVAL(INITZ3,V1,zero,GRND)
The velocity field for the third sub-domain, IZ=3, is
the superposition of the velocity components in two
previous slabs and the cell node coordinates. It
results in 45 degree flow of 10 m/s from north-west
edge of the domain.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
NAMSAT=MOSG
U1AD=GRND
VELAD=U1[,,2]+YG2D
REGION(,NX-1) /IZ.EQ.4
V1AD=GRND
VELAD=V1[,,2]-XG2D
REGION(,,,NY-1)
IF(IZ.EQ.4)
Although, the velocity field at the fourth slab is the
same as for first one, the add-extra-velocity option
is activated as above pointers tell. The extra
velocities added to the main velocity components
alters the convection fluxes to be equal to ones at
third slab. The expected distribution of convected
property H1 should be identical to the one at IZ=3.
Please note the differences in REGION commands. They
are attributed to the staggered nature of velocity
nodes and the usage of either switch or IF command to
limit the Z-direction extent of velocity alterations.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT( Alteration of convection fluxes )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 621
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 20
XULAST =20.
XFRAC(1)=0.05 ;XFRAC(2)=0.1
XFRAC(3)=0.15 ;XFRAC(4)=0.2
XFRAC(5)=0.25 ;XFRAC(6)=0.3
XFRAC(7)=0.35 ;XFRAC(8)=0.4
XFRAC(9)=0.45 ;XFRAC(10)=0.5
XFRAC(11)=0.55 ;XFRAC(12)=0.6
XFRAC(13)=0.65 ;XFRAC(14)=0.7
XFRAC(15)=0.75 ;XFRAC(16)=0.8
XFRAC(17)=0.85 ;XFRAC(18)=0.9
XFRAC(19)=0.95 ;XFRAC(20)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 20
YVLAST =20.
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 = 4
ZWLAST =4.
ZFRAC(1)=0.25 ;ZFRAC(2)=0.5
ZFRAC(3)=0.75 ;ZFRAC(4)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(3)=U1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(14)=H1
NAME(150)=HPOR
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(U1,Y,N,N,N,N,Y)
SOLUTN(V1,Y,N,N,N,N,Y)
SOLUTN(W1,Y,N,N,N,N,Y)
SOLUTN(H1,Y,Y,N,N,N,Y)
SOLUTN(HPOR,Y,N,N,N,N,Y)
EPOR = 0 ;HPOR = 150 ;NPOR = 0 ;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(H1,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
U1AD = GRND ;V1AD = GRND ;W1AD =0.
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-04 ;ENUT =0.
PRNDTL(H1)=1.
PRT(H1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(U1)=-5. ;FIINIT(V1)=5.
FIINIT(W1)=1.0E-10 ;FIINIT(H1)=1.0E-10
FIINIT(HPOR)=0.
PATCH(INITZ2 ,INIVAL, 1, 20, 1, 20, 2, 2, 1, 1)
INIT(INITZ2 ,U1 ,0. , GRND )
INIT(INITZ2 ,V1 ,0. , GRND )
PATCH(INITZ3 ,INIVAL, 1, 20, 1, 20, 3, 3, 1, 1)
INIT(INITZ3 ,U1 ,0. , GRND )
INIT(INITZ3 ,V1 ,0. , GRND )
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(FIXSOR ,CELL , 10, 10, 10, 10, 3, 4, 1, 1)
COVAL(FIXSOR ,H1 , FIXVAL ,1. )
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(H1)=20
ENDIT(H1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(U1,LINRLX,1.)
RELAX(V1,LINRLX,1.)
RELAX(W1,LINRLX,1.)
RELAX(H1,FALSDT,1.0E+09)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(U1)=6.563 ;EX(V1)=6.563
EX(W1)=1.0E-10 ;EX(H1)=0.02162
EX(HPOR)=1.0E-11
************************************************************
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(U1,Y,N,Y,N,N,N)
OUTPUT(V1,Y,N,Y,N,N,N)
OUTPUT(W1,Y,N,Y,N,N,N)
OUTPUT(H1,Y,N,Y,Y,Y,Y)
OUTPUT(HPOR,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 5 ;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 = 1 ;IXPRL = 10000
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
STOP