PHOTON USE
p
parphi
1 4 1
msg Axi-symmetrical turbulent jet with Ys multiplied by 4
msg Velocity vectors:
vec x 1 sh
pause; vec off;red
msg average mixl contours:
con mixl x 1 fi;0.001;pause; con off;red
msg concentration contours:
con conc x 1 fi;0.001; pause; con off;red
msg average f contours:
con avef x 1 fi;0.001; pause; con off;red
msg root-mean-square fluctuation contours:
con mnsq x 1 fi;0.001; pause; con off;red
msg f1 contours:
con f1 x 1 fi;0.001; pause; con off;red
msg f5 contours:
con f5 x 1 fi;0.001; pause; con off;red
msg f10 contours:
con f10 x 1 fi;0.001; pause; con off;red
msg eddy-viscosity (enut) contours:
con enut x 1 fi;0.001 ; pause; con off;red
msg dimensionless length contourscontours:
con diml x 1 fi;0.001 ; pause; con off;red
msg dimensionless viscosity contours:
con dimv x 1 fi;0.001
enduse
l(cls
GROUP 1. Run title and other preliminaries
TEXT(Round Jet; MFM; parabolic
TITLE
The locally-defined variables are as follows:
WJET Jet velocity at the inlet (m/s)
CJET The jet concentration at the inlet (C)
CFREE The concentration of the free stream (C)
REAL(WJET,WFREE,CJET,CFREE)
WJET=10.0; WFREE=0.5*WJET; CJET=1.0; CFREE=0.0
INTEGER(NFLUIDS,NFLR,NFLF)
CHAR(MFMMOD)
NFLF=10; NFLR=1
NFLUIDS=NFLF*NFLR
REAL(RATEXP,VISCON,LENCON,CONMIX)
MFMMOD=MNSQ; LENCON=0.5; VISCON=20.0; CONMIX=5.0 ! model & constants
! see also enut = grnd10
! $MNSQ patch
GROUP 3. X-direction grid specification
CARTES=F;XULAST=0.1
nx=1
grdpwr(x,12,6.28,1)
GROUP 4. Y-direction grid specification
*** Linear grid expansion of DYGDZ
change dimensions and grid
REAL(REALNY,REALNZ)
NY=20; NZ=100
ny=5;nyprin=1
REALNY=NY; REALNZ=NZ
YVLAST = 1.000E-01 ! half-width of the injection nozzle
YFRAC(1) = -REALNY ;YFRAC(2) = 1./REALNY ! uniform y-grid
AZYV=1.0 ! exponent in width-expansion formula.
! 1.0 signifies linear expansion
ZWADD=20.0*YVLAST ! indirectly sets grid-enlargement angle
GROUP 5. Z-direction grid specification
PARAB=T
NZ=100
AZDZ=PROPY ! i.e. grnd2, means dz is proportional to yvlast
DZW1=0.1 ! the proportionality constant
GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1,C1); NAME(C1)=CONC; STORE(ENUT,LEN1)
SOLVE(MIXL)
STORE(EPKE,RATE)
STORE(AVEF); VARMAX(AVEF)=1.0; VARMIN(AVEF)=0.0
STORE(MNSQ); VARMAX(MNSQ)=1.0; VARMIN(MNSQ)=0.0
INTEGER(III)
III=NFLUIDS+1
DO II = 1,NFLUIDS ! DO loop order chosen so that
III = III-1 ! F1 is solved first
SOLVE(F:iii:)
VARMAX(F:iii:)=1.0; VARMIN(F:iii:)=0.0; RELAX(F:iii:,LINRLX,0.5)
ENDDO
inform7begin
(stored var diml is mixl/yvlast )
(stored var dimm is mnsq/w1 )
(stored var dimv is enut/(yvlast*(w1[1&1&]-:wfree: ) ) )
inform7end
GROUP 8. Terms (in differential equations) & devices
DIFCUT=0.0; TERMS(CONC,N,Y,Y,Y,Y,Y)
GROUP 9. Properties of the medium (or media)
ENUT=GRND10 ! this makes turbulent viscosity viscon*mnsq*mixl
(property var len1 is mixl)
GROUP 11 initial values
fiinit(len1)= yvlast ! it is important to initialise to
fiinit(mixl)= yvlast ! reasonable values
PATCH(START2,INIVAL,1,nx,NY/2+1,NY,1,1,1,1)
COVAL(START2,F:NFLUIDS:,0.0,1.0)
PATCH(START1,INIVAL,1,1,1,NY/2,1,1,1,1)
COVAL(START1,F1,0.0,1.0)
GROUP 13. Boundary conditions and special sources
1. Outer Boundary-- free stream
PATCH(FREE,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(FREE,P1,1.0E4,0.0);COVAL(FREE,W1,ONLYMS,WFREE)
COVAL(FREE,V1,ONLYMS,0.0);COVAL(FREE,CONC,ONLYMS,CFREE)
DO II = 1,NFLUIDS
COVAL(FREE,F:II:,ONLYMS,0.0)
ENDDO
COVAL(FREE,F:NFLUIDS:,ONLYMS,1.0)
2. Inlet Boundary-- nozzle
PATCH(NOZZLE,LOW,1,nx,1,NY/2,1,1,1,1)
COVAL(NOZZLE,P1,FIXFLU,RHO1*WJET);COVAL(NOZZLE,W1,ONLYMS,WJET)
COVAL(NOZZLE,CONC,ONLYMS,CJET)
COVAL(NOZZLE,F1,FIXVAL,1.0)
DO II = 1,NFLUIDS
COVAL(NOZZLE,F:II:,ONLYMS,0.0)
ENDDO
3. Inlet Boundary-- free stream
PATCH(UPSTFREE,LOW,1,NX,NY/2+1,NY,1,1,1,1)
COVAL(UPSTFREE,P1,FIXFLU,RHO1*WFREE);COVAL(UPSTFREE,W1,ONLYMS,WFREE)
COVAL(UPSTFREE,CONC,ONLYMS,CFREE)
DO II = 1,NFLUIDS
COVAL(UPSTFREE,F:II:,ONLYMS,0.0)
ENDDO
COVAL(FREE,F:NFLUIDS:,ONLYMS,1.0)
the following statements dictate that there is a source of MIXL,
per unit volume, equal to LENCON*MNSQ
PATCH($MNSQ,VOLUME,1,NX,1,NY,1,NZ,1,1)
COVAL($MNSQ, MIXL, LENCON*1.0E-5, 1.E5)
GROUP 14. Downstream pressure for PARAB=T
IPARAB=1 ! pressure is fixed; and, since no PBAR setting is made
! longitudinal pressure gradients are absent
GROUP 16. Termination of iterations
SELREF=T;RESFAC=1.E-2;LITHYD=25
GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,10.0);RELAX(W1,FALSDT,10.0)
relax(u1,falsdt,1.e-1)
RELAX(CONC,FALSDT,10.0)
relax(p1,linrlx,0.1)
GROUP 18. Limits on variables or increments to them
VARMIN(V1)=-1.E3;VARMAX(V1)=1.E3
GROUP 19. Data communicated by SATELLITE to GROUND
SPEDAT(MFM,MFMMOD, C,:MFMMOD:)
SPEDAT(MFM,NFLUIDS,I,:NFLUIDS:)
SPEDAT(MFM,NFLR, I,:NFLR:)
SPEDAT(MFM,NFLF, I,:NFLF:)
SPEDAT(MFM,CONMIX, R,:CONMIX:)
SPEDAT(MFM,VISCON, R,:VISCON:)
SPEDAT(MFM,POPMIN, R,:WFREE:)
SPEDAT(MFM,POPMAX, R,:WJET:)
GROUP 22. Monitor print-out
IZMON=NZ/2;IYMON=1; ITABL=1;NPLT=1;IPLTL=LITHYD
TSTSWP=-5;NYPRIN=ny/5
GROUP 23. Field print-out and plot control
ORSIZ=0.4
PATCH(IZEQNZ,PROFIL,1,1,1,NY,NZ,NZ,1,1) ! final cross-stream
PLOT(IZEQNZ,W1,0.0,0.0); PLOT(IZEQNZ,CONC,0.0,0.0)
PLOT(IZEQNZ,ENUT,0.0,0.0)
NZPRIN=NZ
----------------------------------------- ! longitudinal
PATCH(MIDDLE,PROFIL,1,nx,NY/2,NY/2,1,NZ,1,1)
COVAL(MIDDLE,F1,0.0,0.0); COVAL(MIDDLE,F4,0.0,0.0)
COVAL(MIDDLE,F7,0.0,0.0); COVAL(MIDDLE,F10,0.0,0.0)
CHAR(NAMPROF)
DO II=1,NFLUIDS
NAMPROF=PROF:II:
PATCH(:NAMPROF:,PROFIL,1,1,1,NY,1,NZ,1,1)
COVAL(:NAMPROF:,F:II:,0.0,0.0)
ENDDO
PATCH(FINAL,PROFIL,1,1,1,NY,NZ,NZ,1,1)
COVAL(FINAL,W1,0.0,0.0); COVAL(FINAL,MNSQ,0.0,0.0)
COVAL(FINAL,MIXL,0.0,0.0); COVAL(FINAL,ENUT,0.0,0.0)
ORSIZ=ORSIZ/2
IDISPA=1 ! dump for photon plot at each z step
GROUP 24. Dumps for restarts
late modifications
TSTSWP=-1; LITHYD=10; IYMON=1; IZMON=1; RESFAC=1.E-3
STOP