PHOTON USE
  p
  phi
 
 
 
  msg grid
  gr x 1;gr x m
  gr y 1;gr y m
  gr z 1;gr z m
  msg Press  to continue
  pause
  cl
  vi x
  up z
  scale
  cl
  msg 1st-phase (light phase) velocity vectors on plane ix=1
  set vec ref;4
  vec x 1 sh
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg 2nd-phase (dense phase) velocity vectors on plane ix=1
  set vec ref;4
  set vec comp
  - v2 w2
  vec x 1 sh
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of pressure on plane ix=1
  con p1 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 1st-phase (steam) volume fraction on plane ix=1
  con stea x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of 2nd-phase (water) volume fraction on plane ix=1
  con liq x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of steam enthalpy on plane ix=1
  con h1 x 1 fi;0.002
  gr ou x 1
  msg Press  to continue
  pause
  cl
  msg contour of water enthalpy on plane ix=1
  con h2 x 1 fi;0.002
  gr ou x 1
  msg Press e to END
  enduse
 
    GROUP 1. Run title
TEXT(3D STEAM GENERATOR :                W802
TITLE
MESG(PC486/50 time last reported as appx. 1.min 45.sec
REAL(PORE,PORN,PORH,PORC,ACO,BCO,DCO,ECO)
ACO=-100.0;BCO=-150.0;DCO=-4.0;ECO=-6.0
PORE=-0.4395;PORN=-0.4395;PORH=-0.4395;PORC=-0.4395
    GROUP 3. X-direction grid specification
CARTES=F;GRDPWR(X,6,3.14159,1.0)
    GROUP 4. Y-direction grid specification
NY=7;YVLAST=2.3;YFRAC(1)=0.13;YFRAC(2)=0.261
YFRAC(3)=0.3913;YFRAC(4)=0.5217
YFRAC(5)=0.6522;YFRAC(6)=0.7826;YFRAC(7)=1.0
    GROUP 5. Z-direction grid specification
GRDPWR(Z,12,10.0,1.0)
    GROUP 7. Variables stored, solved & named
ONEPHS=F;SOLVE(P1,U1,U2,V1,V2,W1,W2,R1,R2,H1,H2)
NAME(R1)=STEA;NAME(R2)=LIQ
  ** Provide storage for inter-phase mass transfer
INTMDT=20;NAME(INTMDT)=MDOT;SOLUTN(INTMDT,Y,N,N,N,N,N)
  ** Provide storage for porosities
EPOR=21;NAME(21)=EPOR;NPOR=22;NAME(22)=NPOR
HPOR=23;NAME(23)=HPOR;VPOR=24;NAME(24)=VPOR
STORE(EPOR,NPOR,VPOR,HPOR)
    GROUP 8. Terms (in differential equations) & devices
DIFCUT=0.0
    GROUP 9. Properties of the medium (or media)
PRNDTL(STEA)=1.E10;PRNDTL(LIQ)=1.E10
PRNDTL(U2)=1.E10;PRNDTL(V2)=1.E10
PRNDTL(W2)=1.E10;PRNDTL(H2)=1.E10
RHO1=38.82;RHO2=733.2;ENUL=0.5/RHO1
    GROUP 10. Inter-phase-transfer processes and properties
  **Set a constant inter-phase friction coefficient and
    activate the computation of inter-phase mass transfer
    by means of a heat balance.
CFIPS=1.E3;CMDOT=HEATBL;RLOLIM=0.001
  **Set the inter-phase heat transfer coefficients.
CINT(H1)=3.E4/CFIPS;CINT(H2)=3.E1/CFIPS
  **Set the enthalpy values at the interface.
PHINT(H1)=2.767E1;PHINT(H2)=1.286E1
    GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=1.0;FIINIT(W2)=1.0
FIINIT(LIQ)=0.8;FIINIT(STEA)=0.2
FIINIT(EPOR)=1.0;FIINIT(NPOR)=1.0
FIINIT(HPOR)=1.0;FIINIT(VPOR)=1.0
FIINIT(H1)=PHINT(H1);FIINIT(H2)=PHINT(H2)
FIINIT(MDOT)=0.0
  **Set field values in blocked region
INIADD=T
PATCH(INIT,INIVAL,1,NX,NY,NY,1,9,1,1)
INIT(INIT,W1,0.0,-1.0);INIT(INIT,W2,0.0,-1.0)
INIT(INIT,STEA,0.0,-0.25);INIT(INIT,LIQ,0.0,0.2)
INIT(INIT,H2,0.0,-1.286E1);INIT(INIT,H1,0.0,-2.767E1)
       **** East-Face Porosities ****
PATCH(E00,INIVAL,1,NX,NY,NY,1,9,1,1);INIT(E00,EPOR,0.0,-1.0)
PATCH(E7,INIVAL,1,NX,1,NY-1,1,7,1,1);INIT(E7,EPOR,0.0,PORE)
PATCH(E8,INIVAL,1,NX,1,4,8,8,1,1);INIT(E8,EPOR,0.0,PORE)
PATCH(E9,INIVAL,1,NX,1,2,9,9,1,1);INIT(E9,EPOR,0.0,PORE)
PATCH(E10,INIVAL,1,NX,NY,NY,10,10,1,1)
INIT(E10,EPOR,0.0,-0.775)
PATCH(E11,INIVAL,1,NX,NY,NY,11,11,1,1)
INIT(E11,EPOR,0.0,-0.4)
PATCH(E12,INIVAL,1,NX,NY,NY,12,12,1,1)
INIT(E12,EPOR,0.0,-0.03)
     ***** North Face Porosities *****
PATCH(N00,INIVAL,1,NX,NY-1,NY,1,9,1,1)
INIT(N00,NPOR,0.0,-1.0)
PATCH(N7,INIVAL,1,NX,1,NY-2,1,7,1,1)
INIT(N7,NPOR,0.0,PORN)
PATCH(N8,INIVAL,1,NX,1,3,8,8,1,1)
INIT(N8,NPOR,0.0,PORN)
PATCH(N9,INIVAL,1,NX,1,1,9,9,1,1)
INIT(N9,NPOR,0.0,PORN)
PATCH(N10,INIVAL,1,NX,NY,NY,10,NZ,1,1)
INIT(N10,NPOR,0.0,-1.0)
     **** High Face Porosities *****
PATCH(H00,INIVAL,1,NX,NY,NY,1,9,1,1)
INIT(H00,HPOR,0.0,-1.0)
PATCH(H7,INIVAL,1,NX,1,NY-1,1,7,1,1)
INIT(H7,HPOR,0.0,PORH)
PATCH(H8,INIVAL,1,NX,1,2,8,8,1,1)
INIT(H8,HPOR,0.0,PORH)
PATCH(H10,INIVAL,1,NX,NY,NY,10,10,1,1)
INIT(H10,HPOR,0.0,-0.55)
PATCH(H11,INIVAL,1,NX,NY,NY,11,11,1,1)
INIT(H11,HPOR,0.0,-0.1)
     **** Cell Volume Porosities ****
PATCH(V00,INIVAL,1,NX,NY,NY,1,9,1,1)
INIT(V00,VPOR,0.0,-1.0)
PATCH(V7,INIVAL,1,NX,1,NY-1,1,7,1,1)
INIT(V7,VPOR,0.0,PORC)
PATCH(V8,INIVAL,1,NX,1,4,8,8,1,1)
INIT(V8,VPOR,0.0,PORC)
PATCH(V9,INIVAL,1,NX,1,2,9,9,1,1)
INIT(V9,VPOR,0.0,PORC)
PATCH(V10,INIVAL,1,NX,NY,NY,10,10,1,1)
INIT(V10,VPOR,0.0,-0.775)
PATCH(V11,INIVAL,1,NX,NY,NY,11,11,1,1)
INIT(V11,VPOR,0.0,-0.4)
PATCH(V12,INIVAL,1,NX,NY,NY,12,12,1,1)
INIT(V12,VPOR,0.0,-0.03)
    GROUP 13. Boundary conditions and special sources
     **Inlet of liquid phase at bottom and radially inward
       at IY=NY-1
INLET(INLET,NORTH,1,NX,NY-1,NY-1,1,1,1,1)
VALUE(INLET,P2,500.)
VALUE(INLET,V2,-500./RHO2)
VALUE(INLET,H2,9.729)
  ** Heat Sources
PATCH(HEAT1,VOLUME,1,NX/2,1,6,1,7,1,1)
COVAL(HEAT1,H2,FIXFLU,1.2E3)
PATCH(HEAT2,VOLUME,1,NX/2,1,4,8,8,1,1)
COVAL(HEAT2,H2,FIXFLU,1.2E3)
PATCH(HEAT3,VOLUME,1,NX/2,1,2,9,9,1,1)
COVAL(HEAT3,H2,FIXFLU,1.2E3)
PATCH(HEAT4,VOLUME,NX/2+1,NX,1,6,1,7,1,1)
COVAL(HEAT4,H2,FIXFLU,6.0E2)
PATCH(HEAT5,VOLUME,NX/2+1,NX,1,4,8,8,1,1)
COVAL(HEAT5,H2,FIXFLU,6.0E2)
PATCH(HEAT6,VOLUME,NX/2+1,NX,1,2,9,9,1,1)
COVAL(HEAT6,H2,FIXFLU,6.0E2)
  ** Gravity
PATCH(GRAVITY,PHASEM,1,NX,1,NY,1,NZ,1,1)
COVAL(GRAVITY,W1,FIXFLU,-9.81)
COVAL(GRAVITY,W2,FIXFLU,-9.81)
  ** Outlet
PATCH(OUTLET,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(OUTLET,P2,RHO2,0.0)
COVAL(OUTLET,P1,RHO1,0.0)
COVAL(OUTLET,H1,ONLYMS,2.767E1)
COVAL(OUTLET,H2,ONLYMS,1.286E1)
      ** Resistance of Baffles **
PATCH(RBAF1,FREEH,1,NX,1,NY,2,2,1,1)
COVAL(RBAF1,W1,-0.5*RHO1,0.0)
COVAL(RBAF1,W2,-0.5*RHO2,0.0)
PATCH(RBAF2,FREEH,1,NX,1,NY,7,7,1,1)
COVAL(RBAF2,W1,-0.5*RHO1,0.0)
COVAL(RBAF2,W2,-0.5*RHO2,0.0)
PATCH(RBAF3,FREEH,1,NX,1,NY,11,11,1,1)
COVAL(RBAF3,W1,-0.5*RHO1,0.0)
COVAL(RBAF3,W2,-0.5*RHO2,0.0)
      ** Resistance of Tubes **
PATCH(RTUBE1,PHASEM,1,NX,1,6,1,7,1,1)
COVAL(RTUBE1,U1,ACO,0.0);COVAL(RTUBE1,U2,BCO,0.0)
COVAL(RTUBE1,V1,ACO,0.0);COVAL(RTUBE1,V2,BCO,0.0)
COVAL(RTUBE1,W1,DCO,0.0);COVAL(RTUBE1,W2,ECO,0.0)
PATCH(RTUBE2,PHASEM,1,NX,1,4,8,8,1,1)
COVAL(RTUBE2,U1,ACO,0.0);COVAL(RTUBE2,U2,BCO,0.0)
COVAL(RTUBE2,V1,ACO,0.0);COVAL(RTUBE2,V2,BCO,0.0)
COVAL(RTUBE2,W1,DCO,0.0);COVAL(RTUBE2,W2,ECO,0.0)
PATCH(RTUBE3,PHASEM,1,NX,1,2,9,9,1,1)
COVAL(RTUBE3,U1,ACO,0.0)
COVAL(RTUBE3,U2,BCO,0.0);COVAL(RTUBE3,V1,ACO,0.0)
COVAL(RTUBE3,V2,BCO,0.0);COVAL(RTUBE3,W1,DCO,0.0)
COVAL(RTUBE3,W2,ECO,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=25
RESREF(P1)=1.E-4;RESREF(V1)=1.E-4
RESREF(V2)=1.E-4;RESREF(W1)=1.E-4
RESREF(W2)=1.E-4;RESREF(H1)=1.E-4
RESREF(H2)=1.E-4;RESREF(STEA)=1.E-4
RESREF(LIQ)=1.E-4;RESREF(U1)=1.E-4
RESREF(U2)=1.E-4
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,1.0);RELAX(V1,FALSDT,0.05)
RELAX(V2,FALSDT,0.05);RELAX(W1,FALSDT,0.05)
RELAX(W2,FALSDT,0.05);RELAX(U1,FALSDT,0.05)
RELAX(U2,FALSDT,0.05);RELAX(STEA,LINRLX,0.5)
RELAX(LIQ,LINRLX,0.5);RELAX(MDOT,LINRLX,0.01)
RELAX(H1,FALSDT,1);RELAX(H2,FALSDT,1)
    GROUP 20. Preliminary print-out
INIFLD=F
    GROUP 21. Print-out of variables
OUTPUT(H1,Y,Y,Y,Y,Y,Y)
OUTPUT(H2,Y,Y,Y,Y,Y,Y)
OUTPUT(STEA,Y,Y,Y,Y,Y,Y)
OUTPUT(W1,Y,Y,Y,Y,Y,Y)
OUTPUT(W2,Y,Y,Y,Y,Y,Y)
OUTPUT(MDOT,Y,Y,Y,Y,Y,Y)
    GROUP 22. Spot-value print-out
TSTSWP=1;IYMON=3;IZMON=7;IXMON=3
    GROUP 23. Field print-out and plot control
ICHR=1;NXPRIN=1;IPLTL=LSWEEP;IXPRL=5
NROWCO=36;NCOLCO=28;IPROF=1
PATCH(PLOT,CONTUR,1,1,1,NY,1,NZ,1,1)
PLOT(PLOT,LIQ,0.0,10.0);PLOT(PLOT,W1,0.0,10.0)
PLOT(PLOT,W2,0.0,10.0);PLOT(PLOT,H2,0.0,10.0)
PATCH(COPLOT,CONTUR,NX,NX,1,NY,1,NZ,1,1)
PLOT(COPLOT,LIQ,0.0,10.0);PLOT(COPLOT,W1,0.0,10.0)
PLOT(COPLOT,W2,0.0,10.0);PLOT(COPLOT,H2,0.0,10.0)
PATCH(PRESS,PROFIL,1,1,2,2,1,NZ,1,1)
PLOT(PRESS,P1,0.0,0.0)
tstswp=-1