Abstract
To advance the traditional structural integrity and life estimation approach, this article presents an explanation of damage and develops an alternative approach to damage assessment based on the fundamental laws of thermodynamics. This explanation provides a new way to perform prognostics and health management in terms of the irreversible thermodynamic framework. Specifically, it is proposed and demonstrated that entropy generated due to material degradation can be used as an index of damage. It is claimed that the use of an entropy-based damage index offers a unified and more comprehensive measure of damage. This is because energy dissipation associated with failure mechanisms that lead to material degradation involves irreversible processes, resulting in entropy generation, which is physically measurable as the materials degrade. We therefore conclude that entropy generation can be used to assess the degree of damage, the amount of the life of materials expended, and the extent of the life remaining. As an application of the proposed concept, this article discusses the entropic characterization of the corrosion-fatigue degradation mechanism, in which a standard entropy-based damage prognostics and health management approach is employed in integrity assessment and remaining useful life prediction of aluminum 7075-T651 specimens. Using entropy as a thermodynamic state function for damage characterization is shown to be effective in handling the endurance threshold uncertainties for prediction purposes.
Keywords
Introduction
The definition of damage has traditionally relied on the field variables chosen to describe the anticipated aging or degradation process. For example, for fatigue, damage variables used to describe the expected damage include the ratio of instantaneous to final crack lengths, 1 the linear summation of normalized number of cycles accumulated at different stresses (Miner’s 2 rule), the accumulation of normalized surface layer stress, 3 the cumulative plastic strain ratio, 4 the linear summation of persistent slip band length ratio, 5 the accumulation of the endurance limit ratio,6,7 and the ratio of defect surface area to total cross-sectional area. 8 For wear damage, the local plastic deformation, 9 the wear volume, 10 and delamination of composites 11 have been used. Moreover, the definition of damage can vary at different geometric scales. For example, in the corrosion-fatigue degradation mechanism, transmission electron microscopy and dislocation modeling studies show that at the nano-scale, adsorbed hydrogen ions localize plastic deformation and lower the energy required for fracture, resulting in degradation. 12 At the micro-scale level, it has been shown that weakness of oxide film for extensive plastic deformation relative to the underlying metal may result in the localization of slip bands and a reduction in near-surface plasticity, ultimately leading to the formation of dislocation structures and voids. 13 At the meso-scale level, damage constitutes the growth and coalescence of micro-cracks that initiate a crack. At a macro-scale level, damage is defined as the growth and propagation of macro-cracks that result in a final fracture. 14
Difficulties in proposing a consistent definition of damage from the physical and mechanical points of view have compelled some researchers to look for a microscopically consistent definition in the context of the continuum damage mechanics (CDM).8,15 For example, in CDM, damage,
where
As the measurement of damage in conventional definitions is subjective and dependent on the choices of the observable variables that uniquely correlate with and describe the underlying physical processes leading to damage, entropy generation has been proposed by a number of researchers as a more inclusive and comprehensive measure to characterize damage. The underlying principle of this approach is that all damage processes and failure mechanisms share a common feature at a deeper level: energy dissipation. Energy dissipation is a fundamental measure for irreversibility, which in a thermodynamic treatment of non-equilibrium processes can be quantified by entropy generation. Energy dissipative processes can be directly related to thermodynamic entropy, for example, plasticity and dislocations,16,17 erosion corrosion, 18 wear and fracture, fretting corrosion, 19 thermal degradation, 20 and associated failure of tribological components.21,22
There are two distinct interpretations of entropy for damage characterization. In the statistical mechanics microscopic interpretation of entropy, introduced by Boltzmann,
23
the connection between damage and entropy may be established using equation (2), where the entropy,
where
Feinberg and Widom35,36 made early efforts to describe a thermodynamic-based definition of damage. They developed a thermodynamic characterization of degradation dynamics, which employs the Gibbs free energy as a measure of material or component parameter degradation. According to their model, disorder/damage must increase entropy or reduce thermodynamic energies. The use of Gibbs free energy and its reduction to represent the degradation, however, is only valid when the pressure and temperature of the system remain constant, 37 which constrained their model of degradation to isothermal–isobaric conditions (e.g. Gibbs energy-based damage modeling cannot capture degradations due to thermal losses). Furthermore, they assumed that change in the system characteristic results from a log-time aging behavior versus time. This is true only if the damage accumulates at a fixed rate. 34 In fact, the emerging field of damage mechanics tracks entropy and thermodynamic energies to quantify the behavior of irreversible degradations including tribological processes such as friction and wear38–42 and fretting and fatigue30,25,43–46 that incline to increase entropy generation of the components or structures. Following this trend, Bryant et al. 47 established a more general relationship between the damage of systems undergoing irreversible dissipative processes and the associated entropy generation called degradation-entropy generation (DEG) theorem. In their approach, Bryant et al. linked the traditional damage measure based on the observable markers of degradation (e.g. wear volume, crack size, and corrosion mass loss) to the associated entropy generation through a linear relation. Their actual results focused on tribological systems with a single dissipative process of wear. 47 They concluded that the famous Archard 10 law for measuring wear in fretting wear is the consequence of the DEG theorem. Clearly, connection of the entropy generation to markers of damage requires that the entropy generated exclusively be generated for the specific degradation process that caused the observable marker. While this is theoretically valid, measurement of the entropy that can be exclusively attributed to a marker of damage would be nearly impossible. In another research performed by Sosnovskiy and Sherbakov, 48 it is mathematically confirmed that entropy generation represents damage metric. They used entropy generation to evaluate damage evolution of materials subjected to fatigue under variable temperature.
Imanian and Modarres 49 offered an alternative degradation assessment approach to formally describe the resulting damage. Their approach uses a theorem relating entropy generation to the dissipation process via the generalized thermodynamic forces and thermodynamic fluxes. In this theorem, entropy generation due to all degradation processes, including those that act synergistically, describes a normalized damage index. Failure occurs when the total entropy generation reaches a level called entropic endurance. At this level, the corresponding damage index reaches unity whereby cumulative damage reaches a level beyond which the specimen or component reaches a zone of concern (but not necessarily the point of failure that might be interpreted as an observable damage). Examples of the entropic endurance level include a point where it is judged that a small margin to failure might remain. The following attributes describe key properties of the entropic endurance:
Entropic endurance describes the capacity of a unit to withstand entropy.
Entropic endurance of the same units is equal.
Entropic endurance of different units is different.
Entropic endurance can be measured and involves stochastic variability.
Accordingly, if the entropy generation,
In this definition, the evaluation of damage is relative to the initial damage, which is not necessarily zero. The initial entropic damage can be calculated using the correlation between the rate of damage and damage at different stages of degradation. 50
Using the entropy-based damage approach allows analysts to develop the corresponding structural integrity and assess the remaining useful life (RUL) of such structures. This recognition lets us to introduce an entropy-based prognostics and health management (PHM) framework which provides a strong foundation for failure prediction and RUL assessment of critical components and structures. The approach builds on the beneficial characteristics of the entropy-based damage model and conventional data-driven PHM frameworks. 51
PHM approaches provide information about the performance and RUL of components by modeling degradation propagation. The techniques included in the PHM provide warnings before failures happen, reduce life cycle cost of inspection, and improve qualification tests assisted in design and manufacturing. Generally speaking, PHM models51,52 are developed to combine the knowledge of the past trend and the current state of degradation with an endurance threshold to estimate the RUL. 53 Traditionally, these methods rely on the condition of the performance data, which are evaluated by empirically based physics-of-failure (PoF) models54–58 or data-driven techniques such as machine learning algorithms.59–61 The PoF models utilize the specific knowledge of products, such as failure mechanisms, material properties, loading profile, and geometry.54–58 For example, in the PoF modeling framework, the corrosion-fatigue damage process is divided into four stages: (1) the pit nucleation, (2) growth of the nucleated pit, (3) surface crack initiation and its growth into a through-crack phase, and (4) through-crack growth to a prescribed critical length. For each of these stages, researchers have proposed empirical models for the life estimation.62–69 In these models, the pit growth rate is widely followed by the Kondo 62 model, whereby the pit remains hemispherical in shape and grows at a constant volumetric rate. The linear elastic fracture mechanism (LEFM) and the elastic–plastic fracture mechanism (EPFM) are often applied for analysis of the crack propagation stage and final instability.57,65,70 However, such empirical methods are limited to simple failure mechanisms and are hard to model when multiple competing and common-cause failure mechanisms are involved.
The data-driven methods such as neural networks, 59 decision tree classifiers, 60 and Bayesian techniques71,72 can obtain the complex relationship and the degradation trend in the data without the need for the particular product characteristics such as the degradation mechanism or material properties; however, they cannot identify the differences between various failure modes and mechanisms. In data-driven approaches, RUL is typically assessed by estimating damage trend, performing an appropriate extrapolation to the damage trend, and calculating RUL from the intersection of the extrapolated damage to a defined endurance threshold. In these methods, the threshold, such as the mean value of endurance thresholds, or a conservative lower bound of endurance thresholds, is commonly selected deterministically. However, the specification of the endurance threshold has a critical effect on the RUL prediction performance. 73 For example, setting a high threshold value can cause a component to fail before its scheduled time. However, a conservative low-threshold value can result in the premature abandonment of the component from the operation.
Using entropy generation as an index of damage can remedy some of the shortcomings of PoF and data-driven approaches in several ways. In comparison with PoF models, which traditionally rely on the dominant failure mechanism, entropy generation provides a consistent and far broader interpretation and measure of damage by incorporating multiple competing and common-cause degradation mechanisms. By its nature, entropy “as a measure of uncertainty”28,31,74 comprises the variabilities in the microstructure states (e.g. chemical composition phases, point and planar defects, grain orientations, grain boundaries, dislocations, and configurations of micro-cracks and cavities). Therefore, it can provide a better estimate of life in the absence of detailed microstructural observations. 31 Furthermore, the entropy-based PHM framework utilized in this article enables us to tackle the issue of endurance threshold determination by taking advantage of entropy as a thermodynamic state function independent of the path of the failure (which commonly depends on factors such as geometry, load, and frequency of load) from the initial state to the final failed state of the material, considering a known endurance limit.34,43
The remainder of this article is organized as follows: in section “Entropy generation and damage,” the necessary thermodynamic background for assessment of entropic damage is presented. Section “Entropic assessment of damage in corrosion fatigue” describes the entropy generation function for evaluating corrosion-fatigue degradation. Section “Entropy-based prognostic framework” proposes an entropy-based PHM framework to predict the RUL of critical components and structures. Section “RUL prediction of AL samples subjected to corrosion fatigue” employs the entropy-based PHM method for RUL prediction of AL 7076-T651 samples subjected to corrosion fatigue, followed by conclusions in section “Conclusion.”
Entropy generation and damage
Using entropy as a broad index of damage, in order to assess damage induced in the system, every fundamental dissipative process active in the system should be identified. Using the conservation (balance laws) of energy, mass, and momentum, along with the Gibbs
75
relation and under the hypothesis of local equilibrium,33,37,76,77 the entropy generation,
Examples of different thermodynamic forces and corresponding fluxes in some dissipative processes.
Qualification of the local equilibrium hypothesis for slow degradation mechanisms enables the use of equation (4) and allows us to express the evolutionary trend of the total volumetric entropy generation,
We must note that all degradation mechanisms leading to damage are in non-equilibrium states. Therefore, equilibrium thermodynamics, which is only concerned with equilibrium conditions, whereby a system achieves its maximum entropy with zero entropy generation rate, σ = 0, cannot describe the degradation processes which are commonly in non-equilibrium conditions. The local equilibrium hypothesis, however, enables us to deal with non-equilibrium changes from place to place and over time.34,49 Two consequences of the local equilibrium hypothesis are that entropy “as a damage index” remains a valuable state function even in non-equilibrium situations and that it leads to an explicit expression for the entropy production from which dissipative processes in the bilinear form of thermodynamic forces and fluxes are inferred. However, the hypothesis is not appropriate for describing fast- or high-frequency phenomena such as ultrasound propagation, shock waves, and nuclear interactions. 33
Entropic assessment of damage in corrosion fatigue
Figure 1 shows the setup for an experimental process that involves a component subjected to the corrosion-fatigue degradation mechanism and being electrochemically monitored. The corrosion process occurs on the surface of the specimen under a cyclic load with the following oxidation and reduction reactions for a metallic electrode,
where

Entropy flow in the control volume under corrosion fatigue.
The entropy generation function for the corrosion-fatigue degradation mechanism in terms of the main dissipative processes may be expressed as follows. 49 Consider a specific control volume: the main dissipative processes during the corrosion-fatigue failure mechanism include heat losses (the first term in equation (8)), electrochemical losses induced by electronic current leakage resulting from activation overpotential (the second to fifth term in equation (8)) and Ohmic overpotential (the sixth term in equation (8)), diffusion losses (the seventh and eight terms in equation (8)), and mechanical losses (the ninth and tenth terms in equation (8)). Equation (8) expresses the contribution of these factors to the entropy generation in a specific control volume as
where
Using this experimental setup, corrosion-fatigue tests were performed in the high-cycle fatigue regime under the following assumptions for quantification of associated entropy generation during corrosion fatigue: (1) the corrosion process occurs close to the equilibrium state i.e. open circuit potential (OCP), and the classical theory of irreversible processes 33 does not deal with situations far from equilibrium; (2) entropy flow due to heat exchange is negligible in accordance with Naderi et al. 43 and Ontiveros, 44 whose results showed insignificant entropy generation due to heat conduction inside the solid under high-cycle fatigue. Although the high-cycle fatigue does not result in significant heat conduction, low-cycle fatigue can experience noticeable thermal heat generation and conduction; 79 (3) diffusion losses are eliminated since the solution species are assumed to be well mixed; (4) the effects of diffusivity and concentration of hydrogen at the crack surface are excluded from entropy calculation in accordance with Manson’s 54 study, indicating that above a critical frequency of load cycling (of 0.001 Hz), shorter cycle periods provide less time for diffusion and accumulation of hydrogen; (5) the Ohmic overpotential is ignored due to the placement of the Luggin capillary close to the working electrode; and (6) the corresponding current to OCP records are obtained based on the fit of the sum of two exponential functions to the part of the polarization curve where OCPs vary (i.e. around the second breakdown potential) in the corrosion-fatigue experiments. 49 According to these assumptions, the amount of entropy generated in each cycle of the corrosion-fatigue experiment was obtained in terms of the area of the hysteresis energy loops resulting from mechanical stress–strain 80 and the energy resulting from variation of corrosion potential–electrical charge in each cycle of loading and unloading. 49
Figure 2(a) represents the results of multiple entropic damage tests showing the evolution trend of corrosion-fatigue volumetric entropy for different loading conditions. Markers show the final value of fracture corrosion-fatigue entropy. The results show that the cumulative entropy values at the fracture points remain between 0.7 and 1.5 MJ m−3 K−1. Figure 2(b) presents the box plot of the fracture entropies for each loading condition. The slim distribution of entropy to failure data points justifies the ability of entropy in treating the uncertainties associated with microstate variabilities.
49
Furthermore, it reveals the independence of entropy to the loading condition (i.e. failure path). Existence of the fracture entropy distribution can be interpreted as the presence of uncertainties such as instrumental measurement errors, legitimacy of assumptions considered in entropy evaluation, weak control of the experimental, operational and environmental conditions, and human error. Using

Volumetric entropy generation evolution. In the Figure 2(a), P represents the tensile stress.

Damage evolution trends.
Entropy-based prognostic framework
In this study, we propose an entropy-based PHM framework to estimate the RUL of critical components and structures. This framework, shown in Figure 4, benefits from conventional data-driven techniques and the superior advantages of the use of entropy as a damage index in comparison with common observable markers of damage.

Entropy-based damage prognostic framework.
According to this framework, the health assessment and RUL prediction can be implemented in four steps. First, the dissipative processes and associated data in the critical components under aging are determined. These processes and relevant parameters can be identified by the traditional failure modes, mechanisms, and effects analysis (FMMEA), which identifies the potential failure mechanisms for products under certain environmental and operating conditions. The volumetric entropy as a parameter of damage, which includes all the interactive failure mechanisms, is quantified then. The second step is the extraction of the features (e.g. root mean square, mean, variance, kurtosis, crest factor, peak to peak, auto correlation, and cross correlation) related to the monitored variables (e.g. temperature, electrical potential, corrosion current, mechanical stress, and solution concentration) used to measure dissipative processes and consequent entropy value, detection, and isolation of anomaly states by fault diagnosis approaches. The fault diagnosis approaches can include artificial neural network (ANN),60,81,82 wavelet theory,83,84 relevance vector machine, 85 support vector machine (SVM), 86 order tracking method, 87 and instantaneous power spectrum statistical analysis. 88 The third step includes the use of prognostic techniques (e.g. ANN, 59 fuzzy logic, 89 wavelet theory, 81 relevance vector machine, 85 Bayesian methods (e.g. Kalman filter and particle filter (PF)71,72), regression, 90 demodulation, 91 time series analysis, 92 and fusion prognostics techniques 93 ) for entropy-based damage model parameter estimation and prediction purposes.
The fourth and final step is RUL prediction. As discussed earlier, endurance threshold determination is critical in the performance of RUL prediction.73,94 The entropy-based PHM framework suggested in this study handles the uncertainties related to endurance threshold by taking the advantage of entropy characteristics (e.g. providing the best estimate of degradation evolution and being independent of the failure path as a thermodynamic state function).
RUL prediction of AL samples subjected to corrosion fatigue
This article adapts the entropy-based damage evolution data from Imanian and Modarres’s 49 study to obtain the RUL of the AL7075-T651 coupons subjected to corrosion fatigue. Tests were carried out at the maximum cyclic stress chosen as 87%, 80%, 70%, 63%, and 57% of the corroded specimen yield stress, the stress ratio of 0.01, and the loading frequency of 0.04 Hz. The stop criterion was the failure of specimens. Figure 5 shows the implementation of the entropy-based PHM framework to the corrosion-fatigue degradation mechanism. The health assessment and RUL estimation of aluminum samples are demonstrated in two steps using two well-known methods in PHM. These steps are offline data analysis, including processing of training data such as extraction of features and classification of training data points to healthy and faulty labels, and online data analysis, including health assessment of test data points, estimation of online degradation model parameters, and RUL prediction. The details of each phase are discussed below.

The application of the entropy-based PHM approach.
Offline modeling
Assuming that the corrosion and fatigue are the main dissipative processes, features of strain and corrosion potential variations (i.e. maximum amplitude, root mean square, and skewness) along with energy loop values resulting from stress–strain hysteresis and corrosion potential–electrical charge hysteresis were processed for training samples. Three specimens were selected as training samples from each test condition having a group of five specimens. To evaluate the state of the component’s health, the k-nearest-neighbor (k-NN) 95 approach was employed.
Selection of k-NN among other classification methods such as ANNs, 96 SVMs, 97 and fuzzy classifiers (FCs) 98 is based on the fact that k-NN needs the least number of assumptions regarding the data. For example, ANN and SVM require to be parametrized, especially if certain assumptions are needed to be made (e.g. if the data is linearly separable in SVM). Nevertheless, in a separate research, authors considered to evaluate the performance of k-NN compared with other classifiers mentioned above, as they can outperform each other based on the data set characteristics.99,100 The k-NN approach is a non-parametric classification method in which a new data point is classified based on its vicinity to the k-neighboring data points of known classes. In this method, the Euclidean distance is calculated from a new data point to the centroid of the nearest neighbors from each class, and the new data is classified to the shortest centroid distance. The healthy and faulty classes were selected based on the portion of the specimen life. Data from the initial 10%–30% and 70% lifetime defined the healthy and faulty class, respectively.
In k-NN, the selection of k is subject to the type of data. To obtain a proper value of k at which the classification output is stabilized, sensitivity analysis should be performed. In this study, the impact of the presence of noise in the classification was examined in altering the location of the nearest neighbors’ centroid. In the sensitivity analysis, the change in the distance of a test point to the centroid of the nearest neighbors was calculated as the value of k was increased by one at each step, as shown in Figure 6. k was selected such that enough neighbors were selected to produce a stable neighbor’s centroid.

Selecting k based on sensitivity study for the training data.
Fault level (FL) is determined by taking a moving average of the k-NN binary output over time with a window of 100 data points. This enables the definition of anomaly detection threshold. An anomaly was acknowledged when the FL > 0.7. Based on the faulty portion of the lifetime, this threshold is defined as a minimum operating requirement level for an item (the item may still be operational, but not satisfactory). An example of a classification output plot is depicted in Figure 7.

An example of a k-NN classification output.
Online modeling
By choosing the normalized volumetric entropy parameter as the index of damage to be monitored, the Bayesian updating method can be used to estimate the parameters of the degradation model and to predict age. The Bayesian updating approach provides a general, rigorous method for dynamic state estimation problems. The idea is to build a probability density function (PDF) of the system model states (i.e. the degradation model parameters in this article) based on all available information. PF is a method for implementing a recursive Bayesian filter using Monte Carlo simulations, where in comparison with other approaches (extended Kalman filter, Kalman filter, and unscented Kalman filter), PF provides a better performance when either the system dynamics or noises are nonlinear. Herein, it approximates the damage evolution trend (Figure 3) parameters’ PDF by a set of particles sampled from the normal PDF and a set of associated weights denoting probability masses. 71
As the PF requires models developed from known system behavior, the particles in the PF method are generated and recursively updated by the system states evolution model shown in equation (9), where the random walk model is considered for evolution of degradation trend parameters
where
The goal of PF implementation is then to estimate the probability distribution of the degradation model parameters,
Within the Bayesian framework, the posterior PDF of the system model states,
where
Sampling/importance resampling is a commonly used algorithm to attribute importance weight,
A recursive formulation in order to update the weights can be obtained as
where the importance distribution
RUL prediction
As discussed earlier, given a series of measured entropic damage values,
The estimated PDF of the damage prediction,
Defining a deterministic entropic endurance which corresponds to endurance threshold level as
The PDF of the RULs at cycle
However, if the variability in the endurance threshold is considered independent of the deterioration distribution, Abdel-Hameed
101
showed that the RUL cumulative density function at cycle
Prognostic results
Using the k-NN classification method, anomalies were identified when the test data FL crosses the anomaly detection threshold. Once the anomaly was detected, the PF algorithm was initiated to predict the RULs. The initial values of the degradation evolution trend parameters

(a) An example of RUL prediction and (b) distribution of the degradation trend at the anomaly criteria.
The same procedure was applied to the remaining specimens. The percentage of the mean of the predicted RULs and actual RULs and the minimum and maximum absolute value of error between the predicted RULs at different endurance thresholds and actual RULs are shown in Table 2. The error interval falls in the range of 0.1%–25.5%. This implies that using entropy as a measure of damage can handle the uncertainties related to the endurance threshold to a good extent. The performance of predictions, however, can be enhanced by reducing the uncertainties discussed earlier. Other studies such as Nystad 73 and Fantoni and Nordlund 103 suggest statistical solutions to manipulate the endurance threshold uncertainty in predicting the RUL, or using the opinion of experts with expertise within the relevant field to provide useful information about the threshold probability distribution. With the use of entropy as a damage index, however, this study tried to provide a practically scientific approach to control the endurance threshold uncertainty effects on RUL prediction.
RUL prediction results.
RUL: remaining useful life.
Conclusion
This article discusses an entropy-based definition of damage and its uses in a PHM approach to estimate RUL prediction of critical components and structures within the irreversible thermodynamic framework. The suggested approach not only offers a science-based foundation for PHM methods but also benefits from the superior advantages of the use of entropy as a damage index, in contrast with the common observable markers of degradation. Entropy as a measure of damage provides a unified non-empirical model of degradation evolution that can capture the effect of multiple competing and common-cause failure mechanisms. Entropy as a measure of uncertainty can provide the best estimate of the degradation evolution trend in the absence of microstructural observation. Furthermore, entropy as a thermodynamic state function is independent of the failure path and is capable of handling the uncertainties related to endurance thresholds.
The proposed prognostic approach was demonstrated by modeling the corrosion-fatigue degradation mechanism and tested in a laboratory environment. The thermodynamic entropy generated in the corrosion-fatigue degradation mechanism was then used to predict the RUL of AL 7075-T651 samples. Combining entropy-based damage theory with real experimental data and the standard PHM framework, we showed that the structural integrity of complex failure mechanisms such as corrosion fatigue can be evaluated scientifically, experimentally, and statistically. Furthermore, using the entropy as a degradation measurement enables us to handle the uncertainties associated with the endurance thresholds and to efficiently predict RUL. Although the RUL prediction results support the adequate performance of the entropy-based damage model, the advantages of this method remain to be explored for different materials, various operational conditions (e.g. different loadings, loading frequencies, loading ratios, and activation overpotentials), different geometries, various environmental conditions (e.g. different temperatures, corrosive solution pHs, solution flow rates, and radiations), and other failure mechanisms.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was a part of ongoing research through grant number N000141410005 from the Office of Naval Research (ONR).
