Next Article in Journal
Dynamic Fault Tree Model of Civil Aircraft Avionics Network Transmission Failure Based on Optimized Extended Fuzzy Algorithm
Next Article in Special Issue
The Key Role of Research in Flight Dynamics, Control, and Simulation for Advancing Aeronautical Sciences
Previous Article in Journal
Simulation and Experimental Study of Gas Turbine Blade Tenon-Root Detachment on Spin Test
Previous Article in Special Issue
Investigation of Fluid Dynamics in Various Aircraft Wing Tank Designs Using 1D and CFD Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Control-Oriented System Identification of Turbojet Dynamics

by
Francisco Villarreal-Valderrama
1,
Eduardo Liceaga-Castro
2,
Diana Hernandez-Alcantara
3,
Carlos Santana-Delgado
1,
Selcuk Ekici
4 and
Luis Amezquita-Brooks
1,*
1
Facultad de Ingeniería Mécanica y Eléctrica, Universidad Autónoma de Nuevo León, San Nicolás de los Garza 66455, Mexico
2
Escuela Superior de Ingeniería Mecánica y Eléctrica—Ticomán, Instituto Politécnico Nacional, Ciudad de México 07340, Mexico
3
Escuela de Ingeniería y Tecnologías, Universidad de Monterrey, San Pedro Garza García 66238, Mexico
4
Department of Aviation, Igdir University, Iğdır 76000, Türkiye
*
Author to whom correspondence should be addressed.
Aerospace 2024, 11(8), 630; https://doi.org/10.3390/aerospace11080630
Submission received: 9 June 2024 / Revised: 18 July 2024 / Accepted: 23 July 2024 / Published: 1 August 2024
(This article belongs to the Special Issue Flight Dynamics, Control & Simulation (2nd Edition))

Abstract

:
The autonomous operation of turbojets requires reliable, accurate, and manageable dynamical models for several key processes. This article describes a practical robust method for obtaining turbojet thrust and shaft speed models from experimental data. The proposed methodology combines several data mining tools with the intention of handling typical difficulties present during experimental turbojet modeling, such as high noise levels and uncertainty in the plant dynamics. The resulting shaft speed and thrust models achieved a percentage error of 0.8561% and 3.3081%, respectively, for the whole operating range. The predictive power of the resulting models is also assessed in the frequency domain. The turbojet cut frequencies are experimentally determined and were found to match those predicted by the identified models. Finally, the proposed strategy is systematically tested with respect to popular aeroengine models, outperforming them both in the time and frequency domains. These results allow us to conclude that the proposed modeling method improves current modeling approaches in both manageability and predictive power.

1. Introduction

The propulsion systems of Unmanned Aerial Vehicles (UAVs) require as much performance and reliability (or even more) as traditional manned aircraft. This becomes clear when considering some of the challenges inherent to the operation of these vehicles, such as adverse operating conditions, multiple operating configurations, variations in the effective payload, stringent reaction-time demands, and performance variability due to component degradation [1,2]. In addition, the variables measured during flight are significantly limited [3]. This combination of operational uncertainty and limited sensor information renders the aeroengine modeling task to be particularly important yet complex.
Turbojet modeling can be useful to estimate variables that cannot be measured during flight and to obtain redundant sensor measurements. For instance, turbojet thrust, which is the ultimate purpose of aircraft propulsion systems, is typically not measured during flight. Therefore, thrust models are computed with off-board experiments and dedicated test facilities [4]. The importance of thrust models relies on their usefulness for the decision process related to aircraft speed and attitude [5]. On the other hand, the shaft speed acts as a safety constraint for turbojet operation [6], is the most commonly regulated variable in commercial set-ups [7,8], and plays an important role as a health [9] and thermal state indicator [10]. The importance of turbojet thrust and shaft speed has led to the introduction of many mathematical representations. These representations can be built either theoretically or experimentally:
  • Theoretical models use physical principles to make predictions of turbojet performance [11,12,13,14]. The underlying principles of theoretical models are based on thermodynamic analyses for mass, energy, and momentum conservation [15]. Most theoretical models yield computationally intensive solutions and/or only deal with static conditions. Thus, these models are often useful for design [16,17], assessing environmental impact [18], and long-term health monitoring [19,20].
  • Experimental models are data-based representations calculated with data measured during operation. Experimental model applications range from assessing aeroengine performance and emissions [21,22], component degradation, and environmental impact [23] to diagnosis [24,25]. Nonetheless, these models often require a priori knowledge of the model structure and larger data sets to compute the corresponding parameters. Complex model structures are prone to be more sensitive to parameter variations [26], especially when dealing with high levels of sensor noise and uncertainty in the system dynamics [27,28]. This is particularly true when considering that in practice, aircraft model parameters are computed from noisy flight test data [29,30,31].
In this context, hybrid combinations of theoretical and experimental models can be used to combine the advantages of each approach. The main difference that segregates fully experimental (black-box) and gray-box models is the intention of finding a consistent model structure rather than focusing solely on data fitting. Gray-box models usually possess a structure that combines well-known physical laws with data-based parameter adjustment [32], as depicted in Figure 1. The purpose of using physically-inspired structures is to have a model structure consistent with the real-world process; nonetheless, the recent development of data mining tools has enabled the discovery of consistent patterns and model structures underlying in experimental data [33,34]. Gray-box models obtained through data mining have already yielded good results in energy forecasting applications [35,36], thus presenting an interesting and novel approach to mathematical modeling.
In order to select the best-suited mathematical representation and system identification method, the application and available information play a key role. For turbojet autonomy, which is required to guarantee performance and operational safety [37,38], the preferred models are those that facilitate maintaining the engine within specific operational conditions and dealing with the system’s inherent nonlinear dynamics [39,40]. From a practical standpoint, turbojet autonomy also requires dynamic models that remain manageable enough to allow its evaluation with metrics typically used by aeronautical authorities, such as GSFC-STD-1000 [41] and MIL-HDBK-516B [42].
The previous discussion shows that current aeroengine models can still be improved to better comply with the specific requirements that involve pursuing turbojet autonomy. For instance, simpler linear approximations are commonly used for this application [7,39,40,43,44,45]. However, UAVs operate in wider regimes than manned aircraft [46,47], and the valid operating range of linear approximations is often exceeded. Neural-network-based models can be used to meet the nonlinear system requirements [48,49]; nonetheless, these models are incompatible with certification metrics, provide little insight into the process dynamic characteristics, and may involve a high computational cost for on-board implementation [50,51]. This explains why, even with the availability of highly accurate modeling techniques, computationally lean models are preferred in many recent aerospace applications [52].
Control-oriented models have been of particular interest to improve the reliability and performance of different applied energy systems, including wind turbines [53], variable-geometry turbochargers [54], battery management systems [55,56,57], fuel-cell systems [58], multienergy systems [59], and prediction of performance in internal combustion engines [60,61,62]. In particular, it is desirable for these models to comply with the following key characteristics: (i) accuracy for the whole operating range, (ii) manageability, and (iii) compatibility with typical certification metrics. Notwithstanding the widespread development of control-oriented models in other application areas, a comprehensive bibliographical review showed that the availability of turbojet control-oriented nonlinear models is limited. Thus, this paper introduces a novel system identification approach that takes these elements into consideration for the main variables that describe turbojet operation with respect to the input fuel flow.
The proposed identification method is formulated with the intention of handling the typical difficulties presented during data-based turbojet modeling, such as high noise-to-signal ratio and uncertainty in the plant dynamics. In particular, the method is based on a hybrid time–frequency domain regression with a clustering + classification technique, which allows obtaining improved low-order gray-box models even in these adverse conditions. The main objective of the proposed identification process is for typical linear control structures to be compatible with the resulting gray-box models with only very minor adaptations. This property simplifies the adoption of the model with the most widespread turbojet control techniques.
In addition, due to the lack of standardized testing procedures for turbojet dynamical model identification and validation, a comprehensive set of experimental tests is proposed for this purpose. The resulting models are compared with two other popular aeroengine modeling approaches in order to assess the performance of the proposed method with current methodologies. The results show that the proposed method is able to capture transient, steady state, and key frequency-domain properties of the engine shaft speed and thrust. Overall, the proposed models show improved predictive power while being structurally simpler.

2. Model Structure Definition

Although linearity can be attractive from an engineering point of view, assuming linearity during controller design can significantly reduce the operating range of the resulting schemes, which is required to take full advantage of the propulsion system’s capabilities. In the context of aeronautical systems, this has historically been dealt with by several techniques that extend the operating range of typical structures, such as gain scheduling [63]. Another option is the direct integration of more complex model structures into the controller design stage [64]. A salient characteristic of these approaches is that the resulting controllers are heavily dependent on the mathematical approach used to represent the plant’s nonlinear dynamics. This introduces several limitations, since the use of a particular model structure limits the viable options for control design and analysis.
Block-structured models are separable models assembled by linear and nonlinear subsystems. The most common configurations are the Hammerstein and Wiener structures, which consist in a cascade connection of a memoryless static nonlinear subsystem and a linear time-invariant subsystem; the difference being the arrangement of the subsystems, as shown in Figure 2.
Wiener-structured models are prone to be compatible with aeroengine modeling [65], since they allow to better represent static and dynamical uncertainties. In [66], the precision of both Wiener and Hammerstein representations is compared, concluding that the Wiener structures outperform Hammerstein structures in modeling a nonlinear process with less identification data requirements. Although previous studies have pursued aeroengine modeling through Wiener models, the present paper assesses the system identification problem through a more realistic approach by dealing with real-process uncertainties and noise through an experimental framework. Moreover, the proposed identification method uses a Wiener structure, specifically due to its particular properties regarding controller design and analysis, which are described as follows:
  • Wiener models can retain the stability, robustness, and performance properties of the linear dynamic component for each possible operating condition [67] by exploiting the properties of the linear subsystem [68].
  • Wiener models can be analyzed in terms of classical robustness margins [69].
These properties allow combining Wiener models with a wide variety of control design techniques, such as model predictive control [70], adaptive control [71], classical control [72], gain scheduling schemes [73], feed-forward strategies [74], and optimal control methods, such as H controllers [75]. Thus, the resulting models can be seamlessly integrated into previous research based on linear representations to assess the resulting system’s stability and performance characteristics. As an example, in [38], the structure of gray-box models was exploited to design an optimal direct-thrust controller with the aid of a thrust observer, resulting in reduced fuel consumption. The previous discussion shows how this structure allows complying with the specific goal of achieving model compatibility with the most widespread control design techniques.

3. Novel Data-Based Identification Method

The proposed identification method used for the turbojet model estimation is described in the following section. This is a two-step process that segregates causal and stochastic dynamical components from the system signals and thereafter captures the nonlinearities in the process through a simple memoryless mapping function.

3.1. Linear Subsystem Identification

During this stage, the purpose of the identification algorithm is to capture the different dynamical and frequency components of the signal; however, the resulting parameters are influenced by measurement noise and other types of distortion. In experimental conditions, the measured response of the system encompasses:
  • Causal components (from input variables) that represent the true dynamic behavior of the system (a consistent behavior in the response).
  • Stochastic components can be considered to be contamination external to the process.
This introduces the necessity of segregating these causal and stochastic components. For instance, it has been shown that high levels of non-Gaussian noise can mistakenly be identified as additional zeros–poles, a phenomenon akin to overlearning in other data-based identification methods. The proposed method is designed to overcome the difficulties associated with dynamic subsystem identification. This method departs from a well-known approach to obtain a discrete-time-invariant model, represented by transfer function G(z), as shown next:
G ( z ) = v ( z ) u ( z ) = i = 0 n b b i z i j = 0 n a a j z j
where v is the linear output to an arbitrary input u, and a i and b i are the coefficients of the discrete transfer function. Therefore, the number of poles and zeros of the linear subsystem is given by n a and n b , respectively. Thus, the resulting sampled time response v ( k ) of linear Model (1), presented as a difference equation, is
v ( k ) = a 1 v ( k 1 ) a n a v ( z k n a ) + b 1 u ( k 1 ) + b n b u ( k n b ) + e ( k )
where k is the sample number. The coefficients of the polynomial function can be easily computed through parameter optimization by the minimization of modeling errors σ of n output measurements y m :
σ 2 = k = 1 n ( y m ( k ) v ( k ) ) 2
If θ m i n ( | | σ | | ) , least squares algorithm can used to determine transfer Function (1) parameters:
θ = ( ϕ T ϕ ) 1 ϕ T Y
where θ contains the transfer Function (1) parameters θ = [ a 1 a n a b 1 b n b ] T that correspond to its poles and zeros polynomial coefficients, and ϕ is the regression vector ϕ = [ v ( k 1 ) v ( k n a ) u ( k ) u ( k 1 ) u ( k n b ) ] . This structure associates the known measured values of past states and input conditions with the unknown parameters. Although least squares algorithms can provide a good level of noise rejection allowing the identification of underlying dynamical components, their performance can be limited when dealing with nonlinear systems, time-varying systems, and high levels of non-Gaussian noise [76]. Furthermore, least squares algorithms require a priori knowledge of the zero–pole structure (i.e., n a and n b ). An initial idea of the zero–pole structure can be obtained from theoretical analysis; however, complex models, such as turbojets, require one to account for possible additional unmodeled dynamics. This yields a two-pronged problem with respect to the previous formulation: (1) determining the ideal zero–pole structure and (2) discriminating valid and invalid zeros–poles (i.e., reducing overlearning). This paper introduces a practical and novel approach to these problems using data mining to analyze and classify pattern formations. The proposed method to obtain the linear subsystem consists in the following:
  • Defining a set of j transfer functions, each one obtained from a small portion of the operating range using the aforementioned least squares approach. The complete set of transfer functions must cover the whole operating range. In addition, each one of these transfer functions must be estimated considering a large number of poles and zeros, generating a set of high-order identified linear models. If a large enough zero–pole structure is used, these models must then contain both valid and invalid zeros–poles. Arguably, the valid components must be consistent along all the operating ranges (with small variations due to changes in operating conditions) [77]. On the other hand, the invalid components introduced by noise, overlearning, and other stochastic elements must be of a more random nature, that is, inconsistent along the set of identified linear models. Moreover, the fundamental frequency of the identified components also reveals important information regarding their validity. For instance, turbojet thrust measurements are normally severely contaminated by high-frequency vibration noise. However, it is well-known that thrust dynamics are much slower.
  • Segregating valid and invalid components of the identified models using the following considerations: (1) repeatability and (2) fundamental frequency. This can be easily achieved by superimposing the zeros–poles of identified transfer functions in a polar plane and looking for clusters. The clusters must be defined according to the number of poles–zeros within a radius determined by the experiments. The threshold defined by the circle boundaries allows determining the zero–pole values that yield a low standard deviation with respect to the circle centroid. That is, the estimates of the zero–pole locations of the system are bounded within conditions that ensure consistency and repeatability in the system time and frequency response characteristics. Inside the circumscribed area there must exist a pole–zero for each identified transfer function to ensure repeatability. Note that increasing the circle radius yields more variety in the identified system dynamics characteristics and may integrate stochastic dynamics that are not part of the system dynamics. The resulting plot allows observing dynamics that are being consistently identified in all the estimated models as clusters, while those that lack a cluster pattern correspond to the stochastic components. This pole–zero clustering analysis allows for distinguishing the causal and stochastic components present in the measurements and determining the appropriate pole–zero structure at the same time. Moreover, the frequency of the identified components is also easily assessed by the angle of the identified cluster. In the case of shaft speed and thrust dynamics, clusters near the real axis, being of the lowest frequency, are preferred (i.e., dominant in the steady state).
  • Forming a transfer function with poles–zeros equal to the average of each cluster. That is, each cluster is translated into a single pole–zero. This yields a time-invariant low-order transfer function with the averaged dynamical properties found to be consistent in the complete operating range.
It is clear that the proposed method depends on the clustering formation, which must be adjusted according to the experiment, mainly the size and number of elements are relevant parameters (i.e., dispersion and aggregation conditions). In particular, the number of elements should be equal or close to the number of identified linear models (i.e., operating points) to ensure that the dynamic modes are effectively present in all cases. This aspect defines the main aggregation condition for the clustering. On the other hand, cluster size can be determined experimentally. A smaller cluster size indicates more homogeneity of the dynamic modes along the complete operating range and is desirable. A larger cluster size could indicate either a higher level of nonlinearity in the experimental dynamic response or a higher noise level. That is, the dispersion of the clusters is related to these aspects. In practice, it was observed that the method lends well to a visual adjustment of this variable since a 2D polar plot is always obtained no matter the number of operating points. For instance, if the Euclidean distance is used to define the clusters, the poles and zeros would have to remain within a distance d c circumscribed with respect to the centroid c c . That is, a cluster of poles exists if and only if
c c p j < d c j
where p j belongs to the roots of b 0 + b 1 z 1 + b n b z n b of each j t h identified transfer function. Respectively, a cluster of zeros exists if and only if
c c z j < d c j
where z i belongs to the roots of a 0 + a 1 z 1 + a n a z n a for each j t h identified transfer function. The automated adjustment of cluster properties, such as k-means [78] and hierarchical clustering [79], can also be used, as long as the specific goal of clearly distinguishing between consistent and stochastic system dynamics is met.

3.2. Static Function Nonlinear Approximation

The typical operation regimes, which involve operating at steady state conditions, introduce the necessity of accurate prediction of the static characteristics. If there exists a family of equilibrium conditions parameterized by the output and each of these conditions represents a feasible operating condition, then an operator can be designed for each local time-invariant approximation. Hereby, gain scheduling provides a simple approach to nonlinear analysis by interpolating between the members of each family. Thus, the nonlinear static function output, defined as y = f ( v ) , can be obtained through any nonlinear mapping. Consider a ratio r ( y m , v ) among the identified linear model response v ( k ) and the measurements y m ( k ) in the form of
r ( y m , v ) = y m ( k ) v ( k )
This ratio represents the remaining nonlinearities that were not captured by the linear model alone, ideally with y m ( k ) = r ( y m , v ) v ( k ) . If f ( v ) f ( v ) r ( y m , v ) v ( k ) , it is possible to define and minimize a residual δ for a given experiment length k = 1 n as
δ 2 = k = 1 n ( y m ( k ) f ( v ( k ) ) ) 2
Then, the resulting model estimates y ( k ) = f ( v ( k ) ) are close to the experimental data available. That is, | | δ | | 0 with the nonlinear static gain approximation. On the other hand, the compatibility of the resulting model with specific control design tools depends on the properties of the nonlinear output function.
With this in mind, simpler nonlinear approximations are preferred for this application. The studies presented in [15,80] analyzed several Exhaust Gas Temperature (EGT) models, providing an experimental criterion for model selection based on the application and the available input data. The results suggest that, for this particular application, static microturbojet nonlinearities can be approximated with polynomials. These structures have also presented good results at capturing the stationary properties of the shaft speed of large- and small-scale turbojets [81,82]. In addition, one such approximation was experimentally tested in [83] with the thrust of a small-scale turbojet. Following these results, in this article, polynomial approximations of order i (which is determined through exploratory data analysis) are used to model the shaft speed and thrust nonlinearities:
r ( y m , v ) v ( k ) f ( v ) = β 0 + β 1 v ( k ) β i v ( k ) i
with a vector of polynomial coefficients β = [ β 0 β 1 β i ] T . Since different nonlinear static gain approximations are possible, it is convenient to define the optimization problem in a generalized form. For instance, gradient descent can be applied in this case to compute the coefficients of Approximation (9) with l iterations. That is,
β l : = β l 1 γ k = 1 n f ( v ( k ) ) y m ( k ) v ( k )
where γ is the learning rate. The goal of this algorithm is to find β m i n | | δ | | , which implies fitting the nonlinear static gain Approximation (9) to the Ratio (7). The role of the nonlinear static gain is to narrow the gap between the linear model estimates and the real process behavior.

3.3. Identification Method Summary

Figure 3 presents a flowchart of the novel identification method. Additionally, a summary of this method is presented next.
  • High-order linear dynamic subsystem set identification: A set of high-order transfer functions is obtained. Each transfer function is obtained from a small portion of the operating range. The parameters of each transfer function are determined through a squared error minimization problem.
  • Cluster analysis: The parameter identification results are superimposed in a polar plot, in the form of the poles and zeros of the identified transfer functions.
  • Model simplification: Pole–zero clusters are formed and identified according to the number of elements and their proximity, the dynamics being identified consistently are considered to represent the causal components of the input.
  • Nonlinear approximation: A nonlinear gain is characterized through a polynomial function with the information from the different steady-state values.
  • Final model determination: The time-invariant transfer function and the nonlinear static gain function are assembled in a Wiener structure.
Note that the resulting models are as good as the input data used for identification. Therefore, different techniques of data reconciliation are advised as a preliminary step to perform a consistency check on the measurements, detect errors, and validate that the data are useful.

4. Experimental Set-Up

In this article, the SR-30 gas turbine from Turbine Technologies (Chetek, WI, USA) was used as an experimental case study involving the modeling of a commercial microturbojet. The obtained information allows assessing the performance of the proposed gray-box approach in an experimental framework.

4.1. SR-30 Engine

The microturbojet comprises an inlet, a centrifugal compressor, an annular combustion chamber, a single-stage axial flow turbine, and a subsonic nozzle (Figure 4). The maximum operational rotational speed is 87,000 RPM with a maximum thrust of 178 N. The fuel flow is adjusted by a manual lever, which is connected to a valve that controls the flow. The fuel blend consists of Jet A with a specific energy of 42.80 MJ/kg. The thrust force is measured through a strain gauge cell, while a pressure transducer system is used to measure the fuel flow and a 2-pole tachometer generator measures the rotational shaft speed.

4.2. Data Acquisition System

The platform includes the PCI 4351 data acquisition system; an A/D board with 24-bit resolution for 16 analog inputs with a 60 samples/s capability and a Virtual Bench Logger data acquisition program for real-time monitoring of the measured variables on a host PC. The maximum practical sampling frequency (8.85 Hz) was selected for all the experiments.

4.3. Design of Experiments

Attempting to model with poorly informative data and/or using an inappropriate model structure are the most common errors during system identification [65]. Therefore, the identification and validation data sets must contain enough information regarding the system properties. Although there are well-known testing regimes used for static turbojet modeling, used for example to calculate performance maps [84], this is not the case for dynamic turbojet identification. Therefore, the following identification and validation tests for microturbojet modeling are proposed to help standardize microturbojet dynamic model identification and validation.
Three separate tests are specified. These tests are designed to provide appropriate data sets for model identification and extensive evaluation of the resulting models. The main objectives and practical considerations for each test are described next and presented in Figure 5.
  • Identification regime. For the excitation signal, stepped inputs of varying duration and magnitude are preferred in this application since they allow extracting information of the dynamic and static characteristics of the engine [85] and are required to test the repeatability as stated by the identification method. The suggested identification procedure is presented in Figure 5a. The measured data are later divided according to each step (labeled as A, B, C, D, E, F, and H). The experiment must include the whole operating range and each step should be maintained sufficient time for the system to stabilize. The proposed experimental design weights the mid-lower shaft speed ranges higher than the top range (only one step in the upper range is used), because this is the most common operating region. The experimental design (i.e., number of steps and operating range) can be modified to weight any desired region according to particular applications.
  • Time-domain validation. To validate the resulting models a validation data set must be obtained and the resulting models should perform well under these untrained test conditions. In this case, both the transient and steady characteristics of the turbojet must be excited. The excitation signals suggested for this purpose must contain a variety of steady-state conditions and dynamical behaviors, simulating a wide range of real operating conditions. In this case, a series of fast steps, ramp signals, sinusoidal inputs, and whole-rage step signals with high frequency content are suggested. This combination of signals allows testing for transient, steady state as well as slower thermodynamic and other nonmodeled phenomena. These excitation signals can be applied in separate experiments or combined in a single experiment as the one shown in Figure 5b.
  • Frequency-domain validation. Finally, the frequency-domain properties of the plant are fundamental for the development of model-based controllers and estimators because they are responsible for key stability and robustness properties. The predictive power of the resulting models in the context of the frequency domain can be assessed with a chirp-like input as the one shown in Figure 5c. The frequency range of the input signal depends on the turbojet capabilities, with the upper bound being the most important. However, it is easy to find this upper bound experimentally by increasing the input frequency until the measured output variables reach a negligible amplitude. For instance, for the tested microturbojet the frequency range obtained using this procedure was 0.574 rad/s–13.78 rad/s.
In all cases, relevant environmental variables must be registered during the tests, including temperature, pressure, and relative humidity. The tests include steady flow at the engine inlet (no induced airspeed) and constant altitude (390 m asl) conditions.

5. Applied System Identification

In this section, the identification of a microturbojet considering the model structure and the method described in the previous sections is presented. In addition, two baseline models are incorporated in this section for further assessment of the model performance.

5.1. Shaft Speed Identification

The proposed high-order transfer functions comprise 19 poles and 3 zeros before the cluster analysis, with a total of 7 identified transfer functions, one for each step input of the identification regime. The identified zeros and poles of each transfer function are superimposed in Figure 6a. The clusters were identified according to the amount of elements (one for each identified model to ensure repeatability) and their proximity. The shaft speed pole clusters are located near 0 and ± 25 deg/s. The high-frequency pole clusters near the Nyquist–Shannon frequency limit are not included in the reduced model to reduce the effects of sampling distortions. The rest of the poles and zeros are randomly scattered over the Z plane. Although some elements resemble clusters, the number of elements condition was not satisfied. That is, these were not identified consistently in all operating conditions and are of higher frequency, suggesting that these components are due to over-learning. Considering the previous procedure, the reduced shaft speed model yields
G s ( z ) = 82.70 ( z 0.936 ) ( z 2 1.608 z + 0.692 )
Figure 6b shows the frequency response of the resulting reduced Model (11) in comparison with the high-order models. This figure allows to evaluate the effects of the cluster simplification.
The static function approximation is obtained considering the ratio among the prediction of the reduced linear dynamic model and each steady state from Figure 6a. The shaft speed static function ( f s ) was approximated through the following polynomial:
f s = 2.878 × 10 4 + 4.021 × 10 1 v s + 1.532 × 10 6 v s 2
Figure 7b presents the accuracy of the nonlinear function estimation in relation with the experimental data.
The complete shaft speed nonlinear model is arranged in the structure shown in Figure 2. This representation consists in transfer function G s ( z ) in series with static function f s . Figure 7a presents the performance of the linear and nonlinear models in comparison with the experimental data. This figure shows that although the dynamic components may be well represented by the reduced linear subsystem; the static gain can diverge depending on the operation condition. This is especially true for the shaft speed dynamics, where the nonlinear behavior is more dominant.

5.2. Thrust Identification

Thrust modeling presents different challenges. In contrast to the measured shaft speed data, thrust measurements contain higher levels of noise, while the nonlinear behavior is less dominant. The proposed high-order transfer functions consist in 21 poles and 19 zeros before the cluster analysis. The zeros and poles of all transfer functions are presented in Figure 8a. This figure shows that the main thrust pole-cluster is located near 0 deg/s. The rest of the poles and zeros appear to be randomly scattered over the polar plane, showing that the clustering method allows extracting the main transient characteristics of the turbojet under the uncertainty and noise conditions commonly found in these devices. Again, although some elements resemble clusters, the number of elements condition was not satisfied.
Figure 8b shows that the reduced model, obtained through the cluster analysis, also captures the averaged frequency–response of the high-order models up to the frequency where they are consistent. The resulting transfer function is
G t ( z ) = 0.7962 z 0.936
The nonlinear static function approximation was calculated similarly to the shaft speed presented before, resulting in the nonlinear output Function (14).
f t = 3.881 × 10 1 1.004 v t + 2.356 × 10 2 v t 2
The resulting approximation and the fitted data are presented in Figure 9b. The final thrust model has the structure presented in Figure 2 with a linear time-invariant model G t ( z ) and the memoryless nonlinear model f t . Lastly, Figure 9a presents the simulation results of the complete thrust model compared with experimental data. These results show that the proposed gray-box identification approach allows accurate modeling of the thrust behavior even in the presence of high levels of measurement noise.

5.3. Baseline Models

In general, assessing the performance of a particular model requires comparing the quality of the model to some relevant baseline. With this purpose, Hammerstein–Wiener and Neural Network AutoRegressive with eXogenous input (NNARX) models were obtained using a commercial parameter-computing tool for both the shaft speed and thrust. These models were obtained using the same identification regime data set.

5.3.1. Hammerstein–Wiener Model

This methodology has been used successfully in commercial applications by computing the parameters of each stage from measured data using least-squares estimation [86]. The accessibility of this approach has led to widespread use in different engineering applications with control and estimation purposes, such as modeling borefield thermal dynamics [87] and two-stage shape memory alloys [88]. Using this approach, Models (15) and (16) were computed for the shaft speed and thrust, respectively, using the same data set from the identification regime test and the structure presented in Figure 10.
The same pole–zero configuration obtained in the previous section was maintained. However, in this case, a more classical Hammerstein–Wiener identification method (as implemented in popular commercial toolboxes) was used mainly for comparison purposes. This methodology does not include the high-order clustering, and simplification steps proposed in this article. Therefore, comparing the resulting models will highlight the importance of these steps. The shaft speed Hammerstein–Wiener model is given by
u ^ s = 0.007 u 2 + 0.873 u 0.7745 G s h ( z ) = z ( 0.026 + z ) ( z 0.786 ) ( z 2 1.043 z + 0.725 ) y ^ s h = 1.526 v s h 2 + 1.517 × 10 3 v s h + 2.990 × 10 4
Similarly, the thrust model is
u ^ t = 0.038 u 2 + 0.404 u 1.357 G t h ( z ) = 1 ( z 0.773 ) y ^ t h = 3.118 × 10 1 v t h 2 4.0531 v t h + 34.823

5.3.2. Neural Network AutoRegressive with Exogenous Input

Dynamic neural network models are widespread and have been used in applications such as the control of heat exchangers [89] and hydro-power plants [90]. Although few reports deal with the use of NN to capture the transient characteristics of turbojets, there are several studies where NNs have been implemented to estimate different static conditions in turbojets [5,15,91]. This shows an interesting opportunity to test the neural network capabilities in estimating the transient characteristics of turbojets. For training purposes, an NN with the structure presented in Figure 11a is used. This allows training the network with the true measured output values, achieving a more accurate representation of the turbojet.
The neural network consists of a feed-forward neural network with one hidden layer composed of 10 neurons, as suggested in [15,91]. The inputs of the NN are the fuel flow and three delayed outputs to maintain an agreement with the third-order dynamics identified during the linear model identification from the previous sections. This structure is congruent with the results presented in [49,92]. These results were obtained through a sensitivity analysis of the network structure (i.e., different number of layers, neurons, activation functions, etc.) for turboshaft dynamic modeling. The best-fitting neural network was selected after an iterative process considering different initial conditions with parameters computed through Bayesian regularization. After the network identification, the neural network is manipulated into an Output Error (OE) structure (Figure 11b), where the delayed estimation is operated in a closed loop. The resulting NNARX allows addressing the estimation problem in a practical sense, where the output is not measurable. This is a representative scenario particularly for thrust, since this variable can be measured in ground stations but is not available during flight. The resulting NNs for both shaft speed and thrust are included as a novel baseline for evaluating the developed gray-box models.

5.4. Models with the Identification Data Set

The resulting models obtained through the identification methods presented in this section are shown in Figure 12a,b for shaft speed and thrust, respectively. Regarding the shaft speed, it is clear that the double nonlinearity of the Hammerstein–Wiener model captures the lower steady-state conditions with greater accuracy than the proposed gray-box model with a single nonlinearity; while the NN clearly outperforms both in terms of steady-state conditions. In transient conditions, the gray-box model provided greater accuracy than the Hammerstein–Wiener model, while the NN provided the slowest response. The resulting models during the thrust identification performed similarly to the shaft speed identification, although the NN was more heavily influenced by the measurement noise.

6. Experimental Validation

The models are evaluated using extensive validation conditions, which are divided into time domain and frequency domain. It is important to note that this additional experimental data was not used to refine the existing models. Although the scope of the paper is to develop control-oriented models rather than evaluate particular parameters of the engine, it is important to note that the resulting measurements (i.e., thrust and shaft speed), specific fuel consumption (SFC) and overall efficiency are within the same range of values as those reported in studies using the same engine, that is within 0–87,000 RPM, 20 N–178 N and up to 1.2 lb/(lb h) for each one of these parameters [8,93,94].

6.1. Time Domain

The time-domain validation results are presented in Figure 13. For the shaft speed (Figure 13a), the Hammerstein–Wiener model tends to have large errors during transient conditions. This suggests that while the nonlinear characteristics are accurately identified, the transient information is not well captured. The NN, on the other hand; struggles to provide the same performance when structured in the OE form, which is necessary for some practical applications. This loss of accuracy when restructuring the model to its OE form has been reported in previous studies [95]. The gray-box model provides a response close to the actual measured data during all the test conditions.
The major difficulty presented during the thrust modeling was handling the effects of measurement noise. The performance of the models with the thrust validation data is similar, with the main differences being that the NN shows more consistent results and the Hammerstein–Wiener model filters high-frequency components of the signal. These validation experiments confirm that the gray-box models are able to predict the microturbojet dynamic and static behaviors with a high level of accuracy, reliability, and simplicity.

6.2. Frequency Domain

An innovative framework for the evaluation of the resulting models based on the frequency domain properties of the system is presented in this section. Figure 14 presents the results of an experimental frequency-swept test of the microturbojet. These results can be compared with the theoretical cut-off frequency of the resulting models. The results show that the gray-box model successfully captures the shaft speed and thrust dynamic characteristics. The experimental amplitude decay observed in Figure 14 corresponds to the behavior of the identified models right after the cut-off frequency. For the shaft speed, the NN response does not follow the sine wave accurately, while the Hammerstein–Wiener model overpredicts the shaft speed values and presents an out-of-phase response. For the thrust, both the NN and Hammerstein–Wiener models capture the low-frequency components of the measured response, but the higher frequencies are amplified. This indicates that the frequency-domain characteristics and the order of the system have been identified better by the gray-box models up to the test bench capabilities.
While some modeling approaches are based solely on frequency response identification [93,96], in this article, frequency-domain elements were indirectly considered during the clustering stage of the identification process. In this section, an experimental frequency response analysis is used for validation. This validation is important to ensure that the models can be used properly with frequency domain controller and estimator design techniques [97].
The results are also comparable with other existing literature reports. For instance, in [98] a microturbojet engine was modeled using a single linear third-order state-space representation with good experimental results. However, the use of a fixed linear representation yields a very limited operating range, as noted in [98]. In addition, instead of adjusting the model parameters using a set of three experimental measurements (in this case, fuel flow, shaft speed, and thrust), in [98], very precise and extensive knowledge of the machine parameters and variables is necessary. In particular, the method proposed in [98] requires the measurement, calculation, reconstruction, or previous knowledge of over 30 machine parameters, which can be difficult to achieve in many setups. In another report [10], a microturbojet is modeled as a single linear discrete transfer function calculated from averaging a set of linear models obtained at different operating points. The average model is then used with a combination of robust control design techniques to account for the model uncertainty. While the approach presented in [10] yields good results, it requires special control techniques to account for the model uncertainty due to the use of a single average model. Instead, the proposed approach achieves a low level of modeling error with a single model. In another study [99], the authors use recurrent neural networks to model the dynamic response of the fuel flow of a microturbojet with good results. Although limited to this variable (which in our case is assumed to be measured), the results shown in [99] suggest that input valve fuel flow may be easier to model than the variables studied here. Finally, an older study using the same engine was presented in [8]; in this study, the authors model the shaft speed using a linear autoregressive model and a neural network, with the neural network yielding the best fit. This is in line with our results, since a single linear model fails to capture the complete operating range, and while a neural network model can achieve a better fit, its modeling performance lags behind both the HW model and the proposed gray-box model. Finally, it is important to note that in these literature reports (and in others not discussed here due to lack of space), several features of the proposed method are lacking, mainly (1) a wide set of identification and validation tests (i.e., not being the same kind of maneuvers), (2) validation tests including frequency-domain testing, and (3) a proposal for obtaining an average linear representation based on overfitting, clustering simplification, and frequency-domain considerations instead of simple averaging. Overall, the results presented in this section, particularly the benchmark models (i.e., the HW and the NN models), are in line with the literature reports, thus validating their use for the assessment of the proposed gray-box model.

7. Predictive Power Analysis

The following section presents the performance assessment of the developed models using different accuracy and precision metrics, as well as some important remarks regarding the computed models.

7.1. Models Error Analysis

In this section, a brief analysis of the model estimation error is presented. Figure 15 presents the Box Plot of the modeling error of each model during the time-domain validation. In the case of the gray-box model, the nonlinear shaft speed modeling error is bounded within a band of ± 5 % error, while the thrust nonlinear modeling error is delimited among ± 12 % (this is highly influenced by the noise in the measurements). This benchmark allows for evaluating the propagation error of the sensors’ uncertainty to the model outputs. The NNs follow in precision, with the Hammerstein–Wiener models having the greater error spread with the validation data.
On the other hand, Figure 16 presents the statistical analysis of the modeling errors during the frequency-domain validation. This validation data set is particularly characterized by the absence of steady states and wide dynamic conditions. Thus, it is designed to reveal identification errors in the process dynamics rather than in the static gain. This analysis shows that the gray-box model retains the accuracy and precision properties presented during the time-domain validation (i.e., similar characteristics in the error Box Plot during the time-domain and frequency-domain validations). This contrasts with the NN and Hammerstein–Wiener models, which presented an increased error distribution during this test. The previous observations suggest that the estimation errors of the gray-box are not strongly related to a specific dynamic condition (static or dynamic), implying that the proposed identification method captures the turbojet dynamics more accurately.
Additionally, the accuracy and precision of the models are studied using the Root Mean Squared Percentage Error (RMSPE) and Standard Error (SE), respectively. These metrics are presented in Table 1 with both validation data sets. The RMSPE is a classical error indicator included to analyze the models’ accuracy (i.e., low error percentage) even with high noise-to-signal ratios. In this case, the gray-box models provide the highest accuracy. On the other hand, the SE in combination with Figure 15 allows one to conclude that the gray-box models are also highly precise (i.e., bounded error in a small region with low variability with respect to the experiments). When comparing similar test conditions, the resulting models can match complex estimation architectures in terms of accuracy and robustness [22]. In summary, the models obtained with the proposed identification method provide greater precision and accuracy in untrained conditions when compared with other popular modeling alternatives, as well as being more robust to propagation errors while also having a simpler structure.

7.2. Error Propagation Analysis

In addition to the total error assessment, the error propagation properties of each model are studied next. For this application, the main source of uncertainty for the models is the sensor noise. The experimental sensor noise was determined by operating the engine in a steady-state regime so that the transient dynamics were minimized. Then, 300 data points of each variable were obtained and fitted to a normal distribution by optimizing the likelihood [100].
Figure 17 shows the resulting histogram of the measured data and the fitted distributions. The resulting log-likelihoods were 1359.71, −1688.49, and −679.634 for the fuel flow, shaft speed, and thrust, respectively. These log-likelihoods were compared with several other fitted distributions, including Nakami, uniform, inverse Gaussian, and extreme value. The normal distribution proved to be either the best fit or very close to the best; thus, it was decided to move forward using this distribution due to its better-known properties. The resulting standard deviations of the fitted distributions were 0.00327, 50.7, and 2.08 for fuel flow, shaft speed, and thrust, respectively. The mean values for each variable were 5.43, 73,849.6, and 87.79 for fuel flow, shaft speed, and thrust, respectively. This demonstrates that both shaft speed and fuel flow have a low level of measurement noise, while thrust has the highest level of noise. This is normal since thrust is measured using a load cell, which is well known to suffer from higher levels of noise.
Although error propagation in static models is a well-established topic, in the case of dynamic autoregressive (particularly nonlinear autoregressive models), the error propagation properties can be more involved to determine (see, for instance, [101]). On the other hand, a simple tool that allows numerically assessing the error propagation is a Monte Carlo (MC) simulation [102]. Thus, an MC simulation is presented next to assess the error propagation. The MC simulation was performed by comparing the output of each model to a noiseless constant fuel flow input from a zero initial condition with a noisy fuel flow input for a simulation time of 20 s (enough time so that all models reach the steady state). This allowed for assessing both the transient and stationary properties of each model. For each simulation, the noise added to the input fuel flow was obtained using the experimentally fitted distribution of Figure 17. A set of 100 simulations for each variable and each model were performed using this methodology. The difference between the nominal (i.e., using the noiseless fuel flow input) and the noisy outputs (i.e., using the noisy fuel flow input) was calculated and aggregated for all the simulations to evaluate the error propagation of each model.
Finally, Figure 18 shows the aggregated error histogram of each model together with the best-fit normal distribution for the thrust and the shaft speed. In the case of the shaft speed obtained with the neural-network-based model, the maximum-likelihood normal distribution does not represent the data properly; thus, a logistic and a stable distribution were also fitted, with the stable distribution resulting in the best log-likelihood. This could be the result of the nonlinear activation functions within the neural network, which tend to distort the normal distribution from the noise source.
The resulting standard deviations of the best-fit normal distributions for the shaft speed were 6.035, 14.08, and 26.5 for the gray-box, the HW, and the NN models, respectively. On the other hand, the resulting standard deviations for the thrust were 0.0136, 0.0352, and 0.1064 for the gray-box, the HW, and the NN models, respectively.
One of the main observations regarding the MC simulation is that the gray-box model yields the lowest error propagation among the compared models (as measured by the standard deviation), while the neural network model has the highest. Another important observation is that all the models have a smaller standard deviation than that of the experimentally measured signals (i.e., 50.7 and 2.08 for the shaft speed and thrust, respectively). This can be explained by observing that, although all models amplify the noise from the noisy fuel flow signal, this particular measurement has the lowest noise level of all the experimental measurements (i.e., 0.00327); thus, the error propagation from the noisy fuel flow yields a lower noise level than the experimental measurements of the output variables. It is important to note that this only quantifies the error propagation from the noisy input signal. Other sources of error and uncertainty are not included and are responsible for the total modeling error observed in the previous section (see Figure 15).
As a final note, the resulting error propagation characteristics can be considered a salient feature of the proposed gray-box model. In addition to yielding the lowest total error and the lowest input-signal error propagation, the gray-box model is suitable for thrust-control applications. In many of these applications, thrust measurements suffer from a high level of sensor noise (as noted in the previous paragraphs) and are difficult/expensive to implement. Therefore, a particularly attractive configuration is that of using a state observer for reconstructing the thrust by measuring the fuel flow as the input signal and the shaft speed as the output signal (see [38] for an example of this configuration). This application demands the lowest level of sensor noise in the fuel flow and shaft speed, and a low level of error propagation from the fuel flow to both thrust and shaft speed, which is the particular combination of characteristics that the gray-box model yields.

7.3. Remarks Regarding Microturbojet Nonlinear Behavior

Finally, important remarks regarding the cluster analysis and the static function approximation of the gray-box identification method are presented. During the cluster analysis, the model poles were located within a small region, allowing their location to vary slightly from one identification condition to another. These variations can be induced due to changes in the engine dynamics, which are dependent on the engine’s thermal state (at any given operating condition). The latter effect has been observed previously for larger turbojets [103,104]. On the other hand, this document presents, an experimental evaluation of this nonlinear behavior in microturbojets. It is important to note that the static function obtained through the results of the proposed identification process is similar in structure to those used for larger engines. The importance of this finding lies in the fact that this structural similitude has not been found using thermodynamic modeling approaches. That is, thermodynamic models for large turbojets have struggled to capture microturbojet static and dynamic behavior [105,106], but the present results demonstrate that there are indeed similarities.

8. Conclusions

A novel model for the main variables that describe the operation of turbojets is presented in this article. The resulting models are compatible with typical controller, observer, and estimator topologies. The identification methodology is formulated with the intention of providing an identification approach capable of handling the typical difficulties presented during data-based turbojet modeling, such as high noise-to-signal ratio and uncertainty in the plant dynamics.
Experimental results show that the proposed method is capable of segregating stochastic and causal components present in the measurements of a microturbojet in an experimental setup. This identification method has the following stages (1) high-order linear model identification, (2) cluster analysis, (3) dynamic model simplification, (4) nonlinear static gain approximation, and (5) final model integration using a Wiener structure. The resulting model encompasses information on the whole turbojet operation regime, allowing one to increase the operating range of linear control schemes.
The experimental validation of the models included a wide spectrum of operating regimes to test their accuracy and precision with untrained conditions. The developed models were also systematically tested along with other widespread models to identify the best-fitting aeroengine representation. These widespread models included Hammerstein–Wiener and Neural Network AutoRegressive with Exogenous Input structures. Along the experimental tests, the gray-box models performed remarkably well under several error metrics, including their error propagation characteristics, outperforming the baseline models. Moreover, the main frequency-domain properties of the turbojet were obtained experimentally and were found to closely match those predicted by the proposed gray-box models. This contrasts with the failure of the baseline models in predicting these important features. These results support the conclusion that the proposed gray-box models have successfully captured the main dynamic and nonlinear static features required for high-performance and robust turbojet control while also maintaining the lowest error propagation among the compared models.

Author Contributions

Conceptualization, S.E. and L.A.-B.; Formal analysis, F.V.-V.; Investigation, F.V.-V. and C.S.-D.; Methodology, E.L.-C. and L.A.-B.; Project administration, L.A.-B.; Resources, C.S.-D.; Software, D.H.-A.; Supervision, E.L.-C. and D.H.-A.; Validation, S.E.; Visualization, D.H.-A.; Writing—original draft, L.A.-B.; Writing—review & editing, E.L.-C. and S.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Some or all data, models, or codes that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Pisano, W.; Lawrence, D. Control limitations of small unmanned aerial vehicles in turbulent environments. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, Chicago, IL, USA, 10–13 August 2009; p. 5909. [Google Scholar]
  2. Villarreal-Valderrama, F.; Takano, L.; Liceaga-Castro, E.; Hernandez-Alcantara, D.; Zambrano-Robledo, P.; Amezquita-Brooks, L. An integral approach for aircraft pitch control and instrumentation in a wind-tunnel. Aircr. Eng. Aerosp. Technol. 2020, 92, 1111–1123. [Google Scholar] [CrossRef]
  3. Léonard, O.; Borguet, S.; Dewallef, P. Adaptive estimation algorithm for aircraft engine performance monitoring. J. Propuls. Power 2008, 24, 763–769. [Google Scholar] [CrossRef]
  4. Przysowa, R.; Gawron, B.; Białecki, T.; Łegowik, A.; Merkisz, J.; Jasiński, R. Performance and emissions of a microturbine and turbofan powered by alternative fuels. Aerospace 2021, 8, 25. [Google Scholar] [CrossRef]
  5. Dalkiran, F.Y.; Toraman, M. Predicting thrust of aircraft using artificial neural networks. Aircr. Eng. Aerosp. Technol. 2021, 93, 35–41. [Google Scholar] [CrossRef]
  6. Soares, C. Gas Turbines: A Handbook of Air, Land and Sea Applications; Elsevier: Amsterdam, The Netherlands, 2011. [Google Scholar]
  7. Surendran, S.; Chandrawanshi, R.; Kulkarni, S.; Bhartiya, S.; Nataraj, P.S.; Sampath, S. Model predictive control of a laboratory gas turbine. In Proceedings of the 2016 Indian Control Conference (ICC), Hyderabad, India, 4–6 January 2016; pp. 79–84. [Google Scholar] [CrossRef]
  8. Aly, A.; Atia, I. Neural modeling and predictive control of a small turbojet engine (sr-30). In Proceedings of the 10th International Energy Conversion Engineering Conference, Atlanta, GA, USA, 30 July–1 August 2012. [Google Scholar] [CrossRef]
  9. Jaw, L.C.; Mattingly, J.D. Aircraft Engine Controls Design, System Analysis and Health Monitoring; American Institute of Aeronautics and Astronautics: Reston, VA, USA, 2009. [Google Scholar]
  10. Komjáty, M.; Főző, L.; Andoga, R. Experimental identification of a small turbojet engine with variable exhaust nozzle. In Proceedings of the 2015 16th IEEE International Symposium on Computational Intelligence and Informatics (CINTI), Budapest, Hungary, 19–21 November 2015; pp. 65–69. [Google Scholar] [CrossRef]
  11. Bakalis, D.P.; Stamatis, A.G. Data analysis and performance model calibration of a small turbojet engine. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2012, 226, 1523–1533. [Google Scholar] [CrossRef]
  12. Koruyucu, E.; Ekici, S.; Karakoc, T.H. Performing thermodynamic analysis by simulating the general characteristics of the two-spool turbojet engine suitable for drone and uav propulsion. J. Therm. Anal. Calorim. 2021, 145, 1303–1315. [Google Scholar] [CrossRef]
  13. Coban, K.; Ekici, S.; Colpan, C.O.; Karakoç, T.H. Performance of a microjet using component map scaling. Aircr. Eng. Aerosp. Technol. 2022, 94, 1201–1219. [Google Scholar] [CrossRef]
  14. Erario, M.L.; Giorgi, M.G.D.; Przysowa, R. Model-based dynamic performance simulation of a microturbine using flight test data. Aerospace 2022, 9, 60. [Google Scholar] [CrossRef]
  15. Villarreal-Valderrama, F.; Liceaga-Castro, E.; Zambrano-Robledo, P.; Amezquita-Brooks, L. Experimental evaluation of different micro-turbojet egt modeling approaches. J. Aerosp. Eng. Forthcom. 2021, 34, 04020087. [Google Scholar] [CrossRef]
  16. Guan, Y.; Zhao, G.; Xiao, X. Design and experiments of plasma jet igniter for aeroengine. Propuls. Power Res. 2013, 2, 188–193. [Google Scholar] [CrossRef]
  17. Najjar, Y.S.; Balawneh, I.A. Optimization of gas turbines for sustainable turbojet propulsion. Propuls. Power Res. 2015, 4, 114–121. [Google Scholar] [CrossRef]
  18. Şöhret, Y. A comprehensive approach to understanding irreversibility in a turbojet. Propuls. Power Res. 2018, 7, 129–137. [Google Scholar] [CrossRef]
  19. Cruz-Manzo, S.; Panov, V.; Zhang, Y. Gas path fault and degradation modelling in twin-shaft gas turbines. Machines 2018, 6, 43. [Google Scholar] [CrossRef]
  20. Zhao, H.; Liao, Z.; Liu, J.; Li, M.; Liu, W.; Wang, L.; Song, Z. A highly robust thrust estimation method with dissimilar redundancy framework for gas turbine engine. Energy 2022, 245, 123255. [Google Scholar] [CrossRef]
  21. Kayaalp, K.; Metlek, S.; Ekici, S.; Şöhret, Y. Developing a model for prediction of the combustion performance and emissions of a turboprop engine using the long short-term memory method. Fuel 2021, 302, 121202. [Google Scholar] [CrossRef]
  22. Kilic, U.; Villareal-Valderrama, F.; Ayar, M.; Ekici, S.; Amezquita-Brooks, L.; Karakoc, T.H. Deep Learning-Based Forecasting Modeling of Micro Gas Turbine Performance Projection: An Experimental Approach, Engineering Applications of Artificial Intelligence; Elsevier: Amsterdam, The Netherlands, 2024; Volume 130. [Google Scholar] [CrossRef]
  23. Menga, N.; Mothakani, A.; Giorgi, M.G.D.; Przysowa, R.; Ficarella, A. Extreme learning machine-based diagnostics for component degradation in a microturbine. Energies 2022, 15, 7304. [Google Scholar] [CrossRef]
  24. Ying, Y.; Cao, Y.; Li, S.; Li, J.; Guo, J. Study on gas turbine engine fault diagnostic approach with a hybrid of gray relation theory and gas-path analysis. Adv. Mech. Eng. 2016, 8, 1687814015627769. [Google Scholar] [CrossRef]
  25. González-Castillo, I.; Loboda, I. Analysis of nonlinear gas turbine models using influence coefficients. Ing. Investig. Tecnol. 2021, 22, 1–17. [Google Scholar] [CrossRef]
  26. Ntantis, E.; Li, Y. The impact of measurement noise on gas turbine gpa diagnostics. In Proceedings of the 6th International Conference on Condition Monitoring & Machinery Failure Prevention Technologies, Dublin, Ireland, 23–26 June 2009; pp. 23–25. [Google Scholar]
  27. Hajela, P.; Vittal, S. Optimal design in the presence of modeling uncertainties. J. Aerosp. Eng. 2006, 19, 204–216. [Google Scholar] [CrossRef]
  28. Schoukens, M.; Noël, J. Wiener-hammerstein benchmark with process noise. In Proceedings of the Workshop on Nonlinear System Identification Benchmarks, Lugano, Switzerland, 24–26 April 2016; pp. 15–19. [Google Scholar]
  29. Rubio-Sierra, C.; Delgado, A.; Fernández-Martínez, J.L. Maneuver optimization for simultaneous airspeed calibration and wind estimation. J. Aerosp. Eng. 2022, 35, 04022004. [Google Scholar] [CrossRef]
  30. Hale, L.E.; Patil, M.; Roy, C.J. Aerodynamic parameter identification and uncertainty quantification for small unmanned aircraft. J. Guid. Control Dyn. 2016, 40, 680–691. [Google Scholar] [CrossRef]
  31. Loboda, I.; Zárate, L.A.M.; Yepifanov, S.; Herrera, C.M.; Ruiz, J.L.P. Estimation of gas turbine unmeasured variables for an online monitoring system. Int. J. Turbo Jet-Eng. 2020, 37, 413–428. [Google Scholar] [CrossRef]
  32. Mehrpanahi, A.; Payganeh, G.; Arbabtafti, M. Dynamic modeling of an industrial gas turbine in loading and unloading conditions using a gray box method. Energy 2017, 120, 1012–1024. [Google Scholar] [CrossRef]
  33. Huang, J.; Chen, Y.; Guo, K. Novel approach to multibody system modeling: Cascading and clustering. J. Aerosp. Eng. 2014, 27, 279–290. [Google Scholar] [CrossRef]
  34. Li, X.; Wang, Y.; Zhu, Y.; Yang, G.; Liu, H. Temperature prediction of combustion level of ultra-supercritical unit through data mining and modelling. Energy 2021, 231, 120875. [Google Scholar] [CrossRef]
  35. Paulescu, M.; Brabec, M.; Boata, R.; Badescu, V. Structured, physically inspired (gray box) models versus black box modeling for forecasting the output power of photovoltaic plants. Energy 2017, 121, 792–802. [Google Scholar] [CrossRef]
  36. Villarreal-Valderrama, F.; Juárez-Pérez, P.; García-Pérez, U.; Amezquita-Brooks, L. Applicability of correlational data-mining to small-scale turbojet performance map generation. Int. J. Turbo Jet-Eng. 2024, 40, s67–s75. [Google Scholar] [CrossRef]
  37. Sheng, H.; Zhang, T. Aircraft engine thrust estimator design based on gsa-lssvm. Int. J. Turbo Jet-Eng. 2017, 34, 279–285. [Google Scholar] [CrossRef]
  38. Villarreal-Valderrama, F.; Delgado, C.S.; Zambrano-Robledo, P.; Amezquita-Brooks, L. Turbojet direct-thrust control scheme for full-envelope fuel consumption minimization. Aircr. Eng. Aerosp. Technol. 2021, 93, 437–447. [Google Scholar] [CrossRef]
  39. Andoga, R.; Főző, L.; Kovács, R.; Beneda, K.; Moravec, T.; Schreiner, M. Robust control of small turbojet engines. Machines 2019, 7, 3. [Google Scholar] [CrossRef]
  40. Tang, W.; Wang, L.; Gu, J.; Gu, Y. Single neural adaptive pid control for small uav micro-turbojet engine. Sensors 2020, 20, 345. [Google Scholar] [CrossRef] [PubMed]
  41. NASA. Rules for the Design, Development, Verification, and Operation of Flight Systems; GSFC-STD-1000; NASA: Washington, DC, USA, 2023. [Google Scholar]
  42. Airworthiness Certification Criteria, MIL-HDBK-516B; Department of Defense: Arlington, VA, USA, 2005.
  43. Yazar, I.; Kiyak, E.; Caliskan, F.; Karakoc, T.H. Simulation-based dynamic model and speed controller design of a small-scale turbojet engine. Aircr. Eng. Aerosp. Technol. 2018, 90, 351–358. [Google Scholar] [CrossRef]
  44. Montazeri-Gh, M.; Rasti, A.; Jafari, A.; Ehteshami, M. Design and implementation of mpc for turbofan engine control system. Aerosp. Sci. Technol. 2019, 92, 99–113. [Google Scholar] [CrossRef]
  45. Ahmed, S.; Mohamed, K.; Ashry, M.M. Controller design for micro turbojet engine. In Proceedings of the 2020 12th International Conference on Electrical Engineering (ICEENG), Cairo, Egypt, 7–9 July 2020; pp. 436–440. [Google Scholar]
  46. Xin, M.; Pan, H. Robust control of pvtol aircraft with a nonlinear optimal control solution. J. Aerosp. Eng. 2010, 23, 265–275. [Google Scholar] [CrossRef]
  47. Ramirez, J.C.H.; Nahon, M. Trajectory tracking control of highly maneuverable fixed-wing unmanned aerial vehicles. In Proceedings of the AIAA Scitech 2020 Forum, Orlando, FL, USA, 6–10 January 2020; p. 2074. [Google Scholar]
  48. Yu, B.; Shu, W.; Cao, C. A novel modeling method for aircraft engine using nonlinear autoregressive exogenous (narx) models based on wavelet neural networks. Int. J. Turbo Jet-Eng. 2018, 35, 161–169. [Google Scholar] [CrossRef]
  49. Giorgi, M.G.D.; Strafella, L.; Ficarella, A. Neural nonlinear autoregressive model with exogenous input (narx) for turboshaft aeroengine fuel control unit model. Aerospace 2021, 8, 206. [Google Scholar] [CrossRef]
  50. Zhang, B.; Pi, Y. Enhanced robust fractional order proportional-plus-integral controller based on neural network for velocity control of permanent magnet synchronous motor. ISA Trans. 2013, 52, 510–516. [Google Scholar] [CrossRef] [PubMed]
  51. Perng, J.-W.; Wu, B.-F.; Lee, T.-T. Limit cycle prediction of a neural vehicle control system with gain-phase margin tester. In Proceedings of the The 2006 IEEE International Joint Conference on Neural Network Proceedings, Vancouver, BC, Canada, 16–21 July 2006; pp. 4972–4977. [Google Scholar] [CrossRef]
  52. Amezquita-Brooks, L.; Liceaga-Castro, E.; Gonzalez-Sanchez, M.; Garcia-Salazar, O.; Martinez-Vazquez, D. Towards a standard design model for quad-rotors: A review of current models, their accuracy and a novel simplified model. Prog. Aerosp. Sci. 2017, 95, 1–23. [Google Scholar] [CrossRef]
  53. Qian, G.-W.; Song, Y.-P.; Ishihara, T. A control-oriented large eddy simulation of wind turbine wake considering effects of coriolis force and time-varying wind conditions. Energy 2022, 239, 121876. [Google Scholar] [CrossRef]
  54. Bahiuddin, I.; Mazlan, S.A.; Imaduddin, F. A new control-oriented transient model of variable geometry turbocharger. Energy 2017, 125, 297–312. [Google Scholar] [CrossRef]
  55. Suri, G.; Onori, S. A control-oriented cycle-life model for hybrid electric vehicle lithium-ion batteries. Energy 2016, 96, 644–653. [Google Scholar] [CrossRef]
  56. Gao, Y.; Liu, C.; Chen, S.; Zhang, X.; Fan, G.; Zhu, C. Development and parameterization of a control-oriented electrochemical model of lithium-ion batteries for battery-management-systems applications. Appl. Energy 2022, 309, 118521. [Google Scholar] [CrossRef]
  57. Li, D.; Yang, L.; Li, C. Control-oriented thermal-electrochemical modeling and validation of large size prismatic lithium battery for commercial applications. Energy 2021, 214, 119057. [Google Scholar] [CrossRef]
  58. Zhang, Q.; Tong, Z.; Tong, S.; Cheng, Z.; Lu, L. Control-oriented modeling of purging process and cold start of proton-exchange membrane fuel cell. J. Energy Eng. 2021, 147, 04021031. [Google Scholar] [CrossRef]
  59. Long, S.; Marjanovic, O.; Parisio, A. Generalised control-oriented modelling framework for multi-energy systems. Appl. Energy 2019, 235, 320–331. [Google Scholar] [CrossRef]
  60. Firoozabadi, M.D.; Shahbakhti, M.; Koch, C.; Jazayeri, S. Thermodynamic control-oriented modeling of cycle-to-cycle exhaust gas temperature in an hcci engine. Appl. Energy 2013, 110, 236–243. [Google Scholar] [CrossRef]
  61. Molina, S.; Guardiola, C.; Martín, J.; García-Sarmiento, D. Development of a control-oriented model to optimise fuel consumption and nox emissions in a di diesel engine. Appl. Energy 2014, 119, 405–416. [Google Scholar] [CrossRef]
  62. d’Ambrosio, S.; Finesso, R.; Fu, L.; Mittica, A.; Spessa, E. A control-oriented real-time semi-empirical model for the prediction of nox emissions in diesel engines. Appl. Energy 2014, 130, 265–279. [Google Scholar] [CrossRef]
  63. Poksawat, P.; Wang, L.; Mohamed, A. Gain scheduled attitude control of fixed-wing uav with automatic controller tuning. IEEE Trans. Control Syst. Technol. 2017, 26, 1192–1203. [Google Scholar] [CrossRef]
  64. Grimble, M.J. Robust Industrial Control Systems: Optimal Design Approach for Polynomial Systems; John Wiley & Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  65. Kulikov, G.G.; Thompson, H.A. Dynamic Modelling of Gas Turbines: Identification, Simulation, Condition Monitoring and Optimal Control; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  66. Iqbal, I.M.; Aziz, N. Comparison of various wiener model identification approach in modelling nonlinear process. In Proceedings of the 2011 3rd Conference on Data Mining and Optimization (DMO), Putrajaya, Malaysia, 28–29 June 2011; pp. 134–140. [Google Scholar] [CrossRef]
  67. Shamma, J.S.; Athans, M. Analysis of gain scheduled control for nonlinear plants. IEEE Trans. Autom. Control 1990, 35, 898–907. [Google Scholar] [CrossRef]
  68. Smith, R.S. Model validation for robust control: An experimental process control application. Automatica 1995, 31, 1637–1647. [Google Scholar] [CrossRef]
  69. Cervantes, A.L.; Agamennoni, O.E.; Figueroa, J.L. A nonlinear model predictive control system based on wiener piecewise linear models. J. Process. Control 2003, 13, 655–666. [Google Scholar] [CrossRef]
  70. Norquay, S.J.; Palazoglu, A.; Romagnoli, J. Model predictive control based on wiener models. Chem. Eng. Sci. 1998, 53, 75–84. [Google Scholar] [CrossRef]
  71. Aström, K.J.; Wittenmark, B. Adaptive Control; Courier Corporation: Chelmsford, MA, USA, 2013. [Google Scholar]
  72. Wigren, T.; Brus, L. Reduction of amplitude dependent gain variations in control of non-linear wiener type systems. IFAC Proc. Vol. 2007, 40, 336–341. [Google Scholar] [CrossRef]
  73. Leith, D.; Leithead, W. Equivalence of gain-scheduling and input-output linearisation for a class of commonly occurring plants. In Proceedings of the 36th IEEE Conference on Decision and Control, San Diego, CA, USA, 12 December 1997; Volume 1, pp. 430–431. [Google Scholar] [CrossRef]
  74. Biagiola, S.I.; Agamennoni, O.E.; Figueroa, J.L. Robust control of wiener systems: Application to a ph neutralization process. Braz. J. Chem. Eng. 2016, 33, 145–153. [Google Scholar] [CrossRef]
  75. Kalafatis, A.D.; Wang, L.; Cluett, W.R. Linearizing feedforward–feedback control of ph processes based on the wiener model. J. Process. Control 2005, 15, 103–112. [Google Scholar] [CrossRef]
  76. Hu, Y.; Wu, J.; Zeng, C. Robust adaptive identification of linear time-varying systems under relaxed excitation conditions. IEEE Access 2020, 8, 8268–8274. [Google Scholar] [CrossRef]
  77. Li, Z.; Nikolaidis, T.; Nalianda, D. Recursive least squares for online dynamic identification on gas turbine engines. J. Guid. Control Dyn. 2016, 39, 2594–2601. [Google Scholar] [CrossRef]
  78. Yu, S.-S.; Chu, S.-W.; Wang, C.-M.; Chan, Y.-K.; Chang, T.-C. Two improved k-means algorithms. Appl. Soft Comput. 2018, 68, 747–755. [Google Scholar] [CrossRef]
  79. Cohen-Addad, V.; Kanade, V.; Mallmann-Trenn, F.; Mathieu, C. Hierarchical clustering: Objective functions and algorithms. J. ACM (JACM) 2019, 66, 1–42. [Google Scholar] [CrossRef]
  80. Atasoy, V.E.; Suzer, A.E.; Ekici, S. A comparative analysis of exhaust gas temperature based on machine learning models for aviation applications. J. Energy Resour. Technol. 2022, 144, 082101. [Google Scholar] [CrossRef]
  81. Chiras, N.; Evans, C.; Rees, D. Global nonlinear modeling of gas turbine dynamics using narmax structures. J. Eng. Gas Turbines Power 2002, 124, 817–826. [Google Scholar] [CrossRef]
  82. Kovacs, R.; Fozo, L.; Andoga, R.; Moravec, T. A non-linear model of a complex system using variable parameters. In Proceedings of the 2018 IEEE 16th World Symposium on Applied Machine Intelligence and Informatics (SAMI), Kosice and Herlany, Slovakia, 7–10 February 2018; pp. 000113–000118. [Google Scholar] [CrossRef]
  83. Torres, M.P.; Sosa, G.; Amezquita-Brooks, L.; Liceaga-Castro, E.; Zambrano-Robledo, P.d.C. Identification of the fuel-thrust dynamics of a gas turbo engine. In Proceedings of the 52nd IEEE Conference on Decision and Control, Firenze, Italy, 10–13 December 2013; pp. 4535–4540. [Google Scholar] [CrossRef]
  84. Kong, C.; Lim, S. Inverse generation of gas turbine component performance maps from experimental test data. Int. J. Turbo Jet-Engines 2010, 27, 135–148. [Google Scholar] [CrossRef]
  85. Yu, D.; Zhao, H.; Xu, Z.; Sui, Y.; Liu, J. An approximate non-linear model for aeroengine control. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2011, 225, 1366–1381. [Google Scholar] [CrossRef]
  86. Ljung, L.; Singh, R.; Zhang, Q.; Lindskog, P.; Iouditski, A. Developments in the mathworks system identification toolbox. Ifac Proc. Vol. 2009, 42, 522–527. [Google Scholar] [CrossRef]
  87. Atam, E.; Schulte, D.O.; Arteconi, A.; Sass, I.; Helsen, L. Control-oriented modeling of geothermal borefield thermal dynamics through hammerstein-wiener models. Renew. Energy 2018, 120, 468–477. [Google Scholar] [CrossRef]
  88. Copaci, D.S.; Moreno, L.E.; Blanco, M.D. Two-stage shape memory alloy identification based on hammerstein-wiener model. Front. Robot. AI 2019, 6, 83. [Google Scholar] [CrossRef] [PubMed]
  89. Mulyana, T. Nnarx model structure for the purposes of controller design and optimization of heat exchanger process control training system operation. In AIP Conference Proceedings; AIP Publishing LLC.: Long Island, NY, USA, 2017; Volume 1831, p. 020040. [Google Scholar]
  90. Kishor, N.; Saini, R.; Singh, S. Small hydro power plant identification using nnarx structure. Neural Comput. Appl. 2005, 14, 212–222. [Google Scholar] [CrossRef]
  91. Koleini, I.; Roudbari, A.; Marefat, V. Egt prediction of a micro gas turbine using statistical and artificial intelligence approach. IEEE Aerosp. Electron. Syst. Mag. 2018, 33, 4–13. [Google Scholar] [CrossRef]
  92. Giorgi, M.G.D.; Quarta, M. Hybrid multigene genetic programming-artificial neural networks approach for dynamic performance prediction of an aeroengine. Aerosp. Sci. Technol. 2020, 103, 105902. [Google Scholar] [CrossRef]
  93. Watanabe, A.; Ölçmen, S.M.; Leland, R.P.; Whitaker, K.W.; Trevino, L.C.; Nott, C. Soft computing applications on a sr-30 turbojet engine. Fuzzy Sets Syst. 2006, 157, 3007–3024. [Google Scholar] [CrossRef]
  94. Francisco, V.-V. Analysis and Modeling of Micro Turbojets: A Comprehensive Model Based on Multiphysics Principles. Master’s Thesis, Universidad Autonoma de Nuevo Leon, San Nicolás de los Garza, Mexico, 2019. [Google Scholar]
  95. Tavakolpour-Saleh, A.; Nasib, S.; Sepasyan, A.; Hashemi, S. Parametric and nonparametric system identification of an experimental turbojet engine. Aerosp. Sci. Technol. 2015, 43, 21–29. [Google Scholar] [CrossRef]
  96. Tagashira, T.; Sugiyama, N.; Koh, M. Dynamic characteristic tests of single spool turbojet engine using altitude test facility. In Proceedings of the 43rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, Cincinnati, OH, USA, 8–11 July 2007; p. 5012. [Google Scholar] [CrossRef]
  97. Billings, S.; Tsang, K. Spectral analysis of block structured nonlinear systems. Mech. Syst. Signal Process. 1990, 4, 117–130. [Google Scholar] [CrossRef]
  98. Khaoula, D.; Beneda, K. Linear Dynamic Mathematical Model and Identification of Micro Turbojet Engine for Turbofan Power Ratio Control. Aviation 2019, 23, 54–64. [Google Scholar] [CrossRef]
  99. Salehi, A.; Montazeri-Gh, M. Black box modeling of a turboshaft gas turbine engine fuel control unit based on neural NARX. Proc. Inst. Mech. Eng. Part M J. Eng. Marit. Environ. 2019, 233, 949–956. [Google Scholar] [CrossRef]
  100. Cousineau, D.; Brown, S.; Heathcote, A. Fitting distributions using maximum likelihood: Methods and packages. Behav. Res. Methods Instrum. Comput. 2004, 36, 742–756. [Google Scholar] [CrossRef]
  101. Wei, X.; Liu, X.; Fan, Y.; Tan, L.; Liu, Q. A Unified Test for the AR Error Structure of an Autoregressive Model. Axioms 2022, 11, 690. [Google Scholar] [CrossRef]
  102. Johnson, A.R.; Irwin, M.; Freund, E.J. Probability and Statistics for Engineers, Engineering Applications of Artificial Intelligence; Pearson Education: London, UK, 2000. [Google Scholar]
  103. Evans, C.; Rees, D.; Borrell, A. Identification of aircraft gas turbine dynamics using frequency-domain techniques. Control Eng. Pract. 2000, 8, 457–467. [Google Scholar] [CrossRef]
  104. Ruano, A.E.; Fleming, P.J.; Teixeira, C.; RodrıGuez-Vázquez, K.; Fonseca, C.M. Nonlinear identification of aircraft gas-turbine dynamics. Neurocomputing 2003, 55, 551–579. [Google Scholar] [CrossRef]
  105. Wessley, G.J.J.; Chauhan, S. Investigation on scaling of gas turbine engines for drone propulsion. Int. J. Eng. Technol. Manag. Appl. Sci. 2017, 5, 48–53. [Google Scholar]
  106. Wessley, G.J.J.; Chauhan, S. Parametric analysis of a down-scaled turbo jet engine suitable for drone and uav propulsion. In AIP Conference Proceedings; AIP Publishing: Long Island, NY, USA, 2018; p. 020024. [Google Scholar] [CrossRef]
Figure 1. Gray-box models combine the advantages of theoretical models, which yield certainty in the model structure, and experimental models, which yield certainty in the computation of the model parameters. Data mining can also be used to determine a consistent model structure, resulting in a gray-box model as well.
Figure 1. Gray-box models combine the advantages of theoretical models, which yield certainty in the model structure, and experimental models, which yield certainty in the computation of the model parameters. Data mining can also be used to determine a consistent model structure, resulting in a gray-box model as well.
Aerospace 11 00630 g001
Figure 2. Block diagram of Wiener- and Hammerstein-structured models, respectively.
Figure 2. Block diagram of Wiener- and Hammerstein-structured models, respectively.
Aerospace 11 00630 g002
Figure 3. Proposed gray-box model identification flowchart.
Figure 3. Proposed gray-box model identification flowchart.
Aerospace 11 00630 g003
Figure 4. Microturbojet cutaway view. The SR-30 is a commercial microturbojet used in UAV propulsion.
Figure 4. Microturbojet cutaway view. The SR-30 is a commercial microturbojet used in UAV propulsion.
Aerospace 11 00630 g004
Figure 5. Fuel flow input for the identification, validation, and frequency-swept tests: (a) Identification data set. The first data segment of this test corresponds to the engine warm-up to avoid any nonlinearities caused by thermal effects and the second segment includes the identification data used for each high-order linear model. (b) Validation data set. Excitation signal used for model validation considering different operating conditions. (c) Frequency–response test. A frequency-swept test is performed to assess the predictive power of the identified models with respect to the process dynamics.
Figure 5. Fuel flow input for the identification, validation, and frequency-swept tests: (a) Identification data set. The first data segment of this test corresponds to the engine warm-up to avoid any nonlinearities caused by thermal effects and the second segment includes the identification data used for each high-order linear model. (b) Validation data set. Excitation signal used for model validation considering different operating conditions. (c) Frequency–response test. A frequency-swept test is performed to assess the predictive power of the identified models with respect to the process dynamics.
Aerospace 11 00630 g005
Figure 6. Shaft speed cluster and frequency domain analysis. Each identified model corresponds to an independent step–response test (i.e., A, B, C, D, E, F and H) while the reduced model is G s . (a) Cluster analysis of the shaft speed models. The reduced model is a 3rd order transfer function. (b) Frequency-domain analysis of the identified shaft speed linear models. The reduced model G s captures the behavior of the high-order models up to a frequency of 4 rads/s, rejecting inconsistent higher-frequency components.
Figure 6. Shaft speed cluster and frequency domain analysis. Each identified model corresponds to an independent step–response test (i.e., A, B, C, D, E, F and H) while the reduced model is G s . (a) Cluster analysis of the shaft speed models. The reduced model is a 3rd order transfer function. (b) Frequency-domain analysis of the identified shaft speed linear models. The reduced model G s captures the behavior of the high-order models up to a frequency of 4 rads/s, rejecting inconsistent higher-frequency components.
Aerospace 11 00630 g006
Figure 7. Shaft speed response identification. (a) Comparison of the linear model and the nonlinear model accuracy. (b) Experimental identification of the nonlinear behavior of the shaft speed static gain.
Figure 7. Shaft speed response identification. (a) Comparison of the linear model and the nonlinear model accuracy. (b) Experimental identification of the nonlinear behavior of the shaft speed static gain.
Aerospace 11 00630 g007
Figure 8. Thrust cluster and frequency domain analysis. Each identified model corresponds to an independent step-response test (i.e., A, B, C, D, E, F, and H), while the reduced model is G t . (a) Cluster analysis of the identified and reduced thrust models. The reduced model is a 1st order transfer function. (b) Frequency-domain analysis of the identified thrust linear models. The reduced model G t captures the behavior of the high-order models up to a frequency of 2 rads/s. The effects of a high noise-to-signal ratio are presented in the Bode plot scattering of the identified models. This phenomenon is caused by the influence of the measurement noise and uncertainty in the identified models.
Figure 8. Thrust cluster and frequency domain analysis. Each identified model corresponds to an independent step-response test (i.e., A, B, C, D, E, F, and H), while the reduced model is G t . (a) Cluster analysis of the identified and reduced thrust models. The reduced model is a 1st order transfer function. (b) Frequency-domain analysis of the identified thrust linear models. The reduced model G t captures the behavior of the high-order models up to a frequency of 2 rads/s. The effects of a high noise-to-signal ratio are presented in the Bode plot scattering of the identified models. This phenomenon is caused by the influence of the measurement noise and uncertainty in the identified models.
Aerospace 11 00630 g008
Figure 9. Thrust response identification: (a) Comparison of the thrust linear and nonlinear model accuracy. (b) Experimental identification of the nonlinear behavior of the thrust static gain.
Figure 9. Thrust response identification: (a) Comparison of the thrust linear and nonlinear model accuracy. (b) Experimental identification of the nonlinear behavior of the thrust static gain.
Aerospace 11 00630 g009
Figure 10. Structure of Hammerstein–Wiener models. These models include a nonlinear static function dependent on the input and a nonlinear static function dependent on the linear model G ( z ) output.
Figure 10. Structure of Hammerstein–Wiener models. These models include a nonlinear static function dependent on the input and a nonlinear static function dependent on the linear model G ( z ) output.
Aerospace 11 00630 g010
Figure 11. Neural networks used for training and practical application: (a) Neural network diagram considering the dynamic behavior similarly to a nonlinear difference equation model. (b) Neural Network AutoRegressive with Exogenous Input. This is the Output Error (OE) configuration with delayed predicted outputs instead of measurements.
Figure 11. Neural networks used for training and practical application: (a) Neural network diagram considering the dynamic behavior similarly to a nonlinear difference equation model. (b) Neural Network AutoRegressive with Exogenous Input. This is the Output Error (OE) configuration with delayed predicted outputs instead of measurements.
Aerospace 11 00630 g011
Figure 12. Models performance with the identification data. Overall, all models perform similarly for this test with the Hammerstein–Wiener model (orange) providing a slightly faster response than the proposed gray-box model (blue). Note that the neural network (green), for this test only, is not structured in the OE form: (a) Performance of the resulting models with the shaft speed identification data. (b) Performance of the resulting models with the thrust identification data.
Figure 12. Models performance with the identification data. Overall, all models perform similarly for this test with the Hammerstein–Wiener model (orange) providing a slightly faster response than the proposed gray-box model (blue). Note that the neural network (green), for this test only, is not structured in the OE form: (a) Performance of the resulting models with the shaft speed identification data. (b) Performance of the resulting models with the thrust identification data.
Aerospace 11 00630 g012
Figure 13. Shaft speed and thrust models validation. The model validation experiments included ramp, sinusoidal, and step-like changes in the fuel flow input to the microturbojet: (a) Shaft speed linear and nonlinear models accuracy with the validation data set. (b) Thrust linear and nonlinear models accuracy with the validation data set.
Figure 13. Shaft speed and thrust models validation. The model validation experiments included ramp, sinusoidal, and step-like changes in the fuel flow input to the microturbojet: (a) Shaft speed linear and nonlinear models accuracy with the validation data set. (b) Thrust linear and nonlinear models accuracy with the validation data set.
Aerospace 11 00630 g013
Figure 14. Experimental frequency-swept response and models approximatio: (a) Shaft speed models evaluation. (b) Thrust models evaluation.
Figure 14. Experimental frequency-swept response and models approximatio: (a) Shaft speed models evaluation. (b) Thrust models evaluation.
Aerospace 11 00630 g014
Figure 15. Statistic analysis (Box Plot) of the modeling errors with the time-domain validation (GB for the gray-box model, H-W for the Hammerstein–Wiener model, and NN for the neural network).
Figure 15. Statistic analysis (Box Plot) of the modeling errors with the time-domain validation (GB for the gray-box model, H-W for the Hammerstein–Wiener model, and NN for the neural network).
Aerospace 11 00630 g015
Figure 16. Statistic analysis (Box Plot) of the modeling errors with the frequency-domain validation.
Figure 16. Statistic analysis (Box Plot) of the modeling errors with the frequency-domain validation.
Aerospace 11 00630 g016
Figure 17. Sensor measurements and best fit normal distributions at static regime.
Figure 17. Sensor measurements and best fit normal distributions at static regime.
Aerospace 11 00630 g017
Figure 18. Aggregated error histogram of the MC simulation and best-fit distribution of each model.
Figure 18. Aggregated error histogram of the MC simulation and best-fit distribution of each model.
Aerospace 11 00630 g018
Table 1. Nonlinear models performance metrics.
Table 1. Nonlinear models performance metrics.
VariableModelRMSPE %SERMSPE %SE
Time DomainFrequency Domain
Shaft speedG-B 1.3429 26.613 2.0327 11.641
W-H 5.5167 121.06 8.5607 47.835
NN 3.1928 52.280 3.8578 27.555
ThrustG-B 4.3739 0.0580 6.4324 0.0309
W-H 8.0221 0.1431 14.838 0.0552
NN 7.4271 0.1466 17.136 0.0489
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Villarreal-Valderrama, F.; Liceaga-Castro, E.; Hernandez-Alcantara, D.; Santana-Delgado, C.; Ekici, S.; Amezquita-Brooks, L. Control-Oriented System Identification of Turbojet Dynamics. Aerospace 2024, 11, 630. https://doi.org/10.3390/aerospace11080630

AMA Style

Villarreal-Valderrama F, Liceaga-Castro E, Hernandez-Alcantara D, Santana-Delgado C, Ekici S, Amezquita-Brooks L. Control-Oriented System Identification of Turbojet Dynamics. Aerospace. 2024; 11(8):630. https://doi.org/10.3390/aerospace11080630

Chicago/Turabian Style

Villarreal-Valderrama, Francisco, Eduardo Liceaga-Castro, Diana Hernandez-Alcantara, Carlos Santana-Delgado, Selcuk Ekici, and Luis Amezquita-Brooks. 2024. "Control-Oriented System Identification of Turbojet Dynamics" Aerospace 11, no. 8: 630. https://doi.org/10.3390/aerospace11080630

APA Style

Villarreal-Valderrama, F., Liceaga-Castro, E., Hernandez-Alcantara, D., Santana-Delgado, C., Ekici, S., & Amezquita-Brooks, L. (2024). Control-Oriented System Identification of Turbojet Dynamics. Aerospace, 11(8), 630. https://doi.org/10.3390/aerospace11080630

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop