TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2020, Date 16/06/20 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY Library case 651: Tutorial: Boiling in a pipe (IPSA) ENDDIS ************************************************************ 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.0E+05 ;TEMP0 =273. RHO1 =749.940002 ENUL =1.544E-05 CP1 =1005. * Phase 2 material index is 67 signifying: * WATER at 20. deg C SETPRPS(2, 67) RHO2 =35.549999 DVO1DT =3.41E-03 DVO2DT =1.18E-04 PHINT(H1)=0.1 ;PHINT(H2)=1.513E+06 PRNDTL(H1)=-0.0258 ;PRNDTL(H2)=-0.597 PRT(KE)=0.75 ;PRT(EP)=1.15 ************************************************************ Group 10.Inter-Phase Transfer Processes CFIPS =2500. RLOLIM =0. ;CMDOT = HEATBL CINT(H1)=10. ;CINT(H2)=10. ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(W1)=2. ;FIINIT(W2)=2. FIINIT(R1)=0.999 ;FIINIT(R2)=1.0E-04 FIINIT(H2)=1.513E+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, V2, FIXFLU, GRND1) COVAL(BUOYANCY, W1, FIXFLU, GRND1) COVAL(BUOYANCY, W2, FIXFLU, GRND1) BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81 BUOYE =20. EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 750 RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(R1)=1 LITER(R2)=1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(V1 ,FALSDT,0.1 ) RELAX(V2 ,FALSDT,0.1 ) RELAX(W1 ,FALSDT,0.1 ) RELAX(W2 ,FALSDT,0.1 ) RELAX(R1 ,LINRLX,0.25 ) RELAX(R2 ,LINRLX,0.25 ) RELAX(KE ,LINRLX,0.5 ) RELAX(EP ,LINRLX,0.5 ) RELAX(H1 ,FALSDT,10. ) RELAX(H2 ,FALSDT,10. ) RELAX(CMDO,LINRLX,0.1 ) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(R1)=1. ;VARMIN(R1)=1.0E-06 VARMAX(R2)=1. ;VARMIN(R2)=1.0E-06 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-7.0E+04 VARMAX(H2)=1.0E+10 ;VARMIN(H2)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station GENK = T PARSOL = F IENUTA = 2 ISG62 = 1 SPEDAT(SET,FACETDAT,NUMOBJ,I,4) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=1.256E+03 ;EX(V1)=0.02349 EX(V2)=0.02354 ;EX(W1)=0.1554 EX(W2)=0.1548 ;EX(R1)=0.8894 EX(R2)=0.1106 ;EX(KE)=0.01628 EX(EP)=0.27 ;EX(H1)=3532. EX(H2)=1.513E+06 ;EX(CMDO)=6.808E-07 EX(EPKE)=6.67 ;EX(DEN1)=749.900024 EX(EL1)=3.825E-03 ;EX(ENUT)=1.981E-04 ************************************************************ 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,-0.986196,-0.010995,0.16522) GVIEW(UP,0.165253,-2.182596E-03,0.986249) GVIEW(VDIS,0.237748) GVIEW(CENTRE,2.930111E-03,0.02935,0.1875) > 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 > GRID, RSET_X_1, 1, 1.000000E+00 > GRID, RSET_Y_1, 1, 1.000000E+00 > GRID, RSET_Y_2, 9, 1.000000E+00 > GRID, RSET_Z_1, 8, 1.000000E+00 > GRID, RSET_Z_2, 12, 1.000000E+00 > OBJ, NAME, INLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 5.870000E-03, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, INLET > OBJ, DENSITY, 749.940002 > OBJ, PRESSURE-2, 0. > OBJ, VELOCITY, 0. ,0. ,2. > OBJ, VELOCITY-2, 0. ,0. ,0. > OBJ, TEMPERATURE, -6.0E+04 > OBJ, TEMPERATURE-2,0. > OBJ, TURB-INTENS, 5. > OBJ, NAME, HEATER > OBJ, POSITION, 0.000000E+00, AT_END, 1.500000E-01 > OBJ, SIZE, TO_END, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, PLATE > OBJ, SURF_ENT1, 0. ,1.366E+07 > OBJ, NAME, OUTLET > OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END > OBJ, SIZE, TO_END, 5.870000E-03, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, poldef > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, SAME > OBJ, COEFFICIENT, 20. > OBJ, TEMPERATURE-2, SAME > OBJ, COEFFICIENT-2,1. > OBJ, TURBULENCE, SAME , SAME STOP