TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2025, Date 04/11/24 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Film Boiling in 2D THINC ) ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,2.,4000) * Modify regions RSET(T,1,4000,1.) ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,64,289,1,1.0E-05,1.0E-04,1.0E-03) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names * Solved variables list SOLVE(P1,U1,V1,SURN,VFOL,TEM1) * Stored variables list STORE(SPH1,DEN1,PRPS,MEVA,KAPA,MEVP,SIGM) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices TERMS (SURN,N,N,N,N,Y,N) TERMS (VFOL,N,N,N,N,Y,N) GALA = T RUPLIM =0.999 ;RLOLIM =1.0E-03 ************************************************************ Group 9. Properties PRESS0 =1.01325E+05 RHO1 =200. ENUL =5.0E-04 CP1 =400. ENUT =0. DVO1DT =1.18E-04 PRNDTL(TEM1)=-40. * List of user-defined materials to be read by EARTH MATFLG=T;IMAT=2 * Name *Ind. Dens. Viscos. Spec.heat Conduct. Expans. Compr. *34 5.0 1.0E-3 200.0 1.0 3.33E-3 * 78 200.0 5.0E-4 400.0 40.0 1.18E-4 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(SIGM)=0.1 FIINIT(PRPS)=-1. ;FIINIT(SURN)=1. FIINIT(VFOL)=1. ;FIINIT(TEM1)=0. No PATCHes used for this Group INIADD = F NAMFI ='phida' ************************************************************ Echo save-block settings for Group 11 save11begin REAL(LAMD0, PI, A0, A1) CHAR(FORM1) LAMD0=0.078684409 PI=3.1415927 A0=0.0029191916 !0.0371*LAMD0 A1=0.00059800151 !0.0076*LAMD0 !FORM1=A0+A1*cos(2.0*PI*XG/0.078684409) FORM1=A0+A1*cos(79.85299*XG) (INIT of SURN is 0.0 with if(YG.LE.(:FORM1:))) (INIT of VFOL is 0.0 with if(YG.LE.(:FORM1:))) (INIT of TEM1 is 5.0-5.0*YG/:FORM1: with if(YG.LE.(:FORM1:))) TEMP0=0.0 RSG163=0.01 save11end ************************************************************ 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, 4000) COVAL(BUOYANCY, U1, FIXFLU, GRND2) COVAL(BUOYANCY, V1, FIXFLU, GRND2) PATCH(STENS, PHASEM, 0, 0, 0, 0, 0, 0, 1, 4000) COVAL(STENS, U1, FIXFLU, GRND3) COVAL(STENS, V1, FIXFLU, GRND3) BUOYA =0. ; BUOYB =-9.81 ; BUOYC =0. BUOYD =5. EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 50 RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(SURN)=1 LITER(VFOL)=1 ;LITER(TEM1)=200 ENDIT(P1)=1.0E-06 ;ENDIT(TEM1)=1.0E-06 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,0.5 ) RELAX(U1 ,FALSDT,0.05 ) RELAX(V1 ,FALSDT,0.05 ) RELAX(SURN,LINRLX,1. ) RELAX(VFOL,LINRLX,1. ) RELAX(TEM1,LINRLX,0.1 ) ************************************************************ Group 18. Limits VARMAX(SURN)=1. ;VARMIN(SURN)=0. VARMAX(TEM1)=5. ;VARMIN(TEM1)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station SURF = T PARSOL = F IPRPSA = 78 ;IPRPSB = 34 ISG21 = 5 ISG62 = 1 ISURFA = 5 BTATHC =3.5 ISRFTA = 1 SURFTA =0.1 ;SURFTB=0. ;SURFTC =0. ;SURFTD=0. IEVAP = 2 LATH =1.0E+04 ;TSAT=0. ;EVAPCO=0.1 NONCON = T RSG163 =1000. CSG10 ='q1' SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,SURFTENS,NSMOOTH,I,3) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,VOF,LTSS,L,F) SPEDAT(SET,MAXINC,SURN,R,1.) SPEDAT(SET,MAXINC,VFOL,R,1.) ************************************************************ Group 20. Preliminary Printout ************************************************************ Group 21. Print-out of Variables OUTPUT(SURN,Y,Y,Y,Y,Y,Y) OUTPUT(VFOL,Y,Y,Y,N,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 38 ;IZMON = 1 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts IDISPA = 25 ;IDISPB = 0 ;IDISPC = 0 CSG1 ='A' GVIEW(P,0.,0.,1.) GVIEW(UP,0.,1.,0.) GVIEW(DEPTH,3000.) GVIEW(VDIS,0.117168) GVIEW(CENTRE,-1.714228E-03,0.08579,5.808009E-03) > DOM, SIZE, 3.934200E-02, 1.770390E-01, 1.000000E-02 > DOM, MONIT, 3.073594E-04, 1.311604E-02, 5.000000E-03 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02 > GRID, AUTO, F F T > GRID, RSET_X_1, 64, 1.000000E+00 > GRID, RSET_Y_1, 24, 1.000000E+00 > GRID, RSET_Y_2, 265, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > DOM, T_AMBIENT, 0.000000E+00 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, S-WALL > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, SURF_TEMP, 0. ,5. > OBJ, NAME, N-OUT > OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00 > OBJ, SIZE, TO_END, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, OUTLET > OBJ, PRESSURE, P_AMBIENT > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, EXTERNAL, HEAVY > OBJ, NAME, B5 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 4.317000E-03, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, wirexyz > OBJ, VISIBLE, NO > OBJ, TYPE, NULL ************************************************************ STOP