Integrated Energy and Catalyst Thermal Management for Plug-In Hybrid Electric Vehicles

: With plug-in hybrid electric vehicles (PHEVs), the catalyst temperature is below the light-off temperature due to reduced engine load, extended engine off period, and frequent engine on/off shifting. The conversion efﬁciency of a three-way catalyst (TWC) and tailpipe emissions were proven to depend heavily on the temperature of the catalyst. The existing energy management strategy (EMS) of the PHEVs focuses on the improvement of fuel efﬁciency and emissions based on hot engine characteristics, but neglects the effect of catalyst temperature on tailpipe emissions. This paper presents a new EMS that incorporates a catalyst thermal management method. First, an additional cost is established to implement additional constraints on catalyst temperature, and then the global cost function is created using this additional cost and the fuel consumption. Second, we ﬁnd the global optimal solution using Pontryagin’s minimum principle method, which provides an optimal control policy and state trajectories. Then, based on the analysis of the optimal control policy, an engine on/off ﬁlter (eng on/off ﬁlter) is introduced to command the engine on/off shifting. This ﬁlter plays an important role in adjusting both the energy and catalyst thermal management strategy for PHEVs. Finally, a practical approach based on the eng on/off ﬁlter is developed, and a genetic algorithm is applied to optimize the time constants of this ﬁlter. Simulation results demonstrate that the proposed approach‘s fuel consumption increased slightly, but the tailpipe emissions of HC (hydrocarbons), CO (carbon monoxide) and NOx (nitrogen oxide) signiﬁcantly decreased compared with the standard approach. after eight repeated NEDC cycles, and increased by 2.27% under the UDDS–HWFET–UDDS driving cycle, but the tailpipe emissions of HC, CO, and NOx under the PMP-based EMS with catalyst temperature decreased by 17.9, 25.9, and 24.5%, respectively, under eight repeated NEDC cycles and decreased by 16.65, 24.15, and 23.88%, respectively, under the UDDS–HWFET–UDDS driving cycle. after


Introduction
An energy management strategy (EMS) is a crucial technology for plug-in hybrid electric vehicles (PHEVs) owing to its impact on fuel economy and emissions performance [1,2]. For PHEVs, the most direct and easiest EMS is the charge depletion-charge sustaining (CD-CS)-based energy management strategy [3] that applies a threshold to the state of charge (SOC) to control the mode change. Before the SOC reaches the predetermined threshold, the vehicle is mainly powered by the battery, which is a process called CD mode. After the SOC reaches the threshold, the vehicle is powered by the engine and the battery together, which is a process called CS mode [4,5]. As this strategy is not a blended strategy, it makes the charge deplete in the whole trip. CD-CS is not able to fully exploit the potential of the plug-in hybrid system [6]. Thus, many studies have been performed to improve the energy management of PHEVs. A variety of optimal methods, such as deterministic dynamic programming (DDP) [7][8][9][10], a two-scale dynamic programming (DP) approach [11,12], stochastic

Thermal Model of Engine and Three-Way Catalyst (TWC)
In this paper, the model for the supervisory control approach can be divided into two sub-models: engine thermal model and TWC thermal model. The steady-state engine model only outputs hot engine data, so establishing the engine thermal model that takes the engine's coolant temperature into account to predict cold engine outputs is necessary, especially for engine-out emissions during cold-start. The TWC thermal model includes catalyst temperature dynamics, which are required for computing the conversion efficiency of the TWC.

Thermal Model of Engine and Three-Way Catalyst (TWC)
In this paper, the model for the supervisory control approach can be divided into two sub-models: engine thermal model and TWC thermal model. The steady-state engine model only outputs hot engine data, so establishing the engine thermal model that takes the engine's coolant temperature into account to predict cold engine outputs is necessary, especially for engine-out emissions during cold-start. The TWC thermal model includes catalyst temperature dynamics, which are required for computing the conversion efficiency of the TWC.

Engine Thermal Model
The engine thermal model can be further divided into two sub-blocks, coolant temperature dynamics, and cold-engine correction factor. A lumped-capacitance thermal network model, depicted in Figure 2, was defined for the coolant temperature dynamics model. The heat transfer equations for calculating the coolant temperature are listed in Appendix A.

Engine Thermal Model
The engine thermal model can be further divided into two sub-blocks, coolant temperature dynamics, and cold-engine correction factor. A lumped-capacitance thermal network model, depicted in Figure 2, was defined for the coolant temperature dynamics model. The heat transfer equations for calculating the coolant temperature are listed in Appendix A. For control purposes, a simplified model that predicts cold engine outputs should be built. One approach is to simply multiply hot engine outputs by a cold factor, which is a function of the coolant temperature. This method was investigated by Murrell et al. [31], which can be expressed as follows: (1 )  For control purposes, a simplified model that predicts cold engine outputs should be built. One approach is to simply multiply hot engine outputs by a cold factor, which is a function of the coolant temperature. This method was investigated by Murrell et al. [31], which can be expressed as follows: where κ is the engine's coolant temperature factor, T set is the engine cooling system's thermostat set point, and T cool is the coolant temperature. These temperatures are expressed in degrees centigrade. m no x _hot are the hot engine fuel consumption rate and hot engine outputs rate of HC, CO, and NOx, respectively. These values can be obtained by looking up maps by engine speed and engine torque, shown in Figure 3. From Figure 3a, it can be seen that fuel enriches at high speed or torque. This can be explained by looking at the fundamentals of the internal combustion engine. For the high speed or torque, the only way to increase power is to richen the mixture after the wide-open-throttle. Thus, after the wide-open-throttle, the mixture becomes richer and richer, the air/fuel ratio gets smaller and smaller and fuel efficiency is lower and lower as the required power increases. m no x _cold are the cold engine fuel consumption rate and cold engine outputs rate of HC, CO, and NOx, respectively; and φ 1 , γ 1 , φ 2 , γ 2 , φ 3 , γ 3 , φ 4 , and γ 4 are the curve-fitting parameters. The cold-start test data of fuel consumption and emissions are not available to the author, so above parameters come from Advanced Vehicle Simulator (ADVISOR). Their values are 1, 3.1, 7.4, 3.072, 9.4, 3.21, 0.6 and 7.3, respectively.

TWC Thermal Model
The schematic diagram of the TWC is shown in Figure 4. The TWC can be simplified as the catalyst monolith (A), catalyst internal shell (B), and catalyst external shell (C). The TWC thermal model uses the exhaust gas flow rate ( flow w ) and the exhaust gas temperature ( in T ) at the TWC inlet as its inputs. The value of these inputs can be obtained by looking up maps indexed vertically by engine speed and horizontally by engine torque. The exhaust gas flow rate map, shown in Figure 5a, is available from the fuel map and air-fuel ratio (A/F) according to the following equation: The map of the exhaust gas temperature is shown in Figure 5b, which can also be obtained from the fuel map through Equations (3)- (5): where waste_pwr is the engine's waste heat, m fuel is the mass of fuel from fuel map, hv l is the lower heating value of the fuel, eng w is the speed of the engine with the unit of rad/s, M eng is the engine torque, ex_pwr is the power of the exhaust gas, and f_rac is the fraction of waste heat that goes to the exhaust. This parameter can be estimated by engine speed, it is shown in Figure 5. ex_tmp is the temperature of the exhaust gas, flow m is the exhaust gas flow with unit g/s, flow c is the capacitance of the exhaust gas, and amb T is the ambient temperature. It is known from Figure 6b, the exhaust gas temperature decreases with the increasing engine torque in most regions. According to Equation (5), since flow c and amb T are constants, the exhaust gas temperature is determined by

TWC Thermal Model
The schematic diagram of the TWC is shown in Figure 4. The TWC can be simplified as the catalyst monolith (A), catalyst internal shell (B), and catalyst external shell (C). The TWC thermal model uses the exhaust gas flow rate (w f low ) and the exhaust gas temperature (T in ) at the TWC inlet as its inputs. The value of these inputs can be obtained by looking up maps indexed vertically by engine speed and horizontally by engine torque. The exhaust gas flow rate map, shown in Figure 5a, is available from the fuel map and air-fuel ratio (A/F) according to the following equation: The map of the exhaust gas temperature is shown in Figure 5b, which can also be obtained from the fuel map through Equations (3)- (5): where waste_pwr is the engine's waste heat, m f uel is the mass of fuel from fuel map, l hv is the lower heating value of the fuel, w eng is the speed of the engine with the unit of rad/s, M eng is the engine torque, ex_pwr is the power of the exhaust gas, and f_rac is the fraction of waste heat that goes to the exhaust. This parameter can be estimated by engine speed, it is shown in Figure 5. ex_tmp is the temperature of the exhaust gas, m flow is the exhaust gas flow with unit g/s, c flow is the capacitance of the exhaust gas, and T amb is the ambient temperature. It is known from Figure 6b, the exhaust gas temperature decreases with the increasing engine torque in most regions. According to Equation (5), since c flow and T amb are constants, the exhaust gas temperature is determined by the power of the exhaust gas and the exhaust gas flow. To explain the above phenomenon in Figure 6b, we take 3000 rpm as an example. The change curves of the exhaust gas power, flow and temperature with engine torque when engine speed is 3000 rpm are shown in Figure 7. According to this figure, when the engine torque is below 100 Nm, the exhaust gas flow's increasing rate is greater than the rate of the exhaust gas power; therefore, the exhaust gas temperature decreases with the increasing torque in this torque range. When the engine torque is between 100 Nm and 140 Nm, the line segment AB is nearly parallel to CD, and BE is nearly parallel to DF, these mean that the exhaust gas flow's increasing rate is nearly the same as the rate of the exhaust gas power; thus, the exhaust gas temperature is almost unchanged during this torque range. When the engine torque is above 140 Nm, the exhaust gas power's increasing rate is obvious greater than the rate of the exhaust gas flow, so the exhaust gas temperature increases with the increasing engine torque during this torque range.
Energies 2018, 11, x FOR PEER REVIEW 6 of 30 the power of the exhaust gas and the exhaust gas flow. To explain the above phenomenon in Figure 6b, we take 3000 rpm as an example. The change curves of the exhaust gas power, flow and temperature with engine torque when engine speed is 3000 rpm are shown in Figure 7. According to this figure, when the engine torque is below 100 Nm, the exhaust gas flow's increasing rate is greater than the rate of the exhaust gas power; therefore, the exhaust gas temperature decreases with the increasing torque in this torque range. When the engine torque is between 100 Nm and 140 Nm, the line segment AB is nearly parallel to CD, and BE is nearly parallel to DF, these mean that the exhaust gas flow's increasing rate is nearly the same as the rate of the exhaust gas power; thus, the exhaust gas temperature is almost unchanged during this torque range. When the engine torque is above 140 Nm, the exhaust gas power's increasing rate is obvious greater than the rate of the exhaust gas flow, so the exhaust gas temperature increases with the increasing engine torque during this torque range.     the power of the exhaust gas and the exhaust gas flow. To explain the above phenomenon in Figure 6b, we take 3000 rpm as an example. The change curves of the exhaust gas power, flow and temperature with engine torque when engine speed is 3000 rpm are shown in Figure 7. According to this figure, when the engine torque is below 100 Nm, the exhaust gas flow's increasing rate is greater than the rate of the exhaust gas power; therefore, the exhaust gas temperature decreases with the increasing torque in this torque range. When the engine torque is between 100 Nm and 140 Nm, the line segment AB is nearly parallel to CD, and BE is nearly parallel to DF, these mean that the exhaust gas flow's increasing rate is nearly the same as the rate of the exhaust gas power; thus, the exhaust gas temperature is almost unchanged during this torque range. When the engine torque is above 140 Nm, the exhaust gas power's increasing rate is obvious greater than the rate of the exhaust gas flow, so the exhaust gas temperature increases with the increasing engine torque during this torque range.     the power of the exhaust gas and the exhaust gas flow. To explain the above phenomenon in Figure 6b, we take 3000 rpm as an example. The change curves of the exhaust gas power, flow and temperature with engine torque when engine speed is 3000 rpm are shown in Figure 7. According to this figure, when the engine torque is below 100 Nm, the exhaust gas flow's increasing rate is greater than the rate of the exhaust gas power; therefore, the exhaust gas temperature decreases with the increasing torque in this torque range. When the engine torque is between 100 Nm and 140 Nm, the line segment AB is nearly parallel to CD, and BE is nearly parallel to DF, these mean that the exhaust gas flow's increasing rate is nearly the same as the rate of the exhaust gas power; thus, the exhaust gas temperature is almost unchanged during this torque range. When the engine torque is above 140 Nm, the exhaust gas power's increasing rate is obvious greater than the rate of the exhaust gas flow, so the exhaust gas temperature increases with the increasing engine torque during this torque range.      Neglecting the exhaust gas heat loss to the exhaust manifold and catalyst inlet/outlet pipes, the thermal network representing the catalyst and associated thermal elements is shown in Figure 8. The TWC was modeled through a three-node lumped capacitance model, which includes monolith ( A T ), internal shell ( B T ) and external shell ( C T ). Heat is exchanged from the exhaust gas to the node of the monolith and the node of the internal shell through convection. A part of the monolith's heat transmits to the internal shell via conduction. Heat is exchanged from the node of the internal shell to the node of the external shell via conduction and radiation. Heat transfers from the node of the external shell to the ambient air through convection and radiation.
where k is thermal conductivity, A is the corresponding surface area, x is the representative distance between nodes, h is the convective heat transfer coefficient, h T and l T are the surface temperatures, ε is the emissivity, and σ is the Steffan-Boltzman constant. Neglecting the exhaust gas heat loss to the exhaust manifold and catalyst inlet/outlet pipes, the thermal network representing the catalyst and associated thermal elements is shown in Figure 8. The TWC was modeled through a three-node lumped capacitance model, which includes monolith (T A ), internal shell (T B ) and external shell (T C ). Heat is exchanged from the exhaust gas to the node of the monolith and the node of the internal shell through convection. A part of the monolith's heat transmits to the internal shell via conduction. Heat is exchanged from the node of the internal shell to the node of the external shell via conduction and radiation. Heat transfers from the node of the external shell to the ambient air through convection and radiation. R cn1 and R cn2 are the resistance of thermal conduction, R cv1 and R cv2 are the resistance of convective heat transfer, and R rd1 and R rd2 are the resistance of radiative heat transfer. These thermal resistances can be obtained by the following equations: where k is thermal conductivity, A is the corresponding surface area, x is the representative distance between nodes, h is the convective heat transfer coefficient, T h and T l are the surface temperatures, ε is the emissivity, and σ is the Steffan-Boltzman constant. Then, the equation for the lumped capacitor model is described as: where m A is the mass of the catalyst monolith (ceramic); C A is the lumped thermal capacitance of the catalyst monolith T A is the temperature of the catalyst monolith; h g2A is the convective heat transfer coefficient between exhaust gas and catalyst monolith, which is a function of exhaust gas flow, and the function, from ADVISOR, is expressed as follows.
A A is the inner (honeycomb) surface area of the catalyst monolith; T B is the temperature of the catalyst internal shell; . Q A is the change rate of the catalyst monolith's thermal energy; . Q g2A is the net heat flow from the exhaust gas to the catalyst monolith through convection; . Q A2B is the net heat flow from the catalyst monolith to the catalyst internal shell via conduction; . Q gen is the net heat flow from chemical reactions of the exhaust gas; and X hc , X co , and X no x are the flow rates of HC, CO, and NOx in the exhaust gas, these parameters can be calculated with Equation (1). According to Equation (1), X hc , X co and X no x , also expressed as . m hc_cold , . m co_cold , and . m no x _cold , are the function of engine torque, engine speed, and the engine's coolant temperature, while the engine's coolant temperature is the function of the engine torque and engine speed. Therefore, these parameters are dependent on engine torque and engine speed. M hc , M co , and M no x are the molar masses of HC, CO, and NOx, respectively; ∆H hc , ∆H co , and ∆H no x are the heat production of HC, CO, and NOx, respectively, by chemical reaction. η hc , η co , and η no x are the conversion efficiencies of HC, CO, and NOx, respectively. Catalyst conversion efficiencies are the function of catalyst temperature. Additionally, there is a catalyst efficiency adjustment (decrease), which is made at high exhaust flow rates. The relationship between conversion efficiency and the catalyst temperature can be described by arctan functions. Energies 2018, 11, x FOR PEER REVIEW 8 of 30 Figure 8. Thermal network representing the catalyst and associated thermal elements.
Then, the equation for the lumped capacitor model is described as: where A m is the mass of the catalyst monolith (ceramic); A C is the lumped thermal capacitance of the catalyst monolith A T is the temperature of the catalyst monolith; 2 g A h is the convective heat transfer coefficient between exhaust gas and catalyst monolith, which is a function of exhaust gas flow, and the function, from ADVISOR, is expressed as follows.
A is the inner (honeycomb) surface area of the catalyst monolith; B T is the temperature of the catalyst internal shell; A Q  is the change rate of the catalyst monolith's thermal energy; is the net heat flow from the exhaust gas to the catalyst monolith through convection; Q  A2B is the net heat flow from the catalyst monolith to the catalyst internal shell via conduction; gen Q  is the net heat flow from chemical reactions of the exhaust gas; and hc X , co X , and x X no are the flow rates of HC, CO, where i represents each type of emission; T lighto f f _i is the light off temperature of each type of emissions; K i is a tuning parameter that represent a slope of the efficiency function; ω exh_i is a correction factor for the exhaust gas flow rate, this factor can be approximated by a linear function as shown below [32].
where X i is flow rate of each emission; K i , a i and b i are the parameters estimated by experimental data. Finally, the list of parameters is shown in Table 2. The list of equations used for the catalyst internal shell (B) is described as follows: where m B is the mass of the catalyst internal shell; C B is the lumped thermal capacitance of the catalyst internal shell, h g2B is the convective heat transfer coefficient between the exhaust gas and the catalyst internal shell; A B is the surface area of the catalyst internal shell; T C is the temperature of the catalyst external shell; . Q B is the change rate of the catalyst internal shell's thermal energy; . Q g2B is the net heat flow from the exhaust gas to the catalyst internal shell through convection; . Q B2C is the net heat flow from the catalyst internal shell to the external shell via conduction and radiation; ε is the emissivity; and σ is the Stefan-Boltzmann constant m 2 .
The list of equations used for catalyst external shell (C) is as follows: where m C is the mass of the catalyst external shell; C C is the lumped thermal capacitance of the catalyst external shell, and the value is 460 J/kgK; T am is the ambient temperature; . Q C is the change rate of the catalyst's external shell's thermal energy; . Q C2am is the net heat flow from the catalyst external shell to the ambient air via convection and radiation; and A Cr is the surface area of the catalyst's external shell.
After converting the above equations, the state equation of the catalyst monolith's temperature was obtained: where M e (t) is the engine torque and n(t) is the engine speed.

Parameter Estimation and Model Validation
We focus on parameter estimation and model validation of the engine and TWC thermal model in this section. Some of these thermal models' parameters, shown in Table 3, are from ADVISOR, which was developed by the American National Renewable Energy Laboratory (NREL) for rapid analysis of the performance and fuel economy of conventional, electric, and hybrid vehicles and contains a component data file (its directory is ":\ADVISOR 2002\date"), this file contains several sub files, such as fuel converter data files, exhaust after treatment files, transmission data files, and driving cycle files. We can obtain the parameters in Table 3 by inspecting the fuel converter data files and exhaust after treatment files. The remaining parameters, more sensitive to coolant and catalyst temperature, come from calibration. These remaining parameters are calibrated by comparing the model's temperature response with the test data. The parameters of the engine thermal model are tuned to match the coolant temperature responses of the model to those of the road test data, after the engine thermal model is tuned properly, parameters of the TWC model are then tuned to generate the catalyst temperature responses that match with those of road test data. Note that only limited the coolant and catalyst temperatures test data are available to the author due to the difficulty of experiment set-ups for the engine out emissions and tail out emissions, and thus, the tail-pipe emission of the model response with those of test data is not performed. The process of the real vehicle on road test and parameter estimation are shown as follows. Firstly, we chose a route in our university campus for the plug-in hybrid electric vehicle's road testing. Then, the vehicle starts by the ISG motor. After starting, the operating mode of the vehicle was switched from pure electric mode to only engine driving mode. The signals of the vehicle speed, the engine speed, the engine torque, the coolant temperature and the catalyst temperature are collected during the test. The vehicle speed is shown in Figure 9a, and the engine torque and motor torque are shown in Figure 9b. Finally, using the vehicle speed, the engine speed, and the engine torque test data as inputs to the thermal model of engine and TWC, the parameters of the engine thermal model are tuned to match the coolant temperature responses of the model to those of the road test data, the comparison between model coolant temperature and test data temperature is shown in Figure 9c. The parameters of the TWC thermal model are tuned to match the catalyst temperature responses of the model to those of the road test data, the comparison between model catalyst temperature and test data temperature is shown in Figure 9d. The list of parameters, obtained from the tuning process, is shown in Table 4. The heat production of HC by chemical reaction 2306 (J/mol) ∆H co The heat production of CO by chemical reaction 410 (J/mol) ∆H nox The heat production of NOx by chemical reaction 417 (J/mol) Table 4. Thermal models' parameters determined from model tuning and validation.

Parameters Description Value
Engine thermal model The coefficient of conductance between engine cylinder and engine interior. 500 (W/K) The coefficient of conductance between engine interior and engine exterior. 500 (W/K)

K x2h
The coefficient of conductance between engine exterior and engine hood.
The product of the mass and the lumped thermal capacitance of the catalyst internal shell. 1232 (J/K) The product of the mass and the lumped thermal capacitance of the catalyst external shell. 1048 (J/K) The resistance of thermal conduction between the catalyst monolith and the internal shell.
The resistance of thermal conduction between the catalyst internal shell and external shell. Torque (Nm) Figure 9. Parameter estimation by matching the model's temperature response to the test data. Figure 9. Parameter estimation by matching the model's temperature response to the test data.
After parameter estimation, another road test date is used to validate above engine and TWC thermal model. This road test's vehicle speed is shown in Figure 10a, the engine torque and motor torque are shown in Figure 10b, the comparison between model temperature and test data temperature is shown in Figure 10c,d. From Figure 10c,d, the thermal model of the engine and TWC can predict the coolant and catalyst temperature well, despite some discrepancy existence.

Global Cost Function and Constraints
The energy management strategy based on optimal control theory seeks to minimize the global cost function over the total length of the trip. Most commonly, the global cost function is only the

Global Cost Function and Constraints
The energy management strategy based on optimal control theory seeks to minimize the global cost function over the total length of the trip. Most commonly, the global cost function is only the fuel consumption. However, considering some other costs in the global cost function it is possible to implement additional constraints, for example on emissions [33], drivability [34], or temperature [26]. In this paper, the proposed global cost function includes an additional cost on catalyst temperature evolution in addition to the fuel consumption, as shown in Equation (13). The constraints in this equation are the limitations of the motor, the engine, and the battery. The torque delivered by the engine has to be greater than its minimum torque and within its maximum torque, as is the case for the motor. The SOC window is limited to guarantee the performance and longevity of the battery. This study also required that the final SOC be near 0.25.
where t represents the time; . m f l_cold is the instantaneous fuel consumption that can be calculated by Equation (1); t f is the terminal time; u is the control variable, selecting the output torque of the motor as the control variable; P req is the power required at the inlet of the CVT; M req is the torque required at the inlet of the CVT;v veh is the vehicle speed; T cat is the catalyst temperature; M mmin is motor's minimum torque; M mmax is motor's maximum torque; M emin and M emax are the minimum and maximum torque of the engine, respectively; SOC min and SOC max are the lower and upper limits of SOC, respectively; Φ SOC t f is a function ensuring a solution meeting the final requirement on the battery's SOC; and k(T cat ) is a weighting factor, which creates a trade-off between fuel consumption and the TWC's conversion efficiency. The value of this factor is shown in Figure 11. The factor is set to 0 when the catalyst temperature is between 350 • C and 800 • C because this temperature range is a suitable temperature range for the catalyst, and the conversion efficiency of the catalyst is high at this temperature range. When the catalyst temperature exceeds 800 • C, the value of the factor increases rapidly with increasing catalyst temperature, to prevent the catalyst temperature from rising further, as an excessively catalyst temperature will cause the catalyst to sinter and quicken the catalyst aging. When the catalyst temperature is below 350 • C, the value of the factor is negative, then the additional cost is also negative, which favors catalyst warming.
Energies 2018, 11, x FOR PEER REVIEW 14 of 30 evolution in addition to the fuel consumption, as shown in Equation (13). The constraints in this equation are the limitations of the motor, the engine, and the battery. The torque delivered by the engine has to be greater than its minimum torque and within its maximum torque, as is the case for the motor. The SOC window is limited to guarantee the performance and longevity of the battery. This study also required that the final SOC be near 0.25.
where t represents the time; fl_cold m  is the instantaneous fuel consumption that can be calculated by Equation (1) Figure 11. The factor is set to 0 when the catalyst temperature is between 350 °C and 800 °C because this temperature range is a suitable temperature range for the catalyst, and the conversion efficiency of the catalyst is high at this temperature range. When the catalyst temperature exceeds 800 °C, the value of the factor increases rapidly with increasing catalyst temperature, to prevent the catalyst temperature from rising further, as an excessively catalyst temperature will cause the catalyst to sinter and quicken the catalyst aging. When the catalyst temperature is below 350 °C, the value of the factor is negative, then the additional cost is also negative, which favors catalyst warming.

Establishment of the Hamilton Function
To minimize the global cost function, Pontryagin's minimum principle (PMP) was used, and the Hamiltonian function is:

Establishment of the Hamilton Function
To minimize the global cost function, Pontryagin's minimum principle (PMP) was used, and the Hamiltonian function is: where λ(t) and ρ(t) are the Lagrange factor, and SOC and T cat are the state variables, u i is the candidates of the control variable, u i = M req or u i ∈ [u min (t) : ∆u : u max (t)], i = 1, 2, 3 . . . n + 1. The simulation step size is 1 s. At t moment, T cat (u i , t). Then, Equation (15) is expressed as follows.
The optimal control u * is obtained, when the following conditions are satisfied.
Therefore, implementing a penalty for the catalyst temperature is actually implementing a penalty for the derivative of the catalyst temperature.
The state equation of the battery's SOC is as follows: S . OC(t) = − I(SOC,P m (t)) where V oc is the open circuit voltage of the battery, R int is the battery's internal resistance, Q 0 is the rated capacity of the battery, and P m is the output power of the battery. The dynamic of the co-state on SOC is defined as: .
Ignoring the impact of the battery SOC on the internal resistance and open-circuit voltage of the battery, then . λ(t) = 0 [35][36][37]. Therefore, λ(t) = λ 0 . The dynamic of the second co-state on catalyst temperature is defined as: .
According to Figure 11, ∂T cat is equal to 2.56 × 10 −3 , 0, or 10 −2 , whereas T cat (u) ≤ 5.3 × 10 −5 ≈ 0. From Equation (8), the Map of h g2A can be obtained, and it is shown in Figure 12. According to Figure 12, the maximum value of the 2 g A h is 33.67.
Equation (18) can be translated to the following equation.
According to Equation (22), we know that the Hamiltonian function is made up of three terms: the fuel consumption, equivalent fuel consumption of battery and an additional cost on catalyst temperature. Since the main tasks of the energy management system is coordination between the first two terms. The value of the additional term should be less than any remaining terms.   According to Figure 12, the maximum value of the h g2A is 33.67. Equation (18) can be translated to the following equation.
According to Equation (22), we know that the Hamiltonian function is made up of three terms: the fuel consumption, equivalent fuel consumption of battery and an additional cost on catalyst temperature. Since the main tasks of the energy management system is coordination between the first two terms. The value of the additional term should be less than any remaining terms.

(k(T cat ) + ρ(t))
. Figure 13 shows the change curve of the s . oc(u) and .
T cat (u) with respect to motor torque. It can be seen that the difference between s . oc(u) and . T cat (u) is 4 orders of magnitude. As a consequence, if the inequality (23) is satisfied, the difference between λ(t) and k(T cat ) + ρ(t) is at least four orders of magnitude. λ(t) is obtained when PMP with the only one Lagrange factor λ is simulated, and the value of λ(t) is about 10 3 , Therefore, k(T cat ) + ρ(t) is about 10 −1 . According to Figure 12, the maximum value of the 2 g A h is 33.67.
Equation (18) can be translated to the following equation.
According to Equation (22), we know that the Hamiltonian function is made up of three terms: the fuel consumption, equivalent fuel consumption of battery and an additional cost on catalyst temperature. Since the main tasks of the energy management system is coordination between the first two terms. The value of the additional term should be less than any remaining terms.     T cat (u) with respect to motor torque.
As a consequence, Thus, .

Simulation Results and Discussion
The calculation process of the integrated strategy is shown in Figure 14, where N is the speed of the ISG motor, M req is the torque required at the inlet of the CVT, u(t) is the control variable, in this paper, the control variable is the motor output torque M m , u min (t) is the lower limit of the control variable at t moment, and u max (t) is the upper limit of the control variable at t moment. The operating principle of the proposed strategy is as follows: T cat (t, u i ) and S . OC(t, u i ) are calculated for every candidate control variable u i . The next step is the calculation of the Hamiltonian functions. The last step is comparing and obtaining the optimal command that corresponds to the smallest Hamiltonian function.
In fact, the procedure, shown in Figure 14, is known as the calibration method. It starts from an initial guess for λ 0 and ρ 0 . The solution of the problem is then obtained by replacing at each time the value of M m (motor torque) resulting from the minimization Hamiltonian function. If the final value of the SOC does not match the desired terminal condition SOC t f − 0.25 ≤ 0.01, and the engine's catalyst has not been lighted off in all segments except for the starting stage, the value of λ 0 and ρ 0 are adjusted until above conditions on the state is met. This approach does not yield strictly optimal results, but is more representative of what could be implemented in an actual vehicle.
To prove the significance of integrating the catalyst temperature on the PMP-based EMS of PHEVs, simulations of the integrated strategy (PMP-based with catalyst temperature) and energy-only management strategy (PMP-based without catalyst temperature) were conducted under eight repeated New European Driving Cycle (NEDC) and an urban dynamometer driving schedule/highway/urban dynamometer driving schedule (UDDS-HWFET-UDDS) driving cycle, which represents a typical work-home commute, it starts in a suburban area, characterized by environmental protection agency (EPA) urban dynamometer driving schedule (UDDS), then continues on a highway (HWFET), and finally arrives to downtown urban area, UDDS. The engine output torque and catalyst temperature trajectory are depicted in Figures 15 and 16, respectively. The working points of the engine are shown in Figures 17 and 18, and the fuel consumption and tailpipe emissions information is provided in Table 5.  Figure 15 demonstrates the simulation results of two EMSs, PMP-based without catalyst temperature and PMP-based with catalyst temperature, over eight repeated NEDC driving cycles. The engine output torque and catalyst temperature trajectories for the PMP-based EMS without catalyst temperature are illustrated in Figure 15a. Under this strategy, the engine on/off switch is activated frequently. The light-off temperature of the catalyst is 300 °C, which is obtained from the  The engine output torque and catalyst temperature trajectories for the PMP-based EMS without catalyst temperature are illustrated in Figure 15a. Under this strategy, the engine on/off switch is activated frequently. The light-off temperature of the catalyst is 300 • C, which is obtained from the catalyst's product introduction. Then, as shown in Figure 15a, the catalyst has not been lighted off in the a1-a8 segments, so the tailpipe emissions will be bad. The engine output torque and catalyst temperature trajectories for the PMP-based EMS with catalyst temperature are illustrated in Figure 15b. Compared with the results of the PMP-based strategy without catalyst temperature, the number of times that the engine on/off switch is activated has been reduced significantly under the PMP-based strategy with catalyst temperature. After the engine starts, the output torque of the engine obviously increases under the condition of minimal engine efficiency loss. So, the engine's catalyst temperature increased rapidly, and the catalyst lighted off in all segments except for b1. catalyst's product introduction. Then, as shown in Figure 15a, the catalyst has not been lighted off in the a1-a8 segments, so the tailpipe emissions will be bad. The engine output torque and catalyst temperature trajectories for the PMP-based EMS with catalyst temperature are illustrated in Figure 15b. Compared with the results of the PMP-based strategy without catalyst temperature, the number of times that the engine on/off switch is activated has been reduced significantly under the PMP-based strategy with catalyst temperature. After the engine starts, the output torque of the engine obviously increases under the condition of minimal engine efficiency loss. So, the engine's catalyst temperature increased rapidly, and the catalyst lighted off in all segments except for b1.  Figure 16 demonstrates the simulation results of two EMSs, PMP-based without catalyst temperature and PMP-based with catalyst temperature, over a UDDS-HWFET-UDDS driving cycle. Figure 16a is the UDDS-HWFET-UDDS driving cycle. The engine output torque and catalyst temperature trajectories for the PMP-based EMS without catalyst temperature are illustrated in Figure 16b. These trajectories for the PMP-based EMS with catalyst temperature are illustrated in Figure 16c. Compared with the results of the PMP-based strategy without catalyst temperature, the number of times that the engine on/off switch is activated was reduced significantly under the PMP-based strategy with catalyst temperature. The catalyst under the EMS based on PMP without catalyst temperature did not light off most of the time, whereas the catalyst under EMS based on PMP with catalyst temperature lighted off in most segments except for the starting stages of c1 and c2.
The working points of the engine under eight repeated NEDC cycles and a UDDS-HWFET-UDDS driving cycle are shown in Figures 17 and 18, respectively. Most engine working points of EMSs based on PMP without catalyst temperature and PMP with catalyst temperature are very close to each other; only a few points of the EMS based on PMP with catalyst temperature diverged from the points of the EMS based on PMP without catalyst temperature, but these points still fell within the oil economic zone.  Figure 16 demonstrates the simulation results of two EMSs, PMP-based without catalyst temperature and PMP-based with catalyst temperature, over a UDDS-HWFET-UDDS driving cycle. Figure 16a is the UDDS-HWFET-UDDS driving cycle. The engine output torque and catalyst temperature trajectories for the PMP-based EMS without catalyst temperature are illustrated in Figure 16b. These trajectories for the PMP-based EMS with catalyst temperature are illustrated in Figure 16c. Compared with the results of the PMP-based strategy without catalyst temperature, the number of times that the engine on/off switch is activated was reduced significantly under the PMP-based strategy with catalyst temperature. The catalyst under the EMS based on PMP without catalyst temperature did not light off most of the time, whereas the catalyst under EMS based on PMP with catalyst temperature lighted off in most segments except for the starting stages of c1 and c2.
The working points of the engine under eight repeated NEDC cycles and a UDDS-HWFET-UDDS driving cycle are shown in Figures 17 and 18, respectively. Most engine working points of EMSs based on PMP without catalyst temperature and PMP with catalyst temperature are very close to each other; only a few points of the EMS based on PMP with catalyst temperature diverged from the points of the EMS based on PMP without catalyst temperature, but these points still fell within the oil economic zone.  The fuel consumption and tailpipe emissions of these EMSs are provided in Table 5. The tailpipe emissions are calculated by the following equation: Compared with the results of the PMP-based EMS without catalyst temperature, the fuel consumption under the PMP-based EMS with catalyst temperature increased by 2.2% after eight repeated NEDC cycles, and increased by 2.27% under the UDDS-HWFET-UDDS driving cycle, but the tailpipe emissions of HC, CO, and NOx under the PMP-based EMS with catalyst temperature decreased by 17.9, 25.9, and 24.5%, respectively, under eight repeated NEDC cycles and decreased by 16.65, 24.15, and 23.88%, respectively, under the UDDS-HWFET-UDDS driving cycle. Thus, after   The fuel consumption and tailpipe emissions of these EMSs are provided in Table 5. The tailpipe emissions are calculated by the following equation: Compared with the results of the PMP-based EMS without catalyst temperature, the fuel consumption under the PMP-based EMS with catalyst temperature increased by 2.2% after eight repeated NEDC cycles, and increased by 2.27% under the UDDS-HWFET-UDDS driving cycle, but the tailpipe emissions of HC, CO, and NOx under the PMP-based EMS with catalyst temperature decreased by 17.9, 25.9, and 24.5%, respectively, under eight repeated NEDC cycles and decreased by 16.65, 24.15, and 23.88%, respectively, under the UDDS-HWFET-UDDS driving cycle. Thus, after The fuel consumption and tailpipe emissions of these EMSs are provided in Table 5. The tailpipe emissions are calculated by the following equation: Compared with the results of the PMP-based EMS without catalyst temperature, the fuel consumption under the PMP-based EMS with catalyst temperature increased by 2.2% after eight repeated NEDC cycles, and increased by 2.27% under the UDDS-HWFET-UDDS driving cycle, but the tailpipe emissions of HC, CO, and NOx under the PMP-based EMS with catalyst temperature decreased by 17.9, 25.9, and 24.5%, respectively, under eight repeated NEDC cycles and decreased by 16.65, 24.15, and 23.88%, respectively, under the UDDS-HWFET-UDDS driving cycle. Thus, after integrating catalyst temperature in the PMP-based EMS of PHEVs, the vehicle's tailpipe emissions were considerably reduced with minimal fuel consumption increase. integrating catalyst temperature in the PMP-based EMS of PHEVs, the vehicle's tailpipe emissions were considerably reduced with minimal fuel consumption increase.

Energy Management Strategy (EMS) Based on Pontryagin's Minimal Principle (PMP) with Eng on/off Filter
Although the PMP-based approach could be implemented in the real-time control of a vehicle, the proposed PMP with catalyst temperature strategy requires an extra co-state after integrating catalyst temperature. As PMP-based EMSs are sensitive to the co-state, PMP-based EMSs with only one co-state are difficult to calibrate, let alone PMP-based EMSs with two co-states.
According to the results, as shown in Figures 15 and 16, the proposed PMP-based strategy with catalyst temperature influences the catalyst to light off by reducing the engine on/off shifting times. Based on the analysis of the optimal control policy, an engine on/off filter (eng on/off filter) was introduced to command the engine on/off shifting. The PMP-based strategy with an eng on/off filter can reduce the engine on/off shifting times without adding an extra co-state; therefore, it is easily applicable to a standard optimal online EMS.

Energy Management Strategy (EMS) Based on Pontryagin's Minimal Principle (PMP) with Eng on/off Filter
Although the PMP-based approach could be implemented in the real-time control of a vehicle, the proposed PMP with catalyst temperature strategy requires an extra co-state after integrating catalyst temperature. As PMP-based EMSs are sensitive to the co-state, PMP-based EMSs with only one co-state are difficult to calibrate, let alone PMP-based EMSs with two co-states.
According to the results, as shown in Figures 15 and 16, the proposed PMP-based strategy with catalyst temperature influences the catalyst to light off by reducing the engine on/off shifting times. Based on the analysis of the optimal control policy, an engine on/off filter (eng on/off filter) was introduced to command the engine on/off shifting. The PMP-based strategy with an eng on/off filter The calculation process of the PMP-based EMS with eng on/off filter is shown in Figure 19. The engine on/off filter, without adding any control variable or co-state, can prevent the engine from experiencing frequent starts and stops, and engine on/off requests sent to the filter are basically determined by the engine torque command. If a non-zero engine torque is distributed by the control Figure 19. Calculation process of PMP-based EMS with eng on/off filter.
The calculation process of the PMP-based EMS with eng on/off filter is shown in Figure 19. The engine on/off filter, without adding any control variable or co-state, can prevent the engine from experiencing frequent starts and stops, and engine on/off requests sent to the filter are basically determined by the engine torque command. If a non-zero engine torque is distributed by the control strategy, then an engine on request is prompted to the filter. Otherwise, an engine off request is sent to the filter. The filter has two time constants, t s1 s and t s2 s. The rules of the filter can be expressed as: (1) If the engine is on and the elapsed time of the zero engine torque signal is t s1 s, then the engine off request passes through the engine on/off filter, and the engine off command will be sent to the engine plant in the next time step. (2) If the engine is on and the elapsed time of the zero engine torque signal is shorter than t s1 s, then the engine off request cannot pass through the engine on/off filter, and the engine must remain in an idle state. (3) If the engine is off and the elapsed time of the non-zero engine torque signal is t s2 s, then the engine on request passes through the engine on/off filter, and the engine on command and engine torque command will be sent to the engine plant. (4) If the engine is off and the elapsed time of the non-zero engine torque signal is shorter than t s2 s, then, the engine on request cannot pass through the engine on/off filter. In this situation, the engine torque command cannot be implemented by the engine, which is still off. So, the original torque distribution will be redistributed, and the torque originally assigned to the engine will be transferred to the ISG motor. (5) Except for the situations mentioned above, the engine will remain in its current state.

Optimization of Filter Time Constants Based on Genetic Algorithm
The filter's time constants, t s1 s and t s2 s, considerably affect a vehicle's fuel economy and emissions. The fuel economy is worse because t s1 and t s2 change the original optimal torque distribution from PMP-based EMS, but if the values of t s1 and t s2 are properly optimized, the engine on/off shifting times decrease, the engine works concentrated, the catalyst will warm up rapidly, and vehicle emissions will decrease considerably at the expense of a minimal fuel increase. So, optimizing these time constants is necessary. The time constant optimization problem is highly non-linear, so finding the optimal constants through analytical or numeric methods is difficult. The genetic method was introduced to effectively determine these constants. The fitness function is expressed as: (27) where FC is engine's fuel consumption; HC, CO, and NOx are the tailpipe emissions of HC, CO, and NOx, respectively. FC max , HC max , CO max , and NOx max are the maximum values of the engine's fuel consumption and tailpipe emissions, and the fuel consumption and tailpipe emissions under the CD-CS control strategy for eight repeated NEDC driving cycles are defined as this maximum values; w 1 is the weight of the engine's fuel consumption; and w 2 , w 3 , and w 4 are the weight of tailpipe emissions of HC, CO, and NOx, respectively. In this paper, the parameters of t s1 and t s2 are integers, so we transform the floating point numbers of all the produced individuals to integers and then operate them on the PMP-based on EMS with eng on/off filter. The round function is used for rounding individuals to the nearest integer. The parameters of the algorithm were set as follows. The maximum iteration of the genetic method was 80, the population size was 100, the crossover probability was 0.7, and the mutation probability was 0.01. The optimization result is shown in Figure 20. With the continuous evolution of the population, the fitness function value decreased, and this value converged to 2.5269, and its corresponding best individual was (t s1 ,t s2 ) = (5,2). Since the best individual is optimized for a given driving cycle, this individual may not suit for another velocity profile. In practical applications, we can optimize these parameters of t s1 and t s2 under different driving cycles, driving distances, and initial battery SOC off-line, then, a compromise value can be chosen as the final value.

Simulation Results of PMP-Based EMS with Eng on/off Filter
To validate the effect of the real-time EMS based on PMP with eng on/off filter, three EMSs were simulated under eight repeated NEDC driving cycles and a UDDS-HWFET-UDDS driving cycle. These EMSs included a PMP-based EMS without catalyst temperature, a PMP-based EMS with catalyst temperature and a PMP-based EMS with eng on/off filter. The simulation results are depicted in Figures 21 and 22, and the fuel consumption and tailpipe emissions information are provided in Tables 6 and 7

Simulation Results of PMP-Based EMS with Eng on/off Filter
To validate the effect of the real-time EMS based on PMP with eng on/off filter, three EMSs were simulated under eight repeated NEDC driving cycles and a UDDS-HWFET-UDDS driving cycle. These EMSs included a PMP-based EMS without catalyst temperature, a PMP-based EMS with catalyst temperature and a PMP-based EMS with eng on/off filter. The simulation results are depicted in Figures 21 and 22, and the fuel consumption and tailpipe emissions information are provided in Tables 6 and 7, respectively.

Simulation Results of PMP-Based EMS with Eng on/off Filter
To validate the effect of the real-time EMS based on PMP with eng on/off filter, three EMSs were simulated under eight repeated NEDC driving cycles and a UDDS-HWFET-UDDS driving cycle. These EMSs included a PMP-based EMS without catalyst temperature, a PMP-based EMS with catalyst temperature and a PMP-based EMS with eng on/off filter. The simulation results are depicted in Figures 21 and 22, and the fuel consumption and tailpipe emissions information are provided in Tables 6 and 7 Figure 21 outlines the simulation results under the three control strategies for eight repeated NEDC driving cycles. According to the figure, under PMP-based EMS without catalyst temperature, the engine starts and stops frequently, and the engine's working time is dispersed. The catalyst does not light off for most of the time, so the efficiency of the catalytic converter is low. under the PMP-based EMS with catalyst temperature and the PMP-based EMS with an eng on/off filter, the engine on/off shifting times decreased considerably, the engine's working time concentrated, the engine's catalyst temperature increased rapidly, and the catalyst lighted off in most segments except for the starting stage. Table 6 is the engine's fuel consumption and tailpipe emissions of the three EMSs under eight repeated NEDC driving cycles. Compared with the results of the PMP-based EMS with catalyst temperature, the fuel consumption, and HC, CO, and NOx emissions under the PMP-based EMS with an eng on/off filter increased by 0.8%, 5.33%, 5.45%, and 4.73%, respectively. However, compared with the results of the PMP-based EMS without catalyst temperature, the tailpipe emissions of HC, CO, and NOx under PMP-based EMS with an eng on/off filter decreased by 13.57%, 21.9%, and 20.99%, respectively, with the fuel consumption only increasing by 3%.   Figure 22a-c are the change curves of the engine torque and catalyst temperature under the different EMSs. According to the figure, the catalyst under the PMP-based EMS without catalyst temperature does not light off for the majority of the time, whereas the catalyst under the PMP-based EMS with catalyst temperature and the PMP-based EMS with an eng on/off filter lighted off in most segments except for the starting stage. This is because the engine under the PMP-based EMS without catalyst temperature starts and stops frequently, and the engine's working time is dispersed. However, under the PMP-based EMS with catalyst temperature and the PMP-based EMS with an eng on/off filter, the engine on/off shifting times obviously decreased, the engine's working time concentrated, and the engine's catalyst temperature increased rapidly. Table 7 provides the engine's fuel consumption and tailpipe emissions of the three EMSs under the UDDS-HWFET-UDDS driving cycle. Compared with the results of the PMP-based EMS with catalyst temperature, the fuel consumption, HC emissions, CO emissions and NOx emissions under the PMP-based EMS with an eng on/off filter increased by 1.62%, 6.06%, 6.26%, and 4.84%, respectively. However, compared with the results of the PMP-based EMS without catalyst temperature, the tailpipe emissions of HC, CO, and NOx under the PMP-based EMS with an eng on/off filter, decreased by 11.59%, 19.4%, and 20.19%, respectively, with the fuel consumption only increasing by 3.92%.
From above analysis, we concluded that the HC, CO, and NOx emissions of the proposed real-time PMP-based approach with eng on/off filter significantly decreased with a slight fuel consumption increase, compared with the standard PMP-based EMS without catalyst temperature.

Conclusions
The catalyst temperature has an important effect on the conversion efficiency of a three-way catalyst and vehicle tailpipe emissions. Therefore, the PMP-based EMS with catalyst temperature considered a soft constraint, also called an additional cost, on the catalyst temperature, and added the additional cost to the global cost function to be minimized. The problem was solved using Pontryagin's minimum principle, Our EMS significantly reduced the number of engine starts and stops times, enabled the rapid warmup of the catalyst, and greatly reduced the tailpipe emissions with an only slight fuel consumption increase. However, this EMS adds a co-state after integrating the catalyst temperature, which restricts its application in real-time control.
To solve the above problem, we introduced an engine on/off filter to command the engine on/off shifting. Based on the eng on/off filter, a real-time PMP-based EMS with eng on/off filter was developed, and a genetic algorithm was applied to optimize the time constants of the filter. To validate the effect of the real-time PMP-based EMS with an eng on/off filter, we simulated three EMSs under eight repeated NEDC driving cycles and a UDDS-HWFET-UDDS driving cycle. Simulation results demonstrated that the fuel consumption of our proposed approach slightly increased, but the tailpipe emissions of HC, CO, and NOx significantly decreased compared to the standard PMP-based approach without catalyst temperature.
Currently, the proposed strategy has only been verified through simulations. The next step would be to perform hardware-in-the-loop tests or experimental validations.
Author Contributions: Y.Z. wrote the paper and provided algorithms; Y.C. and C.C. built the simulation model and completed the simulation for different EMSs; and G.K. and W.G. analyzed the simulation results.

Conclusions
The catalyst temperature has an important effect on the conversion efficiency of a three-way catalyst and vehicle tailpipe emissions. Therefore, the PMP-based EMS with catalyst temperature considered a soft constraint, also called an additional cost, on the catalyst temperature, and added the additional cost to the global cost function to be minimized. The problem was solved using Pontryagin's minimum principle, Our EMS significantly reduced the number of engine starts and stops times, enabled the rapid warmup of the catalyst, and greatly reduced the tailpipe emissions with an only slight fuel consumption increase. However, this EMS adds a co-state after integrating the catalyst temperature, which restricts its application in real-time control.
To solve the above problem, we introduced an engine on/off filter to command the engine on/off shifting. Based on the eng on/off filter, a real-time PMP-based EMS with eng on/off filter was developed, and a genetic algorithm was applied to optimize the time constants of the filter. To validate the effect of the real-time PMP-based EMS with an eng on/off filter, we simulated three EMSs under eight repeated NEDC driving cycles and a UDDS-HWFET-UDDS driving cycle. Simulation results demonstrated that the fuel consumption of our proposed approach slightly increased, but the tailpipe emissions of HC, CO, and NOx significantly decreased compared to the standard PMP-based approach without catalyst temperature.
Currently, the proposed strategy has only been verified through simulations. The next step would be to perform hardware-in-the-loop tests or experimental validations.
Author Contributions: Y.Z. wrote the paper and provided algorithms; Y.C. and C.C. built the simulation model and completed the simulation for different EMSs; and G.K. and W.G. analyzed the simulation results.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: m f l_cold · l hv − M e · w e − . m exh c exh (T exh − T amb ) − Q f ric (A1)