TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2022, Date 21/03/22 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY Library case 653: Tutorial: Water flow into a bund (VOF) ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Water pouring into a bund ) ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,1.,1000) * Modify regions ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,28,12,17) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(144)=PRPS ;NAME(145)=VFOL NAME(146)=SURN ;NAME(147)=EPKE NAME(148)=DEN1 ;NAME(149)=EL1 NAME(150) =ENUT * Solved variables list SOLVE(P1,U1,V1,W1,VFOL,SURN) * Stored variables list STORE(ENUT,EL1,DEN1,EPKE,PRPS) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) TURMOD(KECHEN) ************************************************************ 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. RHO1 =998.22998 ENUL =1.006E-06 CP1 =4181.799805 DVO1DT =1.18E-04 PRT(KE)=0.75 ;PRT(EP)=1.15 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(PRPS)=-1. 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, 1000) COVAL(BUOYANCY, U1, FIXFLU, GRND1) COVAL(BUOYANCY, V1, FIXFLU, GRND1) COVAL(BUOYANCY, W1, FIXFLU, GRND1) BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81 EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 20 RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(SURN)=1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(KE ,LINRLX,0.5 ) RELAX(EP ,LINRLX,0.5 ) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(SURN)=1. ;VARMIN(SURN)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station GENK = T SURF = T PARSOL = F CONWIZ = T IENUTA = 2 IPRPSA = 67 ISG62 = 1 ISURFA = 1 SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout ************************************************************ 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 = 22 ;IYMON = 6 ;IZMON = 2 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 = 10 ;IDISPB = 0 ;IDISPC = 0 CSG1 ='N' GVIEW(P,0.176433,0.975371,0.132371) GVIEW(UP,-0.047601,-0.125869,0.990905) GVIEW(VDIS,0.804453) GVIEW(CENTRE,0.5,0.5,0.15) > DOM, SIZE, 1.000000E+00, 1.000000E+00, 3.000000E-01 > DOM, MONIT, 7.555072E-01, 4.642857E-01, 3.190332E-02 > 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, -14, 1.100000E+00,G > GRID, RSET_X_2, -3, 1.100000E+00,G > GRID, RSET_X_3, -11, 1.100000E+00,G > GRID, RSET_Y_1, 2,-1.100000E+00,G > GRID, RSET_Y_2, -4, 1.100000E+00,G > GRID, RSET_Y_3, -4, 1.100000E+00,G > GRID, RSET_Y_4, 2, 1.100000E+00,G > GRID, RSET_Z_1, 3,-1.100000E+00,G > GRID, RSET_Z_2, -2, 1.100000E+00,G > GRID, RSET_Z_3, -12, 1.100000E+00,G > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, BEAM > OBJ, POSITION, 5.000000E-01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, TO_END, 6.000000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 198,Solid with smooth-wall friction > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, WATER-IN > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 5.000000E-01, 9.000000E-02 > OBJ, GEOMETRY, box > OBJ, TYPE, INLET > OBJ, FLUID, HEAVY > OBJ, VELOCITY, 5. ,0. ,0. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURB-INTENS, 5. > OBJ, NAME, OUTLET > OBJ, POSITION, AT_END, 2.000000E-01, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 6.000000E-01, 9.000000E-02 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, OUTLET > OBJ, PRESSURE, P_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, FREE-TOP > OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, OUTLET > OBJ, PRESSURE, P_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, SAME , SAME STOP