Validated Analytical Modeling of Diesel Engines Intake Manifold with a Flexible Crankshaft

: The flexibility of a crankshaft exhibits significant nonlinearities in the analysis of diesel engines performance, particularly at rotational speeds of around 2000 rpm. Given the explainable mathematical trends of the analytical model and the lack of available analytical modeling of the diesel engines intake manifold with a flexible crankshaft, the present study develops and validates such a model. In the present paper, the mass flow rate of air that goes from intake manifold into all the cylinders of the engine with a flexible crankshaft has been analytically modeled. The analytical models of the mass flow rate of air and gas speed dynamics have been validated using case studies and the ORNL and EPA Freeway standard drive cycles showing a relative error of 7.5% and 11%, respectively. Such values of relative error are on average less than those of widely recognized models in this field, such as the GT-Power and the CMEM, respectively. A simplified version for control applications of the developed models has been developed based on a sensitivity analysis. It has been found that the flexibility of a crankshaft decreases the mass flow rate of air that goes into cylinders, resulting in an unfavorable higher rate of exhaust emissions like CO. It has also been found that the pressure of the gas inside the cylinder during the intake stroke has four elements: a driving element (intake manifold pressure) and draining elements (vacuum pressure and flow losses and inertial effect of rotating mass). The element of the least effect amongst these four elements is the vacuum pressure that results from the piston's inertia and acceleration. The element of the largest effect is the pressure drop that takes place in the cylinder because of the air/gas flow losses. These developed models are explainable and widely valid so that they can help in better analyzing the performance of diesel engines.


Introduction
The continual demand for diesel powertrain development is due to the distin- guished operating features the diesel powertrains exhibit [1].As a cost-effective tool for developing diesel powertrains, modeling assumes a key role in such endeavor [2,3].This is particularly true for Intelligent Transportation Systems (ITS) [4,5].Due to its relative location, intake manifold significantly influences the performance of the rest of the elements of diesel powertrain so that its modeling comes first with respect to diesel powertrain.
The research area of mechanical energy modeling and analysis still has been a hot area of research.This is evident in very recent research papers developed at highly esteemed institutions such as [6,7].In addition, very recently the featured paper in the "Energies" journal-Volume 14-Issue 1 in 2021, has been on diesel engines performance.This indicates that this area of research is still a hot topic that is still investigated by highly esteemed European institutions such as [8].
The intake manifold along with the crankshaft, which are interconnected through the engine cylinder assembly, significantly influence the net power of the engine.Therefore, the modeling of their interconnection becomes critical.Because the flexibility of the crankshaft is an intrinsic characteristic of the dynamics of the crankshaft, the influence of such flexibility should also be modeled.
Guzzella et al. [9] reported that there are two categories of modeling the engine dynamics: (1) Mean Value Model which spreads out major processes in the engine for modeling the input-output behavior with good accuracy but low computational cost; (2) Discrete Event Model which explicitly describes the reciprocating behavior of the engine.However, for the cylinder and intake manifold, the mean value models do not show the influence of the crank angle.In addition, Guzzella et al. did not present a discrete event model of the air mass flow rate in intake manifold and engine's cylinder as a function of crank angle analytically.
Acquati et al. [10] proposed a mean value model of the air flow and pressure in the intake manifold of spark ignition engines.However, the mean value model does not show the influence of the crank angle on the mass flow rate of air in the engine's cylinder and on the gas speed dynamics in the engine's cylinder.
Among the key models in this field is the model developed by Ouenou-Gamo et al. [11] that analytically modeled the air intake process in a turbocharger diesel engine using the camshaft angular position.That study showed the influence of the turbocharger on the manifold and on the cylinders' parameters during their filling.However, that model neglects the gases' kinetic energy and the gas flow's momentum conservation.Other researchers tried to explain the mechanism of fuel consumption that influences the actual Air/Fuel ratio based on understanding the collision behavior of fuel drops in the cylinder of diesel engines [12].
Aguilera-Gonza et al. [13], developed a pressure model based on fuzzy logic for regulating the intake manifold pressure in a four-cylinder diesel engine air-path system with Exhaust Gas Recirculation (EGR) and Variable Geometry Turbocharger (VGT).However, that model did not show analytically the influence of the intake manifold's pressure on the dynamics of the crankshaft.
Considering irreversible thermodynamic processes, Mehta and Bharti [14] presented irreversible Otto cycle model for net power generation and efficiency calculation by including bypass heat leakage, frictional losses and internal irreversibility factors associated with compression and expansion process.They showed that the power is independent of heat leakage term but is affected by irreversibility in terms of flow losses and frictional losses.
The technology of a balance shaft with eccentric weight for offsetting the vibrations of the crankshaft has been also modeled.In this line of modeling, Kumhar et al. [15] have identified the two sources of diesel engine's vibrations: (i) vibrations because of rotating and reciprocating inertial forces that act on the crankshaft and the parts connecting with it, (ii) Torsional vibration of the crankshaft due to combustion pulses.Kumhar et al., proposed a computational model of a balancer shaft with eccentric masses to counter balance these two vibrational forces.However, they have not shown analytically the influence of the intake manifold's pressure on the dynamics of the crankshaft.
Chaubey and Tiwari [16] also developed a computational model but of the intake manifold of a four-cylinder engine.They showed the effects of shape of runners on the intake manifold's output.However, that research did not analytically model the mass flow rate of air in the engine's cylinder nor did it analytically model the gas dynamics in the engine's cylinder.
Vojtíšek and Kotek [17] compared real driving measurement of air flow in gasoline engines as well as in diesel engines with a generic analytical model of the mass flow rate based on volumetric efficiency and throttling.They reported that the relative error of such a model exceeds 10%.However, the paper did not take into account the effect of the flexibility of the crankshaft on the mass flow rate of air in the engine's cylinder.
Rakha et al. [18] indicated that a widely valid and microscopic model of diesel engines is yet sought for evaluating the diesel powertrain technologies.Faris et al. [19] therefore modeled analytically the intake manifold of supercharged diesel engines treating the crankshaft of the engine as an inflexible shaft.However, the flexibility of the crankshaft was not taken into consideration in that modeling.
Albeit vehicle analytical modeling describes the physical phenomena associated with engines with explainable mathematical trends, no analytical model has been developed as yet of the supercharged diesel intake manifold taking into account the flexibility of the crankshaft.The significance of this research is indicated in light of recent research that reflects the importance of the line of research of the effect of gas pressure on diesel engines performance to the research community [20].The significance of the flexibility of the crankshaft captured the attention of several researchers, such as Ma et al. [21], but they opted for computational modeling of the crankshaft based on finite element analysis.The flexibility of the crankshaft in diesel engines exhibits significant nonlinearities.Particularly, at crankshaft rotational speeds of around 2000 rpm or higher, which is the near optimal rotational speed for near optimal mechanical efficiency, the influence of such flexibility on diesel engines performance becomes considerable.

Materials and Methods
The present research is quantitative and is based on a theoretical approach.This research employs exploratory and descriptive techniques.Analytical modeling and simulation are employed in this research.A sensitivity analysis is conducted for simplifying the developed models.Field data of two case studies gathered from literature and the ORNL and EPA Freeway standard drive cycles are used for experimental validation.The experimental data gathered from these two case studies are used as a benchmark for validating the results of the analytical models presented in this paper.For checking the consistency of the developed analytical model with the experimental and field data, the range of the crankshaft rotational angle (θ) is equally and symmetrically divided into a few segments [22].The developed model is thus tested at several values of θ.
The research assumptions in this paper include: (i) Because there is no change in phase in the air flow throughout the compressor and the minimum temperature of the air flow throughout the compressor is far above the critical point of atmospheric air, air is considered an ideal gas; (ii) Because air is considered non viscous fluid, the flow in the intake manifold is non viscous so that the gas flow is considered one-dimensional flow; (iii) The difference in pressure head in the air/gas flow due to elevation in the engine is neglected.

Research Idea and Its Novelty
The key idea of this study is to develop and validate an analytical model of supercharged diesel intake manifolds taking the flexibility of the crankshaft of engines into account.The present research analytically models: (a) the rate of mass flow of air that goes into cylinders with a flexible crankshaft, (b) The intake manifold gas speed dynamics.
The following contributions are accomplished in the present study: 1. Analytically modeling the mass flow rate of air goes into cylinders with a flexible crankshaft (Section 4); 2. A simplified version for control applications of the developed models of both the mass flow rate of the air that goes from intake manifold into all the cylinders of the engine and the gas speed dynamics has been developed based on a sensitivity analysis (Section 5); 3. Validating these developed analytical models using case studies and analyzing the simulated results statistically (Section 6); 4. The simulation and statistical analysis presented in Section 6 reveal that the relative errors are of 7.5% and 11% which are on average less than that of well-known models such as the GT-Power and CMEM, respectively; 5.It has been found that the flexibility of the crankshaft decreases the mass flow rate of air that goes into cylinders, resulting in unfavorable higher rate of exhaust emissions such as CO; 6.It has been also found that the pressure of the gas inside the cylinder during the intake stroke has four elements: a driving element (intake manifold pressure), and draining elements (vacuum pressure and flow losses and inertial effect of rotating mass).
The developed models are simplified in order to suit control applications of supercharged diesel engines.The developed models are widely valid with explainable mathematical trends so that their developed simplified version can help in better analysing and optimizing the performance of diesel engines.

Analytical Modeling of the Mass Flow Rate of Air that Goes into Cylinders with a Flexible Crankshaft
The flexibility of the crankshaft of the engine is illustrated in Figure 1.In Figure 1a, the flexibility of the crankshaft is reflected in a deflection in the crankshaft as indicated by the crank-webs.Figure 1b shows the stressed crank of the crankshaft before and after such a deflection.The present research analyses the intake manifold of four-cylinder supercharged diesel engines equipped with an Electronic Throttle Control (ETC).This set-up configuration is illustrated in Figure 2. Figure 2 depicts an idealized abstraction/pictorial model of the real supercharged diesel engine in an endeavour to analyze the real diesel engine, which is a widely accepted approach [23].Because pressure drop occurs instantaneously, PCyl during the exhaust stroke becomes equalling to the ambient pressure, PAtm [24,25].During the intake stroke, the pressure of the gas inside the cylinder during the intake stroke, PCI, becomes equalling to the intake manifold pressure, Pi.However, the pressure in the cylinder drops during the intake stroke [26].Thus, because of the effect of accelerating the reciprocating piston, a slight pressure drop takes place in the cylinder because of the generated vacuum pressure that results from the piston's inertia and acceleration, PVac, representing the inertial effect of the reciprocating mass [27].In addition, further slight pressure drop takes place in the cylinder due to the air/gas flow losses, PFlow [28].Moreover, the pressure influence due to the inertial effect of the rotating mass, PRot, further affects the in-cylinder pressure [29].It is noteworthy in the intake stroke regarding the PRot that in the range of the crankshaft rotational angle of 360° to 450° the PRot boosts the effect of Pi, as would be indicated in the experimental validation "Section 6".However, in the range of the crankshaft rotational angle of 451° to 495° the PRot due to the centrifugal effect of the rotating mass counteracts the effect of Pi, as would be indicated in the experimental validation "Section 6".Thus, given the variation of the sign of the PRot with the crankshaft rotational angle, in the analytical modelling of the effect of PRot on PCyl, the absolute value of PRot is used.Therefore, the pressure of the gas inside the cylinder during the intake stroke, PCI, can be evaluated as follows: Near the end of the intake stroke, pressure equalizes throughout the runners of the intake manifold so that PCyl becomes equal to Pi.Hence, the in-cylinder pressure, PCyl, can be formulated as follows: The PVac in the diesel engine cylinder can be formulated based on the definition of pressure: where: F ɑ is the force generated because of the piston's acceleration.
B is Cylinder (Bore) diameter.
The Fa indicated in Equation ( 3) can be evaluated following from the definition of force: where: mP is the Piston's mass.mC is the mass of the gas in the clearance volume of the engine cylinder that is equal to the total mass of the gas in the engine cylinder before compression based on the principle of conservation of mass.P a is the piston's instantaneous acceleration.
However, due to the fact that the mP is much greater than the mC with several orders of magnitude, the influence of the mC on the Fa in Equation ( 4), is negligible.The flexibility of the crankshaft is represented in the distance "d" indicated in Figure 3.In an endeavour to evaluate P a indicated in Equation ( 4), it follows from the definitions of dis- placement of the piston, Sd, instantaneous linear velocity of the piston, SP, and acceleration of the piston, P a , indicated in Figure 3 that where: θ is the crankshaft rotational angle.l is the length of connecting rod.The instantaneous linear velocity of the piston, SP, can be thus formulated by differentiation using Equation ( 5): However, the crankshaft rotational speed in radian per second, w, in terms of the crankshaft rotational speed in rpm, Nm, can be formulated as: Hence, combining Equations ( 6) and ( 7) leads to:  (8) The acceleration of the piston, P a , can be therefore formulated by differentiation following from Equation ( 8):  (9) By combining Equations ( 7) and ( 9):  (10) The mass of the piston, mP, can be formulated following from the definition of mass as: where: ρP is the piston' s material density, which is usually aluminium alloy.VP is the volume of the piston.Therefore, the force generated due to the piston's acceleration, Fa, can be evaluated by combining Equations ( 4), (10), and ( 11):  (12) where: Y is the piston's characteristic height.Substituting Equation (12) in Equation ( 3):  (13) The pressure drop due to the inertial effect of the rotating mass, PRot, can be analytically expressed based on the definition of pressure and the geometry indicated in Figure 3 as follows [29]: (14) where: mCP is the mass of the protruded crank portion.mCR is the mass of the connecting rod.
In addition, the pressure drop that takes place in the cylinder due to the air/gas flow losses, PFlow, can be analytically formulated following from the first principles of physics.At 300 K and higher temperature up to 1000 K, the compressibility factor of air remains equaling "1" as long as the pressure of air remains below 55 bar [30].The cylinder's pressure during the intake stroke does not exceed 55 bar [26].Thus, we can treat the air flow inside the intake system of the engine as incompressible fluid.
Let us consider two arbitrary points one and two on the stream of the air flow, one of them is in the intake manifold where the fluid medium is air and the other point is in the cylinder of the engine where the fluid medium is a mixture of air and atomized fuel.This is particularly true in light of the fact that early fuel injection is recommended for improved combustion.The pressure drop between these two points is analytically formulated following from the principle of conservation of energy and Bernoulli's equation that relates the changes in fluid, i.e., air, speed (c) and fluid, i.e., air, pressure (P) along the streamline [31]: Because one of these two points is considered stagnant relative to the other, this equation becomes: where: is the arithmetic mean of fluid density between these two points, averaging the density of air ( A ρ ) and the density of atomized fuel ( F ρ ).
Because the geometry of the intake manifold plenum and runners and the shape of the intake valve opening together affect the fluid flow, this effect is taken into consideration in terms of the loss coefficient, KL, that is introduced into Equation ( 16) as follows: The conservative value of the loss coefficient KL in this regard is about "1.95" for Tees Branch flow of fluid from the intake manifold plenum into the cylinders of the engine [32].Thus, by combining Equations ( 1), ( 2), ( 13), (14), and ( 17) the pressure of the air that goes into cylinder PCyl can be formulated as follows:  (18) where: θIVC is the rotational angle of the crankshaft at which the intake valve closes.θIVO is the rotational angle of the crankshaft at which the intake valve opens.
The mass flow rate of air that goes into one cylinder of the engine with a rigid crankshaft, Cyl m • , is analytically formulated in [33]: where: ηV is the volumetric efficiency of the engine.
Vd is the displaced volume of the engine.n is the number of crank revolutions for each power stroke per cylinder (2 for four-stroke cycles; 1 for two-stroke cycles).
NCyl is the number of cylinders in the internal combustion engine.R is the gas constant.Ti is the temperature of the air of the intake manifold.
Therefore, by combining Equations ( 18) and ( 19), the mass flow rate of the air that goes from intake manifold into all the cylinders of the engine with a flexible crankshaft, AC m • , can be evaluated as:  (20) In order to simplify this developed model, a sensitivity analysis of which will be conducted in the next section.

Sensitivity Analysis of the Developed Models
A sensitivity analysis of the developed models in this study is conducted in this section for simplifying the developed models in order to suit control applications.The sensitivity analysis is based on applying the sensitivity ratio, RS, to each of the independent variables.While evaluating RS for each independent variable.Meanwhile the values of the remaining independent variables remain unchanged from the baseline values.The RS is formulated as follows [34]: N is 1550 rpm [26].The corresponding operating conditions are: Cyl P 0.15 × 10 6 N/m 2 which is equivalent to a boosting pressure at rotational speed 1550 rpm and compressor's adiabatic efficiency of 50% in comparison with the naturally aspirated engines [35,36], V η 0.9, i T 298 K, and Vd 0.0032 m 3 [33].The sensitivity analysis of the analytical model of AC m • at the crank angle of opening the intake valve θIVO is presented in Table 1 following from Equations ( 20) and ( 21).It is noteworthy that the term "N/A" in Table 1 technically means "Not Applicable", that is the cell labeled with "N/A" cannot logically carry a feasible value of the specific variable indicated in the respective horizontal line that simultaneously represents the header of the respective column of this cell.The term "N/A" thus does not mean missing data or unavailability of data.Therefore, the dependent variable, AC m • , is sensitive to the changes in the value of the independent variables.Thus, the simplified analytical model of AC m • remains the same as Equation (20).The simplified analytical model of AC m • is hence formulated in Equation (20).

Sensitivity Analysis of the Developed Model of im c
• The intake manifold air speed dynamics, i.e., acceleration, in diesel engines, im c • in μm/s 2 , is analytically modeled in another research paper of the authors [19]: where: A ρ is the density of air inside the intake manifold (kg/m 3 ).

L D
is the ratio of the diameter of plenum of intake manifold to the length of intake manifold.
ζ is the flow friction coefficient.
The key explanatory variable in this model is c .The baseline value of c is 30 m/s.
The corresponding operating conditions are: i P 0.1 × 10 6 N/m 2 , ζ 0.55, V η 0.9, i T 298 K, D/L 0.2, and Vd 0.00272 m 3 [19,35,36].The sensitivity analysis of the analytical model of im c • is presented in Table 2 following from Equations ( 22) and ( 21).It is noteworthy that the term "N/A" in Table 2 technically means "Not Applicable", that is the cell labeled with "N/A" cannot logically carry a feasible value of the specific variable indicated in the respective horizontal line that simultaneously represents the header of the respective column of this cell.Thus, The gas speed c indicated in Equation ( 23) is formulated as: The analytical model of AC m • developed in Equation ( 20) and the analytical model of presented in Equation ( 23) are validated in this section.The validation would be conducted through experimental case studies and the results would be analyzed statistically.20) is investigated in this section.A case study has been conducted based on the four-cylinder 3.2 L Mitsubishi 4M4 engine [37].Field data that represent operating points of rotational speed of 1550 rpm and compressor's adiabatic efficiency of 50% are shown in Table 3 [36].The engine speed of 1550 rpm is a recommended engine speed during the test for intake stroke [26].This rotational speed is approximately the average speed of such four-cylinder diesel engines [15,38].The length of connecting rod (l) is 0.148 m and the crank length (α) is 0.048 m [39].The piston's material density (ρP) is 2700 kg/m 3 and the piston's characteristic height (Y) is 0.055 m [39].Based on the torsional deformation that reflects the amplitude of the rigid body motion of the center of the crank shaft, the flexibility of the crankshaft (D) is 2 × 10 −4 m [40].This very small value of the (d) corresponds to the reported very small angular deformation of the crankshaft of about 0.2° in average [40].The cylinder (bore) diameter (B) is 0.098 m [37].The summation of the mass of piston, mass of connecting rod and mass of crank pin collectively is 2.95 kg [40].The molar weight of air is 0.0289 kg/mol.

Statistically Analyzed Experimental Validation of the
Mechanical superchargers can operate up to the compressor's impeller speed of 100,000 rpm [41].The supercharger boosting pressure is not constant, since it depends on the speed of the supercharger impeller which in turn is driven by the crankshaft through a belt [42].The faster the supercharger's impeller rotates, the higher the boosting pressure would be.The supercharger efficiency becomes lower during operation, because the crankshaft rotational speed fluctuates and engine's accessories such as the supercharger have moments of inertia, so that belts slip on pulleys.Thus, this results in wear, power loss and reduced efficiency.Thus, in the intake stroke we aim at almost 50% boosting pressure (i.e., Pi = 1.5 bar) at relatively high efficiency of 67% [36,43].
For checking the consistency of the developed analytical model with the experimental and field data, the range of the crankshaft rotational angle (θ) would be equally and symmetrically divided into a few segments [22].Thus, the developed model would be tested at the following values of θ: (1) At θ1 = 360° when the intake valve is opened [46,47], (2) At θ2 = 405°, (3) At θ3 = 450°, (4) At θ4 = 495° when the intake valve is closed [46,47].For a volumetric flow rate of air in the intake manifold per cylinder, im Q , of 40 × 10 -3 m 3 /s the gas speed, c, becomes 30 m/s according to equation (24).Therefore, following from Equation (20) and Table 3, for the four-stroke thermodynamic cycle, the mass flow rate of the air that goes from intake manifold into all the cylinders of the engine, AC m • , would be evaluated at each of these values of θ.At θ1 = 360° when the intake valve is opened: ) (298 K) mol K 0.0243 kg/s = At θ3 = 450°: Within the range of the crankshaft rotational angles in the intake stroke, this crankshaft rotational angle is a singular and inflection point at which the fluid that flows from the manifold plenum into the cylinders of the engine switches from Tees Branch flow to Tees Line flow.This results in even further reduction of the AC m • inside the cylinders of the engine.Thus, the coefficient of flow loss, KL, becomes 0.25 rather than 1.95 [32].
)(298K) mol K 0.019 kg/s = At θ4 = 495° when the intake valve is closed: N m 60 (8.3145 ) (298 K) mol K 0.0246 kg/s = In order to validate this analytical model, the modeling results are compared with the corresponding values of the standard ORNL and EPA Freeway LOS A-C drive cycles.This comparison is illustrated in Table 4.The Air/Fuel ratio is 15 [33].Based on the four-cylinder 3.2 L Mitsubishi 4M4 engine, the fuel mass flow rate as read from field data is 0.00165 kg/s [33,37].Thus, the mass flow rate of air that goes into all cylinders of the engine ORNL and EPA Freeway High Speed drive cycle (Average Speed 101 km/h) [32] The Air/Fuel ratio is 15 [33].Based on the four-cylinder 3.2 L Mitsubishi 4M4 engine, the fuel mass flow rate as read from field data is 0.0018 kg/s The average value of the modeled AC m • is 0.024 kg/s.[33,37].Thus, the mass flow rate of air that goes into all cylinders of the engine AC m • is (0.0018 kg/s × 15) that is 0.027 kg/s.
The statistical measure in this research work is the relative error of the model, R ε , that can be evaluated as follows [48]: where: xi is the interval variable which is the analytically expected value.yi is the interval variable that is the measured value, i.e., field dataset.zn is the number of records in the sample.
The R ε is utilized in this study as a quantifiable measure of validation.The com- parison of the data of the case study and standard cycles with the results of the developed model of AC m • is statistically analyzed in Table 5. where:

−
x is the average value of the interval variable which is the analytically expected value.
− y is the average value of the interval variable that is the measured value, i.e., field dataset.

Statistically Analyzed Experimental Validation of the Developed Model of
The experimental validation of the analytical model of presented in Equation ( 23) is investigated in this section.A case study has been conducted based on the 200 cc Four-Stroke Single Cylinder Go Kart ATV Jeep 2016 [49].Field data that represent an operating point of rotational speed of 2500 rpm and compressor's adiabatic efficiency of 50% are shown in Table 6.This rotational speed is the average speed of such diesel engines.Because the diesel engine is usually operated at wide-open throttle, the corresponding air flow's acceleration in the intake manifold ranges from slow to medium acceleration.ηCA 50% [35,36] where: ηCA is the compressor's adiabatic efficiency.An intercooler is often used at the outlet of the supercharging compressor for maintaining the air temperature drawn into the intake manifold to be on average below 35 °C under severe driving conditions [51].In average the corresponding compressor's pressure ratio at the operating point in the intake manifold becomes about 1.5.
Based on Equations ( 23) and ( 24) and Table 6, the analytical value of

Discussion
The analytical model presented in the present paper favourably facilitates the endeavours of improving the performance of diesel engines towards having environmentally friendly diesel engines.Such an analytical approach to modeling is favourably adopted in automotive and environmental engineering [52].The present paper has thus presented the analytical modeling of the diesel engine intake manifold with a flexible crankshaft.The analytical models presented in this paper have followed from the first principles of physics as a way of further validating the developed models.The paper has started with the analytical modeling of the mass flow rate of air goes into cylinders with a flexible crankshaft indicated in Equation (20).The distance "d" indicated in Figure 3 represents the flexibility of the crankshaft.The nonlinearities brought to the analysis of diesel engine intake manifold due to the flexibility of the crankshaft in diesel engines are apparent in the analytical model of the piston's instantaneous linear velocity, SP, shown in Equation (8).The analytical model of the intake manifold gas speed acceleration, The analytical model developed in Equation (1) indicates that the pressure of the gas inside the cylinder during the intake stroke has a driving element and draining elements.The driving element is the intake manifold pressure.The draining elements are: (i) the vacuum pressure that results from the piston's inertia and acceleration, (ii) the pressure drop that takes place in the cylinder due to the air/gas flow losses, (iii) the pressure influence due to the inertial effect of the rotating mass.The experimental validation "Section 6.1" indicates that the element of the least effect amongst these four elements is the vacuum pressure that results from the piston's inertia and acceleration.The experimental validation section "Section 6.1" also indicates that the elements of the largest effect amongst these four elements is the pressure drop that takes place in the cylinder due to the air/gas flow losses.
The developed model of AC m • has been validated as well through a case study and the ORNL and EPA Freeway standard drive cycles and the R ε on AC m • has been 7.5% as indicated in Table 5 which is less than that of widely recognized models in this field such as the CMEM [18] and and the GT-Power [53].Table 5 and Equation (20) indicate that the influence of PRot on PCyl is larger than the influence of PVac on PCyl.
As proven in Table 5 through the crank angle of the crankshaft, the flexibility of the crankshaft decreases AC m • , resulting in unfavorable higher rate of exhaust emissions.This is proven is several senses.As mentioned in the couple of lines succeeding Equation ( 4), the flexibility of the crankshaft is represented in the distance "d" indicated in Figure 3.As can be gathered from Figure 3, the distance "d" indicated in Figure 3 vanishes at θ = θIVO = 360° (i.e., = 0°), so that the influence of flexibility does not appear.As θ increases therefrom, the distance "d" indicated in Figure 3 increases, so that the influence of flexibility becomes apparent.Thus, following from Section 6.1, by comparing the mass flow rate of air that goes into cylinders in the cases of varying θ, the flexibility of the crankshaft decreases the mass flow rate of air that goes into cylinders.This decrease in the mass flow rate of air negatively influences the actual Air/Fuel ratio increasing the Air/Fuel ratio.This decrease in the mass flow rate of air is in the range of up to 10%, as can be gathered from the results in Section 6.1.Thus, this decrease in the Air/Fuel ratio reduces the output power and exhibits unfavorable higher rate of exhaust emissions, such as CO exhaust emission [54,55].Carbon monoxide (CO) emissions results from incomplete combustion due to incomplete oxidation.The CO concentration is basically dependent on air/fuel mixture and is highest where the air/fuel ratio is less than the stoichiometric concentration, which is classified as "rich mixture" [54].In the rich mixture, due to air deficiency and reactant concentration, not all Carbon atoms can form CO2 molecules, so that the concentration of the CO molecules becomes noticeable.The CO emissions are caused especially at the times of starting and instantaneous acceleration of engine where the rich mixture is required [54].This makes sense and is in accord with the present results since under the transient condition, the operational vibrations become higher which result in higher deflection in the crankshaft which in turn results in reduction in the air/fuel ratio which eventually causes a higher rate of CO emissions.The increase in CO exhaust emission occurs also with bio-diesel fuel [56].The lower value of PCyl due to the crankshaft flexibility when the crank rotational angle becomes larger than θIVO further explains the decrease in the mass flow rate of air.The influence of the flexibility of the crankshaft on the mass flow rate of gas in cylinders, would be even more pronounced in case of bio-fuel/bio-diesel since the organic oil viscosity is several times higher than that of diesel fuel [57].Given the result that the flexibility of the crankshaft decreases the mass flow rate of air that goes into cylinders, resulting in unfavorable higher rate of exhaust emission, strengthening the microstructure of the crankshaft with nano-I-beam is expected to

Figure 1 .
Figure 1.(a) Deflected crankshaft as indicated by the crank-webs, (b) crankshaft pin before and after deflection due to pressure.

Figure 2 .
Figure 2. Intake manifold of a four-cylinder supercharged diesel engine [19,23].In order to evaluate the mass flow rate of the air that goes from intake manifold to all the cylinders of the engine,

Figure 3 .
Figure 3. Abstract geometrical illustration of the cylinder, piston, connecting rod and flexible crankshaft.

ν•
is the value of the dependent variable in the analytical model with the incre- ment/decrement of change.DB ν is the baseline value of the dependent variable in the analytical model.IC ν is the value of the independent variable in the analytical model with the in- crement/decrement of change.IB ν is the baseline value of the independent variable in the analytical model.The Sensitivity analysis is conducted firstly on the developed model of AC m • at the crank angle of opening the intake valve θIVO.The Sensitivity analysis on the developed model of • im c then follows.5.1.Sensitivity Analysis of the Developed Model of AC mThe mass flow rate of air goes into all cylinders of the engine, AC m• , is analytically modeled in Equation(20).The key explanatory variable in this model is m N .As a rec- ommended engine speed during the test for intake stroke, the baseline value of m im c • is insensitive to the changes in the value of ζ .Hence, ζ should be considered a constant in the formulation of im c , sensitive to the changes in the value of i P , A ρ , L, and D. Hence, i P , A ρ , L, and D should not be considered constants in the formulation of im c • .The simplified analytical model of the acceleration of air in the intake manifold of diesel engines, im c • , can be therefore formulated as follows: seen the simplified formulation of the developed model of • im c , let us now investigate the experimental validation of the developed model through a case study and the statistical analysis thereon.6. Statistically Analyzed Experimental Validation of the Developed Models of AC m

•
Developed Model of AC m The experimental validation of the analytical model of AC m • developed in Equation (

•
is (0.00165 kg/s × 15) that is 0.025 kg/s.The average value of the modeledAC m • is 0.024 kg/s.
in Table7.The positive sign of • im c indicates accelerated flow.The comparison of the data of the case study with the results of the developed model of • im c is statistically analyzed in Table7based on Equation(25).

,
has been simplified in Equation(23) based on the sensitivity analysis.The • im c has been formulated for the critical moment at which the largest moment is exerted on the crankshaft.

Table 1 .
Sensitivity analysis of the analytical model of AC m • at the crank angle of opening the intake valve θIVO.

Table 2 .
Sensitivity analysis of the analytical model of im

Table 3 .
Field data on a case study conducted with respect to AC

Table 4 .
Comparison of the modeling results and the results of standard cycles on AC

Table 5 .
Summary of the statistical analysis on AC m • .

Table 6 .
Field data on a case study conducted with respect to

Table 7 .
Summary of the statistical analysis on