[ Article ]
Transactions of the Korean Hydrogen and New Energy Society - Vol. 32, No. 6, pp.470-476
ISSN: 1738-7264 (Print) 2288-7407 (Online)
Print publication date 30 Dec 2021
Received 01 Dec 2021 Revised 16 Dec 2021 Accepted 20 Dec 2021

# Parametric Study on High Power SOEC System

TUANANH BUI1, 2 ; YOUNG SANG KIM1, 2, ; VAN-TIEN GIAP1 ; DONG KEUN LEE1 ; KOOK YOUNG AHN1, 2
1Department of Clean Fuel and Power Generation, Korea Institute of Machinery & Materials (KIMM), 156 Gajeongbuk-ro, Yuseong-gu, Daejeon 34103, Korea
2Department of Environment and Energy Mechanical Engineering, University of Science and Technology (UST), 156 Gajeongbuk-ro, Yuseong-gu, Daejeon 34113, Korea
고출력 SOEC 시스템의 매개변수 연구
뚜안앵1, 2 ; 김영상1, 2, ; 잡반티엔1 ; 이동근1 ; 안국영1, 2
1한국기계연구원 청정연료발전연구실
2과학기술연합대학원대학교 환경에너지기계공학과

Correspondence to: yskim@kimm.re.kr

## Abstract

In the near future, with the urgent requirement of environmental protection, hydrogen based energy system is essential. However, at the present time, most of the hydrogen is produced by reforming, which still produces carbon dioxide. This study proposes a high-power electrolytic hydrogen production system based on solid oxide electrolysis cell with no harmful emissions to the environment. Besides that, the parametric study and optimization are also carried to examine the effect of individual parameter and their combination on system efficiency. The result shows that the increase in steam conversion rate and hydrogen molar fraction in incoming stream reduces system efficiency because of the fuel heater power increase. Besides, the higher Faraday efficiency does not always result a higher system efficiency.

## Keywords:

Solid oxide electrolysis cell, Parametric study, Faraday efficiency, Steam conversion rate

## 키워드:

고체 산화물 전해전지, 매개변수 연구, 파라데이 효율, 스팀 전환율

## 1. Introduction

Currently, the pressure on fossil fuel reduction is increasing year by year due to global warming and climate change. This scenario opens a very promising future for hydrogen production. Not only playing as a clean fuel, hydrogen is also an alternative energy storage chemical in which electricity can be converted into hydrogen to be stored. Most of the current hydrogen produced is by reforming technologies because these technologies are mature at the time and low production cost, between 1-2 dollars1). However, producing hydrogen by reforming also generates carbon dioxide as by product while the future hydrogen society definitely requires a less polluted or even zero-carbon emitted technology. Electrolysis hydrogen production has recently been the center of focus because it can effectively utilize the excess electricity generated from renewable sources such as wind energy, solar energy, etc. The most challenging barrier of electrolysis hydrogen production is its high cost, normally from 3 to 6 $/kg H2 in case of polymer exchange membrane electrolysis cell2), or 2.8 to 5.8$/kg H2 in case of solid oxide electrolysis cell (SOEC)3,4), a more efficient cell. In order to reduce production cost, there are several ways. Among them, increasing stack power to produce more hydrogen with the same capital cost attracting attention from researchers1,5). This study introduces a high power SOEC system, then parametric study will be used to find optimal point.

## 2. High power SOEC system and its modeling by simulation

### 2.1 Description of analyzed system

Fig. 1 shows a flow diagram of the analyzed high power SOEC system. Liquid water is fed to pre-vaporizer where it is heated up by the high temperature flow out from fuel heat exchanger (F-HEX). After pre-vaporizer, hot water comes to vaporizer, receiving heat from heat source such as nuclear reactor, concentrated solar panel, etc. vaporizes into steam. This steam recuperates the heat from extremely hot flue gas stream out from SOEC at the F-HEX, before going through an electrical heater (F-heater) which increases the steam temperature to SOEC operating temperature. In the remaining side of SOEC, air is blown to an air heat exchanger (A-HEX) where it receives heat from the outgoing air stream from SOEC, then once again the air stream is heated up by an electrical heater (A-heater) before being supplied to SOEC. Through the electrochemical reaction of SOEC, the supplied electricity and steam are converted into oxygen and hydrogen. Then, mixes with the unconverted steam to be H2-rich hot stream, exiting SOEC to F-HEX. After transferring heat to cold steam at F-HEX and pre-vaporizer, H2-rich stream is split into two streams; one stream is recirculated to vaporizer to recover steam and also keep H2 molar fraction in inlet steam higher than certain value6). The remaining flows through condenser where steam condenses, and is removed from gaseous mixture which is finally separated to get pure hydrogen.

Flow diagram of the high power SOEC system

### 2.2 Assumption used in the system simulation

Table 1 summarizes the assumptions used in system modeling and simulation7-11). The current density is assumed to be 1 A/cm2 higher than almost normal SOEC5,12,13). Beside, operating condition of cell is chosen at 715˚C and 1 atm. Furthermore, the 36 kW stack includes 380 cells. Other essential coefficients such as inverter efficiency, Faraday efficiency, heat loss are set at widely used values: 92%, 92%, and 5%, respectively10). Lastly, the heat supplied to SOEC system is assumed from a 300˚C waste steam stream with mass flow rate of 15 kg/h.

Assumptions and parameters used in system simulations

### 2.3 Genetic algorithm optimization

The analyzed cycle was modeled using the EBSILON® (Steag, Zwingenberg, Germany) Professional commercial software package14). Fig. 2 shows a screenshot of the EBSILON software used for the analyzed cycle. During optimization process genetic algorithm is employed to find optimal point. Though genetic algorithm does not guarantee the global optimal point, this result suggests a parameter set for further investigation or analysis.

Screenshot of the EBSILON software used for the simulation

## 3. Performance analysis

To examine the performance of SOEC system, parametric analysis method is conducted with external steam temperature, heat exchanger effectiveness, steam conversion rate, hydrogen molar fraction, Faraday efficiency, and cell voltage. Table 2 shows the range of each variable in parametric study and optimization6,9,11,15). In each examination, only studied parameter varies, others are kept at design condition. On the other hand, in the final step all listed parameters with their ranges are brought into an optimization to find the optimal condition. In both parametric studies and optimization, system efficiency is measured for comparison. Eq. (1) describes how system efficiency is calculated.

 ${\eta }_{sys}=\frac{{\stackrel{˙}{m}}_{{H}_{2}out}×{LHV}_{{H}_{2}}}{{E}_{in}}$ (1)
• where:
• ηsys : system efficiency
• ${\stackrel{˙}{m}}_{{H}_{2}out}$ : mass flow rate of produced hydrogen
• LHVH2 : low heating value of hydrogen at 25˚C

Variable range in parametric study and optimization

## 4. Result and discussion

### 4.1 Effect of external steam temperature and heat exchanger effectiveness

The figures show the system efficiency as function of each parameter in parametric study. In Fig. 3, external steam temperature is proportional to the amount of heat fed to the system, so the system efficiency increases when steam temperature increases, as expected. In the same way, the HEX effectiveness also affects positively to the system efficiency because it performs the heat recuperation of heat exchanger as shown in Fig. 4. Moreover, the dependence of system efficiency to HEX effectiveness is very high; about 10 percent points increases in system efficiency when HEX effectiveness changes from 0.7 to 0.9.

System efficiency depends on the external steam temperature

System efficiency depends on the HEX effectiveness

### 4.2 Effect of steam conversion rate and hydrogen molar fraction

Fig. 5 shows the system dependence on steam conversion. When steam conversion rate increases, the H2 content in the output of SOEC stack is also higher, leading to the lower recirculation flow, and then less recovered heat from fuel off-gas. As a result, the fuel heater power increases, causing slightly decrease in system efficiency. In Fig. 6, the effect of hydrogen molar fraction is quite similar to that of steam conversion rate. A higher hydrogen molar fraction in stack inlet stream requires a smaller recirculation flow, leading to less recuperated heat amount. Consequently, the F-heater has to provide more heat to achieve deign temperature, reducing total system efficiency. However, the amplitude is bigger than that in Fig. 5 due to the higher recycle blower power consumption.

System efficiency depends on the steam conversion

System efficiency depends on the H2 molar fraction

### 4.3 Effect of Faraday efficiency and cell voltage

Fig. 7 shows the dependency of system efficiency on Faraday efficiency. The SOEC efficiency increases with Faraday efficiency. However, it affects the heat release from or absorption into SOEC stack, consequently the whole system efficiency does not always increase with Faraday efficiency’s increase. The maximum system efficiency is 89.7% when Faraday efficiency is 0.96. In Fig. 8, system efficiency decreases when cell voltage increases because SOEC works at low faraday efficiency (0.92) that is in exothermic condition. Therefore, the higher voltage makes the higher loss.

System efficiency depends on faraday efficiency

System efficiency depends on the cell voltage

### 4.4 Optimization

Because Faraday efficiency effect is not monotonous, it is set at discrete values from 0.92 to 0.97 in optimization. Besides, the cell voltage also chosen at 1.285 V where the stack is almost at thermal balance point, neither absorb nor release heat9). In addition, the parametric study in the previous section also pointed out the system efficiency is maximum with cell voltage value of 1.285 V.

Table 3 shows the optimization results with Faraday efficiency from 0.92 to 0.97. In all the cases, the optimal results are found at the boundary value of the external steam temperature, 700˚C. However, the result also shows that, though the effect of all parameters are monotonous, the optimal point is not always achieved when the examined variables at the boundary. For example, when Faraday efficiency is 0.97, optimal point is achieved at HEX effectiveness of 0.787, steam conversion rate of 0.57 and hydrogen molar fraction of 0.200. This reveals the combination and correlation effect between the parameters.

Optimization result

Though genetic algorithm does not guarantee the global optimal point, this result suggests a parameter set for further investigation or analysis.

## 5. Conclusions

This study proposes a high power SOEC system and then, examined the effect each parameter including external steam temperature, heat exchanger effectiveness, steam conversion rate, hydrogen molar fraction, Faraday efficiency, and cell voltage to the system efficiency. Finally, the optimization based on genetic algorithm was conducted to see the combining effect of all listed parameters on system performance. In conclusion, there are several points can be summarized as following:

Higher steam conversion rate leads to high concentration of H2 in the recirculation flow, hence the fuel heater power increases. As a result, the system efficiency drops as fuel utilization increases.

Faraday efficiency increase can make SOEC changes from exothermic to endothermic condition, thus the system efficiency trend and the optimal points of the operation changed too.

The genetic optimization method showed the optimal point in which endothermic operation is better option.

The concentration of H2 in the feed fuel, steam conversion rate, and Faraday efficiency should be simultaneously considered in optimizing system operation.

## Acknowledgments

This work was supported by the Korea Institute of Energy Technology Evaluation and Planning (KETEP) and the Ministry of Trade, Industry & Energy (MOTIE) of the Republic of Korea (No. 20213030040110).

## References

• F. O. Ayodele, S. I. Mustapa, B. V. Ayodele, and N. Mohammad, “An overview of economic analysis and environmental impacts of natural gas conversion technologies”, Sustainability, Vol. 12, No. 23, 2020, pp. 10148. [https://doi.org/10.3390/su122310148]
• D. Peterson, D. A. DeSantis, and M. Hamdan, “DOE hydrogen and fuel cells program record 20004: cost of electrolytic hydrogen production with existing technology”, 2020.
• D. Peterson, E. L. Miller, A. Brisse, J. Hartvigsen, R. J. Petri, G. G. Tao, and S. Satyapal, “DOE Hydrogen and fuel cells program record 16014: hydrogen production cost from solid oxide electrolysis”, Hydrogen, 2016. Retrieved from https://www.hydrogen.energy.gov/pdfs/16014_h2_production_cost_solid_oxide_electrolysis.pdf, .
• V. T. Giap, Y. D. Lee, Y. S. Kim, and K. Y. Ahn, “Techno-economic analysis of reversible solid oxide fuel cell system couple with waste steam”, Trans Korean Hydrogen New Energy Soc, Vol. 30, No. 1, 2019, pp. 21-28.
• Y. Zheng, J. Wang, B. Yu, W. Zhang, J. Chen, J. Qiao, and J. Zhang, “A review of high temperature co-electrolysis of H2O and CO2 to produce sustainable fuels using solid oxide electrolysis cells (SOECs): advanced materials and technology”, Chemical Society Reviews, Vol. 46, No. 5, 2017, pp. 1427-1463. [https://doi.org/10.1039/C6CS00403B]
• Q. P. Fang, L. Blum, R. Peters, M. Peksen, P. Batfalsky, and D. Stolten, “SOFC stack performance under high fuel utilization”, International Journal of Hydrogen Energy, Vol. 40, No. 2, 2015, pp. 1128-1136. [https://doi.org/10.1016/j.ijhydene.2014.11.094]
• R. Elder, D. Cumming, and M. B. Mogensen, “Chapter 11-high temperature electrolysis”, Carbon Dioxide Utilisation, 2015, pp. 183-209. [https://doi.org/10.1016/B978-0-444-62746-9.00011-6]
• M. B. Mogensen, M. Chen, H. L. Frandsen, C. Graves, J. B. Hansen, K. V. Hansen, A. Hauch, T. Jacobsen, S. H. Jensen, T. L. Skafte, and X. Sun, “Reversible solid-oxide cells for clean and sustainable energy”, Clean Energy, Vol. 3, No. 3, 2019, pp. 175-201. [https://doi.org/10.1093/ce/zkz023]
• V. Saarinen, J. Pennanen, M. Kotisaari, O. Thomann, O. Himanen, S. Di Iorio, P. Hanoux, J. Aicart, K. Couturier, X. Sun, M. Chen, and B. R. Sudireddy, “Design, manufacturing, and operation of movable 2 × 10 kW size rSOC system”, Fuel Cells, Vol. 21, No. 5, 2021, pp. 477-487. [https://doi.org/10.1002/fuce.202100021]
• Y. S. Kim, Y. D. Lee, K. Y. Ahn, D. K. Lee, S. M. Lee, and E. J. Choi, “Operation characteristics according to steam temperature and effectivenss of external steam-related SOEC system”, Trans Korean Hydrogen New Energy Soc, Vol. 31, No. 6, 2020. pp. 596-604. [https://doi.org/10.7316/KHNES.2020.31.6.596]
• J. Schefold, A. Brisse, and H. Poepke, “23,000 h steam electrolysis with an electrolyte supported solid oxide cell”, International Journal of Hydrogen Energy, Vol. 42. No. 19, 2017, pp. 13415-13426. [https://doi.org/10.1016/j.ijhydene.2017.01.072]
• J. Schefold, A. Brisse, A. Surrey, and C. Walter, “80,000 current on/off cycles in a one year long steam electrolysis test with a solid oxide cell”, International Journal of Hydrogen Energy, Vol. 45, No. 8, 2020, pp. 5143-5154. [https://doi.org/10.1016/j.ijhydene.2019.05.124]
• G. Schiller, M. Lang, N. Monnerie, H. v. Storch, J. Reinhold, and P. Sundarraj, “Solar heat integrated solid oxide steam electrolysis for highly efficient hydrogen production”, Journal of Power Sources, Vol. 416, 2019, pp. 72-78. [https://doi.org/10.1016/j.jpowsour.2019.01.059]

### Fig. 1.

Flow diagram of the high power SOEC system

### Fig. 2.

Screenshot of the EBSILON software used for the simulation

### Fig. 3.

System efficiency depends on the external steam temperature

### Fig. 4.

System efficiency depends on the HEX effectiveness

### Fig. 5.

System efficiency depends on the steam conversion

### Fig. 6.

System efficiency depends on the H2 molar fraction

### Fig. 7.

System efficiency depends on faraday efficiency

### Fig. 8.

System efficiency depends on the cell voltage

### Table 1.

Assumptions and parameters used in system simulations

Parameter Value Unit
SOEC stack
Operating pressure 1 atm
Temperature 7157) °C
Current density 18) A/cm2
Cell voltage 1.2859) V
Cell number 350 cells
Heat loss 510) %
Input heat (external steam)
Temperature 30011) °C
Mass flow rate 15 kg/h

### Table 2.

Variable range in parametric study and optimization

Variable name Unit Range
HEX eff - 0.7-0.9
Steam temperature ˚C 200-70011)
Steam conversion rate - 0.5-0.815)
H2 molar fraction (xH2) - 0.2-0.56)
Cell voltage V 1.285-1.3059)

### Table 3.

Optimization result

effectiveness
External steam
temperature
Steam conversion
rate
H2 molar
fraction inlet
System
efficiency
0.92 1.285 0.9 700 0.58 0.280 87.167
0.93 0.9 700 0.68 0.251 88.833
0.94 0.9 700 0.58 0.207 91.382
0.95 0.9 700 0.75 0.285 91.437
0.96 0.9 700 0.58 0.201 92.133
0.97 0.787 700 0.57 0.200 92.434