Next Article in Journal
Prediction of Transient Pressure Fluctuations within a Low-Pressure Turbine Cascade Using a Lanczos-Filtered Harmonic Balance Method
Previous Article in Journal
Modelling the Condensation Phenomena within the Radial Turbine of a Fuel Cell Turbocharger
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pulsed Flow Turbine Design Recommendations †

by
Florian Hermet
*,
Nicolas Binder
,
Jérémie Gressier
and
Gonzalo Sáez-Mischlich
ISAE-SUPAERO, Université de Toulouse, 31400 Toulouse, France
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in the Proceedings of the 14th European Turbomachinery Conference, Gdansk, Poland, 12–16 April 2021.
Int. J. Turbomach. Propuls. Power 2021, 6(3), 24; https://0-doi-org.brum.beds.ac.uk/10.3390/ijtpp6030024
Submission received: 18 May 2021 / Revised: 1 July 2021 / Accepted: 6 July 2021 / Published: 8 July 2021

Abstract

:
A preliminary analysis of turbine design, fit for pulsed flow, is proposed in this paper. It focuses on an academic 2D configuration using inviscid flows, since pressure loads due to wave propagation are several orders of magnitude higher than friction and viscous effects do not significantly impinge on the inviscid part, as previously shown by Hermet, 2021. As such, a large parametric study was carried out using the design of experiments methodology. A performance indicator adapted to unsteady environment is carefully defined before detailing the factors chosen for the design of experiments. Since the number of factors is substantial, a screening design to identify the factors influence on the output is first established. The non-influential factors are then omitted in a more quantitative study of the output law. The surface response calculation allows determining the factor level favouring the best output. Consequently, the main trends in the turbine design driven by a pulsed flow can be stated.

1. Introduction

In many applications, the turbine is subject to temporal variations of its inlet conditions. The most extreme case is found when the turbine is involved in a thermodynamic cycle including the isochoric combustion process. During the last decade, the rising interest in aeronautics for isochoric combustion cycles has stimulated research on axial turbines supplied by severe unsteadiness. However, major contributions to this question are credited to the automotive turbocharger community. In Baines [1], a detailed summary of nearly 20 years of research on pulsed flow in radial turbines is given. This review highlights the progress and contradictions in the scientific community, particularly regarding the pulsed flow influence on turbine performance. The complexity of this problem comes from the coupling between different phenomenon involving time scales that are not easily separable (pulse time scale, propagative time scale, advective time scale, etc.). Moreover, the definition of a turbine performance indicator is made complex by this severe flow unsteadiness. This partly explains why there is no consensus on the pulsed flow influence in turbines, and more specifically whether or not the performance could benefit from the unsteadiness. The necessity to adapt the design process to the unsteady environment is increasing in the literature [2].
The complexity of the geometries involved in the different studies of the literature is also an additional difficulty to enlighten the comprehension. Recent numerical work [3] focused on that question through a simplified cascade approach. It has been shown that instantaneous loading differs from the quasi-steady approach in response to a rapid increase in turbine inlet pressure. The relevant time scale of the inlet perturbation associated with the additional work extracted by the turbine was also identified by the authors. In short, it appears that the curvature intensifies the different pressure effects during the wave propagation. Thus, a compression wave propagating in the direction of the flow overloads the blades, compared with a quasi-steady transformation. Other combinations of the nature of waves (compression or expansion) and the direction (upstream or downstream) have different consequences. For example, a streamwise expansion underloads the blades. In pulsating flows, a succession of waves of different natures and directions of propagation occurs. The net benefit of unsteady effects over a complete cycle thus needs to be quantified. It is one of the targets of this paper. However, the main objective is to give a clear statement on which geometrical parameters are likely to promote unsteady performance.
The present work is thus a direct continuation of the work of Hermet et al. [3]. The physics of pulsed flows is terribly complex. Therefore, the analysis is conducted by means of parametric studies for which the numerical design of experiments (DOE) is built. The adequate performance indicator is firstly discussed before describing the selected factors. Because of the large number of factors involved, an initial screening phase is conducted. This allows identifying the influence of each factor on the response. Noninfluential factors are then discarded for a more quantitative study of the output law. The quantitative prediction is built by calculating the response surface. As a conclusion, the main trends observed in 2D are exploited to propose first recommendations for the design of turbines able to benefit from the large unsteadiness of a pulsating flow.

2. Methodology

The parametric study focused on the simulation of stabilized pulsed flow through simplified cascade approaches of the stator and rotor with an in-house IC 3 solver, forked from the CharLES X solver [4]. IC 3 is based on the resolution of the compressible formulation of the Navier–Stokes equations in their conservative form, spatially filtered, on an unstructured mesh using a finite volume method. An explicit third-order Runge–Kutta scheme was used for time advancement while an essentially non-oscillatory (ENO) second-order shock-capturing scheme was applied to compute the flux. The nonconformal and sliding interface, which develops between the rotor and stator domains due to the relative movement of the rotor and stator parts, was treated using a sliding mesh method based on the approach proposed by Bohbot et al. [5]. Such a methodology relies on the computation of geometrical intersections between the faces/edges located at the nonconformal interface and requires the reassessment of the parallel rank connectivity involved in the computation. The proposed sliding mesh method enables flux conservation across the interface and ensures accurate interpolations from cell to face averages by using information from cells and faces at each side of the interface. Moreover, the periodic properties of the solution were taken into account by duplicating the sliding interface geometry and solution data for each Runge–Kutta substep. Additionally, the grid motion was taken into account by means of the Arbitrary Lagrangian-Eulerian (ALE) [6] formulation, which allows solving conservation laws with moving grids in the absolute frame by adding specific flux terms to the discretized Navier–Stokes equations. The common face fluxes were computed following the Riemann HLLC solver adapted for moving grids [7]. Regarding the mesh grid density, a mesh convergence was performed on the first simulation of Table 1 to determine a reference density grid. For the other configurations, the number of points in the axial and normal direction was adapted in order to obtain the same density grid as the reference ones. Meshes were automatically generated by ICEM-CFD [8] using a script written in the Tool Command Language (TCL).
A stabilized pulsed flow is a perfectly periodic flow, similar to flows caused by the cyclic opening and closing of a valve separating the combustion chamber from the turbine. The choice of the performance criteria selected for the response of the DOE is first discussed. The selected factors of the design of experiments are then presented.

2.1. Response of the Design of Experiments

Most of the performance indicators used in the turbocharger literature do not take into account the full complexity of flow unsteadiness, as shown in [9]. In order to design a turbine under pulsed flow, the performance indicators of the turbine must be clearly defined. For this, it is appropriate to mention the mass conservation equation, as well as the first and second thermodynamic principles of an adiabatic unsteady system ( Q ˙ = 0 ) , Equation (1).
m ˙ i n ( t ) m ˙ o u t ( t ) = d d t V ρ ( t ) d V m ˙ i n ( t ) h t , i n ( t ) m ˙ o u t ( t ) h t , o u t ( t ) = d d t V ρ ( t ) e t ( t ) d V + W ˙ m ˙ i n ( t ) s i n ( t ) m ˙ o u t ( t ) s o u t ( t ) + V S ˙ g e n ( t ) d V = d d t V ρ ( t ) s ( t ) d V
For stabilized pulsed flow (perfectly periodic), Equation (1) can be simplified by integrating over the inlet valve cycle duration. Actually, the temporal variations of the storage effects are cancelled out over the cycle. The system of Equation (1) then reduces to Equation (2).
Δ t c y c l e m ˙ i n d t Δ t c y c l e m ˙ o u t d t = 0 Δ t c y c l e m ˙ i n h t , i n d t Δ t c y c l e m ˙ o u t h t , o u t d t = Δ t c y c l e W ˙ d t Δ t c y c l e m ˙ i n s i n d t Δ t c y c l e m ˙ o u t s o u t d t + Δ t c y c l e V S ˙ g e n d V d t = 0
From Equation (2), turbine performance indicators over a valve cycle can be defined. The aim of the paper is to give design trends that extract the maximum energy from the flow. It is then necessary to compare the energy extracted during a cycle to the energy injected into the turbine in a cycle; see Equation (3). When η c y c l e = 1 , the turbine extracts all the energy from the flow. This efficiency can be negative when flow waves cause a reversal in the direction of the aerodynamic force applied to the blades. In this case, the system operates as a compressor.
η c y c l e = Δ t c y c l e W ˙ ( t ) d t Δ t c y c l e m ˙ i n ( t ) c p T t , i n ( t ) d t
The indicator Equation (3) is not an efficiency as usually found in the literature and should not be considered as such. It corresponds to the output for the various factorial combinations considered in this paper, allowing a fair comparison among different pulsating profiles. It could be interpreted as a recovery coefficient. Now, the factors that may influence this output are described.

2.2. Factors of the Design of Experiments

The system design is not only based on the turbine geometry, but also on the inlet and outlet boundary conditions. In fact, the geometry allows extracting energy from the flow while the inlet and outlet boundary conditions (this is also the case of the geometry since it modifies the flow physics) contribute to creating favourable flow conditions for the work extraction.

2.2.1. Geometric Parameters

Several assumptions were made to reduce the number of geometric factors. First of all, the stator and rotor cascades were taken as thin and two-dimensional. Cascades were built with Equation (4). Angles at the leading edge of cascades in relation to the direction of the inlet flow were taken to be zero. Moreover, the axial lengths of the stator and the rotor blades were considered identical.
y ( x ) = 4 tan α o u t π 2 1 + sin π 2 sin π 2 x x c
In addition to the blades’ geometrical parameters, the lengths between boundary conditions and the blades are significant for pulsed flows. They influence the location of the waves’ interaction in the domain. The length between the rotor trailing edge and the outlet boundary condition was considered very large to eliminate this parameter. The DOE geometric factors are shown in Figure 1.

2.2.2. Dynamic Parameters

The inlet valve cycle was chosen as a square wave. The opening and closing of the valve was instantaneous. During the opening phase, total pressure and total temperature were prescribed. A wall boundary condition was imposed during the closing phase. Under these assumptions, it took 4 factors to set up the valve cycle: P t , o p e n , T t , o p e n , Δ t c y c l e and Δ t o p e n . To set up the simulation, an outlet pressure condition is applied: P s , o u t l e t . The calculation information at the interface between the stator and rotor cascade was exchanged via a sliding mesh model, and rotor speed must be given: U r o t o r . The rotor speed was considered constant in the numerical simulations. The assumption was made that the valve cycles are much too fast for the rotor to adapt to the flow changes. All the other boundary conditions were supposed periodic.

2.2.3. Fluid Parameters

The analysed fluid was air. Thanks to the ideal gas law, the flow density is known. Only inviscid simulations were performed in this paper. Comparisons of large eddy simulations with inviscid simulations of transient flow within a 2D turbine cascade, carried out at the department [10], show that viscosity has no influence on the work prediction during the transient regime. No fluid parameter was added in the DOE factors.

2.2.4. Dimensionless Factors

The system behaviour law involves n = 13 -dimensional parameters and k = 3 fundamental units. Thanks to Vaschy–Buckingham’s theorem, the behaviour law can be determined with n k = 10 dimensionless parameters. The dimensionless factors of the design of experiments are listed in Table 2.
Most of the dimensionless factors are relatively easy to understand since they are usually used in steady-state flow turbines. However, some parameters, specific to pulsed flow, need to be specified. Π 10 represents the cycle ratio, while Π 5 and Π 6 set, partly, the interaction wave locations. Π 9 gives an indication about the distance travelled by a wave during Δ t c y c l e . Indeed, r T t , o p e n is similar to the sound speed; thus, Π 9 = 1 means that a wave propagating at the sound speed travels x c during Δ t c y c l e .
The factor range of variation is also presented in Table 2. These ranges are large since the regions of the experiment domain that led to the best performance were unknown at the beginning of the study. The ranges were centred and reduced between 1 and + 1 .

3. Results

3.1. Most Influential Factors

The design of experiments was built with the JMP® [11] software. To determine the factors’ influence on the output, a two-stage screening design was selected. This kind of design makes it possible to sort out from a large sample of factors those that have a non-negligible role on the system output and those that have little influence. These designs lead to a reduction in the dimensions of the factor space and a simplification of the problem. In this screening design, the level of each factor corresponds to the boundaries of the variation range, i.e., + 1 and 1 as a coded variable. Consequently, the predicted output law is linear. The experimental design selected is a fractional factorial design. These designs are subsets of the full factorial designs. Fractional factorial designs are a good choice when resources are limited and the number of factors in the design is large; see [12] for more information. The simulations’ number of fractional designs was equal, here 2 n p = 2 10 6 = 16 . The screening design is listed in Table 1.
Screening design results are shown in Table 1. Two additional responses besides the recovery efficiency η were proposed. Since the rotational speed is fixed, the strong decrease of mass-flow can make the compressor operating mode appear. The relative proportion of time in which this happens is quoted in %. Furthermore, the standard deviation of the work signal is reported in order to quantify the loading fluctuations. The recovery efficiency of many simulations is close to η c y c l e = 0 ; the turbine extracts almost no energy from the flow over a cycle. The system alternates between turbine and compressor phases. The balance between these modes is particularly visible on the temporal evolution of W ˙ ( t ) for the simulation n 3 in Figure 2. Compressor modes are driven by waves that generate a greater force on the suction side of the blade than on the pressure side. The temporal evolution of W ˙ ( t ) shows behaviours in agreement with the results of [3]. The valve opening causes a shock wave propagation, which generates a work increase. The valve closing generates an expansion wave, which causes a work decrease. An instantaneous visualization of the density gradient is also shown in Figure 2. This instantaneous visualization highlights the flow complexity.
The linear prediction extracted from the screening phase is given in Figure 3. It shows a fair accuracy R 2 = 0.98 . Moreover, the distribution of the measured values in the observed performance range is relatively homogeneous, which gives confidence in the behavioural law prediction and legitimises the conclusions regarding the true influence of the different factors.
For the numerical experimental designs, the natural variability of the output can be considered as negligible. As a result, all factors in the design of experiments are statistically significant on the system output. However, some factors may be neglected by comparing the sensitivity values of the prediction model. A factor was considered to be noninfluential on the response when its sensitivity was less than 5 % of the maximum sensitivity of the model.
The sensitivity analysis showed that only six factors have a significant role in stage performance; see Table 3. The influence of these factors was further investigated by calculating the response surface. Understanding in detail why factors ( α o u t , r o t o r , L i n l e t , s t a t o r / x c , L s t a t o r , r o t o r / x c and Δ t c y c l e r T t , o p e n ) do not influence the turbine recovery efficiency is difficult with such a flow complexity. However, the spatial distribution of the time-averaged energy recovered shows no sensitivity to these parameters. This result will not be presented here.

3.2. Surface Response

The reduction of the number of factors makes it possible to carry out a more quantitative study of the system output thanks to the response surface methodology [12]. The response surface modelling allows determining the factor that optimizes the output. For this design of experiments, the behaviour law prediction was performed by a quadratic function (see Equation (5)) with x i each selected factors.
y = y ¯ + i λ i x i + γ i x i 2 + i j , j i β i , j x i x j
In order to achieve a quadratic prediction of the output law, three levels for each factor must be considered in the simulations. In addition to the high + 1 and low levels 1 , the central level 0 was added. For the calculation, the stator and rotor channels must be multiples of each other. The central value is then adjusted for h s t a t o r / x c and h r o t o r / x c . The level of factors is indicated in Table 4. In the same way as the screening design, the experimental design was constructed using [11]. In order to minimize the simulations number, a D-optimality criterion was adopted for the design of experiments; see [12]. The surface response design was based on twenty-eight simulations for six factors. The levels of the noninfluential factors on the output, namely Π 4 = α o u t , Π 5 = L i n l e t , s t a t o r / x c , Π 6 = L s t a t o r , r o t o r / x c , and Π 9 = Δ t c y c l e r T t , o p e n / x c , were taken at their high level + 1 to carry out the simulations.
The results of the experimental design are provided in Table 5. The response surface accuracy can be examined by investigating Figure 4. It shows that the output law is very well modelled by the response surface R 2 = 0.99 . The spatial distribution of the measured values in the observed performance range is relatively homogeneous, which ensures that the predicted output law can be trusted. The model quality is also assessed on the residual value ϵ i = y i y ^ i (Figure 4). The residuals are small and have no outliers.
It is interesting to focus on a few particular simulations showing original behaviours. Simulations n 17 or n 18 revealed that it is possible, with the wrong set of parameters, to design a system operating exclusively in compressor mode, whereas the target was to design a turbine. n 25 shows a fairly marked compressor mode over a cycle (13.5%), whereas this simulation leads to the highest recovery efficiency observed.
The factor levels that optimize the response are shown in Table 6. The maximum efficiency in the experimental domain is η c y c l e = 30.1 % . The result on the coded variables reveals that the real optimum of this system is outside the experimental domain. Indeed, except for α o u t , s t a t o r , the optimal value of the factors is located on the boundary of the experimental domain. This is promising since it means that it is possible with this system to extract more than 30 % of the injected energy into the turbine over a valve cycle, if the experimental domain is redefined in order to find the global optimum. For comparison, the η c y c l e value calculated for the usual turbofan is about 0.27 for the high-pressure turbine stage ( T t , i n l e t H P 1800 K, T t , o u t l e t H P 1300 K).
The position of the optimum gives the main trends for the design of a turbine fed by a pulsed flow, as far as 2D analysis can be applied. It reinforces the usefulness of the stator for this system. This result was not evidence since the instantaneous performance is dictated by the waves’ propagation, and not by the usual steady analysis expressed by Euler’s theorem. The stator must be composed of blades with a high solidity ( h s t a t o r / x c ) , which forces a large deviation ( α o u t , s t a t o r ) to the flow. This kind of stator generates much more intense wave reflections and diffraction than if h s t a t o r / x c and α o u t , s t a t o r were kept at their lowest levels. The reflections’ amplification causes an inlet energy flux reduction during the opening phases. The diffraction intensification, especially at the trailing edge, causes a large reduction in the transmitted waves’ intensity to the rotor. The rotor is then less sensitive to the flow unsteadiness. This behaviour can be seen by comparing the standard deviation of the work signal between Simulations n 7 and n 27 , where only stator geometrical parameters are modified.
The rotor solidity h r o t o r / x c must also be high (Table 6). The sensitivity estimation of the response surface shows that this factor is the most important for the turbine design. Observations of the DOE results (Table 5) prove that it is impossible to reach high efficiency when h r o t o r / x c is low. It is possible that this is due to the large recirculation zones that take place at the leading edge on the suction side of the rotor blades, and that persist more or less over time, when the stator geometrical factors are at their highest levels. When the size of these zones is of the order of magnitude of the pitch, pressure profiles on both sides of the rotor blade tend to be the same. Therefore, h r o t o r / x c must be high enough to overcome this difficulty.
Table 6 reveals that the turbine is more efficient to extract energy when high-pressure ratios are applied over a short ratio cycle. In addition to the stator design, low cycle ratios attenuate the energy injected into the turbine over a cycle, while high-pressure ratios increase it. This may seem contradictory. In addition, the unsteadiness related to the shock wave propagation during the valve opening is maximal for these cycle features. Indeed, for low cycle ratios, the upstream flow of the shock wave is mostly at rest in the stator. The downstream shock state is set by the inlet boundary condition, and the shock wave intensity is then maximal for this cycle feature (this explains the compressor phases on Simulation n 25 ). As explained, the stator reduces the unsteadiness, which propagates to the rotor. A compromise must therefore be found among h s t a t o r / x c , α o u t , s t a t o r , P t , o p e n / P s , o u t l e t and Δ t o p e n / Δ t c y c l e in order to reduce the energy flux through the turbine while maximizing the unsteadiness benefits within the rotor.
Finally, in order to minimize the drawbacks related to unsteadiness, more specifically compressor operation during the valve closing phase, Π 8 must be as low as possible. The smaller Π 8 , the lower the rotor speed compared to the characteristic bulk velocity.

4. Conclusions

Preliminary 2D design recommendations for a turbine driven by a pulsed flow were given in this paper. In order to catch the main trends for the turbine design, a parametric study was carried out. Recommendations were given on geometric and dynamic parameters that maximize a recovery ratio based on the amount of energy entering the control volume defined by the stage. Ten dimensionless parameters were selected as the input of the degree of the experiment. Response surface modelling was carried out on the influential factors following the screening design. The response surface calculation allowed showing the turbine characteristics maximizing η c y c l e ; see Figure 5. The set of factors α o u t , s t a t o r , h s t a t o r / x c , P t , o p e n / P s , o u t l e t , Δ t o p e n / Δ t c y c l e must allow finding a balance between minimizing the energy injected into the turbine and amplifying the benefits linked to the unsteadiness caused by the shock wave propagation during the valve opening. In addition, the rotor solidity h r o t o r / x c needs to be high so that the pressure distribution on both sides of the blade is not only controlled by the recirculation zones. Finally, the rotor speed must be low compared to the average bulk velocity in order to avoid the occurrence of compressor mode.
Now that a better knowledge of the factor influence and range has been identified, a gradual increase in geometrical complexity is scheduled, in order to integrate the 3D effects and thickness distribution.

Author Contributions

Conceptualization, F.H., N.B. and J.G.; software, F.H., J.G. and G.S.-M.; formal analysis, F.H.; investigation, F.H.; visualization, F.H.; data curation, F.H.; writing—original draft preparation, F.H.; writing—review and editing, F.H., N.B. and J.G.; supervision, N.B. and J.G.; project administration, N.B. and J.G.; funding acquisition, N.B. and J.G. All authors read and agreed to the published version of the manuscript.

Funding

This research was funded by the Direction Générale de l’Armement (DGA).

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

The authors would like to thank the DGA for their financial support and for permitting the publication of the research. This work was performed using HPC resources from GENCI-IDRIS (Grant 2020-100431), GENCI-CINES (Grant 2020-A0082A07178) and CALMIP (Grant 2020-p1425).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; nor in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
α o u t Blade metal angle
LDistance
Δ t c y c l e Cycle time
P s , o u t l e t Static outlet pressure
Δ t o p e n Opening phase time
P t , o p e n Total pressure during the opening phase
m ˙ Mass flow rate
rGas constant
sEntropy
W ˙ Work turbine
T t Total temperature
σ Standard deviation
x c Axial chord
h r o t o r Rotor pitch
h s t a t o r Stator pitch
h t Total enthalpy
Subscripts
inInlet
outOutlet
HPHigh pressure

References

  1. Baines, N.C. Turbocharger Turbine Pulse Flow Performance and Modelling—25 Years on; Engineering Conferences Online; EIT, Inc.: Egg Harbor Township, NJ, USA, 2010. [Google Scholar]
  2. Liu, Z.; Copeland, C. Optimization of a Radial Turbine for Pulsating Flows. J. Eng. Gas Turbines Power 2020, 142, 051009. [Google Scholar] [CrossRef]
  3. Hermet, F.; Binder, N.; Gressier, J. Transient flow in infinitely thin airfoil cascade. In Proceedings of the 13th European Conference on Turbomachinery Fluid Dynamics & Thermodynamics ETC13, Lausanne, Switzerland, 8–12 April 2019. [Google Scholar]
  4. Brès, G.A.; Ham, F.E.; Nichols, J.W.; Lele, S.K. Unstructured large-eddy simulations of supersonic jets. AIAA J. 2017. [Google Scholar] [CrossRef]
  5. Bohbot, J.; Grondin, G.; Corjon, A.; Darracq, D. A parallel multigrid conservative patched/sliding mesh algorithm for turbulent flow computation of 3D complex aircraft configurations. In Proceedings of the 39th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 8–11 January 2001. [Google Scholar]
  6. Hirt, C.; Amsden, A.; Cook, J. An arbitrary Lagrangian-Eulerian computing method for all flow speeds. J. Comput. Phys. 1974, 14, 227–253. [Google Scholar] [CrossRef]
  7. Luo, H.; Baum, J.D.; Löhner, R. On the computation of multi-material flows using ALE formulation. J. Comput. Phys. 2004, 194, 304–328. [Google Scholar] [CrossRef]
  8. Ansys ICEM CFD User’s Manual; Version 15.0; ANSYS, Inc.: Canonsburg, PA, USA, 2013.
  9. Lee, J.; Tan, C.S.; Sirakov, B.T.; Im, H.S.; Babak, M.; Tisserant, D.; Wilkins, C. Performance metric for turbine stage under unsteady pulsating flow environment. J. Eng. Gas Turbines Power 2017, 139, 072606. [Google Scholar] [CrossRef]
  10. Hermet, F. Simulation des transitoires violents et des écoulements pulsés dans les turbines. Ph.D. Thesis, ISAE-Supaero, Toulouse, France, 2021. [Google Scholar]
  11. JMP®; Version 15.0; SAS Institute Inc.: Cary, NC, USA, 1989–2019; Available online: https://www.jmp.com (accessed on 20 January 2020).
  12. Goupy, J. Introduction aux plans d’expériences, 2nd ed.; Technique et ingénierie Série Conception; Dunod: Paris, France, 2001. [Google Scholar]
Figure 1. Geometric parameters.
Figure 1. Geometric parameters.
Ijtpp 06 00024 g001
Figure 2. Above: Temporal evolution of W ˙ ( t ) . W ˙ < 0 corresponds to compressor mode. Below: Instantaneous visualization of the density gradient. Simulation n 3 .
Figure 2. Above: Temporal evolution of W ˙ ( t ) . W ˙ < 0 corresponds to compressor mode. Below: Instantaneous visualization of the density gradient. Simulation n 3 .
Ijtpp 06 00024 g002
Figure 3. Prediction of the screening design. R 2 = 0.98 .
Figure 3. Prediction of the screening design. R 2 = 0.98 .
Ijtpp 06 00024 g003
Figure 4. Predictions of the surface response design. R 2 = 0.998 .
Figure 4. Predictions of the surface response design. R 2 = 0.998 .
Ijtpp 06 00024 g004
Figure 5. Design recommendation for a turbine driven by a pulsed flow. Recommendations are identified by blue lines. The gray dashed lines show the experimental domain. It should be noted that the speed triangle shown is a pictorial, albeit very inadequate, way of characterizing the rotor speed in relation to a semblance of flow speed.
Figure 5. Design recommendation for a turbine driven by a pulsed flow. Recommendations are identified by blue lines. The gray dashed lines show the experimental domain. It should be noted that the speed triangle shown is a pictorial, albeit very inadequate, way of characterizing the rotor speed in relation to a semblance of flow speed.
Ijtpp 06 00024 g005
Table 1. Screening design and associated results.
Table 1. Screening design and associated results.
N Π 1 Π 2 Π 3 Π 4 Π 5 Π 6 Π 7 Π 8 Π 9 Π 10 η   ( % ) σ W ˙ ( t ) < W ˙ ( t ) > Compressor   ( % )
1 1 + 1 1 1 + 1 1 + 1 1 + 1 1 2.9 4.0 50.0
2 + 1 + 1 + 1 + 1 + 1 + 1 + 1 1 1 1 10.9 0.3 0.0
3 + 1 1 1 1 1 + 1 + 1 + 1 + 1 1 0.02 2.7 70.0
4 1 + 1 1 + 1 1 + 1 1 + 1 1 + 1 2.8 0.5 3.0
5 + 1 + 1 + 1 1 1 1 1 + 1 + 1 + 1 3.0 0.3 0.0
6 1 1 + 1 + 1 1 1 + 1 + 1 1 1 2.1 1.5 25.0
7 + 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 1 19.0 1.5 28.0
8 + 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 1 + 1 22.0 0.2 0.0
9 + 1 + 1 + 1 + 1 + 1 + 1 + 1 1 + 1 + 1 13.0 0.1 0.0
10 + 1 + 1 + 1 + 1 + 1 + 1 1 + 1 + 1 + 1 2.8 0.3 0.0
11 + 1 + 1 + 1 + 1 + 1 1 + 1 + 1 + 1 + 1 24.1 0.1 0.0
12 + 1 + 1 + 1 + 1 1 + 1 + 1 + 1 + 1 + 1 24.2 0.1 0.0
13 + 1 + 1 + 1 1 + 1 + 1 + 1 + 1 + 1 + 1 24.0 0.3 0.0
14 + 1 + 1 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 14.3 0.5 0.0
15 + 1 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 3.0 0.2 0.0
16 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 + 1 24.2 0.5 0.0
Table 2. Factors and experimental domain of the screening design.
Table 2. Factors and experimental domain of the screening design.
Π 1 Π 2 Π 3 Π 4 Π 5
h s t a t o r / x c h r o t o r / x c α o u t , s t a t o r α o u t , r o t o r L i n l e t / s t a t o r / x c
Real [ 0.1 ; 0.8 ] [ 0.1 ; 0.8 ] [ 40 ; 80 ] [ 40 ; 80 ] [ 0.1 ; 2.0 ]
Coded [ 1 ; + 1 ] [ 1 ; + 1 ] [ 1 ; + 1 ] [ 1 ; + 1 ] [ 1 ; + 1 ]
Π 6 Π 7 Π 8 1 Π 9 Π 10
L s t a t o r / r o t o r / x c P t , o p e n / P s , o u t l e t U r o t o r / V θ Δ t c y c l e r T t , o p e n / x c Δ t o p e n / Δ t c y c l e
Real [ 0.1 ; 1.0 ] [ 1.2 ; 3.0 ] [ 0.25 ; 0.75 ] [ 0.5 ; 5.0 ] [ 0.1 ; 0.9 ]
Coded [ 1 ; + 1 ] [ 1 ; + 1 ] [ 1 ; + 1 ] [ 1 ; + 1 ] [ 1 ; + 1 ]
1 V θ = f ( Π 7 , Π 3 , r T t , o p e n ) corresponds to the tangential velocity at the trailing edge of the stator during the opening phase, assuming that all the expansion takes place in this stage.
Table 3. Factors’ influence on the output. ✓ is associated with an influential factor, while ✕ is related to a noninfluential factor.
Table 3. Factors’ influence on the output. ✓ is associated with an influential factor, while ✕ is related to a noninfluential factor.
Π 1 Π 2 Π 3 Π 4 Π 5
h s t a t o r / x c h r o t o r / x c α o u t , s t a t o r α o u t , r o t o r L i n l e t / s t a t o r / x c
Π 6 Π 7 Π 8 Π 9 Π 10
L s t a t o r / r o t o r / x c P t , o p e n / P s , o u t l e t U r o t o r / V θ Δ t c y c l e r T t , o p e n / x c Δ t o p e n / Δ t c y c l e
Table 4. Factors and experimental domain of the surface response design.
Table 4. Factors and experimental domain of the surface response design.
Π 1 Π 2 Π 3
h s t a t o r / x c h r o t o r / x c α o u t , s t a t o r
Coded [ 1 , 0.14 , + 1 ] [ 1 , 0.14 , + 1 ] [ 1 , 0 , + 1 ]
Real [ 0.1 , 0.4 , 0.8 ] [ 0.1 , 0.4 , 0.8 ] [ 40 , 60 , 80 ]
Π 7 Π 8 Π 10
P t , o p e n / P s , o u t l e t U r o t o r / V θ Δ t o p e n / Δ t c y c l e
Coded [ 1 , 0 , + 1 ] [ 1 , 0 , + 1 ] [ 1 , 0 , + 1 ]
Real [ 1.2 , 2.1 , 3.0 ] [ 0.25 , 0.5 , 0.75 ] [ 0.1 , 0.5 , 0.9 ]
Table 5. Surface response design selected and associated results.
Table 5. Surface response design selected and associated results.
N Π 1 Π 2 Π 3 Π 7 Π 8 Π 10 η   ( % ) σ W ˙ ( t ) < W ˙ ( t ) > Compressor   ( % )
1−1−1−1−1+10−0.031.666.0
2−1−1−1+1−1−1−1.043.150.0
3−1−10−1−1+1−0.420.3100
4−1−1+1−10−1−0.084.958.0
5−1−1+1+1+1+1−0.320.2100
6−1+1−1−1−1−11.4614.048.0
7−1+1−1+1+1+15.20.63.4
8−1+1+1−1+1+11.960.10.0
9−1+1+1+1−1+16.470.10.0
10 0.14 0.14 −100+10.0113.349
11 0.14 0.14 00+1−10.321.049.8
12 0.14 0.14 0+1000.0113.050.0
13 0.14 0.14 +10−10−0.032550.0
14−1+1000016.50.81.8
15+1−10−1−1−1−0.0129.041.0
16+1−1−1+1+1+1−0.530.497
17+1−1+1−1+1+1−1.750.4100
18+1−1+1+1−1+1−3.90.1100
19+1−100000.0150.050.0
20+1+1−1−1+1+111.90.53.4
21+1+1−1+1−1+123.20.31.8
22+1+1−1+1+1−13.413.555.0
23+1+1+1−1−1+116.20.20.0
24+1+1+1−1+1−12.923.052.0
25+1+1+1+1−1−129.41.313.5
26+1+1+1+1+1−119.12.335.0
27+1+1+1+1+1+124.30.20.0
28+1−1+1+1+1+13.00.30.0
29+1+1−1+1+1+114.30.30.0
30+1+1+1−1+1+19.00.30.0
Table 6. Factors optimizing the turbine performance indicator in the experimental domain η c y c l e , m a x = 30.1 % .
Table 6. Factors optimizing the turbine performance indicator in the experimental domain η c y c l e , m a x = 30.1 % .
Π 1 Π 2 Π 3 Π 7 Π 8 Π 10
h stator / x c h rotor / x c α out , stator P t , open / P s , outlet U rotor / V θ Δ t open / Δ t cycle
Coded + 1 + 1 + 0.75 + 1 1 1
Real 0.8 0.8 75 3.0 0.25 0.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hermet, F.; Binder, N.; Gressier, J.; Sáez-Mischlich, G. Pulsed Flow Turbine Design Recommendations. Int. J. Turbomach. Propuls. Power 2021, 6, 24. https://0-doi-org.brum.beds.ac.uk/10.3390/ijtpp6030024

AMA Style

Hermet F, Binder N, Gressier J, Sáez-Mischlich G. Pulsed Flow Turbine Design Recommendations. International Journal of Turbomachinery, Propulsion and Power. 2021; 6(3):24. https://0-doi-org.brum.beds.ac.uk/10.3390/ijtpp6030024

Chicago/Turabian Style

Hermet, Florian, Nicolas Binder, Jérémie Gressier, and Gonzalo Sáez-Mischlich. 2021. "Pulsed Flow Turbine Design Recommendations" International Journal of Turbomachinery, Propulsion and Power 6, no. 3: 24. https://0-doi-org.brum.beds.ac.uk/10.3390/ijtpp6030024

Article Metrics

Back to TopTop