`Advances in Mechanical EngineeringVolume 2012 (2012), Article ID 931584, 6 pageshttp://dx.doi.org/10.1155/2012/931584`
Research Article

## Temperature Estimation in the Combustion Chamber of an Internal Combustion Engine

Department of Mathematics, Bu Ali Sina University, Hamedan, Iran

Received 1 December 2011; Accepted 1 June 2012

Copyright © 2012 Gholam Reza Safakish. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The mathematical model of heat transfer phenomena is considered at the combustion chamber wall in an internal combustion (IC) engine. The mathematical model of proposed phenomena is established with respect to the crank angle. An inverse heat conduction problem is derived at the cylinder wall, and this problem is investigated numerically using Alifanov's regularization method. This problem studied as an optimization problem in which a squared residual functional is minimized with the conjugate gradient method. To show the ability of the proposed method, some test problems are considered.

#### 1. Introduction

The solution of inverse heat transfer problems is becoming an essential part in the development of several relevant applications in engineering, such as design of thermal equipment, systems, and instruments. The use of inverse methods represent a new research direction, where the results obtained from numerical simulation and from experiments are not simply compared a posteriori, but a closed synergism exist between experimental and theoretical researches during a course of study, in order to obtain the maximum of information regarding the physical problem under consideration. Most of the methods for the solutions of inverse problems, which are currently in common use, were formalized in the last four decades in terms of their capabilities to treat ill-posed unstable problems. The basis of such formal methods resides on the idea of reformulating an inverse problem in terms of an approximate well-posed problem, by utilizing some kind of regularization techniques [110].

In the inverse analysis, most studies employed the nonlinear least-squares method [1012] to determine the inverse problem. This method minimizes the formulation from the sum of the squares of the difference between the experimental measurements and the calculated response of the system. Based on the nonlinear least-squares method, various researchers have put their efforts in the field of inverse problems. In solving the problems, different algorithms have been adopted such as the conjugate gradient method, the Davidson-Fletcher-Powell method, the Monte-Carlo technique, the covariance analysis, and the dynamic programming. More sophisticated methods also have been developed such as the nonlinear least-squares formulation modified by the addition of a regularization term, the sequential estimation approach, and the adjoint equation approach coupled to the conjugate gradient method [17].

Motivation for the problem investigates here arises from the area of modeling of heat transfer phenomenon at the combustion chamber wall in an internal combustion (IC) engine. Heat transfer between the working fluid and the combustion chamber in an IC engine is one of the most important parameters for cycle simulation and analysis. High temperatures request to combust the fuel, so it is necessary to keep the temperature at a controllable level in order to operate the engine safely. Once the temperature in the engine has reached intolerable values, the engine blocks and components may suffer damage. Therefore it is essential to have a heat removal process which will maintain the engine at a safe operating condition. Heat transfer measurement by thermal sensors such as thermocouples lead to poor bandwidths and large uncertainties [14]. After the compression stroke and during combustion stroke, there is heat transfer to the surroundings from the hot gas through the cylinder walls. The peak gas temperatures of combustion are the order of 3000 (k) and that is why the cylinder walls of the chamber overheat. The only way in which energy can be transferred away from the combustion chamber is through convection and conduction. The heat transfer phenomena in internal combustion engines has been investigated by many authors [16].

In comparison with previous works, The heat transfer problem in the cylinder wall is studied by using the crank angle information. Heat transfer and cylinder pressure are two closed problems in IC engines. Cylinder pressure can be analyzed as a function of crank angle, so for this reason, the amount of temperature in the cylinder chamber can be state as a function of crank angle [1, 2, 7]. The goal of this study is estimation of internal temperature as a function of crank angle by solving an inverse heat conduction problem (IHCP). It is assumed that internal temperature are unknown function of crank angle and at the external surface heat flux are known. In addition we suppose that temperature at the some angles are available as overspecified data.

The inverse problem considered in this work is solved by using a function estimation approach [46], where no information is a priori available regarding the forms of the unknown function, except for the functional space that they belong to. It is assumed that the unknown function belong to Hilbert space of square integrable function in the spatial domain of interest ().

The solution of inverse problems by using the conjugate gradient method with adjoint problem for function estimation consists of the following basic steps: (i) direct problem formulation, (ii) inverse problem formulation, (iii) sensitivity problem formulation, (iv) adjoint problem formulation, (v) gradient equations, (vi) iterative solution procedure, (vii) iterative procedure stopping criterion, and (viii) computational algorithm. Highlights of such steps are presented below as applied to the inverse problem of interest.

#### 2. Direct Problem Formulation

The physical model is presented in Figure 1. The following hypotheses have been taken into account. (i)A cylinder of internal radius and external radius is considered. (ii)The thermophysical properties are constant, and heat transfer is one-dimensional. (iii)Thermal and physical properties are considered constant.(iv)The cylinder wall is treated as thin rectangular sheet.

Figure 1: A simple view of cylinder chamber of an IC engine.

Under these conditions, the mathematical model of heat transfer process in cylinder wall is given by the following system of equations: where is the specific heat, is density, is hermal conductivity,

is the initial temperature, is the temperature of internal surface, and is the heat flux at the external surface of the combustion chamber wall. In an combustion chamber let show the crank angle, TDC is top dead center (piston position), BDC is bottom dead center (piston position), DC is dead center, and is combustion chamber volume. One may see that the crank angle varies with time; that is, , so we have

In the engine is where is the engine angular speed, and is the engine speed Therefore, for the one cycle, the problem (1) can be rewritten as follows:

Now, suppose Then, the problem (4) is reduced to the following standard heat conduction problem where

#### 3. Inverse Problem

In the problem (6)–(9) suppose that and are known functions and is an unknown function. In addition as an overspecified data, suppose that the transient measurements of taken at the location , and are available as .

As mentioned before, the peak gas temperatures of combustion are the order of 3000 (k) and for this resason, the exact measurement of internal wall temperature is impossible. To this end, for the estimation of internal temperature, that is, , we make use of a minimization procedure involving the following objective functional: where is the temperature at and computed from the solution of the direct problem (6)–(9), by using the estimated values for .

This inverse heat conduction problem (IHCP) is recast as an optimum control problem of finding the unknown control function such that it minimizes the functional (11) [3, 811].

#### 4. Sensitivity Problem

The sensitivity function, solution of the sensitivity problem, is defined as the directional derivative of in the direction of the perturbation of the unknown function [12]. The sensitivity problem is obtained by assuming that the dependent variable is perturbed by when the function is perturbed by . Here, is a real number. The sensitivity problem for is then obtained by applying the following limiting process: where and are the direct problem formulations written in operator form for perturbed and unperturbed quantities, respectively. The application of the limiting process given by (12) results in the following sensitivity problem:

A Lagrange multiplier is utilized in the minimization of the functional (11), because the estimated dependent variable appearing in such functional needs to satisfy a constraint, which is the solution of the direct problem. Such Lagrange multiplier, needed for the computation of the gradient equations (as will be apparent below), is obtained through the solution of problems adjoint to the sensitivity problem (13) [2, 3]. In order to derive the adjoint problem, the governing equation of the direct problem (6) is multiplied by the Lagrange multiplier , integrated in the space and time domains of interest and added to the original functional (11). The following extended functional is obtained:

Directional derivatives of in the directions of perturbations in is defined by where denote the extended functional (14) written for perturbed . After performing some lengthy but straightforward manipulations and letting the directional derivatives of go to zero, which is a necessary condition for the minimization of the extended functional (14), the following adjoint problem for the Lagrange multiplier is obtained:

In this section, we consider that the unknown function is approximated in the form of cubic B-spline as follows: where is the number of approximation parameters, are unknown approximation parameters, and are given basis cubic B-spline.

As a result the IHCP is reduced to the estimation of a vector of parameters . Now, we consider the conjugate gradient iterative procedure to obtain the estimated that minimizes the functional (11). Unknown function is considered as an element of space.

The iterative procedure of the conjugate gradient method is written as where is the iteration index, is descent parameter and a positive scaler, and is the descent direction that can be found as follows: where the parameter is computed as follows: where and respectively are the usual scaler product and norm in the space of approximation parameters. The realization of the iterative procedure (18) is based on the computing the vector at each iteration. This vector is determined by the relationship for the residual functional variation: where is the residual functional gradient in space of approximation parameters and is the gradient in space of parameterized functions. The vector is determined by the following equation: where is the Gram's matrix for basis functions This matrix is positive definite and symmetric. The most effective method for calculating the adjoint in space is based on introducing an adjoint problem [8, 11].

During the limiting processes used to obtain the adjoint problem, applied to the directional derivatives of in the directions of perturbations in , the following integral term is obtained:

The gradient of the residual functional is determined by solving adjoint problem. The following expression for the gradient components can be derived: Determination of descent parameter in each iteration deals with solution of sensitivity problem.

Descent parameter is determined minimizing the residual functional (11) with respect to . By using a Taylor series expansion, can be derived as follows: The iterative regularization method is used here to obtain stable solution of the proposed inverse heat transfer problem. The main idea is to terminate the iterative procedure with the residual criterion: where is defined by where is an estimate of the standard deviation for temperature measured at points .

#### 7. Results and Discussion

In order to examine the ability and accuracy of this model and the inverse analysis presented here for the estimation of the internal temperature, in problem (6)–(9) we consider that

Here we try to recover the function as follows

Since experimental data were not available we generate simulated transient temperature data , by using the computed exact temperature , as follows: where is a random number.

Figures 1, 2, 3, and 4 show the results after 30 iterations when we use , , and basis functions for estimating as an unknown function. Figure 5 shows the behavior of the functional when the number of basis functions is increased.

Figure 2: The comparison between exact and approximate solutions using cubic B-spline basis functions.
Figure 3: The comparison between exact and approximate solutions using cubic B-spline basis functions.
Figure 4: The comparison between exact and approximate solutions using cubic B-spline basis functions.
Figure 5: The behavior of the functional with respect to the number of basis functions.

#### 8. Conclusion

This paper presents a new differential model of heat conduction in the cylinder wall of internal combustion engine based on crank angle. Since a lot of information about cylinders problem such as in-cylinder pressure and temperature amount could be found by using crank angle, it is convenient to use crank angle to estimate temperature at the internal surface of cylinder wall. This model can be used to estimate heat flux and temperature at the other parts of engine that deal with combustion chamber as a function of crank angle.

#### References

1. J. K. Ball, R. R. Raine, and C. R. Stone, “Combustion analysis and cycle-by-cycle variations in spark ignition engine combustion. Part 1: an evaluation of combustion analysis routines by reference to model data,” Proceedings of the Institution of Mechanical Engineers D, vol. 212, no. D5, pp. 381–399, 1988.
2. C. Luo, B. Moghtaderi, S. Hands, and A. Page, “Determining the thermal capacitance, conductivity and the convective heat transfer coefficient of a brick wall by annually monitored temperatures and total heat fluxes,” Energy and Buildings, vol. 43, no. 2-3, pp. 379–385, 2011.
3. J. B. Heywood, Internal Combustion Engine Fundamentals, McGraw-Hill, New York, NY, USA, 1988.
4. J. Su and A. J. S. Neto, “Two-dimensional inverse heat conduction problem of source strength estimation in cylindrical rods,” Applied Mathematical Modelling, vol. 25, no. 10, pp. 861–872, 2001.
5. K. S. Chana, T. S. Wilson, and T. V. Jones, “High bandwidth heat transfer measurements in an internal combustion engine under low load and motored conditions,” in Proceedings of the RTO AVT Symposium on Advanced Flow Management. Part B: Heat Transfer and Cooling in Propulsion and Power Systems, Loen, Norway, May 2001.
6. A. C. Alkidas, Puzinauskas, P. V. Peterson, and R. C. Combustion, “and heat transfer studies in a spark-ignited multivalve optical engine,” SAE Paper 900353, 1990.
7. J. Kou, K. P. Judd, and J. R. Saylor, “The temperature statistics of a surfactant-covered air/water interface during mixed convection heat transfer and evaporation,” International Journal of Heat and Mass Transfer, vol. 54, no. 15-16, pp. 3394–3405, 2011.
8. E. C. Trujillo, F. J. Jiménez-Espadafor, J. A. B. Villanueva, and M. T. García, “Methodology for the estimation of cylinder inner surface temperature in an air-cooled engine,” Applied Thermal Engineering, vol. 31, no. 8-9, pp. 1474–1481, 2011.
9. M. N. Özisik, Heat Conduction, John Wiley and Sons, New York, NY, USA, 1993.
10. H. M. Park and W. S. Jung, “A recursive algorithm for multidimensional inverse heat conduction problems by means of mode reduction,” Chemical Engineering Science, vol. 55, no. 21, pp. 5115–5124, 2000.
11. K. L. Hoag, “Measurements and analysis of the effect of wall temperatures on instantaneous heat flux,” SAE Paper 860312, 1986.
12. O. M. Alifanov, E. A. Artyukhin, and S. V. Rumyanstev, Extreme Method for Solving Ill-Posed Problems with Application to Inverse Heat Transfer Problems, Begell House, New York, NY, USA, 1995.