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 p126 Zalesak Disk ENDDIS VRV USE SLICE Z PROBE STYLE SPHERE PROBE SIZE 1.000000E+00 PROBE COLOUR 0 0 255 VARIABLE SURN; VARIABLE RANGE -1.000000E-05 1.000100E+00 CONTOUR ON PAUSE ANIMATE ENDUSE ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Test Zalesak disk CICSAM ) ************************************************************ 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,100,100,1) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(147)=PRPS ;NAME(148)=DEN1 NAME(149)=VFOL ;NAME(150)=SURN * Solved variables list SOLVE(P1,U1,V1,VFOL,SURN) * Stored variables list STORE(DEN1,PRPS) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) ************************************************************ Echo save-block settings for Group 7 save7begin save7end ************************************************************ 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 ENUT =0. DVO1DT =1.18E-04 ************************************************************ 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 ************************************************************ Echo save-block settings for Group 11 save11begin PATCH(IINIT,INIVAL,1,NX,1,NY,1,NZ,1,1) (INITIAL of U1 at IINIT is -6.28*(YG-0.5)) (INITIAL of V1 at IINIT is 6.28*(XG-0.5)) save11end ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources No PATCHes used for this Group EGWF = T ************************************************************ Echo save-block settings for Group 13 save13begin PATCH(ROTVUS, VOLUME, 1, NX, 1, NY, 1, NZ, 1, LSTEP) (source of P1 at ROTVUS is 0.0 with FIXVAL) (source of U1 at ROTVUS is -6.28*(YG-0.5) with FIXVAL) (source of V1 at ROTVUS is 6.28*(XG-0.5) with FIXVAL) save13end ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 40 RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(SURN)=1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) 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 CONWIZ = T IPRPSA = 67 ISG62 = 1 ISURFA = 1 SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,GXMONI,CLASSIC,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,N,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 51 ;IYMON = 50 ;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 = 10 ;IDISPB = 0 ;IDISPC = 0 CSG1 ='A' GVIEW(P,0.,0.,1.) GVIEW(UP,0.,1.,0.) GVIEW(DEPTH,3000.) GVIEW(VDIS,0.831384) GVIEW(CENTRE,0.5,0.5,0.5) > DOM, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, MONIT, 5.000000E-01, 4.950000E-01, 5.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, RSET_X_1, 38, 1.000000E+00 > GRID, RSET_X_2, 10, 1.000000E+00 > GRID, RSET_X_3, 5, 1.000000E+00 > GRID, RSET_X_4, 10, 1.000000E+00 > GRID, RSET_X_5, 37, 1.000000E+00 > GRID, RSET_Y_1, 50, 1.000000E+00 > GRID, RSET_Y_2, 20, 1.000000E+00 > GRID, RSET_Y_3, 5, 1.000000E+00 > GRID, RSET_Y_4, 25, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, OUTL1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 0.000000E+00, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube12t > OBJ, VISIBLE, NO > OBJ, TYPE, OUTLET > OBJ, COLOR-MODE, FROM_FIL > OBJ, PRESSURE, P_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, EXTERNAL, HEAVY > OBJ, NAME, BLOC2 > OBJ, POSITION, 3.750000E-01, 5.000000E-01, 0.000000E+00 > OBJ, SIZE, 2.500000E-01, 2.500000E-01, 1.000000E+00 > OBJ, GEOMETRY, cylinder > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 212 > OBJ, OPAQUE, 20 > OBJ, MATERIAL, 67, WATER at 20. deg C (HEAVY) > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, BLOC3 > OBJ, POSITION, 4.750000E-01, 5.000000E-01, 0.000000E+00 > OBJ, SIZE, 5.000000E-02, 2.000000E-01, 1.000000E+00 > OBJ, GEOMETRY, cubet > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 2 > OBJ, OPAQUE, 98 > OBJ, MATERIAL, 0, Air at 20 deg C, 1 atm, treated as incompressible (LIGHT) > OBJ, TIME_LIMITS, ALWAYS_ACTIVE STOP