The Effects of Ultra-Low Viscosity Engine Oil on Mechanical Efﬁciency and Fuel Economy

: The development of a sustainable powertrain requires improved thermal efﬁciency. Reducing frictional power losses through the use of ultra-low viscosity oil is one of the most effective and economical ways. To assess the potential for efﬁciency enhancement in a new generation of future engines using low-viscosity oils, a technical analysis was conducted based on numerical simulation and theoretical analysis. This study proposes a numerical method coupling the whole multi-dynamics model and lubrication model under mixed lubrication regimes. Then, load distribution was calculated numerically and veriﬁed experimentally. Finally, this paper compares the bearing load and frictional energy loss of the main bearings when using The Society of Automotive Engineers (SAE) 15W40 and SAE 0W20 oil. The results indicate that the application of ultralow-viscosity lubricant can reduce the hydraulic friction loss up to 24%, but the asperity friction loss would increase due to the reduction in load capacity. As a result, the design of a new generation of high efﬁciency internal combustion engines requires careful calculation and design to balance the trade-off relations between hydraulic friction and asperity friction.


Introduction
As the population increases the global demand for energy saving and emission reduction increases. The transport sector is one of the major consumers of energy and reducing fuel consumption of internal combustion engines (ICE) is the key to cut energy consumption and reduce CO 2 emissions [1][2][3][4][5]. Holmberg and Erdemir [6] also reports that around 20% (103 EJ) of the world's total energy is used to overcome friction between moving parts. By taking advantage of the tribology design and method like lubrication technologies for friction reduction in vehicles, energy losses caused by friction can be reduced by 40% in the next 15 years and 1460 MtCO 2 emissions can also be reduced globally. Since most of the friction losses in the ICE occur at the cylinder inner surface with piston and piston ring mating and crankshaft bearing friction, the exploration of friction reduction for these sliding friction pairs is mainly in the areas of lubrication additives, application of advanced coating and surface technologies, precise lubrication theory and simulation techniques. Other works performed an overview of the combustion technologies that have been tested able to improve the CO 2 and the NOx-Soot trade-offs such as specific bowl design, innovative fuel injection systems [7][8][9][10] and clean combustion technologies [11,12]. Studies have shown that reducing oil viscosity can be effective in improving fuel efficiency, but the problem that accompanies the application of lower viscosity oils is the accentuation of wear problems caused by the inadequate lubrication carrying capacity of lower viscosity oils [13]. The use of low viscosity lubricants in the automotive sector has been studied for a long time; however, increasingly stringent legislation is forcing the automotive industry to The above Stribeck curve shows that the lubrication state gradually transits from boundary lubrication to hybrid lubrication and dynamic lubrication with the increase of dimensionless film thickness. As we known that multiple engine main crankshaft bearings usually support one crankshaft under the condition of being heated by the combustion chamber; therefore, the position and structure of bearings and bearing seats are strictly limited owing to the structure and working condition of the ICE. In addition, the flexibility of the main components of the internal combustion engine, such as crankshaft and engine block, should be considered under instantaneous high explosive pressure load working condition. Additionally, the difference of tribological properties between the bearings cannot be ignored considering the uneven stresses and the interaction force of them. Thus, a better design and analysis of bearings, including the lubrication and interaction between the friction surfaces of multiple bearings, are of critical importance.
The load capacity, friction loss and wear condition of the bearings are determined by the lubricant and the viscosity of which plays a leading role in the process of hydrodynamic lubrication. For the two surfaces of the friction pair separating by the lubricating oil film, the thickness of the oil film and the viscosity of the lubricating oil play a leading role in the friction loss and under the condition of mixed lubrication or even boundary lubrication, the nominal gap between the two surfaces is very small and the friction caused by the contact of rough peaks between the two sliding surfaces is more significant. The FEA simulation of internal combustion engine main bearing based on low viscosity oil needs to consider the following aspects: 1. as the internal combustion engine develops requiring higher explosion pressure and lighter weight, the influence of oil film local pressure on elastic deformation increases; 2. the viscosity-pressure effect and shear thinning effect on lubrication friction cannot be ignored for the local force of bearing oil film getting higher with high-speed shear force applying to extremely thin oil film; 3. mixed lubrication and even metal contact (side load) effect of main bearing exists caused by the dynamic load and complex working conditions of the internal combustion engine; 4. the influence of the surface quality and texture direction of the friction pair on the lubrication flow cannot be ignored for the oil film thickness and the roughness of the rough surface of the friction pair are being in the same order of magnitude; The above Stribeck curve shows that the lubrication state gradually transits from boundary lubrication to hybrid lubrication and dynamic lubrication with the increase of dimensionless film thickness. As we known that multiple engine main crankshaft bearings usually support one crankshaft under the condition of being heated by the combustion chamber; therefore, the position and structure of bearings and bearing seats are strictly limited owing to the structure and working condition of the ICE. In addition, the flexibility of the main components of the internal combustion engine, such as crankshaft and engine block, should be considered under instantaneous high explosive pressure load working condition. Additionally, the difference of tribological properties between the bearings cannot be ignored considering the uneven stresses and the interaction force of them. Thus, a better design and analysis of bearings, including the lubrication and interaction between the friction surfaces of multiple bearings, are of critical importance.
The load capacity, friction loss and wear condition of the bearings are determined by the lubricant and the viscosity of which plays a leading role in the process of hydrodynamic lubrication. For the two surfaces of the friction pair separating by the lubricating oil film, the thickness of the oil film and the viscosity of the lubricating oil play a leading role in the friction loss and under the condition of mixed lubrication or even boundary lubrication, the nominal gap between the two surfaces is very small and the friction caused by the contact of rough peaks between the two sliding surfaces is more significant. The FEA simulation of internal combustion engine main bearing based on low viscosity oil needs to consider the following aspects: as the internal combustion engine develops requiring higher explosion pressure and lighter weight, the influence of oil film local pressure on elastic deformation increases; 2.
the viscosity-pressure effect and shear thinning effect on lubrication friction cannot be ignored for the local force of bearing oil film getting higher with high-speed shear force applying to extremely thin oil film; 3.
mixed lubrication and even metal contact (side load) effect of main bearing exists caused by the dynamic load and complex working conditions of the internal combustion engine; 4.
the influence of the surface quality and texture direction of the friction pair on the lubrication flow cannot be ignored for the oil film thickness and the roughness of the rough surface of the friction pair are being in the same order of magnitude; Energies 2021, 14, 2320 4 of 20 5. the distribution of lubricating oil film and the lubrication friction performance will be affected by which the main bearing will tilt in the main bearing for its own flexible deformation when the combustion pressure in cylinder breaking out; 6. the application of low viscosity oil will have a nonlinear effect on the combination of the above factors.
Therefore, a simulation method for the main bearings of ICE under the application scenario of low viscosity oil, with a multibody dynamicsmodel of crankshaft and multiple main bearings, considering the influence of asperity contact model and rough surface morphology characteristics must be established to evaluate the tribological characteristics of the ICE accurately. This work can be useful for the design of a new generation of high efficiency ICE and evaluation of the energy and economic benefits of the reduced friction losses from the use of ultra-low viscosity oil.

Test Procedure
This section aims to implement an engine bench test under properly selected working conditions. The data-such as the cylinder pressure, oil supply pressure and rotational speed, together with the ambient temperature in a large range of working conditionscomprising typical loads were collected and were used to gain the input boundary conditions of the simulation, since the accuracy of the kinetic model calculations strongly dependent on the dynamic boundary conditions.

The Object Engine and Overall Arrangement
For the experiment, a six-cylinder diesel engine is selected as the research object and its the main technical parameters are shown in Table 1. The test bench and the overall arrangement of the equipment are shown as in Figure 2: A dynamometer is connected to the rear end of the engine; signals of the cylinder pressure, the piston top dead point position, the oil supply pressure, the rotational speed and torque of the crankshaft are transmitted to the data acquisition system, and the test bench are manipulated by a numerical control test system. The properties of the devices are shown in Table 2.

Test Procedures and Operational Conditions
The bench test was conducted under ambient circumstances and the selected typical operating conditions are shown in Table 3; the test procedure was also conducted accordingly. First, inject the fuel through the numerical control system to the combustion chamber to start the engine and drive the dynamometer, examine the output data of the engine speed, the torque and the temperature of the outlet of the engine cooling water; ensure that the data displayed are within the reasonable range and stable, then start the testing process.
Adjust engine speed to the rated speed (1900 r/min), the combustion peak pressure to 25 MPa; record when the data are stabilized. In addition, the lower the engine speed to 1000 r/min, the minimum speed on full load; record when the data are stabilized, too. Other subsequent working conditions are shown in the order in Table 3. The recorded combustion curve with the peak pressure of 25 MPa and the engine calibration was drawn to curves.

NO.
Peak Pressure (MPa) The Rotational Speed (r/min) 1 25 1900 A dynamometer is connected to the rear end of the engine; signals of the cylinder pressure, the piston top dead point position, the oil supply pressure, the rotational speed and torque of the crankshaft are transmitted to the data acquisition system, and the test bench are manipulated by a numerical control test system. The properties of the devices are shown in Table 2.

Test Procedures and Operational Conditions
The bench test was conducted under ambient circumstances and the selected typical operating conditions are shown in Table 3; the test procedure was also conducted accordingly. First, inject the fuel through the numerical control system to the combustion chamber to start the engine and drive the dynamometer, examine the output data of the engine speed, the torque and the temperature of the outlet of the engine cooling water; ensure that the data displayed are within the reasonable range and stable, then start the testing process. Adjust engine speed to the rated speed (1900 r/min), the combustion peak pressure to 25 MPa; record when the data are stabilized. In addition, the lower the engine speed to 1000 r/min, the minimum speed on full load; record when the data are stabilized, too. Other subsequent working conditions are shown in the order in Table 3. The recorded combustion curve with the peak pressure of 25 MPa and the engine calibration was drawn to curves.

Results
During the endurance test, the measured cylinder pressure curves are depicted in Figure 3. Other critical parameters, such as torque, power and oil supply, pressures are shown in Table 4. Engine calibration 2150

Results
During the endurance test, the measured cylinder pressure curves are depicted in Figure 3. Other critical parameters, such as torque, power and oil supply, pressures are shown in Table 4.

Governing Equations
The dynamic equation of flexible multi-body considering both translation and rotation is described as follows: where the generalized coordinates are

Governing Equations
The dynamic equation of flexible multi-body considering both translation and rotation is described as follows: where the generalized coordinates are q = (q 1 , q 2 , · · · · · · , q n ) , f * is the generalized force of the interaction between components; f a is the external load; p * is nonlinear inertial force and moment and its expression is as follows: where f g is the force and moment produced by the acceleration of the rigid body; f is the centrifugal gyro force and moment.
In the process of modeling, due to the complexity of the actual structure of cylinder liner piston connecting rod crankshaft and engine block, if each node is calculated after the finite element meshing, it is required to establish a huge equation group, which will lead to a huge amount of calculation. So, to simplify the equations, it is necessary to reduce the degrees of freedom of some nodes which have little influence on the calculation results. The commonly used CMS method and the substructure mode synthesis method are used to reduce the model. The transformation involved in the reduction process is as follows: where q is the reduced front node DOF vector; q a is the reduced degree of freedom; G f a is the transformation matrix. By substituting (5) into (1), the reduced form of the kinetic equation is The above formula can be simplified as Simplify the Navier-Stokes equations in the direction of oil film thickness; the Reynolds equation is used in the circumferential and axial direction of the bearing: The lubricating oil film expands along the circumference from the vertex of the bearing into a long square area, as shown in Figure   The two interacting surfaces of sliding bearing pair are rough in the micro level after machining and the oil film thickness can reach the magnitude of the pair surface roughness under heavy load. In order to consider the influence of rough peaks and valleys on the flow in lubricating film, the Patir and Cheng average Reynolds equation of elastohydrodynamic lubrication model considering micro dynamic lubrication effect is adopted, its expression is as follows: In the above expression, lubrication model considering micro dynamic lubrication effect is adopted, its expression is as follows: In the above expression, Here x, z are the circumferential and axial coordinates; is film pressure and film thickness; φ x , φ y is the pressure flow factor along the x, z direction; φ s is the shear flow factor, used to characterize the effect of rough surface on flow; θ = θ(x, z, t) is the filling rate of lubricating oil; u 1 and u 2 are the circumferential velocity of shaft and bearing; t is time; h T is the mathematical expectation of actual oil film thickness under normal distribution. The oil film pressure at both ends of the bearing is considered equal to the gas pressure of the crankcase and the oil film pressure at the oil hole position is considered equal to the oil supply pressure. The Jakobsson-Floberg-Olsson (JFO) method is used to consider capitation. The boundary condition at the oil film rupture is The boundary condition at the oil film formation is The axial boundary condition is Periodic boundary condition is where p c is the hole pressure; V n is the normal velocity of the interface; p e is environmental pressure.

Oil film thickness relation expression
In the calculation of actual oil film thickness, the effects of surface roughness, elastic deformation and wear are considered: where h is the nominal oil film thickness, which is the distance between the nominal dimensions of two rough surfaces; δ T is thermal deformation; δ is the film thickness caused by roughness; δ p is the elastic deformation of the node on the surface. According to the displacement translation equation and the principle of minimum potential energy in elasticity, the node deformation is established to solve the elastic displacement: where F is the resultant force vector of the node external force; λ is the node displacement array; E is the modulus of elasticity; D is a shape function that can be transformed into nodes; Ω is the finite element region.

Viscosity model of lubricating oil
Because the working condition of internal combustion engine changes frequently, the viscosity of lubricating oil is sensitive to temperature to adapt to different working conditions, such as lubrication and heat dissipation function. The Vogel equation is often used as viscosity temperature equation of lubricating oil [25]. The parameters are shown in Table 5, according to the model: The pressure in the oil film increases and the intermolecular distance decreases; the intermolecular force of the lubricating oil increases; this leads to the increase of viscosity. In the future, high speed and heavy load development tendency will increases the local oil film pressure to be as high as 200 MPa and the viscosity pressure effect of the lubricating oil will be very significant at the same time. The Barus equation is often used to express the viscosity pressure relationship of lubricating oil: Under the dynamic load and high running speed, the shear rate in the bearing may exceed 2 × 10 7 1/s. With the increase of the shear rate, the shear thinning effect will occur, resulting in the obvious decrease of the lubricant viscosity. The Cross equation can be used to fit test data: For the properties' parameters of standard oil SAE 0W20, refer to the Table 6 below [26], "HTHS" represents high temperature, high sheer rate. When only viscosity temperature relationship is considered, the comparison of viscosity temperature relationship between SAE 15W40 and SAE 0W20 standard oil depicts in Figure 5, using data from the standard oil database in AVL Excite software. When only viscosity temperature relationship is considered, the comparison of viscosity temperature relationship between SAE 15W40 and SAE 0W20 standard oil depicts in Figure 5, using data from the standard oil database in AVL Excite software.  The comparison of viscosity pressure and shear thinning effect of SAE 0W20 standard oil is shown in Figure 6: The comparison of viscosity pressure and shear thinning effect of SAE 0W20 standard oil is shown in Figure 6: The comparison of viscosity temperature relationship between SAE 15W40 and SAE 0W20 standard oil depicts in Figure 5, using data from the standard oil database in AVL Excite software.

Rough contact model
After multiple machining processes, the surface of sliding bearing pair cannot be completely smooth and the dynamic pressure oil film thickness is so small under heavy load conditions and even reaches the order of roughness (several microns), that the rough contact of rough surface cannot be ignored. The contact pressure between rough peaks can be calculated by Greenwood/Tripp model [8]: The comparison of viscosity temperature relationship between SAE 15W40 and SAE 0W20 standard oil depicts in Figure 5, using data from the standard oil database in AVL Excite software.

Rough contact model
After multiple machining processes, the surface of sliding bearing pair cannot be completely smooth and the dynamic pressure oil film thickness is so small under heavy load conditions and even reaches the order of roughness (several microns), that the rough contact of rough surface cannot be ignored. The contact pressure between rough peaks can be calculated by Greenwood/Tripp model [8]: where , ν k (k = 1, 2) and η s represents the elastic modulus, Poisson's ratio and the number of roughness peaks in the normal contact area; the number represents the two surfaces of the friction pair; β is the density of the asperity peak; β Curvature radius at the top of asperity. For general engineering surfaces, η s βσ s = 0.04, σ s β = 0.001, σ s is the comprehensive RMS (Root Mean Square) roughness between two rough surfaces, which can be calculated by the following formula: The ratio of rough contact area in lubrication area is

Friction loss in mixed lubrication condition
The friction loss under mixed lubrication can be obtained by integrating the oil film shear stress and rough peak shear stress along the whole bearing surface [27]: According to the shear flow factor and pressure flow factor introduced by Equation (6), the oil film shear stress can be expressed as where "+" and "−" represent the surface of bearing bush and bearing surface respectively and the contact shear stress of asperity is: The friction coefficient of boundary lubrication is selected according to reference [10] µ Bound = 0.02.

Engine parameters, material properties and boundary conditions
In the modeling process, the model parameters, material properties and dynamic boundary conditions of the research object are shown in Table 1.
Material parameters of bearing pair can be seen in Table 7. From the free end to the flywheel end, there is No. 1 to No. 7 bearings. To help understand the structure of the research object, the mesh models were established by the Hypermesh software of the engine block, cylinder liners, bearing seats, main bearings and the crankshaft are shown in Figure 7. Bearing clearance (mm) 0.051

Calculation condition, measured loading and oil supply condition
For the crankshaft bearing nonlinear dynamic model requires high boundary conditions, the dynamic parameters measured by laboratory bench test are mostly used as the boundary conditions under typical working conditions, as shown in Table 8. Solving the above equation, the Reynolds equation, energy equation and the heat conduction equation of bearing bush simultaneously, the lubrication between the crankshaft main journal and main bearing experience hydrodynamic lubrication, mixed lubrication and boundary lubrication. According to the ratio of the distance between the nodes of the shaft bearing surface and the comprehensive roughness of the two surfaces (h/σ), the different states of the shaft bearing are distinguished and calculated by the corresponding model. The Reynolds equation is solved by the finite difference method and the elastic displacement equation and the energy equation are solved by the finite element method. The whole calculation process is shown in Figure 8 and the calculation control accuracy of each parameter is shown in Table 9.

Calculation condition, measured loading and oil supply condition
For the crankshaft bearing nonlinear dynamic model requires high boundary conditions, the dynamic parameters measured by laboratory bench test are mostly used as the boundary conditions under typical working conditions, as shown in Table 8. Solving the above equation, the Reynolds equation, energy equation and the heat conduction equation of bearing bush simultaneously, the lubrication between the crankshaft main journal and main bearing experience hydrodynamic lubrication, mixed lubrication and boundary lubrication. According to the ratio of the distance between the nodes of the shaft bearing surface and the comprehensive roughness of the two surfaces (h/σ), the different states of the shaft bearing are distinguished and calculated by the corresponding model. The Reynolds equation is solved by the finite difference method and the elastic displacement equation and the energy equation are solved by the finite element method. The whole calculation process is shown in Figure 8 and the calculation control accuracy of each parameter is shown in Table 9.

Model Accuracies Analyses
In the lubrication numerical calculation, the number of grids plays an important role in the simulation results. In order to obtain the grid independence results and verify the grid independence, the selected grids are (number of axial grids × number of circumferential grids) 12 × 120, 18 × 120, 24 × 120, 24 × 240, 48 × 240, and 24 × 480; the maximum oil film pressure and minimum oil film thickness calculated by each meshing method are shown in Figure 9.
Load, geometry and material properties and initial boundary condition input

Model Accuracies Analyses
In the lubrication numerical calculation, the number of grids plays an important role in the simulation results. In order to obtain the grid independence results and verify the grid independence, the selected grids are (number of axial grids × number of circumferential grids) 12 × 120, 18 × 120, 24 × 120, 24 × 240, 48 × 240, and 24 × 480; the maximum oil film pressure and minimum oil film thickness calculated by each meshing method are shown in Figure 9. It can be seen from the result that with the increase of mesh density, the calculation results tend to be stable and the error is less than 1% after 24 × 240 groups. Considering the calculation accuracy and efficiency, the 24 × 240 mesh a method is adopted in the further simulation.
The simulation model of this work is established by coupling the kinetic model of the engine and lubrication model of the lubrication oil film. The accuracy of the whole model is dependent on the accuracy of each model. For the engine kinetic model, the calculated results are strongly influenced by the dynamic boundary conditions, such as the rotating speed and cylinder pressure, so the measured data for the boundary conditions from a real engine bench test under typical working conditions is used. In the experiment, the speed signal of the crankshaft is measured and controlled by the closed-loop method by the dynamometer with an accuracy of ≤±1 r/min as shown in Table 2; the dynamic load transmitted from the combustion cylinder to the main bearing calculated by the kinetic simulation model is verified by the bearing interaction force obtained by the quasi-static simple beam theory of structural mechanics, the comparison of the two curves in Figure  10 shows that the dynamic load curve and the theoretical results have similar trends, except that the dynamic simulation results fluctuate slightly due to explosive pressures from It can be seen from the result that with the increase of mesh density, the calculation results tend to be stable and the error is less than 1% after 24 × 240 groups. Considering the calculation accuracy and efficiency, the 24 × 240 mesh a method is adopted in the further simulation.
The simulation model of this work is established by coupling the kinetic model of the engine and lubrication model of the lubrication oil film. The accuracy of the whole model is dependent on the accuracy of each model. For the engine kinetic model, the calculated results are strongly influenced by the dynamic boundary conditions, such as the rotating speed and cylinder pressure, so the measured data for the boundary conditions from a real engine bench test under typical working conditions is used. In the experiment, the speed signal of the crankshaft is measured and controlled by the closed-loop method by the dynamometer with an accuracy of ≤±1 r/min as shown in Table 2; the dynamic load transmitted from the combustion cylinder to the main bearing calculated by the kinetic simulation model is verified by the bearing interaction force obtained by the quasi-static simple beam theory of structural mechanics, the comparison of the two curves in Figure 10 shows that the dynamic load curve and the theoretical results have similar trends, except that the dynamic simulation results fluctuate slightly due to explosive pressures from other cylinders of this multi-cylinder engine. The reliability of the simulation results is obtained from the consistency of the two curves above. applicability of which has been confirmed by many scholars, but the accuracy of the method is strongly related to the mesh size of the oil film, so we conducted the grid independence analysis for several mesh patterns, as shown in Figure 9. Furthermore, the load boundary condition of the Reynolds equation uses the force of the bearing joint obtained by the above kinetic model as the input force boundary condition for the oil film solution, the initial pressure of the solution is obtained from the measured oil supply pressure in the bench test to ensure the validity of the entire coupled simulation model.

Results and Discussion
This section contains simulation results analysis of the crankshaft bearings and economic analysis of profit gained by using low viscosity oil.

Simulation Data Analysis
As shown in Figure 11, the hydraulic friction loss of lower viscosity oil is lower than that of normal viscosity oil in each typical working condition and the decreasing amplitude at 2150 rpm is more obvious with the increase of rotating speed, which reach up to 24%: nearly a quarter. While the low viscosity oil is applied, the bearing capacity of oil film drops, the asperity friction power loss of typical working conditions is multiplied several times, even dozens of times high as when using normal viscosity oil, which may lead to the decrease of hydrodynamic friction power loss by using low viscosity oil is partly or even fully neutralized by asperity friction power loss and the rough contact between two surfaces can also lead to bearing wear, which not only changes the lubrication and friction characteristics, but also causes severe failure of other components. For the fluid lubrication model, the generalized Reynolds equation is adopted, the applicability of which has been confirmed by many scholars, but the accuracy of the method is strongly related to the mesh size of the oil film, so we conducted the grid independence analysis for several mesh patterns, as shown in Figure 9. Furthermore, the load boundary condition of the Reynolds equation uses the force of the bearing joint obtained by the above kinetic model as the input force boundary condition for the oil film solution, the initial pressure of the solution is obtained from the measured oil supply pressure in the bench test to ensure the validity of the entire coupled simulation model.

Results and Discussion
This section contains simulation results analysis of the crankshaft bearings and economic analysis of profit gained by using low viscosity oil.

Simulation Data Analysis
As shown in Figure 11, the hydraulic friction loss of lower viscosity oil is lower than that of normal viscosity oil in each typical working condition and the decreasing amplitude at 2150 rpm is more obvious with the increase of rotating speed, which reach up to 24%: nearly a quarter. While the low viscosity oil is applied, the bearing capacity of oil film drops, the asperity friction power loss of typical working conditions is multiplied several times, even dozens of times high as when using normal viscosity oil, which may lead to the decrease of hydrodynamic friction power loss by using low viscosity oil is partly or even fully neutralized by asperity friction power loss and the rough contact between two surfaces can also lead to bearing wear, which not only changes the lubrication and friction characteristics, but also causes severe failure of other components.
To link the typical working conditions of this engine to the endurance cycle, here reports an example of a statistical average benefit of the engine during service. The reduction in hydraulic friction loss at the lower viscosity oil was as much as 24% at 2150 rpm (maximum idle speed) operating conditions, but the benefits are less under other operating conditions. In order to assess the friction reduction effect in a more practical way, the average friction reduction effect is assessed using the endurance cycle. The endurance cycle load spectrum is based on the company's internal design criteria and experience with this model, which involves a statistical value for the percentage of operating time for each typical operating condition and can be used as an equivalent to the actual service process. The friction reduction effect and time share of each typical operating condition in one cycle are shown in the Table 10 below. To link the typical working conditions of this engine to the endurance cycle, here reports an example of a statistical average benefit of the engine during service. The reduction in hydraulic friction loss at the lower viscosity oil was as much as 24% at 2150 rpm (maximum idle speed) operating conditions, but the benefits are less under other operating conditions. In order to assess the friction reduction effect in a more practical way, the average friction reduction effect is assessed using the endurance cycle. The endurance cycle load spectrum is based on the company's internal design criteria and experience with this model, which involves a statistical value for the percentage of operating time for each typical operating condition and can be used as an equivalent to the actual service process. The friction reduction effect and time share of each typical operating condition in one cycle are shown in the Table 10 below. The increase in friction loss due to rough peak friction under the endurance cycle is hydraulic friction reduced multiplied by the corresponding time share under each working condition. To conclude the example, the overall hydrodynamic friction reduction effect of the internal combustion engine according to its endurance cycle is about 20%. Other standard cycles can be applied to assess the applicable benefits given the corresponding tested boundary conditions.
The minimum oil film thickness is the key parameter to characterize the lubrication performance. Figure 12a-c show the variation of the minimum oil film thickness in a period with rotating speed in each cycle of the main bearing when using the lubricant with different viscosity. After comprehensively comparing the minimum oil film thickness in a period of each main bearing under typical conditions when, respectively, using SAE 15W40 ( ) t η , SAE 0W20 ( ) t η and SAE 0W20 ( , , ) t p η γ  , it can be seen that whether  The increase in friction loss due to rough peak friction under the endurance cycle is hydraulic friction reduced multiplied by the corresponding time share under each working condition. To conclude the example, the overall hydrodynamic friction reduction effect of the internal combustion engine according to its endurance cycle is about 20%. Other standard cycles can be applied to assess the applicable benefits given the corresponding tested boundary conditions. The minimum oil film thickness is the key parameter to characterize the lubrication performance. Figure 12a-c show the variation of the minimum oil film thickness in a period with rotating speed in each cycle of the main bearing when using the lubricant with different viscosity. After comprehensively comparing the minimum oil film thickness in a period of each main bearing under typical conditions when, respectively, using SAE 15W40 η(t), SAE 0W20 η(t) and SAE 0W20 η(t, p, . γ), it can be seen that whether considering the pressure-viscosity and shear thinning effects or not, the use of lower viscosity oils (SAE 0W20) causes obvious decrease of minimum oil film thickness in a period under typical working conditions with respect to when using normal viscosity oil (SAE 15W40). When SAE 0W20 oil is used, two conditions of whether the pressure-viscosity and shear thinning effects are included in the viscosity model or not are compared. In addition, it can be concluded that, for the same kind of lubricant, the pressure-viscosity effect caused by the local increase of oil film pressure induces the increase of both the oil viscosity and bearing capacity of the main bearing. The shear thinning effect caused by local strong shear reduces the oil viscosity and the bearing capacity of the main bearing. After applying the new viscosity model which has considered the two effects above, the minimum oil film thickness in a period of main bearing shows a slight increasing trend on the whole, which explains that the bearing capacity of a main bearing is promoted after adopting the new model. It follows that the increase of oil viscosity caused by the pressure-viscosity effect after using low viscosity oil plays a more important role than that of shear thinning effects. local strong shear reduces the oil viscosity and the bearing capacity of the main bearing. After applying the new viscosity model which has considered the two effects above, the minimum oil film thickness in a period of main bearing shows a slight increasing trend on the whole, which explains that the bearing capacity of a main bearing is promoted after adopting the new model. It follows that the increase of oil viscosity caused by the pressure-viscosity effect after using low viscosity oil plays a more important role than that of shear thinning effects.
There are seven main bearings in a 6-cylinder engine, interacting with seven main shaft segments in the same crankshaft. MB1 refers to No.1 main bearing and MB2 refers to No.2 main bearing and so on. To simplify the expression, the name of each main bearing in the subsequent text is represented by an abbreviation "MB", plus the corresponding number. According to the oil film thickness, it can be seen that the main bearings remain in the mixed lubrication state under the rated speed and rated high torque conditions. The asperity contact area ratio between the two surfaces of the friction pair of main bearing under each typical working condition is extracted, as shown in Figure 12. Figure 13a-c, respectively, show the variation of asperity contact percentage in a period of main bearing with rotating speed when adopting different viscosity. After comprehensively comparing the asperity contact percentage in a period of each main bearing under typical conditions when, respectively, using SAE 15W40 ( ) t η , SAE 0W20 ( ) t η and SAE 0W20 ( , , ) t p η γ  , it can be seen that whether or not considering the oil stick after pressure effect and shear thinned effect, after using lower viscosity oil (SAE 0W20), the asperity contact percentage of main bearing under different typical working conditions significantly increase compared to when using normal viscosity oil (SAE 15W40); among them, the asperity contact percentage of 3# and 5#main bearings under typical conditions increases to more than four times. When the SAE 0W20 oil is used, two conditions of whether the pressure-viscosity and shear thinning effects are included in the viscosity model or not are compared. In addition, it can be concluded that, for the same kind of lubricant, the pressure-viscosity effect caused by the local increase of oil film pressure induces the increase of both the oil viscosity and bearing capacity of the main bearing. The shear thinning effect caused by local strong shear reduces the oil viscosity and the bearing There are seven main bearings in a 6-cylinder engine, interacting with seven main shaft segments in the same crankshaft. MB1 refers to No.1 main bearing and MB2 refers to No.2 main bearing and so on. To simplify the expression, the name of each main bearing in the subsequent text is represented by an abbreviation "MB", plus the corresponding number.
According to the oil film thickness, it can be seen that the main bearings remain in the mixed lubrication state under the rated speed and rated high torque conditions. The asperity contact area ratio between the two surfaces of the friction pair of main bearing under each typical working condition is extracted, as shown in Figure 12. Figure 13a-c, respectively, show the variation of asperity contact percentage in a period of main bearing with rotating speed when adopting different viscosity. After comprehensively comparing the asperity contact percentage in a period of each main bearing under typical conditions when, respectively, using SAE 15W40 η(t), SAE 0W20 η(t) and SAE 0W20 η(t, p, . γ), it can be seen that whether or not considering the oil stick after pressure effect and shear thinned effect, after using lower viscosity oil (SAE 0W20), the asperity contact percentage of main bearing under different typical working conditions significantly increase compared to when using normal viscosity oil (SAE 15W40); among them, the asperity contact percentage of 3# and 5#main bearings under typical conditions increases to more than four times. When the SAE 0W20 oil is used, two conditions of whether the pressure-viscosity and shear thinning effects are included in the viscosity model or not are compared. In addition, it can be concluded that, for the same kind of lubricant, the pressure-viscosity effect caused by the local increase of oil film pressure induces the increase of both the oil viscosity and bearing capacity of the main bearing. The shear thinning effect caused by local strong shear reduces the oil viscosity and the bearing capacity of the main bearing. After applying the new viscosity model which has considered the two effects above, the asperity contact percentage in a period of main bearing shows a slight decreasing trend on the whole, which explains that the bearing capacity of main bearing is promoted after adopting the new model. It follows that the increase of oil viscosity caused by the pressure-viscosity effect after using low viscosity oil plays a more important role than that of shear thinning effects, which are consistent with the influence law of minimum oil film thickness in a period. the rear end, the torque output of a multi-cylinder internal combustion engine transmits the torque and accumulates from the front end of the crankshaft to the rear end. Moreover, two thrust bearings are mounted on both sides of the MB6 to limit the axial displacement of the crankshaft. These are likely to result in a greater load on the main bearing No. 6 compared to other bearings and resulting in a relatively lower load carrying capacity for the same oil supply conditions; thus, its cyclic minimum oil film thickness is less than other bearings in all typical operating conditions. To synthesize the specific performance of the low viscosity oil and the mechanism, load capacity parameters, oil film thickness and the asperity contact ratio under typical working conditions were investigated, the pressure-viscosity relation is dominant when low viscosity oil is adopted. A drop in the cyclic minimum film thickness and an increase in the asperity contact ratio indicates that lower viscosity in engine bearing may induce lubrication deterioration and increase friction caused by rough peaks of the joint surfaces.

Economic Analysis
To estimate the economic gain by using ultra low viscosity oil in ICEs, for example, regarding the SAE 0W20 oil, up to 24% of frictional energy consumption will be saved. At present, there are about 1040 million road vehicles in the transportation sector; the energy used in road cars is 83 EJ annually [6] and, of that total, about 3%-4% used to overcome engine bearing friction; about 0.96% of the total energy will be saved, which is 0.8 EJ annually.
The economic impact of friction energy loss and wear related issue can be estimated based on the energy cost of overcoming friction, the cost of the energy required to manufacture wear-related replacement parts and spare equipment and the cost of maintenance To simplify the expression, the name of each main bearing in the subsequent text is represented by an abbreviation "MB" plus the corresponding bearing number, too.
For the six-cylinder engine, No. 6 main bearing (MB6) suffered the worst lubrication condition. Different from the other main bearings, the minimum film thickness drops to near zero, which means asperity contact occurs between the bearing and journal surface. To analyze the reason behind this phenomenon, load characters of the difference for different cylinders. The peak explosive pressure value in each cylinder is almost the same, with a phase difference of 120 • CA (Crank angle). The 7 main bearings supporting the same crankshaft system. With a pulley at the front end of the crankshaft and a flywheel at the rear end, the torque output of a multi-cylinder internal combustion engine transmits the torque and accumulates from the front end of the crankshaft to the rear end. Moreover, two thrust bearings are mounted on both sides of the MB6 to limit the axial displacement of the crankshaft. These are likely to result in a greater load on the main bearing No. 6 compared to other bearings and resulting in a relatively lower load carrying capacity for the same oil supply conditions; thus, its cyclic minimum oil film thickness is less than other bearings in all typical operating conditions.
To synthesize the specific performance of the low viscosity oil and the mechanism, load capacity parameters, oil film thickness and the asperity contact ratio under typical working conditions were investigated, the pressure-viscosity relation is dominant when low viscosity oil is adopted. A drop in the cyclic minimum film thickness and an increase in the asperity contact ratio indicates that lower viscosity in engine bearing may induce lubrication deterioration and increase friction caused by rough peaks of the joint surfaces.

Economic Analysis
To estimate the economic gain by using ultra low viscosity oil in ICEs, for example, regarding the SAE 0W20 oil, up to 24% of frictional energy consumption will be saved. At present, there are about 1040 million road vehicles in the transportation sector; the energy used in road cars is 83 EJ annually [6] and, of that total, about 3-4% used to overcome engine bearing friction; about 0.96% of the total energy will be saved, which is 0.8 EJ annually.
The economic impact of friction energy loss and wear related issue can be estimated based on the energy cost of overcoming friction, the cost of the energy required to manufacture wear-related replacement parts and spare equipment and the cost of maintenance and repair work. The global average price of 1 GJ of energy is 18 EUR used [28], for all the road vehicle the amount of 14,343 million EUR cost will be saved by using SAE 0W20 oil. The energy consumption of 1 PJ results in 0.0683 MtCO 2 carbon dioxide emissions on average global level according to research [6]; by reducing the 0.8 EJ friction energy consumption, the CO 2 emissions can be cut down by 54.6 Mt, which means that up to 3.7% of the estimated carbon dioxide emission in the next 8 years-according to Holmberg's investigation of 1460 MtCO 2 of emissions-can be reduced globally [6]. There are additional potential energy savings if new anti-wear material and coating techniques are applied to Energies 2021, 14, 2320 19 of 20 the surface of the friction pairs to reduce the friction loss caused by asperity friction due to low viscosity lubricants in the future.

Conclusions
In this paper, the working and lubricating properties of the main bearings of ICE with demands of the application of low viscosity oil are analyzed, the economical profits of using low viscosity to improve the efficiency of road cars were estimated as well. The results show that the application of low viscosity oil makes the hydrodynamic friction power loss decrease up to 24%, but the increase of friction loss caused by the aggravation of asperity contact between surfaces may counteract the reducing effect of hydrodynamic friction power loss and the aggravation of asperity contact may also lead to bearing wear.
For the economic profit, the application of low viscosity oil in future high efficiency engine can bring about a saving of up to 14,343 million EUR, cutting down 3.7% of the estimated carbon dioxide emission in the next 15 years. In this study, the method of main bearing lubrication simulation based on low viscosity oil can provide support for the study of the application of low viscosity oil in the internal combustion engine to improve efficiency and can provide references for the technology policy development.
In this paper, we consider rough surface and the flexibility of the crankshaft and the main bearings for low viscosity oil. In the future, we will include the cavitation in the model and analyze all these factors' effects on friction consumption gain, such as oil leakage. More specifically, we will also establish a test rig to verify the simulation results.