1. Introduction
Turbocharged engines are very popular today. Since the turbocharger generates extra power, it allows smaller turbocharged engines to equal larger non-turbocharged ones in acceleration (Chen, 2012; Bell, 1997; Holmbom, Liang, & Eriksson, 2017; Miller, 2008). But because the turbocharged engines are smaller, they are also generally capable of better fuel economy (Abdelmadjid, Mohamed, & Boussad, 2013; Bontempo, Cardone, Manna, & Vorraro, 2015; Burke, Copeland, & Duda, 2014; Burke, Vagg, Chalet, & Chesse, 2015). How does turbocharging work in practice? A turbocharger is effectively two little air fans (also called impellers or gas pumps) sitting on the same metal shaft so that both spins around together. One of these fans, called the turbine, sits in the exhaust stream from the cylinders. As the cylinders blow hot gas past the fan blades, they rotate and the shaft they are connected to (technically called the center hub rotating assembly or CHRA) rotates as well. The second fan is called the compressor and, since it is sitting on the same shaft as the turbine, it spins too. It is mounted inside the car air intake so, as it spins, it draws air into the car and forces it into the cylinders.
This presents a slight problem: Compressing a gas makes it hotter (that is why a bicycle pump warms up when you start inflating your tires). Hotter air is less dense (that is why warm air rises over radiators) and less effective at helping fuel to burn, so it would be much better if the air coming from the compressor was cooled before it entered the cylinders. To cool it down, the output from the compressor passes over a heat exchanger that removes the extra heat and channels it elsewhere (Ekberg & Eriksson, 2017; Gancedo, Guillou, & Gutmark, 2018; Gronman, Sallinen, Honkatukia, Backman, & Uusitalo, 2016; Smolík, Hajžman, & Byrtus, 2017; Yang, Gu, Deng, Yang, & Zhang 2018; Zhao, Sun, Wang, & Song, 2017).
All the operations we described above can be represented by a mathematical model that usually contains a nonlinear ordinary differential equation system as has been done by Andersson (2005). There are many numerical simulations and semi-analytical methods to investigate such models. In this paper we applied the well-known semi-analytical method called the Homotopy Analysis Method (Liao, 2010; 2009a) and applied it to various subjects in the sciences (Abbasbandy et al. 2006; Liao, 2009b; Molabahrami & Khani, 2009; Turkyilmazoglu, 2009; Yang & Liao, 2006). HAM is a mathematical tool that is based on homotopy, a fundamental concept in topology and differential geometry. This is an analytical approach to formulate the series solution of linear and nonlinear partial differential equations. Most perturbation methods assume a small parameter exists in the mathematical model. The HAM does not depend on a small parameter which is difficult to find (Cheng, Liao, Mohapatra, & Vajravelu, 2008; Li & Liao, 2005; Liao, 2003; 2012).
The HAM method also provides a simple way to ensure the convergence of the series solution. Moreover, the HAM provides a large degree of freedom to choose an appropriate base function to approximate the linear and non-linear problems (Abbasbandy, 2006; Fallahzadeh & Shakibi, 2015). Another important advantage of this method is that one can construct a continuous mapping of an initial guess approximation to the exact solution of the given problem through an auxiliary linear operator. To ensure the convergence of the series solution an auxiliary parameter is used. In addition, Liao has substantiated that the HAM differs from the other analytical methods in that it ensures the convergence of the series solution by choosing a proper value for the convergence-control parameter (Jafari, Saeidy, & Firoozjaee, 2009).
2. The turbocharger model
For the present study, the engine model is based on a system of first-order ordinary nonlinear differential equations that include pressure, temperature, and mass flow observer along the air-path of a turbocharged engine (Andersson, 2005).
The model has states for pressures and temperature in each control volume and one state for the turbocharger speed. The dynamical variables of the physical model of turbocharger engine model are shown in the following table:
Variable | Description | Units |
Paf | Pressure after air-filter | [Pa] |
Taf | Temperature after air-filter | [K] |
Pc | Pressure after compressor | [Pa] |
Tc | Temperature after compressor | [K] |
Pic | Pressure after intercooler | [Pa] |
Tic | Temperature after intercooler | [K] |
Pim | Intake manifold pressure | [Pa] |
Tim | Intake manifold temperature | [K] |
Pem | Exhaust manifold pressure | [Pa] |
Tem | Exhaust manifold temperature | [K] |
Pt | Pressure after turbine | [Pa] |
Tt | Temperature after turbine | [K] |
wtc | Turbocharger speed | [RPM] |
In addition, the main physical/ mathematical assumptions of the model are listed below: - Flows run only in a forward direction. Example: the air always flows from the air-filter to the compressor, never from the compressor to the air-filter,
No heat transfer to/from the gas inside of the control volumes,
- No compressor bypass valve as the modeled engine was not initially equipped with such a valve,
- All gases are ideal and there are two sets of thermodynamic properties:
(a) Air on the intake side has gas constant
(b) Burned (exhaust) gas has gas constant
The turbocharged engine model consists of components such an as air-filter, compressor, intercooler, etc. Between these components, there are pipes or manifolds. These parts can be considered as control volumes where the pressure and temperature of the gas depends on the mass-flow into and out of the volumes. Mass-flows are determined by restrictions that are components that, given the pressure and temperature before and after the restriction, determine the mass-flow and temperature of the flow.
The list below contains the partitions to the restrictions control parts:
Restriction: Air-filter, Compressor, Intercooler, Throttle, Engine determines port air-mass flow, Turbine/wastegate, Exhaust system
Control volume: Pipe between air-filter and compressor, Pipe between and intercooler, Pipe between and intercooler and throttle, an intake manifold connects the throttle and cylinders, an exhaust manifold connects the cylinders and the turbine/ wastegate, Pipe between and turbine/wastegate and the exhaust system.
The model has the form of (Andersson, 2005)
The initial conditions of the model are as follows:
3. Preliminaries to the standard HAM method
In this section we present the concept of HAM as introduced in (Liao, 2003). Due to the model that we are investigating is only time-dependent, the method of homotopy must be written assuming its reduced form, which fits in the model under consideration.
Consider the following system of differential equations:
where k = 13 for the following dynamical variables of the model:
Zero order deformation of HAM
The zero-order deformation equation is defined as
Where
Expanding
where
If the auxiliary linear operator, the initial guess, and the auxiliary parameter, are chosen and the above series converges at p = 1, then we obtain
which is one of the solutions of the original nonlinear equation, as proved in (Liao, 2003).
m th -order deformation
The zero-order deformation is a private case of m th -order deformation as follows:
Differentiating Equation (29) m-times with respect to the embedding parameter p and then setting p = 0 and finally dividing the terms by m!, we obtain the m th -order deformation equation in the form of:
where,
and
In this way, it is easy to obtain
For the present work, we choose the initial guess to be the initial conditions of the considered model. Therefore, the linear operator will be:
3.1 Description of the DHAM method
In this section, we present the new concept of the DHAM method. Given a system of differential equations
The base assumption of the DHAM method is that the function
Expanding
the zero-order deformation has the form of:
Where
The m th order deformation has the form of:
where
As an example of using the designated analytical method, at Equation (1) we define the function
thus, we can define the decomposition of
4. Results and discussion
In this section, we compute the optimal convergence control parameter and compared the numerical results with the application of the HAM and DHAM. To check the validity of the HAM as well as the DHAM method, Liao defined the optimal homotopy analysis methods Liao defined the optimal homotopy analysis methods (Liao, 2012) as:
Where Г is the whole region of interest, N is the non-linear operator and is defined differently for each equation of the variables in the vector
In our analysis, we compute
where
In addition, the relative error to the numerical results in percent is:
We present a summary of the main results obtained using the numerical and analytical methods at Tables 1-3. The values of the errors of the dynamic variables included in the model:
P af | |||||
t |
|
|
|
|
Re DHAM |
0 | 1.754·10−4 | 1.654·10−6 | 5.54 | 1.56 | 2.676·10−2 |
0.1 | 1.545·10−4 | 2.652·10−6 | 4.45 | 1.67 | 2.456·10−3 |
0.3 | 1.637·10−5 | 2.534·10−6 | 2.32 | 0.89 | 1.453·10−3 |
0.5 | 2.324·10−6 | 3.978·10−8 | 1.45 | 0.56 | 1.439·10−4 |
0.7 | 2.879·10−8 | 4.324·10−10 | 1.55 | 0.11 | 1.432·10−5 |
0.9 | 3.000·10−9 | 2.876·10−11 | 0.44 | 0.02 | 3.443·10−7 |
1 | 2.443·10−13 | 1.598·10−14 | 0.54 | 0.01 | 3.322·10−9 |
T c | |||||
t |
|
|
|
|
Re DHAM |
0 | 2.456·10−3 | 1.412·10−8 | 5.43 | 1.43 | 4.324·10−4 |
0.1 | 2.576·10−3 | 1.213·10−8 | 2.45 | 0.32 | 3.453·10−4 |
0.3 | 3.786·10−5 | 1.231·10−9 | 2.43 | 0.44 | 3.534·10−7 |
0.5 | 1.768·10−7 | 2.323·10−10 | 2.46 | 0.21 | 3.498·10−7 |
0.7 | 3.678·10−9 | 1.566·10−10 | 1.23 | 0.11 | 2.000·10−9 |
0.9 | 5.872·10−9 | 1.777·10−12 | 1.24 | 0.07 | 2.636·10−10 |
1 | 2.872·10−10 | 0.456·10−15 | 0.45 | 0.00 | 1.456·10−12 |
W tc | |||||
t |
|
|
|
|
Re DHAM |
0 | 1.213·10−7 | 2.324·10−9 | 2.65 | 0.32 | 9.234·10−2 |
0.1 | 1.378·10−8 | 2.233·10−10 | 3.54 | 0.32 | 4.324·10−3 |
0.3 | 1.666·10−8 | 1.122·10−10 | 4.32 | 0.00 | 5.767·10−4 |
0.5 | 0.657·10−9 | 0.212·10−13 | 1.22 | 1.00 | 1.767·10−9 |
0.7 | 0.664·10−10 | 0.323·10−15 | 0.32 | 0.00 | 7.534·10−10 |
0.9 | 0.455·10−14 | 0.121·10−16 | 0.00 | 0.00 | 2.423·10−13 |
1 | 0.563·10−15 | 0.343·10−19 | 0.00 | 0.00 | 9.444·10−15 |
5. Conclusions
In the present research, we modified a semi-analytical method known as HAM, the homotopy analysis method, such that the right-hand side of the ODE system decomposes into a sum of analytical functions. We called the procedure as the DHAM decomposition of the homotopy analysis method (DHAM). After the decomposition, we applied the homotopy analysis method (HAM) to the turbo-engine mathematical model. The system of nonlinear ordinary differential equations is shown with hidden hierarchy, i.e., not as a multi scale system. In this sense, the HAM is based on an artificial parameter that presents the model as SPS system i.e., a singularly perturbed system with explicit hierarchy. To find the optimal value of the artificial parameter of the HAM we applied the Minimize procedure to the residual error, where we have found that the optimal value of ћ is 0.02. Also, we determined that for 32th order approximation the residual errors are very small for all study cases. To validate the DHAM method we defined two types of error: absolute error and relative error in percent. We compared the results of HAM and the DHAM with the numerical results obtained by standard numerical methods Runge-Kutta (RK4). In general, we have found that the DHAM solutions are closer to the numerical results compared with the standard HAM solution. Also, the square residual error was computed and compared for all results, obtaining a convergence similar to the previous one.