Next Article in Journal
Cooperation between Sea Ports and Carriers in the Logistics Chain
Next Article in Special Issue
Marine Power Systems
Previous Article in Journal
Formation Control of Unmanned Vessels with Saturation Constraint and Extended State Observation
Previous Article in Special Issue
Research on Improving the Working Efficiency of Hydraulic Jet Submarine Cable Laying Machine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy, Economic and Environmental Effects of the Marine Diesel Engine Trigeneration Energy Systems

1
Department of Maritime Sciences, University of Zadar, Mihovila Pavlinovića 1, 23000 Zadar, Croatia
2
Faculty of Engineering, University of Rijeka, Vukovarska 58, 51000 Rijeka, Croatia
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(7), 773; https://doi.org/10.3390/jmse9070773
Submission received: 15 June 2021 / Revised: 11 July 2021 / Accepted: 13 July 2021 / Published: 16 July 2021
(This article belongs to the Special Issue Marine Power Systems)

Abstract

:
The paper discusses the possibility of applying the trigeneration energy concept (cogeneration + absorption cooling) on diesel-powered refrigerated ships, based on systematic analyses of variable energy loads during the estimated life of the ship on a predefined navigation route. From a methodological point of view, mathematical modeling of predictable energy interactions of a ship with a realistic environment yields corresponding models of simultaneously occurring energy loads (propulsion, electrical and thermal), as well as the preferred trigenerational thermal effect (cooling and heating). Special emphasis is placed on the assessment of the upcoming total heat loads (refrigeration and heating) in live cargo air conditioning systems (unfrozen fruits and vegetables) as in ship accommodations. The obtained results indicate beneficiary energy, economic and environmental effects of the application of diesel engine trigeneration systems on ships intended for cargo transport whose storage temperatures range from −25 to 15 °C. Further analysis of trigeneration system application to the passenger ship air conditioning system indicates even greater achievable savings.

1. Introduction

Marine diesel engine trigeneration energy systems (MDETES) represent coupling diesel engine cogeneration systems with absorption cooling systems, which allow the use of a co-generation effect to balance the occurring overall heating loads on marine motor ships intended for the transport of moderately cooled cargoes, [1]. Modern marine four stroke engines’ energy efficiency rate is in the range of about 49% to 50% at its nominal load, [2,3]. Two stroke engines’ energy efficiency is slightly higher, at the gas mode about 53% [4,5]; however, on the diesel mode, energy efficiency is lower at about 52%. The calculation for the two-stroke engine is based on the LHV of the MDO 42.7 MJ/kg, and LHV of a typical LNG 48.0 MJ/kg, [6,7]. As the marine engines cannot convert all heat energy to power, the remaining exhaust heat energy may be utilized for the boiling process of the absorption cooling unit. The typical process utilization overview was given in [8,9]. The exhaust gas heat energy of marine engines [10,11] is not the only source of heat energy; waste cooling energy of the marine engines may be utilized for that purpose as well. The waste heat capacity potential of the typical stationery and marine engines is analyzed in [12,13], where it is concluded that this part of the energy also has potential for trigeneration purposes. The improvement of the waste heat recovery from marine diesel engines with the best fulfilment of the vessel needs in terms of mechanical, electric and thermal energies is analyzed in [14,15], where various solutions are proposed, with combined diesel engine, steam and gas turbines recovering part of the thermal energy of the diesel engine exhaust gas. The similar approach of feedwater regeneration to the boiler from the marine engines was presented in [16], which also concluded that waste heat will beneficially contribute to the efficiency of waste heat recovery. In marine diesel engine cogeneration energy systems (MDECES), the propulsion of the main diesel engine (MDE) directly or indirectly drives the propeller (P), and the shaft electric generator (SEG), balancing both occurring propulsion and electrical load during the voyage, while the occurring heating load is balanced by the waste heat recovery steam generator (WHRSG). During both maneuver and standstill, the ship’s occurring electric load is balanced by the diesel aggregates (DA), and the resulting heating load of the ship is balanced mainly by the fired steam generator (FSG) due to the insufficient availability of waste heat contained in the exhaust gases of the DA. The overview of the energy balance for the cruise ship and the chemical tanker is given in [17,18]. As the vessel is sailing, it is changing the local microclimate’s environmental condition. The engine is not always operating in the steady condition but is changing load, which affects the absorption process. In order to reduce the negative effect on the absorption process, [19] carried out where proposed a solution to the amortization of the heat load changes with the small diameter tube bundle heat exchangers with large specific surface area. The changes of the local microclimate environment should not always affect the absorption process in the negative direction, according to the study under [20]; off design performance of the LiBr-H2O, particularly in the lower generator or higher condenser temperature conditions, generates both higher COP and exergy efficiency. However, in this particular study the absorption process in the complex rolling and pitching environment of the ship is not considered.
Due to economic considerations, the ship’s diesel engines are supplied with cheaper, heavy fuel (HFO, IFO380) [21,22]; therefore, the combustion gases must not substantially cool down because of the presence of sulfur oxides in them, which consequently limits the cogeneration plant heating effect. During navigation, unsteadiness of the cogeneration thermal effect is emphasized due to its dependence on emerging propulsion load (which is reflected through fluctuation of power and engine speed) and on the temperature of the surrounding air, which was discussed in [23,24,25]. On the other hand, the heating load of the ship has been determined by its purpose, respectively prescribed microclimate of the contained commodities compartments, and of the navigation route, so that it is continuously changing according to which climate zone the ship is in [26]. Apart from the need to balance the heating load, it is necessary to create an appropriate cooling effect for balancing the cooling load [27], which is in almost all modern ships balanced by the compression refrigeration plant (CRP); this significantly increases the electric load, while the heating effect of the installed cogeneration plant is being considerably unexploited [28]. By the use of an absorption refrigeration unit (ARU) with thermodynamic mixtures of water-lithium bromide (H2O-LiBr) [29,30], and ammonia water (NH3-H2O) [31,32], and with negligible electrical load, the utilization of the cogeneration thermal output for a production appropriate cooling effect is allowed. The attainable lower temperature limit in the evaporator ammonium ARU (AARU) restricts its application on both accommodation air conditioning systems and commercial cargo that is being transported and stored in a moderate temperature range (−25 to 15 °C). On the other hand, an attainable lower temperature limit in the evaporator lithium bromide ARU (LBARU) restricts its application on the accommodation air conditioning systems. By coupling the DECES and ARU, the (MDETES) is created, which enables cost-effective, energy-efficient, and environmentally friendly balancing of the occurring heating loads of both air conditioning systems and moderate commercial refrigerated cargo by the use of the cogeneration heating effect. Reducing the electrical load automatically reduces the mechanical loads of the MDE during navigation, as well as of the DA during maneuvering and the ship’s standstill in the terminals, which implies substantial fuel savings and consequently lower operating costs with the reduction of the majority of ecologically harmful effects (primarily CO2 and SO2 emissions) according to [33,34].
This paper continues with efforts in the rational usage of the waste heat from the marine diesel engines. It presents a comprehensive mathematic model of the proposed trigeneration system for reefer use. The concept interrelates the trajectory of a vessel and the wave conditions with a detailed described quasi-static model of the energetic, economic and environmental performance of a marine trigeneration unit installed onboard.

2. Defining Technically Possible Application Fields ARU’s

Technically, the obtainable upper temperature of the absorption chiller cooker is ϑcc ~ 170 °C (which is determined by the parameters of steam from the cogeneration plant), and the cooling sea temperature ϑsd ~ 32 °C. These two parameters determine the applicability of the absorption cooling with ammonia mixture on those cooling systems where the temperature of ammonia primary refrigerator in evaporator does not fall below ϑse ~ −35 °C and considering the usual arrangement of the cooling system with CaCl2 brine as secondary refrigerant, which corresponds to the storage temperature limited at the ϑca ~ −25 °C. Such achievable lower temperature storage facilitates the transporting of a wide range of palletized (packed bulk) and liquid cargo, where the majority of commercial bulk palletized cargos comprise foodstuffs, transporting either frozen (dead) or unfrozen (live). The liquid cargos include nutritional liquids such as various fruit juices (e.g., orange juice), and industrial organic liquids such as some alkanes (paraffin), alkenes (olefins), alkyne (acetylene), alkanes (aldehydes), etc.
Among the frozen commodities, meat and fish dominate, and are generally stored at equal low temperature ϑFM ~ −18 °C, although the meat is often carried unfrozen in vacuum packs in the temperature range ϑMV ~ (−2 to 10) °C depending on the type of meat and the duration of storage [35].
Micro-climate storage of the unfrozen (living) perishable products, among which prevail bananas, citruses, deciduous fruits and vegetables, is characterized by a relatively high relative humidity (φlp ~ 0.8 to 0.95), moderately low temperature ϑlp ~ −2 to 15 °C and a relatively small proportion of fresh air (a small refreshment g0 1.5 ~ 2.5%). Nutritional liquids are refrigerated and held in separate tanks in an inert gas atmosphere, with the temperature depending on their respective sugar contents, for example, orange juice is held at ϑOJ ~ −7 °C, [36].
Typical industrial organic liquids are stored at appropriate saturation temperature ϑs (pa) at normal atmospheric pressure pa, and are as follows: paraffins (C4H10 n-butane −0.5 °C, isobutane −11.7 °C, cyclobutane 10 °C), olefins (C4H8 butene −6.6 °C, C4H6 cyclobutene 2.2 °C), acetylenes (propyne −23.2 °C, butyne 7.7 °C and aldehydes) methanal −19 °C, ethanal 20.2 °C, etc., [37].
With cooling sea water of ϑsd ~ 32 °C and with a heating sink of ϑ ˜ s ~ 120 °C (low-pressure saturation steam p ˜ s ~ 2 bar), attainable temperature from mixture extracted water as primary refrigerator in the evaporator of lithium bromide ARU amounts to ϑse ~ 5 °C; hence, by the usual arrangement of the cooling system with water as secondary refrigerator in ships’ accommodation air-conditioned systems, the conditioning air can be cooled down to a temperature of ϑac ~ 13 °C, and consequently the application of this ARU to create the desired level of comfort is particularly suitable in air conditioning systems of passenger ships. Generally, there is a possibility that the same ship simultaneously carries different cargo stored in separate compartments with the appropriate prescribed microclimate. Taking into account ships’ accommodation, balancing of the appearing cooling load of such a ship can be achieved either with two apart single-stage AARUs, one for “dead” and another for “live” products (which is expensive), or with a double-stage ARU for the both “dead” and “live” products, and with one lithium bromide LBARU for ships’ accommodation.
In the case of the balancing of cooling load from “dead” product compartments, the overall occurring cooling load is absorbed by the circulation air, which is then cooled in the pertaining air coolers (AC) by the CaCl2 brine, and furthermore, it is refrigerated in evaporator AARUs by the very rich ammonia water mixture. Transportation of “live” product is characterized, in addition to moderate temperature storage, by the removal of the products of their own metabolism followed by fresh air, and with continuous control of relative humidity air in compartments. Regulating the moisture content while balancing the ongoing cooling load and the continuous (or periodic) supply of fresh air, in addition to cooling it, also requires its heating in the air heater.
In the secondary cooling circuit of both cases, the secondary refrigerator CaCl2 brine circulates by an electric motor driven cooling brine pump (CBP), while the compartment air is recirculated by an electric motor-driven recirculation air fan (RAF), and in transportation of “living” perishable goods the extraction of a part of the air saturated with products of their metabolism is ensured by a special extraction fan (EF). Although many times, depending on market conditions, the reefers simultaneously carry two or more different moderate refrigerated commercial cargoes to separate compartments, in order to define the appropriate cooling systems, it is assumed that during a voyage cycle reefer transport only one type of commercial cargo is transported. Here is an unambiguously defined sailing route. The transport of five different types of cargo is considered: dead-frozen meat (FM) and meat in vacuum packs (MVP), live perishable bananas (B), citruses (C—citrus fruits, e.g., oranges, lemons, grapefruits) and deciduous fruits (DF—kiwi, pears and apples). There is a possibility that, on the one hand, during the regime of navigation that the majority of occurring unsteady cooling and heating loads are balanced by MDETES, while on the other hand, during standstill reefer at final terminal destination, balancing of the corresponding generated cooling and heating loads is carried out by either CRP and FSG, respectively, or by installed absorption chillers that are in this case powered by heat flow produced in FSG. ARU cookers and HFO final heather (HFOFH) together with other ship’s steam consumers (OSSC) determine the overall total unsteady heating load Φ h l ( t ) , which is balanced by the dry-saturated steam of ps ~ 8 bar produced in a single-pressure cogeneration system characterized with medium-pressure evaporator (MPE) and steam drum (MPD). In the case of the extended form of the cogeneration plant, whose application is based on the using of corrosion-resistant materials, additional main items are included, including a low-pressure steam drum (LPD) and low-pressure evaporator (LPE) that utilizes the remained waste heat from diesel motor exhausted gas after its passage through MPE. Simplified functional diagram Figure 1 illustrates the plain original MDETES, the application of which could be found on ships for transport of moderate refrigerated cargoes. Based on the submitted schemes, the following can be concluded. The reference design version of the ship’s heating and cooling system is characterized by the required heating effect being produced by the combined steam generator (CSG = FSG + WHRSG), while the cooling effect is produced by a CRP. The introduction of the trigeneration concept offers a trigeneration design version characterized by the fact that during the entire transport cycle the occurring cooling load is balanced by the appropriate ARUs, while by the WHRSG unbalanced amount of the overall unsteady heating load is balanced by the FSG. In order to achieve sufficient redundancy, it is recommended in the case of application of the trigeneration concept to install the CRP of the corresponding cooling effect.

Description of the Trigeneration Scheme

The enclosed Figure 1 illustrates a simplified scheme of generic MDETES. The functional-interactive connection of the involved subsystems with the flows of the involved media—engine exhaust, water vapor, heavy fuel oil, cooling freshwater, cooling brine, cooling sea, and air of accommodation air conditioning system (AACS) and air of cargo storage air conditioning systems (CACS)—does not illustrates interactive links within the building subsystems of the involved absorption devices (AARU and LBARU).
During navigation, the two-pressure WHRSG is driven by the exhaust gases of the MDE, which simultaneously produces medium-pressure and low-pressure water vapor to balance the corresponding unsteady heat loads—HFOFH and AARU with medium-pressure steam, LBARU, AHa, AHc, and OSSC (Other Ship’s Steam Consumers, such as heaters in HFO S&TS (HFO Storage and Treatment System), as well in HPW (Hot Portable Water)) with low-pressure steam—whereas in cases of insufficient cogeneration effect, the unbalanced part of the occurring heat load (including the heat load of the heavy fuel heater of the fired boiler (FHFSG)) is balanced with medium-pressure steam generated by the fired boiler. In both cases, the medium-pressure and low-pressure condensate are transferred to the hot well (HW) by the appropriate condensate lines. In the air conditioning systems of housing and storage, the extraction of the precisely prescribed part of the metabolic products of saturated air is carried out by extraction fans (EFa and EFc). The amount of the fresh air replacing extraction saturated air and that is mixing with the recirculated saturated is regulated by the regulation flap (RF), immediately in front of recirculation air fans (RAFa and RAFc). The mass flow of this mixture, its temperature, and relative humidity at the outlet of the associated conditioning units, must lie in the prescribed range of comfort parameters to ensure the balance both of humidity and the resulting overall heat load of the air-conditioned space. That is achieved by an appropriate distribution system and by appropriate cooling and/or heating, as well as possible humidification of the same. The unsteady heat loads in the air coolers (ACa and ACc) are balanced by the secondary coolers (fresh water in AACS and brine in CACS), while unsteady heat loads in the associated air heaters (AHa and AHc) are balanced by the low-pressure steam. The secondary coolers are cooled in the evaporators of the respective ARUs. The total heat loads there are increased by the thermal gains of the pertaining conductor systems. The unsteady cooling effect’s generation in ARU evaporators causes the occurrence of unsteady heat loads of contained heat exchangers: cookers that are balanced by water vapor from the cogeneration system, as well absorbers and condensers that are balanced with the cooling sea supplied by the electrically driven cooling brine pump (CBP) of low power rate. In practice, by applying the trigeneration concept, the overall thermal load of the AACS and CACS is balanced by generating the appropriate cogeneration heat output. When transporting “dead” products, there is no need for air conditioning of the storage air, so despite the significantly lower storage temperatures than for “living” products, the thermal load of the storage of “living” products is much higher. This is primarily due to the appearance of increased infiltration heat gains inherent to refreshing. Practically, when maintaining the prescribed micro-climatic conditions of the storage space, a slight change in the intensity of refreshing significantly changes the occurring heat load of the air conditioning system, which may lead to the need for a short-term (but permissible), slight reduction of refreshing in conditions of temporary insufficiency of cogeneration. When modeling MDETES, this scenario was also taken into account, which is confirmed by the results of the simulations of the occurring heat load during navigation for the set refreshing intensities.
The Figure 1 symbols are; AF—Air Filter, AC—Air Cooler, AH—Air Heater, CBP—Cooling Brine Pump, CSP—Cooling Sea Pump, H—Humidifier, HFOFH—Heavy Fuel Oil Final Heater, OSSC—Other ship’s Steam Consumer, RF—Regulation Flap, RAF—Recirculation Air Fan, FP—Feed Pump, O—Orifice, WHRSG—Waste Heat Recovery Steam Generator, FSG—Fire Steam Generator, EF—Extraction Fan, CFWP—Cooling Fresh Water Pump.

3. Mathematical Modeling Methodology of Energy Interactions of Ship and Environment during Navigation

In order to calculate the performance of the proposed system, models of each respective component are built. Modeling and simulation of the complex physical system are made in Wolfram Mathematica environment.
During the navigation on the predefined navigation route, the ship’s energy system balances the occurring unsteady loads which are fluctuating in a wide range of values due to the action of unsteady and intermittent excitation of the occurring characteristic scalar and vector quantities of the environment, as illustrated in Figure 2. The contained sizes are: vectors— v s c ( φ r , t Y ) and v w ( φ r , t Y ) are sea current and wind velocity, respectively (these vectors are aligned in the tangential plane of the sphere), q s i ( φ r , t Y ) is vector of the solar irradiance; scalar— ϑ a Y ( φ r , t Y ) , Δ ϑ a d ( φ r , t Y ) and Δ ϑ a m ( φ r , t Y ) are yearly average daily air temperatures (DAT), yearly average temperature differences between max and min DAT, monthly average temperature differences between yearly average DAT, respectively; χ s ( μ r , t Y ) —the saturation level of the moist air; ϑ s d ( φ r , t Y ) —sea temperature at the depth d measured from the free surface; Γ c ( φ r , t Y ) and Ψ c ( φ r , t Y ) —sky coverage with the clouds and the attenuation of the sun’s radiation due to clouds, respectively; H s ( φ r , t Y ) and T s ( φ r , t Y ) —significant wave height and period, respectively.
By available reliably statistically processed data of the occurring environment values over a defined navigation route [38], it is possible to create appropriate mathematical formulations that express the dependence of these quantities on calendar time ty and geographical position with given either latitude φ or longitude μ. The mathematical model for each defined navigation route is a precise mathematical explicit approximate functional dependence between latitude and longitude—φ(μ) or μ(φ), respectively—either in the form of a polynomial or a trigonometric order with sine terms, as follows:
μ r φ r = i n a i φ r i   or   μ r φ r = i n a i sin i φ r + ε i ,
while in the case of sailing on an orthodrome, μrr) is defined by the following formula:
μ r φ r = μ A + arcsin sin μ B μ A cot φ B φ A tan φ r ( t ) φ A .
Based on that, any emerging j-th scalar or vector item of the environment that interacts with the ship’s power system on the navigation route can be represented in the form of the product of trigonometric series of sine functions depending on tγ and φr as follows:
S j t Y , φ r = i t n t i φ n φ A i t sin i t ω t t Y + ε i t A i φ sin i φ φ r + ε i φ
Thus, for example, based on the obtained statistical data for the characteristic ambient air temperatures, it is possible to create a mathematical model for the atmospheric air temperature on a certain navigation route, as follows:
ϑ a φ r , t Y = ϑ a y φ r , t Y + 1 2 Δ ϑ a m φ r , t Y sin ω s m t Y + γ s m φ r , t Y 1 2 Δ ϑ a d φ r , t Y cos ω s d t Y + γ s d φ r , t Y
For the fluctuation period of the ϑ a m φ r , t Y , it seems appropriate to take half an average synodic month τ ¯ s m = 29.55 τ ¯ s d , where the length of the mean solar day is τ ¯ s d = 24 h , from which it follows for an unsteady circular frequency ω s m = π / τ ¯ s m . The fluctuation period for Δ ϑ a d ( φ r , t Y ) can be obtained by the fact that the daily temperature minimum occurs just after dawn, whose unsteady phase shift is γ s d ( φ r , t Y ) , and which possesses the apparent circular frequency ω s d = 2 π / τ ¯ s d . Using this methodology, models of characteristic air temperatures were developed, which are illustrated in the Figure 3. Similarly, Figure 4 illustrates wind speed and significant wave height.
Ignoring the geoid shape due to the flatness of the sphere, the navigation route can be defined with sufficient accuracy by a position vector which, with a defined dependence μ r ϕ r , taking into account the radius of the sphere re = 6376.5 km, becomes:
r r φ r = r e cos φ r cos μ r φ r i + cos φ r sin μ r φ r j + sin φ r k ,
According to, the sailing speed can be defined as follows:
v s t = d r r φ r t d t = v s φ r t , φ ˙ r t ; φ ˙ r t = d φ r t d t ,
The overall navigational resistance increased by the amount of propulsion reduction for a given ship’s loading condition T lo = R t o t / 1 t ^ , which is charged to the propeller, except for the navigation speed v s and the continuous increase in resistance due to hull pollution and fouling, depends on the occurring vector quantities: wind v w , sea currents v s c , and waves v s w [39,40], which according to (3) are dependent on t Y = t + τ 0 and φ r . The propeller-generated thrust force T p required to balance this load and the required torque of the propeller Q p , in addition to its speed n or wing pitch p, depend on the average sea inflow speed into the propeller disk v a = v s 1 w , which depends on the modulus of sailing speed v s and the mean coefficient of wake w, whose functional dependence as well as ship’s thrust deduction coefficient t ^ are defined by the expression:
w = w v s , v w , v s c , v s w t ^ = t ^ v s , v w , v s c , v s w
Taking into account the above mathematical formulations, the quasi-static characteristics of thrust and torque take on generalized functional dependencies:
T p = T p n , v a = T p n , φ r t , φ ˙ r t , t , τ 0 Q p = Q p n , v a = Q p n , φ r t , φ ˙ r t , t , τ 0
The generalized functional dependence of the overall quasi-static thrust load (taking into account the total time spent in the service after the docking τ s = τ 0 ) is formulated by the expression:
T l o = T l o v s , v w , v s c , v s w , t Y , τ s = T l o φ r t , φ ˙ r t , t , τ 0 ,
The actual torque Q d e = Q d e k F , n of the MDE which balances the oncoming propulsion load increased by the counter friction torque of the propeller shafting Q c s = Q c s n , except for n, depends on the motor load factor k F which is on the fixed geometry of the propeller (FPP, p = const.), the only control variable which ensures (with constant fluctuation of n and v s ) safe, reliable and seakeeping acceptable navigation.
Ignoring the emerging capacitances within the DEPS (diesel engine propulsion system), to conduct credible navigation simulations through a realistic environment, it is necessary to set up an appropriate system of dynamic equilibrium of thrust and torque equations as follows:
m s = d v s φ r t , φ ˙ r t d t = T p φ r t , φ ˙ r t , t , τ 0 T l o φ r t , φ ˙ r t , t , τ 0 ,
I p = π 30 d n ( t ) d t = Q d e k F t , n t Q c s n t Q p n , φ r t , φ ˙ r t , t , τ 0 ,
where: m s is the mass of the loaded ship increased by the mass of the surrounding water affected by the movement of the ship (about 10% of the mass of the ship), and I p is the polar moment of inertia of all rotating masses of the propulsion system increased by the inertia of the propelled sea mass.
This system of nonlinear differential equations (2nd order φ r ( t ) and 1st order n ( t ) ), solvable exclusively numerically based on the prescribed initial conditions φ r ( 0 ) , φ ˙ r ( 0 ) , n ( 0 ) , gives unknown time dependencies φ r A B ( t ) and n AB ( t ) during navigation between destinations A and B, sailing duration τ A B . Furthermore, by the numerical processing of φ r A B ( t ) , φ ˙ r A B ( t ) is obtained and according to (6) unsteady sailing speed v s A B ( t ) is calculated; additionally, based on both the n A B ( t ) and SMCR (Service Maximum Continuous Rating) of MDE, the time dependency control variable k F A B ( t ) is also obtained. In principle, the simulation of navigation from B to A is the same, except that the prescribed initial parameters and conditions are different, and thus the functional dependences of unsteady quantities on time φ r A B ( t ) , n B A ( t ) , v s A B ( t ) , k F B A ( t ) are obtained with the duration of navigation τ B A ( t ) , which are illustrated in Figure 5.
With the estimated time intervals of stay in the final destinations τ A and τ B , the total duration of the kth transport cycle is: τ t c k = τ A k + τ A B k + τ B k + τ B A k , so when simulating navigation in the (k + 1)’s transport cycle, the time parameter τ 0 is obtained:
τ 0 ( A ) k = k = 1 n k 1 τ t c ( k 1 ) , τ 0 ( A B ) k = τ 0 ( A ) k + τ 0 ( A B ) k , τ 0 ( B ) k = τ 0 ( A ) k + τ 0 ( A B ) k , τ 0 ( B A ) k = τ 0 ( A ) k + τ A k + τ A B k + τ B k
By using models for air and sea temperature (according to Equation (3)) and ship sailing speeds during the transport cycle (obtained from simulated navigation through realistic surroundings) the unsteady profiles of air and sea temperature are obtained as illustrated in Figure 6.

3.1. Quasi-Static Effect of Cogeneration System Scheme

Taking into account predictable values of mass flow m ˙ e g and temperature ϑ e g ( t ) of MDE exhaust gases at the SMCR (Service Maximum Continues Rating), and by estimating techno-economically acceptable minimum temperature difference Δ ϑ min at the WHRSG outlet, the main cogeneration design parameters can be determined, such as exhaust gas outlet temperature and heat exchanger area A. Its required area A is determined for the selected configuration, which is interactively related to the overall heat transfer coefficient k.
During navigation, the fluctuation m ˙ e g and ϑ e g ( t ) at the WHRSG inlet, with the determined temperature of the generated steam ϑ s and the selected configuration of the overall flow surface A, leads to volatility k ( t ) , so based on the energy balance and the applicable mathematical formulation for heat exchange via mean logarithmic temperature difference (MLTD- Δ ϑ m ), unsteady cogeneration heat effect is obtained:
Φ c o t = C e g ( t ) ϑ e g ( t ) ϑ s 1 exp K ( t ) C e g ( t ) ,
where: C e g ( t ) = m ˙ e g ( t ) c p ( t ) is the heat capacity of exhaust gases, and K ( t ) = k i ( t ) A i is the exchanger heat transfer capability.
The total unsteady overall heat transfer coefficient k i ( t ) corresponding to Ai is defined by the expression:
k i t = 1 α s t + r i λ t ln r e r i + r i r e 1 α e g ( t ) 1 ,
where the appropriate heat transfer coefficients are α s ( t ) and α e g ( t ) on both the steam and exhaust side, and α e g ( t ) is markedly dependent on the oncoming mean temperature ϑ ¯ e g ( t ) and on the mean flow rate through the exchanger section v ¯ e g ( t ) .
Taking into account the functional dependence of the emerging operating parameters of the diesel engine k F ( t ) and n ( t ) , as well as the ambient air temperature ϑ a ( t ) during navigation, their involvement in (10) easily results in the emerging thermal effects of the cogeneration system Φ c o A B ( t ) and Φ c o B A ( t ) .
When using a two-pressure cogeneration system, the steam production process takes place in such a way that even when Φ h l ( t ) < Φ c o ( t ) < Φ H L t , all available flue gas heat flow is used to produce medium-pressure steam, which results in the flue gas temperature at the inlet to the low-pressure evaporator being exactly equal to the flue gas temperature outlet of the medium-pressure evaporator ϑ ¯ e g i ( t ) = ϑ ¯ e g o ( t ) :
ϑ e g t = ϑ s + ϑ e g ( t ) ϑ s exp K ( t ) C e g ( t ) ,
Thus the available heat output of the low-pressure evaporator is obtained:
Φ ˜ c o t = C e g t 1 exp K ˜ ( t ) C e g ( t ) ϑ s ϑ ˜ s + ϑ e g t ϑ s 1 exp K ( t ) C e g ( t ) ,
where K ˜ t = k ˜ i t A ˜ i is the heat transfer capability of the low-pressure evaporator, the flow surface A ˜ i and the total heat transfer coefficient k ˜ i t . Using mathematical models expressed in formulas (13) and (16), the unsteady MP cogeneration effects Φ c o A B ( t ) and Φ c o B A ( t ) (during navigation; from A to B and from B to A, respectively), as well unsteady LP cogeneration effects Φ ˜ c o A B ( t ) and Φ ˜ c o B A ( t ) (during navigation; from A to B and from B to A, respectively) for 1st and 2nd transport cycles were calculated and shown in Figure 7.

3.2. Quasi-Static Heating Loads Scheme

The total unsteady heat load of the ship’s power system during navigation between the final destinations A and B Φ H L ( t ) = Φ h l ( t ) + Φ ˜ h l ( t ) consists of an MP heat load balanced with MP steam Φ h l t and an LP heat load balanced with LP steam Φ ˜ h l t , defined by the corresponding terms as follows:
Φ h l ( t ) = Φ C A A ( t ) + Φ F H H F O ( t ) ,
Φ ˜ h l ( t ) = Φ ˜ C L B A ( t ) + Φ ˜ A H a ( t ) + Φ ˜ A H c ( t ) + Φ ˜ O S S C ( t ) ,
where unsteady heat loads are contained as follows: Φ C A A t represents AARU cookers; Φ F H H F O t represents HFOFH; Φ ˜ C L B A t represents LBARU cookers; Φ ˜ A H a ( t ) and Φ ˜ A H c ( t ) represent AH in the AACS and CACS, respectively; and Φ ˜ O S S C ( t ) represents OSSC heaters. Unsteady heat loads of AARU and LBARU cookers are defined by the terms:
Φ C A A ( t ) = Φ C E c ( t ) C O P A A ( t ) , Φ C L B A ( t ) = Φ C E a ( t ) C O P L B A ( t ) ,
where Φ C E c t and Φ C E a t are cooling effects of AARU-a and LBARU, respectively, and C O P A A t and C O P L B A t are their coefficients of performance (COP). C O P A A t and C O P L B A t are functionally dependent on the equilibrium pressures in the evaporator p C E t , and on the unsteady temperature of the cooling sea water ϑ c s = ϑ s d ( t ) , which limits the evaporation range for a certain mixture when there is the steady heating source of the cooker (water vapor from the cogeneration system). As ϑ c s is lower, the possible evaporating range is higher, which increases COP, as illustrated in Figure 8.
Unsteady cooling effects corresponding to the occurring heat loads of AARU and LBARU evaporators are:
Φ C E c t = i = 1 n s Φ A C s i t + P c b p + Φ b c s t ,
Φ C E a t = Φ A C a t + P c w p + Φ w c s t .
where: Φ A C s i t and Φ A C s i t are the heat loads of air coolers in climatization systems of cargo spaces and accommodation, respectively; P c b p and P c f w p are the mechanical powers of the cooling brain pumps (for cargo space) and cooling fresh water pump (for accommodation), respectively; and Φ f w c s t and Φ b c s t are the heating incomes of the conducting systems from pertaining surroundings.
The heat loads of both climatization system (storage space and accommodation) depend on both predefined microclimatic comfort parameters (mainly air temperature and relative humidity), and total sensible and infiltration heating loads. In principle there are two emerging characteristics operating scenarios, either humidification or dehumidification of the wet air mixture as a result of mixing of the recirculation air from climatization space and fresh air from surroundings (Figure 9).
Consequently, heat loads of air coolers are defined according to the corresponding formulas for cases of humidification and dehumidification as follows:
Humidification and cooling, if x c s ( t ) x M c ( t ) > 0 , and Φ A C ( t ) ( 1 ) < 0 ,
Φ A C ( t ) ( 1 ) = m ˙ s a h c s ( t ) h M c ( t ) x c s ( t ) x M c ( t ) c w ϑ w ( t ) + Φ h l s ( t ) P R F Φ d s ( t ) ,
Dehumidification and heating, if x a c ( t ) x M c ( t ) < 0 , x a c ( t ) = x c s ( t ) m ˙ w e / m ˙ s a ,
Φ A C ( t ) ( 2 ) = m ˙ s a h d p ( t ) h M c ( t ) P R F Φ d s ( t ) ,
where moisture contents are x ( t ) and specific enthalpies are h ( t ) (per kg dry air); x c s ( t ) and h c s ( t ) are for air of air conditioned space; x M c ( t ) , x M h ( t ) and h M c ( t ) , h M h ( t ) are for wet air mixture obtained by mixing of the recirculated air and fresh environmental air, in the corresponding cases when cooling or heating is carried out; x d p ( t ) and h d p ( t ) are for saturated wet air after excessive moisture extraction (from the condition that is x d p ( t ) = x a c ( t ) ); and ϑ w ( t ) is temperature of water for humidification.
Specific enthalpies of wet air when φ 1 are defined by the following formula:
h i ( t ) = c p a ϑ i ( t ) + x i ( t ) r 0 + c p s ϑ i ( t ) ,
where: c p a , c p s and c w are specific heat capacities at p = const. of the dry air, water steam and water, respectively; r 0 is specific evaporation heat of water steam at 0 ° C ; and ϑ i ( t ) is wet air temperature.
Further specific enthalpies and moisture contents of the wet air mixtures are defined by the corresponding expressions as follows:
h M ( t ) = g f a h f a ( t ) + g r a h r a ( t ) ,
x M ( t ) = g f a x f a ( t ) + g r a x r a ( t ) ,   where   x r a ( t ) = x c s ( t )
and according to (24), (25) and (26) the temperature of the mixture is defined as follows:
ϑ M ( t ) = h M ( t ) r 0 x M ( t ) c p a + c p s x M ( t ) = g f a h f a ( t ) + g r a h r a ( t ) r 0 g f a x f a ( t ) + g r a x r a ( t ) c p a + c p s g f a x f a ( t ) + g r a x r a ( t ) ,
where: g f a = m ˙ f a / m ˙ s a and g r a = m ˙ r a / m ˙ s a are mass fractions of the fresh air and recirculated air, respectively; and m ˙ f a , m ˙ r a , m ˙ s a = m ˙ f a + m ˙ r a are mass flows of the fresh air, recirculated air and air mixture, respectively.
Finally, there are: P R A F —power of recirculating fan electromotor, Φ d s ( t ) —heat gains of the air distribution ducting system, and Φ h l s ( t ) —sensible heating loads of the climatizated space.
In the occurring air conditioning processes of both air conditioning systems (for cargo space and accommodation), it is necessary that the air is heated, apart from in the case of humidification plus cooling (Figure 9a); hence, heating loads of air heaters are defined as follows:
Humidification and heating, if x c s ( t ) x M h ( t ) > 0 , and Φ A H ( t ) ( 1 ) > 0 ),
Φ A H ( t ) ( 1 ) = m ˙ s a h c s ( t ) h M h ( t ) x c s ( t ) x M h ( t ) c w ϑ w ( t ) + Φ h l s ( t ) P R F Φ d s ( t ) ,
Dehumidification and heating, if x d c s ( t ) x M c s ( t ) < 0 and Φ A H ( t ) ( 2 ) > 0 ,
Φ A H ( t ) ( 2 ) = m ˙ s a h c s ( t ) h d p ( t ) + Φ h l s ( t ) ,

3.3. The Sensible Heating Load

This load involves heat flows exchanging between the air conditioned space and surroundings Φ s h ( t ) , and internal heat gains Φ i h ( t ) including the respiration heat flow Φ r ( t ) , either from live products (cargo hold space) or people (accommodation), as well from contained energized equipment such as lighting, etc. For defining the Φ s h ( t ) , a concept of the quasi-static thermal network is used for both cargo hold space and accommodation, illustrated for the latter in Figure 10. By implementation of this concept, sensible heat load of a ship’s air conditioned space is defined by the formula:
Φ s h ( t ) = k = 1 n a ϑ a ( t ) ϑ s a ( t ) R k o a k A k q s i k ( t ) R k o + R k i + p = 1 n p ϑ a ( t ) ϑ a p ( t ) R T p + s w = 1 n s ϑ s w ( t ) ϑ s a ( t ) R T s w ,
where contained heat resistances are as follows:
R k O = 1 A k 1 α k O + i = 1 n k i δ i λ i , R k I = 1 A k α k I , R T p = 1 A p 1 α p I + i = 1 n p i δ i λ i + 1 α p O , R T s w = 1 A s w 1 α s w I + i = 1 n s w i δ i λ i + 1 α s w O .
Involved items in formulas (30) and (31) are: α k O , α k I , α p O , α p I and α s w O , α s w I , the convective heat transfer coefficients on the outer (indices o) and inner (indices I) flat surfaces of ship’s enclosures that are in interacting with; environmental air (indices k), surrounding air non air conditioned ship’s compartments (indices p), and sea (indices sw); A k , A p and A s w are areas of the ship’s flat surfaces of the enclosures indexed by k, p and sw, respectively; and δ i and λ i are thicknesses and heat conductivity of the involved multilayer enclosures ( n k i , n p i and n s w i are pertaining layers’ number enclosures indexed, respectively, by k, p and sw), respectively.
Further, ϑ a p ( t ) is temperature of the pth non-air conditioned ship’s compartment and ϑ s w ( t ) is sea temperature on the place corresponding to sea depth at the center of swth flat surface wetting by sea water, and a k is absorbance of the kth ship’s flat surface exposed to solar irradiation, the intensity of which is displayed as q s i k ( t ) .

3.4. Solar Irradiation

For any kth ship’s flat surface exposed to solar radiation intensity of the overall acting solar irradiation, Φ s i k ( t ) is defined by the following formula:
Φ s i k ( t ) = A k q s i k ( t ) = A k q d ( t ) cos ϕ k ( t ) + q d i f ( t ) ( 1 + cos η k ) / 2 + q r ( t ) ( 1 cos η k ) / 2 ,
where: q d ( t ) , q d i f ( t ) and q r ( t ) are intensities of the solar irradiation directing normally on the sphere tangential surface (Figure 11), diffuse sky irradiation and overall reflected solar irradiation by the sea surface, respectively. These components of irradiation are defined by the corresponding expressions as follows:
q s i k ( t ) = q d e ( t ) exp Ψ c ( t ) B ( t ) csc Θ ( t ) ,
q d i f ( t ) = q d ( t ) D ( t Y ) , .
q r ( t ) = r s s ( t ) q d ( t ) sin Θ ( t ) + q d i f ( t ) ,
where D ( t Y ) = 0.0904 0.04116 cos ω G ( t Y t o ) is the diffuse irradiation parameter depending on changeable amounts of the moisture and dust particles in the atmosphere [41], q d i f ( t ) and q r ( t ) are the intensities, and r s ( t ) = 1 0.95 sin Θ ( t ) is the reflectivity of the sea surface.
By using the solar constant q S = 1.373 kW / m 2 , the direct extraterrestrial irradiation is defined by approximated expression q ˜ d e = q S [ 1 + 0.033 cos ( ω Y t Y ) ] [42], or more precisely by the formula:
q d e ( t ) = σ T s 4 r ¯ S 2 [ R S E ( t ) r ¯ S r ¯ a ] 2 ,
where r ¯ S = 6.923 10 5 km is the mean Sun chromosphere radius, r ¯ a = 6467.5 km is the mean radius of the Earth’s atmosphere, σ T s 4 r ¯ S / R ¯ S E 2 = q S , σ = 5.67 10 8 W / ( m 2 K 4 ) is the Stefan–Boltzmann constant, T S 5791 K is the absolute chromosphere temperature, and module of radii vector of the Earth’s path around of the Sun is:
R S E ( t ) = a ( 1 2 ) 1 + cos Ω ( t Y , t 0 ) ,
where a = 149.5 10 6 km are the big semi-axes and is the = 0.017 eccentricity of the elliptical path, and finally the revolution angle of the Earth around of Sun is defined, according to [43], by the formula:
Ω ( t Y , t 0 ) = ω Y ( t Y t o ) + i = 1 2 k i sin 2 ω Y ( t Y t o ) ,
where k 1 = 0.033985 and k 2 = 3.61128 10 4 are contained constants, ω Y = 2 π / τ Y is the mean angular velocity of Earth’s revolution around the Sun, and τ Y 8766 h is the revolution period.
The atmosphere attenuation factor is defined by B ( t Y ) = 0.17164 0.034686 cos ω G ( t Y t o ) [40], while irradiation attenuation produced by the cloudiness depends on the factor Ψ c ( φ r , t Y ) .
That is simulated by applying formula (3) on the obtainable data from isonephs maps (line connecting the places with equal mean cloudiness), which is illustrated in Figure 12a and by use results for navigation trough realistic surroundings, the unsteady attenuation factor while transport cycle is obtained as shown in Figure 12b.
Sun high sin Θ ( t ) is defined by the scalar product of unity vectors, the normal of the tangential surface of the sphere n 0 ( t ) and acting Sun radiation q 0 ( t ) (Figure 11):
sin Θ ( t ) = n 0 ( t ) q 0 ( t ) ,
where unity vectors are involved as follows:
n 0 ( t ) = cos φ ( t ) cos μ ( t ) i + cos φ ( t ) sin μ ( t ) j + sin φ ( t ) k ,
q 0 ( t ) = cos δ ( t ) i + sin δ ( t ) k ,
where δ ( t Y ) is the Sun’s inclination according to:
δ ( t Y ) = δ 0 cos Ω ( t Y , t 0 ) Ω δ 0 = δ 0 cos ω Y ( t Y t 0 τ ˜ δ 0 ) + i = 1 2 k i sin i ω Y ( t Y t 0 ) .
Further, cos ϕ k ( t ) is defined by formula:
cos ϕ k ( t ) = n k 0 ( t ) q 0 ( t ) ,
where n k 0 ( t ) is the unity normal of kth arbitrarily placed ship’s flat rectangular surface, defined by vector’s product of the pertaining unity base vectors as follows:
n k 0 ( t ) = r β 0 ( t ) × r η 0 ( t ) , r β 0 ( t ) = cos β k r μ 0 ( t ) sin β k r φ 0 ( t ) , r η 0 ( t ) = cos η k sin β k r μ 0 ( t ) + cos η k cos β k r φ 0 ( t ) + sin η k n 0 ( t )
In this formula, contained orts are:
r μ 0 ( t ) = 1 r μ n 0 ( t ) μ , r φ 0 ( t ) = 1 r φ n 0 ( t ) φ .
Further, contained angles β k and η k present azimuth and elevation of a kth ship’s flat surfaces, respectively, as shown in Figure 11b.
Other parameters involved in above expressions are: δ 0 = 23.45 0 —min. inclination for north hemisphere, t 0 = 68 h —time shift for perihelion, τ ˜ δ 0 = 230 h —time shift for winter solstice to which corresponds angle shift Ω δ 0 = Ω δ 0 ( τ ˜ δ 0 ) = 12.66 0 ; finally, of the course in modelling through a transport sequence, ts timing relations t Y = t + τ 0 t s must be taken into account. Figure 13 illustrates overall solar irradiations for the (a) navigation A to B, and (b) navigation B to A.
Finally, by applying the above-mentioned for characteristic quasi-static ship’s energy items during the navigation through any transport cycle, the example illustrated by Figure 14 for unsteady heat load is obtained.
Taking into account that all design parameters of the involved heat exchangers are predetermined on the extreme base values of pertaining heat loads, such as heat exchange areas, overall heat transfer coefficients, mean logarithmic temperature, etc., and assuming that heat capacities of the involved cooling (heating) media (fresh water and brain), as well as heat transfer capabilities of exchangers, are unchangeable during navigation, it can be possible to determine the corresponding quasi-static balanced temperatures on the side cooling media contained in LBARU and AARU.
As an example, Figure 15 illustrates design and unsteady temperatures of the involved cooling media in the air conditioned system of the accommodation, while expressions (46) and (47) define corresponding heat transfer capability and capacities of conditioned air and cooling fresh water:
K e x = k e x A e x = Φ e x max Δ ϑ m d = Φ e x max δ 1 d δ 2 d ln δ 1 d δ 2 d ,
C a c ( t ) = Φ A C ( t ) ϑ a c ( t ) ϑ a c ( t ) , C c w = Φ A C ( t ) ϑ s w ( t ) ϑ s w ( t ) ,
where δ 1 d and δ 2 d are design temperature differences on the ends of the counter flow exchanger.
Further, for characteristics of unsteady balancing temperatures the corresponding expressions are determined:
The inlet temperature of the cooling fresh water in AC is:
ϑ c w ( t ) = ϑ a c ( t ) + Φ A C ( t ) C c w 1 exp K A C 1 C c w 1 C a c ( t ) 1 ,
The outlet temperature of the cooling fresh water in AC is:
ϑ c w ( t ) = ϑ a c ( t ) + Φ A C ( t ) C c w exp K A C 1 C c w 1 C a c ( t ) 1 exp K A C 1 C c w 1 C a c ( t ) 1 ,
The temperature of the cooling media in the cooling evaporator is:
ϑ s ( t ) = ϑ ˜ c w ( t ) Φ C E ( t ) C c w exp K C E C c w 1 exp K C E C c w 1 ,
The temperature of the heating media (LP water steam) in the air heater is:
ϑ s ( t ) = ϑ a h ( t ) + Φ A H ( t ) C s a 1 exp K A h C s a 1 .
By adding equal amounts of the heat gains of the cooling media conducting system Φ C S ( t ) = P c p + Φ c s t on the inlet and outlet of the cooling evaporator, ϑ ˜ c w ( t ) = ϑ c w ( t ) + 0.5 Φ C S ( t )   ϑ ˜ c w ( t ) = ϑ c w ( t ) + 0.5 Φ C S ( t ) are obtained for the inlet and outlet temperatures of the secondary cooling media (fresh water or brine) in the cooling evaporator. By using above mentioned quasi-static balanced temperature of the involved media in AACS are obtained, as illustrated in Figure 16.

4. Preliminary Economic, Energy, and Environmental Indicators of the Application of the Trigeneration Concept

By applying the appropriate original mathematical model for unsteady thermal loads of the ship and the unsteady cogeneration heating effect, the techno-energy applicability of the trigeneration concept on the ships intended for the transport of moderate refrigerated bulk “live” and “dead” commodities is examined. In this sense, a one-year reefer service through a simplified orthodromic navigation route between terminal A (corresponding to Marcaibo, Venezuela) and terminal B (corresponding to an estuary of the river Laba, Hamburg) for reefer which main particulars are given in Table 1 is simulated. The simulation includes 14 transport cycles; both navigations loaded reefer and steaming in the ballast, where the continuous change of energy effects and energy loads that are balanced are taking the place. The application of the double-stage cogeneration system that balances both medium-pressure (MP) heating load (thermal load resulting from AARU cooker and the HFOFH) and low-pressure (LP) heating load (thermal load that results from all other consumers of heat energy including LBARU cooker) is considered for the typical reefer as follows:

4.1. Energy Sufficiency

Combined display of total unsteady heating load of a loaded ship with characteristic bulk cargo (Banana—B, Citrus—C, deciduous fruit—DF, frozen meat—FM), and unsteady cogeneration heating effect, for two characteristic transport cycles, are illustrated in Figure 17. It is shown that the MP cogeneration heating effect is sufficient to balance the generated unsteady MP heating loads in the transport of all above-mentioned commodities.
Possible, short-term insufficiency of the double-stage cogeneration effect in balancing the generated heating load can easily be remedied by the short-term changing of microclimate parameters (relatively small increase in relative humidity, or a small reduction in the proportion of fresh air g f a ), as illustrated in Figure 18.
The total LP unsteady heating load of the reefer transporting characteristic com-modities for certain transport cycles is obtained by summation of the contained basic LP heating loads, as illustrated in Figure 19.
The corresponding total (low-pressure and high-pressure) heating loads for the ballast voyage when maintained storage temperature of cargo was predetermined by the product to be transported in the opposite direction, for the most demanding heating loads (frozen meat—FM, deciduous fruit—DF, orange—O, and bananas—B), is shown in Figure 20. In addition, Figure 5 shows a common view unsteady heating effect of the medium-pressure cogeneration system and the characteristic load caused by the overall thermal load of the ship, which shows energy sufficiency of the high-pressure cogeneration system in balancing the overall heating load of the ship. During practically any condition in the ballast voyage, there is no need for activating the low-pressure cogeneration system.

4.2. Economic Impacts

Furthermore, in brief, the simplified economic review of the positive economic effects of applying DETES based on two-pressure cogeneration systems for reefer ships is given.
In this regard, primarily referring to the estimated scenario of the commercial engagements of a ship, we determine fuel savings compared to competitive design solutions with CRPs in the CACS and in the AACS. Fuel saving corresponds to the fuel consumption of the CRP in the CACS during the steaming in ballast for several characteristic transport cycles, illustrated in Figure 21, while saving fuel for the CRP in AACS for several characteristic transport cycles, illustrated Figure 22.
Fuel consumption accounting for the operation of the compressor plant to balance the occurring storage cooling load during navigation of a loaded ship for two characteristic transport cycles, through a one-year period, is illustrated in Figure 23.
By the numerical integration of fuel consumption (mass flow) over the navigation period, the following results are obtained: average fuel savings for the AACS over one transport cycle and one year of service which consists of 14 transport cycles, as well as the total savings for the 25-year ship lifetime period (SLT), were taken into account alongside the yearly trend price of low sulfur heavy fuel oil (LSHFO, containing sulfur s 0 . 5   % ,   m a s . ), p ¯ 2021 500 US $ / t [44], giving a total saving of FSAC = 1.3146 mil. US$. In this case, in which the reefer in one direction navigates loaded with typical product, and the second direction navigates in ballast, the obtained results of the possible combinations of transport over one year of service, by the numerical integration of the curves of mass fuel consumption, are listed in Table 2.
In the economically more favorable case, which is the case of a complete engagement of the ship in transporting of moderate refrigerated commodities, the various possible combinations of transport are illustrated in Table 3.
Taking into account the current specific price [45], CRP, ARU and the CP, which are p C R P = 137 , 000 US $ / MWce , p A R U = 411 , 000 US $ / MWce , and p C P = 120 , 000 US $ / MWhe   , respectively, in the economically most unfavorable case when the 100% redundancy cooling-heating system of the ship is required, both the CRP and ARU rated cooling effects are installed, and for the same reasons FSG has an extra heating effect installed, the following net cost-benefit equation is obtained:
P E B = p L S H F O m F S I A A I L B A I F S G ,
where, IAA, ILBA and IFSG are the present investment values of AARU, and LBARU and FSG, respectively. Contained investment values are defined by the following terms:
I A R U = p A R U Φ C E , Φ C E = Φ A A c e + Φ L B A c e ,
I F S G = p F S G Φ T C Φ C R P ,
where are the nominal cooling effects of the AARU and LBARU; Φ A A c e = 1 . 85   MW and Φ L B A c e = 0.35 MW , correspondingly, while the nominal cogeneration heating effects are; for trigeneration energy concept (DECES+ARU) Φ T C = 3 MW , and for conventional energy concept (DECES+CRP) Φ C R P = 1 MW .
Following the above, the total, present, additional investment cost of ITOT = 1.045 mil. US$ is obtained, which is less than the economic fuel savings achieved by applying LBARU in the AACS (FSAC = 1.3146 mil. US$). Accordingly, for economic gains, the fuel savings m F S obtained by applying AARU in CACS for various transport combinations through the 25-year economic life of the ship, alongside the same savings that are contained in Table 2 and Table 3, are illustrated in Figure 24.
Based on the above it can be concluded that the economic effects of the application of trigeneration energy systems on motor ships intended for the transport of moderately refrigerated commodities are appreciable. By the extrapolating of techno-economic parameters (the price of heavy fuel and reefers’ hold capacity), the value of preliminary economic benefits would be increased. For example, for a reefer twice the capacity of that considered, at the same price of fuel, profits would be nearly doubled, while at a price of fuel at least twice that of the current (which is realistic to expect after a five-year period), economic profit for the considered ship would be double, but for a ship with double the capacity it would increase nearly fourfold.

4.3. Environmental Impact

Most environmentally harmful substances are manifested through the emission of environmentally harmful combustion gases: CO2, CO, SO2, SO3, and various nitrogen oxides NOx. By reducing fuel consumption, these emissions are almost proportionally reduced, but also in the assumed complete combustion of fuel reduced mass amounts of CO2 and SO2 will be easily achieved based on of the following terms:
Δ m C O 2 = c M C O 2 M C m F τ S L T ,             Δ m S O 2 = s M S O 2 M S m F τ S L T
where c = 0.84 kgc/kgf; s = 0.005 kgs/kgf are the average concentrations of carbon and sulfur in heavy fuel, respectively; MCO2 = 44, MC = 12, MS = 32 and MSO2 = 64 kg/kmol are appropriate molar mass; mF is annual fuel savings; and τSLT is the ship’s economic lifetime.
The calculation results for the 25-year ship’s economic lifetime are given in the following Table 4. The application of trigeneration energy systems on ships intended for transporting moderately cooled products considerably contributes to reducing environmentally negative effects.
This confirms the adequacy of diesel engine trigeneration energy systems on motor reefers in balancing the occurring overall heating load during the voyage, as well as the respectable positive economic, energy, and environmental effects of its application. Therefore, the conclusion is that the application of diesel engine trigeneration energy systems on ships intended for the transport of moderately refrigerated cargo is appreciably economically beneficial, and environmentally acceptable. In addition to the considered types of ships, trigeneration energy systems can also be applied on modern passenger ships, where the necessary energy for ensuring a high level of comfort is very high and often exceeds the power of the ship’s propulsion engines.

5. Conclusions

By applying mathematical models of characteristic environment sizes (based on available WMO data), which interact with the ship system and with developed mathematical models of quasi-static characteristics of energy components involved with the diesel engine trigeneration energy system, we created models of unsteady energy balancing during the ship’s characteristic operating intervals. The concept interrelates the trajectory of a vessel and the wave conditions with a detailed described quasi-static model of the energetic, economic and environmental performance of a marine trigeneration unit installed on-board. The developed mathematical models of quasi-static characteristics for single-stage AARUs and LBARUs, enable modelling of unsteady thermal loads of their cookers, absorbers and condensers. Furthermore, the models enable the management of absorption cooling processes in the most energy-efficient way within the emerging unsteady environment.
Within the defined technically possible area of application of absorption cooling, with the application of LBARU for balancing the refrigeration load in the AACS, the application of AARU in CACS is tested for several characteristic moderately refrigerated cargoes: frozen meat, not frozen vacuum packed meat, bananas, citrus and deciduous fruits, where it is shown that on the increasing of the thermal loads, the most influenced variables are: higher intensity of refreshing, lower relative humidity and to a lesser extent the respiratory heat flux of the transported commodity. This implies that the total refrigeration loads when transporting live cargo at the prescribed storage temperature are almost an order of magnitude larger than the corresponding sensible heat loads. Applying these models over a one-year period of service of the ship, which includes 14 repeated transport cycles, it was shown that the trigeneration system is energy sufficient to balance the overall unsteady heat load of ships intended for the transport of moderately refrigerated “dead” and “live” cargo. Quasi-static characteristics of the (CRP) enable the determination of the appropriate fuel consumption that falls on its drive when it balances the unsteady refrigeration load during the navigation.
By integrating the quasi-static fuel consumption curves of the CRP over the intervals of navigation routes for all contained 14 transport cycles, and summing the obtained values, one-year fuel consumption is obtained, i.e., one-year fuel savings. Based on one-year fuel savings, for the estimated economic life of the ship, taking into account the current fuel price as well as additional investment costs of trigeneration and cogeneration plant, significant economic savings were obtained for certain transport combinations, which fully confirm the economically positive effects of the trigeneration concept on ships intended for the transport of moderately refrigerated cargo.
Finally, using the mean reference values for the composition of the LSHFO, and assuming complete combustion, no negligible amounts of reductions in greenhouse gas emissions, especially CO2, were obtained. In conclusion, according to the above, the energy, economic and environmental effects of the application of diesel engine trigeneration energy systems on ships intended for the transport of moderately refrigerated cargo are significantly positive. In addition, the trigeneration energy concept can be applied on modern passenger ships, where the energy amounts required to ensure a high degree of comfort are significant, and often exceed the installation power of propulsion diesel engines.

Author Contributions

Conceptualization, I.G. (Ivan Gospić) and I.P.; methodology, I.G. (Ivan Gospić); software, I.G. (Ivica Glavan); validation, I.G. (Ivan Gospić), I.P. and V.M.; formal analysis, I.G. (Ivan Gospić); investigation, I.G. (Ivan Gospić); resources, I.G. (Ivica Glavan); data curation, I.G. (Ivica Glavan); writing—original draft preparation, I.G. (Ivan Gospić); writing—review and editing, I.P.; visualization, I.G. (Ivica Glavan); supervision, V.M.; project administration, I.P.; funding acquisition, V.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Nomenclature
Latin symbols
Aarea, m2
Batmosphere attenuation factor
Cflow stream capacity, W/K
Ddiffusely irradiation parameter,
Ippolar moment inertia of rotating mass, kg m2
Iinvestment cost, US$
Kheat transfer capability, W/K
Ppower, W
Rtotal thermal resistance, m2K/W
RSEmodule of radii vector between Sun and Earth, m
Qppropeller torque, Nm
Tabsolute temperature, K
Tppropeller trust, N
abig semi-axes of Earth eliptic path, km
cspecific heat capacity, J/(kgK)
gmass fractions a component in the mixture
hspecific enthalpy, J/kg
eccentricity of Earth elliptical path
koverall heat transfer coefficient, W/m2K
kFdiesel engine load factor
mmass, kg
mass flow rate, kg/h
nrotational speed, rev/min
rradius, m
rospecific evaporation heat of water steam at 0 °C, J/kg
qheat flow intensity, W/m2
ttime, s
t ^ trust deduction coefficient
xmoisture content in wet air,
vSpeed, velocity, m/s
wwake coefficient
Greek symbols
Γ sky coverageness with the clouds
Δ ϑ m mean logaritmic tempearture difference, K
Θ Sun hight angle, °
Φ heat flow, W
Ψ cloudiness attenuation factor
Ω revolution angle of the Earth around of Sun, °
α convective heat transfer coefficient, W/(m2K)
β azimuth of ship’s flat surface, °
δ Sun inclination, °
δ i thickness of the ith layer of the multilayer wall, m
η elevation of ship’s flat surface, °
ϑ temperature, ° C
λ heat conductivity, W/(mK)
μ r geographic longitude (on navigation route), °
φrgeographic latitude (on navigation route), °
φrelative humidity of wet air
ϕ angle between Solar irradiation and normal vector of a flat ship’s surface, °
χ the saturation level of the moist air
ψ Sun azimuth, °
σ Stefan-Boltzmann constant, W/(m2K4)
τ time interval, s
ω angular speed, rad/s
Abbreviations
AACSaccommodation air conditioning system
AARUammonia absorption refrigeration unit
ARUabsorption refrigeration unit
ACair cooler
AFair filter
AHair heater
Bbananas
Ccitruses (citrus fruits)
CACScargo storage air conditioning system
CBPcooling brine pump
CFWPcooling fresh water pump
COPcoefficient of performance
CRPcompression refrigeration plant
CSPcooling sea pump
DAdiesel aggregate
DATdaily air temperatures
DEPSdiesel engine propulsion system
DFdeciduous fruits
EFextraction fan
FMfrozen meat
FPPfixed pitch propeller
FSACfuel savings of the accommodation climatization
FSGfired steam generator
HFOheavy fuel oil
HFOFHheavy fuel oil final heater
LBARUlithium-bromide absorption refrigeration unit
LHVlower heating value
LNGliquefied natural gas
LPlow pressure
LSHFOlow sulphur heavy fuel oil
MCRmaximum continuous rating
MDEmain diesel engine
MDECESmarine diesel engine cogeneration energy system
MDETESmarine diesel engine trigeneration energy system
MPmedium pressure
MVPmeat in vacuum packs
OSSCother ship’s steam consumers
Ppropeller
RAFrecirculation air fan
RFregulation flap
SEGshaft electric generator
SMCRservice maximum continuous rating
WHRSGwaste heat recovery steam generator
t.c.transport cycle
Subscripts
aaccommodation, atmospheric air
acair conditioned
ccargo space
cacargo space air
cecooling effect
cbpcooling brine pump
cfwpcooling freshwater pump
cocogeneration
cwcooling water
ddirect solar irradiation
difdiffusely solar irradiation
dsducting system
dedirect extraterrestrial irradiation
exheat exchanger
fafresh air
fcfuel consumption
fwcsfreshwater conducting system
hlheating load
lplive product
rreflected irradiation
rarecirculating air
ssaturation, sailing
scsea current
sdseawater on the sea depth
seevaporator cooling steam
sisolar irradiation
swseawater, seawaves
wwind
AAammonia ARU
ABnavigation from A to B
ACair cooler
AHair heather
BAnavigation from B to A
CEcooling evaporator, total cooling effect
CSconducting system
CRPcompressor refrigeration plant
HLtotal heating load
FSfuel saving
FSGfired steam generator
LBAlithium-bromide ARU
LPlow pressure
Mmixture
MPmedium pressure
SLTship lifetime
TCtrigeneration concept
Yyear, timing
Superscripts
heat exchanger inlet
heat exchanger outlet
~physical items related on the low-pressure system
mean value

References

  1. Gospić, I.; Donjerković, P. The Application of the Tri-generation Systems on Offshore Vessels. In Proceedings of the 4th International Symposium Power and Process Plant, Dubrovnik, Croatia, 22 May 2000. [Google Scholar]
  2. Improving Efficiency. Available online: https://www.wartsila.com/sustainability/innovating-for-sustainability/improving-efficiency (accessed on 3 May 2021).
  3. Casisi, M.; Pinamonti, P.; Reini, M. Increasing the Energy Efficiency of an Internal Combustion Engine for Ship Propulsion with Bottom ORCs. Appl. Sci. 2020, 10, 6919. [Google Scholar] [CrossRef]
  4. Low-Speed Engines 2021. Available online: https://www.wingd.com/en/documents/general/brochures/wingd-low-speed-engines-booklet-2021/ (accessed on 3 May 2021).
  5. Efficiency of MAN B&W Two-Stroke Engines. Available online: https://man-es.com/docs/default-source/energy-storage/5510-0200-02_man-bw-two-stroke-engines_v4_lt_web.pdf?sfvrsn=c0d11a61_10 (accessed on 3 May 2021).
  6. Zamiatina, N. Comparative Overview of Marine Fuel Quality on Diesel Engine Operation. Procedia Eng. 2016, 134, 157–164. [Google Scholar] [CrossRef] [Green Version]
  7. Tana, R.; Durub, O. Prapisala Thepsithard, Assessment of relative fuel cost for dual fuel marine engines along major Asian container shipping routes. Transp. Res. E 2020, 140, 102004. [Google Scholar] [CrossRef]
  8. Moussawi, H.A.; Fardouna, F. Hasna Louahliab, Selection based on differences between cogeneration and trigeneration in. various prime mover technologies. Renew. Sustain. Energy Rev. 2017, 74, 491–511. [Google Scholar] [CrossRef]
  9. Al-Sulaiman, F.A.; Hamdullahpur, F.; Dincer, I. Trigeneration: A comprehensive review based on prime movers. Int. J. Energy Res. 2011, 35, 233–258. [Google Scholar] [CrossRef]
  10. Yao, Z.M.; Qian, Z.Q.; Li, R.E.H. Energy efficiency analysis of marine high-powered medium-speed diesel engine base on energy balance and exergy. Energy J. 2019, 176, 991–1006. [Google Scholar] [CrossRef]
  11. Dere, C.; Deniz, C. Effect analysis on energy efficiency enhancement of controlled cylinder liner temperatures in marine diesel engines with model based approach. Energy Convers. Manag. 2020, 220, 113015. [Google Scholar] [CrossRef]
  12. Cavalcanti, E.J.C. Energy, exergy and exergoenvironmental analyses on gas-diesel fuel marine engine used for trigeneration system. Appl. Therm. Eng. 2021, 184, 116211. [Google Scholar] [CrossRef]
  13. Marty, P.; Hétet, J.-F.; Chalet, D.; Corrignan, P. Exergy Analysis of Complex Ship Energy Systems. Entropy 2016, 18, 127. [Google Scholar] [CrossRef] [Green Version]
  14. Benvenuto, G.; Trucco, A.; Campora, U. Optimization of waste heat recovery from the exhaust gas of marine diesel engines. Proceedings of the Institution of Mechanical Engineers. Part M J. Eng. Marit. Environ. 2016, 230, 83–94. [Google Scholar] [CrossRef]
  15. Altosole, M.; Benvenuto, G.; Zaccone, R.; Campora, U. Comparison of Saturated and Superheated Steam Plants for Waste-Heat Recovery of Dual-Fuel Marine Engines. Energies 2020, 13, 985. [Google Scholar] [CrossRef] [Green Version]
  16. Theotokatos, G.; Livanos, G. Techno-economical analysis of single pressure exhaust gas waste heat recovery systems in marine propulsion plants. Proceedings of the Institution of Mechanical Engineers. Part M J. Eng. Marit. Environ. 2013, 227, 83–97. [Google Scholar]
  17. Baldi, F.; Ahlgren, F.; Nguyen, T.-V.; Thern, M.; Andersson, K. Energy and Exergy Analysis of a Cruise Ship. Energies 2018, 11, 2508. [Google Scholar] [CrossRef] [Green Version]
  18. Baldi, F.; Johnson, H.; Gabrielii, C.; Andersson, K. Energy and exergy analysis of ship energy systems—The case study of a chemical tanker. Int. J. Thermodyn. 2015, 18, 82–93. [Google Scholar] [CrossRef] [Green Version]
  19. Du, S.; Wang, R.Z.; Chen, X. Development and experimental study of an ammonia water absorption refrigeration prototype driven by diesel engine exhaust heat. Energy 2017, 130, 420–432. [Google Scholar] [CrossRef]
  20. Ren, J.; Qian, Z.; Yao, Z.; Gan, N.; Zhang, Y. Thermodynamic Evaluation of LiCl-H2O and LiBr-H2O Absorption Refrigeration Systems Based on a Novel Model and Algorithm. Energies 2019, 12, 3037. [Google Scholar] [CrossRef] [Green Version]
  21. Oil Monster. Available online: https://www.oilmonster.com/ (accessed on 3 May 2021).
  22. BIX Bunker Index. Available online: https://www.bunkerindex.com/ (accessed on 3 May 2021).
  23. Vallianou, V.A.; Frangopoulos, C.A. Dynamic Operation Optimization of a Trigeneration System. Int. J. Thermodyn. 2012, 15, 239–247. [Google Scholar] [CrossRef] [Green Version]
  24. Glavan, I.; Prelec, Z. The analysis of trigeneration energy systems and selection of the best option based on criteria of ghg emission, cost and efficiency. Eng. Rev. 2012, 32, 131–139. [Google Scholar]
  25. Frangopoulos, C.A. Developments, Trends, and Challenges in Optimization of Ship Energy Systems. Appl. Sci. 2020, 10, 4639. [Google Scholar] [CrossRef]
  26. Gospić, I.; Boras, I.; Mravak, Z. Low-Temperature Ship Operations, Brodogradnja: Teorija i Praksa Brodogradnje i Pomorske Tehnike. 2011, Volume 62. Available online: https://hrcak.srce.hr/75662 (accessed on 3 May 2021).
  27. Lugo-Villalba, R.A.; Guerra, M.Á.; López, B.S. Calculation of marine air conditioning systems based on energy savings. Ciencia Y Tecnología Buques 2017, 11, 103–117. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, S.; Vargas, J.; Ordonez, J. Ship HVAC System Analysis and Optimization Tool. In Proceedings of the 2019 IEEE Electric Ship Technologies Symposium, Arlington, VA, USA, 14–16 August 2019. [Google Scholar] [CrossRef]
  29. Mărcuș, G.; Lungu, C.I. Partial load efficiency analysis of a CCHP plant with RICE and H2O-LiBr absorption chiller. In Proceedings of the CLIMA 2019 Congress, Bucharest, Romania, 26–29 May 2019. [Google Scholar] [CrossRef]
  30. Han, B.; Li, W.; Li, M.; Liu, L.; Song, J. Study on Libr/H2O absorption cooling system based on enhanced geothermal system for data center. Energy Rep. 2020, 6, 1090–1098. [Google Scholar] [CrossRef]
  31. Lima, A.A.S.; Leite, G.d.N.P.; Ochoa, A.A.V.; Santos, C.A.C.d.; Costa, J.A.P.d.; Michima, P.S.A.; Caldas, A.M.A. Absorption Refrigeration Systems Based on Ammonia as Refrigerant Using Different Absorbents: Review and Applications. Energies 2021, 14, 48. [Google Scholar] [CrossRef]
  32. Ouadha, A.; El-Gotni, Y. Integration of an ammonia-water absorption refrigeration system with a marine Diesel engine: A thermodynamic study. Proc. Comput. Sci. 2013, 19, 754–761. [Google Scholar] [CrossRef] [Green Version]
  33. IMO. Prevention of Air Pollution from Ships. International Maritime Organization (IMO). 2017. Available online: https://www.imo.org/en/OurWork/Environment/Pages/Air-Pollution.aspx (accessed on 3 May 2021).
  34. Ammar, N.R.; Seddiek, I.S. Eco-environmental analysis of ship emission control methods: Case study RO-RO cargo vessel. Ocean Eng. 2017, 137, 166–173. [Google Scholar] [CrossRef]
  35. Food Handling, Storage and Preparation. Maritime Labour Convention, 2006 (MLC, 2006). Available online: https://webcache.googleusercontent.com/search?q=cache:Nbx7BJ33CHIJ:https://www.ilo.org/dyn/normlex/en/f%3Fp%3D1000:53:::NO:53:P53_FILE_ID:3132751+&cd=1&hl=en&ct=clnk&gl=kr (accessed on 3 May 2021).
  36. Panama Maritime Authority Department of Seafarers, Instructions for Handling, Storage and Preparation of Food Under the Maritime Labour Convention, 2006 (Rule 3.2 Rule 3.2). Available online: http://webcache.googleusercontent.com/search?q=cache:udhN1MNh1ggJ:www.ilo.org/dyn/normlex/en/f%3Fp%3D1000:53::::53:P53_FILE_ID:3130264+&cd=1&hl=en&ct=clnk&gl=kr (accessed on 3 May 2021).
  37. Lemmon, E.W.; Bell, I.H.; Huber, M.L.; McLinden, M.O. REFerence fluid PROPerties 10, Applied Chemicals and Materials Division National Institute of Standards and Technology: Boulder, CO, USA.
  38. WMO-World Meteorological Organization. Available online: https://public.wmo.int/en (accessed on 3 May 2021).
  39. Martić, I.; Degiuli, N.; Farkas, A.; Gospić, I. Evaluation of the Effect of Container Ship Characteristics on Added Resistance in Waves. J. Mar. Sci. Eng. 2020, 8, 696. [Google Scholar] [CrossRef]
  40. Farkas, A.; Degiuli, N.; Martić, I. Assessment of hydrodynamic characteristics of a full-scale ship at different draughts. Ocean Eng. 2018, 156, 135–152. [Google Scholar] [CrossRef]
  41. Thomas, A.P.; Thekaekara, M.P. Experimental and Theoretical Studies on Solar Energy for Energy Conversion. In Proceedings of the Joint Conference, Winnipeg, MB, Canada, 15–20 August 1976. [Google Scholar]
  42. Kreider, J.F. Rabi A: Heating and Cooling of Buildings, Design for Efficiency; McGraw-Hill, Inc.: New York, NY, USA, 1994. [Google Scholar]
  43. Gospić, I. Modeliranje Brodskih Dizelmotornih Trigeneracijskih Energetskih Sustava. Ph.D. Thesis, University of Zagreb, Zagreb, Croatia, 2008. [Google Scholar]
  44. Global Average Bunker Price. Available online: https://shipandbunker.com/prices/av/global/av-glb-global-average-bunker-price#MGO (accessed on 7 June 2021).
  45. Summerheat. EU Intelligent Energy Europe Programme EIE-06-194. April 2009. Available online: https://www.euroheat.org/wp-content/uploads/2016/04/SUMMERHEAT_Report.pdf (accessed on 7 June 2021).
Figure 1. The trigeneration scheme.
Figure 1. The trigeneration scheme.
Jmse 09 00773 g001
Figure 2. The navigation route through realistic surroundings.
Figure 2. The navigation route through realistic surroundings.
Jmse 09 00773 g002
Figure 3. Air temperature models en route: (a) ϑ a Y φ r , t Y , (b) Δ ϑ a m φ r , t Y , (c) Δ ϑ a d φ r , t Y , (d) ϑ a d ( φ r , t Y ) .
Figure 3. Air temperature models en route: (a) ϑ a Y φ r , t Y , (b) Δ ϑ a m φ r , t Y , (c) Δ ϑ a d φ r , t Y , (d) ϑ a d ( φ r , t Y ) .
Jmse 09 00773 g003
Figure 4. (a) wind velocity modules v w φ r , t Y and (b) significant wave high H s φ r , t Y en route.
Figure 4. (a) wind velocity modules v w φ r , t Y and (b) significant wave high H s φ r , t Y en route.
Jmse 09 00773 g004
Figure 5. Unsteady balanced items of MDE for navigation in 1st and 2nd transport cycle.
Figure 5. Unsteady balanced items of MDE for navigation in 1st and 2nd transport cycle.
Jmse 09 00773 g005
Figure 6. (a) Air ϑ a t and (b) sea temperature ϑ s d t through one transport cycle.
Figure 6. (a) Air ϑ a t and (b) sea temperature ϑ s d t through one transport cycle.
Jmse 09 00773 g006
Figure 7. Unsteady cogeneration effects for navigation in first and second transport cycle.
Figure 7. Unsteady cogeneration effects for navigation in first and second transport cycle.
Jmse 09 00773 g007
Figure 8. Quasi-static characteristics (a) C O P A A   p C E , ϑ c s and (b) C O P L B A   p C E , ϑ c s .
Figure 8. Quasi-static characteristics (a) C O P A A   p C E , ϑ c s and (b) C O P L B A   p C E , ϑ c s .
Jmse 09 00773 g008
Figure 9. The air conditioning processes for the ship’s spaces in a Mollier h-x diagram.
Figure 9. The air conditioning processes for the ship’s spaces in a Mollier h-x diagram.
Jmse 09 00773 g009
Figure 10. Simplified quasi-static thermal network of the accommodation.
Figure 10. Simplified quasi-static thermal network of the accommodation.
Jmse 09 00773 g010
Figure 11. Main parameters of the solar geometry on the sailing route; (a) for tangential sphere surface and (b) for arbitrary place kth ship’s flat surface.
Figure 11. Main parameters of the solar geometry on the sailing route; (a) for tangential sphere surface and (b) for arbitrary place kth ship’s flat surface.
Jmse 09 00773 g011
Figure 12. Equivalent attenuation of the irradiation by the cloudiness: (a) Ψ c ( φ r , t Y ) and (b) Ψ c ( t ) 1 . t c .
Figure 12. Equivalent attenuation of the irradiation by the cloudiness: (a) Ψ c ( φ r , t Y ) and (b) Ψ c ( t ) 1 . t c .
Jmse 09 00773 g012
Figure 13. Overall solar irradiation for the ship’s accommodation enclosures during navigation through first transport cycle.
Figure 13. Overall solar irradiation for the ship’s accommodation enclosures during navigation through first transport cycle.
Jmse 09 00773 g013
Figure 14. Characteristics of heat loads for ship’s accommodation during navigation: (a) from A to B through 1st transport cycle and (b) during navigation from B to A through 10th transport cycle.
Figure 14. Characteristics of heat loads for ship’s accommodation during navigation: (a) from A to B through 1st transport cycle and (b) during navigation from B to A through 10th transport cycle.
Jmse 09 00773 g014
Figure 15. Quasi-static balanced temperatures for (a) conditioned air and cooling water in AC and (b) cooling water and cooling steam in evaporator of LBARU-a.
Figure 15. Quasi-static balanced temperatures for (a) conditioned air and cooling water in AC and (b) cooling water and cooling steam in evaporator of LBARU-a.
Jmse 09 00773 g015
Figure 16. The characteristics of quasi-static balanced temperatures of involved process’s media in the ship’s accommodation climatization system; (a) for 1st transport cycle and (b) for the 4th transport cycle.
Figure 16. The characteristics of quasi-static balanced temperatures of involved process’s media in the ship’s accommodation climatization system; (a) for 1st transport cycle and (b) for the 4th transport cycle.
Jmse 09 00773 g016
Figure 17. MP heating load Φ h l A B ( t ) and heating effect of MP cogeneration Φ c o A B ( t ) for various commodities during 1st and 12th transport cycle for navigation (a) from A to B, and (b) B to A.
Figure 17. MP heating load Φ h l A B ( t ) and heating effect of MP cogeneration Φ c o A B ( t ) for various commodities during 1st and 12th transport cycle for navigation (a) from A to B, and (b) B to A.
Jmse 09 00773 g017
Figure 18. Cogeneration effect Φ c o A B ( t ) for different evaporator area and heating loads Φ h l A B ( t ) C for different citrus storage microclimatic parameters φ a and g f a , during 10th transport cycle.
Figure 18. Cogeneration effect Φ c o A B ( t ) for different evaporator area and heating loads Φ h l A B ( t ) C for different citrus storage microclimatic parameters φ a and g f a , during 10th transport cycle.
Jmse 09 00773 g018
Figure 19. Heating effect of LP cogeneration Φ ˜ c o A B ( t ) and LP heating load Φ ˜ h l A B ( t ) for various commodities during 1st and 11th transport cycle for navigation from A to B for (a) 1st and (b) 11th transport cycle.
Figure 19. Heating effect of LP cogeneration Φ ˜ c o A B ( t ) and LP heating load Φ ˜ h l A B ( t ) for various commodities during 1st and 11th transport cycle for navigation from A to B for (a) 1st and (b) 11th transport cycle.
Jmse 09 00773 g019
Figure 20. Unsteady cogeneration heating effect and total unsteady low-pressure heating loads for products B, C, DF and FM during the ballast voyage, from B to A for 10th transport cycle.
Figure 20. Unsteady cogeneration heating effect and total unsteady low-pressure heating loads for products B, C, DF and FM during the ballast voyage, from B to A for 10th transport cycle.
Jmse 09 00773 g020
Figure 21. Fuel savings (mass flow kg/h) in the CACS for products B, C, DF and FM during navigation of unloaded reefer from B to A for 10th transport cycles.
Figure 21. Fuel savings (mass flow kg/h) in the CACS for products B, C, DF and FM during navigation of unloaded reefer from B to A for 10th transport cycles.
Jmse 09 00773 g021
Figure 22. Fuel savings in air conditioning systems of the ship’s accommodation for 1th and 10th transport cycles during navigation of loaded ship from A to B for; (a) 1st and (b) 10th transport cycle.
Figure 22. Fuel savings in air conditioning systems of the ship’s accommodation for 1th and 10th transport cycles during navigation of loaded ship from A to B for; (a) 1st and (b) 10th transport cycle.
Jmse 09 00773 g022
Figure 23. Fuel savings in maintenance of the prescribed microclimate in the storage space for various products during navigation of loaded reefer from A to B for; (a) 4th and (b) 10th transport cycles.
Figure 23. Fuel savings in maintenance of the prescribed microclimate in the storage space for various products during navigation of loaded reefer from A to B for; (a) 4th and (b) 10th transport cycles.
Jmse 09 00773 g023
Figure 24. Estimated economic benefits for the characteristic combination of transport over the 25-year ship’s lifetime.
Figure 24. Estimated economic benefits for the characteristic combination of transport over the 25-year ship’s lifetime.
Jmse 09 00773 g024
Table 1. The model ship principal particulars.
Table 1. The model ship principal particulars.
Principal Particulars
Length, overall162.3 m
Length, between perpendiculars150 m
Breadth, middle23.4 m
Depth, middle13.2 m
Draught, scantling middle9 m
Draught, ballast condition6.3 m
Deadweight13,390 tones
Hold Capacity16,999.1 m3
Main Engine MAN B&W S60MC-C
Output MCR13.56 MW at 105 min−1
Service speed17 knots
Table 2. Fuel savings for a combination laden voyage from A to B and a ballast voyage (BL) from B to A, where: t.c.—transport cycle, B—bananas, C—citruses, DC—deciduous fruits, FM—frozen meat.
Table 2. Fuel savings for a combination laden voyage from A to B and a ballast voyage (BL) from B to A, where: t.c.—transport cycle, B—bananas, C—citruses, DC—deciduous fruits, FM—frozen meat.
t.c.B + BLC + BLDF + BLFM + BL
1. t. c. A-B + B-A19.406 + 0.85528.298 + 1.93231.373 + 2.8356.500 + 5.977
2. t. c. A-B + B-A16.275 + 0.87529.400 + 1.98831.923 + 2.9176.606 + 6.035
3. t. c. A-B + B-A17.912 + 0.92530.600 + 2.00533.410 + 2.9227.021 + 6.340
4. t. c. A-B + B-A17.545 + 0.96228.237 + 2.03731.840 + 2.9426.615 + 6.062
5. t. c. A-B + B-A17.418 + 1.21227.33 + 2.23730.949 + 3.2426.609 + 6.662
6. t. c. A-B + B-A21.299 + 1.45431.513 + 2.65534.052 + 3.4826.797 + 6.934
7. t. c. A-B + B-A24.680 + 1.69533.776 + 2.85536.721 + 3.8617.125 + 7.132
8. t. c. A-B + B-A23.225 + 1.70533.507 + 2.91536.306 + 3.9027.222 + 7.144
9. t. c. A-B + B-A23.449 + 1.69533.740 + 2.89536.369 + 3.8757.309 + 7.115
10. t. c. A-B + B-A23.459 + 1.66333.781 + 2.82536.920 + 3.8307.555 + 7.109
11. t. c. A-B + B-A22.979 + 0.97233.223 + 2.10736.528 + 2.9927.486 + 6.102
12. t. c. A-B + B-A23.065 + 0.91234.043 + 2.01137.744 + 2.9007.795 + 6.058
13. t. c. A-B + B-A19.778 + 0.88530.785 + 1.99633.561 + 2.9116.916 + 6.095
14. t. c. A-B + B-A19.350 + 0.86528.220 + 1.94431.651 + 2.8566.702 + 5.985
Σ = 306.515 tonesΣ = 468.855 tonesΣ = 524.814 tonesΣ = 189.068 tones
TSLT = 25 years7662.875 tones11,721.34 tones1312.35 tones4736.7 tones
p = 500 US$/ton
TSLT = 25 years3.832 mil. US$5.861 mil. US$6.560 mil. US$2.368 mil. US$
Table 3. Savings for various combinations of cargo in the event of total transport engagement.
Table 3. Savings for various combinations of cargo in the event of total transport engagement.
t.c.BCDFFM
1. t. c. A-B19.40628.29831.3736.500
2. t. c. A-B16.27529.40031.9236.606
3. t. c. A-B17.91230.60033.4107.021
4. t. c. A-B17.54528.23731.8406.615
5. t. c. A-B17.41827.3330.9496.609
6. t. c. A-B21.29931.51334.0526.797
7. t. c. A-B24.68033.77636.7217.125
8. t. c. A-B23.22533.50736.3067.222
9. t. c. A-B23.44933.74036.3697.309
10. t. c. A-B23.45933.78136.9207.555
11. t. c. A-B22.97933.22336.5287.486
12. t. c. A-B23.06534.04337.7447.795
13. t. c. A-B19.77830.78533.5616.916
14. t. c. A-B19.35028.22031.6516.702
Σ = 289.840 tones Σ = 436.453 tones Σ = 479.340 tonesΣ = 98.318 tones
p = 500 US$/ton 1 yearTLP 25 year
Combination 1B + C346.146 US$ 9.079 mil. US$Average savings AS = 8.15 mil. US$
Combination 2B + DF384.594 US$9.615 mil. US$
Combination 3B + FM194.079 US$4.852 mil. US$
Combination 4C + DF457.901 US$11.448 mil. US$
Combination 5C + FM267.386 US$6.685 mil. US$
Combination 6FM + DF288.832 US$7.221 mil. US$
Table 4. Emissions reduction for various combinations of cargo in the event of total transport engagement.
Table 4. Emissions reduction for various combinations of cargo in the event of total transport engagement.
Transport CombinationmF [tones]ΔmCO2 [tones]ΔmSO2 [tones]
1. B + BL306.51523,601.65576.629
2. C + BL468.85536,101.835117.214
3. DF + BL524.81440,410.678131.204
4. FM + BL189.06814,558.23647.267
5. B + C726.29355,924.561181.573
6. B + DF769.18759,227.399192.297
7. B + FM.388.15829,888.16697.039
8. C + DF915.80070,516.600228.950
9. C + FM534.77141,177.367133.693
10. DF + FM577.66544,480.205144.416
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gospić, I.; Glavan, I.; Poljak, I.; Mrzljak, V. Energy, Economic and Environmental Effects of the Marine Diesel Engine Trigeneration Energy Systems. J. Mar. Sci. Eng. 2021, 9, 773. https://doi.org/10.3390/jmse9070773

AMA Style

Gospić I, Glavan I, Poljak I, Mrzljak V. Energy, Economic and Environmental Effects of the Marine Diesel Engine Trigeneration Energy Systems. Journal of Marine Science and Engineering. 2021; 9(7):773. https://doi.org/10.3390/jmse9070773

Chicago/Turabian Style

Gospić, Ivan, Ivica Glavan, Igor Poljak, and Vedran Mrzljak. 2021. "Energy, Economic and Environmental Effects of the Marine Diesel Engine Trigeneration Energy Systems" Journal of Marine Science and Engineering 9, no. 7: 773. https://doi.org/10.3390/jmse9070773

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