Next Article in Journal
Numerical Simulation and Experimental Study on the Aerodynamics of Propulsive Wing for a Novel Electric Vertical Take-Off and Landing Aircraft
Previous Article in Journal
Gaze Movements of Helicopter Pilots during Real and Simulated Take-Off and Landing Maneuvers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anti-Icing System Performance Prediction Using POD and PSO-BP Neural Networks

1
Environmental Control and Oxygen System Department, COMAC Shanghai Aircraft Design and Research Institute, Shanghai 201210, China
2
Laboratory of Fundamental Science on Ergonomics and Environmental Control, School of Aeronautic Science and Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Aerospace 2024, 11(6), 430; https://doi.org/10.3390/aerospace11060430
Submission received: 12 March 2024 / Revised: 22 May 2024 / Accepted: 23 May 2024 / Published: 26 May 2024

Abstract

:
The anti-icing system is important for ice protection and flight safety. Rapid prediction of the anti-icing system’s performance is critical to reducing the design time and increasing efficiency. The paper proposes a method to quickly predict the anti-icing performance of the hot air anti-icing system. The method is based on Proper Orthogonal Decomposition (POD) and Back Propagation (BP) neural networks improved with the Particle Swarm Optimization (PSO) algorithm to construct the PSO-BP neural network. POD is utilized for data compression and feature extraction for the skin temperature and runback water obtained by numerical calculation. A lower-dimensional approximation is derived from the projection subspace, which consists of a set of basis modes. The PSO-BP neural network establishes the mapping relationship between the flight condition parameters (including flight height, atmospheric temperature, flight speed, median volume diameter, and liquid water content) and the characteristic coefficients. The results show that the average absolute errors of prediction with the PSO-BP neural network model on skin temperature and runback water thickness are 3.87 K and 0.93 μm, respectively. The method can provide an effective tool for iteratively optimizing hot air anti-icing system design.

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.

2. Methodology

2.1. Proper Orthogonal Decomposition

The Proper Orthogonal Decomposition method allows a series of basis modes to be obtained based on the observed complex physical field. For a sample of the target information field U = { U 1 , , U i , , U M } , Ui represents a set of data for the information field. The sample U is a snapshot matrix of order N × M, where N represents the dimension of the data sample Ui, i.e., the number of the grid nodes in Computational Fluid Dynamics (CFD). The essence of the POD algorithm is to find a set of optimal orthogonal basis ψ such that the information field of a sample can be represented by its linear representation as follows [32]:
U i = j = 1 M α j i ψ j ,
where α j i is the orthogonal basis coefficient, and ψ j is the orthogonal basis. All ψ j form an optimal orthogonal basis matrix, and the optimal orthogonal basis satisfies that the average value of the projections of all samples on the optimal basis is maximized. Since the optimal orthogonal basis matrix ψ is of the same order as the sample matrix U, it can be expressed as follows:
ψ = U V ,
where the matrix V is the matrix consisting of the eigenvectors of the matrix C, and it satisfies the following:
C V = Λ V .
Λ is the eigenvalue of the matrix C that can be computed from the sample matrix as follows:
C i j = 1 M k = 1 N U k i U k j .
Then the sample matrix can be expressed as a linear superposition of the basis modes. The eigenvalue reflects the proportion of energy in the corresponding basis mode in the sample [31]. Sorting the eigenvectors of matrix C in descending order of the corresponding eigenvalues, the energy of the flow field in the direction of the jth optimal orthogonal basis as a proportion of the total energy is as follows:
E n / E t o t a l = j = 1 n λ j / j = 1 N λ j ,
where λ j stands for the corresponding eigenvalue of the jth basis mode. It can realize the sample dimensionality reduction by ignoring some basis modes with smaller energy and selecting the reserved basis modes according to the energy proportion.

2.2. PSO-BP Neural Network

2.2.1. BP Neural Network

The BP neural network is a multilayer forward neural network trained with the error backpropagation algorithm [35], which consists of an input layer, hidden layers, and an output layer connected sequentially. A typical BP neural network is shown in Figure 1. x1, x2, …, xm are the inputs of the neural network; υ h i denotes the connection weight between the input layer and the hidden layer; a1, a2, …, ag, …, aq are the inputs of the neurons in the hidden layer; b1, b2, …, bg, …, bq are the outputs of the neurons in the hidden layer; w k h is the connection weight between the hidden layer and the output layer; β1, …, βn are the inputs of the neurons in the output layer; and y1, …, yn are the outputs of the neural network.
The training process of a BP neural network includes forward propagation of the work signal and back propagation of errors. The main steps are described as follows:
Step 1: Initially, the connection weights and thresholds of the BP neural network are randomly initialized in the range (0, 1).
Step 2: Input the training datasets X and compute the input ah and the output bh of each neuron in the hidden layer, as well as the input β k and output y k of each neuron in the output layer.
a h = i = 1 m v h i x i b h = f ( a h γ h ) ,         h = 1 , 2 , , g , , q
β k = h = 1 q w k h b h y k = f ( β k θ k ) ,         k = 1 , 2 , , n
where γ h is the threshold of the hth neuron in the hidden layer; q represents the number of neurons in the hidden layer; θk is the threshold of the kth neuron in the output layer; n stands for the number of neurons in the output layer; and f is the activation function.
Step 3: Calculate the mean square error e of a neural network. If the mean square error e is less than or equal to the target computational accuracy of the neural network ε, the neural network training will be completed; conversely, training needs to continue.
e = 1 n k = 1 n ( y k y ^ k ) 2
Step 4: The error signal is back propagated. The connection weights and thresholds of the output layer are adjusted according to the output error signal e.
e w k h = e y k y k β k β k w k h = y k 1 y k y ^ k y k b h
e δ k = e y k y k δ k = y k 1 y k y ^ k y k
w k h = w k h + Δ w k h = w k h η e w k h
θ k = θ k + Δ θ k = θ k η e θ k
Δ θ k shows the variation of the connection weights between the kth neuron in the output layer and the hth neuron in the hidden layer.
Step 5: Similar to adjusting the connection weights and thresholds of the output layer, the error signal is back-propagated. The connection weights and thresholds of the hidden layer are adjusted based on the output error signal e.
e v h i = e y k y k β k β k b h b h a h a h v h i = y k 1 y k y ^ k y k b h 1 b h x i k = 1 n w k h
e γ h = e y k y k β k β k b h b h γ h = y k 1 y k y ^ k y k b h 1 b h k = 1 n w k h
v h i = v h i η e v h i
γ h = γ h η e γ h
Step 6: Finally, update the connection weights and thresholds of the neural network and repeat steps 2 to 5 until the accuracy requirement is satisfied.

2.2.2. PSO Optimization Algorithm

The PSO algorithm was first introduced by Eberhart and Kennedy and used for the optimization of continuous nonlinear functions [36]. The PSO algorithm is a population-based random optimization technique that simulates the process of bird and fish feeding in nature and searches for the global optimal solution to the problem through the collaboration of the population. The particles in the PSO algorithm have only two attributes, including velocity and position. Each particle searches for the optimal solution independently in the search space, which is recorded as the current individual optimal position ppbest. The value of the optimization objective function is set to fp. The particle shares its individual extreme value with the other particles in the whole particle swarm. The optimal individual extreme value is taken as the current global optimal position of the whole particle swarm pgbest. The value of the optimization objective function is denoted as fg. All particles in the particle swarm dynamically adjust their velocities and positions according to the current individual extreme values and the current global optimal solution shared by the whole swarm [37]. The particles update their positions and velocities according to the following equation:
v i d ( k + 1 ) = w v i d ( k ) + c 1 r 1 p pbest , i d ( k ) x i d ( k ) + c 2 r 2 p gbest , d ( k ) x i d ( k )   and
x i d ( k + 1 ) = x i d ( k ) + v i d ( k + 1 ) ,
where i means the particle serial number; d is the particle dimension number; k is the number of iterations; w is the inertia weight; c1 is the individual learning factor, and c2 is the population learning factor; and r1 and r2 are the random numbers in the range (0, 1), which are used to increase the randomness of the search. vid(k) denotes the velocity vector of the dth dimension for the particle i during the kth iteration process. xid(k) presents the position vector of the dth dimension for the particle i during the kth iteration process. Moreover, ppbest,id is the individual optimal solution of the dth dimension for the particle i during the kth iteration process, and pgbest,d is the population optimal solution of the dth dimension for the population during the kth iteration process.

2.2.3. PSO-BP Neural Network Construction

The biggest challenges for designing and training neural networks are the selection and determination of various hyperparameters, including the number of hidden layers, the number of neurons in each layer, the activation function, the training algorithm, the number of iterations, and the methods to prevent overfitting. In order to obtain the optimal hyperparameters of the PSO-BP neural network model conveniently and realize the adaptive adjustment and optimization of the network structure, the neuron number in the hidden layer, the iterations, the dropout regularization rate, and the training batch size of the BP neural network are optimized with the PSO algorithm in this paper. The construction process of the PSO-BP neural network is shown in Figure 2. The PSO-BP neural network is constructed based on the Google artificial intelligence system TensorFlow (version 2.10.0) in Python (version 3.9.16).

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:
t Ω W a d Ω + Ω ( F a , c F a , v ) d S = 0 .
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 m ˙ w can then be calculated as follows:
β = α u d n   and
m ˙ w = L W C U β ,
where u d 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 U 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).
Q ˙ V + Q ˙ F = Q ˙ imp + Q ˙ rad + Q ˙ evap + Q ˙ ice + Q ˙ conv + Q ˙ cond
This indicates that energy change in the control body is caused by the heat flow of the impinging water Q ˙ imp the radiation heat flow Q ˙ rad , the water evaporative heat flow Q ˙ evap , the heat flow of the frozen water Q ˙ ice , the heat flow of runback water Q ˙ conv , and the heat flow of solid wall thermal conductivity Q ˙ cond . Equation (22) is expressed in local differential form as follows:
ρ w h f c p , w T ˜ t + d i v u ¯ f h f c p , w T ˜ = c p , w T ˜ d , + u ¯ d 2 2 × U L W C β + σ ε T 4 T 4 L evap + L sub + c p , w T ˜ m ˙ evap + L fus c p , cee T ˜ m ˙ ice + Q ˙ conv + Q ˙ cond
where Q ˙ cond is obtained by thermal conductivity calculation; ρ , c f , c s , σ , ε , L evap , L fus are the physical parameters of the fluid and solid wall; T and U are the temperature and velocity of the far-field air flow. m ˙ evap is the evaporated water mass flow rate, which is attained from Q ˙ conv . m ˙ ice 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 V i = [ x i 1 , x i 2 , , x i m ] ,   i = 1 , 2 , , n 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 T i = [ T i 1 , T i 2 , , T i p ] and F i = [ F i 1 , F i 2 , , F i p ] , 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 ψ T = ψ T 1 , ψ T 2 , , ψ T q and ψ F = ψ F 1 , ψ F 2 , , ψ F q , which can reflect the characteristics of the skin temperature and runback water thickness distribution, and q-dimensional skin temperature samples T POD i = [ T POD , i 1 , T POD , i 2 , , T POD , i q ] and runback water thickness samples F POD i = [ F POD , i 1 , F POD , i 2 , , F POD , i q ] of characteristic coefficients with a number of n.
Step 4: The PSO-BP neural network is utilized to separately establish the mapping relationship for V i and T POD i , as well as V i and F POD i , 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 V Predict can be fed into each PSO-BP neural network model. Then, we can achieve the characteristic coefficients of the skin temperature distribution T POD Predict and the runback water thickness F POD Predict . 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 MAET and MAEF as follows:
MAE T = i = 1 p T i T ^ l p MAE F = i = 1 p f i f ^ l p
where Ti and fi are the skin temperature and runback water thickness obtained by numerical simulation at the ith grid point. T ^ l and f ^ l 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:
Error T = j = 1 n MAE T , j n Error F = j = 1 n MAE F , j n ,
where ErrorT and ErrorF are the fitting error or prediction error of the skin temperature and runback water thickness for all samples, respectively. MAET,j and MAEF,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.

4. Results and Discussion

4.1. POD Fitting Results

The basis modes, eigenvalues, and occupied energy of the skin temperature distribution and the runback water thickness distribution are obtained by downscaling the sample data using the aforementioned POD method. The results are displayed in Table 4 and Table 5. The variance contribution of a basis mode reflects the degree of the basis mode containing the sample characteristics. The cumulative variance contribution of the first 10 basis modes in the skin temperature datasets reaches 99.26%, while the cumulative variance contribution of the first 10 basis modes in the runback water thickness datasets is calculated to be 99.08%. The MAE of the skin temperature (MAET,POD) and runback water thickness MAEF,POD with the POD fitting method are illustrated in Figure 5 and Figure 6. As the characteristics of the skin temperature distribution on the anti-icing surface are more easily captured, the average absolute error ErrorT,POD and average relative error between the simulation datasets and the skin temperature data after the POD and inverse POD processes are 1.65 K and 0.36% when selecting the first 10 basis modes, respectively. The distribution of runback water on the anti-icing surface is more complex than the temperature distribution. In the case of higher surface temperatures, the distribution area of runback water is very small or non-existent, making it more difficult to capture these characteristics. The average absolute error ErrorF,POD and the average relative error between the simulation datasets and the runback water thickness data after the POD and inverse POD processes are 0.42 μm and 6.28%, respectively, when selecting the first 10 basis modes.
The first 10 orders of basis modes are used to fit the distribution of skin temperature and runback water thickness. The comparison of the POD fitting results with the CFD simulation results and the prediction errors for two of the training samples are shown in Figure 7 and Figure 8. It can be seen that the POD fitting results and CFD simulation results are in good agreement, and the POD fitting results can accurately extract the characteristics of the skin temperature distribution and the thickness distribution of runback water. Figure 7b also demonstrates the power of the POD method to predict the runback water distribution when there is no water film on the skin surface.

4.2. PSO-BP Neural Network Prediction Results

According to the analysis in Section 4.1, the anti-icing performance of the hot air anti-icing system can be accurately predicted using the first 10 orders of POD basis modes. Therefore, PSO-BP neural network models with a number of 10 are established to rapidly predict the skin temperature and runback water thickness, i.e., setting up V i and X POD , i 1 , X POD , i 2 ,…, X POD , i 1 , respectively. After training the PSO-BP neural network models, the characteristic coefficients of the skin temperature distribution T POD Predict and the runback water thickness F POD Predict of 100 test samples can be predicted by the corresponding POS-BP models. The details of part of the 100 test samples for skin temperature prediction and runback water thickness prediction are displayed in Table 6 and Table 7, respectively. The reverse Proper Orthogonal Decomposition (POD) method is used to predict the skin temperature and runback water thickness on 2580 grid nodes of the skin surface, enabling rapid prediction of wing anti-icing performance.
Figure 9 and Figure 10 demonstrate the mean absolute error MAE of skin temperature (MAET,PSO-BP) and runback water thickness (MAEF,PSO-BP) obtained by comparing the PSO-BP neural network model and CFD simulation results for 100 test samples. The average absolute error ErrorT,PSO-BP of the PSO-BP neural network model in predicting skin temperature is 3.87 K for all test samples. The average absolute error ErrorF,PSO-BP of the PSO-BP neural network model to predict runback water thickness is 0.93 μm for all test samples. The prediction results for the anti-icing performance are of the same order of magnitude as those of the POD fitting results, indicating that the established fast prediction PSO-BP model is able to effectively predict the skin temperature and runback water thickness distribution characteristics on the skin surface of the hot air anti-icing system.
Figure 11 and Figure 12 display the prediction results of the PSO-BP neural network model and the CFD simulation results for two of the training samples, along with a comparison between these two sets of results. By synthetically analyzing the prediction effect of the fast prediction model, the established PSO-BP neural network can accurately predict the skin temperature and runback water thickness of the hot air anti-icing system. The PSO-BP neural network model can judge the distribution of runback water on the anti-icing surface. When the skin temperature is low and the thickness of the runback water is large, the model can accurately capture the distribution characteristics of the runback water. When there is only a small amount of runback water on the anti-icing surface or no runback water, the distribution characteristics are more difficult to capture. Furthermore, there is a possibility of a discrepancy in the prediction of the PSO-BP neural network model, which leads to a large prediction error.
The model training and prediction process is conducted on a computer equipped with an Intel(R) Core (TM) i7-9700 3.00 GHz CPU (Intell®, Santa Clara, CA, USA). The results indicate that the time required for one model to train is around 10 min. The cost time of skin temperature prediction or runback water thickness prediction with the PSO-BP neural network and the process for output as a text file for a sample is 2.1255 s on average. Comparing with the average 3.5 h required for the CFD numerical simulation, the PSO-BP neural network significantly improves prediction efficiency and effectively resolves the issues of high cost and long duration associated with the CFD numerical simulation of the hot air anti-icing system. This is achieved by reducing calculation complexity and the consumption of computing resources.

5. Conclusions

In this paper, a fast prediction method for the performance of a hot air anti-icing system based on the POD method and the PSO-BP neural network is established. The flight condition parameters, including flight height, atmospheric temperature, flight speed, MVD, and LWC, are used as input parameters for the fast prediction model. A dataset with a number of 1434 obtained by FENSP-ICE is applied to train and test the POS-BP neural network model. The prediction effect for rapidly predicting the skin temperature and runback water thickness distribution on the anti-icing surface is evaluated for a hot air anti-icing system. The conclusions are as follows:
  • The POD method can accurately extract the skin temperature and runback water distribution characteristics of the anti-icing surface. The average absolute error ErrorT,POD and ErrorF,POD for the skin temperature and runback water thickness are 1.65 K and 0.42 μm when selecting the first 10 basis modes, respectively. It means that downscaling the sample data to obtain the basis modes, eigenvalues, and their occupied energies makes the fast prediction strategy feasible.
  • The prediction results with the PSO-BP neural network show that the established model can more accurately predict the skin temperature and runback water thickness distribution characteristics of the hot air anti-icing system. The average absolute error ErrorT,PSO-BP and ErrorF,PSO-BP of the PSO-BP neural network model to predict skin temperature and runback water thickness are 3.87 K and 0.93 μm for all 100 test samples. The distribution of runback water thickness becomes more challenging to capture when there is only a small amount of runback water or none at all on the anti-icing surface. This leads to a more difficult prediction of the runback water thickness distribution.
  • For cases requiring extensive performance analysis of a hot air anti-icing system, the model based on the POD method and the PSO-BP neural network can realize fast and accurate predictions of the performance of the hot air anti-icing system. This model holds high engineering application value and development potential due to its computational efficiency and accuracy.
  • In the future, different neural networks, such as the convolutional neural network (CNN), will be implemented to construct the neural network framework for improved training results. Meanwhile, the wing anti-icing experiment data obtained from flight tests or wind tunnel tests will be fed into the training models as a more convincing training and testing dataset with high accuracy.

Author Contributions

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

Funding

This research received no external funding.

Data Availability Statement

The dataset is available upon request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ANNartificial neural network
BPback propagation
CFDComputational Fluid Dynamics
CNNconvolutional neural network
IoTInternet of Things
LWCliquid water content
MAEmean absolute error
MLmachine learning
MVDmedian volume diameter
PODProper Orthogonal Decomposition
PSOParticle Swarm Optimization
WPTwavelet packet transform
XGBoostExtreme Gradient Boosting

References

  1. Lynch, F.T.; Khodadoust, A. Effects of ice accretions on aircraft aerodynamics. Prog. Aerosp. Sci. 2001, 37, 669–767. [Google Scholar] [CrossRef]
  2. Hamid, H.A.; Yong, W.K.; Yusoff, H.; Ismail, M. CFD Impingement Flow Study on Temperature Profile of Concave Plate. Adv. Res. Fluid Mech. Therm. Sci. 2022, 95, 1–16. [Google Scholar] [CrossRef]
  3. Fukudome, K.; Tomita, Y.; Uranai, S.; Mamori, H.; Yamamoto, M. Evaluation of Anti-Icing Performance for an NACA0012 Airfoil with an Asymmetric Heating Surface. Aerospace 2021, 8, 294. [Google Scholar] [CrossRef]
  4. Green, S. A study of US inflight icing accidents and incidents, 1978 to 2002. In Proceedings of the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 9 January 2006. [Google Scholar]
  5. Papadakis, M.; Wong, S.H.; Yeong, H.W.; Wong, S.C.; Vu, G.T. Experimental investigation of a bleed air ice protection system. SAE Trans. 2007, 116, 643–663. [Google Scholar]
  6. Brown, J.M.; Raghunathan, S.; Watterson, J.K.; Linton, A.J.; Riordon, D. Heat transfer correlation for anti-icing systems. J. Aircr. 2002, 39, 65–70. [Google Scholar] [CrossRef]
  7. Myers, T.G. Extension to the Messinger model for aircraft icing. AIAA J. 2001, 39, 211–218. [Google Scholar] [CrossRef]
  8. Lee, J.; Rigby, D.; Wright, W.; Choo, Y. Analysis of Thermal Ice Protection System (TIPS) with Piccolo Tube Using State-of-the-Art Software. In Proceedings of the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 9 January 2006. [Google Scholar]
  9. Wang, H.; Tran, P.; Habashi, W.G.; Chen, Y.C.; Zhang, M.; Feng, L.J. Anti-icing simulation in wet air of a piccolo system using FENSAP-ICE. SAE Trans. 2007, 116, 715–723. [Google Scholar]
  10. Thomas, S.K.; Cassoni, R.P.; MacArthur, C.D. Aircraft anti-icing and de-icing techniques and modeling. J. Aircr. 1996, 33, 841–854. [Google Scholar] [CrossRef]
  11. Hinton, G.E.; Osindero, S.; Teh, Y.W. A fast learning algorithm for deep belief nets. Neural Comput. 2006, 18, 1527–1554. [Google Scholar] [CrossRef]
  12. Johnstone, C.; Sulungu, E.D. Application of neural network in prediction of temperature: A review. Neural Comput. Appl. 2021, 33, 11487–11498. [Google Scholar] [CrossRef]
  13. Smith, B.A.; McClendon, R.W.; Hoogenboom, G. An enhanced artificial neural network for air temperature prediction. Int. J. Comput. Inf. Eng. 2007, 1, 2166–2171. [Google Scholar]
  14. Patil, K.; Deo, M.C. Prediction of daily sea surface temperature using efficient neural networks. Ocean Dyn. 2017, 67, 357–368. [Google Scholar] [CrossRef]
  15. Sarkar, P.P.; Janardhan, P.; Roy, P. Prediction of sea surface temperatures using deep learning neural networks. SN Appl. Sci. 2020, 2, 1458. [Google Scholar] [CrossRef]
  16. Wei, L.; Guan, L.; Qu, L.Q. Prediction of sea surface temperature in the South China Sea by artificial neural networks. IEEE Geosci. Remote Sens. Lett. 2019, 17, 558–562. [Google Scholar] [CrossRef]
  17. Piotrowski, A.P.; Napiorkowski, M.J.; Napiorkowski, J.J.; Osuch, M. Comparing various artificial neural network types for water temperature prediction in rivers. J. Hydrol. 2015, 529, 302–315. [Google Scholar] [CrossRef]
  18. Ma, Y.L.; Geng, Y.F.; Chen, X.H.; Lu, Y.X. Prediction for asphalt pavement water film thickness based on artificial neural network. J. Southeast Univ. Engl. Ed 2017, 33, 490–495. [Google Scholar]
  19. Wang, B.; Chen, B.; Wang, G.; Li, R.; Wen, J.; Lu, C.; Tian, R.; Deng, J. Back propagation (BP) neural network prediction and chaotic characteristics analysis of free falling liquid film fluctuation on corrugated plate wall. Ann. Nucl. Energy 2020, 148, 107711. [Google Scholar] [CrossRef]
  20. Ogretim, E.; Huebsch, W.; Shinn, A. Aircraft ice accretion prediction based on neural networks. J. Aircr. 2006, 43, 233–240. [Google Scholar] [CrossRef]
  21. Chang, S.N.; Leng, M.Y.; Wu, H.W.; Thompson, J. Aircraft ice accretion prediction using neural network and wavelet packet transform. Aircr. Eng. Aerosp. Technol. Int. J. 2016, 88, 128–136. [Google Scholar] [CrossRef]
  22. Yi, X.; Wang, Q.; Chai, C.C.; Guo, L. Prediction Model of Aircraft Icing Based on Deep Neural Network. Trans. Nanjing Univ. Aeronaut. Astronaut. 2021, 38, 535–544. [Google Scholar]
  23. Strijhak, S.; Ryazanov, D.; Koshelev, K.; Ivanov, A. Neural network prediction for ice shapes on airfoils using icefoam simulations. Aerospace 2022, 9, 96. [Google Scholar] [CrossRef]
  24. Dong, Y.Q. An application of deep neural networks to the in-flight parameter identification for detection and characterization of aircraft icing. Aerosp. Sci. Technol. 2018, 77, 34–49. [Google Scholar] [CrossRef]
  25. Caliskan, F.; Aykan, R.; Hajiyev, C. Aircraft icing detection, identification, and reconfigurable control based on Kalman filtering and neural networks. Aerosp. Eng. 2008, 21, 51–60. [Google Scholar] [CrossRef]
  26. Aykan, R.; Caliskan, F.; Hajiyev, C. F16 icing identification based on neural networks. IFAC Proc. Vol. 2004, 37, 201–206. [Google Scholar] [CrossRef]
  27. Li, S.B.; Qin, J.K.; He, M.; Paoli, R. Fast evaluation of aircraft icing severity using machine learning based on XGBoost. Aerospace 2020, 7, 36. [Google Scholar] [CrossRef]
  28. Cheng, X.; Shi, F.; Zhao, M.; Li, G.Y.; Zhang, H.X.; Chen, S.Y. Temporal attention convolutional neural network for estimation of icing probability on wind turbine blades. IEEE Trans. Ind. Electron. 2021, 69, 6371–6380. [Google Scholar] [CrossRef]
  29. Abdelghany, E.S.; Farghaly, M.B.; Almalki, M.M.; Sarhan, H.H.; Essa, M.E.S.M. Machine learning and iot trends for intelligent prediction of aircraft wing anti-icing system temperature. Aerospace 2023, 10, 676. [Google Scholar] [CrossRef]
  30. Yang, Q.; Guo, X.F.; Li, Q.; Dong, W. Hot air anti-icing performance estimation method based on POD and surrogate model. Acta Aeronaut. Astronaut. Sin. 2023, 44, 626992. (In Chinese) [Google Scholar]
  31. Ostrowski, Z.; Białecki, R.A.; Kassab, A.J. Estimation of constant thermal conductivity by use of proper orthogonal decomposition. Comput. Mech. 2005, 37, 52–59. [Google Scholar] [CrossRef]
  32. Cazemier, W.; Verstappen, R.; Veldman, A.E.P. Proper orthogonal decomposition and low-dimensional models for driven cavity flows. Phys. Fluids 1998, 10, 1685–1699. [Google Scholar] [CrossRef]
  33. Lappo, V.; Habashi, W.G. POD-Kriging Approximations of Multi-Disciplinary CFD Simulation with Application to In-Flight Icing. In Proceedings of the 17th Annual Conference of the CFD Society of Canada, Ottawa, ON, Canada, 3–5 May 2009. [Google Scholar]
  34. Jung, S.K.; Shin, S.; Myong, R.S.; Cho, T.H. An efficient CFD-based method for aircraft icing simulation using a reduced order model. Mech. Sci. Technol. 2011, 25, 703–711. [Google Scholar] [CrossRef]
  35. Rumelhart, D.E.; Hinton, G.E.; Williams, R.J. Learning representations by back-propagating errors. Nature 1986, 323, 533–536. [Google Scholar] [CrossRef]
  36. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95-International Conference on Neural Networks, Perth, WA, Australia, 27 November–1 December 1995; IEEE: Piscataway, NJ, USA, 1995; Volume 4, pp. 1942–1948. [Google Scholar]
  37. Xu, Y.S.; Zeng, B.; Yin, X.W. BP Neural Network and its applications based on improved PSO. Comput. Eng. Appl. 2009, 45, 233–235. [Google Scholar]
  38. Beaugendre, H.; Morency, F.; Habashi, W.G. FENSAP-ICE’s three-dimensional in-flight ice accretion module: ICE3D. J. Aircr. 2003, 40, 239–247. [Google Scholar] [CrossRef]
  39. Bourgault, Y.; Habashi, W.G.; Dompierre, J.; Baruzzi, G.S. A finite element method study of Eulerian droplets Impingement models. Int. J. Numer. Methods Fluids 1999, 29, 429–449. [Google Scholar] [CrossRef]
Figure 1. Typical BP neural network structure.
Figure 1. Typical BP neural network structure.
Aerospace 11 00430 g001
Figure 2. The construction process of the PSO-BP neural network.
Figure 2. The construction process of the PSO-BP neural network.
Aerospace 11 00430 g002
Figure 3. The model of the fast prediction for a wing’s hot air anti-icing performance.
Figure 3. The model of the fast prediction for a wing’s hot air anti-icing performance.
Aerospace 11 00430 g003
Figure 4. Temperature difference between the numerical simulations and experiments.
Figure 4. Temperature difference between the numerical simulations and experiments.
Aerospace 11 00430 g004
Figure 5. MAE of skin temperature with the POD fitting method.
Figure 5. MAE of skin temperature with the POD fitting method.
Aerospace 11 00430 g005
Figure 6. MAE of runback water thickness with the POD fitting method.
Figure 6. MAE of runback water thickness with the POD fitting method.
Aerospace 11 00430 g006
Figure 7. Comparison of the POD fitting results with the CFD simulation results (case 1): (a) skin temperature; (b) runback water thickness (no water film).
Figure 7. Comparison of the POD fitting results with the CFD simulation results (case 1): (a) skin temperature; (b) runback water thickness (no water film).
Aerospace 11 00430 g007
Figure 8. Comparison of the POD fitting results with the CFD simulation results (case 2): (a) skin temperature; (b) runback water thickness.
Figure 8. Comparison of the POD fitting results with the CFD simulation results (case 2): (a) skin temperature; (b) runback water thickness.
Aerospace 11 00430 g008
Figure 9. MAE of skin temperature with the PSO-BP neural network.
Figure 9. MAE of skin temperature with the PSO-BP neural network.
Aerospace 11 00430 g009
Figure 10. MAE of runback water thickness with the PSO-BP neural network.
Figure 10. MAE of runback water thickness with the PSO-BP neural network.
Aerospace 11 00430 g010
Figure 11. Comparison of the PSO-BP prediction results with the CFD simulation results (case 1): (a) skin temperature; (b) runback water thickness.
Figure 11. Comparison of the PSO-BP prediction results with the CFD simulation results (case 1): (a) skin temperature; (b) runback water thickness.
Aerospace 11 00430 g011aAerospace 11 00430 g011b
Figure 12. Comparison of the PSO-BP prediction results with the CFD simulation results (case 2): (a) skin temperature; (b) runback water thickness.
Figure 12. Comparison of the PSO-BP prediction results with the CFD simulation results (case 2): (a) skin temperature; (b) runback water thickness.
Aerospace 11 00430 g012
Table 1. Cases for anti-icing simulation model calibration.
Table 1. Cases for anti-icing simulation model calibration.
CaseFlight Height (ft)T (°C)March NumberMVD (μm)LWC (g/m3)
110,473−6.90.4413.90.288
210,773−7.50.4215.60.375
311,000−8.10.4116.80.400
413,115−10.80.4515.00.298
Table 2. Flight icing conditions for some cases.
Table 2. Flight icing conditions for some cases.
CaseFlight Height (ft)T (°C)March NumberMVD (μm)LWC (g/m3)
1000.18150.80
1000−150.25300.17
50015,000−300.60200.14
143431,000−400.82500.05
The case sequence number increases from 1 to 1434 in increments of 1. The symbol of “⋮” means the conditions for some cases are omitted and not shown here.
Table 3. The optimal values of hyperparameters for the PSO-BP neural network.
Table 3. The optimal values of hyperparameters for the PSO-BP neural network.
HyperparameterEigenvalues λj
Number of the hidden layers 4
Number of neurons in the input layer5
Number of neurons in the first hidden layer30
Number of neurons in the second hidden layer64
Number of neurons in the third hidden layer64
Number of neurons in the fourth hidden layer30
Number of neurons in the output layer10
Activation functionReLU
Loss functionMSE
Optimization algorithmAdam Optimizer
Learning rate0.0001
Table 4. Eigenvalues and occupied energy for POD downscaling of skin temperature.
Table 4. Eigenvalues and occupied energy for POD downscaling of skin temperature.
OrderEigenvalues λjOccupied Energy Ej
1702,557.54950.7182
2165,618.12130.1693
333,996.05860.0348
424,281.47960.0248
518,064.11570.0185
611,198.06210.0114
75171.19980.0053
84163.14180.0042
93474.41300.0036
102385.20140.0024
Table 5. Eigenvalues and occupied energy for POD downscaling of runback water thickness.
Table 5. Eigenvalues and occupied energy for POD downscaling of runback water thickness.
OrderEigenvalues λjOccupied Energy Ej
1157,450.94760.8940
210,544.69520.0599
32160.22160.0123
41523.90890.0086
51054.44960.0060
6492.18540.0028
7437.27160.0025
8347.11640.0020
9251.62070.0014
10216.22440.0012
Table 6. Test samples of a fast prediction model for skin temperature distribution.
Table 6. Test samples of a fast prediction model for skin temperature distribution.
CaseFlight Height (ft)T (°C)March NumberMVD (μm)LWC (g/m3)
10−7.50.35400.42
104000−50.35251.55
4019,000−300.50450.16
10032,000−300.60350.05
The case sequence number increases from 1 to 100 in increments of 1. The symbol of “⋮” means the conditions for some cases are omitted and not shown here.
Table 7. Test samples of a fast prediction model for runback water thickness.
Table 7. Test samples of a fast prediction model for runback water thickness.
CaseFlight Height (ft)T (°C)March NumberMVD (μm)LWC (g/m3)
1000.18400.10
3014,500−100.45250.30
4019,000−300.50450.16
10032,000-300.82350.05
The case sequence number increases from 1 to 100 in increments of 1. The symbol of “⋮” means the conditions for some cases are omitted and not shown here.
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

Mao, H.; Lin, X.; Li, Z.; Shen, X.; Zhao, W. Anti-Icing System Performance Prediction Using POD and PSO-BP Neural Networks. Aerospace 2024, 11, 430. https://doi.org/10.3390/aerospace11060430

AMA Style

Mao H, Lin X, Li Z, Shen X, Zhao W. Anti-Icing System Performance Prediction Using POD and PSO-BP Neural Networks. Aerospace. 2024; 11(6):430. https://doi.org/10.3390/aerospace11060430

Chicago/Turabian Style

Mao, Handong, Xiaodan Lin, Zhimao Li, Xiaobin Shen, and Wenzhao Zhao. 2024. "Anti-Icing System Performance Prediction Using POD and PSO-BP Neural Networks" Aerospace 11, no. 6: 430. https://doi.org/10.3390/aerospace11060430

APA Style

Mao, H., Lin, X., Li, Z., Shen, X., & Zhao, W. (2024). Anti-Icing System Performance Prediction Using POD and PSO-BP Neural Networks. Aerospace, 11(6), 430. https://doi.org/10.3390/aerospace11060430

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