PHOTON USE
p
0.20443E+04 0.15633E+04 CR
@
0.23455E+04 0.24296E+04 CR
Geometry outline@
use patgeo
upause 2
@
0.14483E+04 0.17000E+04 CR
Particle trajectories@
use genuse
@
0.14996E+04 0.11268E+04 CR
Press RETURN to end.@
pause
end
ENDUSE
GROUP 1. Run title and other preliminaries
TEXT(Test for tracers in BFC: G704
TITLE
GROUP 6. Body-fitted coordinates or grid distortion
REAL(DIAM,RBEND,ABEND,LPIP1,LPIP2)
INTEGER(NCI,NCJ,NCK1,NCK2,NCK3)
DIAM=0.1;RBEND=5*DIAM;ABEND=3.14159/4.0;LPIP1=1;LPIP2=1
NCI=10;NCJ=10;NCK1=10;NCK2=10;NCK3=10;NCI=5;NCJ=5;NCK1=5
NCK2=10;NCK3=5;REAL(WIN);WIN=0.5;REAL(KEIN);KEIN=(0.1*WIN)**2
REAL(EPIN);EPIN=0.1643*KEIN**1.5/(0.1*DIAM);REAL(RADIUS,ROOT)
RADIUS=DIAM/2;ROOT=RADIUS*0.707;BFC=T
GSET(P,A,-ROOT,-ROOT,0);GSET(P,B,ROOT,-ROOT,0)
GSET(P,C,ROOT,ROOT,0);GSET(P,D,-ROOT,ROOT,0)
GSET(L,AB,A,B,NCI,1.0,ARC,0,-RADIUS,0)
GSET(L,BC,B,C,NCJ,1.0,ARC,RADIUS,0,0)
GSET(L,CD,C,D,NCI,1.0,ARC,0,RADIUS,0)
GSET(L,AD,A,D,NCJ,1.0,ARC,-RADIUS,0,0)
GSET(F,ABCD,A,-,B,-,C,-,D,-)
GSET(D,NCI,NCJ,NCK1+NCK2+NCK3)
GSET(M,ABCD,+I+J,1,1,1,LAP10)
GSET(C,K11,F,K1,1,NCI,1,NCJ,+,0,0,LPIP1,INC,1.0)
GSET(C,K21,F,K11,1,NCI,1,NCJ,RY,ABEND,RBEND,LPIP1,INC,1)
GSET(C,K31,F,K21,1,NCI,1,NCJ,+,LPIP2*0.707,0,LPIP2*0.707,INC,1)
GSET(C,K6,F,K1,1,NCI,1,NCJ,+,0,0,LPIP1,INC,1.0)
GSET(C,K16,F,K6,1,NCI,1,NCJ,RY,ABEND,RBEND,LPIP1,INC,1)
GSET(C,K21,F,K16,1,NCI,1,NCJ,+,LPIP2*0.707,0,LPIP2*0.707,INC,1)
NONORT=T
GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1,W1);SOLUTN(P1,Y,Y,Y,P,P,P);TURMOD(KEMODL);KELIN=1
DIFCUT=0.0
GROUP 9. Properties of the medium (or media)
RHO1=1000.0;ENUL=1.002E-06
GROUP 11. Initialization of variable or porosity fields
INIADD=F
FIINIT(W1)=WIN;FIINIT(P1)=RHO1*WIN/1000.0
FIINIT(KE)=KEIN;FIINIT(EP)=EPIN
RESTRT(ALL);NAMFI=F704
GROUP 13. Boundary conditions and special sources
INLET(GXINLET,LOW,1,NX,1,NY,1,1,1,LSTEP)
VALUE(GXINLET,P1,RHO1*WIN);VALUE(GXINLET,W1,WIN)
COVAL(GXINLET,U1,ONLYMS,0.0);COVAL(GXINLET,V1,ONLYMS,0.0)
COVAL(GXINLET,KE,ONLYMS,KEIN);COVAL(GXINLET,EP,ONLYMS,EPIN)
OUTLET(GXIT,HIGH,1,NX,1,NY,NZ,NZ,1,LSTEP)
COVAL(GXIT,P1,1000.0,0.0);COVAL(GXIT,U1,ONLYMS,0.0)
COVAL(GXIT,V1,ONLYMS,0.0);COVAL(GXIT,W1,ONLYMS,0.0)
COVAL(GXIT,KE,ONLYMS,SAME);COVAL(GXIT,EP,ONLYMS,SAME)
WALL(IX=1,WEST,1,1,1,NY,1,NZ,1,LSTEP)
WALL(IX=NX,EAST,NX,NX,1,NY,1,NZ,1,LSTEP)
WALL(IY=1,SOUTH,1,NY,1,1,1,NZ,1,LSTEP)
WALL(IY=NY,NORTH,1,NX,NY,NY,1,NZ,1,LSTEP)
GROUP 15. Termination of sweeps
LSWEEP=40
GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,0.5);RELAX(U1,FALSDT,2.0/WIN/30)
RELAX(V1,FALSDT,2.0/WIN/30);RELAX(W1,FALSDT,2.0/WIN/30)
RELAX(W1,FALSDT,2.0/WIN/30);RELAX(W1,FALSDT,2.0/WIN/30)
GROUP 24. Dumps for restarts
#GENTRA
*-------------------------------------------------------
* GENTRA GROUP 1: Particle physics
*-------------------------------------------------------
* Particle type - 10
GPTYPE=10
* Gravity components/buoyancy
GGRAX=0;GGRAY=0;GGRAZ=0;GBUOYA=F
* Stochastic model of turbulence
GSTOCH=F
*-------------------------------------------------------
* GENTRA GROUP 2: Boundary conditions for particles
*-------------------------------------------------------
* Inlet-data file-name
GINFIL=Q1
-----------------------------------
* X Y Z
0.01E+00 0.00E+00 0.00E+00
0.02E+00 0.00E+00 0.00E+00
0.03E+00 0.00E+00 0.00E+00
-0.01E+00 0.00E+00 0.00E+00
-0.02E+00 0.00E+00 0.00E+00
-0.03E+00 0.00E+00 0.00E+00
-0.04E+00 0.00E+00 0.00E+00
0.00E+00 0.01E+00 0.00E+00
0.00E+00 0.02E+00 0.00E+00
0.00E+00 0.03E+00 0.00E+00
0.00E+00 0.04E+00 0.00E+00
0.00E+00 -0.01E+00 0.00E+00
0.00E+00 -0.02E+00 0.00E+00
0.00E+00 -0.03E+00 0.00E+00
0.00E+00 -0.04E+00 0.00E+00
------------------------------------
* Wall treatment, and rest coefficient if appropriate
GWALLC=1;GWREST=1.0
* Porosity threshold
GPOROS=0
*-------------------------------------------------------
* GENTRA GROUP 3: Numerical controls
*-------------------------------------------------------
* 1st GENTRA sweep; frequency of calls
GSWEP1=40;GSWEPF=1
* Maximum Lagrangian time-step
GDTMAX=1
* Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=500;GTIMMX=5
*-------------------------------------------------------
* 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
*-------------------------------------------------------
* PIL provisions for the GENTRA-Earth run
*-------------------------------------------------------
* Storage of Cartesian components in BFC cases
STORE(UCRT,VCRT,WCRT)
NCRT=1
* Storage for interphase sources of momentum
* Storage for cell residence-time
STORE(REST)
* 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(58) are reserved
RG(51)= 7.05E+01;RG(52)= 3.92E+03;RG(53)= 2.72E-01
RG(54)= 1.26E+01;RG(55)= 3.26E-02;RG(56)= 3.12E+01
RG(57)= 0.00E+00;RG(58)= 0.00E+00
EX(P1)=3.056E+01;EX(U1)=7.769E-03;EX(V1)=2.797E-03;
EX(W1)=4.979E-01;EX(KE)=1.891E-03;EX(EP)=3.707E-03