Experimental Thermal Analysis of Diesel Engine Piston and Cylinder Wall

Knowledge of piston and cylinder wall temperature is necessary to estimate the thermal stresses at different points; this gives an idea to the designer to take care of weaker cross section area. Along with that, this temperature also allows the calculation of heat losses through piston and cylinder wall. The proposed methodology has been successfully applied to a water-cooled four-stroke direct-injection diesel engine and it allows the estimation of the piston and cylinder wall temperature. The methodology described here combines numerical simulations based on FEM models and experimental procedures based on the use of thermocouples. Purposes of this investigation are to measure the distortion in the piston, temperature, and radial thermal stresses after thermal loading. To check the validity of the heat transfermodel, measure the temperature through directmeasurement using thermocouple wire at several points on the piston and cylinder wall. In order to prevent thermocouple wire entanglement, a suitable pathway was designed. Appropriate averaged thermal boundary conditions such as heat transfer coefficients were set on different surfaces for FE model. The study includes the effects of the thermal conductivity of the material of piston, piston rings, and combustion chamber wall. Results show variation of temperature, stresses, and deformation at various points on the piston.


Introduction
For proper functioning of the internal combustion diesel engine, accurate piston temperature distribution is required because piston temperature has an important influence on ignition process of engine, ignition time delay, rate of burning, thermal efficiency, and production of pollutants. Knowledge of heat transfer in internal combustion engines is important to understand such systems [1,2]. It contributes to engine development and design, processes simulation, and emissions reduction. In engine piston experiences high forces due to combustion chamber pressure and thermal load, which come from combustion process and from the huge temperature gradient between intake and exhaust gas flows [3][4][5][6] so, it is important to guarantee the durability of engine components like piston, piston rings, valves, and cylinder wall, to avoid engine body distortions and to improve engine design related to weight and auxiliary energy consumption. In the case of the engine piston and cylinder, such knowledge is indispensable to have a thorough understanding of heat flux, temperature, and the distribution of these parameters. A common procedure used by some authors is to approximate the mean temperature of the distribution with one or more (very few) local measurements obtained by thermocouples [1,7]. These approaches implicitly assume errors that can be acceptable for heat balances but can lead to uncertainties in simulation cycles or temperature distribution analyses. Many of these models include the gas-side wall temperature as a variable to obtain the heat flux through the cylinder walls [8]. Numerous mathematical models have been proposed including correlations based on dimensional analysis, which are widely accepted. Although the models propose different 2 Journal of Engineering heat fluxes, their evolution over the cycle is similar. In addition, finite element method (FEM) codes used for heat transfer simulations require the estimation of temperature to provide boundary conditions where convergence is attained through an iterative process [7]. A finite element model of gasoline spark engine is successfully developed and simulated and had analyzed heat transfer during combustion process and obtains temperature distribution across the major engine component [9]. Furthermore, thermal analyses require the gas-side wall temperature to evaluate temperature distribution and the thermomechanical behaviour of components with the use of thermal barrier coating [10][11][12][13]. Other researchers, having identified one correlation for one small-scale air-cooled engine [14], acknowledge that these correlation parameters are not valid for another smallscale air-cooled engine with similar characteristics [15]. In most heat transfer analyses the external surface temperatures from which heat is extracted are not measured. This is the case for water-cooled engines, where this temperature is often assumed to be equal to the coolant temperature or is calculated from hypotheses that are specific for water-cooled engines [16] and normally assume constant temperature for all operating points. This approach cannot be extrapolated to air-cooled engines because the temperature field in the inner surface changes with the operating conditions [17]. Some studies concerning heat transfer in air-cooled engines have been published and most of them consider two-stroke and/or spark-ignition engines [18][19][20].
In this paper a methodology is proposed for the estimation of the temperatures in piston and cylinder wall, piston body distortions, and radial stresses of a water-cooled fourstroke single cylinder direct injection diesel engine. These temperatures were obtained by simulation and experimental setup with the help of thermocouples and other sensors and other calculations that are particularized to that engine.

Engine Description and Experimental Procedure
The experimental investigation was conducted on a fourstroke single-cylinder direct-injection (DI) compression ignition engine. The main characteristics of this engine are given in Table 1 and a cutaway view of the piston and cylinder is shown in Figure 1. The experimental measurements covered four different loading conditions, namely, no load, half load, three-quarters load, and full engine load. The temperature of the piston can be estimated by direct measurement using thermocouple as well as with the help of numerical method, which consists of finite element analysis method. The purpose is to obtain a generalized method (FEM) for analyzing temperature field, piston distortion, and corresponding thermal stresses so simulated temperatures are to be verified by direct measured temperatures. In this work, seven thermocouples were used, in which four thermocouples were mounted on the piston inner surface and three were mounted on the cylinder wall as shown Figure 2. The coordinates of these seven nodal points where thermocouples were attached are shown in Table 2. In the experimental setup engine, piston is made of aluminium alloys of MSFC-388-T5 having a thermal conductivity equal to 136 W/mK.
For direct measurement of the temperature at piston inner surface a four-bar mechanism was designed, which is shown in Figure 3. With the help of this four-bar mechanism, thermocouple wires safely come out from the engine chamber. Complete experimental setup has been shown in Figure 4. By installing thermocouple at seven points on piston and cylinder wall, variation in temperature of piston and cylinder wall at no load, half load, three-quarters load, and full load conditions was determined. For analysis of stresses and distortion of the piston, two sections 1-1 and 2-2 have been chosen as shown in Figure 1

Thermomechanical FE Analysis
The method is based on the subdivision of the structures into elements with mathematically defined characteristics. Characters of the complex structure are then solved with the aid of computer using matrix algebra. The input data consist of loading conditions on subdivided structures or mesh constraints with physical properties of the material. Thermal loading involves the initial calculations of temperatures using boundary conditions on gas side, coolant side, and air side of piston cylinder. The analysis presented in this paper is divided into two sections, the temperature field distribution and the thermal stresses. The finite element technique with triangular element is used to reduce the variational formulation to a set of algebraic equations. The expressions to calculate nodal temperatures and the corresponding thermal stresses at every element are derived. The construction of finite element approach starts from the variational statement of the problem and then using proper shape function a number of algebraic equations are developed which are equal to the number of nodal elements in the problem domain. Then by minimizing the approximate function a set of governing 2D view of axis symmetric quarter piston equation is developed for the piston and cylinder assembly. These equations are solved by using computer. Computer algorithm and a FORTRAN program code are developed to solve these equations in order to find the unknown parameters, that is, temperature at different nodal points of the piston. Computer program is based on heat transfer through conduction, convection, matrix multiplication, matrix inversion, heat flow, and stiffness. By using these subroutine and main program, temperatures and heat flow field were calculated. Mathematical modeling of conduction equation, convection equation, and contact heat transfer equation has been done which is shown below. The generalized governing differential equation for heat conduction can be represented as [7,21] where is thermal conductivity in radial ( ) and axial ( ) direction, respectively. is heat conduction per unit volume. is density of the material. is heat capacity of the material. is temperature. is time.
The variational formulation for conductive boundary can be represented as [7,21]  ] where [ ] ( ) = stiffness matrix. The generalized governing differential equation for contact boundary can be represented as [7,21] = ℎ ( − ) . ( The variational formulation for contact boundary between 2 elements ( ) and ( ) can be written as On solving it further in similar fashion as done in convective boundary, it has been found that contact boundary variational integral of heat transfer after differentiation with respect to temperature of contact surface yields a set of linear equations as a contribution to the global set of equation. Consider The generalized governing differential equation for heat convection can be represented as [7,21] − ( ) = ℎ ( − ∞ ) .
The variational formulation for convective boundary can be represented as where bconv.
From (2), (5), and (9) variational integral of heat transfer globally is developed and represented as follows: This equation can be written in its popular form as Here [ ] is conduction matrix, [ ] is convection matrix, { } is unknown temperature column vector at all nodal points, and {ℎ } is a column vector of known quantity. All above matrixes are global matrix of size ( × ), where is the number of nodes. The column vector {ℎ } is of size ( × 1).
From (11), temperatures at all the nodes of the piston are to be found out, which is to be represented in Figures 7-10. After prediction of the temperature at all the nodes of the piston, radial thermal stress will be analyzed with the help of (12), (13), and (14) to show the radial strain, angular strain, and axial strain, respectively, as follows: Thus stress strain relationship can be represented in matrix form as shown in where [ ] is the position matrix, = radial stress, = angular stress, and = axial stress. After inversion of position matrix, put the value of this in the above equation and find the stresses at all the nodes of the piston. Consider where = thermal stress, thermal strain { 0 } = ( − 0 ), and simple strain { } = { }; V is the Poisson ratio and it is 0.33 and From (12), (14), and (19), radial displacement, axial displacement, and radial thermal stresses at all the nodes of the piston have been calculated and are represented in Figures  11-16.

Finite Element Formulation for Checking Heat Balance.
This formulation can be used to check the heat balance of the problem. The amount of heat balanced is checked for accuracy and satisfied with result by observing that the amount of heat supplied at gas side of the piston is equal to the amount of heat loss to both water and air side of the piston.
Consider one side of element having two nodal points , , which faced the boundary where heat is either supplied or rejected. The formulation of the equation is as follows. The generalized heat transfer equation from connective surfaces is represented as follows: where where (24) By using (24), heat transfer through the different surfaces can be determined easily. Heat transfer from piston through combustion chamber side, air side, and water side is represented by , , and , respectively.

Development of a Geometrical and a FE Model
A geometrical model of the piston was developed based on geometry of the actual object. Properties of the geometrical objects allowed determining such parameters as the piston diameter, the diameter of the piston pin hole, and dimensions of the piston ring grooves. A final CAD model is presented in Figure 1. In this model, a certain geometric simplification was assumed, including missed bends with a small radius on the edge of crown and lateral surface of the head of the piston. Finally, the geometrical model was discretized into tetrahedron finite elements. Such elements had to be applied due to a complex shape of the piston. The size of finite elements was different in respective sections of the piston; larger elements were used for the piston crown and skirt, whereas the smaller ones were employed close to the oil channels. The total number of nodes and elements in the quarter portion of piston FE model was  equal to 311 and 272, respectively. Figure 5 shows the FE discrete model of the piston.

Thermal Boundary Conditions
The thermal boundary conditions consist of applying a convection heat transfer coefficient and the bulk temperature, and they are applied to the piston crown, piston ring land sides, piston ring groove lands, and piston under crown surfaces. The temperature and heat transfer coefficients in the combustion chamber in all the loading condition were identified based on data provided in previous research paper [7,13] which are presented in Table 3. The adopted heat transfer coefficient on the contact surfaces are (heat transfer coefficient at piston under crown surface) = 174.3 w/m 2 k, 1 (heat transfer coefficient at ring lands and piston skirt upper and lower side) = 2905.4 w/m 2 k, 2 (heat transfer coefficient at ring lands and piston skirt contact surfaces) = 20 w/m 2 k, 3 (heat transfer coefficient between piston rings and cylinder wall contact surfaces) = 38346 w/m 2 k, 4 (heat transfer coefficient between piston and cylinder wall contact surfaces) = 2324 w/m 2 k, (heat transfer coefficient through cylinder wall to water) = 1859.2 w/m 2 k, and temperature on water side ( ) was 120 ∘ C and on crank case side ( ) was 80 ∘ C.

Results
The temperature distribution and heat flow field were studied with the help of finite element analysis using thermal boundary conditions. To check the validity of the heat transfer model, the heat balance approached was adopted. According to the principle of conversion of energy, at steady state condition the heat entering into the piston from gas side is equal to the heat rejected to water and heat rejected to air. Figure 6 shows the variation of heat received ( ) by the piston from the hot gases, heat rejected to water ( ), and heat lost to air ( ) at the four different thermal loading conditions. It seems that heat received from the hot gas is increased with increase in the engine combustion temperature ( ). Similarly heat rejected to water and heat rejected to air increase with engine combustion temperature. Here in the present analysis the heat balance equation is satisfied for all the different loading conditions which were shown in Figure 6. In Figure 6, it can be observed that the error is very small between heat supplied ( ) by the combustion gasses and heat rejected to water ( ) and air sides ( ). It represents a uniform percentage throughout all the tests, shown in Table 4. As expected, for steady state conditions the heat transfer rate increases with engine loading condition, with a maximum observed at full load. Figures 7-10 show the temperature distribution. Through the analysis authors found that the maximum temperature occurs at crown of the piston, because it is exposed to the hot gases in the combustion chamber, while the minimum temperature occurs at the carter end of the piston, which is exposed to the air. These variations in temperature are mainly responsible for the development of temperature stress, causing appearance of crack in the body of the piston. Experimental setup established and measured the temperatures through temperature sensors at seven testing points (nodal points) where temperature sensors were mounted, shown in Table 5. As expected, for all load conditions the simulated and experimental measured temperature increased with engine loading condition, with a maximum observed at full load. Three-quarters load The nonuniform variation in temperature is solely responsible for the development of thermal strain on piston body. Initially the piston is at 25 ∘ C and it is assumed that initially there is no thermal strain and stresses are present. Here in the analysis first the displacement for each element is obtained and the compatibility equation is found to be satisfied. Being CST the strain in the whole element is constant and assumed to be located at the center gravity of the element. Thus an element either expands or contracts owing to its temperature variation to maintain the constant strain in whole element. The thermal strains are related to stresses by means of Hook's law of linear isothermal elastically. Thus the thermal stress developed depends directly on the mean temperature. The present analyses on two different sections on piston body are chosen to study the behavior of the stress. One section on the piston body is chosen near the piston crown (section 1) having 7 nodes and the other at bottom of the body (section 2) having 7 nodes. Figures 11 and 12 show the radial displacements of the nodal points with respect to the piston radius under four thermal loadings at section 1 and section 2, respectively. It was found that the radial displacements of the nodes of these sections are varied in increasing order from centre to the circumference of the piston in all four loading conditions. And it is observed Three-quarters load that the displacement of section 1 is more than that of the section 2 due to high temperature. Figures 13 and 14 show the vertical displacement of the nodal points with respect to the piston radius under four thermal loadings at sections 1 and 2, respectively. Figures 13 and 14 indicate a sharp variation in magnitude and direction of vertical displacement of nodal elements. This is because every element tends to maintain the constant strain property. If a node expands for a certain element then it must compress for another element. Figures 15 and 16 show the variation of radial stresses with respect to the piston radius under four thermal loadings at sections 1 and 2, respectively. These stresses are introduced in the body due to high temperature. Figures 15 and 16 show a sharp variation in magnitude and direction of a radial stress. As discussed it is because every element tends to maintain the constant strain property. Thus it undergoes a simultaneous expansion or compression between two adjacent elements. However by considering a very large number of elements these sharp variations in stress can be eliminated and a more continuous curve can be obtained. Generally at high temperature piston is distorted slightly in outward direction, and this was shown in Figure 17. In this analysis, distortion in the piston body was observed at full loading condition. The maximum value of the dislocation is approximately 8 × 10  meter and it was found nearer to the piston ring land sides, due to low strength at this part of the piston.

Conclusions
As engine load increases, temperature of the piston and cylinder wall increases exponentially and has a positive relationship. The piston temperature for every engine load condition tested was estimated and good agreement was obtained with the expected results. These results are also consistent with those described in the technical literature.
The proposed methodology can also be extended to determine temperatures of other components of the combustion chamber, of course taking into account its particularities. These temperatures, together with experimental measurements and the estimations calculated, can be used to obtain