GROUP 1. Run title and other preliminaries
TEXT(0D UNSTEADY BOILING IN VESSEL:      W410
TITLE
  DISPLAY
  The boiling of water in an enclosed vessel is considered.
  Heat is added at a fixed rate to the liquid, and the liquid
  gradually boils away. The latent heat is taken as constant,
  which in view of the enormous increase in pressure with time
  is physically incorrect. This defect is made good in the
  following case: comparison of the two cases is interesting.
  Only one control cell is used, ie. no velocities are solved.
  This is done to focus attention on the formulation of
  interphase transport in EARTH, & especially the deduction of
  the rate of evaporation from the net heat transferred to the
  interface between the phases.
  ENDDIS
    GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,80,15.0,1.5)
    GROUP 7. Variables stored, solved & named
ONEPHS=F;SOLVE(P1,R1,R2,H1,H2);NAME(R1)=GAS;NAME(R2)=LIQ
STORE(RHO1,MDOT)
    GROUP 8. Terms (in differential equations) & devices
HUNIT=1.0E-05
    GROUP 9. Properties of the medium (or media)
RHO2=1000.0;DRH1DP=COMPRESS;ENUL=1.0E-05
RHO1=COMPRESS;RHO1A=9.22E-06;RHO1B=0.959;PRESS0=1.0E+05;RHO1C=0.0
    GROUP 10. Inter-phase-transfer processes and properties
CFIPS=1.0E+00;CMDOT=HEATBL;CINT(H1)=1.0E+02;CINT(H2)=1.0E+02
PHINT(H2)=5.0;PHINT(H1)=2.6E+01
    GROUP 11. Initialization of variable or porosity fields
FIINIT(GAS)=0.5;FIINIT(LIQ)=0.5
FIINIT(P1)=0.0E+00;FIINIT(H1)=26.0;FIINIT(H2)=5.0
    GROUP 13. Boundary conditions and special sources
   ** Heat input....
PATCH(HEAT,CELL,1,1,1,1,1,1,1,LSTEP);COVAL(HEAT,H2,FIXFLU,5.0E+02)
    GROUP 15. Termination of sweeps
LSWEEP=10;VARMIN(DEN1)=0.1;VARMIN(P1)=0.0
RELAX(MDOT,LINRLX,1.0);RELAX(H1,FALSDT,1.0);RELAX(H2,FALSDT,1.0)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 23. Field print-out and plot control
PATCH(PLOT,PROFIL,1,1,1,1,1,1,1,LSTEP)
PLOT(PLOT,P1,0.0,0.0);PLOT(PLOT,MDOT,0.0,0.0)
PLOT(PLOT,LIQ,0.0,1.0);PLOT(PLOT,GAS,0.0,1.0);NTPRIN=LSTEP/5
TSTSWP=-1;nprmnt=5