Digital Commons@Becker Digital Commons@Becker Interring gas dynamic analysis of piston in a diesel engine Interring gas dynamic analysis of piston in a diesel engine considering the thermal effect considering the thermal effect Interring Gas Dynamic Analysis of Piston in a Diesel Engine considering the Thermal Effect

Understanding the interaction between ring dynamics and gas transport in ring pack systems is crucial and needs to be imperatively studied. The present work features detailed interring gas dynamics of piston ring pack behavior in internal combustion engines. The model is developed for a ring pack with four rings. The dynamics of ring pack are simulated. Due to the fact that small changes in geometry of the grooves and lands would have a significant impact on the interring gas dynamics, the thermal deformation of piston has been considered during the ring pack motion analysis in this study. In order to get the temperature distribution of piston head more quickly and accurately, an efficient method utilizing the concept of inverse heat conduction is presented. Moreover, a sensitive analysis based on the analysis of partial regression coefficients is presented to investigate the effect of groove parameters on blowby.


Introduction
The piston ring performs its function as a seal of high and low pressure sides in a cylinder and a medium of heat transfer from piston to cylinder walls. The operation of piston ring packs influences the performance, efficiency, durability, and emissions of engines in terms of friction, wear, oil consumption, and gas blowby. For those reasons, the performance of piston rings in combustion engines has been a topic of research for many years. Ting and Mayer [1,2] developed early models of ring lubrication and blowby. Their models were aimed at the eventual prediction of cylinder wear. The work of Dowson et al. [3] and Ruddy et al. [4,5] has resulted in the development of more integrated simulations, which include effects of ring dynamics. Rhode [6] has incorporated into ring lubrication calculations a model of mixed lubrication, which allows more realistic calculation of ring friction; it has been so since adopted by a number of other investigators. The effort toward further integration of various submodels into what has become an established scheme of two-dimensional axisymmetric analysis of ring packs has continued in the eighties [7], with more studies involving various design parameters as well as comparison to data [8,9].
Classical theory of lubrication is based on the Reynolds equation, which can be derived from control volume analysis under certain simplifying assumptions. This equation calculates the oil film pressure given the film thickness, squeeze velocity, and the pressure at the boundaries (typically ambient pressure). Then attention had been shifting from ring dynamics and blowby, which can be predicted with reasonable accuracy, to the modeling of oil consumption, the accuracy of the current oil consumption models [9][10][11]. Keribar et al. [12] developed an integrated model of piston ring pack including ring-liner hydrodynamic and boundary lubrication, ring axial, radial, and twist dynamics, interring gas dynamics, and blowby. Blow-by gas flow refers to the undesired gas flow from combustion chamber to the crankcase. It reduces the efficiency of the engine and contaminates oil with the combustion products present in the gases; the ring pack system serves as a channel for blow-by gases. As a result, it is important for engine manufacturers and lubricant suppliers to optimize ring pack system to take it into account in controlling oil consumption resulting from blow-by gases. Such research pursues to decrease the amount of blow-by gas flow while upholding sufficient lubrication and minimum friction.

Mathematical Problems in Engineering
Most of the existing models above are based on the axisymmetric assumption, called 2D models. With the help of these 2D models, a great deal of detailed analysis of ring lubrication has been performed and some conclusions have been used successfully in guiding the practical application [13]. The efforts to model the ring dynamics in 3D manner can also be found in the publications. Tian et al. [14] presented ring dynamics and gas flow model to study ring/groove contact, blowby, and the influence of ring static twist, keystone ring/groove configurations, and other piston ring parameters. Ejakov et al. [15] simulated the dynamic twist angle of the piston ring pack along the circumference by using finite element analysis. In [16], a combination of a detailed elastic analysis and a complementarity method was used to evaluate the ring/bore conformability in an accurate way. Ma [17,18] implemented a flow-continuity algorithm in the hydrodynamic lubrication submodel, which ensured an accurate oil transport simulation. Liu [19] developed a piston ring pack model considering nonaxisymmetric characteristics of power cylinder system considering oil transport along the liner. The oil transport along the ring pack system is driven by the motion of the piston rings [20]. Sometimes, oil is entrained in blow-by gases, transported within the ring pack. The demonstration of the interaction between ring performance and gas transport can help to understand oil transport mechanism and lead to effective ways of reducing the oil consumption. Li [21] took the ring geometry, assembly load, and mechanical and thermal properties into consideration in their simulation of piston ring behavior.
In order to know the ring pack motion, the most important part is to analyse the interring gas dynamic. Due to the fact that small changes in geometry of the grooves, lands, and bore would have a significant impact on the interring gas dynamics, the thermal deformation should be taken into consideration during ring pack motion analysis. In order to determine the thermal deformation of bore and piston, a 2D finite element (FE) model is used in this study. In the thermal analysis part, an efficient method utilizing the concept of inverse heat conduction is presented to get the temperature distribution of piston. After getting the temperature distribution of piston, a thermostructure analysis is presented to get the deformation of piston lands and grooves and bore, as the input data for ring pack dynamics and tribological analysis. The second part of this study focuses on the interring gas dynamics of piston ring pack behavior in internal combustion engines. The model is developed for a ring pack with four rings. The dynamics of the top three rings and the gas pressure in the regions about the oil control ring are simulated.

Modeling of Piston Ring Pack Dynamics
The motion of piston ring within the piston groove can be described by the axial rotational (toroidal twist) and radial motions in the three respective degrees of freedom. The ring motion in the circumference direction is neglected in this study. As shown in Figure 1, a small section for the ring at a where cg and ℎ cg are the axial and radial position of the center of gravity of the ring and is the twist angle. , , and are ring mass, moment of inertia for toroidal rotation, and cross-sectional torsional stiffness.
is the radial tension stiffness, ℎ the reduction in ring radius at installation, and ℎ 0 the minimum ring-bore oil film thickness ( tension = (ℎ + ℎ 0 )). 's and 's are forces and moments acting on the ring cross-section, as shown in the free body diagrams of Figures 1 and 2. In terms of forces and moments in (1) the first subscript indicates location on the ring ( = groove, = cylinder, and = land) and the second subscript describes the source of the force or moment (oil = oil pressure, asp = normal pressure due to asperity contact, pre = gas pressure, and frc = hydrodynamic or boundary friction). The torsional moment is calculated as = . ( For a complete ring with rectangular cross-section, the cross-sectional torsional stiffness ( ) is given by where is the modulus of elasticity of piston ring, is the axial height of piston ring section, is the inner diameter, Mathematical Problems in Engineering and is the outer diameter of piston ring. The axial motion of a ring relative to piston can be expressed as where is the axial position of the piston and is the axial position of the ring related to the piston.
Forces and moment associated with land and groove gas pressure are calculated using pressure solutions from gas dynamics submodel. Those associated with oil pressure are obtained from the lubrication submodel, and those associated with local asperity contact pressure are calculated by asperity deformation model. Ring axial position and twist influence gas flow paths and the forces at the ring groove interface. Ring twist also affects the effective profile presented by the ring face to the cylinder bore and thus oil film thickness and ring friction. The lubrication model, asperity deformation model, and oil film squeezing model at ring side and groove interface have specific introduction in [12,14].

Interring Gas Dynamics
The gas flow between the rings and groove flanks and the flow behind the ring are calculated based on the actual interring volumes due to actual clearances and ring positions. The real mass flow rate is determined from the ideal flow by use of discharge coefficients. The purpose of the interring gas dynamics model is the calculation of the rate of gas blowby through the leakage paths between series of gas volumes. The gas volumes in ring grooves as well as ring-lands are considered and calculated from the piston, groove, and ring geometries and positions, as shown in Figure 2. The leakage paths are the gaps between rings and grooves (due to areas that become available during ring motions in grooves), the gaps between the ring face and cylinder bore (due to ring lift), and the ring-end gaps. The blowby is assumed to be an unsteady adiabatic flow satisfying the perfect gas law. It is also assumed that the combustion chamber pressure versus crank angle relation remains unaffected despite the gas leakage. The continuity equations for the land and groove gas volumes are written as where is the mass of gas in volume anḋdenotes the mass flow rate through the various flow paths identified by subscripts "leg"/"geg" (ring-end gaps connected to groove volume and connected to land volume), "ggt"/"ggb" (ring groove gaps above and below ring), and "rlg" (ring-left gaps). For the first land volume, the continuity equation is wherėl in,1 is the mass flow rate through the gap between piston and bore at the piston top position. The gas flow rates through the flow paths are calculated by the following relationship.
The gas equation of state is applied to land volumes and the groove volumes: In this paper, the whole system includes three piston rings and five gas zones that are formed by rings, piston, and liner as sketched in Figure 3. Gas temperatures in the land volumes and groove volumes are assumed to be equal to the instantaneous area-weighted average of the temperature of surfaces surrounding each volume. And those temperatures and volumes would be gotten by the thermomechanical analysis of piston. In this study, an efficient method utilizing the concept of inverse heat conduction based on the FEA is presented for the thermal analysis of piston.

Thermomechanical Analysis of Piston Based on Using the Inverse Heat Conduction Method
According to the structure of piston head and bore, shown in Figure 3, in order to get the gas temperatures in land volumes and groove volumes, it is important to know the temperature distribution of piston head and bore. In recent years, finite element method has been used to calculate piston temperature [21][22][23][24][25][26][27][28][29]. It should be noticed that after having identified one correlation for one engine, we acknowledge that those correlation parameters are not valid for another engine with similar characteristics. This is because the convection heat transfer coefficient has many influence factors, and it is hard to give an accurate correlation to describe the distribution of HTC. In this paper, an efficient method utilizing the concept of inverse heat conduction based on the FEA is presented for the thermal analysis of pistons [27]. And using the inverse heat conduction method, the HTC can be effectively determined without giving the correlation.

Inverse Heat Conduction
Method. The relationship between temperature and HTC (ℎ ) can be given by Thus, the inverse heat conduction problem can be described to find where is the real temperature (testing temperature), related to (ℎ ) (predicted temperature) and is the number of conduction segments. Due to the specific piston, the structure of this piston is known, and the HTC is the only unknown parameter. According to (11), if the error between calculated temperature and test temperature can be satisfied, the calculation requirement, the HTC ℎ can be tried as the real HTC during conduction in real situation. The governing equation of the steady state thermal conduction in the piston can be expressed as [27] where and represent the temperature and thermal conductivity, respectively. Similar to earlier studies on the thermal analysis of the piston [30], the convection boundary condition can be applied to the piston surface according to the equation as follows: where , ℎ, and ∞ represent the outward normal heat flux, the HTC, and the surrounding temperature at the surface, respectively. The surrounding temperature ∞ can be defined as the combustion gas temperature gas at the combustion-side top surface of the piston head, the cooling Mathematical Problems in Engineering water temperature at the side, and the oil temperature oil at the underside and cooling gallery [31]. In order to find an accurate solution to the inverse heat conduction problem, the finite element (FE) method is employed in this study. The 2D FE model of piston is built by quadratic solid 55 elements in order to analyze the axisymmetric problem (as shown in Figure 4(a)). And the head of piston is divided into numerous boundary segments ( = 23) of which HTCs are a piecewise constant (as shown in Figure 4(b)). The distribution of testing point is shown in Figure 4(c). Assume the temperature of node is (ℎ ), and the initial magnitude is ℎ 0 = (ℎ 10 , ℎ 20 , . . . , ℎ 0 ). Thus using Taylor's formula, (ℎ ) can be expressed as If (ℎ ) = , (14) can be written as Thus, those equations can compose a linear equation system relating to (ℎ 1 , ℎ 2 , . . . , ℎ ), and it can be written as where the HTC matrix h c = (ℎ 1 , ℎ 2 , . . . , ℎ ) , the equation , and each element can be calculated by And = { 1 , 2 , . . . , }; vector element in this matrix can be given by If the partial derivative of each HTC is known, the solution is the optimum solution. However, it is hard to know the linear correlation between HTC and temperature and the following equation can be used to get the better approximate solution: where { } is the predicted temperature which is unknown and [ ] is the function of HTC and heat conductivity coefficient (HCC); it can be given as = ( , ℎ ). There are a number of partial derivatives in (19). In order to improve the solving accuracy, the first-order method of the ANSYS design optimization module is used in this study.

Numerical Simulation Procedure.
The inverse heat conduction problem of the equation above can be solved using the optimization technique; the procedure to solve the problem is described in Figure 5. The HTC varies depending on the location of the piston surface. The HTC at the surface of head, excluding the piston head, is determined by Lu et al. [31], and the initial value of all design variables, that is, the HTC, is set uniformly as 100 W/m 2 K. An inverse heat conduction problem to determine ℎ can be formulated as an optimization problem as finding [32] ℎ ≡ {ℎ 1 , . . . , ℎ } (20) to minimize The numerical implementation to solve the optimization problem is performed by ANSYS. The first-order method of the ANSYS design optimization module is used, the gradients being calculated with a design variable increment of 0.05 percent of the difference between the upper and lower

Temperature Distribution and Thermal Deformation of
Piston Lands and Grooves. The optimization problem for the piston is analyzed. Less than 9 iterations are taken to obtain the optimum solutions; the history of the objective function and testing points' temperature during iterations is shown in Figure 6. The initial optimum value of the objective function is very large and the optimum value of objective function is very close to zero. As an inverse problem, the temperatures of testing points become close to the test value (shown in Figure 6) very fast. And the temperature distributions calculated by initial and optimum value of the HTC are shown in Figure 7. The analyzed temperature using the optimum HTC is compared with the measured temperature in

Performance of Ring
Pack. The geometry of the piston rings and structure of piston head used in this study are shown in Figures 8 and 9. Each ring is modeled as a single mass; twisting (including pretwist angles) is considered. The equations of motion, which consider equilibrium condition of moments and forces for each ring, are solved. The dynamic components of the ring motion are calculated by means of time integration methods of explicit type. The mechanism of piston ring sealing is equivalent to a labyrinth seal, where the gap clearances are determined by the actual position of the rings in the groove in consideration of the global movement and tilting of the piston [33]. The law for isotropic flow of ideal gas is used for the calculation of the gas flows, as (9) showed. The calculation of the updated gas pressures due  to changes in mass and vessel volume is described by the isotherm equation of state. The particular piston temperature and thermal deformation are given by FEA, listed in Tables  2 and 3. The Reynolds equations are solved iteratively in each time step to determine the hydrodynamic pressure distribution between ring running surface and liner. The last part of the simulation procedure deals with the evaluation of oil consumption to the combustion gas [34] and the friction loss [35,36]. All relevant mechanisms are considered; the calculation is done quasistatically per time step, and the program code is designed by Fortran.

The Pressure Distribution in Ring Grooves and Piston
Lands. The calculated axial position of rings in grooves is  shown in Figure 10. Wear of the parallel surfaces in piston ring grooves, commonly called ring groove wear, occurs mainly in the top groove. The main reason for the wear is the combined effect of gas forces and radial motion of the ring, and the wear process is accelerated by poor lubrication and a high temperature. The reasons for the radial motion of the ring are the cylinder distortion, the secondary movement of the piston, and piston tilt allowed by the piston/cylinder clearance. Mass forces, friction forces, axial ring movement, and ring rotation increase the ring groove wear. Instationary gas pressure and gas blowby may cause radial vibrations in the ring, which accelerates the ring groove wear at the ring groove contact areas [37]. And the friction loss and gas blowby at the piston rings are summarized in Table 4. It is clear that the ratios of the first ring in friction loss and inverse blowby are majority (57.14% and 97.44%, resp.), and this means the design of the first ring groove is significant for the piston design.

Piston Groove Parameters Sensitive Analysis.
Blow-by gas flow refers to the undesired gas flow from the combustion chamber to the crankcase. It reduces the efficiency of engine and contaminates oil with the combustion products present    in the gases. The ring pack system serves as a channel for blow-by gases. As a result, it is important for engine manufacturers and lubricant suppliers to optimize ring pack system to take it into account in controlling oil consumption resulting from blow-by gases. In order to investigate the effect of groove parameters on blowby, a sensitive analysis based on the analysis of partial regression coefficients and an analysis of variance (ANOVA) table for the simple linear regression model can be given in Table 5 [38]. During the multiple linear regression (MLR) analysis, the blowby of first ring is the dependent variable (DV), and the parameters of groove, such as Δℎ1( ), Δℎ2( ), Δ 1( ), and Δ 2( ) ( is the number of grooves) shown in Table 2, are the independent variables (IVs). In this study, when one IV increased by one unit (1 mm), all the other IVs are held constant (as shown in Table 3). The analysis results of multiple linear regression related to the groove parameters are listed in Tables 6-9, respectively. Based on those tables, it is obvious that Δℎ1(1) and Δℎ2(1) have significant influence on blow-by gas flow and other parameters without any obvious impact on it. Due to the thermal effect on Δℎ1(1) and Δℎ2 (1), the interring gas dynamic analysis should consider the thermal effect.

Conclusions
The present work features detailed interring gas dynamics of piston ring pack behavior in internal combustion engines. The model is developed for a ring pack with four rings. The dynamics of ring pack are simulated. Due to the fact that small changes in geometry of the grooves and lands would have a significant impact on the interring gas dynamics, the thermal deformation of piston and bore has been considered during the ring pack motion analysis in this study. In order to get the temperature distribution of piston head more quickly and accurately, an efficient method utilizing the concept of inverse heat conduction is presented. The result of the temperature analyzed using the optimum HTC is compared  with the measured temperature, and reasonable agreement is obtained. Moreover, a sensitive analysis based on the analysis of partial regression coefficients is presented to investigate the effect of groove parameters on blowby. And the result shows that Δℎ1(1) and Δℎ2(1) (shown in Table 2) have significant influence on blow-by gas flow, while other parameters have no obvious impact on it. As the thermal effect on Δℎ1(1) and Δℎ2 (1), the interring gas dynamic analysis should consider the thermal effect.