```
TALK=T;RUN(1,1)
PHOTON USE
AUTOPLOT
file
phi 5

da 1
p1
screen
msg pressure profile
pl 1
msg press  to continue
pause
cl
da 1
1u 1 19
screen
msg gas velocity
pl 1
msg press  to continue
pause
cl
da 1
gas
screen
msg gas volume fraction
pl 1
msg press  to continue
pause
cl
da 1
h2
screen
msg liquid enthalpy
pl 1
msg press e to END
enduse
DISPLAY
A liquid flows into the bottom of a vertical pipe at the
saturation temperature. Heat is added in the upper half of the
pipe,causing steam to be formed. There is a small amount of
interphase friction. The treatment is one-dimensional.

Despite the simplicity of its formulation, this is a difficult
phenomenon to simulate. For cases of this type, a steady-state
solution may not exist for all conditions; unstable or
oscillatory boiling is often observed in reality.
It is interesting to explore how sensitive the solution is to the
various relaxation parameters, and see whether the same values
are satisfactory also for different values of heat flux, density
ratio, gravitational acceleration, interphase-friction law, etc
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************

TEXT(1D STEADY BOILING IN PIPE:          W350)

************************************************************
************************************************************

IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 20
XULAST =1.
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 = 1
YVLAST =1.
YFRAC(1)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 1
ZWLAST =1.
ZFRAC(1)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = F
EQUVEL = F
NAME(1)=P1 ;NAME(3)=1U
NAME(4)=2U ;NAME(9)=GAS
NAME(10)=LIQ ;NAME(14)=H1
NAME(15)=H2 ;NAME(150)=MDOT
* 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(1U,Y,Y,N,N,N,Y)
SOLUTN(2U,Y,Y,N,N,N,Y)
SOLUTN(GAS,Y,Y,N,N,N,Y)
SOLUTN(LIQ,Y,Y,N,N,N,Y)
SOLUTN(H1,Y,Y,N,N,N,Y)
SOLUTN(H2,Y,Y,N,N,N,Y)
SOLUTN(MDOT,Y,N,N,N,N,Y)
INTMDT = 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(1U,Y,Y,N,Y,Y,Y)
TERMS(2U,Y,Y,N,Y,N,Y)
TERMS(GAS,Y,Y,N,Y,Y,Y)
TERMS(LIQ,Y,Y,N,Y,N,Y)
TERMS(H1,N,Y,N,Y,Y,Y)
TERMS(H2,N,Y,N,Y,N,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
EQDVDP = F
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.
RHO2 =1000. ;TMP2 =0.
EL2 =0.
DVO2DT =0. ;DRH2DP =0.
ENUL =1.0E-05 ;ENUT =0.
PHINT(1U)=-2.022E+04 ;PHINT(2U)=-2.022E+04
PHINT(GAS)=-2.022E+04 ;PHINT(LIQ)=-2.022E+04
PHINT(H1)=2.5E+06 ;PHINT(H2)=4.0E+05
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
CFIPS =-2500.
RLOLIM =1.0E-02 ;CMDOT = HEATBL
CINT(P1)=1. ;CINT(1U)=1.
CINT(2U)=1. ;CINT(GAS)=1.
CINT(LIQ)=1. ;CINT(H1)=1.0E-02
CINT(H2)=1.
PHINT(P1)=-2.022E+04 ;PHINT(1U)=-2.022E+04
PHINT(2U)=-2.022E+04 ;PHINT(GAS)=-2.022E+04
PHINT(LIQ)=-2.022E+04 ;PHINT(H1)=2.5E+06
PHINT(H2)=4.0E+05
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(1U)=1.0E-02
FIINIT(2U)=1.0E-02 ;FIINIT(GAS)=1.0E-02
FIINIT(LIQ)=0.99 ;FIINIT(H1)=2.5E+06
FIINIT(H2)=4.0E+05 ;FIINIT(MDOT)=1.0E-10
No PATCHes yet used for this Group
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(INLET ,WEST , 1, 1, 1, 1, 1, 1, 1, 1)
COVAL(INLET ,1U ,0. ,0. )
COVAL(INLET ,2U ,0. ,1.0E-02 )
COVAL(INLET ,GAS , FIXFLU ,0. )
COVAL(INLET ,LIQ , FIXFLU ,10. )
COVAL(INLET ,H1 ,0. ,0. )
COVAL(INLET ,H2 ,0. ,4.0E+05 )

PATCH(OUTLET ,CELL , 20, 20, 1, 1, 1, 1, 1, 1)
COVAL(OUTLET ,1U ,0. ,0. )
COVAL(OUTLET ,2U ,0. ,0. )
COVAL(OUTLET ,GAS ,1. ,0. )
COVAL(OUTLET ,LIQ ,1000. ,0. )
COVAL(OUTLET ,H1 ,0. , SAME )
COVAL(OUTLET ,H2 ,0. , SAME )

PATCH(GRAVITY ,PHASEM, 1, 20, 1, 1, 1, 1, 1, 1)
COVAL(GRAVITY ,1U , FIXFLU ,-9.81 )
COVAL(GRAVITY ,2U , FIXFLU ,-9.81 )

PATCH(HEATER ,PHASEM, 6, 20, 1, 1, 1, 1, 1, 1)
COVAL(HEATER ,H2 , FIXFLU ,1.0E+05 )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 200 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
ISWR1 = 1 ;ISWR2 = 10000
SELREF = T
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(1U)=10
LITER(2U)=10 ;LITER(GAS)=1
LITER(LIQ)=1 ;LITER(H1)=20
LITER(H2)=20
ENDIT(P1)=1.0E-03 ;ENDIT(1U)=1.0E-03
ENDIT(2U)=1.0E-03 ;ENDIT(GAS)=1.0E-03
ENDIT(LIQ)=1.0E-03 ;ENDIT(H1)=1.0E-03
ENDIT(H2)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(1U,FALSDT,0.1)
RELAX(2U,FALSDT,0.1)
RELAX(GAS,LINRLX,0.1)
RELAX(LIQ,LINRLX,0.1)
RELAX(H1,FALSDT,1.)
RELAX(H2,FALSDT,1.)
RELAX(MDOT,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(1U)=1.0E+06 ;VARMIN(1U)=-1.0E+06
VARMAX(2U)=1.0E+06 ;VARMIN(2U)=-1.0E+06
VARMAX(GAS)=1. ;VARMIN(GAS)=1.0E-06
VARMAX(LIQ)=1. ;VARMIN(LIQ)=1.0E-06
VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
VARMAX(H2)=1.0E+10 ;VARMIN(H2)=-1.0E+10
VARMAX(MDOT)=1.0E+10 ;VARMIN(MDOT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG62 = 1
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,N,N,N,Y,Y,Y)
OUTPUT(1U,N,N,N,Y,Y,Y)
OUTPUT(2U,N,N,N,Y,Y,Y)
OUTPUT(GAS,N,N,N,Y,Y,Y)
OUTPUT(LIQ,N,N,N,Y,Y,Y)
OUTPUT(H1,N,N,N,Y,N,N)
OUTPUT(H2,N,N,N,Y,N,N)
OUTPUT(MDOT,N,N,N,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 10 ;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
NXPRIN = 1 ;IXPRF = 1 ;IXPRL = 10000
IPLTF = 1 ;IPLTL = -1 ;NPLT = 10
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 2 ;IPROF = 3
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20

PATCH(LONGPLOT,PROFIL, 1, 20, 1, 1, 1, 1, 1, 1)
PLOT(LONGPLOT,1U ,0. ,2. )
PLOT(LONGPLOT,2U ,0. ,2. )
PLOT(LONGPLOT,LIQ ,0. ,1. )
PLOT(LONGPLOT,MDOT,0. ,0.3 )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP
```