1. Introduction
Aircraft icing occurs when supercooled droplets that may be present in clouds at air temperatures below freezing are caught on the windward surfaces of components during flight, such as the wings and engine inlet. Ice accumulation has great damage to the aerodynamic characteristics of aircraft, resulting in reduced lift by 30% and increased drag by 40% [
1,
2,
3], and even causing flight accidents [
4]. Therefore, anti-icing systems or de-icing systems are essential for protecting aircraft from ice accumulation for a safe flight. The hot air anti-icing system is widely used in the aircraft industry [
5,
6], where hot bleed airflow induced by the engine is discharged from the holes of the piccolo tube installed at the wing leading edge to heat the skin surface. However, it will result in increasing engine fuel consumption, reducing engine thrust, and wasting large amounts of energy. It is necessary to carry out iterative optimization of the hot air anti-icing system design to enhance the design efficiency and minimize energy consumption.
Advanced numerical simulation of thermal anti-icing systems has been highly anticipated as a supplementary design and certification tool, in addition to anti-icing experiments and flight tests [
7,
8,
9]. Modeling and numerical simulation for hot air anti-icing system design is a complex problem involving multiphase flow, conjugate heat, and mass transfer with phase change [
10]. The calculation of the flow field and temperature field with high accuracy requires a significant amount of computational time, which increases as the grid number increases. Aircraft icing conditions vary vastly for anti-icing system design. The factors mentioned above result in a significant increase in computing time for anti-icing system performance. Therefore, it is essential to establish an effective method to obtain skin temperature and the runback water distribution of other icing conditions from the existing computational or experimental data quickly.
Recently, artificial neural networks (ANNs) have greatly advanced research in extracting data features with strong adaptivity, learning ability, and fault tolerance [
11]. It can accurately approximate the complex nonlinear function mapping relationships between data and learn the essential characteristics of data from a limited sample set. ANN has been widely used in air temperature prediction as part of weather prediction [
12,
13], sea surface temperature prediction [
14,
15,
16], and water temperature prediction [
17]. Meanwhile, the prediction of asphalt pavement water film thickness with ANN [
18] and the thickness of liquid films on corrugated plate walls with the back propagation (BP) neural network [
19] also demonstrate the power of neural networks. As to aircraft wing anti-icing and de-icing applications, ANN have been widely used for fast prediction of ice shape [
20,
21,
22,
23], in-flight parameters for icing detection [
24,
25,
26], aircraft icing severity [
27], icing probability [
28], skin temperature [
29], and runback water flow [
30]. Ogretim E. [
20] proposed a method of ice shape fast prediction using a neural network with icing conditions as input and ice shape parameters as output. The experimental data on ice accretion in the NACA 0012 wing model were used as the training sample sets for the neural network. Chang et al. [
21] presented a new technique that combines wavelet packet transform (WPT) and ANN to predict ice accretion on the surface of an airfoil. Results showed an advantage of WPT in performing the analysis of ice accretion information, and the prediction accuracy was improved as well. Strijhak et al. [
23] discussed the procedure and method for the ice accretion prediction for different airfoils using ANNs, which were based on the results of the numerical experiments and performed well. Yiqun Dong [
24] applied a deep neural network to identify and characterize aircraft icing for in-flight parameter detection. In [
27], the authors introduced a purely data-driven approach to finding the complex pattern between different flight conditions and aircraft icing severity prediction using machine learning based on the Extreme Gradient Boosting (XGBoost) algorithm. Abdelghany et al. [
29] presented a novel approach based on machine learning (ML) and the Internet of Things (IoT) to predict the thermal performance characteristics of a partial span wing anti-icing system constructed using the NACA 23014 airfoil section. A high-precision computation of skin temperature and runback water thickness relies on a large grid number. Therefore, using anti-icing performance datasets with higher dimensionality as the training samples of the neural network will result in significantly increased computing time.
Proper Orthogonal Decomposition (POD) is a powerful method for order reduction and data compression. It offers an efficient way to capture the dominant features of a system with multiple degrees of freedom and represent the desired precision using a relevant set of modes, thereby reducing the order of the system [
31,
32]. Skin temperature, runback water thickness, and other anti-icing performance parameters can be decomposed into the basis modes that express data characteristics with POD. Then, the linear fitting coefficients needed to reconstruct the anti-icing performance are obtained using the basis modes [
30]. As a result, the high-dimensional datasets with anti-icing performance can be downgraded to low-dimensional samples of fitting coefficients, leading to a significant reduction in computing time and data storage space. Habashi et al. [
33] established a fast prediction model of aircraft 3D icing ice shape based on POD, and the results showed that the POD method can improve prediction accuracy by 600~800 times. SungKi Jung et al. [
34] used POD for dimensionality reduction and integrated it with a neural network to accurately predict the collection efficiency and ice accretion shapes on an airfoil.
In this paper, a fast prediction model for the performance of a hot air anti-icing system based on POD and PSO-BP neural networks that combine BP neural networks with Particle Swarm Optimization (PSO) is proposed. The anti-icing performance, including skin temperature and runback water thickness, is obtained through numerical simulation using FENSAP-ICE as the original datasets. The high-dimensional icing performance data are then order-reduced by POD to attain the basis modes and characteristic coefficients. Finally, the PSO-BP neural network is used to establish the mapping relationship between the flight condition parameters, including flight height, atmospheric temperature, flight speed, median volume diameter (MVD), and liquid water content (LWC), and the characteristic coefficient above, which realizes a fast prediction of the hot air anti-icing system performance under various flight icing conditions.
3. Prediction of Anti-Icing Performance with the PSO-BP Neural Network
3.1. Dataset Preparation
3.1.1. Model and Cases for Anti-Icing Simulation
The parameters affecting the simulation results of the hot air anti-icing system can be categorized into external icing conditions and internal hot air anti-icing parameters. Flight tests have shown that aircraft icing conditions depend on several factors, including meteorological parameters like atmospheric temperature, cloud extent, MVD, and LWC; the flight status of the aircraft, such as flight height, flight speed, angle of attack, etc.; and the factors that determine the heat transfer characteristics of the flow in the hot air anti-icing cavity are mainly the flow, pressure, temperature of the bleed air jetted from the piccolo tube, and so on. The flight condition parameters of flight height, atmospheric temperature, and flight speed, represented by the Mach number later, MVD, and LWC, are taken as the input parameters of the fast prediction model, as shown in
Figure 3. The angle of attack and the flow, pressure, and temperature of the bleed air jetted from the piccolo tube are the parameters directly related to flight speed and flight height according to the actual situation, which can be determined by interpolation calculation. The output parameters of the PSO-BP neural network include skin temperature and runback water distribution.
In the study, a part of the wing skin for a hot air anti-icing system is selected as the target research object, which is shown in
Figure 3. The length of the wing chord and spanwise are 1.08 m and 0.08 m, and the thickness of the wing skin is 1.8 mm. The meshing result of the target wing skin (solid domain) turns out to be 82,831 grid nodes. The reliability of the meshing result for CFD simulation has been verified by previous engineering projects. The target skin surface has 2580 grid nodes, i.e., the dimension
p of a single simulation sample is 2580. The model of the hot air anti-icing system for CFD simulation is calibrated with the test data. Some of the test cases are displayed in
Table 1. The numerical simulation is conducted under the same conditions as the test cases with FENSAP-ICE (version 19.2). FENSAP-ICE is used to simulate the anti-icing performance of a hot air anti-icing system, which is a commercial CFD software for icing and anti-icing calculation [
38]. The difference between the simulation results and test data for surface temperature along the wing chordwise direction is illustrated in
Figure 4. The values along the
x-axis direction in
Figure 4 indicate the distance between the target point and the stationary point, and the point where
x = 0 means the stationary point position. The position where
x > 0 is located on the upper surface of the airfoil, and the points with
x < 0 lie on the lower surface. The results show that the difference between the CFD simulation data and the test data for the upper surface temperature is within 10 K, and the comparison results for the lower surface temperature are within 15 K. The maximum of the average temperature for all four cases in the direction of wing span is counted to be 9.98 K, which meets the CFD calculation accuracy requirements for engineering applications.
Homogeneous sampling is conducted within the parameter space of flight conditions, which is determined by the intersection of the limiting icing envelope, the relationship between icing meteorological conditions, and the flight envelope. To obtain a more representative dataset, appropriate encryption is applied when sampling the limiting icing envelope and the flight envelope. Then, 1434 uniformly distributed samples are obtained for model training and testing, and details of some cases are shown in
Table 2.
3.1.2. Dataset Preparation with FENSAP-ICE
FENSAP-ICE is the second generation of icing and anti-icing analysis software, which applies modular thinking to separate and combine various steps of icing and anti-icing simulation to obtain different target results. Firstly, the external air flow field and internal flow field in the anti-icing cavity can be computed with the FENSAP-ICE module named FENSAP. The classical compressible Navier–Stokes equations are employed as the governing equation, which can be written in the following integral and conservative forms:
Then the droplet impingement properties are analyzed by the FENSAP-ICE module named DROP3D. The governing equations of the droplet impingement are based on the Eulerian model proposed by Bourgault [
39]. This is essentially a two-fluid model consisting of a set of Navier–Stokes equations augmented by droplet-related continuity and momentum equations. The local collection efficiency
β and the mass flow rate of impact water
can then be calculated as follows:
where
stands for the droplet velocity vector;
α denotes the volume fraction, i.e., the proportion of volume occupied by water droplets in the control volume; and
means the velocity of air flow in the far field. The shallow-water model is used in the FENSAP-ICE module named ICE3D to simulate surface water flow and heat transfer.
For the simulation calculation of a hot air anti-icing system, the FENSAP-ICE module named CHT3D is applied and adopts a loosely coupled method to exchange the data of the external flow field, the water film motion, the solid heat conduction, and the internal flow field. Then, the distribution of the skin temperature and runback water thickness after convergence can be obtained. The CHT3D module takes into account the energy balance relationship in the simulation process of anti-icing, as shown in Equation (22).
This indicates that energy change in the control body is caused by the heat flow of the impinging water
the radiation heat flow
, the water evaporative heat flow
, the heat flow of the frozen water
, the heat flow of runback water
, and the heat flow of solid wall thermal conductivity
. Equation (22) is expressed in local differential form as follows:
where
is obtained by thermal conductivity calculation;
are the physical parameters of the fluid and solid wall;
T∞ and
U∞ are the temperature and velocity of the far-field air flow.
is the evaporated water mass flow rate, which is attained from
.
is the icing mass flow rate.
hf is the water film thickness, and
T is the wall equilibrium temperature. The average time on CFD simulation for a case is about 3.5 h.
3.2. Fast Prediction Model of Hot Air Anti-Icing System Performance Based on POD and PSO-BP Neural Network
In this paper, a fast prediction model for the performance of a hot air anti-icing system is established based on the abovementioned POD method and PSO-BP neural network. The whole process is as follows:
Step 1: Determine the input parameters with a number of m and obtain m-dimensional vectors of flight condition parameters by data sampling homogeneously in the flight condition parameter space.
Step 2: Obtain the performance of the hot air anti-icing system by the simulation method in the target cases with the input parameters identified in the first step. If the grid number of the skin model is p, the skin temperature and runback water distribution samples of p dimensions with a number of n can be attained and recorded as and , respectively.
Step 3: Reduce the dimensionality of the anti-icing performance parameter samples with the POD method. Then, we can get the first q-order basis modes and , which can reflect the characteristics of the skin temperature and runback water thickness distribution, and q-dimensional skin temperature samples and runback water thickness samples of characteristic coefficients with a number of n.
Step 4: The PSO-BP neural network is utilized to separately establish the mapping relationship for and , as well as and , respectively. The PSO-BP neural network models are established for each dimension of the characteristic coefficients. The skin temperature prediction models and runback water thickness prediction models can both be built with a number of q.
Step 5: After the training of the PSO-BP neural network models is completed, the flight condition parameters to be predicted can be fed into each PSO-BP neural network model. Then, we can achieve the characteristic coefficients of the skin temperature distribution and the runback water thickness . Then, the inverse POD method can be used to obtain the skin temperature and runback water thickness, which allows for fast prediction of the anti-icing performance of the hot air anti-icing system.
In the study, the dimension of the skin temperature and runback water distribution samples
p is 2580, and the number of datasets
n is 1434. The number of input parameters
m is set to 5, which is the number of neurons in the input layer. The number of neurons in the output layer corresponds to the order of basis modes
q after the POD process, which turns out to be 10. The optimal values of hyperparameters in the multilayer neural network are shown in
Table 3.
3.3. Error Analysis
In order to effectively evaluate the fitting effect of POD and the prediction effect of the PSO-BP neural network model, the mean absolute error (MAE) is used as the evaluation indices. The smaller the MAE is, the more accurate the model is. In this paper, the skin temperature and runback water thickness are the output values that we focus on. Therefore, the MAE of skin temperature and runback water thickness for a single sample are denoted as MAE
T and MAE
F as follows:
where
Ti and
fi are the skin temperature and runback water thickness obtained by numerical simulation at the
ith grid point.
and
are the skin temperature and runback water thickness attained from the POD fitting results or PSO-BP neural network prediction at the
ith grid point, respectively.
p means the dimension of the single sample.
The error of the POD fitting model and the PSO-BP neural network prediction model is defined as the average error of all samples and calculated as follows:
where Error
T and Error
F are the fitting error or prediction error of the skin temperature and runback water thickness for all samples, respectively. MAE
T,j and MAE
F,j are the fitting error or prediction error of the skin temperature and runback water thickness for the
jth sample.
n means the number of all samples.