AUTOPLOT USE
file
ghis 2
data
time diam
pl
red
@
0.97408E+03 0.14786E+04 CR
Particle diameter vers time@
@
0.31273E+04 0.16935E+03 CR
Time (sec.)@
@
0.12176E+03 0.28855E+04 CR
Diameter@
@
0.19161E+04 0.84676E+03 CR
Press RETURN to continue@
pause
cl
file
phi 5
data
vapo
pl
red
@
0.09638E+04 0.15613E+04 CR
Vapour concentration along particle
trajectory.
Press e to END.@
ENDUSE
GROUP 1. Run title and other preliminaries
TEXT(PARTICLE EVAPORATING IN PIPE: G502
TITLE
DISPLAY
Particle evaporates in a pipe with constant gas
temperature and particle velocity, but changing
vapour concentraton.
Particle temperature is also kept constant.
ENDDIS
GROUP 2. Transience; time-step specification
** Activate unsteady term in the solved equations
STEADY=T;GRDPWR(T,20,0.04,1.0)
GROUP 3. X-direction grid specification
CARTES=T
GROUP 5. Z-direction grid specification
RSET(D,CHAM,2.0,2.0,10.0,1,0)
RSET(M,1,1,10)
GROUP 7. Variables stored, solved & named
** Solve for three extra variables, A, B AND C, treated as
temperatures with different Prandtl numbers.
RESTRT(ALL)
NAME(C6)=VAPO;NAME(50)=REST
SOLVE(P1);STORE(H1)
SOLVE(VAPO,W1)
STORE(DEN1);PRESS0=1.E+05
STORE(MOMZ,HEAT,MASS)
TERMS(H1,N,Y,Y,Y,Y,N)
GROUP 9. Properties of the medium (or media)
ENUL=1.E-5;PRNDTL(H1)=0.7
TMP1A=0.0;TMP1B=1.0
DRH1DP=GRND3;RHO1A=0.0;RHO1B=1.0
RHO1=GRND
GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=10.0
FIINIT(P1)=0.0;FIINIT(H1)=400.0;FIINIT(DEN1)=3.125E-02
FIINIT(VAPO)=0.0
VARMIN(VAPO)=0.0;VARMAX(VAPO)=1.0
GROUP 13. Boundary conditions and special sources
1. Wall at north boundary
WALL(PIPEWALL,NORTH,1,NX,NY,NY,1,NZ,1,LSTEP)
COVAL(PIPEWALL,W1,1.0,0.0)
COVAL(PIPEWALL,VAPO,0.0,0.0)
COVAL(PIPEWALL,H1,1./PRNDTL(H1),7000.0)
2. Inlet-- uniform flow
PATCH(UNIFORM,LOW,1,NX,1,NY,1,1,1,LSTEP)
COVAL(UNIFORM,P1,FIXFLU,1.715E-01*10.0)
COVAL(UNIFORM,W1,ONLYMS,10.0)
COVAL(UNIFORM,VAPO,ONLYMS,0.0)
COVAL(UNIFORM,H1,ONLYMS,7000.0)
2. Outlet-- fixed pressure
PATCH(GXOUT,HIGH,1,NX,1,NY,NZ,NZ,1,LSTEP)
COVAL(GXOUT,P1,FIXVAL,0.0);COVAL(GXOUT,U1,ONLYMS,0.0)
COVAL(GXOUT,H1,ONLYMS,SAME)
COVAL(GXOUT,VAPO,ONLYMS,SAME)
GROUP 15. Termination of sweeps
LSWEEP=20
RESREF(P1)=1.E-16
RESREF(H1)=1.E-16;RESREF(W1)=1.E-16
RESREF(VAPO)=1.E-16
GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,0.7);RELAX(W1,FALSDT,0.1)
RELAX(VAPO,LINRLX,0.1)
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=1;NZPRIN=1
NOWIPE=T
#GENTRA
*-------------------------------------------------------
* GENTRA GROUP 1: Particle physics
*-------------------------------------------------------
* Particle type - 22 - STUBBON PARTICLE
GPTYPE=22
* 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;GNUSS=GRND1
GKONV=GRND1;GCPCON=1.0;GCPVAP=1.0
GMWCON=10.0;GMWVAP=50.0;GLIQST=GRND1
GCPLIQ=GRND1;GHLIQD=GRND1;GVAPST=GRND1
GLATVP=100000;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
*_ZP_WP_DI_LDEN_MDOT_TEM_(NUM)_(STR)_(END)
0.0 100.0 0.05 1.0 0.1 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=1;GSWEPF=1
* Maximum Lagrangian time-step; time step size multplier
GDTMAX=1.e-04;GRTFRC=0.1
* Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=2000;GTIMMX=100
*-------------------------------------------------------
* 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=NONE;GT1STC=NONE
* The first, last particles and the interval for
writing history and trajectory files
NGWSTR=1;NGWEND=20;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
* 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,W1,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=0.1
RELAX(MOMZ,LINRLX,GLNRLX)
RELAX(HEAT,LINRLX,GLNRLX)
RELAX(MASS,LINRLX,GLNRLX)
* Activation of GENTRA GROUND-station
L(GENSET
* ----------------------------------------------------------
* GENTRA GROUP 7: GENTRA test
* ----------------------------------------------------------
* GENTRA auto test RG(51)--RG(58) are reserved
RG(51)= 9.89E-02;RG(52)= 1.06E+03;RG(53)= 9.38E-05
RG(54)= 1.00E+00;RG(55)= 1.00E+00;RG(56)= 9.89E+00
RG(57)= 3.00E-01;RG(58)= 0.00E+00
EX(P1 )= 2.487E+02; EX(W1 )= 5.753; EX(H1)= 4.0E+02
EX(MOMZ)= 3.466E+02; EX(DEN1)= 3.004E-01
LIBREF= 502