Research Article  Open Access
The OffDesign Performance Simulation of Marine Gas Turbine Based on Optimum Scheduling of Variable Stator Vanes
Abstract
As one of the antisurge techniques, the adjusting scheme of VSV under offdesign conditions has a significant impact on the performance of gas turbines. In this paper, the onedimensional characteristic of the compressor calculation program is embedded into the zerodimensional overall gas turbine model, which replaces the original compressor characteristic module. Based on the assembling relationship of the actual components of the marine gas turbine, the architecture of the modular model library is designed, and an integrated simulation platform of marine gas turbine is developed by using MATLAB/GUI software. The influence of the first 3 rows of variable stator vanes of the 9stage axial compressor working alone on the performance of the compressor at different speeds and different angles was analyzed by the HARIKA compressor characteristic calculation program. Taking the economics and stability of the gas turbine as the optimization objective, the optimization of the first threestage stator vanes regulation schemes under different working conditions was carried out. The steadystate performance parameters under each working condition of gas turbine of power generation with or without variable stator vane mode were calculated. The study results can provide references for the adjusting scheme of VSV under gas turbine offdesign conditions operating process.
1. Introduction
As an important part of the gas turbine, the compressor should have high performance and enough margin and have higher efficiency [1]. Under the offdesign condition, to avoid entering the unstable working state and ensure enough operating performance requirements, the variablegeometry compressor technology has become the mainstream. The matching and regulation of the variablegeometry compressor will directly determine the working performance of the gas turbine. In order to improve the operating efficiency of the gas turbine and to improve the matching relationship between the compressor and the gas turbine, how to adjust the stator blade angle will become a problem that researchers need to face. Compared with the experimental method to determine the regulation scheme of adjustable stationary blades, using the simulation method not only can save a lot of manpower and material resources but also can determine a reasonable adjustable stationary blades adjustment scheme in a short time.
Today, gas turbines are developing in the direction of complexity, scale, and diversification. The simulation process of a gas turbine becomes more complex [2]. Many researchers at home and abroad have carried out research in this area. In the late 1990s, National Aeronautics and Space Administration (NASA), united airlines industries, universities, and government institutions proposed the Numerical Propulsion System Simulation (NPSS) platform. This platform, using the simulation algorithm and the discipline coupling technology which can freely zoom and change dimensions, deeply reveals the complex flow inside the combustion engine and the mutual relations among the components and realizes the detailed simulation of the whole machine [3, 4]. Simulation flow chart of high precision component modeling based on NPSS framework is shown in Figure 1. NPSS mainly contains five key technologies: () standard data interface; () using a flexible, modular, objectoriented program structure; () simulation and analysis of the specific physical process of the aeroengine according to the needs of the specific physical process; () multidisciplinary integrated technology; and () high performance parallel and distributed computing technology [5].
The Russian Central AeroEngine Research Institute (CIAM) developed the Computer Gas Turbine Engine Simulator (CGTES) to support the CT3 (Turbine Engine Test Technology) program [6]. CGTES is an independent computing system designed for the development of aeroengines. The simulation system can calculate the steadystate and transientstate parameters. The normal operation of the program is controlled by the special task input system. However, the precision of this kind of simulation technology is mainly affected by the establishment of a loss model and the selection of an empirical coefficient.
The offdesign performance of marine gas turbines has great significance for the economy and stability of ships. As one of the surgepreventing measures of the compressor, variable stator vane (VSV) can be used to change the installation angle of the cascade with the working conditions changing to improve the compressor performance in the offdesign conditions [7–11].
VSV technology has a long history and has been used in various applications of gas turbines. In the early 1940s, Germany’s first turbojet engine Jumo004 applied the technology of VSV [12]. Figure 2 shows the stage numbers of VSV and the total stage numbers in several generations of aeroengine compressors.
In the field of ground power generation and marine gas turbines, VSV technology is also widely used [13, 14]. US GE company’s productions, 9FAtype gas turbine and PG6541type singleaxis gas turbine, both use inlet guide vanes, while the Swiss Sulzer company uses VSV technology in its AV series axial compressor. As a typical case in marine gas turbines, LM2500 gas turbine’s top seven stages of its inlet guide vane and 16stage stator blades are variable stator vanes [15, 16]. Hu et al. established a method for calculating the characteristics of a multistage axialflow compressor based on the inclusion of axialflow compressor flow loss and lag angle model combining stepwise primitive blade model [17–20]. And based on the genetic evolutionary algorithm of the tournament, a twoobjective genetic optimization algorithm is obtained. Zhang and Ren studied the influence of the adjustment of VSV on the compressor performance through the method of experiment. In order to improve the stage matching relationship and improve the performance of the compressor, an experimental procedure of adjusting the installation angle is introduced [21].
Through the review of the research status, we can know that the current research on VSV mainly focuses on the effect of VSV on compressor performance. The research on the influence of VSV on the performance of the whole gas turbine is relatively scarce. The research on the adjustment regulation of VSV mainly focuses on the performance of the compressor, and the research on the performance of the whole gas turbine is relatively less. Under the background of gas turbine integrated simulation, the method of simulating research on a certain part of the gas turbine from different perspectives is gradually mature by using the integrated modeling resources. Aiming at the problem mentioned above, this paper takes one threeshaft marine gas turbine for generation as the research object, adopts the modular modeling idea, establishes the general simulation model of (01D) variable dimension gas turbine considering a variablegeometry compressor, and optimizes the adjustment scheme of the top three stages of the variable stator vanes of the lowpressure compressor, and the influence of the variable stator vanes on the steadystate performance of the gas turbine was studied emphatically. The main contents of this paper are as follows.
Design an integrated simulation platform for marine gas turbines and build a simulation model library including a compressor, combustion chamber, turbine, rotor, and volume. According to the characteristics of the variablegeometry compressor, build the new compressor characteristics modules and interface modules.
A variable dimension turbine simulation model considering the variablegeometry compressor is established.
From the angle of the performance of the whole gas turbine, taking the top three grades of the variable stator vanes of the lowpressure compressor as the optimization variable, the economic efficiency and stability of the gas turbine under various working conditions are optimized. The results show that the optimized variable stator vane adjustment scheme can improve the performance of the gas turbine.
2. Methodology
2.1. Collaborative Running Interface Module
As the core part of the variablegeometry compressor model, the compressor interface module is the key to ensure the data transfer between the zerodimensional models and onedimensional characteristic calculation program. The whole interface module was built by Sfunction/MATLAB. The function of the compressor interface module mainly contains two parts: () to call the onedimensional compressor performance calculation program and () to complete the update and data transfer of onedimensional computing program’s input files and output files.
The process of collaborative running is shown in Figure 3. When the simulation begins, the rotating speed and pressure ratio of the compressor can be obtained through the rotor module and the volume module. The flow coefficient and static blade opening adjustment value will also be given.
When Simulink is working, it first enters the simulation initialization phase. In this phase, the solver determines all the parameters of the module through the callback function and passes the signal width, data type, and state number. Meanwhile, the execution order, initial value, and the sampling time of each module are also determined.
Finally, other tasks in the initialization phase are performed before the simulation cycle phase. After entering the simulation cycle, the Simulink solver calculates the output of each module according to the predetermined order and uses the callback function to calculate the module status in the current sampling time, update the discrete state value, and calculate the differential value and the output value of the corresponding module.
After the simulation cycle phase, the Simulink solver calls the callback function to release the system storage space and performs other tasks that need to be done after the simulation step is completed.
The input files of HARIKA include the following: import conditions of variablegeometry compressors (i.e., import total pressure, import total temperature, and import airflow angle), relative speed, the main geometric size at all levels, and parameters such as the opening of variable static blade.
The output file mainly involves the compressor characteristic parameters including pressure ratio, flow rate, speed, and efficiency, which are calculated by the calculation program of compressor characteristics.
2.2. Mathematical Model of Marine Gas Turbine with VariableGeometry Compressor
2.2.1. 1D VariableGeometry Compressor Characteristic Model
When multiple variable stator vane angles change, the characteristics of the compressor will be changed in different angles. So, using only the zerodimensional simulation cannot meet the need of modeling. Onedimensional simulation is mainly related to the component loss model, backward angle model of rotating parts, and some empirical coefficients. Therefore, the onedimensional compressor performance calculation program can accurately reflect the variablegeometry compressor characteristics.
In this paper, HARIKA algorithm was used to calculate the onedimensional compressor characteristic. Based on a large number of test data, HARIKA algorithm uses the statistical relationship of the stage synthesis characteristics to determine the stage parameters under offdesign condition. The algorithm has been empirically modified along the vane height, and it has been empirically modified for various characteristic parameters of vane type as well. After obtaining the total pressure ratio and efficiency of one stage of the compressor, the velocity triangle at the outlet of this stage is further obtained. Using the velocity triangle as the inlet parameter of the next stage of the compressor, the parameters at the outlet of the next stage are obtained by the same method. The characteristics of the whole compressor can be obtained by repeating the process. The algorithm is more suitable as the calculation program of the characteristics of variablegeometry compressor components.
Figure 4 shows the compressor characteristics of experiment and calculation. The results show that the maximum error of the experiment and calculation is 4.8%.
2.2.2. Compressor Performance Calculation
According to the compressor inlet temperature and compressor pressure ratio, we can determine the logarithmic ratio of the compressor outlet. Then, we determine the compressor exit enthalpy, and we can calculate the outlet temperature and compressor power.
The specific calculation process is as follows:
2.2.3. Combustion Chamber Model
During the process of establishing the model of combustion, we always consider the combustion chamber as a uniform field.
for the outlet of the combustion chamber pressure and and for inlet and outlet temperature of the combustion chamber; , , and are the flow at the inlet and outlet of the combustion chamber and fuel injection quantity; and are the combustion chamber volume and index adiabatic of combustion chamber; and are the specific enthalpy at the inlet and outlet of the combustion chamber; is the low calorific value of fuel; is specific heat capacity at constant pressure of the combustion chamber volume.
2.2.4. Turbine Model
The turbine module is similar to the compressor module, and the mathematical model is as follows:
where is the turbine inlet temperature; and are the inlet and outlet turbo pressure; is the turbine corrected speed; is the turbine rotor speed; is the ratio of turbine expansion; is the turbine corrected mass flow; is turbine efficiency.
2.2.5. Rotor Model
From the moment of momentum theorem, the product of the angular velocity of the rotor and the moment of inertia is equal to all the external torque on the rotor.
where is the rotor inertia; is the turbine torque; is the compressor torque; is the mechanical loss torque; is the load torque; is the turbine power; is the power consumption of the compressor; is the mechanical loss power; is the power loss load.
2.2.6. Volume Module Model
With the connection of the main components of the gas turbine, the flow space of the volume module can be regarded as global. In the volume, the flow rate of the working fluid is fast and the retention time is short, so we assume that the medium of heat transfer and import and export pressure is zero and assume central pressure instead of the average pressure of the working fluid inside the volume. According to the principle of conservation of mass, there is the following relationship:where is the volume of inlet temperature; is the volume internal average pressure; and are the import and export flows of volume.
3. Application, Results, and Analysis
This paper takes the triaxial gas turbine of a certain type as the physical model (as it is shown in Figure 5). Performance specifications of the gas turbine at standard sea level static ISA condition is shown in Table 1. The comparison results of the various parameters of the gas turbine engine used in the marine power generation are shown in Figures 12–15 and Table 3.

3.1. The Gas Turbine Simulation Model
In the environment of MATLAB/SIMULINK, a simulation model of variable specific heat capacity of three shaft gas turbine was built, which includes compressor, combustion chamber, turbine, volume, rotor and engine controller, which is shown in Figure 6.
3.2. Optimized Scheme of VSV for Marine Gas Turbine
Based on the gas turbine’s overall performance, this part of the paper has considered variablegeometry compressor’s (01) variable dimensions of the gas turbine overall simulation model, aiming at optimizing the gas turbine’s economics by taking angles of variable stator vane as the optimized variables. In the process of selecting optimized variables, the adjustable angles of the variable stator vane have been taken as optimized variables to improve the design of stator vane adjustment scheme. The optimization target is the economy of the gas turbine, whose specific measurement index is the fuel consumption rate and efficiency during the operation of the gas turbine. The gas turbine efficiency is defined as the formula below:where is the fuel mass flow; is the fuel low calorific value; is the turbine power.
At the same time, the surge margin should be considered in the optimization target of the optimization scheme as well, and the surge margin is defined as follows:
There are 3 stages of variable stator vanes in this paper. Considering the combination number of three optimization variables and the selection rules of optimization variables, this paper adopts the optimization scheme of traverse tree algorithm (Figure 7) to optimize the overall performance of the gas turbine under different conditions.
In the optimization scheme shown in Figure 7, the optimization process of the multistage variable stator vane opening degree combination can be seen by the process of the traversing tree. When the zerostage variable stator vanes take a value within its limits, the firststage variable stator vanes will take all the opening degree values in its range; at the same time, when the opening degree of the firststage variable stator vanes is determined, the thirdstage variable stator vanes can be selected by traversal as well. And the number of multistage variable stator vanes opening degree combinations in the final optimization scheme is the number of all leaf nodes in the “tree.”
Due to the scope of the optimization variables and the limit of selection rules, the total number of all opening degree combinations is still at an acceptable level. At the same time, considering the complexity of the relationship between the optimization variable and the optimization target, it can be found that the optimization problem has obvious unimodality in the optimization process. In each variable stator vane stage, there is always an optimization target corresponding to the opening degree in the range of its opening degree, which is the optimal value, showing a unimodal phenomenon. According to the law, we can greatly reduce the number of variable stator vane opening degree combinations and simplify the optimization process.
3.3. Optimization Results and Analysis
The lowpressure compressor variable stator blade adjustable regulation law has been optimized at different offdesign points of the gas turbine. Table 2 shows the results of the optimization.


It can be calculated in the lowpressure compressor characteristic curves of each condition by inputting optimized adjustment schemes of variable stator vane into the model, as shown in Figures 8 and 9. It can be seen from Figure 8 that the compressor surge margin that has been optimized seems more gentle compared to the original surge margin. In Figure 9, the efficiency of the lowpressure compressor has significantly improved compared to the original scheme.
The changing curves of gas turbine efficiency and fuel quantity in various conditions by applying the optimized variable stator blade adjustment scheme to the simulation model are shown in Figures 1011. It can be seen from Figures 10 and 11 that the overall efficiency was significantly higher than the original regulation scheme above 50% conditions, and the fuel mass flows are lower than the original scheme.
3.4. Results and Analysis of the SteadyState Performance of Marine Gas Turbine
The steadystate calculation of a threeshaft gas turbine for a certain type of ship’s power generation is carried out, which is under the original adjustable stator blade adjustment scheme, the optimal regulation scheme, and the no adjustment scheme. The comparison results of the various parameters of the gas turbine engine used in the marine power generation are shown in Figures 12–15.
Specific analysis of simulation results is as follows.
From Figures 1213, when the adjustable stator blade maintains the design position, the lowpressure rotor speed is low and the highpressure rotor speed is higher in each steady state. At the same time, compared to the original adjustment scheme, the optimized stator blade adjustment scheme is more different than the original method. It can reduce the consumption of fuel flow and improve the economic performance of the gas turbine by adjusting the adjustable stator blade, which is shown in Figure 14. In Figure 15, in the use of adjustable stationary blade adjustment mode, lowpressure compressor surge margin under different operating conditions is at a safe level. Through the adjustment of the adjustable stationary blades, the working condition of the compressor is in good coordination and it can enhance the stability of the gas turbine and scope of work of the compressor. As can be seen from Figure 16, the efficiency of the whole machine is greatly improved by adjusting the adjustable stator blade in each steady state of the combustion engine.
4. Conclusions
In consideration of the characteristics of the variablegeometry compressor, a characteristic calculation program of a onedimensional compressor was used in the 0D simulation model of the gas turbine through scaling simulation. In the end, the simulation model of the gas turbine with variable dimension was built. The threestage VSV of the lowpressure compressor is used as the optimization variable, the economy of the gas turbine under different conditions is used as the optimization goal, and the best adjustment angles of the threestage VSV were got under different working conditions. From the optimization results, we can find that the optimized control scheme can improve the efficiency of the gas turbine and reduce fuel consumption. Although the surge margin is less than the original program’s surge margin under high working conditions, in the low and medium conditions, the optimization scheme has more advantages. Through the reasonable adjustable stationary blades adjustment scheme, we can improve the gas turbine compressor surge margin under the working condition of low and medium. Meanwhile, we can expand the operating range of the compressor, and the gas turbine can get a better operation state. The operating performance of the gas turbine under different operating conditions is improved. It is believed that the reasonable adjustment scheme of the stator has a positive effect on the overall performance of the gas turbine.
Notation
AcronymsVSV:  Variable stator vane 
LPC:  Lowpressure compressor 
HPC:  Highpressure compressor. 
:  Enthalpy 
:  Temperature 
:  Power 
:  Mass flow 
:  Pressure 
:  Specific heat at constant pressure 
:  Gas constant 
:  Low calorific value of fuel 
:  Volume 
:  Rotor inertia 
:  Torque 
:  Power. 
:  Pressure ratio 
:  Efficiency 
:  Index adiabatic of combustion chamber 
:  Angular speed. 
in:  Inlet of compressor 
out:  Outlet of compressor 
:  Compressor 
:  Isentropy 
:  Combustion 
in:  Inlet 
out:  Outlet 
:  Turbine 
:  Fuel 
:  Load. 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The present work is supported by the Fundamental Research Funds for the Central Universities of China (no. HEUCFM170301).
References
 J. Zhang, Model Research of Engine Performance Based on Level Superposition of Compressor, Beihang University, Beijing, China, 2010.
 W. Deng, Research on Simulation Technology of Gas Turbine, Nanjing University of Aeronautics and Astronautics, Nanjing, China, 2012.
 D. P. Petters and J. L. Felder, “Engine system performance of pulse detonation concepts using the NPSS program,” in Proceedings of the 38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, July 2002. View at: Google Scholar
 R. Sampath, R. Plybon, C. Meyers, R. Irani, and M. Balasubramaniam, “High fidelity system simulation of aerospace vehicles using NPSS,” in Proceedings of the 42nd AIAA Aerospace Sciences Meeting and Exhibit, AIAA, Reno, Nevada. View at: Publisher Site  Google Scholar
 M. G. Turner, R. Ryder, J. A. Reed, and J. P. Veres, “Multifidelity simulation of a turbofan engine with results zoomed into minimaps for a ZeroD cycle simulation,” in Proceedings of the 2004 ASME Turbo Expo, pp. 219–230, ASME, Vienna, Austria, June 2004. View at: Google Scholar
 N. D. Semkin, A. V. Piyakov, K. E. Voronov, N. L. Bogoyavlenskii, and D. V. Goryunov, “A linear accelerator for simulating micrometeorites,” Instruments and Experimental Techniques, vol. 50, no. 2, pp. 275–281, 2007. View at: Publisher Site  Google Scholar
 M. T. Schobeiri, “Active aerodynamic control of multistage axial compressor instability and surge by dynamically adjusting the stator blades,” in Proceedings of the ASME Turbo Expo 2001: Power for Land, Sea, and Air, GT 2001, ASME, New Orleans, La, USA, June 2001. View at: Publisher Site  Google Scholar
 V. Cyrus, “Aerodynamic performance of an axial compressor stage with variable rotor blades and variable inlet guide vanes,” in Proceedings of the ASME International Gas Turbine and Aeroengine Congress and Exhibition, GT 1998, ASME, Stockholm, Sweden, June 1998. View at: Publisher Site  Google Scholar
 T. R. Camp and I. J. Day, “A study of spike and modal stall phenomena in a lowspeed axial compressor,” Journal of Turbomachinery, vol. 120, no. 3, pp. 393–401, 1998. View at: Publisher Site  Google Scholar
 J. E. RoyAikins, “Considerations for the use of variable geometry in gas turbines,” in Proceedings of the ASME International Gas Turbine and Aeroengine Congress and Exposition, ASME, Brussels, Belgium. View at: Publisher Site  Google Scholar
 C. Celis, P. D. M. R. Pinto, R. S. Barbosa, and S. B. Ferreira, “Modeling of variable inlet guide vanes affects on a one shaft industrial gas turbine used in a combined cycle application,” in Proceedings of the ASME Turbo Expo, pp. 1–6, ASME, Berlin, Germany, June 2008. View at: Publisher Site  Google Scholar
 A. Wiedermann, D. Frank, U. Orth, and M. Beukenberg, “Computational and experimental analysis of an industrial gas turbine compressor,” in Proceedings of the ASME Turbo Expo: Turbine Technical Conference and Exposition, GT2011, pp. 319–329, ASME, Vancouver, British Columbia, Canada, June 2011. View at: Publisher Site  Google Scholar
 Y. Hu, Study of DSP Control Technology for Inlet Guide Vane/Stator in Multistage Axial Compressor [Dissertation, thesis], School of Chinese Academt of Sciences, China, 2009.
 Y. Hu and C. Nie, “Exploration of variable working condition of guide vane/stator vane/the rotational alignment in axial compressor,” Chinese Science:Science and Technology, vol. 7, pp. 765–771, 2010. View at: Google Scholar
 J. F. Klapproth, M. L. Miller, and D. E. Parker, “Aerodynamic development and performance of the CF66/LM2500 compressor,” AIAA, 1979, 797030. View at: Google Scholar
 A. R. Wadia, D. P. Wolf, and F. G. Haaser, “Aerodynamic design and testing of an axial flow compressor with pressure ratio of 23.3:1 for the LM2500+ gas turbine,” Journal of Turbomachinery, vol. 124, no. 3, pp. 331–340, 2002. View at: Publisher Site  Google Scholar
 W. Hu, M. Kai, and S. Na, “Optimization of variable geoemetry and performance constrains of multistage axial compressor,” Aircraft Engine, vol. 38, article 15706, pp. 11–15, 2012. View at: Google Scholar
 W. Hu, “Multi objective optimization of multi stage axial compressor with variable geometry,” Journal of Engineering Thermal Physics, vol. 3512, pp. 2363–2366, 2014. View at: Google Scholar
 W. Hu and S. Na, “Optimization analysis of multi blade row adjustment and expansion of multi stage axial compressor,” Mechanics and Practice, vol. 3202, pp. 26–31, 2010. View at: Google Scholar
 W. Hu and S. Na, “Sun Na.Research on variable geometry optimization method of multi stage axial compressor,” Journal of Aerospace Power, vol. 2411, pp. 2558–2563, 2009. View at: Google Scholar
 J. Zhang and M. Ren, “Experimental study on the influence of static blade angle on compressor performance,” Journal of Aerospace Power, vol. 1, pp. 27–30, 2000. View at: Google Scholar
Copyright
Copyright © 2017 Zhitao Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.