Next Article in Journal
TDE-Based Adaptive Integral Sliding Mode Control of Space Manipulator for Space-Debris Active Removal
Previous Article in Journal
Atmospheric Disturbance Modelling for a Piloted Flight Simulation Study of Airplane Safety Envelope over Complex Terrain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameter-Matching Algorithm and Optimization of Integrated Thermal Management System of Aircraft

1
School of Aeronautic Science and Engineering, Beihang University, Beijing 100191, China
2
Technical Institute of Physics and Chemistry, Chinese Academy of Sciences, Beijing 100190, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(2), 104; https://doi.org/10.3390/aerospace9020104
Submission received: 4 November 2021 / Revised: 28 January 2022 / Accepted: 12 February 2022 / Published: 15 February 2022
(This article belongs to the Section Aeronautics)

Abstract

:
The integrated thermal management system of aircraft is essential to maintain a suitable environment for the cabin crew and devices. The system is composed of the air-cycle refrigeration subsystem, the vapor-compression refrigeration subsystem, the liquid-cooling subsystem and the fuel-cycle subsystem, which are coupled with each other through heat exchangers. Due to the complex structure and large number of components in the system, it is necessary to design a corresponding parameter-matching algorithm for its special structure and to select the appropriate optimization design method. In this paper, the structure of an integrated thermal management system is analyzed in depth. A hierarchical matching algorithm of system parameters was designed and realized. Meanwhile, a sensitivity analysis of the system was performed, where key parameters were selected. Besides, a variety of optimization algorithms was used to optimize the design calculations. The results show that the particle swarm optimization and genetic algorithm could effectively find the global optimal solution when taking the fuel penalty as the objective function. Furthermore, the particle swarm optimization method took less time.

1. Introduction

With the development of aviation technology, the aircraft is required to carry high-power equipment, electronic warfare, radar and other advanced avionics components. The overall performance of the electromechanical system needs to be greatly improved and the power-to-weight ratio of the system also needs to be increased [1]. The design of the heat dissipation and heat management systems puts forward more stringent requirements [2], where the mismatch of energy consumption and heat sink is one of the most severe contradictions [3,4]. On the one hand, the aircraft is equipped with a large amount of electromechanical equipment, covering energy conversion, transmission, utilization, etc. During the process, the energy losses are converted into heat loads, which substantially increases the heat-dissipation requirements [5,6]. On the other hand, aircrafts mainly use ram air and fuel as heat sinks [7,8], but technical measures such as the extensive use of composite materials in the aircraft structure, the restriction of the ram-air inlet for stealth and the removal of large energy-consuming bleed air further reduce the available heat sink [9], causing a huge gap in the thermal management of the whole aircraft.
In order to design a better integrated thermal management system for the whole aircraft, it is necessary to analyze various heat-generating components on the aircraft and the elements of the thermal management system. Generally, energy-generating components and executive equipment are distributed in different parts of the aircraft. The environment, surface heat-dissipation characteristics, operation conditions and mission profiles must be considered in the thermal load modeling and temperature prediction control [10]. The aircraft electromechanical system includes relatively independent sub-systems, such as electrical system, hydraulic system, fuel system and environmental control system; their transmission and utilization are closely related if analyzed from the perspective of energy and heat flow. Therefore, it is necessary to change the design concept from the optimization of individual subsystems to the optimization of the entire system. Walters et al. [11] established the interface control and specification document for advanced modeling and simulation and realized the analysis and optimization of the system by integrating subsystem models. Further, Donovan [12] used the tip-to-tail model to simulate the overall thermal management system of the whole aircraft, which realizes the optimal design of the system and meets the heat-dissipation requirements of the integrated airborne electromechanical system.
For a simple thermal system, all the characteristic equations of each component can be simultaneously solved to obtain the air state parameters at each point of the system. However, for the integrated thermal management system, the structure is more complicated. The system contains tens of components and state points. Each state point contains temperature, pressure, humidity, specific enthalpy and other parameters [13]. At the same time, the system includes calculation equations for physical parameters of air and coolant, module and component characteristic equations, connection equations and other constrained equations. Some of them are nonlinear equations and the solution is complicated. Therefore, it is difficult to solve them directly by the method of simultaneous equations. In order to ensure the correctness of the calculation and reduce the possibility of errors, it is necessary to design the corresponding calculation process or framework for the special structure of the integrated thermal management system. Researchers have sought various methods to solve this problem. For the thermal management system of electric aircrafts, Freeman et al. [14] designed a process for calculating the total take-off weight and developed a corresponding optimization program, which mainly focuses on the thermal management system and its integration within the propulsion system. Doman [15] determined the optimal thermal cruise conditions for jet aircrafts by deriving and solving sets of ordinary differential equations.
Besides, the thermal management system involves many variables. If the variables are directly used to optimize the design, it leads to excessive calculations and lower calculation efficiency [16]. In addition, it is easier to cause the cost function to fall into a local extremum point during optimization. In this case, it is hard to find the global optimal solution [17]. Jiang et al. [18] carried out a sensitivity analysis for the ACS, obtained key parameters and simplified the optimization design process. The process can greatly reduce the calculation time. Therefore, for integrated thermal management systems, it is necessary to analyze the parameters that have a greater impact during the process based on the parameter-matching calculation. The analysis process is helpful in reducing the complexity of optimization calculation.
To address the difficulty in the parameter matching of the system and to speed up the process of optimization, this paper proposes an optimization process for the aircraft thermal management system with three steps. First, the structure of the integrated thermal management system was analyzed in depth and the multi-level matching calculation of the system parameters was designed and realized. Then, the sensitivity analysis of the system parameters was carried out, where key parameters were chosen to reduce the number of optimization variables. Lastly, a variety of optimization algorithms were used for optimization design calculations, where different numbers of optimization variables were selected to compare the results.

2. Structure of the Integrated Thermal Management System

The integrated environmental control system mainly includes the air-cycle refrigeration subsystem (ACS), the vapor-compression refrigeration subsystem (VCS), the liquid-cooling-cycle subsystem and the fuel-cycle subsystem. The structure of the integrated environmental control system is shown in Figure 1.
The cabin of aircrafts requires air ventilation; thus, it uses an ACS to provide suitable pressure and temperature air. Meanwhile, low-power electronic equipment also uses air as the cooling medium. For high-power electronic equipment, the low-temperature liquid-cooling cycle is used to transfer the heat to the VCS and then to the high-temperature liquid-cooling cycle; the fuel is used as the final heat sink.
These five subsystems are coupled with each other through heat exchangers to achieve the goal of thermal management. The functions of each subsystem are detailed below.
  • The air-cycle refrigeration subsystem
    Due to the relatively small thermal load generated by the cabin, crew and low-power electronic equipment, the ACS can meet the needs of regulating temperature and pressure. The type of ACS is two-wheel bootstrap [19]. The high-temperature air bleeding from the engine compressor enters the subsystem, passes through the primary air heat exchanger and enters the compressor. The air with high temperature and pressure from the compressor enters the secondary heat exchanger and the heat exchanger, which are coupled with the high-temperature liquid-cooling-cycle subsystem, then expands in the turbine. Finally, the air passes through the heat exchanger coupled with the low-temperature liquid-cooling-cycle subsystem and, finally, flows to the cabin at a suitable temperature and pressure.
  • The vapor compression refrigeration subsystem
    The onboard VCS uses R134a as the refrigerant [20]. The refrigerant R134a evaporates in the evaporator, absorbing the heat of the liquid-cooling working fluid in the low-temperature liquid-cooling cycle. The gaseous R134a from the evaporator enters the compressor. By receiving the work of the compressor, the temperature and pressure of the refrigerant increase. After that, the refrigerant enters the condenser. In the condenser, the R134a with high temperature transfers heat to the working fluid of the high-temperature liquid-cooling-cycle subsystem. Subsequently, the temperature of R134a decreases and cools down into a liquid state. Later, the refrigerant flows through the throttle valve, which can be described as an isenthalpic process with its pressure and temperature decreasing. After that, the refrigerant flows back into the inlet of compressor.
  • The low-temperature liquid-cooling-cycle subsystem
    The cooling capacity of the equipment can be generated by the VCS, but it also needs a liquid cycle to transfer the heat between equipment and the VCS. Liquid cooling is usually adopted for the cooling of the circuit board of the high-power integrated electronic equipment cabin. In this subsystem, high-power electronic equipment is directly connected to the low-temperature liquid-cooling cycle. Driven by the pump, the low-temperature liquid-cooling cycle absorbs the heat load from the electronic cabin and transfers the heat to the VCS in the evaporator.
  • The high-temperature liquid-cooling-cycle subsystem
    In the high-temperature liquid-cooling-cycle subsystem, the working fluid is driven by a pump and passes through multiple heat exchangers. The condenser of the VCS, the secondary heat exchanger of the ACS and the liquid-cooled heat exchanger of the fuel cycle are connected in sequence to transfer heat to the fuel-cycle subsystem. The fuel is used as a heat sink to take the heat away.
  • The fuel-cycle subsystem
    In the integrated thermal management system, the fuel-cycle subsystem mainly involves (1) the heat exchanger between the fuel and the high-temperature liquid-cooling cycle; (2) the heat exchanger between the fuel and the thermal load of the secondary power system (SPS), including lubricating oil cycle and hydraulic oil cycle; and (3) the heat exchanger between fuel and ram air. Among them, the heat exchanger No.2 is simplified to take the thermal load as a design variable.

3. System Parameter Matching

3.1. Hierarchical Parameter-Matching Algorithm

The system structure of the integrated thermal management system is typically hierarchical. According to the structural characteristics of the system, the parameter-matching process can be divided into four levels, namely, system level, subsystem level, component level and thermodynamic-property level. The general planning of the parameter-matching calculation process of the integrated thermal management system is shown in Figure 2.
  • System level
    At the system level, the total weight of the system can be calculated according to the system structure configuration. Then, the fuel penalty can be calculated from the temperature, pressure, humidity, specific enthalpy and other parameter values of each state point, combined with engine performance parameters and flight conditions.
  • Subsystem level
    The integrated thermal management system consists of subsystems such as the ACS, the VCS, the liquid-cooling-cycle subsystem and the fuel-cycle subsystem. As the working fluids of each subsystem are not the same, matching calculations need to be performed according to the characteristics of their respective subsystems. After the component calculation is completed, results are returned to the upper system parameter calculation.
    For the ACS, the core part is the turbine-compressor bootstrap component. Therefore, the calculation of this subsystem mainly covers turbines, compressors and a series of heat exchangers related to them, including secondary ram-air heat exchanger and the heat exchanger between the ACS and high-temperature liquid-cooling cycle. For the VCS, the calculation of the high-temperature liquid-cooling cycle and the low-temperature liquid-cooling cycle determine the heat exchange of its evaporator and condenser during system parameter-matching calculations. Then, the capacity of performance (COP) of the VCS can be obtained.
  • Component level
    At the component level, the parameters of air, coolant and fuel at the inlet and outlet can be calculated through the characteristic equations of each component. After the calculation of each component, the results are returned to the subsystem level. In particular, the calculation of components, such as the primary heat exchanger, does not need iteration related to other components or subsystems. The results can be directly passed to the system level after the calculation.
  • Thermodynamic-property level
    In the thermal management system, the parameters of most major components can be directly calculated and obtained without iterative calculations. Iterative calculations are necessary only when the phase-change process is involved and the saturation state needs to be determined. The state points where calculations are required are located in the VCS. At the thermodynamic-property level, only the enthalpy of the fluid is known, while the pressure, temperature and specific entropy are coupled with each other. Therefore, the results of the parameters can be solved by constructing a set of equations among the three parameters.

3.2. Conditions and Assumptions

The system parameter-matching design model can be constructed based on the structure of the integrated thermal management system and analysis of various components of the system. Furthermore, the fuel penalty of the system is calculated as the objective function of the following optimization: For parameter matching during system design, the efficiency of each component is given. During this process, the parameters of air, coolant and fuel at each key point are calculated. Finally, a complete design scheme is obtained.
When designing an integrated thermal management system, the following parameters are usually considered as known:
  • Engine performance parameters, including outlet temperature and pressure parameters of the compressor, bleed air conditions, etc.;
  • Flight parameters, including flight Mach number, time and altitude, to determine the ambient temperature and pressure;
  • Cabin-air-supply parameters, including air-supply flow rate, temperature and pressure;
  • Engine-fuel-supply parameters, including fuel-supply flow rate and fuel temperature;
  • Thermal load parameters, including heat generation of electronic equipment, heat generation of secondary power systems, etc.
In order to simplify the parameter-matching process, the following assumptions were made:
  • Ignoring the humidity of air, dry-air properties were used for calculations [21];
  • Ignoring the dynamic process of fuel tanks and pipes, only steady-state conditions were considered [21];
  • The heat generated by liquid-cooling pumps and fuel pumps were ignored;
  • The resistance of the pipeline was considered in neighboring components;
  • We assumed that the pressure drop ratio of the heat exchanger was constant [22];
  • The thermodynamic properties of the working fluid were considered to be constant [21].

3.3. Parameter-Matching Calculation Process

As analyzed above, the structure of a thermal management system is typically hierarchical. The parameters are matched according to the thermodynamic-property–component–subsystem–system scheme. At the thermodynamic-property calculation level, the corresponding temperature and specific enthalpy are calculated according to the specific heat capacity of each point. At the component calculation level, the calculation of the heat-exchanger efficiency and heat-transfer balance processes are the main concern. At the subsystem calculation level, the parameters are calculated according to the characteristics of each subsystem. For the ACS, the main calculation includes the turbine-compressor component; for the VCS, the calculation of the refrigeration cycle and energy efficiency calculations under a given cooling capacity is important. At the system calculation level, according to the parameters of each key point and component, system-level parameters such as the total mass and fuel penalty of the system can be calculated.
According to the structure of the thermal management system model and the characteristic equations of each component, the framework of the parameter-matching algorithm was designed according to the structural composition and the degree of correlation, which was deemed convenient for the following analysis and the design of the specific algorithm flow. The algorithm framework was divided into several parts and each part was solved independently, but, at the same time, each part has some connection with the outside components, as shown in Figure 3.
As can be seen from Figure 3, despite the complex structure of the thermal management system, it can be divided into several parts through analysis. Each part is composed of one or more components, which are coupled by transferring temperature and heat flow. The data transfer among various parts might be unidirectional or bidirectional. Therefore, the sequence of calculations needs to be considered. For the interior of each part, it is necessary to construct an iterative process to solve the parameters according to the characteristics of the components of each part.
Furthermore, according to the parameter-matching algorithm framework of the above-mentioned preliminary designed thermal management system and the specific characteristic equations of each component, an algorithm flow that is suitable for computer programming was designed as shown in Figure 4.
According to the system parameter-matching calculation process, the specific calculation steps are as follows:
  • p 16 is calculated according to the pressure of the cabin supply air p 17 .
    p 16 = p 17 1 Δ p h
  • The outlet temperature and pressure of the heat exchanger HX1r1 ( T 12 and p 12 ) are calculated according to the temperature and pressure of the bleed air [23].
    T 12 = T 11 η H X 1 T 11 T r a m
    p 12 = p 11 1 Δ p h
  • The outlet temperature of the pump in the low-temperature liquid-cooling cycle is assumed as T 31 .
    (a)
    Ignoring the temperature change of the fluid in the pump, T 34 = T 31 .
    (b)
    The circulating liquid absorbs the heat from the electronic equipment and the temperature rises.
    T 32 = T 31 + Q d e v q m 3 c p 3
    (c)
    According to the energy balance and the formula of heat-exchanger efficiency, the air-side inlet temperature T 16 of the heat exchanger HX13 is calculated.
    T 16 = T 17 η H X 13 T 32 1 η H X 13
    (d)
    Assuming the turbine pressure ratio π t , we perform the following iterative calculation:
    • Calculate the turbine inlet pressure and temperature, p 15 and T 15 [24].
      p 15 = p 16 · π t
      T 15 = T 16 1 η t 1 π t 0.286
    • Calculate the compressor outlet temperature T 13 , according to the energy balance.
      c p ( T 13 T 12 ) = c p ( T 15 T 16 )
    • Calculate the pressure ratio of compressor and recalculate the turbine expansion ratio [24].
      π c = T 13 T 12 1 η c + 1 3.5
      π t = p 12 p 16 π c 1 Δ p h 2
    • Compare the result with the assumed turbine pressure ratio; return to recalculate if not matched.
    (e)
    Calculate the air-side outlet temperature T 14 of the heat exchanger HX1r2.
    T 14 = T 13 η H X 1 r 2 T 13 T r a m
    (f)
    Calculate the liquid-side outlet temperature of the heat exchanger HX14 ( T 43 and T 44 ).
    T 43 = T 14 T 14 T 15 η H X 14
    T 44 = T 43 + q m c p 1 q m c p 4 T 14 T 15
    (g)
    Ignoring the temperature change of the liquid-cooling pump, T 41 = T 44 .
    (h)
    Calculate the fuel temperature T 52 before absorbing the heat generated by the secondary energy system.
    T 52 = T 54 Q S P S q m c p 5
    (i)
    Calculate the fuel temperature after the ram-air heat exchanger ( T 55 ).
    T 55 = T 54 η H X 5 r T 54 T r a m q m c p 5 , r a m q m c p 5
    (j)
    Recalculate the inlet temperature of the coolant side of the heat exchanger HX45 ( T 41 ), according to the conservation of energy.
    T 41 = q m c p 5 T 52 T 51 q m c p 4 + T 42
    (k)
    Compare T 41 and T 41 ; return to assume T 31 again if not matched.
  • Calculate the heat exchange of the evaporator.
    Q c d s = q m c p 4 T 43 T 42
  • Calculate the heat exchange of the condenser.
    Q e v p = q m c p 3 T 33 T 34
  • Calculate the COP of the VCS.
    C O P = 1 Q c d s Q e v p 1
  • Assume the outlet temperature T 21 of the evaporator.
    (a)
    Calculate the specific enthalpy h 21 , specific entropy s 21 and pressure p 21 of point 21, according to the thermodynamic properties of the refrigerant.
    (b)
    Calculate h 22 and h 23 according to the heat exchange of evaporator and condenser.
    h 23 = h 21 Q e v p / q m 2
    h 22 = h 23 + Q c d s / q m 2
    (c)
    Calculate T 22 and p 22 according to the thermodynamic properties of the refrigerant.
    (d)
    Let p 23 = p 22 ; recalculate h 23 according to p 23
    (e)
    Compare h 23 and h 23 ; return and assume T 21 again if not matched.
  • Calculate the efficiency of the evaporator.
    η e v p = h 21 h 24 h 21 h 33 · q m 2 q m , m i n
  • Calculate the efficiency of the condenser.
    η c d s = h 22 h 23 h 22 h 42 · q m 2 q m , m i n
In the calculation process described above, in order to ensure the completeness of the parameter calculation process, the thermodynamic properties of the working fluids need to be known in advance.
Besides, the following 11 design parameters need to be selected:
  • η H X 1 r 1 : Efficiency of ACS’s primary HX
  • η H X 1 r 2 : Efficiency of ACS’s secondary HX
  • η H X 13 : Efficiency of air/low-temperature liquid-cooling HX
  • η H X 14 : Efficiency of air/high-temperature liquid-cooling HX
  • η H X 45 : Efficiency of fuel/high-temperature liquid-cooling HX
  • η 5 r : Efficiency of fuel/ram-air HX
  • η c : Efficiency of ACS’s compressor
  • η t : Efficiency of ACS’s turbine
  • q m 3 : Flow rate of low-temperature liquid-cooling cycle
  • q m 4 : Flow rate of high-temperature liquid-cooling cycle
  • q m 5 : Flow rate of fuel cycle
According to the results of parameter-matching calculations, other system parameters can be calculated and the fuel penalty of the whole system can be obtained after parameter matching is finished. The results can provide data for sensitivity analyses and optimization calculations.

3.4. Results of Parameter Matching

The parameter-matching process was realized using MATLAB code.
For a certain aircraft, there are several states during flight. Usually, the ground maximum state corresponds to the harshest conditions and the cruise may last the longest time. Therefore, to estimate the performance of the integrated thermal management system, we should pay attention to these two conditions. The parameters and flight conditions are shown in Table 1.
For comparison, we also developed codes in Engineering Equation Solver (EES), which is commercial software for solving common equations. The pressure and temperature results for the two conditions are shown in Table 2. We can see from the table that the results calculated from MATLAB and EES are very close. However, EES required a good initial guess for every variable.
For better intuition, the temperature results were also marked in the system scheme diagrams, as shown in Figure 5 and Figure 6.
Under cruise conditions, the supply air pressure is lower; thus, it needs a higher pressure ratio of the turbine, while the cooling capacity of ACS is higher than that under ground conditions. Therefore, the flow rate of the liquid-cooling cycle and the power of VCS are also lower.
The results show that the algorithm could perform parameter matching with a wide range of input data, covering different stages of flight.

3.5. Fuel Penalty Calculation

The total weight of the aircraft increases after it is equipped with a thermal management system. Besides, ram air, bleed air and power extracted from the engine also cause energy loss. To quantitatively describe the impact of these factors on aircraft flight performance, they can be converted into the increment of fuel consumption, which is also called the fuel penalty method. In this paper, fuel penalty is used as the objective function for further optimization.
The fuel penalty m f is comprised of four parts, i.e., system device mass m f , d e v , ram-air resistance m f , r a m , engine power extraction m f , P and engine bleed air m f , b l .
m f = m f , d e v + m f , r a m + m f , P + m f , b l
The detailed calculation process of each part is shown in Appendix A.

3.6. Sensitivity Analysis of Fuel Penalty

According to the parameter-matching calculation process and the thermal management system parameter-matching algorithm framework (Figure 3), the analysis of variance (ANOVA) was used to calculate the influence of the design variables on fuel consumption for the sensitivity analysis, as shown in Figure 7.
The results show that the ACS’s primary heat-exchanger efficiency, fuel-cycle flow rate, ACS’s secondary heat-exchanger efficiency, high-temperature liquid-cooling-cycle flow rate, low-temperature liquid-cooling-cycle flow rate and fuel/ram-air heat-exchanger efficiency had a greater impact on the fuel penalty. The total contribution rate of the parameters exceeds 0.95.
The heat exchangers contribute a large part of the system weight and their efficiency is important to the weight; thus, these parameters are more sensitive. Besides, since the mass flow rate of the subsystems is important to the heat flow of the heat exchangers, they also play a key role in weight calculation.

4. Optimal Design Calculation

4.1. Optimization Algorithm

Due to the high non-linearity of the parameter-matching calculations of the thermal management system, its fuel penalty value presents multi-peak or multi-valley characteristics. Therefore, the classical gradient-based optimization algorithm is easily affected by the initial value and falls into a local extreme point rather than the global optimal value. Under different initial conditions, based on the gradient-based method, there are large deviations during the optimization calculation process. In order to avoid the excessive influence of the initial value and make the optimization calculation results more accurate, this paper adopted the pattern search algorithm, genetic algorithm and particle swarm algorithm to optimize the design of the integrated thermal management system.

4.1.1. Pattern Search Algorithm

The calculation process of the pattern search algorithm is composed of two moving processes, namely, exploratory search and pattern move. Among them, the exploratory search is used to search for the general direction of the decline in the objective function to improve the accuracy of the search direction; the pattern move is to accelerate the movement according to the direction of the function change and find a better point [25]. The mathematical description of the searching process [26] is shown below.
For an n dimension function f ( x ) , set the initial point x ( 0 ) = x 1 ( 0 ) , x 2 ( 0 ) , , x n ( 0 ) , step size α i > 0 , base vector e i and the minimum calculation error ϵ . For the exploratory move, at the k-th step, x ^ ( k , i ) means the point after i-th exploration; therefore,
x ^ ( k , i ) = x ( k ) i = 0 a ( k , i ) = x ^ ( k , i 1 ) + α i e i f ( a ( i , i ) ) < f ( x ^ ( k , i 1 ) ) b ( k , i ) = x ^ ( k , i 1 ) α i e i f ( a ( i , i ) ) < f ( x ^ ( k , i 1 ) ) f ( b ( k , i ) ) x ^ ( k , i 1 ) otherwise
After an exploratory move, if x ^ ( k , n ) = x ^ ( k , 0 ) , the step size α should be reduced and the process repeated again until x ^ ( k , n ) x ^ ( k , 0 ) or α < ϵ .
For the pattern move with step λ , the direction can be written as D ( k ) = x ^ ( k , n ) x ( k ) ; then, the new point y ^ ( k + 1 ) can be calculated as
y ^ ( k + 1 ) = x ^ ( k , n ) + λ D ( k )
If f ( y ^ ( k + 1 ) ) < f ( x ^ ( k , n ) ) , this round of iteration is finished and x ( k + 1 ) = y ^ ( k + 1 ) . Otherwise, the step size should be reduced and the process repeated.
By combining a series of exploratory searches and pattern moves, the algorithm can realize the optimization calculation of the objective function. The pattern search algorithm is relatively simple to implement and has high computational efficiency [27].

4.1.2. Genetic Algorithm

The genetic algorithm is a random search algorithm designed according to the natural selection and genetic mechanisms of the biological world.
Unlike traditional search algorithms, genetic algorithms start the search process from a set of randomly generated initial solutions, called a population. Each individual in the population is a solution, called a chromosome. These chromosomes evolve in subsequent iterations, called genetic. The algorithm optimizes the calculations by simulating the crossover, mutation and reproduction phenomena in natural selection and genetic process [28].
Crossover or mutation operations generate the next generation of chromosomes, called progeny. The goodness of chromosomes is measured by fitness. A certain number of individuals are selected from the previous generation and the offspring according to the fitness, as the next generation population; then, the evolution continues. After several generations, the algorithm converges on the best chromosome, which is likely to be the optimal solution.
The concept of fitness is used in the genetic algorithm to measure the degree to which each individual in the population is likely to reach the optimal solution in the optimization computation. The definition of the fitness function is generally related to the specific solution problem.
In each round of calculations, a set of candidate solutions are retained; then, better individuals are selected from them and these individuals are combined using operations to form a new generation of candidate program groups. The above process has to be repeated until the convergence condition of the solution is met. The genetic algorithm can realize distributed information collection and exploration in the entire solution space and can avoid the solution process falling into the local extreme point. The flow diagram of the genetic algorithm [29] is shown in Figure 8.

4.1.3. Particle Swarm Algorithm

The particle swarm algorithm is a kind of random algorithm that simulates the group behavior of flocks or groups of fish in nature. These populations adopt a cooperative mode when searching for food and each particle (individual) adjusts its search strategy according to its own and surrounding experience [30], so that the entire population can reach the goal faster.
The algorithm is population-based and moves individuals to better regions based on their fitness to the environment. However, instead of using evolutionary operators for individuals, it treats each individual as a volumeless particle (point) in a N-dimensional search space, flying at a certain speed in the search space, which is dynamically adjusted according to its own flight experience and the flight experience of its peers. The i-th particle is denoted as X i = x i 1 , x i 2 , , x i N and the best position it has experienced (with the best adaptation value) is denoted as P i = p i 1 , p i 2 , , p i N , also known as pbest. The index number of the best position experienced by all particles in the population is denoted by the symbol g, i.e., P g , also known as gbest. The velocity of particle i is denoted by V i = v i 1 , v i 2 , , v i N . For each generation, its ( n + 1 ) -st dimension ( 1 n + 1 N ) varies according to the following equation:
v i , n + 1 = w · v i , n + c 1 · rand ( ) · ( p i , n x i , n ) + c 2 · rand ( ) · ( p g , n x i , n )
x i , n + 1 = x i , n + v i , n + 1
where w is the inertia weight, c 1 and c 2 are acceleration constants and rand() generates random values varying in the range [ 0 , 1 ] .
In addition, the velocity V i of the particle is limited by a maximum velocity V m a x . If the current acceleration of the particle causes its velocity v i , n in a dimension to exceed the maximum velocity v m a x , n , the velocity is limited to it.
In Equation (27), the first part is the inertia of the particle’s previous behavior; the second part is the “cognition” part, which represents the thinking of the particle itself; and the third part is the “social” part, which represents the information sharing and mutual cooperation between the particles.
The particle swarm algorithm uses the following psychological assumption: In the process of seeking consistent cognition, individuals tend to remember their own beliefs and consider the beliefs of their colleagues at the same time. When they perceive that their colleagues’ beliefs are better, they adapt.
The algorithmic process of the standard particle swarm algorithm is shown as follows:
  • Initialize a group of particles (population size, m), including random positions and velocities;
  • Evaluate the fitness of each particle;
  • For each particle, compare its adaptation value with that of the best position pbest it has experienced and, if it is better, use it as the current best position pbest;
  • For each particle, compare its adaptation value with that of the best position gbest experienced globally and reset the index number of gbest if it is better;
  • Vary the velocity and position of the particles according to Equations (27) and (28);
  • If the end condition is not reached (usually a sufficiently good adaptation value or reaching a preset maximum algebra G m a x ), return to step 2.

4.2. Selection of Optimal Variables

The variables were chosen to optimize the value of fuel penalty, which is based on the parameter-matching calculation of the thermal management system. The ranges of the variables are shown in Table 3.
The optimization was performed twice within the scope of the design, as follows:
  • All 11 design variables were taken as optimized design variables to perform optimization calculations;
  • According to the calculation results of the sensitivity analysis, six variables, i.e., (1) efficiency of ACS’s primary HX, (2) flow rate of fuel cycle, (3) efficiency of ACS’s secondary HX, (4) flow rate of high-temperature liquid-cooling cycle, (5) flow rate of low-temperature liquid-cooling cycle and (6) efficiency of fuel/ram-air HX, affected the fuel penalty significantly. Therefore, these variables were set as optimization design variables to perform optimization calculations. The remaining five variables were set as constants.

4.3. Optimization Calculation Results

Using the global optimization toolbox in MATLAB, we developed the code for the optimization algorithms. According to the maximum and minimum values of each variable defined in Table 3, pattern search, genetic algorithm and particle swarm algorithm were used to optimize the 11-variable scheme and the simplified 6-variable scheme. The ground conditions were chosen as the design conditions. The initial optimization variables and the optimization results obtained by the three optimization algorithms are summarized in Table 4. The optimization results and calculation time for different variables and optimization algorithms are shown in Figure 9. The 6-variable optimization calculations consumed less time than those of the 11-variable optimization. Except for the 11-variable pattern search algorithm, the optimized fuel penalty results are not significantly different.

4.4. Discussion

The particle swarm algorithm and genetic algorithm are intelligent optimization algorithms, which avoid the multi-peak or valley phenomena of the objective function. They performed well and could find the global optimal value. The results show that the difference between the minimum fuel penalty is small. Furthermore, the results of the optimized variables are relatively close.
Pattern search is an improved gradient-based optimization algorithm. The algorithm overcame the shortcomings of traditional optimization algorithms to a large extent, which performed well in the optimization process of 6 variables. Besides, the algorithm basically achieved the global optimal results. However, during the process of optimizing 11 variables, the variables still fell into a local minimum, where no global optimal solutions were found. Therefore, in the case where the number of variables was higher, there were still certain limitations related to the pattern search optimization algorithm.
The calculation results of the 6 optimized variables selected by the variable sensitivity analysis are generally close to the calculation results of all 11 variables. Besides, the drop in the number of variables greatly reduced the calculation time.
In general, the pattern search algorithm obtained better results with the least run time in the six-variable calculations. Although the particle swarm algorithm took a longer time, it performed well in the optimization calculation process of both the 11 variables and 6 variables. Additionally, the optimal value of each variable could be obtained during the process. Among the three optimization algorithms, the genetic algorithm took the longest time.
Therefore, it is recommended to comprehensively consider the trade-off of consumed time and calculation accuracy when optimizing the integrated thermal management system. The six-variable pattern search or particle swarm optimization algorithm is the first choice.

5. Conclusions

This paper studies the integrated thermal management system. We developed a whole process for optimizing the system, including the hierarchical parameter-matching algorithm, the sensitivity analysis of the system parameters and optimization design calculations. The main conclusions are reported below.
  • The composition characteristics of the integrated thermal management system were analyzed in depth. Based on the analysis, a four-level parameter-matching calculation method of system, subsystems, components and thermodynamic properties was designed.
  • The parameter sensitivity analysis was conducted based on parameter-matching calculations. Through the parameter-matching process and system thermodynamic analysis, 11 design variables were initially selected. The 11 design variables were reduced to 6 variables by the sensitivity analysis. The six variables are (1) efficiency of ACS’s primary HX, (2) flow rate of fuel cycle, (3) efficiency of ACS’s secondary HX, (4) flow rate of high-temperature liquid-cooling cycle, (5) flow rate of low-temperature liquid-cooling cycle and (6) efficiency of fuel/ram-air HX.
  • The optimization of the integrated thermal management system is complicated because of the multi-extremum problem. The variables optimized by traditional gradient-based optimization algorithms easily fall into the local extremum point. Therefore, the pattern search, genetic algorithm and particle swarm algorithm were selected for optimization calculation and comparative analyses. The results show that the particle swarm optimization and genetic algorithm could find the global optimal value more easily. Among them, the particle swarm algorithm took the lowest amount of time.
  • The algorithm in this paper is mainly aimed at the structural design of the integrated thermal management system of aircrafts. However, it could also be applied to other fields. For thermal systems with similar hierarchical structures, the same idea can be used for parameter matching, sensitivity analyses and optimization design calculations.

Author Contributions

Conceptualization, R.W., S.D. and H.J.; methodology, H.J.; validation, P.L., H.Z. and H.J.; formal analysis, H.J.; investigation, R.W.; resources, S.D.; data curation, P.L. and H.Z.; writing—original draft preparation, R.W. and P.L.; writing—review and editing, R.W. and P.L.; visualization, H.J.; supervision, S.D.; project administration, S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research study received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ACSAir-cycle refrigeration (sub)system
COPCoefficient of performance
HXHeat exchanger
NTUNumber of transfer units
SPSSecondary power (sub)system
TMSThermal management system
VCSVapor-compression refrigeration (sub)system

Nomenclature

The following variables are used in this manuscript:
C e Thrust-specific fuel consumption
c p Specific heat capacity
DDiameter
fFuel
gGravity
hEnthalpy
H v Fuel calorific value
KLift-to-drag ratio
lLength
mMass
n d Rotational speed
pPressure
Δ p Pressure drop/pressure difference
PPower
QHeat
q m Mass flow rate
q V Volumetric flow rate
TTemperature
uCircumferential velocity
vVelocity
x 0 Velocity ratio of turbine
ϵ c Coefficient of combustion completeness
ζ Ratio of the heat capacity rate
η Efficiency
π Pressure ratio
τ 0 Time of flight
Ψ H X Coefficient of mass of heat exchanger

Appendix A. Fuel Penalty Calculation

The fuel penalty is comprised of four parts, system device mass, ram-air resistance, engine power extraction and engine bleed air. In a certain flight time τ 0 , each part can be calculated as shown in the following subsections [22].

Appendix A.1. Fuel Penalty of System Device Mass

The fuel penalty of the system device mass can be calculated as follows:
m f , E = m E exp C e τ 0 g K 1
where m E is the total mass of system devices, C e is the thrust-specific fuel consumption, g is gravity and K is the lift-to-drag ratio.
For some devices, their masses change little under different design conditions. However, for heat exchangers and turbo machines, the masses are related to their efficiencies, powers, pressure differences or other conditions, which can be calculated as follows [22]:
m H X = C H X Ψ H X x 1 N T U x 2 a + ζ l 1 l 2 x 3 · q m
m t u r b = C t u r b · D t u r b 3
where C H X , x 1 , x 2 , x 3 and C t u r b are coefficients; l 1 and l 2 are the lengths of the heat exchanger; N T U is the number of transfer units; Ψ H X is a variable related to inlet pressure and temperature; ζ is the ratio of the heat capacity rate; D t u r b is the diameter of the turbomachine. For an aluminum crossflow heat exchanger onboard, C H X = 0.77 , x 1 = 0.317 , x 2 = 1.317 , x 3 = 0.635 , C t u r b = 21.5 . Other coefficients can be calculated as follows:
Ψ H X = p i n / 101   kPa 2 T i n / 288   K · Δ p p i n
ζ = q m c p max q m c p min
D t u r b = u π n d = x 0 2 c p T i n 1 π t u r b 0.286 π n d
where u is the circumferential velocity of the turbine, n d is the rotational speed and x 0 is the velocity ratio.

Appendix A.2. Fuel Penalty of Ram Air

The fuel penalty of ram air, which is mainly due to momentum resistance, can be calculated as follows:
m f , r a m = q m v K g exp C e τ 0 g K 1
where q m is the mass flow rate of ram air and v is the velocity of the aircraft.

Appendix A.3. Fuel Penalty of Engine Power Extraction

The fuel penalty of engine power extraction can be calculated as follows:
m f , P = q m , f , P K C e g exp C e τ 0 g K 1
where q m , f , P is the additional fuel consumption. For turbojet engines, it can be calculated as
q m , f , P = 2.94 c p , f T 3 H v ϵ c T 1 π e n g i n e 0.286 1 · P
where T 1 and T 3 are the inlet temperatures of engine compressor and turbine; H v is the fuel calorific value; ϵ c is the coefficient of combustion completeness; π e n g i n e is the pressure ratio of the engine compressor; and P is the power extracted from the engine.
The power extracted from the engine is converted into electric power and provided to various devices. Among these devices, avionics can be considered as rated power consumption. However, the power consumed by the pumps should be calculated, since it varies under different conditions.
P p u m p = q V Δ p η p u m p
where q V is the volumetric flow rate, Δ p is the pressure difference and η p u m p is the efficiency of the pump.

Appendix A.4. Fuel Penalty of Engine Bleed Air

The fuel penalty of engine bleed air can be calculated as follows:
m f , b l = q m , f , b l K C e g exp C e τ 0 g K 1
where q m , f , b l is the additional fuel consumption, which can be calculated as
q m , f , b l = D · C e + q m , b l c p , f T 3 H v ϵ c
where D is the aerodynamic resistance and q m , b l is the mass flow rate of bleed air.

References

  1. Dooley, M.; Lui, N.; Newman, R.; Lui, C. Aircraft Thermal Management-Heat Sink Challenge; Technical Report 2014-01-2193, SAE Technical Paper; SAE: Warrendale, PA, USA, 2014. [Google Scholar]
  2. Newman, R.W.; Dooley, M.; Lui, C. Efficient propulsion, power, and thermal management integration. In Proceedings of the 49th AIAA/ASME/SAE/ASEE Joint Propulsion Conference, San Jose, CA, USA, 14–17 July 2013; p. 3681. [Google Scholar]
  3. Nuzum, S.R. Aircraft Thermal Management using Liquefied Natural Gas. Master’s Thesis, Ohio Northern University, Ada, OH, USA, 2014. [Google Scholar]
  4. Morioka, N.; Saito, H.; Takahashi, N.; Seta, M.; Oyori, H. Thermal Management System Concept with an Autonomous Air-Cooled System; Technical Report, SAE Technical Paper; SAE: Warrendale, PA, USA, 2014. [Google Scholar]
  5. Wang, J.; Li, Y.; Liu, X.; Shen, C.; Zhang, H.; Xiong, K. Recent active thermal management technologies for the development of energy-optimized aerospace vehicles in China. Chin. J. Aeronaut. 2020, 34, 1–27. [Google Scholar] [CrossRef]
  6. Mehta, J.; Charneski, J.; Wells, P. Unmanned aerial systems (UAS) thermal management needs, current status, and future innovations. In Proceedings of the 10th International Energy Conversion Engineering Conference, Atlanta, GA, USA, 30 July–1 August 2012; p. 4051. [Google Scholar]
  7. Kellermann, H.; Lüdemann, M.; Pohl, M.; Hornung, M. Design and Optimization of Ram Air–Based Thermal Management Systems for Hybrid-Electric Aircraft. Aerospace 2021, 8, 3. [Google Scholar] [CrossRef]
  8. Rheaume, J.; Lents, C.E. Design and simulation of a commercial hybrid electric aircraft thermal management system. In Proceedings of the 2018 AIAA/IEEE Electric Aircraft Technologies Symposium, Cincinnati, OH, USA, 9–11 July 2018; p. 4994. [Google Scholar]
  9. Miller, C.; Maser, A.; Garcia, E.; Mavris, D. INVENT surrogate modeling and optimization of transient thermal responses. In Proceedings of the 50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Nashville, TN, USA, 9–12 January 2012; p. 1123. [Google Scholar]
  10. Moir, I.; Seabridge, A. Environmental Control Systems. In Aircraft Systems: Mechanical, Electrical, and Avionics Subsystems Integration; John Wiley & Sons: Hoboken, NJ, USA, 2011; Chapter 7; pp. 259–263. [Google Scholar]
  11. Walters, E.; Amrhein, M.; O’Connell, T.; Iden, S.; Lamm, P.; Yerkes, K.; Wolff, M.; McCarthy, K.; Raczkowski, B.; Wells, J.; et al. INVENT modeling, simulation, analysis and optimization. In Proceedings of the 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 4–7 January 2010; p. 287. [Google Scholar]
  12. Donovan, A. Vehicle Level Transient Aircraft Thermal Management Modeling and Simulation. Ph.D. Thesis, Wright State University, Dayton, OH, USA, 2016. [Google Scholar]
  13. Chen, L.; Zhang, X.; Wang, C.; Yang, C. A novel environmental control system facilitating humidification for commercial aircraft. Build. Environ. 2017, 126, 34–41. [Google Scholar] [CrossRef]
  14. Freeman, J.; Osterkamp, P.; Green, M.; Gibson, A.; Schiltgen, B. Challenges and opportunities for electric aircraft thermal management. Aircr. Eng. Aerosp. Technol. Int. J. 2014, 86, 519–524. [Google Scholar] [CrossRef]
  15. Doman, D.B. Optimal cruise altitude for aircraft thermal management. J. Guid. Control Dyn. 2015, 38, 2084–2095. [Google Scholar] [CrossRef]
  16. Pérez-Grande, I.; Leo, T.J. Optimization of a commercial aircraft environmental control system. Appl. Therm. Eng. 2002, 22, 1885–1904. [Google Scholar] [CrossRef]
  17. Vargas, J.V.; Bejan, A. Thermodynamic optimization of finned crossflow heat exchangers for aircraft environmental control systems. Int. J. Heat Fluid Flow 2001, 22, 657–665. [Google Scholar] [CrossRef]
  18. Jiang, H.; Sujun, D.; Zhang, H.; Fengming, A.; Zhang, Z.; Jun, W. Optimization on conventional and electric air-cycle refrigeration systems of aircraft: A short-cut method and analysis. Chin. J. Aeronaut. 2020, 33, 1877–1888. [Google Scholar] [CrossRef]
  19. Ashford, R.; Brown, S. F-22 Environmental Control System/Thermal Management System (ECS/TMS) Flight Test Program-Downloadable Constants, an Innovative Approach; Technical Report 2000-01-2265, SAE Technical Paper; SAE: Warrendale, PA, USA, 2000. [Google Scholar]
  20. Ghanekar, M. Vapor Cycle System for the F-22 Raptor; Technical Report 2000-01-2265, SAE Technical Paper; SAE: Warrendale, PA, USA, 2000. [Google Scholar]
  21. Yang, H.; Yang, C.; Zhang, X.; Xiugan, Y. Influences of Different Architectures on the Thermodynamic Performance and Network Structure of Aircraft Environmental Control System. Entropy 2021, 23, 855. [Google Scholar] [CrossRef] [PubMed]
  22. Shou, R.; He, H. Aircraft Environmental Control; Beihang University Press: Beijing, China, 2003. [Google Scholar]
  23. Yu, J.; Gao, H.; Xie, Y. Calculation of Heat Exchanger Heat Transfer and Resistance. In Heat Exchanger Principle and Design; Beihang University Press: Beijing, China, 2019; Chapter 2; pp. 9–37. [Google Scholar]
  24. Jennions, I.; Ali, F.; Miguez, M.E.; Escobar, I.C. Simulation of an aircraft environmental control system. Appl. Therm. Eng. 2020, 172, 114925. [Google Scholar] [CrossRef]
  25. Chinneck, J.W. Pattern Search for Unconstrained NLP. In Practical Optimization: A Gentle Introduction; Carleton University: Ottawa, ON, Canada, 2015; Chapter 17; pp. 1–4. [Google Scholar]
  26. Tang, C.; Mo, H.; Liu, S. The application of the pattern search method for slope stability analysis. J. South China Univ. Technol. (Nat. Sci. Ed.) 2000, 28, 42–46. [Google Scholar]
  27. Al-Sumait, J.; Al-Othman, A.; Sykulski, J. Application of pattern search method to power system valve-point economic load dispatch. Int. J. Electr. Power Energy Syst. 2007, 29, 720–730. [Google Scholar] [CrossRef] [Green Version]
  28. Ma, Y.; Yun, W. Research progress of genetic algorithm. Appl. Res. Comput. 2012, 29, 1201–1206. [Google Scholar]
  29. Pachori, S. Survey Report on Stock Selection and Portfolio Optimization; Technical Report 12BCE055; Institute Of Technology, Nirma University: Ahmedabad, India, 2015. [Google Scholar]
  30. Wang, D.; Tan, D.; Liu, L. Particle swarm optimization algorithm: An overview. Soft Comput. 2018, 22, 387–408. [Google Scholar] [CrossRef]
Figure 1. Structure of the integrated thermal management system.
Figure 1. Structure of the integrated thermal management system.
Aerospace 09 00104 g001
Figure 2. Levels of parameter matching of integrated thermal management system.
Figure 2. Levels of parameter matching of integrated thermal management system.
Aerospace 09 00104 g002
Figure 3. Structure of parameter matching algorithm of integrated thermal management system.
Figure 3. Structure of parameter matching algorithm of integrated thermal management system.
Aerospace 09 00104 g003
Figure 4. Flow chart of parameter matching of integrated thermal management system.
Figure 4. Flow chart of parameter matching of integrated thermal management system.
Aerospace 09 00104 g004
Figure 5. Temperature results under ground conditions.
Figure 5. Temperature results under ground conditions.
Aerospace 09 00104 g005
Figure 6. Temperature results under cruise conditions.
Figure 6. Temperature results under cruise conditions.
Aerospace 09 00104 g006
Figure 7. Analysis of variance of main influencing factors on fuel consumption.
Figure 7. Analysis of variance of main influencing factors on fuel consumption.
Aerospace 09 00104 g007
Figure 8. Flow diagram of the genetic algorithm.
Figure 8. Flow diagram of the genetic algorithm.
Aerospace 09 00104 g008
Figure 9. Comparison of optimization results and run times.
Figure 9. Comparison of optimization results and run times.
Aerospace 09 00104 g009
Table 1. Design parameters of system.
Table 1. Design parameters of system.
ParametersGroundCruise
Supply air temperature of the cabin266 K266 K
Supply air pressure101 kPa80 kPa
Flow rate of cabin supply air0.2 kg/s0.2 kg/s
Heat of electronic device50 kW50 kW
Heat of secondary energy system40 kW40 kW
Supply fuel temperature of engine343 K343 K
Flow rate of supply fuel of engine0.3 kg/s0.2 kg/s
Ram-air total temperature288 K267 K
Flow rate of low-temperature liquid-cooling cycle1.25 kg/s0.8 kg/s
Flow rate of high-temperature liquid-cooling cycle2.22 kg/s1.67 kg/s
Table 2. Results of parameter matching.
Table 2. Results of parameter matching.
Ground StateCruise State
MATLABEESDifferenceMATLABEESDifference
p 11 350.00350.000.000%350.00350.000.000%
p 12 348.30348.250.014%348.30348.250.014%
p 13 661.10661.62−0.078%636.00636.53−0.083%
p 14 657.70658.31−0.093%632.90633.35−0.071%
p 15 654.50655.02−0.079%629.70630.18−0.076%
p 16 101.50101.51−0.007%80.4080.40−0.002%
p 17 101.00101.000.000%80.0080.000.000%
T 11 600.00600.000.000%600.00600.000.000%
T 12 350.40350.400.000%423.00423.000.000%
T 13 444.40444.42−0.005%529.00529.06−0.012%
T 14 350.60350.570.009%387.50387.53−0.008%
T 15 337.00337.02−0.007%353.10353.070.009%
T 16 232.60232.550.021%235.20235.22−0.009%
T 17 266.00266.000.000%266.00266.000.000%
T 31 276.50276.54−0.015%274.70274.76−0.023%
T 32 288.30288.300.001%286.50286.52−0.007%
T 33 286.70286.72−0.006%285.00285.07−0.023%
T 41 328.30328.33−0.010%331.20331.25−0.014%
T 42 320.10320.070.010%317.10317.12−0.005%
T 43 328.00327.990.002%330.10330.090.002%
T 51 276.70276.670.011%265.00265.01−0.005%
T 52 318.00318.000.000%318.00318.000.000%
T 54 343.00343.000.000%343.00343.000.000%
Pressure unit, kPa; temperature unit, K.
Table 3. Ranges of thermal management system optimization variables.
Table 3. Ranges of thermal management system optimization variables.
VariablesNameMinMax
η H X 1 r 1 Efficiency of ACS’s primary HX0.450.9
η H X 1 r 2 Efficiency of ACS’s secondary HX0.450.9
η H X 13 Efficiency of air/low-temperature liquid-cooling HX0.450.85
η H X 14 Efficiency of air/high-temperature liquid-cooling HX0.450.85
η H X 45 Efficiency of fuel/high-temperature liquid-cooling HX0.450.85
η H X 5 r Efficiency of fuel/ram-air HX0.450.85
η c Efficiency of ACS’s compressor0.60.8
η t Efficiency of ACS’s turbine0.60.8
q m 3 Flow rate of low-temperature liquid-cooling cycle0.62.5
q m 4 Flow rate of high-temperature liquid-cooling cycle0.62.5
q m 5 Flow rate of fuel cycle0.62.5
Table 4. Comparison of two optimization designs’ results of the integrated thermal management system.
Table 4. Comparison of two optimization designs’ results of the integrated thermal management system.
Variables11 Variables6 Variables
SearchGAParticleSearchGAParticle
η H X 1 r 1 0.61370.53290.53720.54980.53700.5372
q m 5 0.62150.60890.60090.60680.60570.6013
η H X 1 r 2 0.79920.72400.72870.71220.72860.7294
q m 3 1.88730.72520.70180.71910.70230.6998
q m 4 0.67690.60010.60000.60000.60000.6000
η H X 5 r 0.72540.62030.61700.62120.61730.6170
η H X 45 0.55420.49740.49540.5 (fixed)0.5 (fixed)0.5 (fixed)
η H X 14 0.45000.45040.45000.45 (fixed)0.45 (fixed)0.45 (fixed)
η t 0.80000.80000.80000.8 (fixed)0.8 (fixed)0.8 (fixed)
η c 0.75050.80000.80000.8 (fixed)0.8 (fixed)0.8 (fixed)
η H X 13 0.59970.45050.45000.45 (fixed)0.45 (fixed)0.45 (fixed)
m f 123.4105107.7315105.9917106.0776108.0401105.9921
Search, pattern search algorithm; GA, genetic algorithm; Particle, particle swarm algorithm.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, R.; Dong, S.; Jiang, H.; Li, P.; Zhang, H. Parameter-Matching Algorithm and Optimization of Integrated Thermal Management System of Aircraft. Aerospace 2022, 9, 104. https://doi.org/10.3390/aerospace9020104

AMA Style

Wang R, Dong S, Jiang H, Li P, Zhang H. Parameter-Matching Algorithm and Optimization of Integrated Thermal Management System of Aircraft. Aerospace. 2022; 9(2):104. https://doi.org/10.3390/aerospace9020104

Chicago/Turabian Style

Wang, Ri, Sujun Dong, Hongsheng Jiang, Peiru Li, and Hainan Zhang. 2022. "Parameter-Matching Algorithm and Optimization of Integrated Thermal Management System of Aircraft" Aerospace 9, no. 2: 104. https://doi.org/10.3390/aerospace9020104

APA Style

Wang, R., Dong, S., Jiang, H., Li, P., & Zhang, H. (2022). Parameter-Matching Algorithm and Optimization of Integrated Thermal Management System of Aircraft. Aerospace, 9(2), 104. https://doi.org/10.3390/aerospace9020104

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