AUTOPLOT USE
 
  file
  h00001 2
  h00002
  h00003
  h00004
  h00005
 
  data 1
  time sfrc
  data 2
  time sfrc
  data 3
  time sfrc
  data 44
  4
  time sfrc
  data 5
  time sfrc
  pl
  red
  @
   0.11856E+04 0.23840E+04 CR
  Particle solid fraction vers time for
  five particles introduced each 0.002
  second.
 
 
 
 
 
                 Press e to END.@
  ENDUSE
TEXT(SOLIDIFYING,1D,TRANSIENT,M=3.0:G401
TITLE
  DISPLAY
  Constant gas temperature, no slip.
  Constant properties.
  The particle temperature as a function of time within
  the solidus-liquidus temperature range is given by:
 
  (Cp + 3.L(Tl-Tg)**2/DTLS**3).log((Tg-Tp)/(Tg-T0)) +
  3L/DTLS**3(2Tl-Tg-0.5(Tp+T0))(T0-Tp) = alpha.t/mass
 
  where: alpha - particle heat transfer coefficient
         DTLS  - solidus-liquidus temperature range
         L     - latent heat of solidification
         mass  - particle mass
         t     - time (time=0.0 when Tp=Tl)
         Tg    - continuous-phase temperature
         Tl    - liquidus temperature
         Tp    - particle temperature
         T0    - particle initial temperature
  ENDDIS
REAL(WIN);WIN=100.0
STEADY=F;GRDPWR(T,5,1.0E-02,1.0)
GRDPWR(X,5,1.0,1.0)
SOLVE(P1,U1);STORE(TMP1)
TMP1=GRND1;TMP1A=500.0
FIINIT(U1)=WIN;FIINIT(P1)=RHO1*WIN;FIINIT(TMP1)=500.0
PATCH(GXIN,WEST,1,1,1,1,1,1,1,LSTEP)
COVAL(GXIN,P1,FIXFLU,RHO1*WIN);COVAL(GXIN,U1,ONLYMS,WIN)
PATCH(GXOUT,EAST,NX,NX,1,1,1,1,1,LSTEP)
COVAL(GXOUT,P1,FIXP,0.0)
LSWEEP=1
TRACE=F
#GENTRA
  *-------------------------------------------------------
  * GENTRA GROUP 1: Particle physics
  *-------------------------------------------------------
  * Particle type - 50 - MELTING/SOLIDIFYING_PARTICLES
GPTYPE=50
  * Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0;GGRAY=0;GGRAZ=0;GBUOYA=F
  * Stochastic model of turbulence
GSTOCH=F
  * Data for melting/solidifying particles
GDRAG=GRND1;GKONC=3.3400E-02;GNUSS=1.9060E+01
GLHEAS=2.5000E+05;GSOLIN=3;GSOLST=1450
GLIQST=1600;GCPLIQ=1000;GHLIQD=GRND1
GCPSOL=1000
  *-------------------------------------------------------
  * 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_SDEN_(NUM)_(STR)_(END)
   0.1 100 5e-5 8000 1e-6 1800 8000
  
  * 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=50;GSTEMX=0;GTIMMX=0
  *-------------------------------------------------------
  * 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=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)
  * Interphase sources of heat and momentum
PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENPAT,U1,FIXFLU,GRND)
  * Relaxation of interphase sources
GLNRLX=1
RELAX(MOMX,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)= 2.90E-02;RG(52)= 7.38E+02;RG(53)= 5.90E-04
RG(54)= 8.40E+00;RG(55)= 7.50E+00;RG(56)= 7.50E+00
RG(57)= 1.74E+01;RG(58)= 7.00E-01
 EX(P1)=1.000E+02;EX(U1)=1.000E+02;EX(TMP1)=5.0E+02