Online Model-Based Remaining-Useful-Life Prognostics for Aircraft Cooling Units Using Time-Warping Degradation Clustering

: Remaining-useful-life prognostics for aircraft components are central for efﬁcient and robust aircraft maintenance. In this paper, we propose an end-to-end approach to obtain online, model-based remaining-useful-life prognostics by learning from clusters of components with similar degradation trends. Time-series degradation measurements are ﬁrst clustered using dynamic time-warping. For each cluster, a degradation model and a corresponding failure threshold are proposed. These cluster-speciﬁc degradation models, together with a particle ﬁltering algorithm, are further used to obtain online remaining-useful-life prognostics. As a case study, we consider the operational data of several cooling units originating from a ﬂeet of aircraft. The cooling units are clustered based on their degradation trends and remaining-useful-life prognostics are obtained in an online manner. In general, this approach provides support for intelligent aircraft maintenance where the analysis of cluster-speciﬁc component degradation models is integrated into the predictive maintenance process.


Introduction
Globally, aircraft maintenance, repair and overhaul (MRO) costs account for 9% of the total airline operational costs [1]. To reduce these maintenance costs and, in particular, to reduce the costs of the maintenance needed in the case of an unexpected failure, MROs benefit from predictive maintenance. Here, the analysis of sensor health-monitoring data for aircraft components is central. By analyzing this data, the aim is to make predictions about the remaining-useful-life (RUL) of the components.
In most such studies, an inherent assumption is made that there is only one degradation mechanism or trend that governs the failure of components [14]. However, given that systems are often complex, with multiple interconnected components, and given that each component is subject to various operational and environmental factors, it is expected that the degradation of components follows more than one failure mechanism or trend [14,15]. As such, several degradation models should be integrated into the RUL prediction methodology, in order to capture the various degradation trends toward component failure.
Only a limited amount of studies developed RUL predictions using multiple degradation models. Data-driven approaches with multiple degradation models were considered in [16,17]. In [16], a Long Short-Term Neural Network was developed such that degradation patterns under multiple operational conditions and belonging to several fault modes were detected. With this approach, the RUL of aircraft engines was estimated.
In [17], a deep learning approach was combined with a genetic algorithm to predict the RUL of aircraft engines operating under multiple conditions and failing under several fault modes. Machine-learning approaches, however, are often seen by aircraft maintenance practitioners as "black-boxes" [6]. This is an obstacle for the implementation of purely datadriven prognostic methods in practice [18]. In this paper, we, therefore, focus on a modelbased prognostic approach that considers multiple models to capture the degradation trends of components.
Model-based RUL prognostic approaches that incorporate multiple degradation models and fault modes were considered in [14,[19][20][21][22][23]. In [19], four physics-based failure models of subsea pipelines were incorporated into one dynamic Bayesian network. In [14], the authors considered multiple resistance degradations for batteries, each with its own degradation trend. With these degradation models and a particle filtering algorithm, the authors predicted the RUL of batteries. A switching Kalman filter to model multiple degradation trends within one degradation process was used in [22,23] for bearings and in [20] for aircraft engines. In contrast to these studies, where multiple degradation models are integrated in one comprehensive degradation model, we first cluster the degradation trends of components and propose cluster-specific degradation models. In this way, different components are associated with different degradation models.
For model-based RUL prognostic methods, two frequently considered models are the exponential and the linear degradation models [2,24]. Linear degradation models were developed in [25][26][27][28]. In [25], a linear model with Brownian drift and random shocks, together with a particle filtering algorithm, was used to estimate the RUL of milling machines. A linear degradation model was also used in [26] for batteries, in [27] for aircraft engines, and in [28] for engine bleed valves.
Exponential degradation models were used in [29][30][31][32][33]]. An exponential model was used together with a particle filtering algorithm to predict the RUL of bearings in [29], and to predict the RUL of batteries in [30]. An exponential model for the degradation of bearings was also considered in [31], for batteries in [32], and for railway turnout systems in [33]. Similarly, in this paper, we consider clusters of components that degrade according to a linear or an exponential model.
In this paper, we propose an end-to-end approach to obtain online, model-based RUL prognostics for aircraft components by exploiting the knowledge obtained from clusters of component degradation trends. First, using sensor monitoring data, we construct a health indicator, which is used to diagnose components as being healthy or unhealthy. As soon as a component is diagnosed as unhealthy, a cluster-specific degradation model is selected for this component based on a dynamic time-warping clustering of a library of health indicators. These degradation models, together with a particle filtering algorithm, are further used to obtain RUL prognostics.
We illustrate our approach for the case of several aircraft Cooling Units, originating from a fleet of aircraft. We consider operational aircraft data, i.e., the sensor measurements were collected during actual operation of the aircraft. The results show that our proposed cluster-based RUL prognostics approach was robust and consistent across multiple components.
The main contributions of this paper are as follows: (i) We propose an end-to-end methodology from sensor measurements to diagnostics of a component as healthy or unhealthy to RUL prognostics. Using a health indicator, a component is diagnosed as healthy or unhealthy. Once a component reaches the unhealthy stage, a degradation model is selected based on the similarity between the degradation trend of this component and clusters of a library of health indicators. This approach exploits the potential learning from degradation trends of other components.
(ii) We consider multiple degradation models. A generic clustering method of the degradation of components is proposed. Each obtained cluster is associated with one degradation model and a corresponding failure threshold. (iii) Our proposed approach is illustrated using operational aircraft data, i.e., the sensor measurements were collected during the actual operation of the aircraft. In general, case studies are often conducted with simulated data or sensor data generated in a laboratory, which may not be as noisy as the aircraft operational data.
The remainder of this paper is organized as follows. In Section 2, we provide a methodology for online, model-based RUL prognostics using dynamic time-warping and a particle filtering algorithm. For time-warping, we consider several, same-type components, which are clustered based on the similarity between their degradation trends. In Section 3, we illustrate our methodology for several aircraft cooling units. Section 4 provides our conclusions and recommendations for further research.

Online, Model-Based Rul Prognostics for Multiple Components
In this section, we provide an approach to obtain online RUL prognostics for aircraft components. We first construct a health indicator from the sensor monitoring data. Based on this health indicator, we online diagnose a component as being healthy or unhealthy (step 1). Once a component is diagnosed as being unhealthy, a degradation model and corresponding failure threshold are selected for this component using dynamic timewarping data clustering of a library of series of health indicators (steps 2 and 3). Having selected a degradation model, RUL prognostics are obtained using a particle filtering algorithm (step 4). This process is illustrated in Figure 1.  Definition 1 (Component failure). For a component that degrades according to {X t , t ≥ 0}, we say that this component has failed if X t ≥ D, t > 0. Definition 2 (Component RUL). Given a current time f cur , the RUL of a component is defined as follows: RUL = inf{t : X f cur +t ≥ D|X 0 , . . . , X f cur }.

Step 1: Constructing a Health Indicator and Defining the Health Stage of the Component
Let I i f be a health indicator of component i after flight f . We assume that I i f is a function of the sensor measurements y f , y f −1 , . . . , y f −N , i.e., Here, a measurement y f depends on the degradation, y f = g(X i f , ψ f ) with X i f as the actual degradation of component i after flight f and ψ f , the noise.
We diagnose the component as being healthy or unhealthy by analyzing whether I i f exceeds an alarm threshold T alarm defined by the Chebyshev's inequality [34][35][36], which specifies that, for any probability distribution with a specified mean µ and standard deviation σ, at most 1 k 2 percent of the values from this distribution fall outside the µ ± kσ interval, k > 0. This implies that where µ is the mean and σ is the standard deviation of the health indicator while it is healthy. We, thus, use the following alarm threshold T alarm : As soon as I i f > T alarm , a component is diagnosed as unhealthy (see Figure 1). Let f i,alarm denote the first flight when a component i is diagnosed as unhealthy. Let f i, f ail denote the flight when component i fails.

Step 2: Selecting a Degradation Model for a Component
As soon as a component is diagnosed as being unhealthy, we select a degradation model and a corresponding failure threshold D based on a similarity analysis of clusters in a library of health indicators. As a last step, we use this selected degradation model together with a particle filtering algorithm to obtain RUL prognostics.
We consider an offline library of n health indicators. Each health indicator is constructed starting m flights before the component is diagnosed as unhealthy (see Step 1) until the component fails.
Let the series of indicators We are interested in clustering the set of time series I i , i ∈ {1, 2, . . . , n} into clusters where the underlying deterioration of the health indicators follows a similar path toward failure (similar degradation paths).
We propose (i) a Dynamic Time Warping (DTW) algorithm [37,38] to determine the minimum distance between any two condition series . . , f j, f ail ]} of component i and component j, respectively, and, (ii) we cluster the n series of health indicators based on this minimum distance metric.

(i) Dynamic Time-Warping for an Offline Library with Time-Series of Health Indicators
We define the following Euclidean distance between any two points I i f and I j g of the health indicators time series Here, we see d(I i f , I j g ) as the distance or the dissimilarity between points (health indicator values) I i f and I j g . The larger the distance, the more dissimilar the health indicator values are.
We next find a |I i | × |I j | warping path of minimum distance between the two health indicator series I i and I j , where |I i | is the number of flights in the health indicator timeseries of component i, and |I j | the number of flights in the health indicator time-series of component j. We define a |I i | × |I j | warping path between I i and I j as follows: Condition (a) in Definition 3 implies that the first and last elements of health indicator series I i and I j are warped (aligned), respectively. Condition (b) ensures that the time moments for both series are monotonically increasing. Last, condition (c) ensures that no element of {I i } and {I j } can be omitted or repeated in the warping (alignment), i.e., the elements in the sequence need to be contiguous. Figure 2 shows an example of a warping path for two time-series of health indicators, Let p be a warping path between degradation series I i and I j . Let d(I i f , I j g ), as introduced in Equation (4), be a similarity metric between the degradation points I i f ∈ I i and I j g ∈ I j . Let D(p, I i , I j ) denote a similarity metric between two series of health indicators I i and I j with warping path p, where We are interested in a warping path P between the two series I i and I j with a minimum cumulative Euclidean distance W(I i , I j ) = D(p, I i , I j ), i.e, we are interested in obtaining: Lastly, we obtain a warping path of minimum cumulative distance W(I i , I j ) using dynamic programming as follows [38]. Let C( f , g ) = W(Ĩ i ,Ĩ j ) be the minimum distance between a part of the time series I i and a part of the time series I j .
We can now solve the recursive equation to find the minimum distance between two health indicator time-series as:

(ii) Clustering Series of Health Indicators of a Library of n Components
Having determined the distance W(I i , I j ) between the components i, j ∈ {1, 2, . . . , n}, i = j, we construct a graph G(V, E, δ), where a node v i ∈ V corresponds to a series of health indicators for component i, i ∈ {1, 2, . . . , n}, and thus |V| = n. We construct an edge between a node v i to the closest δ nodes v j , i.e., the shortest δ distances between a component i and the rest of the components j ∈ {1, 2, . . . , n}, i = j. Here, we normalize the distance with the maximum number of health indicator values in I i and I j . The connected components of this graph are considered to be the clusters of health indicator series (or equivalently, clusters of components with similar underlying degradation trends).
We assume that all components in the same cluster have a degradation trend according to a stochastic process that has the same functional form, and that all components in the same cluster have the same failure threshold D. These cluster-specific degradation models and failure thresholds are ultimately used to estimate the RUL for a component using a particle filtering algorithm.

Step 3: Online Clustering of (Non-Failed) Components
Let component i be monitored online. Let an offline library consist of n series of health indicators {I j , j ∈ {1, 2, . . . , n}, j = i. Let C be the set of clusters in the offline libraries, which are obtained using dynamic time-warping (see Section 2.2.2).
As soon as the health indicator of component i exceeds an alarm threshold T alarm (i.e., the component is diagnosed as unhealthy), the online monitored component i is assigned to a cluster in the set C using dynamic time-warping (Section 2.  Figure 3 illustratesĨ i corresponding to an online monitored component i. The partial health indicator seriesĨ j for component j from the offline library is defined as: where |Ĩ i | is the length of seriesĨ i . Figure 3 shows an example for an online monitored component i and a library with component j ∈ {1, 2}. Lastly, we assign component i to a clusterc ∈ C such that the average minimum distance W (Ĩ i ,Ĩ j ) to the series of health indicators of components j in clusterc is minimised, i.e.,: Now, component i assumes the same degradation model and failure threshold as the degradation model and a failure threshold specific to clusterc. This degradation model is further used to estimate the RUL for component i.

Step 4: RUL Prognostics
Together with the degradation model and failure threshold D identified in Step 3, we use a particle filtering algorithm [39] to estimate the RUL of an online monitored component. Let us consider the following recurrence function: where X f −1 is the actual degradation level of a component after flight f − 1 and ω f are i.i.d. Gaussian variables representing the noise in the degradation process. In general, X f is assumed to have a certain functional form, which depends on the type of component and/or the type of degradation that the component is expected to undergo. In our case, we consider cluster-specific functional forms for all components in the same cluster (see Steps 2 and 3). The health indicator I i f is a function of X f (see Equation (1)): To obtain RUL prognostics, the particle filtering algorithm has four steps [39]: (i) prediction, (ii) updating, (iii) resampling, and (iv) prognostics.

(i) Prediction
We are interested in the conditional probabilities: where Equation (12) is the conditional probability density function of the degradation level of the component after flight f − 1, given the recorded degradation measurements I f −1 , . . . , I 1 , whereas Equation (11) is the transition probability density function to reach future degradation state X f , given the current degradation state X f −1 .
(ii) Updating As soon as new measurements are available, the state probability density function is updated, using Bayes' theorem, as: . . , I 1 ) .

(iii) Resampling
We approximate Equation (13), the probability density function for the state degradation, numerically using Importance Sampling [40] as follows. First, we sample M particles from a probability density functionp X f |X f −1 ,I f ,...,I 1 (X f |X f −1 , I f , . . . , I 1 ). Then, using importance sampling, the probability density function of the degradation state of the component after flight f is approximated as: where δ is a Dirac function, X j f is the estimated degradation level of particle j of the component after flight f , and w j f is the weight of the jth particle, j ∈ {1, 2, . . . , M} after flight f , which is updated and normalized as follows: In every re-sampling cycle, the particles are re-sampled proportionally to their weight.
A new particle set {x j f }, j ∈ {1, 2, . . . , M} is generated by re-sampling M new particles, where the probability to resample a particle j equalsŵ j f . (iv) The RUL prognostic Lastly, we consider a threshold D > 0 and define the RUL as the stopping time:

Numerical Case Study-Multiple Cooling Units from a Fleet Of Aircraft
To illustrate our approach from Section 2, we considered eight cooling units (CUs) originating from a fleet of aircraft operated by a large European airline. For each CU, a time series of measurements was recorded during the operation of these aircraft until the components failed (run-to-failure data).
The CU is a vapor cycle refrigeration unit consisting of a compressor, a condenser, an evaporator and a flash tank (see Figure 4). After many days of operations, the CU becomes clogged with burned oil, moisture, and sludge from the compressor. This accelerates the compressor wear. Long-term exposure to these conditions negatively affects the health of the CU, which, in time, leads to failure [41].

Health Indicator For CUs
For each CU, run-to failure measurements were considered. These measurements were available after some initial usage of the component. The run-to-failure data series consists of sequences of 50 flights (short series) up to 390 flights (long series). For each CU, there were nine sensors S = {1, 2, . . . , 9}, each generating a measurement every 10 s during each flight. For the purpose of our analysis, the data sets are anonymized, and thus it is unknown the type of measurement each sensor generates. We observed, in the sensor data, an increasing trend in the mean and maximum sensor measurements toward failure. As an example, Figure 5 shows the mean and maximum sensor measurements per day until the moment of failure for one CU and for each of the nine available sensors. Let y i,s f ,b denote the b th measurement during flight f for CU i generated by sensor s. We first normalize the sensor measurements as follows: with the min s and max s the available minimum and maximum measurements generated by sensor s, respectively, and B i f as the total number of measurements recorded during flight f of CU i.
With these normalized sensor measurements, we considered several health indicators, as described in [2,42]. We constructed each health-indicator for the measurements of each sensor and calculated, for each sensor, the correlation coefficient (trendability) between the health indicator and the time to failure for the last 50 flights before failure (length of the shortest serie). Table 1 provides an overview of the considered health-indicators and the corresponding correlation coefficients. An extensive description of the health indicators can be found in [42].
The highest correlation coefficients were obtained for the Root Mean Square (RMS) health indicator. We, thus, constructed a health indicator based on the RMS of the measurements [10,42,43]. A health indicator based on the RMS of measurements was often employed in the literature [2,42], for example for the health monitoring of gearboxes [44], turbine-cutting tools [43], and rolling element bearings [45]. The RMS of sensor s of CU i during flight f is: A health indicator I i f for CU i during flight f was obtained as the moving average of the maximum RMS obtained by the sensors s ∈ S , as follows: where N = 10 and where S ⊆ S is the set of sensors for which the moving average of the RMS of the measurements obtained in the last 50 flights before failure has a correlation coefficient with the time of failure of at least 0.70. We, therefore, included sensors 4, 5, 7, and 8 in the health-indicator. This final health indicator emphasizes the increase in the mean and the maximum measurements toward failure, which we observe in Figure 5. Table 1. Overview of the considered health indicators. For each sensor, the correlation coefficient between the health indicator and the time to failure is given. The highest absolute correlation coefficient is given in bold.
Health-Indicator Sensors 1 2 3 4 5 6 7 8 9  Figure 6 shows the health indicator and the alarm threshold T alarm when considering the run-to-failure data of the eight CUs. For the alarm threshold, we assume k = 2 (see Equation (3)), obtaining a threshold of T alarm = 11.54. When analyzing all CUs, the alarm threshold is reached between 2 to 40 flights before the actual failure time. Figure 6 also shows that the degradation trends toward failure, which differed across components. The failure thresholds differed among the eight components. As such, in the next section, we cluster the series of health indicators to identify the main degradation models and the corresponding failure thresholds.

Clustering the Health Indicators and Determining Clustering-Specific Degradation Models
We performed dynamic time-warping (see Section 2.2, Step 2) to cluster the n = 8 time series of health indicators and to identify cluster-specific degradation models for the CUs. For this, we used the run-to-failure data of the eight CUs, starting m = 10 flights before the CU reached the alarm threshold until failure time. Figure 7 shows the minimum distance W(I i , I j ), i, j ∈ {1, 2, . . . , 8}, between the series of health indicators of the eight CUs, normalized with the maximum length of the series I i and I j ). Here, the minimum distances W(I i , I j ) between the degradation of any two components varies from 0.44 (high similarity between the degradation trends) to 4.24 (large dissimilarities between the degradation trends).
Using W(I i , I j ), i, j ∈ {1, 2, . . . , 8}, between the eight series of health indicators, we constructed the graph G(V, E), with |V| = 8. Here, a node v i ∈ V in the graph corresponds to the series of health indicators for component i, i ∈ {1, 2, . . . , 8}. For each CU i, and a given δ, an edge was constructed between those δ nodes (CUs) that have the smallest distance to node v i (CU i). For example, if δ = 2 and we consider CU i = 4, an edge was constructed between node v 4 and node v 6 and between node v 4 and v 8 , since the distances between node v 4 and nodes v 6 and v 8 are the smallest (see Figure 7).   In Figure 8, we set δ = 1. This results in two clusters of CUs: {4, 6, 8}, and CU {1, 2, 3, 5, 7}. The health indicator in Figure 6 shows that the CUs 4, 6, and 8 failed soon after the health indicator reached the alarm threshold-when the health indicator reached a value of approximately 15. The degradation trend around failure showed a monotonic increase. In contrast, CUs 1,2,3,4,5, and 7 failed when the health indicator reached a value of 20. Furthermore, the degradation of these CUs often accelerated toward failure.
For the two clusters in Equations (21) and (22), we now define the following clusterspecific degradation models.

Cluster 1-Linear Degradation Model
The health indicators of CUs 4, 6, and 8 in Cluster 1 are given in Figure 6d,f,h. All these CUs fail when their health indicator reached a value of approximately 15. Thus, we consider a failure threshold D = 15 for Cluster 1. Once the health indicator of the CUs in Cluster 1 reaches the alarm threshold T alarm , the indicator had a sharply, but monotonic increasing, trend toward failure. We, thus, consider the following linear degradation model for the CUs in Cluster 1 [8]: where X i t is the degradation of component i at time t, α i is the initial degradation, β i is a model parameter, and ω i is the noise in the degradation process. As in [30], we ignore the process noise because it is handled through the uncertainty in the model parameters and in the measurements [46]. Linear degradation models are often considered for prognostics [24], for example for milling machines [25], batteries [26], aircraft engines [27], and engine bleed valves [28].
Following the degradation model in Equation (23), we consider the derivative: Rewriting the above equation in the form of Equations (10) and (1) with dt = 1 flight, we obtain: where  2). The parameter distributions are initialized such that sample impoverishment, the degeneracy problem and non-convergence of the particle filtering algorithm are avoided [46].

Cluster 2-Exponential Degradation Model
The health indicators of CUs 1,2,3,5, and 7 in Cluster 2 are given in Figure 6. All these CUs failed when their health indicator reached a value of approximately 20. Thus, we consider a failure threshold D = 20 for Cluster 2.
Moreover, the increments of the health indicator steadily increased toward failure for these CUs. We, thus, consider the following exponential degradation model for the CUs in Cluster 2: where X i t is the degradation of component i at time t, α i is the initial degradation, β i is a model parameter, and ω i is the noise in the degradation process. As before, we ignore the process noise [30]. Exponential degradation models have been considered in many prognostic studies [24], for bearings [29,31], batteries [30,32], and railway turnout systems [33].
Taking the logarithm of Equation (27): As in [8], we assume that α i = 0. Following the degradation model in Equation (28), we consider the derivative: Rewriting the above equation in the form of Equations (10) and (1) with dt = 1 flight, we obtain: where  2). The parameter distributions are initialized such that sample impoverishment, the degeneracy problem, and the non-convergence of the particle filtering algorithm are avoided [46].
Lastly, with the degradation models introduced above, we apply the particle filtering algorithm (see Section 2, Step 4) to estimate the RUL of the CUs.

Rul Estimation Results
In this section, we estimate the RUL of the eight CUs using leave-one-out cross validation as follows. We select one CU i ∈ {1, 2, . . . , 8}. We consider the (partial) health indicator of this CU i up to the moment of generating a RUL prediction. We determine the minimum distance between this (partial) health indicator and the (partial) health indicators of the CUs in Cluster 1 and Cluster 2 (see Equations (21) and (22)), which do not include CU i. We assign CU i to that cluster for which the average distance between the cluster and CU i is the minimum (see Step 3, Section 2). Then, CU i assumes the degradation model and the failure threshold of this cluster. Lastly, we estimate the RUL of CU i using particle filtering with 10,000 particles. Table 2 shows for each CU the cluster it is assigned to as soon as T alarm is reached, i.e., the CU is diagnosed as unhealthy, and there is a corresponding RUL prediction. CU 1 was diagnosed as unhealthy, and an alarm was triggered at an early stage at 40 flights before the actual failure time. CUs 2, 3, 4, 5, 6, and 7 were diagnosed as unhealthy up to six flights before the actual failure. For CU 8, an alarm was triggered only two flights before the actual failure.
CUs 1, 2, 5, and 7 were assigned to Cluster 2, which assumed an exponential degradation model, whereas CUs 3, 4, 6, and 8 were assigned to Cluster 1, which assumed a linear degradation model. Only CU 3 was assigned to a different cluster (namely cluster 1) than in the offline clustering, which led to underestimation of the RUL. These results show that, even when considering clustering with partial health indicators, i.e., up to the moment when a prediction is made, our approach can well identify the underlying linear or exponential degradation trend and the corresponding failure threshold. Table 2 shows the RUL estimates 10 flights before the time of failure. CUs 4, 6, and 8 failed within less than 10 flights as soon as they were diagnosed as unhealthy, thus, no RUL predictions are made in this case. CUs 1, 2, 3, 5, and 7 were all assigned to the same exponential cluster as in the offline clustering. Table 2 also shows that the RUL of the CUs was well estimated as soon as these CUs were diagnosed as unhealthy. The Root Mean Square Error (RMSE) was low as well. The fact that we were able to select appropriate degradation models and failure thresholds based on cluster-specific trends contributes to the accuracy of the RUL predictions.

Conclusions
We proposed an online, model-based RUL estimation approach for aircraft components using a clustering of component degradation trends. Using clustering, we determined cluster-specific degradation models and failure thresholds. Together with a particle filtering algorithm, these degradation models were used to estimate the RUL of multiple, same-type aircraft components.
As an illustration, we considered several aircraft Cooling Units for which sensor measurements were recorded during the operation of the aircraft. The results show that our proposed methodology was able to identify the degradation models of components and estimate their RUL. From a practical point of view, our RUL estimation results have the potential to support aircraft maintenance stakeholders with maintenance task scheduling.
For future work, we plan to extend our methodology for other aircraft components. At the same time, we plan to propose maintenance task scheduling algorithms to evaluate the impact of prognostic-driven maintenance on the maintenance costs and aircraft availability.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:

CU
Cooling unit RMS Root mean square RUL Remaining-useful-life