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