PHOTON USE
p
0.20443E+04 0.15633E+04 CR
use patgeo
vec x 1
@
0.86514E+03 0.84676E+03 CR
Velocity vectors of the continuous phase.
Press RETURN to continue.
pause
red
use genuse
@
0.72415E+03 0.89887E+03 CR
Particles bounced on the wall and
the central blockage.
Press e to exit.
ENDUSE
GROUP 1. Run title
TEXT(Particles in 2D curved duct: G209
TITLE
DISPLAY
This case deals with conjugate heat transfer between a metal
solid which stands in the middle of a curved 2-D duct, and a
cold air stream which flows through the duct,passing the
'obstacle'. The geometry is represented by means of body-fitted
coordinates
ENDDIS
GROUP 2. Transience; time-step specification
GROUP 3. X-direction grid specification
GROUP 4. Y-direction grid specification
NY=10;YVLAST=.05
GROUP 5. Z-direction grid specification
NZ=10;ZWLAST=.20
GROUP 6. Body-fitted coordinates or grid distortion
BFC=T
* Set points
XPO= 0.0000E+00;YPO= 0.0000E+00;ZPO= 0.0000E+00;GSET(P,P1 )
XPO= 0.0000E+00;YPO= 1.3500E-02;ZPO= 5.0000E-02;GSET(P,P2 )
XPO= 0.0000E+00;YPO= 0.0000E+00;ZPO= 1.0000E-01;GSET(P,P3 )
XPO= 0.0000E+00;YPO= 4.3000E-02;ZPO=-2.5000E-02;GSET(P,P4 )
XPO= 0.0000E+00;YPO= 6.3500E-02;ZPO= 5.0000E-02;GSET(P,P5 )
XPO= 0.0000E+00;YPO= 4.3000E-02;ZPO= 1.2500E-01;GSET(P,P6 )
* Set lines/arcs
GSET(L,A1,P1,P3,10,1.0,ARC,P2)
GSET(L,A2,P4,P6,10,1.0,ARC,P5)
GSET(L,L3,P1,P4,10,1.0)
GSET(L,L4,P3,P6,10,1.0)
* Set frames
GSET(F,F1,P1,-,P3,-,P6,-,P4,-)
* Match a grid mesh
GSET(M,F1,+K+J,1,1,1,TRANS)
* Copy/Transfer/Block grid planes
GSET(C,I2,F,I1,1,10,1,10,+,2.0000E-01,0,0,INC,1)
GROUP 7. Variables stored, solved & named
SOLVE(P1,V1,W1,TEM1)
SOLUTN(P1,P,P,Y,P,P,P)
SOLUTN(TEM1,P,P,Y,P,P,P)
STORE(PRPS,BLOK,DEN1,VISL,VCRT,WCRT)
GROUP 8. Terms (in differential equations) & devices
IVARBK=50;ISOLBK=1
TERMS(TEM1,N,Y,Y,Y,P,P)
GROUP 9. Properties of the medium (or media)
SETPRPS(1,0)
GROUP 10. Inter-phase-transfer processes and properties
GROUP 11. Initialization of variable or porosity fields
FIINIT(W1)=1.0;FIINIT(BLOK)=1.0
CONPOR(SOLID,-1,CELL,1,1,-5,-6,-5,-6)
INIT(SOLID,PRPS,0.0,103)
INIT(SOLID,BLOK,0.0,2.0)
INIT(SOLID,TEM1,0.0,5.0)
GROUP 12. Unused
GROUP 13. Boundary conditions and special sources
INLET(GXINL,LOW,1,NX,1,NY,1,1,1,1)
VALUE(GXINL,P1,1.29);VALUE(GXINL,W1,1.00)
VALUE(GXINL,TEM1,5.00)
PATCH(HOTSOLID,CELL,1,NX,5,6,5,6,1,1)
COVAL(HOTSOLID,TEM1,FIXFLU,10.0)
PATCH(GXOUT,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(GXOUT,P1,1.0E3,0.0)
WALL(WALL1,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(WALL1,W1,1.0,0.0)
WALL(WALL2,SOUTH,1,NX,1,1,1,NZ,1,1)
COVAL(WALL2,W1,1.0,0.0)
GROUP 14. Downstream pressure for PARAB=.TRUE.
GROUP 15. Termination of sweeps
LSWEEP=500
RELAX(P1,LINRLX,0.2);RELAX(V1,FALSDT,0.001)
RELAX(W1,FALSDT,0.001);RELAX(TEM1,FALSDT,1E8)
GROUP 16. Termination of iterations
GROUP 17. Under-relaxation devices
GROUP 18. Limits on variables or increments to them
GROUP 19. Data communicated by satellite to GROUND
GROUP 20. Preliminary print-out
GROUP 21. Print-out of variables
GROUP 22. Spot-value print-out
IYMON=8;IZMON=7;TSTSWP=-1
GROUP 23. Field print-out and plot control
NYPRIN=0;NZPRIN=0
#GENTRA
*-------------------------------------------------------
* GENTRA GROUP 1: Particle physics
*-------------------------------------------------------
* Particle type - 30 - ISOTHERMAL_PARTICLES
GPTYPE=30
* Gravity components in GENTRA Cartesian system/buoyancy
GGRAX=0;GGRAY=-9.8;GGRAZ=0;GBUOYA=F
* Stochastic model of turbulence
GSTOCH=F
* Data for isothermal particles
GDRAG=GRND1
*-------------------------------------------------------
* GENTRA GROUP 2: Boundary conditions for particles
*-------------------------------------------------------
* Inlet-data file-name
GINFIL=Q1
* Use this space for inlet-data table if GINFIL=Q1
*__YP__ZP__VP__WP__DIAM_LIQDEN_FRATE_(NUMB)
0.0001 0.0001 10. 2.0 0.0001 1000.0 0.001
0.0001 0.0001 8.0 2.0 0.0001 1000.0 0.001
0.0001 0.0001 7.0 2.0 0.0001 1000.0 0.001
0.0001 0.0001 6.0 2.0 0.0001 1000.0 0.001
0.0001 0.0001 5.0 2.0 0.0001 1000.0 0.001
0.0001 0.0001 4.0 2.0 0.0001 1000.0 0.001
0.0001 0.0001 3.0 2.0 0.0001 1000.0 0.001
0.0001 0.0001 2.7 2.0 0.0001 1000.0 0.001
0.0001 0.0001 2.5 2.0 0.0001 1000.0 0.001
0.0001 0.0001 2.2 2.0 0.0001 1000.0 0.001
0.0001 0.0001 2.0 2.0 0.0001 1000.0 0.001
0.0001 0.0001 1.8 2.0 0.0001 1000.0 0.001
* 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=20
* Maximum Lagrangian time-step; time step size multplier
GDTMAX=0.001;GRTFRC=1.000
* Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=1000;GTIMMX=5
*-------------------------------------------------------
* GENTRA GROUP 4: Output controls
*-------------------------------------------------------
* Restart-file, history-file 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 of Cartesian components in BFC cases
NCRT=1
* Storage for interphase sources of momentum
STORE(MOMY)
STORE(MOMZ)
* Interphase sources of heat and momentum
PATCH(GENPAT,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GENPAT,V1,FIXFLU,GRND)
COVAL(GENPAT,W1,FIXFLU,GRND)
* Relaxation of interphase sources
GLNRLX=0.5
RELAX(MOMY,LINRLX,GLNRLX)
RELAX(MOMZ,LINRLX,GLNRLX)
* Activation of GENTRA GROUND-station
L(GENSET
* ----------------------------------------------------------
* GENTRA GROUP 7: GENTRA test
* ----------------------------------------------------------
* GENTRA auto test RG(51)--RG(56) are reserved
RG(51)=1.433E+00;RG(52)=1.895E+03;RG(53)=9.163E-03
RG(54)=6.000E+00;RG(55)=2.523E-01;RG(56)=1.346E+00
RG(57)=1.200E+17;RG(58)=1.200E+00