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