Optimal Planning of Integrated Nuclear-Renewable Energy System for Marine Ships Using Artificial Intelligence Algorithm

Ocean-going ships are one of the primary sources of Greenhouse Gas (GHG) emissions. Several actions are being taken to reduce the GHG emissions from maritime vessels, and integration of Renewable Energy Sources (RESs) is one of them. Ocean-going marine ships need a large amount of reliable energy to support the propulsive load. Intermittency is one of the drawbacks of RESs, and penetration of RESs in maritime vessels is limited by the cargo carrying capacity and usable area of that ship. Other types of reliable energy sources need to be incorporated in ships to overcome these shortcomings of RESs. Some researchers proposed to integrate fossil fuel-based generators like diesel generators and renewable energy in marine vessels to reduce GHG emissions. As the penetration of RESs in marine ships is limited, fossil fuel-based generators provide most of the energy. Therefore, renewable and fossil fuel-based hybrid energy systems in maritime vessels can not reduce GHG emissions to the desired level. Fossil fuel-based generators need to be replaced by emissions-free energy sources to make marine ships free from emissions. Nuclear energy is emissions-free energy, and small-scale nuclear reactors like Microreactors (MRs) are competent to replace fossil fuel-based generators. In this paper, the technical, environmental, and economic competitiveness of Nuclear-Renewable Hybrid Energy Systems (N-R HES) in marine ships are assessed. The lifecycle cost of MR, reliability of the proposed system, and limitations of integrating renewable energy in maritime vessels are considered in this study. The proposed N-R HES is compared with three different energy systems, namely ‘Standalone Fossil Fuel-based Energy Systems’, ‘Renewable and Fossil Fuel-based Hybrid Energy Systems’, and ‘Standalone Nuclear Energy System’. The cost modeling of each energy system is carried out in MATLAB simulator. Each energy system is optimized by using the Differential Evolution Algorithm (DEA), an artificial intelligence algorithm, to find out the optimal configuration of the system components in terms of Net Present Cost (NPC). The results determine that N-R HES has the lowest NPC compared to the other three energy systems. The performance of the DE algorithm is compared with another widely accepted artificial intelligence optimization technique called ‘Particle Swarm Optimization (PSO)’ to validate the findings of the DE algorithm. The impact of control parameters in the DE algorithm is assessed by employing the Adaptive Differential Evolution (ADE) algorithm. A sensitivity analysis is carried out to assess the impact of different system parameters on this study’s findings.


Introduction
If global shipping were a country, it would be considered as the sixth-largest producer of CO 2 after the United States, China, Russia, India, and Japan. Ocean-going shipping is responsible for more than 3% of global CO 2 emissions [1]. Although international shipping is not a part of the Paris Agreement, the International Maritime Organization (IMO) is

System Modeling
An oil tanker called 'Baltic Sunrise' is studied in this paper. Automatic Identification System (AIS) data of this ship is collected from a vessel data collection company named 'FleetMon'. There are MRs, wind turbines, solar PV panels, energy storage systems, controllers, and converters in the proposed system. Estimation of the energy demand of the ship is calculated by using the AIS data. The system's project lifetime is 40 years with the inflation rate and the discount rate as 2% and 8%, respectively. Renewable energy resource data is collected from the NASA Surface meteorology and Solar Energy database. The study considers the lifecycle cost of MR and diesel generators to make the analysis more practical. This research also addresses the system's reliability and constraints of implementing energy systems in marine vessels due to the limited area and cargo carrying capacity of the ship. Four energy systems, namely 'Nuclear Renewable Hybrid Energy Systems (N-R HES)', 'Standalone Fossil Fuel-based Energy System', 'Fossil Fuel and Renewable Hybrid Energy System', and 'Standalone Nuclear Energy System', are studied here. Matlab software is used to model the cost of each energy system. All the systems are then optimized by using the Differential Evolution Algorithm (DEA) to find out the optimal system configuration with minimal NPC. The study compares all the optimized systems to identify the best energy system for maritime transportation. The Particle Swarm Optimization (PSO) is used to measure the performance of the Differential Evolution Algorithm (DEA), and validate the findings of DE optimization technique. To analyze the

System Modeling
An oil tanker called 'Baltic Sunrise' is studied in this paper. Automatic Identification System (AIS) data of this ship is collected from a vessel data collection company named 'FleetMon'. There are MRs, wind turbines, solar PV panels, energy storage systems, controllers, and converters in the proposed system. Estimation of the energy demand of the ship is calculated by using the AIS data. The system's project lifetime is 40 years with the inflation rate and the discount rate as 2% and 8%, respectively. Renewable energy resource data is collected from the NASA Surface meteorology and Solar Energy database. The study considers the lifecycle cost of MR and diesel generators to make the analysis more practical. This research also addresses the system's reliability and constraints of implementing energy systems in marine vessels due to the limited area and cargo carrying capacity of the ship. Four energy systems, namely 'Nuclear Renewable Hybrid Energy Systems (N-R HES)', 'Standalone Fossil Fuel-based Energy System', 'Fossil Fuel and Renewable Hybrid Energy System', and 'Standalone Nuclear Energy System', are studied here. Matlab software is used to model the cost of each energy system. All the systems are then optimized by using the Differential Evolution Algorithm (DEA) to find out the optimal system configuration with minimal NPC. The study compares all the optimized systems to identify the best energy system for maritime transportation. The Particle Swarm Optimization (PSO) is used to measure the performance of the Differential Evolution Algorithm (DEA), and validate the findings of DE optimization technique. To analyze the impact of the control parameters of the DE, Adaptive Differential Evolution (ADE) is employed in this study. A sensitivity analysis is accomplished to see the impact of different system parameters.  [15].
The effective ship propulsive power is calculated in each position using the speed of the ship in that position, draught, breadth, residual resistance, and frictional resistance. ITTC-57 method and Gertler series data chart are used in this calculation. Appendix A summarizes the parameters and assumptions that are considered during the calculation of the estimation of ship electric power. Detail calculation of estimating the energy demand of the ship can be found in [16][17][18].
The maximum rating of propulsive electric power of this ship is 24.64 MW. Another type of electric power, namely 'Auxiliary Load', is considered along with propulsive power. Auxiliary load emulates the electrical load demand for light, cooling, feed systems, pumps, and air compressors. Grzegorz et al. [19] assessed the auxiliary load to propulsive load ratio for different types of maritime vessels by collecting the data from ships register, logbook record and interviewing the captains of ships, chief engineers, and stations pilot. For oil tankers, they found the ratio as 0.211. Appendix B shows the yearly propulsive electric power (MW) and yearly total load demand with auxiliary load (MW).

Diesel Generator
Fossil Fuel-based Generator (FFG) is considered in this study for comparing the proposed N-R HES with the conventional energy systems of ocean-going marine ships. A marine vessel's traditional energy system consists of a main engine, auxiliary engine, heat recovery steam generator, and shaft generator [5]. The cost of a diesel generator can be categorized broadly into three major categories-'Capital Cost', 'Operating Cost', and 'Decommissioning Cost'. Since the decommissioning cost of a diesel generator is minimal, it can be ignored.

Solar Energy
Solar PV systems convert sunlight into electricity. The generation of solar power from solar PV depends on the ambient temperature, solar irradiance, and solar PV area. The output power from solar PV can be calculated by using the following equation [21].
where, P SPV (t) is the output power from the solar PV (W), Num SPV is the number of solar PV, n SPV (t) is the instantaneous efficiency of the solar PV, PV A is the area occupied by the solar PV (5 m 2 ), and R t (t) is the solar irradiance (W/m 2 ). The instantaneous efficiency of the solar PV can be obtained by the following equation [22]. n SPV (t) = n pv−re f × n MPPT × 1 + β T c (t) − T c−re f ∀t ∈ t total (2) where, n SPV (t) is the instantaneous efficiency of the solar PV (%), n pv−re f is the reference efficiency of the solar PV (%), n MPPT is the efficiency of the Maximum Power Point Tracking (MPPT) unit which is 1 in this study, T c (t) is the temperature of the PV cell ( • C) at time step t, T c−re f is the reference temperature of the PV cell, and β is the temperature coefficient of efficiency ( • C −1 ). The PV cell temperature can be obtained by the proposed PV model by Markvar [23].
The hourly data for one year holds a total of 8760 data points. The solar radiation differs in the shipping route with longitude and latitude, time, and date. To keep the calculation simple, the data solar irradiance and temperature are considered based on the monthly position of the ship in the adjacent ports. Appendix C represents the hourly solar irradiance and temperature for one year, respectively. The usable area of the marine ship for solar PV installation is limited. Solar PV systems must not hinder the transfer area and must not conceal any area with financial impact [28]. This study addresses this constraint while planning the optimal configuration of the system components.

Wind Power
The wind is a great source of energy and has the potential to generate electricity without any fuel. Several models can be used to estimate the output power of the wind turbine, like the model based on Weibull parameters [29,30], the linear model [31,32], and the quadric model. In this study, the quadric model is used. To calculate the output power of the wind turbine, the wind speed at hub height needs to be obtained by the following equation.
where, V whub (t) is the wind speed at the hub height (m/s) at time step t, V wre f (t) is the wind speed at the anemometer height (m/s) at time step t, H wre f is the anemometer height (m), H whub is the hub height (m), and µ is the power-law coefficient of the wind speed. After measuring the wind speed at the hub height, it is used to calculate the wind turbine power by using the following equations [33].
where, P W (t) is the output power of the wind turbine (kW) at time step t, P r is the rated power of the wind turbine (kW), V r is the rated wind speed of the wind turbine (m/s), V cin is the cut-in speed of the wind turbine (m/s), V cout is the cut-out wind speed (m/s), and V whub is the measured wind speed at the hub height (m/s  [34], [35], [36], [37], [38], [39], [40].

vSMR/MR
This study evaluates the lifecycle cost of vSMR and focuses on the economic model of it. SMR is a fourth-generation nuclear reactor having power equivalent to 300 MWe or less. Modular fabrication, passive safety systems, and enhanced protection against man-made or natural hazards are some advanced features of SMR. A subcategory of SMR is vSMR, which power rating is less than 15 MWe [41]. Another subcategory of SMR is Microreactor (MR), which typically ranges from 1 MWe to 50 MWe [42].
Based on the lifecycle phases of SMR and nuclear fuel, the cost of vSMR can be divided into three major categories that are capital costs, annual operating costs, and backend costs [20]. The first deployment cost of any technology is always higher compared to the subsequent deployment. The rate of cost reduction of the next deployment depends on the lesson learned from the prior deployment. The cost reduction depends on the learning rate of the factories. The learning rate refers to the cost reduction of a factory that can be achieved by experience and learning how to execute a process and use different tools to deliver a product. The learning rate of SMR could be 8% to 10%. The relation between learning rate and cost reduction can be expressed by the following equation [16].
where, R is the rate cost reduction, and LR is the learning rate. In this study, the learning rate is considered as 10% [13]. The relation between the vSMR capital cost and the rate of cost reduction can be expressed by the following equation.
where, vSMRC u is the unit cost of N th vSMR ($), vSMRC 1st is the 1st vSMR cost ($), and R is the rate of cost reduction. vSMR is capable to perform both baseload operation and load-following operation. However, the SMR vendor outlines that load following feature of SMR affects the main steam supply and coolant system components and frequent replacement might be required. During the Extended Reduced Power Operation (ERPO), reactors can operate as low as 30% of its' rated power for several days to several months. In this study, the baseload operation of vSMR is considered for all the time except the time of loading and unloading the cargo/liquid (oil) in the port when ERPO (30% of the rated power) is considered. Although there are no CO 2 emissions during the operation of SMR, CO 2 emits in the lifecycle of SMR, mainly in the construction stage due to concrete production and in the fuel processing stage due to the high demand for electricity and fossil fuel. The Reactor size, Plant lifetime, Overnight capital cost, Fixed O&M cost, Fuel cost, Refueling cost of fuel module, Core lifetime, Decommissioning cost, Capacity factor, Plant efficiency, and CO 2 emission are 1000 (kWe), 40 (Years), 15,000 ($/kWe), 350 ($/kWe), 10 ($/MWh), 20 (million $), 10 (Years), 5 ($/MWh), 95 (%), 40 (%), and 4.55 (kg/MWh), respectively [20], [43].

Electrochemical Energy Storage
The EES is considered in this study to store the excess energy in the form of chemicals and supply when it requires. The rating of the EES depends on the rating of the load, autonomy hour, depth of discharge of the battery, the efficiency of the battery, and the inverter. The rating (kWh) of the battery can be expressed by the following equation [44].
where, P load is the rating (kW) of the load, AH is the autonomy hours, N bat is the number of battery, DOD bat is the depth of discharge of the battery, n inv is the efficiency of the inverter, and n bat is the efficiency of the battery. In this study, the load is considered as 70 kW [45]. The autonomy hour refers to the required amount of time to serve the load continuously in an emergency. The autonomy hour is between 12 h to 18 h [46]. In this paper, autonomy hour is taken as 18 h. The DOD bat , n inv , and n bat is 80% [46], 95% [47], and 85% [48], respectively. The state of charge (SOC) of the battery bank depends on the DOD and rating of the battery. The maximum and minimum SOC of the battery can be calculated by the following equations.
BAT SOC,max = BAT rating × DOD bat (10) where, BAT SOC,max is the maximum SOC of the battery bank, and BAT SOC,min is the minimum SOC of the battery bank. The charging and discharging of the battery bank at any time step depend on the available energy that can be given to the battery bank, the amount of energy that is needed from the battery bank, and the SOC of the battery bank in the previous time step. BAT SOC,min is considered as the initial SOC of the battery in this paper. The energy that can be given to the battery and the energy that is needed from the battery bank depend on the total energy generation and total load demand of the system. The below equations are used to calculate the amount of energy that can be given and the amount of energy that is needed from the battery bank.
P gen (t) = P pv (t) + P w (t) + P GNRT (t) ∀t ∈ t total (14) where, BAT inm (t), BAT outm (t), P genm (t), P L (t), P pv (t), P w (t), P vSMR (t), P GNRT (t) is the amount of energy that can be given to the battery bank, the amount of energy that is needed from the battery bank, the total electric energy generation, the total electric load, the total electric energy from the solar panel, the total electric energy from the wind turbine, the total electric energy from MR, and the total electric energy from diesel generator at the time step (t), respectively. The amount of energy that the battery bank can take and give at any time step can be formulated by the below equations.
where, BATe(t − 1) is the battery energy at time step (t − 1), BAT take (t) is the amount of energy that can be taken by the battery bank at the time step (t), and BAT give (t) is the amount of energy that can be given by the battery bank at the time step (t). The battery capacity after charging and discharging at time step (t) depends on the load demand, total energy generation, maximum SOC, minimum SOC, and battery capacity at time step (t − 1). The charging and discharging of the battery bank can be formulated by using the below equations BAT dischrg (t) = BAT SOC,min ; BAT e (t) = BATe(t − 1); P gen (t) = P L (t) ∀t ∈ t total (22) where, BAT chrg (t), BAT dischrg (t), and BAT e (t) refer to the battery capacity after charging, battery capacity after discharging, and battery capacity when energy generation and demand are exactly the same, respectively at time step (t). If there is any excess energy that will be dumped. The Capital Cost, Replacement Cost, O&M Cost, Lifetime, Efficiency, Depth of Discharge, and Inverter Efficiency are 398 ($/kWh), 398 ($/kWh), 10 ($/kWh/Year), 5 (Years), 85 (%), 80 (%), and 95 (%), respectively.

Key Performance Indicators
The KPIs are used to compare different energy systems to evaluate the feasibility. This study employs technical, economical, and environmental KPIs as follows.

Net Present Cost (NPC)/Net Present Value (NPV)
The Net Present Cost (NPC) and Net Present Value (NPV) are similar concepts and they only differ in sign. The NPV is the present value of all future cash flows (positive and negative) of an investment with a discount rate [49]. The lower NPC implies a higher profit to the investors. The following equation can be used to calculate the NPC. As NPV and NPC differ only in sign, NPC and NPV are related as follows.

Net Present Cost
where, rld, f, i, and t denote the real discount rate (%), nominal discount rate (8%), inflation rate (2%), and time periods number, respectively. In the nominal discount rate/interest rate, the inflation rate is covered. To factor out the impact of the inflation rate, the real discount rate is considered in this study. The equation for calculating the real interest rate is also known as Fisher equation [50]. For calculating the NPC for a longer project lifetime with multiple cash flows, the formula can be modified as per the below equation.
where R t and n refer to net cash inflow-outflow in a unit time period, and project lifetime, respectively.

Levelized Cost of Energy (LCOE)
The LCOE is the ratio of the total annualized system cost ($) to the total electricity (kWh) delivered by the system annually [51]. The lower LCOE refers to the lower unit price of the energy the system can serve which implies a higher profit to the investors. The LCOE is expressed by the following equation [13].
where, NPC total is the total NPC of the system ($). P gen (t) is the energy generation from the system at time step t, n is the project lifetime (40 years), and ∆t is the time step considered in the analysis.

Generation Reliability Factor (GRF)
GRF is one of the reliability factors that indicate the amount of demand of a system is accomplished by the energy system [14]. The mathematical formula of GRF for the demand of the system can be expressed by the following equation [13].
where, GRF e is the generation reliability factor. P gen (t) and P L (t) represent the power generation and electric load demand at time step t. ∆t shows the time step considered in the calculation. The higher value of GRF implies the system is more reliable.

Loss of Power Supply Probability (LPSP)
The LPSP is another reliability factor that determines the amount of unserved energy of a system when the system demand is higher than the generation. To ensure the reliability of the system, it is important to maintain the LPSP within a boundary. Therefore, this KPI is considered as one of the constraints in the optimization of the system. The LPSP is determined for each time step and for the total time step to assure the utmost reliability of the system. The following equations are used to determine the LPSP [13].
where, LPSP e represents the loss of power probability for the total time step, LPSP e (t) shows the loss of power probability at each time step. The maximum and minimum value of LPSP is 1 and 0, respectively. The lower LPSP infers the higher reliability of the system.

CO 2 Gas Emissions
Several pollutants are emitted from the energy generators namely Carbon Monoxide, Particulate matter, Nitrogen Oxides, Carbon Dioxide, Sulfur Dioxide, and Unburned Hydrocarbons (UHC) in the lifecycle of the generators. In this study, only CO 2 gas emission is considered. The amount of CO 2 emissions from any generators is calculated by the following equation [52].
where, TAEG is the total annual electricity generation from the generators. The emissions factor for diesel generators and SMR is 700 kg/MWh and 4.55 kg/MWh, respectively [20], [53]. The penalty for CO 2 emissions is considered during the calculation of NPC. The annual penalty for CO 2 emissions can be calculated by using the following equation.
where, CCE is the cost associated with CO 2 emissions penalty ($), ACE is the annual CO 2 emissions (tonne), and CEP is the CO 2 emissions penalty ($/tonne). Although nation-wise carbon tax varies a lot, IMF calculated large CO 2 emitter countries should be charged between 50 USD/ton to 100 USD/ton by 2030 to fulfill their commitments to lessen the carbon emissions [54]. For the shipping industry, IMF proposed 30 USD/ton CO 2 emissions as a penalty [55]. In this study 30 USD/ton CO 2 emissions is considered as CO 2 emissions penalty to show the impact on NPC and LCOE.

Power to Weight Ratio (PWR)
The PWR is an important criterion to choose the energy system for the marine ship. The PWR measures the performance of any power source or system. The heavyweight of power system of marine ships results in higher displacement and lower speed.
The higher PWR allows the marine ship to carry more load while maintaining the high speed [56]. The PWR of a power system can be calculated by using the following equation [57].

PWR =
Power rating o f the system (kW) Weight o f the system (kg) (33) PWR is considered as one of the constraints in the optimization problem. Appendix D shows the PWR of the different components of the system.

Problem Formation
This section covers the formation of the optimization problem including objective function, constraints, decision variables, and steps of DE, ADE, and PSO algorithm implementation.

Objective Function
The objective of the optimization problem is to identify the optimal configuration of the HES to achieve the lowest NPC. Along with the N-R HES, other three energy systems, namely 'Standalone Fossil Fuel-based Energy Sytems', 'Renewable and Fossil Fuel-based Hybrid Energy System', and 'Standalone Nuclear Energy System' are also optimized to compare with the proposed N-R HES. The objective function represents the economic KPI and the constraints of the optimization problem covering the technical and other KPIs, such as LPSP and PWR. The fitness function of the optimization problem is the total NPC of all the system components. The fitness function of the optimization problem can be formulated as below.
where, K is the system component set, and NPC j refers to the NPC of the jth component. NPC j = C capc,j + C repc,j + C O&Mc,j + C f uelc,j − C salvc,j ∀i ∈ G where, C capc,j , C repc,j , C O&Mc,j , C f uelc,j , and C salvc,j represent the present value of capital cost, replacement cost, O&M cost, fuel cost, and the salvage value of the jth component, respectively. The salvage value of a component refers to the remaining value of a component at the end of the project lifetime where the component is used [58]. The capital cost of any component is calculated at the beginning of the project. The total capital cost involves the capital cost of each component and the number of components. The capital cost can be calculated by using the following formula.
where, C capc,j , N com,j , and C capc,unit(j) refer to the capital cost, the number of components, and the unit cost of the jth component, respectively.
The capital cost of MRs is calculated in a different formula as the cost of reduction is included in the capital cost. The rate of cost reduction of MRs is related to the learning rate, which is discussed earlier. The capital cost of MR can be calculated by using the formula as follows [59].
where, C capc,MR , C capc,MR(1st) , N MR , and R imply the total capital cost of MR, the unit price of the 1st MR, the number of MR, and the rate of cost reduction, respectively. The O&M cost of a component occurs each year and continues till the end of the project lifetime. The O&M cost of a component is considered the same each year and is calculated by the following equation [60].
where, C O&Mc,j , N com,j , and C O&M,yearly(j) indicate the present worth of the total O&M cost, number of components, and yearly O&M cost of the jth component, respectively.
The replacement of any component happens at the end of its' lifetime. The number of replacements depends on the component lifetime and the lifetime of the project. The present value of the replacement cost of the component is calculated as follows [60].
where, C repc,j is the present value of the replacement cost of the jth component, NR is the number of required replacement, CLT j is the lifetime of the jth component, and C rep,unit(j) is the per-unit replacement cost of the jth component. ceil(X) is a function that rounds X to the nearest integer equal or greater than X. The fuel cost is considered for MRs and fossil fuel-based generators. There is no fuel cost associated with renewable energy sources like solar PV and wind turbines. Also, there is no fuel cost for energy storage devices. The fuel cost is calculated yearly, and the present value of the fuel cost is formulated by using the below equation [60].
where, C f uelc,j is the present worth of the fuel cost ($) of the jth component, E yearly,j is the yearly energy generation (MWh) from the jth component, CU f uel,j fuel price per unit energy generation of the jth component, and RAT j is the rating of the j 4 component. The salvage value of a component is the remaining useful value of a component at the end of the project lifetime. The linear depreciation is considered in the calculation of the salvage value. Linear depreciation refers that the salvage value of the component is directly proportional to the remaining life of the component. The present worth of the salvage value is calculated as follows [58].
where, C salv,j is the present worth of the salvage value of the jth component, and CLT rem,j is the remaining year of the jth component at the end of the project lifetime. f loor (X) is a function that rounds X to the nearest integer equal or less than X. The decommissioning cost and refueling cost are considered only for MRs. Although The decommissioning cost of MRs occurs at the end of the project, the annual distribution of decommissioning cost is considered in this study [43].
where, C decom,MR is the present value of the total decommissioning cost ($) of MR, E yearly,MR is the yearly energy generation (MWh) from the MR, and CU decom,MR is the decommissioning cost per unit ($/MWh). The lifetime of the fuel module is considered 10 years in this study. Hence, the refueling of the MR occurs every ten years. The refueling cost of MRs does not include the cost of fuel. The refueling cost includes the cost associated with transportation of the fuel, labor cost, etc. The present value of the refueling cost is calculated by the following equations [60].
where, C re f ueling,MR is the present value of the refueling cost of the MR, C re f ueling,MR(unit) is the refueling cost in every 10 years, LT f b(MR) is the lifetime of the fuel bundle, and MR re f uel is the number of the required refueling in the project lifetime.

Constraints
Constraints are used in the optimization problem to ensure the reliability, resiliency, and performance of the N-R HES. Some constraints are used to address some real-world limitations to implement the N-R HES in the marine ship.
The energy generation from any energy source must be equal to or less than the maximum capability of that energy source. This can be represented as follows.
where, P j gen (t) is the electric generation of the j th component at time step t and P j gen,max is the maximum electric generation capacity of the j th component at time step t.
To ensure the reliability and resiliency of the systems, the total generation of the system must be greater or equal to the total requirement of the system. The following equations are used to represent the constraints of the energy management of the system.
where, P z gen (t) is the generation of electric power (kW) at time t of the year z. P z EL (t) refers to the power demand at time step t of the year z. Num PV , Num WT , and Num MR indicate the number of solar PV, the number of wind turbine, and the number of MR, respectively. P pv (t), P w (t), and P MR (t) refer to the power generation from solar PV, wind turbine, and MR at the time t, sequentially.
The power generation for Renewable and Fossil Fuel-based Hybrid Energy Systems (F-R HES), Standalone Nuclear Energy System, and Standalone Fossil Fuel-based Energy System can be formulated by the Equation (54), (55), and (56), respectively.
The energy storage device is used to store excess energy and supply energy during generation shortage. The following constraints are considered in this study to ensure the proper energy management of the system and maintain the desired SOC of the storage device.
where, SOC z EES (t) is the SOC (%) of the energy storage device at the time step t in the year z. SOC EES,min and SOC EES,max refer to the minimum and maximum SOC of the energy storage component. This constraint will ensure that the usage of energy storage systems will be within the upper and lower limit which is important for the longevity of the energy storage systems (battery).
To ensure the reliability, resiliency, and performance of the energy system, the reliability parameter, LPSP, is considered as a constraint of the optimization problem. The LPSP is considered for each time step and also for the total time step to ensure the utmost reliability of the system. The lower value of the LPSP refers to the higher reliability of the system. The reliability constraints can be formulated as follows.
where, LPSP e,max is the maximum limit of the LPSP for the total time step and LPSP e(t),max represents the maximum limit of the LPSP in each time step. The maximum limit of LPSP e and LPSP e(t) is considered as 8% in this study, which is a typically acceptable margin for HES [48].
To ensure the cargo-carrying capacity of the marine ship, the ratio between the total weight of the energy system and the total cargo-carrying capacity of the ship is inserted as one of the constraints in the optimization problem. This constraint can be formulated by the following equation.
SYS wr ≤ SYS wr,max where, SYS wr is the ratio between the total weight of the energy system and the total cargocarrying capacity of the ship, and SYS wr,max is the maximum limit of the ratio. In this study, the SYS wr,max is taken as 0.005 considering the weight of the current energy system that is being used in the ship to ensure that the proposed system will not reduce the load-carrying capacity of the marine vessel. The total weight of the components of the N-R HES, F-R HES, Standalone Nuclear Energy System, and Standalone Fossil Fuel-Based Energy System can be calculated by using the Equation (61), (62), (63), and (64), respectively.
where, TW MR , TW PV , TW Wind , TW DG , and TW BAT represent the total weight of the MR, solar PV, wind turbine, diesel generator, and battery, respectively. The weight of the MR, diesel generator, solar PV, wind turbine, and the battery can be formulated by the Equation (65), (66), (67), (68), and (69), respectively.
where, PWR MR , PWR DG , PWR PV , PWR Wind , and PWR BAT refer to the power to weight ratio of MR, fossil fuel-based generator, solar PV, wind turbine, and battery, respectively. The Rat MR , Rat DG , Rat PV , Rat Wind , and Rat BAT refer to the rating of the MR, fossil fuelbased generator, solar PV, wind turbine, and battery, sequentially. Unlike land-based system, marine ship has a limited area to accommodate the energy system components. To ensure that the system components can be accommodated with the available useable area of the ship, another constraint is addressed in the optimization problem. Solar PV is installed on the upper surface of the ship and requires a considerable amount of area compared to other energy sources [61]. Therefore, this constraint is used in the optimization problem to confirm that the solar PV can be implemented within the area of the ship. This constraint can be formulated as below.
where, Solar ratio is the ratio between solar PV area and ship area, and Solar ratio,max is the maximum limit of the ratio between solar PV and ship area. The area of the solar PV can be calculated as follows.
where, Area PV is the area of the solar PV (m 2 ), and U A PV is the area requirement of solar PV for per unit power (m 2 /kW). The U A PV is considered as 4.952 (m 2 /kW) in this study [62]. Considering the area requirement for wind turbine installation and all the area of the marine ship is not practically available for solar PV, the Solar ratio,max is considered as 0.4 which is within the reasonable margin [22].

Decision Variables
The target of this optimization problem is to identify the optimal configuration of the energy system by minimizing the NPC of the systems. Four different energy systems are optimized in this research and compared with each other to find out the most feasible energy system for the marine ship. The decision variables of this optimization problem are the number of the required MRs, fossil fuel-based generators, solar PV, wind turbine, and battery. These decision variables are used to identify the size of the MR, fossil fuel-based generator, solar PV, wind turbine, and battery. The decision variables can be expressed as follows.
where, Num MR,max , Num DG,max , Num PV,max , Num WT,max , and Num BAT,max refer to the maximum number of MR, fossil fuel-based generator (Diesel Generator), solar PV, wind turbine, and battery, respectively.

Implementation of Optimization Algorithm: Differential Evolution
The differential evolution algorithm was first proposed by Kenneth Price and Rainer Store [63]. It is a popular optimization technique, which uses a population of individual solutions to solve multidimensional real-valued functions. The principle idea of the DE algorithm is to produce temporary individual best based on the individual differences within the population and then rebuild the population evolutionary. Because of its robustness, the suitability of different numerical optimization, and better global convergence, it has been used extensively in many problems like image classification, linear array, neural network, and monopoles antenna [64]. To overcome the low search ability and premature convergence of quantum inspired differential evolution (QDE), an enhanced MSIQDE algorithm was proposed in [65]. For solving Minimum Sum-of-Squares Clustering problem, Memetic Differential Evolution (MDE) was employed in [66]. The authors offered four algorithm variants that differ on the matching approach. They compared their method with HG-MEANS and found competitive performance based on numerical experiments. Sunita et al. employed DE with dynamic mutation technique to solve multi-objective optimization model of a microgrid consists of solar and wind energy. They found that the value of crossover probability does not affect the optimization result significantly. The authors concluded that for getting the best result, the crossover probability should be 0.5 [67]. To overcome the limitations of Harmony Search (HS) algorithm, an improved differential HS was employed in [68]. The authors used twenty-two benchmark functions to assess the performance of eight optimization algorithms. For setting the parameters of the optimization algorithms, they utilized the recommendations of several authors. They set the crossover probability and scaling factor of DE algorithm as 0.6 and 0.6. For solving the slow convergence speed, low solution efficiency, and lack of diversity of DE algorithm, Quantum Differential Evolution (QDE) was employed in [69]. The author set the value of parameters based on the literature review. They set the value of crossover probability and scaling factor for different DE algorithms ranging from 0.5 to 0.9. However, to evaluate the impact of crossover probability and scaling factor on optimization result, Adaptive Differential Evolution (ADE) algorithm is employed which is discussed in Section 4.5. The differential evolution has four major stages that are initialization, mutation, recombination/crossover, and selection.
In the initialization stage, random solutions are generated within a boundary. Then, each solution undergoes mutation and recombination phases, and lastly, the best solution is selected. This process takes place until it meets any stopping criteria. The program needs to generate the initial solutions randomly within the minimum and maximum value of each variable. Any random solution that is generated in the initialization phase can be formulated as follows [63].
where, N refers to the Nth solution to be optimized and d represents the parameter variables to be optimized. S (min,d) and S (max,d) represent the minimum and maximum range of the d parameter variables. The rand is a random number between 0 and 1.
In the mutation stage, the random solutions (target solutions) which are generated in the initialization stage are optimized by mutating the solution. This solution is called the donor solution. The new solution after the mutation phase can be formulated as follows [63].
where, r1, r2, and r3 are random solutions selected from 1 to N p, and these solutions are mutually exclusive. F represents the scaling factor, which is between 0 and 2. The target solution does not participate in the mutation phase. Hence, in the differential evolution, the minimum population size is four. After the mutation phase, the solution undergoes a crossover/recombination phase to add diversity to the solution. The solution following the crossover phase is called the trial solution. In this phase, the trial solution can be obtained either from the target solution or the donor solution based on the following condition [63].
where, C p is the crossover probability, and ϕ is the random number between 0 and 1. λ is a randomly selected variable location from d. The N = λ condition ensures that at least one solution in the crossover phase is taken from the donor solution. After getting the trial solution, its' boundary needs to be checked. If the trial solution is greater than the maximum limit then the solution is bounded back to the maximum value, and if the trial solution is less than the minimum value then it is bounded back to the minimum value.
In the selection stage, the objective function values are computed for all the solutions. The better solution, between the target solution and the trial solution, is selected for the next iteration; this is called the greedy selection method. After completing all the iteration, the best objective function value is obtained, and the associated solution is the global optimal solution. Differential Evolution implementation steps are discussed in detail as follows.
Step 1: Read the following input data of the proposed energy systems: A.
Load the electrical demand data (hourly propulsive electric and auxiliary load data for one year). B.
Load the climatological data (hourly wind speed, solar irradiance, and ambient temperature) to calculate the renewable energy where applicable. C.
Load the system component's features (energy generation from MR, fossil fuel-based generator, solar PV, wind turbine, and charging and discharg-ing from energy storage device). D.
Load economic and other parameters of the system and system compo-nents, such as project lifetime, inflation rate, discount rate, capital cost, O&M cost, replacement cost, fuel cost, refueling cost, decommissioning cost, and the salvage value.
Step 2: Initialize the parameters of DE and required system components: A. Set the population size, N p = 50, and the maximum number of iterations as 200. B.
Set the scaling factor, F = 0.5 C.
Set the crossover probability, C p = 0.7 Step 3: Generate the initial solutions randomly within the boundary of the decision variables. The number of the randomly generated solution is equal to the population size. These solutions are called target solutions.
Step 4: Generate the donor solutions by mutating the target solutions.
Step 5: Generate the trail solutions by recombing the donor the solutions (crossover).
Step 6: Bound the trial solutions within the boundary of the decision variables.
Step 7: Evaluate the objective function for the target solutions and trial solutions. The solution that gives the lower value (minimization problem) of the objective function is selected as the target solution for the next iteration. Store the best cost value (lowest value) in every iteration.
Step 8: When the number of iteration reaches the maximum limit of the iteration, then stop. Otherwise, continue from Step 3 to Step 7. After completing all the iteration, the best objective function value is obtained, and the associated solution is the global optimal solution.

Adaptive Differential Evolution
Currently, the adjustment of control parameters (scaling factor and crossover probability) of DE is considered a global behavior. However, these parameters affect the convergence speed and may cause a premature phenomenon. An adaptive adjustment method of the control parameters was discussed in [64]. The authors used sine and cosine functions with a value of (−1, 0), and (0, 1), respectively. The scaling factor (F) and the crossover probability (C p ) can be written as follows [64].
where, F (t) is the scaling factor, C p(t) is the crossover probability, t is the iteration number, and maxit is the maximum number of iterations. Y and β are constants, and the value of Y and β is 0.8 and 0.75, respectively, in this study. This study compares the performance of standard DE and ADE, and evaluates the impact of control parameters of the DE algorithm on finding the optimized result.

Implementation of Optimization Algorithm: Particle Swarm Optimization (PSO)
The PSO is another metaheuristic optimization algorithm, proposed by Kennedy and Eberhart in 1995. The main idea of the PSO algorithm is communication among the population called the swarm and, evolutionary computation. In this technique, two values are compared in each position of a particle. The first value is the best value obtained by that particle so far, called the personal best and the 2nd value is the best value among all the particles, called the global best. The optimizer compares these two values for each particle and stores the best value till it reaches a stopping criteria (number of iterations or pre-defined value of the fitness function). Each particle moves towards the personal best and global best based on its current position and velocity. The position of a particle can be expressed as per the below equation [70].
where, Z is the position of the particle, i is the number of the particle, V is the velocity of the particle, and T is the number of iterations. The velocity of the particle can be obtained by the following equation [71].
where, W is the constriction coefficient, p i is the personal best position, p g is the global best position, c 1 is the acceleration coefficient (individual), c 2 is the acceleration coefficient (social), and r 1 and r 2 are random numbers between 0 and 1. The constriction coefficient can be calculated as the following equation.
The typical value of k, c 1 , and c 2 is 1, 2.05, and 2.05, respectively. The implementation steps of PSO are discussed below.
Step 1: Read the following input data of the proposed energy systems: A.
Load the electrical demand data (hourly propulsive electric and auxiliary load data for one year). B.
Load the climatological data (hourly wind speed, solar irradiance, and ambient temperature) to calculate the renewable energy where applicable. C.
Load the system component's features (energy generation from MR, fossil fuel-based generator, solar PV, wind turbine, and charging and discharging from energy storage device). D.
Load economic and other parameters of the system and system components, such as project lifetime, inflation rate, discount rate, capital cost, O&M cost, replacement cost, fuel cost, refueling cost, decommissioning cost, and the salvage value.
Step 2: Initialize the parameters of DE and required system components: Step 3: Generate the initial solutions randomly within the boundary of the decision variables.
Step 4: Use the individual particle position to determine the objective function value.
Step 5: Update the best position of the individual particle by comparing it with all other populations.
Step 6: Compare the personal best with the global best and update the global best with the individual particle's personal best that has the minimum value of the objective function.
Step 7: Update the velocity of the individual particle.
Step 8: Update the position of the individual particle.
Step 9: Repeat step 3 to step 8 till all particles are evaluated in each iteration. Store the best cost value.
Step 10: Stop the simulation if it reaches the maximum number of iterations.

Results
The findings of this study are discussed in this section. Section 5.1 covers the comparison among the proposed energy systems in terms of technical, economic, and environmental KPIs. The performance of the DE algorithm is evaluated and validated in Section 5.2.

Comparison among the Proposed Hybrid Energy Systems
The Differential Evolution (DE) algorithm is used to determine the lowest Net Present Cost of the four energy systems proposed for marine transportation. The four energy systems are defined as Case-01: Standalone Fossil Fuel-based Energy Systems, Case-2: Renewable and Fossil Fuel-based Energy Systems, Case-3: Standalone Nuclear Energy Systems, and Case-04: Nuclear Renewable Hybrid Energy Systems.
Case-04 shows the lowest NPC (532.11 million $) among all the energy systems, while Case-01 has the highest NPC (877.61 million $). After introducing renewable energy with fossil fuel in Case-02, the NPC (875.90 million $) reduces around 1.7 million USD compared with Case-01. Case-03 exhibits lower NPC (538.05 million $) than Case-01 and Case-02 but a little bit higher than Case-04. The LCOE of all four cases follows a similar trend of NPC. The NPC and LCOE of all four cases are shown in Figure 2. Step 6: Compare the personal best with the global best and update the global best with the individual particle's personal best that has the minimum value of the objective function.
Step 7: Update the velocity of the individual particle.
Step 8: Update the position of the individual particle.
Step 9: Repeat step 3 to step 8 till all particles are evaluated in each iteration. Store the best cost value.
Step 10: Stop the simulation if it reaches the maximum number of iterations.

Results
The findings of this study are discussed in this section. Section 5.1 covers the comparison among the proposed energy systems in terms of technical, economic, and environmental KPIs. The performance of the DE algorithm is evaluated and validated in Section 5.2.

Comparison among the Proposed Hybrid Energy Systems
The Differential Evolution (DE) algorithm is used to determine the lowest Net Present Cost of the four energy systems proposed for marine transportation. The four energy systems are defined as Case-01: Standalone Fossil Fuel-based Energy Systems, Case-2: Renewable and Fossil Fuel-based Energy Systems, Case-3: Standalone Nuclear Energy Systems, and Case-04: Nuclear Renewable Hybrid Energy Systems.
Case-04 shows the lowest NPC (532.11 million $) among all the energy systems, while Case-01 has the highest NPC (877.61 million $). After introducing renewable energy with fossil fuel in Case-02, the NPC (875.90 million $) reduces around 1.7 million USD compared with Case-01. Case-03 exhibits lower NPC (538.05 million $) than Case-01 and Case-02 but a little bit higher than Case-04. The LCOE of all four cases follows a similar trend of NPC. The NPC and LCOE of all four cases are shown in Figure 2.  Table 1 shows the optimal rating of each component of all four energy systems. For all the energy systems, most of the energy is delivered by the diesel generator or MR based on the type of the energy system. The renewable energy penetration is higher in Case-04 compared to Case-02. The renewable energy penetration in Case-04 is about 1.61% whereas the renewable energy penetration in Case-02 is around 1.39%. Case-02 shows the highest dependency on the energy storage device, whereas Case-03 exhibits the lowest dependency on the energy storage device for the optimal energy configuration. The ratings of the energy storage device for Case-01,   Table 1 shows the optimal rating of each component of all four energy systems. For all the energy systems, most of the energy is delivered by the diesel generator or MR based on the type of the energy system. The renewable energy penetration is higher in Case-04 compared to Case-02. The renewable energy penetration in Case-04 is about 1.61% whereas the renewable energy penetration in Case-02 is around 1.39%. Case-02 shows the highest dependency on the energy storage device, whereas Case-03 exhibits the lowest dependency on the energy storage device for the optimal energy configuration. The ratings of the energy storage device for Case-01, Case-02, Case-03, and Case-04 are 2341.38 kWh, 3907.92 kWh, 637.17 kWh, and 1351.92 kWh, respectively. The Loss of Power Probability (LPSP) of all energy systems is within the defined range (less or equal to 8%) that ensures the utmost reliability while optimizing the systems. The GRF is another reliability factor that ensures the system's reliability in terms of energy generation. For a reliable energy system, usually, the GRF should be at least 100%. Among the four cases, Case-03 has the highest value of GRF (145.38%), and Case-01 has the lowest GRF (141.31%). The GRF of Case-04 is 144.82%, which is very close to Case-03. Therefore, Case-04 performs better than the other three energy systems in the ship in terms of economic (NPC, LCOE) and reliability (LPSP, GRF) KPIs.
The CO 2 emissions is an environmental KPI that determines the impact of energy systems on the environment. In terms of environmental aspects, Case-04 shows the least impact on the environment. Only 935.57 tons of CO 2 is emitted from N-R HES whereas, 144,714.33 tons of CO 2 emits each year from Case-01. Hence, the proposed N-R HES will reduce around 99% of CO 2 emissions compared to fossil fuel-based energy systems (Case-01). The penalty for emitting CO 2 is the highest for Case-01 and the lowest for Case-04. The weight of the energy systems and W ratio (the ratio between the total weight of the energy system and the total cargo-carrying capacity of the ship) are two important factors that determine the cargo-carrying performance of the ship. The energy system weight for Case-03 and Case-04 are lower compared to Case-01 and Case-02. The energy system weight of Case-01, Case-02, Case-03, and Case-04 is 1,546,864.99 kg, 1,545,210.58 kg, 445,816.43 kg, and 779,799.66 kg, respectively. Compared to Case-01, the energy system weight of Case-03 and Case-04 is around 3 and 2 times lower, respectively. By using the proposed N-R HES in the ship, the ship operator will be able to carry 767 tons more cargo, which will increase the revenue of the ship. The lower value of W ratio is expected to carry more cargo by ship. The energy system weight and W ratio of four energy systems are presented in Figure 3.  The cost distribution of different energy systems is shown in Figure 4. The capital cost of Case-01 and Case-02 is very low compared to Cae-03 and Case-04. Case-03 has the highest capital cost (313.23 million $), and Case-02 has the lowest capital cost (25.88 million $). After introducing renewable energy, the capital cost is reduced for Case-02 (25.88 million $) compared to Case-01 (26.08 million $,). Also, for Case-03 the capital cost is 313.23 million $ which is reduced to 308.56 million $ after renewable energy addition in Case-04. The operation and maintenance cost follows the same trend as the capital cost. Case-02 has the lowest O&M cost (16.31 million $), whereas Case-03 has the highest O&M cost (154.91 million $). The fuel cost shows that Case-01 and Case-02 have a higher cost associated with fuel compared to Case-03 and Case-04. Case-01 has the highest fuel cost (637.77 million $) since the price of diesel fuel (per MWh) is much higher than nuclear fuel. After incorporating renewable energy in Case-02, the fuel cost (637.31 million $) is reduced a little bit. With nuclear energy, the fuel cost is reduced to a great extent. The fuel cost of Case-03 is 32.48 million $, and the fuel cost is reached the lowest value in Case-04 (31.40 million $) by introducing renewable energy with nuclear energy. The replacement cost of Case-01 is the maximum (131.03 million $) among all the energy systems. The replacement cost reduced a bit in Case-02 (129.93 million $). Unlike other costs, Case-04 has a higher replacement cost (3.71 million $) compared to Case-03 (0.65 million $). In Case-03, only battery replacement cost contributes to the total replacement cost, but replacement cost of renewable energy along with the battery replacement cost added to the total replacement cost in Case-04.  The cost distribution of different energy systems is shown in Figure 4. The capital cost of Case-01 and Case-02 is very low compared to Cae-03 and Case-04. Case-03 has the highest capital cost (313.23 million $), and Case-02 has the lowest capital cost (25.88 million $). After introducing renewable energy, the capital cost is reduced for Case-02 (25.88 million $) compared to Case-01 (26.08 million $). Also, for Case-03 the capital cost is 313.23 million $ which is reduced to 308.56 million $ after renewable energy addition in Case-04. The operation and maintenance cost follows the same trend as the capital cost. Case-02 has the lowest O&M cost (16.31 million $), whereas Case-03 has the highest O&M cost (154.91 million $). The fuel cost shows that Case-01 and Case-02 have a higher cost associated with fuel compared to Case-03 and Case-04. Case-01 has the highest fuel cost (637.77 million $) since the price of diesel fuel (per MWh) is much higher than nuclear fuel. After incorporating renewable energy in Case-02, the fuel cost (637.31 million $) is reduced a little bit. With nuclear energy, the fuel cost is reduced to a great extent. The fuel cost of Case-03 is 32.48 million $, and the fuel cost is reached the lowest value in Case-04 (31.40 million $) by introducing renewable energy with nuclear energy. The replacement cost of Case-01 is the maximum (131.03 million $) among all the energy systems. The replacement cost reduced a bit in Case-02 (129.93 million $). Unlike other costs, Case-04 has a higher replacement cost (3.71 million $) compared to Case-03 (0.65 million $). In Case-03, only battery replacement cost contributes to the total replacement cost, but replacement cost of renewable energy along with the battery replacement cost added to the total replacement cost in Case-04.  The cost distribution of different energy systems is shown in Figure 4. The capital cost of Case-01 and Case-02 is very low compared to Cae-03 and Case-04. Case-03 has the highest capital cost (313.23 million $), and Case-02 has the lowest capital cost (25.88 million $). After introducing renewable energy, the capital cost is reduced for Case-02 (25.88 million $) compared to Case-01 (26.08 million $,). Also, for Case-03 the capital cost is 313.23 million $ which is reduced to 308.56 million $ after renewable energy addition in Case-04. The operation and maintenance cost follows the same trend as the capital cost. Case-02 has the lowest O&M cost (16.31 million $), whereas Case-03 has the highest O&M cost (154.91 million $). The fuel cost shows that Case-01 and Case-02 have a higher cost associated with fuel compared to Case-03 and Case-04. Case-01 has the highest fuel cost (637.77 million $) since the price of diesel fuel (per MWh) is much higher than nuclear fuel. After incorporating renewable energy in Case-02, the fuel cost (637.31 million $) is reduced a little bit. With nuclear energy, the fuel cost is reduced to a great extent. The fuel cost of Case-03 is 32.48 million $, and the fuel cost is reached the lowest value in Case-04 (31.40 million $) by introducing renewable energy with nuclear energy. The replacement cost of Case-01 is the maximum (131.03 million $) among all the energy systems. The replacement cost reduced a bit in Case-02 (129.93 million $). Unlike other costs, Case-04 has a higher replacement cost (3.71 million $) compared to Case-03 (0.65 million $). In Case-03, only battery replacement cost contributes to the total replacement cost, but replacement cost of renewable energy along with the battery replacement cost added to the total replacement cost in Case-04.   Figure 5 shows the convergence plot for Case-01. The static penalty is added to the NPC till the 8th iteration. From the 9th iteration, the NPC is calculated without any penalty, and the optimization moves towards the minimum NPC. Figure 6 illustrates the cost breakdown of Case-01. More than two-thirds of the total cost of Case-01 is associated with the fuel cost (72.67%), and the O&M cost has the smallest portion (1.87%). The replacement cost has the second-largest share (14.93%) of the total cost. The capital cost is accounted for 2.97% of the cost while the carbon penalty made up 7.56% of the total cost.   Figure 5 shows the convergence plot for Case-01. The static penalty is added to the NPC till the 8th iteration. From the 9th iteration, the NPC is calculated without any penalty, and the optimization moves towards the minimum NPC. Figure 6 illustrates the cost breakdown of Case-01. More than two-thirds of the total cost of Case-01 is associated with the fuel cost (72.67%), and the O&M cost has the smallest portion (1.87%). The replacement cost has the second-largest share (14.93%) of the total cost. The capital cost is accounted for 2.97% of the cost while the carbon penalty made up 7.56% of the total cost.  The convergence plot of Case-02 is presented in Figure 7. The static penalty is added up to the 8th iteration and reaches the optimized value within the 200th iteration. The cost distribution for Case-02 is similar to Case-01 as shown in Figure 8. The fuel cost makes up the largest portion (72.78%) of the total cost. The share of the replacement cost for Case-02 (14.83%) is a little bit lower compared to Case-01 (14.93%). The capital cost and the O&M cost comprise only 2.95% and 1.86%, respectively, of the total cost. The CO2 emissions penalty contributes 7.57% of the total cost of Case-02.    Figure 5 shows the convergence plot for Case-01. The static penalty is added to the NPC till the 8th iteration. From the 9th iteration, the NPC is calculated without any penalty, and the optimization moves towards the minimum NPC. Figure 6 illustrates the cost breakdown of Case-01. More than two-thirds of the total cost of Case-01 is associated with the fuel cost (72.67%), and the O&M cost has the smallest portion (1.87%). The replacement cost has the second-largest share (14.93%) of the total cost. The capital cost is accounted for 2.97% of the cost while the carbon penalty made up 7.56% of the total cost.  The convergence plot of Case-02 is presented in Figure 7. The static penalty is added up to the 8th iteration and reaches the optimized value within the 200th iteration. The cost distribution for Case-02 is similar to Case-01 as shown in Figure 8. The fuel cost makes up the largest portion (72.78%) of the total cost. The share of the replacement cost for Case-02 (14.83%) is a little bit lower compared to Case-01 (14.93%). The capital cost and the O&M cost comprise only 2.95% and 1.86%, respectively, of the total cost. The CO2 emissions penalty contributes 7.57% of the total cost of Case-02.  The convergence plot of Case-02 is presented in Figure 7. The static penalty is added up to the 8th iteration and reaches the optimized value within the 200th iteration. The cost distribution for Case-02 is similar to Case-01 as shown in Figure 8. The fuel cost makes up the largest portion (72.78%) of the total cost. The share of the replacement cost for Case-02 (14.83%) is a little bit lower compared to Case-01 (14.93%). The capital cost and the O&M cost comprise only 2.95% and 1.86%, respectively, of the total cost. The CO 2 emissions penalty contributes 7.57% of the total cost of Case-02.   Figure 9 represents the convergence plot of Case-03. The algorithm finds the NPC without any static penalty from the 1st iteration. The algorithm moves to the optimal NPC range very quickly. The cost breakdown of Case-03 is shown in Figure 10. Unlike Case-01 and Case-02, the capital cost contributes more than half of the total cost (58.22%). The O&M cost makes up the second-largest (28.79%) portion of the total cost. The fuel cost, refueling cost, and decommissioning cost comprise 6.04%, 3.73%, and 3.02% cost, respectively. The shares of the replacement cost and the CO2 emissions penalty cost are very low (less than 1%) for Case-03.   Figure 9 represents the convergence plot of Case-03. The algorithm finds the NPC without any static penalty from the 1st iteration. The algorithm moves to the optimal NPC range very quickly. The cost breakdown of Case-03 is shown in Figure 10. Unlike Case-01 and Case-02, the capital cost contributes more than half of the total cost (58.22%). The O&M cost makes up the second-largest (28.79%) portion of the total cost. The fuel cost, refueling cost, and decommissioning cost comprise 6.04%, 3.73%, and 3.02% cost, respectively. The shares of the replacement cost and the CO2 emissions penalty cost are very low (less than 1%) for Case-03.  Figure 9 represents the convergence plot of Case-03. The algorithm finds the NPC without any static penalty from the 1st iteration. The algorithm moves to the optimal NPC range very quickly. The cost breakdown of Case-03 is shown in Figure 10. Unlike Case-01 and Case-02, the capital cost contributes more than half of the total cost (58.22%). The O&M cost makes up the second-largest (28.79%) portion of the total cost. The fuel cost, refueling cost, and decommissioning cost comprise 6.04%, 3.73%, and 3.02% cost, respectively. The shares of the replacement cost and the CO 2 emissions penalty cost are very low (less than 1%) for Case-03.  The convergence plot of Case-04 is presented in Figure 11. The cost distribution of Case-04 is illustrated in Figure 12. Like Case-03, the capital cost is accounted for more than half of the total cost (57.97%). The O&M cost comprises 28.63% of the cost. The fuel cost, decommissioning cost, and refueling cost make up 5.90%, 2.95%, and 3.77%, respectively, of the total cost. The replacement cost and the cost associated with CO2 emissions comprise only 0.70% and 0.08% cost, sequentially. 58   The convergence plot of Case-04 is presented in Figure 11. The cost distribution of Case-04 is illustrated in Figure 12. Like Case-03, the capital cost is accounted for more than half of the total cost (57.97%). The O&M cost comprises 28.63% of the cost. The fuel cost, decommissioning cost, and refueling cost make up 5.90%, 2.95%, and 3.77%, respectively, of the total cost. The replacement cost and the cost associated with CO2 emissions comprise only 0.70% and 0.08% cost, sequentially. 58  The convergence plot of Case-04 is presented in Figure 11. The cost distribution of Case-04 is illustrated in Figure 12. Like Case-03, the capital cost is accounted for more than half of the total cost (57.97%). The O&M cost comprises 28.63% of the cost. The fuel cost, decommissioning cost, and refueling cost make up 5.90%, 2.95%, and 3.77%, respectively, of the total cost. The replacement cost and the cost associated with CO 2 emissions comprise only 0.70% and 0.08% cost, sequentially.

Comparison and Validation of Optimization Techniques
From the findings of Section 5.1, it is evident that Case-04 is best among all the four energy systems for marine transportation in terms of technical, economic, and environmental KPIs. The Differential Evolution (DE) optimization technique is used to optimize the systems. Hence, it is important to compare the findings of DE optimization with other widely used optimization techniques to validate and reinforce the results. In this section, the NPC of Nuclear-Renewable Hybrid Energy Systems (N-R HESs) obtained by DE are compared with PSO and ADE algorithm. Figure 13 shows the NPC of Case-04 for DE, ADE, and PSO. The lowest NPC ($532,112,182.13) is obtained by using the DE optimization technique although the NPC (532,114,720.56) determined by the ADE is very close to DE. The difference between DE and ADE is only 0.0005%. This suggests that the impact of the control parameters of DE optimization is very minimal on finding the lowest value 57

Comparison and Validation of Optimization Techniques
From the findings of Section 5.1, it is evident that Case-04 is best among all the four energy systems for marine transportation in terms of technical, economic, and environmental KPIs. The Differential Evolution (DE) optimization technique is used to optimize the systems. Hence, it is important to compare the findings of DE optimization with other widely used optimization techniques to validate and reinforce the results. In this section, the NPC of Nuclear-Renewable Hybrid Energy Systems (N-R HESs) obtained by DE are compared with PSO and ADE algorithm. Figure 13 shows the NPC of Case-04 for DE, ADE, and PSO. The lowest NPC ($532,112,182.13) is obtained by using the DE optimization technique although the NPC (532,114,720.56) determined by the ADE is very close to DE. The difference between DE and ADE is only 0.0005%. This suggests that the impact of the control parameters of DE optimization is very minimal on finding the lowest value 57

Comparison and Validation of Optimization Techniques
From the findings of Section 5.1, it is evident that Case-04 is best among all the four energy systems for marine transportation in terms of technical, economic, and environmental KPIs. The Differential Evolution (DE) optimization technique is used to optimize the systems. Hence, it is important to compare the findings of DE optimization with other widely used optimization techniques to validate and reinforce the results. In this section, the NPC of Nuclear-Renewable Hybrid Energy Systems (N-R HESs) obtained by DE are compared with PSO and ADE algorithm. Figure 13 shows the NPC of Case-04 for DE, ADE, and PSO. The lowest NPC ($532,112,182.13) is obtained by using the DE optimization technique although the NPC (532,114,720.56) determined by the ADE is very close to DE. The difference between DE and ADE is only 0.0005%. This suggests that the impact of the control parameters of DE optimization is very minimal on finding the lowest value

Sensitivity Analysis
The N-R HES for marine transportation employs several economic, technical, and environmental variables, such as water temperature, water density, solar irradiance, wind speed, air temperature, discount rate, inflation rate, project lifetime, fuel cost, replacement cost, capital cost, operation and maintenance cost, CO2 penalty, Wratio, and LPSP. Hence, a sensitivity analysis is crucial to assess the impact of these variables on the findings of the base model. The sensitivity analysis reinforces the results. The sensitivity analysis is also important to determine how sensitive the system to different parameters. The following subsections will assess the impact of different economic, technical, and environmental variables on the NPC of the four energy systems. The energy systems are optimized separately to get the optimized value of NPC for each parameter change in the sensitivity analysis.

Sensitivity Assessment of Discount Rate on NPC
In the base case, the discount rate is considered as 8%. In the sensitivity analysis, the discount rate is varied from 5% to 10% to determine the impact on NPC. Figure 14 shows the NPC of the four energy systems for different discount rates varied from 5% to 10%. From the figure, it can be observed that Case-04 always demonstrates the lowest NPC for the marine vessels irrespective of the discount rate. Also, the NPC of all the energy systems decreases with the increment of the discount rate.

Sensitivity Assessment of Inflation Rate on NPC
The inflation rate was considered 2% in the base case of the analysis. In the sensitivity analysis, the inflation rate is varied from 1% to 5% to identify the impact of the inflation rate on NPC of different energy systems. Figure 15 shows the NPC of different energy systems for different inflation rates. The NPC increases with the increase of the inflation rate with a varying range. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the inflation rate.

Sensitivity Assessment of Project Lifetime on NPC
The project life is considered as 40 years in the base case analysis. The project lifetime can be different based on the requirement of the ship operator and lifetime of the ship.

Sensitivity Analysis
The N-R HES for marine transportation employs several economic, technical, and environmental variables, such as water temperature, water density, solar irradiance, wind speed, air temperature, discount rate, inflation rate, project lifetime, fuel cost, replacement cost, capital cost, operation and maintenance cost, CO 2 penalty, W ratio , and LPSP. Hence, a sensitivity analysis is crucial to assess the impact of these variables on the findings of the base model. The sensitivity analysis reinforces the results. The sensitivity analysis is also important to determine how sensitive the system to different parameters. The following subsections will assess the impact of different economic, technical, and environmental variables on the NPC of the four energy systems. The energy systems are optimized separately to get the optimized value of NPC for each parameter change in the sensitivity analysis.

Sensitivity Assessment of Discount Rate on NPC
In the base case, the discount rate is considered as 8%. In the sensitivity analysis, the discount rate is varied from 5% to 10% to determine the impact on NPC. Figure 14 shows the NPC of the four energy systems for different discount rates varied from 5% to 10%. From the figure, it can be observed that Case-04 always demonstrates the lowest NPC for the marine vessels irrespective of the discount rate. Also, the NPC of all the energy systems decreases with the increment of the discount rate.

Sensitivity Assessment of Inflation Rate on NPC
The inflation rate was considered 2% in the base case of the analysis. In the sensitivity analysis, the inflation rate is varied from 1% to 5% to identify the impact of the inflation rate on NPC of different energy systems. Figure 15 shows the NPC of different energy systems for different inflation rates. The NPC increases with the increase of the inflation rate with a varying range. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the inflation rate.

Sensitivity Assessment of Project Lifetime on NPC
The project life is considered as 40 years in the base case analysis. The project lifetime can be different based on the requirement of the ship operator and lifetime of the ship. Hence, a sensitivity analysis is carried out for different project lifetimes. The project lifetime is varied from 20 years to 60 years with 10 years interval. The NPC of different energy systems for different project lifetimes is presented in Figure 16. The figure illustrates that with the increase of project lifetime NPC increases for all the cases. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the project lifetime.

Sensitivity Assessment of Electrical Power Requirement on NPC
This study employs several parameters, like seawater temperature, seawater density, length and width of the ship, average speed of the ship, and draught of the ship, to estimate the ship power. These parameters change with the position of the ship, weather, and type of the ship. Any change in these parameters will result in different electrical power requirements. Hence, it is important to assess the performance of the energy systems for different electrical power requirements. In this section, the electrical power of the ship is varied from −20% to +20% with 10% increment. Figure 17 shows the NPC of different energy systems with different electrical power. The figure tells that the NPC increases accordingly with the increase of electrical power. However, Case-04 illustrates the lowest NPC irrespective of the electrical power.
Hence, a sensitivity analysis is carried out for different project lifetimes. The project lifetime is varied from 20 years to 60 years with 10 years interval. The NPC of different energy systems for different project lifetimes is presented in Figure 16. The figure illustrates that with the increase of project lifetime NPC increases for all the cases. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the project lifetime.

Sensitivity Assessment of Electrical Power Requirement on NPC
This study employs several parameters, like seawater temperature, seawater density, length and width of the ship, average speed of the ship, and draught of the ship, to estimate the ship power. These parameters change with the position of the ship, weather, and type of the ship. Any change in these parameters will result in different electrical power requirements. Hence, it is important to assess the performance of the energy systems for different electrical power requirements. In this section, the electrical power of the ship is varied from −20% to +20% with 10% increment. Figure 17 shows the NPC of different energy systems with different electrical power. The figure tells that the NPC increases accordingly with the increase of electrical power. However, Case-04 illustrates the lowest NPC irrespective of the electrical power.    Hence, a sensitivity analysis is carried out for different project lifetimes. The project lifetime is varied from 20 years to 60 years with 10 years interval. The NPC of different energy systems for different project lifetimes is presented in Figure 16. The figure illustrates that with the increase of project lifetime NPC increases for all the cases. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the project lifetime.

Sensitivity Assessment of Electrical Power Requirement on NPC
This study employs several parameters, like seawater temperature, seawater density, length and width of the ship, average speed of the ship, and draught of the ship, to estimate the ship power. These parameters change with the position of the ship, weather, and type of the ship. Any change in these parameters will result in different electrical power requirements. Hence, it is important to assess the performance of the energy systems for different electrical power requirements. In this section, the electrical power of the ship is varied from −20% to +20% with 10% increment. Figure 17 shows the NPC of different energy systems with different electrical power. The figure tells that the NPC increases accordingly with the increase of electrical power. However, Case-04 illustrates the lowest NPC irrespective of the electrical power.

Sensitivity Assessment of Solar Irradiance on NPC
In this study, the voyage route of the 'Baltic Sunrise' is considered from Iraq to Singapore. Renewable energy sources depend on the position of the ship, time, and direction. The solar irradiance will not be the same all the time, and the ship will not use the same route each time. Hence, a sensitivity assessment of solar irradiance is necessary to determine the impact on NPC. Since Case-01 and Case-03 do not have any solar PV, any change in the solar irradiance will not affect the NPC of these two cases. Hence, only Case-02 and Case-04 are considered for the assessment of solar irradiance on NPC. In this study, solar irradiance is varied from −20% to +20% with a 10% increment. Figure 18 shows the NPC of Case-02 and Case-04 for different solar irradiance. The NPC of Case-02 decreases slightly with the increase of solar irradiance while no impact on the NPC of Case-04 is observed with the change of solar irradiance. It should be noted that the penetration of

Sensitivity Assessment of Solar Irradiance on NPC
In this study, the voyage route of the 'Baltic Sunrise' is considered from Iraq to Singapore. Renewable energy sources depend on the position of the ship, time, and direction. The solar irradiance will not be the same all the time, and the ship will not use the same route each time. Hence, a sensitivity assessment of solar irradiance is necessary to determine the impact on NPC. Since Case-01 and Case-03 do not have any solar PV, any change in the solar irradiance will not affect the NPC of these two cases. Hence, only Case-02 and Case-04 are considered for the assessment of solar irradiance on NPC. In this study, solar irradiance is varied from −20% to +20% with a 10% increment. Figure 18 shows the NPC of Case-02 and Case-04 for different solar irradiance. The NPC of Case-02 decreases slightly with the increase of solar irradiance while no impact on the NPC of Case-04 is observed with the change of solar irradiance. It should be noted that the penetration of

Sensitivity Assessment of Solar Irradiance on NPC
In this study, the voyage route of the 'Baltic Sunrise' is considered from Iraq to Singapore. Renewable energy sources depend on the position of the ship, time, and direction. The solar irradiance will not be the same all the time, and the ship will not use the same route each time. Hence, a sensitivity assessment of solar irradiance is necessary to determine the impact on NPC. Since Case-01 and Case-03 do not have any solar PV, any change in the solar irradiance will not affect the NPC of these two cases. Hence, only Case-02 and Case-04 are considered for the assessment of solar irradiance on NPC. In this study, solar irradiance is varied from −20% to +20% with a 10% increment. Figure 18 shows the NPC of Case-02 and Case-04 for different solar irradiance. The NPC of Case-02 decreases slightly with the increase of solar irradiance while no impact on the NPC of Case-04 is observed with the change of solar irradiance. It should be noted that the penetration of solar PV in Case-04 is very minimal. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the change in solar irradiance.

Sensitivity Assessment of Wind Speed on NPC
Like solar irradiance, wind speed also varies with the time and the route of the ship. Hence, a sensitivity assessment is carried out to identify the impact of wind speed on NPC on different energy systems. As there is no wind turbine for Casse-01 and Case-03, there is no impact with the change of wind speed on these cases. In this assessment, wind speed is varied from −20% to +20% with a 10% increment. Figure 19 shows the NPC of Case-02 and Case-04 for different wind speeds. The NPC of Case-04 decreases with the increase of wind speed while no impact on the NPC of Case-02 is observed with the change of wind speed. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the change in solar irradiance.

Sensitivity Assessment of Wind Speed on NPC
Like solar irradiance, wind speed also varies with the time and the route of the ship. Hence, a sensitivity assessment is carried out to identify the impact of wind speed on NPC on different energy systems. As there is no wind turbine for Casse-01 and Case-03, there is no impact with the change of wind speed on these cases. In this assessment, wind speed is varied from −20% to +20% with a 10% increment. Figure 19 shows the NPC of Case-02 and Case-04 for different wind speeds. The NPC of Case-04 decreases with the increase of wind speed while no impact on the NPC of Case-02 is observed with the change of wind speed. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the change in solar irradiance.

Conclusions, Contribution and Future Scope of Work
Several studies have been conducted on stand-alone nuclear energy systems and fossil fuel-renewable hybrid energy systems for ocean-going marine ships to reduce GHG

Sensitivity Assessment of Wind Speed on NPC
Like solar irradiance, wind speed also varies with the time and the route of the ship. Hence, a sensitivity assessment is carried out to identify the impact of wind speed on NPC on different energy systems. As there is no wind turbine for Casse-01 and Case-03, there is no impact with the change of wind speed on these cases. In this assessment, wind speed is varied from −20% to +20% with a 10% increment. Figure 19 shows the NPC of Case-02 and Case-04 for different wind speeds. The NPC of Case-04 decreases with the increase of wind speed while no impact on the NPC of Case-02 is observed with the change of wind speed. However, Case-04 shows the lowest NPC among all the energy systems irrespective of the change in solar irradiance.

Conclusions, Contribution and Future Scope of Work
Several studies have been conducted on stand-alone nuclear energy systems and fossil fuel-renewable hybrid energy systems for ocean-going marine ships to reduce GHG

Conclusions, Contribution and Future Scope of Work
Several studies have been conducted on stand-alone nuclear energy systems and fossil fuel-renewable hybrid energy systems for ocean-going marine ships to reduce GHG emissions. However, both systems have some limitations in terms of economic, regulatory, and technical points of view. Hence, a nuclear-renewable hybrid energy system on marine ships has been proposed in this paper to overcome these limitations. The energy demand of a ship is estimated by using actual data to simulate the real-world scenario. The mathemat-ical modeling and simulation are carried out in MATLAB environment by implementing a Differential Evolution (DE) algorithm. The findings of this paper suggest that an integrated nuclear-renewable hybrid energy system could be the best solution to make ocean-going ships free from emissions, and it is technically and economically feasible. Among the four energy systems, N-R HES shows the lowest NPC (532.11 million $) compared to Standalone Fossil Fuel-based Energy System (877.61 million $), Renewable Energy and Fossil Fuel-based Hybrid Energy System (875.90 million $), and Stand-alone Nuclear Energy System (538.05 million $). The proposed N-R HES emits the lowest CO 2 gas compared to other energy systems in the entire project life. The study identifies the 'Capital Cost' and 'O&M Cost' as the major contributors of NPC for Case-03 and Case-04, whereas 'Fuel Cost' and 'Replacement Cost' have the major contribution towards NPC for Case-01 and Case-02. This study also ensures the reliability of the proposed with adequate margins during the optimization of the systems. The research determines MRs as a competent candidate to replace the diesel generator from maritime transportation.
The optimized result of Case-04 obtained by the Differential Evolution (DE) algorithm is compared with PSO and ADE to validate the results. Among all the optimization techniques, DE performs better than the other two algorithms. The NPC (N-R HES) obtained by DE algorithm ($532,112,182) is the lowest compared to ADE ($532,114,720), and PSO ($532,118,322). In the latter section, sensitivity analysis is carried out to reinforce the findings. The sensitivity assessment is conducted by varying the discount rate, inflation rate, project lifetime, electrical power requirement, solar irradiance, and wind speed. The sensitivity analysis shows NPC increases with the increase of inflation rate, project lifetime, and electrical power demand for energy systems. The NPC decreases with the increase of discount rate. The sensitivity analysis concludes Case-04 always shows the lowest NPC among all the energy systems irrespective of any change in the parameter.
The main contribution of this study is summarized below A comprehensive literature review on the impact of maritime transportation on the environment, renewable and fossil fuel-based hybrid energy systems, nuclear propulsion, development of N-R HES in land-based applications to identify the problem and gaps. Estimate the energy requirement of an oil tanker with the collaboration of industry (FleetMon) to address a real-world problem in this study. Four energy systems are introduced to determine the most feasible energy system for maritime transportations. The feasibility is assessed based on technical, economical, and environmental KPIs. To ensure the reliability of the energy system, certain constraints are added.
The economical model of all the system components is developed mathematically. This mathematical model is then optimized in the popular and versatile simulator MATLAB. The Differential Evolution (DE) algorithm is used to optimize the energy systems. The performance of the DE algorithm is compared to another optimization technique, PSO. Also, the impact of the control parameters of the DE algorithm is measured by the Adaptive Differential Algorithm (ADE).
To address the variability of the parameters like discount rate, inflation rate, project lifetime, electrical power requirements, solar irradiance, and wind speed, a sensitivity analysis is conducted to reinforce the findings.
The future scope of work is summarized below: In this study, average solar irradiance, temperature, and wind speed are used while optimizing the systems. In the future, the actual solar irradiance, temperature, and wind speed can be used to make the study more realistic.
In this study, a specific ship route is considered. However, analysis can be extended to different major international shipping routes to assess the most feasible route for N-R HES.