TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2015, Date 20/04/16 CPVNAM=VDI; SPPNAM=Core ************************************************************ Echo DISPLAY / USE settings vrv use PROBE 9.712500E-02 8.750001E-02 1.602857E-01; PROBE ON SLICE X SLICE OUTLINE ON * View and up directions VIEW 9.456690E-02 7.978453E-01 -5.953989E-01 UP 9.955179E-01 -7.642952E-02 5.570072E-02 * View centre VIEW CENTRE 2.500000E-01 2.500000E-01 2.500000E-01 * View size VIEW SIZE 4.156922E-01 * View perspective VIEW DEPTH 3.000000E+00;VIEW TILT 0.8 NEARPLANE 1.000000E-03 VARIABLE C1 SURFACE ON; SURFACE VALUE 5.000000E-01 msg Press OK to animate the results PAUSE ANIMATE enduse ************************************************************ IRUNN = 1 ;LIBREF = 111 ************************************************************ Group 1. Run Title TEXT(chemical reaction between stirred fluids) ************************************************************ Echo save-block settings for Group 1 save1begin BOOLEAN(MOVOB) MOVOB=T save1end ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,1.,20) * Modify regions ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,20,20,20) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(144)=DEN1 ;NAME(145)=OBID NAME(146)=FUEL ;NAME(147)=TMP1 NAME(148)=PROD ;NAME(149)=OXID NAME(150) =MIXF * Solved variables list SOLVE(P1,U1,V1,W1,H1,C1,MIXF) * Stored variables list STORE(OXID,PROD,TMP1,FUEL,OBID,DEN1) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(H1,Y,Y,Y,N,N,Y) ************************************************************ Group 8. Terms & Devices TERMS (H1 ,N,Y,Y,Y,Y,Y) DIFCUT =0. HUNIT =1. ************************************************************ Group 9. Properties PRESS0 =1.0E+05 ;TEMP0 =273. * Domain material index is 0 signifying: * Air at 20 deg C, 1 atm, treated as incompressible SETPRPS(1, 0) RHO1 =1.189 RHO1A =29. ;RHO1B =29. RHO1C =28. CP1 = GRND10 CP1A =1005. ;CP1B =1005. CP1C =1005. ;CP1D =0. TMP1 = GRND7 TMP2A =0.5 ;TMP2B =0. TMP2C =0. ENUT =3.088E-03 DRH1DP =1.0E-05 DVO1DT =3.41E-03 PRNDTL(H1)=-0.0258 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(P1)=0. ;FIINIT(U1)=0. FIINIT(V1)=0. ;FIINIT(W1)=0. FIINIT(H1)=3.015E+05 ;FIINIT(C1)=0. FIINIT(DEN1)=1.2 ;FIINIT(TMP1)=27. FIINIT(MIXF)=0. This patch is attached to object FLU1IN PATCH(INIC1, INIVAL, -1, 0, 0, 0, 0, 0, 1, 20) INIT(INIC1 , C1, 0., 1.) INIT(INIC1 , MIXF, 0., 1.) INIADD = F ************************************************************ Group 12. Convection and diffusion adjustments No PATCHes used for this Group ************************************************************ Group 13. Boundary & Special Sources PATCH(CHSO, VOLUME, 0, 0, 0, 0, 0, 0, 1, 20) EGWF = T ************************************************************ Echo save-block settings for Group 13 SAVE13BEGIN IF(MOVOB) THEN (MOVOB of CHAM is OFFSET(0&0&0)) (MOVOB of SHAFT is OFFSET(0.15&0.25&0.25) with PARENT=CHAM) (MOVOB of SHAFT is POS(0&0&0&360*tim&0&0)) (MOVOB of ROD is OFFSET(0&0&0) with PARENT=SHAFT) (MOVOB of PADDLE1 is OFFSET(0&0&0) with PARENT=ROD) (MOVOB of PADDLE2 is OFFSET(0&0&0) with PARENT=PADDLE1) ENDIF SAVE13END ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 20 RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(H1 ,FALSDT,1. ) RELAX(C1 ,FALSDT,1. ) RELAX(MIXF,FALSDT,1. ) ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-5.0E+04 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=0. VARMAX(C1)=1. ;VARMIN(C1)=0. VARMAX(MIXF)=1. ;VARMIN(MIXF)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station ISG62 = 0 SPEDAT(SET,GXMONI,TRANSIENT,L,F) SPEDAT(SET,MOFOR,MOFFILE,C,NOTSET) ************************************************************ Echo save-block settings for Group 19 save19begin IF(MOVOB) THEN SPEDAT(SET,MOFOR,mofFILE,C,NOTSET) ELSE SPEDAT(SET,MOFOR,mofFILE,C,AGITATED.mof) ENDIF save19end ************************************************************ Group 20. Preliminary Printout ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 4 ;IYMON = 4 ;IZMON = 6 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 10000 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0 CSG1 ='M' GVIEW(P,0.318169,0.776036,-0.544552) GVIEW(UP,0.9472,-0.236126,0.216926) GVIEW(VDIS,0.415692) GVIEW(CENTRE,0.25,0.25,0.25) > DOM, SIZE, 5.000000E-01, 5.000000E-01, 5.000000E-01 > DOM, MONIT, 9.712500E-02, 8.750001E-02, 1.602857E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, RSET_X_1, 4, 1.000000E+00 > GRID, RSET_X_2, 1, 1.000000E+00 > GRID, RSET_X_3, 1, 1.000000E+00 > GRID, RSET_X_4, 1, 1.000000E+00 > GRID, RSET_X_5, 1, 1.000000E+00 > GRID, RSET_X_6, 1, 1.000000E+00 > GRID, RSET_X_7, 1, 1.000000E+00 > GRID, RSET_X_8, 10, 1.000000E+00 > GRID, RSET_Y_1, 1, 1.000000E+00 > GRID, RSET_Y_2, 4, 1.000000E+00 > GRID, RSET_Y_3, 4, 1.000000E+00 > GRID, RSET_Y_4, 2, 1.000000E+00 > GRID, RSET_Y_5, 4, 1.000000E+00 > GRID, RSET_Y_6, 4, 1.000000E+00 > GRID, RSET_Y_7, 1, 1.000000E+00 > GRID, RSET_Z_1, 7, 1.000000E+00 > GRID, RSET_Z_2, 1, 1.000000E+00 > GRID, RSET_Z_3, 1, 1.000000E+00 > GRID, RSET_Z_4, 2, 1.000000E+00 > GRID, RSET_Z_5, 1, 1.000000E+00 > GRID, RSET_Z_6, 1, 1.000000E+00 > GRID, RSET_Z_7, 7, 1.000000E+00 > OBJ, NAME, SHAFT > OBJ, POSITION, 1.500000E-01, 2.250000E-01, 2.250000E-01 > OBJ, SIZE, 3.500000E-01, 5.000000E-02, 5.000000E-02 > OBJ, GEOMETRY, cylinder > OBJ, ROTATION24, 9 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, ROD > OBJ, POSITION, 1.250000E-01, 1.250000E-01, 2.250000E-01 > OBJ, SIZE, 5.000000E-02, 2.500000E-01, 5.000000E-02 > OBJ, GEOMETRY, cylinder > OBJ, ROTATION24, 5 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, PADDLE1 > OBJ, POSITION, 1.110000E-01, 3.750000E-01, 2.870000E-01 > OBJ, SIZE, 1.000000E-01, 1.000000E-01, 2.500000E-02 > OBJ, ROT-ANGLE, 0.000000E+00, 4.500000E+01, 0.000000E+00 > OBJ, GEOMETRY, cubet > OBJ, TYPE, BLOCKAGE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 7 > OBJ, MATERIAL, DOMAIN > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, PADDLE2 > OBJ, POSITION, 1.320000E-01, 2.500000E-02, 2.040000E-01 > OBJ, SIZE, 1.000000E-01, 1.000000E-01, 2.500000E-02 > OBJ, ROT-ANGLE, 0.000000E+00,-4.500000E+01, 0.000000E+00 > OBJ, GEOMETRY, cubet > OBJ, TYPE, BLOCKAGE > OBJ, COLOR-MODE, USER > OBJ, COLOR-VAL, 7 > OBJ, MATERIAL, DOMAIN > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, FLU1IN > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 2.000000E-01, 5.000000E-01, 5.000000E-01 > OBJ, GEOMETRY, cubet > OBJ, TYPE, USER_DEFINED > OBJ, PATCHES, INIC1 > OBJ, NAME, B6 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 2.775000E-02, 2.500000E-02, 2.914286E-02 > OBJ, GEOMETRY, cubet > OBJ, TYPE, PRESSURE_RELIEF > OBJ, PRES_RELIEF, 1000. ,0. STOP