Next Article in Journal
Surface Morphology of Three-Dimensionally Printed Replicas of Upper Dental Arches
Next Article in Special Issue
Benchmark Sea Trials on a 6-Meter Boat Powered by Kite
Previous Article in Journal
Antihypertensive Effect of Amaranth Hydrolysate Is Comparable to the Effect of Low-Intensity Physical Activity
Previous Article in Special Issue
Method to Increase the Accuracy of Large Crankshaft Geometry Measurements Using Counterweights to Minimize Elastic Deformations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance In-Live of Marine Engines: A Tool for Its Evaluation

Department of Engineering and Architecture, University of Trieste, Via A.Valerio 10, 34127 Trieste, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(16), 5707; https://doi.org/10.3390/app10165707
Submission received: 9 July 2020 / Revised: 7 August 2020 / Accepted: 13 August 2020 / Published: 17 August 2020
(This article belongs to the Special Issue Ship Energy Systems)

Abstract

:
Currently, most ships use internal combustion engines (ICEs) either as propulsion engines or generator sets. The growing concern in environmental protection along with the consequent international rule framework motivated shipowners and designers to replace conventional power systems in order to mitigate pollutant emissions. Therefore, manufacturers have made available on the market many technological solutions to use alternative fuels (Liquefied Natural Gas or LNG, methanol, etc.). However, the main energy source is still fossil fuel, so almost all the ICEs are made up of turbocharged diesel engines (TDEs). TDEs have still the potential to improve their efficiency and reduce fuel consumption and pollutant emissions. In particular, the interpretation of Industry 4.0 given by manufacturers enabled the installation of a robust network of sensors on TDEs, which is able to allow reliable power management systems and make ships much more efficient regarding operating costs (fuel consumption and maintenance) and environmental footprint. In this paper, a software tool that is capable of processing the in-live performance of TDEs is described. The great novelty consists in the ability to process all the information detected by the sensor network in-live and dynamically optimize TDEs’ operation, whereas the common practice involves the collection of performance data and their off-line processing.

1. Introduction

On board all operative ships, internal combustion engines (ICEs) are installed with the function of both propulsion engines and auxiliaries in electric energy generators. Specifically in large ships, ICEs are mainly 2-stroke and 4-stroke turbocharged diesel engines. Two-stroke engines are used exclusively for the propulsion of merchant ships such as very large container ships, tankers, and ro-ro due to the required great power. Moreover, the layout of the general arrangements of these types of ships does not create particular constraints on the installation of these mammoth-sized engines. On the other hand, 4-stroke engines, coupled with suitable generators, are used for the production of electric energy practically on board all ships. Moreover, thanks to the high specific power, they are mainly employed in the propulsion of ferries, naval ships (up to frigates), pleasure crafts, fishing vessels, and support vessels for offshore operations. It is still worth underlining that 4-stroke engines are also used in passenger ship gensets, which, as known, adopt an electric propulsion system.
For both types of engines, the power they can deliver ranges from 10 to 100 MW; therefore, they must be fed with low-cost fuels that make the shipping of goods and passengers profitable for shipowners. In fact, as reported in [1], global shipping is responsible for moving about 90% of the world trade, while cruise ships in 2017 allowed nearly 25 million of passengers to sail the oceans. These significant results were obtained thanks to the great flexibility offered by shipping and the low cost per ton transported of this mode of transportation. Despite maritime transport still being the most fuel-efficient means of good transit, it should be noted that ships have, as mentioned before, big ICEs installed on board; as a consequence, shipping greatly contributes to the emission of air pollutants (PM, SOX, NOX), sea pollutants (noise, thermal, discharges), and greenhouse gases (GHG).
To be honest, in order to reduce all forms of pollution from ships, since 1973, the International Maritime Organization (IMO) has already issued a whole series of regulations. Moreover, the IMO has drawn up a roadmap for reducing emissions from ships in order to provide shipowners and shipbuilders with a clear long-term strategy to guide their research and development activities. With particular reference to air emissions, in 2006, MARPOL Annex 6 [2] was issued introducing the definition of Emission Control Areas, which are the areas of the planet with precise limits of air pollutant from ships allowed. Given the size of the world fleet and the partially retroactive effect of these regulations, several steps and implementation dates have been identified. On the same occasion, the compulsory adoption of the Ship Energy-Efficiency Management Plan (SEEMP) was introduced in order to limit waste as much as possible. Subsequently, IMO also ruled in favor of reducing GHG emissions [3], setting the 50% reduction of CO2 produced by ships as a target in 2050 compared to that produced in 2011 (year of entry into force of [2]).
On this premise, from 2006 onwards, countless studies have been carried out to evaluate reliable, economic, and efficient technological solutions capable of achieving these objectives. In particular, three solutions have been identified [4]: reducing the hull resistance; increasing the overall efficiency of the propulsion and power generation systems; and using alternative fuels with low or zero pollutant contents. It is easy to understand that the first strategy is really effective if applied on the construction of new ships, while the other two ones have found multiple applications also in the retrofitting of existing ships. In particular, the use of Liquefied Natural Gas (LNG) as a fuel is now a consolidated reality; additionally, many studies are being carried out in order to evaluate the use of other fuels such as hydrogen, methanol, and ammonia [5,6,7]. However, given the effectiveness of the application of exhaust gas treatment systems such as SCR (Selective Catalytic Reduction) and Scrubber, the main energy source in marine industry is still Heavy Fuel Oil (HFO). In reality, this choice is also a consequence of the fact that shipowners must guarantee an adequate return on investment, and the major contributors in the life cycle cost of a ship are precisely maintenance and fuel. Therefore, HFO is preferred because besides having the lowest cost among fuels, it is very available on the market, has a low sulfur content variant, and has a consolidated and branched distribution network. However, it should also be noted that the standards currently achieved by the main marine diesel engine manufacturers ensure that the conventional marine turbocharged diesel engines (TDEs) have the potential to further reduce their fuel consumption and their emissions, improving their energy efficiency.
Therefore, monitoring and setting the performance of TDEs is vital: indeed, the efficiency of any onboard machinery is directly related to its performance in order to get an efficient combustion. Shipowners know very well that ensuring an efficient combustion reduces not only pollutant emissions but also fuel and maintenance costs. Conventionally, there are several ways to monitor the performance of a TDE: measuring the peak pressure by mechanical peak pressure gauge; indicator card measurement; digital pressure monitoring; intelligent combustion monitoring; monitoring of engine control parameters; engine parameters; log book monitoring; and engine emission evaluation. All these techniques, which are actually very effective, currently have the limit of providing for the offline processing of the data obtained in the monitoring and subsequent adjustment of the engine. However, numerous studies [8,9,10,11,12] have demonstrated the possibility of simulating engine operation with remarkable precision.
The adoption of the Industry 4.0 paradigm by engine manufacturers has led the development of a new generation of products equipped with advanced automation systems [13]. These automation systems are made up of a dense network of sensors that allows the detection of performance data of the TDEs and a series of actuators that are able to adjust the functional parameters of most of the auxiliary devices, in order to make the TDE always work at a maximum of its possibilities. Unfortunately, the common current practice foresees that, especially when new TDEs are installed, the manufacturer controls the monitoring functions, and the data collected are exploited only for diagnostic and maintenance purposes.
In this paper, a software tool that is capable of processing the in-live performance of TDEs is described. In particular, with reference to the main international standards, the tool instantly processes the data collected by the network of sensors installed on board, obtaining the engine performance data in-live and allowing the dynamic optimization of TDE operation. Particular attention was paid to writing the routine for the evaluation of air polluting emissions, because the latter can be significantly reduced in the first instance by acting on the combustion parameters and subsequently on the actuators that regulate the operation of both the SCR and Scrubber. Given the importance of determining the emissions, the tool calculates them by using two methods, the Air Nozzle and the Carbon Balance methods, both widely described below.

2. Rule Framework

Regarding the software’s mathematical steps required for the in-live performance evaluation, the following regulations have been analyzed:
  • Regulation ISO 5167-1:2003 [14]. In this document, the Technical Committee establishes the general principles for methods of measurement and computation of the flowrate of fluid flowing in a conduit by means of pressure differential devices. Thus, Regulation ISO 5167-3:2003 [15] provides the procedure for the evaluation of the incoming air mass flowrate to the compressor m ˙ a i r by means of an Air Nozzle. Via the calculation of an ideal air mass flow through an Air Nozzle with Bernoulli’s equation and the implementation of discharge coefficient C D and expansion coefficient ϵ, the real air mass flow can be rated;
  • Resolution MEPC.177 (58) (NOx Technical Code 2008) [16]. This technical code, besides providing mandatory procedures for the testing, survey, and certification of marine diesel engines, gives a method to estimate exhaust mass flowrate m ˙ e x h by means of Carbon Balance. The method is demonstrated in ISO 8178-1 [17] through the equivalence between carbon entering the engine versus carbon leaving the engine;
  • Regulation ISO 8217:2010 [18]. Within this regulation, standard values for fuels (in terms of density, viscosity, flash point, and much more) are provided. These are used if correct values from tests on in-use fuels are not available;
  • CIMAC Number 27|2007 [19]. The overall efficiency of a turbocharger η T l can be evaluated by using the general formula provided in this document. Anyway, some manufactures (ABB and Mitsubishi) allow using simplified formulas, derived from the general one, which are reported in manufacturer’s technical notes;
  • Resolution MEPC.176(58)|2008 [20]. For air pollution emission in atmosphere (in terms principally of nitrogen oxides NOX and sulfur oxides SOX emissions), the worldwide thresholds are fixed in this International Convention. This resolution embeds the previous 1997 Protocol that sets limits on sulfur oxide and nitrogen oxide emissions from ship exhausts and prohibits deliberate emissions of ozone-depleting substances. As regards the European Union, the Directive EU 2016/802 [21], commonly known as “EU Sulphur Directive”, concerning the maximum amount of sulfur dioxide emissions in the atmosphere, and the Regulation C(2016)8381 [22], concerning the general and technical requirements for emission limits and type approval for non-road ICEs, are issued.

3. Tool Description

3.1. General Principles

The software described in this paper was developed using as environment LABVIEW (Laboratory of Virtual Instrumentation and Engineering Workbench), which is system engineering software for applications that requires tests, measurements, and control. The programming language is called Graphic Language (G-Language) and differs from the traditional ones because the syntax is not written but graphic. This language is able to perform multithreading spontaneously without explicit control from the programmer; consequently, it is an excellent tool for handling all the processes in parallel. In order to define a digital configuration environment that is capable of perfectly adapting to the real on-board configuration, the number of independent choices to perform is very high. As a result, the software requires several processes.
The calculation procedure of the performance values starts with the evaluation of the incoming air mass flowrate to the engine m ˙ a i r   e n g .
This quantity, expressed in kg / s , is evaluated through two different procedures:
  • by means of the evaluation of incoming air mass flow to the compressor m ˙ a i r , which is measured if the Air Nozzle is installed;
  • by means of the evaluation of exhaust gas mass flow m ˙ e x h , calculated using the Carbon Balance method if CO2, CO and THC (Total HydroCarbon) detectors are present.
As regards the logical functioning of the software, four consequential phases can be identified (Figure 1).
During the first phase of Configuration Set Up, the software requires the user to make all the choices previously mentioned in order to generate a Set Up Environment. The latter will be used to communicate to the actual engine of the software, which is called Main, the calculations to perform and the constants and variables necessary for them; furthermore, it provides an association between the variables and the sensors detected in the network.
The next phase requires running the parsing—that is, the association between the variables used in the Calculation Routine and the numerical values (i.e., the physical sensors) detected in the Sensor Network (SN). The user has to perform this phase manually every time.
Once the parsing is completed, the software starts with the Data Processing phase. The Main determines the calculations to perform and the input to use, which is based on available in-use configuration information, manually inserted constant values, and variables collected by the SN. Indeed, the analysis of both classic bibliography and the Rules in Performance Calculation area allowed identifying the functions required to define the interest value in Performance Evaluation. A huge number of them can be calculated only if the Border Conditions provided by the Set Up Environment are fulfilled; otherwise, the software does not perform the calculation of those functions and proceeds to verify all the others.
The last phase concerns the real calculation and the in-live visualization of the performance through a continuous update of the input values detected by the SN.

3.2. Configuration Set Up

As regards to Configuration Set Up phase (Figure 2), three principal Routines can be defined. Each of them refers to an .xml database file that is able to contain all the choices made in that specific Routine by the user during the Input Configuration phase (if saved). These databases will be used to speed up the future operations of the Configuration Set Up. The Routines are the following ones:
  • Routine Test Bed: routine used to insert all the information concerning the Test Bed instrumentation installed; it is composed of five SubRoutines: SubRoutine Air Nozzle, SubRoutine Engine Torque System, SubRoutine Fuel Analysis and Measurement, SubRoutine Emissions, and SubRoutine Humidity.
    Furthermore, the SubRoutine Fuel Analysis and Measurement refers to its own .xml database file, containing the chemical characteristics of the fuel used by the engine. In this way, the opportunity of monitoring the variation in terms of carbon, hydrogen, nitrogen, oxygen, and sulfur concentration of a fuel compared to another one is provided.
  • Routine Engine Characteristic: routine used to insert all the information concerning the engine;
  • Routine Turbocharger: routine used to insert all the information concerning the turbocharger.
With these Routines, the user must make a series of choices necessary to define a Set Up Environment that is able to identify uniquely both the sensors and instrumentations list installed and the Border Conditions required for the determination of computability of the following mathematical functions.
Among the choices, the user has also to answer about the presence on the SN of every sensor. If the sensor is not online, then the software requires manually inserting a constant value; otherwise, the variable name becomes available in the parsing phase for the matching with the relevant sensor.

3.3. Output

The Routines allow the calculation of several outputs: in particular, the following quantities, summarized in Table 1, represent the most important parameters for Performance Evaluation.
Turbocharger efficiencies in terms of compressor efficiency η S V , turbine efficiency η T , and overall efficiency η T l can be determined by both formulations contained in the literature (e.g., CIMAC [19]) and formulations recommended by manufacturer (e.g. ABB and Mitsubishi).
The compressor pressure ratio (total to total) Π V , turbine pressure ratio (total to partial) Π T , and corrected turbocharger speed n 298 represent the inputs for the most common in-use compressor and turbine maps. Indeed, the implementation of in-live visualization of the operating point of the compressor and the turbine inside the relative map is still under study.
The Brake Specific Fuel Consumption (evaluated for each contribution of fuel, corrected for L H V I S O (ISO Lower Heating Value) of the fuel in use and total) and the Effective Compression Ratio are necessary to evaluate the engine consumption and its thermal efficiency.
Eventually, the main parameters for Performance Evaluation are NOX, CO, CO2, THC, O2, and SO2 concentrations, which are expressed in g / kWh . This unit of measurement is very important because is the one used to verify the compliance with the standards regarding air pollution established by MARPOL 2017 [2].
Once these outputs have been calculated, to optimize the TDE’s functioning (in terms of both increasing turbocharger efficiency and decreasing fuel consumption and emissions, given in g/kWh), the engine automation, in the forthcoming future, will be directly connected with the tool by means of a dedicated Programmable Logic Controller (PLC). Thus, through the automatic modification of some electromechanical processes (such as wastegate opening and inlet and exhaust valve closure), the self-acting improving of Performance values will be achieved.

4. Tool Application to Evaluate Pollutant Emissions

For marine engines, pollution production is strictly related to their working. Indeed, pollutant substances such as NOX, CO, CO2, THC, O2, and SO2 are a natural consequence of the combustion process. However, their production can be mitigated through more efficient and high-performing engines and combustion process. In this regard, the use of a lower quantity of fuel, without affecting the output power, is one of the first measures to apply. Furthermore, also a cleaner combustion process ensures a significant reduction of pollutant emissions. As a result, in order to achieve the best performance from a TDE, it’s fundamental to keep a continuous monitoring and control of the engine while running.
As regards the paper’s subject, the evaluation of emissions will be thoroughly examined hereafter in terms of both the formulation’s explanation and Emission Routine’s logic. As a reference for the thermodynamic quantities, the WIT (Wärtsilä Italia) Convention (in Appendix A) is used.
To define the polluting values in terms of g / kWh , the exhaust mass flowrate must be calculated. As mentioned above, two methods have been identified:
  • By means of the Air Nozzle approach, from which the incoming air mass flowrate to the compressor can be measured, and then the exhaust mass flowrate can be calculated by simple mathematical steps; and
  • By means of the Carbon Balance method, from which the exhaust mass flowrate can be directly evaluated.

4.1. Air Nozzle Approach

The Air Nozzle is a device consisting of a convergent inlet connected to a cylindrical section installed at the compressor inlet. The reference Regulation regarding the Air Nozzle approach is the above-mentioned ISO 5167-3:2003 [15].
In this approach, the formula for the air flowrate calculation is the same for both two-stroke and four-stroke engines. It derives from the Bernoulli’s equation for ideal fluids passing through a section shrinkage (Figure 3). Indeed, for incompressible fluids, the following equation is valid:
p 1   +   1 2 ρ V 1 2   =   p 2   +   1 2 ρ V 2 2
in which subscripts “1” and “2” mark the fluid situation before/upstream and after/downstream the shrinkage, respectively. In both cases, the fluid situation must be measured at a point at which the flow can be considered stationary.
The continuity equation gives the following equation:
V ˙ i d e a l   =   A 1 V 1     = A 2 V 2 V 1   =   V ˙ i d e a l A 1   and   V 2     =   V ˙ i d e a l A 2
where:
  • V ˙ i d e a l   [ m 3 / s ] is the volumetric flowrate;
  • A 1   and   A 2   [ m 2 ] are Section 1 and Section 2 areas, respectively;
  • V 1   and   V 2   [ m / s ] are Section 1 and Section 2 speeds, respectively.
Eventually, through the following equations, the ideal mass flowrate m ˙ i d e a l   [ kg / s ] can be calculated:
V ˙ i d e a l =   m ˙ i d e a l ρ
in which ρ   [ kg / m 3 ] is the density of the fluid.
The previous equations allowed calculating the ideal flowrate m ˙ i d e a l . For the real flowrate, the following two coefficients are necessary:
V 1 = V ˙ i d e a l A 1       and       V 2 = V ˙ i d e a l A 2 p 1 p 2 = 1 2 ρ [ ( V ˙ i d e a l A 2 ) 2 ( V ˙ i d e a l A 1 ) 2 ] V ˙ i d e a l = m ˙ i d e a l ρ }       m ˙ i d e a l = A 2 2 ρ | p 1 p 2 | 1 ( A 2 A 1 ) 2
  • Coefficient of discharge C D (non-dimensional coefficient): it takes into account friction, viscosity, and turbulence effects;
  • Expansion factor ϵ (non-dimensional coefficient): it takes into account fluid compressibility.

4.1.1. Coefficient of Discharge Calculation

The coefficient of discharge C D can be calculated through formulations based on the following equations:
β = D 2 D 1 ; R e D = V 1 D 1 v = 4 m ˙ i d e a l π μ D 1 ( Reynolds   number   for   shrinkage ) .
The formula for calculating C D depends on the type of the Air Nozzle [15]. In particular, we have the following different formulations:
ISA   1932   Nozzle : C D = 0.9900 0.2662 β 4.1 [ 0.00175 β 2 0.0033 β 4.15 ( 10 6 R e D ) 1.15 ]
Long   radius   Nozzle : C D = 0.9965 0.00653 β 0.5 ( 10 6 R e D ) 0.5 .
Eventually, the formula for calculating the real mass flowrate m ˙ r e a l is the following:
m ˙ r e a l = C D A 2 2 ρ | p 1 p 2 | 1 ( A 2 A 1 ) 2 = C D 1 β 4 π D 2 2 4 2 ρ | p 1 p 2 | = C π D 2 2 4 2 ρ | p 1 p 2 |
in which C D 1 β 4 is defined as the flow shrinkage coefficient C (non-dimensional).

4.1.2. Expansion Factor calculation

In order to take into account the fluid compressibility, the expansion factor ϵ must be introduced as follows:
ϵ   =   ( k τ 2 / k k 1 ) ( 1 β 4 1 β 4 τ 2 / k ) ( 1 τ ( k 1 ) / k 1 τ )                   if     p 2 p 1     0.75
where:
  • τ   =   p 2 p 1 ;
  • k is the adiabatic expansion coefficient, which is equal to the constant-pressure specific heat and constant-volume specific heat ratio ( k   =   c p / c v ) .

4.1.3. Real Flowrate Calculation

The previous equations can be combined to obtain the formulation for calculating the air flowrate as follows m ˙ n o z z l e :
m ˙ nozzle = ϵ m ˙ real = C ϵ π D 2 2 4 2 ρ air | p 1 p 2 | = 2 C π D 2 2 4 ρ air | Δ p air   nozzle | = = K ρ air | Δ p air   nozzle |
in which:
  • K = 2 C π D 2 2 4 [ m 2 ] is the Air Nozzle coefficient;
  • ρ a i r   [ kg / m 3 ] is the air density, equal to ρ a i r = 100 p 0 R a i r ( t 1 + 273.15 ) , where R a i r   =   287.04   J / kgK ;
  • | Δ p a i r   n o z z l e |   [ mbar ] is the orifice plate differential pressure. Usually, values are between 100 and 250   mmH 2 O .
In order to take into account the measurement units used for | Δ p a i r   n o z z l e |   [ mbar ] , it is necessary to introduce a small change in the previous formula as follows:
m ˙ n o z z l e = K 100 ρ a i r | Δ p a i r   n o z z l e | [ Pa ] .
Eventually, it is worth noting that the air flowrate that effectively reaches the engine is lower, due to compressor losses. Moreover, a situation with multiple inlets to the turbocharger ( n ° T C ) could be the real one. As a result, the following equation is the most correct one to evaluate the real air flowrate m ˙ a i r   e n g :
m ˙ a i r   e n g = ( 1 T C s e a l i n g   a i r 100 ) n ° T C m ˙ n o z z l e
in which T C s e a l i n g   a i r   [ % ] represents the air loss within the compressor seal. Its approximated value is usually equal to 1.5%.
Finally, the exhaust mass flowrate m ˙ e x h can be calculated as follows:
m ˙ e x h = 3600 m ˙ a i r   e n g + 10 3 P e n g B S F C t o t a l [ kg / h ]
where B S F C t o t a l   [ g / kWh ]   i s the total brake-specific fuel consumption, which is defined below. Indeed, for the flow equivalence principle, the exhaust mass flowrate is equal to the air mass flowrate plus the burnt-fuel flowrate.

4.2. Carbon Balance Method

The Carbon Balance method allows calculating the exhaust flowrate by means of the concentration of some components of the exhaust gases and the fuel concentration and consumption. The approach is regulated within the revision of MEPC (Marine Environment Protection Committee) Resolution 176(58)|2008 [20], of the NOx Technical Code 2008 [16] and partially within the EU 2016/802 [21].
In order to apply this method, both fuel analyses and sensors able to evaluate exhaust gases and calculate engine consumptions are necessary.
In the following sections, thermodynamics quantities will be distinguished on the basis of their water content in two categories: wet and dry quantities. This is important because some sensors could be damaged while measuring quantities (e.g., oxygen detectors).
Another important necessary distinction is related to dual-fuel engines: indeed, for some quantities, the fuel gas duration injection will be reported. In particular, a fuel gas duration injection < 1 µs implies an engine working with fuel oil; in such a way, the transition between fuel and gas will not be considered.

4.2.1. Demonstration of the Carbon Balance Method

The basic assumption of the Carbon Balance method implies that the exhaust gas components are in the form of CO, CO2, and HC (HydroCarbon). Consequently, the flowrate of wet exhaust gases is defined by the following equation [16]:
m ˙ e x h ,   w e t   =   m ˙ f u e l [ 1.4 F u e l m i x   c o n t c F u e l m i x   c o n t c 1 1.293 ( 1.4 F u e l m i x c o n t C f C + 0.08936 F u e l m i x   c o n t H 1 ) + F S C f c f c + 0.08936 F u e l m i x   cont H 1 ] ( 1 + H a 1000 ) + 1 [ kg / h ] .
This formula is based on the principle for which the carbon-inlet quantity is equal to the carbon-outlet quantity: Carbon in = Carbon out. The following sections will clearly define the meaning of the two quantities.
Carbon in
The carbon-inlet quantity is calculated through the following formula:
C a r b o n i n = m ˙ f u e l   F u e l m i x   c o n t   c
where m ˙ f u e l   [ kg / h ] is the net fuel flowrate and F u e l m i x   c o n t   C   [ %   mass ] is the carbon content in the fuel mixture, which is calculated as follows:
F u e l m i x   c o n t   C   = =   { F u e l o i l   c o n t   C [ B S F C F u e l   o i l l o s s   +   B S F C P i l o t   o i l ] B S F C t o t a l if   Fuel   gas   duration   injection     <   1   μ s F u e l o i l   c o n t   C B S F C P i l o t   o i l   +   F u e l g a s   c o n t   C B S F C F u e l   g a s B S F C t o t a l if   Fuel   gas   duration   injection     1   μ s
in which:
  • F u e l o i l   c o n t   C   [ %   mass ] is the carbon content in the primary fuel, provided by fuel tests;
  • B S F C F u e l   o i l l o s s   [ g / kWh ] is the brake-specific fuel consumption minus losses and is calculated through the following formula:
    B S F C F u e l   o i l l o s s = 10 3 m ˙ f u e l P e n g i n e = 10 3 m ˙ f u e l , t o t a l 60 m f u e l , l o s s t f u e l   l o s s P e n g i n e
where:
  • m ˙ f u e l   [ kg / h ] is the net fuel flowrate;
  • m ˙ f u e l , t o t a l   [ kg / h ] is the fuel flowrate, losses included;
  • m f u e l , l o s s   [ g ] is the mass of fuel lost;
  • P e n g   [ kW ] is the engine power;
  • t f u e l   l o s s   [ min ] is the losses measurement time.
  • B S F C P i l o t   o i l   [ g / kWh ] is the brake-specific pilot fuel consumption, which is calculated with the following formula:
    B S F C P i l o t   o i l   =   m ˙ P i l o t   o i l P e n g
in which m ˙ P i l o t   o i l   [ g / h ] is the pilot fuel flowrate and P e n g   [ kW ] is the engine power.
F u e l g a s   c o n t   C   [ %   mass ] : fuel gas carbon content;
B S F C F u e l   g a s   [ g / kWh ] : brake-specific fuel gas consumption, equal to B S F C F u e l   g a s   =   10 3 m ˙ g a s P e n g i n e where m ˙ g a s   [ kg / h ] is the gas flowrate and P e n g   [ kW ] is the engine power.
B S F C t o t a l   [ g / kWh ] is the total brake-specific fuel consumption, which is calculated through the following formulae:
B S F C t o t a l   =   { B S F C F u e l   o i l l o s s   +   B S F C P i l o t   o i l if   Fuel   gas   duration   injection     <   1 μ s B S F C F u e l   g a s   +   B S F C P i l o t   o i l if   Fuel   gas   duration   injection     1   μ s
Carbon out
The carbon-outlet quantity is calculated through the following formula:
C a r b o n o u t   =   C a t o m i c   m a s s   [ C O 2   e m i s s i o n s C O 2   m o l a r   m a s s   +   C O e m i s s i o n s   C O m o l a r   m a s s   +   H C e m i s s i o n s H C m o l a r   m a s s ]
in which:
  • C a t o m i c   m a s s   =   12.011   g is the relative carbon mass;
  • C O 2   e m i s s i o n s , C O e m i s s i o n s and   H C e m i s s i o n s   [ g / h ] are the substances flowrate emissions;
  • C O 2   molar   mass ,   C O m o l a r   m a s s and   H C m o l a r   m a s s [ g / mol ] are the molar masses of pollutant.
Emissions can be calculated through the following general equation:
P o l l u t a n t e m i s s i o n s = P o l l u t a n t m o l a r   m a s s P o l l u t a n t m o l a r   v o l u m e ρ e x h P o l l u t a n t c o n c e n t r a t i o n m ˙ e x h
where P o l l u t a n t m o l a r   v o l u m e   [ m 3 / mol ] is the molar mass and the molar volume of pollutant CO2, CO, or HC, respectively. It is worth noting that the molar mass/molar volume ratio of a pollutant is equal to its density, as shown below:
ρ p o l l u t a n t   =   P o l l u t a n t m o l a r   m a s s   P o l l u t a n t m o l a r   v o l u m e .
As a result, by substituting the density in the carbon out formulation, the following equation is obtained:
C a r b o n o u t     =   C a t o m i c   m a s s m ˙ e x h ρ e x h   [ C O 2   c o n c e n t r a t i o n   C O 2   m o l a r   v o l u m e     +   C O c o n c e n t r a t i o n   C O m o l a r   v o l u m e   +   H C c o n c e n t r a t i o n H C m o l a r   v o l u m e ]
Exhaust gas density calculation
The exhaust gas density is equal to the ratio between mass flowrate m ˙ e x h   [ kg / h   ] and volumetric flowrate V ˙ e x h   [ m 3 / h ] :
ρ e x h   =   m ˙ e x h V ˙ e x h .
The volumetric flowrate V ˙ e x h can be calculated by summing up the volumetric flowrate of all the exhaust gas components; however, since they are dependent mainly on CO, CO2, and HC content, the formula can be simplified as follows:
V ˙ e x h   =   V ˙ C O   +   V ˙ N O   +   V ˙ N O 2   +   V ˙ H C   +   V ˙ H 2 O   +   V ˙ C O 2   +   V ˙ O 2   +   V ˙ N 2   +   V ˙ S O 2 V ˙ e x h     V ˙ C O   +   V ˙ H C   +   V ˙ H 2 O   +   V ˙ C O 2 .
The last expression of V ˙ e x h implies that the exhaust gas volumetric flowrate is calculated by taking into account also the water content (i.e., it is a wet quantity). To obtain the dry exhaust gas volumetric flowrate, the water contribution V ˙ H 2 O must be deducted as follows:
V ˙ e x h , d r y   =   V ˙ e x h     V ˙ H 2 O .
The volumetric flowrate for CO and HC can be evaluated through the single components concentrations as follows:
V ˙ C O   =   V ˙ e x h C O c o n c e n t r a t i o n     ;       V ˙ H C   =   V ˙ e x h H C c o n c e n t r a t i o n .
On the other hand, as regards the volumetric flowrate for H2O and CO2, the following formulae must be used:
V ˙ H 2 O = m ˙ c o m b   a i r , w e t C o m b   a i r c o n t   H 2 O H 2 O m o l a r   v o l u m e H 2 O m o l a r   m a s s + m ˙ f u e l F u e l m i x   c o n t   H 2 O H 2 O m o l a r   v o l u m e 2 H a t o m i c   m a s s
V ˙ C O 2 = m ˙ c o m b   a i r , w e t C o m b   a i r c o n t   C O 2 C O 2   m o l a r   v o l u m e C O 2   m o l a r   m a s s + m ˙ f u e l F u e l m i x   c o n t   C   C O 2   m o l a r   v o l u m e C a t o m i c   m a s s
These equations are based on the principle for which water and carbon dioxide are contained also in the combustion air; therefore, to consider it, a first addend, whose factors are described below, is necessary:
  • m ˙ c o m b   a i r , w e t   [ kg / h ] is the wet mass flowrate of the combustion air;
  • C o m b   a i r c o n t   H 2 O   and   C o m b   a i r c o n t   C O 2 are the concentration of water and carbon dioxide in the combustion air, respectively;

4.2.2. Carbon Balance Quantities

The quantities in the Carbon Balance method are the following:
  • m ˙ f u e l   [ kg / h ] : fuel flowrate;
  • f C : emission carbon factor. It is a coefficient that allows converting activity data into GHG emissions and is calculated as follows:
    f C = 0.5441 ( C O 2   c o n c e n t r a t i o n C O 2   a m b   a i r   c o n c e n t r a t i o n ) + C O concentration 18522 + H C c o n c e n t r a t i o n 17355
    where the dry quantities are the C O 2   a m b   a i r   c o n c e n t r a t i o n [%] (concentration of CO2 in the ambient air), C O 2   c o n c e n t r a t i o n   [ % ]   and   C O c o n c e n t r a t i o n   [ ppm ] (concentration of CO2 and CO in the exhaust gas, respectively), while H C c o n c e n t r a t i o n   [ ppm ] (concentration of HC in the exhaust gas) is a wet quantity. Since the carbon dioxide can be found in larger quantities than the other components, it is usually measured in percentage;
  • H a [ g H 2 O / kg air , dry ] : absolute moisture in the aspiration air, defined as the ratio between steam density ρ s and dry air density ρ a i r , d r y at the same temperature;
  • F u e l m i x   c o n t   C   [ %   mass ] : carbon content in the fuel mixture;
  • F u e l m i x   c o n t   H   [ %   mass ] : hydrogen content in the fuel mixture;
  • F S C : fuel dry specific constant for exhaust gases, calculated as follows:
    F S C = 0.055593 F u e l m i x   c o n t   H + 0.008002 F u e l m i x   c o n t   N + 0.0070046 F u e l m i x   c o n t   O
The basic assumption of the Carbon Balance method is flawed: indeed, since it considers only the emissions based on carbon (CO, CO2 and HC), it eliminates all the other components based on nitrogen and oxygen (e.g., nitrogen oxides NOX and the oxygen itself). Consequently, corrective coefficients able to consider all the other fuel components are necessary:
  • F u e l m i x   c o n t   N   [ %   mass ] : nitrogen content in the fuel mixture;
  • F u e l m i x   c o n t   O   [ %   mass ] : oxygen content in the fuel mixture.

4.3. Emissions Calculation

The values of pollutant concentration measured by the emission sensors are the basis for the calculation of emissions, which are calculated through the following formula:
E m i s s i o n [ g kWh ] = 10 3 m ˙ e x h E m i s s i o n c o n c e t r a t i o n , d r y μ k d t w P e n g
Here, two new quantities are introduced:
  • k d t w : dry to wet coefficient, used to convert a dry quantity to a wet one by means of the following formulation provided in [16]:
    k d t w =   ( 1   +   0.005   ×   11.9164 F u e l m i x   c o n t   H F u e l m i x   c o n t   C ( C O 2   c o n c e t r a t i o n   +   C O c o n c e n t r a t i o n ) + 1.608 H a 1000 + 1.608 H a + 0.00794 0.01 0.5 × 11.9164 F u e l m i x   c o n t   H F u e l m i x   c o n t   C C O c o n c e t r a t i o n ( C O 2   c o n c e t r a t i o n + C O c o n c e t r a t i o n ) ( 3 C O 2   c o n c e t r a t i o n + C O c o n c e t r a t i o n ) ) 1
  • μ : component specific factor, equal to the ratio between gas component density (Table 2) and exhaust gas density. The NOX Technical Code [16] provides specific tables containing μ values (Table 3).
Formula (32) for the calculation of emissions is applicable to all the pollutant emissions except for nitrogen oxides, due to their dependency on ambient air conditions. A coefficient that takes into account ambient air temperature and moisture is introduced as follows:
k h u m N o x = 1 1 0.012 ( H a 10.71 ) 0.00275 ( t 1 24.85 ) + 0.00285 ( T 3 T 3 n o m )
where T 3   n o m is the nominal air receiver temperature. The coefficient k h u m   N O X applies to the moisture range 0 25 g H 2 O / kg a i r , d r y .
As a result, the NOX emissions can be calculated as dry quantities with the formula:
NO X [ g kWh ] = 10 3 N O X   c o n c e n t r a t i o n k d t w k h u m   N O X m ˙ e x h μ N O X P e n g
The NO X emissions calculated as wet quantities can be obtained omitting the coefficient k d t w in the previous Formula (35).

4.4. Emission SubRoutine

The Emission SubRoutine (Figure 4) is the part of the code dedicated to managing the user’s choices in terms of modules for emissions detection within the Test Bed. The current analyses have identified three modules for the emissions detection:
  • CO/CO2/O2 module;
  • NOX module;
  • THC module.
As mentioned previously, in this Routine, the user is called to answer the presence of every above-listed module in the SN (in Figure 4, the red boxes represent Online data, which are available for the parsing, while the blue ones represent Offline data, to insert as a constant value).
Finally, as regards the NOX and THC modules, the user is called to answer another important question about the method used for acquiring emissions values. Indeed, the measurement of these pollutants can be performed either with or without the water contribution; in the first case, the software refers to wet measurements, whereas in the second one, the software refers to dry measurements. For practical purposes, the main difference between these two cases lies in the k d t w coefficient that is present or not in Formula (32).

5. Tool–First Validation

In order to make a first validation about the correct working of the software, a comparison between the results obtained from both the tool and the WIT offline procedure, based on the same engine configuration, has been carried out.
The configuration used for this comparison (Table 4) is about a Test Bed with an in-line four-stroke engine W6L50DF, a two-stage ABB turbocharger and, as equipment, an Air Nozzle and emission machine with all the pollution detection modules.
The analyzed operating point is at 110% of the engine nominal load in only gas mode.
To perform this comparison, the pilot contribution to the combustion and the nitrogen and oxygen values in the fuel gas were not considered due to their negligible values.
Once the Set Up Environment has been properly generated (Table 5) and the parsing has been performed, the software determined which functions had to be calculated, verifying whether the Border Conditions were fulfilled or not. On the basis of the selected engine configuration, the software established 111 calculable functions out of 584.
Table 6, Table 7 and Table 8 presents the output of the first validation of the tool determined by both Air Nozzle and Carbon Balance methods, compared with the same features evaluated by the WIT offline procedure.
The root mean square error of all features evaluated with the tool and the WIT offline procedure is equal to 1.46%, which is a more than acceptable result given the extent of the analyzed quantities.

6. Conclusions

The paper presents a tool for the evaluation of the performance in-live of turbocharged diesel engines for marine application. The tool has been developed in G-Language within the LABVIEW environment; its management is very simple and rapid, even if it requires many formulas in order to reliably describe the thermodynamic phenomena resulting from supercharged combustion inside the engine. One of the main challenges fought has concerned a process of synthesis and selection of the input variables, which are in large part to be acquired from the network of sensors installed on the engine. In order to fully describe the engine performance, 584 different formulas were used in the different routines and subroutines. In the present work, the routine for evaluating air pollutant emissions in terms of NOX, SOX, COX and HC has been completely described. In particular, the two methods that are used (i.e., the Air Nozzle approach and Carbon Balance approach) have been thoroughly presented.
The great novelty offered by this tool is the ability of processing all the data detected by the sensor network in-live in order to monitor and dynamically optimize the TDEs’ operation, whereas the common practice involves the collection of performance data and then their off-line processing. The connection between the presented tool and the automation of the TDEs is currently under development. In this way, the software, in the forthcoming future, will be able to automatically act on the electromechanical processes of TDEs from the perspective of Performance improvement. All this has been made possible thanks to the commitment of the engine manufacturers in the implementation of the principles of Industry 4.0. In particular, in order to offer the market increasingly cheaper products (in terms of fuel consumption) and maintenance costs, the new engines are equipped with a dense network of sensors that is capable of measuring all the quantities related to the definition of the engine performance. The processing of this large amount of in-live data and the further connection with the various actuators able to perform all the adjustments on the engine makes the TDEs of the latest generation very competitive. In recent years, the competition among the various global marine engine manufactures is in fact being fought not only from an economic point of view, but also in the continuous research for new technological solutions that are capable of satisfying all the international regulations issued in order to preserve the global environment.
In order to reduce the release of polluting substances into the atmosphere, many new technologies have been invented and applied in ship propulsion systems. The use of alternative fuels to HFO has been studied, and eventually, LNG can be considered a valid option (even if in this specific case, the bunkering logistics network is not adequately functional at a global level). Hydrogen and ammonia seem very promising, but currently, there is no availability of the large quantities needed to refuel ships and the cost is orders of magnitude greater than that of HFO; moreover, ship propulsion systems would become even more complex than they are.
In this scenario, shipowners need to make adequate profits in order to be able to carry out their service continuously. They cannot operate their ships unless they comply with the international regulations in force. Therefore, the use of the tool described in this paper will certainly benefit them because TDEs working near the efficiency peak will reduce fuel consumption and maintenance, reducing the periods of inactivity of the ship. Future developments of this work (already under consideration) will concern the live visualization of the operating point of the compressor and the turbine within the relative map and the addition of further components for the detection of PM and lightning. Furthermore, the software will be equipped with a digital twin routine: by comparing the real behavior of the engine with the expected one, it will be possible to predictively identify any wear and tear of components inside the engine and set an artificial intelligence to autonomously modify the operating parameters with a view to optimizing consumption, performance, energy, and maintenance.

Author Contributions

Conceptualization, M.D. and V.B.; methodology, M.D. and V.B.; software, M.D.; validation, M.D. and S.B.; formal analysis, M.D. and S.B.; investigation, M.D.; resources, A.M.; data curation, M.D.; writing—original draft preparation, M.D. and S.B.; writing—review and editing, A.M. and V.B.; visualization, M.D., S.B., A.M. and V.B.; supervision, A.M. and V.B.; project administration, V.B.; funding acquisition, A.M. and V.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by 2014–2020 Interreg V-A, Italy-Croatia CBC Programme, by means of the 2017 Standard call research project “METRO—Maritime Environment-friendly TRanspOrt systems”, Priority Axis: Maritime transport.

Acknowledgments

The authors gratefully acknowledge Wärtsilä Italia S.p.a. for the high skilled technical support received in order to carry out the studies described in this work.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

m ˙ a i r Incoming air mass flowrate to the compressor [kg/s]
C D Discharge coefficient
ϵ Expansion coefficient
m ˙ e x h Exhaust gas mass flowrate [kg/s]
m ˙ a i r   e n g Incoming air mass flowrate to the engine [kg/s]
η S V Compressor efficiency [%]
η T Turbine efficiency [%]
η T l Overall efficiency of turbocharger [%]
Π V Compressor pressure ratio
Π T Turbine pressure ratio
n 298 Corrected turbocharger speed [rpm]
B S F C   Brake-Specific Fuel Consumption [g/kWh]
B S F C c o r r ISO correction Brake-Specific Fuel Consumption [g/kWh]
B S F C t o t a l Total Brake-Specific Fuel Consumption [g/kWh]
m ˙ i d e a l Ideal mass flowrate [kg/s]
β Section diameters ratio
R e D Reynolds number for shrinkage
m ˙ r e a l Real mass flowrate [kg/s]
C Flow shrinkage coefficient
τ Pressure ratio
k Adiabatic expansion coefficient
K Air Nozzle coefficient [m2]
ρ a i r Air density [kg/m3]
| Δ p   a i r   n o z z l e | Orifice plate differential pressure [mbar]
m ˙ n o z z l e Air Nozzle incoming mass flowrate to the compressor [kg/s]
T C   s e a l i n g   a i r Air loss within the compressor seal [%]
n ° T C Number of turbochargers
P e n g Engine power [kW]
C a r b o n i n   Carbon-inlet quantity [kg/h]
m ˙ f u e l Net fuel flowrate [kg/h]
F u e l m i x   c o n t   C Carbon content in fuel mixture [% mass]
F u e l o i l   c o n t   C Carbon content in the primary fuel [% mass]
B S F C F u e l   o i l l o s s Brake-Specific Fuel Consumption minus losses [g/kWh]
m ˙ f u e l ,   t o t a l Fuel flowrate, losses included [kg/h]
m f u e l   l o s s Mass of fuel lost [g]
t f u e l   l o s s Losses measurement time [min]
B S F C P i l o t   o i l Pilot Brake-Specific Fuel Consumption [g/kWh]
m ˙ P i l o t   o i l Pilot fuel flowrate [g/h]
F u e l g a s   c o n t   C Fuel gas carbon content [% mass]
B S F C F u e l   g a s   Brake-Specific Fuel gas Consumption [g/kWh]
m ˙ g a s Fuel gas flowrate [kg/h]
C a r b o n o u t Carbon-outlet quantity [kg/h]
C a t o m i c   m a s s Carbon relative atomic mass [g]
P o l l u t a n t e m i s s i o n s Pollutant emissions [g/h]
P o l l u t a n t m o l a r   m a s s Pollutant molar mass [g/mol]
P o l l u t a n t m o l a r   v o l u m e Pollutant molar volume [m3/mol]
P o l l u t a n t c o n c e n t r a t i o n Pollutant concentration
ρ e x h Exhaust gas density [kg/m3]
V ˙ e x h Wet exhaust gas volumetric flowrate [m3/h]
V ˙ e x h ,   d r y Dry exhaust gas volumetric flowrate [m3/h]
m ˙ c o m b   a i r , w e t Wet mass flowrate of combustion air [kg/h]
C o m b   a i r c o n t   H 2 O Water content in combustion air [% mass]
C o m b   a i r c o n t   C O 2 Carbon dioxide content in combustion air [% mass]
ρ H 2 O Water density [kg/m3]
ρ C O 2 Carbon dioxide density [kg/m3]
f C Emission carbon factor
H a Absolute moisture in aspiration air [ g H 2 O / kg air   dry ]
F u e l m i x   c o n t   H Hydrogen content in fuel mixture [% mass]
F S C Fuel dry specific constant for exhaust gas
F u e l m i x   c o n t   N Nitrogen content in fuel mixture [% mass]
F u e l m i x   c o n t   O Oxygen content in fuel mixture [% mass]
k d t w Dry to wet coefficient
μ Component specific factor
k h u m   N O x Ambient-air temperature and moisture coefficient
T 3   n o m Air receiver temperature [°C]

Abbreviations

ABBAsea Brown Boveri
CIMACInternational Council on Combustion Engines
COCarbon Monoxide
CO2Carbon Dioxide
ECREffective Compression Ratio
EUEuropean Union
GHGGreenHouse Gases
HCHydroCarbon
HFOHeavy Fuel Oil
ICEInternal Combustion Engine
IMOInternational Maritime Organization
ISOInternational Organization for Standardization
LABVIEWLaboratory of Virtual Instrumentation and Engineering Workbench
LHVISOISO Lower Heating Value
LNGLiquefied Natural Gas
MARPOLInternational Convention for the Prevention of Marine Pollution
MEPCMarine Environment Protection Committee
NOXNitrogen Oxides
O2Oxygen
PLCProgrammable Logic Controller
SCRSelective Catalytic Reduction
SEEMPShip Energy Efficiency Management Plan
SNSensor Network
SO2Sulphur Dioxides
TCTurbocharger
TDETurbocharged Diesel Engine
THCTotal HydroCarbon
WITWärtsilä Italia

Appendix A

In the paper, the following WIT Convention is used to standardize the terms/words related to thermodynamic quantities (i.e., pressure and temperature [°C]) within performance documents. It was chosen for the following reasons:
  • The WIT standard simplifies reading formulae in internal documentation;
  • On test engines, the sensors for thermodynamic quantities (e.g., pressure gauges and thermocouples) are installed on the basis of the WIT convention, and then additional calculations are not necessary;
  • The WIT standard does not distinguish between four-stroke and two-stroke engines. The only difference in terms of number of sensors refers to single-stage and double-stage cases.
Single-stage case
As regards the single-stage case (single turbocharger), with reference to Figure A1, the WIT convention is the following:
  • Subscript “0” marks all the room/atmospheric/ambient-condition quantities;
  • Subscript “1” marks all the compressor-inlet quantities;
  • Subscript “2” marks all the compressor-outlet quantities;
  • Subscript “3” marks all the air-receiver quantities;
  • Subscript “4” marks the exhaust-gas manifold temperature, calculated as the average of all the exhaust-gas manifold temperatures (each single temperature is marked through an additional subscript that specifies the cylinder number). It is worth noting that the exhaust-gas pressure is not calculated due to its significant fluctuation;
  • Subscript “5” marks all the turbine-inlet quantities. It is worth noting that, in this case, the pressure value is stationary and can be measured;
  • Subscript “6” marks all the turbine-outlet quantities.
Figure A1. WIT convention: single-stage case.
Figure A1. WIT convention: single-stage case.
Applsci 10 05707 g0a1
As regards measurement units, the WIT uses the following:
  • Room/atmospheric pressure p 0 [mbar];
  • Compressor-inlet pressure p 1 [mbar];
  • Compressor-outlet pressure p 2 [bar];
  • Receiver pressure p 3 [bar];
  • Turbine-inlet pressure p 4 [bar];
  • Turbine-outlet pressure p 5 [bar];
  • Temperatures [°C].
The standard for measurement units simplifies the measurement process of quantities, but it also has some disadvantages. Indeed, since the SI unit for pressure is Pascal [Pa] and for temperature is Kelvin [K], proper transformations are needed.
Double-stage case
As regards the double-stage case (two turbochargers in series), the WIT convention is the following, and it is sketched in Figure A2:
  • Subscript “0” marks all the room/atmospheric/ambient-condition quantities;
  • Subscript “1LP” marks all the low-pressure compressor-inlet quantities;
  • Subscript “2LP” marks all the low-pressure compressor-outlet quantities;
  • Subscript “1HP” marks all the high-pressure compressor-inlet quantities;
  • Subscript “2HP” marks all the high-pressure compressor-outlet quantities;
  • Subscript “3” marks all the air-receiver quantities;
  • Subscript “4” marks the exhaust-gas manifold temperature, calculated as the average of all the exhaust-gas manifold temperatures;
  • Subscript “5HP” marks all the high-pressure turbine-inlet quantities;
  • Subscript “6HP” marks all the high-pressure turbine-outlet quantities;
  • Subscript “5LP” marks all the low-pressure turbine-inlet quantities;
  • Subscript “6LP” marks all the low-pressure turbine-outlet quantities.
Figure A2. WIT convention: Double-stage case.
Figure A2. WIT convention: Double-stage case.
Applsci 10 05707 g0a2
As regards measurement units, the WIT uses the following:
  • Room/atmospheric pressure p 0 [mbar];
  • Compressor-inlet pressure p 1 L P and p 1 H P [mbar];
  • Compressor-outlet pressure p 2 L P and p 2 H P [bar];
  • Receiver pressure p 3 [bar];
  • High-pressure turbine-inlet pressure p 5 H P [bar];
  • High-pressure turbine-outlet pressure p 6 H P , low-pressure turbine-inlet pressure p 5 L P and low-pressure turbine-outlet pressure p 6 L P [bar];
  • Temperatures [°C].
In order to simplify calculations, the previous convention was slightly modified, in particular with reference to inlet and outlet quantities. Indeed, the double-stage case can be equated to a single-stage case having the following characteristics:
  • Compressor inlet = low-pressure compressor inlet;
  • Compressor outlet = high-pressure compressor outlet;
  • Turbine inlet = high-pressure turbine inlet;
  • Turbine outlet = low-pressure turbine outlet.
In addition, thermodynamic quantities regarding the high-pressure turbine outlet and the low-pressure turbine inlet can be simplified. Since these values are very similar, they are considered equal and marked with the subscript “5LP” (the subscript of low-pressure turbine-inlet quantities).
In conclusion, the modified WIT convention for double-stage cases is presented in Figure A3.
Figure A3. Modified WIT convention: Double-stage case.
Figure A3. Modified WIT convention: Double-stage case.
Applsci 10 05707 g0a3

References

  1. Brenna, M.; Bucci, V.; Falvo, M.C.; Foiadelli, F.; Ruvio, A.; Sulligoi, G.; Vicenzutti, A. A review on energy efficiency in three transportation sectors: Railways, electrical vehicles and marine. Energies 2020, 13, 2378. [Google Scholar] [CrossRef]
  2. International Maritime Organization. International Convention for the Prevention of Pollution from Ships-MARPOL 2017 Consolidated Edition. Available online: http://www.imo.org/en/Publications/Documents/Supplements%20and%20CDs/English/QQE520E_022019.pdf (accessed on 14 February 2020).
  3. International Maritime Organization. Third IMO Greenhouse Gas Study. Available online: http://www.imo.org/en/OurWork/Environment/PollutionPrevention/AirPollution/Documents/Third%20Greenhouse%20Gas%20Study/GHG3%20Executive%20Summary.pdf (accessed on 14 February 2020).
  4. Altosole, M.; Benvenuto, G.; Campora, U.; Laviola, M.; Zaccone, R. Simulation and performance comparison between diesel and natural gas engines for marine applications. Proc. Inst. Mech. Eng. Part M J. Eng. Marit. Environ. 2017, 231, 690–704. [Google Scholar]
  5. Taccani, R.; Ustolin, F.; Zuliani, N.; Pinamonti, P.; Pietra, A. Fuel Cells and Shipping Emissions Mitigation, Technology and Science for Ships of the Future. Proceedings of NAV 2018: 19th International Conference on Ship & Maritime Research, Trieste, Italy, 20–22 June 2018; pp. 885–892. [Google Scholar]
  6. Kim, K.; Roh, G.; Kim, W.; Chun, K. A preliminary study on an alternative ship propulsion system fueled by ammonia: Environmental and economic assessments. J. Mar. Sci. Eng. 2020, 8, 183. [Google Scholar] [CrossRef] [Green Version]
  7. Brynolf, S.; Fridell, E.; Andersson, K. Environmental assessment of marine fuels: Liquefied natural gas, liquefied biogas, methanol and bio-methanol. J. Clean. Prod. 2014, 74, 86–95. [Google Scholar] [CrossRef]
  8. Jaurola, M.; Hedin, A.; Tikkanen, S.; Huhtala, K. Optimising design and power management in energy-efficient marine vessel power systems: A literature review. J. Mar. Eng. Technol. 2019, 18, 92–101. [Google Scholar] [CrossRef] [Green Version]
  9. Muše, A.; Jurić, Z.; Račić, N.; Radica, G. Modelling, performance improvement and emission reduction of large two-stroke diesel engine using multi-zone combustion model. J. Therm. Anal. Calorim. 2020, 141, 337–350. [Google Scholar]
  10. Stoumpos, S.; Theotokatos, G.; Mavrelos, C.; Boulougouris, E. Towards marine dual fuel engines digital twins-integrated modelling of thermodynamic processes and control system functions. J. Mar. Sci. Eng. 2020, 8, 200. [Google Scholar] [CrossRef] [Green Version]
  11. Ntonas, K.; Aretakis, N.; Roumeliotis, I.; Pariotis, E.; Paraskevopoulos, Y.; Zannis, T. Integrated simulation framework for assessing turbocharger fault effects on diesel-engine performance and operability. J. Energy Eng. 2020, 146, 04020023. [Google Scholar] [CrossRef]
  12. Altosole, M.; Campora, U.; Figari, M.; Laviola, M.; Martelli, M. A diesel engine modelling approach for ship propulsion real-time simulators. J. Mar. Sci. Eng. 2019, 7, 138. [Google Scholar] [CrossRef] [Green Version]
  13. Aiello, G.; Giallanza, A.; Vacante, S.; Fasoli, S.; Mascarella, G. Propulsion monitoring system for digitized ship management: Preliminary results from a case study. Procedia Manuf. 2020, 42, 16–23. [Google Scholar] [CrossRef]
  14. ISO 5167-1:2003. In Measurement of Fluid Flow by Means of Pressure Differential Devices Inserted in Circular Cross-Section Conduits Running Full-Part 1: General Principles and Requirements; International Organization for Standardization (ISO): Geneva, Switzerland, 2003.
  15. ISO 5167-3:2003. In Measurement of Fluid Flow by Means of Pressure Differential Devices Inserted in Circular Cross-Section Conduits Running Full-Part 3: Nozzles and Venturi Nozzles; International Organization for Standardization (ISO): Geneva, Switzerland, 2003.
  16. Marine Environment Protection Committee. Amendments to the Technical Code on Control of Emission of Nitrogen Oxides from Marine Diesel Engines (NOx Technical Code 2008). Resolution MEPC 177(58). Available online: http://www.imo.org/en/KnowledgeCentre/IndexofIMOResolutions/Marine-Environment-Protection-Committee-(MEPC)/Documents/MEPC.177(58).pdf#search=Amendments%20to%20the%20technical%20code%20on%20control%20of%20emission%20of%20Nitrogen%20Oxides%20from%20marine%20Diesel%20Engines%20%28NOx%20Technical%20Code%202008%29%2E%20Resolution%20MEPC%20177%2858%29 (accessed on 14 February 2020).
  17. ISO 8178-1:1996. In Reciprocating Internal Combustion Engines—Exhaust Emission Measurement—Part 1: Test-Bed Measurement of Gaseous and Particulate Exhaust Emissions; International Organization for Standardization (ISO): Geneva, Switzerland, 1996.
  18. ISO 8217-4:2010. In Petroleum Products—Fuel (Class F). Specification of Marine Fuels; International Organization for Standardization (ISO): Geneva, Switzerland, 2010.
  19. International Council of Combustion Engines. Turbocharging Efficiency. Definitions and Guidelines for Measurement and Calculation. Available online: https://www.cimac.com/cms/upload/Publication_Press/Recommendations/Recommendation_27_rev_081007.pdf (accessed on 14 February 2020).
  20. Marine Environment Protection Committee. Amendments to the Annex of the Protocol of 1997 to Amend the International Convention for the Prevention of Pollution from Ship, 1973, as Modified by the Protocol of 1978 Relating Thereto. Resolution MEPC 176(58) 2008. Available online: http://www.imo.org/en/OurWork/Environment/PollutionPrevention/AirPollution/Documents/176(58).pdf#search=176%2058 (accessed on 14 February 2020).
  21. European Commision. Directive (EU) 2016/802 of the European Parliament and of the Council of 11 May 2016 Relating to a Reduction in the Sulphur Content of Certain Liquid Fuels. EU 2016/802. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/?qid=1597385743484&uri=CELEX:32016L0802 (accessed on 14 February 2020).
  22. European Commission. Annex 7: Technical and General Requirements Relating to Emission Limits and Type-Approval for Internal Combustion Engines for Non-Road Mobile Machinery C(2016)-8381. Available online: https://ec.europa.eu/transparency/regdoc/rep/3/2016/EN/C-2016-8381-F1-EN-ANNEX-1-PART-1.PDF (accessed on 14 February 2020).
Figure 1. Main block diagram of the software.
Figure 1. Main block diagram of the software.
Applsci 10 05707 g001
Figure 2. Configuration Set Up phase.
Figure 2. Configuration Set Up phase.
Applsci 10 05707 g002
Figure 3. Section shrinkage.
Figure 3. Section shrinkage.
Applsci 10 05707 g003
Figure 4. Emission SubRoutine.
Figure 4. Emission SubRoutine.
Applsci 10 05707 g004
Table 1. Identified output.
Table 1. Identified output.
Field of InterestOutputSymbolUnit
EngineEffective Compression RatioECR/
TurbochargerCompressor efficiency η s V /
Turbine efficiency η T /
Overall efficiency η T I /
Compressor pressure ratio (total to total) Π V /
Turbine pressure ratio (total to total) Π T /
Corrected turbocharger speed n 298 rpm
FuelBrake Specific Fuel Consumption B S F C g/kWh
ISO Correction Brake Specific Fuel Consumption B S F C c o r r g/kWh
Total Brake Specific Fuel Consumption B S F C t o t a l g/kWh
EmissionsNOX emissions NO X g/kWh
CO emissions CO g/kWh
CO2 emissions CO 2 g/kWh
THC emissions THC g/kWh
O2 emissions O 2 g/kWh
SO2 emissions SO 2 g/kWh
Table 2. Gas component density.
Table 2. Gas component density.
N O X C O H C C O 2 O 2 C H 4
ρ c o m p o n e n t
[ kg / m 3 ]
2.0531.250-1.96361.42770.716
Table 3. Component specific factor μ .
Table 3. Component specific factor μ .
ρ e x h a u s t [ kg / m 3 ] N O X C O H C C O 2 O 2 C H 4
Fuel oil1.29430.0015860.0009660.0004820.0015170.0011030.000553
Ethanol
(ED95)
1.27680.0016090.0009800.0007800.0015390.0011190.000561
Natural
Gas
1.26610.0016210.0009870.0005280.0015510.0011280.000565
Propane1.28050.0016030.0009760.0005120.0015330.0011150.000559
Butane1.28320.0016000.0009740.0005050.0015300.0011130.000558
LPG1.28110.0016020.0009760.0005100.0015330.0011150.000559
Gasoline
(E10)
1.29310.0015870.0009660.0004990.0015180.0011040.000553
Ethanol
(E85)
1.27970.0016040.0009770.0007300.0015340.0011160.000559
Table 4. Engine configuration.
Table 4. Engine configuration.
EngineW6L50DF
TurbochargerTwo-stage ABB
Test bed equipmentAir Nozzle
Emission machine:
-
CO, CO2, O2 module;
-
THC module (Wet);
-
NOX module (Wet).

No Wastegate detected
Load109.4%
FuelOnly gas
Table 5. Input values.
Table 5. Input values.
t 0 [°C]35.000 E N G I N E   S P E E D [rpm]599.5 m G A S   N E T [kg/h]1277.900
t 1 [°C]35.100 Z [-]6 C O 2   D R Y   C O N C [%]5.360
t 2   L P [°C]200.200 B O R E [m]0.500 C O D R Y   C O N C [%]163.470
t 1   H P [°C]61.500 S T R O K E [m]0.580 O 2   D R Y   C O N C [%]11.080
t 2 [°C]149.000 C R [-]13 N O X   W E T   C O N C [ppm]131.690
t 3 [°C]55.800 R O D   T O   C R A N K   R A T I O [-]4.167 T H C W E T   C O N C [ppm]631.990
t 4 [°C]532.100 T I V C [rad]−60 t R E L A T I V E   H U M I D I T Y [°C]20.100
t 5 [°C]564.200 k [-]0.630 H R [gH2O/kgairdry]37.400
t 5   L P [°C]479.500 Δ p N [mbar]34.900 C O M P   L E A K L P [%]1.500
t 6 [°C]324.700 n T C [-]1 C O M P   L E A K H P [%]1.000
p 0 [mbar]1013.300 D 2   L P [m]0.500 P E N G [kW]8530
p 1 [mbar]25.980 D 2   H P [m]0.263 L H V G A S [MJ/kg]49.820
p 2   L P [bar]2.670 D 5   L P [m]0.500 A F R G A S [-]16.970
p 1   H P [mbar]2670.000 D 5   H P [m]0.262 T U R B   S P E E D L P [rpm]18011
p 2 [bar]6.280 C G A S   C O N T E N T [%]75.200 T U R B   S P E E D H P [rpm]16393
p 3 [bar]6.220 H G A S   C O N T E N T [%]24.700 t R E C E I V E R   N O M I N A L [°C]50.000
p 5 [bar]3.830 N G A S   C O N T E N T   [%]0.000
p 5   L P [bar]2.040 O G A S   C O N T E N T [%]0.000
p 6 [bar]0.00843 S G A S   C O N T E N T [%]0.000
Table 6. Comparison between the output values–Running parameters. BSFC: Brake-Specific Fuel Consumption, WIT: Wärtsilä Italia.
Table 6. Comparison between the output values–Running parameters. BSFC: Brake-Specific Fuel Consumption, WIT: Wärtsilä Italia.
SymbolStageUnitIn-Live ToolOffline WIT ProcedurePercentage Error
ECR /10.54610.5500.038%
B S F C g/kWh149.812149.8200.005%
B S F C c o r r g/kWh149.273149.2900.011%
n 298 LPrpm295.225295.000−0.076%
HPrpm257.887258.0000.044%
Table 7. Comparison between the output values–Air Nozzle method results.
Table 7. Comparison between the output values–Air Nozzle method results.
SymbolStageUnitIn-Live ToolOffline WIT ProcedurePercentage Error
η s V LP%81.54185.2764.381%
HP%79.75881.9432.740%
η T LP%78.95776.963−2.525%
HP%80.03176.175−4.818%
η T I LP%64.38265.6311.940%
HP%63.83162.420−2.211%
Π V LP/3.7753.7780.079%
HP/1.9921.973−0.954%
Π T LP/3.0043.0070.100%
HP/1.6341.590−2.693%
NO X g/kWh1.1151.109−0.538%
CO g/kWh330.817346.0294.598%
CO 2 g/kWh0.6420.6724.673%
THC g/kWh1.8331.823−0.546%
O 2 g/kWh497.347520.2174.598%
SO 2 g/kWh0.0000.0000.000%
Table 8. Comparison between the output values–Carbon Balance method results.
Table 8. Comparison between the output values–Carbon Balance method results.
SymbolStageUnitIn-Live ToolOffline WIT ProcedurePercentage Error
η s V LP%81.50780.418−1.336%
HP%79.72881.7392.522%
η T LP%79.00779.2970.367%
HP%80.15377.379−3.461%
η T I LP%64.39663.769−0.974%
HP%63.90463.249−1.025%
Π V LP/3.7703.584−4.934%
HP/1.9901.970−1.005%
Π T LP/3.0043.0070.100%
HP/1.6321.630−0.123%
NO X g/kWh1.0881.0930.460%
CO g/kWh322.808340.9365.616%
CO 2 g/kWh0.6270.6625.581%
THC g/kWh1.7881.7960.447%
O 2 g/kWh485.307512.5615.616%
SO 2 g/kWh0.0000.0000.000%

Share and Cite

MDPI and ACS Style

Dodero, M.; Bertagna, S.; Marino’, A.; Bucci, V. Performance In-Live of Marine Engines: A Tool for Its Evaluation. Appl. Sci. 2020, 10, 5707. https://doi.org/10.3390/app10165707

AMA Style

Dodero M, Bertagna S, Marino’ A, Bucci V. Performance In-Live of Marine Engines: A Tool for Its Evaluation. Applied Sciences. 2020; 10(16):5707. https://doi.org/10.3390/app10165707

Chicago/Turabian Style

Dodero, Matteo, Serena Bertagna, Alberto Marino’, and Vittorio Bucci. 2020. "Performance In-Live of Marine Engines: A Tool for Its Evaluation" Applied Sciences 10, no. 16: 5707. https://doi.org/10.3390/app10165707

APA Style

Dodero, M., Bertagna, S., Marino’, A., & Bucci, V. (2020). Performance In-Live of Marine Engines: A Tool for Its Evaluation. Applied Sciences, 10(16), 5707. https://doi.org/10.3390/app10165707

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