Bruno Jeanneret
email: bruno.jeanneret@univ-eiffel.fr

Alice Guille

Des Buttes

Jérémy Pelluet
email: jeremy.pelluet@u-bourgogne.fr

Alan Keromnes
email: alan.keromnes@unilim.fr

Serge Pélissier
email: serge.pelissier@univ-eiffel.fr

Luis Le Moyne
email: luis.le-moyne_isat@u-bourgogne.fr

Optimal Control of a Spark Ignition Engine Including Cold Start Operations for Consumption/Emissions Compromises

Keywords: 0D engine model, consumption, pollutant emissions, optimal control problem, cold start

This study presents a semi-empirical modeling approach based on an extensive parametric study using a spark-ignition port-injection engine. The experimental results are used to derive engine-out emission models for each regulated pollutant (CO, HC, NO x ) as a function of engine operating parameters. Such parameters include engine speed, intake manifold pressure, equivalence ratio, and spark advance. The proposed models provide accurate predictions over a large range of engine operating conditions. The adequate accuracy and low computational burden of the models are promising in the context of optimal control theory. Dynamic programming is applied in order to find the best operating parameters that define trade-off between fuel consumption and emissions over driving cycles.

Introduction

Transportation is one of the major contributors to pollutant emissions and greenhouse gas emissions [START_REF]Emissions of Air Pollutants from Transport[END_REF][START_REF] Crippa | Fossil CO 2 and GHG Emissions of All World Countries: 2019 Report[END_REF]. Therefore, more stringent regulations have been introduced in order to reduce the negative impact of transportation systems [3]. Regulated pollutants include unburned hydrocarbons (HC), nitrogen oxides (NO x ), carbon monoxide (CO) and particulate matters. More and more efficient and expensive after-treatment systems have been introduced in order to reduce pollutant emissions from internal combustion engines such as the three way catalyst (TWC), the Diesel particulate filter (DPF) and the selective catalytic reduction (SCR). These devices are highly efficient but require specific operating conditions, as well as long and complex calibration procedures which induce a very high development cost of the internal combustion engine [START_REF] Castagné | Comparison of Engine Calibration Methods Based on Design of Experiments (DoE)[END_REF]. In the past, only new vehicles had to go through homologation procedures in a vehicle test bench over a predefined speed cycle with controlled environments. Now, Real Driving Emissions (RDE) measurements are also performed and this makes calibration procedures even more complex and expensive.

Although internal combustion engine (ICE) vehicles are expected to be replaced in the future by hybrid and zero emissions vehicles, optimization of fuel consumption and pollutant emissions is still a challenging issue. Coming regulations become highly stringent and make mandatory for next generation ICE vehicles to combine a low level of CO 2 and very low pollutant emissions. The same applies to hybrid vehicles where ICE will run in a large number of driving phases. This work is aiming to provide tools for effective and simultaneous reduction of fuel consumption and pollutant emissions for spark ignition ICE during complete drive cycles. It takes into account aspects related to catalyst start-up and heating phases, and explores how engine parameters variations from standard steady-state calibration can lead, over a full cycle, to better consumption vs. pollution trades. To do so, we develop a spark ignition engine-out emission model to be included in a dynamic vehicle model. Since three way catalyst efficiency highly depends on temperature [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF], we include a thermal model of the TWC [START_REF] Pandey | A simplified thermal model for the three way catalytic converter[END_REF] to calculate both temperature and efficiency. This is mandatory to estimate tailpipe emissions. We adopt an optimal control theory in order to find the best operating parameters of the engine. Dynamic programming is explored [START_REF] Scordia | Global Approach for Hybrid Vehicle Optimal Control[END_REF][START_REF] Vinot | Model simulation, validation and case study of the 2004 THS of Toyota Prius[END_REF]. It is a robust optimization method even if it is quite computationally intensive when the number of states and decision variables increases.

To optimize fuel consumption while respecting the standard emissions thresholds, models of pollutants have been developed. This topic is widely reported in the literature. Some are based on chemical kinetics mechanisms such as the extended Zeldovitch mechanism [START_REF] Lavoie | Experimental and Theoretical Study of Nitric Oxide Formation in Internal Combustion Engines[END_REF] estimating the production of nitric oxides (NO x ) emissions. CO and HC emissions mainly depend on the physical development of the flame front during the combustion process. The accurate prediction of pollutant emissions requires the coupling of 3D CFD model with chemical kinetics mechanisms in order to fully describe their production inside the combustion chamber on a crank-angle basis [START_REF] Kwon | Combustion and emission modelling of a direct-injection spark-ignition engine by combining flamelet models for premixed and diffusion flames[END_REF]. Such models lead to very high computation costs, unacceptable with a multi disciplinary approach such as vehicle design. Even 0D/1D models reduce the complexity of the problem but still have a computing time that cannot fit with optimization process over an integrated approach [START_REF] Karvountzis-Kontakiotis | Improvement of NO and CO predictions for a homogeneous combustion SI engine using a novel emissions model[END_REF].

Other models are entirely driven by physical assumption and experimental data. They allow defining mean value model for control oriented modeling, i.e., they model input/output behavior of the system with a reasonable precision but low computation complexity [START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF]. The model proposed by Abida and Claude [START_REF] Abida | Spark ignition engines and pollution emission: New approaches in modelling and control[END_REF] in 1994 calculates pollutant raw emissions of a spark ignition engine and takes into account numerous parameters such as throttle angle, spark advance and fuel/air ratio. They minimize step by step a criteria taking into account fuel consumption and raw pollutant emissions over the first elementary profile of speed in the European driving cycle. Successive studies have then explored this topic in the literature, studying the influence of other variables, reducing the complexity of models or using techniques such as artificial neural networks to formulate accurate models from a limited number of experimental data [START_REF] Sayin | Performance and exhaust emissions of a gasoline engine using artificial neural network[END_REF][START_REF] Nikzadfar | An extended mean value model (EMVM) for control-oriented modeling of diesel engines transient performance and emissions[END_REF]. In 2012, Andrianov et al. [START_REF] Andrianov | A Physics-Based Integrated Model of a Spark Ignition Engine and a Three-Way Catalyst[END_REF] developed a model that took into account the same variables as Abida's, but introduced the valve timing (synchronization of cylinder opening and closing with piston movement) and a three way catalyst model. The integrated model comprised 65 state variables and many non linearities. A reduction of the model made it possible to neglect some of these variables depending on the pollutant considered. According to this work, at least 7 state variables were mandatory in order to simulate exhaust pollutants with a good precision. However, this is far from being compatible with optimal control algorithm. Consequently, another approach has to be considered.

The main disadvantage of empirical hypothesis is that the models obtained are valid for a specific engine. Indeed, they are difficult to adapt to other engines because they do not take into account specificity such as combustion chamber parameters (geometry, surface/volume ratio, turbulence intensity). It is, therefore, necessary to calibrate the models for the engine under study [START_REF] Shayler | Generic functions for fuel consumption and engine-out emissions of HC, CO and NOx of spark-ignition engines[END_REF].

The control parameters selected for our study are intake manifold pressure, fuel/air ratio, spark advance and engine speed. Our simulation framework will include the following "main" parts:

•

A simple backward model of the drivetrain (i.e., from drive cycle to engine) as described by [START_REF] Horrein | Forward and Backward simulations of a power propulsion system[END_REF] to calculate at each time step the engine speed and torque.

•

A Mean Value Engine Model inspired from a Willans description [START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF] and adapted from Asus et al. [START_REF] Asus | Simple method of estimating consumption of internal combustion engine for hybrid application[END_REF] to include the effect of control parameters over overall efficiency.

•

Simple forms to describe the behavior of emissions with a semi empirical model based on experimental data.

• A catalytic converter model describing its thermal and efficiency behaviors with only one spatial discretization as proposed by Pandey et al. [START_REF] Pandey | A simplified thermal model for the three way catalytic converter[END_REF].

•

A dynamic programming algorithm to calculate the optimal control parameters of the engine over an entire drive cycle by minimizing a weighted sum of pollutant emissions and fuel consumption. This implementation was derived from [START_REF] Guille Des Buttes | Energy management strategy to reduce pollutant emissions during the catalyst light-off of parallel hybrid vehicles[END_REF].

The main objective of this work is to explore trade-offs between fuel consumption and pollutant emissions while adopting a systemic approach. By that, we mean considering the entire vehicle as a whole system, starting from a drive cycle up to the after treatment system, and taking into account drive train, engine and so on. The models are described in Section 2, followed by the results in Section 3. The originality of this approach will be further discussed in Section 4.

Models and Methods

Experimental Setup and Procedure

The DRIVE laboratory's test bench (see Figure 1) is equipped with the Renault D4F port fuel injection (PFI) engine (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014). It allows measuring the exhaust gas temperature at the engine outlet with a type K thermocouple. A Fourier Transform Infrared spectrometer (FTIR) collects gas samples before the catalyst to measure the concentrations of carbon monoxide (CO), unburned hydrocarbons (HC) and nitric oxide (NO) at the engine outlet. A balance of consumption deduces the fuel mass flow rate used to feed the engine. The control parameters close to the engine (spark timing, air/fuel ratio of the mixture and the position of the throttle valve controlling the intake pressure) are adjustable by means of the open engine control unit (ECU). The measurements were performed over a wide range of operating conditions as described in Table 1. The minimum spark advance value (i.e., spark retard compared to optimal which optimizes engine torque) corresponds to a temperature of 1000 K which should not be exceeded to avoid damage to the systems. The maximum value of spark advance is 25 degrees before the optimum advance. Depending on the chosen intake pressure (which determines the maximum pressure and temperature in the cylinder), some advances cannot be achieved because they may cause knocking. This phenomenon corresponds to a spontaneous-ignition of a part of the end-gas ahead of the flame. This ignition rapidly releases a large energy and causes high-frequency pressure oscillations which may damage the engine [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF]. Almost 500 points were recorded in this campaign.

Each operating point was stabilized during at least one minute before measuring. The data presented in this paper are available in the supplementary materials (see https://www. mdpi.com/2076-3417/11/3/971/s1). The fuel is oxidized during combustion producing mainly CO, particularly abundant around the flame front [START_REF] Bowman | Kinetics of pollutant formation and destruction in combustion[END_REF]. This species is then itself oxidized in a second step to give carbon dioxide if the temperature and the quantity of oxygen are sufficient. Andrianov expresses CO concentration only as a function of the equivalence ratio that controls the amount of oxygen present in the cylinder [START_REF] Andrianov | A Physics-Based Integrated Model of a Spark Ignition Engine and a Three-Way Catalyst[END_REF]. This approach is consistent with the experimental results obtained as shown in Figure 2. Indeed, the curves representing the series where the speed, spark advance and pressure vary, are very close to horizontal straight lines. Whatever the value chosen for the variable on the abscissa, CO concentration is constant. Figure 3 is a correlogram between intake pressure, CO concentration, fuel/air ratio φ, engine speed N and Spark Advance SA. It allows analyzing the relationship between each pair of numeric variables in a dataset. We notice that neither spark advance, nor engine speed or intake pressure show high correlation factor with CO concentration. Andrianov [START_REF] Andrianov | A Physics-Based Integrated Model of a Spark Ignition Engine and a Three-Way Catalyst[END_REF] proposes to use a polynomial of degree 5 to model the dependence of the CO concentration on fuel/air ratio. Abida and Claude [START_REF] Abida | Spark ignition engines and pollution emission: New approaches in modelling and control[END_REF] use the root of a polynomial of degree 2, which gives a suitable shape by reducing the number of coefficients to be adjusted. We have therefore chosen to use the form proposed by Abida and Claude (Equation ( 1)) and to adjust the coefficients with the least squares method. In Equation (1), φ represents the fuel/air equivalence ratio, co 1 and co 2 are the adjusting parameters. Symbols used in the document are described in the nomenclature section and values of the coefficients are listed in Table A1. [CO] = co 1 φ -1 + (φ -1) 2 + co 2 (1)

Figure 4 shows that the proposed fit is in good agreement with the measurements. A large number of observations are concentrated around the origin of the graph which corresponds to low values of CO concentration. Indeed, a major part of the data have a fuel/air ratio lower than 1 (lean fuel mixture), so there is enough oxygen to convert CO into CO 2 . The coefficient of determination is 0.99 and the mean error between the data and the predictions is 31%.

HC Model

Unburned hydrocarbons in gasoline engines mainly come from 4 mechanisms: flame quenching at the cylinder walls, entrapment of the unburned mixture in dead volumes (crevices), adsorption-desorption of the fuel in the oil film that covers the walls, and incomplete combustion or misfire [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF].

Equivalence ratio has an influence on the emissions of unburned hydrocarbons. Indeed, if the mixture is very low in fuel, the flame speed is reduced which can lead to misfire and higher cyclic variability leading to significant increase in HC concentration. With a rich mixture, some of the fuel remains after combustion when all the oxygen has been used to produce carbon monoxide or dioxide and water. Furthermore, increasing φ results in the raise of fuel vapor concentration in the bulk gas which will increase the absorption and hence the desorption from the oil layer [START_REF] Lakshminarayanan | Hydrocarbon Emissions from Spark Ignition Engines[END_REF]. Therefore the concentration of HC increases in both rich and lean mixes. This corresponds well to the behavior observed in Figure 5. Thus, we decided to introduce an augmented variable which is modeled by a second degree function around a minimum φ value of 0.9. The spark advance changes the HC emissions due to the storage of fresh gases in the crevices volumes and in the oil film on the walls [START_REF] Zhao | The effects of EGR and ignition timing on emissions of GDI engine[END_REF]. This is because trapped hydrocarbons can be released and oxidized during the expansion phase if the temperature and amount of oxygen are sufficient. When the spark advance is high, combustion ends earlier, which does not allow efficient oxidation and increases HC concentrations. Figure 5 shows this behavior, which we have modeled by a linear function.

The speed has an influence on the internal aerodynamics of the engine. Increasing the speed improves the homogeneity of the mixture and promotes turbulence during combustion. The flame speed increases which improves combustion efficiency and reduces HC emissions [START_REF] Kaiser | Engine-Out Emissions from a Direct-Injection Spark-Ignition (DISI) Engine[END_REF] . Figure 5 shows this effect, which we modeled using a linear function.

When the mixture is lean, the fuel concentration is lower and the average free path of the molecules increases. On the other hand, when the intake pressure decreases, the fraction of residual gas in the cylinder increases which also contributes to dilute the fresh gas. Simultaneously, a decrease in the wall temperature is observed [START_REF] Lakshminarayanan | Hydrocarbon Emissions from Spark Ignition Engines[END_REF]. Low intake pressure can thus lead to incomplete combustion and increased HC emissions, especially in combination with lean burning. We have chosen to model this influence with a linear function.

As shown in Figure 6, HC emissions are positively correlated with Spark Advance SA and the augmented variable (φ -0.9) 2 , and negatively correlated with intake pressure P intake and engine rotation speed N. The complete HC model expresses the causalities with the variables where hc1 to hc4 are the adjusting parameters. It is a linear model saturated to avoid negative concentrations that are not physically correct. Its expression is given by the Equation ( 2). [HC] = max(0,

hc 1 • N + hc 2 • (φ -0.9) 2 + hc 3 • SA + hc 4 • P intake + hc 5 ) (2)
The values for the different parameters are listed in Table A1. Figure 7 shows that the proposed adjustment correlates with the observations. The coefficient of determination is 0.74 and the mean error between the measurements and the estimate is 27%. About ten points have measured HC emissions higher than 840 ppm and are slightly off-centre with respect to the fit line between the model and the data. These points correspond to low equivalence ratio (φ < 0.8) and low or intermediate intake pressures (P intake ≤ 0.5 bar) except for one point (φ = 0.65, P intake = 1 bar).

NO Model

The oxides of nitrogen NO, N 2 O and NO 2 are all present in the exhaust of a spark ignition engine. However, nitrous oxide is not significant as it represents less than 1% of NO emitted. Nitrogen dioxide can be more present with the maximum occuring at wide open throttle and a rich mixture. Nevertheless, its concentration is less than a few percent of NO [START_REF] Raine | Modeling of nitric oxide formation in spark ignition engines with a multizone burned gas[END_REF]. This is confirmed by Bowman et al. [START_REF] Bowman | Kinetics of pollutant formation and destruction in combustion[END_REF], which states that: "NO is the predominant nitrogen oxide emitted by combustion devices". In this context, the present work will focus on NO species. It can be produced by three distinct mechanisms:

•

Early NO are formed at the flame front and in fuel-rich areas. They result from the combination of nitrogen N 2 with hydrocarbon radicals (C, C2, CH 2 or CH 2 ) to give the intermediate products HN, HCN, CN or CNH 2 . These can recombine with oxygen to give nitrogen oxide NO [START_REF] Bowman | Kinetics of pollutant formation and destruction in combustion[END_REF].

•

The NO-fuel are derived from the oxidation of the nitrogen atoms present in the fuel [START_REF] Fenimore | Formation of nitric oxide in premixed hydrocarbon flames[END_REF].

•
Thermal NO is formed by the direct reaction of nitrogen and oxygen present in the air when the temperature and oxygen concentration are sufficient. This process is an important contributor to the total NO x emissions and is known by the extended Zeldovich mechanism [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF][START_REF] Lavoie | Experimental and Theoretical Study of Nitric Oxide Formation in Internal Combustion Engines[END_REF][START_REF] Raimbault | Benefit of Air Intake Optimization for New Turbocharged Gasoline Engine[END_REF].

The rate of NO formation depends strongly on temperature as described by the Zeldovich mechanism. The temperatures necessary for NO formation are reached by the intense heat released during combustion (see Equation (3)).

d[NO] dt = 6 • 10 6 √ T • exp ( -69090 T ) • [O 2 ] 1/2 • [N 2 ] (3)
The heat release varies depending on the mass of fuel injected. In gasoline engines, this is mainly controlled by the intake pressure, which regulates the amount of mixture introduced into the cylinder. Figure 8 shows that NO emissions increase with intake pressure. This behavior will be modeled by a linear function.

The spark advance also influences the maximum temperature reached during combustion. Indeed, the timing of the ignition with the compression phase of the gases by the piston influences the temperature and pressure at the beginning of combustion and during the propagation of the flame [START_REF] Yang | Variation of Exhaust Gas Temperature with the Change of Spark Timing and Exhaust Valve Timing During Cold Start Operation of an SI Engine[END_REF]. Higher pressure increases the volume concentration of the mixture and reduces the average free path of the molecules. Combustion is faster and the heat release is more intense, which increases the maximum combustion temperature and NO emissions. Figure 8 thus shows that NO emissions increase with spark advance. This influence will be expressed by a linear function. Equivalence ratio φ also plays an important role in the formation of NO. Under slightly lean fuel conditions (0.9 ≤ φ ≤ 0.95), the fuel consumed allows the temperature of the cylinder to increase sufficiently and since oxygen is present in significant quantities, this favors the formation of nitrogen oxides. Leaner fuel conditions reduce the burned mass and therefore the maximum combustion temperature, which reduces NO emissions. Finally, a fuel-rich mixture reduces NO emissions because oxygen is consumed to produce CO or CO 2 via reactions that require less energy. The concentration of NO is therefore at a maximum and decreases in both rich and lean mixes. This corresponds well to the behavior observed in Figure 8 that we have modeled by a second degree function (φ -0.9) 2 .

The correlogram (see Figure 9) also show a good correlation with (φ -0.9) 2 . The engine speed does not have a strong influence on NO emissions as shown in Figure 8. We have therefore chosen not to include it in our model. Here again, we get a linear model with the augmented variable (φ -0.9) 2 , where the four adjusting parameters are no 1 , no 2 , no 3 and no 4 . Figure 10 shows that the proposed adjustment corresponds well to the measurements with a coefficient of determination of 0.78. To illustrate this part of our work, Figure 11 presents the concentration at (N = 2000 rpm, φ = 0.9 Bar) for various spark advance values. 14)).

[NO] = max([NO] min , no 1 • (φ -0.9) 2 + no 2 • SA + no 3 • P intake + no 4 ) (4)

Internal Combustion Engine Model

We need to make explicit the relationship between our control parameters and the engine model. This later comprises three parts: fuel consumption, brake torque and exhaust gas temperature. Fuel flow and torque models are mainly based on the literature, with forms inspired by those proposed by other authors. The temperature model rely on measurements made in the DRIVE laboratory with coefficients adjusted by the method of least squares.

Fuel Flow

The air flow rate is proportional to the speed and intake pressure. The fuel flow rate ṁ f uel is itself proportional to the air flow rate, with the exception that it depends to the chemical composition of the mixture.

To increase the accuracy of the model, the volumetric efficiency η vol is introduced. It expresses the real efficiency of the engine compared to a perfect volumetric pump and therefore its capacity to introduce fresh gases (Equation ( 5)).

ṁ f uel = φ φ + (m air /m f uel ) st • η vol • P intake × 10 5 • V d • N 60R c • M gaz R gaz • T adm (5)
The fill rate is highly dependent on the speed that controls the gas flow rate at the inlet. The flow is limited by flow friction, especially around the valves, which increases with speed (and fluid velocity). Moreover, at very high speed, a critical flow rate is reached and limits the flow around the valve during at least part of the intake phase (chocking phenomenon) [START_REF] Bordjane | Modélisation et Caractérisation Dynamique des Circuits d'admission et d'échappement des Moteurs à Combustion Interne[END_REF].

Conversely, filling at high speed can be improved by delaying the closing of the intake valve to take advantage of the inertia of the air column. This type of adjustment can result in intake backflow that degrades efficiency at low speed [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF]. Finally, filling can be maximized for certain engine speeds thanks to an acoustic tuning that increases intake pressure by wave propagation phenomena [START_REF] Bordjane | Modélisation et Caractérisation Dynamique des Circuits d'admission et d'échappement des Moteurs à Combustion Interne[END_REF].

The speed-dependent filling therefore has a maximum and decreases at both low and high speed. The position of the maximum depends on design choices and in particular on the geometry of the nozzles [START_REF] Raimbault | Benefit of Air Intake Optimization for New Turbocharged Gasoline Engine[END_REF][START_REF] Malkhede | Maximizing Volumetric Efficiency of IC Engine through Intake Manifold Tuning[END_REF]. According to the studies, the filling rate evolves in a range of about 0.1 to 0.3 points and the corresponding maximum speed varies from 1250 to 4000 rpm [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF][START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF][START_REF] Bordjane | Modélisation et Caractérisation Dynamique des Circuits d'admission et d'échappement des Moteurs à Combustion Interne[END_REF][START_REF] Malkhede | Maximizing Volumetric Efficiency of IC Engine through Intake Manifold Tuning[END_REF].

We have chosen to model the effect of engine speed on fill rate by a quadratic function which has an amplitude of 0.2 points and a maximum for an engine speed of 2250 rpm (Equation ( 6) and Figure 12).

The filling rate is also influenced by the equivalence ratio of the mixture. This is because the fuel is injected as a liquid and vaporizes in the intake pipes. The heat consumed lowers the temperature and increases the density of the gases and the filling rate of the cylinder [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF].

We chose to model this effect by varying the maximum value of the filling ratio (Equation [START_REF] Scordia | Global Approach for Hybrid Vehicle Optimal Control[END_REF] and Figure 12). This value is equal to 0.9 for a filling rate less than 1 (excess air) and increases for excess fuel (filling rate greater than 1). 2 (6)

η vol = η vol,φ -1.4 × 10 -8 • (N -2250)
η vol,φ = 0.9 φ ≤ 1 0.5 • φ + 0.4 1 < φ ≤ 1.2 (7)
Finally, the engine speed has an influence on the intake temperature [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF]. Indeed, at low speed, the residence time of gases in the intake system is more important, which favors thermal exchanges with the hot walls and increases their temperature [START_REF] Borman | Internal-combustion engine heat transfer[END_REF]. At high speed, the gases do not have time to warm up and remain at room temperature.

We have modeled this influence by a linear function of the speed, which has the largest temperature increase observed during the measurement campaign (16 °C) for the idle speed and decreases to room temperature for the maximum speed of 6000 rpm (Equation ( 8) and Figure 12).

T intake = T amb -0.003

• N + 18 (8)
We can see in Figure 13 that the model has a good fit to the measurements, with a determination coefficient of 0.96 and a mean error of less than 10%. Parameters values are listed in Table A3.

Brake Torque

By reversing the consumption model formulated by Asus et al. [START_REF] Asus | Simple method of estimating consumption of internal combustion engine for hybrid application[END_REF], we obtain a relationship between brake torque Γ ICE and fuel flow ṁ f uel (Equation ( 9)). In this expression, brake torque is proportional to the ratio between flow rate and engine speed N, with an extra term that estimates friction torque. The coefficient of proportionality corresponds to the engine efficiency, which depends on the parameters of close control of the engine (speed, equivalence ratio and spark advance), based on the shapes proposed by Guzzella and Onder [START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF].

The evolution of combustion efficiency with speed is taken from the model proposed by Asus et al. [START_REF] Asus | Simple method of estimating consumption of internal combustion engine for hybrid application[END_REF] (Equation [START_REF] Kwon | Combustion and emission modelling of a direct-injection spark-ignition engine by combining flamelet models for premixed and diffusion flames[END_REF] and Figure 14). The shape of the friction torque Γ f is based on a Heywood model [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF] (Equation ( 11)) .

Γ ICE = η f i • η comb • η φ • η SA • ṁ f uel • LHV πN/30 -Γ f ( 9
)
η comb = η c0 - A B + N (10)
Γ f = (97 + 0.015 • N + 5 × 10 -6 • N 2 ) • V d 2πR c (11)
Efficiency as a function of equivalence ratio decreases when the mixture is rich. This is because the extra fuel is not converted into energy due to a lack of oxygen. In addition, oxygen depletion can lead to incomplete combustion and reduce thermodynamic efficiency. On the contrary, in a lean mixture, all the fuel introduced can be consumed and the efficiency keeps a value of 1 over a wide range [START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF].

The model of Guzzella and Onder takes into account these two phenomena. They thus propose a curve proportional to the excess air, a transition then a step at 1. However, it is only applicable for equivalence ratio greater than 0.77. For fuel contents below this limit, combustion misfires may occur due to the low fuel concentration [START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF]. Performance deteriorates thus quickly until combustion becomes impossible (lean flammability limit) [START_REF] Liu | Hydrogen effect on lean flammability limits and burning characteristics of an isooctane-air mixture[END_REF].

We have chosen to model this behavior by a portion proportional to the excess air with a slightly larger directing coefficient than that of the rich mixture. The complete model is given by Equation ( 12) and plotted in Figure 14.

η φ =            1.373/φ -0.373 0.7 ≤ 1/φ < 0.95 0.9313 + (1 -0.9313) • sin 1/φ-0.95 1-0.95 0.95 ≤ 1/φ < 0.95 + π 2 (1 -0.95) 1 0.95 + π 2 (1 -0.95) ≤ 1/φ < 1.3 -1.4286/φ + 2.8571 1.3 ≤ 1/φ < 1.54 (12)
Spark advance has a significant impact on the conversion of thermal energy into mechanical energy. For a given speed, fuel/air ratio and fuel flow, there is an optimal advance that maximizes the torque and thus the engine's efficiency [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF][START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF][START_REF] Lumsden | Characterizing Engine Emissions with Spark Efficiency Curves[END_REF].

A relative spark advance (symbol ∆SA) can therefore be defined, which is the difference between the selected setting and the optimum spark advance for a given operating point. For the D4F engine used for this study, the optimal spark advance as a function of engine speed and intake pressure was determined during previous test campaigns (optimal values are reported in Table A2). We thus determine the relative advance for the various measurements carried out on the bench using the table of optimal advances (according to Equation ( 14)). The shape of the efficiency with regards to relative advance is inspired by Guzzella and Onder [START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF] (Equation ( 13) et Figure 14).

η SA = 1 -0.04 × 10 -2 • ∆SA 2 (13)
∆ SA = SA -SA opti (N, P intake ) (14)
We notice from Equation ( 14) that a negative value of ∆ SA brings the ignition closer to top dead center, and a positive value of ∆ SA increases the spark advance. We can see in Figure 15 that the model has a good fit to the measurements, with a determination coefficient of 0.92 and a mean error of 30%.

Exhaust Temperature

Last but not least exhaust temperature model has to be defined with regards to its action on catalytic converter behavior. Figure 16 represents the measurements. In the bottom left graphic, the exhaust temperature is described versus the relative spark advance. It should be noted that this was not the case for emissions plots presented above. The temperature rise at the engine outlet depends on the proportion of energy involved in combustion that remains in the exhaust and the thermal capacity of the gases. The proportion of energy still present in the exhaust depends on the engine's close control parameters. The exhaust temperature depends on the temperature rise and the intake temperature (Equation ( 15)). The latter is calculated using the same formula as for the flow model (Equation ( 8)).

T exh = T intake + k φ • k SA • k N • k Pintake • LHV C p,gaz • (1 + (m air /m f uel ) st φ ) (15)
Engine speed and intake pressure control the flow of gases into the engine. The higher the flow rate, the less time available for heat exchange with the walls relative to the quantities involved [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF][START_REF] Guzzella | Introduction to Modeling and Control of Internal Combustion Engine Systems[END_REF]. Heat losses at the walls decrease and the exhaust temperature increases. This behavior is apparent in Figure 16, and we have chosen to model the influence of these parameters with two linear functions (Equations ( 16) and ( 17)).

k N = a N + b N • N (16)
k Pintake = a Pintake + b Pintake • P intake [START_REF] Shayler | Generic functions for fuel consumption and engine-out emissions of HC, CO and NOx of spark-ignition engines[END_REF] The influence of equivalence ratio on exhaust temperature is similar to its effect on maximum combustion temperature (see Section 2.2.3) [START_REF] Heywood | Internal Combustion Engine Fundamentals[END_REF]. Both excess fuel and excess air form an inert mass that lowers the temperature in the cylinder and at the exhaust. The latter has therefore a maximum around the stoichiometric air/fuel ratio as shown in Figure 16. As the temperature model already depends on the equivalence ratio by the multiplicative factor in the denominator, we have added a linear function of the equivalence ratio to the numerator (Equation ( 18)) which allows obtaiingn a parabolic shape for the global influence of this parameter.

k φ = a φ + b φ • φ (18)
Finally, the later the ignition occurs (negative relative spark advance, i.e., spark retard), the more the combustion shifts towards the opening of the exhaust valve and the more the temperature increases, as shown in Figure 16. We have chosen to model this effect with a parabola (Equation ( 19)).

k SA = a SA + b SA • (∆ SA -∆ SA,min ) 2 (19)
The 9 model fitting parameters were determined using the method of least squares (2 constants for speed, pressure and fuel/air ratio, and 3 for spark advance). The Figure 17 shows the shapes of the coefficients which express the portion of residual energy at the exhaust. The global contribution of the wealth (in orange) corresponds to the (scaled) ratio between k φ and the coefficient depending on φ in the denominator. The model corresponds well to the measurements as can be seen in Figure 18, with a determination coefficient of 0.94 and a mean error of 5.6%.

Catalyst Model

The catalytic converter used in this study conforms to the EURO 6 emission standard. It is light and closely coupled to the engine for an accelerated light-off. The model (that is efficiency as a function of temperature and equivalence ratio, thermal behavior) has been widely described in reference [START_REF] Guille Des Buttes | Energy management strategy to reduce pollutant emissions during the catalyst light-off of parallel hybrid vehicles[END_REF] and the reader can refer to this article for further informations. Figure 19

Vehicle Dynamics

Vehicle dynamics and resulting power demand are modeled using a backward approach. The velocity as a function of time is known in advance and this driving cycle corresponds to a necessary torque at the wheels (Γ wheel ) to overcome the inertia of the vehicle as well as the resistive forces acting on it. The mechanical transmission components have constant efficiencies and the engine speed is saturated to its idle speed to model the clutch behavior.

Γ wheel = F res R tire + J veh R tire dv dt (20)
F res = m veh gk rol + 1 2 ρS F C x v 2 (21)
J veh = m veh R 2 tire + 4J wheel ( 22
)
Γ ICE = Γ wheel η e DF η e GB k DF k BV + J ICE ωICE (23)
e = -1 Γ wheel > 0 (traction) 1 Γ wheel < 0 (brake) ( 24)

ω ICE = max ω idle , v R tire k DF k GB (25)

Problem Formulation and Discretization Values

Dynamic programming is based on Bellman's principle of optimality. State variables are discretized temporally and on a scale of values. A graph whose vertices are defined by the value taken by the state variable at a given instant is thus constructed. The path from one vertex to another is associated with a value of the command and a cost from the point of view of the objective function. The path that minimizes the sum of the costs in the graph corresponds to the sequence of the optimal values of the command.

The associated optimal control problem is given by Equation ( 26), which is the weighted sum of pollutant emissions and fuel consumption.

J = ∑ cycle ṁ f uel ṁ f uel,re f δt + α 3 ṁCO ṁCO,re f + ṁHC ṁHC,re f + ṁNO ṁNO,re f δt ( 26
)
Both fuel consumption and pollutant emissions show extremely different orders of magnitude. To normalize them we consider the simulation of the reference vehicle for which mixture is stoichiometric and spark advance is optimal.

There are three control variables in the optimal control problem, intake pressure, relative spark advance and fuel/air equivalence ratio. By mixing Equations ( 5) and ( 9), we calculate the necessary intake pressure that satisfies Equation ( 23) and the required engine torque at each time step, thus reducing the problem complexity with two control variables only.

The state variable is the catalyst temperature as it drives the catalyst overall efficiency and the associated emitted pollutants. The problem is subject to its first-order dynamics expressed in discrete form. It can vary freely during the driving cycle, its initial temperature is ambient temperature and its final value is unconstrained.

The catalyst monolith temperature varies continuously despite the discrete control variables. To obtain a temperature mesh, it is rounded as expressed by Equation [START_REF] Raimbault | Benefit of Air Intake Optimization for New Turbocharged Gasoline Engine[END_REF]. The impact of the mesh precision will be discussed in the following section.

∆T cata = round δt m cata C p ( Qamb + Qgaz + Qreac ) • 1 δT δT ( 27
)
Table 2 summarizes the discretization steps for the state and the control parameters.

Results

The reference vehicle used in the present study is a conventional one whose characteristics are given in Table 3. The numerical parameters are given in Table A4. Figure 20 represents the vehicle speed and gearbox ratio during the NEDC cycle that was used to obtain the first results. The value of the weighting factor α between pollutant emissions and fuel consumption allows the definition of several scenarios. In the first part of this sec-tion, we present a pollution centered scenario which is compared to the reference strategy where nominal parameters are retrieved (i.e., equivalence ratio equals to 1 and optimal spark advance).

Pollution Centered Scenario

One can see the performance of the pollution centered scenario in Figure 21 (in this case, the weighting parameter, α, is equal to 5) compared to the reference strategy. Fuel consumption increases by 15%, but we observe a great reduction in pollutant emissions: 73% for CO, 57% for HC and 84% for NO. Numerical results are presented in Table 4. As it can be seen in Figures 22 and23, the catalyst light-off (instant where catalyst efficiency reaches 50% which corresponds to a temperature of 280 °C) is reached at 140 s. At this time, more than 80% of total CO have already been emitted, while almost 70% of total HC are produced. 50% of total tailpipe NO have already been dispersed in the environment. The behavior of the engine is far different with the pollution centered scenario (see Figure 24): the light-off of the catalyst is quicker, as it happens after only 60 s. Furthermore, a very small amount of pollutant is emitted during this period as illustrated in Figure 25: CO are reduced by a factor of 15, HC by 4, and NO x emissions before light-off are 11 times lower than in the nominal configuration. Figure 26 shows the positions of the actuators during the pollution centered scenario: We can distinguish between two different behaviors:

•

Catalyst light-off: during this period, we observe a late spark advance (10 to 15°c ompared to optimal ignition). This is combined with a lean mixture, Φ = 0.9. These two parameters increase the exhaust gas temperature. For constant engine torque, the reduced fuel conversion efficiency associated with late combustion requires greater fuel and air flow rate (see intake pressure in Figure 26). As a consequence, exhaust enthalpy and convective effects are high. This shortens the period where catalyst efficiency is low. • warm catalyst period: the strategy adopts a stoichiometric combustion which represents a good compromise to oxidize CO and HC while reducing NO x . In the same time, we observe a relative delay in spark advance to lower exhaust gas temperature and reduce NO (and HC to a lesser extent) formation. The emission of CO depends only on the fuel/air ratio (see Equation ( 1)) of the mixture, and stoechiometry is a good compromise for this pollutant. This delay increases with the brake power as it can be seen in the extra urban part of the NEDC cycle, again to lower NO emissions.

Trade-off between Fuel Consumption and Pollutant Emissions

In the previous section, we presented a pollution centered scenario where the weighting parameter α gives priority to pollutant emissions. In this section, we made several runs with different values of α. This balances the fuel economy compared to pollutant emissions. Results are presented in Figure 27. First, it can be noticed that a low weighting parameter (α = 0.1) of the objective function makes it possible to reduce pollutant emissions to a large extent without altering fuel consumption. In this case, the heating strategy of the catalyst is already implemented. With a hot catalyst, combustion is stoichiometric and a slight delay is observed compared to optimal spark advance (2°in this particular case).

When the value of α increases, NO emission can be further reduced by lowering exhaust gas temperature. As we can see in Figure 27, a good trade-off can be obtained by choosing a value of α between 1 and 2. Table 5 and Figure 29 show a big improvement obtained with an appropriate strategy (gains ranging from 50 to almost 80% depending on the pollutant) while maintaining a good performance in fuel economy (17% increase in fuel consumption). Of course, here again, several trade-off between fuel consumption and emissions can be adopted. First we notice in Figure 30 that both strategies present short light-off time: 75 s for the reference simulation, and only 19 s for the pollution centered scenario. These short times are the consequence of the more agressive dynamics of the drive cycle compared to NEDC. The strategy for the pollution centered scenario remains the same: significant deterioration of the spark advance (20 to 30°BTDC compared to optimal SA), and lean combustion (Φ = 0.9). Throttle is full open. This induces an important exhaust flow at a high temperature.

Worldwide Harmonized Light-Duty Vehicles Test Cycle

When the catalyst is initiated, here again the strategies are similar to what happened during NEDC cycle: stoichiometric combustion, and slight delay of spark ignition in order to lower CO

Discussion

Originality of the Work

There is a great number of publications about catalyst behavior to satisfy more and more stringent emissions targets. Some authors focus on the three way catalyst model ( [START_REF] Pandey | A simplified thermal model for the three way catalytic converter[END_REF][START_REF] Brandt | Three-way catalyst diagnostics for advanced emissions control systems[END_REF][START_REF] Sanketi | A simplified catalytic converter model for automotive coldstart control applications[END_REF][START_REF] Tomforde | A Post-Catalyst Control Strategy Based on Oxygen Storage Dynamics[END_REF]). Other authors work on minimizing one species during the cold start period ( [START_REF] Shaw | A simplified coldstart catalyst thermal model to reduce hydrocarbon emissions[END_REF][START_REF] Mozaffari | A Nonlinear Model Predictive Controller With Multiagent Online Optimizer for Automotive Cold-Start Hydrocarbon Emission Reduction[END_REF]). However, there is little literature that proposes trade-off between fuel consumption and pollutant emission during a whole drive cycle, without any empirical assumption of the duration of the cold start period nor arbitrary values for spark advance or lean/rich combustion.

Engine calibration is a critical step in the vehicle development process. Castagné et al. [START_REF] Castagné | Comparison of Engine Calibration Methods Based on Design of Experiments (DoE)[END_REF] make a comprehensive survey of different methods used to tune all calibration parameters. Every method is quite intensive and not completely satisfying. Most of the existing experimentally derived maps have not been optimized to be used for real-world applications. The model-based development appears to be a good process since it allows reducing the number of experiments with an intensive use of simulations. Developing engine and vehicle models that can be mathematically assisted with optimization methods can bring a solution to calibrate parameters. Moreover, rule-based maps can only represent the features of a specific engine, and hence, by a slight change in the ICE's configuration, tedious design-of-experiment efforts should be repeated to revise the previous cold-start engine input maps. In this paper, we present a methodology that can be applied to every vehicle configuration.

In this sense, the tool presented in this paper can participate to the calibration task in two ways:

• defining the maps to be used by real time control strategies for both hot and cold catalyst • serving as reference values for comparison with online methods

Discussion of the Method

An optimal control method allows calculating the best strategy, regarding the weighting parameter, with a steady-state formulation of the problem. This is an important fact since the optimization is made possible at the cost of simplifications such as neglecting certain state variables and their dynamics: we can cite the temperature rise of the cylinder head, which can induce high emissions of HC during a cold start period. We also neglect the wall wetting phenomena, typical in port fuel injection engines [START_REF] Yildiz | Spark ignition engine fuel-to-air ratio control: An adaptative control approach[END_REF][START_REF] Ye | Model predictive control on wall wetting effect using Markov Chain Monte Carlo[END_REF]. Certain portion of the fuel droplet would condense on the wall of the intake manifold, forming the fuel puddle, which has a significant impact on the mixture. Detrimental air to fuel ratio excursion could occur, giving rise to catalyst efficiency problems, particularly during transients.

Another important simplification is that all measurements have been made under steady-state operating conditions to feed the control model. However, any drive cycle includes transients during which steady-state engine behavior is not reached, leading to spikes in Lambda value. One can argue that part of these deviations can be partially mitigated by O 2 storage in the catalyst [START_REF] Muske | Estimating the oxygen storage level of a three way automotive catalyst[END_REF].

Challenges and Limitations

The driving cycle has to be known in advance to establish the optimal energy management strategy with dynamic programming, which cannot be implemented in real vehicle. Other work focuses on the prediction of real-life driving cycle in order to adapt dynamic programming solutions to real-time control strategies [START_REF] Joud | Smart Energy Management for Series Hybrid Electric Vehicles Based on Driver Habits Recognition and Prediction[END_REF]. Other approaches such as minimum Pontryaguin principle could be investigated and adapted in real time. However, one purpose of this algorithm is to ascertain the potential of reduction of pollutant emissions for conventional vehicles. It provides a reference of comparison to evaluate other sub optimal real-time strategies.

The generic functions developed in this paper have been derived from a specific data set from a particular engine. Nevertheless, for similar engine architectures, one can expect trends to be the same, although quantitatively there may be variations. The general optimization methodology and results presented here should therefore apply for SI engines with equivalent features. The units considered for pollutant models are concentrations, so an increase in the mass flow will naturally result in pollutants increase. Engine models introduce a displacement volume parameter, which can be used for scaling purposes. Next project will include the latest technologies for direct injection spark ignition engines (DISI) and corresponding control parameters in the optimization process (turbo charging, variable valve lifting, . . . ). The models will be evolving with these new technologies.

A new objective will consist of adapting the algorithm to hybrid vehicles. Indeed, one can wonder if the additional degree of freedom brought by hybridization would change the engine management. From this perspective, a new state variable that represents the battery state of charge has to be introduced in the algorithm. We get three dimensions for the states (time, state of charge and catalyst temperature) and three dimensions for the control variables (intake pressure, spark advance and equivalence ratio). This is very challenging regarding to the use of computer's memory and calculation time.

Conclusions

The original hypothesis of this work is that optimal control theory can provide information on how to drive the actuators of an internal combustion engine both during cold start and warmed catalyst. Indeed, the simulation results show that some trade-off can be found between fuel consumption and pollutant emissions. This is achieved in two ways: retarding the spark advance combined with a lean mixture to enhance catalyst ligth-off; keeping a slightly late ignition with a stoichiometric combustion as soon as catalyst is hot.

These findings are determined using two-dimensional dynamic programming with a weighted objective function. The entire model comprises eight sub-models and 35 tuning parameters. Each sub-model is based on semi-empirical assumptions. We take into account the thermal dynamics of the catalyst monolith, simple models for consumption, emissions and exhaust temperature. This decreases the computational burden while retaining the associated physical trends.

Figure 2 .

 2 Figure 2. CO concentration as a function of control parameters: (a) engine speed, (b) fuel/air equivalence ratio, (c) spark advance and (d) intake pressure. Each curve represents a series of measurements where only the parameter in abscissa changes.

Figure 3 .

 3 Figure 3. Correlogram of CO emission with governing parameters.

Figure 4 .

 4 Figure 4. Comparison between [CO] measurement and [CO] prediction based on [CO] model.

Figure 5 .

 5 Figure 5. HC concentration as a function of control parameters (each curve represents a series of measurements where only the parameter in abscissa changes).

Figure 6 .

 6 Figure 6. Correlogram of HC emission with governing parameters.

Figure 7 .

 7 Figure 7. Comparison between [HC] measurement and [HC] prediction based on [HC] model.

Figure 8 .

 8 Figure 8. NO concentration as a function of control parameters (each curve represents a series of measurements where only the parameter in abscissa changes).

Figure 9 .

 9 Figure 9. Correlogram of NO emission with governing parameters.

Figure 10 .

 10 Figure 10. Comparison between [NO] measurement and [NO] prediction based on [NO] model.

Figure 11 .

 11 Figure11. Emission model at (N = 2000 rpm, φ = 0.9 Bar)-relative spark advance ∆ SA varies from -20°to 20°( see Equation (14)).

Figure 12 .

 12 Figure 12. Influence of parameters on fuel flow model.

Figure 13 .

 13 Figure 13. Comparison between ṁ f uel measurement and ṁ f uel prediction based on ṁ f uel model.

Figure 14 .

 14 Figure 14. Influence of parameters on torque model.

Figure 15 .

 15 Figure 15. Comparison between Γ ICE measurement and Γ ICE prediction based on Γ ICE model.

Figure 16 .

 16 Figure 16. Exhaust temperature as a function of control parameters (each curve represents a series of measurements where only the parameter in abscissa changes).

Figure 17 .

 17 Figure 17. Influence of parameters on exhaust temperature model.

Figure 18 .

 18 Figure 18. Comparison between T exh measurement and T exh prediction based on T exh model.

Figure 19 .

 19 Figure 19. Catalyst temperature during WLTC cycle (model vs. measurements).

Figure 20 .

 20 Figure 20. NEDC driving cycle speed and gearbox ratio.

Figure 21 .

 21 Figure 21. Optimal results over NEDC cycle.

Figure 22 .Figure 23 .

 2223 Figure 22. Thermal behavior of the catalyst (reference strategy); red star: TWC light-off.

Figure 24 .

 24 Figure 24. Thermal behavior of the catalyst (pollution centered scenario); red star: TWC light-off.

Figure 25 .

 25 Figure 25. Cumulated emissions (pollution centered scenario); red star: TWC light-off.

Figure 26 .

 26 Figure 26. Optimal parameters to control the engine.

Figure 27 .

 27 Figure 27. Trade-off between fuel consumption and pollutant emissions (NEDC).

Figure 28 Figure 28 .

 2828 Figure28shows the vehicle speed and gearbox ratio during the WLTC that was used to validate the results on a more realistic and dynamic driving cycle.

Figure 30 .

 30 Figure 30. Cumulated emissions (WLTC); yellow star: TWC light-off for the reference strategy-purple star: TWC light-off for the pollution centered strategy.

∆T 1 ω 2 J 1 N

 121 cata variation of catalyst temperature over a step time in K ∆ SA,min adjusting parameter (effect of spark advance on exhaust temperature) ∆ SA relative spark advance in CAD BTDC δt discretization step for the time in s Qamb thermal flux of catalyst with ambient in W Qgaz thermal flux of catalyst with exhaust gas in W Qreac thermal flux of catalyst corresponding to chemical reaction in W η φ fuel/air eq ration's contribution to engine efficiency η c0 reference combustion efficiency η comb combustion efficiency η DF differential efficiency η f i indicated fuel efficiency η GB gear box efficiency η SA spark advance's contribution to engine efficiency η vol,φ Equivalence ratio's contribution to volumetric efficiency η vol volumetric efficiency of the engine Γ f friction torque in N•m Γ ICE brake torque of internal combustion engine in N•.m Γ wheel Wheel torque in N•m ω ICE engine speed in rad•s -idle idle speed in rad•s -1 φ fuel/air equivalence ratio ρ air density in kg•m -3 a φ adjusting parameter (effect of fuel/air eq ratio on exhaust temperature) a N adjusting parameter (effect of engine speed on exhaust temperature) a Pintake adjusting parameter (effect of intake pressure on exhaust temperature) a SA adjusting parameter (effect of spark advance on exhaust temperature) b φ adjusting parameter (effect of fuel/air eq ratio on exhaust temperature) b N adjusting parameter (effect of engine speed on exhaust temperature)) b Pintake adjusting parameter (effect of intake pressure on exhaust temperature) b SA adjusting parameter (effect of spark advance on exhaust temperature)C p,gaz thermal capacity of exhaust gas in J•K -1 •kg -1 C p thermal capacity of the catalyst in J•K -1 •kg -1 C x drag coefficient co X Adjusting parameters for the [CO] model F res rolling resistance force in N g acceleration of gravity in m•s -2 hc X adjusting parameters for the [HC] model J ICE Engine inertia in kg•m 2 J veh vehicle inertia in kg•m wheel Inertia of a wheel in kg•m 2 k DF differential gear ratio k GB gear box ratio k N effect of engine speed on exhaust temperature k phi effect of fuel/air eq. ratio on exhaust temperature k Pintake effect of intake pressure on exhaust temperaturet k rol rolling coefficient k SA effect of spark advance on exhaust temperature LHV low heating value of the fuel in W • s • g -1m cata catalyst mass in kg ṁ f ,re f average fuel mass flow of the reference vehicle running the NEDC test cycle in g•s -1 ṁ f uel fuel mass flow in g•s -1 m veh vehicle mass in kg ṁX,re f average mass flow of pollutant species X on the NEDC test cycle in g•s -1 ṁX mass flow of each pollutant X in g•s -engine speed in rpm no X adjusting parameters for the [NO] model P intake intake manifold pressure in bar R c number of revolutions per combustion cycle (2 for four stroke engine) R gaz perfect gas constant in J•mol -1 •K -1 R tire tire radius in m S F frontal area of the vehicle en m 2 SA spark Advance (CAD BTDC) SA opti optimal spark advance (maximum torque) in CAD BTDC T temperature of the gas in K t time in s T amb ambiant temperature in K T exh exhaust temperature in K T intake intake temperature in K v vehicle speed in m•s -1 V d displacement volume in liter (m air /m f uel ) st air to fuel ratio at stoichiometry Appendix A. Model Parameters Appendix A.1. Engine Parameters

Table 1 .

 1 Operating conditions.

							Intake manifold pressure (bar)			0.3 0.5 0.7 0.9 1
							Engine speed (rpm)			800 1000 1500 2000 3000 4000
								fuel/air eq. ratio (-)			0.7 0.8 0.85 0.9 0.95 1 1.1 1.2
							Spark Advance (°BTDC)						-10 : 5 : 50
						2.2. Emissions Model
						2.2.1. CO Model
							Figure 2 is a graphical illustration of the measurements of CO concentration. Each
						graph represents the CO concentration as a function of one control parameter, i.e., (a) engine
						speed, (b) fuel/air equivalence ratio, (c) spark advance and (d) intake pressure. Each curve
						represents a series of measurements where only the parameter in abscissa changes.
		1.5									8	10 4
	CO (ppm)	0.5 1				(a)				CO (ppm)	6 4					(b)
												2
			500 0 8	1000	1500	2000 Engine speed (rpm) 2500 3000	3500	4000	4500	0 0.6 8	10 4	0.7		0.8 Equivalence ratio 0.9	1 (-)	1.1	1.2
			6									6
	CO (ppm)	4				(c)				CO (ppm)	4					(d)
			2									2
			0									0
			-20	-10	0	10	20	30	40	50	60	0.2		0.3	0.4	0.5	0.6	0.7	0.8	0.9	1
						Spark advance (°BTDC)								Intake pressure (bar)

Table 2 .

 2 Discretization step size used in this paper.

	Variable	Steps	Min. Value	Max. Value
	catalyst temperature	0.1 K	not bounded
	Fuel/air eq. ratio	0.01	0.7	1.1
	Relative spark advance	2°-30°10°S
	tep time	1 s	NA

Table 3 .

 3 Main vehicle characteristics and types of model used in this study.

		Component		Size		Type of Model
	Internal combustion engine		1.8 L		Mean Value Engine Model
	Three way catalytic converter	EURO 6 compliant		0D model
	Chassis (Peugeot 308 SW-2009)	1470 kg		longitudinal forces
						5
						4
	Vehicle speed in km/h					2 3	Gearbox (-)
						1
						0
	200	400	600	800	1000	1200
		Time (s)

Table 4 .

 4 Performance of the pollution centered scenarios compared to the reference strategy (NEDC).

		Ref. Strategy SA opti φ = 1.0	α = 5 SA var φ var	Relative Variation
	Consumption (L/100 km)	6.57	7.59	+15%
	CO (g/km)	0.53	0.144	-73%
	HC (g/km)	0.056	0.024	-57%
	NO (g/km)	0.212	0.035	-84%

Table 5 .

 5 Performance of the pollution centered scenarios compared to the reference strategy (WLTC).

			Ref. Strategy SA opti φ = 1.0	α = 5 SA var φ var	Relative Variation
		Consumption (L/100 km)	6.86	8.02	+17%
		CO (g/km)	0.2554	0.128	-50%
		HC (g/km)	0.029	0.013	-55%
		NO (g/km)	0.179	0.040	-77%
		1.4
				Consumption
		1.2		CO HC
				NO
		1
	Performance	0.6 0.8
		0.4
		0.2
		0
		= 1 ; SA opt	Pollution-centered
		Figure 29. Optimal results over WLTC cycle.

  and NO concentrations.

						Cumulative emissions
		CO (g)	T >50 %
			200	400	600	800	1000	1200	1400	1600	1800
	HC (g)	0.5
			200	400	600	800	1000	1200	1400	1600	1800
		NO (g)	Reference strategy Pollution centered strategy
			200	400	600	800	1000	1200	1400	1600	1800
						Time (s)

Table A1 .

 A1 Parameters values for emission models.

	co 1	164.073	hc 1	-0.0866	no 1	4.666
	co 2	0.00217	hc 2	4976	no 2	-73.687
			hc 3	5.864	no 3	105.7
			hc 4	-0.1995	no 4	-2669
			hc 5	480.4	[NO] min	300

Table A2 .

 A2 Optimal values for spark advance as a function of intake pressure and engine speed.

									N (rpm)
			500	750	1000	1500	2000	2500	3000	3500	3750	4000	4500	5000	5500 6000
		200	29	30	33	36	38	40	40	40	44.1	40	40	40	40
		300	30	30	30	33	36	40.1	40.1	38.1	42.1	40.1	40	40	40
		400	28	28	30	32	33.3	35.1	35.1	34.1	34.1	34.1	34	34	34
	P intake (mBar)	500 600 700	22 20 19	25 23 20	28 25 20	30 27 21	31.1 28.1 23.1	33.1 29.1 27.1	31.1 28.1 28.1	29.1 27.1 26.1	32.1 29.1 26.1	31.1 30.1 26.1	31 29 26	31 29 26	31 29 26
		800	21	20	20	24	21.1	26.1	28.1	24.1	25.1	24.1	25	25	25
		900	20	22	20	22	21.1	24.1	24.1	24.1	24.1	24.1	24	24	24
		1000	20	20	20	17	19.1	21.1	21.1	20.1	23.1	24.1	21	21	21

Table A3 .

 A3 Parameters values for fuel flow, brake torque and exhaust temperature models.

	R c	2	V cyl	1.8
	η c0	0.98	η f i	0.43
	A	300	R gaz	8.314
	B	2000	LHV	44,000
	a SA	5.507 • 10 -5	(m air /m f uel ) st	14.7
	b SA C p,gaz	0.311 1.2 • 10 3	∆ SA,min a N	20 1.468 • 10 -4
	a φ	-0.771	b N	0.251
	b φ	1.503	a Pintake	0.00108
			b Pintake	0.805

Funding: This research received no external funding.

Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.

Data Availability Statement:

The data presented in this study are available in https://www.mdpi. com/2076-3417/11/3/971/s1.

Supplementary Materials:

The following are available online at https://www.mdpi.com/2076-341 7/11/3/971/s1, Data collection for emissions and engine models.


Carole Pelissier

Josephine Lemaitre

Hubert Marotte

Vincent Bonneterre

Barbara Charbotel

Frédéric Dutheil

Paul Vercherin

Luc Fontana

Factors associated with back complaints among French medical students

Keywords: Odd Ratio, CI: Confidence Interval, ORadj: adjusted Odd Ratio P-value NS : non significant * p-value<0.05, ** p-value<0.01, *

Factors associated with back complaints among French medical students

Carole PELISSIER (a,b,*) , Joséphine LEMAITRE (a) , Hubert MAROTTE (c) , Vincent BONNETERRE (d) , Barbara CHARBOTEL (e) , Frédéric DUTHEIL (f) , Paul VERCHERIN (g) Musculoskeletal complaints are now a growing health problem affecting both adults and young people (1) . Medical students are more likely to be affected by increased workload in wards of hospitals, stress and long study hours leading sedentary life, during the period of their medical training (2) . Besides they are increasingly using computers in teaching and learning (2) . These events are considered modifiable musculoskeletal complaints risk factors that may increase the prevalence of musculoskeletal complaints among medical students. The study sought to assess the prevalence of back complaints among French medical students and to explore the relationship between back complaints and occupational factors in order to provide prevention measures.

A descriptive, cross-sectional multicentre study was conducted between November 2016 and February 2017 in France among medical students in 4th or 5th years who have clinical training in hospital and medical teaching in faculty. A self-questionnaire comprising 25 questions, covering the medical and occupational factors was developed for this study. Musculoskeletal complaints were assessed on the Nordic Musculoskeletal Questionnaire, a valid instrument for epidemiological survey into musculoskeletal disorders (ache, pain, discomfort) (3) . A visual analogue scale (VAS), validated for stress assessment in clinical practice, was used (4) . A cut-off point of 7 was proposed for clinical discrimination between high stress and moderate-low stress [START_REF] Fx | Validity of occupational stress assessment using a visual analogue scale[END_REF] . Data collection and analysis were approved by the national data protection authority and by the local 'Ethics Committee' (IRBN642016/CHUSTE).

A total of 786 students (62.6% female) with an average age of 22 years (standard deviation 2.06 years) were enrolled in the study (participation rate of 54%). During the last seven days, the prevalence of neck, upper back and low back were respectively 47.7%, 33.5%, and 49.4%. During the last twelve months, the prevalence of neck, upper back and low back were respectively 59.0%,46.6%, and 70.1%.

Results from univariate analysis show that, among occupational factors, pace of work with clinical training half-day, a high level of global stress, the daily workload more than 8 hours per day, were associated with back complaints during the last seven days and during the last twelve months (Table 1). However, these associations do not remain significant after multivariate analysis.

We observed a high frequency of neck and lower back complaints and a significant association between back complaints and modifiable occupational factors, such as the pace of work and the level of global stress. Our findings regarding prevalence are consistent with those in the literature. According to Nordic classification medical students were, lifetime prevalence (73.4%), 1-year prevalence (46.1%-59.9%) of low back pain (6) . Active coping is considered a better way to deal with stress than avoidant coping strategies (7,8) . Shearer et al. suggest that a short mindfulness training program can be a powerful yet simple way to help students cope with stress and anxiety (9) .

In conclusion, to prevent back disorders among medical students, multidisciplinary actions that include changes in the work environmental factors along with stress management programs that teach people how to cope better with stressful events could be developed in faculties.

  Luc FONTANA (a,b) a : Service de Santé au travail CHU de Saint-Etienne b : Univ Lyon, Univ Lyon 1, Univ St Etienne, IFSTTAR, UMRESTTE, UMR_T9405, 42005 SAINT ETIENNE, France c : Service de Rhumatologie, CHU de Saint-Etienne SAINBIOSE INSERM U1059, University of Lyon, Saint-Etienne CIC 1408, University of Lyon, Saint-Etienne d : Equipe EPSP, laboratoire TIMC-IMAG, Université Joseph Fourier-UFR de Médecine, Domaine de la Merci, La Tronche, Cedex 38706, France. e : Univ Lyon, Univ Lyon 1, IFSTTAR, UMRESTTE, UMR_T9405, F-69373 LYON, France f: University Hospital of Clermont-Ferrand (CHU), Preventive and Occupational Medicine, Clermont-Ferrand, France. g: Santé Publique et Information Médicale, CHU de Saint-Etienne *Corresponding author: Carole Pélissier, Email adress : pelissier.carole.chu@gmail.com

Funding

The authors received no funding for this work.

Contributions of authors.

All authors have role in conducting the research and writing the manuscript. The final version of this manuscript is approved by all authors.

Disclosure of interest

The authors declare that they are no competing interest. 0.9 [START_REF] Lesage | Clinical stress assessment using a visual analogue scale[END_REF]