Use of LNG Cold Potential in the Cogeneration Cycle of Ship Power Plants

: This paper presents the results of a numerical study on the parameters that a ﬀ ect the e ﬃ ciency of the cogeneration cycle of a ship’s power plant. The e ﬃ ciency was assessed based on the excess power (N gen. ) of a free turbine, operated with the inﬂow of gaseous nitrogen, which was used to generate electricity. A mathematical model and simulation of the regenerative cycle were created and adjusted to operate with a dual-fuel (diesel-liquid natural gas (LNG)) six-cylinder four-stroke engine, where the energy of the exhaust gas was converted into mechanical work of the regenerative cycle turbine. The most signiﬁcant factors for N gen . were identiﬁed by parametrical analysis of the cogeneration cycle: in the presence of an ‘external’ unlimited cold potential of the LNG, N gen. determines an exhaust gas temperature Te g of power plant; the pressure of the turbo unit and nitrogen ﬂow are directly proportional to N gen . When selecting the technological units for cycle realization, it is rational to use high ﬂow and average π T pressure (~3.0–3.5 units) turbo unit with a high adiabatic e ﬃ ciency turbine. The e ﬀ ect of the selected heat exchangers with an e ﬃ ciency of 0.9–1.0 on N gen. did not exceed 10%. With LNG for ‘internal’ use in a ship as a fuel, the lowest possible temperature of N 2 is necessary, because each 10 K increment in N 2 entering the compressor decreases N gen. by 5–8 kW, i.e., 5–6%.


Ship Propulsion Plant Environmental Regulations
In 2015, 193 countries adopted the 2030 Agenda for Sustainable Development and its 17 Sustainable Development Goals (SDGs) because greenhouse gas emissions are attaining their highest levels in history [1]. Global emissions have more than doubled from the early seventies and increased by approximately 40% since 2000 [2].
The transport sector is one of the largest consumers of energy from burning petroleum products and is one of the largest polluters of the environment with exhaust gas. Virtually all (95%) of the world's transportation energy comes from petroleum-based fuels, largely gasoline and diesel [3].
In 2016, global transport accounted for a quarter of total emissions at approximately 8 GT of CO 2 , which was 71% higher than in 1990 [2]. In 2016, the transport sector contributed 27% of total EU-28 The EEDI requires a minimum energy efficiency level per capacity mile (e.g., tonne mile) for different types of ship and size segments.
In summary, these measures have the common goal of reducing CO 2 emissions and other pollutants from ships, but they do not specify concrete methods of achieving this. Such methods are also not reflected in the EEDI accounting reports. The increasing amount of ecological and energy efficiency requirements for marine engines and the rigidity of standards causes shipbuilders and engine manufacturers to seek cost-effective methods to satisfy them [19][20][21].
MEPC plans to adjust the minimum energy efficiency level every five years. This can cause technological difficulties for ship design companies and engine manufacturers in the future or it can become a technically unsolvable problem for ships powered by conventional fuels [16].

Ship Propulsion Decarbonization Technologies
Currently, seagoing ships use relatively cheap heavy fuel with a sulfur content not exceeding 3.5% or petroleum distillates [22,23]. Switching to more expensive low-sulfur fuel or installing scrubbers to treat the exhaust gas of an internal combustion engine involves significant operating costs [24,25].
The most effective method to solve the compliance and cost dilemma is to use gaseous and alternative fuels [26,27]. Among the most prospective fuels are natural gas, hydrogen, methane-based biogas, and liquefied propane and butane gases.
For seagoing ships, which often have power plants of several tens of megawatts, require reduction of emissions by several times, and require bunkering volumes associated with the supply of alternative fuels, the most rational marine fuel would be liquefied natural gas (LNG) [28][29][30].
In a deeper examination of studies, research results on the single-cylinder marine engine Wärtsilä 20DF were observed as particularly significant [36] to implementing Tier III requirements in shipping. The authors of the paper presented principles on rationally organizing the indication process of the internal combustion engine to achieve acceptable ecological standards without using secondary technologies for waste gas cleaning.
The determining value of the spatiotemporal distribution of a highly reactive pilot portion of injected liquid fuel in the volume of the combustion chamber is justified by the comprehensive experimental studies and calculation methods of mathematical modeling, as a method to control and regulate the dual-fuel fuel combustion process [36]. As a result, a decrease in NO X concentration in exhaust gas from 3000 to 500 mg/Nm 3 was achieved, which was even lower than Tier III normative within a significantly wide operational engine range (P mc = 2.1-1.3 MPa).
Wärtsilä, as a global leader of 4-stroke medium-speed marine engine solutions for the marine and energy markets, maximizes the environmental and economic performance of vessels running on a wide range of liquid and gaseous fuels. Wärtsilä engines can operate on heavy fuel oil (HFO), marine diesel oil (MDO), low viscosity or low-sulfur fuels, LNG, ethane gas (LEG), or petroleum gas (LPG) [37,38]. Moreover, the Wärtsilä 31 engine was awarded a Guinness World Records title for the most efficient 4-stroke diesel engine [39]: the engine consumes, on average, 8-10 g/kWh less fuel compared with the closest competitor across its entire load range, maintenance costs are reduced by approximately 20%, remote access to operational data reduces unscheduled maintenance visits on board, etc. Wärtsilä supplies dual-fuel engines such as the 31DF, 34DF, 46DF, and 50DF that satisfy MARPOL Tier III requirements in the dual-fuel mode and generate power in the range of 1110-17,550 kW.
The widely recognized MAN Diesel & Turbo SE company is a global leader of 2-stroke low-speed marine engines that supplies a variety of engines in terms of power and speed. Their dual-fuel engines also satisfy MARPOL 73/78 Annex VI Tier III ecological requirements. The 2-stroke high power series ME modified engines S70-ME-C8-GI (16,350-26,160 kW), S65-ME-C8-GI (14,350-22,960 kW), and S60-ME-C8-GI (11,900-19,040 kW) can operate on liquid and gaseous fuels and spirits such as methanol, ethanol, and dimethyl ether [40].
In parallel to the dynamically developing segment of highly efficient dual-fuel engines, the number of LNG-powered ships and LNG carriers has been constantly increasing. At the beginning of 2019, the LNG-powered fleet increased globally from 118 LNG-powered vessels in operation in 2017 to 143 LNG-powered vessels in operation-with a further 135 on order and 135 LNG-ready ships either in operation or on order. By March 2019, 18 cruise vessels under construction were LNG-powered [41].
At the end of 2018, 495 LNG carriers above a capacity of 100,000 m 3 (excluding floating storage regasification unit (FSRU), floating storage unit (FSU), and floating liquefied natural gas (FLNG) vessels) were in service, with 55 units having been delivered during the year [42].
Additionally, the increase in orders for LNG powered ships is based on the intensive development of on-shore and floating LNG terminals, LNG storage facilities, specialized infrastructure for supply, storage, and delivery to vessels [42,43].
These facts indicate that adopting LNG as the fuel of choice for ships is a dynamically improving and prospective field.
The use of LNG aids in reducing CO 2 emissions represented by calculating the EEDI and the Energy Efficiency Operational Indicator (EEOI). The result of the EEDI calculation partly depends on the main engine power, auxiliary engine power, implemented innovative energy-saving technologies for electricity generation on shipboards such as for the main engine exhaust gas cogeneration technology, and innovative energy-saving technologies for propulsion improvement such as wind power (Magnus effect, for instance) or solar power.
Cogeneration is a procedure that uses primary energy introduced in a process to produce electric power and thermal energy from secondary energy sources, i.e., exhaust gas and cooling circuits, in the same technological system [44]. Compared with other possibilities of improving the EEDI, secondary utilization of the main engine energy is crucial to increasing energy efficiency [45][46][47].
With an effective efficiency of 50% or more [39,40,48], implemented in modern engines, the reserve for increasing the total efficiency of a ship power plant is 40-35% of fuel energy. Meanwhile, the use of the diesel engine as the main engine of a ship power plant is the best solution for cogeneration systems in terms of energy efficiency. The authors [45] compared various cogeneration ship systems that generate 25 MW of electrical energy: with steam generators and steam turbines (with backpressure and condensation), with a gas turbine, with a diesel engine, and combined systems. They concluded that, for most types of ships with the ratio of the mechanical energy of propulsion to thermal energy as ψ = 2:10, the best typical values of the energy conversion coefficient and efficiency of the ship's power plant are realized in the system with a diesel engine. The exception is oil tankers, where for low ψ caused by the requirement to heat fuel tanks, cogeneration systems based on steam and gas turbines are more rational.

Use of Secondary Propulsion Plant Heat Sources
The structure of cogeneration systems used on ships is determined by the level of technical excellence of the main engine, requirement for electric propulsion, capability to recover energy in a power turbine, types of fuel used, etc.
Thus, the authors [47], evaluating the prospects of increasing the overall efficiency of a power plant in the 90s of the last century, classified ship cogeneration systems in the following categories: • Deep recovery (utilization) system for exhaust gases in a steam turbine for the production of electrical energy and a charge air cooling system for auxiliary requirements of a main power plant (heating of heavy fuel, coolant, lubricating oils); • Heat recovery system for exhaust gases of 4-stroke powerful engines (power higher than 15,000 kW) with a shaft generator, providing electric movement operation mode; • A 'simplified' heat recovery system for exhaust gases of highly efficient 2-stroke diesel engines (characterized by a low temperature of exhaust gases of ≤275 • C) with a shaft generator and exhaust gas power recovery turbine.
Progress in further improving the fuel economy of modern main marine engines [46,49,50] has contributed to the further development of efficient heat recovery systems. In particular, owing to the low temperature of the exhaust gases, which is insufficient to generate steam to the turbogenerator, the Rankine system with low boiling liquids (R141b, R123, R245f, etc.) [49,50], and heat pumps [51,52] have received further development.
In a review [53], a comparative analysis, discussion on features and further prospects of the shipboard application of the following modern heat recovery systems for the main engine, to achieve a lower EEDI, was conducted: thermoelectric generator (TEG), organic Rankine cycle (OCR), closed Brayton cycle and their combinations, turbo compound systems.
The thermoelectric power generator converts thermal energy into electrical energy. The TEG generates a current as heat is applied to a circuit at the junction of two different conductors (Seebeck effect). Advantages of the TEG include maintenance-free, silent operation, and high reliability with no moving and complex mechanical parts [54]. Currently, the TEG is more often used for automotive waste heat recovery owing to relatively low generated energy values (up to 600 W) and a high-temperature (>1000 • C) waste heat demand [55,56]. According to experts, TEGs operating on lower-temperature waste heat from marine engines are to be developed in the future. Another of the most popular solutions of the Rankine cycle is the replacement of water steam by OCRs with low boiling temperatures, 30-50 • C of the heat carrier [49,50], in combination with TEGs to increase the recuperation efficiency of heat exchangers [57].
The marine industry's solution provider Wärtsilä with Turboden of Italy jointly developed the Wärtsilä marine engine combined cycle (ECC) based on ORC technology. Adopted in a cogeneration power plant of a ship, this solution saves 17.1% of the total energy [58].
Another solution by Wärtsilä for high-efficiency waste heat recovery associated with Wärtsilä RT-flex common-rail low-speed marine engines enables up to 12% of the main engine shaft power to be recovered as electrical power for use as additional ship-propulsion power or for shipboard services of container carriers [59]. Exhaust energy is applied in both a steam turbine and an exhaust-gas power turbine to generate electrical power. The generated electrical power can be used either in a shaft motor to aid propulsion or in supplying shipboard services. Additionally, Wärtsilä's waste heat recovery system is one of the most promising methods to improve EEDI. Electric turbo compounding (ETC), turbo compounding-mechanical, and turbo-generator systems regenerate exhaust gas energy in power turbines connected with electric energy generators. In transient engine operating processes, typical for land transport, the power turbine and generator operating as an electric motor provide the necessary acceleration characteristics for the compressor of the main engine [59,60].
In recuperation systems, including those applied in water transport, a closed gas turbine Brayton cycle is used, where the heat released in the combustion chamber and release of energy into the environment as the exhaust gas are replaced by heat exchangers to supply or extract heat. The energy generated by the cycle is determined by the difference between the work produced during expansion in the turbine and compression in the compressor. Therefore, the improvement of the recovery system is more promising if the working fluids have low specific volumes, moreover, if the system works efficiently using exhaust gas low-potential heat of the modern ship engines.
One of the intensely refined solutions is the supercritical CO 2 (S-CO 2 ) Brayton cycle, which offers high plant efficiencies and beneficial economics for a variety of heat sources [61]. For example, using the S-CO 2 cycle on a US navy vessel exhibited an increase in the key engine power indicators by 20-24% [62].
The Brayton cycle, because of its small size and mass and relatively easy integration into transport systems, is increasingly being used as a versatile tool to improve energy efficiency and eco-efficiency and to practically reduce the CO 2 emissions of internal combustion engines of all purposes [63][64][65][66][67].
The thermal performance of the waste heat recovery systems of offshore oil production facilities, including the components of diesel engines, thermal boilers, and waste heat boilers, were analyzed in Reference [63]. The research provided guidelines on improving the waste heat recovery at offshore oil production facilities by applying a waste heat recovery boiler.
An integrated system based on the closed recuperative Brayton cycle (CRBC) for power generation and engine cooling was presented in Reference [64]. In CRBC applied to hypersonic vehicles-next-generation aircraft and spacecraft with broad applications-the combustion heat dissipation is transferred by liquid metal Na and partly converted into electric power.
In automotive transport, the Brayton cycle is used for either traditional oil-based fuel-driven transport means (ICEs) [67] or electric transport [66]. Because of multi-variability, the parametric optimization of the cycle is performed using mathematical modeling methods. A prerequisite is that the engine generating secondary heat must be included in the system's simulation [65,68,69].
For natural gas (NG)-powered transport, even more opportunities to increase the energy efficiency of recovery systems exist. This is primarily applicable to maritime transport because of large engine rooms, which have sufficient space to install the systems of Rankine and Brayton cycles [65][66][67][68].
For NG-fuelled vehicles, gas chemical energy increases by adding the physical energy of liquefaction condition of the NG, which is equal to the difference between the mechanical energy of different pressures in the fuel tank of LNG and compressed gaseous fuel in the engine supply system.
For LNG-powered vehicles, the LNG cold potential, which consists of the amount of heat energy required to vaporize the LNG and heat the NG from storage temperature to engine-supply temperature, is added. In the cogeneration systems used, the physical NG energy is not utilized, i.e., the mechanical energy generated by throttling compressed natural gas (CNG) and the cold potential released by heating the gas to the supply temperature [46].
Known solutions have been developed to harness the potential of the LNG cold potential. For example, the phase production of dry ice or solid CO 2 based on CO 2 re-sublimation, which results in significant reductions of greenhouse gas [70].
The solution to exploit the LNG cold potential and thus extend the temperature range in Brayton closed-loop recuperation systems designed to utilize the secondary heat generated by modern marine engines is seemingly more rational because the exhaust gas temperature of these engines, as well as for duel-fuel engines, is low. For this purpose, the References [71][72][73] reviewed the structural features of marine propulsion NG fuel supply systems; the structural solutions of the engines of the gas carriers MAN Diesel 2 Turbo and Wärtsilä were reviewed [39,40].
Unlike the classic steam turbine propulsion complex of operating LNG carriers, 85% of newly built fleets of LNG carriers will be equipped with a propulsion system with gas diesel engines. This is based on minimizing cargo vapor losses (boil-off gas (BOG)) and reducing air emissions from ships. About 30% of the LNG carrier fleet is already equipped with the two-stroke dual-fuel engines MAN B&W ME-GI [71]. The efficiency of gaseous fuel supply to the engine system depends on engine load regime, freight conditions, BOG intensiveness, gaseous fuel composition, etc. [71]. Therefore, a constant supply of gaseous fuel to the engine can be ensured by using combined fuel supply systems, which consist of the LNG tank BOG high-pressure (150-300 bar) compressor circuit and the LNG compression to a 300-bar circuit, followed by regasification and heating to 45 • C [74].
For BOG compression, a 5-stage and 6-cylinder reciprocating machine, LABY ® -GI, with the intercooler and bypass, manufactured by Burckhardt Compression, is widely used. The compressor is equipped with a sophisticated automatic control system to provide the necessary gas pressure in different floating modes, with the additional gas liquefaction circuit developed by Hamworthy Gas Systems [74].
If the entire BOG is not consumed as a fuel, the remainder, after two steps of compression (stage 1-2), is directed to a heat exchanger where it is liquefied with nitrogen and returned to the LNG tanks after separation.
If the ship's power plant is powered by liquid fuel only, the entire BOG is liquefied. A similar system with BOG re-liquefaction has been developed by Wärtsilä Oil & Gas Systems and is installed on gas tankers [72].
Thus, dual-fuel vessels can be equipped with recuperative engine systems consisting of cryogenic heat exchangers (BOG liquefaction and LNG regasification), enabling the cold potential of LNG to be exploited in secondary heat energy recuperation systems.
The analysis of various research results proved that the Brayton closed gas turbine is the most promising secondary main engine heat recovery system. Its relatively lightweight and overall size make it easy to integrate into a ship's power plant. By analogy with modern gas fuel supply systems, the rational working agent to use in this system is liquid nitrogen.
The partners of the Shanghai 2017 'Science and Technology Innovation Action Plan' Intergovernmental International Science and Technology Cooperation Project 'The Study of Small LNG Ship and Cascade Use of Cold Energy' developed an efficient gas-to-gas heat exchanger to achieve a nitrogen closed-cycle of cold energy exchange.
In the system presented in this paper, the exhaust gas heat recovery of the main engine is realized in a closed Brayton cycle. Gaseous nitrogen circulating in the loop is heated using acetyl glycol, which in turn is heated by the exhaust gas of the engine. Two regenerator heat exchangers are provided in the cycle to increase energy efficiency. The N 2 -LNG heat exchanger serves to decrease the temperature of N 2 entering the compressor; in the N 2 -N 2 heat exchanger, the compressed nitrogen is heated by the leaving the turbine nitrogen heat and conveyed to the engine's heat exchanger. The electricity in the generator is generated by the difference in power output between the turbine and compressor.
The goals of this study were focused on two main aspects of this technology research: • To create a mathematical model to study the cogeneration cycle technology of an LNG ship power plant and to conduct their adaptation in various calculations. • From the mathematical model, to create an electronic simulator and perform parametric analysis of the regeneration cycle to increase the energy efficiency indicators by rationally selecting the components of the cogeneration system (turbine, compressor, and heat exchangers).

Methodological Aspects
The boundary conditions for the analytical decisions comprised the following key basic regulations: • Energy transformation processes in the system's equipment are reversible (aerodynamic pressure losses are zero).

•
Turbine and compressor analytical connection of parameters η C.ad = f (G, π C ) and π C = f (G, n TC ) are not desirable: the parameters T iC = T iT are independent of η C.ad and η T.ad .

•
The LNG cold potential is external and unlimited (the LNG flow exceeds the demand of supply energy for the engine).
Additionally, two alternative scenarios are discussed: heat transfer in heat exchangers is determined by the thermal efficiency η HE , or LNG flow exclusively used as a fuel for the engine. The analytical description of the system's function consists of three equations (1, 2, and 3) of axial machines and the theory of heat machine equipment [75,76]: where G d and G d are the gas flow rates through the turbine and compressor, respectively, C p and C p are the isobaric specific heats of gas flow through the turbine and compressor, respectively, and K d is the gas adiabatic indicator; the pressure P i and temperature T i marking and location are based on the scheme (Figure 1). The algorithm of the analytical solutions is shown in Figure 2.  The thermal energy, generated by the exhaust gases of the marine engine (DE), is transferred to the cogeneration cycle through a heat exchanger (HE3). After the heat exchanger HE3, nitrogen circulating in the circuit with a high temperature 3 enters the turbine T, where part of its internal energy during expansion is transformed into mechanical energy of the turbine.
After passing through the HE1 heat exchanger, the nitrogen is cooled by LNG in the HE2 heat exchanger and enters the compressor inlet (C) at a low temperature ( 1 ), where its temperature rises to as a result of compression. In the cogeneration heat exchanger HE1, nitrogen is additionally  The thermal energy, generated by the exhaust gases of the marine engine (DE), is transferred to the cogeneration cycle through a heat exchanger (HE3). After the heat exchanger HE3, nitrogen circulating in the circuit with a high temperature 3 enters the turbine T, where part of its internal energy during expansion is transformed into mechanical energy of the turbine.
After passing through the HE1 heat exchanger, the nitrogen is cooled by LNG in the HE2 heat exchanger and enters the compressor inlet (C) at a low temperature ( 1 ), where its temperature rises to 2 as a result of compression. In the cogeneration heat exchanger HE1, nitrogen is additionally The thermal energy, generated by the exhaust gases of the marine engine (DE), is transferred to the cogeneration cycle through a heat exchanger (HE3). After the heat exchanger HE3, nitrogen circulating in the circuit with a high temperature T 3 enters the turbine T, where part of its internal energy during expansion is transformed into mechanical energy of the turbine.
After passing through the HE1 heat exchanger, the nitrogen is cooled by LNG in the HE2 heat exchanger and enters the compressor inlet (C) at a low temperature (T 1 ), where its temperature rises to T 2 as a result of compression. In the cogeneration heat exchanger HE1, nitrogen is additionally heated by the flow of nitrogen leaving the turbine, and at temperature T 5 enters the heat exchanger of the engine exhaust gases HE3. The cogeneration heat exchanger improves the efficiency of the cycle by expanding its operating temperature range (T 3 − T 1 ).
The algorithm is primarily linear; however, one of the key aspects that affect the efficiency is solved using the iteration method: • After entering the technical and operational parameters (π C , G N2 , etc.), the nitrogen temperatures are calculated at the characteristic angles of the co-generation cycle (T 1 ; T 2 ; T 3 ; T 4 ; T 5 ). It is assumed that there are no hydrodynamic losses of nitrogen flow in the cycle heat exchangers channels. As a result, P 1 = P 4 = P 6 and P 2 = P 3 = P 5 . Thus, the algorithm realizes the equality of the degree of pressure increase in the compressor and the degree of pressure decrease in the turbine π C = π T .

•
The outgoing temperature T 4 from the turbine, including thermal efficiency η HE1 of the heat exchanger, directly affect T 5 , which in turn affects the temperatures T 3 and T 4 after increase in heat exchanger HE3.

•
The nitrogen preheat in heat exchanger HE3 is limited by the engine's exhaust gas heat transfer potential to nitrogen. Depending on the parameters used in the system, the nitrogen temperature T 3 is adjusted according to condition Q eg ≥ Q N 2 HE3 . • The temperature T 1 for an alternative supply for diesel engine use was determined based on the potential of the difference between saturated and overheated LNG steam phases.
The effect of the cogeneration cycle, converting excess energy to electric energy, is calculated using the analytical equation For further flexible development, the algorithm of the cogeneration cycle was realized in the form of an electronic simulator using the software LabView (National Instruments) (Figure 3). A medium-speed six-cylinder four-stroke diesel engine (bore/stroke 30/50) was used to check and adapt the algorithm and its simulator to real maritime conditions. Table 1 represents the indicators of engine operating in propulsion characteristics.

Research Results and Discussion
In the first step of this research, a task was set and accomplished: based on the cycle calculations (see Section 2), to create a nomogram to determine the energy parameters of the cogeneration cycle. Graphical evaluations are convenient for creating technological parameters for actual cycles. The foundation of calculations was the first equation of the turbo unit that consisted of the turbine and compressor [75]: The calculation conducted with the relative form of the parameter ∆ , i.e., one unit (1.0) is the enthalpy ( 2 • 2 ) of the circulating nitrogen: where for simplicity, 2 = 1.06 / • was used for the range of 2 = 1.04 ÷ 1.08 / • , which is characteristic of the temperature range of 200-700 ℃; the energy pulse factor = 1 because the system was isobaric.
As a result, a graphical form (Figure 4) implementing the determination of the parameter values is shown in Figure 4.
The use of nomograms does not require an analytical evaluation of the effect of different impacting factors on the relative parameter ∆ because, with the help of a nomogram, it is realized for every evaluated variant.
Note that at this research stage, all key factors were independent, except for the equality = that determined the structure of the closed circulating nitrogen cycle. In this nomogram, technological (  . ,  . ,  ) and energy ( , 2 ) parameters of an operation were set and regarded as independent. In practice, most of them are interrelated:  Remark: P e is the effective power, P C is the charge air pressure, T eg is the exhaust gas temperature, G eg is the exhaust gas flow rate, and λ is the excess air ratio.
The engine work mode was set using simulator tools that imitated the engine control system (Figure 3). The parameters T eg , G eg , were the characteristics of this mode and used in cogeneration cycle calculations. A wide range of tasks, related to cogeneration cycle efficiency and parameter values harmonization, were solved using the simulator: A part of this research is represented in Section 3 below.

Research Results and Discussion
In the first step of this research, a task was set and accomplished: based on the cycle calculations (see Section 2), to create a nomogram to determine the energy parameters of the cogeneration cycle. Graphical evaluations are convenient for creating technological parameters for actual cycles. The foundation of calculations was the first equation of the turbo unit that consisted of the turbine and compressor [75]: The calculation conducted with the relative form of the parameter ∆N TK , i.e., one unit (1.0) is the enthalpy G N 2 ·c pN 2 n of the circulating nitrogen: where for simplicity, c pN 2 = 1.06 kJ/kg·K was used for the range of c pN 2 = 1.04 ÷ 1.08 kJ/kg·K, which is characteristic of the temperature range of 200-700 • C the energy pulse factor β = 1 because the system was isobaric. As a result, a graphical form (Figure 4) implementing the determination of the parameter ∆N TC G N 2 ·c pN 2 well revealed the structure of analytical expression; this was the second part of Equation (6)  The solutions for the gas engine that used natural gas 4 are presented below.

Limitations of 3 Compared with
The parameter was independent of the 2 cycle and determined by the load regime of a particular engine. The condition 3 = was for an ideal scenario: ƞ 3 = 1.0, and maintaining the Note that at this research stage, all key factors were independent, except for the equality π C = π T that determined the structure of the closed circulating nitrogen cycle. In this nomogram, technological (η C.ad , η T.ad , η HE ) and energy T eg , G N 2 parameters of an operation were set and regarded as independent. In practice, most of them are interrelated: As a result, a graphical form (Figure 4) implementing the determination of the parameter   Figure 4.
The use of nomograms does not require an analytical evaluation of the effect of different impacting factors on the relative parameter ∆ because, with the help of a nomogram, it is realized for every evaluated variant.
Note that at this research stage, all key factors were independent, except for the equality = that determined the structure of the closed circulating nitrogen cycle. In this nomogram, technological (  . ,  . ,  ) and energy ( , 2 ) parameters of an operation were set and regarded as independent. In practice, most of them are interrelated: ✓ 3 depends on exhaust gas temperature , exhaust gas flow , circulating 2 flow, , , and efficiency factors of a heat exchanger ( ); T 3 depends on exhaust gas temperature T eg , exhaust gas flow G d , circulating N 2 flow, π k , π T , and efficiency factors of a heat exchanger (η HE ); which is characteristic of the temperature range of 200-700 ℃; the energy pulse factor = 1 because the system was isobaric.
As a result, a graphical form (Figure 4) implementing the determination of the parameter ∆ 2 • 2 well revealed the structure of analytical expression; this was the second part of Equation (6) values is shown in Figure 4.
The use of nomograms does not require an analytical evaluation of the effect of different impacting factors on the relative parameter ∆ because, with the help of a nomogram, it is realized for every evaluated variant.
Note that at this research stage, all key factors were independent, except for the equality = that determined the structure of the closed circulating nitrogen cycle. In this nomogram, technological (  . ,  . ,  ) and energy ( , 2 ) parameters of an operation were set and regarded as independent. In practice, most of them are interrelated: ✓ 3 depends on exhaust gas temperature , exhaust gas flow , circulating 2 flow, , , and efficiency factors of a heat exchanger ( ); T 1 depends on the heat transfer balance in the heat exchanger HE2 (which are the parameters G CH 4 , G N 2 , and T 6 ) Therefore, in solving the practical tasks of the nomogram, it was rational to describe the limitations of the existing relations between the parameters.
The solutions for the gas engine that used natural gas G CH 4 are presented below.

Limitations of T 3 Compared with T eg
The parameter T eg was independent of the N 2 cycle and determined by the load regime of a particular engine. The condition T 3 = T eg was for an ideal scenario: η HE3 = 1.0, and maintaining the following conditions: the heat is removed from the exhaust gas up to the 'dew point'~132 • C (normal conditions), which is sufficient to warm up the nitrogen G N 2 up to temperature T eg . Heat balance equation is expressed as: G N 2 ·c pN 2 ·∆T ≤ G eg ·c p.eg T eg − 405 (7) or G N 2 ·c pN 2 ·∆T ≤ G eg ·c p.eg T eg − 405 (8) where c p.eg is the specific heat of exhaust gases. If the inequality is true, the parameter ∆T is an increase in N 2 temperature from T 5 due to heat exchange with engine exhaust gas up to T eg . In the opposite scenario, T 3 will be lower than T g by the value ∆T f .
Engine technical documentation consists of fuel consumption data G f [kg/h]; at certain load mode, a coefficient of excess air α and T eg are known. Subsequently, G eg = (αL 0 + 1)G f and Equation (7) is represented by the following equation: At different α and temperatures, c pN 2 and c p.eg are determined by interpolation from standard references. Table 2 represents c p.eg values for approximate evaluation when the engine operates on diesel fuel and natural gas (CH 4 ). For a preliminary evaluation (modes comparable to nominal power, because α intensively increasing at low load modes) λ = 2.0 value is accepted; in a result G eg = 35G f . If the G N 2 value is lower than G eg c p.eg (Teg−405) ∆T·c pN 2 , temperature T 3 will attain T eg . In an opposite scenario, reduce temperature T 3 in comparison to T g : where For the calculation of temperatures ∆T f or T 3 values, the concretization of ∆T is necessary using (T eg − T 4 ). For T 3 to be equal to T eg (at η HE3 = 1, 0), the inequality must be: However, according to an accepted condition η HE3 = 1, 0, T 5 = T 4 (i.e., N 2 attains the temperature after expansion in a turbine).
Finally, for the condition execution: Equality of T 3 = T g is ensured if a left side of Equation (12) is higher than or equal to the right side. Equation (12) enables the evaluation of the G N 2 boundary values for different π T values (or π C based on the condition π C = π T ) for certain operation modes of specific engines: Example of Analytical Calculation Application for Medium-Speed Six-Cylinder Four-Stroke Engine The evaluation of a six-cylinder four-stroke engine (bore/stroke 30/50) was conducted in nominal and part load modes of propulsion characteristics. Nominal mode: P e =3500 kW; G eg = 23,950 kg/h = 6.65 kg/s; T eg = 775 K; c p.eg ∼ 1.15 kJ/kg K; c pN 2 ∼ 1.06 kJ/kg K; η T.ad = 0.8. Part load mode: P e = 400 kW; G eg = 7500 kg/h = 2.08 kg/s; T eg = 575 • K; η T.ad = 0.8. The results of boundary nitrogen flow boundary G N 2 bound are shown in Table 3. Entering the received data into the cogeneration cycle nomogram enables a rational field for combinations of parameters for the selected engine to be described.

Determination of Temperature T 1
If a simulated engine uses natural gas (CH 4 ) from a supply line, a cold potential of CH 4 assessed as unlimited, i.e., T 1 extends to −161 • C (LNG liquid phase temperature). A scenario is represented below where CH 4 flow is used only as a fuel to feed an engine, i.e., all flow G N 2 is used for engine feed.
The warming of N 2 from T 6 to T 1 (exiting the heat exchanger HE2) is revealed by Equation (15): (at P = 0.5-0.6 MPa), where c ev is the evaporation heat of liquid CH 4 , (293 − 112) K − CH 4 vapor warming up to +20 • C for delivery to an engine.
For a more accurate calculation: where i i is the enthalpy of CH 4 superheated steam at P = 0.5-0.6 MPa, and i 2 is the enthalpy of liquid CH 4 . According to Reference [77] , nitrogen will be chilled in heat exchanger HE2 up to 112 K. On the assumption that the efficiency of the heat exchanger HE1 is η HE1 = 1.0, a flow of nitrogen will be cooled from temperature T 4 (after expansion in the turbine) up to the nitrogen temperature T 2 (exiting from compressor). That means that the N 2 flow in heat exchanger HE2 is cooled from T 2 to T 1 .
As a result, Equation (16) is transformed into Equation (17): For a boundary condition to attain T 1 = 112 K (CH 4 liquid phase temperature), Equation (18) is used: Finally, As an example in the nomogram (Figure 4), a nitrogen temperature T 1 = 112 K is ensured by drawing a limiting curve for a six-cylinder four-stroke diesel engine (bore/stroke 30/50) according G f /G N 2 ) and π C parameters. A specific π C exists (assumption η C.ad = idem) for each G f /G N ratio when T 1 = 112 K is ensured. An increase in temperature T 1 from -112 K by ∆T 1 is determined according to Equation (14): where indexes '1' and '2' are assigned to the 'limit' variant to ensure T 1 = 112 K and T 1 increases above 112 K, respectively. After the solution of the equation system: ∆T 1 changes from limit value 112 K if any of impacting factors G f , G N 2 , π C, and possibly η C.ad are changed. Based on Equation (10), boundary values of G N 2 and G N 2 bound. parameters are determined (Table 4) for a six-cylinder four-stroke diesel engine (bore/stroke 30/50) operating in the nominal load (P e = 3500 kW) and part load (P e = 400 kW) modes. The obtained results ( Figure 4) justify the significant effect of this aspect on the functioning of the cogeneration cycle to ensure high energy consumption rates. Analytical evaluations and nomograms conducted for technology operating in an ideal condition when the efficiency of all heat exchangers is EC = 1.0. The effect of deviations from the ideal scenarios of heat exchangers useful performance factor on electricity generation assumed by decreasing HE from 1.0 to 0.9. The results of ∆N iC change are shown in Figure 4. The evaluation was performed with two limit options of T eg = 200 and 500 • C, and π C = π T changed in range of 2-5 units. Figure 4 shows the effect of an increase in the turbo compressor and turbine efficiency on the final effect of the system's ∆N CT functioning (for a base version accepted condition η C.ad = η T.ad = 0.8). Furthermore, the cogeneration cycle nomogram aids in determining the interconnection of parameters ensuring the efficiency of a turbine and compressor operation to produce electricity.
Thus, point 1, corresponding to selected π T (e.g., 3), in the first quadrant of the nomogram is determined for a selected temperature T 1 (e.g., 273 • C). The projection of this point on the axis N C forms the corresponding point 1. The specific compressor output N C c PN 2 ·G N 2 (point 1, corresponding to 110 K) is determined for a selected compressor efficiency η C.ad value (e.g., 0.9). At a boundary condition N C = N T , a corresponding result on the N T axis must be obtained from projections of the T eg − π T relationship. Subsequently, the relationship between the parameters T eg and π C is defined as the intersection point 2 in the second quadrant using the corresponding values of the projections (e.g., selecting η T.ad =0.9; point 2).

Combinations of Optimization
To ensure an effective operation of the system during electricity production (i.e., N T > N C ), the values of parameters T eg and π T are necessary from the second quadrant field (at T eg and π T values higher than at point 2). Changing π T , η C.ad , and η T.ad , preliminary conceptual decisions (equipment functioning and construction parameters) of the system are concretized with the aid of the nomogram ( Figure 5) for aggregate assembly.
The task could be solved in a more general context. The 2nd quadrant of the nomogram defines the field of positive values of ∆N TC , where the power of a turbine N T is higher than the compressor power N C for any ratio π T − T eg − G N 2 .
In the 1st quadrant, 4-5 dots are fixed on the selected line T 1 at different values of π C from the predetermined range. The projections determined from selected π C to the selected value of C.ad line and further to N C / c PN ·G N 2 and axes. The points obtained from the axis N T / c PN 2 ·G N 2 designed on the selected 3rd quadrant line η T.ad and further into the 2nd quadrant, where the values cross the corresponding π C values. The obtained intersection points are joined by a line, which separates the field N T ≥ N C in the 2nd quadrant.

Factors Defining the Energy Efficiency of the System
The factors with the highest effect on the energy efficiency of the cogeneration system were identified by performing various simulations, the results of which are shown in Figures 5 and 6  Ngen, kN πc=πt N2 temperature T1=112 K; N2 temperature T1=150 K; N2 temperature T1=190 K; N2 temperature T1=230 K; Figure 6. Effect of technological system parameter π C = π T on N gen. . Figure 5a shows the data of the engine exhaust temperature T eg for a generated electricity (N gen ) at different π C = π T values. Increasing T eg by 220 K approximately doubled the energy efficiency of the system operation: beginning from π C > 2.0, the parameter N gen increased by 70-85% or 3.2-3.9% for 10 K of T eg . The approximation dependencies N gen = f π C , T eg var exhibited a non-linear nature: the increase in pressure of a turbocharger unit above π C (π T ) > 3.0 decreased the effect on N gen . For example, an increase of π C (π T ) for every 0.5 units up to 3.0 determined an increase in N gen by 7-10 kWand only by 4-5 kWin range π C (π T ) = 3.0-5.0 units. In contrast, the value of N gen was directly proportional to the flow of N 2 circulating in the closed loop (G N2 = 1.0 kg/s in the variation calculation). Thus, the realization of a cogenerating cycle is rational when π C (π T ) ≤ 3.0 at a maximal realized G N2 . Furthermore, a realization of a high pressure in a TK is known to be related to the difficulty in ensuring high efficiency in a wide range of TK operating modes [78,79]. Figure 5b shows the impact evaluation results of adiabatic efficiency of the compressor η C.ad and turbine η T.ad . The effect of η C.ad increasing from 0.8 to 0.9 on N gen did not exceed 6 kW(regardless of T eg ). The effect of a turbine's η T.ad in the same range was significantly higher: an increase of N gen was from 9 kWat T eg = 555 K and up to 15 kWat T eg = 775 K or~13% (Figure 5c). Thus, the primary accents of the T C model selection should be given to select the following combination of parameters: high productivity T eg and medium pressure π C (π T )~3.0 units at a maximal efficiency η T.ad of the turbine. Figure 5d shows the results of the effect of the efficiency of heat exchangers (η HE ) on the evaluation. The results ranged from an ideal scenario (η HE =1.0) for thermal efficiency of heat exchangers up to an actual scenario. In practice, the value of~0.9 does not have a significant effect on N gen in the range of π C (π T ) ≤ 3.0. A decrease in the energy efficiency of the system one by one in η HE system results only in 2-3 kW.

Effect of the Efficiency of Heat Exchangers
A decrease from 1.0 to 0.9 of the general η HE in existing heat exchange system results in a decrease in system efficiency by 12 kWat π C (π T ) = 3.0, which is~11%, and limiting values of π C (π T ) = 5.0 efficiency results in a decrease of 21 kWor~15%. Generally, the effect on N gen is related to an efficiency change of LNG heat exchanger N 2 (Figure 5d). When the η HE of LNG-N 2 heat exchanger decreases, the temperature T 1 of entering nitrogen increases by 14-17 K from entry-level 112 K.
The increase in temperature T eg is determined by the change in the efficiency of the LNG-N 2 heat exchanger. Earlier variation calculations are based on an 'external' unlimited LNG cold potential which in itself ensures the value T eg ≈ 112 K η HE = 1.0 of exiting flow from the exchanger. System operation variants are related to an LNG 'internal' use, i.e., with direct use only in an engine, and when LNG's cold potential is lacking; this results in an intensive increase in temperature T 1 , particularly at high N 2 flow rates. Figure 6 shows the results of the evaluation of the effect of temperature T 1 from the LNG-N 2 heat exchanges on N gen . The evaluation was conducted for a T 1 range of 125-230 K at actual values of HE = 0.9 and maximal temperature T eg = 775 K. The effect of change in T 1 on N gen was comparable to linear dependency. Each increment of T 1 by 10 K resulted in a lower value of N gen bỹ 5 kWat π C (π T ) = 3.0. At higher TK pressure π C (π T ) = 5.0, the effect of change in T 1 increased up to 8 kW. Hence, when T 1 increased, the effect of pressure change decreased. When T 1 attained a range of 220-230 K, N gen became independent from π C (π T ). Thus, an effective application of the cogeneration system in a ship can be ensured by making an alignment between G N 2 flow rate and parameters of the LNG engine's fuel consumption and energy efficiency of the exhaust gas.

Conclusions
Based on the theoretical foundations of turbo machines and heat exchangers, a mathematical model of the closed-loop nitrogen cogeneration cycle and its electronic simulation (using LabView [National Instruments]) for a ship power plant were created.
For higher efficiency, for a practical selection of rational parameters of a power plant cogeneration cycle, a nomogram of an interconnection of cogeneration cycle parameters was created and experimental data were adapted. With the aid of the nomogram, in the function of selected cycle technological parameters (turbo unit pressure and efficiency rate, exhaust gas temperature, flow rate of circulating nitrogen, and efficiency rate of heat exchangers), the amount N gen of electricity produced in the cogeneration cycle was determined.
1. Using a variable simulation the key parameters determining the energy efficiency was investigated for a six-cylinder four-stroke engine fueled by LNG, and technological principles of energy efficiency (N gen ) improvement were formulated:

•
In the judgment of the authors, the main advantage for the practical use of the completed development, in contrast to a number of similar ones, is an open algorithm that provides for the expansion of the model in accordance with the technological features of the heat recovery systems under study, as well as operating mainly with thermodynamic parameters without concrete definition of the constructive nodes. The latter, including on the basis of the created electric simulator, expands the possibilities of a variant search and justification of rational ways to increase the indicators of the energy efficiency of cogeneration systems.

•
With an unlimited 'external' cold potential of LNG (for example, an LNG carrier power plant), the energy efficiency of a cogeneration cycle determined by the exhaust gas temperature T eg of the power plant, turbo unit pressure π C (π T ) (for practical implementation in a rational range up to 3.0-3.5 units), and flow rate G N 2 of circulating nitrogen is directly proportional to the implemented N gen . A turbo unit with high output and average pressure π C (π T )~3.0-3.5 with a high adiabatic efficiency turbine is a more rational selection for a cogeneration cycle assembling by technological units. The effect of selected thermal efficiency rates of the heat exchangers (in range of change 1-0.9) on N gen does not exceed~10%.

•
With a limited potential of LNG cold flow and LNG 'internal' use only as fuel for ship engines, a harmonization between the flow rate of the circulating nitrogen and engine load modes is necessary. The temperature T 1 of the exiting nitrogen from an LNG-N 2 heat exchanger has a significant effect on the energy efficiency N gen . For each temperature increase of 10 K, the average decrease of N gen value is 5-8 kW or 5-6%.
2. The created mathematical model of an electronic simulator, using accepted assumptions, addressed the alternative search of rational and constructive solutions for the cogeneration cycle. A further aspect of mathematical model improvement provided for optimization of operation parameters:

•
To implement a simulation of the TK compressor and turbine parameters based on the real characteristics of aggregates G N 2 = (π C , η C.ad , n TC ); G N 2 = (π T , T 3 , n TC ).

•
To apply a mathematical model to determine the operating process parameters of an LNG engine and to ensure its operation together with the cogeneration cycle model.

Conflicts of Interest:
The authors declare no conflict of interest. Heat transfer in heat exchanger number 3 T eg

Nomenclature
Exhaust gases temperature C PN 2 Nitrogen specific isobaric heat C Peg Exhaust gases specific isobaric heat β Energy pulse factor G N 2 Nitrogen flow G eg Exhaust gas flow G f Fuel consumption π C Degree of pressure increase of the compressor π T Degree of pressure decrease of the compressor k Adiabatic coefficient G N 2 bound Boundary nitrogen flow N C Compressor power N T Turbine power n Engine maximum operating speed A Coefficient of nitrogen flow impulses effect λ excess air ratio lT turbine operation lC compressor operation nTC compressor speed ψ outflow leakage function