A high-fidelity real-time simulation code based on a lumped, nonlinear representation of gas turbine components is presented. The code is a general-purpose simulation software environment useful for setting up and testing control equipments. The mathematical model and the numerical procedure are specially developed in order to efficiently solve the set of algebraic and ordinary differential equations that describe the dynamic behavior of gas turbine engines. For high-fidelity purposes, the mathematical model takes into account the actual composition of the working gases and the variation of the specific heats with the temperature, including a stage-by-stage model of the air-cooled expansion. The paper presents the model and the adopted solver procedure. The code, developed in Matlab-Simulink using an object-oriented approach, is flexible and can be easily adapted to any kind of plant configuration. Simulation tests of the transients after load rejection have been carried out for a single-shaft heavy-duty gas turbine and a double-shaft aero-derivative industrial engine. Time plots of the main variables that describe the gas turbine dynamic behavior are shown and the results regarding the computational time per time step are discussed.

A method for predicting the onset of acoustically driven combustion instabilities in gas turbine combustor is examined. The basic idea is that the governing equations of the acoustic waves can be coupled with a flame heat release model and solved in the frequency domain. The paper shows that a complex eigenvalue problem is obtained that can be solved numerically by implementing the governing equations in a finite element code. This procedure allows one to identify the frequencies at which thermo-acoustic instabilities are expected and the growth rate of the pressure oscillations, at the onset of instability, when the hypothesis of linear behavior of the acoustic waves can be applied. The method can be applied virtually to any three-dimensional geometry, provided the necessary computational resources that are, anyway, much less than those required by computational fluid dynamics methods proposed for analyzing the combustion chamber under instability condition. Furthermore, in comparison with the “lumped” approach that characterizes popular acoustics networks, the proposed method allows one for much more flexibility in defining the geometry of the combustion chamber. The paper shows that different types of heat release laws, for instance, heat release concentrated in a flame sheet, as well as distributed in a larger domain, can be adopted. Moreover, experimentally or numerically determined flame transfer functions, giving the response of heat release to acoustic velocity fluctuations, can be incorporated in the model. To establish proof of concept, the method is validated at the beginning against simple test cases taken from literature. Over the frequency range considered, frequencies and growth rates both of stable and unstable eigenmodes are accurately evaluated. Then the method is applied to a much more complex annular combustor geometry in order to evaluate frequencies and growth rates of the unstable modes and to show how the variation in the parameters of the heat release law can influence the transition to instability.

No abstract

No abstract

Nowadays Pumps working as Turbines (PaT) are devices widely used to perform energy recovery in hydraulic grids, thus improving their overall efficiency, and to build small hydropower plants. In this work, a centrifugal pump has been numerically investigated in turbine operating mode by means of the open-source CFD code OpenFOAM with emphasis on the flow field at the runner outlet. Due to the reduced number of blades in a PaT, the mean outlet relative velocity angle differs from the blade angle. In order to account for this phenomenon, the slip factor is introduced. The slip factor is investigated and its application to a 1D model is shown in order to highlight the improvement in predicting the characteristic curve of a centrifugal pump used in reverse mode as a turbine (PaT) especially at its part-load.