AUTOPLOT USE
  file
  ghis 2
 
  data
  time tpar
  pl
  red
  @
   0.27492E+04 0.19541E+03 CR
  Time (sec.)@
  @
   0.17303E+03 0.24230E+04  T
  T(part.)@
  @
   0.17175E+04 0.17261E+04 CR
  Particle temperature vers time.@
  @
   0.25634E+04 0.10226E+04  P
  Press e to END
  @
  ENDUSE
    GROUP 1. Run title and other preliminaries
TEXT(Particle heating in pipe : G301
TITLE
  DISPLAY
  The problem describes transient particle heating in a pipe
  with constant gas temperature
  The result is compared against analytical solution:
   PART. TEMPERATURE=370.0-70.0*EXP(-TIME*0.8722375)
  ENDDIS
    GROUP 2. Transience; time-step specification
  ** Activate unsteady term in the solved equations
STEADY=F
GRDPWR(T,50,10.,1.0)
    GROUP 3. X-direction grid specification
CARTES=F
    GROUP 4. Y-direction grid specification
GRDPWR(Y,2,1.0,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,10,50.,1.0)
RSET(X,1,4,1)
RSET(D,CHAM,6.2832E+00,1.0,50,1,0)
RSET(M,4,2,10)
    GROUP 7. Variables stored, solved & named
  ** Solve for three extra variables, A, B AND C, treated as
     temperatures with different Prandtl numbers.
STORE(P1,U1,V1,W1,H1)
    GROUP 9. Properties of the medium (or media)
ENUL=1.E-5;PRNDTL(H1)=0.7
TMP1A=0.0;TMP1B=1.0
    GROUP 11. Initialization of variable or porosity fields
FIINIT(U1)=0.0;FIINIT(V1)=0.0;FIINIT(W1)=0.0;FIINIT(P1)=0.0
FIINIT(H1)=370.0
    GROUP 13. Boundary conditions and special sources
    1. Wall at north boundary
WALL(PIPEWALL,NORTH,1,1,NY,NY,1,NZ,1,LSTEP)
COVAL(PIPEWALL,W1,1.0,0.0);COVAL(PIPEWALL,H1,1./PRNDTL(H1),1.0)
    2. Inlet-- uniform flow
INLET(UNIFORM,LOW,1,NX,1,NY,1,1,1,LSTEP)
VALUE(UNIFORM,P1,0.1);VALUE(UNIFORM,W1,0.1)
VALUE(UNIFORM,V1,0.0);VALUE(UNIFORM,H1,1.0)
    2. Outlet-- fixed pressure
PATCH(GXOUT,HIGH,1,1,1,NY,NZ,NZ,1,LSTEP)
COVAL(GXOUT,P1,FIXVAL,0.0)
COVAL(GXOUT,V1,ONLYMS,0.0);COVAL(GXOUT,W1,ONLYMS,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=2
RESREF(P1)=1.E-6;RESREF(V1)=1.E-6
RESREF(H1)=1.E-6;RESREF(W1)=1.E-6
    GROUP 17. Under-relaxation devices
RELAX(H1,LINRLX,0.0)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 22. Spot-value print-out
NPRMON=LSWEEP;IYMON=5
    GROUP 22. Spot-value print-out
IYMON=NY/2;IZMON=NZ/2;IPLTL=30
    GROUP 23. Field print-out and plot control
NTPRIN=LSTEP;NZPRIN=2
NOWIPE=T
TRACE=F
#GENTRA
  *-------------------------------------------------------
  * GENTRA GROUP 1: Particle physics
  *-------------------------------------------------------
  * Particle type - 21 - STUBBORN_PARTICLES with heat transfer
GPTYPE=21
  * Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0;GGRAY=0;GGRAZ=0;GBUOYA=F
  * Stochastic model of turbulence
GSTOCH=F
  * Data for particles with heat transfer
GDRAG=GRND1;GKONC=2.6300E-02;GNUSS=GRND1
GCPCON=1;GCPLIQ=1000.
  *-------------------------------------------------------
  * GENTRA GROUP 2: Boundary conditions for particles
  *-------------------------------------------------------
  * The times at which injection starts and ends
  * Inlet-data file-name
GINFIL=Q1
  * Use this space for inlet-data table if GINFIL=Q1
  
  *__XP__YP__ZP__UP__VP__WP__DIAM_LIQDEN_FRATE_TEMP_(NUMB)
   0.1 0.05 0.0 0.0 0.0 1.0 0.1 1.0 1.0 300 0.0 1.e-03
  
  * 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=0.1;GRTFRC=1.0
  * Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=100;GTIMMX=7
  *-------------------------------------------------------
  * 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
  * Particle number for residence-time calculation
GRESTI=0
  *-------------------------------------------------------
  * GENTRA GROUP 5: Provisions for the GENTRA-Earth run
  *-------------------------------------------------------
  * Storage for interphase sources of momentum
  * Interphase sources of heat and momentum
  * Relaxation of interphase sources
  * Activation of GENTRA GROUND-station
L(GENSET
  * ----------------------------------------------------------
  * GENTRA GROUP 7: GENTRA test
  * ----------------------------------------------------------
   * GENTRA auto test  RG(51)--RG(56) are reserved
RG(51)= 1.00E+01;RG(52)= 1.20E+02;RG(53)= 4.33E+00
RG(54)= 5.25E+01;RG(55)= 5.02E+01;RG(56)= 2.55E+02
RG(57)= 1.81E+01;RG(58)= 5.00E+03
 EX(H1  )=   3.700E+02