Journal of Combustion

Volume 2016 (2016), Article ID 8593523, 9 pages

http://dx.doi.org/10.1155/2016/8593523

## Review of Sensing Methodologies for Estimation of Combustion Metrics

Michigan Technological University, 1400 Townsend Drive, Houghton, MI 49931, USA

Received 3 July 2015; Accepted 10 January 2016

Academic Editor: Sergey M. Frolov

Copyright © 2016 Libin Jia et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

For reduction of engine-out emissions and improvement of fuel economy, closed-loop control of the combustion process has been explored and documented by many researchers. In the closed-loop control, the engine control parameters are optimized according to the estimated instantaneous combustion metrics provided by the combustion sensing process. Combustion sensing process is primarily composed of two aspects: combustion response signal acquisition and response signal processing. As a number of different signals have been employed as the response signal and the signal processing techniques can be different, this paper did a review work concerning the two aspects: combustion response signals and signal processing techniques. In-cylinder pressure signal was not investigated as one of the response signals in this paper since it has been studied and documented in many publications and also due to its high cost and inconvenience in the application.

#### 1. Introduction

Determination of combustion metrics for an internal combustion engine has the potential of providing feedback for closed-loop combustion phasing control to meet current and upcoming emission and fuel consumption regulations. Closed-loop control of the combustion process has been a focus for engine research and development [1–3]. Open-loop operation based on calibration maps which are conservatively set based upon laboratory operation can give a quick response and is relatively easy to control. However, the open loop cannot adapt to the changes caused by the condition variations such as injector aging and fuel quality [2].

In comparison, closed-loop control considers the condition changes in the control mechanism and enables operation closer to the optimum fuel consumption and emissions target. Combustion metrics detected or estimated through a sensor, which is referred to as combustion sensing output, provides feedback information to control the combustion process. In this paper, the combustion metrics being addressed involves SOC (start of combustion), peak pressure location, peak apparent heat release rate location, CA50 (crank-angle location for 50% fuel burnt), and so forth. The feedback information provided by the combustion sensing process also refers to the abnormal combustion phenomena including knock and misfire events. In this paper, a review was performed for combustion sensing methodologies concerning combustion metrics estimation with respect to combustion sensing response signals and signal processing techniques.

#### 2. Response Signal for Combustion Event

In-cylinder pressure waveform is the most commonly used signal which provides the information for engine combustion control [4–6]. In-cylinder pressure signal has been historically used in the laboratory and more recently in series production to derive the combustion metrics and provide feedback for combustion phasing control [7–9]. However, the measurement of the in-cylinder pressure is typically obtained with intrusive sensors that require a special mounting process and engine structure modification. Also the in-cylinder pressure transducer has a high cost for mass production for diesel engines [10]. So the response signals discussed in this paper do not include the in-cylinder pressure signal.

##### 2.1. Crank-Shaft Speed Fluctuation

Due to the variations of the in-cylinder pressure waveform during a combustion cycle, the crank-shaft speed fluctuation varies in a complex way which depends on the engine parameters. How the speed fluctuation varies with the engine in-cylinder pressure changes has been explored by many researchers so as to develop a good alternative to the direct intrusive in-cylinder pressure measurement [11–19].

Based on a model that relates the crank-shaft speed and the in-cylinder pressure, the in-cylinder pressure can be estimated with the input of instantaneous speed signal measured by a crank-shaft speed sensor. The sensor can be an optical encoder or a magnetic pickup transducer which are easy to mount and low in cost.

Moro et al. [11] proposed a linear dependency between the in-cylinder pressure and the engine speed signal and experimentally verified it for 38 different engine running conditions. The equation that can represent this linear dependency is given as [11] where is piston area in , is the crank radius in , is the in-cylinder pressure in case of misfire (bar), is the synthetic engine acceleration in , is the synthetic engine acceleration in case of misfire , is the moment of initial (), and is the crank-slider kinematics function.

A frequency response function between the in-cylinder pressure and the engine speed can be obtained by converting (1) into frequency domain. However, this FRF is sensitive to engine running conditions and the FRF obtained based on one condition does not lead to the estimated in-cylinder pressure with high accuracy when condition varies. So a FRF mapping was created in this paper based on 38 different steady-state engine conditions with the engine speed and the manifold pressure as the condition parameters to distinguish different test conditions. For the conditions falling into the FRF mapping, interpolation technique was used for both real and imaginary harmonic components to obtain the estimated FRF. The pressure recovery results for low speed low load, high speed low load, and low speed high load conditions were shown in this paper.

Connolly and Yagle modeled the cylinder combustion pressure via the crank-shaft velocity from a statistical point of view [12]. The model involves three sequent components. First, by replacing the time domain independent variables with crank-angle variables, a nonlinear differential model between the crank’s shaft speed and the in-cylinder pressure signal can be simplified. Secondly, the in-cylinder pressure signal was parameterized by the sample modeling sequence based on a stochastic model which uses the sum of the deterministic waveform and an amplitude-modulated cosine window. Third, an estimation of the in-cylinder pressure based on the crank-shaft speed signal was achieved through a state-space deconvolution process which utilized a Kalman filter. Moreover, signal to noise ratio effects to the in-cylinder pressure estimation were also evaluated in this paper. Results showed that for low to moderate noise level conditions the reasonable deconvolution can be reached.

Shiao and Moskwa [13] employed a sliding observer to estimate the in-cylinder pressure and combustion heat release for an SI engine. To estimate the in-cylinder pressure with high accuracy, the error between the measured and the estimated crank-shaft speed was taken as the feedback to reduce the dynamic error of the estimated in-cylinder pressure. The unobservability problem arises for the pressure estimation around the top-dead center and thus introduces significant estimation error. This problem was partly solved by adapting the parameters of the observer. Then the estimated in-cylinder pressure was used to compute the cylinder heat release. Also, detection of misfire or abnormal combustion events was achieved through the estimated heat release.

Additional investigators examining combustion metrics analysis based on the crank-shaft speed fluctuation can be found in [14–16]. In addition to estimating the in-cylinder pressure waveform and the heat release, engine crank-shaft speed was also used to recover the engine torque [17–19]. For most cases, the crank-shaft speed was fed into an engine model which was simplified based on assumptions to estimate the engine torque.

##### 2.2. In-Cylinder Ion Current

Ion current in the combustion chamber is measured via the spark plug. After the high-voltage discharge, the ion current across the spark plug gap is obtained by applying a DC voltage across the gap and measuring the resulting current. The ion current is affected by gas flow, geometry of flame, electric potential, ion density, and the angle between the flame and electrode [20]. The ionization of gases in the cylinder occurs in two phases. When the fuel reacts with the oxygen during combustion, the first phase ionization occurs which can be considered as chemical phase. The second phase, defined as thermal phase, occurs when the burnt gases are compressed by the increased in-cylinder pressure [21]. The most consistent dependency between the ion current and in-cylinder pressure occurs on the peaks of the two signals for both amplitude and the location perspectives. This has been verified by the researches in [20–25].

Martychenko et al. [22] detected the breakdown voltage across the spark plug gap and modeled the relationship between the peak of the voltage and the peak of the in-cylinder pressure based on second-order polynomial function. The coefficients of the second-order polynomial function for the conditions with varied engine speed are different. However, the coefficients can be curve-fitted by a linear function of engine speeds. Hellring and Holmberg [21] proposed least squares fit method to estimate the in-cylinder pressure peak position of spark ignited engines based on the ion current signal. This method was proved to have a better robustness and accuracy than multilayer perceptron and Gaussian curve fit methods for peak in-cylinder pressure estimation.

Gazis et al. [28] explored the possibility of estimating in-cylinder characteristics based on the ion current with one simple and computationally inexpensive neural network, adaptive linear type of network. Thirteen extracted characteristics of the ion current were taken as the input and four characteristics of the in-cylinder pressure (peak pressure position, peak pressure magnitude, the width of curve at half of its height, and the area of the curve between inlet valve closing (IVC) and exhaust valve opening (EVO)) as the output to train the network with the purpose of predicting the four characteristics of in-cylinder pressure. Also, based on the same neural network structure but with the whole ion current signal (time domain sampled between IVC and EVO) as the input and the whole in-cylinder pressure signal (time synced with ion current signal) as the output, the in-cylinder pressure curve rather than just some characteristics of the in-cylinder pressure signal can be estimated. The peak pressure location, as one of the most important in-cylinder pressure characteristics, was predicted with the mean error at 0.062 degrees and standard deviation at 2.55 degrees.

Ion current was also used to detect engine knock, misfire, or incomplete combustion [24, 29–31]. Kumar et al. [24] applied a band-pass filter on the ion current signal and the filtered output indicates the engine knock. Danne et al. [29] compared the ion current based knock detection with the conventional methods of pressure based and accelerometer-based knock detection on a large-displacement, air-cooled, V-twin motorcycle engine. It was found that the ion current based method can detect the inaudible knock more accurately and reject the mechanical noise more effectively than the other two conventional methods. Zhu et al. [30] found that the in-cylinder ion current can detect misfire or incomplete combustion. Also, the ion current signal can be used to compute minimum spark advance for best torque (MBT) to measure the combustion stability [30, 32]. However, the results are only limited to a fixed load over a narrow speed range.

##### 2.3. Accelerometer Signal

Accelerometers are mounted externally on the engine block or the engine head to detect the combustion events by measuring the vibration signals which are transmitted from the in-cylinder oscillation to the engine outer surface. However, as the accelerometer detects the vibrations from the sources in addition to the cylinder oscillation including the valve dynamics and piston slaps, the signal may vary from cylinder to cylinder and over operating conditions. So the utilization of the accelerometer signal for combustion metrics detection relies on the signal processing technique which can eliminate the effects from other sources.

Naber et al. [33] evaluated the effectiveness and accuracy of accelerometer-based knock detection. The distributions of the accelerometer-based knock intensity metrics for various operation conditions including varied speeds, loads, cam timings, and knock levels were measured and fitted by a log-norm distribution. The log-norm model was verified to provide a good fit of the distributions and the distribution characteristics including skewness and peakedness. In addition, a good correlation can be seen between the cylinder pressure based knock intensity metrics and the accelerometer-based knock intensity metrics. Guillemin et al. [34] estimated the instantaneous engine knock by fitting the accelerometer signals with Gaussian function on a 2.2 L HCCI engine and measure the start of combustion where the knock level is out of the user-defined threshold.

Characteristics of the accelerometer signal which are related to the characteristics of in-cylinder pressure or apparent heat release rate were investigated and extracted [35–37]. Some characteristics of the in-cylinder pressure or apparent heat release including start of combustion, CA50, and peak pressure crank-angle location are closely related to the combustion process and thus can be used as the feedback to control the combustion process. Arnone et al. [35] band-pass filtered the in-cylinder pressure signal and the accelerometer signal within 650–1000 Hz and found that the accelerometer signal can locate the sudden rise of the in-cylinder pressure signal (so as to denote the start of combustion), diffusive combustion process, and the peak of the in-cylinder pressure on a water cooled Lombardini LDW442CRS direct injection common rail diesel engine. Chiavola et al. [36] computed the cumulative heat release based on the measured in-cylinder pressure and investigated the relationship between the accelerometer signal and the cumulative heat release on a two-cylinder diesel engine equipped with a common rail injection system. By superimposing the filtered accelerometer signal to the cumulative heat release, it was found that the filtered accelerometer signal can locate the start of combustion, the beginning of main combustion, and MFB50 (50% of the burnt fuel mass). Taglialatela et al. [37] investigated the correlation between the in-cylinder pressure signal and the features derived from the accelerometer signal on a 4 L single cylinder SI engine. Time-frequency spectrogram method was utilized to analyze the accelerometer signal to present more features of the accelerometer signal than the analysis in time domain. The result indicates a direct correlation between the peak pressure location and the maximum amplitude of accelerometer signal in time-frequency domain for all the engine operating conditions conducted in this paper. So the maximum amplitude location of the accelerometer signal can be used as the feedback for a closed-loop control system of spark advance.

Polonowski et al. [38] explored the potential of accelerometers to recover the in-cylinder pressure curve on a 1.9 L four-cylinder, turbocharged, HPCR, direct injection diesel engine. In this paper, standard signal processing techniques including Fast Fourier Transform (FFT) and coherence were employed and results showed that a strong coherence presented between the in-cylinder pressure signal and the accelerometer signal within frequency band of 0.5 kHz to 4 kHz with the coherence value over 0.9. Also, this research found that the accelerometer location did result in a varied coherence value between the in-cylinder pressure signal and the accelerometer signal. The optimal locations for the accelerometer placement were determined based on both offline and online coherence analysis. In his later work [39], frequency response function (FRF) was used to quantify the relationship between accelerometer response and AHR and the relationship between accelerometer response and in-cylinder pressure. A technique termed as spectrum weighting was utilized to combine FRFs from all conducted test conditions into a single FRF by weighting the FRF magnitude and phase information for each frequency by the coherent output power at that specific frequency.

This obtained FRF was used to estimate in-cylinder pressure and AHR on a cycle-cycle basis. The maximum pressure gradient determination was shown to have a root mean square error (RMSE) accuracy of 15% of actual maximum pressure gradient. The location based metrics had the RMSE as small as 0.29° and more than 80% of the estimated peak apparent heat locations were within 1° crank-angle.

#### 3. Signal Processing Techniques to Correlate the Combustion Metrics and the Response Signal

##### 3.1. Frequency Response Function (FRF)

The pressure curve recovery depends on the transfer path modeling between the acquired source signal (in-cylinder pressure signal) and the response signal (e.g., vibration signal). Frequency response function which represents the frequency domain relationship between the in-cylinder pressure signal, , and the accelerometer signal, , can be presented by is the transfer path in frequency domain and is defined as the frequency response function (FRF). With the FRF obtained by measuring both and based on the representative operating conditions, the in-cylinder pressure signal of any other condition can be recovered with the measured and the obtained byHowever, as the response signal is sensitive to not only the in-cylinder pressure oscillation but the rotating crank-slider and vibration from other mechanical parts including piston slaps and valve dynamics which varies with the engine operating conditions, the FRF with assumption of linear dependency in frequency domain between the source signal and the response signal does not have a good robustness over engine operating conditions. This conclusion has been confirmed by researchers [39–42]. Gao and Randall [40] applied the FRF computed from 2400 rpm full load condition to reconstruct the in-cylinder pressure with the accelerometer signal measured at 3600 rpm and full load condition as the input. The significant recovery error for the in-cylinder pressure waveform indicated that the transfer path modeled by FRF cannot be considered consistent over engine operating conditions. Morello et al. [41] attempted to overcome the drawback of FRF application for the heat release recovery based on the accelerometer signal by optimizing the time domain window applied to the accelerometer signal. Also, a Vold-Kalman order tracking filter was employed to eliminate the abnormal harmonics of the singular value decomposition results of both the accelerometer signal and in-cylinder pressure signal. However, no significant heat release estimation accuracy improvement can be seen for start of injection (SOI) sweep test conditions. Polonowski [39] also tried to improve the FRF performance by adding a weighting function to the FRF. However, the robustness improvement for FRF is still limited.

Gao and Randall [40] explained why the variation of the FRF results in an error of source estimation with the Laplace transform. It can be seen that the variation of FRF introduces incomplete cancellation of the poles and zeros in the dominator with the accelerometer signal as the numerator. The incomplete cancellation of non-minimum-phase zeros will make the extra poles or zeros of FRF present and make the inverse filtering unstable. As a result errors will be introduced to the estimation result.

##### 3.2. Cepstrum Analysis

Complex cepstral analysis is a nonlinear homomorphic signal process which is being utilized in many areas including machine diagnostics, image processing, speech, and radar signal processing. A cepstrum is reached by taking the inverse Fourier transform of the logarithm of a signal spectrum. The complex cepstrum can be expressed as is the signal in time domain, represents the Fourier transform algorithm, and denotes inverse Fourier transform.

Equation (5) can be obtained by applying logarithm to (2): After computing the inverses Fourier transformAccording to (4), the cepstrum of FRF can be obtained as The advantage of this method is that the convolution process is converted to an addition process in cepstrum domain. El-Ghamry et al. [42] applied the complex cepstrum analysis on the root mean square acoustic emission signal. The complex cepstrum of FRF, , was evaluated for four complete combustion cycles at 1280 rpm 30 Nm condition. By inserting and into (7), the cylinder pressure signal can be identified asHowever, this method only gave good estimation of the in-cylinder pressure signal for the same engine operating conditions based on which was computed.

Another important application of complex cepstral analysis is for signal smoothing. The complex cepstrum was utilized to improve the robustness of the transfer path by smoothing both the source signal and the response signal [34, 40, 43]. Smoothing the FRF actually reduces the variations of the FRF associated with different operating conditions. The schematic illustration for smoothing the FRF is shown in Figure 1.