PHOTON USE p gr out z 1 text 1 1 Velocity vectors 0.22368E+04 0.24817E+04 CR set vec ref 40 vec z 1 text 3 Pressto continue 0.20368E+04 0.44265E+03 CR text cl pause red text 1 5 Particle trajectories 0.22368E+04 0.23514E+04 CR GEOMETRY READ T00001 GEOMETRY READ T00002 GEOMETRY READ T00003 GEOMETRY READ T00004 GEOMETRY READ T00005 GEOMETRY READ T00006 GEOMETRY READ T00007 GEOMETRY READ T00008 text 3 Press e to end 0.20368E+04 0.44265E+03 CR enduse TEXT (Spray dryer (BFC=T) : G201 TITLE * Declare reals defining problem specification REAL(VIN,HIN,CPGAS,TIN) CPGAS=1005;TIN=293;HIN=TIN*CPGAS;VIN=20.0 GROUP 3. X-direction grid specification NX=18 GROUP 4. Y-direction grid specification NY=26 GROUP 5. Z-direction grid specification NZ=1 GROUP 6. Body-fitted coordinates or grid distortion BFC=T;NONORT=T GSET(D,18,26,1,1.0000,1.0000,1.0000) GSET(P,P01,0.0010,0.0000,0.0000) GSET(P,P02,0.0010,2.5000,0.0000) GSET(P,P03,0.0010,3.2000,0.0000) GSET(P,P04,0.0010,3.3000,0.0000) GSET(P,P05,0.0010,3.6000,0.0000) GSET(P,P06,0.0010,4.1000,0.0000) GSET(P,P07,1.0000,3.9500,0.0000) GSET(P,P08,1.5000,3.7000,0.0000) GSET(P,P09,1.5000,3.6000,0.0000) GSET(P,P10,1.5000,3.3000,0.0000) GSET(P,P11,1.5000,3.2000,0.0000) GSET(P,P12,1.5000,2.5000,0.0000) GSET(P,P13,0.5000,0.0000,0.0000) GSET(L,L01,P01,P02,8,1.0) GSET(L,L02,P02,P03,6,-1.5) GSET(L,L03,P03,P04,4,1.0) GSET(L,L04,P04,P05,4,1.5) GSET(L,L05,P05,P06,4,1.0) GSET(L,L06A,P06,P08,18,1.2,ARC,P07) GSET(L,L07,P08,P09,4,1.0) GSET(L,L08,P09,P10,4,-1.5) GSET(L,L09,P10,P11,4,1.0) GSET(L,L10,P11,P12,6,1.5) GSET(L,L11,P12,P13,8,1.0) GSET(L,L12,P13,P01,18,-1.5) GSET(L,L13,P05,P09,18,1.5) GSET(L,L14,P02,P12,18,1.5) GSET(F,F01,P01,-,P13,-,P12,-,P02,-) GSET(F,F02,P02,-,P12,P11.P10,P09,-,P05,P04.P03) GSET(F,F03,P05,-,P09,-,P08,-,P06,-) GSET(M,F01,+I+J,1,1,1,TRANS) GSET(M,F02,+I+J,1,9,1,TRANS) GSET(M,F03,+I+J,1,23,1,TRANS) GSET(C,K2,F,K1,1,18,1,26,RY,-0.1,0.0,0.0) GROUP 7. Variables stored, solved & named STORE(UCRT,VCRT,WCRT,TMP1);NCRT=1 SOLVE(P1,U1,V1,H1) GROUP 8. Terms (in differential equations) & devices GROUP 9. Properties of the medium (or media) * Air Density RHO1=1.24 * Air Laminar Kinematic Viscosity ENUL=1.8E-5/1.24 * Air Turbulent Kinematic Viscosity ENUT=1000*ENUL * Temperature Derivation TMP1=GRND2;TMP1A=0.0;CP1=CPGAS * Air and Vapour Prandtl Numbers PRNDTL(H1)=0.71 GROUP 11. Initialization of variable or porosity fields CONPOR(0.0,CELL,1,3,16,16,1,NZ) GROUP 13. Boundary conditions and special sources * Inlet Boundary Condition PATCH(INLET,NORTH,1,3,NY,NY,1,NZ,1,LSTEP) COVAL(INLET,P1,FIXFLU,RHO1*VIN) COVAL(INLET,V1,ONLYMS,-VIN) COVAL(INLET,H1,ONLYMS,HIN) COVAL(INLET,U1,ONLYMS,0.0) COVAL(INLET,C6,ONLYMS,0.0) * Outlet Boundary Condition PATCH(GXIT,SOUTH,1,7,1,1,1,NZ,1,LSTEP) COVAL(GXIT,P1,1.E+4,0.0);COVAL(GXIT,H1,ONLYMS,SAME) COVAL(GXIT,V1,ONLYMS,SAME);COVAL(GXIT,U1,ONLYMS,SAME) GROUP 15. Termination of sweeps * Number of Sweeps LSWEEP=20 * Number of Iterations LITER(P1)=10;LITER(U1)=1;LITER(V1)=1;LITER(W1)=1 VARMIN(C6)=1.E-10;VARMAX(C6)=1. GROUP 16. Termination of iterations GROUP 17. Under-relaxation devices RELAX(P1,LINRLX,.6);RELAX(V1,FALSDT,10/VIN/NY*25) RELAX(U1,FALSDT,10/VIN/NY*25);RELAX(C6,LINRLX,0.2) GROUP 22. Spot-value print-out IXMON=NX/2;IYMON=NY/2;IZMON=NZ/2 GROUP 23. Field print-out and plot control NPRINT=LSWEEP #GENTRA *------------------------------------------------------- * GENTRA GROUP 1: Particle physics *------------------------------------------------------- * Particle type - PARTICLES_IN_ISOTHERMAL_FLOW GPTYPE=30 * Gravity components/buoyancy GGRAX=0;GGRAY=-9.8;GGRAZ=0.0;GBUOYA=F * Stochastic model of turbulence GSTOCH=F * Data for particles in isothermal flow GDRAG=-10120 *------------------------------------------------------- * GENTRA GROUP 2: Boundary conditions for particles *------------------------------------------------------- * Inlet data -------------------------------------- *X Y U V DIA RHO MASFLOW 0.2 3.20E+00 9.04E+00 -1.06E-01 2.E-04 1000. 1.61E-06 0.2 3.20E+00 4.68E+00 -1.64E+00 2.E-04 1000. 8.85E-07 0.2 3.20E+00 6.60E+00 -6.47E-01 2.E-04 1000. 1.18E-06 0.2 3.20E+00 7.36E+00 -3.18E-01 2.E-04 1000. 1.32E-06 0.2 3.20E+00 9.04E+00 1.06E-01 2.E-04 1000. 1.61E-06 0.2 3.20E+00 7.36E+00 3.18E-01 2.E-04 1000. 1.32E-06 0.2 3.20E+00 6.60E+00 6.47E-01 2.E-04 1000. 1.18E-06 0.2 3.20E+00 4.68E+00 1.64E+00 2.E-04 1000. 8.85E-07 --------------------------------------- * Wall treatment, and rest coefficient if appropriate GWALLC=3;GWREST=1 * Porosity threshold GPOROS=0 *------------------------------------------------------- * GENTRA GROUP 3: Numerical controls *------------------------------------------------------- * 1st GENTRA sweep; frequency of calls GSWEP1=20;GSWEPF=1 * Maximum Lagrangian time-step GDTMAX=1 * Min # of t-steps per cell; max # of t-steps; timeout GLAGTS=5;GSTEMX=1000;GTIMMX=10 *------------------------------------------------------- * GENTRA GROUP 4: Output controls *------------------------------------------------------- * Trajectory-file identifier, history-file identifier and frequency for output GT1STC=T;GSWOUT=10;NGWSTR=1;NGWEND=20;NGWINT=1 GOUTFR=1 *------------------------------------------------------- * PIL provisions for the GENTRA-Earth run *------------------------------------------------------- * Storage for interphase sources of momentum STORE(MOMX);STORE(MOMY);store(prps) * Interphase sources of heat and momentum PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP) COVAL(GENPAT,3,FIXFLU,GRND) COVAL(GENPAT,5,FIXFLU,GRND) * Relaxation of interphase sources RELAX(MOMX,LINRLX,7.0000E-01) RELAX(MOMY,LINRLX,7.0000E-01) L(GENSET * ---------------------------------------------------------- * GENTRA GROUP 7: GENTRA test * ---------------------------------------------------------- * GENTRA auto test RG(51)--RG(56) are reserved RG(51)= 2.03E+01;RG(52)= 1.57E+03;RG(53)= 1.03E-01 RG(54)= 3.47E-01;RG(55)=-1.37E-04;RG(56)= 4.00E+00 RG(57)= 8.00E+16;RG(58)= 1.60E+00 EX(P1)=3.403E+1;EX(U1)=9.810E-1;EX(V1)=3.389;EX(H1)=2.926E+5 EX(TMP1)=2.911E+02 EX(PRPS)=2.269E+00 CP1 = 1.005000E+03