PHOTON USE
 
  p
 
 
 
   0.20443E+04 0.15633E+04 CR
  @
   0.17623E+04 0.24817E+04  @
  Geometry outline@
  use patgeo
  upause 2
  use genuse
  @
   0.41655E+03 0.11073E+04 CR
  Particle trajectories@
  @
   0.21725E+04 0.59273E+03 CR
  Press e to end.@
  ENDUSE
    GROUP 1. Run title
TEXT(Beams in reaction turbine: G706
TITLE
  DISPLAY
   Particle phase:
   This GENTRA Library Case tests the stubborn-particles option
  in a BFC geometry. The beams are introduced at a central
  point in the computational domain, and radiate in all
  directions.
   The geometry is that of PHOENICS Library
  Case 525 (Flow in a reaction turbine); the number of sweeps has
  been reduced to 10, since, although a converged solution is not
  needed for the gas, it is intended to test the beam option with
  fluid flow.
  Continuous phase:
   The case represents the flow between two turbine blades in a
  cascade.  The blade profile information was supplied in an
  arbitrary coordinate system and the coding here transforms
  this into two surfaces set at the appropriate pitch and stagger
  angle.  The non-orthogonal sheared grid has a portion upstream
  of the blades inclined at the gas inlet angle and a downstream
  portion inclined at the estimated outlet angle. Cyclic boundary
  conditions are used in these regions.
  ENDDIS
      K1       grid location of leading edge
      K2       location of trailing edge
      SC       scaling factor
      TA1      tangent of inlet angle (alpha1) also used for grid.
      TA2      tangent of outlet grid angle.
      SL       sine of angle of blade coordinate system to Z axis
      RH       leading edge radius
      RB       trailing edge radius
      XH,YH    leading edge coordinates in 'arbitrary' system
      XB,YB    trailing edge coordinates
      XX       blade surface ordinate
      YYL,YYU  upper and lower blade surface abscissae
INTEGER(K1,K2,NI,NJ,NK,KB)
REAL(CL,SL,XX,YYU,YYL,ZZC,XXC,RH,RB,YH,YB,XH,XB,PITCH,CHORD,SC)
SL=0.242;CL=(1.-SL*SL)**0.5;RH=1.12;RB=0.8;YH=3.1;YB=7.42;XH=-14.12
XB=17.2;PITCH=29.0;CHORD=33.6;K1=6;K2=23
SC=.01;CHORD=CHORD*SC*CL
REAL(TA1,SA1,TA2,CA2);TA1=.615;TA2=1.4
SA1=(1.0-1.0/(1.0+TA1*TA1))**0.5
REAL(WIN,RHOIN)
REAL(PII,YYC,CRLNK,CRLNK0,PRLNI,XTF,ZTF);INTEGER(III,KKK)
   ** Declaration and setting of the arrays which define the upper
      and lower edges of the turbine blade
ARRAY(YY,REAL,2,16)
YY(1,1)=4.25;YY(1,2)=3.36;YY(1,3)=2.60;YY(1,4)=2.00
YY(1,5)=1.56;YY(1,6)=1.28;YY(1,7)=1.18;YY(1,8)=1.18
YY(1,9)=1.33;YY(1,10)=1.68;YY(1,11)=2.15;YY(1,12)=2.80
YY(1,13)=3.60;YY(1,14)=4.64;YY(1,15)=5.90;YY(1,16)=7.40
YY(2,1)=1.10;YY(2,2)=-1.12;YY(2,3)=-2.78;YY(2,4)=-3.98
YY(2,5)=-4.70;YY(2,6)=-5.05;YY(2,7)=-5.10;YY(2,8)=-4.86
YY(2,9)=-4.39;YY(2,10)=-3.61;YY(2,11)=-2.63;YY(2,12)=-1.46
YY(2,13)=-0.10;YY(2,14)=1.48;YY(2,15)=3.20;YY(2,16)=5.06
    GROUP 3. X-direction grid specification
NX=9
    GROUP 5. Z-direction grid specification
NZ=28
NI=NX+1;NJ=NY+1;NK=NZ+1
    GROUP 6. Body-fitted coordinates or grid distortion
BFC=T;NONORT=T
XCYIZ(1,K1-1,T);XCYIZ(K2+1,NZ,T)
   ** Set origin (profile nose)
SETPT(1,1,K1,-SC*PITCH,0.,0.);SETPT(1,2,K1,-SC*PITCH,1.0,0.)
SETPT(NI,1,K1,0.,0.,0.);SETPT(NI,2,K1,0.,1.0,0.)
   ** Blade surface points
DO II=1,2
+ III=1+NX*(II-1);PII=PITCH*(II-2)
+ DO KK=1,16
+   XX=2*(KK-8)-XH;KKK=K1+KK
+   YYC=YY(II,KK)-YH
+   ZZC=(RH+XX*CL-YYC*SL)*SC;XXC=(XX*SL+YYC*CL+PII)*SC
+   DO JJ=1,2
+     XC(III,JJ,KKK)=XXC;ZC(III,JJ,KKK)=ZZC
+   ENDDO
+ ENDDO
ENDDO
XXC=((XB-XH)*SL+(YB-YH)*CL+RB)*SC
SETPT(NI,1,K2,XXC,0.,CHORD);SETPT(NI,2,K2,XXC,1.0,CHORD)
XXC=XXC-SC*PITCH
SETPT(1,1,K2,XXC,0.,CHORD);SETPT(1,2,K2,XXC,1.0,CHORD)
DOMAIN(1,NI,1,NJ,1,K1)
CRLNK0=CHORD/(K1-1);PRLNI=PITCH*SC/(NI-1)
DO KK=1,K1
+ CRLNK=CRLNK0*(KK-K1)
+ DO II=1,NI
+   XXC=PRLNI*(II-NI)-TA1*CRLNK;ZZC=CRLNK
+   DO JJ=1,NJ
+     XC(II,JJ,KK)=XXC;ZC(II,JJ,KK)=ZZC
+   ENDDO
+ ENDDO
ENDDO
DOMAIN(1,NI,1,NJ,K2,NK)
CRLNK0=CHORD/(NK-K2);PRLNI=PITCH*SC/(NI-1)
ZTF=ZC(1,1,K2);XTF=XC(1,1,K2)
DO KK=K2,NK
+ CRLNK=CRLNK0*(KK-K2)
+ DO II=1,NI
+   XXC=PRLNI*(II-1)+TA2*CRLNK+XTF;ZZC=CRLNK+ZTF
+   DO JJ=1,NJ
+     XC(II,JJ,KK)=XXC;ZC(II,JJ,KK)=ZZC
+   ENDDO
+ ENDDO
ENDDO
DOMAIN(1,NI,1,NJ,K1,K2);MAGIC(T)
    GROUP 7. Variables stored, solved & named
   ** Solve for pressure (whole-field) and velocity.
SOLVE(P1,U1,W1);SOLUTN(P1,Y,Y,Y,N,N,N);STORE(RHO1)
    GROUP 9. Properties of the medium (or media)
RHO1=GRND3;DRH1DP=GRND3;PRESS0=1.0E5;ENUL=1.0E-10
RHO1B=1./1.4;RHO1A=1./PRESS0**RHO1B
    GROUP 11. Initialization of variable or porosity fields
WIN=190.0;RHOIN=1.35;FIINIT(W1)=WIN;FIINIT(RHO1)=RHOIN
    GROUP 13. Boundary conditions and special sources
   ** Inlet
PATCH(INLET,LOW,1,NX,1,NY,1,1,1,1)
COVAL(INLET,P1,FIXFLU,WIN*RHOIN);COVAL(INLET,W1,ONLYMS,WIN)
COVAL(INLET,U1,ONLYMS,-SA1*WIN)
   ** Outlet
PATCH(OUTLET,HIGH,1,NX,1,NY,NZ,NZ,1,1);COVAL(OUTLET,P1,1000.,0.0)
COVAL(OUTLET,U1,ONLYMS,0.0);COVAL(OUTLET,W1,ONLYMS,0.0)
    GROUP 15. Termination of sweeps
LSWEEP=10
    GROUP 17. Under-relaxation devices
RELAX(P1,LINRLX,.3);RELAX(RHO1,LINRLX,.3)
RELAX(U1,FALSDT,0.05*SC);RELAX(W1,FALSDT,0.05*SC)
OUTPUT(DEN1,N,N,N,N,N,N)
    GROUP 22. Spot-value print-out
ITABL=3;IXMON=3;IZMON=10
    GROUP 23. Field print-out and plot control
NXPRIN=3;IZPRF=4;IZPRL=24
PATCH(XZ,CONTUR,1,NX,1,1,1,NZ,1,1)
PLOT(XZ,P1,0.0,20.0);PLOT(XZ,W1,0.0,20.0)
PATCH(PRESSURE,PROFIL,NX,NX,1,1,1,NZ,1,1)
PLOT(PRESSURE,P1,0.0,0.0)
PATCH(SUCTION,PROFIL,1,1,1,1,1,NZ,1,1)
PLOT(SUCTION,P1,0.0,0.0)
#GENTRA
  *-------------------------------------------------------
  * GENTRA GROUP 1: Particle physics
  *-------------------------------------------------------
  * Particle type - 20 - STUBBORN_PARTICLES
GPTYPE=20
  * Gravity components in GENTRA Cartesian system/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
  * Use this space for inlet-data table if GINFIL=Q1
  
  *__XP__ZP__UP__WP
  -0.145 0.165 0.0 1.0
  -0.145 0.165 0.0 -1.0
  -0.145 0.165 1.0 0.0
  -0.145 0.165 -1.0 0.0
  -0.145 0.165 1.0 1.0
  -0.145 0.165 1.0 -1.0
  -0.145 0.165 -1.0 1.0
  -0.145 0.165 -1.0 -1.0
  -0.145 0.165 2.0 1.0
  -0.145 0.165 -2.0 1.0
  -0.145 0.165 2.0 -1.0
  -0.145 0.165 -2.0 -1.0
  -0.145 0.165 1.0 2.0
  -0.145 0.165 1.0 -2.0
  -0.145 0.165 -1.0 2.0
  -0.145 0.165 -1.0 -2.0
  -0.145 0.165 3.0 1.0
  -0.145 0.165 -3.0 1.0
  -0.145 0.165 3.0 -1.0
  -0.145 0.165 -3.0 -1.0
  -0.145 0.165 1.0 3.0
  -0.145 0.165 1.0 -3.0
  -0.145 0.165 -1.0 3.0
  -0.145 0.165 -1.0 -3.0
  
  * Wall treatment, and rest coefficient if appropriate
GWALLC=1
  * Porosity threshold
GPOROS=0
  *-------------------------------------------------------
  * GENTRA GROUP 3: Numerical controls
  *-------------------------------------------------------
  * 1st GENTRA sweep; frequency of calls
GSWEP1=10;GSWEPF=1
  * Maximum Lagrangian time-step
GDTMAX=1
  * Min # of t-steps per cell; max # of t-steps; timeout
GLAGTS=5;GSTEMX=1000;GTIMMX=0
  *-------------------------------------------------------
  * GENTRA GROUP 4: Output controls
  *-------------------------------------------------------
  * Trajectory-file identifier, history-file identifier
    and frequency for output
  * Particle number for residence-time calculation
GT1STC=T;GSWOUT=10;NGWSTR=1;NGWEND=20;NGWINT=1;GRESTI=0
  *-------------------------------------------------------
  * GENTRA GROUP 5: 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
  * 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)= 2.92E+00;RG(52)= 8.44E+02;RG(53)= 7.61E-02
RG(54)=-2.82E+00;RG(55)= 1.20E+01;RG(56)= 3.36E+00
RG(57)= 2.40E+17;RG(58)= 2.40E+23
EX(P1)=1.354E+04;EX(U1)=6.839E+01;EX(W1)=1.970E+02;EX(RHO1)=1.006