Next Article in Journal
Influence of Solid Lubricants on the Tribological Performance of Photocurable Resins for Vat Photopolymerization
Previous Article in Journal
Critical Shear Rate of Polymer-Enhanced Hydraulic Fluids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Methodology for Determining the Energy Losses in Auxiliary Systems and Friction Processes Applied to Low Displacement Diesel Engines

by
Sofia Orjuela Abril
1,
Jhan Piero Rojas
2 and
Eder Norberto Flórez
3,*
1
Programa de Administración de Empresas, Universidad Francisco de Paula Santander, Avenida Gran Colombia No. 12E-96 Barrio Colsag, San José de Cúcuta 540001, Colombia
2
Programa de Ingeniería Civil, Universidad Francisco de Paula Santander, Avenida Gran Colombia No. 12E-96 Barrio Colsag, San José de Cúcuta 540001, Colombia
3
Facultad de Ingeniería, Universidad Francisco de Paula Santander, Sede el Algodonal Vía Acolsure, Ocaña-Norte de Santander 546551, Colombia
*
Author to whom correspondence should be addressed.
Lubricants 2020, 8(12), 103; https://doi.org/10.3390/lubricants8120103
Submission received: 29 October 2020 / Revised: 2 December 2020 / Accepted: 2 December 2020 / Published: 3 December 2020

Abstract

:
The problem of climate change and the reduction of fossil fuels has motivated the development of research focused on improving the efficiency of internal combustion engines. This research proposes a methodology based on mathematical models to determine the energy losses caused by auxiliary systems and friction processes in the engine. Therefore, models are proposed for calculating the energy losses in fuel injection, lubrication, and cooling system. In the same way, models are proposed for the energy losses due to friction in the piston, valve train, and bearings. Experimental tests are carried out on a single-cylinder diesel engine under different operating conditions to validate the proposed models. The results showed that the energy losses of the fuel injection, lubrication, and coolant system are equal to 0.61%, 0.30%, and 0.31% of the chemical energy of the injected fuel. In the case of the energy losses by friction processes, the piston, valve train, and bearings represent 5.47%, 1.34%, and 1.85% of the fuel energy, respectively. Additionally, the proposed model allows estimating the minimum lubrication film present in the piston, valve train, and bearings, which in the particular case of the present study were 0.63 µm, 0.10 µm, and 0.57 µm, respectively. In general, the methodology developed in the present work stands as a robust tool to evaluate the modifications and/or designs of auxiliary systems and friction processes to reduce the energy losses and protect the system from wear caused by lubrication problems. Additionally, the methodology allows evaluating the effect of different types of fuels on the lubrication conditions of the piston and the crankshaft bearings.

Graphical Abstract

1. Introduction

Currently, internal combustion engines (ICE) are used in various economic activities, such as industrial processes, agriculture, electricity generation, and transportation [1,2,3,4]. Power generation through ICEs is of utmost importance for interconnected areas and to maintain a reliable operation that accounts for energy fluctuations and shortage on a grid line. The use of generators powered by ICE increases progressively every year [5]. Despite the advantages of ICEs, the activities that involve their use are responsible for between 15% and 18% of global CO2 emissions [6,7]. This percentage tends to increase due to the globalized growth of economic activities [8]. In order to minimize the environmental impact, international organizations and governments have agreed to set strict levels of emission control for ICE [9]. In order to reduce ICE emissions, the implementation of different post-treatment systems has been investigated, among which it is possible to highlight particle filters, selective catalytic reduction, and the catalytic converter [10,11,12]. However, the structure of the engines used for electricity generation limits the use of systems such as catalytic converters and exhaust gas recirculation [5]. Additionally, the installation of these systems causes high back pressures in the engine exhaust system. Therefore, greater pumping work is required, which produces an increase in fuel consumption [13,14]. As a result, it is necessary to implement alternative strategies to reduce the emissions of ICEs for electric power generation.
The reduction of pollutants in ICEs can be achieved by reducing fuel consumption, which is directly related to the energy efficiency of the engine. Efficiency in ICEs is reduced by mechanical losses, which significantly reduces the useful power of the engine. Mechanical losses are the consequence of energy losses in the pumping system, friction processes, and the energy necessary for the movement of engine components [15]. Therefore, one way to improve ICE efficiency is by reducing energy losses associated with the friction processes.
Friction losses account for approximately 12% of fuel energy. The sources of friction losses are located in the interactions between the cylinder liner and the piston skirt, the camshaft, and the bearings. Based on the prominent contribution of friction losses, different investigations have been carried out to minimize energy losses through different strategies such as oil system heating [16], low viscosity oil [17,18], smooth contact surfaces [19], and coating implementation [20,21].
Another alternative proposed for reducing friction losses is the development of mathematical models that allow quantifying and determining the origin of the energy loss associated with friction [22]. This type of model can be useful to improve the design of critical components such as the piston, the compression rings, the camshaft, and the bearings [23]. Dowson et al. [24] developed a one-dimensional model to estimate friction-associated losses in engines. The model used considers the movement of the elements to calculate the friction produced in the bearings. Dubois and Ocvirk [25] built a model to quantify friction losses in the piston mechanism. Taraza and Henein [26] implemented a mathematical model based on the Stribeck diagram.
Different tribological investigations indicate the importance of mechanical losses caused by friction processes and auxiliary systems. Furuhama and Sasaki [27] evaluated the lubrication performance in the contact between the piston rings and the cylinder liner. Tian [28] pointed out that tribological behavior plays an important role in the wear of the components and the variation of the lubrication film.
Over the years, the complexity of the models has increased since they are used to quantify the losses in specific areas of the engine. Liu et al. [29] studied the influence of surface texture on the lubrication performance between compression rings and cylinder liners. The study was carried out by applying a mixed lubrication model to evaluate the tribological behavior with the presence of spherical dimples in the cylinder liner. Similarly, Avan et al. [30] developed a tribological model to study the influence of the design of the piston rings on the energy losses produced by the contact in the cylinder liner and the piston. The model was used to estimate the thickness of the lubrication film and friction under different lubrication conditions. Jocsak et al. [31] evaluated the influence of different surface parameters on the cylinder liner. The analysis was performed through the implementation of non-Gaussian probability density functions, which were unified with the Greenwood and Tripp roughness contact model and the standard piston-ring friction model. Allmaier et al. [32] studied the thermodynamic behavior of the bearings using an isothermal elastohydrodynamic simulation model. The model was used to calculate the frictional power losses in the sleeve bearings. Likewise, Sander et al. [33] analyzed the mixed lubrication in plain bearings using experimental and numerical approaches. The study focused on the influence of the properties of non-Newtonian lubricants and the piezoviscous effect on power losses due to friction processes in bearings for different dynamic load conditions.
The ICE valve train is another mechanism that has been investigated to reduce friction losses. The strategies used involve the use of different surface textures between the cam-tappet and cam-roller connections [34]. In general, the estimation of the friction forces for this type of mechanism is carried out using the Hertz contact model. Currently, there are different models focused on calculating the thickness of the lubrication film between the follower and the cam lobe [35].
Apart from the individual analysis for the dynamic and lubrication characteristics in the engine’s systems and components, different alternatives have been evaluated in relation to the type of fuel used. Specifically, the study centers on the exploration of biodiesel blends from edible and non-edible vegetable oil, microalgae, water emulsions, and the incorporation of additives such as hydrogen, metallic compounds, and hydroxy. In general, the main focus is to obtain a fuel mixture that guarantees optimal performance while reducing emissions. However, the effect of this type of alternative fuels on the lubrication characteristics is normally not taken into consideration.
Generally, the analysis of the mechanical losses in the engine is performed in specific sections of the engine [36,37], such as the piston assembly [38], bearings [18], and valve train [39]. This makes it difficult to apply in a specific engine study at different operating conditions. Additionally, the identification of potential damaging friction conditions, the quantification of global mechanical losses, and the risk of direct contact between surfaces are more complex when different engine subsystems are considered together.
The objective of this study is to determine the energy losses in the auxiliary systems and the engine friction processes through the development of a compound methodology that integrates the effects of fuel injection, lubrication, and cooling systems. Also, the proposed methodology accounts for the energy dissipated by friction in the piston, valve train, and bearings. The series of mathematical models allow determining the thickness of the lubrication film located on the piston, valve train, and crankshaft bearings. This work serves as a further step to close the knowledge gap regarding energy losses in ICEs while providing a robust tool to determine the impact of modifications in engine components and systems or even the type of fuel in the overall performance of the engine.

2. Energy Loss Models for Auxiliary Systems

The energy losses of an engine can be classified into two groups. The first group corresponds to losses due to friction processes, which occur due to the interactions present in the valve train, piston, and bearings. The second group is a consequence of the auxiliary systems of the engine due to the energy required for its operation. The latter corresponds to the fuel injection, lubrication, and cooling systems. The models used to calculate the energy consumption of the auxiliary systems are described below.

2.1. Fuel Injection System

The power of the pump used in fuel injection systems depends on the pressure drop in the common rail and the volumetric flow. Therefore, the energy consumed from the pump ( E f u e l ) is calculated by means of Equation (1).
E f u e l = P c r · V ˙ f η f
where V ˙ f is the volumetric flow of fuel, P c r the common rail pressure (i.e., is considered as the pressure drop of the fuel injection system) and η f the efficiency of the fuel pump. The volumetric flow of fuel is defined by Equation (2).
V ˙ f = m ˙ f ρ f
where m ˙ f and ρ f are the mass flow and density of the fuel, respectively. Replacing Equation (2) in (1) determines the equation used to calculate the energy consumed by the fuel injection system, which is shown in Equation (3).
E f u e l = P c r · m ˙ f ρ f · η f .

2.2. Lubrication System

The energy consumption of the oil pump is determined by Equation (4).
E ˙ l u b = Δ P l u b · V ˙ l u b η l u b
where P l u b is the pressure drop across the lubrication pipe, V ˙ l u b is the volumetric flow of the lubricant and η l u b is the efficiency of the pump. The pressure drop in the lubrication pipe is determined by the Darcy–Weisbach equation, as shown in Equation (5) [40].
Δ P l u b = 8 · f p i p e g · π 2 · L p i p e D p i p e 2 · V ˙ l u b 2
where D p i p e , L p i p e , g ,   and f p i p e are the diameter of the pipe, the length of the pipe, gravity, and the friction factor of the pipe, respectively. f p i p e is calculated from the Moody formula [41], as shown in Equation (6):
f p i p e = 1.375 × 10 3 · 1 + 800 · V ˙ l u b · ρ l u b · r p i p e + π · D p i p e · μ l u b × 10 6 4 · V ˙ l u b · ρ l u b 1 / 3
where μ l u b , ρ l u b , and r p i p e are the dynamic viscosity of the lubricant, the density of the lubricant, and the roughness of the pipe, respectively.

2.3. Cooling System

In the particular case of this study, the engine cooling system is carried out by means of a forced ventilation system caused by the presence of a fan. The calculation of the energy losses in this type of system can be defined as shown in Equation (7) [42].
E ˙ c = τ · N = k · N 2 · N
where N is the rotational speed of the engine, and k is the fan factor, which depends on the characteristic curve and the design of the blades.
A more direct alternative to calculating the energy losses in these systems is by measuring the electrical operating characteristics of the fan, as shown in Equation (8).
E ˙ c = V s · I
where V s is the supply voltage and I is the electric current. In the case of having cooling systems using centrifugal pumps, it is possible to estimate the energy losses by applying the methodology shown in Section 2.2.

3. Energy Loss Models for Friction Processes

The models used to calculate the energy losses caused by friction processes, which occur due to the interactions present in the valve train, piston, and bearings, are described below.

3.1. Energy Losses in the Valve Train

The interaction between the follower and the cam produces contact friction F c f , v , which is a consequence of the viscosity force due to the lubricant shear and the asperity [43,44]. Therefore, contact friction is defined by Equation (9).
F c f , v = F a c , v + F v f , v
where F a c , v and F v f , v are the asperity contact force and viscous friction force, respectively. The asperity contact force ( F a c , v ) is determined by the Equation (10):
F a c , v = c p l c + τ s s A a
where c p ,   l c ,   A a , and τ s s are the coefficient of the shear strength, load carried by the asperities, asperity area, and Eyring shear stress, respectively.
The load carried by the asperities l c and asperity area ( A a ) are determined by Equations (11) and (12) [45].
l c = 2 · r s r a · 16 · π 15 ρ a · r a · r s · E e · H c · F 5 / 2 2
A a = H c · F 2 · π · ρ a · r a · r s 2
where r s ,   r a ,   ρ a ,   E e , and H c are the composite surface roughness parameter, asperities radius of curvature, asperity density, effective elastic modulus, and Hertzian contact area, respectively.
The Hertzian contact area H c is estimated by considering the contact between the follower and the cam as two cylinders. Therefore, H c is calculated, as shown in Equation (13) [46]:
H c = 24 · r e · F n , v · w c a m 2 π · E e 1 2
where r e is the equivalent radius of curvature, w c a m is the cam width and F n , v is the normal force. The effective elastic modulus E e is calculated by Equation (14).
1 E e = 1 ϑ f o l l 2 2 · E f o l l + 1 ϑ c a m 2 2 · E c a m
where ϑ f o l l and ϑ c a m are the Poisson’s ratios of the follower and cam. E f o l l and E c a m are young’s modulus of the follower and cam.
Variables F 5 / 2 and F 2 are statistical functions dependent on the separation parameter s p , as shown in Equations (15) and (16) [43].
F 5 / 2 = 0.93 + 0.61 s p 0.143 s p 2 + 1.117 × 10 2 s p 3 + 1.47 e s p
F 2 = 1.64 + 1.172 s p 0.31 s p 2 + 0.3 × 10 1 s p 3 + 2.26 e s p
The viscous friction force is determined from the model proposed by Goksem and Hargreaves [47], as shown in Equation (17).
F v f , v = 4.318 · r c c · w c a m c p v · c p v · μ · v e r c c 0.658 · F n , v E e · r c c · w c a m 0.0126
where c p v , r c c , v e and F n , v are the pressure viscosity coefficient, combined radius of curvature, entrainment velocity, and normal force, respectively.
The minimum thickness of the lubrication film h m in the valve train is determined by means of Equation (18) [48].
h m = 2.65 · r e · μ · v e E e · r e 0.7 c p v · E e 0.54 F n , v E e · r e 0.13
Considering the above equations, the energy losses by contact friction in the valve train can be calculated by Equation (19) [49].
E ˙ F c f , v = n i n d · 0 2 π F a c , v θ · v c θ   d θ + n e x h · 0 2 π F v f , v θ · v c θ   d θ
where v c is the contact velocity, n i n d is the number of intake valves and n e x h is the number of exhaust valves.

3.2. Energy Losses in the Piston

To determine frictional energy losses in the piston, a configuration with three piston rings is considered, which consists of the top compression ring, the middle ring, and the lubrication control. This consideration is used since it is the most common configuration for internal combustion engines.
The friction forces in the piston rings depend on the interaction between the lining and the ring surface. Therefore, friction is closely related to the properties, thickness, and height of the lubrication film. Thus, the calculation of the friction force is made by considering the mechanics of the lubrication fluid. When constructing the model, a series of considerations are made:
  • The inner diameter of the ring is considered rigid, so radial inertial forces are ignored.
  • Flow processes are considered isothermal.
  • The radial and tangential friction forces are calculated based on the Stribeck function.
  • The Reynolds equation is used to determine the damping force.
  • The flow of gas pressures is considered stationary.
To determine the relationship between the width, height, and thickness of the lubrication film on the surface of the piston rings, the Reynolds equation is used, which is described by Equation (20) [50].
y h 3 · p y + x h 3 · p x = 6 · η a · h x + ν · h y + 12 · μ · d h d t .
In the present study, the lubrication conditions between the cylinder liner and ring surface are determined as a function of the piston position. Therefore, the Reynolds equation can be reduced, as shown in Equation (21) [51]:
x h x 3 · p x = 6 · η υ · h x x + 12 · μ · d h x d t
where v is the piston velocity, μ is the dynamic viscosity of the lubricant, p lubrication film pressure, and h is the thickness of the lubrication film, which is a function of the position on the piston.
The diagram of the forces applied to the piston rings is described in Figure 1.
F f r , p represents the frictional force produced between the ring face and the cylinder liner and is defined by Equation (22).
F f r , p , i = F f r , s , i + F f r , v , i
where the subscript i can take the values of 1,2, and 3, which represent the top compression ring, the middle ring, and the lubrication control ring, respectively. F f r , s is the asperity friction force caused by the contact between the ring and the cylinder liner, and F f r , v the viscous friction force caused by the viscosity of the lubricant.
The viscous friction force F f r , v is calculated by considering the lubricant as a Newtonian fluid, as shown in Equation (23) [52].
F f r , v , i = a i b i υ · μ h + h 2 · p x d x .
The dynamic viscosity of the lubricant μ is calculated by means of Equation (24) [53].
μ = μ o · e l n μ o μ · 1 + p p a t m c p z T 138 T o 138 S o 1  
where Z is the piezoelectric viscosity, S o is the thermo-viscosity indices and μ o is the dynamic viscosity of the lubricant under atmospheric conditions. The parameters Z and S o are determined from Equations (25) and (26) [50].
Z = c p · α o l n μ 0 μ 1
S o = T o 138 · β o · l n μ 0 μ 1
where β o and α o are thermo-viscosity coefficients and atmospheric piezo-viscosity. μ and c p are constants defined via Equations (27) and (28).
μ = 6.31 × 10 5   Pa · s
c p = 1.98 × 10 8   Pa .
The asperity friction force F f r , s is determined by Equation (29).
F f r , s , i = a i b i c f r · P c   d x
where c f r is the metal-metal friction coefficient and P c is the asperity contact pressure, which is determined by means of the Greenwood and Tripp model [54], as shown in Equation (30):
P c = F 5 / 2 λ · 16 · 2 · π 15 · ρ s 2 · σ c 2 + σ p 2 5 / 4 r c 3 / 2 · E
where r c is the radius of curvature of asperities, ρ s is the asperity density of the surface, E is the young’s modulus between cylinder liner and ring. σ c and σ p are the surface roughness of the cylinder liner and piston ring.
The effective elastic modulus between the cylinder liner and ring E e , c , is calculated by the Equation (31).
1 ϑ 1 2 E 1 + 1 ϑ 2 2 E 2 = 2 E e , c
where ϑ is the Poisson ratio, and E is the elastic modulus. Subscripts 1 and 2 represent the surface area of the cylinder liner and piston ring, respectively.
The term F 5 / 2 λ represents the Gaussian roughness distribution, which is calculated by Equation (32) [55].
F 5 / 2 λ = 0.0046 λ 5 + 0.057 λ 4 0.296 λ 3 + 0.785 λ 2 1.078 λ + 0.62 ,   λ = h σ p 2 + σ r 2   2.23   0   ,   λ = h σ p 2 + σ r 2   > 2.23     .
Finally, the energy losses to friction in the piston are calculated by Equation (33).
E f , p = i = 1 3 0 2 π F f r , v , i · v   d θ + 0 2 π F f r , s , i · v   d θ .

3.3. Energy Losses in Bearings

To describe the energy losses in the bearings, the mobility model is taken as a basis, which considers the trajectory inside the bearing, the location of the center of the journal, and the minimum thickness of the lubrication film [56]. Figure 2 describes the bearing geometry used for the model.
The friction force applied to the bearings is determined from Equation (34) [56].
F f r , b = 2 · v b · F l o a d , b · sin γ + ω · e r · c j b · F l o a d , b · sin ϕ ω · D b + μ · ω · D b 2 · L b · J 1 00 4 · c j b
where μ is the dynamic oil viscosity, ω is the angular speed (this speed is considered to be the same for the bearing and journal), D b is the bearing diameter, L b is the bearing length, c j b is the clearance between journal and bearing, v b is the speed of the bearing centre displacement, ϕ is the angle between F l o a d , b and the centre line and γ is the angle between F l o a d , b and v b . The eccentricity ratio e r is calculated by Equation (35).
e r = e b c j b
where e b is the eccentricity between the bearing centre and journal. F l o a d , b is the instantaneous load applied to the bearing, and its value depends on the location (i.e., crankshaft or connecting rod). The applied load on the connecting rod bearing (bearing 1) is determined directly from the dynamic analysis of the piston mechanism. In the case of crankshaft bearings (bearing 2 and bearing 3), the load they support is considered to be half the force exerted on bearing 1. This consideration is in accordance with the results recorded in the literature [57]. Figure 3 shows the forces applied to the previously described bearings.
The thickness of the lubrication film on the bearings is determined by Equation (36).
h = c r c · 1 c j b · cos ϕ .
where c r c is the radial clearance in the bearings.
The parameter J 1 00 present in Equation (34) describes the behavior of the extension and the change in thickness of the lubrication film around the bearing. To calculate this parameter, the equation proposed by Taylor is used [56], which is shown in Equation (37).
J 1 00 = 0 2 π d θ e r · cos θ + 1 = 2 π 1 e r 2 .
The model used in the present work is considered a quasi-stable condition, which implies that there is no translational movement in the bearing. Therefore, the speed of the displacement of the bearing center is zero v b = 0 [26]. Hence, Equation (34) can be expressed, as shown in Equation (38) [56]:
F f r , b = e r · c j b · F l o a d , b · sin ϕ D b + μ · ω · D b 2 · π · L b 2 · c j b 1 e r 2 .
The angle between F l o a d , b and the centre line ϕ , is determined from Equation (39).
ϕ = tan 1 π 4   · 1 e r 2 e r 2 .
Solving Equations (35), (38) and (39) together determine the energy losses by the friction process in the bearings, which is calculated by Equation (40).
E f , b = i = 1 n = i 0 2 π D b , i ω · F f r , b , i 2 d θ
where n is the total number of bearings to consider.

4. Numerical Methodology

The methodology used in this research is divided into two parts, which consist of an experimental part for the compilation of test data that allows determining the energy losses mentioned in Section 2 and a numerical part for the solution of the equations shown in Section 3. The experimental test bench and the measuring instruments used for the study are described in Section 5.
The equations described in Section 3 are solved using the solver tool ode45 of the Matlab® software. Figure 4 describes the diagram for the development of the experimental and numerical methodology.
Based on the characteristics of the engine and the operating conditions, the losses energy calculations caused by the auxiliary systems (Section 2) and the friction processes (Section 3) are made. Verification of the model results is carried out through experimental validation, which is discussed in Section 6.1. Finally, the results of the model are used to evaluate the contribution of the energy losses of each system and process, the influence of the engine operating conditions, and the analysis of the lubrication film on the components.

5. Experimental Methodology

The experimental tests were carried out on a test bench made up of a single-cylinder diesel engine. The technical characteristics of the engine are shown in Table 1.
The properties of the lubricant used are shown in Table 2. In the case of fuel, commercial diesel was used with an addition of 0.075 LPM of hydrogen, which entered through the air intake system of the engine. The properties of diesel and hydrogen are described in Table 3.
For the analysis of the energy losses of the engine, nine operating modes were established, as shown in Figure 5. These modes consist of variations in torque and rotation speed between 4–6 Nm and 3400–3800 rpm, respectively. In this way, it is possible to evaluate the energy losses that occur in the engine for low, medium, and high load conditions, thus covering the different conditions in which the engine can operate.
The control of the torque and rotation speed is carried out by means of a dynamometer, which is coupled to the engine shaft. To calculate the energy losses in the fuel injection system, a gravimetric meter is used to determine the fuel flow (OHAUS PA313, OHAUS Corporation, Parsippany, NJ, USA). Additionally, the common rail pressure measurement is performed by means of a pressure sensor (KISTLER Type 4067-E, KISTLER, Winterthur, Switzerland). The measurement of the lubrication flow is carried out by using a turbine-type flow sensor, which is used to calculate the losses in the lubrication system. Finally, the electrical characteristics of voltage and current of the fan in charge of cooling the engine are recorded. For this, a digital multimeter (Fluke 15b-F15b, Fluke Corporation, Everett, WA, USA) is used. The engine test bench diagram is shown in Figure 6.
To calculate the energy losses for friction processes, it is necessary to determine the piston’s position and the pressure in the combustion chamber since they are input data for the solution of the equations proposed. The combustion pressure measurement is carried out utilizing a piezoelectric sensor (KISTLER type 7063-A, KISTLER, Winterthur, Switzerland) installed inside the chamber. Synchronization between chamber pressure and the crankshaft rotation angle is done using an angle encoder (BECK ARNLEY 180-0420, Beck / Arnley Company, Smyrna, TN, USA). Table 4 shows the characteristics of the measurement instruments used.

6. Results and Discussions

6.1. Experimental Validation

To verify the prediction capacity of the models proposed in Section 2 and Section 3, the total energy losses caused by engine friction in the different operating modes are experimentally determined. The evaluation of friction losses is determined by calculating the indicated mean effective pressure (IMEP) and brake mean effective pressure (BMEP), as shown in Equation (41).
F M E P = I M E P B M E P
where FMEP (Friction Mean Effective Pressure) represents the power required to overcome engine friction. Equation (41) can be rewritten as:
F M E P = p d V v d 4 · π · T v d
where T is the measured torque of the engine and v d is the displaced volume in the chamber.
In the case of auxiliary systems, the characteristic curve of the fuel and lubricant pumps supplied by the manufacturer is used to determine the power consumption associated with each quantity of flow. For the cooling system, the power consumed is calculated directly by the current and voltage conditions of the engine fan. Figure 7 shows the comparison between the results of the model and the experimental ones.
Figure 7 shows that the proposed model maintains the same trend when compared with the experimental results. An average error of 6% was obtained between the model estimates and the experimental data.
In general, low and medium load conditions have the lowest error rates. As higher rotational speeds and torques are operated, the error in the model estimates tends to increase. The foregoing can be the consequence of greater operational instability of the engine under high load conditions, which can be attributed to variations in chamber pressure throughout the combustion cycle, greater vibrations, and changes in the characteristics of the lubricant. However, the maximum error recorded is less than 10%.
In the following sections (Section 6.2, Section 6.3, Section 6.4 and Section 6.5), the analysis of the calculations obtained from the auxiliary systems and the friction processes in the engine is performed.

6.2. Analysis of Energy Losses in the Engine

Figure 8 shows the energy losses of the auxiliary engine systems for the different operating conditions. Notice that the bar graph shown in Figure 8 displays the contribution of each system represented with different patterns as a portion of the total energy losses.
According to the results, the greatest energy losses occur in the fuel injection system. On average, an energy loss of 22.9, 21.8, and 44.1 W was observed for the lubrication, cooling, and fuel injection systems, respectively. The maximum values of energy losses in each of the above systems are 26 W, 27.4 W, and 54 W, which is equivalent to 0.40, 0.39, and 0.73% of the chemical energy of the fuel injected into the engine.
For the fuel injection system, it was observed that the energy loss tends to increase with the increase in the rotational speed and mechanical power of the engine. The foregoing is a consequence of the greater need for fuel pumping to reach the highest levels of load and engine speed. In the range of test conditions, the energy loss was of the fuel system varies between 35–54 W.
In the case of the cooling system, it is shown that the pumping power is increased mainly by the increase in torque in the engine. In general, increases in rotation speed do not cause a considerable change in the consumption of the cooling pump. However, by increasing the torque of the engine by 2 Nm, the energy consumed increases by 24%. The energy consumed by the lubrication system is considerably increased for rotational speeds above 3600 rpm. This is a consequence of the greater need for lubrication flow for components subjected to high accelerations. For the test conditions, the consumption of the lubrication pump was observed between 17–26 W, which is similar to the energy loss developed by the cooling pump.
Figure 9 shows the energy losses produced by the friction processes in the bearings, valve train, and piston.
The results indicate that the main friction energy loss occurs in the piston. In general, the energy loss to friction in the piston is three to four times greater than in the bearings and valve train. On average, a loss of 120.4, 86.8, and 354.72 W was observed for the bearings, valve train, and piston, which is equivalent to 1.85, 1.34, and 5.47% of the total energy available from fuel injection.
In addition, it can be observed that the energy loss in the piston increases with the torque and the rotational speed of the engine. This can be a consequence of the friction coefficient, which is directly related to the speed of the bodies. For the energy loss in the valve train and the bearings, it is noticed that the greatest losses occur at the highest torque. In contrast, the increase in the rotation speed does not produce a significant change in the increase in friction losses.

6.3. Analysis of Energy Loss Distributions

Figure 10 shows the average distribution of the total energy losses of the engine due to friction processes and auxiliary systems.
The results show that friction processes are the main contributors to engine energy losses. Among this type, the piston occupies more than half of the total losses (55%), followed by the bearings (19%) and the valve train (13%). In the case of auxiliary systems, it was found that the lubrication and cooling systems present a similar share of energy loss, which corresponds to between 3% and 4%, respectively. Of the auxiliary systems, fuel injection is the one with the highest percentage of total losses, reaching approximately 6.6%.

6.4. Lubrication Film Analysis

6.4.1. Piston

In Figure 11, the minimum thickness of the lubrication film between the piston skirt and the cylinder liner is shown.
The results show that the thickness of the lubrication film decreases with increasing rotational speed and engine torque. Increasing these parameters causes an increase in the temperature of the lubricant, which reduces the viscosity of the oil. The foregoing allows explaining the reduction in the lubrication film. The thickness of the lubrication film is recorded during the compression and combustion stage, which can be an effect of the increase in pressure in the piston rings.
The minimum lubrication thickness recorded for the operating conditions corresponds to 0.63 μ m . Specifically, an increase of 400 rpm in the rotational speed of the engine causes a reduction of the thickness by 11%, while an increase of 4 Nm in torque reduces approximately 25%.

6.4.2. Valve Train

The valve train is composed of different elements such as bearings, valve stem, and camshaft. These elements are subjected to friction, which leads to energy losses. However, the contact between the follower and the cam represents the largest contribution to the energy losses of the valve train system. Figure 12 shows the energy loss produced by the follower-cam contact throughout the combustion cycle.
The results in Figure 12 indicate that the intake stage is where the energy losses of the valve train are concentrated. This can be explained by considering that in this location, the lobe of the cam exerts high pressure on the roller, which is necessary to force the movement of the valve. Therefore, the friction force reaches its maximum value in this condition. This conclusion can be validated by observing the thickness of the lubrication film between the follower and the cam, as shown in Figure 13. Here, the thickness remains relatively constant with a value of 2 μ m during the compression stage, combustion, and escape. However, during the intake stage, the thickness of the film decreases until it reaches a minimum value of 0.10 μ m .

6.4.3. Bearings

Figure 14 shows the minimum lubrication film thickness for bearings 1, 2, and 3 (see Figure 3).
The big end bearing of the connecting rod (bearing 1) is in charge of supporting the highest load conditions because the pressures resulting from the combustion cylinder are directly applied in this location. Bearings 2 and 3 handle a lesser load. This behavior can be seen when comparing the minimum thicknesses of the lubrication film achieved by the three bearings. For an operating condition of 3800 rpm–8 Nm in the engine, it was observed that the minimum thickness in the lubrication film is 0.57, 1.56, and 1.13 μ m for bearing 1, 2, and 3, respectively. The lower thickness located in bearing 1 is an indication that it is subjected to a greater load.
The lower thicknesses of the lubrication film for all the bearings are obtained during the compression and combustion stage, which is an effect of the high pressures in the combustion chamber developed in these stages. The results indicate that the thickness of bearings 1, 2, and 3 is reduced by 10.57, 5.43, and 8.57 μ m , respectively, throughout the combustion cycle.

6.5. Friction Force Analysis

Figure 15 shows the change in the friction force on the piston during the combustion cycle.
The results show that the maximum levels of friction force occur at the top dead center and during the combustion stage. This can be the consequence of a mixed lubrication condition, in which there is greater interaction between the surface thickness of the rings and the cylinder liner. For the operational condition with maximum torque (8 Nm) and rotational speed (3800 rpm), the highest value of friction force (23.82 N) was obtained.
In general, the friction force is lower during the intake and exhaust stages. This can be attributed to the reduction of friction interactions due to the low combustion pressures during these stages. This behavior is in accordance with the predictions reported in the literature [58].
Figure 16 shows the friction force in the valve train for the different operating conditions.
The results show that the friction force in the valve train varies between 1.91 N and 3.62 N for the tested operating conditions. On average, increasing the rotational speed by 200 rpm and the torque by 2 Nm enlarge the friction force by 9% and 26%, respectively. This is associated with the beginning of a mixed lubrication regime between the follower and cam.
Figure 17 shows the friction force acting on the main bearing of the crankshaft (bearing 1).
The results of Figure 17 show a behavior similar to that described in Figure 16. It was observed that the friction force has a strong influence on the rotation speed and the torque of the engine. For the operating conditions tested, the friction force is maintained in a range of 1.07 N to 4.81 N.

7. Conclusions

In this research, a numerical methodology is developed for estimating the combined effects of energy losses in the auxiliary systems (fuel injection system, lubrication system, and cooling system) and the friction processes of the engine components (piston, valve train, and bearings). The methodology incorporated mathematical models in which the dynamic characteristics and tribological conditions present in the piston, valve train, and bearings are considered. Also, the model is capable of determining the consumption of fuel, lubrication, and cooling pumps, as well as the lubrication film conditions on the piston, valve train, and bearings.
The analysis of the auxiliary systems demonstrated that the fuel injection system features the greatest energy losses. Indeed, these losses double those obtained in the lubrication and cooling systems. In the case of the latter two, it was observed that both systems present an average consumption of 22.35 W for the tested conditions (4–8 Nm and 3400–3800 rpm). The energy losses of the fuel injection system, lubrication, and coolant are equal to 0.61, 0.30, and 0.31% of the total chemical energy of the injected fuel, respectively.
For energy losses due to friction processes, it was observed that the main losses are located in the piston, which reaches a maximum value of 391 W. In the case of the valve train and the bearings, energy losses were found between 59–140 W for the operating range. On average, the energy losses by friction processes in the piston, valve train, and bearings represent 5.47, 1.34, and 1.85% of the total fuel energy, respectively.
In general, friction processes cover 87% of total mechanical losses, of which 55% is produced by the interaction of the pistons. The fuel injection system represents 6.5% of total mechanical losses.
The analysis of the lubrication film in the piston showed a decreasing trend as the torque and engine rotation speed increase, reaching a minimum thickness of 0.63 μ m in the highest load condition. For the valve train, it is shown that its energy loss occurs mainly during the inlet and exhaust stage, in which the minimum lubrication film thickness is 0.10 μ m . By analyzing the three crankshaft bearings, it was shown that the connecting rod bearing is subjected to a higher load, which causes a reduction in the thickness of the lubrication film.
In general, the model developed in this work demonstrated to be a useful tool to estimate the energy losses produced in the auxiliary systems and the components subjected to friction. This can be applied to different operating conditions, including the type of fuel and geometric characteristics of a specific engine. Therefore, the numerical methodology implemented helps to analyze the influence of modifications and/or design of engine components on the global behavior of energy losses, which contributes to identifying critical lubrication conditions that minimize the energy wasted in internal combustion engines.

Author Contributions

Conceptualization, S.O.A., and E.N.F.; Methodology, E.N.F., and J.P.R.; Software, E.N.F., and J.P.R.; Validation, E.N.F.; Formal analysis, E.N.F.; Investigation, E.N.F.; Resources, E.N.F., Writing—Original Draft Preparation, E.N.F.; Writing—Review and Editing, S.O.A., and J.P.R.; Funding acquisition, E.N.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to acknowledge the Universidad Francisco de Paula Santander and Sphere Energy company for their support in the development of this investigation.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
E Energy losses
V ˙ Volumetric flow of fuel
P c r Common rail pressure
m ˙ Mass flow
Δ P Pressure drop
D Diameter
L Length
f Friction factor
g Gravity
r Roughness
N Rotational speed
k Fan factor
V s Supply voltage
I Electric current
F c f , v Contact friction between the follower and the cam
F a c , v Asperity contact force between the follower and the cam
F v f , v Viscous friction force between the follower and the cam
c p Coefficient of the shear strength
l c Load carried by the asperities
A a Asperity area
τ s s Eyring shear stress
r s Composite surface roughness parameter
r a Asperities radius of curvature
ρ a Asperity density
E e Effective elastic modulus
H c Hertzian contact area
r e Equivalent radius of curvature
w c a m Cam width
F n , v Force normal
ϑ f o l l Poisson ratio of the follower
ϑ c a m Poisson ratio of the cam
E f o l l Young module of the follower and cam
E c a m Young module of the cam
F 2 Statistical function
s p Separation parameter
h Thickness of the lubrication film
c p v Pressure viscosity coefficient
r c c Combined radius of curvature
v e Entrainment velocity
v c Contact velocity
n i n d Number of intake valves
n e x h Number of exhaust valves
v Piston velocity
p Lubrication film pressure
F f r , p Frictional force produced between the ring face and the cylinder liner
F f r , s Asperity friction force between the ring face and the cylinder liner
F f r , v Viscous friction force between the ring face and the cylinder liner
Z Piezoelectric viscosity
S o Thermo-viscosity indice
β o Thermo-viscosity coefficients
α o Atmospheric piezo-viscosity
c f r Metal-metal friction coefficient
P c Asperity contact pressure
r c Radius of curvature of asperities
ρ s Asperity density of the surface
σ c Surface roughness of the cylinder liner
σ p Surface roughness of the piston ring
F 5 / 2 Gaussian roughness distribution
ω Angular speed
D b Bearing length
L b Bearing length
c j b Clearance between journal and bearing
v b Speed of the bearing centre displacement
e r Eccentricity ratio
e b Eccentricity between the bearing centre and journal
F l o a d , b Instantaneous load applied to the bearing
c r c Radial clearance
Greek Letters
η Efficiency
μ Dynamic viscosity
ρ Density
Subscripts
f Fuel injection system
l u b Oil pump
o Atmospheric conditions
c Cooling system

References

  1. Zhang, Q.; Chen, G.; Zheng, Z.; Liu, H.; Xu, J.; Yao, M. Combustion and emissions of 2,5-dimethylfuran addition on a diesel engine with low temperature combustion. Fuel 2013, 103, 730–735. [Google Scholar] [CrossRef]
  2. Xiao, H.; Guo, F.; Li, S.; Wang, R.; Yang, X. Combustion performance and emission characteristics of a diesel engine burning biodiesel blended with n-butanol. Fuel 2019, 258, 115887. [Google Scholar] [CrossRef]
  3. Aydin, H.; Ilkılıç, C. Effect of ethanol blending with biodiesel on engine performance and exhaust emissions in a CI engine. Appl. Therm. Eng. 2010, 30, 1199–1204. [Google Scholar] [CrossRef]
  4. Xiao, H.; Guo, F.; Wang, R.; Yang, X.; Li, S.; Xiao, H. Combustion performance and emission characteristics of diesel engine fueled with iso-butanol/biodiesel blends. Fuel 2020, 268, 117387. [Google Scholar] [CrossRef]
  5. Erdoğan, S.; Balki, M.K.; Aydın, S.; Sayın, C. Performance, emission and combustion characteristic assessment of biodiesels derived from beef bone marrow in a diesel generator. Energy 2020, 207, 118300. [Google Scholar] [CrossRef]
  6. Mejía, A.; Leiva, M.; Rincón-Montenegro, A.; Gonzalez-Quiroga, A.; Duarte, J. Experimental assessment of emissions maps of a single-cylinder compression ignition engine powered by diesel and palm oil biodiesel-diesel fuel blends. Case Stud. Therm. Eng. 2020, 19, 100613. [Google Scholar] [CrossRef]
  7. Ochoa, G.V.; Acevedo, C.; Duarte, J. Thermo-Economic Assessment of a Gas Microturbine-Absorption Chiller Trigeneration System under Different Compressor Inlet Air Temperatures. Energies 2019, 12, 4643. [Google Scholar] [CrossRef] [Green Version]
  8. Singh, D.; Sharma, D.; Soni, S.; Inda, C.S.; Sharma, S.; Sharma, P.K.; Jhalani, A. A comprehensive review of physicochemical properties, production process, performance and emissions characteristics of 2nd generation biodiesel feedstock: Jatropha curcas. Fuel 2021, 285, 119110. [Google Scholar] [CrossRef]
  9. European Commission. Setting emisison performace standards for new passenger cars as part of the community’s integrated approach to reduce CO2 emissions from light-duty vehicles. Off. J. Eur. Commun. 2009, 140, 5–6. [Google Scholar]
  10. Duarte, J.; Garcia, J.; Jimenez, J.; Sanjuan, M.E.; Bula, A.; Gonzalez, J. Auto-ignition control in spark-ignition engines using internal model control structure. J. Energy Res. Technol. 2017, 139, 022201. [Google Scholar] [CrossRef]
  11. Pietikäinen, M.; Väliheikki, A.; Oravisjärvi, K.; Kolli, T.; Huuhtanen, M.; Niemi, S.; Virtanen, S.; Karhu, T.; Keiski, R.L. Particle and NOx emissions of a non-road diesel engine with an SCR unit: The effect of fuel. Renew. Energy 2015, 77, 377–385. [Google Scholar] [CrossRef]
  12. Beatrice, C.; Di Iorio, S.; Guido, C.; Napolitano, P. Detailed characterization of particulate emissions of an automotive catalyzed DPF using actual regeneration strategies. Exp. Therm. Fluid Sci. 2012, 39, 45–53. [Google Scholar] [CrossRef]
  13. Diaz, G.A.; Duarte, J.; Rincon, A.; Fontalvo, A.; Bula, A.; Padilla, R.V.; Orozco, W. Characteristics of Auto-Ignition in Internal Combustion Engines Operated With Gaseous Fuels of Variable Methane Number. J. Energy Resour. Technol. 2017, 139, 042205. [Google Scholar] [CrossRef]
  14. Diaz, G.A.; Duarte, J.; Garcia, J.; Rincon, A.; Fontalvo, A.; Bula, A.; Padilla, R.V. Maximum Power From Fluid Flow by Applying the First and Second Laws of Thermodynamics. J. Energy Resour. Technol. 2017, 139, 032903. [Google Scholar] [CrossRef]
  15. Ochoa, G.V.; Isaza-Roldan, C.; Duarte, J. Economic and Exergo-Advance Analysis of a Waste Heat Recovery System Based on Regenerative Organic Rankine Cycle under Organic Fluids with Low Global Warming Potential. Energies 2020, 13, 1317. [Google Scholar] [CrossRef] [Green Version]
  16. Forero, J.D.; Ochoa, G.V.; Alvarado, W.P. Study of the Piston Secondary Movement on the Tribological Performance of a Single Cylinder Low-Displacement Diesel Engine. Lubricants 2020, 8, 97. [Google Scholar] [CrossRef]
  17. Duarte Forero, J.; Valencia Ochoa, G.; Piero Rojas, J. Effect of the Geometric Profile of Top Ring on the Tribological Characteristics of a Low-Displacement Diesel Engine. Lubricants 2020, 8, 83. [Google Scholar] [CrossRef]
  18. Knauder, C.; Allmaier, H.; Sander, D.E.; Salhofer, S.; Reich, F.M.; Sams, T. Analysis of the Journal Bearing Friction Losses in a Heavy-Duty Diesel Engine. Lubricants 2015, 3, 142–154. [Google Scholar] [CrossRef]
  19. Morawitz, U.; Mehring, J.; Schramm, L. Benefits of Thermal Spray Coatings in Internal Combustion Engines, with Specific View on Friction Reduction and Thermal Management. SAE Tech. Pap. Ser. 2013. [Google Scholar] [CrossRef]
  20. Tormos, B.; Ramírez, L.; Johansson, J.; Björling, M.; Larsson, R. Fuel consumption and friction benefits of low viscosity engine oils for heavy duty applications. Tribol. Int. 2017, 110, 23–34. [Google Scholar] [CrossRef]
  21. Tomanik, E. Friction and wear bench tests of different engine liner surface finishes. Tribol. Int. 2008, 41, 1032–1038. [Google Scholar] [CrossRef]
  22. Consuegra, F.; Bula, A.; Guillín, W.; Sánchez, J.; Duarte, J. Instantaneous in-Cylinder Volume Considering Deformation and Clearance due to Lubricating Film in Reciprocating Internal Combustion Engines. Energies 2019, 12, 1437. [Google Scholar] [CrossRef] [Green Version]
  23. Hoshi, M. Reducing friction losses in automobile engines. Tribol. Int. 1984, 17, 185–189. [Google Scholar] [CrossRef]
  24. Dowson, D.; Taylor, C.M.; Yang, L. Friction modelling for internal combustion engines. In Tribology Series; Elsevier: Amsterdam, The Netherlands, 1996; pp. 301–318. [Google Scholar]
  25. DuBois, G.B.; Ocvirk, F.W. Analytical Derivation and Experimental Evaluation of Short-Bearing Approximation for Full Journal Bearings. Natl. Advis. Comm. Aeronaut. 1953, 1157, 1–32. [Google Scholar]
  26. Taraza, D.; A Henein, N.; Bryzik, W. Friction Losses in Multi-Cylinder Diesel Engines; SAE International: Warrendale, PA, USA, 2000. [Google Scholar]
  27. Furuhama, S.; Sasaki, S. Effect of oil properties on piston frictional forces. Int. J. Veh. Des. 1986. [Google Scholar] [CrossRef]
  28. Tian, T. Dynamic behaviours of piston rings and their practical impact. Part 2: Oil transport, friction and wear of ring/liner interface and the effects of piston and ring dynamics. Proc. Inst. Mech. Eng. Part J: J. Eng. Tribol. 2002, 216, 229–248. [Google Scholar] [CrossRef]
  29. Liu, C.; Lu, Y.; Zhang, Y.-F.; Li, S.; Müller, N. Numerical study on the lubrication performance of compression ring-cylinder liner system with spherical dimples. PLoS ONE 2017, 12, e0181574. [Google Scholar] [CrossRef]
  30. Avan, E.Y.; Spencer, A.; Dwyer-Joyce, R.S.; Almqvist, A.; Larsson, R. Experimental and numerical investigations of oil film formation and friction in a piston ring–liner contact. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2012, 227, 126–140. [Google Scholar] [CrossRef]
  31. Jocsak, J.; Wong, V.W.; Tian, T. The Effects of Cylinder Liner Finish on Piston Ring-Pack Friction. In Proceedings of the ASME 2004 Internal Combustion Engine Division Fall Technical Conference, Long Beach, CA, USA, 24–27 October 2004; pp. 841–849. [Google Scholar] [CrossRef]
  32. Allmaier, H.; Sander, D.E.; Reich, F. Simulating Friction Power Losses in Automotive Journal Bearings. Procedia Eng. 2013, 68, 49–55. [Google Scholar] [CrossRef] [Green Version]
  33. Sander, D.; Allmaier, H.; Priebsch, H.; Reich, F.; Witt, M.; Füllenbach, T.; Skiadas, A.; Brouwer, L.; Schwarze, H. Impact of high pressure and shear thinning on journal bearing friction. Tribol. Int. 2015, 81, 29–37. [Google Scholar] [CrossRef]
  34. Turturro, A.; Rahmani, R.K.; Rahnejat, H.; Delprete, C.; Magro, L. Assessment of Friction for Cam-Roller Follower Valve Train System Subjected to Mixed Non-Newtonian Regime of Lubrication. In ASME 2012 Internal Combustion Engine Division Spring Technical Conference; ASME International: New York, NY, USA, 2012; pp. 917–923. [Google Scholar]
  35. Masjedi, M.; Khonsari, M. Film Thickness and Asperity Load Formulas for Line-Contact Elastohydrodynamic Lubrication With Provision for Surface Roughness. J. Tribol. 2012, 134, 011503. [Google Scholar] [CrossRef]
  36. Priest, M.; Taylor, C. Automobile engine tribology—Approaching the surface. Wear 2000, 241, 193–203. [Google Scholar] [CrossRef]
  37. Offner, G. Friction Power Loss Simulation of Internal Combustion Engines Considering Mixed Lubricated Radial Slider, Axial Slider and Piston to Liner Contacts. Tribol. Trans. 2013, 56, 503–515. [Google Scholar] [CrossRef]
  38. Westerfield, Z.; Totaro, P.; Kim, D.; Tian, T. An Experimental Study of Piston Skirt Roughness and Profiles on Piston Friction Using the Floating Liner Engine; SAE International: Warrendale, PA, USA, 2016. [Google Scholar]
  39. Mufti, R.A.; Priest, M. Experimental and Theoretical Study of Instantaneous Engine Valve Train Friction. J. Tribol. 2003, 125, 628–637. [Google Scholar] [CrossRef]
  40. Marušić-Paloka, E.; Pažanin, I. Effects of boundary roughness and inertia on the fluid flow through a corrugated pipe and the formula for the Darcy-Weisbach friction coefficient. Int. J. Eng. Sci. 2020, 152, 103293. [Google Scholar] [CrossRef]
  41. Moody, L.F. Approximate formula for pipe friction factors. Mech. Eng. 1947, 69, 1005–1011. [Google Scholar]
  42. Wang, T.T.; Jagarwal, A.; Wagner, J.R.; Fadel, G. Optimization of an automotive radiator fan array operation to reduce power consumption. IEEE ASME Trans. Mechatronics 2014, 20, 2359–2369. [Google Scholar] [CrossRef]
  43. Teodorescu, M.; Taraza, D.; Henein, N.A.; Bryzik, W. Simplified Elasto-Hydrodynamic Friction Model of the Cam-Tappet Contact. SAE Tech. Pap. Ser. 2003. [Google Scholar] [CrossRef]
  44. Guo, J.; Zhang, W.; Zou, D. Investigation of dynamic characteristics of a valve train system. Mech. Mach. Theory 2011, 46, 1950–1969. [Google Scholar] [CrossRef]
  45. Greenwood, J.A.; Williamson, J.B.P. Contact of nominally flat surfaces. Proc. R. Soc. London. Ser. A Math. Phys. Sci. 1966, 295, 300–319. [Google Scholar] [CrossRef]
  46. Sadeghi, F.; Slack, T. Handbook of Lubrication and Tribology; CRC Press: Boca Raton, FL, USA, 2012; Volume II. [Google Scholar] [CrossRef]
  47. Goksem, P.G.; Hargreaves, R.A. The Effect of Viscous Shear Heating on Both Film Thickness and Rolling Traction in an EHL Line Contact—Part I: Fully Flooded Conditions. J. Lubricants Technol. 1978, 100, 346–352. [Google Scholar] [CrossRef]
  48. Dowson, D.; Higginson, G.R. Elasto-Hydrodynamic Lubrication; Pergamon: Oxford, UK, 1977. [Google Scholar]
  49. Menacer, B.; Bouchetara, M. The compression ring profile influence on hydrodynamic performance of the lubricant in diesel engine. Adv. Mech. Eng. 2020, 12, 168781402093084. [Google Scholar] [CrossRef]
  50. Delprete, C.; Razavykia, A.; Baldissera, P. Detailed analysis of piston secondary motion and tribological performance. Int. J. Engine Res. 2020, 21, 1647–1661. [Google Scholar] [CrossRef]
  51. Jang, S.; Cho, J. Effects of skirt profiles on the piston secondary movements by the lubrication behaviors. Int. J. Automot. Technol. 2004, 5, 23–31. [Google Scholar]
  52. Delprete, C.; Razavykia, A. Piston dynamics, lubrication and tribological performance evaluation: A review. Int. J. Engine Res. 2020, 21, 725–741. [Google Scholar] [CrossRef]
  53. Jeng, Y.-R. Theoretical Analysis of Piston-Ring Lubrication Part I—Fully Flooded Lubrication. Tribol. Trans. 1992, 35, 696–706. [Google Scholar] [CrossRef]
  54. Greenwood, J.A.; Tripp, J.H. The Contact of Two Nominally Flat Rough Surfaces. Proc. Inst. Mech. Eng. 1970, 185, 625–633. [Google Scholar] [CrossRef]
  55. Teodorescu, M.; Balakrishnan, S.; Rahnejat, H. Integrated Tribological Analysis within a Multi-physics Approach to System Dynamics. Tribol. Metal Cutt. 2005, 48, 725–737. [Google Scholar] [CrossRef]
  56. Taylor, C. Engine tribology. Tribol. Int. 1997, 30, 464–465. [Google Scholar] [CrossRef]
  57. Internal combustion engine fundamentals. Choice Rev. Online 1988, 26, 26. [CrossRef]
  58. Nikolakopoulos, P.G.; Mavroudis, S.; Zavos, A. Lubrication Performance of Engine Commercial Oils with Different Performance Levels: The Effect of Engine Synthetic Oil Aging on Piston Ring Tribology under Real Engine Conditions. Lubricants 2018, 6, 90. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Frictional forces applied to piston rings.
Figure 1. Frictional forces applied to piston rings.
Lubricants 08 00103 g001
Figure 2. Bearing internal geometry.
Figure 2. Bearing internal geometry.
Lubricants 08 00103 g002
Figure 3. Load diagram for engine bearings.
Figure 3. Load diagram for engine bearings.
Lubricants 08 00103 g003
Figure 4. Methodology used in this study.
Figure 4. Methodology used in this study.
Lubricants 08 00103 g004
Figure 5. Engine operating modes, A: 4 Nm-3400 rpm, B: 4 Nm-3600 rpm, C: 4 Nm-3800 rpm, D: 6 Nm-3400 rpm, E: 4 Nm-3600 rpm, F: 6 Nm-3800 rpm, G: 8 Nm-3400 rpm, H: 8 Nm-3600 rpm, I: 8 Nm-3800 rpm.
Figure 5. Engine operating modes, A: 4 Nm-3400 rpm, B: 4 Nm-3600 rpm, C: 4 Nm-3800 rpm, D: 6 Nm-3400 rpm, E: 4 Nm-3600 rpm, F: 6 Nm-3800 rpm, G: 8 Nm-3400 rpm, H: 8 Nm-3600 rpm, I: 8 Nm-3800 rpm.
Lubricants 08 00103 g005
Figure 6. Engine test bench. 1. Data acquisition (DAQ), 2. Crank angle encoder, 3. Dynamometer, 4. Resistive bench, 5. Engine, 6. Median variables DAQ, 7. Diesel tank, 8. Fuel inlet valve, 9. Fuel filter, 10. Injection pump, 11. Airflow meter.
Figure 6. Engine test bench. 1. Data acquisition (DAQ), 2. Crank angle encoder, 3. Dynamometer, 4. Resistive bench, 5. Engine, 6. Median variables DAQ, 7. Diesel tank, 8. Fuel inlet valve, 9. Fuel filter, 10. Injection pump, 11. Airflow meter.
Lubricants 08 00103 g006
Figure 7. Energy losses by auxiliary systems and friction processes were obtained experimentally and modeled.
Figure 7. Energy losses by auxiliary systems and friction processes were obtained experimentally and modeled.
Lubricants 08 00103 g007
Figure 8. Energy losses in auxiliary systems.
Figure 8. Energy losses in auxiliary systems.
Lubricants 08 00103 g008
Figure 9. Losses of energy by friction.
Figure 9. Losses of energy by friction.
Lubricants 08 00103 g009
Figure 10. Distribution of energy losses.
Figure 10. Distribution of energy losses.
Lubricants 08 00103 g010
Figure 11. Minimum thickness of the lubrication film on the piston.
Figure 11. Minimum thickness of the lubrication film on the piston.
Lubricants 08 00103 g011
Figure 12. Power losses in valve train (3800 rpm–8 Nm).
Figure 12. Power losses in valve train (3800 rpm–8 Nm).
Lubricants 08 00103 g012
Figure 13. Minimum thickness of the lubrication film between follower-cam (3800 rpm–8 Nm).
Figure 13. Minimum thickness of the lubrication film between follower-cam (3800 rpm–8 Nm).
Lubricants 08 00103 g013
Figure 14. Minimum lubrication film thickness for engine bearings (3800 rpm–8 Nm).
Figure 14. Minimum lubrication film thickness for engine bearings (3800 rpm–8 Nm).
Lubricants 08 00103 g014
Figure 15. Variation of friction force on piston rings.
Figure 15. Variation of friction force on piston rings.
Lubricants 08 00103 g015
Figure 16. Friction forces in the valve train.
Figure 16. Friction forces in the valve train.
Lubricants 08 00103 g016
Figure 17. Friction forces in the main bearing of the crankshaft.
Figure 17. Friction forces in the main bearing of the crankshaft.
Lubricants 08 00103 g017
Table 1. Single cylinder diesel engine.
Table 1. Single cylinder diesel engine.
EngineSOKAN
ModelSK-MDF300
Compression ratio20:1
Number of cylinders1
Fuel injection systemDirect injection
Cylinder stroke/bore [mm]63/78
Intake systemNaturally Aspirated
Volume [cc]300
Cycle4 Strokes
Table 2. Properties of the engine lubricant.
Table 2. Properties of the engine lubricant.
SAE15W-40
Density at 20 °C0.864 g/cm3 (DIN 51757)
Pourpoint−33 °C (ISO 3016)
Viscosity at 40 °C91.76 mm2/s
Flashpoint224 °C (DIN ISO 2592)
Lower explosion limit0.6 vol%
Upper explosion limit6.5 vol%
Table 3. Diesel and hydrogen properties.
Table 3. Diesel and hydrogen properties.
PropertiesUnitsDieselHydroxy Gas
Densitykg/m30.8270.084
Kinematic viscositycSt3.2-
Flash point°C74-
Calorific valueMJ/kg41115
Table 4. Measurement instruments.
Table 4. Measurement instruments.
ParameterInstrumentManufacturerRangeAccuracy
AngleCrankshaft angleBeck Arnley 180-04205–9999 RPM0.03%
Cylinder pressurePiezoelectric transducerKISTLER type 7063-A0–250 bar< ±0.5%
AirflowAir mass sensorBOSCH OE-22680 7J6000–125 g/s1%
Fuel measuringGravimetric meterOHAUS-PA3130–310 g1.5%
TemperatureTemperature sensorType K−200–1370 °C0.1%
PressurePressure sensorKISTLER Type 4067-E0–200 bar1%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abril, S.O.; Piero Rojas, J.; Flórez, E.N. Numerical Methodology for Determining the Energy Losses in Auxiliary Systems and Friction Processes Applied to Low Displacement Diesel Engines. Lubricants 2020, 8, 103. https://doi.org/10.3390/lubricants8120103

AMA Style

Abril SO, Piero Rojas J, Flórez EN. Numerical Methodology for Determining the Energy Losses in Auxiliary Systems and Friction Processes Applied to Low Displacement Diesel Engines. Lubricants. 2020; 8(12):103. https://doi.org/10.3390/lubricants8120103

Chicago/Turabian Style

Abril, Sofia Orjuela, Jhan Piero Rojas, and Eder Norberto Flórez. 2020. "Numerical Methodology for Determining the Energy Losses in Auxiliary Systems and Friction Processes Applied to Low Displacement Diesel Engines" Lubricants 8, no. 12: 103. https://doi.org/10.3390/lubricants8120103

APA Style

Abril, S. O., Piero Rojas, J., & Flórez, E. N. (2020). Numerical Methodology for Determining the Energy Losses in Auxiliary Systems and Friction Processes Applied to Low Displacement Diesel Engines. Lubricants, 8(12), 103. https://doi.org/10.3390/lubricants8120103

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop