TALK=T;RUN( 1, 1)
    GROUP 1. Run title and other preliminaries
TEXT(2D STEADY BOILING IN PIPE CFIPS+ CINT = GRND7)
  This case is a modified w370 for detecting relaxation in e3cmp2
    GROUP 3. X-direction grid specification
load(w370)
  ** Size should be filled by CFIPS=GRND7
store(CFIP,size)
 
    GROUP 9. Properties of the medium (or media)
  ** Switch of dissusion of R1,R2
PRNDTL(R1)=1.E10;PRNDTL(R2)=1.E10
  ** Densities
RHO1=750.;RHO2=35.
    GROUP 10. Interphase-transfer processes and properties
 
  ** Set a non-constant inter-phase friction coefficient and
     activate the computation of inter-phase mass transfer
     by means of a heat balance.
CFIPS=grnd7;CMDOT=HEATBL
  ** particle size
cfipb=5.e-5
  ** minimum slip velocity
cfipa=0.0001
   ** Set the inter-phase heat transfer coefficients.
        CINT(H1)=1.E1;CINT(H2)=1.E1
   ** Heat transfer to spheres
cint(h1)=grnd7;cint(h2)=1.e15
   ** particle size
        cinh2c=5.e-5
 
  ** Set the enthalpy values at the interface
  ** Saturation values < 0 liuid > 1.513E+6 steam
PHINT(H1)=0.;PHINT(H2)=1.513E+6
    GROUP 11. Initialization of variable or porosity fields
  ** Incoming steam and water set to same velocity
FIINIT(R1)=0.8
FIINIT(R2)=0.2
FIINIT(H1)=0.0;FIINIT(H2)=1.513E6;FIINIT(MDOT)=0.0
    GROUP 13. Boundary conditions and special sources
  ** Inlet at bottom
PATCH(INLET,LOW,#1,#1,#1,#NREGY,#1,#1,1,1)
COVAL(INLET,P1,FIXFLU,2.*rho1*0.8)
COVAL(INLET,W1,ONLYMS,2.0)
COVAL(INLET,P2,FIXFLU,2.*rho2*0.2)
COVAL(INLET,W2,ONLYMS,2.0)
  ** enthalpies between saturation values
COVAL(INLET,H1,ONLYMS,0.)
COVAL(INLET,H2,ONLYMS,phint(h2))
 
  ** Outlet
PATCH(OUTLET,HIGH,#1,#1,#1,#NREGY,#NREGZ,#NREGZ,1,1)
COVAL(OUTLET,P1,RHO1*FIXP,0.0);COVAL(OUTLET,P2,RHO2*FIXP,0.0)
COVAL(OUTLET,W1,ONLYMS,0.0);COVAL(OUTLET,W2,ONLYMS,0.0)
COVAL(OUTLET,H1,ONLYMS,SAME);COVAL(OUTLET,H2,ONLYMS,SAME)
    GROUP 15. Termination of sweeps
 
LSWEEP=300
RESREF(P1)=1.E-5;RESREF(V1)=1.E-4
RESREF(V2)=1.E-4;RESREF(W1)=1.E-4
RESREF(W1)=1.E-4;RESREF(H1)=1.E-4
RESREF(H1)=1.E-4;RESREF(R1)=1.E-4
RESREF(R1)=1.E-4
resfac=1e-6
  ** default library case relaxations
 
    GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,0.001);RELAX(V2,FALSDT,0.001)
RELAX(W1,FALSDT,0.001);RELAX(W2,FALSDT,0.001)
RELAX(R1,LINRLX,0.25);RELAX(R2,LINRLX,0.25)
RELAX(MDOT,LINRLX,0.1);RELAX(H1,FALSDT,1.0)
RELAX(H2,FALSDT,1.0)
RELAX(CFIP,LINRLX,0.5)
    GROUP 22. Spot-value print-out
IYMON=5;IZMON=1
    GROUP 23. Field print-out and plot control
output(mdot,p,p,p,p,y,y)
ICHR=1;IPLTL=100;NPLT=2;TSTSWP=-1
PATCH(CONT,CONTUR,1,1,1,NY,1,NZ,1,1)
PLOT(CONT,R1,0.0,10.0);PLOT(CONT,W1,0.0,10.0)
PLOT(CONT,W2,0.0,10.0)
 LIBREF=372
VARMIN(R1)=1E-15; VARMIN(R2)=1E-15
LSWEEP=500