AUTOPLOT USE
 
  file
  h00001 2
  h00002
  h00003
  h00004
  h00005
 
  data 1
  time diam
  data 2
  time diam
  data 3
  time diam
  data 4
  time diam
  data 5
  time diam
  pl
  red
  @
   0.84592E+03 0.15372E+04 CR
  Particle diameter vers time.@
  @
   0.12304E+04 0.22211E+04 CR
  particle 1@
  @
   0.20379E+04 0.21430E+04 CR
  particle 2@
  @
   0.28389E+04 0.20257E+04 CR
  particle 3@
  @
   0.28902E+04 0.23775E+04 CR
  particle 4@
  @
   0.29351E+04 0.26706E+04 CR
  particle 5@
  @
   0.29415E+04 0.14330E+03 CR
  time (second)@
  @
   0.13842E+04 0.82071E+03  P
  Press e to end.@
  ENDUSE
 
TEXT(ISOTHER. EVAPORN, CONST PROP'S:G508
TITLE
  DISPLAY
  Evaporation of a particle with conditions such that
  the particle temperature is constant.
  Constant gas temperature,no slip.
  Particle diameter as a function of time is given by:
 
  d = sqrt(d0**2 - 4.Nu.kg(Tg-Tp).t/rhop/L)
 
  where: d    - particle diameter
         d0   - particle initial diameter
         kg   - thermal conductivity of continuous phase
         L    - latent heat of evaporation
         Nu   - Nusselt number
         t    - time
         rhop - particle liquid density
         Tg   - gas temperature
         Tp   - particle temperature
  ENDDIS
REAL(WIN);WIN=50.0
STEADY=F;GRDPWR(T,5,2.0E-02,1.0)
GRDPWR(X,5,2.0,1.0)
SOLVE(P1,U1);SOLVE(TEM1)
STORE(PRPS);FIINIT(PRPS)=0.0
STORE(RHO1);STORE(ENUL)
PRESS0=1.0E+05
FIINIT(U1)=WIN;FIINIT(P1)=WIN;FIINIT(TEM1)=450.0
PATCH(GXIN,WEST,1,1,1,1,1,1,1,LSTEP)
COVAL(GXIN,P1,FIXFLU,WIN);COVAL(GXIN,U1,ONLYMS,WIN)
COVAL(GXIN,TEM1,ONLYMS,450.0)
PATCH(GXOUT,EAST,NX,NX,1,1,1,1,1,LSTEP)
COVAL(GXOUT,P1,FIXP,0.0)
LSWEEP=1
 
#GENTRA
  *-------------------------------------------------------
  * GENTRA GROUP 1: Particle physics
  *-------------------------------------------------------
  * Particle type - 60 - VAPORISING_DROPLETS
GPTYPE=60
  * Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0;GGRAY=0;GGRAZ=0;GBUOYA=F
  * Stochastic model of turbulence
GSTOCH=F
  * Data for vaporising particles
GDRAG=GRND1;GKONC=3.3400E-02;GNUSS=1.9060E+01
GKONV=2.5000E-02;GCPCON=1000;GCPVAP=2000
GMWCON=28;GMWVAP=18;GLIQST=GRND1
GCPLIQ=4200;GHLIQD=1.0000E+06;GVAPST=373
GLATVP=2.2500E+06;GSTPRE=19614;GDTRCT=0
  *-------------------------------------------------------
  * GENTRA GROUP 2: Boundary conditions for particles
  *-------------------------------------------------------
  * Inlet-data file-name
GINFIL=Q1
  * Use this space for inlet-data table if GINFIL=Q1
  
  *_XP_UP_DI_LDEN_MDOT_TEM_(NUM)_(STR)_(END)
   0.1 50.0 5e-5 1000 1e-6 350
  
  * 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=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=25;GSTEMX=0;GTIMMX=0
  *-------------------------------------------------------
  * GENTRA GROUP 4: Output controls
  *-------------------------------------------------------
  * Restart-file, history-file and frequency for output
GRSFIL=NONE;GHFILE=GHIS;GOUTFR=1;GSWOUT=10
  * The identifier of the individual history and
    trajectory files
GH1STC=H;GT1STC=NONE
  * The first, last particles and the interval for
    writing history and trajectory files
NGWSTR=1;NGWEND=10;NGWINT=1
  * Particle number for residence-time calculation
GRESTI=0
  *-------------------------------------------------------
  * GENTRA GROUP 5: Provisions for the GENTRA-Earth run
  *-------------------------------------------------------
  * Storage for interphase sources of momentum
STORE(MOMX)
  * Storage for interphase source of heat
STORE(HEAT)
  * Storage for interphase source of mass and for
    vapour mass-fraction
STORE(MASS)
  * Prandtl numbers for vapour mass-fraction
SOLVE(VAPO)
  * Interphase sources of heat and momentum
PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENPAT,U1,FIXFLU,GRND)
COVAL(GENPAT,TEM1,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
RELAX(MOMX,LINRLX,GLNRLX)
RELAX(HEAT,LINRLX,GLNRLX)
RELAX(MASS,LINRLX,GLNRLX)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
  * Activation of GENTRA GROUND-station
L(GENSET
  * ----------------------------------------------------------
  * GENTRA GROUP 7: GENTRA test
  * ----------------------------------------------------------
   * GENTRA auto test  RG(51)--RG(58) are reserved
RG(51)= 6.00E-02;RG(52)= 1.82E+02;RG(53)= 4.98E-03
RG(54)= 7.90E+00;RG(55)= 7.50E+00;RG(56)= 7.50E+00
RG(57)= 5.25E+00;RG(58)= 5.55E-01
 EX(P1  )=5.000E+01;EX(U1)=4.205E+01;EX(SPH1)=1.005E+03
 EX(RHO1)=1.189; EX(TEM1)=4.500E+02