PHOTON USE
  p
 
 
 
 
  gr ou y 1
  MSG Velocity vectors
  vec y 1 sh
  msg
  msg Press return to plot pressure contours
  pause
  cont p1 y 1 fil;.01
  msg
  msg Press return to plot streamlines
  pause
  cont cl; red
  stream sh
  y 1
 
  y
  posit
  3
   0.32683E+03 0.18303E+04 CR
   0.34606E+03 0.16414E+04 CR
   0.33965E+03 0.14590E+04 CR
   0.32683E+03 0.12636E+04 CR
   0.32042E+03 0.10878E+04 CR
   0.94845E+03 0.95098E+03  t
  exit
  msg Type e to End
  ENDUSE
 
    GROUP 1. Run title
TEXT(FLOW THROUGH A REACTION TURBINE:    B525
TITLE
  DISPLAY
   This 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=COMPRESS;DRH1DP=COMPRESS;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=120
    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;TSTSWP=-1
SELREF=T; RESFAC=0.01
    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)
    l(30) to activate expert with standard settings
    l(31) to set selref and activate the conjugate-gradient
    solver !! dbs 10.01.09 case 31 no longer exists
#$030
#$031 ! obscure but left, with cngr setting deactivated dbs 21.11.12