Next Article in Journal
SASRT: Semantic-Aware Super-Resolution Transmission for Adaptive Video Streaming over Wireless Multimedia Sensor Networks
Next Article in Special Issue
Spatio-Temporal Synchronization of Cross Section Based Sensors for High Precision Microscopic Traffic Data Reconstruction
Previous Article in Journal
Detailed Urban Land Use Land Cover Classification at the Metropolitan Scale Using a Three-Layer Classification Scheme
Previous Article in Special Issue
An Adaptive Track Segmentation Algorithm for a Railway Intrusion Detection System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Virtual In-Cylinder Pressure Sensor Based on EKF and Frequency-Amplitude-Modulation Fourier-Series Method

1
School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
2
School of Automotive Studies, Tongji University, Shanghai 201804, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(14), 3122; https://doi.org/10.3390/s19143122
Submission received: 29 May 2019 / Revised: 8 July 2019 / Accepted: 10 July 2019 / Published: 15 July 2019
(This article belongs to the Special Issue Intelligent Transportation Related Complex Systems and Sensors)

Abstract

:
As a crucial and critical factor in monitoring the internal state of an engine, cylinder pressure is mainly used to monitor the burning efficiency, to detect engine faults, and to compute engine dynamics. Although the intrusive type cylinder pressure sensor has been greatly improved, it has been criticized by researchers for high cost, low reliability and short life due to severe working environments. Therefore, aimed at low-cost, real-time, non-invasive, and high-accuracy, this paper presents the cylinder pressure identification method also called a virtual cylinder pressure sensor, involving Frequency-Amplitude Modulated Fourier Series (FAMFS) and Extended-Kalman-Filter-optimized (EKF) engine model. This paper establishes an iterative speed model based on burning theory and Law of energy Conservation. Efficiency coefficient is used to represent operating state of engine from fuel to motion. The iterative speed model associated with the throttle opening value and the crankshaft load. The EKF is used to estimate the optimal output of this iteration model. The optimal output of the speed iteration model is utilized to separately compute the frequency and amplitude of the cylinder pressure cycle-to-cycle. A standard engine’s working cycle, identified by the 24th order Fourier series, is determined. Using frequency and amplitude obtained from the iteration model to modulate the Fourier series yields a complete pressure model. A commercial engine (EA211) provided by the China FAW Group corporate R&D center is used to verify the method. Test results show that this novel method possesses high accuracy and real-time capability, with an error percentage for speed below 9.6% and the cumulative error percentage of cylinder pressure less than 1.8% when A/F Ratio coefficient is setup at 0.85. Error percentage for speed below 1.7% and the cumulative error percentage of cylinder pressure no more than 1.4% when A/F Ratio coefficient is setup at 0.95. Thus, the novel method’s accuracy and feasibility are verified.

1. Introduction

The development of cleaner and more efficient engines requires the continuous measurement of in-cylinder pressure which is of fundamental importance for combustion optimization, air-fuel ratio control, noise and pollutant reduction. More specifically, cylinder pressure is a fundamental intermediate variable that indicates engine state and drives models for engine combustion control. High-bandwidth control of ignition and the air-fuel ratio may permit these engines to operate reliably near the lean-burn limit, with significant improvements in efficiency and decreased emission of CO and NOx. Fully utilizing this technology requires a key factor: real-time cylinder pressure [1]. The goal of Rizzoni’s [2] method is to build a stochastic model for the combustion pressure process in the spark-ignition engine. Control of torque balance can improve drivability and can suppress noise from light-duty diesel engines, but as a crucial part of the control algorithm, obtaining pressure data necessitates expensive pressure sensors and demands considerable computational time [3]. The advent of cylinder-pressure transducers seems promising, boasting a better approach to detect the in-cylinder pressure; however, its reliability, complex installation, high cost, and short working lifespan are full of controversy. As new technologies improve, such as computational methods, reconstructing cylinder pressure in a multi-cylinder internal combustion (IC) engines becomes more feasible. Bennett [4] proposed a robust adaptive-gradient descent-trained NARX neural network using both crank velocity and crank acceleration as inputs to reconstruct cylinder pressure in multi-cylinder IC engines. Another work has shown comprehensive results from an experimental assessment of a common rail diesel engine operated with neat diesel fuel and with ethanol substitutions [5]. With different substitutions and different substitution rates, pressure varies rapidly. Researchers have pursued other various approaches to develop a comparatively optimal way to obtain the desired cylinder pressure, directly or indirectly. Rizvi [6] proposed a novel method to detect engine misfire faults based on a hybrid model of the gasoline engine. As this method mainly employs cylinder pressure, closed-loop [7] combustion control [8] becomes feasible and misfire can be detected in multiple ways. However, these approaches vary in cost, reliability, robustness, accuracy, and convenience. Therefore, a low-cost [9] noninvasive soft-pressure sensor with a high level of reliability and accuracy is necessary. Payri [10] presented a step-by-step approach to optimize the signal processing both for offline and online applications based on the characteristics of the signal. Maurya [11] proposed that a method based on standard deviation of pressure, and pressure rise rate is used to find the minimum number of engine cycles to be recorded for averaging to get reasonably accurate pressure data independent of cyclic variability. Bhatti [12] proposed that two robust second-order sliding mode observers are employed that require two state mean value engine models based on inlet manifold pressure and rotational speed dynamics. The design of the second-order sliding mode observer depends largely on the model accuracy. And, more importantly, this method had the problems of chattering and large initial error. As a result, the tracking error of the steady-state speed was only about 4% and required a certain convergence time. Ferrari [13] proposed that pressure of nozzle opening and closing was identified by means of pressure data fitted by bench test. Then the time-frequency analysis was used to obtain the mean instantaneous frequency and adjust the control strategy of the injector. However, this method did not process data effectively, such as filter, before fast Fourier time-frequency transform, so this may affect the accuracy. Eriksson [14] calculated the corresponding relationship between pressure value and crankshaft rotation angle. The error of peak identification based on ion current in (Eriksson, 2003, Figure 4) was still large. Comparing with the corresponding relationship proposed by Eriksson [14], the cylinder pressure time curve is presented in this paper. It is based on burning theory and Law of energy conservation, and then modified by Extended-Kalman-Filter-optimized (EKF) engine model, which has a higher stability.
Observer-based method [12,13,14,15] is proven to be a good way to predict the output of SISO system. Engine cylinder pressure signal can be used to characterize the engine combustion state [16,17]. Due to the influence of multiple variables, it is difficult to establish a precise physical model of cylinder pressure. However, by observing and analyzing the cylinder pressure signal, this signal is cyclical with variable frequency (associated with speed) and variable amplitude (associated with the peak combustion pressure). Meanwhile, the engine cylinder pressure signal has a delay characteristic, so it is necessary to predict the current cylinder pressure value by the input of the current moment.
In this paper, a novel method which can identify engine pressure data accurately will provide a spacious viewpoint and means to process signals, thereby achieving important information for practical engineering. This novel method involving Frequency-Amplitude Modulated Fourier Series (FAMFS) and Extended-Kalman-Filter-optimized (EKF) engine model for identifying the periodic signal with its variable frequency and amplitude is proposed. This cylinder pressure identification method is also called a virtual cylinder pressure sensor.
Figure 1 shows a flow chart of method. Crankshaft speed [18,19,20] plays an important role in identifying the pressure as well as other applications, especially in misfire detection [21]. Instantaneous speed fluctuations [22], crankshaft segment acceleration, and transient rotational speed are the most widely utilized variables to locate misfire. Crank speed can also reveal and identify real-time engine combustion parameters [20], which is adopted in this paper. Speed data is affected by measurement noise from the crankshaft and process noise from the engine. These noise sources result in a huge accumulation error during pressure measurements, and EKF is chosen to optimize prediction and filtering [23,24].
The crankshaft speed in one solution step is related to the speed of the last step, so the speed model is an iterative process. Firstly, a novel speed iteration model based on burning theory and Law of energy Conservation is proposed. This iteration model is associated with the throttle opening value and the load. The EKF is then used to estimate the optimal output of the speed iteration model, and this optimal output can then be used for computing the frequency and amplitude of the cylinder pressure cycle-to-cycle. Secondly, taking Otto cycle of standard four-stroke engine as an example, a 24th order Fourier series is chosen to fit the standard working cycle. Thirdly, associated with the variable frequency and amplitude, a frequency-amplitude modulation method is adopted to modulate the standard pressure model identified in step two. Finally, system performance is evaluated by an actual working engine. Therefore, the main contributions of this paper can be summarized as follows:
(1)
Aimed to identify in-cylinder pressure and its periodic signals with variable frequency and amplitude, a virtual pressure sensor for engine healthy monitoring with low-cost, real-time, non-intrusive, a long-lifespan, and high-reliability in severe working conditions is presented.
(2)
A novel speed model is established according to burning theory, Law of energy Conservation and EKF. The FAMFS method is next developed to fit the periodic signals with variable frequency and variable amplitude.
(3)
The proposed method can be applied to multi-cylinder internal combustion (IC) engines including four-cylinder engines.
In this paper, on the one hand, as the key parameter of cylinder pressure identification, the speed iterative model using the EKF has been predicted with high accuracy. However, the overall deviation came from the calculation of pressure due to the unpredictable process of power stroke. On the other hand, due to only related to the crankshaft speed and the amount of air and fuel injection in engine, this method can be applied to multi-cylinder internal combustion (IC) engines including four-cylinder engines. More specifically, many key parameters of ignition engine such as spark advance angle (SAA) need to be considered in the process of building virtual in-cylinder pressure sensor model. Accordingly, at present, the application of the virtual sensor is only limited to spark ignition engines.

2. Crankshaft Speed Model

2.1. Design of Speed Iteration Model

As generally known, the angular acceleration of crankshaft in one cycle is produced by the output power of the engine, and the output energy per cycle is determined by the amount of air and fuel injected into engine. In internal combustion engines, the amount of air and fuel trapped in the cylinder is relatively influenced by engine speed through the volumetric efficiency, and mainly depends on the engine load. The change of engine load will affect crankshaft speed; hence the amount of air and fuel is deemed relating to the speed. Therefore, it forms a closed-loop iteration process. The crankshaft speed model (CSM) is established in two parts: constant load (Part A) and variable load (Part B). In Part A, the amount of air and fuel is related to the speed of the previous cycle. The power produced through fuel burning compels the crankshaft to rotate; and in Part B, the engine overcomes the load with the price of crankshaft speed reduction. The complete model is then forged by the combination of these two parts.

2.1.1. Constant Load Condition

In the following process, the load is constant. Cylinder pressure is directly related to the amount of fuel and air injection in engine. Thus, the air volume is equal to the engine displacement. The air mass flow of the intake manifold can be calculated as
m ˙ a i r = S e f f C q C m P i m T i m
where, S e f f is the active area of throttle valve; C q is the flow coefficient; C m represents the quality factor, associating with the engine. Once the engine model is specified, these parameters such as S e f f , C q and C m are determined. T i m and P i m are the temperature and pressure of intake manifold, respectively.
According to the actual working conditions, air mass is then obtained by Refs. [25,26,27,28]:
m a i r = A ( ε ) 0 t m ˙ a i r d t π A ( ε ) m ˙ a i r N = π A ( ε ) N S C q C m P i m T i m ( 1 cos θ )
where, the value for ε is the engine volumetric coefficient and considered constant when the crankshaft speed does not dramatically change. The integration of air mass flow will inevitably produce a constant term. The sum of constant term and non-constant term is air mass. Since there must be a proportional relationship between the constant term and the non-constant term, in order to facilitate the subsequent calculation, the proportional relationship is defined as A ( ε ) , so it is considered to be the engine parameter coefficient related to ε . Once the engine model is specified, these two parameters are determined. In addition, t = π N , and t is opening time of intake valve, and N is the average value of speed in one working cycle.
Different fuels and different fuel ratios will strongly affect the combustion state of the engine [5]. The main elements in the fuel are oxygen, hydrogen, and carbon, with other elements being negligible. Thus, the relationship between mass fractions for the oxygen element w o , hydrogen element w h , and carbon element w c in 1 kg of fuel are shown below:
w c + w h + w o = 1
The theoretical air volume for 1 kg of fuel that is entirely burned is shown in Equation (4):
L O = 22.4 0.21 ( w C 12 + w H 4 w O 32 )
The engine output power is determined by the calorific value of the gas mixture, which can be defined as (considering 1 kg of fuel)
Q m i x = h u 22.4 × ( λ L O 22.4 + 1 M τ )
where, λ is the excess-air factor; h u is the low calorific value of the fuel; and M τ is the relative molecular mass of fuel. Therefore, the power output related to the throttle valve is defined as
Q = [ ( L O ρ a i r + 1 ) m a i r 4 ρ f u e l L O ρ a i r ] h u 22.4 × ( λ L O 22.4 + 1 M τ )
where, ρ a i r and ρ f u e l is the density of air and fuel, respectively. By substituting Equation (2) into Equation (6), Equation (6) then becomes
Q = [ ( L O ρ a i r + 1 ) 1 4 ρ f u e l L O ρ a i r π A ( ε ) h u ( λ L O + 22.4 M τ ) S C q C m P i m T i m ] · ( 1 cos θ ) N = B ( 1 cos θ ) N
where,
B = ( L O ρ a i r + 1 ) 1 4 ρ f u e l L O ρ a i r π A ( ε ) h u ( λ L O + 22.4 M τ ) S C q C m P i m T i m
As widely known, the total energy produced through burning cannot be completely converted to useful work. The process of energy transfer is shown in Figure 2.
η t is the theoretical thermal efficiency; η j is the loss coefficient; η m is the mechanical efficiency; and η e is the engine efficiency.
In a four-cylinder four-stroke engine, each cylinder burns once per cycle, and the crankshaft is rotated 720 degrees in one working cycle.
W e = U k i n ( φ ) + U p o t ( φ ) = J 2 N 2
where, U k i n ( φ ) is the kinetic energy of the system; U p o t ( φ ) is the potential energy of the system; J is the rotational inertia parameter (in this paper, it is equal to a flywheel’s rotational inertia parameter [24]); and φ is the crankshaft angle. By combining Equations (7)–(9), the CSM for Part A is shown in the following equation, and the relationship between N k | k and N k | k + 1 is shown in Figure 3.
η e B ( 1 cos θ ) N k | k = J 2 N k | k + 1 2

2.1.2. Variable Load Condition

The change of engine load will affect crankshaft speed. This signal is discretized by computing the mean value of crank speed in power stroke. Since the time step of power stroke is enough small, Equations (11) and (12) are approximately satisfied:
T T L = J d N d t = J N k + 1 | k + 1 N k | k + 1 t c y c N k + 1 | k + 1 = T T L J t c y c + N k | k + 1
t c y c = 1 2 N k | k + 1 60 = 30 N k | k + 1
where, J is the rotational inertia; T L is the load torque; and t c y c is the working cycle time. T is the engine output torque and can be obtained by the empirical model [25], and its coefficients are obtained by calibrating the actual engine.
T = 181.3 + 379.36 m a i r + 21.91 R A / F 0.85 R A / F 2 + 0.26 σ + 0.0028 σ 2 + 0.027 N k | k + 1 0.000107 N k | k + 1 2 + 0.00048 N k | k + 1 σ + 2.55 σ m a i r 0.05 σ 2 m a i r
where, m a i r is the mass of air in cylinder for combustion (g); σ is the spark advance angle (SAA).
In this case, only the load T L is variable, thereby satisfying Equation (14).
T = C + D B 1 cos θ N k | k + 0.027 N k | k + 1 0.000107 N k | k + 1 2 + 0.00048 N k | k + 1 σ
where C = 181.3 + 21.91 R A / F 0.85 R A / F 2 + 0.26 σ + 0.0028 σ 2 , D = 379.36 + 2.55 σ 0.05 σ 2 .
From this analysis, the model for Part B can be described by:
N k + 1 | k + 1 = C + D B 1 cos θ N k | k + 0.027 N k | k + 1 0.000107 N k | k + 1 2 + 0.00048 N k | k + 1 σ T L J · 30 N k | k + 1 + N k | k + 1
By combining Equations (12) and (15), the speed iteration model can then be established:
N k + 1 | k + 1 = C + D B 1 cos θ N k | k + 0.027 2 η e B ( 1 cos θ ) N k | k J 0.000107 2 η e B ( 1 cos θ ) N k | k J + 0.00048 2 η e B ( 1 cos θ ) N k | k J σ T L J · 30 2 η e B ( 1 cos θ ) N k | k J + 2 η e B ( 1 cos θ ) N k | k J
Equation (16) can be expressed as
N k + 1 | k + 1 = g ( N k | k , θ ) N k + 1 = g ( N k , θ )
where, g ( ) is a nonlinear function.
The relationship among N k | k , N k | k + 1 , and N k + 1 | k + 1 is depicted in Figure 3.

2.2. Design of Extended Kalman Filter

Under actual working conditions, the throttle opening value can be decomposed into two parts: an ideal opening value and throttle position signal noise mainly due to EMI. In measuring the crankshaft speed process, speed data is affected by measurement noise from the crankshaft and process noise from the engine. Engine vibration caused by engine knock, etc., is transmitted to the crankshaft. In addition, the disturbance caused by road surface irregularity is also gradually transmitted to crankshaft. According to Kalman filter theory, these fluctuations are considered as process noise and need to be processed. These noise sources result in a huge accumulation error during pressure measurements, and the Kalman filter (KF) [29,30] is the most preferable filtering method for measurement and process noise. The KF is used to predict an optimal output for the current step based on the optimal output of the previous step and the observed value of the current step. However, for nonlinear systems, the KF cannot achieve optimal prediction. Therefore, the Extended Kalman filter (EKF) [31,32,33] is proposed by using the local linear property of nonlinear systems. State prediction is accomplished by using EKF gain to update state and error covariance. After linearization by computing the nonlinear function’s first order Taylor series and Jacobian matrix at the operating point, the EKF algorithm matches the KF and it comprises two parts: time update and state update. A discrete-time state-space model with a zero mean and a Q k variance in the white process noise ( w k 1 ) is shown as
X k = g ( X k 1 ) + w k 1
where, X k is the state vector of the k t h step, and g ( ) is a nonlinear function.
The state observation can be written as
Z k = h ( x k ) + v k
where v k is the white measurement noise with a zero mean; R k is the variance; and h ( ) is a nonlinear function.
The prediction is shown below:
X ^ k | k 1 = g ( X k 1 )
The first-order Taylor series for g ( ) at X ^ k 1 | k 1 and h ( ) at X ^ k | k 1 are obtained as shown in Equations (21) and (22), respectively.
G k = g X | X ^ k 1 | k 1 = g ( X ^ k 1 | k 1 ) + δ g ( X k 1 ) δ X k 1 | X k 1 = X ^ k 1 | k 1 ( X k 1 X ^ k 1 | k 1 )
H k = h X | X ^ k | k 1 = h ( X ^ k | k 1 ) + δ h ( X k ) δ X k | X k = X ^ k | k 1 ( X k X ^ k | k 1 ) + v k
Let:
Δ k = X k X ^ k
Defining P k + 1 as the covariance of Δ k + 1 , the transcendental error covariance can be calculated as
P k + 1 | k = G k P k G k T + Q k
The major purpose of the EKF is to obtain a minimum-variance state estimate. The gain matrix of the EKF is computed subject to minimization of the estimation error covariance.
T r [ P k + 1 x ] = min K k + 1 ! K k + 1 = P k + 1 | k H k T ( H k T P k + 1 | k H k T + R ) 1
State prediction is accomplished by using EKF gain to update state and error covariance.
X ^ k + 1 = X ^ k + 1 | k + K k + 1 ( Z k H k X ^ k + 1 | k )
P k + 1 = ( E K k + 1 H ) P k + 1 | k
The observed vector is written as follows, where the observed matrix H = 1 :
Z k = H k x k + v k

3. Calculation of Cylinder Pressure

In this paper, due to only related to the crankshaft speed and the amount of air and fuel injection in engine, this method can be applied to multi-cylinder internal combustion (IC) engines including four-cylinder engines. More specifically, many key parameters of ignition engine such as spark advance angle (SAA) need to be considered in the process of building virtual in-cylinder pressure sensor model. Accordingly, at present the application of the virtual sensor is only limited to spark ignition engines. In addition, now there are many ways to describe the status of engine working cycle such as Otto-cycle, Diesel-cycle, Sabtache-cycle, and Atkinson-cycle. Among the many cycles, the Otto-cycle optimally depicts the four-stroke ignition combustion engine. Taking Otto cycle of standard four-stroke engine as an example, but the virtual sensor method is not limited to the engine described by Otto cycle.
In this cycle, the crankshaft rotates 720° as one working period, where the 720° can be divided into four sub-cycles each with a crankshaft rotation of 180° [34]. A diagram of the standard working cycle is shown in Figure 4. These four smaller cycles are known as: intake stroke, compression stroke, power stroke, and exhaust stroke. It is the power stroke that involves the burning and expanding process. According to the Otto-cycle, the cylinder pressure during the intake stroke and the exhaust stroke is equal to atmosphere. The compression stroke and the expanding process are considered as isentropic, and the burning process is regarded as constant volume combustion [35].
The pressure at the end of compression stroke is also the initial pressure of the burning process. The peak occurs at the end of the burning process and is known as the peak pressure of the cylinder. This pressure is an important variate to understand the combustion efficiency of the burning process and to monitor the state of the cylinder.

3.1. Cylinder Pressure Peak Time Confirmation

It’s generally accepted that when the spark advance angle (SAA) is 10–15° after top-dead-center (TDC), combustion efficiency becomes maximum, while engine vibration reaches a minimum [35]. In order to obtain the maximum power of engine output at Optimal-Spark-Advance-Angle (OSAA), the OSAA need to be adjusted continuously.
In the case of high speed and low load of engine, the OSAA should be increased. This is due to the delay period of fuel combustion in the cylinder. The faster the speed is, the bigger the OSSA. On the contrary, with the speed increased, the stronger the turbulence formed by combustion gas in the cylinder is, the faster the combustion speed and the lower the advance angle is. combustion turbulence, which positively affects the burning speed, will increase with faster crank speeds. In this case, the corresponding OSAA delay period will be small and because of the delay caused by fuel burning in the power stroke, the OSAA should be slightly smaller. Additionally, in this paper, the pressure identification based on the power flow is adopted. The comparison between identified pressure signal and actual pressure signal can be realized by computing each power. Hence the SAA is deemed as constant. Due to the interaction and method selected, it’s deemed that the OSAA is constant at 15°, and the maximum power output will occur when the crankshaft angle is 10° after the TDC where pressure in the cylinder reaches its peak value.

3.2. Computation of Cylinder Pressure Peak Value

After the burning process, the exhaust gases in cylinders is primarily composed of C O 2 , H 2 O and N 2 .According to Equations (7) and (8), the specific heat capacity at constant volume for the burned gas is
C m i x ( w C 12 w C 12 + w H 8 C C O 2 + 1 2 w H 8 w C 12 + w H 8 C H 2 O ) 2 9 + C N 2 7 9
where, C C O 2 is the specific heat capacity at a constant volume for C O 2 , and C H 2 O is that for H 2 O .
As computed in Equation (30), Cylinder temperature varies with output power.
Δ T = Q C m i x m
where, m is the mass of the burned gas, and according to conservation of mass theory:
m = m a i r + m f u e l
According to the ideal gas state equation,
p V = n R T
Since the burning process is deemed as isovolumetric, the peak pressure in cylinder is written as:
p p e a k = n R ( Q C m i x ( m a i r + m f u e l ) + T c o m p ) V
where, n is the amount of substance which can be gained from the coefficient of the fuel burning equation; V is the volume of cylinder at top dead center; T c o m p is the temperature at the end of the compression stroke; and R is the ideal gas constant.

4. Modeling of in-Cylinder Pressure

As previously mentioned, the cylinder pressure signal is periodic with variating frequency and amplitude, FAMFS is proposed to modulate the cylinder pressure [36]. A high-factorial Fourier series with variating frequency and amplitude is adopted to modulate the cylinder pressure.

4.1. Frequency-Modulated Fourier-Series

If there is a periodic signal f ( x ) , then for every x there exists a positive value L which makes f ( x ) = f ( x + L ) correct, where L is known as the fundamental period and ω 0 = 2 π L is termed the fundamental frequency. The Fourier-series is written as [37]:
f ( x ) = a 0 + n = 1 ( a n cos n x 2 π L + b n sin n x 2 π L )
where, a 0 is the intercept and a n , b n are coefficients of the Fourier-series.
The frequency-modulated method [38] allows the instantaneous frequency of the carrier signal to vary with the change law of the delivery signal. This modulation can be classified as primary or secondary according to the category of the modulation effects.
With a defined carrier signal x c ( t ) = A cos ( 2 π f c t ) and a transmitted signal y ( t ) , the modulating signal can be written as:
x c ( t ) = A cos ( 2 π 0 t [ f c + f Δ y ( τ ) ] d τ )
where, f c is the center frequency of the carrier signal; A is the amplitude; f c + f Δ y ( τ ) is the instantaneous frequency; and f Δ is the frequency deviation gain.
Based on the engineering practice, the 24-order FAMFS is selected. It should be noted that the identification process is completed offline, and the only online requirement is adjustment of the frequency and amplitude of the FAMFS. Thus, the instantaneity of method is well guaranteed. The 24th order FAMFS is shown below:
f ( x ) = A a 0 + A n = 1 24 ( a n cos n ( 2 π t i = 0 t t N n y ( t i ) ) + b n sin n ( 2 π t i = 0 t t N n y ( t i ) ) )
where, t is the total running time; N n is the number of working cycles; and t i is time sequence.

4.2. Pressure Model

The optimal output of the Kalman observer is the speed of each power stroke. For a four-cylinder engine, the cylinders work in a logical sequence, known as firing order. Generally, the firing order is 1-3-4-2. By sampling the optimal output of the EKF, according to the working sequence, the optimal power-stroke speed for each cylinder is obtained:
N n u m 1 = N 1 + k n N n u m 2 = N 2 + k n N n u m i = N i + k n
where, N n u m i is the speed of the i t h cylinder; N i + k n is the optimal speed of ( i + k n ) t h cycle; and k is a natural number.
According to actual working conditions, the cylinder pressure signal can be deemed as a FAMFS [39,40] signal with a center frequency at zero and a unit gain of frequency deviation. Thus, the pressure identification is ultimately written as:
F ( x ) = A k a 0 + P a t m + A k n = 1 24 ( a n cos n ( 2 π t i = 0 t t N n y ( t i ) ) + b n sin n ( 2 π t i = 0 t t N n y ( t i ) )
y ( t i ) = π N k 15 n = 1 k 2 π ω n t i n = 1 k + 1 2 π ω n
A k = n R ( Q k C m i x | k ( m a i r | k + m f u e l | k ) + T c o m p ) V n = 1 k 2 π ω n t i n = 1 k + 1 2 π ω n
where, t is time, and A k is the theoretical pressure peak at the k t h working cycle.

5. Validation and Results

When the engine is in normal combustion state, its air-fuel ratio usually exceeds 0.85. A ratio of 0.7 or 0.8 will make the burning process of inside cylinder unstable and unpredictable, and this condition only happens in a very special condition like startup in an extremely cold environment or engine malfunctioned. No matter in which condition, it’s not a normal working status so the prediction is not that useful. The validity of the proposed method is verified by comparing the identification value of the virtual cylinder pressure sensor with the measured value through the cylinder pressure sensor.
In the following section, data collected from a genuine 2.0-L four-stroke four-cylinder engine was provided by the FAW Group Corporation R&D Center (Changchun, China). The engine model is EA211 and produced by VW (Volkswagen) (Changchun, China). Data was collected over the course of 100 working cycles adopting Kistler Model 6052A Pressure Sensor to test the performance of the proposed method.

5.1. Set Air-Fuel Ratio Coefficient at 0.85

The parameters and setup of the engine are shown in Figure 5 and Table 1.
Figure 6 compares the crankshaft speed outputs of the EKF with the measured results acquired from the genuine engine. As seen in Figure 6 and Figure 7, from 30 cycles to 70 cycles the engine speed measured showed signs of rising and began to decline after reaching the maximum at 70 cycles. The EKF method accurately tracks this trend and its maximum percentage of error is only 9.6%. This shows acceptable agreement between identified and actual values.
Once the optimal speed from the EKF engine model is predicted, the peak pressure of cylinder can be calculated. Furthermore, the 24th order FAMFS is adopted for identification purposes, and its parameters are shown is Table 2:
A comparison of the actual cylinder pressure with the model-identified values is shown in Figure 8. It shows that in early stage of identification process, the identified values are basically in agreement with the measured values, whether in phase or amplitude. However, the phase deviation between measured and identified values begins to emerge and become larger with time. This is due to the cumulative error caused by the fitting error of the 24-order Fourier series and the speed tracking error.
For Figure 8, the area under the curve is multiplied by the piston area to represent the piston impulse. Hence, the cumulative error percentage of piston impulse is the cumulative error percentage of cylinder pressure. Figure 9 shows the curves of cumulative piston impulse including the actual and identification situation within 100 cycles. The two coincide basically. And Figure 10 shows that its cumulative error percentage is less than 1.8%. This agreement as well as the permitted errors, represents good performance from the proposed method.
In general, the engine is a complex and sophisticated system, and small changes in valve or atmospheric temperature, fuel rail pressure, vibration from irregularity, etc., can severely impact the state of the engine; thus, besides a direct comparison of pressure and the cumulative piston impulse, comparison of piston impulse per cycle is also needed. As depicted in Figure 11, the actual and identified IMEP of each cycle one by one are contrasted. Figure 12 shows that these two error percentages are below 5.4%. And then the proposed method is proven to possesses a great ability to track the IMEP calculated by actual pressure.

5.2. Set Air-Fuel Ratio Coefficient at 0.95

The parameters and setup of the engine are shown in Figure 13.
Figure 14 shows the comparing result of speed from EKF and a genuine engine at A/F ratio 0.95. Error percentages of the speed in Figure 15 is below 1.7%. Calculated and actual values have a great uniformity. Results show that proposed method possessed a great performance on tracking speed.
The peak pressure of cylinder can be calculated by using the predicted speed. Furthermore, the 24th order FAMFS is adopted for identification purposes, and its parameters are the same as Table 2. A comparison of the actual cylinder pressure with the model-identified values is shown in Figure 16.
Figure 16 shows that the comparing result of actual and identified cylinder pressure. Figure 17 also shows the curves of cumulative piston impulse including the actual and identification situation but at A/F ratio 0.95 within 100 cycles. The two coincide also basically. As depicted in Figure 18, the actual and identified IMEP of each cycle one by one are contrasted. Compared with Figure 8 and Figure 16, the proposed method possesses a better ability to track the cylinder pressure at A/F ratio 0.95 than at A/F ratio 0.85. The closer the air-fuel ratio coefficient is to the optimal A/F ratio, the more predictable the combustion state of the engine is. This is also proved by the cumulative error percentage of piston impulse is less than 1.4% in Figure 19, and the error percentage of IMEP in each cycle is less than 4.9% in Figure 20. On the contrary, the phase difference of identified cylinder pressure began to emerge and become larger as the air-fuel ratio coefficient getting far away from optimal Air-Fuel ratio coefficient. But no matter how the A/F ratio coefficient changed, if it is within reasonable range, error percentage of cylinder pressure is within tolerance.

6. Conclusions

As a crucial and critical factor in monitoring the internal state of an engine, pressure identification is significantly essential. In this paper, aimed at solving problems associated with reliability and cost of invasive pressure sensors, virtual cylinder pressure identification sensor based on EKF and FAMFS is proposed. This new method employs three key steps:
(1)
An iterative speed model based on burning theory and Law of energy Conservation. Efficiency coefficient is used to represent operating state of engine from fuel to motion. The iterative speed model associated with the throttle opening value and the crankshaft load.
(2)
The EKF is used to estimate the optimal output of this iteration model. The optimal output of the speed iteration model is utilized to separately compute the frequency and amplitude of the cylinder pressure cycle-to-cycle.
(3)
A pressure fitting algorithm is established by a 24th-order FAMFS. With this process an approximate identification method for engine cylinder pressure is developed in a standard engine’s working cycle.
To further verify the validity of the proposed method, data collected from a genuine engine (EA211) provided by the China FAW Group Corporate R&D Center was used. Test results demonstrate that the proposed method exhibits great performance for tracking crank speed and the engine’s real-time cylinder pressure. By comparing the identified pressure outputs with measurement results, the cumulative error percentage for cylinder pressure was below1.8%, the error percentage of IMEP each cycle was no more than 5.4%, and the error associated with speed was less than 9.6%, when A/F Ratio coefficient was set to 0.85. However, the cumulative error percentage for cylinder pressure was below 1.4%, the error percentage of IMEP each cycle was no more than 4.9%, and the error associated with speed was less than 1.7%, when A/F Ratio coefficient was setup at 0.95. The effectiveness and instantaneity are thereby proven and shown to be capable of achieving the desired accuracy despite the uncertainty of the engine’s power stroke. The closer the A/F ratio set to optimal A/F ratio, the more exact the identified pressure will be. Furthermore, it is important to note that the proposed method is valid, in general, and may be applied to more complex situations.

Author Contributions

Study design and investigation, Q.W. and T.S.; Methodology, Q.W.; Resources, Q.W.; Supervision, D.G.; Writing—original draft, Q.W.; Writing—review & editing, Q.W., T.S., Z.L., D.G.

Funding

The author(s) disclose receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (grant numbers 51806143) and sponsored by Shanghai Sailing Program (No. 19YF1434600).

Conflicts of Interest

The author(s) declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

References

  1. Guardiola, C.; Pla, B.; Bares, P.; Stefanopoulou, A. Cylinder charge composition observation based on in-cylinder pressure measurement. Measurement 2019, 131, 559–568. [Google Scholar] [CrossRef]
  2. Rizzoni, G.A. A stochastic model for the indicated pressure process and the dynamics of the internal combustion engine. IEEE Trans. Veh. Technol. 1989, 38, 180–192. [Google Scholar] [CrossRef]
  3. Min, K.; Chung, J.; Sunwoo, M. Torque balance control for light-duty diesel engines using an individual cylinder IMEP estimation model with a single cylinder pressure sensor. Appl. Therm. Eng. 2016, 109, 440–448. [Google Scholar] [CrossRef]
  4. Bennett, C.; Dunne, J.F.; Trimby, S.; Richardson, D. Engine cylinder pressure reconstruction using crank kinematics and recurrently-trained neural networks. Mech. Syst. Signal Process. 2017, 85, 126–145. [Google Scholar] [CrossRef]
  5. Bodisco, T.; Brown, R.J. Inter-cycle variability of in-cylinder pressure parameters in an ethanol fumigated common rail diesel engine. Energy 2013, 52, 55–65. [Google Scholar] [CrossRef] [Green Version]
  6. Rizvi, M.A.; Bhatti, A.I.; Butt, Q.R. Hybrid model of the gasoline engine for misfire detection. IEEE Trans. Ind. Electron. 2011, 58, 3680–3692. [Google Scholar] [CrossRef]
  7. Al-Durra, A.; Canova, M.; Yurkovich, S. A real-time pressure estimation algorithm for closed-loop combustion control. Mech. Syst. Signal Process. 2013, 38, 411–427. [Google Scholar] [CrossRef]
  8. Cook, J.A.; Powell, B.K. Modeling of an internal combustion engine for control analysis. IEEE Control Syst. Mag. 1988, 8, 20–26. [Google Scholar] [CrossRef]
  9. Sellnau, M.C.; Matekunas, F.A.; Battiston, P.A.; Chang, C.-F.; Lancaster, D.R. Cylinder-Pressure-Based Engine Control Using Pressure-Ratio-Management and Low-Cost Non-Intrusive Cylinder Pressure Sensors. SAE Technical Paper Series 2000, 899–918. [Google Scholar]
  10. Payri, F.; Luján, J.M.; Martín, J.; Abbad, A. Digital signal processing of in-cylinder pressure for combustion diagnosis of internal combustion engines. Mech. Syst. Signal Process. 2010, 24, 1767–1784. [Google Scholar] [CrossRef]
  11. Maurya, R.K. Estimation of optimum number of cycles for combustion analysis using measured in-cylinder pressure signal in conventional CI engine. Measurement 2016, 94, 19–25. [Google Scholar] [CrossRef]
  12. Ahmed, Q.; Bhatti, A.I.; Iqbal, M. Virtual sensors for automotive engine sensors fault diagnosis in second-order sliding modes. IEEE Sens. J. 2011, 11, 1832–1840. [Google Scholar] [CrossRef]
  13. Ferrari, A.; Paolicelli, F. A virtual injection sensor by means of time frequency analysis. Mech. Syst. Signal Process. 2019, 116, 832–842. [Google Scholar] [CrossRef]
  14. Eriksson, L.; Nielsen, L. Towards on-board engine calibration with feedback control incorporating combustion models and ion-sense. Automatisierungstechnik 2003, 51, 204–212. [Google Scholar] [CrossRef]
  15. Huo, X.; Ma, L.; Zhao, X.; Zong, G. Observer-based fuzzy adaptive stabilization of uncertain switched stochastic nonlinear systems with input quantization. J. Frankl. Inst. 2019, 356, 1789–1809. [Google Scholar] [CrossRef]
  16. Zhao, X.; Wang, X.; Zhang, S.; Zong, G. Adaptive neural backstepping control design for a class of non-smooth nonlinear systems. IEEE Trans. Syst. Man, Cybern. Syst. 2018, 1–12. [Google Scholar] [CrossRef]
  17. Balyts’kyi, O.I.; Abramek, K.F.; Shtoeck, T.; Osipowicz, T. Diagnostics of degradation of the lock of a sealing ring according to the loss of working gases of an internal combustion engine. Mater. Sci. 2014, 50, 156–159. [Google Scholar] [CrossRef]
  18. Balyts’kyi, O.I.; Abramek, K.F.; Mruzik, M.; Stoeck, T.; Osipowicz, T. Evaluation of the losses of hydrogen-containing gases in the process of wear of pistons of an internal-combustion engine. Mater. Sci. 2017, 53, 289–294. [Google Scholar] [CrossRef]
  19. Desbazeille, M.; Randall, R.B.; Guillet, F.; Badaoui, M.; Hoisnard, C. Model-based diagnosis of large diesel engines based on angular speed variations of the crankshaft. Mech. Syst. Signal Process. 2010, 24, 1529–1541. [Google Scholar] [CrossRef]
  20. Taglialatela, F.; Lavorgna, M.; Mancaruso, E.; Vaglieco, B.M. Determination of combustion parameters using engine crankshaft speed. Mech. Syst. Signal Process. 2013, 38, 628–633. [Google Scholar] [CrossRef]
  21. Naik, S. Advanced misfire detection using adaptive signal processing. Int. J. Adapt. Control Signal Process. 2004, 18, 181–198. [Google Scholar] [CrossRef]
  22. Zhihua, P.L.Y.J.W.; Xinping, Z.Y.Y. Diagnosing Leakage of Valves in Engines by Analyzing Instantaneous Speed Fluctuations. J. Wuhan Transp. Univ. 2000, 1, 017. [Google Scholar]
  23. Seuling, S.; Hamedovic, H.; Fischer, W.; Schuerg, F. Model based engine speed evaluation for single-cylinder engine control. SAE Tech. Pap. 2012. [Google Scholar]
  24. Chang, X.H.; Wang, Y.M. Peak-to-peak filtering for networked nonlinear DC motor systems with quantization. IEEE Trans. Ind. Inform. 2018, 14, 5378–5388. [Google Scholar] [CrossRef]
  25. Gupta, H.N. Fundamentals of Internal Combustion Engines; PHI Learning Pvt. Ltd.: New Delhi, India, 2012. [Google Scholar]
  26. Metghalchi, M.; Keck, J.C. Burning velocities of mixtures of air with methanol, isooctane, and indolene at high pressure and temperature. Combust. Flame 1982, 48, 191–210. [Google Scholar] [CrossRef]
  27. Colin, O.; Benkenida, A.; Angelberger, C. 3D modeling of mixing, ignition and combustion phenomena in highly stratified gasoline engines. Oil Gas Sci. Technol. 2003, 58, 47–62. [Google Scholar] [CrossRef]
  28. Crowl, D.A.; Louvar, J.F. Chemical Process Safety: Fundamentals with Applications; Pearson Education: London, UK, 2001. [Google Scholar]
  29. Zhang, T.; Liao, Y. Attitude measure system based on extended Kalman filter for multi-rotors. Comput. Electron. Agric. 2017, 134, 19–26. [Google Scholar] [CrossRef]
  30. Helm, S.; Kozek, M.; Jakubek, S. Combustion torque estimation and misfire detection for calibration of combustion engines by parametric Kalman filtering. IEEE Trans. Ind. Electron. 2012, 59, 4326–4337. [Google Scholar] [CrossRef]
  31. Zhao, J.B.; Netto, M.; Mili, L. A robust iterated extended Kalman filter for power system dynamic state estimation. IEEE Trans. Power Syst. 2017, 32, 3205–3216. [Google Scholar] [CrossRef]
  32. Al-Durra, A.; Canova, M.; Yurkovich, S. Application of extended Kalman filter to on-line diesel engine cylinder pressure estimation. In Proceedings of the ASME 2009 Dynamic Systems and Control Conference, American Society of Mechanical Engineers, Hollywood, CA, USA, 12–14 October 2009; pp. 541–548. [Google Scholar]
  33. Cordeiro, T.F.K.; Da Costa, J.P.C.; De Sousa Júnior, R.T.; So, H.C.; Borges, G.A. Improved Kalman-based attitude estimation framework for UAVs via an antenna array. Digit. Signal Process. 2016, 59, 49–65. [Google Scholar] [CrossRef]
  34. Department of Theoretical Mechanics of Harbin Institute of Technology. Theoretical Mechanics, 7th ed.; Higher Education Press: Beijing, China, 2009; pp. 288–322. [Google Scholar]
  35. Yu, Z. Automobile Theory; Machinery Industry Press: Beijing, China, 2009. [Google Scholar]
  36. Yang, Y.; Peng, Z.K.; Dong, X.J.; Zhang, W.M.; Meng, G. General parameterized time-frequency transform. IEEE Trans. Signal Process. 2014, 62, 2751–2764. [Google Scholar] [CrossRef]
  37. Phillips, C.L.; Parr, J.M.; Riskin, E.A. Signals, Systems, and Transforms; Prentice Hall: Upper Saddle River, UK, 2003. [Google Scholar]
  38. Waadeland, C.H. Synthesis of asymmetric movement trajectories in timed rhythmic behavior by means of frequency modulation. Hum. Mov. Sci. 2017, 51, 112–124. [Google Scholar] [CrossRef]
  39. Leclere, Q.; Pezerat, C.; Laulagnet, B.; Polac, L. Application of multi-channel spectral analysis to identify the source of a noise amplitude modulation in a diesel engine operating at idle. Appl. Acoust. 2005, 66, 779–798. [Google Scholar] [CrossRef]
  40. Payri, F.; Olmeda, P.; Guardiola, C.; Martín, J. Adaptive determination of cut-off frequencies for filtering the in-cylinder pressure in diesel engines combustion analysis. Appl. Therm. Eng. 2011, 31, 2869–2876. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flow chart of method.
Figure 1. Flow chart of method.
Sensors 19 03122 g001
Figure 2. Flowchart of Engine energy transfer.
Figure 2. Flowchart of Engine energy transfer.
Sensors 19 03122 g002
Figure 3. Relationship of speed in the iteration model.
Figure 3. Relationship of speed in the iteration model.
Sensors 19 03122 g003
Figure 4. Diagram of the standard working cycle.
Figure 4. Diagram of the standard working cycle.
Sensors 19 03122 g004
Figure 5. Parameters and setup of Engine at A/F ratio 0.85.
Figure 5. Parameters and setup of Engine at A/F ratio 0.85.
Sensors 19 03122 g005
Figure 6. Speed tracking at A/F ratio 0.85.
Figure 6. Speed tracking at A/F ratio 0.85.
Sensors 19 03122 g006
Figure 7. Error percentage of speed tracking at A/F ratio 0.85.
Figure 7. Error percentage of speed tracking at A/F ratio 0.85.
Sensors 19 03122 g007
Figure 8. Pressure identification at A/F ratio 0.85.
Figure 8. Pressure identification at A/F ratio 0.85.
Sensors 19 03122 g008
Figure 9. Piston impulse accumulated at A/F ratio 0.85.
Figure 9. Piston impulse accumulated at A/F ratio 0.85.
Sensors 19 03122 g009
Figure 10. Error percentage of piston impulse accumulated.
Figure 10. Error percentage of piston impulse accumulated.
Sensors 19 03122 g010
Figure 11. IMEP at A/F ratio 0.85 per cycle.
Figure 11. IMEP at A/F ratio 0.85 per cycle.
Sensors 19 03122 g011
Figure 12. Error percentage of IMEP at A/F ratio 0.85 per cycle.
Figure 12. Error percentage of IMEP at A/F ratio 0.85 per cycle.
Sensors 19 03122 g012
Figure 13. Parameters and setup of Engine at A/F ratio 0.95.
Figure 13. Parameters and setup of Engine at A/F ratio 0.95.
Sensors 19 03122 g013
Figure 14. Speed tracking at A/F ratio 0.95.
Figure 14. Speed tracking at A/F ratio 0.95.
Sensors 19 03122 g014
Figure 15. Error percentage of speed tracking at A/F ratio 0.95.
Figure 15. Error percentage of speed tracking at A/F ratio 0.95.
Sensors 19 03122 g015
Figure 16. Pressure identification at A/F ratio 0.95.
Figure 16. Pressure identification at A/F ratio 0.95.
Sensors 19 03122 g016
Figure 17. Piston impulse accumulated at A/F ratio 0.95.
Figure 17. Piston impulse accumulated at A/F ratio 0.95.
Sensors 19 03122 g017
Figure 18. IMEP at A/F ratio 0.95 per cycle.
Figure 18. IMEP at A/F ratio 0.95 per cycle.
Sensors 19 03122 g018
Figure 19. Error percentage of piston impulse accumulated.
Figure 19. Error percentage of piston impulse accumulated.
Sensors 19 03122 g019
Figure 20. Error percentage of IMEP at A/F ratio 0.95 per cycle.
Figure 20. Error percentage of IMEP at A/F ratio 0.95 per cycle.
Sensors 19 03122 g020
Table 1. The parameters and setup of the engine.
Table 1. The parameters and setup of the engine.
A/F RatioAtmosphere
Pressure/KPa
Atmosphere
Temperature/°C
Valve Opening Value/°CA
BTDC
Fuel Rail Pressure
/bar
0.851012530270100
0.95101253027090
Table 2. Parameters of the 24th order Fourier series.
Table 2. Parameters of the 24th order Fourier series.
a0~a6a7~a13a14~a20a21~a24, b1~a3b4~b10b11~b17b18~b24
12.06−2.6390.3240.37835.68−1.7750.4077
−11.920.13020.7365−0.3231−5.1360.99670.2995
−3.7681.882−0.9149−0.00951.2670.3929−0.5088
9.395−1.8540.13230.2571.927−1.1460.2242
−4.8130.28430.596513.11−2.7870.73610.202
−1.2151.142−0.5427−11.711.3650.291−0.3184
3.652−1.264−0.004711.3030.8418−0.82350.1462

Share and Cite

MDPI and ACS Style

Wang, Q.; Sun, T.; Lyu, Z.; Gao, D. A Virtual In-Cylinder Pressure Sensor Based on EKF and Frequency-Amplitude-Modulation Fourier-Series Method. Sensors 2019, 19, 3122. https://doi.org/10.3390/s19143122

AMA Style

Wang Q, Sun T, Lyu Z, Gao D. A Virtual In-Cylinder Pressure Sensor Based on EKF and Frequency-Amplitude-Modulation Fourier-Series Method. Sensors. 2019; 19(14):3122. https://doi.org/10.3390/s19143122

Chicago/Turabian Style

Wang, Qiming, Tao Sun, Zhichao Lyu, and Dawei Gao. 2019. "A Virtual In-Cylinder Pressure Sensor Based on EKF and Frequency-Amplitude-Modulation Fourier-Series Method" Sensors 19, no. 14: 3122. https://doi.org/10.3390/s19143122

APA Style

Wang, Q., Sun, T., Lyu, Z., & Gao, D. (2019). A Virtual In-Cylinder Pressure Sensor Based on EKF and Frequency-Amplitude-Modulation Fourier-Series Method. Sensors, 19(14), 3122. https://doi.org/10.3390/s19143122

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