TALK=T;RUN(1,1)
PHOTON USE
* dump xxx0 means no gravity source
* dump xxxg means gravity source linked to calculated salt
* dump xxxh means gravity source linked to calculated salt
* with fixv salt source at outer boundary
* dump xxxi means gravity source linked to calculated salt
* with fixv salt source at patch whole
p
parphi
3 4 1
msg Axi-symmetrical turbulent jet with Ys multiplied by 4
msg Velocity vectors:
vec x 1 sh
upause 1
vec off;red
upause 1
msg average mixl contours:
con mixl x 1 fi;0.1
upause 1
msg concentration contours:
con conc x 1 fi;0.1
upause 1
msg average f contours:
con avef x 1 fi;0.1
upause 1
msg root-mean-square fluctuation contours:
con mnsq x 1 fi;0.1
upause 1
msg f1 contours:
con f1 x 1 fi;0.1
upause 1
msg f5 contours:
con f5 x 1 fi;0.1
upause 1
msg f10 contours:
con f10 x 1 fi;0.1
upause 1
msg eddy-viscosity (enut) contours:
con enut x 1 fi;0.1
enduse
DISPLAY
The development of the wake downstream of an underwater
vessel, the propulsion system of which creates an
axi-symmetrical jet of velocity greater than that of the
vessel.
The jet does not remain axi-symmetrical as distance from
the vessel's stern increases, because it is postulated that
the vessel is moving through density-stratified water. This
causes mixing in the vertical direction to differ
significantly from that in the horizontal direction.
MFM is used for the computation of the turbulence length
scale and the effective viscosity. Longitudinal velocity is
used as the sole PDA.
In-Form is used so as to define the density distribution in
the undisturbed water.
ENDDIS
infrbegin
debug t
stored t
source t
formula t
infrend
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(Round Jet; MFM; parabolic )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =1.
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 20
YVLAST =10.
AZYV =1. ;AZRI =0. ;AZAL =0.
Method of pairs used for grid setting.
YFRAC(2)=0.05
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = T
NZ = 100
ZWADD =200.
ZWLAST =1.
AZDZ = GRND2
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(7)=W1
NAME(16)=CONC ;NAME(134)=F1
NAME(135)=F2 ;NAME(136)=F3
NAME(137)=F4 ;NAME(138)=F5
NAME(139)=F6 ;NAME(140)=F7
NAME(141)=F8 ;NAME(142)=F9
NAME(143)=F10 ;NAME(144)=MNSQ
NAME(145)=AVEF ;NAME(146)=RATE
NAME(147)=SALT ;NAME(148)=MIXL
NAME(149)=RHO1 ;NAME(150)=ENUT
* 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,Y)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(CONC,Y,Y,N,N,N,Y)
SOLUTN(F1,Y,Y,N,N,N,Y)
SOLUTN(F2,Y,Y,N,N,N,Y)
SOLUTN(F3,Y,Y,N,N,N,Y)
SOLUTN(F4,Y,Y,N,N,N,Y)
SOLUTN(F5,Y,Y,N,N,N,Y)
SOLUTN(F6,Y,Y,N,N,N,Y)
SOLUTN(F7,Y,Y,N,N,N,Y)
SOLUTN(F8,Y,Y,N,N,N,Y)
SOLUTN(F9,Y,Y,N,N,N,Y)
SOLUTN(F10,Y,Y,N,N,N,Y)
SOLUTN(MNSQ,Y,N,N,N,N,Y)
SOLUTN(AVEF,Y,N,N,N,N,Y)
SOLUTN(RATE,Y,N,N,N,N,Y)
SOLUTN(SALT,Y,Y,N,N,N,Y)
SOLUTN(MIXL,Y,Y,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
DEN1 = 149
VIST = 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)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(CONC,N,Y,Y,Y,Y,Y)
TERMS(F1,N,Y,Y,Y,Y,Y)
TERMS(F2,N,Y,Y,Y,N,Y)
TERMS(F3,N,Y,Y,Y,Y,Y)
TERMS(F4,N,Y,Y,Y,N,Y)
TERMS(F5,N,Y,Y,Y,Y,Y)
TERMS(F6,N,Y,Y,Y,N,Y)
TERMS(F7,N,Y,Y,Y,Y,Y)
TERMS(F8,N,Y,Y,Y,N,Y)
TERMS(F9,N,Y,Y,Y,Y,Y)
TERMS(F10,N,Y,Y,Y,N,Y)
TERMS(SALT,N,Y,Y,Y,N,Y)
TERMS(MIXL,N,Y,Y,Y,Y,Y)
DIFCUT =0. ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = T
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1025. ;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 = GRND10
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1. ;PRNDTL(CONC)=1.
PRNDTL(F1)=1. ;PRNDTL(F2)=1.
PRNDTL(F3)=1. ;PRNDTL(F4)=1.
PRNDTL(F5)=1. ;PRNDTL(F6)=1.
PRNDTL(F7)=1. ;PRNDTL(F8)=1.
PRNDTL(F9)=1. ;PRNDTL(F10)=1.
PRNDTL(SALT)=1. ;PRNDTL(MIXL)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1. ;PRT(CONC)=1.
PRT(F1)=1. ;PRT(F2)=1.
PRT(F3)=1. ;PRT(F4)=1.
PRT(F5)=1. ;PRT(F6)=1.
PRT(F7)=1. ;PRT(F8)=1.
PRT(F9)=1. ;PRT(F10)=1.
PRT(SALT)=1. ;PRT(MIXL)=1.
CP1 =1. ;CP2 =1.
************************************************************
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(W1)=1.0E-10
FIINIT(CONC)=1.0E-10 ;FIINIT(F1)=1.0E-10
FIINIT(F2)=1.0E-10 ;FIINIT(F3)=1.0E-10
FIINIT(F4)=1.0E-10 ;FIINIT(F5)=1.0E-10
FIINIT(F6)=1.0E-10 ;FIINIT(F7)=1.0E-10
FIINIT(F8)=1.0E-10 ;FIINIT(F9)=1.0E-10
FIINIT(F10)=1.0E-10 ;FIINIT(MNSQ)=1.0E-10
FIINIT(AVEF)=1.0E-10 ;FIINIT(RATE)=1.0E-10
FIINIT(SALT)=1.0E-10 ;FIINIT(MIXL)=10.
FIINIT(RHO1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
PATCH(START2 ,INIVAL, 1, 1, 11, 20, 1, 1, 1, 1)
INIT(START2 ,F10 ,0. ,1. )
PATCH(START1 ,INIVAL, 1, 1, 1, 10, 1, 1, 1, 1)
INIT(START1 ,F1 ,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(FREE ,NORTH , 1, 1, 20, 20, 1, 100, 1, 1)
COVAL(FREE ,P1 , FIXFLU ,384.375031 )
COVAL(FREE ,V1 ,0. ,0. )
COVAL(FREE ,W1 ,0. ,5. )
COVAL(FREE ,CONC,0. ,0. )
COVAL(FREE ,F1 ,0. ,0. )
COVAL(FREE ,F2 ,0. ,0. )
COVAL(FREE ,F3 ,0. ,0. )
COVAL(FREE ,F4 ,0. ,0. )
COVAL(FREE ,F5 ,0. ,0. )
COVAL(FREE ,F6 ,0. ,0. )
COVAL(FREE ,F7 ,0. ,0. )
COVAL(FREE ,F8 ,0. ,0. )
COVAL(FREE ,F9 ,0. ,0. )
COVAL(FREE ,F10 ,0. ,1. )
COVAL(FREE ,SALT,In-Form:source - see Grp 19)
PATCH(UPSTJET ,LOW , 1, 1, 1, 10, 1, 1, 1, 1)
COVAL(UPSTJET ,P1 , FIXFLU ,1.025E+04 )
COVAL(UPSTJET ,W1 ,0. ,10. )
COVAL(UPSTJET ,CONC,0. ,1. )
COVAL(UPSTJET ,F1 ,0. ,1. )
COVAL(UPSTJET ,F2 ,0. ,0. )
COVAL(UPSTJET ,F3 ,0. ,0. )
COVAL(UPSTJET ,F4 ,0. ,0. )
COVAL(UPSTJET ,F5 ,0. ,0. )
COVAL(UPSTJET ,F6 ,0. ,0. )
COVAL(UPSTJET ,F7 ,0. ,0. )
COVAL(UPSTJET ,F8 ,0. ,0. )
COVAL(UPSTJET ,F9 ,0. ,0. )
COVAL(UPSTJET ,F10 ,0. ,0. )
COVAL(UPSTJET ,SALT,0. , SAME )
PATCH(UPSTFREE,LOW , 1, 1, 11, 20, 1, 1, 1, 1)
COVAL(UPSTFREE,P1 , FIXFLU ,5125. )
COVAL(UPSTFREE,W1 ,0. ,5. )
COVAL(UPSTFREE,CONC,0. ,0. )
COVAL(UPSTFREE,F1 ,0. ,0. )
COVAL(UPSTFREE,F2 ,0. ,0. )
COVAL(UPSTFREE,F3 ,0. ,0. )
COVAL(UPSTFREE,F4 ,0. ,0. )
COVAL(UPSTFREE,F5 ,0. ,0. )
COVAL(UPSTFREE,F6 ,0. ,0. )
COVAL(UPSTFREE,F7 ,0. ,0. )
COVAL(UPSTFREE,F8 ,0. ,0. )
COVAL(UPSTFREE,F9 ,0. ,0. )
COVAL(UPSTFREE,F10 ,0. ,0. )
COVAL(UPSTFREE,SALT,0. , SAME )
PATCH(MICROMIX,PHASEM, 1, 1, 1, 20, 1, 100, 1, 1)
PATCH($MNSQ ,PHASEM, 1, 1, 1, 20, 1, 100, 1, 1)
COVAL($MNSQ ,MIXL,5.0E-06 ,1.0E+05 )
PATCH(WHOLE ,PHASEM, 1, 1, 1, 20, 1, 100, 1, 1)
COVAL(WHOLE ,SALT,In-Form:source - see Grp 19)
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
IPARAB = 0
AZPH =0. ;PBAR =0.
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1 ;ISWC1 = 1
LITHYD = 100 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(P1)=100 ;LITER(U1)=10
LITER(V1)=10 ;LITER(W1)=10
LITER(CONC)=20 ;LITER(F1)=20
LITER(F2)=20 ;LITER(F3)=20
LITER(F4)=20 ;LITER(F5)=20
LITER(F6)=20 ;LITER(F7)=20
LITER(F8)=20 ;LITER(F9)=20
LITER(F10)=20 ;LITER(SALT)=20
LITER(MIXL)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(CONC)=1.0E-03 ;ENDIT(F1)=1.0E-03
ENDIT(F2)=1.0E-03 ;ENDIT(F3)=1.0E-03
ENDIT(F4)=1.0E-03 ;ENDIT(F5)=1.0E-03
ENDIT(F6)=1.0E-03 ;ENDIT(F7)=1.0E-03
ENDIT(F8)=1.0E-03 ;ENDIT(F9)=1.0E-03
ENDIT(F10)=1.0E-03 ;ENDIT(SALT)=1.0E-03
ENDIT(MIXL)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.5)
RELAX(U1,FALSDT,1.0E-02)
RELAX(V1,FALSDT,1.)
RELAX(W1,FALSDT,10.)
RELAX(CONC,FALSDT,1.0E+09)
RELAX(F1,LINRLX,0.5)
RELAX(F2,LINRLX,0.5)
RELAX(F3,LINRLX,0.5)
RELAX(F4,LINRLX,0.5)
RELAX(F5,LINRLX,0.5)
RELAX(F6,LINRLX,0.5)
RELAX(F7,LINRLX,0.5)
RELAX(F8,LINRLX,0.5)
RELAX(F9,LINRLX,0.5)
RELAX(F10,LINRLX,0.5)
RELAX(MNSQ,LINRLX,1.)
RELAX(AVEF,LINRLX,1.)
RELAX(RATE,LINRLX,1.)
RELAX(SALT,LINRLX,0.1)
RELAX(MIXL,FALSDT,1.0E+09)
RELAX(RHO1,LINRLX,1.)
RELAX(ENUT,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)=1000. ;VARMIN(V1)=-1000.
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(CONC)=1.0E+10 ;VARMIN(CONC)=-1.0E+10
VARMAX(F1)=1. ;VARMIN(F1)=0.
VARMAX(F2)=1. ;VARMIN(F2)=0.
VARMAX(F3)=1. ;VARMIN(F3)=0.
VARMAX(F4)=1. ;VARMIN(F4)=0.
VARMAX(F5)=1. ;VARMIN(F5)=0.
VARMAX(F6)=1. ;VARMIN(F6)=0.
VARMAX(F7)=1. ;VARMIN(F7)=0.
VARMAX(F8)=1. ;VARMIN(F8)=0.
VARMAX(F9)=1. ;VARMIN(F9)=0.
VARMAX(F10)=1. ;VARMIN(F10)=0.
VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0.
VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10
VARMAX(SALT)=1.0E+10 ;VARMIN(SALT)=-1.0E+10
VARMAX(MIXL)=1.0E+10 ;VARMIN(MIXL)=-1.0E+10
VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
DZW1 =0.1
ISG62 = 1
SPEDAT(SET,PROPERTY,RHO1,C,=SALT)
SPEDAT(SET,INITIAL,SALT,C,=1027.+0.5*(TANH(YG*COS(XG)/YVLAST)+1)*$)
SPEDAT(SET,INITIAL,SALT,C,(1023.-1027.))
SPEDAT(SET,SOURCE,SALT!FREE,C,=1027.+0.5*(TANH(YG*COS(XG)/YVLAST)$)
SPEDAT(SET,SOURCE,SALT!FREE,C,+1)*(1023.-1027.)!FIXVAL)
SPEDAT(SET,SOURCE,SALT!WHOLE,C,=1027.+0.5*(TANH(YG*COS(XG)/YVLAST$)
SPEDAT(SET,SOURCE,SALT!WHOLE,C,)+1)*(1023.-1027.)!FIXVAL)
SPEDAT(SET,MFM,MFMMOD,C,MNSQ)
SPEDAT(SET,MFM,NFLUIDS,I,10)
SPEDAT(SET,MFM,NFLR,I,1)
SPEDAT(SET,MFM,NFLF,I,10)
SPEDAT(SET,MFM,CONMIX,R,5.)
SPEDAT(SET,MFM,VISCON,R,10.)
SPEDAT(SET,MFM,POPMIN,R,5.)
SPEDAT(SET,MFM,POPMAX,R,10.)
SPEDAT(SET,GXMONI,PLOTALL,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(W1,Y,N,Y,Y,Y,Y)
OUTPUT(CONC,Y,N,Y,Y,Y,Y)
OUTPUT(F1,Y,N,Y,Y,Y,Y)
OUTPUT(F2,Y,N,Y,Y,Y,Y)
OUTPUT(F3,Y,N,Y,Y,Y,Y)
OUTPUT(F4,Y,N,Y,Y,Y,Y)
OUTPUT(F5,Y,N,Y,Y,Y,Y)
OUTPUT(F6,Y,N,Y,Y,Y,Y)
OUTPUT(F7,Y,N,Y,Y,Y,Y)
OUTPUT(F8,Y,N,Y,Y,Y,Y)
OUTPUT(F9,Y,N,Y,Y,Y,Y)
OUTPUT(F10,Y,N,Y,Y,Y,Y)
OUTPUT(MNSQ,Y,N,Y,N,N,N)
OUTPUT(AVEF,Y,N,Y,N,N,N)
OUTPUT(RATE,Y,N,Y,N,N,N)
OUTPUT(SALT,Y,N,Y,Y,Y,Y)
OUTPUT(MIXL,Y,N,Y,Y,Y,Y)
OUTPUT(RHO1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;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
NYPRIN = 4 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = 100 ;IZPRF = 1 ;IZPRL = 10000
IPLTF = 1 ;IPLTL = 30 ;NPLT = 1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 1 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.2
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
PATCH(IZEQNZ ,PROFIL, 1, 1, 1, 20, 100, 100, 1, 1)
PLOT(IZEQNZ ,W1 ,0. ,0. )
PLOT(IZEQNZ ,CONC,0. ,0. )
PLOT(IZEQNZ ,ENUT,0. ,0. )
PATCH(MIDDLE ,PROFIL, 1, 1, 10, 10, 1, 100, 1, 1)
PLOT(MIDDLE ,F1 ,0. ,0. )
PLOT(MIDDLE ,F4 ,0. ,0. )
PLOT(MIDDLE ,F7 ,0. ,0. )
PLOT(MIDDLE ,F10 ,0. ,0. )
PATCH(PROF1 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF1 ,F1 ,0. ,0. )
PATCH(PROF2 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF2 ,F2 ,0. ,0. )
PATCH(PROF3 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF3 ,F3 ,0. ,0. )
PATCH(PROF4 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF4 ,F4 ,0. ,0. )
PATCH(PROF5 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF5 ,F5 ,0. ,0. )
PATCH(PROF6 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF6 ,F6 ,0. ,0. )
PATCH(PROF7 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF7 ,F7 ,0. ,0. )
PATCH(PROF8 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF8 ,F8 ,0. ,0. )
PATCH(PROF9 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF9 ,F9 ,0. ,0. )
PATCH(PROF10 ,PROFIL, 1, 1, 1, 20, 1, 100, 1, 1)
PLOT(PROF10 ,F10 ,0. ,0. )
PATCH(FINAL ,PROFIL, 1, 1, 1, 20, 100, 100, 1, 1)
PLOT(FINAL ,W1 ,0. ,0. )
PLOT(FINAL ,MNSQ,0. ,0. )
PLOT(FINAL ,MIXL,0. ,0. )
PLOT(FINAL ,ENUT,0. ,0. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
STOP