Turbine cascade

Download this example

Input file

# Turbine cascade
#
# Run a stationary turbine cascade with prescribed flow angles.
#


# Use a new workdir for each run
workdir: runs/*

# Perfect gas inlet condition
inlet:
  Po: 1e5
  To: 300.
  cp: 1005.
  gamma: 1.4

# Set viscosity from Reynolds number
Re_surf: 4e5

# Turbine cascade design variables
mean_line:
  type: turbine_cascade
  span: [0.01, 0.011]
  Alpha: [40., -65.0]
  Ma2: 0.6
  Yh: 0.029
  htr: 0.99

# Setup annulus using fixed aspect ratios
annulus:
  AR_gap: [1.0, 1.0]
  AR_chord: 2.0

# Specify camber and thickness one one mid-span section
blades:
    - spf: 0.5
      q_thick: [0.05, 0.12, 0.3, 0., 0.00, 0.18]
      q_camber: [10., -2., 0.0]

# Discretise fluid domain using a coarse H-mesh
mesh:
  type: h
  yplus: 30.0
  resolution_factor: 0.5

# Use circulation for number of blades
nblade:
  - Co: 0.7

# Configure the built-in flow solver
solver:
  type: ember
  n_step: 500
  n_step_avg: 100

Log output

*** TURBIGEN v2.3.0 ***
Starting at 2025-05-09T16:45:55
Working directory: /builds/jb753/turbigen/runs/0001
Inlet: PerfectState(P=1.000 bar, T=300.0 K)
MeanLine(
    Po=[1.     0.9927] bar,
    To=[300. 300.] K,
    Ma=[0.311 0.6  ],
    Vx=[81.9 85.1] m/s,
    Vr=[0. 0.] m/s,
    Vt=[  68.7 -182.4] m/s,
    Vt_rel=[  68.7 -182.4] m/s,
    Al=[ 40. -65.] deg,
    Al_rel=[ 40. -65.] deg,
    rpm=[0. 0.],
    mdot=[5.66 5.66] kg/s
    )
Designing annulus...
FixedAR(nrow=1, x=[0.002625], r=[0.995], AR=[2.])
Designing blades...
Nblade: [473]
Tip gaps: [0.]
Re_surf=[4e+05]
Generating mesh...
Making an H-mesh...
ncell/1e6=0.1
Applying 2D guess...
Initialising ember...
Patitioning onto 1 processors...
Starting the main time-stepping loop...
499: tpnps=3.106e-07, remaining=0m0s
  block 0: 6.82e-05 2.46e-02 5.80e-04 2.79e-02 5.98e+00
Elapsed time 16.72s
Average tpnps=3.106e-07
mdot_in/out=5.41/5.61, err=-3.5%
Post-processing...
Variable  Nominal    Actual  Err_abs  Err_rel/%
-----------------------------------------------
Alpha[0]       40      40.1    -0.15      -0.37
Alpha[1]      -65     -63.3     -1.7        2.5
     Ma2      0.6      0.55   0.0496       8.27
      Yh    0.029     0.068   -0.039       -134
     htr     0.99      0.99        0  -1.56e-06
 span[0]     0.01      0.01        0          0
 span[1]    0.011     0.011        0          0
      Ys             0.0681
      RR                  1
 Beta[0]              0.101
 Beta[1]           -0.00785
Efficiency/%: eta_tt=8.9, eta_ts=0.7

Plots