Study on the Simplification Calculation Model of Marine Diesel Engine Exhaust Flow Based on Air-Fuel Ratio

Because of the large discharge of marine diesel engines, direct measurement is dicult and complicated. In order to obtain accurate data of marine diesel engine exhaust ow, the carbon balance method is generally used for calculation, but the carbon balance method has iterative operation, adopts wet concentration calculation, and does not consider sulfur oxide components, thus resulting in errors in the calculation results. In this paper, according to the calculationmodel of air-fuel ratio proposed by S H Chan, the law of conservation of C, H, and O atoms is observed, sulfur oxides in the exhaust are considered, and a new exhaust ow calculationmodel based on air-fuel ratio is obtained.is calculationmodel is the rst to directly use dry gas concentration in the calculation model of air-fuel ratio. en, the simplied calculation model is obtained by ignoring CO and simplifying coecient A∗, the dierence between the carbon balance method, ignoring CO, and the simplied model is compared, and then the exhaust ow value and exhaust density are calculated by the test data. e simplied calculation model can also calculate the mass ow rate of carbon dioxide, unburned hydrocarbon, and carbon monoxide. Finally, a software program is developed for the simplied calculation model, which is convenient for quick calculation and quick acquisition of exhaust ow data, thus laying a solid foundation for the application of the simplied calculation model in real ships.


Introduction
A lot of domestic and foreign scholars have studied the carbon balance method, and relevant references are listed in the following.
Chen et al. [1] constructed a carbon balance pressure index to measure the eco-environment pressure caused by carbon emissions in 77 countries from 2000 to 2015, and the logarithmic mean Divisia index decomposition method was used to identify the key factors related to carbon balance pressure.
e results of Nichiporovich [2] show that the average daily runo values are 0.409 mg C/m (2) in the pelagic zone and 0.647 mg C/m (2) in the littoral zone. e balance of carbon dioxide ux shows the universality of runo over emission.
Chang et al. [3] assume that, under the balanced and optimal economic growth path, the economic growth rate is equal to the consumption growth rate, from which we can obtain the ordinary di erential equation governing the consumption level by solving an optimal control problem.
In the paper by Zhang et al. [4], the hybrid graph is proposed to express the direct and indirect constraint relationship among components. e results of Viglizzo et al. [5] show that grazing lands generate C surpluses that not only could o set rural emissions but also could partially or totally o set the emissions of nonrural sectors. e results of Uri et al. [6]demonstrate that a wellregenerated young Scots pine stand on a former clear-cut area will be able to turn into a C sequestering ecosystem already before ten years after cutting. e results of Park et al. [7] show that CH 4 and carbon ow balance methods can be used to estimate model parameters appropriately and to predict long-term carbon emissions from land lls. e methodology of Crous-Duran et al. [8] is based on comparing the carbon emissions associated with the production of food and the carbon sequestered for that same activity for a particular quantity of food produced over a specific area and over a specific time.
Marin-Muniz and Hernandez [9] outlined a need to conserve and restore wetlands and demonstrated the wetland role as carbon sinks without concerning that they are sources of GHGs. e results of Daniel et al. [10] indicated that the hydromorphological alterations in these sites might convert healthy ecosystems contributing to C sequestration and climate change mitigation into C-emitting ecosystems.
e work of Morales et al. [11] shows that American Petroleum Institute results have proven to be an efficient tool for practitioners and researchers that intend to analyze the greenhouse gas emission of carbon capture utilization and storage systems to estimate, with accuracy, the global warming impact. e analysis of Zou et al. [12] shows that the NH 3 emission is highly correlated with the temperature and relative humidity. e ventilation rate shows a positive correlation with all the three gases.
In the paper of Liu et al. [13], their large-scale terrestrial carbon (C) estimating studies using methods such as atmospheric inversion, biogeochemical modeling, and field inventories have produced different results. e results of Doran-Browne et al. [14] showed that, from 2000 to 2014, Jigsaw Farms reduced its emissions by 48% by sequestering carbon in trees and soil. e study of Ventura et al. [15] examined the combined effect of woody ash and nitrogen (N) fertilization on the productivity, net C balance (NECB), and soil C sink capacity of a poplar SRC plantation established on a former arable land in northern Italy. e study by Rasanen et al. [16] explored carbon dioxide fluxes measured with the eddy covariance method for 3 years at a grazed savanna grassland in Welgegund, South Africa.
Weiler et al. [17] used the DayCent model to predict soil C dynamics for the next 30 years under green harvest management in scenarios with 0, 50, and 75% straw removal and to estimate the C savings of different soil tillage methods.
Many domestic and foreign scholars have studied the calculation of diesel engine exhaust flow; relevant references are listed in the following.
In the paper by Liu et al. [18], the calculation formula of the optimal area ratio (K-opt) of outlet and inlet of exhaust port was derived.
Ye et al. [19] numerically simulated the plume field distribution using an ideal jet model and a semiempirical formulation.
Arif et al. [20] proposed a methodology that will allow more accurate analysis of the effects of exhaust nozzle on the overall performance of a unique approach of analyzing jet exhaust nozzle integrated into aircraft and propulsion system. e experimental and numerical analyses of the standard fume hood features in order to determine the nature of the flow phenomena within the working chamber are presented and studied in the paper of Pietrowicz et al. [21].
In the paper by Pacak et al. [22], a theoretical analysis of heat and mass transfer in a counterflow heat exchanger used for energy recovery in air handling unit (AHU) under subzero operating conditions is presented. e results by Liu et al. [23] show that the far-field sound pressure level of the tailpipe predicted via the 3D numerical simulation method agrees well with the experimental results, indicating that accurate acoustic parameters of loads can effectively improve the source identification accuracy for an IC engine.
All of the above authors' studies used the carbon balance method to calculate, but the carbon balance method has iterative operation, adopts wet concentration calculation, and does not consider sulfur oxide components, thus resulting in errors in the calculation results. e carbon balance method of diesel engine exhaust flow calculation is very complicated, and the air-fuel ratio model is not used in the calculation of diesel engine exhaust flow. Due to the above shortcomings of the carbon balance method and the existing diesel engine exhaust flow calculation, in order to accurately and quickly calculate the marine diesel engine exhaust flow, this paper conducts further research on the airfuel ratio calculation model and finally obtains a set of simplified calculation model.
In this paper, based on the calculation model of air-fuel ratio proposed by S H Chan, a new calculation model of airfuel ratio is obtained by considering the sulfur oxide composition in exhaust gas and a simplified calculation model of exhaust mass flow rate is established.

Modified Calculation Model of the Exhaust Flow Based on Air-Fuel Ratio
We know that the exhaust volume of a marine diesel engine is related to the volume change of the cylinder intake and the volume change before and after fuel spray. In the ideal gas condition, the volume change coefficient between the intake and exhaust volume is the coefficient of molecular weight change of the work medium before and after combustion. According to the coefficient of molecular weight change of the work medium before and after combustion, a modified calculation model for the exhaust flow is established. e specific research ideas are as follows. According to the characteristics and exhaust components of diesel engines, this paper gets the total exhaust flow with the ratio of air to fuel and then gets the flow of each component according to the proportion of each component.
where g c is the amount of carbon contained in a unit mass of fuel; g H is the amount of hydrogen contained in a unit mass of fuel; ϕ is the percentage of carbon expelled as free carbon; φ is the mass fraction of CO in the product; ψ a is the excess air coefficient; and L 0 is the theoretical air-fuel ratio, which is the amount of air required for complete combustion of a unit mass of fuel calculated according to the chemical equation. It is the minimum amount of air required for the complete combustion of 1 kg of fuel. According to the reaction equation of C, H, and O, considering that 1 kg of fuel already contains oxygen (32 kg/mol), as well as the approximate volume percentages of oxygen and nitrogen in the air of 0.21 and 0.79, respectively, the amount of combustion products is obtained as follows: (2) Substituting the expression of the theoretical air-fuel ratio L 0 into the above equation, e volume of the fuel droplets injected into the cylinder is small compared to the volume of the inlet air and is omitted.
e work medium before combustion meets L � ψ a L 0 . In the exhaust composition of low-speed twostroke marine diesel engine, the CO concentration accounts for 2.56% of the harmful gas emissions, and the HC concentration accounts for 7.70%. erefore, the coefficient of molecular weight change of the work medium is δ.

Determining the Air-Fuel Ratio L.
e main components of marine fuel are carbon and hydrogen (about 95% or so), plus small amounts of oxygen, nitrogen, sulfur, and metal compounds. erefore, the molecular formula of fuel oil is represented by "C α H β O c N ξ S ζ ." It is assumed that the main components of air are 78.09% of nitrogen, 20.94% of oxygen, 0.93% of argon, and 0.04% of carbon dioxide. [24]. e equation for combustion (considering air humidity and incomplete combustion and the resulting B) is as follows: From Dalton's law of partial pressures, it is known that, in an ideal gas mixture, the ratio of the moles of the components is equal to the ratio of the partial pressures [25]. e relationship between air humidity and A * is given by the following equation: p is the total pressure of wet air and p v is the water vapor partial pressure in wet air. A * , B, and A are coefficients.

Mathematical Problems in Engineering
Φ is the ratio of p v to p s at the same temperature under the same total pressure state. p s is the water vapor partial pressure in saturated wet air. Φ is the relative humidity. e following equation is obtained: Equations (7) and (8) can be organized as Air-fuel ratio L is expressed as follows: In formula (10), α is the molecular number of C, β is the molecular number of H, c is the molecular number of O, ξ is the molecular number of N, and ζ is the molecular number of S.
To find the air-fuel ratio L, we need to first find A * . e reaction equation of fuel combustion, according to the balance of C, H, O, and N atoms, can be obtained from the following relationship: Normally, a particular instrument is not used to measure water vapor concentration alone. If the water vapor condenses on the instrument, it will cause the instrument to not work properly. e usual method used is to remove the water by condensing the sample gas before it reaches the instrument. e resulting gas concentration is thus the dry gas concentration without water vapor [26]. For a given gas, the dry and wet concentration is related as follows: e dry gas concentration is thus expressed as follows: erefore, equations (11) to (14) can be expressed as follows: 2B + β � 2y 4 + xy 0 5 + 2y 0 where y 0 i is the instrumental measurement and A * , B, y 4 , y 6 , and 11 i�1 v i are unknown. To facilitate the solution, new conditions need to be introduced. e reaction equation for the conversion from water gas is expressed as follows: CO and CO 2 concentrations can be related by the equilibrium constant K.
K is a function of the equilibrium temperature T eq and has values between 3.5 and 3.8 [27]. Usually, the equilibrium temperature is between 400 K and 3200 K, with the following expression: Solving equations (15) and (25) to (18) and (20) together yields Among them,

Mathematical Problems in Engineering
Φp s p − Φp s .
(30) e actual air-fuel ratio can be obtained after formula (22) is brought into formula (10), and the calculation formula for the remaining gases can be obtained. e details are expressed as follows:

A Correction Factor KW from Dry Basis to Wet
Basis Concentration. When calculating the mass flow or volume flow of wet exhaust, if the measured concentration is dry concentration (conc D ), it should be converted to wet concentration (conc W ) according to the following equation: where k w is the dry (wet) concentration correction factor. k w can be calculated. For exhaust, the value of k w can be determined according to the following equation: H a � 6.220 · R a · P a P B − P a · R a · 10 −2 .
(36) e symbol expression is shown in Table 1. For the exhaust gas, K w,r can also be determined by the following equation: For the intake air, the value of K w,a can be determined according to the following equation: where K w,a is the dry (wet) correction factor of the inhaled air.

Determination of the Fuel Consumption per Unit n Mile of the Ship's Main Engine.
ere is a calculated relationship between the fuel consumption of the main engine per hour as g f (kg/h), the effective power of the main engine as P e (kW), and the fuel consumption of the unit effective power output of the main engine per unit time as g fP (kg/kWh) as expressed in the following: Assuming the ship speed V s (kn, n mile/h), the fuel consumption per nautical mile of the main engine range is expressed as follows:

Determining the Formula for Calculating the Exhaust Flow
Rate. If the fuel consumption G f (kg/n mile) per unit mile of the ship's main engine is known, the exhaust composition of the ship's main engine and the concentration of each component can be obtained by gas detection instruments so that the mass flow rate G M (kg/n mile) and volume flow rate of the exhaust can be calculated. e calculation equations are expressed as follows: e expression for the volumetric flow rate of a gas is expressed as follows:

Mathematical Problems in Engineering 5
Another expression for the volumetric flow rate is expressed as follows: In the above equation, where G f is the mass flow rate of fuel consumption; c C , c S , c H , and c N are the mass fractions accounted for by carbon, sulfur, hydrogen, and nitrogen in the fuel, respectively; W denotes the relative atomic mass or relative molecular mass; M represents the molar volume; and ϕ CO , ϕ HC , ϕ O 2 , ϕ NO 2 , and ϕ NO represent the dry volume concentrations of CO, HC, O 2 , NO 2 , and NO in the exhaust, respectively.

Calculation of the Mass Flow Rate of Each Component Gas.
Based on the measured fuel flow and the concentration values of NO X , HC, CO, CO 2 , and O 2 , the mass flow rates of each component gas are calculated. e mass flow rates of the main gas components in the exhaust are calculated as indicated in the following: ϕ HC and ϕ NO x are the measured wet basis concentrations expressed in ppm concentrations, or if the NO X concentration is expressed in dry basis, add K w to (28). ϕ CO is the measured dry basis concentration expressed in ppm.
K HDIES is the correction factor for NO X .
It should be noted that, in the calculation of HC, the C/H ratio is taken as 1 : 1.85. In addition, the coefficients in the NO X equation are obtained by calculating N 2 O, which is the final form of NO X in the atmosphere.
Based on the exhaust mass flow rate obtained from the previous calculation, the mass flow rate of each component gas can be calculated from the following equation, using CO 2 as an example: where G CO 2 is the exhaust mass flow rate of CO 2 and μ takes the value 0.1519 when the exhaust density ρ is 1.293 kg/m 3 . If the exhaust density ρ is not 1.293 kg/m 3 , then the following relationship exists: e parameters μ and ω are taken as shown in Table 2.

Test Data Calculation
e test was conducted with a Shanghai Dongfeng G128ZCo high-speed diesel engine, which has a bore of 135 mm, a stroke of 150 mm, a rated speed of 1500 rpm, and a rated power of 162 kW. Table 3 shows the specified experimental modes and weighting factors. e marine diesel engine was tested in accordance with the E3 mode. Table 4 shows the fuel composition and its content in the test. Table 5 shows the test measurement data.

Calculation of the Air-Fuel Ratio L.
According to the composition of the fuel, it can be known that α: β: c: ξ: ζ is equal to 1 : 1.925 : 0.00087 : 0.001 : 0.0007. X/Z � 1.85, K is taken as 3.5, α � 12 as the base, β � 23.1, c � 0.01, and the relative molecular mass of the fuel is counted as 168.01. Bringing the data into formulas (10) and (25)-(30), respectively, we can find A 1 , A 2 , A 3 , A 4 , A 5 , A, and the value of L, as shown in Table 6.

Calculation of the Fuel Consumption G f per Unit n Mile for
Diesel Engines. Assuming a speed V s (kn) of 10, according to equation (40), G f can be obtained as shown in Table 7.

Symbol
Expression H a Absolute humidity, that is, the mass of water contained in each kg of dry inlet gas (g) R a Relative humidity of the incoming gas (%) P a Saturated vapor pressure of the incoming gas (kPa) P B Total pressure of the incoming gas (kPa) K w,r e correction factor for the dry (wet) concentration of the original exhaust K w2 e volume share of water vapor in the wet inlet gas

Determine the Mass Flow Rate of Exhaust Gas (kg/n Mile).
According to (41) and (45), the mass ow rate values of exhaust gas are shown in Table 8.

Determine the Volume Flow Rate of the Exhaust Gas (m 3 /n Mile).
e gas ow rate υ (m/s) is measured by the gas analyzer, the cross-sectional area of the pipe S (m 3 ) can be calculated to obtain the gas ow rate V (m 3 /s), the ship's speed under stable working condition is V s (n mile/h), and the volumetric ow rate of the gas can be obtained as expressed as follows: e calculated data can be obtained as shown in Table 9. It is found that the average value of the measured exhaust density is 1.2875 kg/m 3 , which is not equal to the air density of 1.293 kg/m 3 , and the deviation was within 1%, which is within the acceptable range.

Simplification of a Calculation Method
e modi ed method of calculation formula discussed above is complex; for the convenience of calculation and application, it is necessary to simplify the calculation formula.
Due to the low value of CO concentration in the diesel engine, if you do not consider the exhaust CO, an analysis of the formula can nd that the product of A 3 A 4 terms can exactly eliminate the terms with CO concentration as the multiplier; namely, Notated as A', the data are recalculated and obtained as shown in Table 10.
It can be found that, compared with the carbon balance method, the e ect of CO concentration on the calculated value of exhaust mass ow is reduced by 1.2% to 1.4%. erefore, CO can be ignored for the calculation. e comparison of the exhaust mass ow rate calculated by ignoring CO is shown in Figure 1.
Moreover, in the process of calculation, it is found that the denominator part of the A * expression is calculated very close to 2, and the numerator part of the c term has a small value, so the expression of A * can be simpli ed as shown as follows: e calculated air-fuel ratio data are shown in Table 11.
A comparison of the exhaust mass ow rate calculated by the simpli ed model is shown in Figure 2.
It can be found that the calculated values of this simpli ed model are small compared with the carbon balance method, and the deviations are 2.8%, 4.3%, 4.9%, and 5.6%, respectively. e deviations of the simpli ed model are all within 10%, which is acceptable due to the requirement of 10% deviation in the NO X emission speci cation issued by IMO.
rough experiments under di erent measurement conditions, it is found that when the simpli ed calculation model needs more gas components, due to the di erent measurement accuracy of various instruments, some gas     measurement is not accurate, which inevitably leads to a little error in the simpli ed calculation model. In addition, because the exhaust gas of the diesel engine contains water vapor, the gas concentration measured by the gas analyzer is usually obtained after the removal of water vapor, so there will be some deviation from the real value in the calculation.

Development of a Software Program for the Simplification Calculation Method
By analyzing, modifying, and simplifying the calculation model of the exhaust mass ow rate of marine diesel engine and programming it with the computer language Visual Basic, the calculation speed can be greatly improved. e   G 4.9 * (X + L)/(X + Y + N) + 9.8 * (Y + W + M)/(X + Y + N) + 5.7-5.
Label9.Caption "" and str(G) 12 Label11.Caption "" and str(K) software program for this calculation model is shown in Table 12. e interface of the designed program and some calculation results are shown in Figures 3 and 4. By running this program, the values of air-fuel ratio and exhaust mass flow rate can be obtained quickly, which is convenient and practical.

Conclusions
e model for calculating the exhaust mass flow of marine diesel engines obtained in this paper has the following advantages compared with the carbon balance calculation methods in MARPOL Convention Annex VI.

Comprehensive Consideration and More Accurate
Calculation.
e calculation model takes into account the sulfur, nitrogen, and oxygen components in the fuel and is applicable to the calculation of common fuel, while the carbon balance method is only applicable to the fuel containing only hydrogen and carbon components. e correction of the two calculation methods in the MARPOL Convention Annex VI caused by changes in the humidity of the intake air (when the absolute humidity of the intake air is greater than or equal to the humidity of the charge air, i.e., when Ha ≥ Hsc, part of the water vapor condensates in the intake air) should be made after the exhaust mass flow rate is derived. e modified method corrects the parameters for the wet intake air before the exhaust flow calculation so that the calculated exhaust mass flow is the corrected mass flow.
e results are more accurate since the correction is not required to account for the amount of fuel consumed.

Avoiding the Iterative Operations.
e carbon balance method requires iterations of the exhaust density or the wet and dry correction factors for the three cases of each exhaust  component concentration: all wet, all dry, or wet and dry, which may lead to problems of nonconvergence of the iterations. e modified calculation model does not need to consider the iterative calculation when the exhaust concentration is all dry; when the measured concentration is both wet and dry, the problem of iterative nonconvergence will not occur through a one-step correction of the concentration.
Although the simplified calculation method in this paper is based on the carbon balance method, it has its own uniqueness. To avoid iterative calculation, dry concentration is directly used in the calculation formula, which is convenient and fast.

Data Availability
All data, models, and code generated or used during the study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.