Abstract
Military technological advancements are crucial for meeting evolving operational needs. Submarines are essential naval assets due to their difficulty in detection. Selecting appropriate equipment for mission efficiency is essential, and automating these tasks with a Decision Support System (DSS) is vital. Air Independent propulsion (AIP) is critical for covert submarine operations, allowing propulsion underwater without surfacing for oxygen. Choosing the right AIP system is challenging because it requires balancing operational efficiency, economic cost-effectiveness, and environmental consequences, all of which entail complex factors that are difficult to quantify or compare. The developed DSS is based on a two-stage Fuzzy Mamdani Model (FMM), applied to perform inference over a submarine AIP. Three reference values were obtained in the first stage to aid operational, economic, and environmental decision-making. In the second stage, these three values create a recommendation value using a percentage system selection. To prove its effectiveness, the DSS was used to evaluate four types of AIPs to establish the level of recommendation and the values of the subcategories in the initial phase. Experts from the Portuguese Navy (PoN) assisted with the system development by defining decision-making rules. The FMM offers advantages over traditional decision-making models by effectively managing vague and qualitative data, enabling reasoning that reflects expert judgment instead of relying on deterministic models. This DSS not only enhances the accuracy and transparency of AIP selection but also provides a scalable and adaptable framework applicable to various complex decision-making scenarios in other engineering applications and fields.
Keywords
Introduction
The Armed Forces (AF) play a key role in technological innovation and scientific development (Chin, 2019). Technological superiority is vital for battlefield advantage, requiring careful equipment selection due to economic and technological constraints. War impacts territories, populations, and the economy, as significant resources divert from essential areas of a country (Liadze et al., 2023).
The underwater realm has perpetually captivated humanity, with the development of improved exploration systems being a longstanding objective (Marx, 1990). Two primary concepts for underwater navigation have emerged (Moorhouse, 2015): submersibles and submarines. Submersibles require a principal vessel for launch and recovery due to their limited power capabilities, whereas submarines can operate independently. The first successful deployment of submersibles in warfare occurred during World War II, conferring a substantial advantage to their users (Jones, 2022). The concept of submarines was first conceived by William Bourne in 1578, although practical testing was not conducted then (Avallone, 2000). The inaugural functional submarine, the Turtle, was constructed in 1776, paving the way for numerous technological advancements to refine this crucial instrument in military operations (Lautenschläger, 1987).
Submarines are efficient platforms with superior stealth, giving them a strategic advantage over other naval vessels (Haffa & Patton, 1991). Conventional submarines need to maximize underwater time for stealth. For instance, diesel-electric submarines must frequently surface to snorkel for fresh air to operate their Diesel Engines (DE) and recharge batteries, compromising their discretion (Brighton et al., 1994). This challenge led to the development of nuclear-powered submarines (Chopra, 2023), despite the high costs and nuclear waste disposal issues (Friedman-Jimenez et al., 2022).
Recent submarine advancements aim to reduce energy use, improve battery storage, and integrate Air Independent Propulsion (AIP) systems (Başhan, 2022). AIP extends a submarine’s submerged operation time and is a cost-effective alternative to diesel generators, which need air to charge batteries. It allows for a longer submerged duration while maintaining stealth (Menon et al., 2020). Standard AIP systems include Fuel Cells (FC), Stirling Engines (SE), closed-cycle Steam Turbines (ST), and closed-cycle DE (Kordesch & Simader, 1996; Menon et al., 2020; Wang et al., 2023).
The increasing complexity of real-world systems and the many variables involved in analyzing them have led to developing methods and frameworks known as Decision Support Systems (DSSs). These DSSs are designed to explore a wide range of variables and, through logical reasoning, offer insights that can assist in a particular decision-making process.
Aristotle founded classical logic, which derives binary true or false conclusions (Gabbay & Woods, 2004). Over time, probability theory and fuzzy sets introduced by (Zadeh, 1965) offered a classification method aligned with human thinking, addressing uncertainties. Fuzzy logic effectively manages imprecise qualitative data and improves decision-making where traditional logic is inadequate (Zadeh, 1988). Initially used in control systems (Stelzer et al., 2007), fuzzy logic applications have expanded across other domains.
The Mamdani model is widely used for deriving conclusions from logical rules (Dadios, 2012). Ebrahim Mamdani initially developed it to control a steam engine and boiler using linguistic control rules from experienced human operators (Mamdani, 1976). The model has since been applied in various areas, such as medical diagnosis systems (Gupta et al., 2022), management (Pourjavad & Shahin, 2018), among many others.
This work develops a DSS to select an AIP system for submarines, enhancing submerged autonomy. A two-stage Fuzzy Mamdani Model (FMM) assigns a recommendation value to each AIP system based on operational, economic, and ambient factors. It aims to choose the most suitable AIP system, reducing human subjectivity and improving decision-making. Experts from the Portuguese Navy (PoN) assisted in defining the decision-making rules to enhance accuracy. The implemented system can also be easily adapted to analyze other systems and cases.
Selecting an appropriate AIP system involves navigating complex trade-offs between qualitative and quantitative factors, often in environments of uncertainty and limited data, which makes conventional DSS approaches less effective. Fuzzy logic is particularly well-suited for this context because it facilitates decision-making based on imprecise, linguistic, and expert-derived information, commonly encountered in defense applications. This is especially true when evaluating submarine propulsion systems, where data may be scarce, uncertain, or subjective. DSS methodologies like Neural Networks (NN), Analytic Hierarchy Processes (AHP), and Bayesian models can also aid in selecting AIP systems. NNs excel in recognizing complex patterns and learning, but usually require extensive, high-quality datasets, which are lacking in this defense decision context (Shafi et al., 2021). AHP provides a structured multi-criteria decision-making framework through pairwise comparisons but struggles with uncertainty and the imprecise data typical in submarine propulsion evaluations (Saaty, 2008). Bayesian models offer a probabilistic approach by integrating prior knowledge with observed data, but require well-defined priors, which may introduce biases and limit adaptability in changing conditions (Koller & Friedman, 2009). Consequently, a fuzzy logic-based DSS was selected for its ability to manage uncertainty in expert assessments while ensuring flexible, interpretable decision-making (Zadeh, 1996). Fuzzy inference enables the use of qualitative inputs, adapting to complex AIP trade-offs and enhancing decision reliability while addressing imprecision in AIP propulsion evaluation.
In the submarine domain, besides some fuzzy logic approaches used to enhance depth control (Ranga et al., 2017), and the design of a fuzzy logic system for optimizing decision-making during submarine attack scenarios (Ranga et al., 2017), a fuzzy DSS for managing repair priorities of shipboard weapon engineering equipment during combat situations (Simões-Marques et al., 2000), or a fuzzy expert system developed to assist submarine commanders in damage control situations (Shi et al., 2009), no submarine or ship DSS has been established to support the selection of AIP systems or related equipment. Selecting the most suitable AIP system is challenging due to the need to balance operational effectiveness, economic cost-effectiveness, and environmental impact, which involve complex and multidimensional factors that are difficult to quantify or compare. The proposed DSS utilizes fuzzy logic to enhance decision-making in the selection of submarine propulsion systems, addressing trade-offs and uncertainties in AIP evaluation. This advancement offers a structured, flexible framework that addresses a crucial gap in decision-support methodologies, especially in the naval sector and, more specifically, submarines.
The article is structured as follows: After introducing the problem and describing the contributions and objectives in this section, Section 2 defines the problem and the adopted methodologies. Section 3 presents and describes the developed DSS. Section 4 describes the developed DSS application and discusses the obtained results. Finally, in Section 5, conclusions are provided, and further research work and development topics are suggested.
Problem Formulation & Methodologies
Defining the problem helps understand it and outline methodologies for applying the framework to other case studies. Section 2.1 illustrates the submarine technology and the importance of the AIP system to its operations, while Section 2.2 describes the method used to develop the DSS based on fuzzy reasoning.
Submarine Technology
Global navies’ use of submarines provides a significant strategic advantage during operations. The design and technical requirements of these submarines should align with their respective nations’ operational and technological needs, influenced by strategic and military demands (Hartley, 2024). Militarily, submarine operations depend heavily on their ability to operate stealthily and independently (Howard, 2011).
Compared to nuclear submarines, the main disadvantage of conventional diesel-electric submarines is their limited submerged endurance (Rains & Mitchell, 2000). They must surface or perform snorkeling to intake air for diesel generators, which makes them vulnerable to detection due to exhaust gases and noise. Advancements in anaerobic energy systems, like AIP, are essential for improving underwater endurance (Hawley et al., 1994).
The endurance limitation is a frequent challenge for traditional submarine operations, particularly for navies lacking nuclear-powered submarines. AIP technologies effectively address this limitation, enabling prolonged submerged missions without surfacing, thereby maintaining stealth and tactical advantages. For many countries, AIP offers a budget-friendly option to enhance underwater operational capabilities without the complexities and high costs associated with nuclear propulsion.
As previously described, AIP systems for submarines typically consist of FC, SE, and ST (Ahmadi et al., 2017; Kerros et al., 1994; Kordesch & Simader, 1996; Menon et al., 2020; Wang et al., 2023). Each of these power sources has its advantages and disadvantages. FC are known for their high efficiency and low acoustic and thermal signatures, providing a quiet power source with water and heat as byproducts. However, they require pure hydrogen, which leads to complex storage and high costs. SE use standard diesel fuel and have a low acoustic signature but face challenges related to size and the need for significant maintenance. Closed-cycle ST deliver high power output, enabling submarines to achieve higher speeds, but they are complex, inefficient, and have higher operational costs. Each AIP method offers unique capabilities, making the propulsion system vital in submarine design and mission success.
Choosing an appropriate AIP system for submarines is crucial for any Navy due to its strategic implications and various factors influenced by technology. National needs and capabilities are essential for successful selection. Selection criteria vary based on the Navy type, mission profile, or analyst, necessitating a careful assessment. This complex process can result in time-consuming and subjective decisions. A clear illustration of the inefficiencies arising from poor AIP system selection is the Spanish S-80 Plus submarine program, which faced substantial delays and redesigns due to difficulties in integrating its complex AIP bioethanol-reforming fuel cell system (Navantia, 2020; Revista Ejércitos, 2024). As described in the next section, fuzzy inference systems can simplify these challenges, reducing decision time and imprecision while aligning results with analyst expertise.
Fuzzy Reasoning
Fuzzy reasoning, typically consists of three stages (Mendel, 1995): Fuzzification (Section 2.2.1), inference mechanism (Section 2.2.2), and defuzzification (Section 2.2.3). First, input data becomes fuzzy sets through fuzzification, changing numbers to linguistic variables using membership functions. Each variable receives membership degrees from 0 to 1, reflecting truth in fuzzy logic. Next, the inference mechanism links input data to knowledge-based rules and merges multiple logical results. Finally, fuzzy values revert to real values based on earlier inferences in defuzzification.
Fuzzification
In the fuzzification stage, we will use a variant of Gaussian membership functions, specifically a membership function based on a mixture of two Gaussian functions, which can be defined as (Zhao et al., 2019):
The Mamdani model (Dadios, 2012; Mamdani, 1976) utilizes max-min composition, where the minimum operation evaluates the degree of truth for antecedents within each rule as an intersection, and the maximum operation combines the consequents from different rules through a union to form a cohesive output (Sakti, 2014). After fuzzification, the inference process applies logical operators to combine input fuzzy sets, defining new subsets from the initial sets and firing any rule where this combination yields a non-zero membership degree. The max-min composition corresponds to the AND and OR logical operators used in this process. The implication method then models the fuzzy output set (consequent) for each rule, with the resulting membership degree truncating the consequent fuzzy set by applying the min function to the maximum defined for that set and the value from the antecedent using the logical AND operator. If a rule is not fired, the consequence is a constant membership degree of zero. All rule consequents aggregate in the final inference step, generating the fuzzy output set as the union of the consequent fuzzy subsets with the logical OR operator.
Defuzzification
Different methods can be utilized for defuzzification, with the centroid method, also referred to as the Center of Gravity (COG) method, being the most common (Takagi & Sugeno, 1985). The centroid method calculates the center of mass of the area under the membership function curve, where the desired value is the coordinate of this center, defined as (Takagi & Sugeno, 1985):
Consider a simple triangular membership function, as illustrated in Figure 1, defined by:

Example of a Triangular Membership Function
Applying the centroid method, we obtain:
As expected, the centroid of this triangular membership function coincides precisely with its midpoint,
The centroid method computes the defuzzified output by identifying the geometric center of the aggregated fuzzy set, resulting in balanced and stable outcomes that demonstrate robustness (Takagi & Sugeno, 1985). In contrast, the bisector method divides the area of the fuzzy set into two halves but can be sensitive to asymmetric membership functions (Ross, 2005). Methods such as the Smallest Of Maximum (SOM) and Largest Of Maximum (LOM) select extreme points linked to maximum membership values, potentially causing abrupt shifts in the output. Conversely, the Mean Of Maximum (MOM) averages all points of maximum membership, which can sometimes produce unintuitive results (Mendel, 2017). Overall, centroid defuzzification typically provides smoother and more reliable outputs, making it particularly beneficial in control and decision-making applications (Mendel, 2017; Ross, 2005).
A two-stage FMM, using fuzzy reasoning, serves as a DSS. First, we evaluate the operational, economic, and environmental criteria with fuzzy rules to generate subcategory scores. Then, these scores are aggregated using higher-level fuzzy rules to provide an overall recommendation. This DSS assigns values to the analyzed AIP systems. The choice of economic, operational, and environmental criteria as the primary focus stems from the necessity of a thorough and well-rounded assessment of AIP systems. Economic criteria encompass the acquisition, maintenance, and operational costs, which are crucial for long-term adoption and budgetary considerations. Operational criteria measure efficiency, power, and acoustic signature, significantly affecting the systems’ tactical and strategic importance. Environmental criteria ensure the consideration of thermal emissions, byproducts, and their respective size. Collectively, these aspects provide a comprehensive perspective that enables informed, responsible, and mission-oriented decisions within defense frameworks. The following sections detail the methodology and analysis.
To develop the DSS, a theoretical study of the various systems that typically constitute an AIP was first conducted, as described in Section 2.1, to identify the key variables that support the decision-making process, as illustrated in Figure 2. After defining and organizing these variables into subcategories, a recommendation level is determined based on rule evaluation to identify the system that best complies.

Defining Key Variables, Organizing them into Subcategories, and Determining the Recommendation Process Level.
After defining the variables, subcategories, and recommendation level, fuzzy sets were created by converting numerical data into linguistic variables and establishing the rules. This process involved analyzing related work and leveraging PoN expertise in submarine systems.
The DSS will explore the relationships between characteristics and subcategories, as well as between subcategories and recommendation levels, to provide a comprehensive profile of each system across three distinct fuzzy subcategories: operational, economic, and environmental. The definitions of variables and subcategories will be outlined in Section 3.1, and the recommendation level will be discussed in Section 3.2.
The systems studied can be installed in modern submarines, with any challenges, like structural stability, expected to be managed by naval engineering and therefore not addressed here. Key characteristics enable a clear comparison among the systems. While safety, autonomy, and consumption might seem non-essential, they are crucial in modern designs. Autonomy largely depends on the system and platform, while consumption reflects operational costs. The selected variables are grouped into three subcategories: operational, economic, and environmental, providing a comprehensive evaluation of technical characteristics, costs, and potential environmental impacts. These include variables like efficiency, power, acoustic signature, acquisition costs, maintenance costs, operating costs, thermal emissions, expelled byproducts, and size, as shown in Figure 3.

The Subcategories and Respective Variables Chosen for Evaluation.
Due to the classified nature of the military systems discussed, access to data is limited. As a result, unavailable data is estimated using qualitative methods, drawing on theoretical studies and insights from PoN submarine experts to establish logical rules. The system’s logical rules were developed based on responses from a questionnaire administered to 21 PoN militaries with submarine specialization. The system is designed to adapt to real-world situations, with the Navy working alongside manufacturers and suppliers to gather precise data. Once obtained, this data modifies the existing system or standardizes new variables based on defined fuzzy set values.
The operational subcategory objectively enhances submarine performance while meeting essential operational requirements. It comprises three variables: efficiency, power, and acoustic signature. Efficiency is the energy or heat input ratio to the work output produced by a system. Power refers to the rate at which energy is converted into work. The acoustic signature indicates the equipment’s noise level at full capacity.
A system with high efficiency offers significant operational advantages, allowing energy supplies to last much longer. In contrast, a system with low efficiency consumes energy faster, resulting in a shorter duration of energy supply. For this variable, four fuzzy sets have been defined based on concrete values ranging from 0 to 100%: {Very Poor (VP), Poor (P), Good (G), Very Good (VG)}. The membership functions for these fuzzy sets are determined using a mixture of two Gaussian functions, as described by Equation 1 and illustrated in Figure 4.

Operation Membership Functions - Efficiency.
Table 1 presents the defined fuzzy sets for efficiency, along with their corresponding mean values
Parameters for Each Efficiency Membership Function.
A system with high power output is advantageous for operational purposes as it enables rapid battery charging. Conversely, a system with a relatively low power output significantly reduces its capacity for fast charging. Four fuzzy sets were defined for the power variable: {VP, P, G, VG}. The membership functions for these sets were established using a mixture of two Gaussian functions, as described in Equation 1 and illustrated in Figure 5. Table 2 presents the defined fuzzy sets and their corresponding parameters

Operation Membership Functions - Power.
Parameters for Each Power Membership Function.
For the acoustic signature variable, four fuzzy sets have been defined: {Very Low (VL), Low (L), High (H), Very High (VH)}. The membership functions for these sets are identical to those for the power variable, as outlined in Equation 1 and illustrated in Figure 6. Table 3 presents the defined fuzzy sets along with their respective parameters
By defining the three previously mentioned variables, efficiency, power, and acoustic signature, we established the first subcategory, the operational. For this subcategory, three fuzzy sets were defined: {Low (L), Medium (M), High (H)}. The membership functions for these sets were created using a combination of two Gaussian functions, as outlined in Equation 1 and illustrated in Figure 7. Table 4 lists the specific parameters

Operation Membership Functions - Acoustic Signature.
As previously stated, the inference method will use the FMM, as outlined in Section 2. With the operational subcategory and its defining variables fully established, creating the rule base is necessary to obtain the required inferences.
The economic subcategory refers to the costs associated with a system, encompassing several essential variables: acquisition, maintenance, and operational costs. The acquisition cost includes all expenses for obtaining the system and any auxiliary systems. Maintenance cost covers the expenses required to keep these systems and their auxiliaries in good working condition. Operational cost pertains to the costs associated with operating the systems, such as consumption and refueling, among other expenses. These three variables are defined by four sets: {VL, L, H, VH}, with respective membership functions aligned with the previously described variables. These membership functions are represented by a mixture of two Gaussian functions, as shown in Equation 1 and illustrated in Figures 8.
Parameters for Each Acoustic Signature Membership Function.
Parameters for Each Acoustic Signature Membership Function.

Subcategory Operation Membership Functions.
Parameters for Each Operational Subcategory Membership Function.

Subcategory Economic: The Same Membership Functions Represent Acquisition, Maintenance, and Operation Costs.
The fuzzy sets defined, along with the parameters

Subcategory Economic Membership Functions.
Parameters for Each Economic Subcategory Membership Function.
The environmental subcategory evaluates the system’s impact on the environment in which it operates. It includes variables such as thermal emission, by-products, and size. Thermal emission reflects the temperatures resulting from the operation of various systems. The byproduct variable refers to the type of byproduct expelled by the systems when operating, indicating whether it is polluting. The size variable pertains to the space or volume occupied by the system.
Obtaining these specific values through theoretical study was possible for the thermal emission variable. For this variable, four fuzzy sets were defined: {VL, L, H, VH}, with respective membership functions represented by a mixture of two Gaussian functions, as shown in Equation 1 and illustrated in Figure 10. Table 6 also indicates the defined fuzzy sets, parameters

Environmental Membership Functions - Thermal Emissions.
Parameters for Each Thermal Emission Membership Function.
Given the inherent subjectivity in assessing the number of byproducts generated by the operation of various systems, this variable is represented by four fuzzy sets: {Low Pollutant (LP), Pollutant (P), Highly Pollutant (HP), Extremely Pollutant (EP)}. The membership functions for these sets are defined by a mixture of two Gaussian functions, as shown in Equation 1 and illustrated in Figure 11. Table 7 lists the parameters

Environmental Membership Functions - Byproducts.
Parameters
As with many other variables, obtaining concrete values for size was impossible, making the fuzzification purely qualitative. This variable is again defined by four fuzzy sets: {VL, L, H, VH}, with respective membership functions defined by a mixture of two Gaussian functions, as shown in Equation 1 and illustrated in Figure 12. Table 8 indicates the fuzzy sets and their respective parameters

Environmental Membership Functions - Size.
Parameters for Each Size Membership Function.
Since the environmental subcategory variable is also composed of the same three fuzzy sets as the operational and economic subcategories {L, M, H}, only Figure 13 is presented, where its membership functions can be visualized. The mean values

Subcategory Environmental Membership Functions.
The recommendation level is derived as a system output variable by combining values obtained from different subcategories and after obtaining the center of mass using the centroid method described in Equation 2. This process allows us to establish the system’s recommendation level. As previously mentioned, the recommendation level is based on rules stemming from insights provided by PoN experts in submarine systems.
To implement these inferences, it is essential to define the output variable. This variable is categorized into five fuzzy sets: {VL, L, M, H, VH}, with respective membership functions defined by a mixture of two Gaussian functions, as shown in Equation 1 and illustrated in Figure 14. Table 9 indicates the fuzzy sets and their respective parameters

Recommendation Level Membership Functions.
Parameters for Each Recommendation Level Membership Function.
To improve the clarity of the final value obtained at this stage, it would be beneficial for the system to provide an output between 0% and 100%. It is important to note that due to the specific nature of the centroid defuzzification method used, normalization was applied to the output values using the following expression:
While the transformation may not be essential since the original values already serve as a ranking, with higher values indicating better recommendations, normalization enhances the clarity of results by presenting them as percentages. In decision-making, users of this system should not only focus on the final value but also compare values across the three subcategories to ensure optimal outcomes. The final value should not be the only determining factor for choosing the best system.
This section presents and justifies the values chosen for each system, based on the variables discussed in Section 4.1, the experimental setup in Section 4.2, and the discussion and analysis of the results obtained, as outlined in Section 4.3.
Data Application
For characteristics for which concrete values could not be obtained, such as power, acoustic signature, acquisition costs, maintenance costs, operation costs, byproducts, and size, a scale ranging from 0 to 10 was established. The values assigned to different systems will be qualitative and based on theoretical study and the expertise of PoN experts in submarine systems. Once again, feedback from the PoN experts was obtained from the 21 questionnaires completed by PoN militaries with submarine specialization.
It is important to note that the scales are designed to reflect real-world conditions, facilitating the system’s future adaptation with specific values. For example, when considering the acoustic signature, the values will be measured in decibels (dB), with lower values indicating better performance. This principle applies to all other variables except power, where higher levels are preferred for better outcomes. Additionally, we will adopt the worst-case scenario values for each input variable when analyzing the different systems. This approach establishes a coherent method for comparison.
Data Application: Efficiency
For the variable of efficiency, FC can achieve efficiencies exceeding 70% (Kordesch & Simader, 1996; Milliken & Ruhl, 2002), which is relatively high compared to the other three systems. SE have an efficiency of about 40% (Menon et al., 2020), ST achieve efficiencies ranging from 25% to 30% (Lus, 2001), and DE have efficiencies around 30% to 35% (Lus, 2001). Therefore, the values to be considered for efficiency are as follows: (i) FC – 70%, (ii) SE – 40%, (iii) ST – 25%, and (iv) DE – 30%.
Data Application: Power
A theoretical analysis of the variable power revealed that the steam turbine system delivers the highest power output among the four systems. In contrast, SE exhibit the lowest power. When comparing the DE and FC systems, it is concluded that DE have a greater power output than FC. Therefore, we can establish the following power hierarchy:
Data Application: Acoustic Signature
According to research conducted and insights provided by experts in PoN submarine systems, it was concluded that the least noisy system is the fuel cell system, which operates solely through electrochemical reactions and produces almost no noise. In contrast, the other three systems generate significant noise due to their moving parts. DE are particularly loud because of the explosions within them. ST also produce considerable noise from the rotation of the turbine itself. However, they are much quieter than DE. Finally, SE generate some noise from their moving parts, mainly the pistons. Still, this noise level is significantly lower than that of ST. Based on this evaluation, the following values will be adopted: (i) FC – 1, (ii) SE – 3, (iii) ST – 7, and (iv) DE – 9.
Data Application: Acquisition Cost
When analyzing the acquisition cost variable, we conclude that FC and ST are complex systems requiring highly specialized infrastructures. Consequently, their acquisition costs are significantly higher, comparable to each other, than those of the other two systems DE and SE. DE are regarded as much less complex than SE, resulting in lower acquisition costs. Based on this assessment, we will adopt the following values: (i) FC – 8, (ii) SE – 5, (iii) ST – 8, and (iv) DE – 3.
Data Application: Maintenance Cost
Based on the findings from the study, it was noted that the maintenance routines for FC are significantly less extensive than those for DE, SE, and ST. However, FC tend to have a shorter lifespan than the other systems, leading to the need for more frequent replacement of various internal components, which can be relatively costly. As a result, fuel cell maintenance costs are comparable to DE. In contrast, the maintenance cycles for DE, SE, and ST are much longer due to exposure to corrosion and contamination. However, ST have a distinct advantage over the other engine types. Additionally, SE, due to their complexity, incur higher maintenance costs than DE. Thus, to compare this variable, the following values will be adopted: (i) FC – 5, (ii) SE – 6, (iii) ST – 8, and (iv) DE – 5.
Data Application: Operating Cost
Regarding operating costs, FC are significantly more expensive than DE and SE. This is primarily due to their need for two essential components, oxygen (
Data Application: Thermal Emission
The theory behind the operation of the different systems made it possible to obtain concrete temperature data. FC operating temperatures range from 50°C to 1000°C (Lus, 2001). However, these values correspond to the Proton Exchange Membrane (PEM) FC installed in the PoN Trident class submarines and will be used for comparison. Assuming this, these cells achieve operating temperatures around 80°C. Regarding SE, they reach operating temperatures around 500°C (Sripakagorn & Srikam, 2011), ST achieve temperatures around 700°C (Lus, 2019), and DE operate around 500°C. Thus, the following values will be adopted: (i) FC – 80°C, (ii) SE – 500°C, (iii) ST – 700°C, and (iv) DE – 500°C.
Data Application: Byproducts
When discussing FC, it is essential to note that the byproduct of their chemical reactions is solely distilled water. This is advantageous for maintaining stability, as it helps regulate weight by compensating for the fuel consumption of hydrogen (
In contrast, the other three systems primarily produce carbon dioxide (
Data Application: Size
Both FC and ST are complex systems that must be implemented appropriately. ST complexity primarily arises from their size, while FC involve additional systems contributing to their complexity. SE also have some complexity, which is lessened by incorporating a specialized module into the submarine’s design. In contrast, DE are the easiest to implement among these options. Based on this assessment, the following values will be adopted: (i) FC – 7, (ii) SE – 5, (iii) ST – 7, and (iv) DE – 3.
Data Application: Variables Resume
Table 10 summarizes the values to be adopted for each variable of the systems mentioned. The accuracy of this definition will directly impact the reliability of the developed DSS. While significant technological advances may gradually influence the values of these variables, the inherent nature of the systems makes abrupt or significant changes unlikely. Nonetheless, a new propulsion system can be integrated into the analysis, ensuring the DSS remains adaptable and up-to-date.
Variables Values in the Considered Subcategories - Summary.
Variables Values in the Considered Subcategories - Summary.
After defining the variables and ensuring the reproducibility of the results, this subsection details the experimental setup for the proposed DSS, as illustrated in Figure 15. Section 4.2.1 will describe the computational tools used, while Section 4.2.2 will outline the fuzzy categories considered based on the specified variables. Section 4.2.3 will summarize the adopted inference framework, and Section 4.2.4 will discuss the defuzzification process and normalization.

Experimental Setup - Simplified Illustration.
The DSS was developed entirely in Python because of its widespread use, user-friendly nature, and strong support for scientific computing and fuzzy logic. The following Python libraries were utilized:
As outlined in Section 3.1, the analyzed variables were categorized into three fuzzy subcategories:
Each variable was modeled using a combination of two Gaussian functions, as detailed in Equation 1. The values assigned to each system were determined following the methodology presented in the preceding subsections.
As described previously, the developed DSS is based on a two-stage FMM, as illustrated in Figure 15, which can be summarized as follows:
The fuzzy inference rules used in both stages are adequately defined in Section 3.
Defuzzification and Normalization
The centroid defuzzification method outlined in Equation 4 converted fuzzy values into numerical outputs in both stages. The final recommendation score was normalized to a percentage scale (0%–100%), as noted in Equation 5, which helps interpret and compare recommendation values.
Discussion of Results
The fuzzy system described in Section 3 has been implemented. In this section, we will analyze the results by presenting the membership degrees of each input across subcategories and their classifications. We will also include the outcomes from both system phases. For a detailed analysis of output values, see Appendix A for the activated rules.
We will defuzzify the fuzzy variables using the centroid calculation method as shown in Equation 2 to complete the process and obtain the final recommendation scores. Alongside the final recommendation score for each system, we will also present its normalization using Equation 5, yielding a score between 0% and 100%.
Discussion of Results: Fuel Cells
Regarding the operational aspects of FC, the data presented in Table 11 reveals that the system rated the efficiency as Good, the power as Poor, and the acoustic signature primarily as Very Low. However, there is also a slight possibility of it being Low, with a weak membership of 0.1.
Degree of Membership and Classification - Operational - FC.
Degree of Membership and Classification - Operational - FC.
These classifications triggered the rules presented in Appendix A, resulting in the area shown in Figure 16. Thus, it can be observed that the operational subcategory scored 79.

Final Fuzzy set of the Operational Subcategory - FC.
Regarding the economic subcategory, the data presented in Table 12 reveals that the system classified the acquisition cost as High and Very High, the maintenance cost as Low and High, and the operating cost as High.
Degree of Membership and Classification - Economic - FC.
Again, after the rules presented in Appendix A were activated, a resulting colored area was obtained, as illustrated in Figure 17, showing that the economic subcategory achieved a score of 42.

Final Fuzzy Set of the Economic Subcategory - FC.
Regarding the environmental subcategory, the data presented in Table 13 reveals that the system classified the thermal emission as primarily Very Low. Still, it could also be considered Low, with a meager membership degree of 0.1. The byproduct is predominantly a Low Pollutant but could also be regarded as a Pollutant, with a membership degree of 0.1 and a High size.
Degree of Membership and Classification - Environmental - FC.
Again, after the rules presented in Appendix A were activated, a resulting area was obtained, as illustrated in Figure 18. It can be observed that the environmental subcategory achieved a score of 79.

Final Fuzzy Set of the Environmental Subcategory - FC.
Considering the scores achieved in each subcategory mentioned before {79, 42, 79}, it is possible to ascertain that each of them received the classifications described in Table 14.
Degree of Membership and Classification - Recommendation - FC.
Again, after the rules presented in Appendix A were activated, resulting in the area illustrated in Figure 19, where the recommendation level for FC reached a High level, with a score of 80.

Final Fuzzy Set of the Recommendation Level - FC.
To better understand the final score of 80, this value was normalized using Equation 5, from which it can be concluded that the fuel cell system achieved a recommendation level of 86%.
Regarding the operational subcategory of SE, based on the data presented in Table 15, it can be observed that the system classified the efficiency as Poor, the power as Poor and also Very Poor, and the acoustic signature as Low.
Degree of Membership and Classification - Operational - SE.
Degree of Membership and Classification - Operational - SE.
Again, after the rules presented in Appendix A were activated, a resulting colored area was obtained in orange, as illustrated in Figure 20, where it is observed that the operational subcategory achieved a score of 50.

Final Fuzzy Set of the Operational Subcategory - SE.
Regarding the economic subcategory, the data presented in Table 16 reveals that the system classified the acquisition cost as both High and Low. The maintenance cost was predominantly High, but it could also be considered Low, with a meager membership of 0.1. The operating cost was classified as Low.
Degree of Membership and Classification - Economic - SE.
Again, after the rules presented in Appendix A were activated, the resulting colored area was obtained, as illustrated in Figure 21, where the economic subcategory scored 50.

Final Fuzzy Set of the Economic Subcategory - SE.
Regarding the environmental subcategory, the data presented in Table 17 reveals that the system classified the thermal emission as predominantly High. Still, it could also be considered Low, with a meager membership of 0.2. Classifies the byproducts as Highly and Extremely Polluting and the size as Low.
Degree of Membership and Classification - Environmental - SE.
Again, after the rules presented in Appendix A were activated, the resulting area was obtained, as illustrated in Figure 22, where the environmental subcategory scored 42.

Final Fuzzy Set of the Environmental Subcategory - SE.
Considering the scores achieved in each of the dimensions above {50, 50, 42}, it is possible to ascertain that each received the classifications described in Table 18.
Degree of Membership and Classification - Recommendation - SE.
Again, after the rules presented in Appendix A were activated, resulting in the area illustrated in Figure 23, where the recommendation level for SE reached a Medium level, with a score of 49. Through normalization, the system achieved a recommendation level of 49%.

Final Fuzzy Set of the Recommendation Level - SE.
Regarding the operational subcategory of ST, based on the data presented in Table 19, it can be observed that the system classified the efficiency as predominantly Poor, but also as Very Poor, with a relatively low degree of membership of 0.1, the power as Good and also Very Good, and the acoustic signature as High.
Degree of Membership and Classification - Operational - ST.
Degree of Membership and Classification - Operational - ST.
Again, after the rules presented in Appendix A were activated, the resulting colored area was obtained, as illustrated in Figure 24, where it can be observed that the operational subcategory achieved a score of 48.

Final Fuzzy Set of the Operational Subcategory - ST.
Regarding the economic subcategory, the data presented in Table 20 reveals that the system classified the acquisition cost as High and Very High, the maintenance cost as High and Very High, and the operating cost as High.
Degree of Membership and Classification - Economic - ST.
Again, after the rules presented in Appendix A were activated, a resulting colored area was obtained, as illustrated in Figure 25, allowing us to observe that the economic subcategory scored 42.

Final Fuzzy Set of the Economic Subcategory - ST.
Regarding the environmental subcategory, the data presented in Table 21 reveals that the system classified the thermal emission as High and Very High, the byproduct as Highly Polluting and Extremely Polluting, and the size as High.
Degree of Membership and Classification - Environmental - ST.
Again, after the rules presented in Appendix A were activated, a resulting area was obtained, as illustrated in Figure 26, where the environmental subcategory scored 16.

Final Fuzzy Set of the Environmental Subcategory - ST.
Considering the scores achieved in each of the dimensions above {48, 42, 16}, it is possible to ascertain that each received the classifications described in Table 22.
Degree of Membership and Classification - Recommendation - ST
Again, after the rules presented in Appendix A were activated, resulting in the area illustrated in Figure 27, where the recommendation level for ST reached a Low level, with a score of 29. The system achieved a recommendation level of 24% after normalization.

Final Fuzzy Set of the Recommendation Level - ST.
Regarding the operational aspects of DE, the data presented in Table 23, it can be observed that the system classified the efficiency as Poor, the power as predominantly Good, but also as Poor, with a meager degree of membership of 0.1, and the acoustic signature as predominantly Very High, but also High, with a meager degree of membership of 0.2.
Degree of Membership and Classification - Operational - DE.
Degree of Membership and Classification - Operational - DE.
Again, after the rules presented in Appendix A were activated, a resulting colored area was obtained, as illustrated in Figure 28, where the operational subcategory scored 18.

Final Fuzzy Set of the Operational Subcategory - DE.
Regarding the economic subcategory, the data presented in Table 24 reveals that the system classified the acquisition cost as Low, the maintenance cost as Low and High, and the operating cost as Low.
Degree of Membership and Classification - Economic - DE.
Again, after the rules presented in Appendix A were activated, a resulting colored area was obtained, as illustrated in Figure 29, where the economic subcategory scored 50.

Final Fuzzy Set of the Economic Subcategory - DE.
Regarding the environmental subcategory, the data presented in Table 25 reveals that the system classified the thermal emission as predominantly High. Still, it could also be considered Low, with a meager membership of 0.2, the byproduct being Highly and Extremely Pollutant, and the size being Low.
Degree of Membership and Classification - Environmental - DE.
Again, after the rules presented in Appendix A were activated, a resulting area was obtained, as illustrated in Figure 30, where the environmental subcategory scored 50.

Final Fuzzy Set of the Environmental Subcategory - DE.
Considering the scores achieved in each of the dimensions as mentioned above {18, 50, 50}, it is possible to ascertain that each received the classifications described in Table 26.
Degree of Membership and Classification - Recommendation - DE.
Again, after the rules presented in Appendix A were activated, resulting in the area illustrated in Figure 31, the recommendation level for DE reached a Low level, with a score of 29. After normalization, the system achieved a recommendation level of 24.

Final Fuzzy Set of the Recommendation Level - DE.
To summarize all the scores across the dimensions and levels of recommendation achieved by the different systems, along with their respective normalization, Table 27 has been created.
Final Obtained Results - Resume.
Final Obtained Results - Resume.
The data shown in Figure 32 indicates that FC outperform other systems. They exhibit significantly higher operational efficiency, produce virtually no noise during operation, and have a minimal environmental impact, generating nearly zero pollution.

Starplot Illustrating the Final Obtained Results.
The final recommendation value is not the only factor in this DSS. Comparing dimension values is essential for decision-makers. Although FC systems are preferred, this analysis is critical as new systems can emerge, complicating decisions. A two-phase system provides a recommendation value while evaluating systems operationally, economically, and environmentally for thorough comparison.
An equation was applied that assigned different weights based on the opinions of PoN submarine systems experts to compare values from the fuzzy system with those from linear and non-fuzzy applications. Normalized values from the first phase were used for the operational subcategory (
Final Recommendation Level Analysis.
These results are further validated by analyzing the most recently developed systems. The FC AIP systems have been integrated into a new generation of non-nuclear submarines, including the Republic of Singapore Navy’s Type 218SG, commissioned in 2020 (Naval News, 2020), the German–Norwegian Type 212CD, which concluded its design phase in 2024 (Army Recognition, 2024a; Naval News, 2024b), and the Turkish Navy’s Type 214TN, with its lead ship commissioned in 2024 (Naval News, 2024c).
This study employed a fuzzy model instead of a linear one, offering the DSS greater flexibility and variability by assigning unlimited membership degrees. It accommodates vague human assessments, providing an effective problem-solving method and intuitively enhancing understanding.
AIP systems boost submarines’ stealth and underwater power. They enable non-nuclear submarines to compete with nuclear ones by reducing the need for battery recharge. Their compact size and low noise make AIP systems strategic assets for countries without a large nuclear fleet. This study confirms the superiority of fuel cells among AIP systems, highlighting the DSS as a key decision-making tool for AIP implementation.
A two-phase FMM approach was developed to create a DSS for selecting the best submarine AIP system. It evaluated four systems: (i) FC, (ii) SE, (iii) closed-cycle ST, and (iv) closed-loop DE. FC received the highest recommendation at about 86%, followed by SE at 49%, and both closed-cycle ST and closed-loop DE at 24%. This systematic method provided clear recommendations, allowing for comparison across key operational, economic, and environmental categories.
The implementation faced challenges in accessing technical data due to military confidentiality. Nevertheless, insights from PoN experts and theoretical studies significantly alleviated this issue. The system’s flexible nature allows for future adjustments as concrete data becomes available through collaboration with manufacturers and suppliers, enabling the expansion of this study to other systems as well. Although the DSS has been validated with naval experts, its applicability remains limited by the availability of high-fidelity operational data, potential refinements in fuzzy rule definitions, and the necessity for integration with real-time propulsion performance monitoring systems.
This method introduced an application of fuzzy logic to AIP selection, filling a void not addressed by the existing DSS frameworks. By incorporating fuzzy logic, the system effectively manages uncertainties and imprecision in expert evaluations, providing a flexible and adaptable solution for assessing AIP systems.
Recent advancements in battery technology yield higher capacities and smaller sizes than traditional lead-acid batteries used in submarines. Future research should focus on meeting Navy requirements through AIP systems and batteries, bypassing traditional generators. This could lead to more efficient submarine designs that align with modern technology and operational needs. Additionally, we should focus on enhancing the system’s adaptability to evolving propulsion technologies and exploring hybrid decision-making techniques to optimize performance further. Integrating machine learning for self-improving decision rules could further enhance the DSS, allowing it to learn from historical data and refine its recommendations over time, thereby improving decision-making accuracy and adaptability in more dynamic environments.
This DSS can provide navies with a robust, data-driven framework for evaluating and selecting AIP technologies. It aligns choices with operational, economic, and environmental demands and priorities while minimizing the risk of technical mismatches, reducing long-term costs, and enhancing the overall effectiveness of submarine missions.
Footnotes
Author Contributions
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by national funds through the Fundação para a Ciência e a Tecnologia (FCT). Nuno Pessanha Santos received support under the projects LA/P/0083/2020, UIDP/50009/2020, and UIDB/50009/2020 (DOI: 10.54499/LA/P/0083/2020, 10.54499/UIDP/50009/2020, and 10.54499/UIDB/50009/2020) – Laboratory of Robotics and Engineering Systems (LARSyS). Ricardo Moura received support under the projects UIDB/00297/2020 and UIDP/00297/2020 (DOI: 10.54499/UIDB/00297/2020 and 10.54499/UIDP/00297/2020) – Center for Mathematics and Applications.
Declaration of Conflicting Interests
The authors declare that there are no known competing financial interests or personal relationships that may have influenced the work reported in this paper.
