Hydrodynamic Simulation of Green Hydrogen Catamaran Operating in Lisbon, Portugal

: Similar to other industries, the maritime industry is also facing increasing restrictions on ships regarding pollution control. The research presented in this paper is aimed at studying the pros and cons of alternative fuels followed by a detailed analysis on hydrogen fuel cells (PEMFC) for a particular ship operating in Lisbon, Portugal. Dynamic forces acting on the ship have been studied for a year. Assessing various scenarios based on these results aids ship operators in making informed decisions regarding the future course of action for their existing vessels. These different cases are ﬁrst: business as usual (diesel engine), second: replacing the diesel engine with a hydrogen hybrid system and, third: replacement of the ship with a new hydrogen hybrid ship. The study is based on the simulation of numerical equations and CFD simulation results. As the result, the second scenario is best suited in both aspects; namely, environmental and economic.


Introduction
It is estimated that around 80-90% of international trade is facilitated by maritime vessels, which is largely attributed to its cost efficiency in transportation [1].There is a direct relationship between world trade and world economic growth [2,3].Therefore, with the increasing world economy and shipping industry, carbon emissions are also expected to swell if business continues as usual.To address this challenge, the International Maritime Organization (IMO) has made several efforts to control these emissions by formulating emissions related to policies and regulations.There are a range of measures that an operator can adopt to achieve these reductions like regulatory requirements (energy efficiency design index, ship energy efficiency management plan, etc.), advancements in technology (improved engines, alternative fuels, etc.), adjustments related to a ship's operation (slow steaming, low emission fuels, etc.), and market-based instruments (cap-and-trade systems, rebate mechanism, etc.) [4].On top of the IMO's regulations, local governments are also actively taking part in policy making.The recent example of such a policy is an agreement related to emission trading systems (EU ETS) made by the EU's legislative bodies in January 2023, and to be enforced in January 2024, which penalizes the ship operators based on the amount of carbon emitted.
It is evident that shipping significantly impacts both the economy and the environment.Therefore, it is important to study available options to reduce the environmental impact while the shipping industry keeps booming.Pros and cons of available alternative fuels are studied to boil down to only one which is than studied on the selected ship.A SIMULINK model is designed for the catamaran ferry "Almadense" (the selected ship) operating on Tagus River.The following objectives are studied:

•
To study an alternative fuel technology (the hydrogen fuel cell) for the ship considered in this study.

•
To develop a MATLAB SIMULINK model of a power system and subsequently reproduce year-long simulation outcomes by utilizing year-long data inputs.

•
To perform lifetime economic and environmental analyses for all the three considered cases.

•
To suggest the optimum size of hybrid system components.

Literature Review
Carbon emissions contribution via shipping increased from 1.8% in 1996 to 2.9% in 2018 [5,6].Numerous studies have been conducted to evaluate the potential application of different technologies across a spectrum of ship sizes and travel routes.The economic factor proved to be the most influencing factor for the implementation of potential technologies [5].
Another key parameter for the selection of new technology is the required power and distance for travel.This is due to the size and storage restrictions along with efficiencies, route predictability, present technology capabilities, and up-front costs [7].
Although very low sulphur fuel oil (VLSFO) is currently in use, because of new IMO regulation enforcements that went into effect in 2020, there are a variety of alternate fuels such as liquefied natural gas (LNG), methanol, liquid petroleum gas (LPG), ethanol, hydrogen, etc., that have better GHG emission reduction potentials.Amongst these, currently LNG is the most studied alternative fuel.On the other hand, green hydrogen provides a cleaner and efficient fuel, but it is limited to short-range shipping [5]; these fuels have been compared with conventional marine gas oils (MGO) in terms of their physical parameters [8], as shown in Table 1.These parameters are especially important when considering their storage requirements.
Table 1.Properties of different fuels that could be used as alternative fuels for power generation [8].In terms of transportation, liquid or gaseous hydrogen transportation is more efficient compared to ammonia [9].

Parameter
A hybrid combination of hydrogen and green ammonia is suggested in [10] for the best environmental performance whereas dual fuel engine with LNG is considered to have the best net present value (NPV) performance.These calculations also include carbon taxes, which are not applicable for the case study presented in this paper.Ammonia is also toxic for humans and it could incur a negative NPV due to its high operation costs [10,11].
For short-distance transits of small-to medium-sized ships, compressed hydrogen emerged as the dominant choice [12].In a hybrid system, fuel cells (FCs) have a higher energy density with a considerably lower power density, whereas batteries exhibit a higher power density, which complements the characteristics of FCs (as shown in Figure 1), summing up to make a hybrid system.
Electrochemical energy is harvested and converted into electrical energy in a hydrogen FC.Fundamentally, there are different types of FC technologies available today, amongst which the PEMFC (proton exchange membrane fuel cell) is the best suited for ships with lower power demands like ferries, boats, and fishing vessels.The PEMFC is amongst the lightest and most compact FC with high power density and efficiency.Additionally, PEMFCs have a higher cold start potential with better transient performance [13].amongst the lightest and most compact FC with high power density and efficiency.Additionally, PEMFCs have a higher cold start potential with better transient performance [13].
Figure 1.Ragone plot adapted from [14] for various types of emerging green technologies to power the energy system based on their power and energy densities.
The study presented in [15] ranks FC technologies on different parameters; namely, safety, emissions, efficiency, cost, lifetime, power output, fuel type, and size on a scale of five.The Molthen Carbonate Fuel Cell (MCFC) operating on diesel oil (diesel oil undergo reforming process to extract hydrogen) scores the highest, although it is not an environment friendly alternative.The PEMFC scores the second highest.
A robust and fairly accurate numerical simulation of the PEMFC was performed by [16] on a marine vessel.The study concentrated on assessing the dynamic response of fuel cells within a hybrid system subjected to step load changes.This highlights the possibility of achieving a good fidelity using dynamic simulations.Another study [17] quantified the emission reduction of a ferry in Denmark that used a hybrid hydrogen system; the authors determined a total decline of 65% in greenhouse gases.This massive reduction highlights the potential of such technology in achieving the desired goal of emissions reduction.But this reduction in emission comes with an increased economic cost.The adoption of the PEMFC as the primary power source for marine ships results in a heightened cost, ranging between 35% and 38% higher than conventional diesel propulsion and with this price surge, a notable reduction of 73-78% in GHGs emissions can be achieved [18].
Further insights into this decline have been examined in reference [19], where a liquified natural gas fuel cell (LNGFC) demonstrated a superior performance compared to a PEMFC.However, LNGFCs comes with a cost that is nearly three times the price of PEM-FCs.AFC, SOFC, and low-temperature PEMFCs exhibit efficiencies in the range of 50% to 60%, surpassing all the other FCs.Quantitatively, a low-temperature PEMFC consumes the least hydrogen fuel (approximately 60 g/kWh) [19].Apart from all these advantages of PEMFCs over other FC types, PEMFCs are easily available in the market and have extended maintenance intervals [19], making it a better choice overall.These are potential figures and need to be verified in actual scenarios.There are 3123 ferries which are equipped with FCs onboard [18], some of which only utilize FCs for auxiliary power demands.This suggests that the implementation of FCs in the marine sector is still in the R&D phase [20].Therefore, a significant economic investment along with technological development is required.[14] for various types of emerging green technologies to power the energy system based on their power and energy densities.

Identified Research Gaps
The study presented in [15] ranks FC technologies on different parameters; namely, safety, emissions, efficiency, cost, lifetime, power output, fuel type, and size on a scale of five.The Molthen Carbonate Fuel Cell (MCFC) operating on diesel oil (diesel oil undergo reforming process to extract hydrogen) scores the highest, although it is not an environment friendly alternative.The PEMFC scores the second highest.
A robust and fairly accurate numerical simulation of the PEMFC was performed by [16] on a marine vessel.The study concentrated on assessing the dynamic response of fuel cells within a hybrid system subjected to step load changes.This highlights the possibility of achieving a good fidelity using dynamic simulations.Another study [17] quantified the emission reduction of a ferry in Denmark that used a hybrid hydrogen system; the authors determined a total decline of 65% in greenhouse gases.This massive reduction highlights the potential of such technology in achieving the desired goal of emissions reduction.But this reduction in emission comes with an increased economic cost.The adoption of the PEMFC as the primary power source for marine ships results in a heightened cost, ranging between 35% and 38% higher than conventional diesel propulsion and with this price surge, a notable reduction of 73-78% in GHGs emissions can be achieved [18].
Further insights into this decline have been examined in reference [19], where a liquified natural gas fuel cell (LNGFC) demonstrated a superior performance compared to a PEMFC.However, LNGFCs comes with a cost that is nearly three times the price of PEMFCs.AFC, SOFC, and low-temperature PEMFCs exhibit efficiencies in the range of 50% to 60%, surpassing all the other FCs.Quantitatively, a low-temperature PEMFC consumes the least hydrogen fuel (approximately 60 g/kWh) [19].Apart from all these advantages of PEMFCs over other FC types, PEMFCs are easily available in the market and have extended maintenance intervals [19], making it a better choice overall.These are potential figures and need to be verified in actual scenarios.There are 3123 ferries which are equipped with FCs onboard [18], some of which only utilize FCs for auxiliary power demands.This suggests that the implementation of FCs in the marine sector is still in the R&D phase [20].Therefore, a significant economic investment along with technological development is required.

•
There is a lack of studies that simulate the model under real environmental conditions.Most of the studies were performed while assuming the average condition of a voyage using simplified numerical methods.

•
There are no studies, to the best of the authors' knowledge, that analyzes the lifetime economic cost of running the ship (with new technology) and the environmental emissions for an already commissioned ship under different possible scenarios.Such an analysis would consider the scenario of the ship's continuous operation on diesel for its remaining lifespan, or with a transition to alternate technology, or through decommissioning and replacement with a new alternate fuel-powered ship.

•
According to the literature review, it is clear that a unique analysis of the selected technology for a particular ship on a specified route is required.

•
There is a lack of studies that calculate the power system component sizing based on the real environmental condition simulation that balances both economic and environmental concerns.

Methodology
This study focuses on the analysis of the catamaran ferry "Almadense", currently in operational service in Lisbon whose specifications are as per Table 2. Operating along the River Tagus, this diesel-powered catamaran connects three key ports: Porto Brandão, Belém, and Trafaria, as illustrated in Figure 2. The vessel boasts a catamaran-style hull design renowned for its high efficiency, primarily attributed to minimized hull drag forces.

Case I-Hydrogen Hybrid System Modelling
To have an optimal design of an energy system, a hybrid power system is designed in MATLAB SIMULINK with the principle as shown in Figure 3.The design simply works on the principle that in order to propel the ship at a certain advance speed, it has to overcome all the resistance that it would encounter which includes hull, waves, and air re- The ship starts its operation in the morning at 05:30 a.m. and the last voyage finishes at 09:30 p.m. when it arrives at Porto Brandão, concluding a daily voyage composed of 55 trips.The block coefficient "C B " of the ship has been calculated numerically using Equation (1) where Fn is the Froude number.This C B is essential to calculating the hydrodynamic forces (wave and hull resistance).There are many ways of estimating it; one that has been used in this study includes Fn (Froude number) according to [22].The study limited the range of Fn between 0.27 and 0.36 for passenger ships [22].This range is used as a verification method to verify the C B .The ship considered in this study satisfies this range with Fn of 0.28.This acts as additional yet simpler approach for C B verification:

Case I-Hydrogen Hybrid System Modelling
To have an optimal design of an energy system, a hybrid power system is designed in MATLAB SIMULINK with the principle as shown in Figure 3.The design simply works on the principle that in order to propel the ship at a certain advance speed, it has to overcome all the resistance that it would encounter which includes hull, waves, and air resistances.If the power supplied overcomes these resistances, the ship will start moving at the desired speed.

Case I-Hydrogen Hybrid System Modelling
To have an optimal design of an energy system, a hybrid power system is desig in MATLAB SIMULINK with the principle as shown in Figure 3.The design simply w on the principle that in order to propel the ship at a certain advance speed, it has to o come all the resistance that it would encounter which includes hull, waves, and ai sistances.If the power supplied overcomes these resistances, the ship will start movin the desired speed.The speed profile of the ship for an average day is used as an input to simulate designed model.The velocity demand profile is kept constant (as shown in Figure 4) The speed profile of the ship for an average day is used as an input to simulate the designed model.The velocity demand profile is kept constant (as shown in Figure 4), but environmental conditions of the ship's operation vary.This variation changes the power demand even for the same velocity demand.
The power needed by the ship is influenced by various external forces as previously mentioned.An outlined representation of the power required is depicted in Figure 5, incorporating factors such as the instantaneous wave height (ξ wave ), water-specific weight (σ wave ), wind speed (v air ), actual wind direction (θ air ), air pressure (P air ), air density (ρ air ), resulting ship speed (v ship ), and direction (θ ship ).To provide this power demand, a modular FC with a unit capacity of 200 kW inspired by Ballard's FcWave module [23] is used in this study.The characteristic and efficiency curves of this FC are assumed as shown in Figures 6 and 7, respectively.The power needed by the ship is influenced by various external forces as previously mentioned.An outlined representation of the power required is depicted in Figure 5, incorporating factors such as the instantaneous wave height (ξwave), water-specific weight (σwave), wind speed (vair), actual wind direction (θair), air pressure (Pair), air density (ρair), resulting ship speed (vship), and direction (θship).To provide this power demand, a modular FC with a unit capacity of 200 kW inspired by Ballard's FcWave module [23] is used in this study.The characteristic and efficiency curves of this FC are assumed as shown in Figures 6 and 7, respectively.environmental conditions of the ship's operation vary.This variation changes the power demand even for the same velocity demand.The power needed by the ship is influenced by various external forces as previously mentioned.An outlined representation of the power required is depicted in Figure 5, incorporating factors such as the instantaneous wave height (ξwave), water-specific weight (σwave), wind speed (vair), actual wind direction (θair), air pressure (Pair), air density (ρair), resulting ship speed (vship), and direction (θship).To provide this power demand, a modular FC with a unit capacity of 200 kW inspired by Ballard's FcWave module [23] is used in this study.The characteristic and efficiency curves of this FC are assumed as shown in Figures 6 and 7, respectively.The FC module generates varying voltage and current outputs, necessitating the use of a DC-DC boost converter to maintain a constant voltage for the propulsion motor.Multiple simulations were conducted, altering the number of FC modules to determine the optimal module numbers.A crucial determining factor is the annual hydrogen consumption, directly influencing the operating expense (OPEX).To calculate the hydrogen used, the hydrogen mass flow is required, which is calculated during the simulation for every simulation step using Equation (2): where  * is the power delivered by the FC, η is the FC efficiency for the corresponding load, and  is the lower heating value of hydrogen.The DC-DC boost converter stabilizes the voltage output to 815 V corresponding to the voltage requirements of the DC propulsion motor (because of high torque).The motor used here is an ABB DC motor with Cat.No. 3BSM003050-XVJ with a power rating of 1355 kW at 970 rpm.
Apart from the power requirement of the propulsion, the power system in this study is also designed to supply energy to the ship's auxiliary demands.The auxiliary system includes all the necessary equipment like communication, navigation, fire management systems, etc., along with all the devices required for comfortable movements of crew and passengers like air conditioning system, lighting, power ramps, etc.This study assumes that this power is fixed at a constant value of 50 kW.The power supplied while in port is sourced from a battery power pack.In the designed hybrid power system, if the number The FC module generates varying voltage and current outputs, necessitating the use of a DC-DC boost converter to maintain a constant voltage for the propulsion motor.Multiple simulations were conducted, altering the number of FC modules to determine the optimal module numbers.A crucial determining factor is the annual hydrogen consumption, directly influencing the operating expense (OPEX).To calculate the hydrogen used, the hydrogen mass flow is required, which is calculated during the simulation for every simulation step using Equation (2): .
where P f c * is the power delivered by the FC, η f c is the FC efficiency for the corresponding load, and LHV H 2 is the lower heating value of hydrogen.The DC-DC boost converter stabilizes the voltage output to 815 V corresponding to the voltage requirements of the DC propulsion motor (because of high torque).The motor used here is an ABB DC motor with Cat.No. 3BSM003050-XVJ with a power rating of 1355 kW at 970 rpm.Apart from the power requirement of the propulsion, the power system in this study is also designed to supply energy to the ship's auxiliary demands.The auxiliary system includes all the necessary equipment like communication, navigation, fire management systems, etc., along with all the devices required for comfortable movements of crew and passengers like air conditioning system, lighting, power ramps, etc.This study assumes that this power is fixed at a constant value of 50 kW.The power supplied while in port is sourced from a battery power pack.In the designed hybrid power system, if the number of FC modules is less (while using a simulation to track different module numbers) but there is a higher power demand exceeding FC capacity, the battery pack supplements the FC system for propulsion.This leads to an increase in the battery size to accommodate the combined power needs.Apart from this, the battery also provides interim power until the FC catches up.Equation (3) represents the total power available in the bus bar where I BATTERY and I BFC are currently drawn from the battery and the fuel cell modules, respectively, and V is the bus bar voltage: The PMS evaluates real-time conditions acting on the vessel, and their effects are compared with required power, i.e., the desired speed.After evaluation, as a result, its output is the net required power for the ferry to continue to propel at the set speed, which instructs the energy producers either to increase or decrease the power.This is accomplished in simultaneous stages by performing simultaneous calculations.The PMS calculates the difference between the delivered power and the power opposing the propulsion ship due to environmental forces (namely, wave and wind forces and air resistance).
The PMS senses the difference in power and guides the energy providers either to increase or decrease power delivery.These decisions are defined as per Equations ( 4)-( 7): P req (t) = P req * (t) + P aux 0 ≤ P * req (t) − P aux ; and P * req (t) − P aux ≤ P max motor (7) where P error (t) is the generated difference, P e f (t) is the total power opposing the propulsion, P del (t) is the power delivered by the propulsion motor, P req * (t) is the power requirement except for auxiliary power, and K p , K d , and K i are the controller constants which are obtained after repeated trials of many combinations.Lastly, F wave (t), F air (t), and F hull (t) are the external environment forces.The total forces are calculated in SIMULINK as shown in Figure 8; P req (t) is the total required power estimated by the PMS and v(t) is the ship's surge velocity.The power is limited by maximum motor power (P max motor ).Here, P req * (t) is estimated by using a PID controller.This controller varies the amount of hydrogen and oxygen (or battery power in the case of battery operation) using the standard PID Equation (5).Auxiliary power demand is added to this, making the total demand power P req (t), whereas the opposing forces demands P e f (t).Fundamentally, this P e f (t) has to be satisfied.After all the system losses, P del (t) is the power that has been delivered by the propeller which has to satisfy P e f (t) need.If there is a mismatch, P error (t) appears, resulting in an increase or decrease in P req * (t).This process of repetitive iterations will continue until P error (t) reaches zero (or to a minimum error limit).The wave resistance is calculated via Equation ( 8) [24].The wave height (ξ(t)) in this equation is dynamic and is different for every month.
Here, σ is the specific weight of water, L is the length of the ship, and B corresponds to the breadth of the ship.For the simulation performed in this study, the average of each month is used for every second of the ship's operational day.Similarly, the air resistance is governed via Equation ( 9) where the density of air varies each month of the year.This density depends on ambient temperature and pressure conditions: The wave resistance is calculated via Equation ( 8) [24].The wave height (ξ(t)) in this equation is dynamic and is different for every month.
Here, σ is the specific weight of water, L is the length of the ship, and B corresponds to the breadth of the ship.For the simulation performed in this study, the average of each month is used for every second of the ship's operational day.Similarly, the air resistance is governed via Equation (9) where the density of air varies each month of the year.This density depends on ambient temperature and pressure conditions: where C d is the drag coefficient of the ship that depends on the angle of attack of the incoming wind; ρ air is the density of air; A proj. is the projected sectional area of the ship for incoming wind that depends on the angle of the incoming wind.Lastly, v attack is the resultant velocity of incoming wind.This velocity is relative to the actual wind speed and direction in addition to the ship's speed and direction.This resultant wind velocity and direction is calculated as a vector product of both the velocities as represented in Figure 9.The wave resistance is calculated via Equation ( 8) [24].The wave height (ξ(t)) in this equation is dynamic and is different for every month.
Here, σ is the specific weight of water, L is the length of the ship, and B corresponds to the breadth of the ship.For the simulation performed in this study, the average of each month is used for every second of the ship's operational day.Similarly, the air resistance is governed via Equation (9) where the density of air varies each month of the year.This density depends on ambient temperature and pressure conditions: where  is the drag coefficient of the ship that depends on the angle of attack of the incoming wind;  is the density of air;  . is the projected sectional area of the ship for incoming wind that depends on the angle of the incoming wind.Lastly,  is the resultant velocity of incoming wind.This velocity is relative to the actual wind speed and direction in addition to the ship's speed and direction.This resultant wind velocity and direction is calculated as a vector product of both the velocities as represented in Figure 9. Lastly, the hull force (the force exerted by water) is calculated in the simulation.This force depends on the draft of the ship as the higher the draft, the higher the hull force.Therefore, this force is simulated by varying the draft using a randomizer block of the SIMULINK library depicting change in the ship's draft (passenger + cars).This randomizer block generates a random signal for each voyage that changes the number of passengers and vehicles transported.The number of passengers and vehicles changes the displacement of the ship and hence affects the draft according to Equations ( 10) and (11).A CFD analysis of the designed model satisfying the physical parameters was conducted using the CFD simulation ANSYS Fluent software [25] to calculate the force exerted by moving water at various drafts.The results are input into the SIMULINK model as a 2D lookup table (see Figure 10).During the simulations, this 2D lookup table provides the force corresponding to the calculated draft.
T = T light weight + T added (10) Here, T is the calculated draft, T light weight is the light weight draft of the ship, T added is the draft added due to the added weight in meters, T added is the added draft in centimeters, W added is the total weight added, and TPC is the parameter of the ship that represents the weight needed to change the draft by 1 cm.The study involves the meticulous consideration of monthly data, wherein environmental parameters were meticulously recorded every second throughout the ship's operational duration.The process involved calculating the average environmental parameters for each specific time instance (e.g., 6:31 a.m.) across the entire month.This method was consistently applied across various months within the year, enabling the representation of average environmental parameters for each month.
Here, T is the calculated draft,  is the light weight draft of the ship,  is the draft added due to the added weight in meters,  is the added draft in centimeters,  is the total weight added, and  is the parameter of the ship that represents the weight needed to change the draft by 1 cm.The study involves the meticulous consideration of monthly data, wherein environmental parameters were meticulously recorded every second throughout the ship's operational duration.The process involved calculating the average environmental parameters for each specific time instance (e.g., 6:31 AM) across the entire month.This method was consistently applied across various months within the year, enabling the representation of average environmental parameters for each month.
The sum of all these forces results in the total force, and consequently the product of this total force and the actual velocity of the ship provides the opposing power  ().

Energy Management System (EMS)
The EMS is responsible for the decision to distribute the power requirement between both the energy providers (FCs and the battery system).To accomplish that, firstly the system must check if the vessel is in port or sailing since the FC system is not designed in The sum of all these forces results in the total force, and consequently the product of this total force and the actual velocity of the ship provides the opposing power P e f (t).

Energy Management System (EMS)
The EMS is responsible for the decision to distribute the power requirement between both the energy providers (FCs and the battery system).To accomplish that, firstly the system must check if the vessel is in port or sailing since the FC system is not designed in this study to operate in port and switch to ideal conditions (air purging).The EMS is responsible for deciding the energy distribution as shown in the decision flow diagram Figure 11.The FC capacity depends on the number of modules that change for each simulation.

Case II-Diesel Engine Modelling
It is important to compare the hybrid system with a conventional diesel engine-powered system (current system) in order to produce a holistic analysis.To perform this analysis, the diesel engine system is also modelled in SIMULINK and simulated under exactly the same conditions.In this study, a MAN Diesel & Turbo 8L23/30A 4-stroke IMO Tier II engine [26] is

Case II-Diesel Engine Modelling
It is important to compare the hybrid system with a conventional diesel enginepowered system (current system) in order to produce a holistic analysis.To perform this analysis, the diesel engine system is also modelled in SIMULINK and simulated under exactly the same conditions.In this study, a MAN Diesel & Turbo 8L23/30A 4-stroke IMO Tier II engine [26] is considered with a power rating of 1280 kW.Also, as there is only one energy provider (a diesel engine), the EMS is omitted.The only variables left to calculate are the efficiency and the fuel consumption.For fuel consumption calculations, the specific fuel oil consumption (SFOC) outlined under the reference conditions has been employed, resulting in the tabulated consumptions presented in Table 3: But in actual conditions, the SFOC changes with the change in ambient conditions (temperature and pressure).These conditions are measured, and a factor named the fuel consumption factor (γ) is calculated using Equations ( 12) and ( 13) which is then multiplied by the reference load SFOC: where, t r , t bar , and p r are the reference ambient temperature, charge air temperature before cylinder, and ambient air pressure, respectively, as provided in the manufacturer's manual [26], whereas t x , t bax , and p x are the ambient temperature, charge air temperature before cylinder, and ambient air pressure, respectively, corresponding to the conditions of the simulation.Lastly, SFOC r is the reference SFOC as given in Table 3 and SFOC x is the calculated SFOC engine under the simulated environmental conditions.The schematic diagram of the designed SIMULINK model is shown in Figure 12.The PMS of Case II (Diesel Engine Modelling) remains the same as the one considered in Case I (hydrogen hybrid modelling).The remaining modelling is majorly based on the SFOC of the engine under consideration.
J. Mar.Sci.Eng.2023, 11, x FOR PEER REVIEW 12 of 23 [26], whereas  ,  , and  are the ambient temperature, charge air temperature before cylinder, and ambient air pressure, respectively, corresponding to the conditions of the simulation.Lastly,  is the reference SFOC as given in Table 3 and  is the calculated SFOC engine under the simulated environmental conditions.The schematic diagram of the designed SIMULINK model is shown in Figure 12.The PMS of Case II (Diesel Engine Modelling) remains the same as the one considered in Case I (hydrogen hybrid modelling).The remaining modelling is majorly based on the SFOC of the engine under consideration.To facilitate a comparative analysis between Case I and Case II in terms of economic and environmental impacts, it is imperative to measure the quantity of the diesel fuel utilized.The calculation of the diesel flow rate involves the utilization of the power output and specific fuel oil consumption (SFOC), as determined via Equation ( 14): To facilitate a comparative analysis between Case I and Case II in terms of economic and environmental impacts, it is imperative to measure the quantity of the diesel fuel utilized.The calculation of the diesel flow rate involves the utilization of the power output and specific fuel oil consumption (SFOC), as determined via Equation ( 14): .
The overall efficiency of the diesel engine is calculated as the ratio between the power delivered by the engine and the power potential of the fuel used to deliver that power.This is computed using Equations ( 15) and ( 16):

Case I-Hydrogen Hybrid System Modelling
Energy distribution amongst the FC modules and the battery largely depends on the size of the FC.The battery acts as a buffer until either the FC system takes over or the FC does not have enough potential.
The hydrogen consumption over a year follows the "U" shaped characteristic as can be seen in Figure 13.This trend emerges due to the decreased external forces during the summer months, resulting in diminished air and water densities alongside reduced wave heights.Hydrogen consumption reaches its peak in the month of November.The influence of the external forces on the ship dynamics reveals a noteworthy dominance of wave and hull forces over air drag.While air drag occasionally aids propulsion rather than opposing it, its impact remains notably minimal.The predominance between wave and hull forces fluctuates seasonally, with the wave force exhibiting significant variation due to seasonal fluctuations in wave heights, whereas the hull force remains relatively stable (see Figure 14).drag.While air drag occasionally aids propulsion rather than opposing it, its impact remains notably minimal.The predominance between wave and hull forces fluctuates seasonally, with the wave force exhibiting significant variation due to seasonal fluctuations in wave heights, whereas the hull force remains relatively stable (see Figure 14).Figure 13 represents the hydrogen consumption for different FC modules.The total consumption per year increases with the increase in the number of FC modules until six and then experiences a dip as shown in Table 4.This is because as the number of FC modules increases, the load on each module decreases and FCs have higher electrical efficiency with low loads (see Figure 7).Figure 13 represents the hydrogen consumption for different FC modules.The total consumption per year increases with the increase in the number of FC modules until six and then experiences a dip as shown in Table 4.This is because as the number of FC modules increases, the load on each module decreases and FCs have higher electrical efficiency with low loads (see Figure 7).The simulations are also helpful in deciding the capacity of the battery to be installed.The power required at every instance remains the same regardless of the FC module number.As an example, if a 1000 kW power is required at a given time but the FC modules only have the capacity of delivering 800 kW (in case of 4 FC modules), the remaining power will be compensated for by the battery.Consequently, the annual total energy provided by the battery is observed to be higher in instances with fewer FC modules, as demonstrated in Table 5.
With the increase in FC modules, the energy consumption of the battery decreases.This is simply because the required energy is primarily delivered through the FCs.But Table 5 shows that after seven modules, the consumption does not decrease further.The explanation behind this fact is that after seven modules, the FCs have enough potential to match the required power and the battery is only essential while in port.

Case II-Diesel Engine Modelling
After simulation, Case II also follows the similar "U" shaped fuel consumption characteristic.This effect can also be seen in the engine efficiency for different months.The maximum fuel consumption occurs in the month of November, similarly to the results obtained for Case I.In total, the diesel engine consumes 472,136 kg of diesel in a year.

Module Selection and Economic Viability Study
Apart from the technical analysis, an economic viability study is equally important in order to make the project's selection eligible for investors.Also, according to a study conducted by [5] (as referred to in the Literature Review section), the economic factor is the most powerful influence in the selection of technology.Therefore, the FC module with the maximum return has been selected.The economic viability study is based on various market parameters, such as inflation, discount rate, etc., that affect the return of the project.These parameters are calculated by using the assumptions established in Table 6.Linear depreciation of equipment is considered.To calculate the lifetime returns of a project, the components' lifetimes, capital expenditure (CAPEX), and OPEX are required.Table 7 shows these values for all the components used in this study.The net income of a project depends on revenues generated and on the expenses.For the ship Almadense, revenues are generated by ticket sales of vehicles and passengers along with some complimentary revenues [32].Expenses like crew salary [32], ship maintenance costs, and fuel charges affect the annual net income.Diesel fuel prices are constantly fluctuating in the market and therefore prices before 2022 have been sourced from Transtejo Soflusa (TTSL) (ship operator, Lisbon, Portugal) reports [23,26].For prices beyond 2022, they are assumed and calculated based on forecasts made by Det Norske Veritas (DNV) in their maritime forecast to 2050 as shown in Figure 15.
The net income of a project depends on revenues generated and on the expenses.Fo the ship Almadense, revenues are generated by ticket sales of vehicles and passenger along with some complimentary revenues [32].Expenses like crew salary [32], ship maintenance costs, and fuel charges affect the annual net income.Diesel fuel prices ar constantly fluctuating in the market and therefore prices before 2022 have been sourced from Transtejo Soflusa (TTSL) (ship operator, Lisbon, Portugal) reports [23,26].For price beyond 2022, they are assumed and calculated based on forecasts made by Det Norsk Veritas (DNV) in their maritime forecast to 2050 as shown in Figure 15.For Case I, green hydrogen was used as it is accepted to be the most environmentally friendly.The ship operates in Portugal; therefore, the prices considered are also associated with Portugal according to the research results of PwC until the year 2050 [36] (see Figur 16).For years beyond 2050, linear forecasting was used.For Case I, green hydrogen was used as it is accepted to be the most environmentally friendly.The ship operates in Portugal; therefore, the prices considered are also associated with Portugal according to the research results of PwC until the year 2050 [36] (see Figure 16).For years beyond 2050, linear forecasting was used.Grid electricity is needed to re-charge the battery while in port and it will be during the off-time operation, i.e., during the nights.Since the requirement of el begins in the year 2026, it requires forecasts.In this study, forecasts prepared by Brainpool in their research for the average baseload price for EU27 countries are ered [37].Beyond 2050, the last trend is followed and kept until 2055 (see Figure Annual profit and loss accounts, annual balance sheet, and statement cash fl generated for all the three case scenarios for the ship's lifetime of 30 years.Fuel e after the calculations emerge as the key factor that affects the economic analysis ap CAPEX and OPEX.Grid electricity is needed to re-charge the battery while in port and it will be charged during the off-time operation, i.e., during the nights.Since the requirement of electricity begins in the year 2026, it requires forecasts.In this study, forecasts prepared by Energy Brainpool in their research for the average baseload price for EU27 countries are considered [37].Beyond 2050, the last trend is followed and kept until 2055 (see Figure 17).
Annual profit and loss accounts, annual balance sheet, and sta generated for all the three case scenarios for the ship's lifetime of 3 after the calculations emerge as the key factor that affects the econom CAPEX and OPEX.To make a sensible conclusion, three base scenarios or mode first scenario considers the ship to continue to use a diesel engine its remaining life.The second and third scenarios incorporate the source.The difference between the latter two is that, in the second m is retrofitted with a FC that completely replaces the diesel engine, new ship is assumed to be operating on FC power in the third m construction involves a substantial timeline.To enhance the practic presumed that the ship order is commissioned towards the end of accounts for all the construction-related formalities (involving act Model 2, the ship is assumed to run as usual until 2026.Both the hybrid FC system; the calculations are carried out by changing FC Annual profit and loss accounts, annual balance sheet, and statement cash flows are generated for all the three case scenarios for the ship's lifetime of 30 years.Fuel expenses after the calculations emerge as the key factor that affects the economic analysis apart from CAPEX and OPEX. To make a sensible conclusion, three base scenarios or models are considered.The first scenario considers the ship to continue to use a diesel engine (business as usual) for its remaining life.The second and third scenarios incorporate the FC as the main power source.The difference between the latter two is that, in the second model, the current ship is retrofitted with a FC that completely replaces the diesel engine, whereas a completely new ship is assumed to be operating on FC power in the third model.In practice, ship construction involves a substantial timeline.To enhance the practicality of this study, it is presumed that the ship order is commissioned towards the end of 2023.A gap of 2 years accounts for all the construction-related formalities (involving actual construction).For Model 2, the ship is assumed to run as usual until 2026.Both the models consider the hybrid FC system; the calculations are carried out by changing FC modules to decide the optimal size.To select the economically viable solution, a balanced approach of the internal rate of return (IRR) and NPV is taken.These two approaches are considered to be the key performance indicators (KPI) which are calculated using Equations ( 17) and (18).
where F VAR is the financing cash flow for the year a, d is the discount rate, NPV is the net present value, and IRR is the internal rate of return of the project.F VAR is calculated using three factors: the financing cash flow, operating cash flow, and investment cash flow.The financing and investment cash flows are the outflows of cash whereas the operating cash flow is the inflow of cash into the company.Figure 18 represents the NPV of Models 2 and 3 for varying FC modules, where it can be seen that the maximum value of the NPV for Model 2 is achieved at the 7 FC module configuration whereas, for Model 3, 8 FC modules lead to a higher NPV (although with a very small margin between 7 and 8 modules); 8 FC modules have an edge but with only a 0.03% increase in the NPV with regards to 7 FC modules.On the other hand, investment increases by 2.35% with regards to the 7-module configuration.
net present value, and  is the internal rate of return of the project. is calculated using three factors: the financing cash flow, operating cash flow, and investment cash flow.The financing and investment cash flows are the outflows of cash whereas the operating cash flow is the inflow of cash into the company.
Figure 18 represents the NPV of Models 2 and 3 for varying FC modules, where it can be seen that the maximum value of the NPV for Model 2 is achieved at the 7 FC module configuration whereas, for Model 3, 8 FC modules lead to a higher NPV (although with a very small margin between 7 and 8 modules); 8 FC modules have an edge but with only a 0.03% increase in the NPV with regards to 7 FC modules.On the other hand, investment increases by 2.35% with regards to the 7-module configuration.For IRR, similar conclusions are obtained, i.e., the 7 FC module configuration leads to the best result.The change in IRR is also similar to the previous case as shown in Figure 19.For IRR, similar conclusions are obtained, i.e., the 7 FC module configuration leads to the best result.The change in IRR is also similar to the previous case as shown in Figure 19.

Selected Configuration Results
With the simulation results, it is validated that, when the ship encounters situations where it fails to reach the specified desired advance speed due to power limitations of the system, the speed is constrained to the threshold where the system attains maximum power under the prevailing conditions.
The hydrogen storage tank is based on the value obtained through simulation and it is designed for the highest consumption.A total of 440.1 kg-H2 is consumed during an average day in November, making this value the highest hydrogen consumption in a year.The tank is oversized by 20%, with a value of 5.50 kg mass.Figure 20 shows the consumption of hydrogen from the starting of its operation day (5:31 a.m.) until the end (09:59 p.m.).

Selected Configuration Results
With the simulation results, it is validated that, when the ship encounters situations where it fails to reach the specified desired advance speed due to power limitations of the system, the speed is constrained to the threshold where the system attains maximum power under the prevailing conditions.
The hydrogen storage tank is based on the value obtained through simulation and it is designed for the highest consumption.A total of 440.1 kg-H 2 is consumed during an average day in November, making this value the highest hydrogen consumption in a year.The tank is oversized by 20%, with a value of 5.50 kg mass.Figure 20 shows the consumption of hydrogen from the starting of its operation day (5:31 a.m.) until the end (09:59 p.m.).

Selected Configuration Results
With the simulation results, it is validated that, when the ship encounters situations where it fails to reach the specified desired advance speed due to power limitations of the system, the speed is constrained to the threshold where the system attains maximum power under the prevailing conditions.
The hydrogen storage tank is based on the value obtained through simulation and it is designed for the highest consumption.A total of 440.1 kg-H2 is consumed during an average day in November, making this value the highest hydrogen consumption in a year.The tank is oversized by 20%, with a value of 5.50 kg mass.Figure 20 shows the consumption of hydrogen from the starting of its operation day (5:31 a.m.) until the end (09:59 p.m.).The battery capacity selection is also carried out in a similar way and the highest value achieved is 369.8 kWh.The battery capacity is obtained using the Equation (19) given below: where  is the installed battery capacity,  is the energy consumed (369.8 kWh),  is assumed DoD of battery (50%), and  is the assumed total battery efficiency.In this case, the battery is designed to operate between 90% and 40% of its total capacity.This range of operation increases the lifetime of a battery and keeps it healthier The battery capacity selection is also carried out in a similar way and the highest value achieved is 369.8 kWh.The battery capacity is obtained using the Equation ( 19) given below: where C BAT is the installed battery capacity, E BAT CON is the energy consumed (369.8 kWh), k DoD is assumed DoD of battery (50%), and η b is the assumed total battery efficiency.In this case, the battery is designed to operate between 90% and 40% of its total capacity.This range of operation increases the lifetime of a battery and keeps it healthier for a longer period of time.Figure 21 shows the variation in battery energy during its operation day for an average day of each month of the year.for a longer period of time.Figure 21 shows the variation in battery energy during its operation day for an average day of each month of the year.

Environmental Pollution
Environmental protection is the main topic of this study.The combustion of each tonnage of oil in the engine produces 3.206 tons of CO2, whereas FCs do not produce any emissions during their operation.Emissions generated by the grid electricity to charge the battery depend on the energy mix of the grid.More than 67% of the electricity came from

Environmental Pollution
Environmental protection is the main topic of this study.The combustion of each tonnage of oil in the engine produces 3.206 tons of CO 2 , whereas FCs do not produce any emissions during their operation.Emissions generated by the grid electricity to charge the battery depend on the energy mix of the grid.More than 67% of the electricity came from renewable sources between 1 January and 30 September 2023, in Portugal as seen in Figure 22 [38].

Environmental Pollution
Environmental protection is the main topic of this study.The combustion of each tonnage of oil in the engine produces 3.206 tons of CO2, whereas FCs do not produce any emissions during their operation.Emissions generated by the grid electricity to charge the battery depend on the energy mix of the grid.More than 67% of the electricity came from renewable sources between 1 January and 30 September 2023, in Portugal as seen in Figure 22 [38].Portugal's government has made a firm commitment to achieve net-zero grid electricity carbon emissions by 2045.Therefore, as a result, a linear reduction in carbon emissions from 2020 to 2045 has been employed for emission calculations (see Figure 23).Portugal's government has made a firm commitment to achieve net-zero grid electricity carbon emissions by 2045.Therefore, as a result, a linear reduction in carbon emissions from 2020 to 2045 has been employed for emission calculations (see Figure 23).Through simulations of Cases I and II for Models 1, 2, and 3, fuel consumptions were recorded for every month of a year.To compare the lifetime emissions of each model, yearly emissions were considered, and models were compared by their lifetime emissions as shown in Table 8.Model 3 evidently reduced emissions to almost 97% whereas Model 2 reduced emissions by nearly 47% as compared to Model 1.  Through simulations of Cases I and II for Models 1, 2, and 3, fuel consumptions were recorded for every month of a year.To compare the lifetime emissions of each model, yearly emissions were considered, and models were compared by their lifetime emissions as shown in Table 8.Model 3 evidently reduced emissions to almost 97% whereas Model 2 reduced emissions by nearly 47% as compared to Model 1.

Conclusions and Future Work
Fuel cells represent an innovative technology aiding the decarbonization of the shipping sector.This paper studied the best possibilities of achieving the FC-driven ferryboat ship on a specified route by examining three options.The simulation outcomes of the chosen FC setup revealed that amongst the external resistive forces affecting the ship's power demand, wave and air drag exert the most significant influence.These forces interchange their dominant roles according to varying months.During the summer months, marked by low wave heights, the hull resistance prevails, while, in the winter months, the wave force becomes dominant.Consequently, the model produces diverse outcomes when tested across different seasons, illustrating the seasonal impact on component sizing.
The simulation performance shows the ability of a hybrid FC model to reduce carbon emissions to just 3% for Model 3 and 53% for Model 2 as compared with Model 1. Economically, the CAPEX of the battery is one of the key influencing factors (apart from fuel costs) as it is a recurring investment where it contributes to 3% and 7.7% for Models 2 and 3, respectively.The upfront FC hybrid system cost along with hydrogen fuel and electricity costs must plunge down to give preference over conventional engines.An increase in FC modules lead to an increase in capital cost, but this cost factor alone does not solely dictate the selection of a technology.The cost of hydrogen as a fuel decreases with the increase in FC modules.This affects the configuration selection, since this cost denotes around 61% of all expenses.With the economic selection criteria, 7 FC modules with a battery size of 770 kWh for Model 3 exhibits the best result with an IRR of 15.63% and an NPV of EUR 16.1 million.On the other hand, Model 2 results makes more sense, as replacing the relatively new ship is a difficult decision.The best way to achieve these results in the real world largely depends on government policies, technology advancements, and where investors interests are, with investors' interest being directly related to government policies.
Future work includes the enhancement of calculations in order to match the Energy Efficiency eXisting ship Index (EEXI) requirements, considering the degradation of the equipment involved over the period of time and the analysis of possible on-board renewable re-charging facilities.The designed simulation model can also be adapted to analyze different ships operating in different locations, with changing only a few parameters and blocks.In order to perform a similar economic analysis, carbon taxes are needed to be accounted for when considering bigger ships.The study can be enhanced to obtain the overall emissions via fuel that includes well-to-pump emissions (emissions that have already been emitted through the production and transportation stages).

Figure 1 .
Figure1.Ragone plot adapted from[14] for various types of emerging green technologies to power the energy system based on their power and energy densities.

Figure 2 .
Figure 2. The ship's route in Tagus River (the red line indicates the followed regular route whereas the yellow line represents a less frequently travelled route).

Figure 2 .
Figure 2. The ship's route in Tagus River (the red line indicates the followed regular route whereas the yellow line represents a less frequently travelled route).

Figure 2 .
Figure 2. The ship's route in Tagus River (the red line indicates the followed regular route wh the yellow line represents a less frequently travelled route).

Figure 3 .
Figure 3.The hybrid power system designed in MATLAB SIMULINK ("*" suffix denotes en providers output before converter) for the ferryboat Almadense (constructed by Navalria ship Grupo Martifer-Aveiro, Portugal) operating en route to Porto Brandão, Belém, and Trafaria in bon.

J
. Mar. Sci.Eng.2023, 11, x FOR PEER REVIEW 6 of 23 environmental conditions of the ship's operation vary.This variation changes the power demand even for the same velocity demand.

Figure 4 .
Figure 4.The speed-time variation of the ship Almadense for an average day of operating en route to Porto Brandão, Belém, and Trafaria in Lisbon.

Figure 4 .
Figure 4.The speed-time variation of the ship Almadense for an average day of operating en route to Porto Brandão, Belém, and Trafaria in Lisbon.

Figure 4 .
Figure 4.The speed-time variation of the ship Almadense for an average day of operating en route to Porto Brandão, Belém, and Trafaria in Lisbon.

Figure 5 .
Figure 5. Ship power demand profile.Figure 5. Ship power demand profile.

Figure 6 .
Figure 6.Assumed FC characteristic curve to be used in this study.Figure 6. Assumed FC characteristic curve to be used in this study.

Figure 6 .
Figure 6.Assumed FC characteristic curve to be used in this study.

Figure 7 .
Figure 7. Assumed FC module efficiency curve used in this study.

Figure 7 .
Figure 7. Assumed FC module efficiency curve used in this study.

Figure 8 .
Figure 8. Individual calculation blocks of all the forces acting on the ship's body; namely, in SIM-ULINK.

Figure 8 .
Figure 8. Individual calculation blocks of all the forces acting on the ship's body; namely, in SIMULINK.

Figure 8 .
Figure 8. Individual calculation blocks of all the forces acting on the ship's body; namely, in SIM-ULINK.

Figure 9 .
Figure 9.A vector representation showing the resultant vector in relation to north.

Figure 10 .
Figure 10.Hull resistances for different ship speeds and drafts.

Figure 10 .
Figure 10.Hull resistances for different ship speeds and drafts.
this study to operate in port and switch to ideal conditions (air purging).The EMS is responsible for deciding the energy distribution as shown in the decision flow diagram Figure 11.The FC capacity depends on the number of modules that change for each simulation.

Figure 11 .
Figure 11.EMS decision flow diagram for the selected case.

Figure 11 .
Figure 11.EMS decision flow diagram for the selected case.

Figure 12 .
Figure 12.Case 2 principle schematic diagram representing the flow of signals and variables required for the simulation.

Figure 12 .
Figure 12.Case 2 principle schematic diagram representing the flow of signals and variables required for the simulation.

Figure 13 .
Figure 13.Hydrogen consumption in kg for an average operating day of the month throughout the year with varying FC module numbers (200 kW unity).Figure 13.Hydrogen consumption in kg for an average operating day of the month throughout the year with varying FC module numbers (200 kW unity).

Figure 13 .
Figure 13.Hydrogen consumption in kg for an average operating day of the month throughout the year with varying FC module numbers (200 kW unity).Figure 13.Hydrogen consumption in kg for an average operating day of the month throughout the year with varying FC module numbers (200 kW unity).

Figure 13 .Figure 14 .
Figure 13.Hydrogen consumption in kg for an average operating day of the month throughout the year with varying FC module numbers (200 kW unity).

Figure 14 .
Figure 14.The force profile of all the external forces acting on the ship for the month (a) of July and (b) of November.Figure 14.The force profile of all the external forces acting on the ship for the month (a) of July and (b) of November.

Figure 15 .
Figure 15.The MGO price forecast by DNV from 2011 to 2050 in EUR per liter [35].

Figure 15 .
Figure 15.The MGO price forecast by DNV from 2011 to 2050 in EUR per liter [35].

Figure 16 .
Figure 16.The green hydrogen price forecast from 2026 to 2055 in EUR/kg for Portugal [36

Figure 16 .
Figure 16.The green hydrogen price forecast from 2026 to 2055 in EUR/kg for Portugal [36].

Figure 18 .
Figure 18.(a) The Model 2 NPV of business over its lifetime for different configurations of FC modules; (b) the Model 3 NPV of business over its lifetime for different configurations of FC modules.

Figure 18 .
Figure 18.(a) The Model 2 NPV of business over its lifetime for different configurations of FC modules; (b) the Model 3 NPV of business over its lifetime for different configurations of FC modules.

Figure 19 .
Figure 19.Models 2 and 3 IRR for the different FC configurations used to select the best solution.

Figure 19 .
Figure 19.Models 2 and 3 IRR for the different FC configurations used to select the best solution.

Figure 20 .
Figure 20.Hydrogen tank variation contour graph for the average operation day of each month of the year.

Figure 20 .
Figure 20.Hydrogen tank variation contour graph for the average operation day of each month of the year.

Figure 21 .
Figure 21.Battery capacity variation contour graph for the average operation day of each month of the year.

Figure 21 .
Figure 21.Battery capacity variation contour graph for the average operation day of each month of the year.

Figure 21 .
Figure 21.Battery capacity variation contour graph for the average operation day of each month of the year.

Table 4 .
Hydrogen fuel consumption per year for different FC module configurations in kg.

Table 5 .
Energy consumed by the battery in a year with varying FC module configurations.

Table 6 .
Market parameter values for economic viability analysis.

Table 8 .
The CO2 emission results of all models considered in this study over the lifetime of their operation in tons.

Table 8 .
The CO 2 emission results of all models considered in this study over the lifetime of their operation in tons.