Model-Based Optimization of
Model-Based Optimization of
Turesson, Gabriel
2018
Document Version:
Publisher's PDF, also known as Version of record
Link to publication
General rights
Unless other specific re-use rights are stated the following general rights apply:
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors
and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the
legal requirements associated with these rights.
• Users may download and print one copy of any publication from the public portal for the purpose of private study
or research.
• You may not further distribute the material or use it for any profit-making activity or commercial gain
• You may freely distribute the URL identifying the publication in the public portal
LUNDUNI
VERSI
TY
PO Box117
22100L und
+4646-2220000
Model-Based Optimization of
Combustion-Engine Control
Gabriel Turesson
The work presented in this thesis is motivated by the need to reliably operate
a compression-ignition engine in a partially premixed combustion (PPC) mode.
Partially premixed combustion is a low temperature combustion concept, where
the ignition delay is prolonged to enhance fuel-air mixing in the combustion
chamber before the start of combustion. A premixed combustion process, in
combination with high levels of exhaust-gas recirculation (EGR), gives low com-
bustion temperatures, which decrease NOx and soot formation. Lowered com-
bustion temperatures also reduce heat-transfer losses which increase the ther-
modynamic engine efficiency. The ignition delay is, however, determined by che-
mical reactions rates, which are sensitive to temperature, gas-mixture composi-
tion, fuel properties and fuel-injection timing. This sensitivity makes PPC more
challenging to operate compared to conventional diesel combustion. Challenges
related to PPC include an increased sensitivity to operating conditions, decre-
ased combustion-timing controllability, high pressure-rise rates, and low com-
bustion efficiency at low engine loads. These challenges put high demands on
the engine control system that needs to be able to adjust fuel-injection timings
and durations to compensate for the combustion sensitivity.
Therefore, this thesis investigates closed-loop combustion control for relia-
ble PPC operation. The feedback loop from pressure-sensor measurement to
fuel-injection actuation is studied in particular. A common theme for the con-
trollers presented is the use of models in the controller design. Either to evaluate
controller performance in simulation, or to optimize engine performance online.
The principle of model predictive control is used for its ability to incorporate mo-
deled system behavior in the controller design, and to control multi-variable sys-
tems with input and output constraints.
The problem of tuning robust and noise insensitive combustion-timing con-
trollers, and its dependence on fuel reactivity is studied in simulation. Simulation
results reveal a nonlinear relation between start of injection and combustion ti-
ming that depends on both load and fuel reactivity. Optimization is used to find
robust and noise-insensitive controller gains. Guidelines for combustion-timing
controller tuning are also presented.
5
Low-order autoignition models are evaluated and compared for the purpose
of model-based controller design. The comparison shows that a simple autoigni-
tion model is sufficient for control of the ignition delay when the cylinder-charge
properties are varied. This model is then used by a model predictive controller to
simultaneously control ignition delay and combustion timing in transient engine
operation, using both gas-exchange and fuel-injection actuation.
The effects of pilot injection on the combustion processes are characterized
experimentally. Experimental results show that a pilot injection can decrease the
main-injection ignition delay and maintain the pressure-rise rate at an accepta-
ble level. This is utilized by a presented fuel-injection controller that manages to
control both combustion timing and pressure-rise rate.
Strategies for improving the low-load performance of PPC are studied expe-
rimentally, where results show that the selection of injection timings and the use
of a pilot injection are important when maximizing the combustion efficiency.
The suggested low-load controller demonstrated a 9 % efficiency increase during
transient engine operation.
This thesis also investigates the design of a controller that utilizes the de-
grees of freedom enabled by multiple injections to efficiently fulfill constraints
on cylinder pressure, NOx emissions and exhaust temperature. A simulation
study shows a potential 2 - 4 % indicated efficiency increase when two injec-
tions are used instead of one. These findings motivated the design of a hybrid
multiple-injection controller that changes the number of injections depending
on operating conditions. The controller designed was capable of reproducing the
found efficiency increase experimentally with respect to constraints on pressure
and NOx emissions.
A model-predictive pressure controller is also introduced. The controller pre-
dicts how the cylinder pressure varies with fuel injection by taking advantage of
the estimated heat-release rate and a cylinder-pressure model. This feature was
used to adjust fuel-injection timings, durations, and number of injections, for
efficient constraint fulfillment in transient engine operation. Experimental re-
sults demonstrate that the pressure controller can also be used for tracking of
cycle-resolved in-cylinder pressure trajectories, as well as finding the most effi-
cient combustion timing.
Heat-release analysis is an essential component in the pressure-sensor feed-
back loop. Methods for calibrating heat-release model parameters with the use
of engine data, and methods for detecting combustion timings are therefore dis-
cussed in the thesis.
The experimental results presented were conducted on a heavy-duty Scania
D13 engine with a modified gas-exchange system. The fuel used was a mixture
(by volume) of 80 % gasoline and 20 % n-heptane, to elevate the fuel octane num-
ber and allow for longer ignition delays.
6
Acknowledgments
I would like to thank my supervisor, Rolf Johansson, for his consistent support
and advice, for sharing his knowledge of scientific work, and for setting high stan-
dards in the project. My co-supervisor Per Tunestål has been very helpful with
discussions and suggestions on theoretical and practical matters related both to
combustion engines and control. Bo Bernhardsson gave interesting courses in
control and optimization that contributed to the project.
The experimental work was done in collaboration with my colleague, Lian-
hao Yin. His knowledge about, and experience with laboratory work and pro-
gramming made it possible for us to quickly build the test rig and conduct the
experiments presented in the thesis. The technicians at the Division of Combus-
tion Engines were very helpful in keeping the engine running and re-building the
experimental platform when necessary. I am also grateful to Ola Stenlåås for ta-
king a great interest in our project and for giving constructive feedback on our
progress. I wish to express my gratitude to the colleagues who took their time
and helped me with the proofreading of this thesis.
I also extend my acknowledgments to colleagues at the Department of Auto-
matic Control and the Department of Energy Sciences for interesting discussions
on research and related subjects, and to the technical and administrative staff at
the departments for assistance on various subjects related to my research. Spe-
cial thanks goes to Leif Andersson for his support in the finalization of this thesis.
Finally, I would like to thank my wife Desire, my family, and my friends for
their encouragement and support.
Financial Support
The author would like to acknowledge the Competence Center for Combus-
tion Processes (KCFP) and the Swedish Energy Agency for the financial support
(project number 22485-3), Scania for supplying the experimental engine, and the
KFCP PPC Control reference group for consistent feedback on the work. The au-
thor is also a member of the LCCC Linnaeus Center and the ELLIIT Excellence
Center at Lund University.
7
Contents
Nomenclature 13
1. Introduction 17
1.1 Combustion-Engine Challenges . . . . . . . . . . . . . . . . . . . 17
1.2 Fundamental Engine Principles . . . . . . . . . . . . . . . . . . . 18
1.3 Low Temperature Combustion . . . . . . . . . . . . . . . . . . . . 23
1.4 Engine Control . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
1.5 Outline and Contributions . . . . . . . . . . . . . . . . . . . . . . 31
2. Modeling 37
2.1 Control-Oriented Modeling . . . . . . . . . . . . . . . . . . . . . . 37
2.2 In-Cylinder Model . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.3 Gas-Exchange System . . . . . . . . . . . . . . . . . . . . . . . . . 45
2.4 Ignition-Delay Modeling . . . . . . . . . . . . . . . . . . . . . . . 50
2.5 NOx -formation Model . . . . . . . . . . . . . . . . . . . . . . . . . 66
3. Control and Estimation Methods 70
3.1 Control Concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
3.2 State Estimation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
4. Heat-Release Analysis 77
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
4.2 Problem Formulation . . . . . . . . . . . . . . . . . . . . . . . . . 81
4.3 Filter Configurations . . . . . . . . . . . . . . . . . . . . . . . . . . 84
4.4 Simulation Results . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
4.5 Experimental Results . . . . . . . . . . . . . . . . . . . . . . . . . 91
4.6 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
4.7 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
5. Experimental Setup 96
5.1 The Scania D13 Engine . . . . . . . . . . . . . . . . . . . . . . . . 96
5.2 Instrumentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
5.3 Control-System Architecture . . . . . . . . . . . . . . . . . . . . . 100
9
Contents
10
Contents
11
Nomenclature
Symbol Descriptions
The table below summarizes the most frequently used notation in the thesis.
Notation Description
A Cylinder area
B Cylinder bore
C1 & C2 Heat-transfer coefficients
cp Specific heat at constant pressure
cv Specific heat at constant volume
γ Ratio of specific heats
d p max Maximum cylinder pressure derivative
Ea Activation energy
η GIE Gross indicated efficiency
η NIE Net indicated efficiency
η th Thermodynamic efficiency
λ Relative air/fuel ratio
φ Fuel/air equivalence ratio
ṁ air Air mass flow
ṁ EGR EGR mass flow
ṁ f Fuel mass flow
mf Fuel mass
NOx Oxides of nitrogen
Nspeed Engine speed
p In-cylinder pressure
p ex Exhaust-manifold pressure
p in Intake-manifold pressure
p IMEPg Gross indicated mean effective pressure
p IMEPn Net indicated mean effective pressure
13
Nomenclature
14
Nomenclature
Abbreviation Description
0D Zero dimensional
BC Boundary condition
CAD Crank angle degree
CDC Conventional diesel combustion
CFD Computational fluid dynamics
DOI Duration of injection
EGR Exhaust gas recirculation
EVO Exhaust valve opening
EKF Extended Kalman filter
FTM Fast thermal management
HC Hydrocarbon
HCCI Homogeneous charge compression ignition
IAE Integrated absolute error
ICE Internal combustion engine
IMEP Indicated mean effective pressure
IVC Inlet valve closing
KF Kalman filter
LTC Low temperature combustion
MPC Model predictive control
MVM Mean value model
NTC Negative temperature coefficient
ON Octane number
PF Particle filter
PI Proportional integral
PPC Partially premixed combustion
PM Particulate matter
PRF Primary reference fuel
rpm Revolutions per minute
RMSE Root mean square error
SOC Start of combustion
SOI Start of injection
TDC Top-dead-center
QP Quadratic program
VI Virtual instrument
15
1
Introduction
17
Chapter 1. Introduction
18
1.2 Fundamental Engine Principles
0.4
0.2
0
0 2 4 6 8
NOx [g/kWh]
Figure 1.1 Legislated emission levels for heavy-duty vehicles in the European
Union (Euro I-VI) during the past 20 years [EU, 2007]. The emission goals have
been met by gradual improvements of emission control techniques and improved
fuel quality. Improvements from Euro IV were done with the help of exhaust-gas
after treatment.
due to downward motion of the piston. The air is then compressed during the
compression stoke and fuel is injected as the piston approaches top-dead-center
(TDC). At this point, the temperature is sufficiently high for autoignition to oc-
cur. Combustion leads to a pressure increase that generates work during the ex-
pansion stroke. Combustion products are then scavenged as the piston moves
upward with the exhaust-valves open during the exhaust stroke. The four stro-
kes and corresponding cylinder pressure and volume curves are presented in
Figs. 1.3 and 1.4. The cylinder pressure was obtained from the engine used in
the experimental work presented in the thesis, whereas the volume was compu-
ted from the cylinder geometry. Fuel injection is indicated by the injector-current
pulse located before TDC.
19
Chapter 1. Introduction
injector
intake valve exhaust valve
TDC
combustion
chamber
piston
BDC
connecting rod θ
crankshaft
Figure 1.2 The basic geometry of an engine cylinder. Combustion of fuel and
air in the combustion chamber leads to a pressure increase that generates linear
piston motion. The linear motion is converted to rotational motion of the crank-
shaft. Flow of fuel and air are governed by poppet valves, a fuel injector and the
motion of the piston. The acronyms TDC and BDC stand for top-dead-center
and bottom-dead-center, indicating the top and bottom positions of the piston.
Crank-angle degree is denoted θ.
the lift-off length, the fuel has vaporized completely. Chemical reactions are ini-
tiated all over the jet cross section after further air entrainment. Initial reactions
are followed by rapid, rich, premixed combustion and the resulting temperature
increase leads to formation of soot due to an excess of fuel. Air entrainment con-
tinues as the reacting fuel travels along the spray axis and a quasi-steady diffu-
sion flame is formed along the jet periphery when stoichiometric conditions are
reached. At this stage, the combustion rate is controlled by how fast the injected
fuel is vaporized, mixed with air and supplied to the diffusion flame. This type
of combustion is therefore referred to as mixing controlled. Temperature reaches
its maximum in the vicinity of the flame which causes nitrogen to oxidize and
form harmful NOx emissions. High temperature in combination with availability
of oxygen also result in soot oxidation which gives the characteristic diesel-flame
luminosity. Furthermore, these conditions lead to almost complete oxidation of
20
1.2 Fundamental Engine Principles
21
Chapter 1. Introduction
Cylinder Pressure
100
p [bar] Intake Compression Expansion Exhaust
fuel injection
50
0
−300 −200 −100 0 100 200 300
θ [CAD]
Cylinder Volume
·10−3
2
V [m3 ]
0
−300 −200 −100 0 100 200 300
θ [CAD]
Figure 1.4 Cylinder pressure and volume during the engine cycle. The cylinder
pressure was measured from the engine used for the experimental work presented
in this thesis, whereas the volume was computed from the cylinder geometry. Fuel
injection is indicated by the injector current pulse located before TDC.
The journey of a fuel particle along the diesel-jet axis is indicated by the
orange trajectory in this diagram. The fuel/air mixture starts rich (high φ) with
fairly low temperature after vaporization (1). Premixed rich combustion is then
initiated as the fuel mixes with air (φ ∼ 2 − 4). Temperature increases steeply and
the rich, high-T , soot-formation region is reached (2). The fuel becomes incre-
asingly diluted and the temperature peaks at the diffusion flame boundary at
close to stoichiometric conditions (3). The lean, high-T region is favorable for
thermal NOx formation and the availability of oxygen promotes soot oxidation
as the burned mixture cools down during the expansion stroke (4). The inter-
section of the NOx -formation zone with the soot-oxidation zone leads to a fun-
damental diesel-combustion trade-off between NOx and soot emissions. If NOx
formation is to be avoided, soot oxidation will simultaneously decrease with hig-
her soot-emission levels as a result.
It should be noted, however, that this conceptual model only gives a quali-
tative understanding of diesel-combustion characteristics. In reality, there is a
distribution of φ and T around the trajectory in Fig. 1.6, [Kitamura et al., 2002].
22
1.3 Low Temperature Combustion
Figure 1.5 Conceptual model of conventional diesel combustion. The left dia-
gram shows a fuel jet 5 crank-angle degrees (CAD) after injection, where liquid
fuel has vaporized and mixed with air. Initial premixed combustion is indicated
in purple. The diagram to the right shows the same jet 1.5 CAD later, where a hot
diffusion-flame front (green) has been established. Source: [Musculus, 2006].
HCCI
The interest in LTC emerged with the discovery of homogeneous charge com-
pression ignition (HCCI). A concept where a diluted homogeneous charge is in-
ducted and autoignitioned by compression. It was first studied in two-stroke en-
gines [Onishi et al., 1979] and was later shown to yield low emission levels in com-
bination with high efficiencies in the low-load operating region of a four-stroke
engine [Epping et al., 2002]. The blue line in Fig. 1.6 indicates the path taken
by a fuel element in the case of ideal HCCI combustion. It is assumed that the
fuel is completely mixed with air prior to the start of combustion, φ < 1. The
HCCI trajectory shows lean, premixed combustion that avoids the soot forma-
23
Chapter 1. Introduction
5
25 % Soot Formation
4
Local Equivalence Ratio, φ [-]
(1)
3
1%
(2)
Diesel Combustion
2
500 ppm
PPC
(3)
1
HCCI (4)
5000 ppm
NOx Formation
Soot Oxidation
0
600 1,000 1,400 1,800 2,200 2,600 3,000
Local Temperature, T [K]
tion region. Furthermore, the reduced temperature results in low NOx -emission
levels. In HCCI, the combustion timing is completely determined by chemical
autoignition kinetics as compared to injection-controlled diesel combustion.
Combustion timing can therefore only be controlled by varying the tempera-
ture and mixture composition during the compression stroke. This poses a chal-
lenge with respect to combustion-timing sensitivity and controllability. Alterna-
tive combustion-timing strategies, such as variable-valve timing [Agrell et al.,
2003], variable compression ratio [Haraldsson et al., 2002] and dual-fuel opera-
tion [Olsson et al., 2001] have been proposed to control the combustion timing in
HCCI. Another challenge with HCCI is the limited operating range due to violent
combustion rates during high-load operation [Olsson et al., 2004].
24
1.3 Low Temperature Combustion
PPC
Low temperature combustion can also be obtained by prolonging the ignition
delay of conventional diesel combustion. This can be done by diluting the in-
ducted air charge with EGR, and by injecting fuel earlier during the compression
stroke, or later during the expansion stroke when the temperature is lower [Ta-
keda and Keiichi, 1995; Kimura et al., 1999; Akihama et al., 2001]. With direct in-
jection, the fuel-air mixture obtained is not as homogeneous as in HCCI. On the
other hand, combustion timing is only to a limited degree controlled by chemical
kinetics. Emission levels can therefore be reduced with maintained combustion
controllability. These concepts are commonly referred to as partially premixed,
and has been applied by means of different techniques, leading to a number of
names and abbreviations existing in the literature. Examples are MK, PCCI and
RCCI [Kimura et al., 1999; Kanda et al., 2005; Kokjohn et al., 2011]. This thesis ad-
dresses topics related to a concept called partially premixed combustion (PPC).
The conceptual model developed by Musculus [2006] is summarized below for
the purpose of describing PPC.
In PPC, the fuel takes an intermediate path through the φ − T diagram. This
is illustrated by the green line in Fig. 1.6. Injection occurs earlier during the com-
pression stroke, which increases the fuel-traveling distance prior to vaporization
due to reduced mixture temperature and density. The reduced temperature and
increased EGR dilution extend the ignition delay and make it possible for the
majority of the fuel to vaporize before the start of combustion. Increased mix-
ing leads to a more spatially distributed fuel jet as compared to CDC. Combus-
tion is therefore initiated more uniformly over the jet cross section. The incre-
ased combustion-zone heat capacity reduces combustion temperature, and the
simultaneous reduction of φ and T reduces formation of both NOx and soot. Two
PPC fuel-jet intersections are presented in Fig. 1.7. Here, the fuel starts to react
later (12 CAD after injection) as compared to Fig. 1.5, which allows for complete
fuel vaporization before the start of combustion. The fuel jet also occupies a lar-
ger region during combustion. This leads to an increased dilution with smaller
regions of soot formation.
It has been proven difficult to obtain sufficient ignition delay for premixed
combustion with diesel fuels at high-load conditions [Noehre et al. 2006]. This
could be remedied by increasing the fuel autoignition resistance. Gasoline PPC
was introduced by Hildingsson et al. [2006] who showed that longer ignition de-
lays could be achieved even at high-load conditions. At the author’s engine lab-
oratory, Manente et al. [2010c] showed that gasoline PPC could achieve gross in-
dicated efficiencies between 52 and 55 % from idle to 26 bar (indicated mean
effective pressure), with Euro VI compatible emission levels. This was achieved
with EGR ratios at 50 %, φ at 0.75, and an advanced fuel-injection strategy. Simu-
lation results in [Fridriksson et al., 2011] attributed the high efficiency obtained
to low heat-transfer losses resulting from low combustion temperatures. These
25
Chapter 1. Introduction
Figure 1.7 Conceptual model of low load, single-injection, and EGR-diluted par-
tially premixed combustion (PPC). In PPC, the fuel reacts later compared to con-
ventional diesel combustion, see Fig. 1.5. This allows for complete fuel vapori-
zation before the start of combustion. The fuel jet also occupies a larger region
during combustion, which leads to increased dilution with smaller regions of soot
formation. Source: [Musculus, 2006].
results were the main motivation for the work presented in this thesis.
The increased dependency on chemical kinetics makes PPC more sen-
sitive to operating conditions, as compared to CDC. Variations in tempera-
ture, dilution and fuel reactivity might lead to undesired combustion in the
emission-formation regions in Fig. 1.6. Other PPC challenges include violent
combustion rates, increased cylinder-to-cylinder variation, and misfire if the
fuel does not ignite properly. This thesis therefore studies the problem of con-
trolling PPC. With the objective of advancing the concept from manual operation
in steady state, towards autonomous and transient operation in a multi-cylinder
engine.
26
1.4 Engine Control
Disturbance
r P e u y
Controller Engine
−1
Figure 1.8 A system controlled in closed loop. The controller decides u based on
deviation between the measured system output y and a set point r .
27
Chapter 1. Introduction
80 p
p [bar] 60
40 injector current
20 dQ c /d θ
0
−30 −20 −10 0 10 20 30 40
θ [CAD]
Figure 1.9 This thesis investigates the problem of how to decide timings and
durations of fuel-injection pulses to control the combustion processes and the
in-cylinder pressure p. This is a control problem with several degrees of freedom
and a highly informative measured system output as feedback signal.
red with piezoelectric pressure transducers, mounted in the cylinder head. These
sensors utilize the piezoelectric effect where a charge is generated when a pie-
zoelectric crystal is exposed to a force. The measured in-cylinder pressure can
be used to compute indicated engine work, heat-release rate and NOx forma-
tion. The heat-release rate dQ c /d θ can be used to compute combustion timings
and ignition delays, which are important indicators for efficiency and emission
formation. Although cylinder pressure sensing is a widely used in engine rese-
arch, development, and calibration, in-cylinder pressure sensors have not yet
reached widespread use in production vehicles due to high technical demands
and associated cost. Recent announcements still indicate that cylinder pressure
sensing might be used in future production vehicles [BorgWarner, 2014; Nagatsu
et al., 2017]. There are several reviews describing the potential use for pressure
sensors in engine control and diagnostics, see [Powell, 1993; Iorio et al., 2003;
Eriksson and Thomasson, 2017].
Cycle-resolved input and output signals are presented as a function of θ in
Fig. 1.9. The figure shows injector-current pulses, the measured in-cylinder pres-
sure p and the computed heat-release rate dQ c /d θ.
Favorable cylinder boundary conditions in terms of intake temperature, pres-
sure and composition were in this work obtained by regulating mass flows
through EGR paths and a charge-air cooler. This was done by actuating valve po-
sitions in the gas-exchange system. Sensors in the gas-exchange system measu-
red temperatures, pressures, air-mass flow and exhaust oxygen concentration. A
more detailed description of the experimental setup and the gas-exchange sys-
tem is given in chapters 2 and 5.
28
1.4 Engine Control
Optimal Control
A common theme for the work presented in this thesis is to represent control
problems as mathematical optimization problems with the system input u as the
optimization variable
29
Chapter 1. Introduction
100 100
open loop closed loop
a) b)
80 80
p [bar]
p [bar]
60 60
40 40
20 20
−10 0 10 20 −10 0 10 20
θ [CAD] θ [CAD]
100 100
misfire
c) d)
80 80
p [bar]
p [bar]
60 60
40 40
knock
20 20
−10 0 10 20 −10 0 10 20
θ [CAD] θ [CAD]
Figure 1.10 This thesis investigates the use of closed-loop control for
reliable PPC operation. Presented controller designs manage to: Reduce
cylinder-to-cylinder variation for consistent and efficient work output with suf-
ficient ignition delay; Compare (a) where same injection durations and timing
are actuated to the different cylinder, with (b), where work output and com-
bustion timings are regulated in closed loop; Control the pressure-rise rate with
pilot-injection adjustment to avoid knock and maintain an efficient combustion
timing (c); Improve low-load operation by adjusting intake conditions and fuel
injection to avoid misfire and incomplete combustion (d).
The approach taken here was to let the controller repeatedly solve (1.1) on
a cycle-to-cycle basis with respect to measured y and a time horizon of future
inputs, states and outputs. This optimal-control technique is called model pre-
dictive control (MPC) and has gained attention in several context, for example
process control, automotive applications and combustion-engine control. Mo-
del predictive control will be described in greater detail in Chapter 3.
Chapters 10 to 12 of this thesis investigates the optimization problem of how
pressure-sensor feedback and actuation of a number of fuel injections could
be combined to efficiently fulfill constraints on cylinder pressure, NOx forma-
tion and exhaust-gas temperature. The aim of this investigation was to design a
controller capable of automatically finding efficient fuel-injection timings, du-
rations and number of injections, as a function of the engine operating point.
30
1.5 Outline and Contributions
Measurement
100 Prediction
p
p limit Constraint
80
NOx limit NOx
p [bar]
60
40
20
dQ c /d θ
0
−20 −10 0 10 20 30 40
θ [CAD]
Two types of controllers were studied. The first controller was designed with sim-
ple proportional/integral controller components and heuristic constraint hand-
ling. The second controller employed the MPC principle and utilized methods for
cylinder-pressure approximation and heat-release analysis to predict how the cy-
linder pressure varies with fuel-injection parameters. This principle is illustrated
in Fig. 1.11, where a controller predicts future-cycle fuel-injection adjustment
and cylinder-pressure variation (dashed) based on a pressure model and the me-
asured pressure from the previous cycle (solid). This allows the controller to re-
peatedly optimize fuel-injection timings and durations on a cycle-to-cycle ba-
sis, so that the predicted engine outputs efficiently fulfill specified constraints
(dash-dotted).
31
Chapter 1. Introduction
analysis and the experimental platform. The main contributions of the thesis are
presented in chapters 6 to 12, where chapters 6 to 9 cover engine experiments
and controller designs related to partially premixed combustion. Chapters 10 to
12 are focused on fuel-injection optimization and constraint fulfillment, where
the results are also applicable to conventional compression-ignition engines.
A more detailed description of the chapters are given below along with refe-
rences to publications on which they are based. Preliminary versions of parts of
the research presented in this thesis was published in the licentiate thesis by the
author:
Chapter 2
This chapter presents models used for simulation, state estimation and control-
ler design. It includes control-oriented models of the gas-exchange system and
in-cylinder processes. Model-calibration results are also presented. The main
contribution of this chapter is an evaluation of low-order ignition-delay models
for the purpose of control. The results show that a fairly simple model can be
used to predict the relation between intake conditions and ignition delay. The re-
lation between injection timing and ignition delay when the gain from injection
timing to ignition delay changes sign was, however, not adequately captured by
the models considered.
Related Publication
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2017). “An investigation on
ignition-delay modeling for control”. Int. J. Powertrains 6:3, pp. 282–306.
Chapter 3
The control principles used are presented in this chapter with emphasis on mo-
del predictive control. This chapter also presents the state-estimation methods
used to estimate heat-release model parameters and EGR mass flow in subse-
quent chapters.
Chapter 4
This chapter introduces heat-release analysis methods used to extract combus-
tion information from in-cylinder pressure measurement. A heat-release detec-
tion method for multimodal heat-release feedback is also presented. The pro-
blem of calibrating unknown heat-release model parameters is represented as
32
1.5 Outline and Contributions
Related Publication
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2016). “Evaluation of nonli-
near estimation methods for calibration of a heat-release model”. SAE Int. J.
Engines 9:2, pp. 1191–1200.
Chapter 5
The setup used in the experimental work is presented in this chapter. Engine spe-
cifications, sensors and actuators are presented together with a description of the
control-system architecture.
Chapter 6
Chapter 6 investigates proportional-integral (PI) combustion-timing controller
design as a function of fuel reactivity, indicated by the fuel octane number. The
investigation was done through simulation and describes how PI controllers
should be tuned for maximized disturbance rejection subject to constraints on
robustness and cycle-to-cycle variation. The obtained results present challenges
and limiting factors for combustion-timing controller performance.
Related Publication
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2018). “Proportional–integral
controller design for combustion-timing feedback, from n-heptane to
iso-octane in compression–ignition engines”. J. Dynamic Systems, Mea-
surement, and Control 140:5, p. 054502.
Chapter 7
This chapter covers simultaneous control of ignition delay and combustion
timing through combined actuation of fuel injection and valve positions in
the gas-exchange system. The suggested model predictive controller utilizes
a physics-based ignition-delay model, previously presented and calibrated in
Chapter 2. The controller was evaluated experimentally and was shown capa-
ble of tracking set points with respect to cylinder-individual combustion timings
and the cylinder-mean ignition delay during engine load and speed changes.
Related Publication
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2015). “Simultaneous control
of combustion timing and ignition delay in multi-cylinder partially premixed
combustion”. SAE Int. J. Engines 8:5, pp. 2089–2098.
33
Chapter 1. Introduction
Chapter 8
Chapter 8 investigates the effects of pilot injection for control of the pressure-rise
rate. An experimental evaluation on how pilot injection affects emission levels
and efficiency is first presented. The experimental results are then used to de-
sign a model predictive controller with the objective to fulfill an upper limit on
pressure-rise rate, and to track combustion-timing and engine-load set points.
Experimental controller-evaluation results are also presented.
Related Publication
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2016). “A double-injection
control strategy for partially premixed combustion”. In: Proc. 8th IFAC Sym-
posium on Advances in Automotive Control (AAC 2016). Vol. 49. 11. Norrkö-
ping, Sweden, pp. 353–360.
Chapter 9
Chapter 9 explores control strategies for improved combustion efficiency at
low-load operation. The results presented were based on experimental engine
data. Saturation of the fuel-injection timing and the introduction of a pilot in-
jection increased the indicated efficiency. Gas-exchange actuation for avoidance
of low-efficiency regions in a φ − T diagram was found through simulation of a
calibrated gas-exchange system model.
Related Publication
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2016). “Control of the low-load
region in partially premixed combustion”. In: Proc. J. Physics: Conference Se-
ries. Vol. 744. 1. Southampton, England.
Chapter 10
Chapter 10 investigates potential efficiency improvements with multimodal
heat-release rates, when constraints on maximum cylinder pressure, NOx and
exhaust temperature are imposed. A simulation study showing an efficiency in-
crease with two injections is first presented. The simulation results suggest a
heuristic hybrid fuel-injection controller that varies the number of injections de-
pending on operating conditions. Experimental controller-performance results
in both steady state and transient operation are presented. The experimental
result showed a 4-5 % efficiency improvement with respect to pressure and NOx
constraints, compared to that of a single-injection controller
Related Publications
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2017). “Efficiency optimal,
maximum-pressure control in compression-ignition engines”. In: Proc. Ame-
rican Control Conf. (ACC 2017). Seattle, WA, USA, pp. 4753–4759.
34
1.5 Outline and Contributions
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2018). “Design and evalua-
tion of a multiple-injection controller for efficient fulfillment of NOx and
exhaust-temperature constraints”. submitted to SAE Int. J. Engines.
Chapter 11
This chapter introduces a model-based pressure-prediction method. The novelty
of this method lies in the use of the estimated heat-release rate to predict how the
cylinder pressure varies with injection timing. This is a computationally cheap
alternative to heat-release modeling. This chapter also presents how this me-
thod can be used to find efficiency-optimal injection timings with respect to con-
straints on maximum pressure and pressure-rise rate.
Related Publication
Ingesson, G., L. Yin, R. Johansson, and P. Tunestål (2015). “A model-based
injection-timing strategy for combustion-timing control”. SAE Int. J. Engines
8:3, pp. 1012–1020.
Chapter 12
This chapter presents how the pressure-prediction method in Chapter 11 can
be used in a model predictive control framework to efficiently fulfill constraints
on in-cylinder pressure, NOx and exhaust temperature with multiple injections.
This controller objective was previously studied in Chapter 10, but was revi-
sed here with a model-based approach. An alternative model predictive con-
trol formulation for pressure set-point tracking is also presented. The control-
lers presented rely on a heat-release detection method, capable of separating the
heat-release rate among several injections. Model-based methods for adding and
removing injections are also discussed.
Related Publication
Turesson, G., L. Yin, R. Johansson, and P. Tunestål (2018). “Predictive pressure
control with multiple injections”. submitted to E-CoSM 2018, Changchun,
China.
Chapter 13
The thesis is concluded in Chapter 13, where aspects on future research are also
discussed.
Authors Contributions
The author was the foremost contributor to the publications listed above. The
author was the main contributor to the work related to modeling, control and
experimental evaluation. This includes problem formulation, as well as develo-
ping, implementing, testing, and evaluating controller designs. The author wrote
35
Chapter 1. Introduction
the papers himself with input from the co-authors. All authors jointly determined
the general direction of the research, for example, that focus should be directed
towards control problems related to fuel-injection and its effect on the combus-
tion processes. Lianhao Yin was the main contributor to the design of the expe-
rimental platform used in the thesis. This includes design of the test-cell setup,
control-system and data-acquisition architecture.
36
2
Modeling
The increasing demands for reduced fuel consumption, emissions levels and im-
proved driveability lead to more actuators, sensors and complex control func-
tions. With the increasing engine complexity in mind, a systematic implementa-
tion of the engine-control system requires mathematical models for simulation,
calibration and controller design [Isermann, 2014, Atkinsson, 2009].
Model-based control is a controller-design approach where a mathematical
model describing physical or empirical system knowledge is utilized. A model
can be used off-line to evaluate controller performance in simulation, and in this
way provide suitable controller parameters and reduce the experimental work-
load. A model can also be used on-line to provide the controller with information
about predicted system behavior in real time. Both approaches to model-based
controller design were adopted in this work.
37
Chapter 2. Modeling
38
2.2 In-Cylinder Model
cycle k
| {z }
sampling of cycle k
BC for 0D model
injection → delay → combustion
BC for MVM
MV gas-exchange model
cycle-to-cycle controller
Figure 2.1 The modeling approach adopted in this work was to represent the
in-cylinder processes with a crank-angle resolved zero-dimensional (0D) model
with empirical expressions for chemical reaction rates. Cylinder boundary condi-
tions (BC) are determined by the states of the intake and exhaust manifolds, obtai-
ned from a mean-value gas-exchange model. Controller sampling and actuation
occur in-between consecutive engine cycles.
A mean-value model was used to model states of the intake and exhaust mani-
folds which determine the boundary conditions for the in-cylinder model.
These models will be described in greater detail in the following sections.
First, the crank-angle resolved in-cylinder pressure and temperature model
is presented. The mean-value gas-exchange model is presented subsequently.
Then, ignition-delay models of different complexity are presented and evaluated.
Finally, an NOx -formation model is presented. System identification was con-
ducted in order to obtain unknown model parameters from engine data. Identi-
fication experiments and results are also presented below.
39
Chapter 2. Modeling
dQ ht
dQ c
combustion chamber dW
piston
Figure 2.2 Open system boundary for the combustion chamber. The combus-
tion is modeled as a release of heat dQ c , whereas the cylinder gas performs work
on the piston dW , and heat dQ ht is transferred to the cylinder walls.
where dU is the change in system internal energy, dQ is the heat added to the
P
system, dW is the work done by the system, and i h i d m i is the enthalpy flux
across the system boundary. The combustion is modeled as a release of heat,
which gives that the heat added to the system dQ is the difference between rele-
ased chemical energy dQ c and heat transferred to the cylinder walls
dQ = dQ c − dQ ht (2.2)
With the piston work given by pdV , the first law can be rewritten as
X
dU = dQ c − dQ ht − pdV + h i d m i (2.3)
i
dU = mc v d T + ud m (2.4)
40
2.2 In-Cylinder Model
by assuming a constant R. Now, by inserting (2.4) and (2.6) into (2.3), and substi-
tuting
R
cv = (2.7)
γ−1
γ γ−1 X 1
dp = − pdV + (dQ c − dQ ht + h i d m i ) + (RT − (γ − 1)u)d m (2.8)
V V i V
Cylinder Geometry
The cylinder volume V was modeled as a slider-crank mechanism
µ r ¶
Vd π π
V = Vc + R v + 1 − cos( θ) − R v2 − sin2 ( θ) (2.10)
2 180 180
where Vd and Vc are displacement and clearance volumes, and R v is the ratio of
the connecting-rod length to the crank radius.
41
Chapter 2. Modeling
Temperature
The in-cylinder temperature T was computed from the ideal-gas equation using
the conditions at intake-valve closing (IVC)
pV TIVC
T= (2.11)
p IVC VIVC
where c vin and c vres are specific heats of inducted charge and residual gases, α is
the mass ratio between trapped residuals and inducted gases and χ is a measure
of residual-gas temperature decrease during gas exchange. The residual gas tem-
perature Tres was computed from the exhaust-valve-opening temperature of the
previous cycle.
1 ³ y´
Cx H y + x+ (O2 + 3.773N2 ) →
φn 4
µ ¶
y 1 ³ y´ 1 ³ y´
xCO2 + H2 O + − x+ O2 + 3.773 x + N2 (2.13)
2 φn 4 φn 4
where y/x is the fuel hydrogen to carbon ratio and φn is the molar equivalence
ratio.
Heat-Transfer Model
The convective heat-transfer rate from in-cylinder gas to piston, cylinder head
and walls dQ ht /d θ was modeled using Newton’s law of cooling
dQ ht hc A
= (T − T w ) (2.14)
dθ 60Nspeed
42
2.2 In-Cylinder Model
Vd TIVC
ω = C 1 S p +C 2 (p − p m ) (2.16)
p IVC VIVC
The first term in (2.16) relates gas motion to the mean piston speed S p , and the
second term captures the effect of charge-density variation during combustion,
where p m is the cylinder pressure of a motored cycle. The empirical parame-
ters C 1 and C 2 are engine dependent, where C 2 = 0 before the start of combus-
tion. The problem of estimating C 2 from cylinder-pressure data is investigated in
Chapter 4. Similar global heat-transfer models presented in [Annand, 1963; Ho-
henberg, 1979] are also commonly used to model (2.14).
Cylinder-Wall Temperature
The cylinder-wall surface temperature T w is determined by the heat flux from
in-cylinder gases to the engine coolant with temperature Tc . The cylinder wall,
cylinder head and piston were modeled as a single mass with conductive coeffi-
cient k c , thickness L c , mass m c and specific heat c p . By assuming that the inner
wall temperature is in steady state and that the outer-wall surface on the coo-
lant side has fixed temperature Tc , the following dynamic equation for T w can be
derived
d Tw 2A(h c + k c /L c ) 2Ah c
=− Tw + (T + Tc ) (2.17)
dθ m c c p 60Nspeed m c c p 60Nspeed
During the intake and exhaust strokes, T was assumed to be constant and equal
to intake- and exhaust-manifold temperatures. This model was previously pre-
sented and used in [Roelle et al., 2006; Widd et al., 2008].
Heat-Release Rate
The heat-release rate dQ c /d θ is difficult to model from first principles due
to its dependency on a multitude of factors such as chemical combus-
tion rates, fuel-injection profile and fuel-air mixing rates. The approach ad-
opted in Chapter 11 was to utilize the heat-release rate computed from
cylinder-pressure measurements to predict how p varies with fuel injection.
This method can not be used to simulate p. A Wiebe [1970] expression for the
43
Chapter 2. Modeling
accumulated heat-release
µ ¶
³ θ − θSOC b+1 ´
Q c (θ) 1 − exp − a for θ ≥ θSOC
= ∆θ (2.18)
Q tot
0 otherwise
was therefore used, mainly because of its simplicity, as an input to simulate (2.9)
in Chapter 10. In (2.18), θSOC is the start of combustion, and the parameters ∆θ,
a and b relate to the duration and shape of the heat-release profile.
44
2.3 Gas-Exchange System
air
intercooler θcool
1
θhot EGR
3 4 EGR
compressor
θLP
θHP
turbine
exhaust manifold
intake manifold
θBP
- valve
exhaust
Figure 2.3 Engine gas-exchange system layout. The engine was boosted by a
fixed-geometry turbocharger. Exhaust-gas recirculation was supplied by a high
and a low-pressure path, and the intake temperature was controlled by the gas
flow through a charge-air cooler prior to the intake manifold. The gas-exchange
system was modeled as 5 interconnected adiabatic volumes (indicated by the blue
regions). Gas flows through the EGR paths and the thermal management system
were controlled by actuating valve positions, denoted θx . The engine layout is fur-
ther described in Chapter 5.
45
Chapter 2. Modeling
then either goes through a charge-air cooler or directly to the intake manifold.
The mass-flow ratio between these paths is determined by the position of two
valves. These paths will later be referred to as the fast thermal-management
part (FTM) of the gas-exchange system. After leaving the exhaust manifold,
the exhaust gases either enter the cooled high-pressure EGR path or expands
through the turbine. After the turbine, some of the exhaust goes through a co-
oled low-pressure EGR path, and the remaining exhaust gas passes through a
back-pressure valve to the exhaust pipe. Mass flows through the EGR paths are
determined by the EGR-valve positions.
Gas-Exchange Dynamics
A gas-exchange system model was used in Chapter 9 to design a low-load PPC
controller. The main objective of the model was to describe in-cylinder tempe-
rature and mass in the low-load operating range of the engine, in order to com-
pute optimal valve-position actuation in simulation. The gas-exchange system
was modeled as five adiabatic ideal-gas control volumes, interconnected with
restrictions, denoted 1 to 5 in Fig. 2.3. Under this assumption, the dynamic equa-
tions with respect to pressure p i and temperature Ti in volume Vi are given by
[Eriksson and Nielsen, 2014]
d p i RTi p i d Ti
= (ṁ in − ṁ out ) +
dt Vi Ti d t
(2.22)
d Ti RTi
= (ṁ in c v,i (Tin − Ti ) + R(ṁ in Tin − ṁ out Ti ))
dt Vi p i c v,i
Here, ṁ in and ṁ out are in- and outgoing mass flows, and c v,i is the gas speci-
fic heat at constant volume. States in the intake and exhaust manifolds are also
denoted X in and X ex in this thesis. Mass flow from Vi to V j , ṁ i j , was mode-
led as turbulent compressible flow through a restriction, where ṁ i j is given by
[Heywood, 1988]
µ ¶ s µ µ ¶¶(γi −1)/γi
A i j p i p j 1/γi 2γi pj
ṁ i j =p 1− (2.23)
RTi p i γi − 1 pi
Here, A i j is the effective flow area and γi is the ratio of specific heats which is
different for air and exhaust. If the pressure ratio is too low
µ ¶ γi
pj 2 γi − 1
< (2.24)
pi γi + 1
the flow becomes choked, and the expression for ṁ i j is given by
µ ¶ γi + 1
A i j p i 1/2 2 2(γi − 1)
ṁ i j =p γ (2.25)
RTi i γi + 1
46
2.3 Gas-Exchange System
The effective flow area A i j was in the case of a valve restriction determined by
the valve angle θi j according to an empirical expression
³ ´
−k (θ −θ 0 )
A i j = A 0i j 1 − e Ai j i j i j (2.26)
in
The mass flow from the intake manifold to the cylinders ṁ cyl was modeled using
a volumetric efficiency η v , the intake-manifold density and the engine displace-
ment rate
in p 3 Vd Nspeed
ṁ cyl = ηv · · (2.27)
RT3 120
The volumetric efficiency was approximated by assuming displacement of new
charge by residual gas
µ ¶1/γ
p4
rc −
p3
ηv = (2.28)
rc − 1
where r c is the engine compression ratio. The cylinder-out mass flow was then
obtained by adding the fuel flow ṁ f
out in
ṁ cyl = ṁ cyl + ṁ f (2.29)
where ṁ 45 is the mass flow over the turbine, modeled using (2.23). The compres-
sor power P c was then related to P t through a mechanical efficiency η m and the
static relation
P c = ηm P t (2.32)
The compressor power P c was in turn related to the power of isentropic compres-
sion over the compressor P cs
P cs = η c P c (2.33)
where P cs is a function of the compressor mass flow ṁ 12
õ ¶ !
s p 2 1−1/γ1
P c = c p,1 T1 − 1 ṁ 12 (2.34)
p1
47
Chapter 2. Modeling
η tc = η c η m η t (2.36)
where p im , Tim are measured pressure and temperature data, and p i , Ti are cor-
responding model outputs. The importance of capturing different outputs are
determined by the tuning parameters αi and βi . Cost function (2.39) was mini-
mized subject to
η tc ≤ 1 (2.40)
and measured boundary conditions p atm , Tatm . Model output and measured
data are shown in Fig. 2.4 for a local minimizer of (2.39) ϑ∗ , where ϑ∗ was
found using the MATLAB nonlinear-optimization toolbox. The data presented
in Fig. 2.4 are combined sequences of steady-state data with different θCool and
θHP positions for p IMEP between 1 and 5 bar. The tuning parameters were set to
prioritize model fit with respect to the intake-manifold state. In Fig. 2.4, it can be
seen that the model managed to predict in-cylinder temperature at the start of
injection TθSOI and the relative air-fuel ratio λ, which was the main purpose of
the model when used in Chapter 9.
48
2.3 Gas-Exchange System
1.2 1.6
1.4
p ex [bar]
p in [bar]
1.1
1.2
1
1
0.9
0 500 1,000 1,500 2,000 0 500 1,000 1,500 2,000
cycle cycle
340 800
TθSOI [K]
320 750
Tin [K]
300 700
280 650
0 500 1,000 1,500 2,000 0 500 1,000 1,500 2,000
cycle cycle
1.4
1.3 40
θ HP [deg]
p 2 [bar]
1.2
20
1.1
1 0
0 500 1,000 1,500 2,000 0 500 1,000 1,500 2,000
cycle cycle
100 20
80
θ cool [deg]
15
60
λ [-]
10
40
20 5
0 0
0 500 1,000 1,500 2,000 0 500 1,000 1,500 2,000
cycle cycle
Figure 2.4 Model output (black) and measured data (blue) for a local minimizer
ϑ∗ of (2.39). The light blue color indicates cycle-to-cycle variation in the intake
and exhaust manifolds. The engine data presented are combined sequences of
steady-state data with different θCool and θEGR positions for p IMEPn from 1 to 5
bar. The model managed to capture in-cylinder temperature at the start of injec-
tion TθSOI and the relative air-fuel ratio λ.
49
Chapter 2. Modeling
d p in RTin in p in d Tin
= (ṁ air + ṁ EGR − ṁ cyl )+ (2.42)
dt Vin Tin d t
p ex VEVO
m res = (2.43)
RTex
This method was used to compute [O2 ]IVC online in Chapter 7, where [O2 ]IVC was
used to predict ignition-delay variation.
θ10 − θSOI
τ= (2.44)
0.006Nspeed
50
2.4 Ignition-Delay Modeling
6,000 p
Qc
Q c [J] 4,000
2,000
Figure 2.5 The ignition delay τ was defined as the time between the start of in-
jection θSOI and the crank angle of 10 % heat released θ10 . This indicator for the
start of combustion was computed from the accumulated heat release Q c . The
start of injection was here determined by the rising flank of the injector-current
pulse and a hydraulic injector delay.
51
Chapter 2. Modeling
M1
The first model is given by the following inverted reaction-rate expression
α β
τ = A[O2 ] [Cx Hy ] p̄ ζ S δp e E a /R̃ T̄ (2.48)
52
2.4 Ignition-Delay Modeling
M2
The second model is a crank-angle resolved extension of M1 where it is assumed
that combustion starts when
Z θSOI +τ
A[O2 ]α [Cx Hy ]β p ζ S δp e −E a /R̃T d θ = 1 (2.50)
θSOI
M3
The third model is a two-stage reaction model with empirically determined reac-
tion rates as presented in [Westbrook and Dryer, 1981]
¡x
y ¢ k1 y
Cx Hy + O2 → xCO + H2 O
+
2 4 2
(2.52)
1 k2
CO + O2 → CO2
2
Here, the fuel Cx Hy reacts with O2 to form CO and H2 O in the first reaction,
CO is then oxidized to CO2 in the second reaction. Reaction rates are given by the
following differential equations
d [Cx Hy ]
= −k 1 [Cx Hy ]β [O2 ]α
dt
d [O2 ] ¡x y ¢
= − + k 1 [Cx Hy ]β [O2 ]α − k 2 [CO][O2 ]1/2
dt 2 4
(2.53)
d [CO]
= xk 1 [Cx Hy ]β [O2 ]α − k 2 [CO][O2 ]1/2
dt
d [CO2 ]
= k 2 [CO][O2 ]1/2
dt
Initial conditions for (2.53) are given by the global cylinder oxygen and fuel
concentrations after fuel injection. The reaction-rate parameters k i are empirical
expressions on the form
i
k i = A i S pi e −E a /R̃T
δ
(2.54)
Finally, the heat-release rate of the reactions is computed from
dQ c
= V (Q 1 k 1 [Cx Hy ]β [O2 ]α + k 2Q 2 [CO][O2 ]1/2 ) (2.55)
dt
53
Chapter 2. Modeling
Table 2.1 Investigated operating points for model calibration and evaluation.
where Q 1 and Q 2 are the lower heating values per mole of fuel and CO in (2.53).
Combustion is assumed to start when Q c has reached 10 % of the expected total
heat released Q ctot . For M3, unknown model parameters are given by
ϑM 3 = [A 1 , A 2 , E a1 , E a2 , α, β, δ1 , δ2 ]T (2.56)
Ignition-Delay Experiments
Ignition-delay experiments were conducted at the four operating points obtai-
ned by combining p IMEPg = 5, 10 bar and Nspeed = 1200, 1500 rpm, see Table 2.1.
Suitable θSOI were found at each operating point with EGR valves closed and
both thermal-management valves opened at 45◦ . A layout of the experimen-
tal engine setup showing valve locations was presented in the section covering
gas-exchange modeling above, see Fig. 2.3.
In order to investigate the τ response to different engine inputs, θSOI ,
thermal-management (θcool/hot ) and EGR-valve positions (θHP/LP ) were varied
manually in steps during a total of 12000 cycles at each operating point.
The thermal-management valve positions were changed by setting
54
2.4 Ignition-Delay Modeling
3
τID τVAL 40
Tin [◦ C]
τ [ms]
30
2
20
1 10
0 0.5 1 0 0.5 1
cycle [-] ·104 cycle [-] ·104
−10
8
−20
−30 7
0 0.5 1 0 0.5 1
cycle [-] ·104 cycle [-] ·104
Figure 2.6 Experiment data at operating point L1S1, top left: τ [ms], top right:
Tin [C ], bottom left: θSOI [CAD ATDC] and bottom right: [O2 ]IVC [mol/m3 ].
Parameter Identification
Sum-of-squares model-error cost functions J M x (τID , ϑM x ) were minimized for
each model with respect to ϑM 1−3 and the identification data set τID to identify
suitable model parameters.
N ¡
X ¢2
J M 1 (τID , ϑM 1 ) = ln(τID M1
i ) − ln(τi ) (2.58)
i =1
with respect to ϑM 1 is a linear regression, where the index i denotes sample num-
ber and N is the number of samples. In the following notation
¡ ¢T
y = ln(τ1M 1 ) ... ln(τM 1
N ) (2.59)
55
Chapter 2. Modeling
Tin [◦ C]
τ [ms]
1.2
40
1
0.8 20
0 0.5 1 0 0.5 1
cycle [-] ·104 cycle [-] ·104
11
−5
10
−10
−15 9
0 0.5 1 0 0.5 1
cycle [-] ·104 cycle [-] ·104
1
1 R̃ T̄1
ln([O2 ]1 ) ln([Cx Hy ]1 ) ln(p̄ 1 ) ln(S p 1 )
. .. .. .. .. ..
Φ = .. . . . . . (2.60)
1
1 ln([O2 ]N ) ln([Cx Hy ]N ) ln(p̄ N ) ln(S p N )
R̃ T̄ N
¡ ¢T
the minimizer of (2.58), ϑ∗ = ln(A ∗ ) E a∗ α∗ β∗ ζ∗ δ∗ , is given by
56
2.4 Ignition-Delay Modeling
Tin [◦ C]
τ [ms]
2 30
1.5 20
1 10
0 0.5 1 0 0.5 1
cycle [-] ·104 cycle [-] ·104
9
−10
8
−20
7
0 0.5 1 0 0.5 1
cycle [-] ·104 cycle [-] ·104
This is not a linear regression, without a closed-form expression for the minimi-
zer. The cost function was therefore minimized numerically, using the MATLAB
Optimization Toolbox. In order to avoid local minima, (2.62) was minimized with
respect to different initial parameters.
M3 The reaction-rate parameters of the third model were found using the same
numerical procedure as for M2. The chemical reactions were simulated by app-
roximating the derivatives in (2.53) using the forward-Euler method and a suffi-
ciently small step size, h = 0.01 ms. The model parameters were then found by
minimizing
N ³
X ´2
J M 3 (τID , ϑM 3 ) = M 3 ID
Q c,i (θ10,i ) − 0.1m f ,i Q LHV (2.63)
i =1
M 3 VAL
where Q c,i (θ10,i ) is the modeled accumulated heat-release at θ10 in the
identification-data set, and m f ,i Q LHV is the injected fuel energy, computed from
fuel-flow measurements.
57
Chapter 2. Modeling
Tin [◦ C]
τ [ms]
60
1
40
0.8
20
0 0.5 1 0 0.5 1
cycle [-] ·104 cycle [-] ·104
−20 10
−40 8
−60 6
0 0.5 1 0 0.5 1
cycle [-] ·104 cycle [-] ·104
Figure 2.9 Experiment data at operating point L2S2. Unfortunately, the last 2000
cycles of τID were corrupted due to a malfunctioning thermal-management valve.
This part of τID was therefore replaced with the last 2000 cycles of τVAL . Note the
unintentional decrease in θSOI after cycle 11000. It was decided to be keep these
data points in the sets τID and τVAL .
Model Evaluation
The three models were evaluated by how well they could explain variation in the
validation-data set τVAL in two different ways:
• By how well a model calibrated by τID from one operating point could pre-
dict τVAL from the same operating point. During model calibration, the
speed dependence was removed (δ = 0). The load dependence was also
removed by setting β = 0 in M1-2 and β = 1 in M3. Furthermore, M1 and
M2 were investigated with and without pressure dependence, i.e., with ζ
free and ζ = 0.
• By how well a model calibrated with the complete τID data set could pre-
dict the complete τVAL data set. Now, β and δ were free parameters during
model calibration.
58
2.4 Ignition-Delay Modeling
L1S1 L1S2
0.15 0.15
RMSE [ms] without p with p
RMSE [ms]
0.1 0.1
0.05 0.05
M1 M2 M3 M1 M2 M3
L2S1 L2S2
0.15 0.15
RMSE [ms]
RMSE [ms]
0.1 0.1
0.05 0.05
M1 M2 M3 M1 M2 M3
Figure 2.10 Model RMSE with respect to validation data for M1, M2 and M3 at
each operating point. Including p dependence in M1 and M2 did not give any
improvement. Increasing model complexity from M1 to M2 and M3 did not ne-
cessarily yield a smaller prediction error.
The reason for evaluating the models with respect to these two aspects was to de-
termine if there is an incentive to use multiple models with fewer parameters, in-
stead of having one model that covers all operating points with more parameters.
Model performance was evaluated by computing the root-mean-square error
v
u N
u1 X
RMSEM x = t (τM x − τVAL )2 (2.64)
N i =1 i i
where τM x is the model output and τVAL is measured τ in the validation-data set.
59
Chapter 2. Modeling
L2S2 but not at L1S1 and L1S2. Increasing model complexity from M2 to M3 im-
proved model performance in all operating points but L2S2.
For a more detailed analysis, the model outputs are displayed together with
τVAL in Figs. 2.11-2.14, where the gray dots indicate τVAL and the red, blue and
green lines are model outputs from M1, M2 and M3 without pressure depen-
dence.
In Fig. 2.11, the model outputs did not follow τVAL well around cycles 100 and
1000. This occured when θSOI was delayed at cycles 6100 and 7000 in Fig. 2.8. The
same behavior is found in Fig. 2.12 from cycle 1 to 800, in Fig. 2.13 at cycle 200,
and in Fig. 2.14 at cycles 500 and 1500. As θSOI was delayed close to TDC, the
computed T during τ increased which caused the models to predict a decreased
τ. Instead, the measured τVAL increased. Hence, the models were incapable of
anticipating the point where a delayed θSOI started to increase τ. This behavior
can clearly be seen in Fig. 2.12 at cycle 1400, where an increase in θSOI resulted in
a large increase in τVAL . The models were overall better at predicting τVAL during
variations in TIVC , [O2 ]IVC and θSOI for early θSOI , see cycle 750 in Fig. 2.11.
Similar model outputs for M1, M2 and M3 in Figs. 2.11-2.14 indicate that
controller designs based on the different models would yield comparable per-
formance.
60
2.4 Ignition-Delay Modeling
2.6
2.4
2.2
1.8
1.6
1.4
0 1,000 2,000 3,000 4,000 5,000 6,000
cycle [-]
Figure 2.11 Modeled τ together with τVAL from operating point L1S1. It can be
noted that the models did not follow τVAL well at cycles 100 and 1000. M1 had
overall comparable performance to M2 and M3.
∂τ(X 0 )
τ(k + 1) = τ(k) + ∆TIVC (k)
∂TIVC
∂τ(X 0 ) ∂τ(X 0 )
+ ∆[O2 ]IVC (k) + ∆θSOI (k) (2.66)
∂[O2 ]IVC ∂θSOI
where ∆ is the forward-difference operator and k denotes cycle index. The li-
near model could then easily be incorporated in controller design frameworks
61
Chapter 2. Modeling
1.3
1.2
1.1
0.9
0 1,000 2,000 3,000 4,000 5,000 6,000
cycle [-]
Figure 2.12 Modeled τ together with τVAL from operating point L2S1. The model
mismatch at cycle 1400 is an example of where θSOI was increased and the models
underestimated τVAL . The models were overall better at predicting τVAL during
variations in TIVC , [O2 ]IVC and θSOI for early θSOI .
such as linear quadratic Gaussian (LQG) control [Kalman, 1960] or linear model
predictive control [Maciejowski, 2002]. Linearization of M1, M2 and M3 can be
conducted by approximating the partial derivatives numerically
62
2.4 Ignition-Delay Modeling
2.2
1.8
1.6
1.4
1.2
0 1,000 2,000 3,000 4,000 5,000 6,000
cycle [-]
Figure 2.13 Modeled τ together with τVAL from operating point L1S2.
Discussion
The models could not accurately detect when an increase in θSOI started to give
an increase in τ. If this was caused by errors in the estimated cylinder-gas state,
errors in the heat-release analysis or by the fact that the models did not account
for spatial and physical effects is not known. Model performance was however
more acceptable when θSOI was advanced earlier during the compression stroke
and when TIVC and [O2 ]IVC were varied. This implies that it is easier to obtain sa-
tisfactory controller performance when regulating τ using TIVC and [O2 ]IVC , un-
less θSOI is sufficiently far away from TDC. When evaluating performance for all
operating points simultaneously (see Fig. 2.16), model performance decreased
with complexity. Small improvements could be obtained when using four local
models instead of one global model.
The results were somewhat unambiguous when comparing model perfor-
mance, and it was only worthwhile to increase model complexity in some ca-
ses, (see L1S1, L2S1, Fig. 2.10). When studying Figs. 2.11-2.14, it can be seen
that model performance did not differ significantly in most cases. When taking
into account for the computational cost of linearization, M1 was superior to M2
and M3 (see Table 2.2). Furthermore, M1 had a closed-form expression (2.61) for
63
Chapter 2. Modeling
1.4
τ [ms] Cross-Validation Data
1.2
0.8
Figure 2.14 Modeled τ together with τVAL from operating point L2S2. Note that
the models were able to predict the steep increase in τ as θSOI was decreased
around cycle 5250. The reason for this was that cycles 4000-6000 at L2S2 were
included τID , see Fig. 2.9.
parameter calibration. Based on these findings, M1 was the suitable choice for
model-based controller design, and was therefore used by the model predictive
controller in Chapter 7.
These models did not include effects from cylinder-wall temperature, fuel va-
porization, atomization and fuel-spray interaction with the combustion cham-
ber walls. If carefully modeled, these effects are believed to improve the model
performance. However, information from these effects were not easily accessible
from the sensors available, making the validation of these effects difficult.
M4
The ignition-delay models presented above did not take fuel reactivity into ac-
count. One of the main motivators for closed-loop combustion control is the
possibility to handle fuel variation. The octane-number dependent τ correlation
64
2.4 Ignition-Delay Modeling
0.09
0.08
RMSE [ms]
0.07
0.06
0.05
M1 M2 M3
Figure 2.15 RMSE with respect to validation data. Increasing the model comple-
xity from M1 to M2 and M3 increased RMSE, and using individual operating-point
models gave a slight reduction in RMSE. Increased prediction performance could
therefore be obtained by switching between different models instead of using one
model for all operating points.
was therefore used in Chapter 6 to investigate how the fuel reactivity affects con-
troller design. This model was developed for primary-reference fuel (PRF) blends
from PRF0 (n-heptane) to PRF100 (iso-octane) and was calibrated with data from
constant-volume simulations in Cantera, using a reduced gasoline-surrogate ki-
netic mechanism. Note that the PRF value and octane number are equal per de-
finition. In (2.68), φ denotes equivalence ratio, p denotes pressure, and xO 2 is the
mole fraction of oxygen of the inducted gas. The parameters α, β, ζ and Λ are
temperature-dependent polynomials multiplied with a PRF-dependent expo-
nential expression, yielding roll-off at low temperatures. Delvescovo et al. [2016]
calibrated the polynomial coefficients with respect to simulated τ data at the fol-
lowing conditions
65
Chapter 2. Modeling
M2
3
2.5
1.5
0.5
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4
cycle [-] ·104
Figure 2.16 Modeled τ together with τVAL for all operating points.
to capture a wide range of engine-like operating conditions. With (2.68), the start
of combustion θSOC was computed using a Livengood-Wu integration criterion
[Livengood and Wu, 1955]
Z θSOC 1
dt = 1 (2.69)
θSOI τ
66
2.5 NOx -formation Model
is the two-zone (burned and unburned gas) combustion model coupled with the
Zeldovich mechanism
O + N2 → NO + N
(2.70)
N + O2 → NO + O
that was previously developed for combustion-control purposes in [Egnell, 2001;
Murić et al., 2014]. This model assumes NOx formation through the thermal-
formation path alone. Less dominant formation paths include reactions between
nitrogen and unburnt hydrocarbons and reactions with nitrogen components
in the fuel [Fenimore, 1971]. More detailed NOx -formation models incorporate
fuel-air mixing and the cylinder air-to-fuel ratio distribution, which affect the
combustion temperature. An extensive NOx -modeling review can be found in
[Miller and Bowman, 1989].
The model used here separates the cylinder gases into a burned zone and
an unburned zone. The unburned zone contains fuel, air and EGR, whilst the
burned zone contains combustion products. When fuel and air react at a lo-
cal air-fuel ratio, λlocal , the reaction products are moved from the unburned
zone to the burned zone where combustion chemistry is assumed to be in
equilibrium after combustion. The local air-fuel ratio λlocal was here used as a
model-calibration parameter. Under these assumptions, the mass m bz and tem-
perature Tbz of the burned zone change every crank-angle increment ∆θ accor-
ding to
1 dQ c (θ)
m bz (θ + ∆θ) = m bz (θ) + (1 + λlocal + r EGR )∆θ
Q LHV d θ
µ ¶γ−1 (2.71)
p(θ + ∆θ) γ 1 dQ c (θ)
Tbz (θ + ∆θ) = Tbz (θ) + ∆θ
p(θ) m bz (θ)c p d θ
where dQ c /d θ is the heat-release rate, and r EGR is the EGR ratio. The
unburned-zone temperature Tuz changes due to isentropic compression and
expansion alone
µ ¶γ−1
p(θ + ∆θ) γ
Tuz (θ + ∆θ) = Tuz (θ) (2.72)
p(θ)
The NO-formation rate in the burned zone is given by
³ ³ [NO] ´2 ´
2r 1 1 −
d [NO] [NO]e [NO] dV
= − (2.73)
dt [NO] r 1 V dt
1+
[NO]e r 2
according to the Zeldovich mechanism [Zeldovich et al., 1947]. This expression
assumes radical species in equilibrium where NO and O equilibrium concentra-
tions, [NO]e and [O]e , can be computed from combustion-product composition,
67
Chapter 2. Modeling
Model Calibration
The NOx -model was calibrated with respect to experimental data. Figure 2.17
shows modeled NOx coincidence with measured NOx concentrations, obtained
from engine experiments. The data originate from operating points at p IMEP =
5 and 10 bar with 0 and 30 % EGR, with combustion timings from 0 to 15 CAD,
pilot-injection durations from 0 to 0.4 ms, and with diesel (×) and gasoline (◦)
fuel.
In Fig. 2.17, performance was found to degrade for low NOx concentrations,
especially with high EGR ratios where the model underestimated the measured
NOx emissions. An explanation for this could be that other formation paths be-
come more important at these operating points. Possible improvement could be
obtained by also including the N2 O formation path as described in [Gong and
Rutland, 2013], which improves NOx prediction at lower combustion tempera-
tures. It can also be seen that model errors were larger for gasoline compared to
diesel. The reason for this is unknown and deserves further investigation.
68
2.5 NOx -formation Model
1,500
1,000
500
100
, × EGR = 0 %, Load = 5 bar
, × EGR = 0 %, Load = 10 bar
75
, × EGR = 30 %, Load = 5 bar
, × EGR = 30 %, Load = 10 bar
Error [%]
50
25
0
500 1,000 1,500 2,000 2,500
Measured NOx [ppm]
Figure 2.17 The NOx model was calibrated at p IMEP = 5 and 10 bar and 1200
rpm, with 0 and 30 % EGR, combustion timings from 0 to 15 CAD, pilot-injection
durations from 0 to 0.4 ms, and with diesel (×) and gasoline (◦) fuel.
69
3
Control and Estimation
Methods
This chapter gives a brief overview of the control and state-estimation methods
used in the thesis.
70
3.1 Control Concepts
which consists of output costs J y (k) over a prediction horizon H p and input costs
J u (k) over a control horizon Hu . The output cost at sample k
p
X
J y (y(k)) = ω y j (y rj (k) − y j (k))2 (3.4)
j =1
is the sum of the squared set-point deviations y rj (k) − y j (k) for each output y j
and corresponding set point y rj , scaled with a cost weight ω y j . The input cost at
sample k
Xm
J u (u(k)) = ωu j u j (k)2 + ω∆u j ∆u j (k)2 (3.5)
j =1
is the square sum of both absolute input values u j and changes ∆u j with corres-
ponding weights ωu j and ω∆u j . The positive cost weights ωx are controller design
parameters that determine how the controller prioritizes different output errors
and control actions.
Feedback is introduced by minimizing (3.3) subject to the system initial con-
ditions x(0) = x 0 , obtained from measurement or state estimation. The main fea-
ture of MPC is its ability to incorporate constraints with respect to inputs, outputs
and states in the optimization problem
y y
−²η min + y min ≤ y(k) ≤ y max + ²η max
−²ηxmin + x min ≤ x(k) ≤ x max + ²ηxmax
(3.6)
−²ηumin + u min ≤ u(k) ≤ u max + ²ηumax
−²η∆u ∆u
min + ∆u min ≤ ∆u(k) ≤ ∆u max + ²η max
71
Chapter 3. Control and Estimation Methods
yr
k =0 k = Hp
The positive variable ² and vectors η determine costs for constraint violation.
When η = 0, the solution U ∗ is not allowed to violate the constraint, but when
η 6= 0, the constraint can be violated with additional cost, determined by the cost
weight ρ ² in (3.3). The design choice of having η > 0 was used in Chapter 8 to en-
sure existence of feasible solutions. The principle of MPC is illustrated in Fig. 3.1
for a single-input / single-output system where the objective is to track a set point
y r , subject to constraints on the input (dashed).
If the system dynamics are linear, the optimization problem of minimizing
(3.3) subject to the equality constraints in (3.2), and the inequality constraint in
(3.6) is a quadratic-programming (QP) problem
1 T
minimize x Hx + f T x (3.7)
x 2
subject to Ax ≤ b
A eq x = b eq
72
3.1 Control Concepts
Linearization
Linear model predictive control requires linear system models. Most of the mo-
dels presented in Chapter 2 are however nonlinear. In order to enable linear con-
troller design, one can linearize the system dynamics at the current operating
point. Given a nonlinear discrete-time model on the form
∂ fi ∂ fi
Ai j = (x 0 , u 0 ), Bi j = (x 0 , u 0 )
∂x j ∂u j
(3.10)
∂h i
Ci j = (x 0 , u 0 )
∂x j
where ∆x(k), ∆u(k) and ∆y(k) are deviations from the linearization point x 0 , u 0 ,
y 0 . Linearization was used in Chapter 7 to linearize an ignition-delay model, and
73
Chapter 3. Control and Estimation Methods
Here, the notation x ∼ p(x|y) indicates that x is distributed according to the con-
ditional probability-density function p of x, given y.
and perturbed with additive Gaussian noise with zero mean and covariance ma-
trices Q and R
v (k) ∼ N (0,Q)
(3.14)
e(k) ∼ N (0, R)
the estimation problem has an analytic solution, and the Kalman Filter (KF)
provides an optimal algorithm for iteratively estimating x(k), see Algorithm 1
[Kalman, 1960]. The Kalman filter is used in Chapter 7 to estimate the EGR mass
flow for [O2 ]IVC computation, and in Chapter 8 to filter cycle-to-cycle variation in
pressure-rise rate and combustion timing.
If the system (3.11) is nonlinear on the form
74
3.2 State Estimation
2: while k > 0 do
3: x̂(k|k − 1) = A x̂(k − 1|k − 1) + B u(k − 1)
4: P (k|k − 1) = AP (k − 1|k − 1)A T +Q
5: e(k) = y(k) −C x̂(k|k − 1)
6: S(k) = C P (k|k − 1)C T + R
7: K (k) = P (k|k − 1)C T S −1 (k)
8: x̂(k|k) = x̂(k|k − 1) + K (k)e(k)
9: P (k|k) = (I − K (k)C )P (k|k − 1)
10: end while
the Kalman filter is unfortunately not directly applicable. There are however
well-established methods for solving nonlinear estimation problems. The exten-
ded Kalman filter (EKF) and the particle filter (PF) are examples of such methods.
They were used in Chapter 4 for automatic calibration of the heat-release model
presented in Chapter 2.
2: while k > 0 do
3: x̂(k|k − 1) = f (x̂(k − 1|k − 1)|u(k − 1))
∂f
4: A(k − 1) = ∂x |x̂(k|k−1),u(k−1)
5: P (k|k − 1) = A(k − 1)P (k − 1|k − 1)A(k − 1)T +Q
6: e(k) = y(k) − h(x̂(k|k − 1))
∂h
7: C (k) = ∂x |x̂(k|k−1)
8: S(k) = C (k)P (k|k − 1)C T (k) + R
9: K (k) = P (k|k − 1)C (k)T S −1 (k)
10: x̂(k|k) = x̂(k|k − 1) + K (k)e(k)
11: P (k|k) = (I − K (k)C (k))P (k|k − 1)
12: end while
75
Chapter 3. Control and Estimation Methods
Np
X
p(x(k)|y(k)) ≈ ωi (k)δ(x(k) − x i (k)) (3.16)
i =1
where δ(·) is the Dirac delta function. The PF performs the sampling proce-
dure using a sequential Monte Carlo technique where particles x i (k + 1) are
sampled sequentially, given the old particles x i (k), and a proposal distribution
q(x(k + 1)|x i (k)). After each time step, the weights are updated to represent the
desired probability density function in (3.12). When choosing q(x(k + 1)|x i (k)) =
p(x(k + 1)|x i (k), u(k)), the update rule for the weights becomes ωi (k + 1) =
ωi (k)p(y(k + 1)|x i (k + 1)). To avoid particle depletion, meaning that only a few
weights contribute to (3.16), the particles have to be resampled according to the
weight distribution. This procedure puts more particles into areas of high proba-
bility and discards particles in regions of low probability. The resampling step is
commonly conducted when the ratio of the number of effective particles Neff
1
Neff = PN (3.17)
p
i =1
(ωi )2
76
4
Heat-Release Analysis
4.1 Introduction
Heat-release analysis refers to the use of physical models, such as the pres-
sure model introduced in Chapter 2, to infer information about the combustion
processes from cylinder-pressure measurements. This information is commonly
used for engine diagnostics, research, control and simulation. Pioneering work
on methods for inferring the heat-release rate from measured in-cylinder pres-
sure were presented in [Rassweiler and Withrow, 1938; Krieger and Borman, 1966;
Gatowski et al., 1984]. From a feedback-control perspective, heat-release analy-
sis provides the possibility to regulate combustion timing and ignition delay on
a cycle-to-cycle basis.
Heat-Release Analysis
In this work, the heat-release rate dQ c /d θ is computed from the measured pres-
sure p by rearranging (2.9)
dQ c γ dV 1 d p dQ ht
= p + V + (4.1)
dθ γ − 1 dθ γ − 1 dθ dθ
Crank angles of x% burned θx are commonly used indicators for the timing and
duration of the combustion process. These are obtained from the accumulated
heat release Z θ
dQ c
Q c (θ) = dθ (4.2)
θIVC d θ
77
Chapter 4. Heat-Release Analysis
6,000 p
Qc
4,000
Q c [J]
2,000
θ50
0
dQ c X M dQ i
c
= (4.4)
dθ i =1 d θ
The procedure used for obtaining dQ ci /d θ from dQ c /d θ was to first detect the
M most significant peaks, as described above. With the peaks detected, dQ c /d θ
78
4.1 Introduction
Model Parameters
The heat-release model in (4.1) has a set of unknown parameters that have to
be tuned for satisfactory performance. These parameters can be tuned man-
ually with knowledge of the appearance of physical heat-release rates and the
influence of the different model components. This procedure is, however, time
consuming and has to be redone from time to time.
The development of automatic calibration methods has been an ac-
tive research area during the past decades, where methods for calibrating
pressure-sensor offset [Tunestål et al., 2001; Brunt and Pond, 1997], polytro-
pic coefficients [Manente et al., 2008; Randolph, 1990], volume-curve offset
[Stas, 2004; Tunestål, 2001] and compression ratio [Klein et al., 2006] have been
presented. In [Klein, 2007], an off-line method for calibration of a large set of
parameters simultaneously was presented and studied in detail. It was, however,
concluded by Eriksson [1998], that all model parameters might not be identifia-
ble simultaneously. This indicates that a calibration problem involving a large
set of parameters is not easily solved.
This chapter investigates on-line calibration of a subset of the model parame-
ters in (4.1). The task is first formulated as a nonlinear estimation problem, where
unknown states of a dynamic system are to be estimated given a statistical mo-
del and sensor measurements. The formulated estimation problem is then sol-
ved using the extended Kalman filter (EKF) and the bootstrap particle filter (PF),
79
Chapter 4. Heat-Release Analysis
600 600
dQ c /d θ [J/CAD] 1. 2.
dQ c /d θ [J/CAD]
400 400
200 200
0 0
−20 0 20 −20 0 20
θ [CAD] θ [CAD]
600 600
3. 4.
dQ c /d θ [J/CAD]
dQ c /d θ [J/CAD]
400 400
200 200
0 0
−20 0 20 −20 0 20
θ [CAD] θ [CAD]
described in Chapter 3. These filters lend themselves nicely for real-time appli-
cations because of their sequential processing of measured data. The problem
representation provides a general framework for which any parameter combina-
tion could be estimated as long as the system is observable.
The chapter is outlined as follows: The estimation-problem formulation and
filter configurations are introduced in Secs. 4.2 to 4.3. Filter-performance results
with respect to simulated and experimental data are then given in Secs. 4.4 to 4.5.
Discussion and conclusions are presented in Secs. 4.6 and 4.7.
80
4.2 Problem Formulation
dQ c γ dV 1 d p dQ ht
= p + V + (4.7)
dθ γ − 1 dθ γ − 1 dθ dθ
where ²1 ∼ N (0, σ21 ) and ²2 ∼ N (0, σ22 ) are i.i.d normally distributed noise proces-
ses, with standard deviations σi . They are introduced to represent sensor noise
and unmodeled effects. The mean injected fuel energy Q ctot is assumed to be
known and determined by the injection duration and common-rail pressure. The
variable d ∼ N (0, σ2d ) is a random offset accounting for stochastic variation in the
injected fuel energy.
¡ 0 0
¢T
C2 θ∆TDC T w0 ∼ N (µ0 , Σ0 ) (4.9)
81
Chapter 4. Heat-Release Analysis
6,000
Qc θ EOC
Q ctot
4,000
Q c [J]
2,000
θ SOC
Figure 4.3 A realization of the statistical model in (4.9-4.13). The dashed line re-
presents the expected accumulated heat release before and after combustion wi-
thout stochastic variation. Deviation from this value is due to an offset d and noise
processes ²1 , and ²2 for which σ21 = 0, σ22 = 2500. It is assumed that the injected
fuel energy Q ctot , the start and end of combustion θSOC and θEOC are available.
dQ c γ dV 1 d p dQ ht
= p + V +
dθ γ − 1 dθ γ − 1 dθ dθ
Z (4.13)
θ dQ c
Q c (θ) = dθ
θIVC dθ
82
4.2 Problem Formulation
where a realization of (4.13) with the correct model parameters is shown in Fig.
4.3. Equations (4.9)-(4.13) are on the form
Observability
It is crucial to evaluate the system observability when conducting state-estimation.
Observability is a measure of how well the system state x(k) can be inferred
from knowledge of the system output y(k), and was introduced to linear-system
theory by Kalman [1959]. In order to apply the concept of observability to the
system in (4.10-4.13), a linearization
was conducted at the operating point presented in Table 4.1. In our case, observ-
ability was investigated by evaluating
∂h i (x(k))
C (k)i j = (4.16)
∂x j
83
Chapter 4. Heat-Release Analysis
System Observability
150
∂h/∂θ∆TDC
100 ∂h/∂C 2 × 0.001
∂h/∂T w × 100
50
0
C
−50
−100
−150
−200
−150 −100 −50 0 50 100 150
θ [CAD]
84
4.4 Simulation Results
Convergence
Filter consistency and rate of convergence were evaluated by initializing the EKF
and the PF with correct model parameters according to Table 4.1, apart from the
incorrect initial filter states
¡ ¢T
x 10 = 1, 0.0095, 595
¡ ¢T (4.20)
x 20 = −1, 0.002, 336
85
Chapter 4. Heat-Release Analysis
Table 4.1 Nominal operating point for filter evaluation. The heat-release rate,
dQ c /d θ, was chosen as a Gaussian function with a standard deviation of 5 CAD.
86
4.4 Simulation Results
EKF PF
1 1
θ ∆ TDC [CAD]
θ ∆ TDC [CAD]
0.5 0.5
0 0
−0.5 −0.5
−1 −1
0 20 40 60 80 0 50 100 150 200
cycle [-] cycle [-]
·10−2 ·10−2
1 1
0.8 0.8
C 2 [m/(sK)]
C 2 [m/(sK)]
0.6 0.6
0.4 0.4
0.2 0.2
0 0
20 40 60 80 50 100 150 200
cycle [-] cycle [-]
600 600
T w [K]
T w [K]
400 400
200 200
20 40 60 80 50 100 150 200
cycle [-] cycle [-]
Figure 4.5 Filter convergence with initial conditions x 01 (dashed) and x 02 (solid).
True state values are indicated in red. Note the different scales on the cycle-axes
for the two filters. It can be seen that the filters converged to the correct state
(red, dash-dotted) and that the EKF had a higher convergence rate. Also note that
θ∆TDC estimates had faster initial transients. This indicates a relatively high filter
sensitivity to θ∆TDC errors.
87
Chapter 4. Heat-Release Analysis
EKF PF
0.2 0.2
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
·10−3 ·10−3
8 8
C 2 RMSE [m/(s K)]
4 4
2 2
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
200 200
T w RMSE [K]
T w RMSE [K]
150 150
100 100
50 50
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
Figure 4.6 Filter RMSE with the true model uncertainty Σ1 (black) and an incre-
ased assumed model uncertainty 16Σ1 (red). The latter increased the filter con-
vergence rate and steady state RMSE.
ased exploration of the parameter space. This did, however, increase steady-state
variation. Filter convergence rates decreased when the assumed Q c noise levels
were increased. This was due to the reduced parameter correction for Q c devia-
tion from (2.9), since variation in Q c was more probable. Assumed noise varian-
ces essentially functioned as tuning parameters that determined the trade-off be-
tween speed of convergence and steady-state variation.
Incorporating the stochastic fuel-energy component d in the filters was
shown to be important. This made constant Q c after combustion more proba-
ble and allowed for constant offsets in Q c after θEOC , which improved robustness
88
4.4 Simulation Results
EKF PF
0.2 0.2
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
·10−3 ·10−3
8 8
C 2 RMSE [m/(s K)]
4 4
2 2
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
200 200
T w RMSE [K]
T w RMSE [K]
150 150
100 100
50 50
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
Figure 4.7 Filter RMSE with the correct Q c noise levels σx (black) and increased
assumed noise levels, with 4σx (red). An increased assumed Q c noise level clearly
decreased the convergence rates.
to variation in the burned fuel energy. Without d , the filters attempted to fit the
mean Q c to Q ctot after θEOC , which resulted in non-physical Q c appearances.
89
Chapter 4. Heat-Release Analysis
EKF PF
6 6
2 2
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
·10−2 ·10−2
1.5 1.5
C 2 RMSE [m/(sK)]
C 2 RMSE [m/(sK)]
1 1
0.5 0.5
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
400 400
T w RMSE [K]
T w RMSE [K]
300 300
200 200
100 100
0 0
0 20 40 60 80 100 0 50 100 150 200
cycle [-] cycle [-]
90
4.5 Experimental Results
4.6 Discussion
For the choice of 250 particles, the PF was slower than the EKF (see Fig. 4.5).
The EKF is the minimum-variance unbiased state estimator if the system is linear
and perturbed with Gaussian noise, which could explain the higher convergence
rate in Fig. 4.5. The PF was however less sensitive to model-parameter errors and
had comparable RMSE in stationarity. Here, N p = 250 was chosen so that the
filters would have comparable computation times, and it is possible that the PF
performance would have been improved with an increased N p . The effects of N p
and the re-sampling criterion on PF performance were not addressed here and is
suggested future work.
The suggested parameter-estimation framework could easily be extended
or modified to cover other model parameters. It should however be kept in
mind that observability or identifiability might be degraded or lost when many
91
Chapter 4. Heat-Release Analysis
Table 4.2 Sensitivity to model-parameter errors. The filters were more sensitive
to errors in p in and r c and not as sensitive to other parameter errors. The PF was
overall less sensitive to model errors compared to the EKF. The units for θ∆TDC ,
C 2 and T w are CAD, m/(sK) and K.
EKF PF
λ error, 2 ± 0.2
θ∆TDC = 0 ± 0.07 (0.09) θ∆TDC = 0 ± 0.05 (0.07)
C 2 = 0.0032 ± 0.00045 (0.0004) C 2 = 0.0032 ± 0.0004 (0.0003)
T w = 497 ± 12 (10) T w = 497 ± 7.5 (7.5)
Tc error, 333 ± 20 K
92
4.6 Discussion
Experimental Convergence
6,000
Q initial
EKFfinal
5,000
PFfinal
Q ctot
4,000
3,000
Q c [J]
2,000
1,000
−1,000
−100 −80 −60 −40 −20 0 20 40 60 80 100 120
θ [CAD]
93
Chapter 4. Heat-Release Analysis
Experimental Convergence
4
PF
EKF
θ TDC [CAD]
2
−2
0 50 100 150 200 250 300 350 400
cycle [-]
·10−2
1
0.8
C 2 [m/(sK)]
0.6
0.4
0.2
0
0 50 100 150 200 250 300 350 400
cycle [-]
800
600
T w [K]
400
200
0 50 100 150 200 250 300 350 400
cycle [-]
Figure 4.10 Filter convergence with respect to 400 cycles of experimental data
from the operating point in Table 4.3. The EKF (black) and the PF (red) have com-
parable performance but converged to slightly different state estimates. Both fil-
ters managed to detect a significant top-dead-center offset close to 1 CAD.
94
4.7 Conclusions
4.7 Conclusions
A statistical framework for estimation of unknown heat-release model para-
meters was introduced in this chapter. Within this framework, the EKF and
the PF both seem to be feasible options for on-line estimation. The simulation
results showed that both filters were consistent in converging to the correct pa-
rameter values. The relation between assumed model and heat-release noise
variance determined a trade-off between convergence rate and steady-state
RMSE, and could be used as a tuning parameter. An assumed stochastic
accumulated-heat-release offset showed to be crucial when the injected fuel
energy varied. In reality, such variations are present due to common-rail pres-
sure fluctuations. The model-error-sensitivity results in Table 4.2 indicated that
the filters were more sensitive to intake-pressure and compression-ratio errors,
compared to other parameter errors. The model-error sensitivity was also found
to be dependent on the assumed heat-release noise variance. Furthermore, both
filters showed consistent convergence from different initial states with respect to
experimental data and manually tuned filter parameters, see Figs. 4.9 and 4.10.
95
5
Experimental Setup
Fuel
The fuel used was a mixture of 80 volume % gasoline and 20 volume % n-heptane.
This ratio was chosen based on previous results showing that a fuel octane num-
ber around 80 could be used over a wide range of engine operating points [Ma-
nente, 2010b]. The fuel was mixed together with an Infimeum fuel lubricant to
increase the lifetime of the fuel-injection system which was developed for con-
ventional diesel fuel.
96
5.2 Instrumentation
Figure 5.1 Scania D13 engine, located in the combustion-engine lab at Lund
University. Figure courtesy of Nhut Lahm.
5.2 Instrumentation
Actuation
Fuel Injection The fuel-injection system was a production extra-high
pressure-injection (XPI) common-rail system with solenoid injectors. The
common-rail pressure was regulated with an inlet-metering valve, positioned
prior to a fuel pump used to elevate the pressure in the common rail volume.
Fuel-injection timings and durations were determined by current pulses sent to
the injectors. Current-pulse timings and durations were set from the LabVIEW
control system and actuated with Drivven direct-injection drivers. A more detai-
led description of the fuel-injection system is given in [Källkvist, 2011].
Gas Flow The engine was equipped with two cooled EGR loops, located before
and after the turbine. EGR flows were regulated with two valves. Two valves posi-
tioned prior to the intake manifold were used to regulate the intake temperature
by adjusting the flow over an intercooler. A back-pressure valve positioned after
the tubine was used to create back pressure for sufficient EGR flow. Servo motors
used for valve actuation were controlled from the LabVIEW control system and
actuated with Drivven drivers. Valve locations are marked in Fig. 5.2.
Engine Speed The engine speed was controlled with an ABB M2BA electrical
motor with a rated power of 355 kW. The motor reference speed was adjusted
manually from the engine control room.
97
Chapter 5. Experimental Setup
air
intercooler
θcool ṁ a
T
p, T
p
θhot
EGR p, T
EGR
p p, T p, T
p, T
T p, T
p
p, T
T compressor θLP
p
p, T θHP
T
p
p turbine
p, T
T
exhaust manifold
intake manifold
λ
p
p, T
T p, T
p
p, T
T
ṁ f
CO2 p
θBP
M,N
fuel em
exhaust
- valve
electrical motor
Sensing
Sensor locations are marked in Fig. 5.2.
Sampling Crank-angle based sampling was enabled through a Leine & Linde
encoder emitting 5 pulses every CAD which triggered sampling of cylinder pres-
sure, engine torque and injector current. Other sensor signals were sampled
every engine cycle.
98
5.2 Instrumentation
• exhaust manifold.
Torque sensor A force sensor integrated in the electrical motor was used to me-
asure engine torque.
Engine Speed Engine speed was obtained from the internal speed measure-
ment of the electrical motor.
Fuel flow A Bronkhorst mini CORI-FLOW M15 mass-flow meter mounted prior
to the fuel system was used to measure the fuel-mass flow.
Air flow A Bronkhorst hot-film air-mass flow meter mounted prior to the com-
pressor was used to measure the air-mass flow.
99
Chapter 5. Experimental Setup
λ sensor A broadband λ sensor mounted after the turbine measured the ex-
haust oxygen concentration.
Emissions Intake and exhaust CO2 and exhaust NOx , HC, CO and O2 levels
were measured with an AVL AMA i60 exhaust-measurement system. Soot levels
were measured with an AVL micro soot sensor measurement unit.
Software
The engine control system was programmed in LabVIEW which is a graphical
programming environment developed by National Instruments. The software
was originally developed by Borgquist for his thesis work [Borgquist, 2013].
Real-time heat-release analysis and controller computations were executed
by the real-time PXI system. Computations were done using floating point arith-
metic, and most of them were done in LabVIEW MathScript RT Module nodes
inside timed loops, triggered every engine cycle. PI controllers were implemen-
ted using the LabVIEW PID advanced VI and QPs were solved using the LabVIEW
quadratic programming VI. The quadratic programming VI had functionalities
useful for model predictive control implementation such as initialization, warm
start of active constraints, various stopping criteria and error flags when feasible
solutions were not found. The user interface was also programmed in LabVIEW.
Signal Processing
The in-cylinder pressure was measured with piezo-electric transducers. This
measurement technique has high cut-off frequency, good linearity and handling
of the environment inside the combustion chamber. The signal given by these
sensors is on the form
p meas = kp + ∆p (5.1)
where p meas is the sensor signal, p the actual pressure, k the sensor conversion
factor and ∆p the sensor offset. Methods for determining k and ∆p were pre-
sented by Randolph [1990]. In this work, k were known from sensor calibration
100
5.3 Control-System Architecture
and ∆p was determined by referencing the cylinder pressure at intake valve clo-
sing (IVC) to the measured intake-manifold pressure p in . High-frequency con-
tent in p meas during combustion and expansion was attenuated using a digital
zero-phase filter.
101
6
Proportional-Integral
Combustion-Timing Control
102
Chapter 6. Proportional-Integral Combustion-Timing Control
Disturbance
r
θ50 P e θSOI θ50
Controller Engine
−1
Figure 6.1 The combustion-timing feedback loop. The indicator for combustion
timing θ50 is obtained from the measured in-cylinder pressure and heat-release
r −θ caused by,
analysis. The controller varies θSOI to counteract the error e = θ50 50
e.g., changes in fuel-amount or intake conditions.
where the injection timing of the following cycle θSOI (k + 1) is determined by the
previous-cycle error
r
e(k) = θ50 (k) − θ50 (k) (6.2)
multiplied with the proportional gain k p , and the integral term I (k), which is the
sum of previous errors, scaled with the integral gain k I . The integral term is intro-
duced to bring e to zero in steady state.
When introducing feedback, the controller has to be robustly designed to en-
sure closed-loop stability. It is also important that the controller does not en-
hance stochastic cycle-to-cycle variation. In PI-controller design, this is done by
carefully deciding the controller gains k p and k I . The problem of deciding k p and
k I is investigated in this chapter by evaluating controller performance through
simulation. Simulation allows for evaluation of a large number of gain combina-
tions at different engine loads.
The effect of fuel reactivity on controller design is addressed by evaluating
controller performance for different primary reference fuels (PRF). A PRF is
a mixture of n-heptane and iso-octane, where the PRF number indicates the
iso-octane volume percentage. Primary reference fuels are commonly used as
reference in engine research and are used to determine the octane number of
a fuel. The octane number (ON) is a measure of the fuel resistance to autoigni-
tion which increases with ON, and is defined by the PRF value needed to provide
equivalent autoignition properties.
The presented controller-gain evaluation provides gains that simultaneou-
sly maximize attenuation of θ50 disturbances and fulfill constraints on robust-
103
Chapter 6. Proportional-Integral Combustion-Timing Control
ness and noise sensitivity. The evaluation also investigates the trade-off between
these performance requirements. The optimization-based approach to PI con-
troller design was inspired by works presented in [Hast et al., 2013; Garpinger
and Hägglund, 2015], where understanding and rules of thumb for PI controller
design were found through optimization.
The chapter is outlined as follows: The model used for controller evaluation
is introduced in Sec. 6.1. The steady-state relation between θSOI and θ50 is stu-
died in Sec. 6.2. Section 6.3 presents the criteria used for controller evaluation.
Robust and noise-insensitive controllers, found through simulation of transient
and steady-state operation for different PRFs are presented in Sec. 6.4, together
with an analysis of the results. Finally, conclusions are given in Sec. 6.5.
6.1 Modeling
The in-cylinder state and wall-surface temperature were modeled using the
zero-dimensional model presented in Sec. 2.2. Constant-volume combustion
and static gas-exchange boundary conditions were assumed to speed up com-
putations. A detailed description of this model is given in [Widd et al., 2008]
where it was shown to successfully predict experimental θ50 and p IMEPg . The
ignition-delay τ was computed using the model M4, presented in Sec. 2.4
ζ(T,PRF) Λ(T,PRF)
τ = φα(T,PRF) p β(T,PRF) xO e (6.3)
2
The model of the closed engine cycle was used in closed-loop simulation ex-
periments, for which the model outputs θ50 and p IMEPg were regulated on a
cycle-to-cycle basis using θSOI and the injected fuel energyf Q ctot . Model para-
meters used are presented in Table 6.1.
104
6.2 Steady-State Characteristics
Table 6.1 Model parameters used in simulation. Three different convection co-
efficients were used during compression (compr.), expansion (exp.) and gas ex-
change (gas ex.).
105
Chapter 6. Proportional-Integral Combustion-Timing Control
20
Q ctot = 2000 J
θ 50 [CAD]
10 ← Not Controllable
Q ctot = 6000 J
0
−10
−40 −30 −20 −10 0 10 20
θ SOI [CAD]
Figure 6.2 Steady-state θ50 / θSOI relation for PRF20 and PRF80 with fuel ener-
gies from 2000 to 6000 J. The Q ctot sweep generates a band of θ50 / θSOI relations.
For late θSOI , θ50 is excessively delayed which means that ignition never occurs.
Controllability is reduced when θSOI decreases.
Table 6.2 Relation between PRF and Tin in order to obtain θ50 = 5 CAD with
λ = 2 and Q ctot = 2000 J.
106
6.2 Steady-State Characteristics
20
Q ctot = 2000 J
15
θ 50 [CAD]
10
Q ctot = 6000 J
0
-5
-10
−40 −30 −20 −10 0 10 20
θ SOI [CAD]
Figure 6.3 Steady-state θ50 / θSOI relation with adjusted Tin are presented as
layered bands. The lower red band corresponds to PRF0 and the upper blue band
to PRF100.
107
Chapter 6. Proportional-Integral Combustion-Timing Control
6
τ [ms]
2 PRF0
PRF100
NTC behavior
0
−30 −20 −10 0 10 20 30
θ [CAD]
Figure 6.4 Ignition delay τ used in the autoignition criterion in (2.69) as a func-
tion of θ with adjusted Tin according to Table 6.2. The low-value PRFs exhibit a
NTC behavior during compression and expansion. Note that both temperature
and PRF value are varied for the different curves, and that the NTC behavior is
dependent on both parameters. It was decided to study NTC behavior after Tin
adjustment since this would be necessary prior to engine operation to get satis-
factory operation for the high PRF values.
implies that a PI controller alone is not sufficient. A controller also has to be able
r
to detect misfire, and if θ50 is infeasible to avoid wind up. The controller should
r
then take action by adjusting θ50 or limit θSOI . Design of such a controller is not
covered here.
108
6.3 Controller Evaluation
PRF0
PRF100
5,000
Q ctot [J]
4,000
3,000
Figure 6.5 The region where 0.1 < ∂θ50 /∂θSOI < ∞. In order for the controller
to perform satisfactorily, θSOI should be limited within this region. Note that the
region is narrower for low load and PRF100.
Load Control
The gross indicated mean effective pressure p IMEPg was also controlled during
simulation. This was done using a PI controller that adjusted Q ctot for tracking of
r
a set point p IMEPg . The PI-controller gains were tuned for a response time of 10
cycles. Robustness and measurement-noise sensitivity were therefore evaluated
with respect to the resulting multiple input/output system to account for cross
coupling from Q ctot to θ50 , and from θSOI to p IMEPg .
Disturbance Rejection
r
The θ50 -controller objective is to follow θ50 changes whilst attenuating the effects
of disturbances, see Fig. 6.1. The controller ability to fulfill this objective was me-
asured by computing the discrete-time integrated absolute error (IAE)
N
X
IAE = |e(k)| (6.5)
k=1
during the simulation experiments, where N is the number of cycles. The IAE is a
commonly used measure for controller evaluation [Åström and Hägglund, 2006].
109
Chapter 6. Proportional-Integral Combustion-Timing Control
Here, it will mainly penalize transient control error since the controllers exhibit
error-free tracking in steady state due to integral action.
Robustness
It is important that the controller is stable around the set point. This can be en-
sured by having sufficient stability margins. For linear discrete-time systems, ro-
bustness is captured by the sensitivity function S and the complementary sensi-
tivity function T [Zhou and Doyle, 1998]
Here, P 11 is the transfer function from θSOI to θ50 , P 22 is the transfer function
from Q ctot to p IMEPg and P 12 , P 21 represent the corresponding cross-coupling
dynamics. Furthermore, k idj is a direct gain, k ires
j
models the one-cycle delayed
residual-gas effect, and a wall and k iwall
j
determine the time constant and gain of
the wall-temperature dynamics. Step responses for P are presented in Fig. 6.6, for
PRF0 at p IMEPg = 6 and θ50 = 12.
110
6.3 Controller Evaluation
P 11 P 12
−1
residual gas
1.4 −1.2
θ 50 [CAD]
θ 50 [CAD]
−1.4
1.2 −1.6
wall temperature
−1.8
direct effect
1 −2
0 100 200 300 400 0 100 200 300 400
cycle [-] cycle [-]
P 21 P 22
0 3
p IMEPg [bar]
p IMEPg [bar]
−0.05
2
−0.1
1
−0.15
−0.2 0
0 100 200 300 400 0 200 400 600 800
cycle [-] cycle [-]
Figure 6.6 Step response for (6.8). Input steps of θSOI = 1 CAD and Q ctot = 1000
J are applied at cycle 1. The P 11 subdiagram indicates the direct effect, the
residual-gas effect and the wall-temperature dynamics, respectively.
k Ii
C i (z) = z −1 (k pi + ) (6.10)
z −1
Where z −1 is the delay from measurement to actuation. The reason for desig-
ning a decoupled controller, and for designing C 2 given C 1 , was because of weak
coupling (see P 12 and P 21 in Fig. 6.6) and the convenience of designing a load
controller with respect to other aspects than θ50 .
Noise Sensitivity
It is necessary to avoid excessive θSOI variation due to stochastic cycle-to-cycle
variation. This requirement was formulated as a constraint on the steady-state
111
Chapter 6. Proportional-Integral Combustion-Timing Control
σθSOI ≤ κσ (6.11)
Optimization Problem
In summary, the design problem is described by the following optimization pro-
blem from PRF0 to PRF100
N
X
minimize |e(k)| (6.12)
kp , kI k=1
subject to M s ≤ 1.4
M t ≤ 1.4
σθSOI ≤ 0.25
The optimization problem was solved by simulating the system model and eval-
uating (6.12) over a grid of controller gains. This was not only done to find opti-
mal gains but also to investigate cost function and constraint characteristics as a
function of k I , k p , Q ctot and PRF.
6.4 Results
The optimization problem (6.12) was solved by evaluating a grid of gains k p , k I ∈
{0.05, 0.1, . . . , 1} and PRFs ∈ {0, 10, . . . , 100} during two simulation experiments:
The relative air-fuel ratio λ was set equal to 2 by adjusting p in , whereas Tin was
adjusted as a function of PRF according to Table 6.2.
112
6.4 Results
0
θSOI
−10
PRF0
PRF100
−20
100 200 300 400 500 600 700
cycle [-]
Figure 6.7 Optimal θ50 and θSOI trajectories from PRF0 to PRF100. Changes in
r
Tin were applied at cycles 50, 100, 200, 250, 400, 450, 550 and 600, p IMEPg chan-
ges were applied at cycles 350 (increase) and 700 (decrease), and θ50r was varied at
r
cycles 150, 300, 350, 500, 650 and 750. During the p IMEPg r was also chan-
steps, θ50
ged in order to increase θSOI variation. The symbols indicate the points of linea-
rization where M s and M t where computed, see Table 6.3 for static-gain values.
When comparing PRF performance, the difference in θ50 was not as significant as
the difference in θSOI for the different PRFs. It can be seen that changes in p IMEPg
r resulted in greater θ
and θ50 SOI variation for lower PRF values, and that the θSOI
response was more comparable during Tin changes.
Optimal Gains
Optimal transient θ50 and θSOI trajectories with respect to (6.12) are presented in
Fig. 6.7, where optimal time constants of the closed loops are within 10 cycles for
all disturbances and set-point changes. When comparing different PRFs, it can
be seen that the difference in θ50 was not as significant as the difference in θSOI .
r
Load and θ50 variation gave greater θSOI variation for low-value PRFs, whilst θSOI
variation was comparable for different PRFs during Tin changes.
Linearization points for which robustness was evaluated are indicated by the
symbols 5, 2, ° and ∗ in Fig. 6.7. Computed static gains with respect to these
points are presented in Table 6.3 for PRF0 and PRF100. It can be seen that the P 11
gains were higher for late θ50 and PRF100, and that the interaction from Q ctot to
r
θ50 (P 12 ) was higher at low p IMEPg .
Optimal steady-state performance is presented in Fig. 6.8. The θSOI standard
deviation σθSOI was higher for the low-load operating points (cycles 1-500), σθSOI
was also higher for lower PRF values.
113
Chapter 6. Proportional-Integral Combustion-Timing Control
Table 6.3 Static gains at the linearization points. Units are given by [-], [CAD/kJ],
[bar/CAD] and [bar/J] for P 11 , P 12 , P 21 and P 22 . Note the following trends: higher
P 11 gains at late θ50 and PRF100 and a larger interaction from Q ctot to θ50 at low
p IMEPg .
PRF0
P5 P2
r r r r
p IMEPg = 5 bar, θ50 = 6 CAD p IMEPg = 5 bar, θ50 = 12 CAD
à ! à !
0.35 −0.005 1.35 −0.007
0.0033 0.0029 −0.07 0.003
P∗ P°
r r r r
p IMEPg = 15 bar, θ50 = 6 CAD p IMEPg = 15 bar, θ50 = 12 CAD
à ! à !
0.9 −0.001 1.35 −0.0015
−0.056 0.0029 −0.16 0.0029
PRF100
P5 P2
r r r r
p IMEPg = 5 bar, θ50 = 6 CAD p IMEPg = 5 bar, θ50 = 12 CAD
à ! à !
0.7 −0.0025 1.85 −0.004
−0.0011 0.0028 −0.1 0.0029
P∗ P°
r r r r
p IMEPg = 15 bar, θ50 = 6 CAD p IMEPg = 15 bar, θ50 = 12 CAD
à ! à !
1 −0.0005 1.4 −0.001
−0.065 0.0029 −0.16 0.0028
Finally, optimal controller gains are presented in Fig. 6.9 as a function of PRF.
Both k p and k I decreased with PRF and k p was slightly lower than k I . The M s
value at the 2-linearizing point in Fig. 6.7 was always found to be the constraint
limiting controller-gain magnitudes. At this point, the system gain ∂θ50 /∂θSOI
was highest among the linearization points, see Table 6.3. This partial derivative
also increased with PRF, which explains the trend in Fig. 6.9. These trends could
have been anticipated by studying Fig. 6.3, where ∂θ50 /∂θSOI increased with θ50 ,
PRF value, and decreased with Q ctot . The remainder of this chapter investigates
how the robustness and noise-sensitivity constraints vary with PRF.
114
6.4 Results
10
σθSOI = 0.052
θ 50 , θ SOI [CAD]
σθSOI = 0.123
σθSOI = 0.044
0
σθSOI = 0.136
σθSOI = 0.111
σθSOI = 0.180
−20
100 200 300 400 500 600 700 800 900 1,000
cycle [-]
Figure 6.8 Optimal steady-state performance from PRF0 to PRF100. The θSOI
standard deviation σθSOI was higher for the low p IMEPg points (cycles 1-500), and
for the lower PRFs. The presented data is a part of the steady-state experiment
consisting of 8000 cycles.
0.5
kp
kI
Controller Gains [-]
0.4
0.3
0.2
0 10 20 30 40 50 60 70 80 90 100
PRF [-]
Figure 6.9 Optimal controller gains as a function of PRF. The M s < 1.4 constraint
at the 2 operating point was consistently limiting the controller gains. The sys-
tem gain ∂θ50 /∂θSOI was highest at this point. Furthermore, ∂θ50 /∂θSOI increased
with PRF which explains the trend in this figure.
115
Chapter 6. Proportional-Integral Combustion-Timing Control
1.4
20
100
1. 4
1. 4
1.4
150
0.6
k I [-]
1. 4
150
1. 4
0.4 1.4
1.4
1.4
200
250
1.4
300
1.4
0.2
400
Figure 6.10 IAE (black), M s = 1.4 (solid) and M t = 1.4 (dashed) level curves as
a function of controller gains k p and k I , for PRF0 and PRF100. The IAE score de-
creased as k p and k I increased. The allowed gains for which M s , M t < 1.4 are
encircled by the colored lines, and optimal gains are indicated by circles.
Robustness
To illustrate robustness-constraint characteristics, IAE level curves (black, solid)
are presented in Fig. 6.10, together with level curves of M s = 1.4 (solid) and M t =
1.4 (dashed) as a function of controller gains for PRF0 and PRF100. The IAE score
decreased as k p and k I were simultaneously increased. The stability limits are
indicated by the steep increase in IAE in the upper-left and lower-right regions
of the figure. The constraint on M t was less restrictive than the constraint on M s ,
and the M s and M t constraints were overall less restrictive for PRF0. The limiting
r r
M s value was computed at the late θ50 and low p IMEPg operating point 2. The
PRF0 fuel also had visible M s and M t constraints for low controller gains due to
interaction with the p IMEPg loop. This was a result of the higher P 12 gain for PRF0.
These results can be explained by linear-systems analysis. With the
116
6.4 Results
Nyquist Diagram
0.5
0
Imaginary Axis [-]
−0.5
−1
Figure 6.11 Nyquist curves with optimal controller gains for PRF0 (red) and
PRF100 (blue), for (6.13) (solid), P 1C 1 (dashed), and P 1C 1 without residual and
wall-temperature dynamics (dotted). It can be seen that the simple model has the
smallest stability margin.
¡ P 12C 2 P 21 ¢
P 11 − C1 (6.13)
1 + P 22C 2
Nyquist curves with k I = 0.35, k p = 0.3 for PRF0 (red) and PRF100 (blue) are pre-
sented in Fig. 6.11. This figure presents Nyquist curves for (6.13) (solid), P 1C 1
(dashed) and P 1C 1 without residual and wall-temperature dynamics (dotted). It
can be seen that the simplest open-loop transfer function has the smallest sta-
bility margin. The intuition behind this result is that both residual-gas dynamics
and the p IMEPg controller counteract the θSOI effect on θ50 , which decreases the
open-loop gain. Analysis of the simple system can therefore be used to compute
stability margins that are sufficient for the more detailed loops.
When omitting wall-temperature, residual-gas and p IMEPg -loop dynamics,
r
the linearized closed-loop pulse-transfer function, from θ50 to θ50 is given by
K C 1 (z)
Hcl (z) = (6.14)
1 + K C 1 (z)
117
Chapter 6. Proportional-Integral Combustion-Timing Control
Stability Region
1.4
← K = 0.9
1.2
←K =1
1
← K = 1.1
0.8
k I [-]
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
k p [-]
Figure 6.12 PI-gain stability region, computed using the critera in (6.16). The
stability region shrinks as K increases. The stability limits are similar to the IAE
level curves in Fig. (6.10).
where K = ∂θ50 /∂θSOI . Inserting the pulse transfer function for (6.10) into (6.14)
yields the characteristic polynomial for the closed-loop poles
z 2 + z(K k p − 1) + K (k I − k p ) (6.15)
By applying the Jury stability criterion [Jury, 1964], and assuming positive con-
troller gains, the conditions for stability is obtained
k p < k I /2 + 1/K
(6.16)
k p > k I − 1/K
Noise Sensitivity
Figure 6.13 presents IAE (black) and σSOI = 0.25 level curves for PRF0 and PRF100
as functions of k p and k I . Allowed gains for which σSOI ≤ 0.25 are within the red
and blue lines. The σSOI constraint was more restrictive for PRF0, and the most
restrictive σSOI constraint was found at late θ50 and low p IMEPg for both fuels, see
Fig. 6.8.
118
6.4 Results
0
10
00
5
0.2
4
0
0
30
15
0.8 0
25 10
0
2 00 100
0.25
5
0.2
0.25
0.25
0.6
k I [-]
0
0.2
5 15
0.25
150
0. 25
0.4
200
250
0.25
300
5
0.2 0.2
400
Figure 6.13 IAE (solid, black) and σSOI = 0.25 level curves for PRF0 and PRF100
as a function of k p and k I . The allowed gains for which σSOI ≤ 0.25 are within
the solid colored lines. The most restrictive σSOI = 0.25 constraints, for both fuels
were found at the late θ50 at the low-load operating points, see Fig. 6.8.
The θSOI standard deviation can be evaluated by studying the H2 -norm of the
transfer function S c
119
Chapter 6. Proportional-Integral Combustion-Timing Control
0 0
∂ θ 50 / ∂Q ctot [CAD/J]
∂ θ 50 / ∂Tin [CAD/K] −0.2
−5
−0.4
−0.6 −10
−0.8
5 10 15 5 10 15
θ 50 [CAD] θ 50 [CAD]
0
∂ θ 50 / ∂p in [CAD/bar]
0.6
−2
||S c (z)||2
0.5
−4
0.4
−6
0.3
5 10 15 5 10 15
θ 50 [CAD] θ 50 [CAD]
Figure 6.14 The closed-loop noise sensitivity depends on θ50 partial derivatives
with respect to Tin , p in , Q ctot and the S c norm. This figure presents these quantities
as a function of θ50 for PRF0, PRF100 with p IMEPg = 5 bar (LL) and 15 bar (HL).
Partial derivatives of θ50 and ||S c ||2 with constant controller gains are presen-
ted in Fig. 6.14 for p IMEPg = 5 and 15 bar. The partial-derivate magnitudes were
clearly larger at low p IMEPg and late θ50 . It can also be seen that PRF0 was more
sensitive to Q ctot and p in whilst PRF100 was more sensitive to Tin . The closed-loop
noise sensitivity ||S c ||2 was higher for PRF100 and decreased with θ50 . Similar ex-
perimental PRF trends were presented in [Sjöberg and Dec, 2005], where Tin was
adjusted according to PRF.
The standard deviation σθSOI , computed using (6.18) is presented in Fig. 6.15.
It can be seen that σSOI decreased with p IMEPg and θ50 . Overall, σSOI was also
higher for PRF0 which agrees with the observed trends in Figs. 6.8 and 6.13.
120
6.5 Conclusions
Noise Sensitivity
0.3
PRF100, HL
PRF100, LL
0.25 PRF0, HL
PRF0, LL
σθSOI [CAD]
0.2
0.15
0.1
0.05
4 6 8 10 12 14
θ 50 [CAD]
Rules of Thumb
To summarize, the simulation results provided the following rules of thumb for
PI-controller tuning:
1. Robustness was limited by operating points with large ∂θ50 /∂θSOI . These
were found at late θ50 and low p IMEPg . At these points, ∂θ50 /∂θSOI incre-
ased with PRF value, meaning that robustness constraints became more
restrictive for higher PRF values.
2. Noise sensitivity was higher at low loads and late combustion timings
where gains from disturbances to θ50 were higher. The noise sensitivity was
also higher for lower PRF values due to an increased Q ctot and p in sensitivity.
3. A robust and noise-insensitive controller-gain choice for all fuels and ope-
rating points is given by 0.2 < k p < k I < 0.35.
6.5 Conclusions
This chapter showed that the θSOI interval for which θ50 is controllable is limited
between a low-gain limit for early θSOI and a misfire limit for late θSOI . In order
to obtain early combustion timings for higher PRFs, Tin had to be adjusted as a
121
Chapter 6. Proportional-Integral Combustion-Timing Control
function of PRF. Even after Tin adjustment, the PRFs have different θ50 characte-
ristics, partly due to the varying NTC behavior for the different fuels. The com-
putation of suitable θSOI limits is a part of the controller design where such limits
could be a function of fuel, load and intake conditions.
Optimal PI-controller gains found through simulation were limited by the
high system gain at late θ50 and low p IMEPg . Controller measurement-noise sen-
sitivity was also found to be higher at this operating point due to an increased θ50
sensitivity to load and intake-condition variation. Controller gain requirements
also varied with the PRF value: robustness was lower for higher PRFs due to an
increased system gain, and noise sensitivity was higher for lower PRFs due to a
higher Q ctot and p in sensitivity.
122
7
Model-Based Control of
Combustion Timing and
Ignition Delay
7.1 Introduction
A sufficiently long ignition delay τ is a prerequisite for the fuel / air mixing le-
ading to premixed low-temperature combustion, as described in Sec. 1.3. Mea-
nwhile, it is also important that the combustion timing θ50 (4.3) is sufficiently
well timed for high thermodynamic efficiency. A too early θ50 leads to increased
heat-transfer and inefficient pressure build-up during the compression stroke,
and a late θ50 results in high exhaust losses and a lowered combustion efficiency.
Simultaneous control of τ and θ50 is an important component for a successful
implementation of partially premixed combustion, and when controlling these
two variables, one has to pay attention to their coupling through θSOI and ther-
modynamic in-cylinder conditions.
This chapter studies model-based multi-cylinder control of τ and θ50 with
combined actuation of the gas-exchange and fuel-injection systems. The objec-
tive is to regulate τ and θ50 during load disturbances and set-point changes. This
is an under-determined control problem due to more output than input variables
for the engine setup used.
A model predictive control (MPC) is suggested, see Sec. 3.1. This is a suita-
ble design choice for multiple input/output systems with actuator constraints.
The controller obtains τ and θ50 from in-cylinder pressure measurement and
heat-release analysis, and the ignition-delay model M1 (see Sec. 2.4) is linearized
every engine cycle for model-based prediction. The MPC feedback loop is illus-
trated in Fig. 7.1, where the multiple cylinders are indicated by the (bold type)
vector notation. The MPC computes fuel-injection timings and valve positions
of the two EGR paths and the fast thermal-management (FTM) system. Valve po-
123
Chapter 7. Model-Based Control of Combustion Timing and Ignition Delay
Disturbance
τr , θ50
r
θ SOI τ, θ 50
MPC Engine
θ EGR
θ FTM
Linearization
−1
Figure 7.1 This chapter presents and studies the following MPC feedback loop.
The controller obtains τ and θ 50 from in-cylinder pressure measurement and
heat-release analysis, and then linearizes the ignition-delay model M1, presented
in Sec. 2.4. The controller computes fuel-injection timings θ SOI and valve posi-
tions of the two EGR paths θ EGR and the fast thermal-management system θ FTM .
Multiple cylinders are indicated by the (bold type) vector notation.
7.2 Modeling
Model M1 (see, Sec. 2.4) was used to model τ, whereas a calibrated stat-
ic model was used to model the gains from gas-exchange valve positions to
intake-manifold composition and temperature. This section describes these
models and how they were used in the controller design.
124
7.2 Modeling
125
Chapter 7. Model-Based Control of Combustion Timing and Ignition Delay
310
300
290
10 20 30 40 50 60 70 80 90
θ hot [deg]
θ LP , p IMEP = 5 bar
10 θ LP , p IMEP = 10 bar
0 10 20 30 40 50 60 70 80 90
θ HP,LP [deg]
Figure 7.2 TIVC and [O2 ]IVC as functions of θLP , θHP and θhot at two different
loads and Nspeed = 1200 rpm. These experimentally obtained functions were used
to model the relation between the gas-exchange valve positions and the intake
manifold conditions. Note that the gain from θHP decreased with load. This was
because of the decreased pressure difference over the exhaust and intake mani-
folds due to increased turbocharger boost.
126
7.2 Modeling
∂τ τ([O2 ]IVC + ∆[O2 ]IVC /2) − τ([O2 ]IVC − ∆[O2 ]IVC /2)
≈
∂[O2 ]IVC ∆[O2 ]IVC
µ ¶ ∆θSOIi (k)
∂τi ∂τi ∂τi
τi (k + 1) = τi (k) + i ∆TIVC (k)
∂θSOI ∂TIVC ∂[O2 ]IVC
∆[O2 ]IVC (k)
θ50i (k + 1) = θ50i (k) + ∆θSOIi (k) (7.7)
µ ¶ ∆θSOIi (k)
∂θ ∂τi ∂τi ∂τi
+ i ∆TIVC (k)
∂t ∂θSOI ∂TIVC ∂[O2 ]IVC
∆[O2 ]IVC (k)
∆θSOIi (k)
µ ¶
∂τi ∂τi ∂τi ∂τi ∆θhot (k)
τi (k + 1) = τi (k) +
i
∂θSOI ∂θhot ∂θHP ∂θLP ∆θHP (k)
∆θLP (k)
(7.8)
∆θSOIi (k)
µ ¶
∂θ50i ∂θ50i ∂θ50i ∂θ50i ∆θhot (k)
θ50i (k + 1) = θ50i (k) +
∂θSOIi ∂θhot ∂θHP ∂θLP ∆θHP (k)
∆θLP (k)
127
Chapter 7. Model-Based Control of Combustion Timing and Ignition Delay
The complete linear state-space model can be written on more compact form
∆θ SOI (k)
à ! à !
θ 50 (k + 1) θ 50 (k) ∆θhot (k)
= +B
∆θ (k)
(7.9)
τ(k + 1) τ(k) HP
∆θLP (k)
where
à !
θ 50 (k) ¡ ¢T
= θ501 (k) ... θ506 (k) τ1 (k) ... τ6 (k)
τ(k)
¡ ¢T
∆θ SOI (k) = ∆θSOI,1 (k) ... ∆θSOI,6 (k)
∂θ ∂θ50,1 ∂θ50,1 ∂θ50,1
50,1
··· 0
∂θSOI,1 ∂θhot ∂θHP ∂θLP
.. .. .. .. .. ..
. . . . . .
(7.10)
∂θ50,6 ∂θ50,6 ∂θ50,6 ∂θ50,6
0 ···
∂θSOI,6 ∂θhot ∂θHP ∂θLP
B =
∂τ1 ∂τ1 ∂τ1 ∂τ1
... 0
∂θ ∂θhot ∂θHP ∂θLP
SOI,1
.. .. .. .. .. ..
. . . . . .
∂τ6 ∂τ6 ∂τ6 ∂τ6
0 ···
∂θSOI,6 ∂θhot ∂θHP ∂θLP
The system is not controllable due to more outputs than inputs, which is indica-
ted by the fact that the controllability matrix C
¡ ¢
C= B ... B (7.11)
has less than 12 independent columns. The system outputs τ and θ 50 can there-
fore not be controlled independently and we are left with two options. Combus-
tion timings θ 50 can either be controlled cylinder individually by varying θ SOI
and then letting gas-exchange actuators control τ. The other option is to control
τ cylinder individually and let the gas-exchange actuators control θ 50 . The first of
the two approaches was taken here. The reason for this choice was the nonlinear
relationship between θ SOI and τ (see Fig. 6.4), and the monotone relationship
between τ and TIVC , [O2 ]IVC . It was also shown in Sec. 2.4 that it is difficult to ac-
curately model the sign of ∂τ/∂θSOI when θSOI is close to TDC. Such model errors
could lead to unwanted closed-loop behavior if θ SOI is set to control τ.
The desired controller can be obtained by MPC-weight tuning, which is cove-
red in the following sections.
128
7.3 Model Predictive Control Formulation
Hp ³
X r
minimize ω1 ||θ50 (k) − θ 50 (k)||22 + ω2 ||τr (k) − τ(k)||22
∆θ SOI , ∆θHP , k=1 ´
∆θLP , ∆θhot + ω3 θHP (k)2 + ω4 θLP (k)2
X ³
+ ω5 ||∆θ SOI (k)||22 + ω6 ∆θhot (k)2
k∈k Hc ´
+ ω7 ∆θHP (k)2 + ω8 ∆θLP (k)2 (7.12)
θ SOI
θhot
subject to θ l ≤
θHP ≤ θ
u
θLP
Here, k is the cycle index, and || · ||2 is the Euclidian norm in R6 . Initial condi-
tions at the current cycle k = 0 are obtained from measurements. The first sum
penalizes θ50 and τ set-point error, and the usage of EGR over the prediction ho-
r
rizon H p . Set points θ50 and τr are considered to be precomputed as a function
of the engine operating point. It was decided to penalize EGR-valve opening area
to favor flow over the turbine and to not use more EGR than needed. The terms
in the second sum penalize control action over the control horizon Hc . The cost
function should be minimized subject to absolute constraints on actuators and
the linearized system dynamics in (7.9).
129
Chapter 7. Model-Based Control of Combustion Timing and Ignition Delay
ω1 ω2 ω3 ω4
2 4
10 10 5000 3000
ω5 ω6 ω7 ω8
−2
0.1 10 400 800
3 ≤ θHP ≤ 40 20 ≤ θLP ≤ 50 5 ≤ θhot ≤ 85 −40 ≤ θSOI ≤ 10
zero, see Fig. 7.2. Other actuator bounds were chosen to fulfill physical limita-
tions, and θSOI was limited to avoid misfire, see Fig 6.2.
The weights and constraints used are presented in Table 7.1. Solution trajec-
tories with these weights are presented in Fig. 7.3 for an arbitrary initial condi-
tion. The dashed black lines in the upper part of Fig. 7.3 are set points τr and
r
θ50 . The cylinder numbers are indicated by different colors, cylinder 1 being the
upper cylinder in Fig. 2.3.
The weights were tuned to obtain the desired behavior where θ SOI controls
cylinder-individual θ 50 . This was done by prioritizing θ 50 tracking, and let the
slow gas-system actuators control the mean τ. Note that the grid over which θ SOI
is allowed to change is denser initially than for the gas system valves. This is be-
cause the actuation of θ SOI is much faster than for the valve positions.
Solving (7.12) is a quadratic program (QP) and was solved in LabVIEW using
the QP-solver VI. The solver used the previous-cycle solution and active set as
initial guesses for the next cycle to shorten computation times. Early termina-
tion was also used for the solver to finish within one engine cycle. These are well
known methods for speeding up MPC execution, see [Wang and Boyd, 2010]. The
average computation time for differentiating Eqs. (7.6) and constructing the QP
matrices was 10 ms, while it took 25 ms on average to solve (7.12). These com-
putations were repeated every engine cycle after sampling of the previous-cycle
cylinder pressure.
130
7.4 Experimental Results
2.6 10
θ 50 [CAD]
τ [ms]
2.4
6
2.2 4
2
20 40 20 40
cycle [-] cycle [-]
18 −11
16 −12
θ SOI [CAD]
θ HP [deg]
14 −13
12 −14
10 −15
20 40 20 40
cycle [-] cycle [-]
80 θ cool θ hot 32
θ hot /θ cool [deg]
70 30
θ LP [deg]
60 28
50 26
40 24
20 40 20 40
cycle [-] cycle [-]
Figure 7.3 Solution trajectories of (7.12) using the weights and constraints in
Table 7.1 and an arbitrary initial condition. The inputs were allowed to change
at predefined cycle indices over Hc . The weights were chosen to prioritize θ 50 tra-
cking, and let the slower gas-system actuators control the mean τ. Note that the
samples for which θ SOI was allowed to change was denser initially than for the
gas-system valves. This was because control action of θ SOI was much faster than
for the valve positions.
Set-Point Tracking
Set-point tracking performance was evaluated by keeping p IMEPn and Nspeed
r
constant and varying θ50 and τr . System inputs and outputs during 800 cycles
r
of θ50 step changes are presented in Fig. 7.4. The controller weights ωi (see Table
131
Chapter 7. Model-Based Control of Combustion Timing and Ignition Delay
3) were set so that tracking of θ 50 was done by changing θ SOI , τ was then jointly
controlled by the gas-exchange system actuators. This tuning resulted in regula-
tion of the mean cylinder τ.
As θ 50 was delayed at cycle 950 in Fig. 7.4, τ decreased due to the increased
temperature at θ SOI . This forced θhot to close while θHP,LP opened. The control-
ler tried to find the lowest possible θHP in stationarity for higher turbine mass
flows. Some chattering is visible in the gas-exchange actuators. This was due to
stochastic cycle-to-cycle variation in θ 50 and τ which could have been reduced
by introducing filtering or increasing ω6−8 .
A zoom-in around cycle 950 is presented in Fig. 7.5. Here, it can be seen that
θ 50 reached the new set point in 5 cycles. There was an internal delay of 5 cycles
from set point to the controller, caused by the communication from the user in-
terface and the real-time system. The gas system managed to adjust for the τ de-
crease in 50 cycles where cylinder-to-cylinder variation created a τ-distribution
among the cylinders. This variation could be caused by non-uniform EGR distri-
bution to the different cylinders or different cylinder-wall temperatures. It can
be seen that τ6 was consistently shorter whilst τ1 was the longest. A hypothetical
explanation for this is that cylinder 1 is closer to the high-pressure EGR path, see
Fig. 2.3.
In Fig. 7.6, system inputs and outputs are displayed during 1400 cycles for
which step changes in τr were made. The tracking of τr was realized by varying
θhot and θHP,LP whilst θ SOI was varied to keep θ 50 constant. The high-pressure
EGR valve opened too much initially which gave a slight overshoot in τ. A
zoom-in around cycle 550 is presented in Fig. 7.7. Here, it can be seen that θ SOI
was varied to keep θ 50 within 0.6 CAD whilst τ reached the new set point in 50
cycles.
Load Changes
In Fig. 7.8, system inputs and outputs are displayed during 1000 cycles for which
r
p IMEPn steps were made between 6 and 10 bar, and p IMEPn reached its new set
point in 20 cycles. The ignition delays decreased as p IMEPn was increased due to
increased cylinder temperature and richer cylinder mixtures. This forced θHP,LP
to increase. The FTM system was limited to the cold-flow upper limit during the
higher p IMEPn values. Adjustments in θ SOI managed to keep θ 50 within 1 CAD.
In-cylinder data at p IMEPn = 6 (dashed) and p IMEPn = 10 bar (solid) are pre-
sented in Fig. 7.9. The figure shows pressure and heat-release rates for the diffe-
rent cylinders, and the injection current from cylinder 1. The controller manages
to maintain constant θ 50 and τ despite the p IMEPn difference.
Speed Changes
In Fig. 7.10, system inputs and outputs are displayed during 2000 cycles for
which the engine speed Nspeed was varied between 1200 and 1500 rpm. The igni-
132
7.4 Experimental Results
4
10
θ 50 [CAD]
3.5 8
τ [ms]
3 6
4
2.5
800 1,000 1,200 1,400 800 1,000 1,200 1,400
cycle [-] cycle [-]
−10
60
θ SOI [CAD]
θ LP [deg]
−15
40
−20
20
−25
800 1,000 1,200 1,400 800 1,000 1,200 1,400
cycle [-] cycle [-]
100
40
θ cool /θ hot [deg]
80
θ HP [deg]
60
40 20
20
0 0
800 1,000 1,200 1,400 800 1,000 1,200 1,400
cycle [-] cycle [-]
Figure 7.4 System inputs and outputs during step changes in θ50 r .θ
hot is indica-
ted in red and θcool in blue. The black dashed lines in the upper diagrams are θ50 r
r
and τ , whilst the dash-dotted lines are constraints on θLP and θHP .
tion delays decreased as Nspeed was increased. Probably due to increased engine
temperatures and cylinder-gas turbulence levels [Heywood, 1988]. This forced
θHP,LP and θcool to open. The FTM system was limited to the cold-flow limit at
r
Nspeed = 1500 rpm, and θ SOI managed to keep θ 50 within 1 CAD from θ50 .
133
Chapter 7. Model-Based Control of Combustion Timing and Ignition Delay
4
10
θ 50 [CAD]
3.5 8
τ [ms]
6
3
4
2.5 2
920 940 960 980 1,000 920 940 960 980 1,000
cycle [-] cycle [-]
−12
60
−14
θ SOI [CAD]
θ LP [deg]
−16
40
−18
−20
20
−22
920 940 960 980 1,000 920 940 960 980 1,000
cycle [-] cycle [-]
100
40
θ cool /θ hot [deg]
80
θ HP [deg]
60
40 20
20
0 0
920 940 960 980 1,000 920 940 960 980 1,000
cycle [-] cycle [-]
Figure 7.5 A zoom in around cycle 960 in Fig. 7.4 where θ 50 reached a new set
point in 5 cycles. There is a visible internal delay from the set point to the control-
ler of about 5 cycles. The gas system managed to adjust for the τ decrease in 50
cycles.
134
7.5 Discussion and Conclusions
4 8
θ 50 [CAD]
3.5
τ [ms]
6
3
5
2.5 4
500 1,000 1,500 500 1,000 1,500
cycle [-] cycle [-]
−14
60
−16
θ SOI [CAD]
θ LP [deg]
40 −18
−20
20
−22
500 1,000 1,500 500 1,000 1,500
cycle [-] cycle [-]
100
40
θ cool /θ hot [deg]
80
θ HP [deg]
60
40 20
20
0 0
500 1,000 1,500 500 1,000 1,500
cycle [-] cycle [-]
Figure 7.6 System inputs and outputs during 1400 cycles for which steps in τr
were made. The tracking of τr was realized by varying θhot and θHP,LP whilst θ SOI
kept θ 50 constant.
ler relation between intake conditions and τ, compared to the relation between
θSOI and τ. Additional actuation techniques such as variable valve timings, varia-
ble compression ratio or spark ignition could be applied to control both θ50 and
τ cylinder individually
The suggested controller was successful in tracking τ and θ50 set points with
response times of 50 and 5 cycles, respectively. Both in stationarity and du-
ring load and speed changes. A general observation was that MPC-weight tun-
ing was a trade-off between short response times and closed-loop robustness
135
Chapter 7. Model-Based Control of Combustion Timing and Ignition Delay
4 8
θ 50 [CAD]
3.5
τ [ms]
6
3
5
2.5 4
450 500 550 600 650 450 500 550 600 650
cycle [-] cycle [-]
60 −16
θ SOI [CAD]
−18
θ LR [deg]
40 −20
−22
20 −24
450 500 550 600 650 450 500 550 600 650
cycle [-] cycle [-]
100
40
θ cool /θ hot [deg]
80
θ HP [deg]
60
40 20
20
0 0
450 500 550 600 650 450 500 550 600 650
cycle [-] cycle [-]
Figure 7.7 A zoom in around cycle 550 in Fig. 7.6. The ignition delay reached its
new set point after 50 cycles and θ SOI kept θ 50 within 0.6 CAD during the step
change.
136
7.5 Discussion and Conclusions
10
3.5
θ 50 [CAD]
8
τ [ms]
3
6
2.5
4
200 400 600 800 1,000 200 400 600 800 1,000
cycle [-] cycle [-]
−10
60
−12
θ SOI [CAD]
θ LP [deg]
−14
40
−16
−18
20
−20
200 400 600 800 1,000 200 400 600 800 1,000
cycle [-] cycle [-]
100
40
θ cool /θ hot [deg]
80
θ HP [deg]
60
40 20
20
0 0
200 400 600 800 1,000 200 400 600 800 1,000
cycle [-] cycle [-]
12 1.8
p IMEPn [bar]
1.6
θ DOI [deg]
10
1.4
8
1.2
6
1
200 400 600 800 1,000 200 400 600 800 1,000
cycle [-] cycle [-]
Figure 7.8 System inputs and outputs during 1000 cycles for which steps in
r
p IMEPn were made. In addition to the signals displayed in the previous figures,
p IMEPn is plotted together with its set point in the lower left figure. Injection du-
rations θDOI are presented in the lower-right subdiagram. In-cylinder data from
this experiment are presented in Fig. 7.9.
137
Chapter 7. Model-Based Control of Combustion Timing and Ignition Delay
120
100
80
p [bar]
60
40
20
0
−20 −15 −10 −5 0 5 10 15 20 25
θ [CAD]
Figure 7.9 In-cylinder data from the experiment in Fig. 7.8 at p IMEPn = 6 (da-
shed) and p IMEPn = 10 bar (solid). The figure shows pressures and heat-release
rates for the different cylinders and the injection current from cylinder 1. The con-
troller maintains constant θ 50 and τ despite the p IMEPn difference.
improve performance. This would also allow for the gas-exchange efficiency to
be included in the MPC cost function.
The MPC framework provided a simple way of prioritizing system output be-
havior. It also took interaction effects into account. Input constraints and the
cost of using EGR were also incorporated. Comparable controller performance
could probably be obtained by adopting an approach where θ50 is controlled with
θSOI and cylinder-individual controllers, and then let the mean τ be controlled
by the gas-system valve positions, see [Karlsson et al., 2008]. This approach de-
mands less on-line computations but is less general compared to the framework
presented here, where input and output constraints can be incorporated. The
MPC could also be extended to cover p IMEPn control by adding the injected fuel
amount and its effect on τ to the model. Variation in fuel amount could in this
way also be included in the prediction to reduce set-point deviation during load
transients.
A sufficient τ was here considered to be an indicator for low temperature
combustion with favorable emission properties. An assumption that was shown
to be accurate in [Karlsson et al., 2008; Lewander et al., 2008]. In future work, this
controller could be evaluated with emission measurements to conclude if this
hypothesis holds or if supplementary control actions need to be taken.
138
7.5 Discussion and Conclusions
3.6 10
3.4 9
θ 50 [CAD]
τ [ms]
3.2
8
3
2.8 7
2.6 6
0 500 1,000 1,500 0 500 1,000 1,500
cycle cycle
−14
60
θ HP / θ LP [deg]
−16
θ SOI [CAD]
40 −18
−20
20
−22
0 −24
0 500 1,000 1,500 0 500 1,000 1,500
cycle cycle
100 1,600
θ cool /θ hot [deg]
80
Nspeed [rpm]
60 1,400
40
20 1,200
0
0 500 1,000 1,500 0 500 1,000 1,500
cycle cycle
Figure 7.10 System inputs and outputs during 2000 cycles for which Nspeed was
varied, see the lower-right subdiagram. Both EGR-valve positions are now plotted
together in the mid-left figure.
139
8
Pilot Injection
8.1 Introduction
It was discovered in [Manente et al., 2009] that long ignition delays in
single-injection PPC give rise to high pressure-rises rate due to violent HCCI-like
combustion. A long ignition delay creates in-cylinder mixtures where the major-
ity of the injected fuel reaches high-temperature reactions simultaneously, which
result in high heat-release and pressure-rise rates. An example of this can be seen
at the higher-load operating point in Fig. 7.9. A high pressure-rise rate is an in-
dicator for high audible noise levels and could also lead to mechanical engine
damage. The pressure-rise rate therefore has to be kept below certain levels in or-
der to ensure silent and safe engine operation. Previous research by Tsurushima
et al. [2009] implies that pressure oscillations resulting from violent combus-
tion rates are able to break insulating gas boundary layers in the cylinder. High
pressure-rise rates could therefore result in increased heat-transfer flux to the
cylinder walls. The issue of having high pressure-rise rates is not as severe in
conventional diesel engines where the heat-release rate is limited by the rates of
fuel-injection and fuel-air mixing.
A means of counteracting the pressure-rise rate problem is to introduce a pil-
ot injection, e.g., by having a smaller fuel injection earlier during the compres-
sion stroke and then inject the majority of the fuel amount closer to TDC in a
main injection. Optical OH chemiluminescence experiments and computatio-
nal fluid dynamics (CFD) simulations suggest that the pilot injection provides a
background mixture whose reactions increase in-cylinder temperature and as-
sist autoignition of the the main injection [Tanov et al., 2014; Solsjö, 2014]. The
reduced ignition delay of the main injection increases in-cylinder stratification
during combustion with decreased combustion rates as a result.
In-cylinder data showing the effect of introducing a pilot injection is pre-
sented in Fig. 8.1. The two fuel-injection configurations produce the same load,
p IMEP = 5 bar, and the same combustion timing θ50 = 6 CAD, but with different
pressure-rise rates d p max = 11 bar/CAD (blue) and d p max = 30 bar/CAD. Pilot in-
jections are also used in conventional diesel engines, both to improve low-load
140
8.1 Introduction
120
100
80
60
p [bar]
40
20
−20
−25 −20 −15 −10 −5 0 5 10 15 20 25
θ [CAD]
Figure 8.1 In-cylinder data showing the effect of introducing a pilot injection.
The two fuel-injection configurations produce the same load, p IMEP = 5 bar, and
the same combustion timing θ50 = 6 CAD, with different maximum pressure-rise
rates d p max = 11 bar/CAD (blue) and d p max = 30 bar/CAD (red). Pressure-rise
rates are indicated by the dashed tangent lines.
performance [MacMillan et al., 2009; Osuka et al., 1994] and to decrease emis-
sions and engine-noise levels [Kiencke and Nielsen, 2000].
With additional injections, the amount of calibration work for optimized
engine performance at different operating points grows exponentially [Meyer,
2011]. It would therefore be advantageous to have a fuel-injection controller
that automatically adjusts injection timings and the fuel distribution among
the injections depending on the engine operating point. Previous work on
pilot-injection combustion control in low-temperature combustion concepts
has been investigated in [Ott et al., 2013; Eichmeier et al., 2012; Ekholm et al.,
2008; Kokjohn et al., 2009]. Whereas previous works have focused either on ca-
libration and control in open-loop, or on decentralized PI control of θ50 and
d p max , this chapter introduces an MPC that aims to decouple control of θ50 and
d p max , where control of d p max is formulated as an output constraint. The sug-
gested feedback loop is presented in Fig. 8.2, where a Kalman filter is used to
attenuate stochastic output variation.
This chapter begins with a presentation of experimental engine-performance
characteristics in terms of efficiency, emissions and maximum pressure-rise rate
controllability, with respect to pilot-injection parameters in the low-to-mid load
141
Chapter 8. Pilot Injection
Disturbance, e
c r m e e
d p max , θ50 θSOI d p max , θ50
MPC Engine
rp
d p max , θ50
Kalman filter
Figure 8.2 This chapter studies the following feedback loop where the ratio of
pilot-fuel injected r p and main-injection timing θSOIm are adjusted to keep the
c
pressure-rise rate d p max below an upper limit d p max , while the combustion ti-
r . A Kalman filter was used to attenuate sto-
ming θ50 should follow a set point θ50
chastic output variation.
range of the engine. The experimental results are then used to design the con-
troller in Fig. 8.2. An experimental controller evaluation during engine load and
speed changes is also presented in this chapter.
Controlled Outputs
The maximum pressure rise rate d p max is here defined as
dp
d p max = max (8.1)
θ dθ
Due to the high cycle-to-cycle variation in d p max , filtering is necessary for d p max
e
to be used as a feedback variable. The computed d p max was therefore modeled
as the cycle mean d p max , corrupted with additive Gaussian noise e with standard
deviation σd p max
e
d p max = d p max + e d p max , e d p max ∼ N (0, σ2d p max ) (8.2)
Stochastic variation in θ50 was also modeled as additive Gaussian noise accor-
ding to
e
θ50 = θ50 + e θ50 , e θ50 ∼ N (0, σ2θ50 ) (8.3)
e e
A Kalman filter was then used to recover d p max and θ50 from d p max and θ50 . The
Kalman filter will be presented in Sec. 8.3.
Input Variables
The input variables considered are the fuel-injection timings and durations defi-
ned by the current pulses sent to the injector, see Fig. 8.3. The injector-current ri-
m p
sing flanks indicate the start of the main and pilot injections θSOI , θSOI , that occur
142
8.1 Introduction
Input Variables
2
1.5
Injector Current [a.u]
p
1 θ SOI θm
SOI
0.5
0
p
θh θ
DOI
θh θm
DOI
−0.5
Figure 8.3 Injector-current signal with a pilot (p) and a main (m) injection toge-
x , θ x , and θ .
ther with definitions of θSOI DOI h
m p
after an hydraulic injector delay, θh . The fuel-injection durations θDOI and θDOI
are defined as the difference between the injector-current pulse width and θh .
The delay therefore determines the minimum current-pulse duration for which
fuel is injected into the cylinder. It was here assumed to be constant at 0.25 ms.
m p m p
Instead of studying the effects of θSOI , θSOI , θDOI and θDOI explicitly, the pilot
ratio r p
p
θDOI
rp = p m
(8.4)
θDOI + θDOI
the injection separation d SOI
m p
d SOI = θSOI − θSOI (8.5)
m
the main-injection timing θSOI , and the total injection duration
tot m p
θDOI = θDOI + θDOI (8.6)
were considered.
tot m
These variables were chosen because θDOI and θSOI are determined by the
desired load and θ50 at a given operating point. The objective is then to determine
the pilot-injection variables r p and d SOI . The influence of r p and d SOI on engine
performance is studied in the following sections.
143
Chapter 8. Pilot Injection
Table 8.1 The investigated operating points. The notion x : y : z indicates that
the corresponding parameter was swept from x to z in steps of y. Every combina-
tion of r p and d SOI was tested.
OP 1 OP 2
p IMEPg [bar] 5 10
Nspeed [rpm] 1200 1200
λ [-] 2.5 1.8
r EGR [-] 0.15 0.25
θ50 [CAD] 6 10
r p [-] 0 : 0.125 : 0.5 0 : 0.075 : 0.3
d SOI [CAD] 12.5 : 12.5 : 50 12.5 : 12.5 : 50
p rail [bar] 800 800
It would have been more convenient to use the injected fuel masses as input
variables as opposed to injection durations. Both for increased physical unders-
tanding and the more direct connection between the injected fuel amount and
the combustion processes. The reason for not doing so was due to the lack of
injector characteristics to compute the injection duration for a given demanded
fuel amount at the time of this study.
144
8.2 Experimental Characterization
20 20
40 40
d SOI [CAD]
d SOI [CAD]
30 15 30 15
20 20
10
10
0 0.2 0.4 0 0.1 0.2
r p [-] r p [-]
Figure 8.4 The relation between the pilot-injection variables and d p max at
p IMEPg = 5, 10 bar. It is clear that r p can be used to control d p max since d p max
decreased with r p . The d p max controllability increased when d SOI was reduced.
τm = θ10 − θSOI
m
[ms] (8.7)
which is presented in Fig. 8.5, where it can be seen that τm and d p max correlate.
The pilot effect weakened with increased d SOI . This could be explained by a
m
more dilute pilot mixture at θSOI , with a lowered temperature increase as a result.
An increased d SOI would also result in more diverse fuel-spray targets for the two
injections. Early pilot injections put more fuel in the crevice volume outside of
the piston bowl where the main-injection is targeted. The resulting spatial sepa-
ration of the injections could explain the observed trends. This was suggested by
optical engine data obtained from a similar heavy-duty engine setup, see [Lönn
et al., 2017]. Furthermore, the results in Figs. 8.4 and 8.5 indicate a significant
trade-off between obtainable τ and d p max , which could be problematic if a long
τ is required.
145
Chapter 8. Pilot Injection
2 2
40 40
d SOI [CAD]
d SOI [CAD]
30 1.5 30 1.5
20 20
1 1
Figure 8.5 The relation between the pilot-injection variables and τm at p IMEPg =
5, 10 bar. The main-injection ignition delay correlated with d p max in Fig. 8.4.
Indicated Efficiency
The pilot-injection impact on efficiency was evaluated by computing the gross
indicated efficiency
p IMEPg Vd
η GIE = (8.8)
m f Q LHV
where m f is the injected fuel amount, computed from the average fuel-flow du-
ring the experiment. The results presented in Fig. 8.6 show that η GIE had a shal-
low maximum when d SOI was short. This trend was more significant for OP 1,
where η GIE was more sensitive to the pilot configuration. Another visible trend is
that η GIE decreased when d SOI increased, and that this effect was stronger when
r p was higher. These results show that it could be efficient to have a pilot injec-
p
tion but that this effect is reversed for very early θSOI .
HC
The measured unburned hydrocarbon (HC) emission levels are presented in Fig.
8.7. HC emissions increased steeply as d SOI and r p were simultaneously incre-
ased. The explanation for this could be that the pilot fuel was injected into the
crevice regions outside of the piston bowl and did not burn completely due to
wall-cooling effects and too lean conditions. This effect is clearer at p IMEP = 5
bar where in-cylinder temperatures were lower. The decrease in combustion eff-
iciency indicated by the HC-emission increase could explain the η GIE decrease
observed in Fig. 8.6.
NOx
The measured NOx emission levels are presented in Fig. 8.8, where it can be
seen that NOx emissions mainly depended on r p and decreased with an in-
146
8.2 Experimental Characterization
40 0.52 40 0.52
d SOI [CAD]
d SOI [CAD]
30 0.5 30 0.5
20 20
0.48 0.48
Figure 8.6 The relation between the pilot-injection variables and η GIE at
p IMEPg = 5, 10 bar. The results show that η GIE had a shallow maximum when the
pilot was close to the main injection. Another visible trend is that η GIE decreased
when d SOI increased, and that this effect was stronger for larger r p .
Soot
Figure 8.9 shows that measured soot levels increased with r p and more so when
d SOI was small. An explanation for this is that the ignition delays for the fuel in-
jections were minimized at these operating points. A decreased air/fuel mixing
time results in richer combustion and increased soot formation. The increase in
soot emissions was not as large for longer d SOI , which means that d p max could
be reduced with a lower soot-emission penalty. Similar NOx and soot-emission
trends were presented in [Manente et al., 2009, 2010b].
Combustion-Timing Controllability
m p
The θ50 controllability was investigated by varying θSOI and θSOI in open loop at
m p
the operating points in Table 8.1. For each operating point, θSOI and θSOI were
varied individually in square-wave sequences with an amplitudes of 1 CAD and
m,p
a period of 25 cycles during 500 cycles. The partial derivatives ∂θ50 /∂θSOI were
then computed.
Computed partial derivatives in Figs. 8.10 and 8.11 indicate that θ50 is con-
m p
trolled by θSOI , and that the controllability decreased slightly with r p . The θSOI
p
effect was an order of magnitude smaller, where θSOI affected θ50 more when r p
147
Chapter 8. Pilot Injection
800 800
40 40
d SOI [CAD]
d SOI [CAD]
600 600
30 30
400 400
20 20
200 200
was increased. A hypothetical explanation is that the main injection initiated the
high-temperature reactions with corresponding heat-release rate. Similar trends
were shown in [Hasegawa and Yanagihara, 2003; Manente et al., 2009].
Another observation is that the pilot injection linearizes the relation between
m
θSOI and θ50 , due to its reduction of τ. This argument is supported by the fin-
m
dings in Chapter 6, where τ contributed to a nonlinear relation between θSOI and
θ50 . A pilot injection therefore increases the controllable region in Fig. 6.2, which
facilitates θ50 control.
Summary
The experimental findings show that r p can be used to control d p max . The choice
of r p gives a trade off between obtainable τ and d p max and a trade off between
soot and NOx emissions. A controller-design approach would be to let a fast
cycle-to-cycle controller decide r p to obtain acceptable d p max , and then let the
gas-exchange system set suitable boundary conditions in terms of EGR ratio and
intake temperature that simultaneously allow for low d p max and high τ.
The efficiency η GIE was shown to vary slightly with r p where η GIE decreased or
increased depending on d SOI . This is believed to be linked to the observed trends
in NOx and HC-emission levels as discussed previously. This hypothesis should
be confirmed with optical experiments.
The injection separation d SOI was shown to be a trade-off between high
d p max controllability and η GIE with shorter d SOI , and simultaneously low NOx
148
8.3 Controller Design
40 1,200 40 1,200
d SOI [CAD]
d SOI [CAD]
1,000 1,000
30 30
800 800
20 20
600 600
0 0.2 0.4 0 0.1 0.2
r p [-] r p [-]
Figure 8.8 The relation between the pilot-injection variables and NOx at
p IMEPg = 5, 10 bar. It can be seen that NOx mainly depended on r p and decre-
ased with an increased r p . A hypothetical explanation for this is the decrease in
τm with r p . The resulting decrease in heat-release rate and in-cylinder pressure
with τm indicate lowered peak temperatures and decreased NOx -formation rates.
m
and soot emissions with increased d SOI . The experiments also showed that θSOI
should be used to control θ50 .
149
Chapter 8. Pilot Injection
40 4 40
d SOI [CAD]
d SOI [CAD]
100
30 30
2 50
20 20
Figure 8.9 The relation between the pilot-injection variables and soot at
p IMEPg = 5, 10 bar. Note the difference in z-axis scale. Soot levels increased with r p
and more so when d SOI was small. An explanation for this is that the ignition de-
lays for the fuel injections were minimized at these operating points which decre-
ased air/fuel mixing time and resulted in richer combustion with increased soot
formation.
∂ θ 50 / ∂ θ m
SOI
[-], p IMEPg = 5 bar ∂ θ 50 / ∂ θ m
SOI
[-], p IMEPg = 10 bar
1.4 1.4
40 1.2 40
d SOI [CAD]
d SOI [CAD]
1.2
30 1 30
1
20 20
0.8
0.8
0 0.2 0.4 0 0.1 0.2 0.3
r p [-] r p [-]
was then formulated, where the assumption of a static relation between u and x
150
8.3 Controller Design
p p
∂ θ 50 / ∂ θ SOI [-], p IMEPg = 5 bar ∂ θ 50 / ∂ θ SOI [-], p IMEPg = 10 bar
0.2 0.2
40 0.1 40 0.1
d SOI [CAD]
d SOI [CAD]
30 0 30 0
−0.1 −0.1
20 20
−0.2 −0.2
0 0.2 0.4 0 0.1 0.2 0.3
r p [-] r p [-]
p
Figure 8.11 The computed gain from θSOI to θ50 at p IMEPg = 5, 10 bar. The
gain was low for all injection configurations. When comparing these trends with
m should be used to control θ .
Fig. 8.10, it is clear that θSOI 50
∂θ ∂θ50
50
∂θSOI ∂r p
1 −3
B = ∂d p max ∂d p max = −0.5 −25 (8.11)
∂θ ∂r p 0 1
SOI
0 1
From the experimental results in figures 8.10 and 8.11, it was found that
∂θ50 /∂θSOI = 1. The partial derivative ∂θ50 /∂r p = −3 was extracted from
Fig. 8.5. A decrease in θSOI resulted in an increased τm which in turn increa-
sed d p max , which explains ∂d p max /∂θSOI = −0.5. The d p max data in Fig. 8.4 gave
∂d p max /∂r p = −25. The matrix C in (8.10) is given by I 3x3 , since the first two
states can be computed directly using in-cylinder pressure measurements and
heat-release analysis. In order to incorporate model uncertainty and stochas-
tic cycle-to-cycle variation, zero-mean Gaussian processes v (k) and e(k) were
introduced with covariance matrices Q v and Q e
¡ ¢
Q v = E v (k)v (k)T
¡ ¢ (8.12)
Q e = E e(k)e(k)T
151
Chapter 8. Pilot Injection
ω1 ω2 ω3 ω4 ω5 ρ
100 0.01 25 6000 8000 1
Kalman Filter
State estimation was used to handle cycle-to-cycle variation in d p max and θ50 . In
order to estimate x from y, u and (8.10), a stationary Kalman filter was used. The
Kalman-filter state estimate x̂ was updated recursively according to
x̂(k + 1) = A x̂(k) + B u(k) + K (y(k) −C ŷ(k)) (8.14)
where the Kalman-filter gain K was given by the steady-state solution to Algori-
thm 1 in Chapter 3
K = APC T (Q e +C PC T )−1
(8.15)
P = AP A T +Q v − APC T (Q e +C PC T )−1C P A T
152
8.4 Controller Evaluation
were chosen to get sufficient measurement-noise attenuation. The choice Q v2,2 >
Q v1,1 reflects an increased model uncertainty in d p max . The Kalman filter only
considered the first two state equations in (8.10).
153
Chapter 8. Pilot Injection
d p max [bar/CAD]
14 15
12
θ 50 [CAD]
10 10
8 5
6
4 0
400 600 800 400 600 800
cycle [-] cycle [-]
−6 0.4
θ SOI [CAD]
−8 0.3
r p [-]
−10 0.2
−12 0.1
−14 0
400 600 800 400 600 800
cycle [-] cycle [-]
1.04 7.4
p IMEPg [bar]
1.02 7.2
[ms]
1 7
DOI
0.98 6.8
θ tot
0.96 6.6
0.94 6.4
400 600 800 400 600 800
cycle [-] cycle [-]
Figure 8.12 Input and output data during a sequence of θ50 r step changes. In
154
8.5 Discussion and Conclusions
d p max [bar/CAD]
12 20
15
θ 50 [CAD]
10
10
8
5
6 0
200 400 600 200 400 600
cycle [-] cycle [-]
−9 0.4
θ SOI [CAD]
0.3
−10
r p [-]
0.2
−11 0.1
0
200 400 600 200 400 600
cycle [-] cycle [-]
1.15 7.5
p IMEPg [bar]
[ms]
1.1 7
DOI
1.05 6.5
θ tot
1 6
200 400 600 200 400 600
cycle [-] cycle [-]
c
Figure 8.13 Input and output data during step changes in d p max . The response
c
time of d p max during a negative d p max step change was 2 cycles, while the re-
c
sponse time of d p max during a positive d p max step change was 10 cycles. The
m
main-injection timing θSOI adjusted for variations in θ50 caused by changes in
rp .
155
Chapter 8. Pilot Injection
d p max [bar/CAD]
12 14
10 12
θ 50 [CAD]
8 10
6 8
4 6
2 4
200 400 600 800 200 400 600 800
cycle [-] cycle [-]
−10 0.4
θ SOI [CAD]
−11 0.3
r p [-]
−12 0.2
−13 0.1
−14 0
200 400 600 800 200 400 600 800
cycle [-] cycle [-]
p IMEPg [bar]
1.2 8
DOI
[ms]
1 6
θ tot
0.8 4
200 400 600 800 200 400 600 800
cycle [-] cycle [-]
Figure 8.14 Input and output data during p IMEPg set-point changes. The re-
sponse time for p IMEPg was approximately 20 cycles. The pilot ratio was only
doing minor adjustments in order to keep d p max below the upper limit. This indi-
tot at this operating point.
cates that d p max was not very sensitive to changes in θDOI
Triple-injection strategies has previously been suggested for PPC operation [Ma-
nente et al., 2010a].
The injection separation d SOI was shown to impact both emissions and effi-
ciency, but was in this work held constant for simplicity. A suggested extension
of the controller presented here is therefore to incorporate d SOI as an input va-
riable in the controller design. This was for instance done in [Yang et al., 2017] to
control soot-emission levels.
156
8.5 Discussion and Conclusions
d p max [bar/CAD]
20
10
15
θ 50 [CAD]
8 10
5
6
0
600 700 800 900 600 700 800 900
cycle [-] cycle [-]
−8 0.4
θ SOI [CAD]
−10 0.3
r p [-]
−12 0.2
−14 0.1
−16 0
600 700 800 900 600 700 800 900
cycle [-] cycle [-]
1,600
Nspeed [rpm]
1,500
1,400
1,300
1,200
600 700 800 900
cycle [-]
Figure 8.15 Input and output data during Nspeed set-point changes. The con-
m in order to fulfill d p c
troller increased r p and decreased θSOI r
max and track θ50
when Nspeed was increased. It can also be seen that the d p max -noise level de-
creased with Nspeed .
157
9
Low-Load Control
9.1 Introduction
To obtain reliable low-load operation is one of the major challenges with ga-
soline PPC, both due to high cycle-to-cycle variation and low combustion ef-
ficiency. Advantages with a high octane number (ON) have mainly been obse-
rved at mid-to-high engine load. With a high ON, the in-cylinder conditions at
low-load operation result in too long ignition delays for combustion to complete
during the closed part of the engine cycle [Fieweger et al., 1997]. The issue of
having too long ignition delays was discussed in Chapter 6, where elevated in-
take temperatures were needed to achieve θ50 close to TDC, see Fig. 6.2.
Incomplete combustion and misfire result in reduced engine efficiency and
increased HC and CO emission levels. Homogeneous-reactor simulations have
shown that high levels of HC and CO are a result of lean mixtures at too low tem-
perature [Kim et al., 2008]. Fuel reactivity is also an important factor. Manente et
al. [2010a] showed that as the fuel ON increased from 70 to 100, the PPC low-load
limit, defined as the p IMEP at the fuel-ignitability limit, increased from 3 to 15 bar
at atmospheric intake conditions.
Results in [Kalghatgi et al., 2006] and [Weall and Collings, 2009] showed that
a remedy to the PPC low-load issue is to increase in-cylinder temperature and
equivalence-ratio stratification levels, since combustion initiated in rich regions
with high temperature aids combustion of the overall cylinder charge. Previous
studies presented in [Borgqvist et al., 2012; 2013] aimed to improve the gaso-
line PPC low-load efficiency using variable-valve actuation. These results showed
that negative valve overlap and re-breathing in combination with a split-main in-
jection strategy were able to increase the low-load performance by increasing the
temperature of trapped residual gases. In [Solaka et al., 2012], the low-load limit
of a single-cylinder light-duty engine was extended down to p IMEPg = 2 bar, using
boosted intake air. The absolute intake pressure needed at this load was appro-
ximately 2 bar for the fuels with the highest ON (88.6 and 87.1). Other actuator
options for improved low-load operation include variable compression ratio [Ha-
raldsson et al., 2002] and fast thermal management [Martinez-Frias et al., 2000].
158
9.2 Low-Load Experiments
This chapter investigates how fuel injection and intake conditions should be
chosen for improved PPC efficiency, with the objective of extending the operation
range towards lower loads. Section 9.2 presents experimental data for p IMEPn ∈
[1, 5] bar. This data were then used to suggest a simple controller that acts to
minimize the ignition delay τ during low-load operation, see Sec. 9.3.
Section 9.4 shows how the controller-design choices affect fuel consumption
during transient operation.
Injection Timing
The combustion timing θ50 has to be phased shortly after TDC for work ouput
to be maximized and temperatures to be sufficient for the chemical reactions to
complete. Experimental θ50 data are presented in Fig. 9.1 as a function of θSOI for
three different fuel-injection durations θDOI and p rail = 800 bar.
The gain from θSOI to θ50 was positive for θSOI > −20 CAD. As θSOI decrea-
sed from this point, the increase in τ with θSOI became more significant, and the
gain from θSOI to θ50 decreased. The gain even became slightly negative for very
early θSOI , as the increase in τ exceeded the decrease in θSOI . Similar θ50 / θSOI
characteristics were observed in Chapter 6, Fig. 6.2, with the difference that the τ
model in Chapter 6 was unable to capture the negative-gain region for early θSOI .
A physical explanation for a negative gain was given in [Kalghatgi et al., 2006],
where similar experimental results were presented: An increase in τ gives more
time for fuel-air mixing, which in turn gives leaner mixtures and even longer τ.
The model in Chapter 6 did not include the time-resolved equivalence ratio his-
tory of the fuel/air mixture when computing the accumulated reactivity, and was
therefore not able to capture the negative-gain effect.
The data in Fig. 9.1 provide a lower bound for which θSOI can be used to ef-
fectively control θ50 with a linear controller. With integral action and a too low
infeasible set point, the controller would keep advancing θSOI and instead in-
crease set-point deviation. It was therefore decided to saturate θSOI at -20 CAD
to maintain a positive gain. In this way, the controller could obtain the wanted
set-point or the earliest obtainable θ50 , without risking unnecessarily long τ.
Pilot Injection
It was previously shown in Chapter 8 that a pilot injection reduces τ of the main
injection. A decrease in main-injection τ makes it possible to obtain a more ad-
vanced θ50 . This pilot-injection effect has previously been observed for conven-
tional diesel combustion [Osuka et al., 1994; Macmillan et al., 2009].
159
Chapter 9. Low-Load Control
20
θ DOI = 0.65 [ms]
θ DOI = 0.70 [ms]
θ DOI = 0.75 [ms]
15
θ 50 [CAD]
10
Figure 9.1 Combustion timing θ50 as a function of fuel-injection timing θSOI for
three different fuel-injection durations θDOI . For injections closer to TDC, the gain
between θSOI and θ50 was positive. As θSOI was advanced, the gain decreased and
became slightly negative. The negative gain was more significant for the shorter
injection durations.
p
Experiments with different pilot- and main-fuel injection durations, θDOI and
m p
θDOI , were conducted to investigate the θDOI effect on the gross indicated effi-
ciency
p IMEPg Vd
η GIE = (9.1)
m f Q LHV
The combustion timing θ50 was kept as close to TDC as possible, whilst the pi-
lot injection was positioned 10 CAD prior to the main injection. A small separa-
tion between pilot and main was previously shown to minimize τ, see Fig. 8.5.
m
Level curves of p IMEPn (blue, dashed) and η GIE (red, solid) as a function of θDOI
p
and θDOI are presented in Fig. 9.2. Figure 9.2 shows that for a given p IMEPn , η GIE
p
could be increased by having a longer θDOI . This effect became less significant
for higher p IMEPn . The pilot-injection effect on heat-release rate is presented in
Fig. 9.3, where it can be seen that the pilot injection both reduced τ and increased
the heat-release rate. From these results it was concluded that a pilot injection
should be used to aid main-injection ignition and allow for a more advanced θ50 .
160
9.2 Low-Load Experiments
0.4
0. 6
0.38
2. 5
2
0.
5
3
0.35 0.5
5
0.33 0.45
θ DOI [ms]
0.4
1
1.
2. 5
5
0.3
2
0.35
p
0.28 0. 4
5
0.3
0.
0.25
4
0.5
0.2
1.5
0.23 0.5 p IMEPn
2.5
1
2
ηGIE
0.2
0.6 0.63 0.65 0.68 0.7 0.73 0.75 0.78 0.8
θm
DOI
[ms]
Figure 9.2 Level curves of p IMEPn (blue, dashed) and η GIE (red, solid) as a func-
m and θ p . For a low p
tion of θDOI DOI IMEPn , η GIE could be increased by having a
p
longer θDOI . This effect became less significant for higher p IMEPn .
80
60
p [bar]
40
20
0
−30 −20 −10 0 10 20 30
θ [CAD]
Figure 9.3 In-cylinder data comparing heat-release and pressure with and wi-
thout a pilot injection. With the same main-injection timing, τ was reduced by the
pilot injection, which gave an advanced θ50 and an increased heat-release rate.
161
Chapter 9. Low-Load Control
ηGIE [-]
0.3
0.4
0.25
0.2 0.35
φ [-]
0.15
0.3
0.1
0.25
Intake Conditions
The gas-system valve positions θHP and θcool were varied at the following injec-
tion durations
m
θDOI = {0.65, 0.7, 0.75, 0.85, 1.05} [ms] (9.2)
to investigate how the intake conditions affect η GIE . The combustion timing was
kept in the interval 0-5 CAD and no pilot injection was used. Computed η GIE du-
ring these experiments are presented as a function of global equivalence ratio φ
and temperature at θSOI , TθSOI in Fig. 9.4. The temperature TθSOI was computed
by assuming adiabatic compression
µ ¶
V (θIVC ) γ−1
TθSOI = Tin (9.3)
V (θSOI )
In Fig. 9.4, the efficiency η GIE decreased steeply when φ and TθSOI were simulta-
neously reduced. It can also be seen that η GIE was more sensitive to TθSOI at low
φ, and could be increased significantly by increasing TθSOI . The effect of increa-
sing Tin can also be seen in Fig. 9.5, where the increase in Tin decreased τ and
increased the heat-release rate.
162
9.2 Low-Load Experiments
80
Tin = 317 K
Tin = 293 K
p [bar] 60
40
20
0
−30 −20 −10 0 10 20 30
θ [CAD]
Figure 9.5 The effect of elevating the intake temperature. This was done by
opening the high-pressure EGR valve and closing the cool thermal-management
valve. An increased temperature gave a decrease in τ as well as an increased
m
heat-release rate, which indicates an improved combustion efficiency since θDOI
was kept constant.
Varying θHP and θcool also affected pumping losses p PMEP , which is the in-
dicated mean-effective pressure during the gas-exchange strokes. The relation
between p PMEP and the intake and exhaust manifold pressures, p in , p ex , is pre-
sented in Fig. 9.6, where the symbols ° and 4 relate the data in Figs. 9.4 and 9.6.
Pumping losses correlated with p in − p ex as expected. Moreover, the symbols in-
dicate that p PMEP in Fig. 9.6 correlated with TSOI in Fig. 9.4. This was because ope-
ning θHP not only heated the intake charge, it also elevated the intake-manifold
pressure, and in that way reduced pumping losses. The data presented in Figs. 9.4
and 9.6 were used for efficiency optimization through simulation in the following
section.
163
Chapter 9. Low-Load Control
p PMEP [bar]
−0.24
1.28
−0.26
1.26
−0.28
1.24 −0.3
p ex [bar]
1.22 −0.32
−0.34
1.2
−0.36
1.18
−0.38
1.16
0.96 0.98 1 1.02 1.04 1.06 1.08 1.1 1.12 1.14 1.16 1.18
p in [bar]
∆θcool
ẋ = f (x, θcool , θHP ), (2.22) − (2.37)
η GIE = g 1 (φ, TθSOI )
p PMEP = g 2 (p 1 , p 2 )
164
9.2 Low-Load Experiments
sented by the first cost term in (9.4). The other two terms penalize deviation from
suitable valve positions at mid-to-high engine load where intake-manifold tem-
perature should be kept low to reduce heat-transfer losses and the low-pressure
EGR path is preferred over the high-pressure EGR path. Valve positions θcool and
θHP are here normalized from 0 to 1, where 1 denotes fully open. The cost func-
tion is defined over a prediction horizon of H p engine cycles, where k denotes
cycle index. The cost function was minimized subject to the dynamics of the
gas-exchange system model in Sec. 2.3. Here, the state x include the pressures
and temperatures of the volumes in Fig. 2.3. The functions g 1 and g 2 were obtai-
ned by interpolating the data in Figs. 9.4 and 9.6.
The gas-exchange model in Sec. 2.3 was then used to simulate the MPC in
(9.4). During this experiment, p IMEPn was set to follow a set-point trajectory using
a PI controller. Combustion timing θ50 was also regulated through θSOI adjust-
ments. No combustion model was used since p IMEPn was given directly by g 1 and
g 2 . An ignition-delay model was however used to compute TθSOI . The nonlinear
MPC problem (9.4) was solved every simulated engine cycle using the MATLAB
nonlinear-optimization toolbox together with the ode23s solver to compute the
gas-system model output. Global solutions of (9.4) can not be guaranteed due to
the nonlinearity of the optimization problem. Solutions obtained were however
justified by confirming that they were consistent for different initial guesses.
Simulation results for three MPCs with different fuel-consumption weights
ω1 = 0, 1, 5, ω2 = ω3 = 10 and H p = 10 are presented in Fig. 9.7. As the cost for fuel
consumption was increased, the controller increased gas-exchange actuation by
closing θcool and opening θHP when p IMEPn was reduced. The controller was in
this way able to avoid the low-efficiency region in the φ − T diagram (see the
lower right subdiagram in Fig. 9.7). This lowered the injected fuel amount m f
needed at low load. Overshoot in p IMEPn was also slightly reduced.
Fuel efficient actuation of θcool and θHP in Fig. 9.7 was approximated with a
static φ-feedback law K (φ) µ ¶ µ ¶
θcool K 1 (φ)
= (9.5)
θHP K 2 (φ)
Figure 9.8 shows K coincidence with the actuated valve positions in Fig. 9.7.
There is some deviation between data and K 1 for higher φ. This deviation was
deliberately chosen to utilize the full range of θcool for reduced intake-manifold
temperatures at higher engine loads. The reason for approximating the MPC be-
havior was due to the inhibiting computation times needed to solve (9.4) online.
165
Chapter 9. Low-Load Control
80 6
p IMEPn [bar]
60
4
m f [mg] 40
2
20
0 0
0 2 4 6 8 0 2 4 6 8
time [s] time [s]
0.3
1
0.2
θ cool [-]
θ HP [-]
0.5
0.1
0 0
0 2 4 6 8 0 2 4 6 8
time [s] time [s]
800 0.5
0.4
TθSOI [K]
750 0.3
φ [-]
0.2
700 0.1
0
0 2 4 6 8 680 700 720 740 760 780
time [s] TSOI [K]
ω1 = 0 ω1 = 1 ω1 = 5
166
9.3 Suggested Controller
90
K 1 (φ)
K 2 (φ)
75 ∗
θcool
∗
θHP
60
θ HP / θ cool [deg]
45
30
15
0
0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 0.28 0.3 0.32
φ [-]
Figure 9.8 The static feedback law K was used to approximate efficient MPC va-
lve positions in Fig. 9.7. It was decided to use the full range of θcool which gave a
∗
deviation between θcool and K 1 for higher φ.
• A pilot injection should be used to reduce τ for advanced θ50 and increased
heat-release rates. It was here decided to use a small pilot, located 10 CAD
prior to the main-injection. The reason the short separation time was to
minimize τ, see Fig. 8.5.
The suggested feedback loop is presented in Fig. 9.9. The controller is fairly sim-
ple, and the decentralized controller design is easy to implement. The controller
would nevertheless be able to maximize efficiency according to the presented ex-
perimental results. The design of a centralized controller structure, that for exam-
ple utilizes ignition-delay and gas-exchange models, as presented in Chapter 7,
is suggested future work. Another possible extension to the controller in Fig. 9.9
is a method that experimentally identifies the θ50 gain, in order to adapt the θSOI
saturation limit.
167
Chapter 9. Low-Load Control
−1
θSOI
r
θ50 P
PI
θ50
Engine
θHP φ
θcool
K (φ)
Figure 9.9 The experimental results suggested the following low-load con-
trol strategy: A range-limited θ50 PI controller in combination with a static
gas-exchange controller with global φ as feedback variable. Load and rail pressure
are controlled separately using PI controllers with feedforward.
For all cases, p IMEPn was controlled by keeping the rail pressure constant and
varying the main-injection duration with PI controllers and feedforward.
Experimental test-cycle results from one cylinder are presented in Fig. 9.10
for cases (1-4). Differences between the cases are more noticeable at low engine
load. The injection timing was not limited in case (1). This led to very early in-
r
jection timings when the earliest obtainable θ50 was larger than θ50 . The early
injection timing in case (1) resulted in a τ increase and a delayed θ50 as a result.
r
When limiting θSOI , as in case (2), θ50 was advanced and the deviation from θ50
was reduced.
168
9.5 Conclusions
When a pilot injection was introduced in case (3), τ decreased, see the θSOI
and θ50 subdiagrams in Fig. 9.10. The reduction in τ allowed the controller to
r
keep θ50 closer to θ50 . In case (4), τ was further decreased as Tin and TθSOI in-
creased due to valve-position actuation. This resulted in smaller θ50 error at low
load.
In-cylinder data showing these trends more clearly are presented in Fig. 9.11
where in-cylinder pressure, injection current and heat-release rate are presented
at cycle 250. It can be seen that the gradual controller adjustments led to decrea-
sed τ. The heat-release rate also differed for the different cases where case 1, with
the longest ignition delay had the lowest heat-release rate.
Accumulated fuel consumption, computed from injection durations is pre-
sented in Fig. 9.12. The fuel-consumption rates differed more clearly at the
low-load operating points. Figure 9.12 shows that fuel consumption decreased
from case (1-4), with reductions from 2 to 9 %. The greatest reduction resulted
from the introduction of a pilot injection.
9.5 Conclusions
This chapter presented a control strategy for improved PPC performance at low
load. The combustion timing θ50 was only controllable with respect to θSOI in a
specific interval, see Fig. 9.1. It was therefore important keep θSOI in this interval
to maintain closed-loop stability with a θ50 -feedback controller.
A pilot injection was shown to increase the combustion efficiency at low load,
see Figs. 9.2 and 9.12. This was due to the decrease in τ, which led to an increased
heat-release rate and the possibility to further advance θ50 .
The problem of maximizing η GIE was formulated as an optimization problem
in the φ-T diagram, where η GIE could be increased by heating the inducted air
charge, see Fig. 9.4. Control of the intake conditions was achieved by varying θHP
and θcool according to a feedback law K (φ), obtained from MPC simulations.
An experimental controller evaluation showed that these findings could re-
duce fuel consumption from 2 to 9 %, see Fig. 9.12.
169
Chapter 9. Low-Load Control
1.2
p IMEPn [bar]
5
θ DOI [ms]
4 1
3 0.8
2 0.6
1 0.4
0
500 1,000 1,500 500 1,000 1,500
cycle [-] cycle [-]
20
−10
θ SOI [CAD]
15
θ 50 [CAD]
10 −20
5
−30
0
500 1,000 1,500 500 1,000 1,500
cycle [-] cycle [-]
25
20 310
θ HP [deg]
Tin [K]
15
300
10
5 290
0
500 1,000 1,500 500 1,000 1,500
cycle [-] cycle [-]
100 0.75
θ cool [deg]
75
0.5
φ [-]
50
0.25
25
0 0
500 1,000 1,500 500 1,000 1,500
cycle [-] cycle [-]
Figure 9.10 Experimental results for the four different controller cases. In case 1,
θSOI was not limited to the positive-gain region. This led to very early θSOI , and as
a result, a delayed θ50 (compare case 1 and case 2). With a pilot injection in case
3, τ decreased, see the θSOI and θ50 subdiagrams. Finally, in case 4, τ decreased
further as θHP was opened and θcool was closed.
170
9.5 Conclusions
p IMEPn = 2 bar
80
Case 1
Case 2
Case 3
60 Case 4
p [bar]
40
20
0
−30 −25 −20 −15 −10 −5 0 5 10 15 20
θ [CAD]
Figure 9.11 In-cylinder data from cycle 950 in Fig. 9.10. The gradual controller
adjustments from case 1 to case 4 resulted in a decreased τ. It can also be seen
that the advanced θSOI in case (1) gave a reduced heat-release rate.
171
Chapter 9. Low-Load Control
·105
3.5
Case 1
3 Case 2 (−2%)
Case 3 (−8%)
Case 4 (−9%)
Fuel Consumption [mg]
2.5
1.5
0.5
0
0 200 400 600 800 1,000 1,200 1,400 1,600 1,800
cycle [-]
Figure 9.12 Accumulated fuel consumption for the four controller cases in
Fig. 9.10 where the fuel-consumption rate differed more at low load. Fuel con-
sumption decreased from case (1-4), with reductions from 2 to 9 % where the gre-
atest reduction came with the introduction of a pilot injection. Here, the injected
fuel mass was computed from fuel-injection durations. The reason for not using
the fuel-mass-flow meter was that the meter was mounted far from the engine,
and was therefore not reliable in transient operation.
172
10
Constraint Handling with
Multiple Injections
10.1 Introduction
Despite its lack of spatial information, the heat-release rate as a function of
crank-angle degree is an important variable when maximizing the thermody-
namic engine efficiency. When assuming a zero-dimensional, adiabatic model
with constant γ, the thermodynamic efficiency is maximized when heat is re-
leased instantaneously at TDC. This heat-release rate minimizes exhaust losses
without generating counterproductive pressure during the compression stroke,
and is equivalent to the ideal Otto cycle with constant-volume combustion at
TDC.
The optimal combustion timing is delayed to after TDC when heat-transfer
to cylinder walls is also accounted for. A delayed combustion timing reduces
in-cylinder temperature and heat-transfer losses. It also reduces peak in-cylinder
pressure which in turn reduces engine friction and heat-losses to crevice volu-
mes. A drawback with delaying the combustion timing is the resulting increase
in exhaust-gas energy. The optimal combustion timing is therefore a compromise
between these losses.
The optimal heat-release rate becomes more involved when constraint ful-
fillment is required. With constraints on maximum cylinder pressure and NOx
formation, motivated by mechanical tolerances and emission regulations, com-
bustion has to be delayed during the expansion stroke to reduce the in-cylinder
pressure and temperature. Furthermore, efficient aftertreatment-system perfor-
mance requires sufficient exhaust-gas temperature, which in turn demands in-
creased exhaust losses [Gieshoff et al., 2000; Katare et al., 2009]. As suggested by
these examples, there are also compromises between constraint fulfillment and
the thermodynamic efficiency.
Optimal heat-release rates with respect to constraints on in-cylinder pres-
sure, NOx formation and knock intensity were computed in [Eriksson and Sivert-
173
Chapter 10. Constraint Handling with Multiple Injections
sson, 2016; Guardiola et al., 2017]. In these studies, optimal heat-release rates
were found to be multimodal to reduce in-cylinder pressure and temperature.
With a single fuel injection and a production heavy-duty fuel-injection system,
however, the heat-release rate controllability is limited. The combustion timing
can be controlled more or less freely by varying the injection timing, but the
heat-release shape depends on rates of fuel injection, mixing and chemical reac-
tions. These rates can only be controlled partially by adjusting cylinder-mixture
properties and the fuel-injection pressure.
The degrees of freedom increase when multiple injections are used. It
was found in [Han et al., 1996] that two injections can be used to provide a
more distributed heat-release rate with reduced peak temperature and NOx
formation. This allows for a more advanced effective combustion timing
without NOx -constraint violation. Okamoto and Uchida [2016] presented a
multiple-injector strategy for heat-release shaping. Experimental results showed
that heat-release shaping could provide a 75 % reduction of NOx emissions
due to suppressed peak average temperature and pressure with maintained in-
dicated efficiency. An indicated specific-fuel consumption reduction of 12 %
was reported in [Dober et al., 2008] with maintained NOx emissions, using a
novel fuel-injection system and multiple injections. The use of post injections
for control of exhaust-gas temperature was demonstrated in [Zheng et al., 2005;
Castellano et al., 2013] among others.
Some insight to the potential efficiency benefit with multiple injections can
be gained by studying the simulated pressure curves in Fig. 10.1. Here, the dif-
ferent heat-release rates and injection configurations generate the same work
output, p IMEPn = 15 bar. Without a constraint on maximum pressure p max , the
optimal combustion timing is close to TDC with the gross indicated efficiency
η GIE = 0.5, see the upper subdiagram. When a p max constraint is imposed (lower
subdiagram), the combustion timing has to be delayed for constraint fulfillment.
This results in a 4.8 % η GIE decrease with one injection (blue). With two injec-
tions (red), however, the combustion timing is allowed to be advanced, and the
decrease in η GIE is only 0.8 %.
An extensive calibration effort is demanded in order to find efficiency-optimal
fuel-injection configurations that fulfill pressure, NOx and exhaust-temperature
constraints over the engine operating range. Optimal configurations are also sen-
sitive to hardware aging, fuel properties and variation in operating conditions.
Moreover, the engine should operate as close to the constraint as possible for
maximized engine efficiency.
This chapter therefore investigates the use of feedback control for automa-
tic fuel-injection adjustment and increased engine efficiency subject to specified
constraints. The suggested controller is a hybrid, multiple-input multiple-output
PI controller that utilizes feedback from in-cylinder pressure-sensor measure-
ments, an NOx -emission model functioning as a virtual sensor, and measured
exhaust temperature. The controller varies the number of injections and adjusts
174
10.1 Introduction
Unconstrained p max
200
100
50
0
−20 −10 0 10 20 30 40 50
θ [CAD]
100
50
0
−20 −10 0 10 20 30 40 50
θ [CAD]
Figure 10.1 Simulated pressure curves with and without constraints on p max .
Without a constraint on p max , the optimal combustion timing is close to TDC with
the gross indicated efficiency η GIE = 0.5, see the upper subdiagram. When a p max
constraint at 125 bar is imposed, the combustion timing has to be delayed (see
lower subdiagram). This results in a 4.8 % η GIE decrease, with one injection (1).
With two injections (2), the combustion timing can be advanced, and the decrease
in η GIE is only 0.8 %.
175
Chapter 10. Constraint Handling with Multiple Injections
with simulation results in Sec. 10.3. The proposed controller design is then intro-
duced in Sec. 10.4. Experimental evaluation results are presented in Sec. 10.5,
where both transient and steady-state operation are evaluated. Controller per-
formance is also compared to that of a simpler single-injection controller. Finally,
discussion and conclusions are given in Secs. 10.6 and 10.7.
p IMEPg Vd
η GIE = (10.1)
m f Q LHV
The efficiency η GIE should be maximized whilst p max , formed NOx emissions and
the exhaust temperature Tex fulfill upper and lower bounds
c
p max ≤ p max
NOx ≤ NOcx (10.2)
c
Tex ≥ Tex
Upper bounds on p, d p/d θ and NOx are motivated by mechanical engine tole-
rances and legislated emission limits. The lower limit for Tex is introduced to gua-
r
rantee after-treatment system performance. The demanded work output p IMEPn
should also be delivered
r
p IMEPn = p IMEPn (10.3)
The case with up to two fuel injections is first considered. The optimization va-
1 2
riables are the injection timings, denoted θSOI , θSOI , and injected fuel masses, m 1f
and m 2f . The optimization problem can be simplified by assuming that the total
fuel mass m tot
f
r
is determined by p IMEPn . The ratio
m 1f
r= (10.4)
m tot
f
can then be optimized instead of m 1f and m 2f . This ratio was here limited to
0.5 ≤ r ≤ 1 (10.5)
1 2
to exclude redundant configurations. Injection durations θDOI and θDOI are de-
tot
termined by r , m f and p rail , with the use of an injector map M inj
¡ 1 ¢T
θDOI 2
θDOI = M inj (r, m tot
f , p rail ) (10.6)
176
10.3 Simulation
The injection pulses are not allowed to overlap in order to ensure that the inj-
ector needle closes in-between injections. This imposes constraints on possible
fuel-injection timings and durations.
The optimization variable
¡ m p ¢T
u = θSOI θSOI r Q tot (10.7)
y = f (u) (10.9)
10.3 Simulation
Simulation experiments with two combustion timings were conducted to find
optimal u with respect to (10.10). The objective was to investigate when multiple
injections are beneficial and how optimal combustion timings are configured as
a function of imposed constraints. This section presents the model used, and
simulation results showing trends for p max , NOx , temperature at exhaust-valve
opening, TEVO , and η GIE as u is varied.
177
Chapter 10. Constraint Handling with Multiple Injections
Model
The model used was the zero-dimensional model presented in Chapter 2
µ ¶
dp γ dV γ − 1 dQ c dQ ht
=− p+ − , p(θIVC ) = p in (10.11)
dθ V dθ V dθ dθ
i
Here, θSOC was related to the i :th injection timing θSOI with the use of an Arrhe-
nius ignition-delay expression
i i
θign = θSOI + τi
(10.13)
τi = Ap −n i e E a /R̃TSOIi
SOI
where ignition delay τi was omitted in (10.13) if combustion had started prior to
injection.
NOx formation was computed with the two-zone, Zeldovich-mechanism mo-
del presented in [Egnell, 2001]. A more detailed description of this model is given
in Chapter 2.
Instead of modeling Tex , it was decided to study the in-cylinder temperature
at exhaust-valve opening TEVO , which correlates with Tex .
Simulated Conditions
Simulation experiments were conducted with two injections and the correspon-
ding accumulated heat-release
where Q c1 and Q c2 are Wiebe expressions on the form of (10.12) and r is defined by
(10.5). The combustion timings of Q c1 and Q c2 , θCT
1 2
and θCT were swept for diffe-
rent r and constant total fuel energy. The reason for limiting the study to two in-
jections and not considering a more general heat-release was that this would re-
semble a realistic scenario that could be realized in the experimental setup used.
Constraints on p max and NOx were evaluated at a higher load compared to TEVO ,
since lower bounds on TEVO are more likely to become active at low load. Model
parameters used are presented in Table. 10.1.
178
10.3 Simulation
Table 10.1 Model parameters used in simulation. Constraints on p max and NOx
were evaluated at a higher load compared to TEVO .
TEVO constraints
All constraints
Simulation Results
1 2
Level curves for p max , NOx , TEVO and η GIE as a function of θCT and θCT are pre-
sented in Figs. 10.2, 10.3 and 10.4. Results with r = 0.5, 0.625, 0.75 and 0.875 are
presented in Fig 10.2, whereas results with r = 0.5 are presented in Figs. 10.3 and
10.4. In these figures, the black lines are η GIE level curves, with the most efficient
combustion timing, marked ×, found close to TDC. The efficiency then decreases
179
Chapter 10. Constraint Handling with Multiple Injections
1 2
with delayed θCT and θCT . The colored lines indicate p max [bar], NOx [ppm] and
TEVO [K] level curves, where the colored symbols indicate η GIE -optimal combus-
tion timings with respect to the different level curves as constraints, for both
1 2 1 2
θCT > θCT (∗) and θCT < θCT (¦). The gray regions indicate combustion-timings
where fuel-injection pulses overlap, i.e., u ∉ U.
Figures 10.2-10.4 show that it is more efficient with two combustion timings
when subject to p max , NOx and TEVO constraints. This is recognized by observing
1 2
that optimal timings occur away from the θCT = θCT line which indicate a unimo-
dal, single-injection heat release rate. This holds even if u ∉ U would have been
allowed.
1 2
The symmetry with respect to θCT = θCT in Figs. 10.3 and 10.4 is due to the
heat-release rates being identical with r = 0.5. The symmetry was altered as r
was varied in Fig. 10.2. In this figure, it can be seen that η GIE was maximized for
1 2 1 2
both θCT > θCT and θCT < θCT . In general, there was no trend favoring any of these
two configurations.
Trade-offs between η GIE and x c for different r are presented in Fig. 10.5. The
p max and NOx trade-offs can be improved by up to 4 % when using two combus-
tion timings instead of one (r = 1), whereas the TEVO trade-off can be improved
c
by up to 2 %. The potential η GIE advantage with two injections increases as p max
c c
and NOx become more conservative. This does not hold for TEVO . The trade-offs
in Fig. 10.5 do not distinguish any clear choice for r < 1.
An attempt to explain the observed trends is made in Fig. 10.6. This figure
presents efficiency-optimal combustion timings for one and two injections with
c
arbitrary p max (upper), NOcx (middle) and TEVO c
(lower). For p max and NOx con-
straints, two injections give a more distributed heat-release. This lowers the
peak pressure and temperature which gives a slower NOx -formation rate. Similar
trends were observed experimentally in [Han et al., 1996]. For the TEVO constraint
in the lower subdiagram, the late injection gives a sufficient contribution to TEVO
for the first combustion timing to be timed optimally. To summarize: the overall
trend in Fig. 10.6 is that two injections allow for an effective or mean combustion
timing closer to TDC, which increases the indicated efficiency η2GIE > η1GIE .
180
10.4 Controller Design
r = 0.875 r = 0.75
150
130
110
130
150
110
20 η = 0.513 20 η = 0.518
η = 0.527
θ 2CT [CAD]
θ 2CT [CAD]
η = 0.531
η = 0.538
η = 0.543
10 10
r = 0.625 r = 0.5
130
110
130
110
20 20
η = 0.523 η = 0.526
θ 2CT [CAD]
θ 2CT [CAD]
η = 0.538
η = 0.537
10 η = 0.542 10 η = 0.541
150 150 η = 0.526
Figure 10.2 Level curves of η GIE (black) for different combustion timings θCT 1 ,
2 , and ratios r . Level curves for p
θCT max are presented in red (110 bar), blue (130
bar) and green (150 bar). The shaded gray areas correspond to infeasible injection
timings where the injection pulses overlap. For each r , the most efficient timings
are marked ×. The colored marks indicate the most efficient feasible points, given
1 > θ 2 and θ 1 <
the different p max constraints for the feasible regions where θCT CT CT
2 . The figure shows that two combustion timings are optimal when subject to
θCT
p max constraints. This is recognized by observing that optimal timings occur away
from the θCT1 = θ 2 line which indicate a unimodal, single-injection heat release
CT
rate. Pressure curves corresponding to the marked high-efficiency points for r =
0.5 are presented in Fig. 10.7.
181
Chapter 10. Constraint Handling with Multiple Injections
30
θ 2CT [CAD]
20
10 0
0
10
1400 η = 0.518
η = 0.545 600
η = 0.535
η = 0.54
0
0 5 10 15 20 25 30 35 40
θ 1CT [CAD]
Figure 10.3 Simulated η GIE level curves (black) as a function of θCT1 and θ 2 .
CT
Level curves for NOx are presented in red (600 ppm), blue (1000 ppm) and green
(1400 ppm). The colored marks indicate η GIE -optimal combustion timings with
respect to the different NOx constraints. NOx -formation curves corresponding to
the marked high-efficiency points are presented in Fig. 10.8.
182
10.4 Controller Design
η = 0.498
30
θ 2CT [CAD]
η = 0.519
20 850
93
0
10
89
0
η = 0.538
0η
0 5 10 15 20 25 30 35 40
θ 1CT [CAD]
Figure 10.4 Simulated η GIE level curves as a function of θCT1 and θ 2 . Level cu-
CT
rves for TEVO are presented in red (930 K), blue (890 K) and green (850 K). It
is optimal to keep one injection close the optimal point ×, and then delay the
latter to fulfill the constraint. Temperature curves corresponding to the marked
high-efficiency points are presented in Fig. 10.9.
The indices 1, 2 and 3 represent the pilot, main and post injection, respectively.
The controller should also deliver the desired load output, which was handled by
adjusting the total fuel mass m tot
f
.
A simple controller design is attractive from an implementation perspective.
Therefore, a hybrid multiple-input multiple-output PI controller
with gains k p and k I was designed to achieve the desired system behavior. The
suggested controller used combustion timings, peak pressure levels, the modeled
NOx -emission level and measured exhaust temperature as feedback signals. This
controller is presented in the remainder of this section.
Combustion Detection
When controlling multiple combustion timings, it is no longer sufficient to use
the crank-angle of 50 % burnt as feedback variable. This quantity might not be re-
lated to a physical combustion timing with a multimodal heat-release rate. The
183
Chapter 10. Constraint Handling with Multiple Injections
0.5 0.5
ηGIE [-]
ηGIE [-]
r = 0.5 r = 0.5
0.45 0.45
r = 0.625 r = 0.625
r = 0.75 r = 0.75
r = 1.0 r =1
0.4 0.4
500 1,000 1,500 2,000 80 100 120 140
NOcx [ppm] cp max [bar]
c
Trade-off Between TEVO & ηGIE
0.55
r = 0.5
r = 0.625
r = 0.75
0.5 r = 1.0
ηGIE [-]
0.45
0.4
900 1,000
cTEVO [K]
184
10.4 Controller Design
200
One Injection
Two Injections
150 η1GIE = 0.476
p [bar]
50
0
−20 −10 0 10 20 30 40 50
θ [CAD]
2,000
η1GIE = 0.49 p
NOx
1,500
η2GIE = 0.50
NOx [ppm]
1,000
500
0
−20 −10 0 10 20 30 40 50
θ [CAD]
2,000
η1GIE = 0.48 p
T
1,500 η2GIE = 0.49
T [K]
1,000
500
0
−20 0 20 40 60 80 100 120 140
θ [CAD]
Figure 10.6 A physical explanation for the difference between one (1) and two (2)
c
injections with arbitrary p max (upper), NOcx (middle) and TEVO
c (lower). For p max
and NOx constraints, two injections give a more distributed heat-release. A di-
stributed heat release lowers the peak pressure and gives a slower NOx -formation
rate, which allows for an effective or mean combustion timing closer to TDC. For
the TEVO constraint, the late injection provides a sufficient TEVO increase for the
first injection to be timed optimally.
185
Chapter 10. Constraint Handling with Multiple Injections
150
p [bar]
100
50
0
−20 −10 0 10 20 30 40 50
θ [CAD]
Small fuel quantities are not always sufficient to generate heat-release. For in-
stance, in Chapter 8, a small pilot was only used to enhance the reactivity of the
main-injection and did not generate a separate heat-release impulse. The detec-
x
tion method therefore only allocated detected θCT to injections with sufficient
fuel mass. Values in the range of 15-30 mg were used as lower limits in the engine
experiments presented below. When a single combustion timing was expected,
θ50 was used instead of θCT as a combustion-timing indicator. A similar proce-
dure was applied for detecting the peak pressure levels generated by the pilot
1 2
and main injection p max , p max which are necessary quantities for the control-
c
ler to fulfill p max . The motored pressure curve was first subtracted from p before
1,2
detecting p max , to not detect the pressure peak generated by compression.
186
10.4 Controller Design
1,500
NOx [ppm]
1,000
500
−20 −10 0 10 20 30 40 50
θ [CAD]
where k is cycle index, e p IMEPn and e θ2 are p IMEPn and θCT errors
CT
r
e θ2 (k) = θCT (k) − θθ2 (k)
CT CT
(10.19)
r
e p IMEPn (k) = p IMEPn (k) − p IMEPn (k)
Pressure Constraint
2 c
The main-injection timing θSOI was adjusted if p max was violated to keep the
2 c
pressure peak corresponding to the main injection, p max , below p max . This be-
havior was obtained by the following controller
( 2 2 c
θ (k) + ∆PICT (e θ2 (k)) if p max (k) < p max
θSOI (k + 1) = SOI
2
2
CT (10.20)
θSOI (k) + ∆PIp max (e p max
2 (k)) otherwise
187
Chapter 10. Constraint Handling with Multiple Injections
1,500
T [K]
1,000
500
where
2 c
e p max
2 (k) = p max (k) − p max (10.21)
and ∆PIx are combustion timing updates from two separate PI controllers on the
c
form of (10.17). A pilot injection was also introduced if p max > p max , in order to
obtain a distributed heat-release rate, see Fig. 10.7. This was done by varying r 1
according to
(
1 r 1 (k) + k r (r ∗1 − r 1 (k)) c
if p max (k) > p max
r (k + 1) = 1 1
(10.22)
r (k) − k r r (k) otherwise
where r ∗1 is a predefined set point. Furthermore, once the pilot injection was in-
1 c
troduced, a pilot-injection controller was used to keep p max below p max
( 1 1 c
1
θSOI (k) + ∆PICT (e θ1 (k)) if p max (k) < p max
θSOI (k + 1) = 1
CT (10.23)
θSOI (k) + ∆PIp max (e p max
1 (k)) otherwise
where
r 1
e θ1 (k) = θCT (k) − θCT (k)
CT
(10.24)
1 c
e p max
1 = p max − p max
188
10.4 Controller Design
NOx Constraint
NOx -constraint fulfillment was handled similarly
(
2
2
θSOI (k) + ∆PICT (e θ2 (k)) if NOx (k) < NOcx
θSOI (k + 1) = CT (10.25)
2
θSOI (k) + ∆PINOx (e NOx (k)) otherwise
where
e NOx (k) = NOx (k) − NOcx (10.26)
A pilot injection was also introduced if NOx (k) > NOcx to obtain a distributed
heat-release rate according to Fig. 10.8
(
1 r 1 (k) + k r (r ∗1 − r 1 (k)) if NOx (k) > NOcx
r (k + 1) = 1 (10.27)
r (k) − k r r 1 (k) otherwise
1
was used to adjust θSOI for constraint fulfillment.
Temperature Constraint
The exhaust temperature Tex was controlled by adjusting the post injection to
obtain the configurations presented in Fig 10.9
(
3 r 3 (k) + k r (r ∗3 − r 3 (k)) c
if Tex (k) < Tex
r (k + 1) = 3
r (k) − k r r 3 (k) otherwise
( (10.29)
3 c
3
θSOI (k) + ∆PICT (e θ3 (k)) if Tex (k) > Tex
θSOI (k + 1) = CT
3
θSOI (k) + ∆PITex (e Tex (k)) otherwise
where
r 3
e θ3 (k) = θCT (k) − θCT (k)
CT
(10.30)
c
e Tex = Tex − Tex
p IMEPn Control
The engine load was controlled by adjusting the total injected fuel mass using the
controller
m tot tot
f (k + 1) = m f (k) + ∆PIp IMEPn (e p IMEPn (k)) (10.31)
189
Chapter 10. Constraint Handling with Multiple Injections
Saturation
Saturation limits were introduced to enforce
u ∈U (10.32)
1 3
This was done by saturating the pilot and post injection timings θSOI , θSOI to
avoid overlap with the main injection. Unfortunately, this could lead to dead-
1 1,2
lock and constraint violation for p max (k) when regulating p max . A remedy to this
1 1 2
problem was to reduce r when p max (k) > p max (k)
1,2,3
Additional minimum and maximum limits on θSOI , m tot
f
and r 1,3 were introdu-
ced as safety margins to ensure feasible injection timings and durations.
Slack Variables
Slack variables ²x were added to the constraint conditions above
c
p max (k) > p max − ²p max
NOx (k) > NOcx − ²NOx (10.34)
c
Tex (k) < Tex + ²Tex
Summary
The controllers presented above are summarized in Algorithm 4. Here, p max - and
NOx -constraint handling were merged using the max function. In this way, the
controller would adjust for the constraint demanding the latest θSOI , since both
p max and NOx decrease with θSOI . The notation
x̃ c = x c ± ²x (10.35)
190
10.5 Experimental Evaluation
Pilot-Injection Controller
c ˜ x then c
1: if p max (k) > p̃ max ∨ NOx (k) > NO
1 1
¡ ¢
2: θSOI (k + 1) = θSOI (k) + max ∆PIp max (e p max
1 (k)), ∆PINOx (e NOx (k))
c 1 2
3: if p max (k) > p̃ max ∧ p max (k) > p max (k) then
4: r 1 (k + 1) = r 1 (k) − k r r 1 (k)
5: else
6: r 1 (k + 1) = r 1 (k) + k r (r ∗1 − r 1 (k))
7: end if
8: else
1 1
9: θSOI (k + 1) = θSOI (k) + ∆PICT (e θ1 (k))
CT
1 1 1
10: r (k + 1) = r (k) − k r r (k)
11: end if
Main-Injection Controller
c ˜ x then c
12: if p max (k) > p̃ max ∨ NOx (k) > NO
2 2
¡ ¢
13: θSOI (k + 1) = θSOI (k) + max ∆PIp max (e p max
2 (k)), ∆PINOx (e NOx (k))
14: else
2 2
15: θSOI (k + 1) = θSOI (k) + ∆PICT (e θ2 (k))
CT
16: end if
Post-Injection Controller
17: if Tex (k) < T̃ex then
3 3
18: θSOI (k + 1) = θSOI (k) + ∆PITex (e Tex (k))
19: r 3 (k + 1) = r 3 (k) + k r (r ∗3 − r 3 (k))
20: else
3 3
21: θSOI (k + 1) = θSOI (k) + ∆PICT (e θ3 (k))
CT
191
Chapter 10. Constraint Handling with Multiple Injections
6: end if
Transient Operation
c
p max c
The ability to handle p max constraints was evaluated with p max = 75. Ex-
r
perimental results with p IMEPn step changes are presented in Fig. 10.10, where
r
p IMEPn increased from 4 to 10 bar in 10 cycles as p IMEPn was increased.
2 r c c
At p IMEPn = 10 bar with θCT = θCT , p max violated p max . To fulfill p max , the con-
1 2
troller increased r and delayed θCT , and reached a new injection configuration
r
in 50 cycles. The controller converged to the initial conditions once p IMEPn was
decreased.
Cycle-resolved data at cycles 1250 and 1450 are presented in Fig. 10.11. At
2 r
cycle 1250, a single injection was used where p IMEPn = 4 bar, and θCT = θCT =8
1
CAD. At cycle 1450, with r = 0.5, the injection timings were adjusted to fulfill
c 1 2 r
p max . The red and blue vertical lines indicate detected θCT and θCT , and θCT is
indicated by the horizontal dashed line.
NOx c Fulfillment of NOcx = 800 ppm is presented in Fig. 10.12. In this figure,
the controller was compared to the single-injection controller in Algorithm 5,
which is indicated in purple. The combustion timings were delayed when NOx
r
was increased due to positive p IMEPn step changes. The double-injection control-
ler introduced a pilot injection when NOcx was violated, and an additional com-
θ θ
k pCT [-] 0.1 k I CT [-] 0.15
NO NO
kp x [CAD/ppm] 10 k I x [CAD/ppm] 15
p p
k p max [CAD/bar] 0.04 k I max [CAD/bar] 0.05
T T
k p ex [CAD/K] 0.3 k I ex [CAD/K] 0.35
mf mf
kp [mg/bar] 1.5 kI [mg/bar] 2
192
10.5 Experimental Evaluation
10
p max [bar]
80
5
60
0
1,200 1,300 1,400 1,500 1,600 1,200 1,300 1,400 1,500 1,600
cycle [-] cycle [-]
0.8
θ SOI / θ CT [CAD]
10
0.6
r 1 [-]
0.4
0
0.2
−10 0
1,200 1,300 1,400 1,500 1,600 1,200 1,300 1,400 1,500 1,600
cycle [-] cycle [-]
bustion timing (red) was detected. The fuel-masses presented are the fuel masses
demanded by the controllers, and not the actual injected fuel mass. Both con-
trollers had comparable p IMEPn response times of 10 cycles and injection-timing
settling times of 25 cycles. The double-injection controller had a slightly larger
NOx overshoot.
c
Tex c
Fulfillment of Tex = 240◦ C is presented in Fig. 10.13. The combustion ti-
ming was delayed for the single-injection controller in order to fulfill the con-
r
straint when Tex was decreased due to negative p IMEPn step changes. A post in-
193
Chapter 10. Constraint Handling with Multiple Injections
p [bar]
50 50
0 0
0 20 0 20
θ [CAD] θ [CAD]
Figure 10.11 Cycle-resolved data belonging to cycle 1250 and 1450 in Fig. 10.10.
1,2 r . The dashed horizontal lines indicate
The vertical lines indicate θCT and θCT
c
p max .
194
10.5 Experimental Evaluation
1,500
Two Injections
NOx [ppm]
1,000
500
One Injection
0
1,000 1,100 1,200 1,300 1,400 1,500 1,600 1,700 1,800 1,900 2,000
cycle [-]
p IMEPn S r
p IMEPn r,S
10 p IMEPn p IMEPn
p IMEPn [bar]
4
1,000 1,100 1,200 1,300 1,400 1,500 1,600 1,700 1,800 1,900 2,000
cycle [-]
m 1f m 2f m Sf
100
m f [mg]
50
0
1,000 1,100 1,200 1,300 1,400 1,500 1,600 1,700 1,800 1,900 2,000
cycle [-]
20
θ SOI / θ CT [CAD]
−20
1,000 1,100 1,200 1,300 1,400 1,500 1,600 1,700 1,800 1,900 2,000
cycle [-]
195
Chapter 10. Constraint Handling with Multiple Injections
320
Tex S
Tex c Tex
300 One Injection Two Injections
Tex [◦ C] 280
260
240
220
200 400 600 800 1,000 1,200 1,400 1,600 1,800 2,000
cycle [-]
12 S r r,S
p IMEPn p IMEPn p IMEPn p IMEPn
10
p IMEPn [bar]
8
6
4
200 400 600 800 1,000 1,200 1,400 1,600 1,800 2,000
cycle [-]
m 2f m 3f m Sf
100
m f [mg]
50
0
200 400 600 800 1,000 1,200 1,400 1,600 1,800 2,000
cycle [-]
20
θ SOI / θ CT [CAD]
−20
200 400 600 800 1,000 1,200 1,400 1,600 1,800 2,000
cycle [-]
196
10.6 Discussion
Pressure Constraint
120
ηGIE = 0.442 [-] ηGIE = 0.468 [-]
100 NOx = 2111 [ppm] NOx = 1860 [ppm]
60
40
20
c
Figure 10.14 Fulfillment of p max = 95 bar with one (blue) and two (red)
injections. Two injections allowed for a combustion timing closer to TDC,
which increased gross-indicated efficiency with 5 %. A disadvantage with the
double-injection strategy was the significant increase in soot emissions.
10.6 Discussion
From the experimental results in Sec. 10.5, it can be concluded that multiple
injections can increase the indicated efficiency when stringent constraints on
p max and NOx are imposed. The efficiency increase was a result of the distribu-
ted heat-release rate with reduced peak in-cylinder pressure and temperature,
which allowed for a more advanced effective combustion timing, see Fig. 10.6.
These results agree with previous simulation and experimental work on optimal
heat-release rates [Eriksson and Sivertsson, 2016; Okamoto and Uchida, 2016;
Guardiola et al., 2017].
The problem of optimally calibrating multiple injection timings and dura-
tions for different engine operating points is both demanding and sensitive
to disturbances. The contribution of the work presented in this chapter is a
197
Chapter 10. Constraint Handling with Multiple Injections
0.49 5
0.48 4.5
Soot [mg/m3 ]
η NIE [-]
0.47 4
0.46 3.5
two injections
one injection
0.45 3
400 600 800 1,000 1,200 400 600 800 1,000 1,200
NOx [ppm] NOx [ppm]
45 340
40 320
HC [ppm]
35 Tex [◦ C] 300
30 280
25 260
400 600 800 1,000 1,200 400 600 800 1,000 1,200
NOx [ppm] NOx [ppm]
Figure 10.15 Steady-state NOx trade-offs with the single and double-injection
controller. The double-injection controller improved the NOx trade-off with η NIE
by 4 %. It also improved the trade-off with HC. The double-injection controller
had a worsened trade off with respect to soot, and a partly worsened trade-off
with Tex .
198
10.6 Discussion
0.48
Soot [mg/m3 ]
0.47
η NIE [-]
0.46
two injections
2
one injection
0.45
210 220 230 240 250 210 220 230 240 250
Tex [◦ C] Tex [◦ C]
200 800
150 700
NOx [ppm]
HC [ppm]
100 600
50 500
0 400
210 220 230 240 250 210 220 230 240 250
Tex [◦ C] Tex [◦ C]
Figure 10.16 Steady-state Tex trade-offs with the single and double-injection
controller. The results show that the double-injection controller worsened the Tex
trade-off with η NIE and NOx for high Tex c , but improved the trade-off with HC
emissions.
199
Chapter 10. Constraint Handling with Multiple Injections
p [bar]
50 50
0 0
0 20 40 0 20 40
θ [CAD] θ [CAD]
p [bar]
50 50
0 0
0 20 40 0 20 40
θ [CAD] θ [CAD]
Figure 10.17 Pressure and heat-release rates for the constraint sweeps presented
in Figs. 10.15 and 10.16.
fillment were 100-200 cycles. These settling times were comparable to those of
the single-injection controller, that adjusted one injection timing to fulfill con-
straints.
Steady-state experiments showed a 4-5 % efficiency increase with respect to
p max and NOx constraints, and a 1.7 % efficiency decrease with respect to Tex
c
constraints. The decrease in efficiency with Tex contradicted the simulation re-
sults in Sec. 10.3. One possible explanation for the efficiency decrease could be
the assumption of a constant heat-release shape with combustion timing. In re-
ality, heat-release rates and combustion efficiency decrease with combustion ti-
ming, which penalize the efficiency of late post injections. It was however not
easy to verify this from the heat-release rates in Fig. 10.17. Despite this fact, the
controller suggested provides a framework for how to introduce post injections
when favorable. Results presented in [Honardar et al., 2011] showed that a post
200
10.7 Conclusions
10.7 Conclusions
A controller was designed and implemented for increased thermodynamic effi-
ciency when constraints on p max , NOx and Tex are imposed. The controller de-
sign was motivated by previous research on multiple injections and a presented
0D-simulation study that showed a 2-4 % efficiency increase when introducing
an additional injection. These results suggested that the controller should ad-
just the number of injections for efficient constraint fulfillment. This was done
by introducing a pilot injection when encountering active p max and NOx con-
straints, whilst exhaust-temperature constraints were handled with a post injec-
tion. A single injection was found optimal when no constraints were active.
The desired controller behavior was obtained with a hybrid, multiple-input
multiple-output PI controller that utilized feedback from in-cylinder pressure
measurements, a NOx -emission model that functioned as a virtual sensor, and
measured exhaust temperature. The suggested controller was experimentally
evaluated, where it showed comparable transient performance to that of a
single-injection controller. The controller exhibited improved p max / η NIE and
NOx / η NIE trade-offs with a 4-5 % increase in η NIE . The controller also showed a
worsened Tex / η NIE trade-off with a 1.7 % decrease in η NIE . Increased soot emis-
sions levels with two injections were also observed.
201
11
Pressure Prediction and
Efficiency Optimization
11.1 Introduction
Combustion timing has traditionally been controlled in open loop by means of
experimentally calibrated injection-timing maps [Guzella and Onder, 2009]. This
approach requires a considerable calibration effort and can be sensitive to va-
riations in hardware and fuel properties, especially in low-temperature combus-
tion modes. Experimental results in the previous chapters of this thesis have
shown that closed-loop combustion-timing control can be used to accurately
track combustion-timing set points and make the combustion timing robust to
disturbances.
Closed-loop combustion-timing controllers can be divided into two sub-
groups, one where the controller tracks a predefined combustion-timing set
point, where set-point optimization is considered to be a separate task. Exam-
ples of such controllers were discussed in chapters 6, 7 and 9, where PI control
and MPC were used for set-point tracking. Then, there are controllers that in-
stead adjust the combustion timing to directly fulfill higher-level performance
targets, such as emission-limit fulfillment and efficiency maximization. An ex-
ample of a high-level performance controller was presented in [Karlsson et al.,
2010], where a dynamical black-box model related the injection timing to p IMEPg ,
d p max and NOx emissions to minimize fuel consumption subject to specified
output constraints. Similar data-driven approaches were presented in [Hafner
et al., 2000; Atkinson et al., 2009]. Extremum-seeking control is another example
of a data-driven controller design that aims to fulfill higher-level specifications.
Extremum-seeking control has previously been used to find efficiency-optimal
combustion timings through set-point perturbation, see [Lewander et al. 2012].
This technique was also investigated by Killingsworth et al. [2009] and Hellström
et al. [2013], in HCCI- and spark-ignition engines, respectively.
202
11.2 Controller Description
Pressure Prediction
First, the controller obtains the previous-cycle cylinder pressure p 0 and compu-
tes the corresponding heat-release rate dQ c0 /d θ using (4.1).
After computing dQ c0 /d θ, the objective is to predict how the pressure chan-
ges with combustion timing. A pressure change due to a change in combustion
203
Chapter 11. Pressure Prediction and Efficiency Optimization
dQ c+ dQ c0 (θ + ∆θ)
=
dθ dθ
(11.1)
dQ c− dQ c0 (θ − ∆θ)
=
dθ dθ
h c ATIVC
µ(p, θ) = (γ − 1) p (11.4)
p IVC VIVC
(see, (2.15)), where α was introduced as a heat-transfer tuning parameter. The re-
ason for linearizing (11.2) is that ∆p can be computed from the solution to (11.3)
Z θ d ∆Q c (ϑ)
∆p(θ) = Φ(θ, ϑ)Γ(ϑ) dϑ (11.6)
θIVC dϑ
where µ Z ¶µ ¶γ
θ ∂µ(p 0 , τ) V (ϑ)
Φ(θ, ϑ) = exp − dτ
ϑ ∂p V (θ)
(11.7)
γ−1
Γ(ϑ) =
V (ϑ)
204
11.2 Controller Description
d ∆Q c+ dQ c+ dQ c0
= −
dθ dθ dθ
(11.8)
d ∆Q c− dQ c− dQ c0
= −
dθ dθ dθ
Injection-Timing Optimization
With p + and p − , variations in quantities such as p IMEPg , p max and d p max can be
computed
Z VEVO Z VEVO
+ 1 1
p IMEPg = p + dV −
p IMEPg = p − dV
Vd VIVC Vd VIVC
+
p max = maxp + −
p max = maxp − (11.9)
θ θ
+
d p max = maxd p + −
d p max = maxd p −
θ θ
A simple model for p IMEPg , p max and d p max in the subsequent engine cycle can
then be formulated with the partial derivatives in (11.10)
∂p IMEPg
p IMEPg p0 ∂∆θ
IMEPg
p max = p 0 + ∂p max ∆θ (11.11)
max
0
∂∆θ
d p max d p max ∂(d p )
max
∂∆θ
205
Chapter 11. Pressure Prediction and Efficiency Optimization
∗ 1 d p IMEPg
∆θuc = (11.13)
2β d ∆θ
c c
Combustion-timing shifts ∆θ that reach the constraint limits p max and d p max are
given by
p c − p max
0
∆θp max
c = max
∂p max /∂∆θ
(11.14)
c 0
d p max − d p max
∆θd p max
c =
∂(d p max )/∂∆θ
Now, if one assumes that p max and d p max are monotonically decreasing with ∆θ,
the solution to (11.12), ∆θ ∗ , is simply given by the largest value among ∆θuc∗
,
c
∆θp max and ∆θd p max , saturated within the limits of ∆θ . Furthermore, the opti-
c c
11.3 Results
This section presents experimental controller results. Open-loop experiments
are first presented for the purpose of evaluating how well the controller predicts
206
11.3 Results
Open-loop Experiments
The injection timing θSOI was swept with a single injection for three different
speed/load combinations to evaluate the pressure-prediction method. Condi-
tions at these three operating points are presented in Table 11.1. Each sweep
consisted of 2000 cycles where θSOI was incremented in steps of one from ap-
proximately -25 to 5 CAD after TDC.
Table 11.1 The data used in the prediction- and controller evaluation were ob-
tained from the following operating points.
Operating Point 1 2 3
207
Chapter 11. Pressure Prediction and Efficiency Optimization
Operating Point 1
0.52
0.5
ηGIE [-]
0.48
0.46
0.44
0.42
−10 −5 0 5 10 15 20
θ 50 [CAD]
Operating Point 2
0.52
0.5
ηGIE [-]
0.48
0.46
0.44
0.42
−10 −5 0 5 10 15 20
θ 50 [CAD]
Operating Point 3
0.52
0.5
ηGIE [-]
0.48
0.46
0.44
0.42
−10 −5 0 5 10 15 20
θ 50 [CAD]
208
11.3 Results
20
OP1 OP2 OP3
15
10
θ 50 [CAD]
−5
−20 −15 −10 −5 0
θ SOI [CAD]
Figure 11.2 The relation between θ50 and θSOI for the θSOI sweeps in Fig. 11.1
as θ50 was increased or decreased outside this interval. The efficiency was sli-
ghtly higher at operating point 1. This is believed to be a result of the difference
in intake temperature Tin , which was due to engine warm up and the order of
which the experiments were conducted.
∂θ50 /∂θSOI Figure 11.2 shows θ50 as a function of θSOI for the three θSOI sweeps
in Fig. 11.1. The assumption of a constant partial derivative ∂θ50 /∂θSOI was ac-
curate close to θ50 = 5. However, at the lower-load operating points, ∂θ50 /∂θSOI
decreased at early θSOI and increased at late θSOI . This can be explained by the in-
crease in ignition delay τ when θSOI was decreased or increased. This effect was
then stronger at low load where τ was longer. Similar trends for ∂θ50 /∂θSOI with
engine load was observed in Chapter 6 and Fig. 6.3.
209
Chapter 11. Pressure Prediction and Efficiency Optimization
140
p0
p−
120
p+
p̂ −
100 p̂ +
p [bar]
80
60
40
−20 −10 0 10 20 30
θ [CAD]
Figure 11.3 The black solid pressure curve is the cycle-averaged pressure p 0 for
0 . The red and blue pressure curves to the left and right of this curve are the
θSOI
0 . The
cycle-averaged pressures p + and p − for θSOI shifted ±1 CAD relative to θSOI
dashed red and blue curves are the predicted pressures p̂ − and p̂ + , computed
with (11.6) and the same θSOI shifts.
restingly, it can also be seen that the pressure deviation resulting from a θSOI shift
closely resembled a heat-release rate.
Figures 11.5-11.7 presents cycle-averaged pressure changes p ± − p 0 toge-
ther with the prediction errors p ± − p̂ ± (as in Fig. 11.4) for all θSOI at operating
points 1-3. Just as in Fig. 11.4, the blue (red) color indicate a pressure change
p + − p 0 (p − − p 0 ) due to a delayed (advanced) θSOI . It can be seen that the pre-
diction error was smaller for θSOI close to TDC for all operating points and that
the error gradually increased when θSOI was advanced or delayed. The prediction
error also changed sign somewhere around TDC. For the red lines, indicating a
delayed θSOI , this meant that p̂ − < p − for late θSOI , and p̂ − > p − for early θSOI .
The opposite trend was found for p̂ + and p + .
The pressure-prediction performance was also evaluated by computing the
coefficient of determination R 2
210
11.3 Results
−5
−10
−15
−20 −10 0 10 20 30
θ [CAD]
Figure 11.4 The pressure differences p + − p 0 (blue, solid) and p − − p 0 (red, so-
lid), together with the prediction errors p + − p̂ + (blue, dashed) and p − − p̂ − (red,
dashed) for the pressure curves in Fig. 11.3. Interestingly, the pressure deviation
closely resembled a heat-release rate.
Operating Point 1
20
10
p [bar]
−10
−20
−10 −5 0 5 10 15 20 25
θ [CAD]
211
Chapter 11. Pressure Prediction and Efficiency Optimization
Operating Point 2
15
10
5
p [bar]
−5
−10
−15
−10 −5 0 5 10 15
θ [CAD]
Operating Point 3
20
10
p [bar]
−10
−20
−4 −2 0 2 4 6 8 10 12
θ [CAD]
212
11.3 Results
Prediction Performance
1
0.9
0.8
R 2 [-]
0.7
0.6
0.5
OP 1 OP 2 OP 3
0.4
−4 −2 0 2 4 6 8 10 12 14 16
θ 50 [CAD]
Figure 11.8 R 2 score as a function of θ50 for the three sweeps in Figs. 11.5-11.7.
The pressure-prediction method worked satisfactorily for θ50 ∈ [0, 6], where R 2 ≥
0.9. The performance then started to degrade outside this interval, and more ste-
eply at the low-load operating points.
load experiment in Fig. 11.8, where τ did not change significantly with an almost
constant ∂∆θ/∂θSOI , see Fig. 11.2.
Closed-loop Experiments
This section demonstrates closed-loop performance. Tuning for best perfor-
mance was not carried out as the experiments were focused on convergence and
parameter sensitivity.
Convergence and β - Sensitivity The controller in Algorithm 6 was evaluated
at the investigated operating points. The parameter β and the initial injection ti-
0
ming θSOI were varied to investigate controller convergence. Convergence results
0
are presented in Figs. 11.9-11.11, where θSOI = {20, 10, 0} [CAD] and β = {0.05, 0.2}.
∗
The controller consistently converged to the same θ50 , independently of the
starting point in Figs. 11.9-11.11. The parameter β clearly influenced the conver-
gence rate, where a larger β gave slower convergence. In Fig. 11.10, the controller
behavior in stationarity became oscillatory for β = 0.05. In the same figure, it also
0
seems as if the controller converged faster for β = 0.2 and θ50 = −2. This was cau-
sed by an unintended active p max constraint. The constraint did not affect the
stationary behavior since the estimated d p IMEPg /d ∆θ was zero in the conver-
gence point.
Figures 11.9-11.11 show that 1/β can be viewed as a controller gain, where
the choice of β is a trade-off between convergence speed and stationary
cycle-to-cycle variation. If β was chosen too small, the derivative d p IMEPg /d ∆θ
213
Chapter 11. Pressure Prediction and Efficiency Optimization
Controller Convergence
β = 0.2
β = 0.05
10
θ 50 [CAD]
Figure 11.9 Controller convergence at operating point 1 where the solid lines
correspond to β = 0.05 and the dashed lines to β = 0.2. The controller converged
∗ = 5.8, θ ∗ = −7, in 50 and 150 cycles, depending on β.
to θ50 SOI
Controller Convergence
15 β = 0.2
β = 0.05
10
θ 50 [CAD]
Figure 11.10 Controller convergence at operating point 2, where the solid lines
correspond to β = 0.05 and the dashed lines to β = 0.2. Here, the point of conver-
∗ = 6.4, θ ∗ = −9.5. In the lower dashed θ trajectory, the fast con-
gence is θ50 SOI 50
vergence was due to an unintended active p max constraint. This did however not
affect the point of convergence, since the computed ∂p IMEPg /∂∆θ was zero in the
convergence point.
214
11.3 Results
Controller Convergence
15 β = 0.2
β = 0.05
10
θ 50 [CAD]
Figure 11.11 Controller convergence at operating point 3, where the solid lines
correspond to β = 0.05 and the dashed lines to β = 0.2. Here, the point of con-
∗ = 5.6, θ ∗ = −9.6, where the controller converged in 5 cycles for
vergence is θ50 SOI
β = 0.05.
caused large controller steps ∆θSOI which led to an oscillatory controller be-
havior. For best performance, β should be increased with load as d p IMEPg /d ∆θ
increases. The point of convergence in Figs. 11.9-11.11 occurred later than the
experimentally found most efficient points in Fig. 11.1. The reason for this will
be discussed in the following section.
∗
Parameter Sensitivity The point of convergence θ50 depends on the parame-
ter values in (11.2), and especially on the parameters of the heat-transfer model.
∗
In order to investigate θ50 sensitivity, the model parameter α (see (11.5)) and the
TDC offset ∆θTDC were varied. These parameters were previously set to α = 5 and
∆θTDC = 0. Convergence results can be viewed in Figs. 11.12 and 11.13 where α
and ∆θTDC were varied stepwise from 1 to 6, and from 2 to -2 CAD, respectively.
Figure 11.12 shows that the magnitude of α affected the convergence point. With
∗
an increased α, θ50 was delayed and the converse was true for a decreased α.
One could view α as a trade-off parameter that weighs efficiency effects from
heat-transfer and exhaust losses, and in that way determines an efficiency opti-
∗
mal θ50 .
∗
In Fig. 11.13, it can be seen that a ∆θTDC of 2 CAD gave a θ50 offset of ap-
proximately 2 CAD. This indicates that the controller accuracy is limited by the
precision of the measured θ, and its synchronization with p.
215
Chapter 11. Pressure Prediction and Efficiency Optimization
α=6
5
α=5
θ 50 [CAD]
α=4
α=3
0
α=2
α=1
100 200 300 400 500
cycle [-]
Figure 11.12 Heat-transfer sensitivity at operating point 1. Here the scaling fac-
tor α was varied stepwise from 1 to 6, which changed the point of convergence.
An α between 3 and 4 would have maximized efficiency according to the data in
Fig. 11.1. In Figs. 11.9-11.11, α = 5 was used.
10 ∆θTDC = 2
θ 50 [CAD]
8 ∆θTDC = 0
6 ∆θTDC = -2
Figure 11.13 TDC-offset sensitivity at operating point 1. Here, ∆θTDC was chan-
ged stepwise from 2 to -2 CAD, which indicates that the controller accuracy is li-
mited by the precision of the measured θ and its synchronization with p.
Pressure-Constraint Handling
Controller performance with respect to constraint fulfillment was also investiga-
c c
ted by varying the constraint limits p max and d p max . Result are presented in Figs.
c c
11.14 and 11.15, where p max was varied between 80 to 120 bar and d p max was va-
ried between 20 to 40 bar/CAD. The controller managed to fulfill the constraints
by initially taking a larger positive step in θSOI and then slowly advancing θ50 to
reach the constraint level from below.
216
11.4 Discussion
120
p max [bar]
100
80
15
θ 50 [CAD]
10
5
0 100 200 300 400 500 600 700 800 900 1,000
cycle [-]
Figure 11.14 In the upper part of the figure, p max (solid) is displayed together
c
with p max (dashed), The corresponding θ50 is displayed in the lower subdiagram.
11.4 Discussion
A desirable combustion timing is a trade-off between exhaust losses, heat trans-
fer, and constraint fulfillment. This chapter introduced a high-level model-based
combustion-timing controller that finds the efficiency-optimal combustion
217
Chapter 11. Pressure Prediction and Efficiency Optimization
50
d p max [bar/CAD]
40
30
20
10
0 100 200 300 400 500 600 700
cycle [-]
20
15
θ 50 [CAD]
10
Figure 11.15 In the upper part of the figure, d p max (solid) is displayed together
c
with d p max (dashed). The corresponding θ50 is displayed in the lower subdia-
gram.
timing online. In this way, the controller solves the problem of deciding a
combustion-timing set point, which is necessary for most combustion-feedback
controllers, see for instance Chapter 6 and [Bengtsson et al., 2004; Chiang and
Stefanopoulou, 2005; Widd et al., 2008]. Furthermore, the model-based approach
presented here allows for faster convergence (5 cycles) than data-based extre-
mum seeking controllers that search for efficiency-optimal θ50 using the com-
puted p IMEPg , see [Killingsworth et al. 2009; Lewander et al. 2012]. The method
presented in this chapter could speed up convergence of such methods by provi-
ding an initial guess.
Feedback was introduced through the estimated heat-release rate, which was
utilized by the controller, together with a linearized 0D model, to predict how
the cylinder pressure varies with injection timing. This is a computationally ef-
ficient alternative to that of modeling the relation between fuel injection and
heat-release rate.
The point of convergence was found to be sensitive to TDC offset and
218
11.5 Conclusions
11.5 Conclusions
A model-based combustion-timing controller was introduced. The controller
utilized the estimated heat-release rate to predict how the in-cylinder pressure
varies with injection timing. The controller converged close to experimentally
found most efficient combustion timing, and was capable of fulfilling constraints
with respect to p max and d p max .
219
12
Predictive Constraint
Handling and Pressure
Tracking
220
12.1 Model Predictive Control Formulation
where m f denotes the fuel mass, θSOI denotes injection timing, and M is the
number of injections. Injector-current pulse durations θDOI are computed from
m f and the common-rail pressure p rail using an injector map
i
θDOI = M inj (m if , p rail ) (12.2)
Constraint Fulfillment
The first MPC problem concerns optimization problem (10.10), which was stu-
died in Chapter 10. Here, (10.10) is reformulated as an MPC problem with a pre-
diction horizon and additional costs on control action and load-tracking error.
221
Chapter 12. Predictive Constraint Handling and Pressure Tracking
222
12.1 Model Predictive Control Formulation
where p IMEP is used to denote p IMEPg in this chapter. The control-action cost
term in (12.3) is given by
Pressure Tracking
The MPC for pressure tracking is given by
Hp
X
minimize J p (k) + J ∆u (k) (12.7)
u(1), . . . , u(H p ) k=1
subject to u(k) ∈ U, k = 1, . . . , H p
θf
X ¡ ¢2
J p (k) = αt r p r (θ j ) − p(θ j , k) (12.8)
θ j =θ0
Input Constraints
Input constraints were enforced through saturation in Chapter 10. Here, the in-
put constraints are instead included in the optimization problem, where con-
straints on u are imposed in order for the solution to be realizable. Injected fuel
masses have to be positive
m if ≥ 0, i = 1, . . . , M (12.9)
223
Chapter 12. Predictive Constraint Handling and Pressure Tracking
140
αv = 1
αv = 0.5
120
αv = 0
100
80
p r [bar]
60
40
20
0
−20 −10 0 10 20 30 40 50
θ [CAD]
Cost-Function Weights
The cost-function weights α, αt r , β and the diagonal matrix
³ ´
1 1 M M
R = diag R m f
R SOI . . . Rm f
R SOI >0 (12.12)
224
12.1 Model Predictive Control Formulation
are parameters that determine the trade-off between control action and fulfill-
ment of the different tracking objectives. Suitable values will be presented toge-
ther with experimental results in Sec. 12.5.
Problem Approximation
We have now obtained two nonlinear-programming problems on the form
minimize J (U ) (12.13)
U
subject to f (U ) ≤ c
U ∈U
where
¡ ¢T
U = u(1) ... u(H p ) (12.14)
1
minimize J (u 0 ) + ∇J (u 0 )T ∆U + ∆U T ∇2 J (u 0 )∆U (12.15)
∆U 2
subject to f (u 0 ) + ∇ f (u 0 )T ∆U ≤ c
u 0 + ∆U ∈ Û
Another motivation for approximating (12.13) was that gradients and Hessians
can be computed using simpler physical models, valid for smaller variations in u.
The optimization problem can then be re-approximated on a cycle-to-cycle ba-
sis when new measurements have been obtained. This means that the gradients
∇ f (u 0 ), ∇J (u 0 ) and Hessian ∇2 J (u 0 ) have to be recomputed as u 0 changes. The
following section presents how to obtain gradients and Hessians with the use of
physical models. A linear approximation of Û was obtained by linearizing M inj
i −1 i −1 i −1
∂M inj (u 0 )
θSOI (0) + ∆θSOI + θDOI (0) + ∆m if −1 + θh ≤ θSOI
i i
(0) + ∆θSOI (12.17)
∂m f
225
Chapter 12. Predictive Constraint Handling and Pressure Tracking
From ∆u to d ∆Q c /d θ
First, the previous cycle dQ c /d θ is obtained from
dQ c γ dV 1 d p dQ ht
= p + V + (12.18)
dθ γ − 1 dθ γ − 1 dθ dθ
dQ c X M dQ i
c
= (12.19)
dθ i =1 d θ
where the bounds l i and d i are determined by the minima between the detec-
ted peaks. The obtained heat-release rates d Q̂ ci /d θ are then smoothed with a
zero-phase filter to obtain more physical heat-release shapes. Finally, d Q̂ ci /d θ
have to be normalized so that (12.19) is fulfilled
dQ ci ³XM d Q̂ i ´−1 d Q̂ i dQ
c c c
= (12.21)
dθ i =1 d θ d θ d θ
The four steps of the detection procedure are illustrated in Fig. 12.2. In the case
of detecting less than M peaks, which is possible for small injections or if com-
bustion from different injections overlap, the detected peaks are allocated to the
largest injections. Furthermore, the detection procedure accounts for the orde-
ring of injections when allocating combustion timings to injections.
226
12.2 Modeling and Heat-Release Detection
600 600
dQ c /d θ [J/CAD] 1. 2.
dQ c /d θ [J/CAD]
400 400
200 200
0 0
−20 0 20 −20 0 20
θ [CAD] θ [CAD]
600 600
3. 4.
dQ c /d θ [J/CAD]
dQ c /d θ [J/CAD]
400 400
200 200
0 0
−20 0 20 −20 0 20
θ [CAD] θ [CAD]
Figure 12.2 The method used for separating dQ c /d θ among the different in-
jections (from 1 to 4). First, the heat-release rate is obtained from the measured
pressure signal (1) and the most significant M peaks are detected (2). dQ c /d θ
is then separated in different intervals according to the peak locations (3). The
heat-release rates d Q̂ ci /d θ in the different intervals are then filtered and normali-
zed (4).
∂ dQ c d 2Q ci
=−
i
∂θSOI dθ d θ2
(12.22)
i
∂θCT
i
=1
∂θSOI
227
Chapter 12. Predictive Constraint Handling and Pressure Tracking
600 600
dQ c /d θ [J/CAD]
dQ c /d θ [J/CAD]
400 400
200 200
0 0
−20 0 20 40 −20 0 20 40
θ [CAD] θ [CAD]
Figure 12.3 The assumed relation between changes in m if (left) and θSOI
i (right)
to changes in dQ c /d θ. An increase in m if results in an increase in accumulated
dQ ci /d θ, and a shift in θSOI
i results in a shift in dQ ci /d θ.
dQ ci /d θ, which gives
∂ dQ c Q LHV dQ ci
= R (12.23)
∂m if d θ dQ ci d θ
These assumptions are illustrated in Fig. 12.3. In both cases, the shape of dQ c /d θ
is preserved. It is also assumed that ∆u j and dQ ci /d θ are decoupled when i 6= j
and that ignition delays remain constant. These are approximations since subse-
quent injections are coupled both through ignition delay and rail pressure. The
shape of dQ c /d θ is also known to change slightly with m f and θSOI . If more accu-
rate combustion models are available, those could be incorporated for potential
controller-performance improvement. The approximations made could still be
motivated for small changes in ∆u, since feedback from subsequent engine cy-
cles corrects for unmodeled effects. Second derivatives, necessary for computing
the Hessians in (12.15) are given by
∂2 dQ c d 3Q ci
=
i
∂(θSOI )2 dθ d θ3
∂2 dQ c Q LHV d 2Q ci
= − R (12.24)
i
∂θSOI ∂m if d θ dQ ci d θ 2
∂2 dQ c
=0
∂(∂m if )2 dθ
228
12.2 Modeling and Heat-Release Detection
µ ¶T
dQ c Q LHV dQ c1 d 2Q c1 Q LHV dQ cM d 2Q cM
∇ = R 1 2
... R
dθ dQ c d θ dθ dQ cM d θ d θ2
Q LHV d 2Q ci
(12.25)
0 − R
dQ dQ i d θ2
2 c c
∇ = diag
dθ
Q LHV d 2Q ci d 3Q ci
−R i d θ2 3
dQ c dθ
From d ∆Q c /d θ to ∆p
The relation between d ∆Q c /d θ and ∆p can now be established with the lineari-
zed pressure model introduced in Chapter 11
µ ¶
d ∆p γ dV d µ(p 0 , θ) γ − 1 d ∆Q c
=− + ∆p + (12.26)
dθ V dθ dp V dθ
µ ¶
d ∆Q c dQ c T
= ∇ ∆u (12.27)
dθ dθ
Z θ d ∆Q c (ϑ)
∆p(θ) = Φ(θ, ϑ)Γ(ϑ) dϑ (12.28)
θIVC dϑ
where
µ Z θ ¶µ ¶γ
d µ(p 0 , τ) V (ϑ)
Φ(θ, ϑ) = exp − dτ
ϑ dp V (θ)
(12.29)
γ−1
Γ(ϑ) =
V (ϑ)
229
Chapter 12. Predictive Constraint Handling and Pressure Tracking
100
p0
p 0 + ∆p
80 p(u 0 + ∆u)
60
p [bar]
40
20
0
−20 −10 0 10 20 30 40 50
θ [CAD]
With (12.28), the gradients with respect to p, p IMEP and d p/d θ are given by
Z θ µ ¶
dQ c (ϑ)
∇p = Φ(θ, ϑ)Γ(ϑ) ∇ dϑ
θIVC dϑ
Z θEVO
1 (12.30)
∇p IMEP = ∇pdV
Vd θIVC
d p d ∇p
∇ =
dθ dθ
Hessians are given by
Z θ µ ¶
2 2 dQ c (ϑ)
∇ p= Φ(θ, ϑ)Γ(ϑ) ∇ dϑ
θIVC dϑ
(12.31)
Z θEVO
2 1 2
∇ p IMEP = ∇ pdV
Vd θIVC
Note that the gradients and Hessians above, except for ∇p IMEP and ∇2 p IMEP ,
are functions of θ.
NOx
The NOx constraint was linearized using the NOx -formation model in Sec 2.5.
Since this model is not easily linearized, partial derivatives of the cylinder-out
230
12.2 Modeling and Heat-Release Detection
10
Data Model
5
Tex [◦ C]
0
−5
−10
0 200 400 600 800 1,000 1,200 1,400 1,600 1,800
cycle [-]
200
TEVO [◦ C]
−200
Figure 12.5 Experimental Tex and TEVO data (solid) together with (12.33) output
(dashed) when varying θSOI . The temperature scales are relative to steady-state
values at p IMEP = 5 bar.
NOx concentration with respect to u was obtained by solving (2.73) and applying
numerical differentiation
In the forward step, NOx (u 0 + ∆u) was computed by solving (2.73) with modified
cylinder pressures, temperatures and heat-release rates p(u 0 + ∆u), T (u 0 + ∆u),
dQ c /d θ(u 0 + ∆u), obtained from the linearized expressions presented above.
With M injections, this amounts to solving (2.73) 2M + 1 times. The most
computationally-demanding part of computing NOx was to compute gas pro-
perties as a function of temperature. To reduce the computational load, it was
decided to also use the gas properties computed at u 0 in the forward steps.
231
Chapter 12. Predictive Constraint Handling and Pressure Tracking
Exhaust Temperature
Examples of lumped-parameter Tex models, suitable for control applications are
presented in [Eriksson and Nielsen, 2014]. This methodology was adopted here
and used to model the relation between the temperature at exhaust-valve ope-
ning (EVO) TEVO and Tex
where ∆Tex (k) and ∆TEVO (k) are temperature deviations from an equilibrium
0 0
point Tex , TEVO , and ΦTex and ΓTex are model parameters. Differentiation with
respect to TEVO
V (θEVO )TIVC
∇TEVO = ∇p(θEVO ) (12.34)
p IVC VIVC
establishes a relation with ∇p(θEVO ), which is given by (12.30).
Experimental Tex and TEVO data (solid) during θSOI step changes are presen-
ted together with (12.33) output (dashed) in Fig. 12.5. The temperature scales are
relative to steady-state values at p IMEP = 5 bar. Values for ΦTex and ΓTex were ob-
tained from engine data and the MATLAB system-identification toolbox.
In order to incorporate the long time constants of Tex into the MPC pro-
blem formulation, exhaust temperature was predicted over a longer horizon
T
H p ex = 500 and with a longer sampling interval, (n Tex = 100 cycles). Moreover, the
model was augmented with a disturbance state d TEVO to keep track of the model
steady-state offset
The Kalman filter presented in Chapter 3 was then used to estimate d TEVO using
(12.35) and Tex measurements.
QP Approximations
Constraint Fulfillment With gradients and Hessians available, the QP approxi-
mation of (12.3) is given by
Hp
X
minimize ˆ ˆ
J θCT (k) + J p IMEP (k) + J ∆u (k) (12.36)
∆U k=1
0 c
subject to p max + ∇p 0 (θ j )∆u(k) ≤ p max θ j = θ0 , . . . , θ f , k = 1, . . . , H p
0 c
d p max + ∇d p 0 (θ j )/d θ∆u(k) ≤ d p max θ j = θ0 , . . . , θ f , k = 1, . . . , H p
NO0x + ∇NO0x ∆u(k) ≤ NOcx k = 1, . . . , H p
0 c T
Tex + ∆Tex (k) ≥ Tex k = n Tex , . . . , H p ex
u 0 + ∆u(k) ∈ Û k = 1, . . . , H p
232
12.3 Adding and Removing Injections
where
r,i i ,0 2
(θCT − θCT )
f θiCT =
−2m if ,0 (θCT
r,i i ,0
− θCT )
(12.38)
r,i i ,0
0 −(θCT − θCT )
Hθi CT = α r,i i ,0
−(θCT − θCT ) m if ,0
Furthermore, the cost function Jˆp IMEP (k) is given by
r
¡ T
¢
Jˆp IMEP (k) = −β(p IMEP − p IMEP,0 ) 2∇p IMEP,0 ∆u(k) + ∆u T (k)∇2 p IMEP,0 ∆u(k)
+ β∆u T (k)∇p IMEP,0 ∇p IMEP,0
T
∆u(k) (12.39)
Index 0 in (12.36) to (12.39) above denotes cycle k = 0.
Pressure Tracking The QP approximation (12.7)
Hp
X
minimize Jˆp (k) + J ∆u (k) (12.40)
∆U k=1
subject to u 0 + ∆u(k) ∈ U, k = 1, . . . , H p
is obtained by computing the approximated pressure-tracking cost
θf
X ¡ ¢
Jˆp (k) = − αt r (p r (θ j ) − p 0 (θ j )) 2∇p 0T (θ j )∆u(k) + ∆u T (k)∇2 p 0 (θ j )∆u(k)
θ j =θ0
233
Chapter 12. Predictive Constraint Handling and Pressure Tracking
100 100
p [bar]
p [bar]
50 50
0 0
−20 0 20 40 −20 0 20 40
θ [CAD] θ [CAD]
Figure 12.6 The optimal number of injections changes depending on the ope-
rating conditions. To the right, a single fuel injection is optimal for tracking the
constant-volume combustion reference pressure (dashed). To the left, it is instead
optimal with three fuel injections for tracking of a limited-pressure cycle (dashed).
Constraint Fulfillment
The simulation results presented in Chapter 10 suggested the following logic for
adding and removing injections when solving (12.3) with up to three injections,
pilot, main and post:
• Add a pilot injection if any of the predicted outputs NOx , p max or d p/d θ
are larger than x c − ²add
x .
c
• Add a post injection if the predicted Tex is smaller than Tex + ²add
Tex .
• Remove the pilot injection if predicted NOx , p max and d p/d θ are smaller
than x c − ²rem 1
x , and the pilot-fuel mass m f is sufficiently low.
c
• Remove the post injection if predicted Tex is larger than Tex + ²rem
Tex , and the
post-fuel mass m 3f is sufficiently low.
2
was therefore introduced. The compensation in (12.42) adjusts θSOI (k + 1) when
a pilot injection is added or removed to compensate for variation in ignition de-
234
12.3 Adding and Removing Injections
100 100
p [bar]
p [bar]
50 50
0 0
−20 0 20 40 −20 0 20 40
θ [CAD] θ [CAD]
Pressure Tracking
A different method was adopted in (12.40) for adding and removing injections.
The method used for adding injections was to extrapolate from already active in-
jections. For example, if injection i is inactive, ∆u with respect to injection i is as-
sumed to correspond to variations of a shifted (∆θ) and rescaled (x i ) heat-release
rate dQ ca /d θ, corresponding to an already active injection
dQ ci dQ ca
= xi (θ − ∆θ) (12.43)
dθ dθ
235
Chapter 12. Predictive Constraint Handling and Pressure Tracking
x 0 + ∇x 0T ∆u(k) ≤ x c + ²x
(12.44)
²x ≥ 0
x 0 + ∇x 0 ∆u(k) ≤ x c + ²x + 2x σ
(12.46)
²x ≥ 0
The constraint margins were here pre-computed from measured output stan-
dard deviations, meaning that the margins provide probabilistic guarantees for
constraint fulfillment. An alternative approach would be to estimate x σ on-line.
The reason for introducing constraint margins was to remediate the problem of
stochastic constraint violation that occured in previous chapters.
Constraint Fulfillment
This section presents experimental results for MPC (12.36). The controller was
r r
tuned for θCT - and p IMEP -tracking response times within 10 engine cycles. Con-
troller gains and input constraints used are presented in Table 12.2. Gain sche-
duling was implemented by increasing R SOI (×4) for the controller to vary θSOI
more cautiously in the vicinity of a constraint limit. A prediction horizon of two
236
12.5 Experimental Results
engine cycles H p = 2 was used. This enabled mean QP-solver computation times
of 2 ms. The QP was constructed below 50 ms, where ≈ 90% of the time was used
to compute the NOx gradients in (12.32). Fulfillment of each constraint was eval-
r
uated separately during p IMEP step changes. Experimental results are presented
in the following sections.
Table 12.2 MPC weights and constraints for (12.36), units for θSOI and m f are
[CAD] and [mg].
237
Chapter 12. Predictive Constraint Handling and Pressure Tracking
20 1
θSOI 2
θSOI 3
θSOI
20
θ SOI [CAD]
θ CT [CAD] 10
0
0
−20
100 200 300 100 200 300
cycle [-] cycle [-]
150 12
10
p IMEP [bar]
100
m f [mg]
8
50
6
0 4
100 200 300 100 200 300
cycle [-] cycle [-]
90
p max [bar]
80
70
60
50 100 150 200 250 300
cycle [-]
Figure 12.8 Evaluation of p max -constraint fulfillment. The controller was set to
r
follow p IMEP step changes from 5 to 10 bar. As p IMEP was increased from 5 to 10
r = 8 CAD could not be maintained without viola-
bar with a single injection, θCT
c
ting p max = 80 bar. The controller acted by increasing m 1f as the constraint was
1 (red) shortly after TDC.
approached. This resulted in an additional detected θCT
2 (blue) was delayed 8 CAD for constraint fulfillment. Crank-angle
Moreover, θSOI
resolved data for cycles 45 and 80 are presented in Fig. 12.9.
Crank-angle resolved data for cycles 45 (left) and 80 (right) are presented in
Fig. 12.9. In-cylinder pressure, injector current and dQ c /d θ are presented toge-
r
ther with vertical lines indicating θCT and θCT . The dashed pressure curves vi-
sible at cycle 45 are predicted cylinder pressures for two subsequent engine cy-
r
cles. The increase in predicted pressure was due to the increase in p IMEP . Two
1 2
θCT are detected at cycle 80 due to the increase in m f and θSOI . The red and
238
12.5 Experimental Results
80 80
cycle 45 cycle 80
60 60
p [bar]
p [bar]
40 40
20 20
0 0
−20 0 20 40 −20 0 20 40
θ [CAD] θ [CAD]
Figure 12.9 In-cycle data for cycles 45 (left) and 80 (right) in Fig. 12.8. In-cylinder
pressure, injector current and dQ c /d θ are presented together with the vertical
r . The dashed pressure curves visible at cycle 45 are
lines indicating θCT and θCT
predicted cylinder pressures for two subsequent engine cycles. The increase in
r
predicted pressure was due to the increase in p IMEP . At cycle 80, it can be seen
how the controller separated dQ c /d θ between the pilot and main injection.
blue heat-release rates indicate how the combustion detection method separa-
ted dQ c /d θ into dQ c1 /d θ and dQ c2 /d θ.
d p/d θ The d p/d θ limit was varied at constant p IMEP in the experiment pre-
1,2
sented in Fig. 12.10. The controller increased m 1f and delayed θSOI c
as d p max was
decreased. The pilot and main fuel burned simultaneously and only one com-
bustion timing was therefore detected, meaning that dQ c /d θ was attributed to
and controlled by the main injection. Disturbances in p IMEP are visible when m 1f
changed. These disturbances could have been attenuated by improving the ca-
adj
libration of ∆m f in (12.42). The d p/d θ limit was occasionally violated due to
cycle-to-cycle variation.
In-cycle data for cycles 15 (left) and 150 (right) are presented in Fig. 12.11. The
c
solid and dashed tangents correspond to d p/d θmax and d p max , respectively. The
1 2 c
pilot fuel amount m f was increased and θCT was delayed to fulfill d p max . The
injected fuel burned simultaneously, which resulted in one combustion timing
detected.
r
NOx The controller was once again set to follow p IMEP step changes to evaluate
NOx -constraint handling where the solution to (2.73) was used as a virtual NOx
r 2
sensor, see Fig. 12.12. As p IMEP increased, θCT was delayed and m 1f increased for
NOx to remain below NOcx = 500 ppm. Constraint violation occurred due to NOx
overshoots as p IMEP was increased. The constraint was however later fulfilled in
steady state. The NOx overshoot could possibly be reduced by adding pilot fuel
239
Chapter 12. Predictive Constraint Handling and Pressure Tracking
20 20 1
θSOI 2
θSOI 3
θSOI
10
θ SOI [CAD]
θ CT [CAD] 10
0
0 −10
−20
100 200 300 100 200 300
cycle [-] cycle [-]
150 10
p IMEP [bar]
100
m f [mg]
50
6
0
100 200 300 100 200 300
cycle [-] cycle [-]
30
d p max [CAD/bar]
20
10
0
50 100 150 200 250 300
cycle [-]
c
Figure 12.10 Evaluation of d p/d θ constraint fulfillment where d p max was va-
1 1,2
ried at constant p IMEP . The pilot fuel amount m f was increased and θSOI were
delayed to fulfill the constraint. The pilot and main fuel burned simultaneously
and only one combustion timing was detected, meanin87g that dQ c /d θ was at-
tributed to and controlled by the main injection. Disturbances in p IMEP are visible
when m 1f changes around cycle 60 although compensation in m 2f was made.
more carefully, allowing for larger θSOI changes, and/or having a longer predic-
tion horizon. Small NOx overshoots could however be acceptable if NOcx was set
to regulate accumulated NOx emissions.
In-cycle data for cycles 25 (left) and 65 (right) from the experiment in
Fig. 12.12 are presented in Fig. 12.13. The solid and dashed purple lines corres-
pond to in-cylinder NOx and NOcx , respectively. The pilot-fuel mass m 1f was in-
2
creased and θCT delayed at cycle 65, in order for NOcx = 500 ppm to be fulfilled.
240
12.5 Experimental Results
80 80
cycle 15 cycle 150
60 60
p [bar]
p [bar]
40 40
20 20
0 0
−20 0 20 40 −20 0 20 40
θ [CAD] θ [CAD]
Figure 12.11 In cycle data for cycles 15 (left) and 150 (right) from the experiment
c
in Fig. 12.10. The solid and dashed tangents correspond to d p/d θmax and d p max
1 2
respectively. The controller increased m f and delayed θCT at cycle 150 to fulfill
c
d p max = 10 bar/CAD.
c
T ex Figure 12.14 shows how the controller managed to keep Tex above Tex =
◦ r
240 C , whilst p IMEP was varied between 5 and 7 bar. The controller introduced a
c
post injection as Tex approached Tex . The post-injection mass m 3f was then used
c
to regulate Tex above Tex . A p IMEP disturbance is visible when the post injection
was introduced at cycles 180 and 750.
In-cycle data from cycles 5 (left) and 350 (right) in Fig. 12.14 are presented
c
in Fig. 12.15. A post injection was introduced at cycle 350 to fulfill Tex = 240◦C
2 r 3
whilst θCT was kept at θCT . The post-injection combustion timing θCT and cor-
responding heat-release rate dQ c3 /d θ (green) were detected by the controller.
Summary
It can be concluded that the controller was able to fulfill the different constraints
as intended. Speed of convergence was higher than for the heuristic control-
ler design in Chapter 10 with 10-20 cycles as compared to 40-50 cycles. The
model-based controller presented here was also better at avoiding constraint
violation. Both because of its predictive capability, and because of its constraint
margins.
Pressure Tracking
This section presents experimental results with MPC (12.40). The controller was
investigated by varying the pressure-reference parameters αv , θSOC and Q cr . Con-
troller parameters used are presented in Table. 12.3.
The ability to follow αv changes with two injections is presented during a
10-cycle transition from p 1r with αv = 0.5 (red) to p 2r with αv = 0.2 (blue) in
241
Chapter 12. Predictive Constraint Handling and Pressure Tracking
20 20 1
θSOI 2
θSOI 3
θSOI
10
θ SOI [CAD]
θ CT [CAD] 10
0
0 −10
−20
200 400 200 400
cycle [-] cycle [-]
150 10
p IMEP [bar]
100 8
m f [mg]
50 6
0 4
200 400 200 400
cycle [-] cycle [-]
800
600
NOx [ppm]
400
200
r
Figure 12.12 Evaluation of NOx -constraint fulfillment. As p IMEP 2
increased, θCT
was delayed and m 1f increased for the constraint to be fulfilled. Constraint viola-
tion occured as p IMEP increased. The constraint was later fulfilled in steady state.
Fig. 12.16. With αv = 0.5, the controller had a relatively large m 1f . With more
constant-pressure combustion, the controller increased m 2f and delayed θSOI
1
2
and θSOI .
Figure 12.17 presents in-cycle data during a transition from p 1r with Q cr = 3000
J (red) to p 2r with Q cr = 6000 J (blue) with one injection. When Q cr was increased,
the controller increased m f whilst θCT was kept constant.
r
A θSOC transition from p 1r with θSOC
r
= 10 CAD (red) to p 2r with θSOC
r
= 5 CAD
r
(blue) is presented in Fig. 12.18. When θSOC was advanced, the controller adjus-
r
ted θSOI whilst m f was kept constant.
242
12.5 Experimental Results
80 80
cycle 25 cycle 65
60 60
p [bar]
p [bar]
40 40
20 20
0 0
0 20 40 60 0 20 40 60
θ [CAD] θ [CAD]
Figure 12.13 In-cycle data for cycles 25 (left) and 65 (right) from the experiment
in Fig. 12.12. The solid and dashed purple lines correspond to in-cylinder NOx
and NOcx respectively. The pilot-fuel mass m 1f was increased and θCT 2 delayed at
Data illustrating the suggested method for introducing injections during pres-
sure tracking is shown in Fig. 12.19. For this experiment, the engine was run with
diesel fuel. In the upper subdiagram, two injections were used to track a p r with
αv = 0.3. As αv was changed to 0.1 and Q cr was increased in the middle subdia-
gram in Fig. 12.19, the controller extrapolated (dashed) from the first detected
heat-release rate to apprehend how a post injection would affect p. The post in-
jection was then introduced, since it would decrease p r error cost, see the lower
subdiagram in Fig. 12.19.
Summary
The pressure-tracking controller was able to adjust fuel injection as p r parame-
ters were varied. The ratio between m 1f and m 2f was changed to adjust the com-
bustion duration as αv was varied. The injected fuel m f and θSOI were adjus-
ted to account for changes in Q cr and θSOC
r
, respectively. Error-free tracking could
not be obtained due to limited controllability and the steep p r increase during
constant-volume combustion. It is believed that improved tracking performance
could be obtained with a smoother p r , and an adjusted γ during the expansion
Table 12.3 Controller weights and constraints for (12.40), units for θSOI and m f
are [CAD] and [mg].
αt r = 1 R θSOI = 0.5 Rm f = 3
−0.5 ≤ ∆θSOI ≤ 0.5 −10 ≤ ∆m f ≤ 10 −25 ≤ θSOI ≤ 20 0 ≤ m f ≤ 120
243
Chapter 12. Predictive Constraint Handling and Pressure Tracking
30 1
θSOI 2
θSOI 3
θSOI
20
θ SOI [CAD]
20
θ CT [CAD]
10 0
0
−20
200 400 600 800 200 400 600 800
cycle [-] cycle [-]
150 8
p IMEP [bar]
100
m f [mg]
6
50
5
0 4
200 400 600 800 200 400 600 800
cycle [-] cycle [-]
350
300
Tex [◦ C]
250
200
150
100 200 300 400 500 600 700 800
cycle [-]
Figure 12.14 Evaluation of Tex -constraint fulfillment. The ability to fulfill Tex
was evaluated by decreasing p IMEP from 7 to 5 bar. A post-injection was intro-
duced when Tex approached Tex c = 240 ◦ C . The controller then regulated T by
ex
adjusting m 3f .
stroke. A fuel-injection system that allows for direct control of the injection rate
would improve controllability further. The controllability could also be improved
by increasing the fuel-injection rate, which would allow for an increased number
of injections in a shorter θ interval.
244
12.5 Experimental Results
80 80
p [bar]
40 40
20 20
0 0
−20 0 20 40 −20 0 20 40
θ [CAD] θ [CAD]
Figure 12.15 In-cycle data for cycles 5 (left) and 350 (right) from the experiment
c = 240◦ C
in Fig. 12.14. At cycle 350, a post injection was introduced to fulfill Tex
2 r 3
whilst θCT was kept at θCT . The post-injection combustion timing θCT and cor-
responding heat-release rate dQ c3 /d θ (green) were detected by the controller.
100
p 1r
80
p 2r
60
p [bar]
40
20
−20 −10 0 10 20 30 40
θ [CAD]
Figure 12.16 The controller’s ability to follow αv changes with two injections.
During this experiment, αv was changed from 0.5 (red) to (blue) 0.2. To obtain
more constant-pressure combustion, the controller increased m 2f and delayed
1 and θ 2 .
θSOI SOI
245
Chapter 12. Predictive Constraint Handling and Pressure Tracking
100
p 2r
80
60
p 1r
p [bar]
40
20
−20 −10 0 10 20 30 40
θ [CAD]
100
p 2r
80
60
p 1r
p [bar]
40
20
−20 −10 0 10 20 30 40
θ [CAD]
r
Figure 12.18 A θSOC transition from p 1r with θSOC
r = 10 CAD (red) to p 2r with
r r
θSOC = 5 CAD (blue). As θSOC was advanced, the controller adjusted θSOI whilst
m f was kept constant.
246
12.5 Experimental Results
100
80
60
p [bar]
40
20
0
−20 −10 0 10 20 30 40
θ [CAD]
100
80
60
p [bar]
40
20
0
−20 −10 0 10 20 30 40
θ [CAD]
100
80
60
p [bar]
40
20
0
−20 −10 0 10 20 30 40
θ [CAD]
Figure 12.19 The suggested method for introducing injections during pressure
tracking. In the upper subdiagram, two injections were used to track a p r with
αv = 0.3. As αv was changed to 0.1 and Q cr was increased in the middle sub-
diagram, the controller extrapolated (dashed) from the first detected heat-release
rate to apprehend how a third post injection would affect p. The post-injection
was introduced since it was found to decrease p r -tracking error, see the lower
subdiagram.
247
Chapter 12. Predictive Constraint Handling and Pressure Tracking
12.6 Discussion
The use of empirical, data-based and mean-valued models to describe
in-cylinder processes has been a common theme in previous works on opti-
mal engine control for constraint fulfillment, see [Hafner et al., 2000; Stewart and
Borelli, 2008; Atkinsson et al., 2009; Karlsson et al., 2010; Grahn et al., 2014].
This chapter combined the heat-release detection and separation method
presented in Chapter 4 with the pressure-prediction method in Chapter 11. These
methods allowed for a cycle-resolved MPC formulation, where the MPC, in con-
trast to previous work, can predict the effect from different injections on the cy-
linder pressure. To the author’s knowledge, this is a novel controller-design fra-
mework that can be used for both constraint fulfillment and pressure tracking.
Constraint Fulfillment
The proposed MPC for constraint handling in (12.36) worked as intended in
the experimental evaluation. It was shown in Chapter 10 that similar transient
behavior could be obtained with a simpler PI-controller design. However, the
model-based approach adopted in this chapter allowed the controller to predict
constraint violation and in that way act beforehand. This, in combination with
constraint margins, resulted in smaller NOx overshoots, and no p max overshoots,
as compared to the PI controller in Chapter 10. Furthermore, the MPC had shor-
ter settling times of 10-20 cycles, as opposed to 40-50 cycles for the PI controller.
The model-based controller accounts for variation in combustion characteris-
tics as a function of operating point which makes the feedback loop more robust.
The centralized MPC design also allowed the controller to take into account for
cross-coupling between control of load and other engine outputs.
Compared to the heuristic PI controller in Chapter 10, the MPC framework
provided a systematic way of handling input and output constraints. All con-
straints were accounted for simultaneously, where it was straightforward to add
or remove constraints as long as meaningful solutions existed. The controller was
also flexible when adding or removing injections and adjusted the size of the op-
timization problems accordingly. The model-based methods for adding and re-
moving injections could however be developed further by, for instance, also op-
timize timings and amounts of the injections introduced.
Even though MPC has its potential benefits, it demands careful tuning so
that desired controller behavior is obtained. Poor tuning could result in patho-
logical behavior, such as control of engine load with injection timing or control
of exhaust temperature with the main-injection duration, with an offset in load
tracking as a result. Further controller development would be to also incorpo-
rate gas-system dynamics and actuators. If more accurate combustion models
are available, those could be incorporated for potential controller-performance
improvement. This may however require a more sophisticated heat-release de-
tection method.
248
12.7 Conclusions
Combustion Detection
The combustion detection method worked well, apart from some exceptions at
cycle 60 in Fig. 12.10 and cycle 50 in Fig. 12.12. These errors occurred when a pilot
injection was introduced and the controller did not detect the pilot heat-release
rate properly. Instead, it detected the main-injection heat-release and a later
heat-release peak. One approach to further develop the heat-release detection
method would be to instead compute the likelihood of a detected peak being a
combustion timing. Such a methodology could make use of a heat-release model
and fuel-injection information.
Pressure Tracking
Although tracking of a cycle-resolved pressure reference is an unconventional
way of controlling an engine, the results presented in this chapter showed that
the presented MPC in (12.40) allowed for this approach. The convergence rate
presented here (10 to 15 engine cycles) was comparable to the results in [Zweigl
et al., 2015]. Control errors presented here were however somewhat larger than
those reported in [Jörg et al., 2015; Zweigl et al., 2015]. With a fuel-injection sys-
tems that allows for additional injections during the engine cycle, and more di-
rect control of the fuel-injection rate, it is possible that the controller presented
in this chapter could exhibit more accurate cylinder-pressure control.
If this controller design is more favorable than conventional cylinder-pressure
feedback controllers that regulate combustion-timing and indicated load re-
mains to be investigated.
12.7 Conclusions
Two model predictive controllers were presented and experimentally evalua-
ted. Both controllers utilized a linearized cylinder-pressure model and a novel
combustion-detection method in order to predict in-cylinder pressure variation
due to fuel-injection changes. Experimental results demonstrated:
In both cases, fuel-injections were added and removed depending on the predic-
ted in-cylinder pressure.
249
13
Conclusions and Future
Research
Demands for reduced emission levels and lowered fuel consumption have cre-
ated a need for accurate engine control. This has been illustrated in this thesis
through an investigation of how model-based closed-loop combustion control
can be used to improve the reliability of a low-emission combustion concept.
This thesis has also investigated how timings and durations of multiple injec-
tions can be decided with feedback control for efficient constraint fulfillment. In
both contexts, feedback control reduces the amount of calibration work needed
for efficient engine operation. It makes the combustion processes more robust
to changes in intake conditions, hardware aging and fuel properties, and lowers
the demands for precise actuators. The main results presented in the thesis are
summarized below, together with suggestions for future research.
PPC
This thesis investigated model-based control for improved operation of partially
premixed combustion (PPC). Designed controllers experimentally demonstrated
control of ignition delay, combustion timing and pressure-rise rate in transient
operation. Gas-exchange and fuel-injection actuations for improved low-load ef-
ficiency were also suggested.
The problems of regulating ignition delay and pressure-rise rate were studied
separately in chapters 7 and 8. Since there is an inherent trade-off between igni-
tion delay and pressure-rise rate, however, it would be interesting to investigate
concurrent control of these variables. The controller objective could be formula-
ted as a set-point tracking problem with respect to ignition-delay with an upper
bound on pressure-rise rate. This behavior could be obtained by combining the
model predictive controllers in (7.12) and (8.13).
The PPC experiments presented in the thesis were mainly limited to the
low-to-mid load operating range of the engine. This suggests that future work
should include a more detailed evaluation of controller performance at higher
250
Chapter 13. Conclusions and Future Research
engine loads. Experiments have, however, showed that it might be difficult to cre-
ate sufficiently long ignition delays at higher loads with the current experimental
setup, and that the combustion characteristics approach those of conventional
diesel combustion at higher loads. Manente et al. [2009] suggested that long ig-
nition delays could be achieved at high load with large injections, located early
during the compression stroke (θSOI < -50 CAD). This strategy was also difficult to
implement in the experimental setup used due to preignition of such injections.
Preignition avoidance is an interesting control problem related to PPC, that has
not been discussed in this thesis. It is possible that hardware adjustments, such
as a decreased compression ratio, decreased swirl, and fuels of even higher ON
would facilitate high-load PPC operation. It is, however, also possible that such
adjustments would make low-load performance more challenging.
Optimal Control
This thesis investigated how multiple injections should be actuated for efficient
fulfillment of constraints on pressure, NOx and exhaust temperature. A hybrid
multivariate PI controller was designed and experimentally evaluated. This con-
troller showed an experimental efficiency improvement of 4-5 % compared to a
single-injection controller, as restrictive constraints on pressure and NOx were
imposed.
A predictive pressure controller was introduced, where a simple pres-
sure model was used to directly control the in-cylinder pressure using model
predictive control. This controller was capable of tracking load and efficient
combustion-timing set points, as well as fulfilling constraints. The controller was
also able to track ideal-cycle pressure curves. Further development of this con-
troller include enhancement of physical model assumptions, design of a more
robust heat-release-detection method, and further investigation of model-based
conditions for varying the number of injections. It could also be interesting to
add the common-rail pressure as a control variable for increased controllability
of the heat-release shape.
Hardware that allows for additional injections and more direct control of the
fuel-injection rate could also increase the usefulness of the controller. In-cycle
control of the cylinder pressure with the use of similar predictive methods could
also be an interesting research topic. An FPGA would be a suitable option for such
an application, as the demand for computational speed increases.
Feedforward
The controllers presented in this thesis only utilized feedback control without
any feedforward action. It is, however, believed that the performance of the
controllers presented could be improved significantly with model-based feed-
forward control, especially with respect to response times in transient operation.
251
Chapter 13. Conclusions and Future Research
Design of model-based feedforward controllers that are compatible with the pre-
sented controllers is therefore an additional suggested research topic.
252
Bibliography
A, B. Kempinski, and J. Rife (1981). Knock in Spark Ignition Engines. SAE Technical
Paper 810147.
Agrell, F., H.-E. Ångström, B. Eriksson, J. Wikander, and J. Linderyd (2003). In-
tegrated Simulation and Engine Test of Closed Loop HCCI Control by aid of
Variable Valve Timings. SAE Technical Paper 2003-01-0748.
Akihama, K., Y. Takatori, K. Inagaki, S. Sasaki, and A. M. Dean (2001). Mechanism
of the Smokeless Rich Diesel Combustion by Reducing Temperature. SAE Te-
chnical Paper 2001-01-0655.
Annand, W. (1963). “Heat transfer in the cylinders of reciprocating internal com-
bustion engines”. Proc. Inst. of Mech. Eng. 177:1, pp. 973–996.
Aoyagi, Y., T. Kamimoto, Y. Matsui, and S. Matsuoka (1980). A Gas Sampling Study
on the Formation Processes of Soot and NO in a DI Diesel Engine. SAE Techni-
cal Paper 800254.
Askin, A. C., G. E. Barter, T. H. West, and D. K. Manley (2015). “The heavy-duty
vehicle future in the United States: a parametric analysis of technology and
policy tradeoffs”. Energy Policy 81, pp. 1–13.
Åström, K. J. and T. Hägglund (2006). Advanced PID Control. ISA – The Instrumen-
tation, Systems, and Automation Society, Research Triangle Park, NC, USA.
Åström, K. J. and R. M. Murray (2010). Feedback Systems: an Introduction for
Scientists and Engineers. Princeton University Press, Princeton, NJ, USA.
Atkinson, C. M., M. Allain, Y. Kalish, and H. Zhang (2009). Model-Based Con-
trol of Diesel Engines for Fuel Efficiency Optimization. SAE Technical Paper
2009-01-0727.
Bemporad, A., M. Morari, V. Dua, and E. N. Pistikopoulos (2002). “The explicit
linear quadratic regulator for constrained systems”. Automatica 38:1, pp. 3–
20.
253
Bibliography
254
Bibliography
255
Bibliography
Eichmeier, J., U. Wagner, and U. Spicher (2012). “Controlling gasoline low tem-
perature combustion by diesel micro pilot injection”. J. Eng. for Gas Turbines
and Power 134:7.
Ekholm, K., M. Karlsson, P. Tunestål, R. Johansson, B. Johansson, and P. Strandh
(2008a). “Ethanol-diesel fumigation in a multi-cylinder engine”. SAE Int. J.
Fuels and Lubricants 1, pp. 26–36.
Ekholm, K., M. Karlsson, P. Tunestål, R. Johansson, B. Johansson, and P. Strandh
(2008b). “Ethanol-diesel fumigation in a multi-cylinder engine”. SAE Int. J.
Fuels Lubricants 1, pp. 26–36.
Enkvist, P., T. Nauclér, and J. Rosander (2007). “A cost curve for greenhouse gas
reduction”. McKinsey Quarterly 1, p. 34.
Epping, K., S. Aceves, R. Bechtold, and J. E. Dec (2002). The Potential of HCCI
Combustion for High Efficiency and Low Emissions. SAE Technical Paper
2002-01-1923.
Eriksson, L. (1998). Requirements for and a Systematic Method for Identifying
Heat-Release Model Parameters. SAE Technical Paper 980626.
Eriksson, L. and L. Nielsen (2014). Modeling and Control of Engines and Driveli-
nes. John Wiley & Sons, Chichester, West Sussex, England.
Eriksson, L. and M. Sivertsson (2016). “Calculation of optimal heat release rates
under constrained conditions”. SAE Int. J. Engines 9:2, pp. 1143–1162.
Eriksson, L. and A. Thomasson (2017). “Cylinder state estimation from measu-
red cylinder pressure traces-a survey”. IFAC-PapersOnLine 50:1, pp. 11029–
11039.
European Commission (2016). A European Strategy for Low-Emission Mobility.
COM(2016) 501. Brussels, Belgium.
European Union (2007). “Cleaner trucks and buses: tighter limits for nitrogen oxi-
des and particulate matter (Euro VI).” Press release, IP/07/1989.
Fenimore, C. (1971). “Formation of nitric oxide in premixed hydrocarbon fla-
mes”. In: Proc. Int. Symposium on Combustion. Vol. 13. 1. Pittsburgh, PA,
pp. 373–380.
Fieweger, K., R. Blumenthal, and G. Adomeit (1997). “Self-ignition of SI engine
model fuels: a shock tube investigation at high pressure”. Combustion and
Flame 109:4, pp. 599–619.
Fiveland, S. B. and D. N. Assanis (2001). Development of a two-zone HCCI com-
bustion model accounting for boundary layer effects. SAE Technical Paper
2001-01-1028.
Fridriksson, H., B. Sundén, S. Hajireza, and M. Tunér (2011). CFD Investigation of
Heat Transfer in a Diesel Engine with Diesel and PPC Combustion Modes. SAE
Technical Paper 2011-01-1838.
256
Bibliography
257
Bibliography
258
Bibliography
IPCC (2014). “Climate change 2014, synthesis report”. IPCC Fifth Assessment Re-
port.
Isermann, R. (2014). Engine Modeling and Control. Springers Berlin Heidelberg,
Berlin, Germany.
Johansson, B. (2006). Förbränningsmotorer. Institutionen för värme-och kraftte-
knik, Lunds Tekniska Högskola, Lund, Sweden.
Johansson, R. (1993). System Modeling and Identification. Prentice Hall, En-
glewood Cliffs, NJ, USA.
Johansson, T. B., P. Kågesson, H. Johansson, L. Jonsson, J. Westin, H. Hejenstedt,
O. Hådell, K. Holmgren, and P. Wollin (2013). “Fossilfrihet på väg”. Ministry of
Enterprise, SOU 2013:84. Stockholm, Sweden.
Jones, J. C. P. and J. Frey (2015). “Threshold optimization and performance eval-
uation of a classical knock controller”. SAE Int. J. Engines 8:3, pp. 1021–1028.
Jörg, C., T. Schnorbus, S. Jarvis, B. Neaves, K. Bandila, and D. Neumann (2015).
“Feedforward control approach for digital combustion rate shaping reali-
zing predefined combustion processes”. SAE Int. J. Engines 8:2015-01-0876,
pp. 1041–1054.
Julier, S. J. and J. K. Uhlmann (2004). “Unscented filtering and nonlinear estima-
tion”. Proc. IEEE 92:3, pp. 401–422.
Jury, E. I. (1964). “Theory and application of the Z-transform method”.
Kalghatgi, G. T., P. Risberg, and H.-E. Ångström (2006). Advantages of Fuels with
High Resistance to Auto-Ignition in Late-Injection, Low-Temperature, Com-
pression Ignition Combustion. SAE Technical Paper 2006-01-3385.
Källkvist, K. (2011). Fuel Pressure Modelling in a Common-Rail Direct Injection
System. MSc Thesis LiTH-ISY-EX–11/4488–SE. Linköping University, Dept. of
Electrical Engineering, Linköping, Sweden.
Kalman, R. (1959). “On the general theory of control systems”. IRE Trans. Auto-
matic Control 4:3, pp. 110–110.
Kalman, R. E. et al. (1960). “Contributions to the theory of optimal control”. Bol.
Soc. Mat. Mexicana 5:2, pp. 102–119.
Kalman, R. E. (1960). “A new approach to linear filtering and prediction pro-
blems”. J. Fluids Eng. 82:1, pp. 35–45.
Kamimoto, T. and M.-h. Bae (1988). High Combustion Temperature for the Reduc-
tion of Particulate in Diesel Engines. SAE Technical Paper 880423.
Kanda, T., T. Hakozaki, T. Uchimoto, J. Hatano, N. Kitayama, and H. Sono (2005).
PCCI Operation with Early Injection of Conventional Diesel Fuel. SAE Techni-
cal Paper 2005-01-0378.
Karlsson, M., K. Ekholm, P. Strandh, R. Johansson, and P. Tunestål (2008). “LQG
control for minimization of emissions in a diesel engine”. In: Proc. IEEE Int.
Conf. Control Applications (CCA 2008). San Antonio, TX, USA, pp. 245–250.
259
Bibliography
260
Bibliography
261
Bibliography
262
Bibliography
Olsson, J.-O., P. Tunestål, and B. Johansson (2004). Boosting for High Load HCCI.
SAE Technical Paper 2004-01-0940.
Onishi, S., S. H. Jo, K. Shoda, P. Do Jo, and S. Kato (1979). Active
Thermo-Atmosphere Combustion (ATAC)-a New Combustion Process for
Internal Combustion Engines. SAE Technical paper, 790501.
Osuka, I., M. Nishimura, Y. Tanaka, and M. Miyaki (1994). Benefits of new Fuel In-
jection System Technology on Cold Startability of Diesel Engines-Improvement
of Cold Startability and White Smoke Reduction by Means of Multi Injection
with Common Rail Fuel System ECD-U2. SAE Technical Paper 940586.
Ott, T., F. Zurbriggen, C. H. Onder, and L. Guzzella (2013). “Cylinder individual
feedback control of combustion in a dual fuel engine”. In: Advances in Auto-
motive Control (ACC 2013). Vol. 7. 1. Tokyo, Japan, pp. 600–605.
Powell, J. D. (1993). “Engine control using cylinder pressure: past, present, and
future”. J. Dynamic Systems, Measurements, and Control 115:2B, pp. 343–350.
Randolph, A. L. (1990). Methods of Processing Cylinder-Pressure Transducer Sig-
nals to Maximize Data Accuracy. SAE Technical Paper 900170.
Rassweiler, G. M. and L. Withrow (1938). Motion Pictures of Engine Flames Corre-
lated with Pressure Cards. SAE Technical Paper 380139.
Roelle, M. J., N. Ravi, A. F. Jungkunz, and J. C. Gerdes (2006). “A dynamic model of
recompression HCCI combustion including cylinder wall temperature”. In:
Proc. IMECE. Chicago, IL, USA, pp. 415–424.
Saracino, R., M. R. Gaballo, S. Mannal, S. Motz, A. Carlucci, and M. Benegiamo
(2015). Cylinder Pressure-Based Closed Loop Combustion Control: A Valid
Support to Fulfill Current and Future Requirements of Diesel Powertrain Sys-
tems. SAE Technical Paper 2015-24-2423.
Shaver, G. M., J. C. Gerdes, and M. Roelle (2004). “Physics-based closed-loop con-
trol of phasing, peak pressure and work output in HCCI engines utilizing va-
riable valve actuation”. In: Proc. American Control Conf. (ACC 2004). Vol. 1.
Boston, MA, USA, pp. 150–155.
Sher, E. (1998). Handbook of Air Pollution from Internal Combustion Engines: Pol-
lution Formation and Control. Academic Press, San Diego, CA, USA.
Shi, Y., H.-W. Ge, and R. D. Reitz (2011). Computational Optimization of Internal
Combustion Engines. Springer Science & Business Media, London, England.
Sjöberg, M. and J. E. Dec (2005). “An investigation into lowest acceptable com-
bustion temperatures for hydrocarbon fuels in HCCI engines”. Proc. Combus-
tion Institute 30:2, pp. 2719–2726.
Solaka, H., U. Aronsson, M. Tunér, and B. Johansson (2012). Investigation of Par-
tially Premixed Combustion Characteristics in Low Load Range with Regards
to Fuel Octane Number in a Light-Duty Diesel Engine. SAE Technical Paper
2012-01-0684.
263
Bibliography
Solsjö, R. (2014). Large Eddy Simulation of Turbulent Combustion in PPC and Die-
sel Engines. PhD thesis TMHP-14/1107, Dept. Energy Sciences. PhD thesis.
Lund University, Lund, Sweden.
Spadaccini, L. J. and J. A. Tevelde (1982). “Autoignition characteristics of
aircraft-type fuels”. Combustion and Flame 46, pp. 283–300.
Sperling, D. and D. Gordon (2009). Two Billion Cars: Driving Toward Sustainabi-
lity. Oxford University Press, New York, NY, USA.
Stas, M. J. (1996). Thermodynamic Determination of TDC in Piston Combustion
Engines. SAE Technical Paper 960610.
Stewart, G. and F. Borrelli (2008). “A model predictive control framework for in-
dustrial turbodiesel engine control”. In: Proc. Conf. Decision and Control.
Cancun, Mexico, pp. 5704–5711.
Szekely, G. A., A. S. Solomon, and P.-H. Tsai (2004). Optimization of the
Stratified-Charge Regime of the Reverse-Tumble Wall-Controlled Gasoline
Direct-Injection Engine. SAE Technical Paper 2004-01-0037.
Takeda, Y. and N. Keiichi (1996). Emission Characteristics of Premixed Lean Diesel
Combustion with Extremely Early Staged Fuel Injection. SAE Technical Paper
961163.
Tanov, S., R. Collin, B. Johansson, and M. Tunér (2014). “Combustion stratifica-
tion with partially premixed combustion, PPC, using NVO and split injection
in a LD-diesel engine”. SAE Int. J. Engines 7:4, pp. 1911–1919.
Tsurushima, T., E. Kunishima, Y. Asaumi, Y. Aoyagi, and Y. Enomoto (2002). The
Effect of Knock on Heat Loss in Homogeneous Charge Compression Ignition
Engines. SAE Technical Paper 2002-01-0108.
Tunestål, P. (2011). “TDC offset estimation from motored cylinder pressure data
based on heat release shaping”. Oil & Gas Science and Technology–Revue
d’IFP Energies nouvelles 66:4, pp. 705–716.
Tunestål, P., J. K. Hedrick, and R. Johansson (2001). “Model-based estimation of
cylinder pressure sensor offset using least-squares methods”. In: Proc. Conf.
Decision and Control. Vol. 4. Orlando, FL, USA, pp. 3740–3745.
Turns, S. R. et al. (1996). An Introduction to Combustion. McGraw-Hill, New York,
NY, USA.
Wahlström, J. and L. Eriksson (2011). “Modelling diesel engines with a
variable-geometry turbocharger and exhaust gas recirculation by optimi-
zation of model parameters for capturing non-linear system dynamics”.
Proc. Inst. Mech. Eng., Part D: J. Automobile Eng. 225:7, pp. 960–986.
Wang, Y. and S. Boyd (2010). “Fast model predictive control using online optimi-
zation”. IEEE Trans. Control Systems Technology 18:2, pp. 267–278.
264
Bibliography
265
Bibliography
Zweigel, R., F. Thelen, D. Abel, and T. Albin (2015). “Iterative learning approach for
diesel combustion control using injection rate shaping”. In: Proc. European
Control Conf. (ECC 2015), pp. 3168–3173.
266
Videos
Videos showing cycle-resolved data for some of the experiments presented are
available in the entry for this thesis in the Lund University Research Portal:
http://portal.research.lu.se/portal/
267