Next Article in Journal
Monitoring Contractile Cardiomyocytes via Impedance Using Multipurpose Thin Film Ruthenium Oxide Electrodes
Previous Article in Journal
Investigation on Wireless Link for Medical Telemetry Including Impedance Matching of Implanted Antennas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Minimization of Torque Deviation of Cylinder Deactivation Engine through 48V Mild-Hybrid Starter-Generator Control

1
Eco-Vehicle Control Design Team, Hyundai KEFICO Corporation, Gunpo 15849, Korea
2
Department of Automotive Engineering, Hanyang University, Seoul 04763, Korea
3
Department of Mechanical and Automotive Engineering, Keimyung University, Daegu 42601, Korea
4
Research & Development Division, Hyundai Motor Company, Hwaseong 18280, Korea
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(4), 1432; https://doi.org/10.3390/s21041432
Submission received: 24 January 2021 / Revised: 14 February 2021 / Accepted: 16 February 2021 / Published: 18 February 2021
(This article belongs to the Section Physical Sensors)

Abstract

:
Cylinder deactivation (CDA) is an effective technique to improve fuel economy in spark ignition (SI) engines. This technique enhances volumetric efficiency and reduces throttling loss. However, practical implementation is restricted due to torque fluctuations between individual cylinders that cause noise, vibration, and harshness (NVH) issues. To ease torque deviation of the CDA, we propose an in-cylinder pressure based 48V mild-hybrid starter-generator (MHSG) control strategy. The target engine realizes CDA with a specialized engine configuration of separated intake manifolds to independently control the airflow into the cylinders. To handle the complexity of the combined CDA and mild-hybrid system, GT-POWER simulation environment was integrated with a SI turbulent combustion model and 48V MHSG model with actual part specifications. The combustion model is essential for in-cylinder pressure-based control; thus, it is calibrated with actual engine experimental data. The modeling results demonstrate the precise accuracy of the engine cylinder pressures and of quantities such as MAF, MAP, BMEP, and IMEP. The proposed control algorithm also showed remarkable control performance, achieved by instantaneous torque calculation and dynamic compensation, with a 99% maximum reduction rate of engine torque deviation under target CDA operations.

1. Introduction

Regulations for energy and exhaust emissions are major challenges in the research and development of powertrains. To overcome such challenges, cylinder deactivation (CDA) technology has drawn attention as a way to make engine displacement flexible; the process involves deactivating a several cylinders under low-load conditions [1,2,3]. The CDA deactivates several cylinders with wider throttle openings for compensating the required torque. Wider throttle opening reduces pumping loss and thus improves fuel economy [4,5,6,7,8,9,10,11].
Hybridization also has been demonstrated as a promising method to improve the fuel efficiency and reduce exhaust emission from the internal combustion engine [12,13,14]. An optimization of several parameters via a genetic algorithm minimizes the brake specific fuel consumption and nitric oxide emission [15]. The CDA combined with a hybridization improves fuel efficiency by activating only three cylinders out of the six-cylinder engine [16].
Aside from the consideration of the fuel efficiency and exhaust emissions, the successful implementation of CDA should suppress adverse effects such as noise, vibration, and harshness (NVH). Such NVH problem results from torque deviation between cylinders when activating CDA because activated cylinders must provide more torque to maintain vehicle driving performance. The absence of combustion in the deactivation sequences reduces, therefore, the number of power strokes per revolution and extends the interval between the maximum torque pulses, making it difficult to deliver consistent and smooth brake torque from the engine [17,18].
Overcoming the NVH problems of CDA has been challenging in terms of cost of additional components, calibration efforts, and limited control accuracy. Such drawbacks hinder the realization of CDA in the commercial market. The conventional approach for mitigating NVH problems of CDA engines was to use advanced hardware systems, such as additional torsional dampers and suitable engine mounts [19,20]. Although this method alleviates vibration of the drivetrain, cost and weight increased due to the additional hardware equipment. Unlike such hardware systems used to cope with NVH problems, the Dynamic Skip Fire (DSF) system improves NVH performance by selecting diverse firing sequences through independent control of individual cylinders [21,22,23,24,25]. The high degree of freedom to activate DSF operation increases calibration efforts. To complement this DSF system, an additional motor generator unit was integrated with a 48V electric hybridization [26]. The combined system presented the ability to smooth torque pulsation based on mild-hybrid torque assist; however, the open loop scheme makes it difficult to guarantee accurate control performance without any feedback on the current engine torque and dynamics compensation.
For application of the spark ignition (SI) engine model for control purposes, the SI combustion model should strike a balance between model fidelity and computational cost [27]. There have been various types of combustion modeling, which can be grouped into non-predictive models, such as the Wiebe function, and predictive models [28]. Non-predictive models impose a burn rate as a function of crank angle, which means that the burn rate will not be affected by factors such as residual fraction or injection timing. A predictive combustion model would be a good choice for control application because the burn rate will respond appropriately to changes in the variable of interest. One predictive combustion model, the SI turbulent model, predicts the burn rate for homogeneous charges and plays a crucial role in the control application. Even in the high nonlinear conditions where the engine equipped the EGR, water injection, and hydrogen enriched control, the SI turbulent model showed a great fidelity with the computational efficiency [29,30]. The SI turbulent model comprises two subparts of laminar flame speed calculation and entrainment/burn-up. The laminar flame speed is calculated as a function of the equivalence ratio, pressure, temperature, fuel type and dilution effect [31]. The rate of entrained and unburned mixture of fuel and air passing into the flame front through the flame area is proportional to the sum of the turbulent and laminar flame speeds. The burn rate is proportional to the amount of unburned mixture behind the flame front divided by a time constant. This model requires calibration of the four dominant factors of dilution, turbulent flame speed, Taylor microscale length, and flame kernel growth [32,33,34].
In this study, we propose an in-cylinder pressure based 48V mild-hybrid starter-generator (MHSG) control strategy to reduce engine torque deviation under CDA operations. We use commercial software and combine the SI turbulent model for engine modeling with the 48V MHSG model with actual part specifications. In contrast to conventional CDA valve-train systems, the target engine implements CDA with a novel structure of separated intake and exhaust manifolds [35]. To cope with the complexity of the specialized CDA and 48V electrification, the target engine model was first derived using the SI turbulent model based on actual computer-aided design (CAD) files and engine experimental data. The four parameters of the SI turbulent model were calibrated with the in-cylinder pressures obtained by the actual engine experiment. The engine experiment was conducted in 50 operating points, and 100 cycles of in-cylinder pressure data were obtained for each operating condition. The SI turbulent engine model was designed with a strong correlation of representative air states, cylinder pressure, and engine torque. With this engine model, the proposed algorithm calculates instantaneous torque profiles of individual cylinders using in-cylinder pressure traces. Based on the dynamic compensation, the torque setpoint of the 48V MHSG is derived to modulate the slow characteristics of the 48V MHSG from the inherent time delay and response. This is a feedforward control to allow precise control performance with a fast response on a cycle-by-cycle basis for practical torque assist with the 48V MHSG. Moreover, the proposed control structure, where the engine torque is calculated in real-time based on the in-cylinder pressure, effectively reduces the effort needed to calibrate the conventional map-based open-loop controller under various types of engine operations. Furthermore, the proposed algorithm effectively reduces engine torque deviation under the main driving range of the target CDA operations.

2. Experimental Setup

2.1. Target Engine Description

The target engine was a 1.6 L gasoline direct injection (GDI) engine, designed with a specialized architecture of separated intake and exhaust manifolds to accommodate the CDA operation [35]. The detailed working principle of this particular CDA is illustrated in Figure 1. Intake manifold 1 is connected to cylinders 2 and 3. In the same manner, intake manifold 2 is linked to cylinders 1 and 4 to form two different closed airflow loops. These two independent loops allow the engine to operate with four or two cylinders through control of the intake manifold valve and the three-way exhaust valve.
As described in Figure 1a, four-cylinder operation is realized by fully opening the intake manifold valve between the separated intake manifolds and closing the three-way exhaust valve. Except for small cylinder-to-cylinder variations from the air flow through an additional throttle, the engine operates as a normal four-cylinder engine. On the other hand, under two-cylinder operation, as shown in Figure 1b, the intake manifold valve is closed to cut off airflow into intake manifold 2 and the three-way exhaust valve is opened to isolate the internal loop of cylinders 1 and 4. This enables deactivation of cylinders 1 and 4, along with the prevention of fuel injection and ignition. Since the closed loop of the deactivated cylinders does not involve a pressure difference between intake and exhaust strokes, the pumping loss can be eliminated simultaneously.
The target engine is equipped with external high-pressure exhaust gas recirculation (HP-EGR), variable valve timing (VVT), an intercooler, and a 48V electric supercharger, as shown in Table 1. The engine is also equipped with a 48V MHSG, which is connected by a pulley to assist in generating engine power; detailed specifications are provided in Table 2.

2.2. Experimental Conditions

The target engine was modeled with actual experimental engine data from 50 steady states. Figure 2 shows the steady states, represented as value of brake mean effective pressure (BMEP) and engine speed. The normal mode of four-cylinder operation includes a total of 25 naturally aspirated regions, whereas the CDA mode contains 15 naturally aspirated and 10 supercharging areas to track the target BMEP values in only two-cylinder operation. To derive a precise combustion model, 100 cycles of in-cylinder pressure data were obtained for each operating condition. The measured pressure signals were sampled every 1 crank angle degree. Depending on the two-cylinder or four-cylinder operation, each spark timing map has been calibrated to produce maximum brake torque to reflect the trapped mass influence caused by the different activated cylinder configuration. The HP-EGR rate, VVT, spark timing, and intercooler temperature set-points were fixed at specific target values for each engine operating point.

3. Target Engine Model and 48V MHSG Model

3.1. Model Overview

To deal with the enormous complexity of cylinder deactivation combined with a 48V mild hybrid system, the target engine model was designed using GT-POWER simulation. In particular, the simulation model was derived based on 3D CAD files to precisely reflect the geometric information of the specialized engine structure. Figure 3 represents the process of importing the actual 3D CAD files into the simulation environment. The discretization process divides a large volume assembly into smaller connected subcomponents. After the ensuing procedure of simplification to eliminate unnecessary obstacles with respect to the model accuracy, the simulation model analyzes the intake and exhaust flow rates based on 1D fluid dynamics. The forward and backward discharge coefficients were addressed based on the unit experimental test for components of valve types, such as throttle, intake manifold valve, and three-way exhaust valve. The 48V electric supercharger and MHSG were also implemented with the specification values, as shown in the overall model schematics in Figure 4.

3.2. SI Turbulent Engine Model

The in-cylinder pressure is essential information indicating the combustion status of individual cylinders in CDA operations; therefore, the engine model was precisely designed using the SI turbulent flame model of GT-POWER to reflect various in-cylinder conditions such as cylinder geometry, spark location, and mixture characteristics. The SI turbulent flame model calculates the entrained mass of mixture ( M e ) from a rate proportional to the sum of the turbulent flame speed ( S T ) and the laminar flame speed ( S L ). The derivative of burned mass ( M b ) is formulated by dividing the difference of M e and M b by the time constant of the burn rate ( τ b ), as follows [28,36]:
d M e d t =   ρ u A e ( S T + S L ) ,
d M b d t =   M e M b τ b ,
where ρ u is the unburned mixture density, and A e is the surface area at the flame front location.
The four dominant phenomena of dilution, turbulent flame speed, Taylor microscale length, and flame kernel growth determine the flame model [32,33,34]. The dilution effect of residuals and external EGR is sophisticatedly calibrated with the dilution multiplier ( C D E ). The C D E scales the dilution effect term ( E D ) with the burned gas fraction in the unburned zone ( f b ). E D is used as a proportional factor to determine S L based on the maximum laminar speed ( B m ), laminar speed roll-off value ( B ϕ ), equivalence ratio ( ϕ ), and equivalence ratio at maximum speed ( ϕ m ), with the following formula:
E D = 1 0.75 · C D E · ( 1 ( 1 0.75 · C D E · f b ) 7 ) ,
S L = ( B m + B ϕ ( ϕ ϕ m ) 2 ) · ( T u T r e f ) α · ( p p r e f ) β · E D ,
where p is the pressure, T u is the unburned gas temperature, α and β are temperature exponent and pressure exponent values, and T r e f and p r e f are constant values to make p and T u dimensionless.
The flame kernel growth multiplier ( C F K G ) and the turbulent flame speed multiplier ( C T F S ) represent magnitude factors of the turbulent flame speed. C F K G has an especially large influence on the initial combustion velocity with the smaller flame radius ( R f ), as follows:
S T = C T F S u ( 1 1 1 + C F K G ( R f L i ) 2 ) ,
where u is the turbulent intensity, and L i is the integral length scale.
The Taylor length scale multiplier ( C T L S ) is a parameter that modulates the propagation speed of the burn rate with respect to the Taylor microscale length (λ). The λ is calculated based on L i and the turbulent Reynolds number ( R e t ). The C T L S factor is multiplied to calibrate the degree of Taylor length scale for the flame propagation:
λ = C T L S L i R e t ,
τ b = λ S L ,
To calibrate the combustion model, in-cylinder pressure measurements of the engine experimental data were utilized. The goal of combustion model calibration is to determine the best set of four scale multipliers C D E , C F K G , C T F S , and C T L S . The optimal set was derived based on a grid search method using in-cylinder pressure measurements for a total of 50 engine operating conditions.

3.3. 48V MHSG Model

From a control point of view, it is important to accurately model the dynamic characteristics to adequately perform the torque assist function of the 48V MHSG. The inherent electromagnetic and mechanical features cause specific response times. Moreover, time delay also occurs depending on certain periods of the internal control modules and communication events. As shown in Figure 5, these dynamic characteristics were implemented using a first-order system with a time delay:
T m ( s ) =   K τ s + 1 e T d s ,
where the parameter K is the DC gain. The time constant ( τ ) and time delay ( T d ) were generated with component specification values using a look-up table (LUT) for the current engine operating conditions. The 48V MHSG model was connected to the engine flywheel with a certain pulley ratio, rated power, and torque specifications as shown in Table 2.

3.4. Modeling Results

Figure 6 provides a comparison of the simulation results and engine experimental measurements under four-cylinder and two-cylinder engine operations. The four figures show the mass air flow (MAF), manifold absolute pressure (MAP), pressure at upstream of throttle (PUT), and brake specific fuel consumption (BSFC). The performance indices of the coefficient of determination ( R 2 ) and root mean squared error (RMSE) are also summarized in Table 3. The representative air states of MAF, MAP, and PUT showed a strong correlation, with R2 values over 0.97. BSFC, a performance indicator of fuel consumption, also presented an accurate trend, with an R2 value over 0.94 and RMSE value less than 2.4 g/kWh.
Figure 7 show modeling results of in-cylinder pressure traces in normal and CDA modes, respectively. These figures show fitted results of in-cylinder pressure from the combustion model over a wide range of engine operations. For quantitative comparison of the combustion models, the combustion characteristics were calculated from the cylinder pressure, as shown in Figure 8 and Table 4. The maximum value of pressure trace (PcylMax) and its crank angle location (CaPcylMax) had values of less than 0.6 bar and 0.3 degree, respectively, showing notable RMSE accuracy. The crank angle location of the mass fraction burned 50% (MFB50) and the burn duration difference between MFB90 and MFB10 (MFB90-MFB10) were also precise, with RMSE values of less than 0.5 degree.
The mean effective pressure values of BMEP, indicated mean effective pressure (IMEP), pumping mean effective pressure (PMEP), and friction mean effective pressure (FMEP) were also analyzed, as shown in Figure 9. BMEP and IMEP, which denote the engine brake torque and indicated combustion torque, showed precise modeling results with R2 values over 0.99 and RMSE values under 0.04 bar, as detailed in Table 5. Both PMEP and FMEP were also calculated with RMSE values below 0.03 bar. However, FMEP at some operating points showed a discrepancy compared to the experiment. It is caused by the limitation of the linear regression model to estimate the FMEP. The regression model was based on a least-square method, called the Chen-Flynn engine friction model [28]. Overall, these results indicate that the engine model accurately simulates the various engine conditions under normal and CDA operations.

4. In-Cylinder Pressure Based 48V MHSG Control Strategy

4.1. Controller Overview

The design objective of the proposed algorithm is to precisely control the 48V MHSG to reduce engine torque deviation under CDA operations. The control algorithm was designed based on the engine model to efficiently handle the high complexity of the target system. Figure 10 shows the overall architecture of the proposed in-cylinder pressure based 48V MHSG controller. The controller was composed of (1) engine torque calculation, (2) 48V MHSG torque profile generation, and (3) dynamic response compensation. The control input of the controller is the in-cylinder pressure in the previous engine cycle, and it can calculate the real-time engine torque in each cylinder. In addition, the target torque profile of the 48V MHSG is generated using the torque deviations derived from the difference between the activated and deactivated cylinders. In more detail, the desired 48V MHSG torque is calculated every 0.5 degree of the crank angle within one engine cycle. The desired target torque profile can minimize the engine torque deviation in an ideal situation. However, the control strategy requires a rapid response to control the 48V MHSG torque in the crank angle domain. Otherwise, the dynamic response of 48V MHSG could distort the torque profile. Therefore, the last step of dynamic compensation finally produces the torque set-point to modulate the time delay and transient response of the 48V MHSG. To alleviate the time delay, the target profile is shifted to match the control instance to the next complete engine cycle. The transient response, which is approximated by the first-order time constant, compensated by a lead compensator which places a dominant pole to the point where has fast response characteristics.

4.2. Engine Torque Calculation

The in-cylinder pressure is important information that allows instantaneous torque calculation for the individual cylinders. As described in Figure 11, the indicated torque of the individual cylinders is calculated using the in-cylinder pressure ( P c y l ) and the geometry. The total indicated torque applied to the crankshaft T Q i n d ( θ ) is formulated for each crank angle ( θ ) as follows:
T Q i n d ( θ ) = n c P c y l ( θ ) S R s i n ( θ ) P c y l ( θ ) S R c o s ( θ ) t a n ( α ) ,
where n c is the number of cylinders, and S , R , and α are the cross-sectional area of the piston, length of the crank, and connecting rod angle, respectively.
The inertia torque T Q i n e r ( θ ) is the part of the engine torque generated from the crankshaft moment of inertia. The formulation is derived as a negative sign with multiplication of the effective inertia moment ( I e f f ) and derivative of the instantaneous engine speed ( N ˙ ):
T Q i n e r ( θ ) = I e f f 2 π 60 N ˙ ( θ ) ,
The friction torque T Q f r i c is approximated as a mean value of one engine cycle. The FMEP is calculated based on the Chen-Flynn engine friction model [28] with a polynomial equation of the PcylMax and a cycle mean value of the engine speed ( N ¯ ), as follows:
FMEP = C 0 + C 1 PcylMax + C 2 N ¯ + C 3 N ¯ 2 ,
T Q f r i c = n c V d 2 π n r FMEP ,
where n r is the number of crankshaft rotations for a complete engine cycle, and V d is the engine displacement volume. The coefficients C 0 , C 1 , C 2 , and C 3 are derived using the least-square method.
As shown in Figure 12, the engine brake torque T Q e n g ( θ ) is calculated as the sum of T Q i n d ( θ ) , T Q i n e r ( θ ) , and T Q f r i c with the following equation:
T Q e n g ( θ ) = T Q i n d ( θ ) + T Q i n e r ( θ ) + T Q f r i c ,

4.3. 48V MHSG Torque Profile Generation

Based on the calculated engine brake torque, the desired torque profile of the 48V MHSG is generated to eliminate engine torque deviation at each crank angle. Since the sequence of activation and deactivation is repeated every 1/4 cycle in a four-cylinder engine, the deviation can be derived from the brake torque trace shifted by 1/4 cycle, as illustrated in Figure 13. Therefore, the desired profile for the torque assist with 48V MHSG is produced as follows:
T Q d e s ( θ ) = 1 2 · ( T Q e n g ( θ ) + T Q e n g ( θ 2 π 4 ) ) T Q e n g ( θ ) = 1 2 · T Q e n g ( θ ) + 1 2 · T Q e n g ( θ 2 π 4 )

4.4. Dynamic Compensation of 48V MHSG

Accurate torque assist requires a rapid response of 48V MHSG on a cycle-by-cycle basis; however, the inherently slow characteristics and time delay of the 48V MHSG deteriorate the control accuracy for tracking the target torque waveform. Figure 14 shows the distorted torque waveform without any compensation of dynamic response. The grey dashed line is desired target torque profile, and the blue line is actual torque trajectory which cannot tracking the target torque profile.
To resolve this problem, the setpoint applied to the 48V MHSG is finally modulated through dynamic compensation of the desired torque profile. The controller already has two dynamic characteristic parameters as a LUT for the current engine operating conditions: the time delay ( T d ) and time constant ( τ ) explained in Equation (8). Figure 15 shows the sequential procedures for the dynamic compensation with the parameters. The signal in Figure 15a represents a desired torque profile of the 48V MHSG. In ideal conditions, the desired torque profile can minimize the engine torque deviation from the CDA. Each signal in Figure 15b,c shows the control instance after the compensation of the T d and τ , respectively. After both compensations, the control instance is applied to the 48V MHSG which has a dynamic response as shown in Figure 15d. Through the compensation, the actual torque profile in Figure 15e which is generated from the 48V MHSG is almost identical to the desired torque profile.
In the case of the T d compensation, the delay effect is mitigated by matching the control instance to the next complete engine cycle. The target profile is further shifted by the difference between the 48V MHSG time delay and the current cycle interval, calculated from the instantaneous engine speed. For the compensation of the τ , the phase-lead compensator is designed to modulate the slow dynamic characteristics. The phase-lead compensator contributes one zero and one pole to the system. The additional zero eliminates the slow response of MHSG, which is approximated as a first-order time constant τ , through the pole-zero cancellation. And the additional pole improves the response characteristic of the system as a new dominant pole. The formula of the phase-lead compensator is as follows:
T c ( s ) = K c τ s + 1 ω · τ s + 1 ,
where the K c is control gain to compensate a steady-state error and parameter ω is an attenuation factor of the lead compensator which has a value between 0 to 1. The compensator contributes one zero and one pole to the system at −1/ τ and −1/ ω τ , respectively. After the initial value search based on the pole-zero cancelation, the K c and ω values fine-tuned to offset the time constant of the 48V MHSG specifications. These values are also generated as a LUT according to the current engine operating condition.

4.5. 48V MHSG Control Results

As shown in Figure 16, the performance of the proposed control algorithm was evaluated according to the reduction rates of the torque deviations, peak-to-peak values of the engine torque traces, and average power of the 48V MHSG without any regeneration. Figure 17a shows the control results for 48V MHSG at constant engine speed of 1000 rpm, BMEP of 2 bar. Compared to the results without 48V MHSG torque assistance, the proposed algorithm reduced the engine torque deviation by 98% and the peak-to-peak value by 26%. In this operating area, the average power of 48V MHSG was 3.54 kW. These results demonstrate that the proposed control algorithm effectively reduces torque fluctuation between the activated and deactivated cylinders under CDA operations. Moreover, they indicate that the specifications of 48V MHSG are sufficient to mitigate torque variations in this engine operating condition, with working range of the 48V MHSG torque and power within rated values of 100 Nm and 10 kW. On the other hand, the results presented in Figure 17b show relatively small reduction rates of 38% and 10% in the deviation and peak-to-peak indices, despite the MHSG consumes a larger average power of 6.9 kW. This is because torque variation was not completely eliminated in this engine operating area at an engine speed of 3000 rpm and BMEP of 2 bar, due to the limited power of the 48V MHSG.
Figure 18 summarize the control results for each engine operating condition. Torque deviations were calculated with a maximum decrease rate of 98% at an engine speed of 1000 rpm and BMEP of 2 bar, and a minimum reduction rate of 0.1% in the engine operation at an engine speed of 4000 rpm and BMEP of 7 bar. For the reduction rate in peak-to-peak values, a maximum of 26% and minimum of 0.2% were derived under operation with an engine speed of 1000 rpm, BMEP of 2 bar, area at engine speed of 4000 rpm, and BMEP of 4 bar. The average power of the 48V MHSG was proportional to engine speed and BMEP in the range of 1.71 kW to 9.67 kW, and it converged to 48V MHSG’s maximum power of 10 kW. These analysis results indicate that the proposed algorithm practically decreases the torque deviations of the CDA engine with 48V MHSG specifications, mainly under low engine speed and load conditions. In other words, the proposed controller provides effective performance under the primary target driving area of CDA operations.

5. Conclusions

In this paper, an in-cylinder pressure based 48V MHSG control strategy was proposed for minimizing engine torque deviation in the CDA engine. To deal with the significant complexity of the integrated system of CDA and mild-hybrid electrification, the proposed control algorithm was designed in a simulation environment.
The SI turbulent engine model was derived based on 3D CAD data to implement the geometric information of the particular structure with a separated intake and exhaust manifolds. The model included a combustion model that was finely calibrated using in-cylinder pressure measurements, and a 48V MHSG model with dynamic characteristics based on detailed specifications. The modeling results of BMEP and IMEP showed a precise accuracy, with RMSE values under 0.04 bar. PcylMax and CaPcylMax presented RMSE values of less than 0.6 bar and 0.3 degree, respectively.
The performance of the proposed 48V MHSG control algorithm was analyzed with respect to the engine operating range. The reduction rate of the torque deviations and peak-to-peak values were calculated and found to have maximum values of 98% and 26%, respectively. At the maximum reduction rate the average power of 48V MHSG was 3.54 kW without any regeneration. These results demonstrate that the proposed strategy effectively reduced engine torque deviations through accurate control of the 48V MHSG, under target CDA operations of low engine speed and load conditions.
As a follow-up research, the proposed 48V MHSG control algorithm which was evaluated through a simulation will be implemented into an actual vehicle for the investigation of NVH issues during the mode transition such as four-cylinder to two-cylinder mode transition, and vice versa. In addition, battery management considering regenerative braking will be considered to analyze fuel efficiency while running specific driving cycles such as NEDC and WLTP.

Author Contributions

Methodology, H.S., D.J., S.H. and D.H.; Writing—original draft, H.S. and D.J.; Writing—review & editing, M.H.; Supervision, M.H.; Conceptualization, S.H. and D.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Industrial Strategy Technology Development Program (No. 10039673, 10060068, 10079961, 10080284), the International Collaborative Research and Development Program (N0001992) under the Ministry of Trade, Industry and Energy (MOTIE Korea), and National Research Foundation of Korea (NRF) grant funded by the Korean government (MEST) (No. 2011-0017495).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Middendorf, H.; Theobald, J.; Lang, L.; Hartel, K. The 1.4-l TSI Gasoline Engine with Cylinder Deactivation. MTZ Worldw. 2012, 73, 4–9. [Google Scholar] [CrossRef]
  2. Michelini, J.; Glugla, C. Control System Design for Steady State Operation and Mode Switching of an Engine with Cylinder Deactivation. In Proceedings of the American Control Conference, Denver, CO, USA, 4–6 June 2003; Volume 4, pp. 3125–3129. [Google Scholar]
  3. Kortwittenborg, T.; Walter, F. Strategy to Control the Cylinder Deactivation. MTZ Worldw. 2013, 74, 18–22. [Google Scholar] [CrossRef]
  4. Bech, A.; Shayler, P.J.; McGhee, M. The Effects of Cylinder Deactivation on the Thermal Behaviour and Performance of a Three Cylinder Spark Ignition Engine. SAE Int. J. Engines 2016, 9, 1999–2009. [Google Scholar] [CrossRef] [Green Version]
  5. Leone, T.G.; Pozar, M. Fuel economy benefit of cylinder deactivation-sensitivity to vehicle application and operating constraints. In Proceedings of the SAE Technical Papers; SAE International: Warrendale, PA, USA, 2001. [Google Scholar]
  6. Boretti, A.; Scalco, J. Piston and valve deactivation for improved part load performances of internal combustion engines. In Proceedings of the SAE 2011 World Congress and Exhibition; SAE International: Warrendale, PA, USA, 2011. [Google Scholar]
  7. Lee, N.; Park, J.; Lee, J.; Park, K.; Choi, M.; Kim, W. Estimation of Fuel Economy Improvement in Gasoline Vehicle Using Cylinder Deactivation. Energies 2018, 11, 3084. [Google Scholar] [CrossRef] [Green Version]
  8. Parker, M.C.; Jiang, C.; Butcher, D.; Spencer, A.; Garner, C.P.; Witt, D. Impact and observations of cylinder deactivation and reactivation in a downsized gasoline turbocharged direct injection engine. Int. J. Engine Res. 2019, 146808741988281. [Google Scholar] [CrossRef] [Green Version]
  9. Souflas, I.; Mason, B.; Cary, M.; Schaal, P. Mode Transition Optimisation for Variable Displacement Engines. In Proceedings of the SAE Technical Papers; SAE International: Warrendale, PA, USA, 2016. [Google Scholar]
  10. Flierl, R.; Lauer, F.; Breuer, M.; Hannibal, W. Cylinder Deactivation with Mechanically Fully Variable Valve Train. SAE Int. J. Engines 2012, 5, 207–215. [Google Scholar] [CrossRef]
  11. Salim, W.S.I.W.; Mahdi, A.A.M.; Ismail, M.I.; Abas, M.A.; Martinez-Botas, R.F.; Rajoo, S. Benefits of spark-ignition engine fuel-saving technologies under transient part load operations. J. Mech. Eng. Sci. 2017, 11, 3027–3037. [Google Scholar] [CrossRef]
  12. An, F.; Santini, D.J. Mass Impacts on Fuel Economies of Conventional vs. Hybrid Electric Vehicles. SAE Trans. 2004, 113, 258–276. [Google Scholar]
  13. Bitsche, O.; Gutmann, G. Systems for hybrid cars. J. Power Sources 2004, 127, 8–15. [Google Scholar] [CrossRef]
  14. Reid, A. Understanding the Importance of an Energy Crisis. Ph.D. Thesis, University of Michigan, Ann Arbor, MI, USA, 2009. [Google Scholar]
  15. Li, Y.; Wang, S.; Duan, X.; Liu, S.; Liu, J.; Hu, S. Multi-objective energy management for Atkinson cycle engine and series hybrid electric vehicle based on evolutionary NSGA-II algorithm using digital twins. Energy Convers. Manag. 2021, 230, 113788. [Google Scholar] [CrossRef]
  16. Mechtenberg, A.R.; Peters, D.L. Control of a 36 Mode Hybrid Vehicle with Driver Option Selection: Incorporating Urban, Suburban, and Highway Driving; American Society of Mechanical Engineers: Fort Lauderdale, FL, USA, 2012; pp. 137–144. [Google Scholar]
  17. Millo, F.; Mirzaeian, M.; Luisi, S.; Doria, V.; Stroppiana, A. Engine displacement modularity for enhancing automotive s.i. engines efficiency at part load. Fuel 2016, 180, 645–652. [Google Scholar] [CrossRef]
  18. McGhee, M.; Wang, Z.; Bech, A.; Shayler, P.J.; Witt, D. The effects of cylinder deactivation on the thermal behaviour and fuel economy of a three-cylinder direct injection spark ignition gasoline engine. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2019, 233, 2838–2849. [Google Scholar] [CrossRef]
  19. Fujiwara, M.; Kumagai, K.; Segawa, M.; Sato, R.; Tamura, Y. Development of a 6-cylinder gasoline engine with new variable cylinder management technology. In Proceedings of the SAE Technical Papers; SAE International: Warrendale, PA, USA, 2008. [Google Scholar]
  20. Falkowski, A.; Mcelwee, M.; Bonne, M. Design and development of the Daimlerchrysler 5.7L HEMI® engine multi-displacement cylinder deactivation system. In Proceedings of the SAE Technical Papers; SAE International: Warrendale, PA, USA, 2004. [Google Scholar]
  21. Chien, L.C.; Younkins, M.; Wilcutts, M. Modeling and Simulation of Airflow Dynamics in a Dynamic Skip Fire Engine. In Proceedings of the SAE Technical Papers; SAE International: Warrendale, PA, USA, 2015; Volume 2015. [Google Scholar]
  22. Eisazadeh-Far, K.; Younkins, M. Fuel Economy Gains through Dynamic-Skip-Fire in Spark Ignition Engines. In Proceedings of the SAE Technical Papers; SAE International: Warrendale, PA, USA, 2016. [Google Scholar]
  23. Ortiz-Soto, E.; Wolk, B.; Chen, H.; Younkins, M. MDSF: Improved fuel efficiency, drivability and vibrations via dynamic skip fire and miller cycle synergies. In Proceedings of the SAE Technical Papers; SAE International: Warrendale, PA, USA, 2019; Volume 2019. [Google Scholar]
  24. Wilcutts, M.; Switkes, J.; Shost, M.; Tripathi, A. Design and Benefits of Dynamic Skip Fire Strategies for Cylinder Deactivated Engines. SAE Int. J. Engines 2013, 6, 278–288. [Google Scholar] [CrossRef] [Green Version]
  25. Serrano, J.; Routledge, G.; Lo, N.; Shost, M.; Srinivasan, V.; Ghosh, B. Methods of Evaluating and Mitigating NVH when Operating an Engine in Dynamic Skip Fire. SAE Int. J. Engines 2014, 7, 1489–1501. [Google Scholar] [CrossRef]
  26. Wilcutts, M.; Nagashima, M.; Eisazadeh-Far, K.; Younkins, M.; Confer, K. Electrified Dynamic Skip Fire (eDSF): Design and Benefits. In Proceedings of the SAE Technical Papers; SAE International: Warrendale, PA, USA, 2018; Volume 2018. [Google Scholar]
  27. Wang, S.; Prucka, R.; Zhu, Q.; Prucka, M.; Dourra, H. A Real-Time Model for Spark Ignition Engine Combustion Phasing Prediction. SAE Int. J. Engines 2016, 9, 1180–1190. [Google Scholar] [CrossRef]
  28. Gamma Technologies. GT-SUITE Engine Performance Application Manual; Gamma Technologies: Westmont, IL, USA, 2018. [Google Scholar]
  29. Wang, J.; Duan, X.; Liu, Y.; Wang, W.; Liu, J.; Lai, M.-C.; Li, Y.; Guo, G. Numerical investigation of water injection quantity and water injection timing on the thermodynamics, combustion and emissions in a hydrogen enriched lean-burn natural gas SI engine. Int. J. Hydrog. Energy 2020, 45, 17935–17952. [Google Scholar] [CrossRef]
  30. Duan, X.; Liu, Y.; Liu, J.; Lai, M.-C.; Jansons, M.; Guo, G.; Zhang, S.; Tang, Q. Experimental and numerical investigation of the effects of low-pressure, high-pressure and internal EGR configurations on the performance, combustion and emission characteristics in a hydrogen-enriched heavy-duty lean-burn natural gas SI engine. Energy Convers. Manag. 2019, 195, 1319–1333. [Google Scholar] [CrossRef]
  31. Heywood, J.B. Internal Combustion Engine Fundamentals; McGraw-Hill: New York, NY, USA, 1988. [Google Scholar]
  32. Hires, S.D.; Tabaczynski, R.J.; Novak, J.M. The Prediction of Ignition Delay and Combustion Intervals for a Homogeneous Charge, Spark Ignition Engine; SAE International: Warrendale, PA, USA, 1978. [Google Scholar]
  33. Blizard, N.C.; Keck, J.C. Experimental and Theoretical Investigation of Turbulent Burning Model for Internal Combustion Engines; SAE International: Warrendale, PA, USA, 1974. [Google Scholar]
  34. Morel, T.; Rackmil, C.I.; Keribar, R.; Jennings, M.J. Model for Heat Transfer and Combustion In Spark Ignited Engines and its Comparison with Experiments; SAE International: Warrendale, PA, USA, 1988. [Google Scholar]
  35. Han, D.; Park, J.; Lim, H.; Hong, S.; Lee, J.; Choi, K.; Oh, H.; Park, Y.; Chu, D.H.; Lee, H. Engine System. US20,190,128,219A1, 2 May 2019. [Google Scholar]
  36. Zhao, J.; Xi, Q.; Wang, S.; Wang, S. Improving the partial-load fuel economy of 4-cylinder SI engines by combining variable valve timing and cylinder-deactivation through double intake manifolds. Appl. Therm. Eng. 2018, 141, 245–256. [Google Scholar] [CrossRef]
Figure 1. Operations of the target engine: (a) Four-cylinder operation; (b) Two-cylinder operation.
Figure 1. Operations of the target engine: (a) Four-cylinder operation; (b) Two-cylinder operation.
Sensors 21 01432 g001
Figure 2. Engine operating conditions: (a) Normal mode; (b) cylinder deactivation (CDA) mode.
Figure 2. Engine operating conditions: (a) Normal mode; (b) cylinder deactivation (CDA) mode.
Sensors 21 01432 g002
Figure 3. Process of importing 3D computer-aided design (CAD) files with discretization.
Figure 3. Process of importing 3D computer-aided design (CAD) files with discretization.
Sensors 21 01432 g003
Figure 4. Schematics of the target engine model.
Figure 4. Schematics of the target engine model.
Sensors 21 01432 g004
Figure 5. Modeling of the 48V mild-hybrid starter-generator (MHSG) with first-order dynamic response.
Figure 5. Modeling of the 48V mild-hybrid starter-generator (MHSG) with first-order dynamic response.
Sensors 21 01432 g005
Figure 6. Modeling results of the air states and brake specific fuel consumption (BSFC).
Figure 6. Modeling results of the air states and brake specific fuel consumption (BSFC).
Sensors 21 01432 g006
Figure 7. Comparison of the in-cylinder pressure traces: (a) four-cylinder operation; (b) two-cylinder operation.
Figure 7. Comparison of the in-cylinder pressure traces: (a) four-cylinder operation; (b) two-cylinder operation.
Sensors 21 01432 g007
Figure 8. Modeling results of the combustion characteristics.
Figure 8. Modeling results of the combustion characteristics.
Sensors 21 01432 g008
Figure 9. Modeling results of the mean effective pressure values.
Figure 9. Modeling results of the mean effective pressure values.
Sensors 21 01432 g009
Figure 10. Overall process of the proposed 48V MHSG control strategy.
Figure 10. Overall process of the proposed 48V MHSG control strategy.
Sensors 21 01432 g010
Figure 11. In-cylinder pressure and cylinder geometry.
Figure 11. In-cylinder pressure and cylinder geometry.
Sensors 21 01432 g011
Figure 12. Calculated engine torques in crank angle domain.
Figure 12. Calculated engine torques in crank angle domain.
Sensors 21 01432 g012
Figure 13. Calculation of the desired 48V MHSG torque profile in crank angle domain: (a) Average torque; (b) Torque variation; (c) Desired 48V MHSG torque.
Figure 13. Calculation of the desired 48V MHSG torque profile in crank angle domain: (a) Average torque; (b) Torque variation; (c) Desired 48V MHSG torque.
Sensors 21 01432 g013
Figure 14. Actual response of the 48V MHSG without dynamic compensation.
Figure 14. Actual response of the 48V MHSG without dynamic compensation.
Sensors 21 01432 g014
Figure 15. Procedures of the dynamic compensation for time delay and time constant: (a) Desired torque; (b) Compensation of time delay; (c) Compensation of lead compensator; (d) Application of dynamic response in 48V MHSG; (e) Actual response of 48V MHSG.
Figure 15. Procedures of the dynamic compensation for time delay and time constant: (a) Desired torque; (b) Compensation of time delay; (c) Compensation of lead compensator; (d) Application of dynamic response in 48V MHSG; (e) Actual response of 48V MHSG.
Sensors 21 01432 g015
Figure 16. Performance index to evaluate the proposed control strategy.
Figure 16. Performance index to evaluate the proposed control strategy.
Sensors 21 01432 g016
Figure 17. 48V MHSG control results: (a) engine speed of 1000 rpm and brake mean effective pressure (BMEP) of 2 bar; (b) engine speed of 3000 rpm and BMEP of 2 bar.
Figure 17. 48V MHSG control results: (a) engine speed of 1000 rpm and brake mean effective pressure (BMEP) of 2 bar; (b) engine speed of 3000 rpm and BMEP of 2 bar.
Sensors 21 01432 g017
Figure 18. Reduction rate with respect to engine operating conditions: (a) the torque deviations; (b) the peak-to-peak values; (c) the average power consumptions.
Figure 18. Reduction rate with respect to engine operating conditions: (a) the torque deviations; (b) the peak-to-peak values; (c) the average power consumptions.
Sensors 21 01432 g018
Table 1. Specifications of target engine.
Table 1. Specifications of target engine.
DescriptionSpecifications
Engine typeGasoline direct injection (GDI), Inline,
Double overhead camshaft (DOHC)
Number of cylinders4
Bore × stroke (mm)72 × 97
Connecting rod (mm)152.2
Displacement volume (L)1.580
Compression ratio12:1
Applied technologiesCylinder deactivation (CDA) with separated manifolds architecture
External high-pressure EGR (HP-EGR)
Intake/exhaust variable valve timing (VVT)
Intercooler
48V electric supercharger
48V mild-hybrid starter-generator (MHSG)
Table 2. Specifications of 48V MHSG.
Table 2. Specifications of 48V MHSG.
DescriptionSpecifications
Pulley ratio2.67
Maximum and minimum power±10 kW
Maximum and minimum torque±100 Nm
Table 3. Performance index of the modeling results for air states and BSFC: the coefficient of determination ( R 2 ) and root mean squared error (RMSE).
Table 3. Performance index of the modeling results for air states and BSFC: the coefficient of determination ( R 2 ) and root mean squared error (RMSE).
Normal Mode
(Four-Cylinder Operation)
Deactivation Mode
(Two-Cylinder Operation)
R 2 RMSE R 2 RMSE
MAF (g/s)0.9990.0190.9950.154
MAP (bar)0.9940.0020.9830.006
PUT (bar)0.9860.0020.9730.004
BSFC (g/kWh)0.9931.7480.9412.363
Table 4. Performance index of the modeling results for combustion characteristics.
Table 4. Performance index of the modeling results for combustion characteristics.
Normal Mode
(Four-Cylinder Operation)
Deactivation Mode
(Two-Cylinder Operation)
R 2 RMSE R 2 RMSE
PcylMax (bar)0.9720.2610.9390.590
CaPcylMax (deg)0.5600.2210.9240.299
MFB50 (deg)0.4240.2530.9010.365
Burn duration (deg)0.5000.4100.5290.484
Table 5. Performance index of the modeling results for mean effective pressure values.
Table 5. Performance index of the modeling results for mean effective pressure values.
Normal Mode
(Four-Cylinder Operation)
Deactivation Mode
(Two-Cylinder Operation)
R 2 RMSE R 2 RMSE
BMEP (bar)0.9990.0050.9980.011
IMEP (bar)0.9990.0140.9940.037
PMEP (bar)0.9490.0080.8610.013
FMEP (bar)0.9730.0040.2600.022
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shin, H.; Jung, D.; Han, M.; Hong, S.; Han, D. Minimization of Torque Deviation of Cylinder Deactivation Engine through 48V Mild-Hybrid Starter-Generator Control. Sensors 2021, 21, 1432. https://doi.org/10.3390/s21041432

AMA Style

Shin H, Jung D, Han M, Hong S, Han D. Minimization of Torque Deviation of Cylinder Deactivation Engine through 48V Mild-Hybrid Starter-Generator Control. Sensors. 2021; 21(4):1432. https://doi.org/10.3390/s21041432

Chicago/Turabian Style

Shin, Hyunki, Donghyuk Jung, Manbae Han, Seungwoo Hong, and Donghee Han. 2021. "Minimization of Torque Deviation of Cylinder Deactivation Engine through 48V Mild-Hybrid Starter-Generator Control" Sensors 21, no. 4: 1432. https://doi.org/10.3390/s21041432

APA Style

Shin, H., Jung, D., Han, M., Hong, S., & Han, D. (2021). Minimization of Torque Deviation of Cylinder Deactivation Engine through 48V Mild-Hybrid Starter-Generator Control. Sensors, 21(4), 1432. https://doi.org/10.3390/s21041432

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