DISPLAY Library case 651: Tutorial: Boiling in a pipe (IPSA) ENDDIS TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 3.6, Date 02/11/04 CPVNAM=VDI;SPPNAM=Core ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Boiling in a pipe (IPSA) ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,1,10,20) * Cylindrical-polar grid CARTES=F ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = F * Non-default variable names NAME(146) =CMDO ; NAME(147) =EPKE NAME(148) =DEN1 ; NAME(149) =EL1 NAME(150) =ENUT * Solved variables list SOLVE(P1 ,V1 ,V2 ,W1 ,W2 ,R1 ,R2 ,H1 ) SOLVE(H2 ) * Stored variables list STORE(ENUT,EL1 ,DEN1,EPKE,CMDO) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,Y) SOLUTN(H1 ,Y,Y,Y,N,N,Y) SOLUTN(H2 ,Y,Y,Y,N,N,Y) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices TERMS (H1 ,N,Y,Y,Y,Y,Y) TERMS (H2 ,N,Y,Y,Y,N,Y) ************************************************************ Group 9. Properties PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02 RHO1 = 7.499400E+02 ENUL = 1.544000E-05 CP1 = 1.005000E+03 * Phase 2 material is: 67 WATER at 20. deg C SETPRPS(2, 67) RHO2 = 3.555000E+01 DVO1DT = 3.410000E-03 DVO2DT = 1.180000E-04 PHINT (H1 ) = 1.000000E-01 ;PHINT (H2 ) = 1.513000E+06 PRNDTL(H1 ) = -2.580000E-02 ;PRNDTL(H2 ) = -5.970000E-01 PRT (KE ) = 7.500000E-01 ;PRT (EP ) = 1.150000E+00 ************************************************************ Group 10.Inter-Phase Transfer Processes CFIPS = 2.500000E+03 RLOLIM = 0.000000E+00 ;CMDOT = HEATBL CINT (H1 ) = 1.000000E+01 ;CINT (H2 ) = 1.000000E+01 ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(W1 ) = 2.000000E+00 ;FIINIT(W2 ) = 2.000000E+00 FIINIT(R1 ) = 9.990000E-01 ;FIINIT(R2 ) = 1.000000E-04 FIINIT(H2 ) = 1.513000E+06 No PATCHes used for this Group INIADD = F ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources PATCH (BUOYANCY,PHASEM,0,0,0,0,0,0,1,1) COVAL (BUOYANCY,V1 , FIXFLU , GRND1 ) COVAL (BUOYANCY,W1 , FIXFLU , GRND1 ) BUOYA = 0.000000E+00 ; BUOYB = 0.000000E+00 BUOYC =-9.810000E+00 BUOYE = 2.000000E+01 EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 300 RESFAC = 1.000000E-03 ************************************************************ Group 16. Terminate Iterations LITER (P1 ) = 200 ;LITER (R1 ) = 1 LITER (R2 ) = 1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) RELAX(V1 ,FALSDT, 1.000000E-01) RELAX(V2 ,FALSDT, 1.000000E-01) RELAX(W1 ,FALSDT, 1.000000E-01) RELAX(W2 ,FALSDT, 1.000000E-01) RELAX(R1 ,LINRLX, 2.500000E-01) RELAX(R2 ,LINRLX, 2.500000E-01) RELAX(KE ,LINRLX, 5.000000E-01) RELAX(EP ,LINRLX, 5.000000E-01) RELAX(H1 ,FALSDT, 1.000000E+01) RELAX(H2 ,FALSDT, 1.000000E+01) RELAX(CMDO,LINRLX, 1.000000E-01) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(V2 ) = 1.000000E+06 ;VARMIN(V2 ) =-1.000000E+06 VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06 VARMAX(W2 ) = 1.000000E+06 ;VARMIN(W2 ) =-1.000000E+06 VARMAX(R1 ) = 1.000000E+00 ;VARMIN(R1 ) = 1.000000E-06 VARMAX(R2 ) = 1.000000E+00 ;VARMIN(R2 ) = 1.000000E-06 VARMAX(H1 ) = 1.000000E+10 ;VARMIN(H1 ) = 1.000000E-10 VARMAX(H2 ) = 1.000000E+10 ;VARMIN(H2 ) = 0.000000E+00 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T GENK = T ASAP = T PARSOL = T IENUTA = 2 SPEDAT(SET,FACETDAT,NUMOBJ,I,4) SPEDAT(SET,KECONST,C1E,R,1.44000E+00) SPEDAT(SET,KECONST,C2E,R,1.92000E+00) ************************************************************ Group 20. Preliminary Printout ECHO = T NULLPR = F NDST = 0 DSTTOL = 1.000000E-02 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 10 ;IZMON = 20 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts GVIEW(P,-9.840119E-01,-5.525740E-02,1.693140E-01) GVIEW(UP,1.686971E-01,1.567861E-02,9.855433E-01) > DOM, SIZE, 1.000000E-01, 5.870000E-02, 3.750000E-01 > DOM, MONIT, 5.000000E-02, 5.576500E-02, 3.656250E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, 5.870000E-03, 0.000000E+00 > OBJ, GEOMETRY, polcu5t > OBJ, ROTATION24, 1 > OBJ, TYPE, INLET > OBJ, DENSITY, 7.499400E+02 > OBJ, PRESSURE-2, 0.000000E+00 > OBJ, VELOCITY, 0.000000E+00, 0.000000E+00, 2.000000E+00 > OBJ, VELOCITY-2, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, TEMPERATURE, -6.000000E+04 > OBJ, TEMPERATURE-2, 0.000000E+00 > OBJ, TURB-INTENS, 5.000000E+00 > OBJ, NAME, HEATER > OBJ, POSITION, 0.000000E+00, 5.870000E-02, 1.500000E-01 > OBJ, SIZE, 1.000000E-01, 0.000000E+00, 2.250000E-01 > OBJ, GEOMETRY, polcu2 > OBJ, ROTATION24, 1 > OBJ, TYPE, PLATE > OBJ, SURF_ENT1, 0.000000E+00, 1.366000E+07 > OBJ, ADIABATI2, 0.000000E+00, 0.000000E+00 > OBJ, NAME, OUTLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 3.750000E-01 > OBJ, SIZE, 1.000000E-01, 5.870000E-03, 0.000000E+00 > OBJ, GEOMETRY, polcubet > OBJ, ROTATION24, 1 > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 2.000000E+01 > OBJ, TEMPERATURE-2, SAME > OBJ, COEFFICIENT-2, 1.000000E+00 > OBJ, TURBULENCE, SAME , SAME STOP