TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2022, Date 05/07/22 CPVNAM=VDI; SPPNAM=Core Echo DISPLAY / USE settings DISPLAY Library case p135 Water and Kerosene in a T Junction ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(VOF Cartes slug Water-kero XY CA 1 ) ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,1.2,8000) * Modify regions RSET(T,1,8000,1.) ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,200,120,1,5.0E-05) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(141)=SURT ;NAME(145)=VPOR NAME(146)=ENUL ;NAME(147)=PRPS NAME(148)=VFOL ;NAME(149)=DEN1 NAME(150) =SURN * Solved variables list SOLVE(P1,U1,V1,VFOL,SURN) * Stored variables list STORE(DEN1,PRPS,ENUL,VPOR,SURT) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(VFOL,Y,Y,N,Y,N,Y) SOLUTN(SURN,Y,Y,N,Y,N,Y) ************************************************************ Group 8. Terms & Devices TERMS (VFOL,N,N,N,N,Y,N) TERMS (SURN,N,N,N,N,Y,N) GALA = T RUPLIM =0.999 ;RLOLIM =1.0E-03 ************************************************************ Group 9. Properties PRESS0 =1.01325E+05 ;TEMP0 =273.149994 RHO1 =998.22998 ENUL =1.006E-06 CP1 =4181.799805 ENUT =0. DVO1DT =1.18E-04 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(VPOR)=1. ;FIINIT(PRPS)=-1. FIINIT(VFOL)=0. ;FIINIT(DEN1)=780. FIINIT(SURN)=0. This patch is attached to object WATLEV PATCH(WATER, INIVAL, -1, 0, 0, 0, 0, 0, 1, 1) INIT(WATER , SURT, 0., 1.) INIT(WATER , VFOL, 0., 1.) INIT(WATER , SURN, 0., 1.) INIADD = F NAMFI ='phida' ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources PATCH(STENS, PHASEM, 0, 0, 0, 0, 0, 0, 1, 8000) COVAL(STENS, U1, FIXFLU, GRND3) COVAL(STENS, V1, FIXFLU, GRND3) EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 40 RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(VFOL)=1 LITER(SURN)=1 ENDIT(P1)=1.0E-05 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(U1 ,FALSDT,0.5 ) RELAX(V1 ,FALSDT,0.5 ) RELAX(VFOL,LINRLX,1. ) RELAX(SURN,LINRLX,1. ) ************************************************************ Group 18. Limits VARMAX(SURN)=1. ;VARMIN(SURN)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station SURF = T PARSOL = F IPRPSA = 67 ;IPRPSB = 58 ISG21 = 5 ISG62 = 1 ISURFA = 5 BTATHC =3.5 ISRFTA = 1 SURFTA =0.045 ;SURFTB=0. ;SURFTC=0. VOFCNT = T CNTANG =1. SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,RESTRT,PBCL,C,PBCL.DAT) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,GXMONI,CLASSIC,L,T) SPEDAT(SET,SOLPHI,P1,I,24) SPEDAT(SET,SOLPHI,U1,I,0) SPEDAT(SET,PREPHI,U1,I,0) SPEDAT(SET,SOLPHI,V1,I,0) SPEDAT(SET,PREPHI,V1,I,0) SPEDAT(SET,SOLPHI,VFOL,I,0) SPEDAT(SET,PREPHI,VFOL,I,0) SPEDAT(SET,SOLPHI,SURN,I,0) SPEDAT(SET,PREPHI,SURN,I,0) SPEDAT(SET,MAXINC,VFOL,R,1.) SPEDAT(SET,MAXINC,SURN,R,1.) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=16.049999 ;EX(U1)=2.757E-03 EX(V1)=1.602E-03 ;EX(SURT)=0.1124 EX(VPOR)=0.2275 ;EX(ENUL)=3.431E-07 EX(PRPS)=194.600006 ;EX(VFOL)=0.1121 EX(DEN1)=202.800003 ;EX(SURN)=0.1121 ************************************************************ Group 21. Print-out of Variables OUTPUT(VFOL,Y,Y,Y,N,Y,Y) OUTPUT(SURN,Y,Y,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 9 ;IYMON = 62 ;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 = 100 ;IDISPB = 0 ;IDISPC = 0 CSG1 ='A' GVIEW(P,0.,0.,1.) GVIEW(UP,0.,1.,0.) GVIEW(DEPTH,21.082172) GVIEW(VDIS,3.728045E-03) GVIEW(CENTRE,4.027397E-03,2.051923E-03,8.398048E-04) > DOM, SIZE, 1.062000E-02, 4.130000E-03, 5.000000E-04 > DOM, MONIT, 2.507501E-04, 2.134412E-03, 2.500000E-04 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02 > GRID, RSET_X_1, 20, 1.000000E+00 > GRID, RSET_X_2, 20, 1.000000E+00 > GRID, RSET_X_3, 160, 1.000000E+00 > GRID, RSET_Y_1, 12, 1.000000E+00 > GRID, RSET_Y_2, 39, 1.000000E+00 > GRID, RSET_Y_3, 17, 1.000000E+00 > GRID, RSET_Y_4, 42, 1.000000E+00 > GRID, RSET_Y_5, 10, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > DOM, T_AMBIENT, 0.000000E+00 > OBJ, NAME, WATLEV > OBJ, POSITION, 0.000000E+00, 3.800000E-03, 0.000000E+00 > OBJ, SIZE, 5.900000E-04, TO_END, 5.000000E-04 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, USER_DEFINED > OBJ, PATCHES, WATER > OBJ, NAME, BLK1 > OBJ, POSITION, 5.900000E-04, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.003000E-02, 1.770000E-03, 5.000000E-04 > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, BLK2 > OBJ, POSITION, 5.900000E-04, AT_END, 0.000000E+00 > OBJ, SIZE, 1.003000E-02, 1.770000E-03, 5.000000E-04 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, INW > OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00 > OBJ, SIZE, 5.900000E-04, 0.000000E+00, 5.000000E-04 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, INLET > OBJ, FLUID, HEAVY > OBJ, VELOCITY, 0. ,-9.4136E-03 ,0. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, INK > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 5.900000E-04, 0.000000E+00, 5.000000E-04 > OBJ, GEOMETRY, box > OBJ, TYPE, INLET > OBJ, VELOCITY, 0. ,9.4136E-03 ,0. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, OUT > OBJ, POSITION, AT_END, 1.770000E-03, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 5.900000E-04, 5.000000E-04 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, OUTLET > OBJ, PRESSURE, 0. > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, EXTERNAL, HEAVY > OBJ, NAME, PLAT7 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, 5.000000E-04 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, NULL8 > OBJ, POSITION, 5.900000E-04, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 5.900000E-04, 4.130000E-04, 5.000000E-04 > OBJ, GEOMETRY, wirexyz > OBJ, TYPE, NULL STOP