Next Article in Journal
Global-Local Heat Demand Development for the Energy Transition Time Frame Up to 2050
Next Article in Special Issue
Identifiability Analysis of Degradation Model Parameters from Transient CO2 Release in Low-Temperature PEM Fuel Cell under Various AST Protocols
Previous Article in Journal
Pore-Scale Modeling of Air–Water Two Phase Flow and Oxygen Transport in Gas Diffusion Layer of Proton Exchange Membrane Fuel Cell
Previous Article in Special Issue
Model-Predictive-Control-Based Reference Governor for Fuel Cells in Automotive Application Compared with Performance from a Real Vehicle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Health-Conscious Optimization of Long-Term Operation for Hybrid PEMFC Ship Propulsion Systems

Department of Engineering and Architecture, University of Trieste, via Valerio 10, 34127 Trieste, Italy
*
Author to whom correspondence should be addressed.
Energies 2021, 14(13), 3813; https://doi.org/10.3390/en14133813
Submission received: 27 May 2021 / Revised: 21 June 2021 / Accepted: 23 June 2021 / Published: 24 June 2021

Abstract

:
The need to decarbonize the shipping sector is leading to a growing interest in fuel cell-based propulsion systems. While Polymer Electrolyte Membrane Fuel Cells (PEMFC) represent one of the most promising and mature technologies for onboard implementation, they are still prone to remarkable degradation. The same problem is also affecting Lithium-ion batteries (LIB), which are usually coupled with PEMFC in hybrid powertrains. By including the combined degradation effects in an optimization strategy, the best compromise between costs and PEMFC/LIB lifetime could be determined. However, this is still a challenging yet crucial aspect, rarely addressed in the literature and rarely yet explored. To fill this gap, a health-conscious optimization is here proposed for the long-term minimization of costs and PEMFC/LIB degradation. Results show that a holistic multi-objective optimization allows a 185% increase of PEMFC/LIB lifetime with respect to a fuel-consumption-minimization-only approach. With the progressive ageing of PEMFC/LIB, the hybrid propulsion system modifies the energy management strategy to limit the increase of the daily operation cost. Comparing the optimization results at the beginning and the end of the plant lifetime, daily operation costs are increased by 73% and hydrogen consumption by 29%. The proposed methodology is believed to be a useful tool, able to give insights into the effective costs involved in the long-term operation of this new type of propulsion system.

1. Introduction

In the last few decades, an increasing number of restrictions has been introduced in the transportation sector in order to reduce greenhouse gases and pollutant emissions. This has posed new technological challenges, and the need to explore new solutions is ever more urgent. In particular, growing efforts have been placed in reducing the emissions of the shipping sector, where the new and upcoming restrictions require a thorough revision of the whole ship building process. The urge to develop and upscale new technological solutions for decarbonizing the maritime transportation sector has promoted the research and development of different strategies. Among them, fuel cell-based propulsion has gained attention, as was demonstrated by several national and international initiatives on this topic [1,2,3,4]. As emerged from a previous study by the authors [3], one of the most promising solutions employs hydrogen-fuelled Polymer Electrolyte Membrane Fuel Cells (PEMFC), hybridized with a Battery Energy Storage System (BESS). In fact, it is well known that hydrogen-fueled PEMFC can potentially guarantee zero-emission propulsion: this feature makes their use fully compliant even to the most stringent regulations. As has already been demonstrated in several analyses, a small to medium sized vessel could benefit from the installation of hybrid PEMFC/BESS propulsion systems. In fact, the technical and regulatory barriers are currently limiting the installed power of PEMFC to few hundred kW, as highlighted in Table 1 and extracted from [3]. As for BESS, Table 1 also points out that the most recent projects use Lithium-Ion Batteries (LIB), as LIB are the commercial product with the best characteristics in terms of specific energy [5].
The main drawback of such hybrid PEMFC/LIB powertrains is that both PEMFC and BESS are prone to degradation, which is responsible for the overall reduction of the plant lifespan. Hence, a correct design and management of hybrid PEMFC/LIB powertrains is fundamental in order to achieve their cost effectiveness [27,28]. The management of hybrid PEMFC/LIB powertrains also includes the energy-efficient power allocation, usually addressed by the definition of an Energy Management Strategy (EMS). The latter can be defined following different approaches, among which rule-based and optimization-based ones are the most investigated in the literature on ship propulsion systems.
As for rule-based EMSs, an example is provided by Taccani et al. [29], who analyzed the implementation of a hybrid PEMFC/LIB power system for the propulsion of a small-sized ferry. The authors proposed a rule-based EMS in which PEMFC work at a constant load in order to cover the average power demand of the vessel, while LIB are used for peak shaving. Moreover, Han et al. [30] proposed a state-based EMS to control the energy flows between the PEMFC and the batteries of the hybrid ship Alsterwasser depending on the vessel power demand and the battery State of Charge (SOC) [23,24,25]. The same vessel was also studied by Bassam et al. [31], who proposed an improved state-based EMS for the hybrid PEMFC/battery system taking into account the PEMFC efficiency as input to the proportional-integral EMS. By adopting the proposed EMS, up to 3.5% of hydrogen daily consumption could be saved in comparison with a previously proposed strategy [30]. The main advantages of rule-based EMSs lie in their simplicity of use, which makes them practical for onboard applications. As a drawback, their operative principle is based on human expertise, and they could easily reflect human errors generated while defining the control strategy.
To overcome this issue, some studies proposed optimization-based EMSs [32,33,34,35,36,37,38,39]. For example, Rivarolo et al. [32,33] developed two computational tools dedicated to the preliminary design and the thermo-economic analysis of different engineering solutions for ship power systems, including PEMFC and Solid Oxide Fuel Cells (SOFC). The proposed methodology provided the best energy system configuration in terms of environmental impact, cost, and overall dimensions. Bassam et al. [34] defined a multi-scheme EMS that minimizes the equivalent fuel consumption over one operation day of the hybrid PEMFC/BESS passenger ship Alsterwasser [23,24,25]. The developed EMS resulted in a 16.7% reduction of the total hydrogen consumption with respect to previously proposed EMSs. Zhang et al. [35] designed a real-time optimization EMS that took into account the degradation of PEMFC and used a high-pass filter to prolong battery life. A similar approach was followed by Chen et al. [36], who employed a whale optimization algorithm to optimize size and frequency for a hybrid battery/supercapacitor fuel cell ferry. Wu et al. [37] included the calculation of PEMFC degradation in the optimization of a hybrid plug-in PEMFC/BESS ferry for coastal navigation, defining the optimal size of the system by considering both greenhouse gas emissions and economic performances. For the same ferry, Wu et al. [38] recently proposed a reinforcement learning approach for the definition of an adaptive EMS. In both [37,38], the BESS degradation rate was set as a constant in each time step, independently from the operating conditions. Pivetta et al. [39] proposed a multi-objective design and operation optimization of three different ferries for coastal navigation, adopting a Mixed-Integer Linear Programming (MILP) approach based on a simultaneous minimization of PEMFC degradation, investment cost and operation cost. In this case, battery degradation was set as a constraint of the optimization model. As a main outcome, it was found that in order to limit both battery and PEMFC degradation, the battery capacity should be increased by up to 136%.
Although EMSs taking into account system degradation (also known as health-conscious EMSs) are of great interest for the shipping sector, the interplay between PEMFC and BESS degradation has been rarely addressed as of yet. Interesting contributions on this point can be found in the recent literature on hybrid PEMFC/BESS powertrains for other types of applications [40,41,42,43,44,45,46,47]. Some remarkable examples have been provided by Ҫinar et al. [43], Liang et al. [44], and Wang et al. [45]. The former group [43] used meta-heuristic algorithms to optimize hydrogen consumption in the aircraft sector. The second one [44] compared two EMSs (composite fuzzy control versus Pontryagin’s minimum principle) for hybrid PEMFC/LIB logistics, while the third one [45] proposed a health-conscious EMS for a hybrid PEMFC/LIB bus. Both PEMFC and LIB degradation were included as the objectives in a cost-optimization problem. In addition, Li et al. [46] studied an online adaptive EMS for a PEMFC hybrid vehicle, which takes into account the degradation of PEMFC, batteries, and supercapacitors. The authors concluded that the degradation of one of the power sources implied an increase in the degradation of the other ones. Moreover, they proved that the proposed EMS guaranteed normal plant operation also in case of plant source degradation. However, a holistic analysis throughout the whole plant lifetime was not addressed, and a quantitative plant lifetime estimation was not provided. Lastly, a comprehensive review on health-conscious EMSs for hybrid fuel cell vehicles was proposed by Yue et al. [47], who concluded that future trends in this research field should address the multi-objective optimization of a hybrid PEMFC/LIB system including degradation modeling and lifetime estimation methods of the power sources, possibly limiting the overall complexity of the problem.
Overall, the literature review pointed out the importance of considering the degradation of both PEMFC and LIB in optimizing the entire hybrid ship-propulsion system. However, it emerged that a holistic optimization procedure that takes into account the progressive worsening of both PEMFC and BESS performances with the progressive ageing over the plant lifetime has not been addressed yet. This indeed would be fundamental for a complete understanding of the costs involved in the operation of hybrid PEMFC/LIB propulsion systems. To fill this gap, the present study proposes an innovative methodology that (i) evaluates the degradation over time of both PEMFC and LIB, (ii) optimizes the operation of a hybrid PEMFC/LIB powertrain over the entire lifetime, and (iii) simultaneously minimizes the fuel consumption and the power source degradation rates. As is widely recognized, the MILP approach is a good candidate for optimizing complex energy systems because it allows a reduction in the computation effort with respect to other optimization techniques [48,49,50,51]. Therefore, in this study, MILP equations have been used to describe the hybrid PEMFC/LIB operation and a multi-objective optimization has been performed.
A small-size ferry for coastal navigation has been taken as case study to perform the optimization of operation over the entire power plant lifespan. The design of the system previously optimized in [39] has been taken as a reference for the operation optimization here proposed. Here, Lithium Iron Phosphate (LFP) chemistry has been chosen as LIB. Indeed, among different LIB chemistries Nickel Manganese Cobalt (NMC) and LFP batteries are addressed as the most mature technologies for shipping, but LFP are slightly better in terms of safety onboard [5].
The first part of the paper describes the methodology proposed for the analysis. Afterwards, the case study is described, and the results are presented and discussed.

2. Methodology

This section shows the main characteristics of the proposed method, with a particular focus on the equations describing the energy conversion and storage units. Given the complexity of the multi-objective optimization, a MILP approach has been adopted to limit the overall computational effort. Therefore, operation and performance degradation of the energy units, i.e., LFP battery and PEMFC, have been described with MILP equations. In the next paragraphs, a description of the procedure followed in the analysis is reported.
At first, the preliminary sizing of the hybrid PEMFC/LFP battery propulsion system has been determined by implementing the design and operation optimization model developed in [39]. The obtained sizes for PEMFC and LFP battery systems have been taken as input parameters of the new optimization approach for the long-term operation.
Starting from time t i = 0, corresponding to a novel installation of both PEMFC and LFP battery, consecutive optimizations have been run in a loop until either PEMFC or LFP battery reached the end of their life. In order to limit the computational effort, it has been assumed that the first day of the month i is representative for the operation of all the days of the month (assuming all months have 30 days). Operation optimization has been performed for the day representative of the entire month i. The cumulative degradation of PEMFC and LIB calculated at the end of month i has been taken as an input parameter for the operation optimization of month i + 1.
The optimization problem has been developed in Python [52] and solved with the MILP optimizer Gurobi Optimization [53]. The general formulation of the MILP optimization problem has been set to find the optimal value of the generic continuous variables x * ( t ) and binary variables δ * ( t ) that maximize or minimize the set objective function(s) Z (Equation (1)). The MILP problem is subject to equality constraints g ( t ) and inequality constraints h ( t ) (Equations (2) and (3)), which include the description of the modeled ship energy system. The continuous and binary variables represent the decision variables of the optimization problem. In particular, binary variables have been implemented to decide about activation/de-activation of each energy unit during operation.
Z = f ( x * ( t ) ,   δ * ( t ) )
g ( x * ( t ) ,   δ * ( t ) ) = 0
h ( x * ( t ) ,   δ * ( t ) ) 0
In the following subsections, the main constraints describing the energy system operation and the degradation models of PEMFC and LFP battery are presented, i.e., the relationships g ( t ) and h ( t ) in Equations (2) and (3). Subsequently, the objective functions of the optimization problem, i.e., the objective functions Z in Equation (1), and the solving procedure are illustrated and discussed.

2.1. Energy System Constraints

It has been assumed that at each time step t the ship energy demand P d e m a n d ( t ) (propulsion and auxiliaries power demand) is fulfilled by the power output of the n fuel cell stacks P F C j ( t ) and by the discharging power of battery P b a t t + ( t ) . P b a t t ( t ) indicates the charging power of battery, which is charged by fuel cells whenever P F C i ( t ) is higher than P d e m a n d ( t ) . The overall power balance has been set as (Equation (4)):
j = 1 n P F C j ( t ) + P b a t t + ( t ) = P d e m a n d ( t ) + P b a t t ( t )
Equations (5) and (6) are the constraints adopted to limit the proposed energy system volume and weight, set to be lower than the values of volume V m a x   and weight w m a x   of the ferry existing energy system (i.e., internal combustion engines powered by marine gasoil). A correction factor ( c o v ) has been introduced to consider the current limited development state of PEMFC and LFP battery technologies.
j = 1 n t i n t f i n ( F F C j ( t ) · d t ) · V H 2 + E b a t t e r y m a x · V b a t t + n · P F C m a x · V F C V m a x · ( 1 + c o v )
j = 1 n t i n t f i n ( F F C j ( t ) · d t ) · w H 2 + E b a t t e r y m a x · w b a t t + n · P F C m a x · w F C w m a x · ( 1 + c o v )
where t i n and t f i n identify the time interval considered for the operation of ferry without refueling. With reference to Equation (5), the volume occupied by hydrogen is expressed by t i n t f i n ( F F C i ( t ) · d t ) · V H 2 , where F F C j ( t ) is the fuel consumption of the stack j at time t and V H 2 is the volume occupied by hydrogen and its storage system. The volume occupied by LIB is expressed as E b a t t e r y m a x · V b a t t , where E b a t t e r y m a x is the energy capacity of the battery and V b a t t is the specific volume. The volume required by PEMFC is determined as n · P F C m a x · V F C , with n being the number of installed PEMFC stacks, P F C m a x the fuel cell stack rated power and V F C the specific volume of a PEMFC stack. Similarly, in Equation (6) the weight of hydrogen is expressed as t i n t f i n ( F F C j ( t ) · d t ) · w H 2 , being w H 2 the weight of bunkered hydrogen and its storage system. Weight of LIB is calculated as E b a t t e r y m a x · w b a t t , where w b a t t is the specific weight of LIB, and weight of PEMFC is n · P F C m a x · w F C , where w F C is the specific weight of PEMFC stack.

2.2. Polymer Electrolyte Membrane Fuel Cells Model

It has been considered that separate PEMFC stacks with fixed rated power are installed on board. Each stack has been assumed to operate independently from the others. Stacks with different rated power or with different performances could be included in the model by modifying the stack operation parameters (see Section 3).
The operation of the PEMFC stack j is described by the following MILP equations (Equations (7)–(11)):
F F C j ( t ) = ( k 1 F · I F C j ( t ) + k 2 F ) · δ F C j ( t ) + δ s t , u p j ( t ) · F s t a r t · P F C m a x
P F C j = ( k 1 P · I F C j ( t ) + k 2 P ) · δ F C j ( t )
I F C m i n · δ F C j ( t ) I F C j ( t ) I F C m a x · δ F C j ( t )
| P F C j ( t ) P F C j ( t 1 ) |   Δ P F C
0 δ F C j   ( t ) δ F C j   ( t + 1 ) + δ s t , u p j ( t )
where I F C j is the current density (in A/cm2), k 1 F , k 2 F ,   k 1 P and k 2 P are the linearization coefficients, δ F C j is the binary variable defining on/off status, δ s t , u p j is the binary variable defining the occurrence of a start-up phase (if δ s t , u p j = 1, it will be counted as one start/stop cycle), F s t a r t is the fuel consumption in the start-up phase, I F C m i n and I F C m a x are the limiting values for the current density, and Δ P F C is the allowed load variation.
The PEMFC stack degradation is a complex phenomenon that affects the performances of several stack components (e.g., electrodes, membranes, bipolar plats, gas diffusion layers). Degradation effects may be different from one cell to another, e.g., near the edges, cells are subject to a faster degradation [54]. Besides, at the component level, the degradation mechanisms are hard to estimate. As reported in [47], the ageing mechanisms in mobility application could be analyzed through three possible approaches: impedance estimation based on electrochemical impedance spectrometry, remaining useful life estimation, and a stack voltage degradation model. Even though the latter approach is less accurate with respect to other degradation models and strongly depends on experimental data, voltage degradation can be expressed by MILP equations, allowing a sensible reduction of the computational effort required by an optimization with respect to, for example, Mixed-Integer Quadratic Programming (MIQP).
Therefore, the stack voltage degradation model has been adopted in this study to describe the PEMFC ageing. PEMFC voltage loss has been considered to depend mainly on three operating conditions: the stack load variation, the start/stop cycle and the idling/high-current operation of the stack [55]. The voltage loss caused by the stack degradation has been expressed as the voltage reduction of a single cell at an equal current density output, as proposed in other studies available in the literature [37,56,57,58].
Performance degradation of the PEMFC stack j at time t is described by Equations (12)–(15).
d V l o a d j ( t + 1 ) = | P F C j ( t ) P F C j ( t + 1 ) | · Δ v l o a d
d V s t , u p j ( t ) = δ s t , u p j ( t ) · Δ v s t , u p
d V P F C j ( t ) = k 1 d v · I F C j ( t ) + k 2 d v
d V j ( t ) = d V l o a d j ( t ) + d V s t , u p j ( t ) + d V P F C j ( t )
where d V l o a d j ( t ) is the voltage reduction due to load variation, Δ v l o a d and Δ v s t , u p are the proportionality constants, d V s t , u p , j ( t ) is the voltage reduction due to start-up, d V P F C , j ( t ) is the voltage reduction depending on the power load, k 1 d V and k 2 d V are the linearization coefficients, d V j ( t ) is the total loss of voltage at time t.
The cumulative voltage degradation evaluated for the single cell modifies the characteristic curves of the entire stack. In particular, the relation between the cumulative loss of voltage (Equation (16)) and power has been approximated by a linear curve, which relates voltage loss to the angular coefficient of the power curve ( k 1 P ), as described in Equation (17). The SOH for PEMFC stack ( S O H P E M F C j ) has been defined as in Equation (18).
d V T O T , j = d V j ( t ) · d t
k 1 P = k 1 d e g · d V T O T , j + k 2 d e g
S O H P E M F C j = V r e f , F C d V T O T , j V r e f , F C
where d V T O T , j is the cumulative loss of voltage, k 1 d e g and k 2 d e g are the linearization coefficients, V r e f , F C is the reference maximum voltage for a novel single cell (voltage at the minimum current density I F C m i n ).
It has been assumed that k 1 F and k 2 F are not affected by the stack degradation d V T O T , j . Therefore, when the power loss caused by degradation increases, the efficiency of the stack tends to decrease. The stack efficiency η F C   has not been included as a decisional variable in the optimization problem, but it has been evaluated as a result (Equation (19)).
η F C = P F C j / F F C j

2.3. Lithium Iron Phosphate Battery Model

At each time step t, the energy stored in the battery E b a t t ( t ) has been defined according to the charging/discharging power of the battery P b a t t ± ( t )   and the relative efficiency η b a t t , assumed constant and equal in the charging/discharging phase (Equation (20)). The SOC of the battery at time t has been calculated as the ratio between E b a t t ( t ) and the energy capacity of the battery E b a t t e r y m a x (Equation (21)). The battery C-rate has been determined as the ratio between the charging/discharging power of the battery and the battery capacity (Equation (22)), and the C-rate has been constrained in order not to exceed the maximum C-rate ( C r a t e m a x ) achievable by the battery (Equation (23)). A constraint has been set on the SOC in the first-time step of the day t 0 to make it correspond to the SOC at the end of the day t f i n (Equation (24)).
E b a t t ( t ) = E b a t t ( t 1 ) + ( η b a t t · P b a t t ( t ) ( 1 / η b a t t ) ·   P b a t t + ( t ) ) · Δ t  
S O C ( t ) = E b a t t ( t ) E b a t t e r y m a x
C r a t e ± ( t ) = P b a t t ± ( t ) E b a t t e r y m a x
| C r a t e ± ( t ) | C r a t e m a x
S O C ( t 0 ) = S O C ( t f i n )
As for battery degradation, the most dominant ageing phenomenon of LIB is the formation of a Solid Electrolyte Interface (SEI) at the electrode/electrolyte interface. The thickness of SEI increases over time, leading to a progressive increase in the impedance of the battery and a consequent decrease in the battery capacity. Different degradation models are available in the literature, based either on electrochemical/equivalent circuit methods (e.g., SEI film thickness models and internal resistance models), or on empirical/semi-empirical methods (e.g., capacity fade models) [47,59]. In this study, it has been chosen to model the battery degradation over time with a capacity fade model. In fact, such an approach allows to limit the overall computational burden with respect to electrochemical models, allowing a comprehensive analysis of the entire system.
Commonly, battery ageing can be classified into two categories: calendar ageing and cycle ageing. Calendar ageing occurs when the battery is at rest condition, i.e., when no current is flowing through the battery. Calendar degradation rates depend on SOC and temperature. Cycle ageing occurs when the battery is charged or discharged, and depends on battery C-rate, SOC, temperature, number of performed equivalent cycles, and Depth of Discharge (DOD) [60,61].
In the proposed study, the hybrid propulsion system has been considered to power a small-size ferry for public transport, which operates 24 h, seven days a week. Therefore, it is reasonable to assume that only cycle ageing occurs on the battery. Nevertheless, to prevent an excessive calendar ageing of the battery in case of a prolonged stop of the vessel it has been decided to set a constraint on the SOC at the end of each day of operation (Equation (25)). The reference S O C c a l has been set equal to 60%, according to the data available in the literature [62].
S O C ( t f i n ) = S O C c a l
Cycle ageing has been evaluated following the experimental studies available in [61] for a single LFP cell. As done for PEMFC (see Section 2.2), in this case it has also been considered in terms of the degradation of a single cell, representative for the whole battery. At each time step t of an operative day, T b a t t ( t ) has been expressed as linear function of C r a t e ( t ) , through the linearization coefficients a 1 T , x and a 2 T , x that depend on the ambient temperature (Equations (26) and (27)). The capacity fade of battery Q l o s s ( t ) has been calculated as a linear function of battery temperature T b a t t ( t ) (Equation (28)).
T b a t t ( t ) = a 1 T , x · C r a t e ( t ) + a 2 T ,   x
a 1 T , x = f ( T a m b , x ) ;     a 2 T , x = f ( T a m b , x )
Q l o s s i ( t ) = a 1 Q T , i · T b a t t ( t ) + a 2 Q T , i
The time dependency of Q l o s s ( t ) in the long-term operation has been taken into account by determining the linearization coefficients a 1 Q T , x and a 2 Q T , x in Equation (27) at each month of the plant lifetime. In particular, a 1 Q T , x and a 2 Q T , x at the month i have been evaluated as a function of the energy throughput A h i at month I, as expressed in Equation (29).
a 1 Q T , i = k 11 , A h · A h i 2 + k 21 , A h · A h i   + k 31 , A h
a 2 Q T , i = k 12 , A h · A h i 2 + k 22 , A h · A h i   + k 32 , A h
where k 11 , A h , k 21 , A h , k 31 , A h , k 12 , A h , k 22 , A h , and k 32 , A h . Note that Equations (29) and (30) have been evaluated at the beginning of the optimization and hence have not been formulated as MILP. Equation (31) expresses the energy throughput A h i at month i.
A h i = n d a y s · N c y c l e i 1 · D O D · A h c e l l
where n d a y s is the number of days per month (assumed all equal to 30 days), N c y c l e i 1 is the number of equivalent cycles performed by the battery during one operation day, representative for the month i − 1. N c y c l e has been determined via a rainflow algorithm. A h c e l l indicates the capacity of the single battery cell. DOD effect on the capacity fade has not been directly taken into account in the evaluation of Q l o s s ( t ) , since its effect on the capacity fade is negligible if compared to the cycling effect [61]. As for the SOC effect, a constraint has been set (Equation (32)) to limit the SOC-window of operation and hence further limit the degradation, as proposed for example in [63].
S O C m i n   S O C ( t ) S O C m a x
At the end of each month i, the SOH of the battery ( S O H b a t t e r y i ) has been determined, as in Equation (32).
S O H b a t t e r y i = E b a t t e r y , m a x Q l o s s i E b a t t e r y , m a x

2.4. Objective Functions

A multi-objective optimization has been performed to determine the optimal long-term operation of the proposed FC/LFP battery hybrid system, taking into account the fuel consumption and the degradation of both PEMFC and LFP battery over their lifetime.
Three objective functions have been specified, namely fuel consumption ( f 1 in Equation (34)), PEMFC performance degradation ( f 2 in Equation (35)), and LFP battery performance degradation ( f 3 in Equation (36)). For the multi-objective optimization, a “Blended Objectives” method has been adopted, considering a linear combination of the objective functions, each with a fixed weight [53]. In this case, a cost-related weight has been assumed for each objective function. The weight of f 1 ( w 1 in Equation (38)) depends on the cost of hydrogen. Weights of f 2 and f 3 ( w 2 and w 3 in Equations (39) and (40)) derive from the cost of components and their lifetime. The multi-objective optimization has been performed by minimizing the linear combination of the three objective functions, defined as f M O (Equation (37)).
A 1% deviation from the optimal value of the f M O function has been allowed.
f 1 = j = 1 n 0 t f i n F F C j ( t )   d t
f 2 = j = 1 n 0 t f i n d V j ( t ) d t
f 3 = 0 t f i n Q l o s s   b a t t e r y ( t ) d t
f M O = M i n i m i z e   ( w 1 · f 1 + w 2 · f 2 + w 3 · f 3 )
w i parameters have been defined as in Equations (38)–(40).
w 1 = c H 2
w 2 = c F C · P F C m a x V o
w 3 = c b a t t · E b a t t e r y , m a x Q l o s s m a x
where c F C is the cost of a PEMFC stack (€/kW), and V o is the total voltage loss permitted for the single cell of the PEMFC stacks (set equal to 20% of V r e f , F C ). Similarly, c b a t t is the cost of LFP battery (€/kWh), and Q l o s s m a x is the maximum capacity fade of the battery (set equal to 20%).

3. Case Study

A typical small size RO-Pax (Roll-on/roll-off Passengers) ferry operating in coastal areas has been chosen as a case study. Table 2 reports the main characteristics of the ferry. The proposed RO-Pax is similar to the represented ferry in Figure 1. Power system volume and weight in Table 2 have been elaborated from data available on diesel storage on board and diesel-electric generator for marine propulsion, scaled by the installed power of the ferry.
The ferry power demands (propulsion and auxiliaries demands) over a typical operation day have been evaluated in [39]. Ferry power demands rated per minute enter as input in the proposed long-term operation optimization. More details on power profile and the typical sailing cycle of the ferry chosen as a case study are reported in Appendix A.
A simplified schematic of the propulsion system for the proposed case study is presented in Figure 2. Resulting from the design optimization performed in [39], the proposed energy system is composed of two 100-kW PEMFC stacks fueled by hydrogen, an energy management system, and a 286-kWh LFP battery. Hydrogen has been assumed to be stored in liquid form in a cryogenic tank, with the aim of reducing the volume and weight for the storage system in comparison to other types of hydrogen storage systems (e.g., compressed hydrogen tanks). It has been assumed that the ferry has enough fuel for a daily operation.
The parameters that have been used for the long-term operation optimization are reported in Table 3. These model parameters represent the main characteristics of the considered ferry and of the current market scenario. In particular, the PEMFC and battery characteristics have been evaluated on the basis of the existing technologies.

4. Results

In this section, the main results of the optimization are presented and discussed. In the first part, a comparison of the plant lifetime obtained with the proposed optimization approach and the one obtained without considering energy units degradation (i.e., non-health-conscious optimization) is presented. Afterwards, the optimal operation profile during one day of operation along the entire lifetime operation is presented and compared with the one obtained through the non-health-conscious optimization. Lastly, the increase/decrease path of different optimization variables during the plant lifetime is analyzed.

4.1. Plant Lifetime Estimation

In this subsection, an evaluation of the energy system lifetime obtained with the proposed long-term optimization approach is compared with a non-health-conscious optimization approach. The latter considers the minimization of fuel consumption only (i.e., the minimization of the objective function f 1 in Equation (34)), and will be later referred to as FM (Fuel consumption Minimization). Figure 3 shows the SOH of both LFP batteries and PEMFC during the plant lifetime, expressed in the number of operative days. The dots indicate the SOH evaluated at the month i. Figure 3a shows the results of the FM optimization, while Figure 3b illustrates the results obtained with the proposed health-conscious multi-objective optimization.
From Figure 3a it can be noticed that with an FM approach the plant lifetime results to be about 210 days. After this time, S O H P E M F C reaches 80%, while S O H b a t t e r y is still about 93%. This is due to the fact that the design of the system is fixed following the methodology explained in [39], where battery degradation was anyway included as a constraint in the methodology. Therefore, if PEMFC degradation is not minimized the battery results to be oversized. This implies low C-rates and hence lower degradation rates of LIB. On the other hand, PEMFC stacks tend to follow the power demand. Such an operation implies that PEMFC operates at a high current density and performs large load variations, causing high degradation rates of PEMFC and hence limited lifetime.
The optimization approach proposed here allows to increase the plant lifetime from 210 to 600 operation days, as shown in Figure 3b. In this case, the first energy unit to reach the threshold value of SOH is the LFP battery, while the S O H P E M F C , at this time is still about 83% (so theoretically PEMFC could work for another 60–90 days). This aspect is also due to the current higher cost of PEMFC with respect to battery (see Table 3). As a consequence, the degradation of PEMFC ( f 2 ) weighs more on the overall objective function ( f M O ), and hence the optimal EMS advantages the minimization of PEMFC degradation instead of the LFP battery one.

4.2. Health-Conscious Optimal Daily Operation over the Entire Lifetime

The proposed optimization approach allows also to determine the optimal daily operation of the ferry, representative for the entire month i, during the whole lifetime of the energy system.
Figure 4 shows the health-conscious optimal operation of the ferry during one day in the first month of operation (Figure 4a) and the last month of operation (Figure 4b). Here, at both the beginning and the end of the plant lifetime, the PEMFC tends to avoid the operation at rated power, as this would result in a higher degradation rate and a lower energy conversion efficiency, and tends also to avoid large load variations and start-up phases, which would also cause higher degradation rates. As for battery, it is possible to see that the operation changes from the beginning till the end of the plant lifetime. As mentioned before, the cost of battery is lower than that of PEMFC, hence the minimization of battery f 3 (Equation (36)) weighs less on the overall objective function than the minimization of PEMFC degradation f 2 (Equation (35)). As a consequence, at the beginning of the plant lifetime when the degradation rate of the battery is low, f 3 slightly affects f M O (Equation (37)), and battery operates to fulfill the power demand of the ferry not covered by PEMFC. Differently, at the last month of operation, the degradation rate of the battery is higher and hence f 3 weighs more on f M O . Consequently, the optimal EMS finds a compromise between the battery degradation and that of the PEMFC. For example, at the beginning of the day in the last month (Figure 4b) the PEMFC stacks operated at a higher load with respect to the first month to avoid high C-rates of the battery. To better understand the advantages of the proposed health-conscious EMS, the results of the long-term operation optimization have been compared with the results of FM optimization. Figure 5 shows the optimal operation of the ferry during one day in the first month of operation (Figure 5a) and the last month of operation (Figure 5b) if the FM approach is followed. It can be noticed that at the first month of operation (Figure 5a) the PEMFC operated to follow the highest possible value of the total energy conversion efficiency. Indeed, storing energy in the LIB would imply a lower overall efficiency, as LIB is subject to a charging/discharging efficiency. In the first month, LIB operates in a range of SOC (about 35% to 65%) smaller than in the health-conscious optimization (about 20% to 65%), and hence the LIB degradation rate is only partially affected by the non-health-conscious EMS at the beginning of the plant lifetime. In the last month of operation (Figure 5b), the situation was different. The non-health-conscious EMS would still advantage the PEMFC operation instead of LIB in order not to decrease the efficiency. However, at the end of the plant lifetime there is an increased use of LIB to overcome the voltage loss experienced by PEMFC stacks.
When focusing attention on H2 consumption, the results show that H2 daily consumptions obtained with the FM approach (155 kgH2/day at the first month and 208 kgH2/day at the last month) are similar to the ones obtained with the health-conscious approach (157 kgH2/day at the first month and 204 kgH2/day at the last month). This indeed demonstrates that energy conversion efficiency has not been affected considerably by limiting the degradation of PEMFC.

4.3. Behaviour of the Energy System Variables during the Plant Lifetime

As seen in Section 4.2, the proposed methodology allows to optimally define the daily operation of the ferry during the entire lifetime, and hence it is possible to analyze the behavior of the energy system variables in the long-term operation. Figure 6 shows the daily H2 consumption (Figure 6a), the PEMFC average efficiency (Figure 6b), the volume and weight of the plant (Figure 6c), and the daily cost (Figure 6d) over time. The PEMFC average efficiency is defined as the average conversion efficiency of PEMFC over an operation day. It can be noticed that the hydrogen consumption at the end of the plant lifetime (204 kg/day) is about 29% higher than the hydrogen consumption at the beginning of the plant lifetime (157 kg/day). While this increase could be partially caused by the change in the optimal operation over time (see Section 4.2), the rise of hydrogen daily consumption is mainly due to the decrease of PEMFC average efficiency (Figure 6b). Indeed, with the progressive degradation of PEMFC, the efficiency is decreased and hence a higher amount of hydrogen is needed to reach the same power output.
An important drawback of the increase of hydrogen daily consumption over time is the increase of volume and weight required by the storage system onboard. Figure 6c reports the increase of volume and the weight of the main components of the system (i.e., LFP battery, PEMFC stacks, hydrogen tank). It is possible to observe that the system volume and weight at the end of the plant lifetime increased by up to 8.9% and 1.7%, respectively, with respect to the volume and weight at the beginning of the plant lifetime. While this does not exceed volume and weight constraints of the plant currently installed on board (Equations (5) and (6)), such aspects should be carefully taken into account when designing hybrid PEMFC propulsion systems, as it could affect the space available for the total payload of the vessel.
Lastly, from Figure 6d it is possible to see that the daily cost of the plant, i.e., the value of the f M O (in Equation (37)), increases from 3183 €/day to 5518 €/day. While this increase is partially due to the increase in hydrogen daily consumption (Figure 6a), it can be mainly related to the ongoing PEMFC and LFP battery degradation over time. The average daily cost over the entire plant lifetime was calculated as being equal to 4273 €/day.

5. Conclusions

The present study proposes a new methodology for the definition of the optimal EMS over the long-term operation of hybrid PEMFC/LFP battery propulsion by developing and applying a holistic MILP optimization methodology to a small-size ferry for coastal navigation.
The proposed EMS allows to increase the overall lifetime of the plant from 210 to 600 typical operation days with respect to a FM optimization. Given the current higher cost of PEMFC in comparison with LFP battery, the first component to reach the end of life is the LFP battery system, as its degradation rate weighs less in the overall objective function of the optimization problem. The results also show that the usage of LIB helps in reducing the degradation rate of PEMFC, and plant operation changes with the passing of time to allow the optimal economic performance of the plant. With the progressive ageing of the components, the decrease in the average PEMFC efficiency causes an increase in hydrogen daily consumption (from 157 kg/day to 204 kg/day at the end of the plant lifetime). This reflects on the need to oversize the hydrogen storage system in order to ensure the daily autonomy of the ferry over the entire lifespan. Such an implication is thus jeopardizing the technical feasibility of directly installing on board the system, as well as on the cost and possible total loss of payload. Moreover, the components ageing also affects the plant daily cost, which was found to increases by up to 73% at the end of lifetime with respect to the daily cost at the beginning of the plant lifetime.
The proposed methodology is a useful tool, which can be easily adapted to other types and designs of hybrid PEMFC/LFP battery propulsion systems, as well as to other types of vessels and to other market scenarios. Moreover, the present study points out that changes in the costs of PEMFC, LFP battery and liquid hydrogen could noticeably affect the optimal long-term operation of the energy system. These results can provide some interesting inputs for future research on this topic, such as the analysis of how the optimization results are affected by uncertainties on hydrogen, PEMC and LIB costs. In addition, the uncertainty on the mission profile, caused by variations on hull resistance (due to, e.g., fouling of ship’s hull) or seasonal variability of power will also be considered in future works.
Future studies will be also addressed to the improvement of hybrid PEMFC/LIB powertrains with the addition of onboard solar energy plants and to the possible connection to a battery charging station when ferries are staying at the quay for a long time. Lastly, it should be stressed that safety and regulatory aspects still play a key role in limiting the onboard installation of such systems and hence further analyses should address these points to provide information about the technical feasibility of the plant.

Author Contributions

Conceptualization, C.D., D.P. and R.T.; methodology, C.D., D.P. and R.T.; software, C.D. and D.P.; validation, C.D., D.P. and R.T.; formal analysis, C.D., D.P. and R.T.; investigation, C.D., D.P. and R.T.; data curation, C.D. and D.P.; writing—original draft preparation, C.D., D.P. and R.T.; writing—review and editing, C.D., D.P. and R.T.; visualization, C.D., D.P. and R.T.; supervision, R.T.; project administration, R.T. 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.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

BESSBattery Energy Storage System
DODDepth of Discharge
EMSEnergy Management Strategy
FMFuel consumption Minimization
H2Hydrogen
LFPLithium Iron Phosphate
LIBLithium-Ion Batteries
MILPMixed-Integer Linear Programming
MIQPMixed-Integer Quadratic Programming
NMCNickel Manganese Cobalt
PEMFCPolymer Electrolyte Membrane Fuel Cells
RO-PaxRoll-On/roll-off Passengers
SEISolid Electrolyte Interface
SOCState of Charge
SOFCSolid Oxide Fuel Cells
SOHState of Health
Z Generic objective function of a MILP optimization problem
x * Generic continuous variable of a MILP optimization problem
δ * Generic binary variable of a MILP optimization problem
g Equality constraints of a generic MILP optimization problem
h Inequality constraints of a generic MILP optimization problem
tTime
P demand Propulsion and auxiliary power demand of the ferry
jIndex for PEMFC stacks
P FC Power output of a PEMFC stack
P batt + Discharging power of LIB
P batt Charging power of LIB
nNumber of PEMFC stacks
t i n Initial time step for one day of operation
t f i n Last time step of one day of operation
F F C Fuel consumption of a PEMFC stack
P F C m a x PEMFC stack rated power
V H 2 Volume occupied by the hydrogen storage
V FC Volume occupied by a PEMFC stack
V b a t t Volume occupied by LIB
w H 2 Weight of bunkered hydrogen
w F C Weight of a PEMFC stack
w b a t t Weight of LIB
c o v Oversizing correction factor
I F C Current density of PEMFC stack
k 1 F , k 2 F Linearization coefficients
k 1 P ,   k 2 P Linearization coefficients
δ F C Binary variable defining on/off status of PEMFC stacks
δ s t .     u p Binary variable defining the occurrence of PEMFC start-up phase
F s t a r t Fuel consumption in a start-up phase
I F C m i n Lower limit of PEMFC current density
I F C m a x Upper limit of PEMFC current density
Δ P F C Allowed load variation of a PEMFC stack
d V l o a d PEMFC voltage loss due to load variation
Δ v l o a d Proportionality constant
Δ v s t , u p Proportionality constant
d V s t , u p PEMFC voltage loss due to start-up
d V P F C PEMFC voltage loss due to power load
k 1 d v , k 2 d v Linearization coefficients
d V Total loss of voltage of PEMFC (not cumulative)
d V T O T Cumulative loss of voltage of PEMFC
k 1 d e g , k 2 d e g Linearization coefficients
V r e f , F C Reference maximum voltage for new PEMFC
η F C PEMFC stack efficiency
E b a t t Energy stored in LIB
η b a t t Charging/discharging efficiency of LIB
E b a t t e r y m a x Energy capacity of LIB
Δ t Time interval between two consecutive time steps
C r a t e ± C-rate of LIB in charging/discharging phase
C r a t e m a x Upper limit on LIB C-rate
S O C cal SOC constraint to limit the calendar ageing
T b a t t Temperature of LIB
a 1 T , x ,   a 2 T , x Linearization coefficients
Q l o s s Capacity fade of LIB
T a m b Ambient temperature
a 1 Q T , a 2 Q T Linearization coefficients
A h Energy throughput of LIB
k 11 , A h ,   k 21 , A h , k 31 , A h Proportionality constants
k 12 , A h ,   k 22 , A h , k 32 , A h Proportionality constants
n d a y s Number of days per month
N c y c l e Number of equivalent cycles of LIB
A h c e l l Capacity of a LIB single cell
S O C m i n   Lower limit on SOC
S O C m a x Upper limit on SOC
w 1 , w 2 , w 3 Weights of the objective functions
f 1 , f 2 , f 3 Objective functions
f MO Blended objective function
c H 2 Cost of hydrogen
c b a t t Cost of LIB
c F C Cost of PEMFC
V o PEMFC loss of voltage allowed before substituting PEMFC
Q l o s s m a x Maximum capacity fade before substituting LIB

Appendix A

In the following, details on the ferry chosen as a case study are reported for convenience, while extensive information is available in [39]. Figure A1 shows the power demand profile of the ferry for a typical day of operation. The considered ferry has been assumed to operate 72 journeys per day, each with a duration of five minutes. The propulsion power demand while sailing has been evaluated as a function of typical sailing speeds obtained by fitting the speed data available in [68]. It has been assumed that the sailing phase can be divided into two phases: a first phase at high power and a second one in which the power demand is reduced thanks to the draft. It has been assumed that only auxiliary power is required while the ship is in port. Auxiliary power has been set to 25 kW, constant during the day. In order to account for possible variations in the sea conditions during the day, it has been assumed that the power demand is higher in the central hours of the day due to increased sea-traffic in the area.
Figure A1. Power demand profile of the ferry chosen as a case study [39].
Figure A1. Power demand profile of the ferry chosen as a case study [39].
Energies 14 03813 g0a1

References

  1. Van Biert, L.; Godjevac, M.; Visser, K.; Aravind, P.V. A review of fuel cell systems for maritime applications. J. Power Source 2016, 327, 345–364. [Google Scholar] [CrossRef] [Green Version]
  2. Tronstad, T.; Åstrand, H.H.; Haugom, G.P.; Langfeldt, L. Study on the Use of Fuel Cells in Shipping; European Maritime Safety Agency: Lisboa, Portugal, 2017.
  3. Dall’Armi, C.; Micheli, D.; Taccani, R. Comparison of different plant layouts and fuel storage solutions for fuel cells utilization on a small ferry. Int. J. Hydrogen Energy 2021, 46, 13878–13897. [Google Scholar] [CrossRef]
  4. Anwar, S.; Zia, M.Y.I.; Rashid, M.; de Rubens, G.Z.; Enevoldsen, P. Towards Ferry Electrification in the Maritime Sector. Energies 2020, 13, 6506. [Google Scholar] [CrossRef]
  5. European Maritime Safety Agency (EMSA). Study on Electrical Energy Storage for Ships; EMSA: Lisboa, Portugal, 2020.
  6. Elektra Project: Ballard Receives PO from BEHALA for 3 × 100 kW Fuel Cell Modules to Power German Push Boat. Available online: https://www.ballard.com/about-ballard/newsroom/news-releases/2019/10/02/ballard-receives-po-from-behala-for-3-x-100kw-fuel-cell-modules-to-power-german-push-boat (accessed on 12 January 2021).
  7. Elektra Project: Brennstoffzellen an Bord von Schiffen. Available online: https://www.e4ships.de/english/inland-shipping/elektra/ (accessed on 12 January 2021).
  8. Energy Observer. Available online: https://www.energy-observer.org/ (accessed on 12 January 2021).
  9. Energy Observer: Energy Balance Nassau—Fort-de-France. Available online: https://www.energy-observer.org/innovations/navigation-nassau-fort-de-france (accessed on 12 January 2021).
  10. Sjø, B. Rapport fase 2 Utviklingskontrakt utslippsfri hurtigbåt. In Doffin 2017-138144 Brødrene Aa Westcon Power & Automation; 2017; Available online: https://www.trondelagfylke.no/contentassets/bd8d4260feb14f6bb7503ddc6360e168/brodrene-aa-sluttrapport-offentlig.pdf (accessed on 19 April 2020).
  11. Brødrene Aa Presents the “AERO” a Newly Developed Hydrogen Fast Ferry Concept—Fuel Cells Works. Available online: https://fuelcellsworks.com/news/brodrene-aa-presents-the-aero-a-newly-developed-hydrogen-fast-ferry-concept/ (accessed on 12 January 2021).
  12. Race for Water—A Foundation to Preserve Water. Available online: https://www.raceforwater.org/en/ (accessed on 12 January 2021).
  13. Future Proof Shipping—Creating a Zero-Emissions Shipping World. Available online: https://www.futureproofshipping.com/ (accessed on 12 January 2021).
  14. Toyota Tsusho Conducts Joint Research with Yanmar on Fuel Cell System for Boats—Toward the Creation of a Hydrogen Society by Expanding the Use of Hydrogen. Available online: https://www.toyota-tsusho.com/english/press/detail/180326_004178.html (accessed on 12 January 2021).
  15. Ballard Fuel Cell Modules Used in Successful Yanmar Boat Test to Support Development of Safety Guidelines in Japan. Available online: https://www.prnewswire.com/news-releases/ballard-fuel-cell-modules-used-in-successful-yanmar-boat-test-to-support-development-of-safety-guidelines-in-japan-678507543.html (accessed on 12 January 2021).
  16. Water-Go-Round. Available online: https://watergoround.com/ (accessed on 12 January 2021).
  17. SW/TCH Maritime. Available online: https://www.switchmaritime.com/ (accessed on 16 March 2021).
  18. Enova Banebrytende Samarbeid Gir Millioner til Utslippsfri Teknologi. Available online: https://www.enova.no/pilot-e/banebrytende-samarbeid-gir-millioner-til-utslippsfri-maritim-teknologi/ (accessed on 12 January 2021).
  19. Enova PILOT-E. Available online: https://www.enova.no/pilot-e/ (accessed on 12 January 2021).
  20. Madsen, R.T.; Klebanoff, L.E.; Caughlan, S.A.M.; Pratt, J.W.; Leach, T.S.; Appelgate, T.B.; Kelety, S.Z.; Wintervoll, H.C.; Haugom, G.P.; Teo, A.T.Y.; et al. Feasibility of the Zero-V: A zero-emissions hydrogen fuel-cell coastal research vessel. Int. J. Hydrogen Energy 2020, 45, 25328–25343. [Google Scholar] [CrossRef]
  21. Sandia Energy Maritime Applications for Hydrogen Fuel Cells. Available online: https://energy.sandia.gov/programs/sustainable-transportation/hydrogen/fuel-cells/maritime-applications/ (accessed on 12 January 2021).
  22. Choi, C.H.; Yu, S.; Han, I.S.; Kho, B.K.; Kang, D.G.; Lee, H.Y.; Seo, M.S.; Kong, J.W.; Kim, G.; Ahn, J.W.; et al. Development and demonstration of PEM fuel-cell-battery hybrid system for propulsion of tourist boat. Int. J. Hydrogen Energy 2016, 41, 3591–3599. [Google Scholar] [CrossRef]
  23. Zemship Project Website. Available online: https://ec.europa.eu/environment/life/project/Projects/index.cfm?fuseaction=search.dspPage&n_proj_id=3081 (accessed on 24 April 2020).
  24. Vogler, F.; Würsig, G. Fuel cells in maritime applications: Challenges, chances and experiences. In Proceedings of the 4th International 527 Conference on Hydrogen Safety (ICHS), San Francisco, CA, USA, 12–14 September 2011. [Google Scholar]
  25. Zemship Project. One Hundred Passengers and Zero Emissions—The First ever Passenger Vessel to Sail Propelled by Fuel Cells, 2013. Available online: https://webgate.ec.europa.eu/life/publicWebsite/index.cfm?fuseaction=search.dspPage&n_proj_id=3081 (accessed on 19 April 2020).
  26. Fuel Cell Boat, Amsterdam, The Netherlands. Available online: http://www.opr-advies.nl/fuelcellboat/efcbboot.html (accessed on 12 January 2021).
  27. He, Y.; Miao, C.; Wu, J.; Zheng, X.; Liu, X.; Liu, X.; Han, F. Research on the Power Distribution Method for Hybrid Power System in the Fuel Cell Vehicle. Energies 2021, 14, 734. [Google Scholar] [CrossRef]
  28. Sorlei, I.S.; Bizon, N.; Thounthong, P.; Varlam, M.; Carcadea, E.; Culcer, M.; Iliescu, M.; Raceanu, M. Fuel Cell Electric Vehicles—A Brief Review of Current Topologies and Energy Management Strategies. Energies 2021, 14, 252. [Google Scholar] [CrossRef]
  29. Taccani, R.; Malabotti, S.; Dall’Armi, C.; Micheli, D. High energy density storage of gaseous marine fuels: An innovative concept and its application to a hydrogen powered ferry. Int. Shipbuild. Prog. 2020, 67, 29–52. [Google Scholar] [CrossRef]
  30. Han, J.; Charpentier, J.-F.; Tang, T. An Energy Management System of a Fuel Cell/Battery Hybrid Boat. Energies 2014, 7, 2799–2820. [Google Scholar] [CrossRef] [Green Version]
  31. Bassam, A.M.; Phillips, A.B.; Turnock, S.R.; Wilson, P.A. An improved energy management strategy for a hybrid fuel cell/battery passenger vessel. Int. J. Hydrogen Energy 2016, 41, 22453–22464. [Google Scholar] [CrossRef] [Green Version]
  32. Rivarolo, M.; Rattazzi, D.; Lamberti, T.; Magistri, L. Clean energy production by PEM fuel cells on tourist ships: A time-dependent analysis. Int. J. Hydrogen Energy 2020, 45, 25747–25757. [Google Scholar] [CrossRef]
  33. Rivarolo, M.; Rattazzi, D.; Magistri, L. Best operative strategy for energy management of a cruise ship employing different distributed generation technologies. Int. J. Hydrogen Energy 2018, 43, 23500–23510. [Google Scholar] [CrossRef]
  34. Bassam, A.M.; Phillips, A.B.; Turnock, S.R.; Wilson, P.A. Development of a multi-scheme energy management strategy for a hybrid fuel cell driven passenger ship. Int. J. Hydrogen Energy 2017, 42, 623–635. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, Z.; Guan, C.; Liu, Z. Real-Time Optimization Energy Management Strategy for Fuel Cell Hybrid Ships Considering Power Sources Degradation. IEEE Access 2020, 8, 87046–87059. [Google Scholar] [CrossRef]
  36. Chen, H.; Zhang, Z.; Guan, C.; Gao, H. Optimization of sizing and frequency control in battery/supercapacitor hybrid energy storage system for fuel cell ship. Energy 2020, 197. [Google Scholar] [CrossRef]
  37. Wu, P.; Bucknall, R. Hybrid fuel cell and battery propulsion system modelling and multi-objective optimisation for a coastal ferry. Int. J. Hydrogen Energy 2020, 45, 3193–3208. [Google Scholar] [CrossRef]
  38. Wu, P.; Partridge, J.; Bucknall, R. Cost-effective reinforcement learning energy management for plug-in hybrid fuel cell and battery ships. Appl. Energy 2020, 275, 115258. [Google Scholar] [CrossRef]
  39. Pivetta, D.; Dall’Armi, C.; Taccani, R. Multi-objective optimization of hybrid PEMFC/Li-ion battery propulsion systems for small and medium size ferries. Int. J. Hydrogen Energy 2021. In Press. [Google Scholar] [CrossRef]
  40. Kandidayeni, M.; Macias, A.; Boulon, L.; Trovão, J.P.F. Online modeling of a fuel cell system for an energy management strategy design. Energies 2020, 13, 3713. [Google Scholar] [CrossRef]
  41. Ghobadpour, A.; Amamou, A.; Kelouwani, S.; Zioui, N.; Zeghmi, L. Impact of powertrain components size and degradation level on the energy management of a hybrid industrial self-guided vehicle. Energies 2020, 13, 5041. [Google Scholar] [CrossRef]
  42. Fletcher, T.; Ebrahimi, K. The Effect of Fuel Cell and Battery Size on Efficiency and Cell Lifetime for an L7e Fuel Cell Hybrid Vehicle. Energies 2020, 13, 5889. [Google Scholar] [CrossRef]
  43. Çınar, H.; Kandemir, I. Active energy management based on meta-heuristic algorithms of fuel cell/battery/supercapacitor energy storage system for aircraft. Aerospace 2021, 8, 85. [Google Scholar] [CrossRef]
  44. Liang, J.; Li, Y.; Jia, W.; Lin, W.; Ma, T. Comparison of Two Energy Management Strategies Considering Power System Durability for PEMFC-LIB Hybrid Logistics Vehicle. Energies 2021, 14, 3262. [Google Scholar] [CrossRef]
  45. Wang, Y.; Moura, S.J.; Advani, S.G.; Prasad, A.K. Power management system for a fuel cell/battery hybrid vehicle incorporating fuel cell and battery degradation. Int. J. Hydrogen Energy 2019, 44, 8479–8492. [Google Scholar] [CrossRef]
  46. Li, H.; Ravey, A.; N’Diaye, A.; Djerdir, A. Online adaptive equivalent consumption minimization strategy for fuel cell hybrid electric vehicle considering power sources degradation. Energy Convers. Manag. 2019, 192, 133–149. [Google Scholar] [CrossRef]
  47. Yue, M.; Jemei, S.; Gouriveau, R. Review on health-conscious energy management strategies for fuel cell hybrid electric vehicles: Degradation models and strategies. Int. J. Hydrogen Energy 2019, 44, 6844–6861. [Google Scholar] [CrossRef]
  48. Ito, K.; Yokoyama, R.; Shiba, T. Optimal operation of a diesel engine cogeneration plant including a heat storage tank. J. Eng. Gas Turbines Power 1992, 114, 687–694. [Google Scholar] [CrossRef]
  49. Rech, S. Smart Energy Systems: Guidelines for Modelling and Optimizing a Fleet of Units of Different Configurations. Energies 2019, 12, 1320. [Google Scholar] [CrossRef] [Green Version]
  50. Pivetta, D.; Rech, S.; Lazzaretto, A. Choice of the Optimal Design and Operation of Multi-Energy Conversion Systems in a Prosecco Wine Cellar. Energies 2020, 13, 6252. [Google Scholar] [CrossRef]
  51. Rech, S.; Lazzaretto, A. Smart rules and thermal, electric and hydro storages for the optimum operation of a renewable energy system. Energy 2018, 147, 742–756. [Google Scholar] [CrossRef]
  52. Python. Available online: https://www.python.org/ (accessed on 24 March 2021).
  53. Gurobi Optimization. Available online: https://www.gurobi.com/ (accessed on 16 March 2021).
  54. Jouin, M.; Bressel, M.; Morando, S.; Gouriveau, R.; Hissel, D.; Péra, M.; Zerhouni, N.; Jemei, S.; Hilairet, M.; Ould, B. Estimating the end-of-life of PEM fuel cells: Guidelines and metrics. Appl. Energy 2016, 177, 87–97. [Google Scholar] [CrossRef] [Green Version]
  55. Lorenzo, C.; Bouquain, D.; Hibon, S.; Hissel, D. Synthesis of degradation mechanisms and of their impacts on degradation rates on proton-exchange membrane fuel cells and lithium-ion nickel—Manganese—Cobalt batteries in hybrid transport applications. Reliab. Eng. Syst. Saf. 2021, 212, 107369. [Google Scholar] [CrossRef]
  56. Balestra, L.; Schjølberg, I. Modelling and simulation of a zero-emission hybrid power plant for a domestic ferry. Int. J. Hydrogen Energy 2021, 46, 10924–10938. [Google Scholar] [CrossRef]
  57. Fletcher, T.; Thring, R.; Watkinson, M. An Energy Management Strategy to concurrently optimise fuel consumption & PEM fuel cell lifetime in a hybrid vehicle. Int. J. Hydrogen Energy 2016, 41, 21503–21515. [Google Scholar] [CrossRef] [Green Version]
  58. Chen, H.; Pei, P.; Song, M. Lifetime prediction and the economic lifetime of proton exchange membrane fuel cells. Appl. Energy 2015, 142, 154–163. [Google Scholar] [CrossRef]
  59. Al-gabalawy, M.; Mahmoud, K.; Darwish, M.M.F.; Dawson, J.A.; Lehtonen, M.; Hosny, N.S. Reliable and robust observer for simultaneously estimating state-of-charge and state-of-health of lifepo4 batteries. Appl. Sci. 2021, 11, 3609. [Google Scholar] [CrossRef]
  60. Redondo-Iglesias, E.; Venet, P.; Pelissier, S. Calendar and cycling ageing combination of batteries in electric vehicles. Microelectron. Reliab. 2018, 88–90, 1212–1215. [Google Scholar] [CrossRef] [Green Version]
  61. Wang, J.; Liu, P.; HIcks-Garner, J.; Sherman, E.; Soukiazan, S.; Verbrugge, M.; Tataria, H.; Musser, J.; Finamore, P. Cycle-life model for graphite-LiFePO4 cells. J. Power Source 2011, 196, 3942–3948. [Google Scholar] [CrossRef]
  62. Dubarry, M.; Qin, N.; Brooker, P. Calendar aging of commercial Li-ion cells of different chemistries—A review. Curr. Opin. Electrochem. 2018, 9, 106–113. [Google Scholar] [CrossRef]
  63. Groot, J. State-of-Health Estimation of Li-Ion Batteries: Ageing Models; Chalmers University of Technology: Gothenburg, Sweden, 2014. [Google Scholar]
  64. Internation Energy Agency (IEA). The Future of Hydrogen; IEA: Paris, France, 2019. [Google Scholar]
  65. International Renewable Energy Agency (IRENA). Green Hydrogen Cost Reduction; IRENA: Abu Dhabi, United Arab Emirates, 2020. [Google Scholar]
  66. Ioannis, T.; Dalius, T.; Natalia, L. Li-Ion Batteries for Mobility and Stationary Storage Applications—Scenarios for Costs and Market Growth; Publications Office of the European Union: Luxembourg, 2018; ISBN 9789279972546. [Google Scholar]
  67. Wang, G.; Huang, F.; Yu, Y.; Wen, S.; Tu, Z. Degradation behavior of a proton exchange membrane fuel cell stack under dynamic cycles between idling and rated condition. Int. J. Hydrogen Energy 2018, 43, 4471–4481. [Google Scholar] [CrossRef]
  68. Marine Traffic. Available online: https://www.marinetraffic.com/ (accessed on 19 April 2020).
Figure 1. A typical RO-Pax small size ferry.
Figure 1. A typical RO-Pax small size ferry.
Energies 14 03813 g001
Figure 2. A simplified schematic of the proposed hybrid PEMFC/LFP battery propulsion system for the small-size ferry taken as a case study. Cryogenic tank provides hydrogen to two PEMFC stacks that convert chemical energy into electricity. The latter is sent to an energy management system that directs the power flows to fulfill the energy demand of the ferry (propulsion and auxiliaries demands) and/or to charge the battery.
Figure 2. A simplified schematic of the proposed hybrid PEMFC/LFP battery propulsion system for the small-size ferry taken as a case study. Cryogenic tank provides hydrogen to two PEMFC stacks that convert chemical energy into electricity. The latter is sent to an energy management system that directs the power flows to fulfill the energy demand of the ferry (propulsion and auxiliaries demands) and/or to charge the battery.
Energies 14 03813 g002
Figure 3. SOH of LFP battery and PEMFC during the plant lifetime for (a) fuel consumption minimization only optimization and (b) newly proposed long-term operation optimization approach. The dashed orange lines indicate the threshold values of SOH. The solid black lines represent the SOH of LFP battery, the dashed black lines the SOH of PEMFC.
Figure 3. SOH of LFP battery and PEMFC during the plant lifetime for (a) fuel consumption minimization only optimization and (b) newly proposed long-term operation optimization approach. The dashed orange lines indicate the threshold values of SOH. The solid black lines represent the SOH of LFP battery, the dashed black lines the SOH of PEMFC.
Energies 14 03813 g003
Figure 4. Health-conscious optimal daily operation of the ferry during the (a) first month and (b) last month of operation. Light blue lines represent the ferry power demand, orange lines the PEMFC output power, the blue lines the SOC of LFP battery.
Figure 4. Health-conscious optimal daily operation of the ferry during the (a) first month and (b) last month of operation. Light blue lines represent the ferry power demand, orange lines the PEMFC output power, the blue lines the SOC of LFP battery.
Energies 14 03813 g004
Figure 5. Non-health conscious daily optimal operation of the ferry during the (a) first month and (b) last month of operation. Light blue lines represent the ferry power demand, orange lines represent the PEMFC output power, the blue lines the SOC of LFP battery.
Figure 5. Non-health conscious daily optimal operation of the ferry during the (a) first month and (b) last month of operation. Light blue lines represent the ferry power demand, orange lines represent the PEMFC output power, the blue lines the SOC of LFP battery.
Energies 14 03813 g005
Figure 6. Energy system variables during the whole lifetime: (a) hydrogen daily consumption; (b) PEMFC average efficiency during one day; (c) volume and weight of the main components of the energy system; (d) total cost during one day of operation (value of f M O in Equation (37)).
Figure 6. Energy system variables during the whole lifetime: (a) hydrogen daily consumption; (b) PEMFC average efficiency during one day; (c) volume and weight of the main components of the energy system; (d) total cost during one day of operation (value of f M O in Equation (37)).
Energies 14 03813 g006
Table 1. List of projects on hybrid PEMFC/BESS propulsion systems. Elaborated from [3].
Table 1. List of projects on hybrid PEMFC/BESS propulsion systems. Elaborated from [3].
Project NameStart DateEnd DatePEMFC PowerBattery Type and CapacityRef.
Elektra and Elektra 2 (e4ships consortium)20172024300 kWLIB, 2.5 kWh capacity[6,7]
Energy observer2017202422 kWLIB, 122 kWh propulsion, 18 kWh daily facilities[8,9]
Aero 42-20232800 kWLIB, 672 kWh[10,11]
Race for water2015202160 kWLIB, 745 kWh[12]
Future Proof Shipping-2021635 kWLIB, 300 kWh[13]
Kamine boat2018202060 kWLIB, 60 kWh[14,15]
Water Go Round Golden Gate Zero Emission-2020360 kW100 kWh, no data on the type[16,17]
Zeff (part of Pilot-E scheme)201820202200 kW50 kWh estimated; type of batteries under evaluation[18,19]
Zero V-20171800 kWLIB provide for 200 kW power[20,21]
Busan tourist boat-201656 kWLIB, 47 kWh[22]
ZEMSHIP FSC Alsterwasser20062013100 kWLead gel, 560 V (7 × 80 V) 360 Ah[23,24,25]
Nemo H22008201160–70 kWLead acid, 55 kW[26]
Table 2. The main characteristics of the RO-Pax ferry proposed as a case study. Power system volume and weight have been elaborated from data available on diesel storage on board and a diesel-electric generator for marine propulsion, scaled by the installed power of the ferry.
Table 2. The main characteristics of the RO-Pax ferry proposed as a case study. Power system volume and weight have been elaborated from data available on diesel storage on board and a diesel-electric generator for marine propulsion, scaled by the installed power of the ferry.
UnitValue
Vessel Type-RO-Pax ferry
Length Overallm42
Breadth Extremem9
Gross Tonnaget280
Propulsion engineskW2 × 206
Auxiliary engineskW2 × 28
Power system volumem315.5
Power system weightkg7185
Table 3. Input parameters of the multi-objective optimization.
Table 3. Input parameters of the multi-objective optimization.
ParametersUnitValueRef.
c H 2 €/kWh0.3[64,65]
c b a t t €/kWh818[5,66]
c F C €/kW3750[64]
c o v %20Assumed
Δ v l o a d µV/kW0.0441[37,55,56,57,58]
Δ P F C %10Assumed
Δ v s t , u p µV/cycle23.91[37,55,56,57,58]
Δ t min1Assumed
η b a t t %90Assumed
F s t a r t %10Assumed
I F C m a x A/cm2600[67]
I F C m i n A/cm275[67]
k 1 d e g µVcm2/A0.243[37,55,56,57,58]
k 2 d e g µV0.159[37,55,56,57,58]
k 1 F kWcm2/A0.354[57,58,67]
k 2 F kW3.837[57,58,67]
k 1 P kWcm2/A0.159[57,58,67]
k 2 P kWcm2/A7.251[57,58,67]
k 11 , A h 1/(KAh2) 7.1 · 10 9 [61,62,63]
k 21 , A h 1/(KAh) 2.2 · 10 4 [61,62,63]
k 31 , A h 1/K0.033[61,62,63]
k 12 , A h 1/(KAh2) 2.1 · 10 6 [61,62,63]
k 22 , A h 1/(KAh)−0.48[61,62,63]
k 32 , A h 1/K−8.885[61,62,63]
S O C m i n %20[63]
S O C m a x %90[63]
V b a t t m3/kWh0.0091[5]
V F C m3/kW0.0312[64]
V H 2 m3/kgH20.025[64,65]
V m a x m315.5Assumed
V o V0.19Assumed
w b a t t kg/kWh8[5]
w F C kg/kW20[64]
w H 2 kg/kgH22.5[64,65]
w m a x kg7185Assumed
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dall’Armi, C.; Pivetta, D.; Taccani, R. Health-Conscious Optimization of Long-Term Operation for Hybrid PEMFC Ship Propulsion Systems. Energies 2021, 14, 3813. https://doi.org/10.3390/en14133813

AMA Style

Dall’Armi C, Pivetta D, Taccani R. Health-Conscious Optimization of Long-Term Operation for Hybrid PEMFC Ship Propulsion Systems. Energies. 2021; 14(13):3813. https://doi.org/10.3390/en14133813

Chicago/Turabian Style

Dall’Armi, Chiara, Davide Pivetta, and Rodolfo Taccani. 2021. "Health-Conscious Optimization of Long-Term Operation for Hybrid PEMFC Ship Propulsion Systems" Energies 14, no. 13: 3813. https://doi.org/10.3390/en14133813

APA Style

Dall’Armi, C., Pivetta, D., & Taccani, R. (2021). Health-Conscious Optimization of Long-Term Operation for Hybrid PEMFC Ship Propulsion Systems. Energies, 14(13), 3813. https://doi.org/10.3390/en14133813

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