Instantaneous in-Cylinder Volume Considering Deformation and Clearance due to Lubricating Film in Reciprocating Internal Combustion Engines

A new methodology for predicting the real instantaneous in-cylinder volume in the combustion chamber of a reciprocating internal combustion engine is implemented. The mathematical model developed as part of this methodology, takes into consideration the deformations due to pressure and inertial forces, via a deformation constant adjusted through ANSYS®, using a high-precision CAD model of a SOKAN SK-MDF300 engine. The deformation constant was obtained from the CAD model using the computational tool ANSYS® and the pressure data was obtained from the engine running at three regimes: 1500, 2500, and 3500 rpm. The results were compared with previous models reported in the literature, showing that the deformation constant obtained has a smaller variation among cycles, which leads to a more precise value of the mechanical deformations. Furthermore, to have a more accurate model of the instantaneous volume variation, a factor taking into consideration the lubricant film behavior is introduced to calculate volumetric variation due to geometrical clearances. The influence of the introduced volumetric variation was evaluated through a process of combustion diagnosis, evidencing the improvement in the predictive capacity of thermodynamic modeling and, therefore, the correct prediction of heat release rate.


Introduction
Internal combustion engines have been studied since their creation and because of their vast number of applications, it is indispensable to have a better understanding of the processes occurring in these machines [1,2].Currently, for combustion analysis, there are computational approaches, such as predictive models [3][4][5][6][7][8] and diagnostic models [9][10][11][12][13], in which the in-cylinder pressure signal is taken, and through the first law of thermodynamics, the rate of heat release is determined.
The input and output data of the model varies considerably depending on the purpose of the model itself.To know the combustion process, a diagnostic model must be used [14].Such a model is intended to obtain the heat release curve inside the combustion chamber.To attain it, this type of modeling requires the chamber pressure during the entire cycle as an input.However, if the chamber pressure is required, it is necessary to use a mirror model to the diagnostic model [14], which is a predictive combustion model.It should be noted that the correct operation of the models described above, requires a previous characterization of the engine and adequate data for the submodels, i.e., the heat transfer submodel.
The engine is characterized using a technique called motoring test, in which the engine is motored to certain engine speed by an electric motor.Under this test, no fuel is injected into the engine.Therefore, no combustion is generated in the chamber, and the heat release curve must be zero [15].One of the models adjusted during the test corresponds to the instantaneous volume of the combustion chamber, which in actual conditions, presents variations respect to the value calculated geometrically.Failure to consider these real variations would generate large differences when estimating the heat release rate [15].Previous works [15] focussed on quantifying the volumetric variations due to the combustion chamber pressure and the inertial loads acting on the piston-rod-crankshaft mechanism, using equivalent approximations.As a consequence of this approach, the volumetric variation is defined as a function of a deformation constant K def , whose value is specific for each engine studied, and should be adjusted based on experimental data obtained from motoring tests.This adjusting process is time-and computational resources consuming due to its iterative nature.Furthermore, there is a high margin of uncertainty when determining coefficient, with variations up to 70% in its magnitude even for the same engine tests.To solve this situation, authors have proposed different methods: (a) the use of alternative ways to characterize the deformations and thus the instantaneous volume [16], and (b) the development of new methodologies to determine the uncertainties of the adjustment constants [12].
Over the years, many publications have developed models to study and analyze internal combustion engines.Sun et al. [17] considered the crankshaft deformations and the influence of the pressure exerted by the oil lubrication film on the bearing surface.They found the influence exerted by the cylinder deformation over the lubrication film thickness, the end leakage flow rate, and the frictional coefficient.Furthermore, the study demonstrates the important role of the lubricating film thickness over the cylinder deformation.Sun et al. [18], applied finite element analysis to study the effect of the lubricating film over the crankshaft stress distribution.
Weinzapfel et al. [19] studied the fatigue life dispersion of a crankshaft, specifically at the running line, considering the material microstructure.Cai et al. [20] analyzed the effect of the dry/wet contact in elastic/plastic materials with rough surfaces, finding the optimum layer parameters.Mo et al. [21] presented a thermodynamic analysis of a self-lubricant bearing developed in Fluent ® for gear pumps, finding the whirling orbit behavior of the temperature.Regarding deformations in internal combustion engines, it was noticed the lubricant film effect needs to be addressed [22,23] because it induces changes in the free height of the piston, hence in the instantaneous volume throughout the thermodynamic cycle [5].
The main objective of this article is the development of a new methodology for the prediction of real instantaneous volume in the combustion chamber of an internal combustion engine.The model will be implemented in ANSYS, to reduce the computational time as well as the computing resources requirements.The model to be developed considers the effects of mechanical deformations by pressure and inertial loads on the instantaneous volume and proposes a submodel for the calculation of volume variation due to changes in the lubricating film, and its effect on the instantaneous global volume.

Model Considerations
The work developed by Duarte et al. [24] was taken into consideration to develop the deformation model of an internal combustion engine.The main considerations are the following:

•
The deformation occurs only in the connecting rod [16].

•
The mass center of the crankshaft does not experience translational movement [16,17].

•
The connecting rod is replaced by a concentrated mass system [16,17] as shown in Figure 1, with values calculated from Equations ( 1) to (3), Therefore, the inertial mass of the system m i is obtained from Energies 2019, 12, x FOR PEER REVIEW 3 of 20 Therefore, the inertial mass of the system  is obtained from Figure 1.Statically equivalent mass system for connecting rod.


The acceleration of the deformation model due to inertial forces, as well as the acceleration of the piston, will be equal [15].


The elastic modulus of the connecting rod is considered as unknown [16,24], but it is related to the standard value of steel elastic modulus [210 GPa] through a deformation constant Kdef, as shown in Equation ( 4),  Deformations due to thermal effect are neglected [16].Besides the previous considerations, for the volume variation due to geometrical clearances, it was assumed that the lubricating film obeys the theory of hydrodynamic lubrication [17].Therefore, a laminar flow for the film was assumed, the viscosity of the fluid is considered constant, and the inertial effects, as well as the direct contact between components, are ignored, leaving only the hydrodynamic behavior of the fluid.Further, to mitigate the thermal effects over the engine pieces, the studied engine has a thermal control system that allows a constant lubricant and refrigerant temperature fluid.In addition to this, all tests done during the development of this work were carried out under steady conditions (constant RPM) during long periods [15].


Vibrational effects in lubrication oil film thickness are neglected [15,25] due to the stabilization time allowed during the tests.This allows that the big initial displacements generated by the engine torsional vibrations were mitigated, creating minimal displacements that can be neglected.

Kinematic Analysis of the Mechanism
A kinematic analysis of the system (piston, connecting rod, and crankshaft pin) is required.The system is considered as a slider-crank mechanism [25].A schematic of the system is shown in Figure 2, in which Element 1 is ground (or the non-moving part of the engine), Element 2 is the crankshaft, and Elements 3 and 4 are the connecting rod and the engine piston, respectively.
Using complex algebra, the velocity and acceleration of the piston in the decentralized slider-

•
The acceleration of the deformation model due to inertial forces, as well as the acceleration of the piston, will be equal [15].

•
The elastic modulus of the connecting rod is considered as unknown [16,24], but it is related to the standard value of steel elastic modulus [210 GPa] through a deformation constant K def , as shown in Equation ( 4), • Deformations due to thermal effect are neglected [16].Besides the previous considerations, for the volume variation due to geometrical clearances, it was assumed that the lubricating film obeys the theory of hydrodynamic lubrication [17].Therefore, a laminar flow for the film was assumed, the viscosity of the fluid is considered constant, and the inertial effects, as well as the direct contact between components, are ignored, leaving only the hydrodynamic behavior of the fluid.Further, to mitigate the thermal effects over the engine pieces, the studied engine has a thermal control system that allows a constant lubricant and refrigerant temperature fluid.In addition to this, all tests done during the development of this work were carried out under steady conditions (constant RPM) during long periods [15].

•
Vibrational effects in lubrication oil film thickness are neglected [15,25] due to the stabilization time allowed during the tests.This allows that the big initial displacements generated by the engine torsional vibrations were mitigated, creating minimal displacements that can be neglected.

Kinematic Analysis of the Mechanism
A kinematic analysis of the system (piston, connecting rod, and crankshaft pin) is required.The system is considered as a slider-crank mechanism [25].A schematic of the system is shown in where  and  are the lengths of the crankshaft arm and connecting rod respectively, while  =  ̇ is the piston velocity with respect to the reference coordinate system.Similarly,  =  ̇ and  / =  ̇, or the derivative of the angle between the connecting rod and the piston displacement axis.This angle can be calculated by the equation By solving this system given in Equations ( 5) and ( 6), the speed and angular acceleration of the connecting rod are obtained, together with the velocity and acceleration of the piston, for each angular position of the mechanism.

Kinetic Analysis of the Mechanism
To determine the forces applied on joints between components (which are necessary for the mathematical clearance model), the system is described as a crank-sliding mechanism, in which each part of the system will be analyzed separately (as shown in Figure 3).In the end, the expressions obtained will be unified in systems of linear equations in order to find the required unknowns.The crankshaft analysis leads to the following Equations Using complex algebra, the velocity and acceleration of the piston in the decentralized slider-crank mechanism can be calculated by solving the Equations [25]: where L 2 and L 3 are the lengths of the crankshaft arm and connecting rod respectively, while v 4 = .R 4y is the piston velocity with respect to the reference coordinate system.Similarly, ω 2 = .
β, or the derivative of the angle between the connecting rod and the piston displacement axis.This angle can be calculated by the equation L 3 (7) By solving this system given in Equations ( 5) and ( 6), the speed and angular acceleration of the connecting rod are obtained, together with the velocity and acceleration of the piston, for each angular position of the mechanism.

Kinetic Analysis of the Mechanism
To determine the forces applied on joints between components (which are necessary for the mathematical clearance model), the system is described as a crank-sliding mechanism, in which each part of the system will be analyzed separately (as shown in Figure 3).In the end, the expressions obtained will be unified in systems of linear equations in order to find the required unknowns.
To determine the forces applied on joints between components (which are necessary for the mathematical clearance model), the system is described as a crank-sliding mechanism, in which each part of the system will be analyzed separately (as shown in Figure 3).In the end, the expressions obtained will be unified in systems of linear equations in order to find the required unknowns.The crankshaft analysis leads to the following Equations The crankshaft analysis leads to the following Equations For the connecting rod, the analysis generates the following Equations While for the piston, taking into account the pressure generated by the gas p gas , the analysis leads to the following Equations By coupling ( 8) to (15), the following matrix is obtained

Instantaneous Geometric Volume Calculation
The volume in the cylinder, when the effects of the instantaneous variations are not considered, is composed of two terms [16]:

•
The minimum combustion chamber volume, which is the free space remaining when the piston reaches the top dead center (TDC).This is not dependent on the angular position of the crankshaft: where r c is the engine compression ratio.

•
The volume displaced by the stroke of the piston in its path, which depends on the angular position of the crankshaft: For centralized geometries (d = 0), R 4y,max = L 2 + L 3 , and (18) becomes: Besides these terms, the pressure of the gases and the inertial forces of the mechanism, generate a variation in the instantaneous volume of the combustion chamber.This variation was studied by previously [25], and it is calculated by (20), with A CR being the area of the slender Section of the connecting rod.

Mathematical Model of Clearance
In its most basic form, a sliding bearing consists of two elements: a cylindrical shaft, a shaft or a journal to which a load and a rotational speed are applied, and a cylindrical sleeve or bearing, which supports the journal as well as the load exerted.When there is no relative motion between these elements or an applied load (Figure 4a), these two elements are concentric to each other, separated by a radial clearance C [26].If a load is applied, or there is a relative movement, these elements are no longer completely concentric.Because of the lubricant film behavior, there is an eccentricity e along the imaginary line joining the centers of the journal and bearing (Figure 4b).This is commonly expressed through the eccentricity ratio ε, which corresponds to the dimensionless relationship between eccentricity and radial clearance.Being defined in this way, the eccentricity ratio can only take values between zero and one.
Energies 2019, 12, x FOR PEER REVIEW 6 of 20 with  being the area of the slender Section of the connecting rod.

Mathematical Model of Clearance
In its most basic form, a sliding bearing consists of two elements: a cylindrical shaft, a shaft or a journal to which a load and a rotational speed are applied, and a cylindrical sleeve or bearing, which supports the journal as well as the load exerted.When there is no relative motion between these elements or an applied load (Figure 4a), these two elements are concentric to each other, separated by a radial clearance C [26].If a load is applied, or there is a relative movement, these elements are no longer completely concentric.Because of the lubricant film behavior, there is an eccentricity e along the imaginary line joining the centers of the journal and bearing (Figure 4b).This is commonly expressed through the eccentricity ratio ε, which corresponds to the dimensionless relationship between eccentricity and radial clearance.Being defined in this way, the eccentricity ratio can only take values between zero and one.
The origin of volume variation due to clearance lies in the behavior of the lubricating film.Therefore, the mathematical modeling of the previous phenomenon is developed from the Reynolds lubrication equation.For pure tangential sliding motion, pressure changes within the lubricant film due to the physical wedge mechanism.In this case, the flow perpendicular to the direction of motion (side-leakage effect) has to be included, and the Reynolds equation takes the form [27]: The origin of volume variation due to clearance lies in the behavior of the lubricating film.Therefore, the mathematical modeling of the previous phenomenon is developed from the Reynolds lubrication equation.For pure tangential sliding motion, pressure changes within the lubricant film due to the physical wedge mechanism.In this case, the flow perpendicular to the direction of motion (side-leakage effect) has to be included, and the Reynolds equation takes the form [27]: where h corresponds to the film height, p to the pressure exerted on the lubricant, D j to the diameter of the bearing, ω j is the angular velocity of the hub, and µ is the dynamic viscosity of the lubricant.Due to the complexity of the solution of ( 22), Vignolo et al. [28] validated an approximate analytical solution developed by Dubois [29], that takes care of side leakage and assumes that the pressure-induced flow in the circumferential direction is small relative to the Couette flow term.This assumption further implies that the Poiseuille flow is less significant in the circumferential direction.This short bearing theory is valid as long as the diameter-to-width ratio is greater than 2, which is fulfilled in this case.
Because of this assumption, the Reynolds equation reduces to: To calculate the film height h, a first approximation is obtained from a geometrical analysis of the bearing with an eccentricity e, as shown in Figure 5.As a result, the film height expression has the following structure [26]: Energies 2019, 12, x FOR PEER REVIEW 7 of 20 Equation ( 24) fulfills the objective of providing an expression for the film height.Nevertheless, it is ambiguous and difficult to implement.Because of this, trigonometric relations from ( 25) are used to obtain (26), In the literature, it is common to find the mathematical modeling of the parameters necessary for the design of a radial short bearing, such as the load capacity that the film can withstand without breaking, and the position of the locus with the smallest film size (which corresponds to the position of the center line).However, these parameters are described under the condition of zero motion for the bearing.That is, the journal is subjected to a defined rotational speed, while the former is held fixed [26].In the event that the above is not met, it is possible to replace the original system-bearing with rotational speed ωj and static journal-by an equivalent system in which the bearing is fixed, and the journal has an opposite rotational speed to the original, in order to maintain relative rotation between the two elements.For this, an additional coordinate system is applied as shown in Figure 6, which has an orientation β concerning the original system.Since this system changes its orientation Equation ( 24) fulfills the objective of providing an expression for the film height.Nevertheless, it is ambiguous and difficult to implement.Because of this, trigonometric relations from ( 25) are used to obtain (26), In the literature, it is common to find the mathematical modeling of the parameters necessary for the design of a radial short bearing, such as the load capacity that the film can withstand without breaking, and the position of the locus with the smallest film size (which corresponds to the position of the center line).However, these parameters are described under the condition of zero motion for the bearing.That is, the journal is subjected to a defined rotational speed, while the former is held fixed [26].In the event that the above is not met, it is possible to replace the original system-bearing with rotational speed ω j and static journal-by an equivalent system in which the bearing is fixed, and the journal has an opposite rotational speed to the original, in order to maintain relative rotation between the two elements.For this, an additional coordinate system is applied as shown in Figure 6, which has an orientation β concerning the original system.Since this system changes its orientation at each instant, it becomes necessary to transform the resulting loads considering the fixed coordinate system to apply the following mathematical analysis.From trigonometric relations, loads of the new system are expressed by [25,29]: From these components, the magnitude of the resulting load and its orientation is known by: Once the above is obtained, the relevant parameters for the free height variation are determined, which are the eccentricity value e (or the eccentricity relation ε), and the location angle  .These variables are represented graphically in Figure 6b.Using the previous complex mathematical analysis, it is shown that the load capacity for a shortwidth journal bearing, is related to the eccentricity ratio, by the expression [25,26]: Assuming a static equilibrium between the journal load and the reaction exerted by the oil film, the value of ε and the angle of inclination of the center line is estimated, which is measured from the line of application of the load, in the direction of rotation of the journal.This is calculated by: Using the definition of lubricant film pressure and circumferential velocity U j , the axis components of the load W j are obtained by the following: Energies 2019, 12, 1437 9 of 21 From these components, the magnitude of the resulting load and its orientation is known by: Once the above is obtained, the relevant parameters for the free height variation are determined, which are the eccentricity value e (or the eccentricity relation ε), and the location angle ζ j .These variables are represented graphically in Figure 6b.
Using the previous complex mathematical analysis, it is shown that the load capacity for a short-width journal bearing, is related to the eccentricity ratio, by the expression [25,26]: Assuming a static equilibrium between the journal load and the reaction exerted by the oil film, the value of ε and the angle of inclination of the center line is estimated, which is measured from the line of application of the load, in the direction of rotation of the journal.This is calculated by: With these two values, it is possible to determine the angle that forms the center line with the x' axis, through two relations: From this, the free height clearance variation is expressed by: However, this is the variation of height related to the coordinate system located along the connecting rod.To bring this to the original coordinate system, the following expression is used: Due to the number of joints present in the mechanism, the individual contribution of each of them over the volume variation should be considered.For this reason, the variation of volume by clearances is expressed by: Energies 2019, 12, 1437 10 of 21

Determination of Instantaneous Volume
From Equations ( 17), ( 18), ( 19) and ( 20), a unique expression can be formulated to describe the cylinder volume, for each angular position of the crankshaft.In this expression, the volume variation due to deformations by pressure and inertia are taken into account, and it is expressed mathematically by: where the constant of mechanical deformations K def , corresponds to a unique value for each engine, which must be adjusted from experimental data.Including the effect of clearances (Equation ( 44)), the above expression becomes: This expression calculates the instantaneous volume, considering all the effects studied (pressure, inertia, and clearance).

Result Methodology
The methodology used covers a series of stages.Initially, a treatment of the pressure data from the combustion chamber measured in the SOKAN SK-MDF300 diesel direct injection engine whose geometric characteristics can be consulted in Table 1, was carried out through the fast Fourier transform (FFT) with a cut-off frequency of 44 kHz.This is an accurate value for the type of analysis developed in this work.This pressure data was measured for three-rotational speed regimes: 1500, 2500, and 3500 rpm, which are common operation values for the studied engine.The filtered pressure data was used for the simulation of the engine in ANSYS ® , from which the piston free height variation due to pressure and inertia deformation was obtained.To do so, the mesh for the connecting rod was generated by the tetra patch conforming method, with a degree of refinement (relevance) of 100, generating tetrahedral elements with an average asymmetry coefficient (skewness) of 0.21.Further, the frictional contact for the simulation was obtained from [30,31], which includes shear stress and viscous effects.The mechanical data of the simulated elements are summarized in Table 2. To verify that the results obtained in the simulation were consistent with experimental data, validation was performed on the universal testing machine, whose characteristics are illustrated in Table 3.The obtained results were compared with the developed model, without taking into account the volume variation due to clearances, as shown in Equation (45), and allowed the validation of the proposed model and the calculation of the mechanical deformation coefficient, K def .Likewise, the volume variation due to clearances was obtained, and its relevance in the instantaneous volume change was emphasized.A sensitivity analysis was carried out using a combustion diagnosis model [15,16].Using the in-cylinder pressure signal, the rate of heat release (ROHR) and the temperature in the combustion chamber are obtained.This will emphasize the importance of the correct prediction of the volume in the combustion chamber and will help to study the influence of the volume variation, which is the main objective of the present research.
The results obtained in the universal test machine were compared with the simulation results in Table 4, validating that the simulations adequately describe the phenomena under study.The validation was carried out at the top dead center, because at this point the load acts at compression, which is a loading mode in which the piston rod assembly can be tested in the universal machine, due to limitations in the fixing of the elements to the machine.This validates the results obtained from the simulation.The variations can be attributed to changes in the simulated geometry used in ANSYS ® , and that affects the deformation values obtained.Likewise, there is uncertainty between the materials simulated and the materials used in the real power train.Despite that, the maximum error in the prediction does not exceed 2%.
In Section 2.6, it was mentioned that pressure in the combustion chamber is the starting point for the calculation of volume variations due to clearances, as the load profile for journal bearings (a mandatory input for the estimation of lubricant film thickness) varies according to the measured pressure in the engine.Given the existence of previous works that evaluate the behavior of lubricant films in conditions similar to that described, the validation was performed using the experimental load profiles measured by Moreau [32] (as shown in Figure 7).These values were used as input in the model and the minimum film thickness was obtained for each angular position of the crankshaft.This procedure is summarized graphically in Figure 8.
between the materials simulated and the materials used in the real power train.Despite that, the maximum error in the prediction does not exceed 2%.
In Section 2.6, it was mentioned that pressure in the combustion chamber is the starting point for the calculation of volume variations due to clearances, as the load profile for journal bearings (a mandatory input for the estimation of lubricant film thickness) varies according to the measured pressure in the engine.Given the existence of previous works that evaluate the behavior of lubricant films in conditions similar to that described, the validation was performed using the experimental load profiles measured by Moreau [32] (as shown in Figure 7).These values were used as input in the model and the minimum film thickness was obtained for each angular position of the crankshaft.This procedure is summarized graphically in Figure 8.In Figure 9, the values calculated by the proposed model were compared with the experimental results obtained by Moreau et al. [32], where deviations are less than 5% when compared to data sets.This graph allowed to verify the capacity of the model to predict the behavior of minimum film thickness, and, therefore, the journal eccentricity originated by the operating conditions.Since Equation (42) establishes the relation between the eccentricity and the instantaneous volume variation in the combustion chamber, it can be inferred that the proposed model can produce a robust prediction of the film effects on the volume changes due to clearances.In Figure 9, the values calculated by the proposed model were compared with the experimental results obtained by Moreau et al. [32], deviations are less than 5% when compared to data sets.This graph allowed to verify the capacity of the model to predict the behavior of minimum film thickness, and, therefore, the journal eccentricity originated by the operating conditions.Since Equation (42) establishes the relation between the eccentricity and the instantaneous volume variation in the combustion chamber, it can be inferred that the proposed model can produce a robust prediction of the film effects on the volume changes due to clearances.results obtained by Moreau et al. [32], where deviations are less than 5% when compared to data sets.This graph allowed to verify the capacity of the model to predict the behavior of minimum film thickness, and, therefore, the journal eccentricity originated by the operating conditions.Since Equation (42) establishes the relation between the eccentricity and the instantaneous volume variation in the combustion chamber, it can be inferred that the proposed model can produce a robust prediction of the film effects on the volume changes due to clearances.

Results and Discussion
As defined above, Equation ( 20) is expressed as a function of an adjustment coefficient Kdef that takes into account the uncertainties of the proposed model and is adjusted so that the rate of heat release is zero in the motored engine, because there is no combustion.Therefore, if the coefficient is not adjusted properly, it causes an error in the variation of instantaneous volume.To estimate the value of the deformation coefficients, the volume variation obtained in ANSYS ® was compared with the calculated value for each angular position of the crankshaft.These results were obtained for each regimen, and after statistical analysis, the mean and standard deviation of the coefficient for each engine speed was obtained.Table 5 shows the values for these statistical parameters.

Results and Discussion
As defined above, Equation ( 20) is expressed as a function of an adjustment coefficient K def that takes into account the uncertainties of the proposed model and is adjusted so that the rate of heat release is zero in the motored engine, because there is no combustion.Therefore, if the coefficient is not adjusted properly, it causes an error in the variation of instantaneous volume.To estimate the value of the deformation coefficients, the volume variation obtained in ANSYS ® was compared with the calculated value for each angular position of the crankshaft.These results were obtained for each regimen, and after statistical analysis, the mean and standard deviation of the coefficient for each engine speed was obtained.Table 5 shows the values for these statistical parameters.This shows that the values obtained by the proposed model meet the homoscedasticity condition, that is, its mean and standard deviation remain almost constant throughout the studied regimes.Additionally, it was noticed that these values are remarkably close to the ideal case.That is, K def is equal to one.From the results of Table 5, the global value of 1.15 was assigned as the coefficient of deformations of the studied engine, which produces adequate results when the adjusted model is compared with the ANSYS ® results, as shown in Figures 10-12 for 1500, 2500, and 3500 rpm.
This shows that the values obtained by the proposed model meet the homoscedasticity condition, that is, its mean and standard deviation remain almost constant throughout the studied regimes.Additionally, it was noticed that these values are remarkably close to the ideal case.That is, Kdef is equal to one.From the results of Table 5, the global value of 1.15 was assigned as the coefficient of deformations of the studied engine, which produces adequate results when the adjusted model is compared with the ANSYS ® results, as shown in Figures 10-12 for 1500, 2500, and 3500 rpm.This shows that the values obtained by the proposed model meet the homoscedasticity condition, that is, its mean and standard deviation remain almost constant throughout the studied regimes.Additionally, it was noticed that these values are remarkably close to the ideal case.That is, Kdef is equal to one.From the results of Table 5, the global value of 1.15 was assigned as the coefficient of deformations of the studied engine, which produces adequate results when the adjusted model is compared with the ANSYS ® results, as shown in Figures 10-12 for 1500, 2500, and 3500 rpm.To demonstrate that the clearance plays a significant role in the instantaneous volume throughout the working cycle, Equation (44) was used, with the help of Octave ® , to obtain the volume variation due to the lubricating film.The clearance model described before was applied to two points of interest: the piston pin-connecting rod and the crankshaft junction, for each regimen, studied (1500, 2500, and 3500 rpm), using the properties of SAE 15W-40 lubricant oil.The geometric properties of To demonstrate that the clearance plays a significant role in the instantaneous volume throughout the working cycle, Equation (44) was used, with the help of Octave ® , to obtain the volume variation due to the lubricating film.The clearance model described before was applied to two points of interest: the piston pin-connecting rod and the crankshaft junction, for each regimen, studied (1500, 2500, and 3500 rpm), using the properties of SAE 15W-40 lubricant oil.The geometric properties of these joints are summarized in Table 6, and the results obtained in Octave ® are shown in Figure 13.To demonstrate that the clearance plays a significant role in the instantaneous volume throughout the working cycle, Equation (44) was used, with the help of Octave ® , to obtain the volume variation due to the lubricating film.The clearance model described before was applied to two points of interest: the piston pin-connecting rod and the crankshaft junction, for each regimen, studied (1500, 2500, and 3500 rpm), using the properties of SAE 15W-40 lubricant oil.The geometric properties of these joints are summarized in Table 6, and the results obtained in Octave ® are shown in Figure 13.Table 6.Geometric characteristics of the joints studied.

Dimension
Piston pin-Connecting Rod Connecting Rod-Crankshaft Length By direct comparison of volume variation due to pressure, inertia and clearances (shown in Figure 14 at 1500 rpm), it is noticed that the volume variation due to the latter represents By direct comparison of volume variation due to pressure, inertia and clearances (shown in Figure 14 at 1500 rpm), it is noticed that the volume variation due to the latter represents approximately 66% of the range comprised by deformations due to pressure and inertia.Furthermore, it was found that clearances model presents an almost constant variation through the thermodynamic cycle, unlike the variations due to pressure and inertia, which are manifested largely near the TDC (approximately about 120 • of rotation of the crankshaft).As a result of the above, it is assumed that the behavior of the lubricating film is a factor required to calculate the overall variation of the free height of the piston, and hence the instantaneous volume variation in the combustion chamber.
Furthermore, it was found that clearances model presents an almost constant variation through the thermodynamic cycle, unlike the variations due to pressure and inertia, which are manifested largely near the TDC (approximately about 120° of rotation of the crankshaft).As a result of the above, it is assumed that the behavior of the lubricating film is a factor required to calculate the overall variation of the free height of the piston, and hence the instantaneous volume variation in the combustion chamber.Taking into account previous research [15,[33][34][35] and the fact that there is a relation between the volume of the combustion chamber and the inside temperature reached, it is evident the influence of the volume variation due to pressure and clearances on the thermodynamic modeling [36].Additionally, diagnostics thermodynamic modeling shows variations in the rate of heat release due to changes in the compression ratio, which in turn is caused by the change of volume in the combustion chamber.Therefore, this investigation shows the influence of volume variation on thermodynamic modeling, using the DIAGNO-DIESEL ® software for 1500, 2500, and 3500 rpm.These angular velocities were chosen to take into account the entire range of operation speed for the engine under study.
After engine characterization, the average temperature in the chamber and rate of heat release (ROHR) are analyzed, considering the crankshaft angle.Figure 15 shows the influence of the volume variation of the combustion chamber on the thermodynamic modeling for the studied regimes, where variations up to 4% in the mean temperature in the chamber are observed.In the same way, the variation in the ROHR for different engine speeds was calculated (and shown in Figure 15) and it is emphasized that the best result is obtained when the model including deformation and clearance in the actual research is used.Taking into account previous research [15,[33][34][35] and the fact that there is a relation between the volume of the combustion chamber and the inside temperature reached, it is evident the influence of the volume variation due to pressure and clearances on the thermodynamic modeling [36].Additionally, diagnostics thermodynamic modeling shows variations in the rate of heat release due to changes in the compression ratio, which in turn is caused by the change of volume in the combustion chamber.Therefore, this investigation shows the influence of volume variation on thermodynamic modeling, using the DIAGNO-DIESEL ® software for 1500, 2500, and 3500 rpm.These angular velocities were chosen to take into account the entire range of operation speed for the engine under study.
After engine characterization, the average temperature in the chamber and rate of heat release (ROHR) are analyzed, considering the crankshaft angle.Figure 15 shows the influence of the volume variation of the combustion chamber on the thermodynamic modeling for the studied regimes, where variations up to 4% in the mean temperature in the chamber are observed.In the same way, the variation in the ROHR for different engine speeds was calculated (and shown in Figure 15) and it is emphasized that the best result is obtained when the model including deformation and clearance proposed in the actual research is used.

Conclusions
The application of the filtering methodology to the supplied pressure data allowed obtaining reliable information to be used in the simulation because it permits the elimination of low-frequency noises (by averaging a certain amount of thermodynamic cycles) and high frequency (by the application of the fast Fourier transform) presented in the original signal.Moreover, its programming simplicity, in tools like Octave ® , makes it a convenient way to isolate the useful information contained in the acquired pressure signal.
On the other hand, the implementation of a simulation package, such as ANSYS ® , using the experimental pressure data in the chamber allowed the direct estimation of the coefficient of mechanical deformations of Equation (45), which was characterized by maintaining a low variability between each studied regimen.This allowed obtaining a unique value Kdef of 1.15 for the studied engine, which produces a reasonable adjustment of the model to the experimental data.From this, the volume variation model due to pressure and inertia is validated as a fast and precise way of characterizing the engine, which leads to the reduction of uncertainties in the process, compared with the iterative process currently used.This leads to an improvement in the calculation of the heat

Conclusions
The application of the filtering methodology to the supplied pressure data allowed obtaining reliable information to be used in the simulation because it permits the elimination of low-frequency noises (by averaging a certain amount of thermodynamic cycles) and high frequency (by the application of the fast Fourier transform) presented in the original signal.Moreover, its programming simplicity, in tools like Octave ® , makes it a convenient way to isolate the useful information contained in the acquired pressure signal.
On the other hand, the implementation of a simulation package, such as ANSYS ® , using the experimental pressure data in the chamber allowed the direct estimation of the coefficient of mechanical deformations of Equation (45), which was characterized by maintaining a low variability between each studied regimen.This allowed obtaining a unique value K def of 1.15 for the studied engine, which produces a reasonable adjustment of the model to the experimental data.From this, the volume variation model due to pressure and inertia is validated as a fast and precise way of characterizing the engine, which leads to the reduction of uncertainties in the process, compared with the iterative process currently used.This leads to an improvement in the calculation of the heat release rate, due to the reduction of uncertainties in one of the sub-models involved in the global balance.
Regarding the volume variation by clearances, the results obtained in Octave ® indicate that the lubricant film has a significant influence on the global variation of instantaneous volume, with variations that represent 66% of the deformation produced by the pressure and inertial load, which causes a significant change in the prediction of the thermodynamic behavior of the engine.For these reasons, it is necessary to include them in thermodynamic modeling and its experimental validation, because of the notable influence on the chamber temperature as seen in Figure 15.
Figure 15 shows an improvement in the prediction of ROHR using the model developed in the present investigation, taking into account that the rate of heat release must be close to zero in the condition of motored test since there is no combustion.The study was limited to this condition to avoid external uncertainties due to the combustion process.When the engine speed was increased, it was verified that the ROHR also increased (Figure 15).However, it is observed that the joint influence of the deformation by pressure/inertia and tolerances is significant since it can generate deviations of up to 300% in the ROHR.
The results show that there is a higher average temperature and ROHR when the volume variation in the combustion chamber as shown by Equation ( 46) is not taken into account.This is because in the previous studies of thermodynamic modeling of diesel engine did not consider volume changes due to clearances.In many cases, the deformation energy due to pressure and inertia forces is not considered, which in the present study was taken into account and caused greater predictability.

Figure 1 .
Figure 1.Statically equivalent mass system for connecting rod.

Figure 2 ,
Figure 2, in which Element 1 is ground (or the non-moving part of the engine), Element 2 is the crankshaft, and Elements 3 and 4 are the connecting rod and the engine piston, respectively.Energies 2019, 12, x FOR PEER REVIEW 4 of 20

Figure 4 .
Figure 4. Schematic of a sliding bearing: (a) with no load applied, and (b) with a load applied to the journal.

Figure 4 .
Figure 4. Schematic of a sliding bearing: (a) with no load applied, and (b) with a load applied to the journal.

Figure 6 .
Figure 6.Graphic representation of the journal bearing parameters: (a) coordinate system transformation angles and (b) resultant load and centerline orientation.

Figure 6 .
Figure 6.Graphic representation of the journal bearing parameters: (a) coordinate system transformation angles and (b) resultant load and centerline orientation.

Figure 7 .
Figure 7. Experimental load profile used for the validation, measured by Moreau et al. at 5500 rpm and full load [32].

Figure 7 .
Figure 7. Experimental load profile used for the validation, measured by Moreau et al. at 5500 rpm and full load [32].Energies 2019, 12, x FOR PEER REVIEW 12 of 20

Figure 8 .
Figure 8. Schematic of the procedure used for the validation of the clearance model.

8 .
Schematic of the procedure used for the validation of the clearance model.

Figure 9 .
Figure 9.Comparison of values obtained for minimum film thickness.

Figure 9 .
Figure 9.Comparison of values obtained for minimum film thickness.

Figure 10 .
Figure 10.Results of ANSYS ® and the adjusted model, at 1500 rpm.

Figure 11 .
Figure 11.Results of ANSYS ® and the adjusted model, at 2500 rpm.

Figure 10 .
Figure 10.Results of ANSYS ® and the adjusted model, at 1500 rpm.

Figure 10 .
Figure 10.Results of ANSYS ® and the adjusted model, at 1500 rpm.

Figure 11 .
Figure 11.Results of ANSYS ® and the adjusted model, at 2500 rpm.Figure 11. Results of ANSYS ® and the adjusted model, at 2500 rpm.

Figure 11 . 20 Figure 12 .
Figure 11.Results of ANSYS ® and the adjusted model, at 2500 rpm.Figure 11. Results of ANSYS ® and the adjusted model, at 2500 rpm.Energies 2019, 12, x FOR PEER REVIEW 14 of 20

Figure 12 .
Figure 12. Results of ANSYS ® and the adjusted model, at 3500 rpm.

Figure 12 .
Figure 12. Results of ANSYS ® and the adjusted model, at 3500 rpm.

Figure 13 .
Figure 13.Volume variation due to clearance for regimens studied.

Figure 13 .
Figure 13.Volume variation due to clearance for regimens studied.

Figure 14 .
Figure 14.Variation of volume by clearance and adjusted model, for 1500 rpm.

Figure 14 .
Figure 14.Variation of volume by clearance and adjusted model, for 1500 rpm.

Figure 15 .
Figure 15.Results obtained for (a, c and e) in-cylinder temperature and (b, d, and f) rate of heat release for 1500 rpm, 2500 rpm, and 3500 rpm, respectively.

Figure 15 .
Figure 15.Results obtained for (a,c,e) in-cylinder temperature and (b,d,f) rate of heat release for 1500 rpm, 2500 rpm, and 3500 rpm, respectively.

Table 2 .
Materials used in the mechanism simulation.

Table 3 .
Characteristics of universal testing machine.

Table 5 .
Values obtained from the coefficient of deformation K def .

Table 6 .
Geometric characteristics of the joints studied.