An Investigation of the Composition of the Flow in and out of a Two-Stroke Diesel Engine and Air Consumption Ratio

The aim of this research was to investigate the mass, substance and energy flow through two-stroke low speed Diesel engines. For this reason, a zero-dimensional model of the combustion in the engine was developed with a calculated amount and composition of exhaust gases. Due to the large amount of oxygen in the exhaust gases, a ratio of real air consumption and stoichiometric amount of air required for combustion of injected fuel was set. The calculated ratio showed that the engine consumed four times more air than needed for combustion in AFRstoich. In this work, this was called the Air Consumption Factor or Ratio, and has not previously been mentioned in scientific literature. The air consumption ratio is defined as a factor of dry or humid air. To be more comprehensive, a modified diagram of the composition of the flow in and out of a two-stroke fuel injection engine and the cylinder was made.


Introduction
Low speed, two-stroke, turbocharged diesel engines are the most common marine propulsion engines used today.These engines are the most efficient among others, exhibiting 50% efficiency.The remaining 50% of the energy released from combustion of fuel is lost to the atmosphere as waste heat.
At the beginning of the development of the internal combustion engine, it was recognized that its improvement would be a lengthy and expensive process.Even by the late 1930s, there were established methods to calculate a large number of physical processes in the engine.In other words, there were calculation methods developed that enabled the analysis of the influential parameters of the working processes in the engine, and also provided data for the development of new and improved engine construction methods.
The theoretical background of engine process calculations rests on the work done by List [1].Alongside this, the development of methods for the design of the high part of the process, taking into account the increasing number of parameters, can be found in [2].
Due to the large impact of changes to the working fluid on the power and efficiency of the engine, developing calculation methods for the low-pressure part of the process are based on the laws of gas dynamics.The simplest one is called the stationary method, which takes into account the processes in the cylinder and distribution devices, ignoring the changes of the gas in front of and behind the valve (channel), based on the work of List [3] and Hasselgruber [4].
The models that are described in the papers of McAulay [5] and Woschni [6] approach the real processes in the engine.The processes in the cylinders are described by differential equations derived from the law of conservation of energy and mass, and also from the equation of the gas state.Moreover, changes in the properties of the gas due to compressibility and dissociation are taken into account.
Boy [7] described in detail the process in a propulsion Diesel engine.The model describes the processes in the engine cylinder, the process in the turbocharger, passage channels, and intake distributor, by the method of "full-empty".Some background information on two-stroke engines can be found in Heywood [8] and Richard Stone [9].The gas exchange processes are comprehensively treated by Sher [10].The composition of the inlet and outlet flows of a two-stroke engine and the cylinder are shown by Sher [10] and Van Basshuysen et al. [11].An experimental study of the flow pattern inside a model cylinder of a uniflow scavenged two-stroke engine is presented by Sher et al. [12].The velocity field as well as the turbulent parameters were mapped under steady-flow conditions with the aid of a hot-wire anemometry technique.Ravi and Marathe [13] in their work showed a multidimensional prediction of the scavenging characteristics of a homogeneous charge uniflow scavenged two-stroke cycle engine.An engine geometry with the same dimensions as a General Motors EMD 710 engine (General Motors, Detroit, MI, USA) was analyzed.A multidimensional program, Cyclic Analysis of Reciprocating Engines (CARE), has been developed for this purpose.The fluid flow problem is subdivided into global and local parts and the two parts are solved simultaneously.Combustion is treated as a stoichiometric heat release phenomenon.The computer program is used to study the effects of port/valve sizes and timings on the scavenging characteristics of the engine, provided that the pressures at the inlet and exhaust ports are held constant.It is observed that a larger inlet port area and early opening-late closing results in a considerable increase in the scavenge ratio, hence resulting in a higher scavenging efficiency and lower trapping efficiency.An increase in the exhaust port area results in an increase in the scavenge ratio, but a decrease in the trapping efficiency; this gives rise to an optimum exhaust area for a given inlet area at which the scavenging efficiency is maximal.Carlucci at al. [14] presented an analytical model for the estimation of the trapping efficiency according to the Oswald diagram to the molar concentration of carbon dioxide and oxygen at the tailpipe and then according to the mass flow of air and fuel.
Larsen et al. [15] investigated two-stroke diesel machinery for ships, with five varying configurations to explore the trade-off between increased NO x emissions and the reduction in fuel consumption.By implementing a waste heat recovery system through the use of an organic Rankine cycle and also a hybrid turbocharger, the fuel consumption and NO x emissions were lowered by up to 9% and 6.5% respectively.Andreadis et al. [16] used a large two-stroke marine diesel engine, operating at its full load to explore pilot injection strategies using simulations of computational fluid dynamics along with an evolutionary algorithm.Guan et al. [17] used a modular zero-dimensional engine model that was built in MatLab and the Simulink environment, to investigate the operation of a large two-stroke marine diesel engine.Engine shop trial values were compared with the derived performance parameters of the engine, which was simulated at steady state conditions.The study purpose of Varbanets et al. [18] was to record the methods by which the ship's diesel process efficiency could be improved.Under the conditions of the fuel equipment and the normal state of the main diesel system, even load distribution between the cylinders was controlled.From previous research, the authors investigated the possibility of increasing efficiencies of a low speed two-stroke turbocharged main diesel engine operating with waste heat recovery through combined heat and power production [19,20].Spahni et al. [21] in particular deals with new electronically controlled, two stroke, low speed marine engines.Carlucci et al. [22] compared different architectures (single turbocharger, double turbocharger, single turbocharger with an electrically-assisted turbocharger, with intercooler or aftercooler) designed to supercharge an aircraft two-stroke Diesel engine for general aviation and unmanned aerial vehicles characterized by very high altitude operation and long fuel distance.A 1D model of the purposely designed engine has been used to compare the performance of different supercharging systems in terms of power, fuel consumption, and their effect on trapping and scavenging efficiency at different altitudes.In [23], Carlucci et al. provided several guidelines about the definition of design and operation parameters for a two-stroke two banks uniflow diesel engine, turbocharged with two sequential turbochargers and an aftercooler per bank, with the goal of either increasing the engine braking power at take-off, or decreasing the engine fuel consumption at cruise conditions.The engine has been modelled with a 0D/1D modelling approach.Wang et al. [24] evaluated scavenge port designs for a boosted uniflow scavenge direct injection gasoline engine by 3D computational fluid dynamics (CFD) simulations.In order to fulfil the potential of the boosted uniflow scavenged direct injection gasoline (BUSDIG) engine, various scavenge ports were designed with different scavenge port numbers, axis inclination angles and swirl orientation angles, and their effects were evaluated by 3D CFD under different intake pressures and engine speeds.The scavenging process was analysed by its delivery ratio, trapping efficiency, scavenging efficiency and charging efficiency.
This investigation focused on how to improve efficiency of the diesel engine plant using waste energy contained in exhaust gases and scavenged air.For this purpose we needed information about the quantity of mass, substance and energy flow through the engine, and a 0D model of the combustion and exchange process in the engine was developed.Our 0D program was directly coupled with the NIST-REFPROP Version 9.0 Standard Reference Database [25] so that all physical, chemical and thermodynamic data for all used fluids could be used in our calculations.The calculated amount and composition of the exhaust gases showed a good correlation with the measured results of Maschinenfabrik Augsburg Nürnberg (MAN) Diesel & Turbo.Since our results showed that the amount of air that passed through the engine was several times greater than required by stoichiometric combustion, it was decided to establish a new term: the Air Consumption Ratio.At the end of this paper, the parameters of substantive changes in two-stroke engines with fuel injection are strictly defined.Instead of using the displacement volume as a reference volume which is common in literature, we used V max as a reference volume, which is more appropriate for two-stroke engines.For the purpose of connecting the air consumption ratio with two-stroke gas performance parameters, we proposed an improved version of the figure-the composition of the flows in and out a two-stroke engine and its cylinder.We also found that air consumption ratio depended on ratio of the excess air in the cylinder, and the trapping efficiency.

Engine Model Description
The propulsion engine model with the main input and output variables is presented in Figure 1.The fluid flow through the engine is steady, while the observer position is stationary with the respect to the control volume surrounding the ship's engine.Engine inlet consists of fuel and air, and additionally of cylinder oil.Engine outlet consists of combustion products and small amount of emissions. .W is shaft power transmitted to the propeller while .Q 1 and .Q 2 are heats removed from the engine at temperatures T 1 and T 2 [26].
Energies 2017, 10, 805 3 of 19 with the goal of either increasing the engine braking power at take-off, or decreasing the engine fuel consumption at cruise conditions.The engine has been modelled with a 0D/1D modelling approach.Wang et al. [24] evaluated scavenge port designs for a boosted uniflow scavenge direct injection gasoline engine by 3D computational fluid dynamics (CFD) simulations.In order to fulfil the potential of the boosted uniflow scavenged direct injection gasoline (BUSDIG) engine, various scavenge ports were designed with different scavenge port numbers, axis inclination angles and swirl orientation angles, and their effects were evaluated by 3D CFD under different intake pressures and engine speeds.The scavenging process was analysed by its delivery ratio, trapping efficiency, scavenging efficiency and charging efficiency.This investigation focused on how to improve efficiency of the diesel engine plant using waste energy contained in exhaust gases and scavenged air.For this purpose we needed information about the quantity of mass, substance and energy flow through the engine, and a 0D model of the combustion and exchange process in the engine was developed.Our 0D program was directly coupled with the NIST-REFPROP Version 9.0 Standard Reference Database [25] so that all physical, chemical and thermodynamic data for all used fluids could be used in our calculations.The calculated amount and composition of the exhaust gases showed a good correlation with the measured results of Maschinenfabrik Augsburg Nürnberg (MAN) Diesel & Turbo.Since our results showed that the amount of air that passed through the engine was several times greater than required by stoichiometric combustion, it was decided to establish a new term: the Air Consumption Ratio.At the end of this paper, the parameters of substantive changes in two-stroke engines with fuel injection are strictly defined.Instead of using the displacement volume as a reference volume which is common in literature, we used Vmax as a reference volume, which is more appropriate for two-stroke engines.For the purpose of connecting the air consumption ratio with two-stroke gas performance parameters, we proposed an improved version of the figure-the composition of the flows in and out a two-stroke engine and its cylinder.We also found that air consumption ratio depended on ratio of the excess air in the cylinder, and the trapping efficiency.

Engine Model Description
The propulsion engine model with the main input and output variables is presented in Figure 1.The fluid flow through the engine is steady, while the observer position is stationary with the respect to the control volume surrounding the ship's engine.Engine inlet consists of fuel and air, and additionally of cylinder oil.Engine outlet consists of combustion products and small amount of emissions.W  is shaft power transmitted to the propeller while 1 Q  and 2 Q  are heats removed from the engine at temperatures T1 and T2 [26].The unused heat is rejected at the environment temperature T0 and in this case the environment is the atmosphere and sea water.The heat from the cylinder liner, cylinder cover, lubricating oil and The unused heat is rejected at the environment temperature T 0 and in this case the environment is the atmosphere and sea water.The heat from the cylinder liner, cylinder cover, lubricating oil and The first law of thermodynamics can be written as: where .
Q is total heat output, .
W is shaft power and .m i and h i are the mass flow rates and enthalpies.The shaft power in Equation ( 1) is actually technical work.Now, Equation (1) can be rewritten as: As a result of heat transfer to the environment .
Now, the first law of flow processes is: .
If the working fluids, control volume and environment are considered to be united as a single isolated system, then according to the Clausius inequality (second law for all processes), the total amount of the change in entropy is: .
If all processes are reversible, the maximum output power is: .
In the Equation (4) a change in enthalpy shows how much more work would be obtained, when there are no heat exchanges with the environment.Rearranging Equation (4): .
or per unit mass flow: where h − T 0 s, is known as the exergy function ζ, and e x is the exergy expended per kilogram of input or output substance.
In the case of the internal combustion engine or a fuel cell, as shown in Figure 1, if the air and fuel enter at atmospheric temperature and pressure, maximum work from the fuel is obtained if the process is in equilibrium and the combustion products leave the system at atmospheric temperature and pressure.If so, the defined state of equilibrium combustion products with the environment, from the Equation ( 5) per kilogram of reactants, are as follows: Energies 2017, 10, 805 5 of 20 Per kilogram of fuel is: T 0 ,p 0 (7) where: ∆H T 0 ,p 0 is the enthalpy of reaction of combustion or the fuel caloric value.The suffix * means a quantity corresponding to one kilogram of fuel, and without suffix * is an amount per kilogram of products.∆H T 0 ,p 0 would be a lower caloric value if H 2 O in the combustion products was in the vapour state, and the higher caloric value if the H 2 O in the combustion products was in the liquid state.
When the liquid hydrocarbon fuel burns with air and forms carbon dioxide and water, under their full balance with the atmosphere, it is considered as their states with the ambient temperature T 0 and the pressures equal to their partial pressures in the atmosphere.In that case, s * products ≈ s * reactants with [∆H] T 0 ,p 0 is equal to fuel a higher caloric value.It follows that for these fuels, w * tehn,max ≈ [∆H] T 0 ,p 0 .For pure carbon, lower and higher caloric values are the same; for hydrogen, w * tehn,max = 0.823[∆H ] T 0 ,p 0 .For gaseous hydrocarbons, w * tehn.max is slightly lower in percentage than [∆H] T 0 ,p 0 .But since in practice, the full amount of w * tehn.max is very difficult to achieve, the exergy efficiency calculated as η ex= w* tehn /H d will be used for the purpose of comparison.

Data MAN B&W CEAS_ERD (Engine Room Dimensioning)
For Newbuilding 447 "Donat" Tankerska plovidba d.d.Zadar, in the shipyard and diesel engine factory SPLIT, a propulsion engine 6S70MC-C7 (Shipyard Split-Division Diesel Engines, Split, Croatia) was chosen with the following characteristics: ´T p H Δ would be a lower caloric value if H2O in the combustion products was in the vapour state, and the higher caloric value if the H2O in the combustion products was in the liquid state.
When the liquid hydrocarbon fuel burns with air and forms carbon dioxide and water, under their full balance with the atmosphere, it is considered as their states with the ambient temperature T0 and the pressures equal to their partial pressures in the atmosphere.In that case, .But since in practice, the full amount of * tehn.maxw is very difficult to achieve, the exergy efficiency calculated as ηex=w*tehn/Hd will be used for the purpose of comparison.

Data MAN B&W CEAS_ERD (Engine Room Dimensioning)
For Newbuilding 447 "Donat" Tankerska plovidba d.d.Zadar, in the shipyard and diesel engine factory SPLIT, a propulsion engine 6S70MC-C7 (Shipyard Split-Division Diesel Engines, Split, Croatia) was chosen with the following characteristics:  The optimising point was 85% SMCR and marked with a light blue point on the propeller curve, since normal continuous rating (NCR), which is also service rating (SCR), was 80% SMCR.
Lines L1-L3 and L2-L4 are isobars and referred to the respective mean effective pressure pme, and lines L1-L2 & L3-L4 are isotachs and related to the constant speed n, of the engine.The blue curve that passes through the point L1 is the heavy propeller curve, and it is expected that the voyage of the The optimising point was 85% SMCR and marked with a light blue point on the propeller curve, since normal continuous rating (NCR), which is also service rating (SCR), was 80% SMCR.
Lines L 1 -L 3 and L 2 -L 4 are isobars and referred to the respective mean effective pressure p me , and lines L 1 -L 2 & L 3 -L 4 are isotachs and related to the constant speed n, of the engine.The blue curve that Energies 2017, 10, 805 6 of 20 passes through the point L1 is the heavy propeller curve, and it is expected that the voyage of the ship will mainly reflect the actual conditions of navigation as hull fouling (covered by shells and algae), and environmental conditions (waves, wind, ocean currents, etc.).
The optimized power is 15,861 kW at 86.2 min −1 , and Service Power or NCR is taken to be 80% SMCR, which is 14,928 kW at 84.5 min −1 .

The Flow of Substances through the Engine
Data for International Organization for Standardization (ISO) ambient condition according to [27] are shown in Tables 1 and 2 for engine loads from 50% to 100%.Specific cylinder oil consumption b e,c.o. is 0.60 g/kWh at all engine loads.From Table 2 is obvious that scavenged air cooler heat had great potential, especially at higher engine loads.Scavenged air amounts, pressure and temperature were very important for the modelling of mass, substance and energy flow during breathing phase of the two stroke marine diesel engine.Scavenged air cooler heat, jacket water cooler heat and data from Table 1 were necessary for the calculation of waste energy which could be used in waste heat recovery (WHR) processes.
Basic data with mass, substances and energy flow through the engine has been shown on Figure 3 showing mass, energy, shaft work and heat flow through the control surface.In this open circuit heat engine, reactants (fuel and air) cross the control surface at inlet and products of combustion (exhaust gases) leave it at exit, while only work, but no heat crosses the control surface.Note, that the exhaust gases, though hot, convey energy but no heat.
. Q 1 and .Q 2 represent heat losses to the environment through cooling water.Cylinder oil .m c.o. was added at engine inlet since most of this oil burns in the cylinder.Air pressure after the turbocharger is 250 mm of the water column and higher than in the air distributor.oil burns in the cylinder.Air pressure after the turbocharger is 250 mm of the water column and higher than in the air distributor.

Calculation of Combustion Based on Fuel, Cylinder Oil and Air Consumption
Although the scavenged air consumption was taken from [27] at ISO standard conditions, the composition of the air was calculated by NIST [25] using our initial values.The fuel composition for MDO and cylinder oil composition were also taken into account.The heat input from the fuel was expressed in terms of its lower calorific value, while the fuel and air mass flow rates must meet the following relation: Specific mass flow through the engine: . .

Calculation of Combustion Based on Fuel, Cylinder Oil and Air Consumption
Although the scavenged air consumption was taken from [27] at ISO standard conditions, the composition of the air was calculated by NIST [25] using our initial values.The fuel composition for MDO and cylinder oil composition were also taken into account.The heat input from the fuel was expressed in terms of its lower calorific value, while the fuel and air mass flow rates must meet the following relation: Specific mass flow through the engine: where subscript a denotes air, f fuel and c.o. cylinder oil.

AIR
Oxygen content at the inlet of the cylinder: where N denotes the specific molar flow rate [kmol/kWh], M is molar mass [kg/kmol], and r is the molar and volume fraction.

FUEL
Fuel energy is: .
m f • ∆H T 0 ,p 0 kJ/s.The amount of carbon, hydrogen and sulfur from fuel at the inlet of the cylinder: .
where g is the mass fraction.

CYLINDER OIL
The amount of carbon and hydrogen from cylinder oil at the inlet of the cylinder: Sulfur content at the inlet was ignored as the size of a lower order.

Calculation of the amount of unburnt fuel and cylinder oil in exhaust gases
The mass of hydrocarbons entering the cylinder: The hydrocarbon (HC) share in the total mixture of the reactants is sm f HC /sm f react .

Emissions and unburnt fuel in the exhaust gases
According N + O → NO The stoichiometric amount of oxygen required for combustion:

Air Consumption Ratio/Factor
Now, the air consumption ratio/factor is: which is a ratio of oxygen entering the engine to oxygen necessary for stoichiometric combustion of fuel.The stoichiometric (dry air) air-fuel ratio is: The mole fraction of the fuel is difficult to calculate because the exact composition of the fuel is unknown.Therefore, the calculation is done by mass fraction.
Since the stoichiometric combustion of 0.1709 kg fuel/kWh required 0.01808288 kmol O 2 /kWh, this will be a set ratio with the amount of O 2 in dry air.
The AFR stoich for fuel combustion (without cylinder oil): The minimum amount of air required for combustion of fuel: . m a,stoich = .
And air consumption ratio can be defined as:

Calculation of Exhaust Gas Composition Using NIST Refprop09
The pressure drop through the exhaust system (spark arrestor, exhaust gas silencer, exhaust gas boiler and piping was given by Aabo [28].
The total quantity kmol/kWh in the exhaust gases: The share of components in the exhaust gases: where "i" is the index of chemical ingredients in the mixture of exhaust gases.

The Composition of the Flow in and out of a Two-Stroke Engine and the Cylinder
The absence of separate induction and exhaust strokes in two-stroke engines is fundamentally different from four-stroke engines.In the two-stroke engine, after a short blowdown exhaust at the end of expansion, the gas exchange process has induction and exhaust processes occurring simultaneously when the piston is near bottom dead centre (BDC).The most efficient scavenging process is the one whereby the products of combustion are completely replaced by a fresh charge, at the charge pressure and charge temperature.In fact, the fresh charge does not replace the products of combustion by a perfect displacement mechanism, and charging efficiency will be always smaller than unity.Another very important feature of the two-stroke exchange gas process is the unavoidable need for supercharging or turbocharging, which is not necessary in four-stroke engines.In Figure 4, the timing diagram for a two-stroke engine is shown.
process is the one whereby the products of combustion are completely replaced by a fresh charge, at the charge pressure and charge temperature.In fact, the fresh charge does not replace the products of combustion by a perfect displacement mechanism, and charging efficiency will be always smaller than unity.Another very important feature of the two-stroke exchange gas process is the unavoidable need for supercharging or turbocharging, which is not necessary in four-stroke engines.In Figure 4, the timing diagram for a two-stroke engine is shown.When defining the parameters of substance changes in the two-stroke engine with fuel injection, the following symbols will be used: p sc = scavenge air pressure at the entrance of scavenging ports T sc = scavenge air temperature at the entrance of scavenging ports T s = gas temperature in cylinder at the end of scavenge process .m sc = mass of air delivered to engine at pressure p sc and temp.T sc per sec m sc,h = mass of humid air delivered to the cylinder/cycle at pressure p sc and temperature T sc m sc,d = mass of dry air delivered to the cylinder/cycle at pressure p sc and temperature T sc m a,H 2 O = mass of water delivered with humid air to the cylinder/cycle at pressure p sc and temperature T sc m ar = mass of fresh air retained in the cylinder/cycle m ar,t = total air mass in the cylinder/cycle after scavenging port closed (SC) and exhaust valve closed (EC) m r = residual mass of gases in the cylinder/cycle m rp = residual (recirculated) mass of products in the cylinder/cycle m tr = m ar + m r = mass of trapped charge (fresh air and residual gas in the cylinder/cycle) after SC and EC m = mass of air that could be caught in the cylinder/cycle at a pressure p sc and a temperature T sc m a.d.stoich = mass of dry air required for stoichiometric combustion in the cylinder/cycle m a.h.stoich = mass of humid air required for stoichiometric combustion in the cylinder/cycle m a,excess = excess air in the cylinder/cycle m a.d.,s-c = mass of dry air, direct blowout from cylinder/cycle m f = the mass of fuel injected into the cylinder/cycle m co = the mass of cylinder oil injected into the cylinder/cycle m a.d.out = the total mass of the dry air at the outlet of the cylinder/cycle m p = mass products of combustion at the exit from the cylinder/cycle (including moisture from the air and emissions) m exh = the total mass of the exhaust gas at the outlet of the cylinder/cycle Figure 5 is very important and it enables proper definition of the parameters of working fluid exchange in the two stroke engines.Such a figure was first suggested by Sher [10] and gradually modified and improved.Figure 5 is significantly improved compared to the same figures presented in [8][9][10][11].In this figure, fuel was introduced as one of the important constituents which passes through the cylinder.Also, humid air intake, which takes part in combustion, is shown as dry air and water.
This compares the actual mass of air for scavenging the engine to the one required by the ideal turbocharging process.The reference mass is defined as the product of the cylinder maximum volume and the density of the ambient air or with a high supercharging engine as the product of the cylinder maximum volume and density of air in the air distributor.
This indicates to what extent the residual (recirculation) gases in the cylinder mr are replaced by a fresh charge mar.
The total mass of charge (trapped cylinder charge) is the sum of the mass of fresh air mar,t, the mass of burnt gas, and the mass of unburnt fuel from previous cycle mrp.This shows how much air comes directly to the exhaust.The parameters of working fluid exchange can be now defined according to Figure 5: (a) Delivery ratio (Scavenge ratio) Λ d : This compares the actual mass of air for scavenging the engine to the one required by the ideal turbocharging process.The reference mass is defined as the product of the cylinder maximum volume and the density of the ambient air or with a high supercharging engine as the product of the cylinder maximum volume and density of air in the air distributor.
(b) Charging efficiency η ch : This shows how successfully the cylinder volume is filled with fresh air.(c) Scavenging efficiency η sc : mass of delivered air retained mass of trraped cylinder charge = m ar m ar + m r = m ar m ar,t + m rp (28) This indicates to what extent the residual (recirculation) gases in the cylinder m r are replaced by a fresh charge m ar .
Energies 2017, 10, 805 13 of 20 The total mass of charge (trapped cylinder charge) is the sum of the mass of fresh air m ar,t , the mass of burnt gas, and the mass of unburnt fuel from previous cycle m rp .
(d) Retaining Efficiency η rt : mass of delivered air retained mass of delivered air / cycle = m ar m sc (29) This shows how much air comes directly to the exhaust.
The charging efficiency can be defined in terms of the delivery ratio and retaining efficiency: At the closing of the scavenging ports and exhaust valve, the retaining efficiency becomes the trapping efficiency.
The factor of relative charge is an indication of the degree of charging, and it is the ratio of charging efficiency and the scavenging efficiency: When the reference mass in the definition of the delivery ratio is the air and gas mass trapped in the cylinder m tr (or a close approximation of it), then: (g) Air-fuel ratio (AFR) for stoichiometric fuel combustion: The minimal quantity of air necessary for stoichiometric burning of injected fuel: . m a.d.stoich = .
m f • AFR stoich refers to the dry air.Now: (h) The excess air in the cylinder, and the relative ratio Λ excess is: Λ excess = 1/Φ where Φ is equivalence ratio.
(i) The air Consumption Factor/Ratio Λ AC : Λ AC = mass of air delivered/cycle mass of air necessary for stoichometric fuel combustion/cycle = m sc m a stoich In the steady state flow, Λ AC,dry = .
m a stoich or Λ AC,h = .

Results and Discussion
Results of our calculations are: % unburnt fuel = 0.9625% % unburnt fuel and cyl oil = sm f HC t unburnt b e + b c.o.
• 100 % unburnt fuel and cylinder oil = 1.2986%Thus, the total amount of unburnt hydrocarbons from fuel and cylinder oil is 1.3% over the adjusted values.
From the point of view of energy balance, emissions do not play a significant role.Depending on the design of the engine, the amount of unburnt fuel can reach more than 3%.However, due to environmental pollution, today much attention is paid to reducing emissions and reducing the CO 2 content in the exhaust gases.As of 14 March 2011, from MAN Diesel & Turbo data, modern 6SME-C diesel engines have emissions of only 0.3 g/kWh CO and 0.4 g/kWh HC.In the engine reported on here, total unburnt hydrocarbons from fuel and cylinder oil in SMCR was 2.227 g/kWh (HC, CO, soot), of which the fuel was only 1.645 g/kWh.Results of our calculations are presented in Tables 3-8.The minimum amount of air required for combustion of fuel: . m a,stoich = .
m f • AFR stoich = 12.94965028 kg air/s at100% SMCR Figure 6 shows AFR actual and stoichiometric air consumption change with main engine loads.Actual (total) AFR was very close to the results presented in the work of Guan et al. [29] for the main engine load from 50 to 100% SMCR.In this work, a new term was defined, the air consumption ratio/factor.This factor can be applied for two-stroke engines, especially for two-stroke fuel injection spark ignition engines, and for twostroke diesel engines.Since the engine runs humid air, A C ,hum id Λ is the definition which should be used, although AC,dry Λ is a more useful definition.Figure 7 shows changes of the air consumption ratio in the function of the main engine load for humid air and calculated dry air.At service power which is usually 80% of SMCR, air consumption was four times larger than necessary for stoichiometric combustion.In this work, a new term was defined, the air consumption ratio/factor.This factor can be applied for two-stroke engines, especially for two-stroke fuel injection spark ignition engines, and for two-stroke diesel engines.Since the engine runs humid air, Λ AC,humid is the definition which should be used, although Λ AC,dry is a more useful definition.Figure 7 shows changes of the air consumption ratio in the function of the main engine load for humid air and calculated dry air.At service power which is usually 80% of SMCR, air consumption was four times larger than necessary for stoichiometric combustion.In this work, a new term was defined, the air consumption ratio/factor.This factor can be applied for two-stroke engines, especially for two-stroke fuel injection spark ignition engines, and for twostroke diesel engines.Since the engine runs humid air, A C ,hum id Λ is the definition which should be used, although AC,dry Λ is a more useful definition.Figure 7 shows changes of the air consumption ratio in the function of the main engine load for humid air and calculated dry air.At service power which is usually 80% of SMCR, air consumption was four times larger than necessary for stoichiometric combustion.In Figure 8, cylinder oil and emissions data were from the MAN source [27].The inlet amount of air was from MAN data, but the composition was calculated using NIST program.The MDO data were data from the diesel factory-the shipyard SPLIT.Exhaust gas composition was calculated by our 0D program in combination with the NIST program.Although emission data was oriented, it was assumed that the figures relating to the emissions of NOx, CO, HC and particles when using HFO M exh.gas = Σr i • M i = 28.90631758kg/kmol R exh.gas = M exh.gas = 287.6305492J/kg K In Figure 8, cylinder oil and emissions data were from the MAN source [27].The inlet amount of air was from MAN data, but the composition was calculated using NIST program.The MDO data were data from the diesel factory-the shipyard SPLIT.Exhaust gas composition was calculated by our 0D program in combination with the NIST program.Although emission data was oriented, it was assumed that the figures relating to the emissions of NO x , CO, HC and particles when using HFO were very close to the actual values of the MDO fuel, for which our calculation was made.This assumes that all the cylinder lubricating oil entered and burned in the combustion chamber, and the small percentage that was scraped by the piston rings and remained below the piston was negligible.The size of the lower order NOx and SOx in the exhaust gases was ignored, while unburnt HC and soot (C) in the calculation was taken over hexane whose calorific value approximately corresponded to the calorific value of MDO.In order to simplify the exhaust gas calculation, SO2 was included in CO2.Results of our calculation for volume/molar composition are shown in Figure 9 for 100% SMCR engine load.Table 9 shows the exhaust gas composition over 50% to 100% of main engine load.The size of the lower order NO x and SO x in the exhaust gases was ignored, while unburnt HC and soot (C) in the calculation was taken over hexane whose calorific value approximately corresponded to the calorific value of MDO.In order to simplify the exhaust gas calculation, SO 2 was included in CO 2 .Results of our calculation for volume/molar composition are shown in Figure 9 for 100% SMCR engine load.Table 9 shows the exhaust gas composition over 50% to 100% of main engine load.Table 9 shows that there was no drastic change in the composition of the exhaust gas in the main engine load ranging from 50 to 100 SMCR.The biggest change was in the composition of CO2 and H2O as expected.The oxygen proportion in the exhaust gas was high and ranged from 15.73% (50% SMCR) to 14.78% (100% SMCR).The main reason was a short circuiting of air during the gas exchange process.The high consumption of air helps to clean exhaust gases in the cylinders, but also lowers  Table 9 shows that there was no drastic change in the composition of the exhaust gas in the main engine load ranging from 50 to 100 SMCR.The biggest change was in the composition of CO 2 and H 2 O as expected.The oxygen proportion in the exhaust gas was high and ranged from 15.73% (50% SMCR) to 14.78% (100% SMCR).The main reason was a short circuiting of air during the gas exchange process.The high consumption of air helps to clean exhaust gases in the cylinders, but also lowers the temperature of the exhaust gases, whose energy could be used in the WHR system.It also requires a higher capacity turbocharger.At the service power (80% SMCR), the engine consumed four times more air than is required for stoichiometric fuel combustion.
Thanks to Dr Fredrik H. Andersen from MAN we obtained the results of performance measurements, exhaust measurements and fuel specification and consumption conducted by MDT for an undisclosed 6S70MC-C7 engine.The data we obtained from Dr Andersen were for validating the model, and not for publication since data are confidential.Even though used MDO was not the same as the fuel used in our calculation, and the tested engine probably had different piston geometry, atomizer and engine lay-out, there were very good correlation between our results and factory test measurements performed.

Conclusions
This work developed a 0D model of combustion in a diesel engine.It was applied on a 6S70MC-C7 low speed engine, and by using the NIST program, thermodynamic parameters of air and fuel at the inlet of the engine and parameters of the exhaust gases were obtained.Our results were in very good correlation with MAN Performance measurements and exhaust measurements conducted on the same engine type.
The obtained data show that the composition of exhaust gases did not change significantly in the range of 50 to 100% SMCR.Also, it is observed that the quantity of oxygen in exhaust gases in the complete observed range was high at a value ca.15% (v/v).Dividing the total quantity of air which enters the engine with the air quantity required for the stoichiometric combustion, the air consumption ratio/factor was obtained and is shown in Figure 7; this parameter has not been previously used in scientific literature.In the main engine service rating, this factor value was at about four, which means that at this engine load it consumed four times more air than was required by stoichiometric combustion.The reason for this is that a great quantity of air flows directly from the air distributor to the exhaust during the exchange of working fluid.This achieves better blowout of cylinders and less residual exhaust gas in the cylinders, but decreases the exhaust gas temperature, therefore it lower efficiency of the WHR process.Further, for such a high quantity of fresh air, a turbo charger at higher capacity was required, which additionally decreases the energy of exhaust gases available for WHR process.Better swirling flow within the cylinder and bigger diameter of exhaust valve could probably decrease the air consumption ratio.
The new term Λ AC is defined as Λ AC,humid and Λ AC,dry .Since, humid air passes through the engine, Λ AC,humid should be used, although Λ AC,dry would be more useful.Although displacement volume V s is used as a reference volume in literature, we used V max for two stroke engines, which is much more appropriate.
Analysis of scavenge parameters of two stroke engines shows: Λ AC = Λ excess /η tr .This shows that air consumption ratio depended on the relative ratio and trapping efficiency.It is necessary to keep high the relative ratio to about 2, which means controlled/lower peak temperature, which is the main parameter governing NO x formation.Better trapping efficiency can be achieved by decreasing the mass of air in the direct blowout cylinder.By measuring and calculating the air consumption ratio and relative ratio, trapping efficiency can be calculated.
Figure 5 shows the flow of working substances in a cylinder and from the cylinder according to the [10,11].Figure 5 is a modified and significantly improved version of the figures shown in the literature [11], particularly in relation to [8][9][10].Two-stroke gas flow performance parameters can be easily defined using the improved diagram in Figure 5.

Figure 1 .
Figure 1.The propulsion engine model depicting main input and output substances and energy flows.

Figure 1 .
Figure 1.The propulsion engine model depicting main input and output substances and energy flows.
Bore = 700 mm Stroke = 2800 mm s/d=4.0c m = 8.5 m/s SMCR L 1 : n = 91 min −1 p me =19.0 bar NMCR = 18,660 kW Using MAN B&W CEAS_ERD (Engine Room Dimensioning) [27] with selection of the specified operating point i.e., Specified Maximum Continuous Rating (SMCR) is 100% Nominal Maximum Continuous Rating (NMCR) at 91 min −1 and marked by point L 1 on Figure 2. of reaction of combustion or the fuel caloric value.The suffix * means a quantity corresponding to one kilogram of fuel, and without suffix * is an amount per kilogram of products.

Figure 2 .
Figure 2. The quad of mean effective pressure and engine speed with chosen nominal continuous rating (NCR) [27].

Figure 2 .
Figure 2. The quad of mean effective pressure and engine speed with chosen nominal continuous rating (NCR) [27].

Figure 3 .
Figure 3. Schematic diagram of the conventional main engine 6S70MC-C7 for ISO conditions and 100% specified maximum continuous rating (SMCR).At the engine intake, a.h.169 600 kg/h = 47.1111kg/s m =  at 100% SMCR and ISO standard conditions (25 °C air temperature at the turbocharger inlet and fresh water cooling temperature, with ambient air pressure 1 bar and relative air humidity Φ=30%).p ws ws p 0.622 0.622 p p x p p p p Φ Φ ⋅ = = − ⋅ − (9) where x = vapor quality, pws = water saturation pressure [Pa], p=humid air pressure [Pa], and pp = partial vapor pressure [Pa].

Figure 3 .
Figure 3. Schematic diagram of the conventional main engine 6S70MC-C7 for ISO conditions and 100% specified maximum continuous rating (SMCR).
EC = Exhaust valve Closed EO = Exhaust valve Opening SC = Scavenging port Closed SO = Scavenging port Opening SP = Scavenging Period EBP = Exhaust Blowdown Period

Figure 4 .
Figure 4.The timing diagram for a two-stroke diesel engine.

Energies 2017, 10 , 805 12 of 19 Figure 5 .
Figure 5.The composition of the flow in and out of a two-stroke fuel injection engine and the cylinder after [10,11].The parameters of working fluid exchange can be now defined according to Figure 5: (a) Delivery ratio (Scavenge ratio) Λd: sc sc d 1 s c mass of delivered air reference mass successfully the cylinder volume is filled with fresh air.(c)Scavenging efficiency η sc: (d) Retaining Efficiency η rt: ar rt sc mass of delivered air retained mass of delivered air / cycle

Figure 5 .
Figure 5.The composition of the flow in and out of a two-stroke fuel injection engine and the cylinder after [10,11].
(e) Trapping efficiency η tr η tr = mass of trapped air charge mass of delivered air / cycle = m ar,t m sc (31) (f) Relative charge (Volumetric efficiency) Λ tr :

Figure 6 .
Figure 6.AFRactual and stoichiometric air consumption in the function of main engine load.

Figure 6 .
Figure 6.AFR actual and stoichiometric air consumption in the function of main engine load.

Figure 7 .
Figure 7. Air consumption ratio in the function of main engine load

Figure 7 .
Figure 7. Air consumption ratio in the function of main engine load.

Figure 8 .
Figure 8.The calculated substance flow through the two-stroke low speed engine at 100% SMCR.

Figure 8 .
Figure 8.The calculated substance flow through the two-stroke low speed engine at 100% SMCR.

Figure 9 .
Figure 9. Exhaust gas composition of the 6S70MC MAN B&W low speed two stroke engine.

Figure 9 .
Figure 9. Exhaust gas composition of the 6S70MC MAN B&W low speed two stroke engine.
removed by cooling water.The reference state of the environment is standard state; 25.0 • C and 1.00 bar.At steady state conditions the continuity equation is .
equal to fuel a higher caloric value.It follows that for these fuels,

Table 3 .
Fuel and cylinder oil consumption and mass composition at 100% SMCR, ISO conditions.

Table 4 .
Air consumption and mass composition at 100% SMCR, ISO conditions.

Table 5 .
Fuel composition and stoichiometric oxygen consumption for fuel and cylinder oil combustion at 100% SMCR, ISO standard conditions.

Table 6 .
Stoichiometric Air Consumption and Composition.

Table 7 .
The share of the components involved in the combustion.

Table 9 .
Exhaust gas composition at various main engine loads.

Table 9 .
Exhaust gas composition at various main engine loads.