TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2018, Date 08/01/18 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY Library case 635: Tutorial: Complete Diaphragm rupture ENDDIS VRV USE * Start of frame * Setting display swiches WIREFRAME ON * Domain scaling factors SCALE 1.000000E+00 1.000000E+00 1.000000E+02 * Settings for current slice PROBE 1.975000E-01 9.750000E-02 9.900001E-04; PROBE ON SLICE Y SLICE OUTLINE ON VIEW 3.230001E-02 5.183482E-01 8.545594E-01 UP -1.225601E-02 8.551467E-01 -5.182412E-01 VIEW CENTRE 1.000000E-01 5.000000E-02 5.000000E-04 VIEW SIZE 1.175755E-01 VIEW DEPTH 3.000000E+00;VIEW TILT 0.8 NEARPLANE 4.000000E-04 VARIABLE Pressure; CON ON CONTOUR BLANK ON CONTOUR AVERAGE ON msg Pressure variations with time in Y plane. (Z axis is time) PAUSE ENDUSE ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Complete Diaphragm Rupture ) ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,1.0E-03,100) * Modify regions ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,40,20,1) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(145)=SPH1 ;NAME(146)=TEM1 NAME(148)=DEN1 ;NAME(150)=PRPS * Solved variables list SOLVE(P1,U1,V1,TEM1) * Stored variables list STORE(PRPS,DEN1,SPH1) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices TERMS (TEM1,Y,Y,Y,Y,Y,Y) NEWRH1 = T NEWENL = T ************************************************************ Group 9. Properties PRESS0 =1.01325E+05 ;TEMP0 =273. * Domain material index is 2 signifying: * Air using Ideal Gas Law, STP SETPRPS(1, 2) ENUT =0. DRH1DP = GRND5 DVO1DT =3.33E-03 PRNDTL(TEM1)=-0.0263 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(TEM1)=23. 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 No PATCHes used for this Group 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 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) ************************************************************ Group 18. Limits VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75 ************************************************************ Group 19. EARTH Calls To GROUND Station PARSOL = F CONWIZ = T ISG62 = 1 SPEDAT(SET,OUTPUT,NOFIELD,L,T) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=4.856E+04 ;EX(U1)=46.02 EX(V1)=2.367E-05 ;EX(SPH1)=1004. EX(TEM1)=23.030001 ;EX(DEN1)=1.781 EX(PRPS)=1. ************************************************************ Group 21. Print-out of Variables OUTPUT(SPH1,N,N,Y,Y,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 40 ;IYMON = 20 ;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 = 1 ;IDISPB = 0 ;IDISPC = 0 GVIEW(P,0.0323,0.518348,0.854559) GVIEW(UP,-0.012256,0.855147,-0.518241) GVIEW(VDIS,0.117576) GVIEW(CENTRE,0.1,0.05,0.05) > DOM, SIZE, 2.000000E-01, 1.000000E-01, 1.000000E-01 > DOM, MONIT, 1.975000E-01, 9.750000E-02, 5.000000E-02 > 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, 20, 1.000000E+00 > GRID, RSET_X_2, 20, 1.000000E+00 > GRID, RSET_Y_1, 20, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > DOM, T_AMBIENT, 2.300000E+01 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, HIGH-P > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, TO_END, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cubet > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, INI_PRESS, 1.0E+05 > OBJ, NAME, PROBE > OBJ, POS-CEN, 1.775000E-01, 4.750000E-02, 5.000000E-02 > OBJ, SIZE, 4.999995E-03, 5.000003E-03, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, default > OBJ, TYPE, POINT_HISTORY > OBJ, VARLIST1, P1 ,U1 ,V1 ,TEM1 > OBJ, NAME, PROBE_3 > OBJ, POS-CEN, 9.750000E-02, 4.750000E-02, 5.000000E-02 > OBJ, SIZE, 5.000003E-03, 5.000003E-03, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, default > OBJ, TYPE, POINT_HISTORY > OBJ, VARLIST1, P1 ,U1 ,V1 ,TEM1 > OBJ, NAME, PROBE_4 > OBJ, POS-CEN, 1.750000E-02, 4.750000E-02, 5.000000E-02 > OBJ, SIZE, 5.000001E-03, 5.000003E-03, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, default > OBJ, TYPE, POINT_HISTORY > OBJ, VARLIST1, P1 ,U1 ,V1 ,TEM1 STOP