TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2017, Date 05/10/17 CPVNAM=VDI; SPPNAM=Core ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT( ) VRV USE * Start of frame SLICE Z * View and up directions VIEW 0 0 -1 UP 0 1 0 * Setting object visibility and painting status OBJECT HIDE TYPE BLOCKAGE VECTOR OFF VARIABLE VFOL; VARIABLE RANGE 0.000000E+00 1.000000E+00 CONTOUR ON CONTOUR FILL ON PAUSE ANIMATE START 10 END 1000 INTERVAL 10 ENDUSE ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,1.,1000) * Modify regions RSET(T,1,1000,1.) ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,64,256,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) ************************************************************ 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 =1.353E+04 ENUL =1.125E-07 CP1 =139.300003 ENUT =0. DVO1DT =1.8E-04 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(PRPS)=-1. FIINIT(VFOL)=1. ;FIINIT(SURN)=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, GRND2) COVAL(BUOYANCY, V1, FIXFLU, GRND2) BUOYA =0. ; BUOYB =-9.81 BUOYC =0. BUOYD =861.150024 EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 50 RESFAC =1.0E-04 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(SURN)=1 ENDIT(P1)=GRND1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(VFOL,FALSDT,10. ) RELAX(SURN,LINRLX,0.5 ) ************************************************************ Group 18. Limits VARMAX(SURN)=1. ;VARMIN(SURN)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station SURF = T PARSOL = F IPRPSA = 66 ;IPRPSB = 59 ISURFA = 1 ISG21 = 5 ISG62 = 1 SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,GXMONI,TRANSIENT,L,F) ************************************************************ 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 = 32 ;IYMON = 148 ;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(NEARPLANE,2.0E-03) GVIEW(VDIS,2.736842) GVIEW(CENTRE,0.5,2.,0.5) > DOM, SIZE, 1.000000E+00, 4.000000E+00, 1.000000E+00 > DOM, MONIT, 4.921875E-01, 2.304688E+00, 5.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, RSET_X_1, 64, 1.000000E+00 > GRID, RSET_Y_1, 96, 1.000000E+00 > GRID, RSET_Y_2, 64, 1.000000E+00 > GRID, RSET_Y_3, 96, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, BLOC6 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 1.500000E+00, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cubet > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 59, Paraffin oil at 25 deg C, 1 atm > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, BLOC1 > OBJ, POSITION, 0.000000E+00, 1.500000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > OBJ, GEOMETRY, sinewave > OBJ, VISIBLE, NO > OBJ, TYPE, BLOCKAGE > OBJ, COLOR-MODE, DEFAULT2 > OBJ, COLOR-VAL, 27 > OBJ, OPAQUE, 25 > OBJ, MATERIAL, 59, Paraffin oil at 25 deg C, 1 atm > OBJ, TIME_LIMITS, ALWAYS_ACTIVE STOP