Project report
Project 1 – TSFS09

Name1 – email1
Name2 – email2

September 10, 2015
This report template includes project 1 in the course Modeling and Control of Engines and Drivelines (TSFS09). The purpose of this template is to assist you as you will present your findings and conclusions. It can be difficult to know how detailed the report should be, therefore some examples are included in this template.

Another major point of this template is that you will have everything at place. It will become easier for us to inform you and help you if you get stuck, but most of all it helps yourselves. What you write in Chapter 2 and what you do in preparatory tasks, is that you will build the rest of the work on. Just remember to bring your report and look at it when you are working on the project!

The template already contains a lot of text. Anything which is written in normal font is text that you can leave to remain. Text written in italics like this, you shall delete as it is only meant to be helpful to you while you write the report.

If you want to write in another word processing tool that is obviously fine, but the different sections and subsections must be well separated and easy to find.

You can choose to write your report in Swedish or English.

There are many constant parameters which are required in the models, it is your task to differentiate between the ones that need to be obtained from the measurements and those which relate to e.g. gas properties.

NOTE: In addition to the measurements required in Project 1, we also need to perform measurements for Project 2. These measurements are added with brief explanatory text in sections 2.8, and 2.9 below, and performed here to avoid duplicate measurement occasions. The star () before the section name indicates that these measurements are for Project 2. The projects 2 measurements are included in projects 1 report only the experiment plans (so you know what to measure and how), and thus not the implementation and validation of the models estimated from these measurements. Implementation and Validation for these components should instead be reported in project report for Project 2.

Contents

1 Introduction
2 Model description and experiment plan
 2.1 Accelerator pedal
  2.1.1 Model description
  2.1.2 Experiment plan
 2.2 Throttle
  2.2.1 Model description
  2.2.2 Experiment plan
 2.3 Intake manifold
  2.3.1 Model description
  2.3.2 Experiment plan
 2.4 Fuel injector
  2.4.1 Model description
  2.4.2 Experiment plan
 2.5 Cylinder
  2.5.1 Model description
  2.5.2 Experiment plan
 2.6 Exhaust system
  2.6.1 Model description
  2.6.2 Experiment plan
 2.7 Lambda sensor
  2.7.1 Model description
  2.7.2 Experiment plan
 2.8 Compressor- and turbine impeller diameters
 2.9 Wastegate
  2.9.1 Model description
  2.9.2 Experiment plan
 2.10 Complete engine
3 Stationary analysis and model validation
 3.1 In-cylinder analysis
  3.1.1 pv diagram
  3.1.2 Otto cycle
  3.1.3 Specific fuel consumption
 3.2 Model validation
  3.2.1 Accelerator pedal
  3.2.2 Throttle
  3.2.3 Intake manifold
  3.2.4 Fuel injector
  3.2.5 Cylinder
  3.2.6 Exhaust system
  3.2.7 Lambda sensor
 3.3 Table of Parameters
4 Implementation of engine model
 4.1 Accelerator pedal
 4.2 Throttle
 4.3 Intake manifold
 4.4 Fuel injector
 4.5 Cylinder
 4.6 Exhaust system
 4.7 Lambda sensor
 4.8 Complete engine model
5 Controller
6 Experiments
 6.1 Acceleration test
 6.2 Test run on part of drive cycle
 6.3 Simulation of complete drive cycle
 6.4 Fuel consumption
 6.5 Emissions
7 Conclusions

Chapter 1
Introduction

This project report deals primarily with a turbocharged engine. A mathematical model for the engine is to be developed and implemented, its model parameters identified from measured data, and the various models are validated. A controller in connected to the fuel injection system, which uses both feedforward and feedback from the discrete lambda sensor. Finally, the engine is coupled with a number of ready-made modules containing driver, clutch and driveline models. This is in order to perform simulations of the complete vehicle. A series of experiments are conducted with the aim to study the emission formation and consumption. The tool that is used for implementation is Matlab combined with Simulink.

Besides working on the engine model, the report addresses a series of studies, whose purpose is to illustrate how the energy stored in fuel is converted into a driving torque during the engine operating cycle, and the losses that may occur.

Chapter 2
Model description and experiment plan

Project 1A.
In this chapter the engine is divided into a number of components, which are then modeled separately. For each component specified input and output signals, mathematical model which describes the behavior and the parameters are determined. It is then described how the model parameters will be identified. An overview of all input signals to the actuators and output signals from the sensors which should be included in project 1 is shown in Figure 2.1 where the red signals and components will be used in project 2 for a turbocharged engine.


PIC

Figure 2.1: An overview of the air path with signals and components to be modeled in project 1 (shown in black) and project 2 (shown in red+black)


Mathematical models are going to be used for the modeling of following engine sub components:

Because the simulations should be done with both warm and cold engine, it is also explored how long it takes for the catalyst to become warm. Sensors that measure air mass flow and pressure in the intake manifold are assumed to have negligible dynamics, while only the delays at lambda measurement are taken into account and exhaust gas mixing dynamics in the exhaust manifold are neglected. The control signals which are the inputs to the engine model are:

The control signals that are going to be mentioned in the experiment plans, must be selected only among these. During measurements in the lab, control signals are the ones which can be controlled from the control room.

The entire engine model block in Simulink should produce the following output signals

For your convenience, a complete model description and experimental plan for the accelerator pedal, and most of model description for the throttle are already included in the report template. Use this as a template when you write model description and experimental plan for other components. You should have replaced every ???? in this template before submitting your report.

2.1 Accelerator pedal

2.1.1 Model description

Model input signal(s):
Accelerator pedal position (βped),
Model output signal(s):
Throttle angle (αth)
Model:
The throttle angle is calculated by the following differential equation:
α˙th = -1(βped - αth)
      τth
(2.1)

where βped and αth can changes continuously between 0 and 1.

Parameter to find (via measurement):
Time constant of the throttle (τth).
Other parameters in the model:
-

2.1.2 Experiment plan

Signal(s) to be measured:
Throttle angle (αth)
Control signal(s):
Accelerator pedal position (βped)
Type of experiment:
Dynamic.
Measurement description:
While engine runs with a constant speed at an operating point many steps (up-down-up-down-...) are performed in accelerator pedal position. The times between steps are selected such that measured signals become stable before the next step is done.
How to find the parameter(s):
Since it is assumed that the system has first order dynamics, the step response expressed in time domain is:
        (      t-t )
αth(t) = 1- e- τth0 (βped,new- βped,old)+ βped,old
(2.2)

where t0 is the time point when the step happens. At time t = τth + t0 is

(     t-t)
 1- e--τt0h  = 1- e-1 ≈ 0.63
(2.3)

Therefore the time constant is determined as the time interval between the instance when the first change in the signal is observed until throttle angle reaches 63% of its final values.

2.2 Throttle

The air mass flow past the throttle depends on throttle angle and ambient pressures and temperatures. In most of the engine operating points it is a significant pressure drop across the throttle which is well described by compressible flow models.

˙m  = √-pus-A  (α  )Ψ(Π),  Π = pds
 at    Rus  eff  th            pus

which is also found in the textbook along with the expression for Ψ(Π). The only thing that is unknown is how the effective area Aeff depends on throttle plate angle αth. Here we use ...

Several models are proposed in the textbook, here a quadratic polynomial in αth is recommended.

2.2.1 Model description

Model input signal(s):
Throttle angle (αth), Pressure before throttle (pbef,thr), Pressure after throttle (is equal to Intake manifold pressure pim).
Model output signal(s):
Air mass flow rate (at)
Model:
Air mass flow rate past the throttle is given by:
at =  pbef,thr
∘R-------
   bef,thrAeff(αth)Ψ(Π), Π =   pim
p-----
 bef,thr (2.4)
Aeff(αth) =???? (2.5)
Ψ =   -------------------
∘   2γ (  -2     γ+1)
   ----- Πγlim - Πliγm
   γ - 1, Πlim = max(   (  2  ) γ-γ1)
  Π, -----
     γ + 1 (2.6)
Parameter to find:
Parameters ???? in the effective area (Aeff) model, temperature before throttle (bef,thr).
Other parameters in the model:
Gas constant (R) and ratio between specific heat capacities, γ = cp∕cv.

2.2.2 Experiment plan

Signal(s) to be measured:
Throttle angle (αth), Air mass flow (at), Intake manifold pressure (pim), Pressure before throttle (pbef,thr)
Control signal(s):
Accelerator pedal position (βped)
Type of experiment:
????
Measurement description:
????
How to find parameters:
Since it is important that the model is well suited    when there is a significant pressure drop across the throttle plate only    points with Π < 0.73 are used in the parameter estimation. ????

2.3 Intake manifold

Pressure dynamics in intake manifold. ????????????????????

2.3.1 Model description

Model input signal(s):
????
Model output signal(s):
????
Model:
????
Parameter to find:
????
Other parameters in the model:
????

2.3.2 Experiment plan

Signal(s) to be measured:
????
Control signal(s):
????
Type of experiment:
????
Measurement description:
????
How to find parameters:
????

2.4 Fuel injector

????????????????????

2.4.1 Model description

Model input signal(s):
????
Model output signal(s):
????
Model:
????
Parameter to find:
????
Other parameters in the model:
????

2.4.2 Experiment plan

Signal(s) to be measured:
????
Control signal(s):
????
Type of experiment:
????
Measurement description:
????
How to find parameters:
????

2.5 Cylinder

In the cylinder model there is one sub-model for the generated torque, one sub-model for the engine temperature and another one for the air/fuel ratio inside the cylinder (λcyl).
tips: Equation 7.55 in the book is recommended for torque modeling, read chapter 7.9.[1,2,3] before doing this part. For friction modeling use the Heywood model on page 191 but you have to find the model constants. For the temperature model, read the project compendium.

2.5.1 Model description

Model input signal(s):
????
Model output signal(s):
????
Model:
????
Parameter to find:
????
Other parameters in the model:
????

2.5.2 Experiment plan

Signal(s) to be measured:
????
Control signal(s):
????
Type of experiment:
????
Measurement description:
????
How to find parameters:
????

2.6 Exhaust system

Exhaust system is comprised of a model for exhaust manifold pressure and also the light-off time for the catalyst.
Considering Figure 2.1, the fact that the Enginemap data is obtained from a turbocharged engine and that a naturally aspirated engine is going to be modeled, which one of Pem or Pes signals in the Enginemap should be used in this task? Why?

tips: The output mass flow from the exhaust manifold can be modeled as an incompressible turbulent restriction, see chapter 7.2 in the book.

????????????????????

2.6.1 Model description

Model input signal(s):
????
Model output signal(s):
Pem
Model:
????
Parameter to find:
????
Other parameters in the model:
????

2.6.2 Experiment plan

Signal(s) to be measured:
????
Control signal(s):
????
Type of experiment:
????
Measurement description:
????
How to find parameters:
????

2.7 Lambda sensor

models for continuous and discrete lambda sensors, before catalyst are required.

2.7.1 Model description

Model input signal(s):
????
Model output signal(s):
????
Model:
????
Parameter to find:
τdbc,τλ
Other parameters in the model:

2.7.2 Experiment plan

Signal(s) to be measured:
????
Control signal(s):
????
Type of experiment:
????
Measurement description:
????
How to find parameter(s):
????

2.8 Compressor- and turbine impeller diameters

As part of the compressor model for Project 2, the outer wheel diameter of compressor impeller (dC2) is needed. This can most easily be measured by using calipers. In the same way, the turbine inlet diameter (dT1) should be measured with calipers.

Parameter to find:
Outer diameter of compressor wheel dC2and outer diameter of turbine wheel dT1.

2.9 Wastegate

Wastegate position of a turbocharged engine is normally controlled by the pressure on one side of a membrane. Pressure on other side of the membrane is normal ambient pressure. The pressure is controlled by using a mixing valve, where two sources with different pressure alternately are connected to the pressure chamber of the wastegate. This mixing valve is controlled by a PWM signal. An overview of the various components included in the wastegate controller on the turbocharged engine lab is shown in Figure 2.2.


PIC

Figure 2.2: Snapshot of the wastegate actuation system. A PWM signal controls the pressure in an actuator chamber, by alternately (and quickly), connecting the chamber to two different pressure sources. The pressure in actuator chamber then actuates a membrane that controls actuation lever. A return spring also serves on the actuator chamber which is why a greater pressure difference between the membrane different sides is required for larger change in position of the arm.


The PWM signal alternates between +5 V and ground, and the relative length difference between the +5 V level and ground-level controls the valve behavior. The frequency of the PWM signal is so high that an almost continuous signal is obtained (ie change between +5 V and ground occurs with the order of hundred Hz, because wastegate system frequency is significantly lower). The position of the wastegate arm can be measured in lab, and the PWM signal to the mixing valve can be controlled if wanted. A first order behavior, as for throttle angle, is assumed here. Model description and experiment plan for the wastegate will therefore largely follow the accelerator pedal.

2.9.1 Model description

Model input signal(s):
PWM signal for wastegate valve (PWMwg)
Model output signal(s):
Position of the wastegate (poswg)
Model:
Similar to the accelerator pedal model the wastegate position can be calculated by ????
Parameter to find:
τwg
Other parameters in the model:
????

2.9.2 Experiment plan

Signal(s) to be measured:
????
Control signal(s):
????
Type of experiment:
????
Measurement description:
????
How to find parameters:
????

2.10 Complete engine

Now that models are described and input-output-control signals are defined for all components, the block diagram of entire engine model can be illustrated as shown in Figure 2.3.


PIC

Figure 2.3: Block diagram of complete engine. Draw a block diagram of the entire engine, in which you can clearly see how the signals (the physical quantities) flow between the various sub-blocks. Having the block diagram in this figure right will later, in project 1C, help you when implementing the models in Simulink.


Chapter 3
Stationary analysis and model validation

Project 1B.
The enginemap can be downloaded from:
www.fs.isy.liu.se/Edu/Courses/TSFS09/Projekt1/EngineTests/

This chapter outlines and discusses the results of projects 1B. The aim in projects 1B is to show how engine torque can be produced from chemically stored energy, and where losses can occur. Approach is described by presenting the utilized equations and the name of the used data file(s). The results of parameter estimations in models of different engine components are also reported and the models are validated against measurements.

This chapter will be submitted as the second submission, i.e. after the project 1B has been implemented. Attach the Matlab code.

Preparatory tasks need not be explained in the report, but be aware that they can provide valuable clues as you e.g. explain figures.

3.1 In-cylinder analysis

3.1.1 pv diagram

3.1.2 Otto cycle

3.1.3 Specific fuel consumption

3.2 Model validation

In this section, the models developed in previous chapter are compared against the measurement data.

3.2.1 Accelerator pedal

The accelerator pedal takes pedal position (βped) as input and provides throttle angle (αth) as output. The model is described by the equation:

       1
α˙th = ---(βped - αth)
      τth
(3.1)

Parameter(s) value:
Time constant for the throttle τth in (3.1) is estimated to be 0.???? sec.
Measurement file:
throttleStep.mat.
Method:
The time constant is estimated as the time that it takes for the step to reach 63% of the final value as explained in the experiment plan 2.1.2.
Validation:
In Figure 3.1 we see that the model fits the measured step response. The rise time certainly depends on how big of a step is made ??but this will not have any significant impact on experiments done in this project.


PIC

Figure 3.1: Estimation of time constant for throttle. Note: This will be replaced by a figure from you.


3.2.2 Throttle

The model of the throttle determines the air mass flow rate into the intake manifold based on throttle angle. Model equation presented in Section 2.2, is repeated here

at =  pamb
√R-----
   ambAeff(αth)Ψ(Π), Π = pim
p---
 amb
Aeff(αth) = ...
Ψ = ∘ -------------------
   2γ  (  2γ     γ+γ1)
  γ --1 Π lim - Π lim, Πlim = max(           γ )
    (  2  )γ-1
 Π,  γ-+-1
What needs to be determined is the parameters of the polynomial for the effective area.
Parameter(s) value:
????
Measurement file:
????
Method:
For each measurement point the effective throttle area is calculated, then the parameters are adjusted with least squares method to the measured data.
Validation:
The model and measurements are compared in Figure ??. The figure shows that the modeled and measured values are in a good agreement in the range ???? but not in ???? range, and ...

3.2.3 Intake manifold

Parameter(s) value:
????
Measurement file:
????
Method:
????
Validation:
????

3.2.4 Fuel injector

Parameter(s) value:
????
Measurement file:
????
Method:
????
Validation:
????

3.2.5 Cylinder

Parameter(s) value:
????
Measurement file:
????
Method:
????
Validation:
????

3.2.6 Exhaust system

Parameter(s) value:
????
Measurement file:
????
Method:
????
Validation:
????

3.2.7 Lambda sensor

Parameter(s) value:
????
Measurement file:
????
Method:
????
Validation:
????

3.3 Table of Parameters

After validation of models, all estimated parameters in different sub models must be presented in this table. This makes it easier for us to correct your reports, and for you to reach and find different values during model implementation in Simulink.

All estimated parameters in the different engine sub models are summarized in table 3.1.


Table 3.1: Estimated parameters in engine sub models




Component Parameter Value






Accelerator pedal τth 0.??? [sec]



Throttle . ...
. ...



Intake manifold . ...
. ...



Fuel injector . ...
. ...



Cylinder . ...
. ...
. ...



Exhaust system . ...
. ...



Lambda sensor . ...

Number of empty rows in front of each component is randomly selected and should not be inferred as number of parameters.

Chapter 4
Implementation of engine model

Project 1C.
This chapter presents how different models for sub components of the engine are implemented in Simulink.

4.1 Accelerator pedal

Implementation:
In Figure 4.1 the implementation of accelerator pedal model in Simulink is shown.


PIC

Figure 4.1: Simulink model of the implemented throttle.


4.2 Throttle

Implementation:

4.3 Intake manifold

Implementation:

4.4 Fuel injector

Implementation:

4.5 Cylinder

Implementation:

4.6 Exhaust system

Implementation:

4.7 Lambda sensor

Implementation:

4.8 Complete engine model

Implementation:
Figure ???? shows the entire engine model implemented in Simulink which is in agreement with the initial design in Figure 2.3.

Chapter 5
Controller

Project 1C, continued.
A fuel controller is connected to the engine fuel injection system. In this chapter, the design of the controller, the way we went about it to determine the controller parameters, and figures showing that we have succeeded are presented. A step in accelerator position is implemented, and the plausibility of the results is assessed.

Implementation:

Chapter 6
Experiments

Project 1C, continued.
This chapter first presents the results of a number of experiments specified in the project compendium. The results are discussed and the plausibility is assessed.

6.1 Acceleration test

6.2 Test run on part of drive cycle

Include figures for: Continuous Lambda, Engine speed, Engine torque, Vehicle Speed VS drive cycle speed

6.3 Simulation of complete drive cycle

Include figures for: Continuous Lambda, Engine speed, Engine torque, Vehicle Speed VS drive cycle speed

6.4 Fuel consumption

6.5 Emissions

Chapter 7
Conclusions

Project 1C, continued.