TALK=T;RUN(1,1) PHOTON USE p gr ou z 1;gr ou y 1 x 17 22 col 15;gr ou y m x 17 22 col 15 msg red colour outlines the Heat Source msg Press RETURN to continue pause msg Velocity vectors for the steam vec z 1 sh msg Pressto continue pause vec cl;red set vec ph 2 msg Velocity vectors for the water droplets vec z 1 sh msg Press to continue pause vec cl;red msg Contours of pressure con p1 z 1 fi;.01 msg Press to continue pause con cl;red msg Contours of mass transfer rate con mdot z 1 fi;.01 con mdot z 1;va;1;0 msg Press to continue pause con cl;red msg Volume fraction of water droplets con r2 z 1 fi;.001 msg Press e to END enduse DISPLAY BOILING DROPLETS IN A DUCT 2-dimensional (x-y), BFC, steady, elliptic simulation In this case the saturated water droplets enter the S-shaped duct, and then boil to produce steam. | | | | heat source | | / | \___|__ \ \ --------\ \ sat. ________/ / --> / water ---------/ enddis ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Boiling Droplets in a Duct ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 14 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 33 XULAST =0.34 ************************************************************ Group 4. Y-Direction Grid Spacing NY = 8 YVLAST =0.03 ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST =0.1 ************************************************************ Group 6. Body-Fitted Coordinates BFC = T ;NONORT = T NCRT = -1 RSTGEO = F ;SAVGEO = F UUP = T ;VUP = F ;WUP = F NGEOM =CHAM NAMXYZ =CHAM ANGMIN =20. DOMAIN( 1, 34, 1, 9, 1, 2) * Set fixed sub-domain FIXDOM( 1, 0, 0, 0, 0, 0, 0) FIXDOM( 2, 0, 0, 0, 0, 0, 0) FIXDOM( 3, 0, 0, 0, 0, 0, 0) FIXDOM( 4, 0, 0, 0, 0, 0, 0) FIXDOM( 5, 0, 0, 0, 0, 0, 0) FIXDOM( 6, 0, 0, 0, 0, 0, 0) FIXDOM( 7, 0, 0, 0, 0, 0, 0) FIXDOM( 8, 0, 0, 0, 0, 0, 0) FIXDOM( 9, 0, 0, 0, 0, 0, 0) FIXDOM( 10, 0, 0, 0, 0, 0, 0) LIJ = F ;LJK = T ;LIK = F ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = F EQUVEL = F NAME(1)=P1 ;NAME(3)=U1 NAME(4)=U2 ;NAME(5)=V1 NAME(6)=V2 ;NAME(9)=R1 NAME(10)=R2 ;NAME(14)=H1 NAME(15)=H2 ;NAME(141)=TMP2 NAME(142)=TMP1 ;NAME(143)=MDOT NAME(144)=CFIP ;NAME(145)=W2CR NAME(146)=V2CR ;NAME(147)=U2CR NAME(148)=WCRT ;NAME(149)=VCRT NAME(150)=UCRT * 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,Y,N,N,N) SOLUTN(U1,Y,Y,N,N,N,Y) SOLUTN(U2,Y,Y,N,N,N,Y) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(V2,Y,Y,N,N,N,Y) SOLUTN(R2,Y,Y,N,N,N,Y) SOLUTN(H1,Y,Y,N,N,N,Y) SOLUTN(H2,Y,Y,N,N,N,Y) SOLUTN(TMP2,Y,N,N,N,N,Y) SOLUTN(TMP1,Y,N,N,N,N,Y) SOLUTN(MDOT,Y,N,N,N,N,Y) SOLUTN(CFIP,Y,N,N,N,N,Y) SOLUTN(W2CR,Y,N,N,N,N,N) SOLUTN(V2CR,Y,N,N,N,N,N) SOLUTN(U2CR,Y,N,N,N,N,N) SOLUTN(WCRT,Y,N,N,N,N,N) SOLUTN(VCRT,Y,N,N,N,N,N) SOLUTN(UCRT,Y,N,N,N,N,N) INTMDT = 143 INTFRC = 144 TEMP1 = 142 TEMP2 = 141 ************************************************************ 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(U1,Y,Y,Y,Y,Y,Y) TERMS(U2,Y,Y,Y,Y,N,Y) TERMS(V1,Y,Y,Y,Y,Y,Y) TERMS(V2,Y,Y,Y,Y,N,Y) TERMS(R2,Y,Y,Y,Y,N,Y) TERMS(H1,Y,Y,Y,Y,Y,Y) TERMS(H2,Y,Y,Y,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 =0.5 ;TMP1 = GRND2 ;EL1 =0. TSURR =0. ;TEMP0 =0. ;PRESS0 =0. DVO1DT =0.5 ;DRH1DP =0. TMP1A =-273. ;TMP1B =0. ;TMP1C =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. RHO2 =950. ;TMP2 = GRND2 ;EL2 =0. DVO2DT =0. ;DRH2DP =0. TMP2A =-273. ;TMP2B =0. ;TMP2C =0. ENUL =1.25E-04 ;ENUT =0.0125 PHINT(U1)=-2.022E+04 ;PHINT(U2)=-2.022E+04 PHINT(V1)=-2.022E+04 ;PHINT(V2)=-2.022E+04 PHINT(R2)=-2.022E+04 ;PHINT(H1)=2.5E+06 PHINT(H2)=4.0E+05 PRNDTL(U1)=1. ;PRNDTL(U2)=1. PRNDTL(V1)=1. ;PRNDTL(V2)=1. PRNDTL(R2)=1. ;PRNDTL(H1)=0.97 PRNDTL(H2)=1.73 PRT(U1)=1. ;PRT(U2)=1. PRT(V1)=1. ;PRT(V2)=1. PRT(R2)=1. ;PRT(H1)=1. PRT(H2)=1. CP1 =6702. ;CP2 =1072. ************************************************************ Group 10.Inter-Phase Transfer Processes CFIPS =2.25E+07 RLOLIM =0. ;CMDOT = HEATBL CINT(P1)=1. ;CINT(U1)=1. CINT(U2)=1. ;CINT(V1)=1. CINT(V2)=1. ;CINT(R2)=1. CINT(H1)=0.687285 ;CINT(H2)=0.385356 PHINT(P1)=-2.022E+04 ;PHINT(U1)=-2.022E+04 PHINT(U2)=-2.022E+04 ;PHINT(V1)=-2.022E+04 PHINT(V2)=-2.022E+04 ;PHINT(R2)=-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(U1)=10. FIINIT(U2)=10. ;FIINIT(V1)=1.0E-10 FIINIT(V2)=1.0E-10 ;FIINIT(R1)=0.95 FIINIT(R2)=0.05 ;FIINIT(H1)=2.5E+06 FIINIT(H2)=3.99E+05 ;FIINIT(TMP2)=1.0E-10 FIINIT(TMP1)=1.0E-10 ;FIINIT(MDOT)=1.0E-10 FIINIT(CFIP)=1.0E-10 ;FIINIT(W2CR)=1.0E-10 FIINIT(V2CR)=1.0E-10 ;FIINIT(U2CR)=1.0E-10 FIINIT(WCRT)=1.0E-10 ;FIINIT(VCRT)=1.0E-10 FIINIT(UCRT)=1.0E-10 No PATCHes yet used for this Group INIADD = F 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, 8, 1, 1, 1, 1) COVAL(INLET ,U1 ,0. ,10. ) COVAL(INLET ,U2 ,0. ,10. ) COVAL(INLET ,V1 ,0. ,0. ) COVAL(INLET ,V2 ,0. ,0. ) COVAL(INLET ,R1 , FIXFLU ,4.75 ) COVAL(INLET ,R2 , FIXFLU ,475.000122 ) COVAL(INLET ,H1 ,0. ,2.5E+06 ) COVAL(INLET ,H2 ,0. ,3.99E+05 ) PATCH(OUTLET ,EAST , 33, 33, 1, 8, 1, 1, 1, 1) COVAL(OUTLET ,V1 ,0. ,0. ) COVAL(OUTLET ,V2 ,0. ,0. ) COVAL(OUTLET ,R1 ,500. ,0. ) COVAL(OUTLET ,R2 ,9.5E+05 ,0. ) COVAL(OUTLET ,H1 ,0. , SAME ) COVAL(OUTLET ,H2 ,0. , SAME ) PATCH(TOP1 ,NWALL , 1, 16, 8, 8, 1, 1, 1, 1) COVAL(TOP1 ,U1 , GRND2 ,0. ) COVAL(TOP1 ,H1 , GRND2 ,2.49E+06 ) PATCH(TOP2 ,NWALL , 23, 33, 8, 8, 1, 1, 1, 1) COVAL(TOP2 ,U1 , GRND2 ,0. ) COVAL(TOP2 ,H1 , GRND2 ,2.4999E+06 ) PATCH(BOT1 ,SWALL , 1, 16, 1, 1, 1, 1, 1, 1) COVAL(BOT1 ,U1 , GRND2 ,0. ) COVAL(BOT1 ,H1 , GRND2 ,2.49E+06 ) PATCH(BOT2 ,SWALL , 23, 33, 1, 1, 1, 1, 1, 1) COVAL(BOT2 ,U1 , GRND2 ,0. ) COVAL(BOT2 ,H1 , GRND2 ,2.4999E+06 ) PATCH(HEAT1 ,NORTH , 17, 22, 8, 8, 1, 1, 1, 1) COVAL(HEAT1 ,H1 , FIXFLU ,3.0E+06 ) PATCH(HEAT2 ,SOUTH , 17, 22, 1, 1, 1, 1, 1, 1) COVAL(HEAT2 ,H1 , FIXFLU ,3.0E+06 ) PATCH(BUOYANCY,PHASEM, 1, 33, 1, 8, 1, 1, 1, 1) COVAL(BUOYANCY,U1 , FIXFLU , GRND2 ) COVAL(BUOYANCY,U2 , FIXFLU , GRND2 ) COVAL(BUOYANCY,V1 , FIXFLU , GRND2 ) COVAL(BUOYANCY,V2 , FIXFLU , GRND2 ) XCYCLE = F EGWF = T WALLCO = GRND2 BUOYA =0. ; BUOYB =-9.81 ; BUOYC =0. BUOYD =0.5 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 150 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 ISWR1 = 1 ;ISWR2 = 10000 SELREF = T RESFAC =1.0E-02 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(U1)=10 LITER(U2)=10 ;LITER(V1)=10 LITER(V2)=10 ;LITER(R2)=1 LITER(H1)=20 ;LITER(H2)=20 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(U2)=1.0E-03 ;ENDIT(V1)=1.0E-03 ENDIT(V2)=1.0E-03 ;ENDIT(R2)=1.0E-03 ENDIT(H1)=1.0E-03 ;ENDIT(H2)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,0.7) RELAX(U1,FALSDT,1.666667E-04) RELAX(U2,FALSDT,1.666667E-04) RELAX(V1,FALSDT,1.666667E-04) RELAX(V2,FALSDT,1.666667E-04) RELAX(R1,LINRLX,0.2) RELAX(R2,LINRLX,0.2) RELAX(H1,FALSDT,1.666667E-04) RELAX(H2,FALSDT,1.666667E-04) RELAX(TMP2,LINRLX,1.) RELAX(TMP1,LINRLX,1.) RELAX(MDOT,LINRLX,0.1) RELAX(CFIP,LINRLX,1.) RELAX(W2CR,LINRLX,1.) RELAX(V2CR,LINRLX,1.) RELAX(U2CR,LINRLX,1.) RELAX(WCRT,LINRLX,1.) RELAX(VCRT,LINRLX,1.) RELAX(UCRT,LINRLX,1.) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06 VARMAX(U2)=1.0E+06 ;VARMIN(U2)=-1.0E+06 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06 VARMAX(V2)=1.0E+06 ;VARMIN(V2)=-1.0E+06 VARMAX(R1)=1. ;VARMIN(R1)=1.0E-06 VARMAX(R2)=1. ;VARMIN(R2)=1.0E-06 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10 VARMAX(H2)=1.0E+10 ;VARMIN(H2)=-1.0E+10 VARMAX(TMP2)=1.0E+10 ;VARMIN(TMP2)=-1.0E+10 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10 VARMAX(MDOT)=1.0E+10 ;VARMIN(MDOT)=-1.0E+10 VARMAX(CFIP)=1.0E+10 ;VARMIN(CFIP)=-1.0E+10 VARMAX(W2CR)=1.0E+10 ;VARMIN(W2CR)=-1.0E+10 VARMAX(V2CR)=1.0E+10 ;VARMIN(V2CR)=-1.0E+10 VARMAX(U2CR)=1.0E+10 ;VARMIN(U2CR)=-1.0E+10 VARMAX(WCRT)=1.0E+10 ;VARMIN(WCRT)=-1.0E+10 VARMAX(VCRT)=1.0E+10 ;VARMIN(VCRT)=-1.0E+10 VARMAX(UCRT)=1.0E+10 ;VARMIN(UCRT)=-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,Y,N,Y,Y,Y,Y) OUTPUT(U1,Y,N,Y,Y,Y,Y) OUTPUT(U2,Y,N,Y,Y,Y,Y) OUTPUT(V1,Y,N,Y,Y,Y,Y) OUTPUT(V2,Y,N,Y,Y,Y,Y) OUTPUT(R1,Y,N,Y,N,N,N) OUTPUT(R2,Y,N,Y,Y,Y,Y) OUTPUT(H1,Y,N,Y,Y,Y,Y) OUTPUT(H2,Y,N,Y,Y,Y,Y) OUTPUT(TMP2,Y,N,Y,N,N,N) OUTPUT(TMP1,Y,N,Y,N,N,N) OUTPUT(MDOT,Y,N,Y,N,N,N) OUTPUT(CFIP,Y,N,Y,N,N,N) OUTPUT(W2CR,Y,N,Y,N,N,N) OUTPUT(V2CR,Y,N,Y,N,N,N) OUTPUT(U2CR,Y,N,Y,N,N,N) OUTPUT(WCRT,Y,N,Y,N,N,N) OUTPUT(VCRT,Y,N,Y,N,N,N) OUTPUT(UCRT,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 24 ;IYMON = 4 ;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 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 No PATCHes yet used for this Group ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP