Next Article in Journal
Forecasting of 10-Second Power Demand of Highly Variable Loads for Microgrid Operation Control
Previous Article in Journal
Performance Characteristics of A Micro Wind Turbine Integrated on A Noise Barrier
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Salah A. M. Elmoselhy
1,
Waleed F. Faris
2,* and
Hesham A. Rakha
3
1
CFisUC, Department of Physics, University of Coimbra, P-3004 516 Coimbra, Portugal
2
Department of Mechanical Engineering, International Islamic University Malaysia, Gombak, Kuala Lumpur 53100, Malaysia
3
Virginia Tech Transportation Institute, Virginia Polytechnic Institute and State University, 3500 Transportation Research Plaza, Blacksburg, VA 24061, USA
*
Author to whom correspondence should be addressed.
Energies 2021, 14(5), 1287; https://doi.org/10.3390/en14051287
Submission received: 17 November 2020 / Revised: 22 January 2021 / Accepted: 25 January 2021 / Published: 26 February 2021

Abstract

:
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.

Article Highlights

  • Validated Analytical model of the mass flow rate of air goes into cylinders with a flexible crankshaft;
  • Flexibility of the crankshaft decreases the mass flow rate of air that goes into cylinders, resulting in unfavorable higher rate of exhaust emission;
  • Cylinder pressure during intake stroke has four elements: a driving element (intake manifold pressure) and draining elements (vacuum pressure, flow losses and inertial effect of rotating mass).

1. Introduction

The continual demand for diesel powertrain development is due to the distinguished 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.

2. 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.

3. 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:
  • Analytically modeling the mass flow rate of air goes into cylinders with a flexible crankshaft (Section 4);
  • 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);
  • Validating these developed analytical models using case studies and analyzing the simulated results statistically (Section 6);
  • 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;
  • 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;
  • 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.

4. 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].
In order to evaluate the mass flow rate of the air that goes from intake manifold to all the cylinders of the engine, m A C , taking into account the flexibility of the crankshaft let us evaluate the cylinder pressure, PCyl, taking into account the flexibility of the crankshaft. 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:
P C I   =   P i   +   P V a c   +   | P R o t |     P F l o w  
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:
P C y l   =   [ P C I θ I V O θ < θ I V C P i θ   =   θ I V C
The PVac in the diesel engine cylinder can be formulated based on the definition of pressure:
P V a c   =   F a π 4   B 2
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:
F a   =   ( m P   + m C )     a P
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.
  • a P 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 a P indicated in Equation (4), it follows from the definitions of displacement of the piston, Sd, instantaneous linear velocity of the piston, SP, and acceleration of the piston, a P , indicated in Figure 3 that
S d   =   a   cos θ   +   l 2     d 2   2   d   a   sin   θ     a 2   sin 2   θ
where:
  • a is the crank length.
  • θ 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):
S P =   a   ( sin θ   +   d   a   cos θ   + a   sin θ   cos θ   l 2     d 2   2   d   a   sin θ   a 2   sin 2 θ )   d θ d t
However, the crankshaft rotational speed in radian per second, w, in terms of the crankshaft rotational speed in rpm, Nm, can be formulated as:
w   =   2   π   ( N m 60 )
Hence, combining Equations (6) and (7) leads to:
S P =   2   a   π   N m 60   ( sin θ   +   d   a   cos θ   + a   sin θ   cos θ   l 2     d 2   2   d   a   sin θ   a 2   sin 2 θ )  
The acceleration of the piston, a P , can be therefore formulated by differentiation following from Equation (8):
a P =     2   a   π   N m 60     ( cos θ     ( d   a   cos θ   + a   sin θ   cos θ )   ( 2   d   a   cos θ   2   a 2   sin θ   cos θ   ) 2   ( l 2     d 2   2   d   a   sin θ   a 2   sin 2 θ ) 3 2   + d   sin θ   a   sin 2   θ   +   a   cos 2   θ   ( l 2   d 2   2   d   a   sin θ   a 2   sin 2   θ ) 1 2 )   d θ d t
By combining Equations (7) and (9):
a P =     a     π 2     N m 2 900     ( cos θ     ( d   a   cos θ   + a   sin θ   cos θ )   ( 2   d   a   cos θ   2   a 2   sin θ   cos θ   ) 2   ( l 2     d 2   2   d   a   sin θ   a 2   sin 2 θ ) 3 2   + d   sin θ   a   sin 2   θ   +   a   cos 2   θ   ( l 2   d 2   2   d   a   sin θ   a 2   sin 2   θ ) 1 2 )
The mass of the piston, mP, can be formulated following from the definition of mass as:
m P   =   ρ P   V P
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):
F a = ρ P π 3 B 2 Y a N m 2 3600   ( cos θ     ( d   a   cos θ   + a   sin θ   cos θ )   ( 2   d   a   cos θ   2   a 2   sin θ   cos θ   ) 2   ( l 2     d 2   2   d   a   sin θ   a 2   sin 2 θ ) 3 2   + d   sin θ   a     sin 2   θ   +   a   cos 2   θ   ( l 2   d 2   2   d   a   sin θ   a 2   sin 2   θ ) 1 2 )
where:
  • Y is the piston’s characteristic height.
Substituting Equation (12) in Equation (3):
P V a c   =   ρ P   π 2   Y   a   N m 2 900     ( cos θ       ( d   a   cos θ   + a   sin θ   cos θ )   ( 2   d   a   cos θ   2   a 2   sin θ   cos θ   ) 2   ( l 2     d 2   2   d   a   sin θ   a 2   sin 2 θ ) 3 2   + d   sin θ   a     sin 2   θ   +   a   cos 2   θ   ( l 2   d 2   2   d   a   sin θ   a 2   sin 2   θ ) 1 2 )
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]:
P R o t   =   ( m C P   +   m C R   )     a     cos θ     ( N m   2 π 60   ) 2 π   ( B / 2 ) 2
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]:
P 1 ρ   +   c 2 1 2   = P 2 ρ   +   c 2 2 2
Because one of these two points is considered stagnant relative to the other, this equation becomes:
P F l o w   =   ρ ¯     c 2 2  
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:
P F l o w   = K L (   ρ ¯     c 2 2 )  
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:
P C y l   =   [ P i   +   ρ P   π 2   Y   a   N m 2 900   ( cos θ     ( d   a   cos θ   + a   sin θ   cos θ )   ( 2   d   a   cos θ   2   a 2   sin θ   cos θ   ) 2   ( l 2     d 2   2   d   a   sin θ   a 2   sin 2 θ ) 3 2   + d   sin θ   a   sin 2   θ   +   a   cos 2   θ   ( l 2   d 2   2   d   a   sin θ   a 2   sin 2   θ ) 1 2 )   θ I V O   <   θ   θ I V C +   | ( m C P   +   m C R   )   a     cos θ     ( N m   2 π 60   ) 2 π   ( B / 2 ) 2 |     K L     ( ρ     c 2   2 ) P i                                                                                          θ   =   θ I V O
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, m C y l , is analytically formulated in [33]:
m C y l   =     P C y l   η V   V d   N m   n 60   N C y l   R   T i  
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, m A C , can be evaluated as:
m A C   =   [ ( P i   +   ρ P   π 2   Y   a   N m 2 900   ( cos θ     ( d   a   cos θ   + a   sin θ   cos θ )   ( 2   d   a   cos θ   2   a 2   sin θ   cos θ   ) 2   ( l 2     d 2   2   d   a   sin θ   a 2   sin 2 θ ) 3 2   + d   sin θ   a     sin 2   θ   +   a   cos 2   θ   ( l 2   d 2   2   d   a   sin θ   a 2   sin 2   θ ) 1 2 )   +   | ( m C P   +   m C R   )     a     cos θ     ( N m   2 π 60   ) 2 π   ( B / 2 ) 2 |     K L     ( ρ     c 2 2 ) )   η V   V d   N m   n 60   R   T i   θ I V O   <   θ θ I V C P i   η V   V d   N m   n 60   R   T i                                                                          θ   =   θ I V O
In order to simplify this developed model, a sensitivity analysis of which will be conducted in the next section.

5. 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]:
R S   =   ν D C     ν D B ν D B ν I C     ν I B ν I B
where:
  • ν D C is the value of the dependent variable in the analytical model with the increment/decrement of change.
  • ν D B is the baseline value of the dependent variable in the analytical model.
  • ν I C is the value of the independent variable in the analytical model with the increment/decrement of change.
  • ν I B is the baseline value of the independent variable in the analytical model.
The Sensitivity analysis is conducted firstly on the developed model of m A C at the crank angle of opening the intake valve θIVO. The Sensitivity analysis on the developed model of c i m then follows.

5.1. Sensitivity Analysis of the Developed Model of m A C

The mass flow rate of air goes into all cylinders of the engine, m A C , is analytically modeled in Equation (20). The key explanatory variable in this model is N m . As a recommended engine speed during the test for intake stroke, the baseline value of N m is 1550 rpm [26]. The corresponding operating conditions are: P C y l 0.15 × 106 N/m2 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, T i 298 K, and Vd 0.0032 m3 [33]. The sensitivity analysis of the analytical model of m A C 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, m A C , is sensitive to the changes in the value of the independent variables. Thus, the simplified analytical model of m A C remains the same as Equation (20). The simplified analytical model of m A C is hence formulated in Equation (20).

5.2. Sensitivity Analysis of the Developed Model of c i m

The intake manifold air speed dynamics, i.e., acceleration, in diesel engines, c i m in μm/s2, is analytically modeled in another research paper of the authors [19]:
c i m =   P i ρ A     L       2   ζ     c 2 D
where:
  • ρ A is the density of air inside the intake manifold (kg/m3).
  • D L 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: P i 0.1 × 106 N/m2, ζ 0.55, ηV 0.9, T i 298 K, D/L 0.2, and Vd 0.00272 m3 [19,35,36]. The sensitivity analysis of the analytical model of c i m 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, c i m is insensitive to the changes in the value of ζ . Hence, ζ should be considered a constant in the formulation of c i m . The c i m is, however, sensitive to the changes in the value of P i , ρ A , L, and D. Hence, P i , ρ A , L, and D should not be considered constants in the formulation of c i m . The simplified analytical model of the acceleration of air in the intake manifold of diesel engines, c i m , can be therefore formulated as follows:
c i m =     P i ρ A     L     1.1     c 2 D  
The gas speed c indicated in Equation (23) is formulated as:
c   =   Q i m π 4   D 2  
Having seen the simplified formulation of the developed model of c i m , 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 m A C and c i m through a Case Study

The analytical model of m A C developed in Equation (20) and the analytical model of c i m 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.

6.1. Statistically Analyzed Experimental Validation of the Developed Model of m A C

The experimental validation of the analytical model of m A C developed in Equation (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/m3 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.
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, Q i m , of 40 × 10-3 m3/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, m A C , would be evaluated at each of these values of θ. At θ1 = 360° when the intake valve is opened:
m A C = ( 0.15 10 6   N m 2 )   ( 0.9 )   ( 0.0032   m 3 )   ( 1550 2 π 60   rad s ) ( 0.0289 kg mol ) ( 2 ) 60   ( 8.3145   N   m mol   K   )   ( 298   K )     =   0.027   kg / s
At θ2 = 405°:
[ ( 0.15 10 6   N m 2 )   + ( 1439   N m 2 )   +   | 349672   N m 2 |     ( 365742   N m 2 ) ] m A C = ( 0.9 )   ( 0.0032   m 3 )   ( 1550 2 π 60   rad s ) ( 0.0289 kg mol ) ( 2 ) 60   ( 8.3145   N   m mol   K   )   ( 298   K )   =   [ 131052   N m 2 ]   ( 0.9 )   ( 0.0032   m 3 )   ( 1550 2 π 60   rad s ) ( 0.0289 kg mol ) ( 2 ) 60   ( 8.3145   N   m mol   K   )   ( 298   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 m A C inside the cylinders of the engine. Thus, the coefficient of flow loss, KL, becomes 0.25 rather than 1.95 [32].
[ ( 0.15 10 6   N m 2 )   + ( 2.8   N m 2 )   +   |   0   N m 2 |     ( 46890   N m 2 ) ] m A C = ( 0.9 )   ( 0.0032   m 3 )   ( 1550 2 π 60   rad s ) ( 0.0289 kg mol ) ( 2 ) 60   ( 8.3145   N   m mol   K   )   ( 298   K )   =   [ 103107   N m 2 ]   ( 0.9 )   ( 0.0032   m 3 )   ( 1550 2 π 60   rad s ) ( 0.0289 kg mol ) ( 2 ) 60   ( 8.3145   N   m mol   K   )   ( 298   K )   =   0.019   kg / s
At θ4 = 495° when the intake valve is closed:
[ ( 0.15 10 6   N m 2 )   + ( 1439   N m 2 )   +   | 349672   N m 2 |     ( 365742   N m 2 ) ] m A C = ( 0.9 )   ( 0.0032   m 3 )   ( 1550 2 π 60   rad s ) ( 0.0289 kg mol ) ( 2 ) 60   ( 8.3145   N   m mol   K   )   ( 298   K )   =   [ 132491   N m 2 ]   ( 0.9 )   ( 0.0032   m 3 )   ( 1550 2 π 60   rad s ) ( 0.0289 kg mol ) ( 2 ) 60   ( 8.3145   N   m mol   K   )   ( 298   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 statistical measure in this research work is the relative error of the model, ε R , that can be evaluated as follows [48]:
ε R   =   i   = 1 z n | y i     x i y i |   100 % z n
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 comparison of the data of the case study and standard cycles with the results of the developed model of m A C 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.

6.2. Statistically Analyzed Experimental Validation of the Developed Model of c i m

The experimental validation of the analytical model of c i m 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.
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 c i m for the engine is shown in Table 7. The positive sign of c i m indicates accelerated flow. The comparison of the data of the case study with the results of the developed model of c i m is statistically analyzed in Table 7 based on Equation (25).

7. 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, c i m , has been simplified in Equation (23) based on the sensitivity analysis. The c i m has been formulated for the critical moment at which the largest moment is exerted on the crankshaft.
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 m A C has been validated as well through a case study and the ORNL and EPA Freeway standard drive cycles and the ε R on m A C 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 m A C , 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 improve the flexural stiffness and strength of the crankshaft. This would result in longer service life along with exhaust emission reduction [58,59].
The developed model of c i m has been validated as well through a case study and the ε R on c i m has been 11% as indicated in Table 7 which is less than that of widely recognized models in this field such as the CMEM [18]. Equation (23) has implied that the influence of P i , ρ A , and L on c i m is as much as about 90%. Equation (23) has also implied that the influence of ζ , c , and D is as little as about 10%. The developed analytical model of c i m improves on corresponding analytical models presented in key references in this field, such as [60]. Such results can help in realizing an accurate prediction of fuel injection and spray characteristics in diesel engines. Such an accurate prediction results in an accurate prediction of pollutants of these engines [61].

8. Conclusions

The present study has presented and validated the analytical modeling of the intake manifold of a supercharged diesel engine taking into account the nonlinearities exhibited due to the flexibility of the crankshaft. The following achievements have been accomplished in this paper: (I) the study has analytically modeled the mass flow rate of the air that goes into cylinders with a flexible crankshaft; (II) a simplified version for control applications of the analytical models of both the mass flow rate of the air that goes from intake manifold into all the cylinders ( m A C ) and gas speed dynamics of the engine ( c i m ) has been developed based on a sensitivity analysis; (III) the relative errors of the developed models are 7.5% and 11%, which are on average less than that of key models such as the GT-Power and the CMEM, respectively; (IV) validating these developed analytical models using case studies and the ORNL and EPA Freeway standard drive cycles and analyzing the simulated results statistically; (V) the paper captures the influence of the flexibility of the crankshaft on the mass flow rate of air that goes into cylinders, resulting in decrease in the mass flow rate of air by as much as 10% which in turn results in unfavorable higher rate of CO exhaust emissions; (VI) the paper has identified the four elements of the pressure of the gas inside the cylinder during the intake stroke—a driving element (intake manifold pressure), and draining elements (vacuum pressure and flow losses and inertial effect of rotating mass); (VII) 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.

9. Limitation and Future Research

The results of the present research are limited to diesel engines. The validation of these results is limited to the presented data set. Due to the fact that the availability of the field data for validating such research is limited, further validation of the present research results can be provided in future research.

Author Contributions

Conceptualization, S.A.M.E., W.F.F. and H.A.R.; methodology, S.A.M.E., W.F.F. and H.A.R.; analysis, S.A.M.E.; resources, S.A.M.E., W.F.F. and H.A.R.; writing—original draft preparation, S.A.M.E.; supervision and review, W.F.F. and H.A.R.; project administration and funding acquisition, W.F.F. and H.A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding from the Center for Sustainable Mobility, Virginia Polytechnic Institute and State University, USA. Also, this research received funding from the IIUM, Malaysia, grant # RMGS 09-10.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The International Islamic University Malaysia (IIUM) and Virginia Polytechnic Institute and State University (Virginia Tech) are thanked for their support. In addition, the support provided by the University of Coimbra and the Foundation for Science and Technology (FCT) is acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Elmoselhy, S.A.M.; Rakha, H.A.; Faris, W.F. Experimentally validated analytical modeling of diesel engine power and in-cylinder gas speed dynamics. J. Mech. Sci. Technol. 2016, 30, 4725–4734. [Google Scholar] [CrossRef]
  2. Faris, W.F.; Rakha, H.A.; Kafafy, R.I.; Idres, M.; Elmoselhy, S. Vehicle fuel consumption and emission modelling: An in-depth literature review. Int. J. Veh. Syst. Model. Test. 2011, 6, 318. [Google Scholar] [CrossRef]
  3. Faris, W.F.; Rakha, H.A.; Elmoselhy, S.A. Analytical model of diesel engines exhaust NOx emission rate. Int. J. Veh. Syst. Model. Test. 2014, 9, 264. [Google Scholar] [CrossRef]
  4. Faris, W.F.; Rakha, H.; Elmoselhy, S.A. Impact of Intelligent Transportation Systems on Vehicle Fuel Consumption and Emission Modeling: An Overview. SAE Int. J. Mater. Manuf. 2014, 7, 129–146. [Google Scholar] [CrossRef]
  5. Faris, W.F.; A Rakha, H.; Elmoselhy, S. Validated analytical modeling of diesel engine regulated exhaust CO emission rate. Adv. Mech. Eng. 2016, 8, 1–15. [Google Scholar] [CrossRef] [Green Version]
  6. Wojtkowiak, D.; Talaśka, K.; Wilczyński, D.; Górecki, J.; Wałęsa, K. Determining the Power Consumption of the Automatic Device for Belt Perforation Based on the Dynamic Model. Energies 2021, 14, 317. [Google Scholar] [CrossRef]
  7. Bondar, V.; Aliukov, S.; Malozemov, A.; Das, A. Mathematical Model of Thermodynamic Processes in the Intake Manifold of a Diesel Engine with Fuel and Water Injection. Energies 2020, 13, 4315. [Google Scholar] [CrossRef]
  8. Palazzo, N.; Zigan, L.; Huber, F.J.T.; Will, S. Impact of Oxygenated Additives on Soot Properties During Diesel Combustion. Energies 2020, 14, 147. [Google Scholar] [CrossRef]
  9. Guzzella, L.; Onder, C.H. Introduction to Modeling and Control of Internal Combustion Engine Systems, 2nd ed.; Springer: Heidelberg, Germany, 2010. [Google Scholar]
  10. Acquati, F.; Battarola, L.; Scattolini, R.; Siviero, C. An Intake Manifold Model for Spark Ignition Engines. IFAC Proc. Vol. 1996, 29, 7945–7950. [Google Scholar] [CrossRef]
  11. Ouenou-Gamo, S.; Ouladsine, M.; Rachid, A. Air intake theoretical model and simulation of a turbocharger diesel engine. In Proceedings of the Twenty-Ninth Southeastern Symposium on System Theory, Cookeville, TN, USA, 9–11 March 1997; IEEE: Piscataway, NJ, USA, 2002; p. 332. [Google Scholar]
  12. Chen, R.-H.; Wang, W.-C.; Chen, Y.-W. Like-drop collisions of biodiesel and emulsion diesel. Eur. J. Mech. B Fluids 2016, 60, 62–69. [Google Scholar] [CrossRef]
  13. Aguilera-Gonza´ Lez, A.; Bosche, J.; El hajjaji, A.; Abidi, I. Feedback design air-path control on a diesel engine based on Takagi-Sugeno fuzzy descriptor systems. In Proceedings of the Mediterranean Conference on Control & Automation, Platanias-Chania, Crete, Greece, 25–28 June 2013; IEEE: Piscataway, NJ, USA, 2013. [Google Scholar]
  14. Mehta, H.B.; Bharti, O.S. Performance analysis of an irreversible Otto cycle using finite time thermodynamics. In Proceedings of the World Congress on Engineering, London, UK, 29 June–1 July 2009; WCE: London, UK, 2009; Volume II. [Google Scholar]
  15. Kumhar, G.K.; Singh, S.K.; Babhale, T.; Machinery, F.E.A. Design and Analysis of Four Cylinder Diesel Engine Balancer Shaft. Int. J. Eng. Res. 2016, 5. [Google Scholar] [CrossRef]
  16. Chaubey, A.; Tiwari, A.C. Design and CFD analysis of the intake manifold for the Suzuki G13bb engine. Int. J. Res. Appl. Sci. Eng. Technol. (IJRASET) 2017, 5, 1258–1276. [Google Scholar]
  17. Vojtíšek, M.; Kotek, M. Estimation of Engine Intake Air Mass Flow using a generic Speed-Density method. J. Middle Eur. Constr. Des. Cars 2014, 12, 7–15. [Google Scholar] [CrossRef]
  18. Rakha, H.A.; Ahn, K.; Faris, W.F.; Moran, K.S. Simple Vehicle Powertrain Model for Modeling Intelligent Vehicle Applications. IEEE Trans. Intell. Transp. Syst. 2012, 13, 770–780. [Google Scholar] [CrossRef]
  19. Faris, W.F.; Rakha, H.A.; Elmoselhy, S.A. Supercharged diesel powertrain intake manifold analytical model. Int. J. Veh. Syst. Model. Test. 2014, 9. [Google Scholar] [CrossRef]
  20. Alhinti, I.; Samhouri, M.; Al-Ghandoor, A.; Sakhrieh, A. The effect of boost pressure on the performance characteristics of a diesel engine: A neuro-fuzzy approach. Appl. Energy 2009, 86, 113–121. [Google Scholar] [CrossRef]
  21. Ma, X.; You, X.; Wen, B. Multi-body dynamics simulation on flexible crankshaft system. In Proceedings of the 12th IFToMM World Congress, Besançon, France, 18–21 June 2007. [Google Scholar]
  22. Kaya, T.; Temiz, V.; Parlar, Z. Optimum connecting rod design for diesel engines. Trans. Motauto World 2016, 1, 96–100. [Google Scholar]
  23. Jung, M.; Glover, K. Control-oriented linear parameter-varying modelling of a turbocharged diesel engine. In Proceedings of the 2000 IEEE International Conference on Control Applications, Taipei, Taiwan, 2–4 September 2004; Conference Proceedings (Cat. No.00CH37162). IEEE: Piscataway, NJ, USA, 2004. [Google Scholar]
  24. Obert, E.F. Internal Combustion Engines and Air Pollution; Harper & Row Publishers, Inc.: New York, NY, USA, 1973. [Google Scholar]
  25. Taylor, C.F.; Glaister, E. The Internal Combustion Engine in Theory and Practice. J. Appl. Mech. 1961, 28, 316–317. [Google Scholar] [CrossRef] [Green Version]
  26. Mamat, R.; Abdullah, N.R.; Xu, H.; Wyszynski, M.L.; Tsolakis, A. Effect of air intake pressure drop on performance and emissions of a diesel engine operating with biodiesel and ultra low sulphur diesel (ULSD). Renew. Energy Power Qual. J. 2009, 1, 787–794. [Google Scholar] [CrossRef]
  27. Rama, D.D.; Durgaiah, D.R. Fluid Mechanics and Machinery; New Age International: New Delhi, India, 2007; p. 784. [Google Scholar]
  28. Xu, J. Flow analysis of engine intake manifold based on computational fluid dynamics. J. Phys. Conf. Ser. 2017, 916, 012043. [Google Scholar] [CrossRef] [Green Version]
  29. Louvigny, Y.; Vanoverschelde, N.; Janssen, G.; Breuer, E.; Duysinx, P. Dynamic analysis and preliminary design of twin cylinder engines for clean propulsion systems. In Proceedings of the Multibody Dynamics 2009 ECCOMAS Thematic Conference, Warsaw, Poland, 29 June–2 July 2009; European Community on Computational Methods in Applied Sciences: Barcelona, Spain, 2009. [Google Scholar]
  30. Green, D.W.; Perry, R.H. Perry’s Chemical Engineers’ Handbook, 6th ed.; McGraw-Hill: New York, NY, USA, 1984; pp. 3–162. [Google Scholar]
  31. Cengel, Y.A.; Cimbala, J.M. Fluid Mechanics; McGraw Hill: New York, NY, USA, 2006. [Google Scholar]
  32. Kudela, H. Lectures in Fluid Mechanics. Department of Numerical Flow Modeling, Politechnika Wrocławska. 2018. Available online: http://fluid.itcmp.pwr.wroc.pl/~znmp/dydaktyka/fundam_FM/Lecture11_12.pdf (accessed on 19 December 2018).
  33. Faris, W.F.; Rakha, H.; Elmoselhy, S.A. Analytical Modelling of Diesel Powertrain Fuel System and Consumption Rate. SAE Int. J. Mater. Manuf. 2015, 8, 139–152. [Google Scholar] [CrossRef]
  34. United States Environmental Protection Agency EPA. Process for Conducting Probabilistic Risk Assessment; Appendix A. RAGS; Environmental Protection Agency EPA: Washington, DC, USA, 2001; Volume 3, Part A. [Google Scholar]
  35. Edelbrock Pro-Flo XT EFI LS1 Chevy Intake Manifold Systems 2021. Available online: https://www.edelbrock.com/pro-flo-xt-efi-intake-manifold-71373.html (accessed on 18 January 2021).
  36. Eaton M90 Supercharger 2021. Available online: http://www.eaton.com/Eaton/ProductsServices/Vehicle/Superchargers/diesel-supercharger/index.htm#tabs-1 (accessed on 19 January 2021).
  37. Mitsubishi Motors. Mitsubishi Motors Develops New 3.2-liter Direct Injection Diesel. Available online: https://www.mitsubishi-motors.com/en/corporate/pressrelease/corporate/detail724.html (accessed on 18 January 2021).
  38. Nikkar, S. Estimation of in-Cylinder Trapped Gas Mass and Composition. Master’s Thesis, Linköping University, Linköping, Sweden, 2017. [Google Scholar]
  39. Lückert, P.; Schommers, J.; Werner, P.; Roth, T. The new four-cylinder diesel engine for the Mercedes-Benz B-Class. MTZ Worldw. eMagazine 2011, 72, 18–25. [Google Scholar] [CrossRef]
  40. Milašinović, A.; Milovanovic, Z.; Knežević, D.; Mujanić, I. Determination of Differential Equations of Motion and Parameters of an Elastic Internal Combustion Engine Crankshaft. Trans. FAMENA 2016, 40, 83–95. [Google Scholar] [CrossRef] [Green Version]
  41. Bosch Automotive Handbook, 9th ed.; Wiley: Hoboken, NJ, USA, 2014.
  42. Kimura, J.; Yamashita, T. Diesel Engine Crankshaft Rotational Speed Fluctuation Analysis. Des. Appl. Perform. Emiss. Mod. Intern. Combust. Engine Syst. Compon. 2003, 263–269. [Google Scholar] [CrossRef]
  43. De Farias, M.S.; Schlosser, J.F.; Martini, A.T.; Dos Santos, G.O.; Estrada, J.S. Air and fuel supercharge in the performance of a diesel cycle engine. Ciência Rural 2017, 47. [Google Scholar] [CrossRef] [Green Version]
  44. Sulaiman, S.A.; Murad, S.H.M.; Ibrahim, I.; Karim, Z.A.A. Study of Flow in Air-Intake System for A Single-Cylinder Go-Kart Engine. Int. J. Automot. Mech. Eng. 2010, 1, 91–104. [Google Scholar] [CrossRef]
  45. Pesic, R.; Davinic, A.; Petkovic, S.; Taranovic, D.; Miloradovic, D. Aspects of volumetric efficiency measurement for reciprocating engines. Therm. Sci. 2013, 17, 35–48. [Google Scholar] [CrossRef]
  46. Haapakoski, J. Medium-Speed Four-Stroke Diesel Engine Cylinder Pressure Effect on Component Dimensioning. Master’s Thesis, University of Oulu, Oulu, Finland, 2016. [Google Scholar]
  47. Idris, A.; Bakar, R.A.; Ismail, A.R. Engine cylinder fluid characteristics of diesel engine converted to CNG engine. Eur. J. Sci. Res. 2009, 26, 443–452. [Google Scholar]
  48. Keller, G. Statistics for Management and Economics, 9th ed.; South Western, Cengage Learning: Mason, OH, USA, 2012. [Google Scholar]
  49. Yongkang Guihou Industry and Trade Co., Ltd. 150cc 4-Stroke Single Cylinder Go Kart ATV Jeep 2016. Available online: https://emtosales.en.made-in-china.com/product/AKwJSHWOXmru/China-150cc-200-Cc-4-Stroke-Single-Cylinder-Air-Cooled-Go-Kart-ATV-Factory-Price-Jeep-2016.html (accessed on 19 January 2021).
  50. Ibrahim, A.; Dunn, P. Effects of temporal flow acceleration on the detachment of microparticles from surfaces. J. Aerosol Sci. 2006, 37, 1258–1266. [Google Scholar] [CrossRef]
  51. Holmgren, A. Mean Value Modelling of the Intake Manifold Temperature. Master’s Thesis, Linkopings Universitet, Linköping, Sweden, 2005. [Google Scholar]
  52. Wortmann, S.; Vilhena, M.T.; Moreira, D.M.; Buske, D. A new analytical approach to simulate the contaminant dispersion in the PBL. Atmos. Environ. 2005, 39, 2171–2178. [Google Scholar] [CrossRef]
  53. Bos, M. Validation Gt-Power Model Cyclops Heavy Duty Diesel Engine. Master’s Thesis, The Technical University of Eindhoven, Eindhoven, The Netherlands, 2007. [Google Scholar]
  54. Reşitoğlu, İ.A.; Altinişik, K.; Keskin, A. The pollutant emissions from diesel-engine vehicles and exhaust aftertreatment systems. Clean Technol. Environ. Policy 2015, 17, 15–27. [Google Scholar] [CrossRef] [Green Version]
  55. Khajepour, A.; Fallah, S.; Goodarzi, A. Electric and Hybrid Vehicles—Technologies, Modeling and Control: A Mechatronic Approach; Wiley-Blackwell: Hoboken, NJ, USA, 2014; ISBN 978-1-118-40310-5. [Google Scholar]
  56. Lapuerta, M.; Hernández, J.J.; Rodríguez-Fernández, J.; Barba, J.; Ramos, A.; Fernández-Rodríguez, D. Emission benefits from the use of n-butanol blends in a Euro 6 diesel engine. Int. J. Engine Res. 2018, 19, 1099–1112. [Google Scholar] [CrossRef]
  57. Malkawi, G.; Mashayek, F.; Yarin, A.L. Breakup mechanisms of electrostatic atomization of corn oil and diesel fuel. J. Appl. Phys. 2010, 108, 064910. [Google Scholar] [CrossRef] [Green Version]
  58. Elmoselhy, S.A.M. Hybrid Organic/Inorganic Nano-I-Beam for Structural Nano-mechanics. Sci. Rep. 2019, 9, 1–13. [Google Scholar] [CrossRef]
  59. Serrano, J.R.; Arnau, F.; Martín, J.; Auñón, Á. Development of a Variable Valve Actuation Control to Improve Diesel Oxidation Catalyst Efficiency and Emissions in a Light Duty Diesel Engine. Energies 2020, 13, 4561. [Google Scholar] [CrossRef]
  60. Heywood, J. Internal Combustion Engine Fundamentals; McGraw Hill: New York, NY, USA, 1988. [Google Scholar]
  61. Elmoselhy, S.; Faris, W.F.; Rakha, H.A. Experimentally validated analytical modeling of diesel exhaust HC emission rate. J. Mech. Sci. Technol. 2014, 28, 4139–4149. [Google Scholar] [CrossRef]
Figure 1. (a) Deflected crankshaft as indicated by the crank-webs, (b) crankshaft pin before and after deflection due to pressure.
Figure 1. (a) Deflected crankshaft as indicated by the crank-webs, (b) crankshaft pin before and after deflection due to pressure.
Energies 14 01287 g001
Figure 2. Intake manifold of a four-cylinder supercharged diesel engine [19,23].
Figure 2. Intake manifold of a four-cylinder supercharged diesel engine [19,23].
Energies 14 01287 g002
Figure 3. Abstract geometrical illustration of the cylinder, piston, connecting rod and flexible crankshaft.
Figure 3. Abstract geometrical illustration of the cylinder, piston, connecting rod and flexible crankshaft.
Energies 14 01287 g003
Table 1. Sensitivity analysis of the analytical model of m A C at the crank angle of opening the intake valve θIVO.
Table 1. Sensitivity analysis of the analytical model of m A C at the crank angle of opening the intake valve θIVO.
Variables−50% Decrement Below Baseline Value−25% Decrement Below Baseline Value Baseline Value+25% Increment Above Baseline Value+50% Increment Above Baseline Value Average
Sensitivity Ratio on Each Independent Variable
P C y l (N/m2)
(Independent Variable)
0.075 × 1060.11 × 1060.15 × 1060.187 × 1060.225 × 106N/A
m A C (kg/s)
(Dependent Variable)
0.01350.0200.0270.03370.040N/A
RS on
P C y l
11N/A111
η V
(Independent Variable)
0.450.670.91.121.35N/A
m A C (kg/s)
(Dependent Variable)
0.01350.0200.0270.03370.040N/A
RS on η V 11N/A111
V d (m3)
(Independent Variable)
0.00160.00240.00320.0040.0048N/A
m A C (kg/s)
(Dependent Variable)
0.01350.0200.0270.03370.040N/A
RS on V d 0.9970.96N/A1.011.041
T i (K)
(Independent Variable)
149223.5298372.5447N/A
m A C (kg/s)
(Dependent Variable)
0.0540.0360.0270.0220.018N/A
RS on T i −0.5−0.78N/A−1−1.4−0.92
Table 2. Sensitivity analysis of the analytical model of c i m .
Table 2. Sensitivity analysis of the analytical model of c i m .
Variables−50% Decrement Below Baseline Value−25% Decrement Below Baseline Value Baseline Value+25% Increment Above Baseline Value+50% Increment Above Baseline Value Average
Sensitivity Ratio on Each Independent Variable
c i m
(Dependent Variable)
0.0580.0950.131 m/s20.1670.204N/A
P i
(Independent Variable)
0.05 × 1060.075 × 1060.1 × 106 N/m20.125 × 1060.15 × 106N/A
RS on P i 1.11.1N/A1.11.11.1
c i m
(Dependent Variable)
0.2760.1790.131 m/s20.1020.083N/A
ρ A
(Independent Variable)
0.851.291.75 kg/m32.192.55N/A
RS on ρ A −2.2−1.47N/A−0.9−0.7−1.3
c i m
(Dependent Variable)
0.2760.1790.131 m/s20.1020.083N/A
L
(Independent Variable)
0.20.30.4 m0.50.6N/A
RS on L−2.2−1.47N/A−0.9−0.7−1.3
c i m
(Dependent Variable)
0.1380.1340.131 m/s20.1270.124N/A
ζ
(Independent Variable)
0.280.410.550.690.83N/A
RS on ζ 00N/A000
c i m
(Dependent Variable)
0.1170.1260.131 m/s20.1340.136N/A
D
(Independent Variable)
0.0400.0600.080 m0.1000.120N/A
RS on D0.210.2N/A0.10.10.2
Table 3. Field data on a case study conducted with respect to m A C .
Table 3. Field data on a case study conducted with respect to m A C .
ParameterAverage Value Extracted from Field Data during the Intake Stroke
P i 0.15 × 106 N/m2[26]
V d 3.2 × 10−3 m3[26,33]
N m 1550 rpm[26]
T i 298 K[26]
ρ A 1.75 kg/m3[36]
ρ F 832 kg/m3[5]
η V 0.9[33,44]
Q i m 40 × 10−3 m3/s[45]
Table 4. Comparison of the modeling results and the results of standard cycles on m A C .
Table 4. Comparison of the modeling results and the results of standard cycles on m A C .
Standard CyclesStandard Cycles ResultAverage Result of the Analytical Model
ORNL and EPA Freeway LOS
A-C drive cycles (Average Speed 95.52 km/h) [33]
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 m A C is (0.00165 kg/s × 15) that is 0.025 kg/s.The average value of the modeled m A C is 0.024 kg/s.
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 [33,37]. Thus, the mass flow rate of air that goes into all cylinders of the engine m A C is (0.0018 kg/s × 15) that is 0.027 kg/s.The average value of the modeled m A C is 0.024 kg/s.
Table 5. Summary of the statistical analysis on m A C .
Table 5. Summary of the statistical analysis on m A C .
Table x y ε R
Table 4, ORNL and EPA Freeway LOS A-C drive cycles 0.024 kg/s. kg/s0.025 kg/s4%
Table 4, ORNL and EPA Freeway High Speed drive cycle 0.024 kg/s. kg/s0.027 kg/s11%
7.5%
Table 6. Field data on a case study conducted with respect to c i m .
Table 6. Field data on a case study conducted with respect to c i m .
ParameterAverage Value Extracted from Field Data During the Intake Stroke
c i m 0.185 m/s2, on average for low to medium flow’s acceleration [50]
P i 0.155 × 106 N/m2[45]
D L 65 mm / 400 mm, i.e., 0.2[35]
Q i m 40 × 10−3 m3/s[45]
ζ 0.55[45]
ρ A 1.8 kg/m3[36]
ηCA50%[35,36]
Table 7. Summary of the statistical analysis on c i m .
Table 7. Summary of the statistical analysis on c i m .
Table x y ε R
Table 40.21 m/s20.185 m/s211%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Elmoselhy, S.A.M.; Faris, W.F.; Rakha, H.A. Validated Analytical Modeling of Diesel Engines Intake Manifold with a Flexible Crankshaft. Energies 2021, 14, 1287. https://doi.org/10.3390/en14051287

AMA Style

Elmoselhy SAM, Faris WF, Rakha HA. Validated Analytical Modeling of Diesel Engines Intake Manifold with a Flexible Crankshaft. Energies. 2021; 14(5):1287. https://doi.org/10.3390/en14051287

Chicago/Turabian Style

Elmoselhy, Salah A. M., Waleed F. Faris, and Hesham A. Rakha. 2021. "Validated Analytical Modeling of Diesel Engines Intake Manifold with a Flexible Crankshaft" Energies 14, no. 5: 1287. https://doi.org/10.3390/en14051287

APA Style

Elmoselhy, S. A. M., Faris, W. F., & Rakha, H. A. (2021). Validated Analytical Modeling of Diesel Engines Intake Manifold with a Flexible Crankshaft. Energies, 14(5), 1287. https://doi.org/10.3390/en14051287

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