Next Article in Journal
Integration of Land Use Potential in Energy System Optimization Models at Regional Scale: The Pantelleria Island Case Study
Next Article in Special Issue
Optimal Configuration of Multi-Energy Storage in an Electric–Thermal–Hydrogen Integrated Energy System Considering Extreme Disaster Scenarios
Previous Article in Journal
An Aspect-Based Review Analysis Using ChatGPT for the Exploration of Hotel Service Failures
Previous Article in Special Issue
Optimizing Water-Light Complementary Systems for the Complex Terrain of the Southwestern China Plateau Region: A Two-Layer Model Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study on Effects of Flow Channel Length on Solid Oxide Fuel Cell-Integrated System Performances

1
Mechanical and Electrical Engineering College, Hainan University, Haikou 570228, China
2
Collaborative Innovation Center of Ecological Civilization, Hainan University, Haikou 570228, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(4), 1643; https://doi.org/10.3390/su16041643
Submission received: 12 January 2024 / Revised: 8 February 2024 / Accepted: 14 February 2024 / Published: 16 February 2024

Abstract

:
The structural dimensions of the SOFC have an important influence on the solid oxide fuel cell (SOFC)-integrated system performance. The paper focuses on analyzing the effect of the flow channel length on the integrated system. The system model includes a 3-D SOFC model, established using COMSOL 6.1, and a 1-D model of the SOFC-integrated system established, using Aspen Plus V11. This analysis was conducted within an operating voltage range from 0.4 V to 0.9 V and flow channel length range from 6 cm to 18 cm for the SOFC-integrated system model. Performance evaluation indicators for integrated systems are conducted, focusing on three aspects: net electrical power, net electrical efficiency, and thermoelectric efficiency. The purpose of the paper is to explore the optimal flow channel length of SOFC in the integrated system. The results indicate that there is inevitably an optimal length in the integrated system at which both the net electrical power and net electrical efficiency reach their maximum values. When considering the heat recycling in the system, the integrated system with a flow channel length of 16 cm achieves the highest thermoelectric efficiency of 65.68% at 0.7 V. Therefore, there is a flow channel length that allows the system to achieve the highest thermoelectric efficiency. This study provides optimization ideas for the production and manufacturing of SOFCs from the perspective of practical engineering applications.

1. Introduction

In recent years, due to rapid environmental changes, researchers have been searching for clean energy that is more environmentally friendly and sustainable. Solid oxide fuel cell (SOFC)-integrated systems have attracted attention from various countries due to their high power and efficiency advantages [1,2]. By integrating thermal cycle components and SOFCs into a single system, fuel utilization efficiency can be effectively improved, and the issue of residual energy in the waste gas after the SOFCs’ reaction can be addressed [3,4]. At the same time, the optimal configuration and operating circumstances of each component in the integrated system have a significant impact on system performance [5]. Currently, the research is mostly focused on increasing fuel cell performance indicators such as power and electrical efficiency by adjusting SOFC structural factors and operating circumstances. In practical applications, the power consumption of system components will also change with the variation in SOFCs’ structural factors and operating circumstances. Therefore, achieving the optimal SOFC performance requires considering their operation in actual systems to better meet the practical needs of industrial production.
Through previous research, it is known that structural parameters are one of the important factors affecting SOFCs [6,7]. In order to analyze the influence of structural dimensions on the performance of SOFCs, numerous scholars have developed simulation models of SOFCs in various states [8,9,10]. Cui et al. developed a two-dimensional model of a flat plate tube segmented tandem solid oxide fuel cell (FTS-SOFC) to analyze the impact of structural dimensions on the uniformity of cell voltage. It was found that modifying the length of the runner greatly improved the voltage uniformity. The 5 mm FTS-SOFC exhibits the highest maximum power density when other cell dimensions were held constant [11]. Liu et al. constructed various runner structure planar direct ammonia SOFC (DA-SOFC) models to investigate the impact of runner structure on the thermoelectric performance of SOFC bipolar plates. It was shown that the cell performance, such as power density, is optimal for SOFCs with an elliptical channel structure. It was found that increasing the width of the channel cross-section while decreasing the rib width can improve the performance of the cell [12]. An et al. fabricated a segmented-in-series solid oxide fuel cell (SIS-SOFC) with varying cell-to-cell distances to investigate the impact of cell-to-cell distance on voltage loss. The study demonstrated that the ohmic resistance decreased slightly with increasing cell-to-cell distance, while the polarization resistance decreased to a greater extent. The open-circuit voltage maximum power output density and open-circuit voltage were highest when the cells were spaced 2 mm apart [13]. Jee et al. developed several models of SOFCs with varying electrolyte layer thicknesses to assess the influence of electrolyte thickness structural parameters on SOFC single cells. They demonstrated that as the thickness of the electrolyte decreases, the cell’s ohmic loss also decreases, leading to an increase in power density. An SOFC stack with an electrolyte thickness of 80 μm demonstrates a high power density. Furthermore, elevated operating temperatures decrease hydrogen consumption [14]. Mushtaq et al. fabricated planar tubular SIS-SOFCs with varying cathode thicknesses to study the impact of cathode thickness on SOFC performance. The study reveals that the optimal thickness of the cathode is 57 mm. Furthermore, it is suggested that a composite cathode with a defined perimeter and geometry may be a potential method to enhance the performance of flat tubular SIS-SOFCs [15]. The above study focuses on the influence of different structural parameters of SOFCs on the output power and efficiency of SOFCs. However, performance analyses from a systemic perspective are lacking.
In practical production, SOFCs are often used as part of an integrated system, rather than being used alone. Many scholars have conducted a parameter analysis of integrated systems for SOFCs [16]. Considering the high construction costs of SOFC-integrated systems, mathematical modeling and simulation methods are generally used [17]. Zhao et al. developed a model of a solid oxide fuel cell–gas turbine (SOFC-GT) system to gain a comprehensive understanding of the optimal integration of the SOFC, gas turbine, and other system components. The model can be enhanced for stability and efficiency by optimizing parameters such as current density, operating temperature, and the fuel utilization factor of the fuel cell [18]. Huang et al. used Simulink and Ebsilon to conduct homeostatic simulations of the SOFC-GT system. The study found the system efficiency increases with operating temperature, while system efficiency may decrease due to the higher fuel utilization [19]. Ameri et al. developed an SOFC-GT model based on the literature from Siemens-Westinghouse company using Aspen Plus. The results of the simulation show that the electricity generation efficiency within the cell tends to increase and then decrease as the fuel utilization rate improves. With minimum fuel input at a 0.85 fuel utilization rate, the efficiency is close to its maximum value [20]. However, they only used a one-dimensional simulation with non-adjustable parameters, which has limitations. In current research on SOFC-integrated systems, most of the analysis focuses on the impact of fuel utilization on the system efficiency and other performance parameters of fuel cell systems. There is relatively little analysis of the impact of SOFCs’ structural parameters on the system. At present, SOFC-integrated systems primarily rely on module modeling using commercial software like Aspen Plus, but this approach has certain limitations [21,22]. The system model created by commercial software is challenging to parameterize for a specific component in SOFC, and it is difficult to simulate certain specific functions. Therefore, it needs to be combined with other simulation models that have adjustable parameters to create more realistic system models.
Therefore, this study establishes an integrated system for SOFCs, taking into account the losses of various subsystems, including the SOFCs. The paper specifically focuses on analyzing the impact of the flow channel length of a single cell on the integrated system. A 3D model of SOFCs with different channel lengths was established based on COMSOL 6.1 [23,24]. The aim was to investigate the optimal flow channel length to achieve the best integrated system performance at different cell voltages. This study provides optimization ideas for the production and manufacturing of SOFCs from the perspective of practical engineering applications.

2. Description of SOFC-Integrated System

The proposed integrated system proposed mainly comprises an SOFC stack, an air compressor, an afterburner, a fuel compressor, a water pump, a mixer, and three heat exchangers, as shown in Figure 1. The SOFC operating pressure is standard atmospheric pressure. The SOFC system’s operational process is as follows.
Air temperature is increased to the operating temperature using heat exchanger 1 after being pressurized by the air compressor, and then passed into the SOFC cathode. Water temperature is increased to the operating temperature using heat exchanger 2 after being pressurized by the water pump. Hydrogen temperature is increased to the operating temperature using heat exchanger 3 after being pressurized by the fuel compressor. The mixer mixes heated and pressurized hydrogen and water vapor. The moist hydrogen is sent to the SOFC anode, where it undergoes a chemical reaction with the oxygen, producing an electric current. The unreacted excess gases at the anode and cathode outlets, such as hydrogen, oxygen, and nitrogen, are combined and directed into a combustion chamber for complete combustion, producing a significant amount of heat. The heat is utilized to preheat the hydrogen, water, and air, thus achieving the SOFC operating temperature.

2.1. Description of the SOFC Model

Utilizing COMSOL 6.1 software, this paper establishes a 3-D, non-isothermal, numerical model for a single SOFC extracted from the SOFC stack, as depicted in Figure 2a. Air and fuel travel in the flow channel in a counter-current fashion. The model includes equations for momentum, charge, energy, mass, and species conservation, and fully considers the electrode kinetics of electrochemical reactions. To simplify the calculations, the metal interconnects were removed, and the heat transfer around the flow channels was replaced with open periodic heat transfer boundary conditions. At the same time, the anode electrode was grounded while the cathode electrode was connected to the cell potential. The model geometric dimensions and certain boundary conditions are shown in Figure 2b.
Based on the relevant literature data [14,25], the geometric dimensions of the SOFC are shown in Table 1. Strontium-doped lanthanum manganese oxide/yttrium-stabilized zirconia (LSM-YSZ) is used as the cathode material, and nickel/yttrium-stabilized zirconia (Ni YSZ) is used as the anode material. The electrolyte material unit’s mass consists of 92% zirconia (ZrO2) and 8% yttrium-stabilized zirconia (YSZ). This model is established on the COMSOL 6.1 platform, which includes a fuel cell electrochemical model, a porous material transport model, and a heat transfer model. In a bid to better analyze its main issues, this study proposes the following methodology [26,27]:
(1)
The SOFC unit cell operates steadily, and the electrochemical reactions reach equilibrium.
(2)
The cathode and anode are constructed from homogeneous and isotropic porous materials.
(3)
In the anodic electrochemical reaction, only hydrogen gas participates.
(4)
SOFCs maintain consistent flow rates at the cathode and anode inlets over different lengths, and the inlet gas pressure is supplied by a compressor.
(5)
Gas pressure at the anode and cathode outlet is equivalent to the standard atmospheric pressure in the SOFC. Furthermore, the model operates at standard atmospheric pressure.
(6)
The SOFC gas inlet temperature and operating temperature are both set at 800 °C.
(7)
Anode and cathode outlet temperatures are computed from the model.
(8)
All gases are ideal gases in SOFC.

2.1.1. Electrochemical Reaction

Unlike the combustion power generation of internal combustion engines, hydrogen, and oxygen generate electricity directly through electrochemical reactions in SOFC. Oxygen gets electrons from the anode and undergoes a reduction reaction at the cathode, transforming into oxygen ions. The oxygen ions undergo an oxidation reaction with hydrogen after passing through the electrolyte layer, producing electrons and water vapor. This process is accompanied by the movement of electrons and ions, which generates a current. When the reaction reaches equilibrium, the cell current density stabilizes at a constant value, and the voltage can be determined accordingly.
In practice, the actual cell potential is less than the theoretical equilibrium potential due to the effects of various types of cell polarization. The three main polarization losses are activation overpotential, concentration overpotential, and ohmic overpotential. The actual output voltage is shown in Equation (1) [26,28,29].
V cell = E r V pol = E r V act V ohm V con
where Er is the SOFC reversible voltage, Vcell is the SOFC actual voltage, Vpol is the total polarization voltage, Vact is the activation voltage, Vcon is the concentration voltage, and Vohm is the ohmic voltage.
In an SOFC, the open circuit voltage (OCV) or the reversible Nernst voltage represent the maximum theoretical voltage that the cell can attain under certain conditions. The reversible voltage of the SOFC is shown in Equations (2) and (3).
E r = E 0 + R g × T n e × F ln p H 2 × p O 2 1 / 2 p H 2 O
E 0 = 1.317 2.769 × 10 4 T
where E0 is the standard electromotive force at atmospheric pressure, Rg is the gas constant (Rg = 8.314 J/(mol·K)), T is the SOFC operating temperature (T = 800 °C), F is the Faraday constant (F = 96485 C/mol), and p is the fractional pressure for each gas component. ne is the number of electrons that are transferred (ne = 2).
Activation polarization is the energy required to overcome the activation energy barrier and allow for an electrochemical reaction to proceed. In the model, activation polarization can be described in relation to the cell current using the Butler–Volmer formula. The Butler–Volmer equation for cathode and anode is shown in Equations (4) and (5) [30,31].
i a n = i 0 , a n exp α a a × n e × F × V a n , a c t R g × T exp ( 1 α a a ) n e × F × V a n , a c t R g × T
i c a = i 0 , c a × exp α c a × n e × F × V c a , a c t R × T exp ( 1 α c a ) n e × F × V c a , a c t R × T
where ian and ica represent the anodic and cathodic current densities. α a a and α c a represent the anodic transfer coefficient of the anode and cathode. Van,act and Vca,act represent anode and cathode activation polarization. i0,an and i0,ca represent the exchange current densities at the anode and cathode. The cell exchange current density is governed by the mass action law and is expressed as follows:
i 0 , a n = R g × T × K a n n e × F exp E a n R g × T
i 0 , c a = R g × T × K c a n e × F exp E c a R g × T
where Ean and Eca represent the anodic and cathodic preindex factor. Kan and Kca represent the anodic and cathodic activation energy.
Since the SOFC has a high current density in this model, the calculation of the activation loss is simplified by approximation with the Tafel formula [32,33].
V a c t = V a n , a c t + V c a , a c t = R g × T α a a × n e × F ( arcsin h ( i a n 2 i 0 , a n ) ) + R g × T α c a × n e × F ( arcsin h ( i c a 2 i 0 , c a ) )
The concentration overvoltage can be determined from the following equation [12,34]:
V c o n = V a n , c o n + V c a , c o n = R g × T 2 F ln ( p H 2 × p H 2 O 0 p H 2 O × p H 2 0 ) + R g × T 4 F ln ( p O 2 p O 2 0 )
where p H 2 , p O 2 and p H 2 O represent the fractional pressure of hydrogen, oxygen, and water in working condition. p H 2 0 , p O 2 0 and p H 2 O 0 represent the fractional pressure of hydrogen, oxygen, and water in the standard state.
During SOFC operation, the flow of electrons in the circuit and the flow of ions in the electrolyte cause ohmic losses, which are calculated as follows [35]:
V o h m = i R o h m
where Rohm represents the approximate total resistance.
Current density is a crucial parameter for evaluating cell performance. The SOFC model current density is described by Faraday’s law [4,36,37].
i = n e × F × Δ v f u e l N c e l l × A c e l l
where Ncell is the amount of SOFC cells in the stack; Δvfuel is electrochemical reaction level (in mol/s); Acell is the active area of a single cell (in m2).
The relationship between electrochemical reaction level and fuel utilization can be expressed by the following equation:
Δ v f u e l = v f u e l × U f u e l
U f u e l = i × N c e l l × A c e l l n e × F × v f u e l
Ufuel is the SOFC fuel utilization rate; vfuel is the hydrogen molar flow rate at the anode (in mol/s).
The single cell output current can be represented as follows:
I = i × A c e l l
The formula for calculating the SOFC stack output power is as follows:
W S O F C = I × V c e l l × N c e l l
As shown in Equations (9)–(13), the expression for the SOFC stack output power is as follows:
W S O F C = V c e l l × U f u e l × v f u e l × n e × F
To facilitate subsequent data analysis, we introduce the concept of electricity generation efficiency. The efficiency of electricity generation is described as the ratio of the SOFC stack power output to the fuel energy consumed by the SOFC over a specific period. The expression for electricity generation efficiency is shown in Equation (17).
ξ e l e = W S O F C Δ v f u e l × Δ h L H V = V c e l l × n e × F Δ h L H V
where ΔhLHV is the hydrogen’s lower heating value (241.8 kJ∙mol−1).
The power generation efficiency of SOFC is represented by Equation (18).
η e l e = W S O F C v f u e l × Δ h L H V = ξ e l e × U f u e l

2.1.2. Model Parameter Settings

The SOFC model consists of a hydrogen fuel cell module, a free and porous media flow channel module, a solid and fluid heat transfer module, and a multi-physics field coupling module. Table 2 lists the boundary conditions used in the process of solving the control equations. Table 3 shows the relevant data for the SOFC model in COMSOL 6.1 [38].
It is essential to consider the gas diffusion and the properties of the anode and cathode substrates in the free and porous media flow module. The electrodes’ substrate parameters are listed in Table 4. The model considers compressible flow (Ma < 0.3) and Darcy flow but does not account for turbulent flow.
In solid and fluid heat transfer models, it is necessary to define the material properties of the anode and cathode substrates. The material property parameters are depicted in Table 5.

2.1.3. SOFC Model Validation

This paper validates the model using experimental data from the literature. The geometric parameters, material characteristics, fuel characteristics, and SOFC model operating conditions used in this paper are consistent with those in the references [25,26]. Figure 3 compares the simulated current–voltage curve with the experimental current–voltage curve. There is a discrepancy in the results of the comparison, and this discrepancy is due to the fact that there is still a gap between the ideal environment of the model setup and the actual experimental conditions. Generally speaking, the error will be larger in the high-current-density region, which is caused by the temperature difference resulting from ohmic heating during the experiment. This pertains to a general trend and does not impact the accuracy of our results. Therefore, the established SOFC model meets the actual operational needs of SOFCs.

2.2. Description of Other Equipment Models in the Integrated System

This paper also describes the properties of other devices in the integrated system. The other equipment mainly comprises heat exchangers, air compressors, fuel compressors, water pumps, mixers, and combustion chambers. To prevent the anode in the SOFC from being oxidized by oxygen, a specific level of hydrogen partial pressure must be maintained at the outlet. Therefore, unreacted hydrogen gas will be present in the anode exhaust gas. The exhaust gases are directed into the combustion chamber to be combusted. The exhaust gas, after combustion, will pass through three heat exchangers consecutively to heat hydrogen, water, and air. This process will fully utilize waste heat, thereby improving the entire system’s energy efficiency. To simplify the simulation process, we assume that the various auxiliary modules in the system operate at an isothermal steady state and that the battery system is fully enclosed.
The three heat exchangers in the system are used to heat air, water, and hydrogen to the specified temperatures. The heat exchangers used in the system are countercurrent heat exchangers, and the amount of total heat loss stays the same and is not affected by changes in load. The expressions for the heat exchanger model are as follows, using the method described in the literature [45]:
ε = M ˙ × C c ( T o u t , c T i n , c ) M ˙ × C min ( T i n , h T i n , c )
where Cc and Cmin represent the gas-specific heat capacity in the cold stream and the gas-specific heat capacity of the smaller fluid between the hot and cold streams at constant pressure (in J/g·K). M ˙ is the gas mass flow rate (in g/s). Tin,c and Tout,c represent the temperatures of the cold inlet stream and outlet stream of the compressor (in °C). Tin,h represents the temperature of the hot inlet stream of the compressor. Additionally, when heating water, the latent heat of vaporization for water is considered.
The heat exchange power is strongly dependent on the performance of the equipment and the gas. In an ideal situation, neglecting heat loss, the energy exchange in the heat exchanger is conserved through the change in the thermal energy of the gas. According to this principle, the heat exchanger power can be determined using Equation (20).
W h = U h × S h × ( T o u t , c T i n , c ) = M ˙ × C c ( T o u t , c T i n , c )
where Sh represents the heat transfer area of the heat exchanger (in m2) and Uh is the convective heat transfer coefficient (in W/(m2∙K)).
The fuel compressor and air compressor use the same compressor model. It is assumed that the compressed gas undergoes isentropic compression in an ideal state. Equation (21) represents the formula used to calculate the outlet gas pressure.
p 2 = p 1 × π C
where p2 represents the compressor outlet pressure (in Pa) and p1 represents the compressor inlet pressure. πC represents compressor pressure ratio.
According to the definition of isentropic efficiency, the actual outlet temperature can be calculated by Equation (22):
η c = T 2 T 1 T 2 T 1
where ηc is the compressor isentropic efficiency, T2* is the ideal outlet temperature of the compressor, T1 is the compressor inlet temperature, and T2 is the compressor outlet temperature. The ideal outlet temperature can be calculated using Equation (23).
T 2 = T 1 × π C γ 1 γ
where γ represents the adiabatic index of the compressor.
The power consumed by the compressor is determined using the following equation:
W c = M ˙ × c p × ( T 2 T 1 )
Due to its high density at standard conditions and low flow rate in the pump, water requires relatively low energy to be pressurized to a specified pressure. The pump power consumption is represented by Equation (25).
W p = M h 20 × v i n p o u t p i n η
where Mh2o is the water mass flow rate (in kg/s); Pout and Pin refer to the outlet and inlet pressures of the water pump (in Pa); η is the water pump efficiency; v is the specific volume of water (in m3/kg).
The combustion chamber is considered to be an ideal combustion chamber, where all the hydrogen gas at the inlet is converted into water. It is also assumed that there is no residual fuel hydrogen at the combustion chamber outlet, and the potential side reactions and oxidative corrosion are not taken into consideration. The combustion chamber is modeled as adiabatic, which means that the enthalpy remains constant at the inlet and outlet, and can be expressed as follows:
H o u t T o u t , n o u t = H i n T i n , n i n
Based on the above conditions, an SOFC-integrated system is constructed with a Vcell between 0.4 V and 0.9 V and an Lcell between 6 cm and 18 cm.

3. System Performance Evaluation Indicators

In our effort to accurately assess the system’s performance, we established several evaluation metrics, including net electrical power, net electrical efficiency, and thermoelectric efficiency. Among them, the net system electrical power is described as the difference between the SOFC stack’s output power and the Balance of Plant (BOP) component’s power consumption (auxiliary equipment power consumption), as expressed in Equations (27) and (28).
W s y s = W S O F C W B O P
W B O P = W p + W c , a i r + W c , f u e l
The system’s net electrical efficiency is described as the ratio of the system’s net electrical output to the low calorific value (LCV) of fuel that enters the system. The representation is shown in Equation (29). This indicator represents the system’s output performance.
η s y s , n e t = W S O F C W B O P v f u e l × Δ h L H V
The thermoelectric efficiency of a system is described as the ratio of the system’s net electricity and the thermal energy recovered by the system to the LCV of the fuel that enters the system. The recovered heat energy is mainly achieved through the use of a heat exchanger. The expression is given by Equations (30) and (31). This indicator mainly reflects the efficient use of energy in the system.
η s y s , t h e r = W S O F C + W h W B O P v f u e l × Δ h L H V
W h = W h , 1 + W h , 2 + W h , 3

4. Results and Analysis

This paper is based on the COMSOL 6.1 and ASPEN software. It employs the “3-D numerical simulation of the SOFC stack + 1-D numerical simulation of the system” method to research the effect of different flow channel lengths on the properties of individual components and integrated systems. The paper establishes seven different lengths of SOFC single-cell models, at an Lcell ranging from 6 to 18 cm, and analyzes their cell performance at Vcell ranging from 0.4 V to 0.9 V. The method of controlling variables was used to conduct analysis and discussion, ensuring the validity of the performance evaluation.

4.1. Influence of Different Flow Channel Lengths on SOFC Stack

Figure 4 illustrates the voltage–current and power density profiles of a single cell for various flow channel lengths. Under the same flow channel length, the voltage and current density outputs show an inverse proportional trend, while the power density exhibits a first-increasing-then-decreasing trend with the increase in current density. These findings are consistent with the existing literature [38]. Moreover, as the flow channel length increases, the polarization phenomenon becomes more pronounced, which is consistent with the single-cell performance characteristics. In a single SOFC, the fuel in the anode flow channel is continuously consumed. As a result, the molar flow rate of the fuel gradually decreases, leading to a reduction in the extent of its reaction and a subsequent decrease in the output power per unit area. Therefore, the power density of the SOFC decreases as the size of the runner increases.
Theoretically, increasing the flow channel length will enlarge the reaction area inside the cell. This can improve fuel utilization and subsequently increase the cell output power, assuming a consistent fuel flow rate. The overall power output curve of the SOFC stack is shown in Figure 5. It can be seen from the curve that the output power of the stack increases with the increase in the flow channel length at the same operating voltage. However, the spacing between the curves indicates a gradual decrease in the rate of increase in output power. This decrease is directly related to fuel utilization during the cell reaction process [46].
Equation (14), which was used for the SOFC output power mentioned above, indicates that when the amount of fuel that is input into the fuel cell remains constant, the SOFC output power is directly proportional to the product of fuel utilization and cell voltage. Therefore, numerical simulations were carried out to analyze the fuel utilization performance at various flow channel lengths, as depicted in Figure 6. The curves indicate that, under the same cell voltage, the fuel utilization rate follows a similar trend to the stack output power as the flow channel length varies. At a working voltage of 0.6 V, the fuel utilization of the SOFC stack increased as the flow channel length increased from 6 cm to 12 cm, rising from 37.45% to 61.87%. This represents a 100% increase in length, while only a 65.2% increase in fuel utilization was recorded. As the flow channel length increases, fuel utilization also increases, but the rate of growth in fuel utilization decreases.

4.2. Influence of Different Flow Channel Lengths on Compressor

The gas pressure in the gas flow path is generated by a fuel compressor and an air compressor, connected to the anode and cathode inlets, respectively. In the SOFC, the flow channel length and the cell operating voltage will affect the gas pressure at the inlets. Through numerical simulation, the air and gas compressors in the integrated system have pressure ratios and power that vary depending on the flow channel length and operating voltage. This is illustrated in Figure 7 and Figure 8.
It is evident from the curves that both the pressure ratio and power consumed by the fuel and air compressors increase as the flow channel length increases, with the cell voltage remaining constant. To ensure a consistent outlet pressure for fuel cells with varying flow channel lengths, the longer cells require higher air and fuel inlet pressures at the entrance to compensate for the channel-generated pressure drop.
Under the same flow channel length, the ratio of air and fuel compressor pressures and power decrease as the operating voltage increases. However, the decrease in pressure ratio and power is slower with shorter flow channel lengths. This is because, when the SOFC stack operates at a lower cell voltage, it is in a higher current operating condition, leading to higher polarization losses in the cell. From Equation (11), it can be seen that the cell current is directly proportional to the fuel consumption, leading to a higher consumption of hydrogen and oxygen. As a result, there is a greater pressure drop, necessitating higher pressure at the inlets provided by the compressors.

4.3. Influence of Different Flow Channel Lengths on Heat Exchangers

The function of heat exchangers in integrated systems is to efficiently harness the residual heat from the exhaust to elevate the temperatures of the fuel and air to the operating level. Heat exchangers 1–3 represent air, water, and hydrogen heat exchangers, respectively. The relationship between power, flow channel lengths, and operating voltages is shown in Figure 9.
From the curve, it can be observed that there is an inverse relationship between the power consumption of each heat exchanger and the flow channel length at the same operating voltage. The reason for this is that the longer the flow channel length in an integrated system, the more work is required by the compressor to compress the fuel and air. Therefore, the fuel and air attain higher enthalpy values and temperatures at the compressor outlet. The relationship between the outlet temperature, various flow channel lengths, and the operating voltages of the gas and air compressors is illustrated in Figure 10. Therefore, the energy needed to heat the fuel gas and air to 800 °C decreases, resulting in a reduction in power consumption for the heat exchanger. Additionally, Figure 9b shows that the power consumption of the water heat exchanger is minimally influenced by the length and voltage of the cell. This is because when water is pressurized by the water pump at 25 °C, the increase in enthalpy and temperature is minimal. As a result, there is little difference in the energy required for the water to be heated to 800 °C after passing through the heat exchanger.

4.4. System Performances Analysis

From the above analysis, we found that the different flow channel lengths of the SOFC have different impacts on the various components of the integrated system. Therefore, we boldly ask whether there is an optimal length dimension that would result in the integrated system having the best performance. We analyzed the system’s optimal performances based on three indicators: net electrical power, net electrical efficiency, and thermoelectric efficiency.

4.4.1. Maximum Net Electrical Power of the System

To minimize the need for external energy inputs, in the integrated system, the SOFC stack supplies electrical energy to other equipment components. The system output power is derived as the difference between the SOFC output power and the power consumed by other energy-consuming equipment in the BOP. The relationship curve between BOP power consumption and different flow channel lengths and operating voltages is shown in Figure 11a, and the relationship curve between the system net electrical power and different flow channel lengths and voltages is shown in Figure 11b. In Figure 11b, when the operating voltage is 0.7 V, the length of the cell ranges from 6 cm to 18 cm, and the system net electrical power shows a gradually increasing trend, but at a slower rate. As the length increases, the SOFC output power also increases. The power consumed by the air compressor also increases simultaneously. However, the growth rate of stack output power is slower than the growth rate of blowout preventer consumption power, leading to a slower growth rate of net system power.
According to Figure 11, we can further determine the correlation between the flow channel length and the system’s net electrical power, as shown in Figure 12. At 0.4 V, the overall system power exhibits a rising and then falling trend. The maximum net electrical power is 471.84 W when the length is 16 cm, with a fuel utilization rate of 88.93%. This means that a longer flow channel length does not necessarily result in a better performance. Furthermore, it can be inferred that there is an optimal length value that maximizes the system’s net electrical power at various voltages.

4.4.2. Maximum Net Electrical Efficiency of the System

On the basis of the previous expressions for the system performance evaluation metrics, the highest electrical efficiency corresponding to the maximum system output power that can be achieved by the SOFC-integrated system is noted when the energy input to the system is kept constant. Therefore, the tendency of the system’s net electrical efficiency mirrors the tendency of the system’s net electrical power, as depicted in Figure 13. As the flow channel length increases, the system’s net electrical efficiency demonstrates a rising but decelerating trend.
Similarly, the changing trend of the system’s net electrical efficiency with flow channel length is further illustrated in Figure 14. At 0.4 V, the net electrical efficiency shows an increasing trend with the increase in flow channel length, followed by a decrease. At a length of 16 cm, the maximum net electrical efficiency is 13.25%. It can be inferred that there is an optimal length value for different voltages that maximizes the net electrical efficiency of the system.

4.4.3. Maximum Thermoelectric Efficiency of the System

From the above discussion, it is evident that, when considering only the power aspect, the system’s electrical efficiency is relatively low. However, when considering the utilization of waste heat from the fuel exhaust gases, the integrated system energy efficiency is improved. The thermoelectric efficiency of the integrated system, obtained through numerical simulation, with the variation in flow channel length, is depicted in Figure 15.
In Figure 15, it can be observed that 0.6 V and 0.7 V are close to the optimal operating voltage. With the increase in flow channel length, the system thermoelectric efficiency tends to increase and then decrease. However, the optimal length of a single battery corresponding to the highest thermal efficiency varies at different battery voltages. At 0.6 V, the thermoelectric efficiency of the integrated system is highest when the length is 12 cm, reaching 64.62%. At 0.7 V, the system thermoelectric efficiency is highest when the length is 16 cm, reaching 65.68%. Additionally, the fuel utilization rate is 65.93%, which aligns with the practical engineering situation of SOFCs. At this point, the system’s net electrical power output is 816.19 W. This is because, as the flow channel length increases, the fuel utilization rate of the SOFC also increases, thereby increasing the SOFC output power. Although the energy in the SOFC exhaust gas decreases at this point, there is still a significant amount of energy available to meet the requirements for fuel and air preheating. Longer stacks require higher air and fuel pressure at the inlet, which inevitably increases the compressor power consumption. As a result, the heat transfer power of the heat exchanger is reduced, leading to a decrease in the energy-utilization efficiency of the exhaust gas. Therefore, there is an optimal length at different voltages that maximizes the system’s thermoelectric efficiency.

5. Summary

This paper investigates an integrated system based on SOFCs. An SOFC was simulated using the multiphysics simulation software COMSOL 6.1 and verified by experimental results. The paper analyzed the impact of the length of a single SOFC cell on system performance, and the results are presented as follows.
Under the same cell voltage, the single flow channel length is directly proportional to the pressure ratio and power consumption of the air and fuel compressors. However, it is inversely correlated with the heat transfer power of each heat exchanger. For an SOFC stack, the variation in single flow channel length affects the SOFC output power and efficiency. Increasing the flow channel length can improve fuel utilization to some extent, thereby enhancing the SOFC output power. For the integrated system, it is inferred that there is an optimal length that maximizes the net electric power and net electric efficiency. When considering heat recycling in the system, an integrated system with a single flow channel length of 16 cm achieves the highest thermoelectric efficiency of 65.68% at 0.7 V. Therefore, the flow channel length is an important factor affecting the performance of the integrated system in real industrial production. This study will provide a theoretical foundation for the further optimization of SOFC power generation systems. In the following study, the proposed numerical model will be combined with a more realistic production-integrated system to design a superior-performing SOFC-integrated system.

Author Contributions

Methodology, software, and writing—original draft preparation, Y.L.; conceptualization, funding acquisition, and supervision, J.L.; investigation, writing—review and editing, L.F.; validation, formal analysis, Q.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Hainan Provincial Natural Science Foundation of China: 520RC540 and 521RC492; and the Collaborative Innovation Center Project of Hainan University: XTCX2022STC16.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumar Yadav, A.; Sinha, S.; Kumar, A. Comprehensive review on performance assessment of solid oxide fuel cell-based hybrid power generation system. Therm. Sci. Eng. Prog. 2023, 46, 102226. [Google Scholar] [CrossRef]
  2. Vinchhi, P.; Khandla, M.; Chaudhary, K.; Pati, R. Recent advances on electrolyte materials for SOFC: A review. Inorg. Chem. Commun. 2023, 152, 110724. [Google Scholar] [CrossRef]
  3. Meng, Q.S.; Han, J.T.; Kong, L.J.; Liu, H.; Zhang, T.; Yu, Z.T. Thermodynamic analysis of combined power generation system based on SOFC/GT and transcritical carbon dioxide cycle. Int. J. Hydrog. Energy 2017, 42, 4673–4678. [Google Scholar] [CrossRef]
  4. Deng, M.L.; Liu, J.Y.; Zhang, X.S.; Li, J.S.; Fu, L.R. Energy and Parameter Analysis of SOFC System for Hydrogen Production from Methane Steam Reforming. J. Therm. Sci. 2022, 31, 2088–2110. [Google Scholar] [CrossRef]
  5. Zeng, Z.Z.; Qian, Y.P.; Zhang, Y.J.; Hao, C.K.; Dan, D.; Zhuge, W.L. A review of heat transfer and thermal management methods for temperature gradient reduction in solid oxide fuel cell (SOFC) stacks. Appl. Energy 2020, 280, 19. [Google Scholar] [CrossRef]
  6. Atif Mahmood, M.; Nawaz Chaudhary, T.; Farooq, M.; Salman Habib, M.; Maka, A.O.M.; Usman, M.; Sultan, M.; Shiung Lam, S.; Chen, B. Sensitivity analysis of performance and thermal impacts of a single hydrogen fueled solid oxide fuel cell to optimize the operational and design parameters. Sustain. Energy Technol. Assess. 2023, 57, 103241. [Google Scholar] [CrossRef]
  7. Zhu, P.; Wu, Z.; Yang, Y.; Wang, H.; Li, R.; Yang, F.; Zhang, Z. The dynamic response of solid oxide fuel cell fueled by syngas during the operating condition variations. Appl. Energy 2023, 349, 121655. [Google Scholar] [CrossRef]
  8. Iliev, I.K.; Gizzatullin, A.R.; Filimonova, A.A.; Chichirova, N.D.; Beloev, I.H. Numerical Simulation of Processes in an Electrochemical Cell Using COMSOL Multiphysics. Energies 2023, 16, 7265. [Google Scholar] [CrossRef]
  9. Cai, W.; Yuan, J.; Zheng, Q.; Yu, W.; Yin, Z.; Zhang, Z.; Pei, Y.; Li, S. Numerical Investigation of Heat/Flow Transfer and Thermal Stress in an Anode-Supported Planar SOFC. Crystals 2022, 12, 1697. [Google Scholar] [CrossRef]
  10. Mueller, M.; Klinsmann, M.; Sauter, U.; Njodzefon, J.-C.; Weber, A. High Temperature Solid Oxide Electrolysis—Technology and Modeling. Chem. Ing. Tech. 2023, 96, 143–166. [Google Scholar] [CrossRef]
  11. Cui, D.A.; Ji, Y.L.; Chang, C.; Wang, Z.; Xiao, X.; Li, Y.T. Influence of structure size on voltage uniformity of flat tubular segmented-in-series solid oxide fuel cell. J. Power Sources 2020, 460, 6. [Google Scholar] [CrossRef]
  12. Liu, Y.; Xu, Y.; Liu, J.; Xiang, M.; Sun, B.; Ya, Y.; Guo, Z.; Cheng, X. Effects of bipolar plate flow channel configuration on thermal-electric performance of direct ammonia solid oxide fuel cell: Part I—Numerical study of channel section geometry. Int. J. Hydrog. Energy 2024, 50, 765–785. [Google Scholar] [CrossRef]
  13. An, Y.-T.; Ji, M.-J.; Hwang, H.-J.; Eugene Park, S.; Choi, B.-H. Effect of cell-to-cell distance in segmented-in-series solid oxide fuel cells. Int. J. Hydrog. Energy 2015, 40, 2320–2325. [Google Scholar] [CrossRef]
  14. Park, J.M.; Kim, D.Y.; Baek, J.D.; Yoon, Y.J.; Su, P.C.; Lee, S.H. Effect of Electrolyte Thickness on Electrochemical Reactions and Thermo-Fluidic Characteristics inside a SOFC Unit Cell. Energies 2018, 11, 473. [Google Scholar] [CrossRef]
  15. Mushtaq, U.; Kim, D.-W.; Yun, U.-J.; Lee, J.-W.; Lee, S.-B.; Park, S.-J.; Song, R.-H.; Kim, G.; Lim, T.-H. Effect of cathode geometry on the electrochemical performance of flat tubular segmented-in-series (SIS) solid oxide fuel cell. Int. J. Hydrog. Energy 2015, 40, 6207–6215. [Google Scholar] [CrossRef]
  16. Sinha, A.A.; Sanjay; Ansari, M.Z.; Shukla, A.K.; Choudhary, T. Comprehensive review on integration strategies and numerical modeling of fuel cell hybrid system for power & heat production. Int. J. Hydrog. Energy 2023, 48, 33669–33704. [Google Scholar]
  17. Peng, J.; Zhao, D.; Xu, Y.; Wu, X.; Li, X. Comprehensive Analysis of Solid Oxide Fuel Cell Performance Degradation Mechanism, Prediction, and Optimization Studies. Energies 2023, 16, 788. [Google Scholar] [CrossRef]
  18. Zhao, Y.R.; Sadhukhan, J.; Lanzini, A.; Brandon, N.; Shah, N. Optimal integration strategies for a syngas fuelled SOFC and gas turbine hybrid. J. Power Sources 2011, 196, 9516–9527. [Google Scholar] [CrossRef]
  19. Huang, S.L.; Yang, C.; Chen, H.; Zhou, N.N.; Tucker, D. Coupling impacts of SOFC operating temperature and fuel utilization on system net efficiency in natural gas hybrid SOFC/GT system. Case Stud. Therm. Eng. 2022, 31, 12. [Google Scholar] [CrossRef]
  20. Ameri, M.; Mohammadi, R. Simulation of an atmospheric SOFC and gas turbine hybrid system using Aspen Plus software. Int. J. Energy Res. 2013, 37, 412–425. [Google Scholar] [CrossRef]
  21. Doherty, W.; Reynolds, A.; Kennedy, D. Computer simulation of a biomass gasification-solid oxide fuel cell power system using Aspen Plus. Energy 2010, 35, 4545–4555. [Google Scholar] [CrossRef]
  22. Ni, T.M.; Si, J.W.; Gong, X.H.; Zhang, K.; Pan, M.Z. Thermodynamic and economic analysis of a novel cascade waste heat recovery system for solid oxide fuel cell. Energy Convers. Manag. 2022, 259, 14. [Google Scholar] [CrossRef]
  23. Hussain, J.; Ali, R.; Akhtar, M.N.; Jaffery, M.H.; Shakir, I.; Raza, R. Modeling and simulation of planar SOFC to study the electrochemical properties. Curr. Appl. Phys. 2020, 20, 660–672. [Google Scholar] [CrossRef]
  24. Qiao, Y.X.; Fan, C.; Sun, K.N. Multiphysics Simulation of a Solid Oxide Fuel Cell Based on COMSOL Method. In Proceedings of the 5th International Conference on Advances in Energy, Environment and Chemical Science (AEECS), Shanghai, China, 26–28 February 2021; EDP Sciences: Les Ulis, France, 2021. [Google Scholar]
  25. Yakabe, H.; Ogiwara, T.; Hishinuma, M.; Yasuda, I. 3-D model calculation for planar SOFC. J. Power Sources 2001, 102, 144–154. [Google Scholar] [CrossRef]
  26. Khazaee, I.; Rava, A. Numerical simulation of the performance of solid oxide fuel cell with different flow channel geometries. Energy 2017, 119, 235–244. [Google Scholar] [CrossRef]
  27. Lee, S.; Park, M.; Kim, H.; Yoon, K.J.; Son, J.W.; Lee, J.H.; Kim, B.K.; Choi, W.; Hong, J. Thermal conditions and heat transfer characteristics of high-temperature solid oxide fuel cells investigated by three-dimensional numerical simulations. Energy 2017, 120, 293–305. [Google Scholar] [CrossRef]
  28. Milewski, J.; Świrski, K.; Santarelli, M.; Leone, P. Advanced Methods of Solid Oxide Fuel Cell Modeling; Springer Science & Business Media: Berlin, Germany, 2011. [Google Scholar]
  29. Patcharavorachot, Y.; Arpornwichanop, A.; Chuachuensuk, A. Electrochemical study of a planar solid oxide fuel cell: Role of support structures. J. Power Sources 2008, 177, 254–261. [Google Scholar] [CrossRef]
  30. Chatrattanawet, N.; Saebea, D.; Authayanun, S.; Arpornwichanop, A.; Patcharavorachot, Y. Performance and environmental study of a biogas-fuelled solid oxide fuel cell with different reforming approaches. Energy 2018, 146, 131–140. [Google Scholar] [CrossRef]
  31. Janardhanan, V.M.; Deutschmann, O. CFD analysis of a solid oxide fuel cell with internal reforming: Coupled interactions of transport, heterogeneous catalysis and electrochemical processes. J. Power Sources 2006, 162, 1192–1202. [Google Scholar] [CrossRef]
  32. Rosner, F.; Rao, A.; Samuelsen, S. Economics of cell design and thermal management in solid oxide fuel cells under SOFC-GT hybrid operating conditions. Energy Convers. Manag. 2020, 220, 112952. [Google Scholar] [CrossRef]
  33. Oryshchyn, D.; Harun, N.F.; Tucker, D.; Bryden, K.M.; Shadle, L. Fuel utilization effects on system efficiency in solid oxide fuel cell gas turbine hybrid systems. Appl. Energy 2018, 228, 1953–1965. [Google Scholar] [CrossRef]
  34. Rosner, F.; Rao, A.; Samuelsen, S. Thermo-economic analyses of solid oxide fuel cell-gas turbine hybrids considering thermal cell gradients. J. Power Sources 2021, 507, 230271. [Google Scholar] [CrossRef]
  35. Sadeghi, M.; Jafari, M.; Hajimolana, Y.S.; Woudstra, T.; Aravind, P.V. Size and exergy assessment of solid oxide fuel cell-based H2-fed power generation system with alternative electrolytes: A comparative study. Energy Convers. Manag. 2021, 228, 113681. [Google Scholar] [CrossRef]
  36. Wu, Z.; Zhang, Z.; Ni, M. Management, Modeling of a novel SOFC-PEMFC hybrid system coupled with thermal swing adsorption for H-2 purification: Parametric and exergy analyses. Energy Convers. Manag. 2018, 174, 802–813. [Google Scholar] [CrossRef]
  37. Akkaya, A.V.; Sahin, B.; Erdem, H.H. Exergetic performance coefficient analysis of a simple fuel cell system. Int. J. Hydrog. Energy 2007, 32, 4600–4609. [Google Scholar] [CrossRef]
  38. Ilbas, M.; Kumuk, B. Numerical modelling of a cathode-supported solid oxide fuel cell (SOFC) in comparison with an electrolyte-supported model. J. Energy Inst. 2019, 92, 682–692. [Google Scholar] [CrossRef]
  39. Ferguson, J.R.; Fiard, J.M.; Herbin, R. Three-dimensional numerical simulation for various geometries of solid oxide fuel cells. J. Power Sources 1996, 58, 109–122. [Google Scholar] [CrossRef]
  40. Celik, A.N. Three-dimensional multiphysics model of a planar solid oxide fuel cell using computational fluid dynamics approach. Int. J. Hydrog. Energy 2018, 43, 19730–19748. [Google Scholar] [CrossRef]
  41. Hajimolana, S.A.; Hussain, M.A.; Daud, W.; Soroush, M.; Shamiri, A. Mathematical modeling of solid oxide fuel cells: A review. Renew. Sust. Energ. Rev. 2011, 15, 1893–1917. [Google Scholar] [CrossRef]
  42. Akhtar, N.; Decent, S.P.; Loghin, D.; Kendall, K. A three-dimensional numerical model of a single-chamber solid oxide fuel cell. Int. J. Hydrog. Energy 2009, 34, 8645–8663. [Google Scholar] [CrossRef]
  43. Andersson, M.; Yuan, J.L.; Sundén, B. SOFC modeling considering electrochemical reactions at the active three phase boundaries. Int. J. Heat Mass Transf. 2012, 55, 773–788. [Google Scholar] [CrossRef]
  44. Kanno, D.; Shikazono, N.; Takagi, N.; Matsuzaki, K.; Kasagi, N. Evaluation of SOFC anode polarization simulation using three-dimensional microstructures reconstructed by FIB tomography. Electrochim. Acta 2011, 56, 4015–4021. [Google Scholar] [CrossRef]
  45. Wang, X.S.; Lv, X.J.; Weng, Y.W. Performance analysis of a biogas-fueled SOFC/GT hybrid system integrated with anode-combustor exhaust gas recirculation loops. Energy 2020, 197, 17. [Google Scholar] [CrossRef]
  46. Campanari, S. Thermodynamic model and parametric analysis of a tubular SOFC module. J. Power Sources 2001, 92, 26–34. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of SOFC-integrated system.
Figure 1. Schematic diagram of SOFC-integrated system.
Sustainability 16 01643 g001
Figure 2. (a) Specific structure of SOFC; (b) dimensions and boundary conditions. The blue arrow at the top indicates the air inlet, and the blue arrow at the bottom indicates the fuel outlet.
Figure 2. (a) Specific structure of SOFC; (b) dimensions and boundary conditions. The blue arrow at the top indicates the air inlet, and the blue arrow at the bottom indicates the fuel outlet.
Sustainability 16 01643 g002
Figure 3. Comparison between experimental and simulation results.
Figure 3. Comparison between experimental and simulation results.
Sustainability 16 01643 g003
Figure 4. Cell voltage–current-power (V-I-P) for all lengths.
Figure 4. Cell voltage–current-power (V-I-P) for all lengths.
Sustainability 16 01643 g004
Figure 5. Effect on flow channel length and cell voltage on power of the SOFC stack.
Figure 5. Effect on flow channel length and cell voltage on power of the SOFC stack.
Sustainability 16 01643 g005
Figure 6. Effect on flow channel length and cell voltage on fuel utilization factor.
Figure 6. Effect on flow channel length and cell voltage on fuel utilization factor.
Sustainability 16 01643 g006
Figure 7. (a) Effect of flow channel length and cell voltage on the pressure ratio of the air compressor; (b) effect of flow channel length and cell voltage on the power of the air compressor.
Figure 7. (a) Effect of flow channel length and cell voltage on the pressure ratio of the air compressor; (b) effect of flow channel length and cell voltage on the power of the air compressor.
Sustainability 16 01643 g007
Figure 8. (a) Effect of flow channel length and cell voltage on the pressure ratio of the fuel compressor; (b) effect of flow channel length and cell voltage on the power of the fuel compressor.
Figure 8. (a) Effect of flow channel length and cell voltage on the pressure ratio of the fuel compressor; (b) effect of flow channel length and cell voltage on the power of the fuel compressor.
Sustainability 16 01643 g008
Figure 9. (a) Effect of flow channel length and cell voltage on the power of the air heat exchanger; (b) Effect of flow channel length and cell voltage on the power of the water heat exchanger; (c) effect of flow channel length and cell voltage on the power of the hydrogen heat exchanger.
Figure 9. (a) Effect of flow channel length and cell voltage on the power of the air heat exchanger; (b) Effect of flow channel length and cell voltage on the power of the water heat exchanger; (c) effect of flow channel length and cell voltage on the power of the hydrogen heat exchanger.
Sustainability 16 01643 g009
Figure 10. (a) Air compressor outlet gas temperature at different flow channel lengths; (b) fuel compressor outlet gas temperature at different flow channel lengths.
Figure 10. (a) Air compressor outlet gas temperature at different flow channel lengths; (b) fuel compressor outlet gas temperature at different flow channel lengths.
Sustainability 16 01643 g010
Figure 11. (a) Effect of flow channel length and cell voltage on BOP power consumption; (b) effect of flow channel length and cell voltage on net electrical power.
Figure 11. (a) Effect of flow channel length and cell voltage on BOP power consumption; (b) effect of flow channel length and cell voltage on net electrical power.
Sustainability 16 01643 g011
Figure 12. Net electrical power.
Figure 12. Net electrical power.
Sustainability 16 01643 g012
Figure 13. Effect of flow channel length and cell voltage on net electrical efficiency.
Figure 13. Effect of flow channel length and cell voltage on net electrical efficiency.
Sustainability 16 01643 g013
Figure 14. Net electrical efficiency.
Figure 14. Net electrical efficiency.
Sustainability 16 01643 g014
Figure 15. Effect of flow channel length and cell voltage on thermoelectric efficiency.
Figure 15. Effect of flow channel length and cell voltage on thermoelectric efficiency.
Sustainability 16 01643 g015
Table 1. SOFC model’s geometric dimensions.
Table 1. SOFC model’s geometric dimensions.
ParameterValueUnitSymbol
Channel width2mmwch
Rib width2mmwrib
Anode thickness0.15mmtan
Electrolyte thickness0.1mmtel
Cathode thickness0.1mmtca
Gas channel height2mmhch
Flow channel length60–180mmLcell
Number of cells1000-Ncell
Table 2. SOFC boundary conditions.
Table 2. SOFC boundary conditions.
ParameterValueUnit
Anode electric potential0V
Cathode electric potential0.4–0.9V
Anode air velocity0.8m/s
Cathode air velocity3m/s
Anode mass fractionH2:H2O = 0.4:0.6-
Cathode mass fractionO2:N2 = 0.15:0.85-
Anode fuel outlet0Pa
Cathode fuel outlet0Pa
Cell operating temperature800
Cell operating pressure101.32kPa
Table 3. SOFC parameters.
Table 3. SOFC parameters.
ParameterValueUnitSymbolReferences
Anode electronic conductivity2149.2S/m σ s , a [39]
Cathode electronic conductivity5093S/m σ s , c [39]
Electrolyte ionic conductivity2.2669S/m σ i o n , l [39]
Anodic transfer coefficient of anode0.5 α a a [40]
Anodic transfer coefficient of cathode3.5 α c a [40]
Anode activation energy6.54 × 10111/(Ω∙m2) K a n [26,41]
Cathode activation energy2.35 × 10111/(Ω∙m2) K c a [26,41]
Exchange current density of anode4637.4A/m2 i 0 , a n [26]
Exchange current density of cathode1166.2A/m2 i 0 , c a [26]
The specific surface area of anode102,5001/m S a n [42]
The specific surface area of cathode102,5001/m S c a [42]
Electrolyte volume fraction0.7- θ [43]
porosity0.4- ε [26,44]
Table 4. Substrate parameters.
Table 4. Substrate parameters.
ParameterValueUnitReferences
Porosity of anode0.4-[26,44]
Porosity of cathode0.4-[26,44]
Permeability of anode1.76 × 10−11m2[26]
Permeability of cathode1.76 × 10−11m2[26]
Table 5. Material properties.
Table 5. Material properties.
ParameterValueUnitReferences
Anode heat capacity450J/(kg∙k)[43]
Cathode heat capacity430J/(kg∙k)[43]
Electrolyte heat capacity470J/(kg∙k)[43]
Anode density3310kg/m2[45]
Cathode density3030kg/m2[45]
Electrolyte density5160kg/m2[45]
Thermal conductivity of the anode11W/(m∙k)[43]
Thermal conductivity of the cathode6W/(m∙k)[43]
Thermal conductivity of electrolyte2.7W/(m∙k)[43]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, Y.; Liu, J.; Fu, L.; Wang, Q. Numerical Study on Effects of Flow Channel Length on Solid Oxide Fuel Cell-Integrated System Performances. Sustainability 2024, 16, 1643. https://doi.org/10.3390/su16041643

AMA Style

Liu Y, Liu J, Fu L, Wang Q. Numerical Study on Effects of Flow Channel Length on Solid Oxide Fuel Cell-Integrated System Performances. Sustainability. 2024; 16(4):1643. https://doi.org/10.3390/su16041643

Chicago/Turabian Style

Liu, Yuhang, Jinyi Liu, Lirong Fu, and Qiao Wang. 2024. "Numerical Study on Effects of Flow Channel Length on Solid Oxide Fuel Cell-Integrated System Performances" Sustainability 16, no. 4: 1643. https://doi.org/10.3390/su16041643

APA Style

Liu, Y., Liu, J., Fu, L., & Wang, Q. (2024). Numerical Study on Effects of Flow Channel Length on Solid Oxide Fuel Cell-Integrated System Performances. Sustainability, 16(4), 1643. https://doi.org/10.3390/su16041643

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop