```

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
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
```