Open Access
Issue
Renew. Energy Environ. Sustain.
Volume 6, 2021
Article Number 18
Number of page(s) 11
DOI https://doi.org/10.1051/rees/2021018
Published online 28 May 2021

© S. Horr et al., Published by EDP Sciences, 2021

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Nowadays, fuel cell technology is a fast growing scientific and technical field. It is in the process of constituting the core of the next industrial revolution. Recently, solid oxide fuel cells studying has attracted the attention of many researchers worldwide for clean energy production [1]. They are regarded as the most advantageous types appropriating renewable energy sources to produce electricity and provide great power density, which can be employed to versatile power generation systems [2]. SOFCs provide several advantages such as: high efficiency [3,4], low emission of pollutants, working in silence [1], fuel flexibility [5], long-term stability [6], use of low-cost catalyst material [7] and the ability of utilization for cogeneration [8,9]. Flow channel design has a vital importance on the performance of SOFC fuel cells. Consequently, various investigations have been carried out to explore the influence of flow channels design. Amongst these, the most often mentioned model given by Achenbach [10] studied with time-dependent the effects of flow manifolding with utilizing differential and finite equations that allow determining heat and mass transfer in SOFCs. The researchers found that the counter-flow design has an impact in improving the performance compared to cross and co-flow. Xu et al. [11] conducted a numerical model to study the planar solid oxide electrolyzer cell (SOEC) with diverse flow forms (co-flow, counter-flow and cross-flow fields). The authors arrived at that the cross-flow design presented the best performance as it produced hydrogen the most compared to the co-flow and counter-flow under similar boundary conditions. Hawkes et al. [12,13] presented 3D simulations on a SOEC stack under cross-flow configuration. They have discussed profiles of activation overpotential, temperature, operating potential, current density, Nernst potential, the gas composition of anode-side and cathode-side and hydrogen production at several deferent operating conditions of the stack using ANSYS FLUENT. Wuxi et al. [14] examined the impact of geometric parameters such as the dimension of the channel, the repeating cell height and the width of the manifold on the flow improvement. The simulation results showed that the ratio of the outlet to the inlet manifold width was affected significantly by the flow uniformity. Valery et al. [15] proposed a modified planar SOFC fuel cell with internal reforming. The study combined parallel flow and a new flow field design. The simulation analysis confirmed that the new modified anode configuration enhanced the performance cell and realized more uniform distribution of species concentration, velocity and current density. Saied et al. [16] compared diverse flow fields including helical, modified parallel design, traditional parallel, single-entry serpentine, double-entry serpentine and triple-entry serpentine. Their results noted that the triple-entry serpentine configuration achieved regular distributions for fuel and oxygen and presented better performance compared to other designs. Huang et al. [17] considered a new design using short guide vanes equally spaced around the feed header of the double-inlet/single-outlet flows of rib-channel interconnects to examine the flow uniformity in several interconnectors and its impact on cell performance. They established that the new design enhanced the degree of flow uniformity in interconnects and achieved a greater value of peak power density. Lin et al. [18] developed analytically the influence of rectangular ribs on the concentration polarization of planar SOFCs. The study showed that the best design was achieved by reducing the overall concentration and Ohmic polarization of the ribs. Qiuwan et al. [19] performed numerical studies examining the performance of SOFC with and without rectangular obstacles. This attempt indicated that the fuel cell with obstacles enhanced the cell performance as it attained better current density, species concentration and temperature distribution. Clearly, the flow field configuration has a vital effect on the performance of the fuel cell and is the key factor in SOFC design optimization. Therefore, further studies on new field designs are recommended to obtain more acute and reliable results on cell performance. The main objective of the present investigation is to enhance the performance of a solid oxide fuel cell (SOFC). This study contributed and help to choose the optimum configuration for the SOFC fuel cell system manufacturers and their applications and commercialization. Thus, a new flow field based on a sinusoidal flow has been proposed and studied.

The SOFC is an electrochemical device as other fuel cells; it is a combination of a reactant (oxygen) and a fuel to produce electricity and heat. The electrochemical reaction occurring in a planar SOFC involves the following individual reactions, at the anode and cathode [20], which could be expressed as follows: At anode H 2 + O 2 H 2 O + 2 e (1) At cathode  1 2 O 2 + 2 e O 2 (2) Overall Cell reaction H 2 + 1 2 O 2 H 2 O + Heat (3)

The cell voltage [2124] could be calculated from: V c e l l = E N e r n s t η O h m η a c t η c o n (4)where ηOhm, ηact and ηcon are the Ohmic, activation and concentration losses over-potential respectively [25].

The Nernst equation is also an open-circuit voltage (OCV); its value is related to the gas composition, the operating pressure, the operating temperature… etc. It is defined by equation (5) [2629]: E N e r n s t = Δ G 2 F + R T 2 F ln ( P H 2 . P O 2 0.5 P H 2 O ) (5)where T (K) is the operating temperature, P (Pa) the partial pressures of reacting species, R (J mol−1 K−1) is the gas constant, F (C mol−1) is Faraday number, and ΔG (J mol−1 K−1) is the Gibbs free energy.

2 Theoretical analysis

Mathematical models of a planar SOFC involve the governing equations of mass, momentum, energy and transport of charges. Moreover, knowledge of the Butler-Volmer equation was necessary for analyzing a model. In this study, a SOFC model was implemented in ANSYS FLUENT 18.1 software, for solving 3D and steady-state governing equations by using the finite volume approach. The applicable equations are:

  • Mass conservation equation [30]:

The steady-state mass conservation equation is given as: . ( ϵ ρ v ) = S m (6)

The source terms, which are applicable at the interface electrode-electrolyte, are obtained from equations (8)(10), calculated as [3134]: S m = S H 2 + S H 2 O + S O 2 (7) S H 2 = ( J a 2 F ) M H 2 (8) S H 2 O = + ( J a 2 F ) M H 2 O (9) S O 2 = ( J c 4 F ) M O 2 (10)where ρ (kg m−3) is the gas density, v (m s−1) the velocity, ε is the porosity and Sm (kg m−3 s−1) the mass source.

  • Momentum equation [30]:

Due to low Reynolds and the steady-state momentum conservation equation can be written as: . ( ϵ ρ v v ) = ϵ p + . [ ϵ μ ( v + ( v ) T ) ] + μ ϵ 2 k g v (11) where kg is the gas phase permeability (m2) and μ (kg m−1 s−1) the gas viscosity.

  • Species conservation [3034]:

The species transport equation inside SOFC is represented by equation (12). The source terms in species equation are earlier defined in equations (8)(10). . ( ρ y i j i n D e f f , i j x j + ρ v y i ) = S i (12)where yi is the mass fraction of species i, Deff (m2 s−1) the effective diffusivity coefficient between species i and j, Si (kg m−3 s−1) the source term for species i and xj the mole fraction for species j.

  • Energy equation [30]:

The energy conservation equation can be formulated as: . ( ϵ ρ c p v T ) = . ( k e f f T ) + S T (13)where keff (W m−1 K−1) is the coefficient of thermal conductivity, cp (J kg−1 K−1) is the specific heat capacity, and ST (W m−3) the heat source.

  • Charge transport [16,3135]:

In order to solve the charge conservation equation, electrons and ions transport should be considered. The electronic charge occurs in the electrodes and interconnectors; while, the ionic charge only exists in the electrodes and electrolyte. Thus, the governing equations for charge by Ohm's law are as follows:

Electronic charge balance: Anode electrode layer  . ( σ a e l ) = J a A V (14) Cathode electrode layer  . ( σ c e l ) = J c A V (15)

Ionic charge balance: Electrolyte layer  . ( σ e l e i 0 ) = 0 (16) Anode electrode layer  . ( σ a i 0 ) = J a A V (17) Cathode electrode layer  . ( σ c i 0 ) = J c A V (18)where σa and σc (s m−1) denote the electrical conductivity for anode and cathode, respectively. While Ja and Jc designating the volumetric current densities of anode and cathode; based on the Butler‑Volmer equation [30] could be computed from equations (19) and (20) as: J a = J 0 , a A V ( exp ( α 2 F η a c t , a R T ) exp ( ( 1 α ) 2 F η a c t , a R T ) ) (19) J c = J 0 , c A V ( exp ( β 4 F η a c t , c R T ) exp ( ( 1 β ) 4 F η a c t , c R T ) ) (20)where J0,a and J0,c (Am−2) are the reference exchange current densities of anode and cathode and AV (m2 m−3) is the reactive surface area per unit volume. The reference exchange current density for anode and cathode could be calculated as [36]: j 0 , a = k a ( P h 2 P r e f ) ( P h 2 o P r e f ) exp ( E a c t , a R T ) (21) j 0 , c = k c ( P O 2 P r e f ) 0.25 exp ( E a c t , c R T ) (22)ka and kc refer to the pre-exponential coefficients. Eact,a and Eact,c are the activation energies of the anode and cathode exchange current densities, respectively. P is the partial pressure of chemical species [36].

3 Methodology

The geometrical SOFC single cell was executed in the commercial software GAMBIT (Version 2.4.6). Figure 1 displays a schematic view of the SOFC with both sinusoidal and single flow field designs. It is constituted of the fuel-flow channel, anode gas diffusion layer (A-GDL), anode catalyst layer (A-CL), electrolyte, cathode catalyst layer (C-CL), cathode gas diffusion layer (C-GDL), air-flow channel, and anode and cathode collectors.

The flow of air and fuel in the investigated cell are counter-flow. Details of the cell dimensions are given in Table 1 [37]. The subsequent step after producing cell geometry was the creation of high-quality mesh. This was done by an ANSYS WORKBENCH MESH. In order to achieve the best quality, the mesh independency was studied by creating three different meshes with 76,800, 194,560 and 304,000 hexahedral cells. The results of mesh independency are presented in Table 2, from which it was noticed that mesh of 304,000 cells led to a reasonable agreement when compared to the numerical results of [37]. For more details about the mesh independency, refer to [34,38,39]. Once the computational mesh was created, it must be imported into the solver ANSYS FLUENT 18.1. After examining the mesh file, setting the parameters, materials and defining the boundary conditions for the two geometries; as listed in Tables 3 and 4, respectively.

thumbnail Fig. 1

Schematic view of SOFC with two different flow designs.

Table 1

Geometrical parameters.

Table 2

Mesh independency study.

Table 3

Specifications of the case study.

Table 4

Inlet boundary conditions for the anode and cathode channels.

4 Model validation

The predicted current density versus voltage (JV) results for our numerical model showed an acceptable accuracy with the data obtained from references [37,46] as shown in Figure 2. Nonetheless, the simulation manifested different values of current density at lower values of voltage. This is most probably due to errors and hypotheses associated with the SOFC model. Conversely, when using default parameters, the results were somewhat far from the other curves. This is why all conditions must be scrutinized and carefully readjusted according to the physical properties, to obtain results comparable to the experimental results.

thumbnail Fig. 2

JV curve comparison between present results and both the experimental and the numerical results obtained by [37,46].

5 Results and discussion

5.1 Distribution of velocity

Figure 3 shows the anode velocity distribution for both flow field designs at 0.6 V. The upmost velocity in anode channel was 10.63 and 10.35 m/s for sinusoidal and single flow field designs, respectively, as displayed in Figure 3a and b. It could be deduced that the distribution near GDL had improved, due to deflection design which increased the anode velocity and augmented the distribution mechanism. As a result, the velocity in the anode channel for the sinusoidal design was higher than the simple single-channel which led to cell performance enhancement.

thumbnail Fig. 3

Velocity distribution contour along the anode channel for two different flow designs at 0.6 V.

5.2 Distribution of hydrogen and water mass fraction

The hydrogen (H2) mass fraction contours on the anode channel for both designs are represented in Figure 4 at 0.6 V. The mass fraction reduction along the z-channel is due to consumption of hydrogen by the electrochemical reaction. It is also found that the mass fraction diminished from 0.97 to 0.38 and 0.43 for the sinusoidal and single flow design respectively. These results signify that the reaction rates at the sinusoidal flow design are greater than those detected in the single flow. The maximum use of hydrogen mass fraction consumed in the sinusoidal and single flow field designs are 60% and 55% respectively for the operating voltage of 0.6 V, as given in Table 5. According to Figure 4 along with Table 5, the sinusoidal form gave a greater level of hydrogen mass fraction. Consequently, this configuration is considered more uniform in diffusing of hydrogen mass fraction because of the best use of the active area in the CL (catalyst layer). On the other hand, increasing hydrogen mass fraction consumption increases the water (H2O) mass fraction produced as displayed in Figure 5. The sinusoidal flow design has a better distribution compared to the single flow field as demonstrated in Figure 5. The highest water mass fraction values produced for the sinusoidal and single flow designs were 42% and 34% respectively for the operating voltage 0.6 V, as shown in Table 6. This uniform distribution leads to regularly generate heat and decrease thermal stresses, suggesting that the sinusoidal cell aids in reducing losses.

thumbnail Fig. 4

Hydrogen mass fraction distribution contour along the anode channel for two different flow designs at 0.6 V.

Table 5

Maximum value consumed of hydrogen mass fraction at the anode/ electrolyte interface.

thumbnail Fig. 5

Water mass fraction distribution contour along the anode channel for two different flow designs at 0.6 V.

Table 6

Maximum value produced of Water mass fraction at the anode/ electrolyte interface.

5.3 Distribution of temperature

Figure 6 displays temperature contours for the sinusoidal and single flow field designs at 0.6 V. An increase in temperature values is clearly observable along the cells for both designs. The sinusoidal flow led to greater values as compared to the single flow, as can be seen in the same figure. This could be explicated by the higher rate of reaction. Specifically, the highest temperatures for the sinusoidal and single flow were 1128.30 K and 1126.95 K respectively for the operating voltage 0.6 V. For all the geometrical configurations, the highest temperature value was always located in the center of the cell which represents the electrolyte due to the counter-flow direction as marked in the literature [47,48].

thumbnail Fig. 6

Temperature distribution contour along the SOFCs cells for two different flow designs at 0.6 V.

5.4 Concentration, activation and ohmic overpotentials

The mathematical model was able to predict the overpotentials inside the SOFC for both sinusoidal and single cells, which is necessary for testing their performance. As presented in the literature, there are three overpotentials, caused by various mechanisms. The concentration overpotential is based on the resistance to transport the fraction of reactants and oxidant. As shown in Figure 7a the sinusoidal flow field led to a reduction in the concentration over potential, whereas, the concentration over potential increased with the current density for both designs. The activation overpotential was slightly higher for the single cell than the sinusoidal as shown in Figure 7b. This is due to increasing the temperature in the sinusoidal cell which contributed in reducing activation polarization. In addition, the sinusoidal cell took part in minimizing the Ohmic overpotential as displayed in Figure 7c. This is due to the improvement of the operation temperature which raises the conductivity of the electrolyte. From the above mentioned arguments, it could be deduced that the sinusoidal cell significantly contributed to the performance enhancement and to minimization of the concentration, activation and Ohmic overpotentials compared to the single cell.

thumbnail Fig. 7

(a) The concentration, (b) the activation, (c) the Ohmic overpotantials curves for two different flow field designs.

5.5 SOFC overall performance

Figure 8 displays the polarization (JV) and the power density curves for the two designs; from which it could be readily noticed that the sinusoidal flow design was more performable than the single flow. Additionally, the maximum power densities of the sinusoidal and single flow designs were 1.43 and 1.35 W/cm2, respectively. This result indicates that the sinusoidal flow design offers the most use of the hydrogen as considered in the previous section. Moreover, the difference in performance between the two designs is well remarked at low voltages due to decreasing in the concentration losses. As a consequence, the sinusoidal flow design has contributed significantly in ameliorating the performance.

thumbnail Fig. 8

The polarization curve (JV) and power density curve for two different flow designs.

6 Conclusions

In this paper, two different flow designs of SOFC have been numerically developed in the software ANSYS FLUENT 18.1. The main motivation behind this research was to ameliorate the performance of SOFC fuel cells. The results obtained from this model were confronted to available research data in the literature. Simulation results including velocity distributions, gas species, temperature, concentration overpotantial, activation overpotantial, Ohmic overpotantial and cell performance for both flow designs have been introduced and discussed. The sinusoidal design gave a better distribution of velocities and temperatures than the single flow, which improved the transport of hydrogen of the cell. As a result, more water was produced at the anode channel and accordingly, the cell performance rose. On the other hand, it was established that the channel design performed an important role in enhancing the performance and decreasing the overpotantials losses of SOFC fuel cells. The CFD results indicated that, when the cell was run at high voltage, the influence of channel design on the cell performance was negligible; whereas, when treated at the low voltage the channel design influence became significant. Lastly, it could be concluded that the sinusoidal design has greater performance than the single flow; and this could be useful in choosing the best configuration for certain applications requiring the greatest achievable performance. Future work will focus on further understanding of the porosity effects of SOFC fuel cells.

Nomenclatures

AV: Reactive surface area per unit volume, m2 m−3

cp: Specific heat capacity, J kg−1 K−1

D: Gas diffusivity, m2 s−1

Eact: Activation energy, J mol−1

ENernest: Nernst potential, V

F: Faraday's constant, C mol−1

J: Current density, A m−2

j: Exchange current density, A m−2

ΔG: Gibbs free energy change, J mol−1 K−1

k: Pre-exponential factor

kg: Gas phase permeability, m2

p: Gas pressure, Pa

R: Universal gas constant, J mol−1 K−1

s: Source term, kg m−3 s−1, W m−3

T: Temperature, K

Vcell: Cell potential, V

x: Mole fraction

y: Mass fraction

Greek letters

0: Standard state

Α: Transfer coefficient in anode

Β: Transfer coefficient in cathode

ϵ: Porosity

∅: Exchange potential (V)

η: Overpotential, V

k: Coefficient of thermal conductivity, W m−1 K−1

σ: Specific conductivity, s m−1

μ: Fluid viscosity, kg m−1 s−1

v: Velocity, m s−1

ρ: Density, kg m−3

Subscripts and superscripts

a: Anode

act: Activation

c: Cathode

conc: Concetration

eff: Effective

el: Electronic

ele: Electrolyte

j: Species j

i: Species i

i0: Ionic

LSM: Cathode material

Ni-YSZ: Anode material

OCV: Open circuit voltage

Ohm: Ohmic

SOFC: Solid oxide fuel cell

TPB: Three phase boundary length

YSZ: Electrolyte material

References

  1. U.S. DOE, National Energy Technology Laboratory, in Fuel Cell Handbook (U.S. DOE, National Energy Technology Laboratory, 2004), 7th ed. [Google Scholar]
  2. A.V. Virkar, J.W. Kim, K. Mehta, K.Z. Fung, Low temperature, high performance, planar solid oxide fuel cells and stacks (1997) [Google Scholar]
  3. C. Sun, U. Stimming, Recent anode advances in solid oxide fuel cells, J. Power Sources 171, 247‑260 (2007) [CrossRef] [Google Scholar]
  4. S. Revankar, P. Majumdar, in Fuel Cells Principles, Design and Analysis (Taylor & Francis Group. LLC, 2014), 1st edn. [Google Scholar]
  5. J. Wu, X. Liu, Recent development of SOFC metallic interconnect, J. Mater. Sci. Technol. 26, 293‑305 (2010) [CrossRef] [Google Scholar]
  6. S. (Supramaniam) Srinivasan, Fuel Cells: From Fundamentals to Applications, Springer (2006) [Google Scholar]
  7. K.H.J. Buschow, Encyclopedia of Materials: Science and Technology, Elsevier (2001) [Google Scholar]
  8. M. Xu, T.S. Li, M. Yang, M. Andersson, I. Fransson, T. Larsson et al., Modeling of an anode supported solid oxide fuel cell focusing on thermal stresses, Int. J. Hydrogen Energy 1, 1–14 (2016) [Google Scholar]
  9. B. Ghorbani, K. Vijayaraghavan, 3D and simplified pseudo-2D modeling of single cell of a high temperature solid oxide fuel cell to be used for online control strategies, Int. J. Hydrogen Energy (2018) [Google Scholar]
  10. E. Achenbach, Three-dimensional and time-dependent simulation of a planar solid oxide fuel cell stack, J. Power Sources 49, 333‑348 (1994) [CrossRef] [Google Scholar]
  11. Z. Xu, X. Zhang, G. Li, G. Xiao, J. Wang, Comparative performance investigation of different gas flow configurations for a planar solid oxide electrolyzer cell, Int. J. Hydrogen Energy (2017) [Google Scholar]
  12. G.L. Hawkes, J. O'Brien, C.M. Stoots, B. Hawkes, 3D CFD model of a multi-cell high-temperature electrolysis stack, Int. J. Hydrogen Energy 34, 4189‑4197 (2009) [Google Scholar]
  13. G.L. Hawkes, J.E. O'Brien, C.M. Stoots, J.S. Herring, M. Shahnam, CFD model of a planar solid oxide electrolysis cell for hydrogen production from nuclear energy, Nucl. Technol. 158, 132‑144 (2007) [CrossRef] [Google Scholar]
  14. W. Bi, D. Chen, Z. Lin, A key geometric parameter for the flow uniformity in planar solid oxide fuel cell stacks, Int. J. Hydrogen Energy 34, 3873‑3884 (2009) [Google Scholar]
  15. V. Danilov, M. Tade, A CFD-based model of a planar SOFC for anode flow field design, Int. J. Hydrogen Energy 34, 8998‑9006 (2009) [Google Scholar]
  16. M. Saied, K. Ahmed, M. Nemat-Alla, M. Ahmed, M. El-Sebaie, Performance study of solid oxide fuel cell with various flow field designs: numerical study, Int. J. Hydrogen Energy (2018) [Google Scholar]
  17. C. Huang, S. Shy, C. Lee, On flow uniformity in arious interconnects and its influence to cell performance of planar SOFC, J. Power Sources 183, 205–213 (2008) [CrossRef] [Google Scholar]
  18. Z. Lin, J.W. Stevenson, M.A. Khaleel, The effect of interconnect rib size on the fuel cell concentration polarization in planar SOFCs, J. Power Sources 117, 92‑97 (2003) [CrossRef] [Google Scholar]
  19. Q. Shen, L. Sun, B. Wang, Numerical simulation of the effects of obstacles in gas flow fields of a solid oxide fuel cell, Int. J. Electrochem. Sci. 14, 1698–1712 (2019) [CrossRef] [Google Scholar]
  20. S.C. Singhal, K. Kendall, High temperature Solid Oxide Fuel Cells: Fundamentals, Design, and Applications, Elsevier, New York (2003) [Google Scholar]
  21. S. Cordiner, M. Feola, V. Mulone, F. Romanelli, Analysis of a SOFC energy generation system fuelled with biomass reformate, Appl. Thermal Eng. 27, 738–747 (2007) [CrossRef] [Google Scholar]
  22. S. Campanari, P. Iora, Definition and sensitivity analysis of a finite volume SOFC model for a tubular cell geometry, J. Power Sources 132, 113‑126 (2004) [CrossRef] [Google Scholar]
  23. A. Luca, Experimental and numerical analysis of a radial flow solid oxide fuel cell, Int. J. Hydrogen Energy 32, 4559–4574 (2007) [CrossRef] [Google Scholar]
  24. D. Sanchez, On the effect of methane internal reforming modelling in solid oxide fuel cells, Int. J. Hydrogen Energy 33, 1834‑1844 (2008) [CrossRef] [Google Scholar]
  25. H. Mahcene, H.B. Moussa, H. Bouguettaia, D. Bechki, S. Babay, M. Salah Meftah, Study of species, temperature distributions and the solid oxide fuel cells performance in a 2-D model, Int. J. Hydrogen Energy 36, 4244‑4252 (2011) [CrossRef] [Google Scholar]
  26. M. Andresson, J. Yuan, Sundén, SOFC modeling considering hydrogen and carbon monoxide as electrochemical reactants, J. Power Sources 232, 42‑54 (2013) [Google Scholar]
  27. S. Yang, Electrochemical analysis of an anode-supported SOFC, Int. J. Electrochem. Sci. 8, 2330–2344 (2013) [Google Scholar]
  28. R.P. O'Hayre, Fuel cells for electrochemical energy conversion, EPJ Web Confer. 148, 00013 (2017) [CrossRef] [Google Scholar]
  29. Y. Sahli, H.B. Moussa, B. Zitouni, Optimization study of the produced electric power by SOFCs, Int. J. Hydrogen Energy (2018) [Google Scholar]
  30. Y. Wang, R. Zhan, Y. Qin, G. Zhang, Q. Du, K. Jiao, Three-dimensional modeling of pressure effect on operating characteristics and performance of solid oxide fuel cell, Int. J. Hydrogen Energy (2018) [Google Scholar]
  31. J. Zhang, L. Lei, D. Liu, F. Zhao, F. Chen, H. Wang, Numerical investigation of solid oxide electrolysis cells for hydrogen production applied with different continuity expressions, Energy Convers. Manag. 149, 646–659 (2017) [CrossRef] [Google Scholar]
  32. M. Ni, Thermo-electrochemical modeling of ammonia-fueled solid oxide fuel cells considering ammonia thermal decomposition in the anode, Int. J. Hydrogen Energy 36, 3153–3166 (2011) [CrossRef] [Google Scholar]
  33. M. Ni, 2D thermal-fluid modeling and parametric analysis of a planar solid oxide fuel cell, Energy Convers. Manag. 51, 714–721 (2010) [CrossRef] [Google Scholar]
  34. M.E. Chelmehsara, J. Mahmoudimehr, Techno-economic comparison of anode-supported, cathode-supported, and electrolyte-supported SOFCs, Int. J. Hydrogen Energy (2018) [Google Scholar]
  35. H. Paradis, M. Andersson, B. Sundén, Modeling of mass and charge transport in a solid oxide fuel cell anode structure by a 3D lattice Boltzmann approach, Heat Mass Transfer (2015) [Google Scholar]
  36. A. Chaisantikulwat, C. Diaz-Goano, E.S. Meadows, Dynamic modelling and control of planar anode-supported solid oxide fuel cell, Comput. Chem. Eng. 32, 2365‑2381 (2008) [CrossRef] [Google Scholar]
  37. Y.M. Barzi, A. Raoufi, N.M. Rasi, S. Davari, Three dimensional simulation of a counter-flow planar solid oxide fuel cell, ECS Trans. 35, 1021‑1033 (2011) [CrossRef] [Google Scholar]
  38. M. Saied et al., Investigations of solid oxide fuel cells with functionally graded electrodes for high performance and safe thermal stress, Int. J. Hydrogen Energy 42, 15887‑15902 (2017) [CrossRef] [Google Scholar]
  39. A. Arvay et al., Convergence criteria establishment for 3D simulation of proton exchange membrane fuel cell, Int. J. Hydrogen Energy 37, 2482–2489 (2011) [CrossRef] [Google Scholar]
  40. M. Chyu, Numerical modeling of transport phenomena in solid oxide fuel cells, In 2005 Taiwan Solid Oxide Fuel Cell Symposium (2005) [Google Scholar]
  41. G. Anandakumar, N. Li, A. Verma, P. Singh, J.H. Kim, Thermal stress and probability of failure analyses of functionally graded solid oxide fuel cells, J. Power Sources 195, 6659‑6670 (2010) [CrossRef] [Google Scholar]
  42. P. Costamagna, K. Honegger, Modeling of solid oxide heat exchanger integrated stacks and simulation at high fuel utilization, J. Electrochem. Soc. 145, 3995–4007 (1998) [CrossRef] [Google Scholar]
  43. J.R. Ferguson, J.M. Fiard, R. Herbin, Three-dimensional numerical simulation for various geometries of solid oxide fuel cells, J. Power Sources 58, 109–122 (1996) [CrossRef] [Google Scholar]
  44. Z. Qu, P. Aravind, S. Boksteen, N. Dekker, A. Janssen, N. Woudstra et al., Three-dimensional computational fluid dynamics modeling of anode-supported planar SOFC, Int. J. Hydrogen Energy 36, 10209‑10220 (2011) [CrossRef] [Google Scholar]
  45. M.M. Hussain, X. Li, I. Dincer, Mathematical modeling of planar solid oxide fuel cells, J. Power Sources 161, 1012‑1022 (2006) [Google Scholar]
  46. K. Keegan, M. Khaleel et al., Analysis of a planar solid oxide fuel cell based automotive auxiliary power unit, SAE Technical Paper Series 2002–01-0413 (2002) [Google Scholar]
  47. B. Zitouni, H. Ben-Moussa et al., Temperature field, H2 and H2O mass transfer in SOFC single cell: electrode and electrolyte thickness effects, Int. J. Hydrogen Energy 34, 5032‑5092 (2009) [CrossRef] [Google Scholar]
  48. B. Zitouni, G.M. Andreadis et al., Two-dimensional numerical study of temperature field in an anode supported planar SOFC: effect of the chemical reaction, Int. J. Hydrogen Energy 36, 4228‑4235 (2001) [CrossRef] [Google Scholar]

Cite this article as: Sabrina Horr, Hocine Mohcene, Hamza Bouguettaia, Hocine Ben Moussa, Performance analysis of AS-SOFC fuel cell combining single and sinusoidal flow field: numerical study, Renew. Energy Environ. Sustain. 6, 18 (2021)

All Tables

Table 1

Geometrical parameters.

Table 2

Mesh independency study.

Table 3

Specifications of the case study.

Table 4

Inlet boundary conditions for the anode and cathode channels.

Table 5

Maximum value consumed of hydrogen mass fraction at the anode/ electrolyte interface.

Table 6

Maximum value produced of Water mass fraction at the anode/ electrolyte interface.

All Figures

thumbnail Fig. 1

Schematic view of SOFC with two different flow designs.

In the text
thumbnail Fig. 2

JV curve comparison between present results and both the experimental and the numerical results obtained by [37,46].

In the text
thumbnail Fig. 3

Velocity distribution contour along the anode channel for two different flow designs at 0.6 V.

In the text
thumbnail Fig. 4

Hydrogen mass fraction distribution contour along the anode channel for two different flow designs at 0.6 V.

In the text
thumbnail Fig. 5

Water mass fraction distribution contour along the anode channel for two different flow designs at 0.6 V.

In the text
thumbnail Fig. 6

Temperature distribution contour along the SOFCs cells for two different flow designs at 0.6 V.

In the text
thumbnail Fig. 7

(a) The concentration, (b) the activation, (c) the Ohmic overpotantials curves for two different flow field designs.

In the text
thumbnail Fig. 8

The polarization curve (JV) and power density curve for two different flow designs.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.