TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2011, Date 23/03/11 CPVNAM=VDI;SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings DISPLAY This case shows the use of InForm to set the inflow boundary for a Scalar-Equation-Model (SEM) case ENDDIS VRV USE SLICE CLEAR PROBE 4.500000E-01 4.750000E-01 5.000000E-01 SLICE Z; SLICE SAVE PROBE 4.500000E-01 4.750000E-01 1.500000E+00 SLICE Z; SLICE SAVE * Settings for current slice PROBE 4.500000E-01 4.750000E-01 2.500000E+00 SLICE Z; SLICE OUTLINE ON * View and up directions VIEW 1.487207E-01 5.626179E-01 8.132296E-01 UP -4.604268E-02 8.254228E-01 -5.626335E-01 * View centre VIEW CENTRE 1.000000E+00 5.000000E-01 1.500000E+00 * View size VIEW SIZE 1.905372E+00 * View perspective VIEW DEPTH 3.000000E+00;VIEW TILT 0.8 NEARPLANE 3.000000E-02 * Setting object visibility and painting status OBJECT WIREFRAME TYPE BLOCKAGE ON OBJECT WIREFRAME TYPE INLET ON OBJECT WIREFRAME TYPE OUTLET ON OBJECT WIREFRAME TYPE PLATE ON VARIABLE VFOL; VARIABLE RANGE 0.000000E+00 1.000000E+00 VECTOR ON; CONTOUR ON MSG Press OK to animate the solution PAUSE ANIMATE ENDUSE ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(SEM Inflow boundary set with InForm ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,1.,10) * Modify regions ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,20,20,3) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(146)=HPOR ;NAME(147)=PRPS NAME(148)=DEN1 ;NAME(149)=VFOL NAME(150) =SURN * Solved variables list SOLVE(P1,U1,V1,W1,VFOL,SURN) * Stored variables list STORE(DEN1,PRPS,HPOR) * 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.6 ;RLOLIM =0.4 ************************************************************ Group 9. Properties PRESS0 =1.01325E+05 ;TEMP0 =273. RHO1 =998.22998 ENUL =1.006E-06 CP1 =4181.799805 ENUT =3.088E-03 DVO1DT =1.18E-04 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(HPOR)=1. 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 ************************************************************ Echo save-block settings for Group 13 save13begin (source of p1 at in2 is 1*998.3 with area) (source of u1 at in2 is 1 with onlyms) (source of vfol at in2 is 1/998.3 with onlyms) (source of surn at in2 is 1 with area) (source of p1 at in3 is coval(fixflu,1*998.3) with area) (source of u1 at in3 is coval(onlyms,1)) (source of vfol at in3 is coval(onlyms,1/998.3)) (source of surn at in3 is coval(fixflu,1) with area) save13end ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 50 RESFAC =1.0E-03 ************************************************************ 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 CONWIZ = T IPRPSA = 67 ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=140.399994 ;EX(U1)=1. EX(V1)=0.1462 ;EX(W1)=0. EX(HPOR)=0.3333 ;EX(PRPS)=55.279999 EX(DEN1)=823.299988 ;EX(VFOL)=0.8244 EX(SURN)=0.853 ************************************************************ 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 = 5 ;IYMON = 10 ;IZMON = 1 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NTPRIN = 1 ;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 CSG1 ='A' GVIEW(P,0.160511,0.524429,0.836188) GVIEW(UP,-0.039454,0.849904,-0.525458) > DOM, SIZE, 2.000000E+00, 1.000000E+00, 3.000000E+00 > DOM, MONIT, 4.500000E-01, 4.750000E-01, 5.000000E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, AUTO, T T F > GRID, RSET_X_1, 20, 1.000000E+00 > GRID, RSET_Y_1, 10, 1.000000E+00 > GRID, RSET_Y_2, 10, 1.000000E+00 > GRID, RSET_Z_1, 1, 1.000000E+00 > GRID, RSET_Z_2, 1, 1.000000E+00 > GRID, RSET_Z_3, 1, 1.000000E+00 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, IN1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube3t > OBJ, TYPE, INLET > OBJ, FLUID, 1. > OBJ, VELOCITY, 1. ,0. ,0. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, OUT1 > OBJ, POSITION, AT_END, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube12t > OBJ, TYPE, OUTLET > OBJ, PRESSURE, P_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, B3 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 5.000000E-01, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cubet > OBJ, TYPE, BLOCKAGE > OBJ, WIREFRAME, YES > OBJ, MATERIAL, 67, WATER at 20. deg C > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, IN2 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube3t > OBJ, TYPE, INLET > OBJ, FLUID, 1. > OBJ, VELOCITY, 1. ,0. ,0. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, IN3 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube3t > OBJ, TYPE, INLET > OBJ, FLUID, 1. > OBJ, VELOCITY, 1. ,0. ,0. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, OUT2 > OBJ, POSITION, AT_END, 0.000000E+00, 1.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube12t > OBJ, TYPE, OUTLET > OBJ, PRESSURE, P_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, OUT3 > OBJ, POSITION, AT_END, 0.000000E+00, 2.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, 1.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube12t > OBJ, TYPE, OUTLET > OBJ, PRESSURE, P_AMBIENT > OBJ, COEFFICIENT, 1000. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, B8 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 1.000000E+00 > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube11 > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, POROSITY, 0. > OBJ, SIDE, BOTH > OBJ, WALL_CO_L, FRICTIONLESS > OBJ, WALL_CO_H, FRICTIONLESS > OBJ, NAME, B9 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 2.000000E+00 > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, cube11 > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, POROSITY, 0. > OBJ, SIDE, BOTH > OBJ, WALL_CO_L, FRICTIONLESS > OBJ, WALL_CO_H, FRICTIONLESS STOP