TALK=T;RUN( 1, 1) ************************************************************ Q1 created by VDI menu, Version 2020, Date 16/06/20 CPVNAM=VDI; SPPNAM=FLAIR ************************************************************ Echo DISPLAY / USE settings ************************************************************ IRUNN = 1 ;LIBREF = 2 ************************************************************ Group 1. Run Title TEXT(Fire-Spray in a compartment ) ************************************************************ Group 2. Transience STEADY=F * Set overall time and no. of steps RSET(U,0.,120.,12) * Modify regions ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,31,21,17) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(128)=SLEN ;NAME(130)=ENUL NAME(137)=MASS ;NAME(138)=HEAT NAME(139)=MOMZ ;NAME(140)=MOMY NAME(141)=MOMX ;NAME(142)=EPKE NAME(143)=DEN1 ;NAME(144)=EL1 NAME(145)=ENUT ;NAME(146)=SPH1 NAME(147)=TEM1 ;NAME(149)=SMOK NAME(150) =MH2O * Solved variables list SOLVE(P1,U1,V1,W1,TEM1,SMOK,MH2O) * Stored variables list STORE(SPH1,ENUT,EL1,DEN1,EPKE,MOMX,MOMY,MOMZ) STORE(HEAT,MASS,ENUL,SLEN) * Additional solver options SOLUTN(P1,Y,Y,Y,N,N,Y) SOLUTN(TEM1,Y,Y,Y,N,N,Y) SOLUTN(SMOK,Y,Y,Y,N,N,Y) SOLUTN(MH2O,Y,Y,Y,N,N,N) TURMOD(KECHEN) ************************************************************ Group 8. Terms & Devices 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) 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(P1)=0. ;FIINIT(U1)=0. FIINIT(V1)=0. ;FIINIT(W1)=0. FIINIT(KE)=1.0E-03 ;FIINIT(EP)=9.0E-03 FIINIT(TEM1)=20. ;FIINIT(SMOK)=0. No PATCHes used for this Group INIADD = F NAMFI ='stdy' ************************************************************ 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, 12) COVAL(BUOYANCY, U1, FIXFLU, GRND2) COVAL(BUOYANCY, V1, FIXFLU, GRND2) COVAL(BUOYANCY, W1, FIXFLU, GRND2) PATCH(GENPAT, CELL, 0, 0, 0, 0, 0, 0, 1, 12) COVAL(GENPAT, U1, FIXFLU, GRND) COVAL(GENPAT, V1, FIXFLU, GRND) COVAL(GENPAT, W1, FIXFLU, GRND) COVAL(GENPAT, TEM1, FIXFLU, GRND) COVAL(GENPAT, MH2O, FIXFLU, GRND) PATCH(GENMAS, CELL, 0, 0, 0, 0, 0, 0, 1, 12) COVAL(GENMAS, P1, FIXFLU, GRND) BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81 BUOYD =1.204938 EGWF = T ************************************************************ Echo save-block settings for Group 13 SAVE13BEGIN SAVE13END ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 100 RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=200 ;LITER(TEM1)=50 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX,1. ) RELAX(KE ,LINRLX,0.5 ) RELAX(EP ,LINRLX,0.5 ) RELAX(MASS,LINRLX,0.3 ) RELAX(HEAT,LINRLX,0.3 ) RELAX(MOMZ,LINRLX,0.3 ) RELAX(MOMY,LINRLX,0.3 ) RELAX(MOMX,LINRLX,0.3 ) RELAX(DEN1,LINRLX,0.3 ) RELAX(TEM1,LINRLX,0.25 ) RELAX(SMOK,LINRLX,0.5 ) RELAX(MH2O,LINRLX,0.3 ) KELIN = 3 ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-7.5E+04 VARMAX(TEM1)=3000. ;VARMIN(TEM1)=8. VARMAX(SMOK)=1. ;VARMIN(SMOK)=0. VARMAX(MH2O)=1. ;VARMIN(MH2O)=0. ************************************************************ Group 19. EARTH Calls To GROUND Station NAMGRD =FLAR GENK = T PARSOL = F CONWIZ = T IENUTA = 2 ISG62 = 1 SPEDAT(SET,FACETDAT,NUMOBJ,I,14) SPEDAT(SET,SLEN,ACON,R,3.) SPEDAT(SET,SLEN,BCON,R,150.) SPEDAT(SET,SLEN,CCON,R,30.) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,IMMERSOL,EMCON,R,5.) SPEDAT(SET,FLAIR,HCOMB,R,1.3E+07) SPEDAT(SET,FLAIR,RADFAC,R,0.) SPEDAT(SET,FLAIR,SYIELD,R,0.02) SPEDAT(SET,FLAIR,STOIC,R,0.992366) SPEDAT(SET,FLAIR,KEXT,R,7600.) SPEDAT(SET,FLAIR,COYIELD,R,0.06) SPEDAT(SET,MAXINC,TEM1,R,10.) SPEDAT(SET,MAXINC,MH2O,R,0.5) L($G001 GENTR = T *------------------------------------------------------ * GENTRA GROUP 1: Particle physics *------------------------------------------------------ * Particle type - 60 GPTYPE = 60 * Gravity components in GENTRA Cartesian system GGRAX =0. ;GGRAY =0. ;GGRAZ =-9.81 * Buoyancy forces GBUOYA = T ;GSURPR = F * Stochastic model of turbulence GSTOCH = F * Data for vaporising particles * Data for melting/solidifying particles GDRAG = GRND1 ;GKONC =0.0263 ;GNUSS = GRND1 GKONV = GRND1 ;GCPCON =1007. ;GCPVAP = GRND1 GMWCON =28.9 ;GMWVAP =18. ;GLIQST =273.149994 GCPLIQ =4131.799805 ;GHLIQD = GRND1 ;GVAPST = GRND1 GLATVP = GRND1 ;GSTPRE = GRND1 ;GDTRCT =0. *------------------------------------------------------- * GENTRA GROUP 2: Boundary conditions for particles *------------------------------------------------------- * Inlet-data file-name GINFIL ='spin' * Wall treatment, and rest coefficient if appropriate GWALLC = 1 * Porosity threshold GPOROS =0. GTSTRT =0. ;GTEND =0. ;GFLENG =0. *------------------------------------------------------- * GENTRA GROUP 3: Numerical controls *------------------------------------------------------- * 1st GENTRA sweep; frequency of calls GSWEP1 = 1 ;GSWEPF = 1 * Maximum Lagrangian time-step; time step size multplier GDTMAX =1. ;GRTFRC =1. * Min # of t-steps per cell; max # of t-steps; timeout GLAGTS = 5 ;GSTEMX = 100 GTIMMX =100. *------------------------------------------------------- * GENTRA GROUP 4: Output controls *------------------------------------------------------- * Restart-file, history-file and frequency for output GRSFIL ='NONE' GHFILE ='GHIS' GOUTFR = 1 * The identifier of the individual history and trajectory files GH1STC ='NONE' GT1STC ='NONE' * The first, last particles and the interval for writing history and trajectory files NGWSTR = 0 ;NGWEND = 0 ;NGWINT = 1 GSWOUT = 1 * Particle number for residence-time calculation GRESTI = 0 L(GENSET ************************************************************ Group 20. Preliminary Printout ************************************************************ Group 21. Print-out of Variables OUTPUT(ENUL,Y,N,N,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 14 ;IYMON = 11 ;IZMON = 8 NPRMON = 100000 NPRMNT = 1 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 NTPRIN = 1 ;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.15279,-0.927805,0.340344) GVIEW(UP,-0.053857,0.336058,0.9403) GVIEW(VDIS,2.946434) GVIEW(CENTRE,2.5,1.4,1.1) > DOM, SIZE, 5.000000E+00, 2.800000E+00, 2.200000E+00 > DOM, MONIT, 1.592500E+00, 1.400000E+00, 9.259092E-01 > DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00 > GRID, MINCELL, 2.5E-03 ,5.0E-03 ,5.0E-03 > GRID, BOUNDS, T F T T T T > GRID, RSET_X_1, -10, 1.200000E+00 > GRID, RSET_X_2, -2, 1.000000E+00 > GRID, RSET_X_3, -1, 1.000000E+00 > GRID, RSET_X_4, -2, 1.000000E+00 > GRID, RSET_X_5, -10, 1.200000E+00 > GRID, RSET_X_6, -1, 1.000000E+00 > GRID, RSET_X_7, 5, 1.500000E+00 > GRID, RSET_Y_1, -7, 1.000000E+00 > GRID, RSET_Y_2, -1, 1.000000E+00 > GRID, RSET_Y_3, -2, 1.000000E+00 > GRID, RSET_Y_4, -1, 1.000000E+00 > GRID, RSET_Y_5, -2, 1.000000E+00 > GRID, RSET_Y_6, -1, 1.000000E+00 > GRID, RSET_Y_7, -7, 1.200000E+00,G > GRID, RSET_Z_1, 3, 1.000000E+00 > GRID, RSET_Z_2, 11, 1.000000E+00 > GRID, RSET_Z_3, 2, 1.000000E+00 > GRID, RSET_Z_4, 1, 1.000000E+00 > DOM, INI_AMB, YES > DOM, INI_BUOY, YES > OBJ, NAME, BLK1 > OBJ, POSITION, 2.800000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, 1.030000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 199,Solid allowing fluid-slip at walls > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, BLK2 > OBJ, POSITION, 2.800000E+00, AT_END, 0.000000E+00 > OBJ, SIZE, 1.000000E-01, 1.050000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 199,Solid allowing fluid-slip at walls > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, BLK3 > OBJ, POSITION, 2.800000E+00, 1.030000E+00, AT_END > OBJ, SIZE, 1.000000E-01, 7.200000E-01, 3.700000E-01 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, BLOCKAGE > OBJ, MATERIAL, 199,Solid allowing fluid-slip at walls > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, NAME, OPN1 > OBJ, POSITION, 2.900000E+00, AT_END, 0.000000E+00 > OBJ, SIZE, TO_END, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, TYPE, OPENING > OBJ, PRESSURE, P_AMBIENT > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 > OBJ, NAME, OPN2 > OBJ, POSITION, 2.900000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, OPENING > OBJ, PRESSURE, P_AMBIENT > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 > OBJ, NAME, OPN3 > OBJ, POSITION, AT_END, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, OPENING > OBJ, PRESSURE, P_AMBIENT > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 > OBJ, NAME, OPN4 > OBJ, POSITION, 2.900000E+00, 0.000000E+00, AT_END > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, OPENING > OBJ, PRESSURE, P_AMBIENT > OBJ, TEMPERATURE, T_AMBIENT > OBJ, COEFFICIENT, 1. > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, TURBULENCE, 1.0E-03 ,9.0E-03 > OBJ, NAME, WALL1 > OBJ, POSITION, 0.000000E+00, AT_END, 0.000000E+00 > OBJ, SIZE, 2.800000E+00, 0.000000E+00, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, SURF_TEMP, 0. ,20. > OBJ, NAME, WALL2 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 2.800000E+00, 0.000000E+00, 2.200000E+00 > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, SURF_TEMP, 0. ,20. > OBJ, NAME, WALL3 > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, 0.000000E+00, TO_END, TO_END > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, SURF_TEMP, 0. ,20. > OBJ, NAME, FLOOR > OBJ, POSITION, 0.000000E+00, 0.000000E+00, 0.000000E+00 > OBJ, SIZE, TO_END, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, SURF_TEMP, 0. ,20. > OBJ, NAME, CEILING > OBJ, POSITION, 0.000000E+00, 0.000000E+00, AT_END > OBJ, SIZE, 2.800000E+00, TO_END, 0.000000E+00 > OBJ, DOMCLIP, NO > OBJ, GEOMETRY, box > OBJ, VISIBLE, NO > OBJ, TYPE, PLATE > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, SURF_TEMP, 0. ,20. > OBJ, NAME, FIRE > OBJ, POSITION, 1.250000E+00, 1.150000E+00, 0.000000E+00 > OBJ, SIZE, 5.000000E-01, 5.000000E-01, 3.000000E-01 > OBJ, GEOMETRY, box > OBJ, TYPE, FIRE > OBJ, INFSRC_P1, (exp(-500.*MH2O)*46.9*tim)*tim*(1+rox)/hcom$ > OBJ, INFSRC_P1, b with whol > OBJ, INFSRC_TEM1, (exp(-500.*MH2O)*46.9*tim)*tim with whol > OBJ, INFSRC_SMOK, 1 with onlyms > OBJ, TIME_LIMITS, ALWAYS_ACTIVE > OBJ, PRE-TEMP, T_AMBIENT > OBJ, MASS-SOURCE, No Source > OBJ, HEAT-SOURCE, Fixed Power > OBJ, FXD-FLU, 0. > OBJ, SCALAR-SOURCE, No Source > OBJ, INLET_SMOK, 1. > OBJ, NAME, SPRAY > OBJ, POSITION, 1.440000E+00, 1.350000E+00, 2.100000E+00 > OBJ, SIZE, 1.000000E-01, 1.000000E-01, 0.000000E+00 > OBJ, GEOMETRY, spray > OBJ, TYPE, SPRAY-HEAD > OBJ, TIME_SPRAY, 120. > OBJ, LINK_T, YES > OBJ, T_ACTIVATE, 50. > OBJ, RTI, 160. > OBJ, AXIS, Z > OBJ, ORIGIN, 1.49 ,1.4 ,2.1 > OBJ, RADIUS, 0.05 > OBJ, PORTS, 8. > OBJ, TOTAL-VOL, 60. > OBJ, TOTAL-VEL, 1. > OBJ, ANGLE, 135. > OBJ, TEMPERATURE, 8. > OBJ, DMEAN, 0.8 > OBJ, NSIZE, 1. > OBJ, DMIN, 0. > OBJ, DMAX, 2.5 > OBJ, SPREAD_EXP, 1.7 STOP