Next Article in Journal
Experimental and Numerical Investigation on Ballast Flight from Perspective of Individual Particles
Next Article in Special Issue
Numerical Study of the Comparison of Symmetrical and Asymmetrical Eddy-Generation Scheme on the Fire Whirl Formulation and Evolution
Previous Article in Journal
Seismic Response Analysis of Multi-Story Steel Frames Using BRB and SCB Hybrid Bracing System
Previous Article in Special Issue
Combustion Process of Canola Oil and n-Hexane Mixtures in Dynamic Diesel Engine Operating Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Similarity Analysis for Time Series-Based 2D Temperature Measurement of Engine Exhaust Gas in TDLAT

1
Research & Business Foundation, Dong-A University, Busan 49315, Korea
2
Technical Center for High-Performance Valves, Dong-A University, Busan 49315, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(1), 285; https://doi.org/10.3390/app10010285
Submission received: 24 October 2019 / Revised: 23 December 2019 / Accepted: 29 December 2019 / Published: 30 December 2019
(This article belongs to the Special Issue Progress in Combustion Diagnostics, Science and Technology)

Abstract

:
As regulations on the emission of pollutants from combustion systems are further tightened, it is necessary to reduce pollutant species and improve combustion efficiency to completely understand the process in the combustion field. Tunable diode laser absorption tomography (TDLAT) is a powerful tool that can analyze two-dimensional (2D) temperature and species concentration with fast-response and non-contact. In this study, stabilized spectra were implemented using the mean periodic signal technique to enable real-time 2D temperature measurement in harsh conditions. A time series statistical-based verification algorithm was introduced to select an optimal spectral cycle to track 2D reconstruction temperature. The statistical-based verification is based on the Two-sample t test, root mean square error, and time-based Mahalanobis distance, which is a technique for similarity analysis between thermocouple and reconstruction temperature of 18 candidate cycles. As a result, it was observed that the statistical-based TDLAT contribute to improving the accuracy of time series-based 2D temperature measurements.

1. Introduction

With global warming prevention and environmental protection emerging as an issue all over the world, emission regulations are gradually being tightened as the effects of emissions from combustion systems on air and human bodies are highlighted. For example, EURO 6 regulations are an effort to reduce emissions such as CO, HC, NOx, and PM by the auto industry. The recent EURO 6d emission regulations require real driving emissions (RDE) tests, which measure pollutants emitted from motor vehicles while they are driven, as an additional approval. As a result, demand for vehicles with low NOx emissions and high engine power has increased, and several studies have been conducted to improve the performance of NOx after treatment systems [1,2,3]. An important element in catalytic after treatment devices is exhaust gas temperature because catalytic performance is a function of temperature. A car engine running at high temperature tends to emit less CO2 and more NOx. On the other hand, a lower engine temperature will result in more CO2 and less NOx. However, the challenge of reducing both CO2 and pollutant emissions simultaneously still remains. Therefore, analysis of temperature and species concentration during combustion is necessary for fundamental understanding of combustion reaction, which may reduce pollutants in the combustion system and improve combustion efficiency. However, the thermocouple (TC), a conventionally used one-point temperature measurement device, has difficulty in analyzing non-uniform flow, such as in engine exhaust. In addition, the techniques such as coherent anti-stokes Raman scattering (CARS) [4,5] and laser induced grating spectroscopy (LIGS) [6], which are recently used in the spectroscopy field, can measure temperature at a fast response speed, but have a limit of point measurement, and filtered Rayleigh scattering (FRS) [7] and planer laser induced fluorescence (PLIF) [8,9] enable two-dimensional temperature, but their measurement systems are expensive and complicated. Tunable diode laser absorption tomography (TDLAT) is a suitable tool for understanding the combustion process by allowing 2D analysis of temperature and species concentration with non-contact fast-response. Therefore, it has been developed with various TDLAT techniques for measuring 2D temperature and species concentration in combustion fields [10,11,12,13,14]. Deguchi et al. demonstrated a 2D temperature distribution measurement by TDLAT on Bunsen-type flame burners and gasoline engines [10]. Ma et al. developed a hyperspectral tomography (HT) method that can simultaneously measure 2D distribution of temperature and H2O concentration with a temporal resolution of 50 kHz at the exhaust plane of a j85 engine [11]. Busa et al. developed a measurement method of combustion efficiency that combines TDLAT and particle image velocimetry (PIV) in the scramjet engine [13]. However, these methods cannot be applied directly for real-time measurements because of excessive computational time and lack of time resolution assessment required for computed tomography (CT) analysis.
This study introduced a spectral cycle parameter improving reconstruction performance with a posteriori statistical technique in conventional absorption tomography to measure the exhaust gas temperature in milliseconds through TDLAT. Signal stabilization technique was applied in various cycle models to obtain stabilized absorption spectra in harsh environments, and statistical verification using TC measurement was performed to evaluate the results of the time series 2D reconstruction temperature distribution.

2. Theory of TDLAS

Tunable diode laser absorption spectroscopy (TDLAS) is a measurement technique commonly used in a variety of applications recently as well as temperature and density of molecules in gas [15,16]. When a laser beam permeates an absorption medium, its intensity will be attenuated by the absorbing gas species. This phenomenon is measured by TDLAS system that continuously scans laser wavelengths. The principle of TDLAS is theoretically expressed in terms of absorption by Beer-Lambert law that describes the relationship between the incident intensity Iv0 and transmitted intensity Iv as follows:
I v   I v 0 = e x p A v   = e x p 0 L α v s d s ,
where A ν denotes the spectral absorbance of representative wavenumber ν . L denotes the optical path length, and α ν is the local mass-weighted absorption coefficient at position s along the laser beam, from 0 to L. This quantity is given by the following expression:
α v = i x i p k T j S i , j T ƒ v v ; v j ,   Δ ν G , Δ ν L   ,
where, k [J/K] denotes Boltzmann’s constant, ƒ v is the Voigt profile, with Lorentzian and Gaussian half width at half maximums (HWHMs), Δ ν L and Δ ν G [17,18,19]. p denotes pressure, assumed to be uniform and xi denotes mole-fraction of absorbing species i, where, x i , p, and T are evaluated locally. Si,j(T) denotes temperature dependent line strength of absorption transition j and can be calculated by a function of temperature as follows:
S i , j T = S i , j T 0 Q T 0 Q T T T 0 e x p h c E i k 1 T 1 T 0 1 e x p h c v 0 , i k T 1 e x p h c v 0 , i k T 0 1 ,
where, h (J·s) and c (cm/s) denote Planck’s constant and the speed of light, v0 (cm−1) denotes the line-center wavenumber, Ei (cm−1) denotes lower-state energy of the transition, T0 (K) is the reference temperature (296 K), and Q(T) denotes the partition function of molecular absorption at a particular temperature. It is usually expressed as a 4th order polynomial function [20]. Although it is difficult to accurately describe the chemical formula of gasoline, iso octane (C8H18) was selected because it is the hydrocarbon component with similar thermal properties to the actual composition. The combustion reaction of iso octane is as follows:
C 8 H 18 +   12.05 O 2 + 3.76 N 2   8 C O 2 +   9 H 2 O   +   12.05 · 3.76 N 2 .
Among the combustion products, nitrogen (N2) is an inert gas and hence, there is no chemical or absorption reaction. Therefore, water-vapor (H2O) was chosen for exhaust gas diagnosis because of its high mole fraction and high absorption strength.
Figure 1 shows the theoretical H2O absorption spectra calculated using the HITRAN2012 molecular spectroscopic database [21]. The patterns of the absorption spectra can be observed as the temperature changes. In this study, three peaks lines (7202.255 cm−1, 7202.909 cm−1, and 7203.891 cm−1) with high line strengths and no interference with other molecules were chosen for calculating reconstruction temperature.

3. Experiment and Tomography Analysis

A schematic diagram of a TDLAS system measuring engine exhaust temperature is shown in Figure 2. The OHC (over head camshaft) gasoline engine (EX17, Subaru Robin, Doylestown, PA, USA) is driven in idle mode to produce emissions. The exhaust pipe installed 5 mm below the measurement cell was 20 mm in diameter and 150 mm in length. The distributed feedback (DFB) laser (NLK1E5GAAA, NTT Electronics, Japan) for H2O gas target is tuned to the near-infrared spectral region at 1388 nm by a laser diode controller (LDC-3900, ILX Lightwave, Bozeman, MT, USA). The modulated laser beam is split into 16 channels by a 1 × 16 fiber optic coupler (SMF-28e, OPNETI, NY, USA). Among them, the 10 beams are irradiated into 10 collimators (50-1310-APC, THORLABS, Newton, NJ, USA) and the transmitted optical signals are detected by 10 photodetectors (G12180-010A, Hamamatsu Photonics, Japan). The received laser signals are amplified and filtered by low pass filters to remove signal noise. The amplified signals are transmitted into a data recorder (8861 Memory High Coda HD Analog16, HIOKI, Japan) in real time at a sampling rate of 500 kHz. Simultaneously, the real temperature for verification is measured by a B-type TC at the sampling rate of 1 kHz at the center of the target area. After absorption, signals collect by the data acquisition card are sent to the computer for reconstruction calculation.
Multiple parallel laser beams pass through the region of interest. The intersected spectral information is used to reconstruct 2D temperature by CT analysis. Figure 3 shows the laser paths and grids in the measurement region. The laser path consists of 10 real beams in a 5 × 5 layout and eight interpolated virtual beams The interpolated virtual beams imply the virtual signal data interpolated from the signal of the actual beam, and in the signal application field a widely introduced cubic spline interpolation is used in this study [22]. The cell length through which the laser passes in the measurement area is 60 mm and reconstruction area is 40 mm, and the spatial region for CT analysis consist of 81 (9 × 9) square grids with a grid size of 2.0 mm × 2.0 mm.
Unexpected error signals due to harsh conditions were removed to reduce reconstruction calculation time for converging. The formula is as follows:
δ t = n = 1 N I ¯ v 0 , i I ¯ v , i t   I ¯ v 0 , i I ¯ v , i t 1 ,
where, the subscript t denotes a cycle at t time and n denotes the data number in a cycle. δ denotes the error rate, which is a ratio of absolute error sum in I ¯ v and I ¯ v 0 of the cycle to that of the previous cycle. Here, the signal cycle was eliminated when δ t is greater than 1.5. This method effectively reduced the error when converting the signal into the spectra. The noises generated by beam steering, window foul, and thermoelectric noises were approximated by Gaussian noise modeling, and the stabilized absorption spectra was obtained by the polynomial reduction method (PRM) averaging a number of absorption signal cycles [10,23,24] Although averaging a significant number of absorption cycles provides more stable signal, it has the drawback of degrading time-resolution and distorting raw signals. Hence, it is necessary to investigate appropriate number of cycles with good reconstruction performance and fast time-resolution simultaneously.
Figure 4 shows the result of mean transmitted signals ( I ¯ v ) and mean absorption spectra using a total of 201 cycles (1 reference cycle, 100 cycles in front and 100 cycle in back) measured by the TDLAS system. Here, one cycle of measured signal consists of 500 samples with a recording time of 1 ms. A total of 5 million signals were received during 10 s of engine operation. Figure 4a is the result of averaging each raw signal acquired from the photodetector at the same time when 10 lasers passed through the target area. To convert these signals into absorption spectra, I ¯ v and the mean incident signal ( I ¯ v 0 ) must be known as shown in Equation (4). I ¯ v 0 , also known as the baseline, was estimated by PRM, which is a least squares fitting using 5th order polynomial curve with I ¯ v [23,25]. To facilitate curve fitting, the saw-tooth waveform had a frequency of 1 kHz and a sampling rate of 2 μs. All measured absorption spectra generated from I ¯ v and I ¯ v 0 is shown in Figure 4b. These spectra consisted of 311 wavenumbers. Laser group 0 and 1 represent a group of vertical and horizontal lasers, respectively. The spectra pattern can be identified by different laser numbers. As shown in Figure 1, the measured absorption spectra are determined by the temperature of the measurement area, so the reconstructed temperature distribution in measurement zone can be traced by a tomographic reconstruction method from the information of each absorption spectrum.
The determination method of cycle and temperature is illustrated in Figure 5. First, raw absorption signals, reference time i and candidate cycle j are input. i denotes analysis time point at engine running time. A total of 990 absorption signals were used in 0.1 s from 1.0 to 9.9 s. j denotes the number of cycles of absorption signal to be averaged. Eighteen candidate cycles were chosen to evaluate the analytical performance in various cycles (31, 51, 71, 91, 111, 151, 201, 251, 301, 351, 401, 451, 501, 601, 701, 801, 901, and 1001 cycles).
In the second step, the PRM and δ value are combined to obtain mean absorption spectra, described in more detail above. This technique was able to remove noise more effectively than the conventional PRM. Reconstruction calculation consisting of 9 × 9 array was calculated using mean absorption spectra. The reconstructed temperature distribution was tracked by repeated calculations until errors in the experimental and theoretical absorption spectra were minimized by the multiplicative algebraic reconstruction technique (MART) algorithm [24].
α v n , m k = α v n , m k 1 × A v L n , m m = 1 M n = 1 N α v n , m k 1 β L n , m ,
where, the superscript k represents the number of iteration and the subscript (n, m) is the gird number. β is the relaxation parameter that affects the rate of convergence and was set to 0.1. The initial values of α v n , m at each grid are an important factor for fast convergence calculation, which were approximated to target values by adopting the multiplication line of sight (MLOS) method [26]. The reconstruction calculation was terminated when difference of α v n , m k and α v n , m k 1 was less than 1 × 10−6.
Figure 6 shows a comparison of the measured absorption spectra with the calculated absorption spectra in the engine running time of 5.0 s. The average root-mean-square error (RMSE) at 311 wavenumber were 9.766 × 10−4, 4.485 × 10−5, 1.273 × 10−3, and 2.144 × 10−3, respectively, in Figure 6a–d. Measured spectra were in good agreement with spectra by the theoretical calculation, but the measured spectra for the 51 and 201 cycles were slightly different from each other, and their reconstruction calculations were also differed slightly. This implies that the number of mean absorption cycles is an important parameter for the reconstruction results.
In the third step, the reconstructed temperatures, TRec(i,j), calculated for each candidate cycle and time are output. To estimate the TRec(i,j) from the finally converged α v n , m at all grids, the theoretical absorption coefficients at the temperature has already been calculated using the HITRAN database for the mole fraction 1.0. The correlation between the theoretical absorption coefficient, α v , t h e , and the reconstructed absorption coefficient, α v , was evaluated to find the value with the most similar pattern, i.e., the nearest correlation coefficient to 1.0, and to obtain TRec(i,j). The correlation coefficient, C n , m , at each grid is described as follows.
C n , m = v = 1 n α v α ¯ α v , t h e α ¯ t h e v = 1 n ( α v α ¯ ) 2 v = 1 n α v , t h e α ¯ t h e 2 ,
where α ¯ and α ¯ t h e denote the average value each of the theoretical and reconstructed absorption coefficient.
Figure 7 shows the results of 2D reconstructed temperature distribution using 51 and 201 cycles at the engine running time of 5.0 s. Comparing Figure 7a,b, the temperature difference at the center point was 12.35 K and the RMSE at all grids was 5.87 K. This is an example of the highest temperature difference compared to other engine running times. The characteristics of the temperature distribution for each time were similar, but the temperature fluctuation was shown differently. From these results, it was found that the number of mean cycles affected the results of CT analysis. Additionally, the output T ¯ Rec was set to the average value of nine grids located in the center of the measurement area as the TC was installed in the center of the exhaust pipe and the temperature was measured at one point. After all, the output TRec consists of a set of elements T ¯ Rec for each i and j. The number of TRec was 17820 in the i × j array (990 × 18).
In the fourth step, the statistical-based analysis was performed using TRec and TThc, which means temperature measured by TC and was used to determine the cycle and TRec most similar to TThc. Further details are described in Section 4.
The computation time of each step on an Intel Core i7 2600K CPU with 16GB RAM was as follows. The CPU time of Process I to convert signal data into mean absorption spectra was 1.38 s and 1.90 s (311 wavenumbers) respectively at 51 and 201 cycles. As the number of mean cycles increased, the CPU time increased linearly. The CPU time of Process II was scored in under 1.79 s (311 wavenumbers), mostly converging on less than 150 iterations.

4. Statistical-Based Mean Difference Test and Similarity Analysis

The statistical-based analysis is performed to identify cycles that represent TRec similar to TThc. The 18 candidate cycles were selected using suitable intervals for analysis from 31 to 1001 cycles. The procedure is shown in Figure 8. Paired t test, which is a mean difference testing method, is conducted to evaluate the similarity between the 18 candidate cycles [27].
The next step is followed by performing a mean difference test and similarity analysis between TThc and TRec for all 18 candidate cycles. The mean difference test uses the representative statistical technique of Two-sample t test. RMSE and time-based Mahalanobis distance (TMD) were used for similarity analysis. RMSE is a measure of generalization of standard deviations that indicates the difference between actual and estimated values [28]. The Mahalanobis distance (MD) represents the distance between data and considers covariance [29].
Finally, the cycles that represent TRec most similar to TThc based on the analysis results were determined.

4.1. [Step 1] Mean Difference Test

TRec derived from 18 candidate cycles consists of 0.1 s intervals from 1 to 9.9 s as shown in Table 1. To test the mean temperature difference in 18 candidate cycles, the analysis was performed on nine TRec(i,j), which is the reconstructed temperature of ith reference time of the j-cycles. The total number of sets comparing TRec(i,j) in 18 candidate cycles to each other for the paired t test were 153. For the actual test, the analysis was performed only at the reference time (1 s, 2 s, 3 s, 4 s, 5 s, 6 s, 7 s, 8 s, and 9 s) for all 153 sets.

4.1.1. Paired t Test

The paired t-test is useful for analyzing the same set of items that were measured under two different conditions, i.e., differences in measurements made on the same subject before and after a treatment, or differences between two treatments given to the same subject.
The paired t-test is based on the test statistic as shown below:
T p t = D ¯ d 0 S D / n ,
where,
D ¯ = k = 1 n T R e c i , A k T R e c i , B k n ,
S D = k = 1 n T R e c i , A k T R e c i , B k D ¯ 2 n 1 ,
where, D ¯ denotes the mean of the differences between the two kinds of cycles.   d 0 denotes a constant that can be set depending on whether the mean of the difference between the two groups is zero or a specific non-zero value. S D denotes the standard deviation of differences between two kinds of cycles. TRec(i,A)k denotes the kth reconstruction temperature at ith time (s) of cycle A, i = 1.0 , 1.1 ,   1.2 ,   ,   9.9 . TRec(i,B)k denotes the kth reconstruction temperature at ith time (s) of cycle B. n denotes the number of TRec in ith time (s).
For paired t test, the hypothesis is as follows:
  • Null hypothesis: H 0 :   D ¯ = d 0 (the population mean of the difference ( D ¯ ) equals the hypothesized mean of the difference ( d 0 = 0 ).
  • Alternative hypothesis: H 1 :   D ¯ d 0 (the population mean of the difference ( D ¯ ) does not equal the hypothesized mean of the difference ( d 0 = 0 ).
Here, the null hypothesis is rejected when the p-value is less than or equal to 0.05, i.e., when there is no difference between the mean of the two groups.

4.1.2. Results of Mean Difference Test in [Step 1]

Table 2 is a summary of the paired t-test results. For 153 sets, 1377 test results would be obtained when each test is performed from 1 to 9 s. If the p-value is less than 0.05, the null hypothesis is rejected. The p-values of 1 s, 3 s, 8 s, and 9 s were less than 0.05, i.e., the means of all cases were not equal to each other. In 2 s, the p-value was less than 0.05 except in four sets (i.e., 251 vs. 301, 451 vs. 501, 451 vs. 601, and 501 vs. 601). In 4 s, only one set (351 vs. 401) showed a p-value greater than 0.05 and the remaining 152 sets had a p-value of less than 0.05. In 5 s, 6 s, and 7 s, the p-value was greater than 0.05, while the remaining 152 sets had a p-value less than 0.05. Therefore, except for eight sets among 1377 sets, the p-value in all sets was less than 0.05 and hence, there was a mean difference between the cycles.

4.2. [Step 2] Similarity Analysis

Table 3 shows the TThc(i) and T ¯ Rec(i,j) of all 18 candidate cycles used in the analysis where each T ¯ Rec(i,j) is the mean value for nine TRec(i,j)k per 0.1 s interval in Table 1. For example, in Table 3, T ¯ Rec(1.0,31) is 406.418, which averaged reconstruction temperatures from TRec(1.0,31)1 = 404.833 to TRec(1.0,31)9 = 404.999 in Table 1. Here, three methods were used for similarity analysis. First, the Two-sample t test was used to test the difference between the means of the two groups. Second, RMSE analysis was performed to address the difference between estimates or models and observed values in real-world settings. RMSE analysis can be used to aggregate the differences between the residuals of two groups into one measure. Finally, in TMD analysis the correlation between two groups with different means was taken into account.

4.2.1. Two-Sample t-Test

Two-sample t-test is a widely used hypothesis test used to compare whether the averages are equal. The assumptions made while doing a Two-sample t-test include: (i) the data are continuous, (ii) the data follows normal probability distribution, (iii) the variances of the two populations are equal, (iv) two-sample is independent, and (v) both samples are simple random samples from their respective populations. The Two-sample t-test is based on the test statistic shown below:
T t w = T ¯ T h c g T = R e c g T S 2 n T g n R g n T g + n R g ,
T S 2 = i = 1 n T g T T h c i T ¯ T h c g 2 + i = 1 n R g T ¯ R e c i , j T = R e c g 2 n T g + n R g 2 ,
where, T S 2 denotes the variance of TThc(i) and TRec(i,j). TThc(i) denotes the TC temperature at ith time (s), i = 1.0 ,   1.1 ,   1.2 ,   ,   9.9 . T ¯ Thcg denotes the mean of all TThc(i) in the gth group, g = 1 , 2 ,   , 9 . T = Recg denotes the mean of all T ¯ Rec(i,j) in the gth group. nTg denotes the number of all TThc(i) in the gth group. nRecg denotes the number of all TRec(i,j) in the gth group.
For the Two-sample t test, the hypothesis is as follows
  • Null hypothesis: H0: = T ¯ Thcg T = Recg = δ 0 (the difference between the population means ( T ¯ Thcg T = Recg) is equal to the hypothesized difference ( δ 0 = 0 )).
  • Alternative hypothesis: H1: = T ¯ Thcg T = Recg δ 0 (the difference between the population means ( T ¯ Thcg T = Recg) is not equal to the hypothesized difference ( δ 0 0 )).
Here, the null hypothesis is rejected when the p-value is less than or equal to 0.05, i.e., when there is no difference between the mean of the two groups. In this study, statistical software version 18 of Minitab was used for the Two-sample t test.

4.2.2. Root Mean Square Error (RMSE)

RMSE is used to measure the difference between values predicted by a model or an estimator and the value observed. The RMSE is a measure of accuracy to compare forecasting errors of different models for a particular dataset and not between the dataset as it is scale-dependent. RMSE is always non-negative and a value of zero would indicate a perfect fit to the data. In general, a lower RMSE is better than a higher one. The formula is as follows:
T R M S E   = 1 m × i = 1 m T T h c i T R e c i , j 2 .

4.2.3. Time-Based Mahalanobis Distance (TMD)

MD was introduced by P. C. Mahalanobis in 1936 to measure distance between a point P and distribution D. It is a multi-dimensional generalization of the idea for measuring the amount of standard deviations P is away from the mean of D [30]. Among various methods for measuring distance (i.e., Euclidean distance, Manhattan distance, etc.), MD is used as it has the advantage of calculating the distance by considering weight according to the magnitude of covariance between variables [31]. Thus, MD is a method that is used to assess similarities between two groups. In this study, TMD analysis is performed to analyze similarities between TThc and T ¯ Rec by reference time. The closer the TMD value is to zero, the more likely TThc and T ¯ Rec are to be similar. The calculation of TMD procedure is as follows:
[Stage 1] Calculate the covariance matrix (S) and TMD using the following formula (14) and (15).
S = V a r T T h c i C o v T T h c i ,   T ¯ R e c i , j C o v T ¯ R e c i , j ,   T T h c i V a r T ¯ R e c i , j ,
M D = X i X ¯ T S 1 X i X ¯ ,
X i X ¯ = T T h c i T ¯ T h c T ¯ R e c i , j T = R e c i , j ,
where, T ¯ Thc denotes the geometric mean of all TThc(i), i = 1.0 ,   1.1 ,   1.2 ,   ,   9.9 . T = Rec(i,j) denotes the geometric mean of all T ¯ Rec(i,j) in jth cycles. S−1 denotes inverse covariance matrix. T indicates that the vector should be transposed.
[Stage 2] The obtained TMD values are calculated using the following formula for each time (s) interval.
T M D =   M D 1 × M D 2 × × M D n 1 n ,   n = 1 , 2 , m .

4.2.4. Results of Similarity Analysis in [Step 2]

Two-sample t test, RMSE, and TMD analysis were performed to select the cycles that indicate a TRec similar to that of TThc. The results of a normality test conducted on TThc and TRec for 18 candidate cycles prior to the Two-sample t test showed that all were satisfactory for normality. However, the equivalence test between TThc and TRec showed that some results did not satisfy the equivalence. The results of the equivalence test are summarized in Appendix A. In general, Welch’s t-test is used if the two groups do not meet equal variances. However, if the very small samples (n < 10), the Two-sample t test performs slightly better than Welch’s t-test [32].
The first Two-sample t-test results are shown in Table 4 and the null hypothesis is rejected when the p-value is less than 0.05 (see Section 4.2.1). The analysis results are explained based on the cycle, all p-values from 1 to 9 s for 31, 51, 91, 111, and 251 cycles were above 0.05. For the 71 cycles, the p-value was close to 1 at 2 s, 5 s, 6 s, and 8 s. In other words, the probability of selecting null hypothesis was high. For 151 and 201 cycles, the p-value was over 0.05 for all seconds except 5 s. On the other hand, the 451 cycles showed p-value less than 0.05 at 5 s and 9 s. The p-value of the 901 cycles was less than 0.05 at 8 s and 9 s. In the 501, 601, 701, and 801 cycles, the p-values were less than 0.05 for 5 s, 8 s, and 9 s. In the case of the 1001 cycles, the p-value of 3 s, 8 s, and 9 s was found to be below 0.05.
The result of RMSE analysis is shown is Table 5. The smaller the RMSE values, the more the similarity between TThc and TRec. First, the analysis results were based on 10 temperatures per group from 1 to 9 s (see Table 3). In terms of time, the RMSE value (1.951) of 401 cycles in 1 s was the smallest. In 2 s, the RMSE value of 1001 cycles was the smallest at 2.270. RMSE values for 601 cycles were the smallest in 3 s and 4 s. The RMSE values for 251, 201, and 701 cycles were the lowest in 5 s, 6 s, and 7 s. RMSE values for 1001 cycles were the smallest in 8 s and 9 s. The analysis using 90 temperatures from 1 to 9.9 s (see Table 3) showed that the RMSE value of 601 cycles was the smallest.
Lastly, the analysis results of TMD are shown is Table 6. The smaller the TMD value, the more the similarity between TThc and TRec. The analysis results were based on 10 temperatures per group from 1 to 9 s (see Table 3). In terms of time, the TMD value (1.626) of 351 cycles in 1 s was the smallest. In 2 s, the TMD value of 111 cycles was the smallest at 0.770. The TMD value for 151, 201, 251, and 801 cycles were the lowest in 3 s, 4 s, 5 s, and 6 s. TMD values for 201 cycles were the smallest in 7 s, 8 s, and 9 s. Next, the analysis using 90 temperature readings from 1 to 9.9 s showed that the TMD value (0.933) of 201 cycles was the smallest.

5. Conclusions

The results of this study are summarized as follows: (i) the paired t test for 18 candidate cycles indicated that, in most cases, TRec was not the same in each cycle. This implies that TRec might vary depending on the cycle set during CT analysis. Therefore, a standard for proper cycle setting is needed. (ii) Three types of similarity analysis yielded different results. First, the Two-sample t test showed that the TRec for 71 cycles was not the much different from that of TThc. Second, the results from RMSE showed that the TRec of the 601 cycles was most similar to that of the TThc. Finally, TMD results show that the TRec for 201 cycles was most similar to that of TThc. (iii) The reason for different results is that all three methodologies have different characteristics. For a more visualized description, a comparison of TThc and TRec was provided based on the similarity analysis result. Figure 9, Figure 10 and Figure 11 show a comparison of the TThc with the TRec for each 71, 601, and 201 cycles, respectively. The movement of TThc and TRec over time shows that TRec for 601 and 201 cycles moved closely with TThc. However, a closer look at the temperature change pattern shows that TRec of 201 cycles was moving similarly to TThc. On the other hand, for 71 cycle, the range of TRec changes over time was observed to be much larger than TThc. (iv) For the Two-sample t-test, it was difficult to determine similarity by sensitively reflecting changes in temperature as the difference in the mean of the overall temperature data in a given period was taken into account. Likewise, the RMSE was purely an indicator of the mean of distance between TThc and TRec, which made it difficult to analyze, including the correlation between changes in TThcand TRec. TMD was a suitable indicator for this study because it considers both the distance and correlation between TThc and TRec. Therefore, if the engine used in this study is in idle mode, setting it to 201 cycles in the CT analysis will increase the accuracy of the analysis while minimizing the analysis time and effort.
Figure 12 shows the 2D reconstruction temperature distribution in 201 cycles in time from 1.0 to 9.0 s. The temperature of the peak area of the 2D distribution in each time was equal to that of the corresponding time at the 201 cycle reconstruction temperature in Figure 11. As illustrated in Figure 12, it can be seen that the maximum temperature was continuously distributed in the central area. For example, Figure 12b,h show a reconstruction temperature of approximately 410 K and 408 K in the central area, and Figure 10 shows also the same temperature at that time. This was the same in other times.
This study presented the statistical analysis procedure and the appropriate cycle to improve TDLAT performance when measuring time series-based 2D temperature in idle mode of engine. To the best of our knowledge, this is the first attempt to present the criteria for cycle selection through the similarity analysis when based on signal stabilization technique. This signal stabilization technique was more effective in error reduction than the conventional polynomial reduction method when converting signals into absorption spectra.
However, limitations of this study include (i) it was only tested in idle mode of a particular engine, so additional experiments are required in operation mode with relatively large temperature variations; and (ii) in idle mode, the temperature variation in the peak area was uniformly measured only for that area, but it is also necessary to check the need for measurements in other areas.
To overcome these limitations, the following studies will be carried out in the future: (i) conduct a time series-based 2D temperature measurement study by applying the latest statistical techniques (i.e., linear Bayesian absorbance inference, Tikhonov regularization, and so on) to provide a fast and accurate case of reconstruction temperature calculation and (ii) typically, the reconstruction temperature is verified by TC measurements. The reliability of verification will depend on the TC measurement method, therefore studies on the TC verification methodology will also be conducted; and (iii) the size of the laser grid will affect the bias of reconstruction temperature. We would like to conduct a study on the selection of optimal grid sizes to minimize bias.

Author Contributions

Conceptualization, H.J., D.C.; Methodology, H.J., D.C.; Modeling, H.J., D.C.; Experiment, D.C.; Analysis, H.J.; validation, H.J., D.C.; Writing—original draft preparation, H.J., D.C.; Writing—Review and Edition, H.J., D.C.; Funding Acquisition, H.J., D.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2018R1C1B5085281, 2019R1G1A1010335).

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The Two-sample t test assumes a normality and equivalence. Since normality was satisfied, only the results of the equivalence test were summarized in Table A1. The equivalence is satisfied when p-value is more than 0.05.
Table A1. Result of the equal variance test between TThc and TRec for 18 candidate cycles.
Table A1. Result of the equal variance test between TThc and TRec for 18 candidate cycles.
CyclesTime (s)
123456789
P-Value
31<0.050.055<0.05<0.050.094<0.05<0.05<0.05<0.05
51<0.05<0.050.175<0.050.0620.105<0.05<0.050.205
71<0.05<0.050.127<0.050.117<0.05<0.05<0.050.071
910.071<0.050.109<0.050.2650.0730.071<0.050.099
1110.1600.0530.142<0.050.4260.1620.3920.0600.207
1510.0990.2130.218<0.050.4010.7060.106<0.050.143
2010.5810.2150.1520.1680.9970.6170.104<0.050.832
2510.8990.1480.1150.2800.6450.6880.3880.3020.678
3010.3490.2760.0660.3550.4660.1730.8010.6350.537
3510.3980.2850.2200.1810.6140.2710.6600.7580.391
401<0.050.2480.1530.3580.499<0.050.5290.5650.062
4510.0540.2430.2840.2270.3480.0640.3580.8670.125
501<0.050.4270.3790.5620.332<0.050.1140.127<0.05
601<0.050.5480.7940.8540.387<0.05<0.05<0.05<0.05
701<0.050.2470.8320.8330.281<0.05<0.05<0.05<0.05
801<0.050.4130.6070.5650.231<0.05<0.05<0.05<0.05
901<0.050.3520.7290.4410.210<0.05<0.05<0.05<0.05
1001<0.050.6870.3430.3830.103<0.05<0.05<0.05<0.05

References

  1. Amin, R.; Galen, B.F.; John, W.H.; Joseph, R.T.; James, D.P.; Christine, K.L. Rapidly pulsed reductants for diesel NOx reduction with lean NOx traps: Effects of pulsing parameters on performance. Appl. Catal. B Environ. 2018, 223, 177–191. [Google Scholar]
  2. Perng, C.C.; Easterling, V.G.; Harold, M.P. Fast lean-rich cycling for enhanced NOx conversion on storage and reduction catalysts. Catal. Today 2014, 231, 125–134. [Google Scholar] [CrossRef]
  3. Lietti, L.; Artioli, N.; Righini, L.; Castoldi, L.; Forzatti, P. Pathways for N2 and N2O formation during the reduction of NOx over Pt–Ba/Al2O3 LNT catalysts investigated by labeling isotopic experiments. Ind. Eng. Chem. Res. 2012, 51, 7597–7605. [Google Scholar] [CrossRef]
  4. Roy, S.; Gord, J.R.; Patnaik, A.K. Recent advances incoherent anti-Stokes Raman scattering spectroscopy: Fundamental developments and applications in reacting flows. Prog. Energy Combust. Sci. 2010, 36, 280–306. [Google Scholar] [CrossRef]
  5. Steuwe, C.; Kaminski, C.F.; Baumberg, J.J.; Mahajan, S. Surface enhanced coherent anti-Stokes Raman scattering on nanostructured gold surfaces. Nano Lett. 2011, 11, 5339–5343. [Google Scholar] [CrossRef] [PubMed]
  6. Williams, B.; Edwards, M.; Stone, R.; Williams, J.; Ewart, P. High precision in-cylinder gas thermometry using laser induced gratings: Quantitative measurement of evaporative cooling with gasoline/alcohol blends in a GDI optical engine. Combust. Flame 2014, 161, 270–279. [Google Scholar] [CrossRef]
  7. Hoffman, D.; Münch, K.U.; Leipertz, A. Two-dimensional temperature determination in sooting flames by filtered Rayleigh scattering. Opt. Lett. 1996, 21, 525–527. [Google Scholar] [CrossRef]
  8. Itani, L.; Bruneaux, G.; Lella, A.; Schulz, C. Two-tracer LIF imaging of preferential evaporation of multi-component gasoline fuel sprays under engine conditions. Proc. Combust. Inst. 2015, 35, 2915–2922. [Google Scholar] [CrossRef]
  9. Cho, K.Y.; Satija, A.; Pourpoint, T.L.; Son, S.F.; Lucht, R.P. High-repetition-rate three-dimensional OH imaging using scanned planar laser-induced fluorescence system for multiphase combustion. Appl. Opt. 2014, 53, 316–326. [Google Scholar] [CrossRef]
  10. Deguchi, Y.; Yasui, D.; Adachi, A. Development of 2D temperature and concentration measurement method using tunable diode laser absorption spectroscopy. J. Mech. Eng. Autom. 2012, 2, 543–549. [Google Scholar]
  11. Ma, L.; Li, X.; Sanders, S.T.; Caswell, A.W.; Roy, S.; Plemmons, D.H.; Gord, J.R. 50-kHz-rate 2D imaging of temperature and H2O concentration at the exhaust plane of a J85 engine using hyperspectral tomography. Opt. Express 2013, 21, 1152–1162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Busa, K.M.; Ellison, E.N.; McGovern, B.J.; McDaniel, J.C.; Diskin, G.S.; DePiro, M.J.; Capriotti, D.P.; Gaffney, R.L. Measurements on NASA Langley Durable Combustor Rig by TDLAT: Preliminary Results. In Proceedings of the 51st AIAA Aerospace Sciences Meeting Including New Horizons Forum and Aerospace Exposition, Grapevine, TX, USA, 7–10 January 2013. [Google Scholar]
  13. Busa, K.M.; Rice, B.E.; McDaniel, J.C.; Goyne, C.P.; Rockwell, R.D.; Fulton, J.A.; Edwards, J.R.; Diskin, G.S. Scramjet combustion efficiency measurement via tomographic absorption spectroscopy and particle image velocimetry. AIAA J. 2016, 54, 2463–2471. [Google Scholar] [CrossRef]
  14. Sun, P.; Zhang, Z.; Li, Z.; Gou, Q.; Dong, F. Study of two dimensional tomography reconstruction of temperature and gas concentration in combustion field using TDLAS. Appl. Sci. 2017, 7, 990. [Google Scholar] [CrossRef] [Green Version]
  15. Zhang, Z.R.; Pang, T.; Yang, Y.; Xia, H.; Cui, X.; Sun, P.; Wu, B.; Wang, Y.; Sigrist, M.W.; Dong, F. Development of a tunable diode laser absorption sensor for online monitoring of industrial gas total emissions based on optical scintillation cross-correlation technique. Opt. Express 2016, 24, 943–955. [Google Scholar] [CrossRef] [PubMed]
  16. Dickheuer, S.; Marchuk, O.; Tsankov, T.V.; Luggenhölscher, D.; Czarnetzki, W.; Ertmer, S.; Kreter, A. Measurement of the magnetic field in a linear magnetized plasma by tunable diode laser absorption spectroscopy. Atoms 2019, 7, 48. [Google Scholar] [CrossRef] [Green Version]
  17. Stritzke, F.; Kely, S.V.D.; Feiling, A.; Dreizler, A.; Wagner, S. TDLAS-based NH3 mole fraction measurement for exhaust diagnostics during selective catalytic reduction using a fiber-coupled 2.2-µm DFB diode laser. Opt. Express 2015, 119, 143–152. [Google Scholar] [CrossRef] [Green Version]
  18. Bolshov, M.A.; Kuritsyn, Y.A.; Romanovskii, Y.V. Tunable diode laser spectroscopy as a technique for combustion diagnostics. Spectrochim. Acta Part B 2015, 106, 45–66. [Google Scholar] [CrossRef]
  19. Xu, L.; Liu, C.; Jing, W.; Cao, Z.; Xue, X.; Lin, Y. Tunable diode laser absorption spectroscopy-based tomography system for on-line monitoring of two-dimensional distributions of temperature and H2O mole fraction. Rev. Sci. Instrum. 2016, 87, 013101. [Google Scholar] [CrossRef] [Green Version]
  20. Gamache, R.R.; Kennedy, S.; Hawkins, R.; Rothman, L.S. Total internal partition sums for molecules in the terrestrial atmosphere. J. Mol. Struct. 2000, 517–518, 407–425. [Google Scholar] [CrossRef]
  21. Rothman, L.S.; Gordon, I.E.; Barbe, A.; Benner, D.C.; Bernath, P.F.; Birk, M.; Bizzocchi, L.; Boudon, V.; Brown, L.R.; Campargue, A.; et al. The HITRAN2012 molecular spectroscopic database. J. Quant. Spectrosc. Radiat. Transf. 2013, 130, 4–50. [Google Scholar] [CrossRef]
  22. Hong, D.H.; Wang, L.; Truong, T.K. Low-complexity direct computation algorithm for cubic-spline interpolation scheme. J. Vis. Commun. Image Represent. 2018, 50, 159–166. [Google Scholar] [CrossRef]
  23. Deguchi, Y.; Noda, M.; Abe, M.; Abe, M. Improvement of combustion control through real-time measurement of O2 and CO concentrations in incinerators using diode laser absorption spectroscopy. Proc. Combust. Inst. 2002, 29, 147–153. [Google Scholar] [CrossRef]
  24. Choi, D.W.; Jeon, M.G.; Cho, G.R.; Kamimoto, T.; Deguchi, Y.; Doh, D.H. Performance improvements in temperature reconstructions of 2-D tunable diode laser absorption spectroscopy (TDLAS). J. Therm. Sci. 2016, 25, 84–89. [Google Scholar] [CrossRef]
  25. Zaatar, Y.; Bechara, J.; Khoury, A.; Zaouk, D.; Charles, J.-P. Diode laser sensor for process control and environmental monitoring. Appl. Energy 2000, 65, 107–113. [Google Scholar] [CrossRef]
  26. Doh, D.H.; Lee, C.J.; Cho, G.R.; Moon, K.R. Performances of Volume-PTV and Tomo-PIV. Open J. Fluid Dyn. 2012, 2, 368–374. [Google Scholar] [CrossRef] [Green Version]
  27. Choi, K.M.; Lee, J.T.; Jeon, S.G.; Hong, T.G.; Paek, J.K.; Han, S.H.; Yim, D.S. A review of fundamentals of statistical concepts in clinical trials. J. Korean Soc. Clin. Pharmacol. Ther. 2012, 20, 109–124. [Google Scholar] [CrossRef] [Green Version]
  28. Cort, J.W.; Kenhi, M. Advantages of the mean absolute error (MAE) over the root mean square error (RMSE) in assessing average model performance. Clim. Res. 2015, 30, 79–82. [Google Scholar]
  29. Seo, M.K.; Yun, W.Y. Clustering based hot strip roughing mill diagnosis using Mahalanobis distance. J. Korean Inst. Ind. Eng. 2017, 43, 298–307. [Google Scholar]
  30. Wikipedia. Available online: https://en.wikipedia.org/wiki/Mahalanobis_distance (accessed on 23 September 2019).
  31. Yu, J.W.; Jang, J.Y.; Yoo, J.Y.; Kim, S.S. Fault detection method for steam boiler tube using Mahalanobis distance. J. Korean Inst. Intell. Syst. 2016, 26, 246–252. [Google Scholar] [CrossRef] [Green Version]
  32. Minitab 19 Support. Available online: https://support.minitab.com/en-us/minitab/18/Assistant_Two_Sample_t.pdf (accessed on 20 November 2019).
Figure 1. Theoretical H2O absorption spectrum.
Figure 1. Theoretical H2O absorption spectrum.
Applsci 10 00285 g001
Figure 2. Schematic diagram of the experimental setup using tunable diode laser absorption spectroscopy (TDLAS).
Figure 2. Schematic diagram of the experimental setup using tunable diode laser absorption spectroscopy (TDLAS).
Applsci 10 00285 g002
Figure 3. Laser paths and measurement cell.
Figure 3. Laser paths and measurement cell.
Applsci 10 00285 g003
Figure 4. Results of mean transmission signals and mean absorption spectra of 10 laser beams: (a) 10 mean transmitted signals using 201 cycles and (b) 10 mean absorption spectra using 201 cycles.
Figure 4. Results of mean transmission signals and mean absorption spectra of 10 laser beams: (a) 10 mean transmitted signals using 201 cycles and (b) 10 mean absorption spectra using 201 cycles.
Applsci 10 00285 g004
Figure 5. Flowchart for tracking cycle and reconstructed temperature using TDLAT.
Figure 5. Flowchart for tracking cycle and reconstructed temperature using TDLAT.
Applsci 10 00285 g005
Figure 6. The spectra of reconstructed and measured absorption data at the engine running time 5.0 s: (a) 51 cycles, boundary; (b) 201 cycles, boundary; (c) 51 cycles, center; and (d) 201 cycles, center.
Figure 6. The spectra of reconstructed and measured absorption data at the engine running time 5.0 s: (a) 51 cycles, boundary; (b) 201 cycles, boundary; (c) 51 cycles, center; and (d) 201 cycles, center.
Applsci 10 00285 g006
Figure 7. 2D reconstructed temperature distributions at the engine running time 5.0 s: (a) 51 cycles and (b) 201 cycles.
Figure 7. 2D reconstructed temperature distributions at the engine running time 5.0 s: (a) 51 cycles and (b) 201 cycles.
Applsci 10 00285 g007
Figure 8. Procedure for mean different test and similarity analysis.
Figure 8. Procedure for mean different test and similarity analysis.
Applsci 10 00285 g008
Figure 9. Comparison of TThc and TRec for 71 cycles.
Figure 9. Comparison of TThc and TRec for 71 cycles.
Applsci 10 00285 g009
Figure 10. Comparison of TThc and TRec for 601 cycles.
Figure 10. Comparison of TThc and TRec for 601 cycles.
Applsci 10 00285 g010
Figure 11. Comparison of TThc and TRec for 201 cycles.
Figure 11. Comparison of TThc and TRec for 201 cycles.
Applsci 10 00285 g011
Figure 12. 2D reconstruction temperature distributions at 201 cycles: (a) 1.0 s; (b) 2.0 s; (c) 3.0 s; (d) 4.0 s; (e) 5.0 s; (f) 6.0 s; (g) 7.0 s; (h) 8.0 s; and (i) 9.0 s.
Figure 12. 2D reconstruction temperature distributions at 201 cycles: (a) 1.0 s; (b) 2.0 s; (c) 3.0 s; (d) 4.0 s; (e) 5.0 s; (f) 6.0 s; (g) 7.0 s; (h) 8.0 s; and (i) 9.0 s.
Applsci 10 00285 g012
Table 1. TRec per 18 candidate cycles used in [Step 1].
Table 1. TRec per 18 candidate cycles used in [Step 1].
No.Time
(s)
TRec
3151 9011001
11.0404.833400.000 405.432404.189
2406.775401.665 407.081405.842
3404.947399.784 405.436404.455
4407.458401.846 407.028405.633
5409.400403.506 408.581407.223
6407.556401.614 406.820405.720
7404.951399.982 405.459404.137
8406.842401.604 406.900405.606
9404.999399.722 405.051405.159
8069.9398.712399.182 400.164400.543
807397.870397.965 402.270402.782
808396.147396.544 401.184401.859
809397.183397.668 401.351401.646
810396.221396.327 403.357403.785
Table 2. Results of paired t-test.
Table 2. Results of paired t-test.
Time (s)Test Setsp-Value
1.0All sets<0.05
2.0251vs.3010.626
451vs.5010.476
451vs.6010.402
501vs.6010.331
3.0All sets<0.05
4.0351vs.4010.371
5.0451vs.5010.208
6.0401vs.4510.633
7.0601vs.7010.221
8.0All sets<0.05
9.0All sets<0.05
Table 3. Temperatures used in [Step 2].
Table 3. Temperatures used in [Step 2].
No.GroupTime
(s)
TThc(i) T ¯
3151 1001
111.0408.000406.418401.080 405.204
21.1408.688402.722401.324 406.271
31.2409.921407.890407.890 406.525
41.3407.798410.285411.842 406.990
51.4411.453402.434411.723 406.453
61.5406.759392.685401.679 406.284
71.6406.005397.013405.508 407.419
81.7405.461418.457416.981 407.537
91.8405.110407.200401.237 406.350
101.9404.842416.607407.599 405.416
8099.0404.658397.217397.462 401.795
819.1405.853394.707393.932 401.638
829.2402.657408.041403.190 400.741
839.3400.688395.484396.351 400.231
849.4401.921402.101398.922 401.280
859.5399.798393.647396.263 401.277
869.6403.453410.205399.067 401.559
879.7400.895390.764398.329 401.343
899.8403.987398.441401.295 400.894
909.9401.752403.739409.931 400.857
Table 4. Results of Two-sample t-test.
Table 4. Results of Two-sample t-test.
CyclesTime (s)
123456789
p-Value
310.6470.4630.8390.2730.6920.4960.9240.9760.169
510.7110.7260.3160.2920.7070.6810.8640.7280.069
710.7740.9850.1470.2920.9820.8960.7270.9970.164
910.6030.8910.2250.4410.3690.9850.8420.6930.362
1110.5310.9830.3900.6670.0920.9010.8280.3520.482
1510.4860.8920.4020.671<0.050.7790.7080.5740.474
2010.3960.9640.3280.381<0.050.9580.6840.4880.177
2510.5100.8580.2710.2890.0740.7430.7910.4670.078
3010.3840.8900.2460.325<0.050.7890.9280.1370.109
3510.4150.9280.1960.411<0.050.7700.9790.1020.073
4010.4060.7900.1770.253<0.050.6550.9840.1120.051
4510.5300.7350.1360.235<0.050.5370.9460.081<0.05
5010.5090.6910.1360.269<0.050.6220.865<0.05<0.05
6010.4780.6200.0970.284<0.050.4950.886<0.05<0.05
7010.4320.6310.0840.367<0.050.4490.911<0.050.054
8010.3570.5500.0580.336<0.050.4300.977<0.05<0.05
9010.3160.5740.0530.4310.0620.4310.871<0.05<0.05
10010.2160.528<0.050.4460.0570.4130.773<0.05<0.05
Table 5. Results of root mean square error (RMSE) analysis between TThc and TRec for 18 candidate cycles.
Table 5. Results of root mean square error (RMSE) analysis between TThc and TRec for 18 candidate cycles.
CyclesTime (s)Total
123456789
RMSERMSERank
318.4878.03910.3258.44112.6227.8697.3205.8875.8878.57818
515.5145.5926.3705.9947.7594.7766.9125.5265.5266.04917
716.2694.4975.4046.3345.3166.0076.0284.4204.4205.53816
914.8793.7754.3315.5593.8945.9215.0733.6203.6204.82215
1114.7023.8604.2825.0463.5805.2654.1563.4963.4964.39314
1513.7553.1294.2164.8333.5642.9574.3144.0654.0653.91813
2013.1772.9783.3493.3762.9002.5293.8103.7363.7363.19912
2513.0243.0633.4963.4002.3383.1663.2432.7552.7553.12011
3012.6032.8843.1592.9612.7492.8442.5412.2672.2672.80910
3512.0562.7332.7992.6862.9553.3332.5292.0992.0992.7039
4011.9512.5832.6412.4512.8493.1492.5462.4562.4562.5844
4512.0472.9092.6622.4032.7853.4392.5122.1562.1562.6417
5012.4332.5802.7032.2782.8943.5432.4332.1402.1402.6438
6012.0962.4532.5612.2652.9473.3682.1822.0392.0392.4911
7012.2012.5912.8992.4262.9253.5441.9051.9871.9872.5553
8012.2092.3183.0552.6292.7713.4921.9071.8701.8702.5492
9012.4382.4993.1742.5682.7353.6932.0541.7151.7152.6246
10012.4342.2703.4382.5762.6933.6062.0411.6601.6602.5975
Table 6. Results of time-based Mahalanobis distance (TMD) analysis between TThc and TRec for 18 candidate cycles.
Table 6. Results of time-based Mahalanobis distance (TMD) analysis between TThc and TRec for 18 candidate cycles.
CyclesTime (s)Total
123456789
TMDTMDRank
311.8601.1241.1881.2081.1741.3830.8890.7310.8801.1214
511.7581.0791.2791.2671.2041.2530.9100.8540.8761.1366
711.8290.9671.2501.3050.9951.3630.9700.7361.0141.1235
911.7800.8731.1701.3481.0061.2990.9150.7731.1481.1123
1111.7770.7701.1551.3071.0221.3040.9210.8350.9621.0822
1511.7270.9361.1331.4561.1431.3031.0520.9971.0711.18110
2012.1320.7981.5520.8881.5490.9920.6840.4290.5060.9331
2511.6751.2771.3911.2730.9651.3611.0160.8121.1631.19014
3011.7191.2711.5061.2521.2581.2911.0500.7550.8721.18512
3511.6261.2801.3941.2611.2911.4810.9960.7450.9111.18813
4011.7081.2971.4771.2731.3631.2631.0460.9110.9931.23718
4511.7721.3751.3121.2921.2051.2580.9720.7881.0331.19416
5011.8711.3421.2641.3281.3091.1530.9450.8111.0561.20017
6011.8801.3731.2621.4391.3200.9810.9970.7261.0981.19115
7011.9561.3551.3411.4581.2580.9900.8820.7541.0531.18211
8012.0121.3201.3471.4791.2990.8490.8140.7271.1471.1658
9012.0611.3191.4111.4001.2461.0550.7790.6931.1381.1779
10012.0581.2331.3821.3931.2431.0060.7880.6551.1531.1547

Share and Cite

MDPI and ACS Style

Jang, H.; Choi, D. Similarity Analysis for Time Series-Based 2D Temperature Measurement of Engine Exhaust Gas in TDLAT. Appl. Sci. 2020, 10, 285. https://doi.org/10.3390/app10010285

AMA Style

Jang H, Choi D. Similarity Analysis for Time Series-Based 2D Temperature Measurement of Engine Exhaust Gas in TDLAT. Applied Sciences. 2020; 10(1):285. https://doi.org/10.3390/app10010285

Chicago/Turabian Style

Jang, Hyeonae, and Doowon Choi. 2020. "Similarity Analysis for Time Series-Based 2D Temperature Measurement of Engine Exhaust Gas in TDLAT" Applied Sciences 10, no. 1: 285. https://doi.org/10.3390/app10010285

APA Style

Jang, H., & Choi, D. (2020). Similarity Analysis for Time Series-Based 2D Temperature Measurement of Engine Exhaust Gas in TDLAT. Applied Sciences, 10(1), 285. https://doi.org/10.3390/app10010285

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