Next Article in Journal
A New Algorithm on Automatic Trimming for Helicopter Rotor Aerodynamic Loads
Previous Article in Journal
Surrogate Aerodynamic Wing Modeling Based on a Multilayer Perceptron
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance-Seeking Control of Aero-Propulsion System Based on Intelligent Optimization and Active Disturbance Rejection Fusion Controller

School of Control Theory and Control Engineering, Dalian University of Technology, Dalian 116024, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(2), 151; https://doi.org/10.3390/aerospace10020151
Submission received: 17 October 2022 / Revised: 1 February 2023 / Accepted: 2 February 2023 / Published: 7 February 2023
(This article belongs to the Section Aeronautics)

Abstract

:
To investigate the performance potential of the aero-propulsion system and the problem of control mode conversion, this paper takes the inlet/engine integrated component-level model as the research object, and a performance-seeking control (PSC) scheme based on the neighborhood-based speciation differential evolution–grey wolf optimization (NSDE-GWO) algorithm is designed and combined with an active disturbance rejection control (ADRC) to establish a multivariable fusion closed-loop control system. The analysis reveals that the NSDE-GWO hybrid algorithm, which takes advantages of the two algorithms, significantly improves the computing efficiency and optimization accuracy, achieving better optimization solutions in three different modes. The intelligent fusion controller is able to achieve a smooth transition of performance modes to ensure that the engine is provided with a stable thrust during operation under supersonic conditions, and the potential for performance optimization is maintained at a reasonable level. Maximum installed thrust mode capable of achieving no thrust loss and a maximum fluctuation rate within 2000 N/s, with the largest variation in thrust during the conversion being less than 0.9% under the minimum turbine temperature mode and the minimum specific fuel consumption mode. This study presents a theoretical foundation and engineering applications for the design of supersonic propulsion system controllers.

1. Introduction

With the rapid development of supersonic aero technology, the propulsion system controller requires increased reliability and performance potential, and the control system object increasingly transfers from engine to comprehensive propulsion system [1,2]. The present control system mostly uses sensor-based control design principles to handle non-measurable parameters, including thrust, temperature, and safety margins indirectly through measurable parameters. Propulsion systems are characterized by complex structures and wide operating envelopes that cannot be fully exploited by conventional controllers, necessitating the use of fusion optimization approaches. performance-seeking control (PSC) [3,4] is an essential mode of intelligent engine control, which is implemented by using optimization algorithms to find the maximum performance potential (maximum thrust, minimum fuel consumption, or minimum turbine temperature) of the model while satisfying various operating limitations. Therefore, research on comprehensive performance-seeking control algorithms for supersonic propulsion systems is required, and the fusion and conversion of performance-seeking controls with sensor-based controllers necessitates concentrated effort.
Scholars have long researched and experimented with PSC, attempting to apply it to the propulsion system [5]. Nobbs. S [6] applied the full envelope PSC to PW1128 powered F-15 aircraft. NASA Dryden Flying Test Center’s Orme J.S [7] conducted the PSC supersonic flight test. Under supersonic conditions, thrust can be enhanced by 9%, fuel consumption can be lowered by 8%, and turbine temperature can be reduced by 48K. In terms of the numerical simulation study, the linear programming (LP) approach is employed as the performance optimization algorithm in the NASA report [8], and Orme, J.S [7] used LP in the airborne adaptive model. Later, numerous academics engaged in the research and expansion of LP, including SLP [9], FSQP [10], and SQCQP [11]. According to Chen HY [12], the FSQP could obtain the optimal operating conditions of the double bypass variable cycle engine in different conditions. With the rapid growth of artificial intelligence, intelligence algorithms have been gradually applied to the PSC, overcoming the initial value dependence of traditional methods. Several academics have used GA and EA to tackle optimization issues [13,14]. PSO provides an ideal optimization effect under maximum thrust optimization condition [15], but it depends on multiple-parameter tweaking, which cannot adapt to the multiple operating states. Other methods, such as the direct method [16], the Beetle Antennae Search algorithm [17], and IA [18] are also proposed. The above study focuses more on the PSC of turbofan engines and less on multivariable propulsion systems; the initial value sensitivity and the local convergence problems of the traditional algorithm will be more prominent due to the increased number of control variables and the harsh constraints of aero-propulsion modeling. The preceding study is primarily concerned with optimizing potential, a lack of quantitative analysis of stability, and the global convergence of intelligent PSC.
Although the PID controller is still employed in actual engines, it has drawbacks while dealing with strongly coupled multivariable objects. The research focus is on the construction of a multivariable controller with a wide flight envelope under supersonic conditions. NASA Glenn Research Center [19] argued that the multivariable control is mainly used for the requirements of engines in specific high-performance engine conditions. The earliest multivariable control study on aero engines was demonstrated by the LQR method [20] but suffered from accuracy dependence on mathematical models. Compared with LQR, LTR [21] has a better response speed and robustness of disturbance suppression but still exposes the defects in practical application conditions [22]. Richter H [23,24] carried out a study of sliding mode control, which is highly robust to parameter fluctuations and disturbances yet is prone to chattering. Active disturbance rejection control (ADRC) was invented by Han JQ [25]. Based on the approximate bandwidth concept, Gao ZQ [26] presented a simpler parameter setting approach for ADRC. Miklosovic R [27] integrated ADRC into turbofan engine control and validated anti-interference performance using NASA’s highly reliable model. Li SQ [28] compared ADRC to PID and determined that ADRC had superior control performance. Zhang HB [29] achieved turbofan engine decoupling control by integrating a decoupling matrix with ADRC. Xue WC [30] utilized ADRC to adjust the engine fuel ratio. It is clear that ADRC can be employed for decoupling and the robust control of turbofan engines, and its anti-interference performance is superior to PID. Previous research has focused on controller performance and robustness analysis rather than the combination of multivariable controllers and PSC. Nevertheless, the process of transitioning control systems to performance modes has received little attention.
To summarize, the performance-seeking control effects of NSDE-GWO and various swarm-intelligence algorithms are compared first, and then the benefits and application value of NSDE-GWO are analyzed. Finally, under typical supersonic conditions, a fusion-control system using the NSDE-GWO optimization algorithm and ADRC was created to accomplish a seamless transition from the PSC mode to quick reaction and low thrust loss.

2. Materials and Methods

The design of the fusion-control system basically contains the inlet/engine integrated model module, the control system module, and PSC modules. Considering the feasibility and cost of validation, the component-level model is commonly used in studies instead of actual engines. In this paper, a military turbofan engine with a variable geometry inlet and small bypass ratio is adopted, and an integrated component-level model of the inlet and engine is developed. The real engine model is commonly replaced by a component-level model in this paper, and a turbofan engine with a variable geometry inlet and small culvert ratio is adopted; the control system module contains an ADRC closed-loop controller and an inlet open-loop controller. The PSC module contains an online model and an intelligent optimization program, among which the online model is replaced by a component-level model, and the optimization program is adopted by the improved NSDE-GWO method. Figure 1 illustrates the operation principle of the fusion controller.

2.1. Supersonic Inlet/Engine Integrated Component-Level Model

Aero-propulsion systems are multivariable, nonlinear, and time-varying, and their models for control generally adopt component-level methods, specifically, models of multiple components based on their characteristic diagrams and aerodynamic-thermodynamic relationships. This paper describes the principle of building an inlet/engine integrated component-level model, which mainly includes two parts: the inlet system and the engine. The main workflow is shown in Figure 2. The main air path components include the external pressure inlet, fan, high-pressure compressor, combustion, high-pressure turbine, low-pressure turbine, mixing chamber, afterburner combustion, and nozzle.
Most of the conventional inlet components are calculated using idealized empirical formulas, which cannot satisfy the requirements of calculation accuracy at supersonic conditions. Considering the shock structure and drag calculation method, quasi one-dimensional flow theory is used to calculate the inlet performance [31]. Combined with basic flight conditions, an external pressure inlet with a combination of “two oblique and one normal” shock waves is selected. Through the geometric relationship of the two-dimensional plane, the structural dimension parameters are determined to seal the shock wave, where the first-stage ramp angle and second-stage ramp angle are adjustable variables and the bleed valve opening is represented by the relative bleed adjustment.
The object is a military small bypass ratio turbofan engine with an operating range of 0–2.5 Mach number interval, in which the fan-guide vane angle, the compressor-guide vane angle, and the nozzle area can be adjusted. The performance parameters of the fan and compressor components are realized using the characteristic interpolation, and the characteristic curves are shown in Figure 2. Simplified linearized characteristics were used for the high-pressure turbine and low-pressure turbine to improve the real-time performance, considering that the turbine mostly works at larger operating conditions and similarity exists in efficiency and flow rate at different speeds.
In addition, the output parameters of each gas path component model are governed by a common set of operating equations. Typical common equations include the pressure balance between the pressurized combustion chamber and nozzle, and the flow balance between the compressors of the fan; the power consumed on the same shaft should form a balance with the power generated. To achieve the inlet/engine matching, the residuals of the inlet and fan flows are taken as the balancing equations, and the ratio of the flow coefficient/total pressure recovery coefficient ( ϕ σ ) is chosen as the starting guess for the iteration. If the number of iterative variables is the same as the balancing equations, the problem is transformed into a nonlinear implicit set of equations with unknown independent variables, and the component-level model is considered to have a reliable solution when the residuals tend to zero. The hybrid iterative algorithm of the damping Newton method [32] is utilized to solve the issue due to the rise in iterative variables. [ x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 ] T denotes the independent assignment variables for the iterative process. The iterative variables (e) are selected, and co-working equations (Eq) are solved in association, expressed as
E q ( 1 ) : f 1 ( x 0 , x 1 , x 2 , , x 8 ) = W i n l e t W f a n , f = e 1 E q ( 2 ) : f 2 ( x 0 , x 1 , x 2 , , x 8 ) = W f a n , a W H P C , f = e 2 E q ( 3 ) : f 3 ( x 0 , x 1 , x 2 , , x 8 ) = I H · N h p r · d N h p r d t L H P T + L H P C L a d d = e 3 E q ( 4 ) : f 4 ( x 0 , x 1 , x 2 , , x 8 ) = W c o m b , a W H P T , f = e 4 E q ( 5 ) : f 5 ( x 0 , x 1 , x 2 , , x 8 ) = W H P T , a W L P T , f = e 5 E q ( 6 ) : f 6 ( x 0 , x 1 , x 2 , , x 8 ) = I L · N l p r · d N l p r d t L L P T + L f a n = e 6 E q ( 7 ) : f 7 ( x 0 , x 1 , x 2 , , x 8 ) = P 6 P 16 = e 7 E q ( 8 ) : f 8 ( x 0 , x 1 , x 2 , , x 8 ) = P 8 P n o z z l e = e 8
[ x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 ] T = [ ϕ σ , N h p r , Z f , N l p r , Z c , T 4 , X W G 41 c , X W G 45 c ] T
W f a n , a = W f a n , f W b y p a s s W c o m b , a = W c o m b , f + W f m
The formula parameters W i n l e t ,   W f a n , f ,   W f a n , a ,   W H P C , f ,   W c o m b , f ,   W c o m b , a ,   W H P T , f , W H P T , a , W L P T , f indicate, respectively, inlet outlet flow, fan inlet flow, fan outlet flow, compressor inlet flow, main combustion chamber inlet flow, main combustion chamber outlet flow, high-pressure turbine inlet flow, high-pressure turbine outlet flow, and low-pressure turbine inlet flow. I H , I L stand for the rotary inertia of high-pressure and low-pressure, and L H P T , L L P T , L H P C , L f a n stand for the power of HPT, HPC, HPC, and fan, respectively. L a d d stands for the additional power of engine starting. N h p r , N l p r indicate the high-pressure and low-pressure rotor speed. Z f ,   Z c ,   T 4 ,   X W G 41 c ,   X W G 45 c indicate the relative pressure ratio coefficients for fan, HPC, HPT, and LPT components, and initial guesses for characteristic interpolation. Due to the introduction of the external characteristics and taking into account the external resistance of the inlet, the installed thrust and ISFC are obtained by the following equation. D i n l e t indicates the external resistance calculated for the inlet tract.
F I T H R = F T H R D i n l e t I S F C = 3600 · w f m F I T H R

2.2. Performance-Seeking Control in Aero-Propulsion System

Under the supersonic conditions of the aero-propulsion system, the inlet/engine matching seriously affect the installed performance. PSC can optimize the performance by adjusting the inlet geometry parameters and compression parts parameters. At a certain flight altitude and Mach number, the model operating state is mainly determined by the control quantity, X, denoted as the state parameters. The optimized control variables(u) available for the advance integrated propulsion system include the following:
X = f ( u ) u = [ w f m , w f a , α f q , δ 1 , δ 2 , A 8 , α f , α c ] T
The formula parameters are as follows: relative bleed adjustment α fq , first-stage ramp angle δ 1 , second-stage ramp angle δ 2 , fan guide vane angle α f , compressor guide vane angle α c , main fuel flow rate wfm , afterburner fuel flow rate wfa , and exhaust nozzle outlet area A 8 . The ideal boundary conditions, as well as the stability of the inlet, must be determined while keeping the feasible area of the optimization space and operating limit in mind. If the Inlet flow coefficient φ in is too low, the inlet will enter an unstable functioning state known as “surge”. The flow coefficient of the surge working state is φ is . Increased inlet distortion may easily cause engine surge and negatively impact engine performance. Formula parameters: low-pressure turbine outlet temperature T 6 , and fan surge margin SMF, and compressor surge margin SMC. To summarize, the following requirements are guaranteed:
M = u min u u max N h p r , min N h p r N h p r , max N l p r , min N l p r N l p r , max T 6 T 6 , max φ i s φ i n φ max S M F min S M F S M C min S M C
PSC optimizes the engine operating point in different flight objectives based on the optimization control mode request, as well as the aircraft attitude, flight circumstances, and other information. It is separated into optimization control modes based on various optimization targets, such as maximum thrust, minimum fuel consumption, and minimum turbine temperature. Formula parameters: min ISFC penalty function coefficient K1, min T45 penalty function coefficient K2, optimized installed thrust F I T H R , o p t , and initial installed thrust F I T H R , o r i . The penalty function is introduced in the objective function to achieve the optimization objective of constant installation thrust. The penalty function parameters need to adopt suitable values; id they are too large, they may reduce the optimization efficiency, and if they are too small they may weaken the constant thrust effect. Considering the different magnitudes of ISFC and T45, K1 and K2 are taken as 0.1 and 100.
(1) Maximum installed thrust control
max   F I T H R Objective   f ( u ) = F I T H R s . t . M , w f m = c o n s t , w f a = c o n s t u = [ α f q , δ 1 , δ 2 , A 8 , α f , α c ] T
(2) Minimum installed specific fuel-consumption control
min   I S F C Objective   f ( u ) = I S F C   +   K 1 · F I T H R , o p t F I T H R , o r i s . t . M , F I T H R = c o n s t u = [ w f m , w f a , α f q , δ 1 , δ 2 , A 8 , α f , α c ] T
(3) Minimum low-pressure turbine inlet temperature control
min   T 45 Objective   f ( u ) = T 45   +   K 2 · F I T H R , o p t F I T H R , o r i s . t . M , F I T H R = c o n s t u = [ w f m , w f a , α f q , δ 1 , δ 2 , A 8 , α f , α c ] T

2.3. The Principle of Improved NSDE-GWO Hybrid Algorithm

The PSC of the aero-propulsion system is confronted with the problem of increasing control variables and the nonlinearity of the objective function, which necessitates that the control algorithm adapt to a wide range of working conditions while maintaining global convergence and computational efficiency. Grey wolf optimization [33] (GWO) is a heuristic algorithm proposed by Mirjalili in 2014. It offers the advantages of fewer adjustment parameters, as well as being simple and efficient. The differential algorithm (DE) has the characteristics of fast convergence and strong robustness. The NSDE-GWO combines the two algorithms to maximize the global search capability and improve the optimization efficiency by the NS evolutionary principles of the PSC.
Based on GWO, the three individuals α , β , γ with the highest fitness value in the historical population are found as the leading class of the wolf pack. By replicating the exploration and surrounding of prey, each generation of the population is updated. The formula for the behavior of encircling prey is as follows:
D = C × X p ( t ) X ( t )
X ( t + 1 ) = X p ( t ) A · D
A = 2 a · r 1 a
C = 2 · r 2
As the number of iterations increases, α decreases linearly from the initial value of 2 to 0. r 1 and r 2 are random vectors between [0,1]. To simulate the hunting behavior of grey wolves, other positions need to be updated according to α , β , γ . The formula for hunting behavior is as follows:
D α = C 1 · X α X D β = C 2 · X β X D δ = C 3 · X δ X
X 1 = X α A 1 · D α X 2 = X β A 2 · D β X 3 = X δ A 3 · D δ
X ( t + 1 ) = X 1 + X 2 + X 3 3
A new generation of population updates based on the grey wolf algorithm is accomplished using the above formula. The GWO population update mostly relies on the historically best three individual values, disregarding other individuals, whereas DE employs a random update technique for all individuals, which allows it to participate in the grey wolf algorithm. This algorithm’s population update must go through the mutation and crossover procedure. The following is the mutation process formula:
v i = x i 1 + F i · x i 2 x i 3 F i = N i 0.5 , 0.5 ,   if   U i 0 , 1 fp δ i ,     otherwise
The formula of the crossover process is as follows:
u i j = v i j ,   if   U j 0 , 1 C R x i j ,   otherwise
The greedy algorithm is used for population selection, and a new generation wolf group of the new generation gray wolf algorithm is constructed based on the DE and the GWO new generation population. The next iteration begins, and so on until the model converges. The greedy algorithm selection procedure has the following formula:
x i = u i ,   if   f u i f x i x i ,   otherwise
Figure 3 depicts the calculation flow of NSDE-GWO. The benefit of DE is that the updating process makes use of all individuals’ information, which is stochastic and global. However, DE optimization efficiency is low in the early stages of iteration, and it is possible to fall into local convergence owing to population aggregation in the late stages of iteration. GWO has high optimization efficiency in the early stages, and the randomness of C ~ is conducive to exploring the global optimal solution in the later stages. The concept of a greedy algorithm is integrated into the population update process, which retains more dominant individuals to engage the following generation, potentially improving the optimization efficiency. Therefore, the NSDE-GWO hybrid algorithm has theoretically better global search capability and optimization efficiency.

2.4. Design of the ADRC Multivariable Controller

In this paper, the ADRC is used as the main closed loop controller in the propulsion system. The main function of ADRC is to use the extended state observer (ESO) to estimate the system disturbance in real time, compensate the disturbance, and improve the robustness. The tracking differentiator (TD) makes the feed-back error change gradually to solve the contradiction between rapidity and over-shoot.
Considering the practical requirements of the engine, the wfa and A8 controlled by sensor-based closed-loop controllers, the turbine total pressure ratio (NPR), and the high-pressure rotor speed ( N h p r ) are chosen as the controlled quantities. The adjustable variable α fq , δ 1 , δ 2 of the supersonic inlet and α f , α c is controlled by an open-loop controller and incorporates a TD to avoid rapid changes in performance parameters. The operating principle of the control system module is shown in Figure 4 [34]. According to the parameter adjustment principle of the ADRC controller proposed by Gao ZQ [35], the parameters are set at ALT = 12.47 km, Ma = 2.0.

3. Results and Discussion

3.1. Comparison of PSC Results of Several Intelligent Algorithms

This research investigates the effects of various intelligent algorithms in the PSC to validate the efficacy of the NSDE-GWO hybrid algorithm. The GWO, DE, DE-GWO, and NSDE-GWO optimization algorithms are all verified. The inlet has a large capture flow under Alt = 12.47 km, Ma = 2.0, but the long inlet size makes the inlet work in the deep subcritical state. The initial population and the number of iteration steps are unchanged in four algorithms to exclude the effects of different parameters. Because of the randomness of the swarm-intelligence algorithm, each algorithm is optimized three times. Table 1 indicates the optimized initial conditions and the range of the optimization search for the control variables.
According to Figure 5a, DE, DE-GWO, and NSDE-GWO may greatly enhance the installed thrust, but GWO is sensitive to local convergence and has an unstable ultimate value of optimization. NSDE-GWO and DE-GWO combine the advantages of two intelligent algorithms and have a faster convergence speed during the iteration process. Furthermore, compared with DE-GWO, NSDE-GWO has faster calculation efficiency and globality, and the installed thrust of the optimization process is increased quickly, allowing for the ideal optimization result to be obtained within a restricted optimization time. To quantify the improvement effect of the installed thrust of different intelligent algorithms, Figure 5b shows the installed thrust optimization potentials of the three intelligent algorithms are similar; both are about 3%. NSDE-GWO has the largest percentage increase in the installed thrust, and the three optimization effects are close and have good stability.
Figure 6a reveals that DE, DE-GWO, and NSDE-GWO reach a lower fuel consumption level in the early stage of the iteration (about 20 steps), but the optimization speed in the later stage of the iteration is significantly reduced due to the influence of the objective function penalty function intervention. The diversity requirement of the population search is high, while the GWO can easily fall into local convergence, and the final value of the optimization is unstable. Figure 6b shows the percentage analysis of the mini-mum installed fuel consumption optimization potential of the three intelligent algorithms. It shows that NSDE-GWO has the highest minimum installed fuel consumption optimization percentage compared with DE-GWO and DE; it can reduce the minimum installed fuel consumption by 0.838%, which is 0.314% and 0.313% higher than DE and DE-GWO, respectively. The analysis shows that NSDE-GWO improves the randomness of the population updating in the later iteration of the model, and it is easier to search for lower fuel-consumption points.
Figure 7a shows the optimization effect of four intelligent algorithms for minimum turbine temperature control. Because the turbine temperature is mainly affected by the working state of the core components, there may be many groups of optimization input variables under the same turbine temperature conditions, which makes the optimization objective function more complex and presents a significant multi-peak state. The four intelligent algorithms do not achieve global convergence in a limited number of iteration steps. However, NSDE-GWO has faster optimization efficiency because it combines the advantages of the other two algorithms. Figure 7b shows the analysis of the lowest turbine temperature potential of the three intelligent algorithms. After 150 steps, the lowest turbine temperature of NSDE-GWO can be reduced by 38.25 k, while DE-GWO and DE can be reduced by 37.75 k and 37.42 k, respectively.
This study carries out a comprehensive analysis of the iterative process in steps 0–100; the algorithm’s iterative process from 0–100 steps is divided into five parts. The data during every 20 steps is used as a data set to draw a box plot. The median and average of each data set are shown in the graph. The 20% quantile is selected for IQR, and the outliers beyond the maximum and minimum observed values are not displayed. Figure 8 shows that the median and average value of GWO is larger than DE in the process of 1–20 steps, indicating that GWO has high optimization efficiency at the beginning of iteration. The volume of the GWO box is smaller than that of DE, which means that the GWO group of individuals gather faster, while that of DE is more scattered. However, the volume of the GWO box was larger and fell into local convergence during the 20–100 iteration steps, whereas the volume of the DE box gradually shrinks and tends to a fixed value, indicating that DE is better than GWO in global convergence. Compared to DE, DE-GWO has a higher median and average value for each data set, particularly in the early iterations. Although the box volume at the beginning of the iteration becomes larger, the average value of optimized thrust increases in 20–100 iteration steps, and the box volume decreases rapidly, demonstrating that the NSDE-GWO has better global convergence and iteration speed. Figure 9 shows that DE, DE-GWO, and NSDE-GWO can gain the ideal low fuel consumption within 1–20 iteration steps. The median, average, and box volume of the NSDE-GWO in steps 1–20 are significantly smaller, resulting in increased efficiency and stability. Figure 10 shows that the volume of the box does not decrease significantly with the increase in the iteration steps, and it is difficult for the four iterative algorithms to achieve global convergence. However, the NSDE-GWO average value is generally lower in 20–100 steps, indicating that the ideal temperature point can be found in limited iterations.
The above analysis demonstrates that a single intelligent algorithm cannot achieve the computational efficiency and accuracy requirement. The NSDE-GWO combines the benefits of the two methods and achieves the best results within all three control modes.

3.2. Analysis of the Performance Mode Transition of the Fusion Controller

In contrast to the previous studies, this study focuses on the analysis of the optimization effect potential brought about by the regulation of the inlet/engine matching at higher Mach number conditions (Ma = 2.0). At the same time, more attention is paid to the smooth switching of the dynamic process fusion controller because the optimization potential varies greatly among engines subject to different constraints and variable ranges. The ADRC closed-loop controller initially simulates the typical acceleration process (Times = 50–60 s) at supersonic conditions and switching from the quasi-steady state to the performance mode (Start time = 100 s). As shown in Figure 11a, the ADRC controller is able to achieve overshoot-free tracking relatively quickly during acceleration. At 100 s, the fusion controller transitions from the quasi-steady state mode to the performance mode, with the implementation of the transition being fed into the ADRC controller via optimized parameters from the PSC module. The maximum variance process during conversion for the dIFN is within 2000 N/s, and thrust smoothing is presented at 103.5 s. This demonstrates that the fusion controller is able to achieve no loss of installed thrust and only minor fluctuations during conversion, as well as obtain maximum thrust performance relatively rapidly. The optimized thrust is increased from 89,658.7 N to 92,147.5 N, about 2.78%, which is slightly below the optimized potential of the steady-state simulation in Section 3.1 and is thought to be caused by the closed-loop controller’s departure from the wfm and A8 dynamic inputs. Figure 11b represents the dynamic change process of multiple open-loop variables; the input process is the change law after the TD module. In Figure 12, X N h p r and XNPR are the controlled variables, wfm and A8 are the control quantities, and Nhc and NPR are the delay curves following the TD module. The TD smooths the target input without creating abrupt changes, and it is tuned using a gradient-optimization approach.
Figure 13 and Figure 14 represent the trends of the performance parameters and optimization variables under the min ISFC mode, and Figure 15 and Figure 16 represent the results under the min T45 mode. Figure 13a and Figure 15a show the results of the dynamic minimum ISFC and minimum T45 modes of the fusion-control system, respectively. Figure 13a reveals that an ISFC reduction of 0.73% can be achieved at 105s, with only a 0.28% reduction in installed thrust, basically achieving the requirement of no loss of thrust performance. The maximum fluctuation in thrust during the transition is within 0.9%, indicating that the fusion controller is able to achieve a smooth transition in the performance mode. Figure 13b indicates that the T45 can drop by 40.11 k, a reduction of 3.01% and a reduction in installed thrust of 0.19%. Compared to the min ISFC mode, the min T45 mode is optimized for high efficiency within Times = 100–101 s. This is due to the reason that the turbine temperature is mainly affected by the adjustment of the fan and compressor guide vane angle, which responds faster than the performance response of the inlet system adjustment. The change pattern of control variables is similar to the max F I T H R mode, which can achieve a faster response and a smoother transition of parameters.
Previous studies on Mach 0.7–1.5 PSC found that the optimization potential gradually decreases as the Mach number increases, with a 5% increase in thrust and a 2% decrease in SFC at Mach 1.4 [36], so the optimization potential of the engine itself will be lower at Mach 2.0 operating conditions, a trend that coincides with the optimization potential results in this paper. The increase in optimization variables leads to greater component potential, and the optimization potential of the lowest turbine temperature is improved compared to previous studies. The experimental study of the supersonic PSC showed a greater optimization potential in terms of thrust and fuel consumption (9%, 8%) due to the elevated amount of enriched fuel and main fuel, while our study kept essentially no increase in the fuel quantity, and the performance optimization potential was achieved mainly through the tuning of accessory components [7].
N f a n , c o r and N c o m , c o r indicate the fan relative corrected speed and the compressor relative corrected speed, respectively. As can be seen from Table 2, the value of N f a n , c o r and N c o m , c o r varies due to the inletmatching and the change of engine operating point. In addition, SMF and SMC are also significantly reduced, and the optimized point works closer to the surge line, suggesting that PSC redistributes the energy of the high- and low-pressure rotors. In the dynamic process Figure 17, the change of engine surge margin also needs to be demonstrated. During the acceleration process, there is a significant decline in the surge margin, but during the performance mode transition process due to the smooth transition of the fusion controller, a trend of rising and then decreasing surge margin is achieved. To sum up, the optimization path of PSC can be realized by simultaneously relying on inlet/engine matching and the guide vane angle adjustment of compression components; the installed thrust optimization can obtain greater optimization potential. The minimum fuel-consumption control and the minimum turbine temperature control can not only rely on the guide vane angle of the compression part but also compensate for the thrust loss caused by reducing the fuel flow through the inlet/engine matching.
In conclusion, the fusion controller is capable of achieving a smooth transition of performance modes to ensure that the engine maintains a stable thrust during the operation and that the dynamic optimization potential is conservative when compared to the ideal steady-state simulation results but has no effect on the overall control effect.

4. Conclusions

To improve the performance potential of the aero-propulsion system controller, the performance-seeking control effects of NSDE-GWO and different swarm-intelligence algorithms are compared, and a fusion-control system using the NSDE-GWO algorithm and ADRC to realize mode conversion is proposed. The following conclusions were obtained:
  • Compared to other intelligent algorithms (DE, GWO, and DE-GWO), the NSDE-GWO hybrid algorithm offers larger potential for optimization, the highest iteration efficiency, and the greatest stability. It can obtain better results in three control modes. Under the conditions of Alt = 12.47 km, Ma = 2.0, the installed thrust is increased by 3.02%, the fuel consumption is reduced by 0.84%, and the temperature after high-pressure turbine is reduced by 3.01%.
  • The fusion controller can produce a seamless transition of performance modes to guarantee that the engine maintains a stable thrust throughout operation; the dynamic optimization potential is conservative when opposed to ideal steady-state simulation outputs. The maximum thrust mode is capable of achieving no thrust drop and a maximum fluctuation rate within 2000 N/s. In the minimum installed fuel consumption mode, the greatest fluctuation in thrust during the changeover is less than 0.9%, and the minimum turbine temperature mode also has smaller volatility and higher optimization efficiency.

Author Contributions

Conceptualization, C.W. and X.D.; methodology, C.W.; validation, C.W.; formal analysis, C.W.; investigation, C.W.; data curation, C.W.; writing—original draft preparation, C.W.; writing—review and editing, C.W. and X.D.; supervision, X.D. and X.S.; project administration, X.S.; funding acquisition, X.S. and X.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by National Natural Science Foundation (NNSF) of China under grants 61890921 and 61903059, and the Aviation Science Foundation of China under Grant 2019ZB063001.

Data Availability Statement

The data used to support the findings of this study are included within the article, which can available from the author [Chen Wang], upon reasonable request.

Conflicts of Interest

We declare that we have no known competing financial and personal relationships with other people or organizations that could have inappropriately influenced our work in this paper. There is no professional or other personal interest of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of, the manuscript entitled.

Abbreviations

The following abbreviations are used in this manuscript:
PSCPerformance-seeking control
ADRCActive disturbance rejection control
ITHRInstalled thrust
SFCSpecific fuel consumption
ISFCInstalled specific fuel consumption
ESOExtended state observer
TDTracking differentiator
GWOGrey wolf optimizer
DEDifferential evolution
NSDENeighborhood-based speciation differential evolution
NPRTurbine total pressure ratio
HPRHigh-pressure rotor
LPRLow-pressure rotor
HPCHigh-pressure compressor
LPTLow-pressure turbine
HPTHigh-pressure turbine
LPLinear programming
SLPSequential linear programming
FSQPFeasible sequential quadratic programming
SQCQPSequential quadratically constrained quadratic programming
GAGenetic algorithm
PSOParticle swarm optimization
EAEvolutionary algorithm
IAInterval analysis

References

  1. Jaw, L.; Mattingly, J. Aircraft Engine Controls; American Institute of Aeronautics and Astronautics: New York, NY, USA, 2009. [Google Scholar]
  2. Montazeri-Gh, M.; Abyaneh, S.; Kazemnejad, S. Hardware-in-the-loop simulation of two-shaft gas turbine engine electronic control unit. Proc. Inst. Mech. Eng. Part I J. Syst. Control. Eng. 2016, 230, 512–521. [Google Scholar]
  3. Wurth, S.; Mahone, T.; Hart, J.; Baxter, J. X-35B integrated flight propulsion control fault tolerance development. In Proceedings of the 2002 Biennial International Powered Lift Conference and Exhibit, Williamsburg, VA, USA, 5–7 November 2002; p. 6019. [Google Scholar]
  4. Landy, R.; Yonke, W.; Stewart, J. Development of HIDEC adaptive engine control systems. In Proceedings of the Turbo Expo: Power for Land, Sea, and Air. American Society of Mechanical Engineers, Dusseldorf, Germany, 8–12 June 1986; Volume 79290, p. V002T02A015. [Google Scholar]
  5. Smith, R.H.; Chisholm, J.D.; Stewart, J.F. Optimizing Aircraft Performance With Adaptive, Integrated Flight/Propulsion Control. J. Eng. Gas Turbines Power 1991, 113, 87–94. [Google Scholar] [CrossRef]
  6. Nobbs, S.; Jacobs, S.; Donahue, D. Development of the full-envelope performance-seeking control algorithm. In Proceedings of the 28th Joint Propulsion Conference and Exhibit, Nashville, TN, USA, 6–8 July 1992. [Google Scholar]
  7. ORME, J.; GILYARD, G. Preliminary supersonic flight test evaluation of performance seeking control. In Proceedings of the 29th Joint Propulsion Conference and Exhibit, Monterey, CA, USA, 28–30 June 1993; p. 1821. [Google Scholar]
  8. Nobbs, S.G. PSC algorithm description. In NASA. Dryden Flight Research Center, An Electronic Workshop on the Performance-Seeking Control and Propulsion Controlled Aircraft Results of the F-15 Highly Integrated Digital Electronic Control Flight Research Program; NASA: Washington, DC, USA, 1995; Volume 95, p. 33013. [Google Scholar]
  9. Ren, X.Y.; Yang, Y.W.; Fan, S.Q. Study on propulsion system overall performance-seeking control algorithm. J. Propuls. Technol. 2010, 31, 62–64. [Google Scholar]
  10. Wang, J.; Zhang, H.; Sun, J.; Huang, X. Experimental verification of aero-engine performance-seeking control based on the hybrid model and FSQP algorithm. J. Propuls. Technol. 2012, 33, 579–590. [Google Scholar]
  11. NIE, Y.; LI, Q.; WANG, Y.; GU, S. Variable cycle engine performance-seeking control based on SQCQP algorithm. J. Beijing Univ. Aeronaut. Astronaut. 2017, 43, 2564. [Google Scholar]
  12. Chen, H.; Zheng, Q.; Gao, Y.; Zhang, H. performance-seeking control of minimum infrared characteristic on double bypass variable cycle engine. Aerosp. Sci. Technol. 2021, 108, 106359. [Google Scholar] [CrossRef]
  13. Li, J.; Fan, D.; Sreeram, V. Optimization of aero engine acceleration control in combat state based on genetic algorithms. Int. J. Turbo Jet Engines 2012, 29, 29–36. [Google Scholar] [CrossRef]
  14. Li, Y.; Ghafir, M.; Wang, L.; Singh, R.; Huang, K.; Feng, X. Nonlinear multiple points gas turbine off-design performance adaptation using a genetic algorithm. J. Eng. Gas Turbines Power 2011, 133, 071701. [Google Scholar] [CrossRef]
  15. Li, Y.; Yuan, H.Q.; Yang, S. Optimization on mistuning blades arrangement of vibration absorption based on genetic particle swarm algorithm in aero-engine. In Proceedings of the Advanced Materials Research. Trans Tech Publ, Xi’an, China, 27–28 July 2013; Volume 655, pp. 481–485. [Google Scholar]
  16. Sun, F.; Miao, L.; Zhang, H. A Study on the Installed performance-seeking control for Aero-Propulsion under Supersonic State. Int. J. Turbo Jet Engines 2015, 2015, 341–351. [Google Scholar] [CrossRef]
  17. Zheng, Q.; Wang, Y.; Sun, F.; Fang, J.; Zhang, H.; Hu, Z. Research on aero-engine steady model based on an improved compact propulsion system model. Proc. Inst. Mech. Eng. Part I J. Syst. Control. Eng. 2020, 234, 726–733. [Google Scholar] [CrossRef]
  18. Zheng, Q.; Fu, D.; Wang, Y.; Chen, H.; Zhang, H. A study on global optimization and deep neural network modeling method in performance-seeking control. Proc. Inst. Mech. Eng. Part I J. Syst. Control. Eng. 2020, 234, 46–59. [Google Scholar] [CrossRef]
  19. Frederick, D.K.; Garg, S.; Adibhatla, S. Turbofan engine control design using robust multivariable control technologies. IEEE Trans. Control. Syst. Technol. 2000, 8, 961–970. [Google Scholar] [CrossRef]
  20. Lehtinen, B.; DeHofft, R.; Hackney, R. Multivariable control altitude demonstration on the FlOO turbofan engine. J. Guid. Control. 1981, 4, 50–58. [Google Scholar] [CrossRef]
  21. Athans, M.; Kapasouris, P.; Kappos, E.; Spang III, H. Linear-quadratic Gaussian with loop-transfer recovery methodology for the F-100 engine. J. Guid. Control. Dyn. 1986, 9, 45–52. [Google Scholar] [CrossRef]
  22. Yang, G.; Sun, J.; Yao, H.; Zang, J.; Liu, A.; Yi, X.; Zhang, L. Experimental verification of H/LTR method for aeroengine control systems. Acta Aeronaut. Astronaut.-Sin.-Ser. 2006, 27, 773. [Google Scholar]
  23. Richter, H. A multi-regulator sliding mode control strategy for output-constrained systems. Automatica 2011, 47, 2251–2259. [Google Scholar] [CrossRef]
  24. Richter, H. Multiple sliding modes with override logic: Limit management in aircraft engine controls. J. Guid. Control. Dyn. 2012, 35, 1132–1142. [Google Scholar] [CrossRef]
  25. Han, J.Q. Auto disturbance rejection controller and it’s applications. Control. Decis. 1998, 13, 19–23. [Google Scholar]
  26. Gao, Z. Scaling and bandwidth-parameterization based controller tuning. In Proceedings of the ACC. Citeseer, Borovets, Bulgaria, 10–12 September 2003; pp. 4989–4996. [Google Scholar]
  27. Miklosovic, R.; Gao, Z. A dynamic decoupling method for controlling high performance turbofan engines. IFAC Proc. Vol. 2005, 38, 532–537. [Google Scholar] [CrossRef]
  28. Li, S.Q.; Zhang, S.X.; Liu, Y.N.; Zhou, S. Simple design and application of auto disturbance rejection controller for aeroengine. Aeroengine 2012, 38, 46–48. [Google Scholar]
  29. Zhang, H.B.; Wang, J.K.; Wang, R.X.; Sun, J.G. Design of an active disturbance rejection decoupling multivariable control scheme for aero-engine. J. Propuls. Technol. 2012, 33, 78–83. [Google Scholar]
  30. Xue, W.; Bai, W.; Yang, S.; Song, K.; Huang, Y.; Xie, H. ADRC with adaptive extended state observer and its application to air–fuel ratio control in gasoline engines. IEEE Trans. Ind. Electron. 2015, 62, 5847–5857. [Google Scholar] [CrossRef]
  31. Jia, L.Y.; Chen, Y.C.; Xie, J.Q.; Cai, F.C. A Simplified Method to Simulate Supersonic Inlet Installed Performance in Terms of Engine and Inlet Matching. J. Propuls. Technol. 2017, 38, 510–518. [Google Scholar]
  32. Wang, C.; Sun, X.; Du, X. The Aero-Engine Component-Level Modelling Research Based on NSDE Hybrid Damping Newton Method. Int. J. Aerosp. Eng. 2022, 13, 1–13. [Google Scholar] [CrossRef]
  33. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar]
  34. Wang, C.; Ai, L.; Du, X.; Sun, X. Multivariable Control and Robustness Analysis of Supersonic Aero-Engine Based on ADRC. In Proceedings of the 5th China Aeronautical Science and Technology Conference, Qingdao, China, 24–26 August 2022; Springer: Berlin/Heidelberg, Germany, 2022; pp. 885–901. [Google Scholar]
  35. Gao, Z. Active disturbance rejection control: A paradigm shift in feedback control system design. In Proceedings of the 2006 American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; IEEE: Piscataway, NJ, USA, 2006; p. 7. [Google Scholar]
  36. Zheng, Q.; Xiang, D.; Chen, C.; Zhang, H. Research on aero-engine performance-seeking control based on the NN-PSM on-board model. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2022, 236, 3505–3517. [Google Scholar]
Figure 1. Intelligent optimization and ADRC fusion controller.
Figure 1. Intelligent optimization and ADRC fusion controller.
Aerospace 10 00151 g001
Figure 2. The main workflow of intet/engine integrated component-level model.
Figure 2. The main workflow of intet/engine integrated component-level model.
Aerospace 10 00151 g002
Figure 3. Flow chart of NSDE-GWO hybrid algorithm.
Figure 3. Flow chart of NSDE-GWO hybrid algorithm.
Aerospace 10 00151 g003
Figure 4. Schematic diagram of the engine’s ADRC control system.
Figure 4. Schematic diagram of the engine’s ADRC control system.
Aerospace 10 00151 g004
Figure 5. (a) Optimization effect of maximum installed thrust control of four intelligent algorithms (DE, GWO, DE-GWO, and NSDE-GWO). (b) Histogram of installed thrust-increase percentage.
Figure 5. (a) Optimization effect of maximum installed thrust control of four intelligent algorithms (DE, GWO, DE-GWO, and NSDE-GWO). (b) Histogram of installed thrust-increase percentage.
Aerospace 10 00151 g005
Figure 6. (a) Optimization effect of minimum installed fuel-consumption control of four intelligent algorithms (DE, GWO, DE-GWO, NSDE-GWO). (b) Histogram of ISFC increase percentage.
Figure 6. (a) Optimization effect of minimum installed fuel-consumption control of four intelligent algorithms (DE, GWO, DE-GWO, NSDE-GWO). (b) Histogram of ISFC increase percentage.
Aerospace 10 00151 g006
Figure 7. (a) Optimization effect of minimum turbine temperature control of four intelligent algorithms (DE, GWO, DE-GWO, and NSDE-GWO). (b) Histogram of T45 increase percentage.
Figure 7. (a) Optimization effect of minimum turbine temperature control of four intelligent algorithms (DE, GWO, DE-GWO, and NSDE-GWO). (b) Histogram of T45 increase percentage.
Aerospace 10 00151 g007
Figure 8. Box plots of different iterative algorithm optimization process in maximum installed thrust control.
Figure 8. Box plots of different iterative algorithm optimization process in maximum installed thrust control.
Aerospace 10 00151 g008
Figure 9. Box plots of different iterative algorithm optimization process in minimum installed specific fuel-consumption control.
Figure 9. Box plots of different iterative algorithm optimization process in minimum installed specific fuel-consumption control.
Aerospace 10 00151 g009
Figure 10. Box plots of different iterative algorithm optimization process in minimum low-turbine inlet temperature control.
Figure 10. Box plots of different iterative algorithm optimization process in minimum low-turbine inlet temperature control.
Aerospace 10 00151 g010
Figure 11. (a) Simulation analysis of incorporating PSC Max F I T H R module and ADRC controller. (b) Dynamic analysis of the open-loop controlled variables α f q , δ 1 , δ 2 , A 8 , α f , α c under Max F I T H R mode.
Figure 11. (a) Simulation analysis of incorporating PSC Max F I T H R module and ADRC controller. (b) Dynamic analysis of the open-loop controlled variables α f q , δ 1 , δ 2 , A 8 , α f , α c under Max F I T H R mode.
Aerospace 10 00151 g011
Figure 12. (a) Dynamic analysis of the controlled variables wfm and N h p r under Max F I T H R mode. (b) Dynamic analysis of the controlled variables A8 and NPR under Max F I T H R mode.
Figure 12. (a) Dynamic analysis of the controlled variables wfm and N h p r under Max F I T H R mode. (b) Dynamic analysis of the controlled variables A8 and NPR under Max F I T H R mode.
Aerospace 10 00151 g012
Figure 13. (a) Simulation analysis of incorporating PSC Min ISFC module and ADRC controller. (b) Dynamic analysis of the open-loop controlled variables α f q , δ 1 , δ 2 , A 8 , α f , α c under Min ISFC mode.
Figure 13. (a) Simulation analysis of incorporating PSC Min ISFC module and ADRC controller. (b) Dynamic analysis of the open-loop controlled variables α f q , δ 1 , δ 2 , A 8 , α f , α c under Min ISFC mode.
Aerospace 10 00151 g013
Figure 14. (a) Dynamic analysis of the controlled variables wfm and N h p r under Min ISFC mode. (b) Dynamic analysis of the controlled variables A8 and NPR under Min ISFC mode.
Figure 14. (a) Dynamic analysis of the controlled variables wfm and N h p r under Min ISFC mode. (b) Dynamic analysis of the controlled variables A8 and NPR under Min ISFC mode.
Aerospace 10 00151 g014
Figure 15. (a) Simulation analysis of incorporating PSC Min T45 module and ADRC controller. (b) Dynamic analysis of the open-loop controlled variables α f q , δ 1 , δ 2 , A 8 , α f , α c under Min T45 mode.
Figure 15. (a) Simulation analysis of incorporating PSC Min T45 module and ADRC controller. (b) Dynamic analysis of the open-loop controlled variables α f q , δ 1 , δ 2 , A 8 , α f , α c under Min T45 mode.
Aerospace 10 00151 g015
Figure 16. (a) Dynamic analysis of the controlled variables wfm and N h p r under Min T45 mode. (b) Dynamic analysis of the controlled variables A8 and NPR under Min T45 mode.
Figure 16. (a) Dynamic analysis of the controlled variables wfm and N h p r under Min T45 mode. (b) Dynamic analysis of the controlled variables A8 and NPR under Min T45 mode.
Aerospace 10 00151 g016
Figure 17. (a) Dynamic analysis of SMF under three modes. (b) Dynamic analysis of SMC under three modes.
Figure 17. (a) Dynamic analysis of SMF under three modes. (b) Dynamic analysis of SMC under three modes.
Aerospace 10 00151 g017
Table 1. Optimization of initial conditions and parameter ranges.
Table 1. Optimization of initial conditions and parameter ranges.
Optimization VariablesInitial DataOptimization Range
Main fuel flow (kg/h)4125.6[3713, 4126]
Afterburner fuel flow (kg/h)16,873.2[15,187, 16,874]
Nozzle area (m2)0.6685[0.50, 0.67]
Fan guide vane angle (deg)0[−10, 10]
Compressor guide vane angle (deg)0[−10, 10]
Relative bleed adjustment (%)0[−10, 20]
First-stage ramp angle (deg)2[1, 5]
Second-stage ramp angle (deg)1[1, 5]
Table 2. Comparison of main parameters of PSC under Alt = 12.47 km Ma = 2.0.
Table 2. Comparison of main parameters of PSC under Alt = 12.47 km Ma = 2.0.
Computational TermOriginal DataOptimized DataOptimization Percentage (%)
Maximum installed thrust control
FITHR(N)89,658.792,147.52.78
SMF0.4970.444−10.66
SMC0.2050.133−35.12
Nfan,cor0.8820.877−0.57
Ncom,cor0.9730.957−1.64
Minimum installed specific fuel-consumption control
FITHR(N)89,658.789,410.63−0.28
ISFC (kg/N.h)2.3432.326−0.73
SMF0.4970.378−23.94
SMC0.2050.144−29.76
Nfan,cor0.8820.863−2.15
Ncom,cor0.9730.952−2.16
Minimum low-pressure turbine inlet temperature control
FITHR(N)89,658.789,492.21−0.19
T45 (K)1332.831292.72−3.01
SMF0.4970.401−19.32
SMC0.2050.126−38.54
Nfan,cor0.8820.824−6.58
Ncom,cor0.9730.954−1.95
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

Wang, C.; Du, X.; Sun, X. Performance-Seeking Control of Aero-Propulsion System Based on Intelligent Optimization and Active Disturbance Rejection Fusion Controller. Aerospace 2023, 10, 151. https://doi.org/10.3390/aerospace10020151

AMA Style

Wang C, Du X, Sun X. Performance-Seeking Control of Aero-Propulsion System Based on Intelligent Optimization and Active Disturbance Rejection Fusion Controller. Aerospace. 2023; 10(2):151. https://doi.org/10.3390/aerospace10020151

Chicago/Turabian Style

Wang, Chen, Xian Du, and Ximing Sun. 2023. "Performance-Seeking Control of Aero-Propulsion System Based on Intelligent Optimization and Active Disturbance Rejection Fusion Controller" Aerospace 10, no. 2: 151. https://doi.org/10.3390/aerospace10020151

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