TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2020, Date 16/06/20 CPVNAM=VDI; SPPNAM=FLAIR ************************************************************ Echo DISPLAY / USE settings DISPLAY Library case i214: Tutorial: Fire and Smoke Modelling No jetfans ENDDIS ************************************************************ IRUNN = 1 ;LIBREF = 214 ************************************************************ Group 1. Run Title TEXT(Fire and Smoke Modelling ) ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,300.,30) * Modify regions ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,28,28,22) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(138)=SPH1 ;NAME(139)=SMOK NAME(140)=SLN2 ;NAME(141)=SLEN NAME(142)=COPM ;NAME(143)=PPPM NAME(144)=PPM ;NAME(145)=SDEN NAME(146)=TEM1 ;NAME(147)=EPKE NAME(148)=DEN1 ;NAME(149)=EL1 NAME(150) =ENUT * Solved variables list SOLVE(P1,U1,V1,W1,SMOK,TEM1) * Stored variables list STORE(ENUT,EL1,DEN1,EPKE,SDEN,PPM,PPPM,COPM) STORE(SLEN,SLN2,SPH1) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(SMOK,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices NEWRH1 = T NEWENL = T ************************************************************ Group 9. Properties PRESS0 =1.0E+05 ;TEMP0 =273. * Domain material index is 2 signifying: * Air using Ideal Gas Law, STP SETPRPS(1, 2) DRH1DP = GRND5 DVO1DT =3.33E-03 PRNDTL(TEM1)=-0.0263 PRT(KE)=0.75 ;PRT(EP)=1.15 ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(SMOK)=0. ;FIINIT(TEM1)=20. 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, 30) COVAL(BUOYANCY, U1, FIXFLU, GRND2) COVAL(BUOYANCY, V1, FIXFLU, GRND2) COVAL(BUOYANCY, W1, FIXFLU, GRND2) BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81 BUOYD =1.189 EGWF = T ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 50 RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(KE)=10 LITER(EP)=10 ;LITER(SMOK)=10 LITER(TEM1)=50 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,0.5 ) RELAX(U1 ,LINRLX,0.5 ) RELAX(V1 ,LINRLX,0.5 ) RELAX(W1 ,LINRLX,0.5 ) RELAX(KE ,LINRLX,0.5 ) RELAX(EP ,LINRLX,0.5 ) RELAX(SMOK,LINRLX,0.5 ) RELAX(TEM1,LINRLX,0.25 ) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-7.5E+04 VARMAX(SMOK)=1. ;VARMIN(SMOK)=0. VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.75 ************************************************************ Group 19. EARTH Calls To GROUND Station NAMGRD =FLAR GENK = T PARSOL = F CONWIZ = T IENUTA = 2 ISG62 = 1 SPEDAT(SET,SLEN,ACON,R,3.) SPEDAT(SET,SLEN,CCON,R,30.) SPEDAT(SET,SLN2,ACON,R,8.) SPEDAT(SET,SLN2,CCON,R,30.) SPEDAT(SET,KECONST,C2E,R,1.92) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,KECONST,C1E,R,1.44) SPEDAT(SET,IMMERSOL,EMCON,R,5.) SPEDAT(SET,FLAIR,HCOMB,R,1.3E+07) SPEDAT(SET,FLAIR,RADFAC,R,0.333333) SPEDAT(SET,FLAIR,SYIELD,R,0.022) SPEDAT(SET,FLAIR,STOIC,R,0.992366) SPEDAT(SET,FLAIR,KEXT,R,7600.) SPEDAT(SET,FLAIR,COYIELD,R,0.02) SPEDAT(SET,MAXINC,TEM1,R,10.) ************************************************************ Group 20. Preliminary Printout ************************************************************ Group 21. Print-out of Variables ************************************************************ Group 22. Monitor Print-Out IXMON = 14 ;IYMON = 14 ;IZMON = 19 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.75052,0.511061,0.418971) GVIEW(UP,-0.363409,-0.210348,0.907572) GVIEW(VDIS,6.93928) GVIEW(CENTRE,5.,5.,1.5) > DOM, SIZE, 1.000000E+01, 1.000000E+01, 3.000000E+00 > DOM, MONIT, 4.869048E+00, 4.869048E+00, 2.070362E+00 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > DOM, INCREMENT, 1.000000E-02, 1.000000E-02, 1.000000E-02 > GRID, AUTO, T T T > GRID, MINCELL, 0.02 ,0.02 ,0.02 > GRID, RSET_X_1, -12, 1.100000E+00,G > GRID, RSET_X_2, -4, 1.100000E+00,G > GRID, RSET_X_3, -12, 1.100000E+00,G > GRID, RSET_Y_1, -12, 1.100000E+00,G > GRID, RSET_Y_2, -4, 1.100000E+00,G > GRID, RSET_Y_3, -12, 1.100000E+00,G > GRID, RSET_Z_1, -7, 1.100000E+00,G > GRID, RSET_Z_2, 15, 1.100000E+00,G > OBJ, NAME, FIRE > OBJ, POSITION, 4.500000E+00, 4.500000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+00, 1.000000E+00, 5.000000E-01 > OBJ, GEOMETRY, box > OBJ, TYPE, FIRE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, PRE-TEMP, 20. > OBJ, MASS-SOURCE, Heat related > OBJ, HEAT-SOURCE, Power of time > OBJ, COEFFS, 46.900002 ,2. ,0. > OBJ, QMAX, 2.0E+07 > OBJ, SCALAR-SOURCE, Heat Related > OBJ, INLET_SMOK, 1. > OBJ, NAME, FLOOR > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+01, 1.000000E+01, 0.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, SURF_TEMP, 0. ,20. > OBJ, NAME, OPEN1 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E+01, 0.000000E+00, 3.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, OPENING > OBJ, WIREFRAME, YES > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, 20. > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 > OBJ, NAME, OPEN2 > OBJ, POSITION, 1.000000E+01, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.000000E+01, 3.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, OPENING > OBJ, WIREFRAME, YES > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, 20. > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 > OBJ, NAME, OPEN3 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, 1.000000E+01, 3.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, OPENING > OBJ, WIREFRAME, YES > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, 20. > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 > OBJ, NAME, OPEN4 > OBJ, POSITION, 0.000000E+00, 1.000000E+01, 0.000000E+00 > OBJ, SIZE, 1.000000E+01, 0.000000E+00, 3.000000E+00 > OBJ, GEOMETRY, box > OBJ, TYPE, OPENING > OBJ, WIREFRAME, YES > OBJ, PRESSURE, 0. > OBJ, TEMPERATURE, 20. > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 STOP