TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 3.6, Date 23/09/04 CPVNAM=VDI;SPPNAM=FLAIR ************************************************************ Echo DISPLAY / USE settings PHOTON USE P;;;; msg Geometry of the problem use patgeo msg msg Velocity vectors vec x 1 sh; vec x m sh vec z m sh; vec z 1 sh; vec y m sh msg Press return to continue pause vec cl; red msg Surface of smoke concentration of 0.1 surf tem1 x col 1;500.35 surf tem1 y col 1;500.35 surf tem1 z col 1;500.35 surf c1 x col 6;.3 surf c1 y col 6;.3 surf c1 z col 6;.3 msg Type E to end ENDUSE DISPLAY Library case Y101: SINGLE ROOM WITH A FIRE This case models a simple three-dimensional room, consisting of four open windows, an open door, a partial wall in the middle and a fire in one corner. Turbulence is modelled by means of the k-e model. The problem is set up with a coarse grid; for more accurate results a finer grid and more sweeps are needed. The results show smoke and temperature distribution for steady-state conditions. (Use PHOTON to display results. Enter USE Q1 in PHOTON command mode, or create a file called U, which contains the one line: USE Q1 before typing RUNPHO.) ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 0 ************************************************************ Group 1. Run Title TEXT(Single room with a fire: I101 ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,21,10,10) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Non-default variable names NAME( 41) =EPOR ; NAME( 44) =ABSV NAME( 46) =BLOK ; NAME( 47) =PRPS NAME( 49) =DEN1 ; NAME( 50) =TEM1 * Solved variables list SOLVE(P1 ,U1 ,V1 ,W1 ,C1 ,TEM1) * Stored variables list STORE(DEN1,PRPS,BLOK,ABSV,EPOR) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,N) SOLUTN(TEM1,Y,Y,Y,N,N,Y) IVARBK = -1 ;ISOLBK = 1 TURMOD(KEMODL) ************************************************************ Group 8. Terms & Devices NEWRH1 = T NEWENL = T ISOLX = 0 ;ISOLY = 0 ;ISOLZ = 0 DENPCO = T ************************************************************ Group 9. Properties PRESS0 = 1.000000E+05 ;TEMP0 = 2.730000E+02 * Domain material is: 2 Air using Ideal Gas Law, STP SETPRPS(1, 2) DRH1DP = GRND5 DVO1DT = 3.330000E-03 PRNDTL(TEM1) = -2.630000E-02 * Non-default SOLUTN setting for KE SOLUTN(KE ,Y,Y,N,N,N,Y) * Non-default SOLUTN setting for EP SOLUTN(EP ,Y,Y,N,N,N,Y) PRT (EP ) = 1.314000E+00 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(KE ) = 1.000000E-02 ;FIINIT(EP ) = 1.000000E-03 FIINIT(C1 ) = 1.001000E-10 ;FIINIT(EPOR) = 1.000000E+00 FIINIT(ABSV) = 1.001000E-10 ;FIINIT(BLOK) = 1.000000E+00 FIINIT(PRPS) = -1.000000E+00 ;FIINIT(DEN1) = 1.001000E-10 FIINIT(TEM1) = 2.000000E+01 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,1) COVAL (BUOYANCY,U1 , FIXFLU , GRND2 ) COVAL (BUOYANCY,V1 , FIXFLU , GRND2 ) COVAL (BUOYANCY,W1 , FIXFLU , GRND2 ) BUOYA = 0.000000E+00 ; BUOYB =-9.810000E+00 BUOYC = 0.000000E+00 BUOYD = 1.189000E+00 EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 550 RESFAC = 1.000000E-02 ************************************************************ Group 16. Terminate Iterations LITER (P1 ) = 50 ;LITER (KE ) = 10 LITER (EP ) = 10 ;LITER (TEM1) = 50 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000000E+00) RELAX(KE ,LINRLX, 5.000000E-01) RELAX(EP ,LINRLX, 5.000000E-01) RELAX(DEN1,LINRLX, 3.000000E-01) KELIN = 1 ************************************************************ Group 18. Limits VARMAX(U1 ) = 1.000000E+06 ;VARMIN(U1 ) =-1.000000E+06 VARMAX(V1 ) = 1.000000E+06 ;VARMIN(V1 ) =-1.000000E+06 VARMAX(W1 ) = 1.000000E+06 ;VARMIN(W1 ) =-1.000000E+06 ************************************************************ Group 19. EARTH Calls To GROUND Station USEGRD = T ;USEGRX = T NAMGRD =FLAR GENK = T ASAP = T PARSOL = T LSG57 = T SPEDAT(SET,FLAIR,CLOTHING,R,9.30000E-02) SPEDAT(SET,FLAIR,RELHUMID,R,5.00000E+01) SPEDAT(SET,FLAIR,TRAD,R,2.50000E+01) SPEDAT(SET,FACETDAT,NUMOBJ,I,8) SPEDAT(SET,IMMERSOL,EMCON,R,5.00000E+00) ************************************************************ Group 20. Preliminary Printout ECHO = T NULLPR = F NDST = 0 DSTTOL = 1.000000E-02 ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 20 ;IYMON = 9 ;IZMON = 9 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NXPRIN = 4 NYPRIN = 1 NZPRIN = 1 NPLT = 1 ISWPRF = 1 ;ISWPRL = 100000 No PATCHes used for this Group ************************************************************ Group 24. Dumps For Restarts GVIEW(P,5.773503E-01,5.773503E-01,5.773503E-01) GVIEW(UP,-4.082483E-01,8.164966E-01,-4.082483E-01) > DOM, SIZE, 1.200000E+01, 2.500000E+00, 6.000000E+00 > DOM, MONIT, 1.112500E+01, 2.125000E+00, 5.125000E+00 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, SNAPSIZE, 1.000000E-02 > GRID, RSET_X_1, 2, 1.000000E+00 > GRID, RSET_X_2, 3, 1.000000E+00 > GRID, RSET_X_3, 3, 1.000000E+00 > GRID, RSET_X_4, 3, 1.000000E+00 > GRID, RSET_X_5, 3, 1.000000E+00 > GRID, RSET_X_6, 6, 1.000000E+00 > GRID, RSET_X_7, 1, 1.000000E+00 > GRID, RSET_Y_1, 2, 1.000000E+00 > GRID, RSET_Y_2, 2, 1.000000E+00 > GRID, RSET_Y_3, 3, 1.000000E+00 > GRID, RSET_Y_4, 1, 1.000000E+00 > GRID, RSET_Y_5, 2, 1.000000E+00 > GRID, RSET_Z_1, 2, 1.000000E+00 > GRID, RSET_Z_2, 2, 1.000000E+00 > GRID, RSET_Z_3, 1, 1.000000E+00 > GRID, RSET_Z_4, 2, 1.000000E+00 > GRID, RSET_Z_5, 2, 1.000000E+00 > GRID, RSET_Z_6, 1, 1.000000E+00 > OBJ, NAME, PART1 > OBJ, POSITION, 7.000000E+00, 0.000000E+00, 3.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.800000E+00, 3.000000E+00 > OBJ, GEOMETRY, cube11 > OBJ, ROTATION24, 1 > OBJ, TYPE, PLATE > OBJ, POROSITY, 0.000000E+00 > OBJ, SIDE, BOTH > OBJ, NAME, FIRE > OBJ, POSITION, 1.150000E+01, 0.000000E+00, 5.500000E+00 > OBJ, SIZE, 5.000000E-01, 5.000000E-01, 5.000000E-01 > OBJ, GEOMETRY, cubet1 > OBJ, ROTATION24, 1 > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, DOMAIN > OBJ, HEAT_FLVX, 0.000000E+00, 2.000000E+06 > OBJ, SCAL_FIXF, 2.000000E+00 > OBJ, FIX_VAL_C1, FIXVAL , 1.000000E+00 > OBJ, INI_C1, 0.000000E+00 > OBJ, NAME, WIN1 > OBJ, POSITION, 1.200000E+01, 1.000000E+00, 1.000000E+00 > OBJ, SIZE, 0.000000E+00, 8.000000E-01, 1.500000E+00 > OBJ, GEOMETRY, cube12t > OBJ, ROTATION24, 1 > OBJ, TYPE, OPENING > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, 2.000000E+01 > OBJ, COEFFICIENT, 1.000000E+03 > OBJ, OUTLET_C1, SAME > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, WIN2 > OBJ, POSITION, 1.200000E+01, 1.000000E+00, 4.000000E+00 > OBJ, SIZE, 0.000000E+00, 8.000000E-01, 1.500000E+00 > OBJ, GEOMETRY, cube12t > OBJ, ROTATION24, 1 > OBJ, TYPE, OPENING > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, 2.000000E+01 > OBJ, COEFFICIENT, 1.000000E+03 > OBJ, OUTLET_C1, SAME > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, WIN3 > OBJ, POSITION, 1.000000E+00, 1.000000E+00, 6.000000E+00 > OBJ, SIZE, 1.500000E+00, 8.000000E-01, 0.000000E+00 > OBJ, GEOMETRY, cube12t > OBJ, ROTATION24, 1 > OBJ, TYPE, OPENING > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, 2.000000E+01 > OBJ, COEFFICIENT, 1.000000E+03 > OBJ, OUTLET_C1, SAME > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, WIN4 > OBJ, POSITION, 4.000000E+00, 1.000000E+00, 6.000000E+00 > OBJ, SIZE, 1.500000E+00, 8.000000E-01, 0.000000E+00 > OBJ, GEOMETRY, cube12t > OBJ, ROTATION24, 1 > OBJ, TYPE, OPENING > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, 2.000000E+01 > OBJ, COEFFICIENT, 1.000000E+03 > OBJ, OUTLET_C1, SAME > OBJ, TURBULENCE, SAME , SAME > OBJ, NAME, DOOR > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 2.000000E+00, 1.000000E+00 > OBJ, GEOMETRY, cube12t > OBJ, ROTATION24, 1 > OBJ, TYPE, OPENING > OBJ, PRESSURE, 0.000000E+00 > OBJ, TEMPERATURE, 2.400000E+01 > OBJ, COEFFICIENT, 1.000000E+03 > OBJ, OUTLET_C1, SAME > OBJ, TURBULENCE, SAME , SAME STOP