TALK=T;RUN( 1, 1)
 ** LOAD(x505) from the x Input Library
 
  PHOTON USE
  p
 
 
 
 
  use patgeo
  text
  1
  1
  Vapour mass fraction
   0.22368E+04 0.26119E+04 CR
  con vapo z 1 fill
  0.01
  text
  3
 
  Press  to continue
   0.20368E+04 0.44265E+03 CR
  pause
  con cl
  text cl
  red
  text
  1
  1
  Velocity vectors
   0.22368E+04 0.24817E+04 CR
  set vec ref
  40
  vec z 1
  text
  1
  5
  Particle trajectories
   0.22368E+04 0.23514E+04 CR
  use genuse
  text cl
  vec off
  con off
  text
  3
 
  Press  to continue
   0.20368E+04 0.44265E+03 CR
  pause
  red
  text
  1
  5
  Particle trajectories
   0.22368e+04 0.24817E+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  to continue
   0.20368E+04 0.44265E+03 CR
  pause
  text cl
  red
  text
  1
  1
  Mass deposited by particles
   0.21368E+04 0.23514E+04 CR
  con mass z 1 fill
  0.005
  text
  3
 
  Press e to end
   0.22368E+04 0.44265E+03 CR
  enduse
 
TEXT(Evapor. particl.in spray dryer:G505
TITLE
  DISPLAY
  Evaporating particles in spray dryer (BFC)
  ENDDIS
     * Declare reals defining problem specification
REAL(VIN,HIN,CPGAS,TIN)
CPGAS=1005.0;TIN=473;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
NAME(C6)=VAPO
STORE(UCRT,VCRT,WCRT,DEN1);NCRT=1
SOLVE(P1,U1,V1,H1,VAPO)
STORE(HEAT,MASS,TMP1,SPH1)
  Echo InForm settings for Group  7
  inform7begin
 (STORED VAR TDEC IS TMP1-273.15)
  inform7end
    GROUP 9. Properties of the medium (or media)
     * Air Density
PRESS0=1.01325E+05
     * Air Laminar Kinematic Viscosity
ENUL=1.8E-5/1.24
     * Air Turbulent Kinematic Viscosity
ENUT=1000*ENUL
     * Temperature Derivation
TMP1=GRND2;CP1=CPGAS
     * Air and Vapour Prandtl Numbers
PRNDTL(H1)=0.71
    GROUP 10. Inter-phase-transfer processes and properties
    GROUP 11. Initialization of variable or porosity fields
CONPOR(0.0,CELL,1,3,16,16,1,NZ)
     * Nozzle Blockage
     * Initialization of Enthalpy and Temperature (Commented Out)
FIINIT(H1)=HIN;FIINIT(DEN1)=7.374E-01;FIINIT(MASS)=1.e-10
FIINIT(TMP1)=TIN
FIINIT(VAPO)=0.0
RHO1=GRND5;DRH1DP=GRND5;RHO1C=0.714
RHO1A   = 0.000000E+00 ;RHO1B  = 3.530020E-03
  inform9begin
 (PROPERTY RHO1 is (P1+:PRESS0:)/((175.2*VAPO+286.7)*TMP1) with IMA$
T<100)                                                              
 (PROPERTY CP1 is (1910.*VAPO+:CPGAS:*(1-VAPO)) with IMAT<100)
 (PROPERTY TMP1 is (H1/(1910.*VAPO+:CPGAS:*(1.-VAPO))) with IMAT<10$
0)                                                                  
  inform9end
    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,7.374E-01*VIN)
COVAL(INLET,V1,ONLYMS,-VIN);COVAL(INLET,H1,ONLYMS,HIN)
COVAL(INLET,U1,ONLYMS,0.0);COVAL(INLET,VAPO,ONLYMS,0.0)
     * Outlet Boundary Condition
PATCH(GXOUT,SOUTH,1,7,1,1,1,NZ,1,LSTEP)
COVAL(GXOUT,P1,1.,0.0);COVAL(GXOUT,H1,ONLYMS,SAME)
COVAL(GXOUT,VAPO,ONLYMS,SAME);COVAL(GXOUT,V1,ONLYMS,SAME)
COVAL(GXOUT,U1,ONLYMS,SAME)
    GROUP 15. Termination of sweeps
     * Number of Sweeps
LSWEEP=500
     * Number of Iterations
LITER(P1)=20;LITER(U1)=1;LITER(V1)=1;LITER(W1)=1
VARMIN(VAPO)=0.0;VARMAX(VAPO)=1.
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,.3)
RELAX(V1,FALSDT,.1);RELAX(U1,FALSDT,.1)
RELAX(H1,FALSDT,0.1)
RELAX(VAPO,FALSDT,0.5)
    GROUP 22. Spot-value print-out
IXMON=NX/2;IYMON=NY/2;IZMON=NZ/2
TSTSWP=-1
OUTPUT(TMP1,Y,N,Y,Y,Y,Y)
    GROUP 23. Field print-out and plot control
  * ----------------------------------------------------------
  * Activate GENTRA
GENTR=T
  * ----------------------------------------------------------
  * GENTRA GROUP 0: PIL variables and default values
  * ----------------------------------------------------------
INTEGER(IX1,IX2,IX3,IX4,IX5,IX6,IX7,IX8,IX9)
CHAR(CX1,CX2,CX3);REAL(RX1);BOOLEAN(LX1)
CHAR(CDRAG,CKONC,CNUSS,CCPLIQ,CLHEAS,CSOLIN,CSOLST)
CHAR(CKONV,CCPCON,CCPVAP,CMWCON,CMWVAP,CLIQST,CLATVP)
CHAR(CVAPST,CSTPRE,CDTRCT,CHLIQD,CCPSOL,CHSOLD)
INTEGER(GPTYPE,GWALLC);GPTYPE=30;GWALLC=3
REAL(GWREST,GPOROS);GWREST=1.0;GPOROS=0.0
REAL(GGRAX,GGRAY,GGRAZ);GGRAX=0.0;GGRAY=0.0;GGRAZ=0.0
BOOLEAN(GBUOYA,GSTOCH,GSURPR);GBUOYA=F;GSURPR=F;GSTOCH=F
REAL(GDRAG,GNUSS);GDRAG=GRND1;GNUSS=GRND1
REAL(GLHEAS,GSOLIN);GLHEAS=3.335E+05;GSOLIN=1
REAL(GSOLST,GLIQST);GSOLST=273.15;GLIQST=273.15
REAL(GKONC,GKONV,GCPCON,GCPVAP,GCPLIQ,GMWCON,GMWVAP,GHLIQD)
REAL(GCPSOL,GDBUG1)
GKONC=0.0263;GKONV=GRND1;GCPCON=1007;GCPVAP=GRND1;GDBUG1=0.0
GHLIQD=GRND1;GCPLIQ=4131.8;GMWCON=28.9;GMWVAP=18;GCPSOL=4131.8
REAL(GVAPST,GSTPRE,GLATVP);GVAPST=GRND1;GSTPRE=GRND1;GLATVP=GRND1
REAL(GDTRCT);GDTRCT=0.0
CHAR(GINFIL);GINFIL='Q1'
INTEGER(GINSYS);GINSYS=12
INTEGER(GSWEP1,GSWEPF,GLAGTS);GSWEP1=1;GSWEPF=1;GLAGTS=5
REAL(GTIMMX);GTIMMX=0.0;INTEGER(GSTEMX);GSTEMX=0
REAL(GDTMAX,GRTFRC);GDTMAX=1.0;GRTFRC=1.0
REAL(GLNRLX);GLNRLX=1.0
CHAR(GRSFIL,GHFILE);GRSFIL=NONE
GHFILE='GHIS'
INTEGER(GOUTFR,GRESTI);GOUTFR=1;GRESTI=0
INTEGER(NGWSTR,NGWEND,NGWINT,GSWOUT);CHAR(GT1STC,GH1STC)
NGWSTR=0;NGWEND=0;NGWINT=1;GSWOUT=1
GT1STC=NONE;GH1STC=NONE
REAL(GFLENG,GTSTRT,GTEND)
CHAR(GENSET);GENSET=$G002
 ** LOAD(x  1) from the x Input Library
 ** LOAD(x  1) from the x Input Library
  *-------------------------------------------------------
  * GENTRA GROUP 1: Particle physics
  *-------------------------------------------------------
  * Particle type - 60 - VAPORISING_DROPLETS
GPTYPE=60
  * Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0;GGRAY=-9.8000E+00;GGRAZ=0;GBUOYA=F
  * Stochastic model of turbulence
GSTOCH=F
  * Data for vaporising particles
GDRAG=GRND1;GNUSS=GRND1
GKONV=GRND1;GCPCON=CPGAS;GCPVAP=GRND1
GMWCON=29;GMWVAP=18;GLIQST=GRND1
GCPLIQ=GRND1;GHLIQD=GRND1;GVAPST=GRND1
GLATVP=GRND1;GDTRCT=2.E-07
  *-------------------------------------------------------
  * GENTRA GROUP 2: Boundary conditions for particles
  *-------------------------------------------------------
  * Inlet-data file-name
GINFIL=Q1
  * Use this space for inlet-data table if GINFIL=Q1
  
  *__XP__YP__UP__VP__DIAM_LIQDEN_FRATE_TEMP_(NUMB)
   0.2  3.20E+00  9.04E+00  8.06E+00  6.0E-05  1000.  5.E-05   300.
   0.2  3.20E+00  7.36E+00  3.18E+00  1.1E-04  1000.  5.E-05   300.
   0.2  3.20E+00  6.60E+00  6.47E+00  1.1E-04  1000.  5.E-05   300.
   0.2  3.20E+00  4.68E+00  5.64E+00  1.1E-04  1000.  5.E-05   300.
   0.2  3.20E+00  9.04E+00 -1.06E+00  0.7E-04  1000.  5.E-05   300.
   0.2  3.20E+00  4.68E+00 -1.64E+00  1.1E-04  1000.  5.E-05   300.
   0.2  3.20E+00  6.60E+00 -6.47E+00  1.1E-04  1000.  5.E-05   300.
   0.2  3.20E+00  7.36E+00 -3.18E+00  1.1E-04  1000.  5.E-05   300.
  
  * 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=50;GSWEPF=1
  * Maximum Lagrangian time-step; time step size multplier
GDTMAX=4.e-4;GRTFRC=1.0
  * Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=2500;GTIMMX=5
  *-------------------------------------------------------
  * 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=10;GH1STC=H
  * Particle number for residence-time calculation
GRESTI=0
  *-------------------------------------------------------
  * GENTRA GROUP 5: Provisions for the GENTRA-Earth run
  *-------------------------------------------------------
  * Storage of Cartesian components in BFC cases
  * Storage for interphase sources of momentum
STORE(MOMX)
STORE(MOMY)
store(vpor,epor,npor)
  * Storage for interphase source of mass and for
    vapour mass-fraction
PRNDTL(VAPO)=1;PRT(VAPO)=1
  * Interphase sources of heat and momentum
PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENPAT,U1,fixflu,GRND)
COVAL(GENPAT,V1,fixflu,GRND)
COVAL(GENPAT,H1,fixflu,GRND)
COVAL(GENPAT,VAPO,fixflu,GRND)
  * Interphase sources of mass
PATCH(GENMAS,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENMAS,P1,FIXFLU,GRND)
  * Relaxation of interphase sources
GLNRLX=1.0000E-01
RELAX(MOMX,LINRLX,GLNRLX)
RELAX(MOMY,LINRLX,GLNRLX)
RELAX(HEAT,LINRLX,GLNRLX)
RELAX(MASS,LINRLX,GLNRLX)
  * Activation of GENTRA GROUND-station
  * ----------------------------------------------------------
  * GENTRA GROUP 6: Transmission of data to GENTRA-Earth
  * ----------------------------------------------------------
IG(20)=GPTYPE;IG(19)=GWALLC;IG(18)=GSWEP1;IG(17)=GSWEPF
IG(16)=GLAGTS;IG(15)=GSTEMX;IG(14)=GOUTFR;IG(13)=GRESTI
IG(12)=NGWSTR;IG(11)=NGWEND;IG(10)=NGWINT;IG(9)=GSWOUT
RG(100)=GGRAX;RG(99)=GGRAY;RG(98)=GGRAZ;RG(97)=GDRAG
RG(96)=GKONC;RG(95)=GNUSS;RG(94)=GLHEAS;RG(93)=GSOLIN
RG(92)=GSOLST;RG(91)=GLIQST;RG(90)=GKONV;RG(89)=GCPCON
RG(88)=GCPVAP;RG(87)=GMWCON;RG(86)=GMWVAP;RG(85)=GDBUG1
RG(84)=GWREST
RG(83)=GTIMMX;RG(82)=GCPLIQ;RG(81)=GPOROS;RG(80)=GDTMAX
RG(77)=GRTFRC;RG(76)=GHLIQD;RG(75)=GVAPST;RG(74)=GLATVP
RG(73)=GDTRCT;RG(72)=GCPSOL;RG(71)=GSTPRE
RG(70)=GTEND;RG(69)=GTSTRT;RG(68)=GFLENG
LG(20)=GBUOYA;LG(19)=GSTOCH;LG(18)=GSURPR
CG(10)=':GT1STC:'
CG(9) =':GH1STC:'
CG(8) =':GINFIL:'
CG(7) =':GRSFIL:'
CG(6) =':GHFILE:'
NOCOPY=T
 ** LOAD(x  2) from the x Input Library
 ** LOAD(x  2) from the x Input Library
  * ----------------------------------------------------------
  * GENTRA GROUP 7: GENTRA test
  * ----------------------------------------------------------
   * GENTRA auto test  RG(51)--RG(58) are reserved
RG(51)= 5.00E+00;RG(52)= 5.00E+03;RG(53)= 8.00E-03
RG(54)= 4.52E+00;RG(55)= 5.20E+00;RG(56)= 4.00E+00
RG(57)= 3.02E+00;RG(58)= 0.00E+00
DISTIL=T
EX(P1  )=3.469E+01;EX(U1  )=8.482E-01;EX(V1  )=5.525E+00
EX(H1  )=4.321E+05;EX(VAPO)=1.922E-02;EX(SPH1)=1.016E+03
EX(NPOR)=9.936E-01;EX(EPOR)=9.936E-01;EX(MOMY)=1.780E-05
EX(MOMX)=6.394E-06;EX(VPOR)=9.936E-01;EX(TMP1)=4.227E+02
EX(MASS)=8.816E-07;EX(HEAT)=1.512E+00;EX(DEN1)=8.180E-01
EX(WCRT)=4.218E-02;EX(VCRT)=5.470E+00;EX(UCRT)=8.428E-01
EX(TDEC)=1.513E+02
 LIBREF  =         505
STOP