硅酸盐学报

相场法模拟固体氧化物燃料电池阳极内Ni-YSZ界面处YSZ相变引发界面形貌演化

PDF下载

程凯 , 焦震钧

硅酸盐学报 | 固体氧化物电池表界面专题—研究论文 2025,53(10): 2995-3002

收起

硅酸盐学报 |固体氧化物电池表界面专题—研究论文 2025 , 53 (10) : 2995 -3002

相场法模拟固体氧化物燃料电池阳极内Ni-YSZ界面处YSZ相变引发界面形貌演化

全屏

[Author(id=1322570447028339212, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, orderNo=0, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=24s058089@stu.hit.edu.cn, emailSecond=null, emailThird=null, correspondingAuthor=0, authorType=1, ext={EN=AuthorExt(id=1322570447133196814, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, authorId=1322570447028339212, language=EN, stringName=Kai CHENG, firstName=Kai, middleName=null, lastName=CHENG, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=School of Science, Harbin Institute of Technology (Shenzhen), Shenzhen 518055, Guangdong, China, bio={"content":"

CHENG Kai (2002-), male, Master candidate. E-mail: 24s058089@stu.hit.edu.cn

"}, bioImg=null, bioContent=

CHENG Kai (2002-), male, Master candidate. E-mail: 24s058089@stu.hit.edu.cn

, aboutCorrespAuthor=null), CN=AuthorExt(id=1322570447208694287, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, authorId=1322570447028339212, language=CN, stringName=程凯, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=哈尔滨工业大学（深圳）理学院，广东深圳 518055, bio={"content":"

程凯（2002—），男，硕士研究生。

"}, bioImg=null, bioContent=

程凯（2002—），男，硕士研究生。

, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1322570446927675912, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, xref=null, ext=[AuthorCompanyExt(id=1322570446944453129, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, companyId=1322570446927675912, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=School of Science, Harbin Institute of Technology (Shenzhen), Shenzhen 518055, Guangdong, China), AuthorCompanyExt(id=1322570446952841738, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, companyId=1322570446927675912, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=哈尔滨工业大学（深圳）理学院，广东深圳 518055)])]), Author(id=1322570447275803153, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, orderNo=1, firstName=null, middleName=null, lastName=null, nameCn=null, orcid=null, stid=null, country=null, authorPic=null, dead=0, email=jiaozhenjun@hit.edu.cn, emailSecond=null, emailThird=null, correspondingAuthor=1, authorType=1, ext={EN=AuthorExt(id=1322570447410020883, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, authorId=1322570447275803153, language=EN, stringName=Zhenjun JIAO, firstName=Zhenjun, middleName=null, lastName=JIAO, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=School of Science, Harbin Institute of Technology (Shenzhen), Shenzhen 518055, Guangdong, China, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null), CN=AuthorExt(id=1322570447493906964, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, authorId=1322570447275803153, language=CN, stringName=焦震钧, firstName=null, middleName=null, lastName=null, prefix=null, suffix=null, authorComment=null, nameInitials=null, affiliation=null, department=null, xref=null, address=哈尔滨工业大学（深圳）理学院，广东深圳 518055, bio=null, bioImg=null, bioContent=null, aboutCorrespAuthor=null)}, companyList=[AuthorCompany(id=1322570446927675912, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, xref=null, ext=[AuthorCompanyExt(id=1322570446944453129, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, companyId=1322570446927675912, language=EN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=School of Science, Harbin Institute of Technology (Shenzhen), Shenzhen 518055, Guangdong, China), AuthorCompanyExt(id=1322570446952841738, tenantId=1146029695717560320, journalId=1263187385517883426, articleId=1319726362705904017, companyId=1322570446927675912, language=CN, country=null, province=null, city=null, postcode=null, companyName=null, departmentName=null, remark=哈尔滨工业大学（深圳）理学院，广东深圳 518055)])])]

程凯, 焦震钧

作者信息

哈尔滨工业大学（深圳）理学院，广东深圳 518055

通讯作者:

焦震钧（1982—），男，博士，教授。

作者简介:

程凯（2002—），男，硕士研究生。

Simulation of Evolution of Interface Morphology Induced by Yttria-Stabilized Zirconia (YSZ) Phase Transformation at Ni-YSZ Interface in Solid Oxide Fuel Cells Anode Using Phase Field Method

Kai CHENG, Zhenjun JIAO

Affiliations

School of Science, Harbin Institute of Technology (Shenzhen), Shenzhen 518055, Guangdong, China

Corresponding Author:

JIAO Zhenjun (1982-), male, Ph.D., Professor. E-mail: jiaozhenjun@hit.edu.cn

CHENG Kai (2002-), male, Master candidate. E-mail: 24s058089@stu.hit.edu.cn

出版时间: 2025-10-20 doi: 10.14062/j.issn.0454-5648.20250194

文章导航

摘要

收起

固体氧化物燃料电池（SOFC）的电解质材料氧化钇稳定氧化锆（YSZ）在长期运行中会发生由立方晶向四方晶（c-t）相变所导致的材料离子电导率下降与力学性能退化，进而加剧全电池性能衰减。针对现有研究对YSZ相变动力学定量分析不足的问题，本工作通过构建微弹性相场模型，对SOFC阳极内部Ni-YSZ界面处的YSZ内部的c-t晶体相变过程进行了定量研究。模拟结果准确预测了实验过程中观察到的界面处YSZ形貌演化，并解释了Ni-YSZ界面附近YSZ内部孔洞的生成机制。研究通过理论建模与实验的相互验证，为提升SOFC阳极内部Ni-YSZ界面稳定性提供了理论支持。

关键词

相场法 / 固体氧化物燃料电池 / 镍-氧化钇稳定氧化锆界面 / 立方晶-四方晶相变 / 形貌演化

Abstract

收起

Introduction

Solid oxide fuel cells (SOFCs) are considered as a leading candidate for clean energy conversion due to their high efficiency and flexibility. Yttria-stabilized zirconia (YSZ) serves as an electrolyte material in SOFCs due to its high ionic conductivity and mechanical stability at a high operating temperature. However, the long-term stability of YSZ can be impaired via the cubic-tetragonal (c-t) phase transformation. This transformation leads to a decrease in ionic conductivity and mechanical integrity, which contributes to an overall degradation of SOFC performance. Despite its importance, most of the previous investigations primarily focused on electrochemical degradation mechanisms, while less attention was paid to the material stability of YSZ, especially at the Ni-YSZ interface. For instance, the atomic diffusion of Ni into YSZ during sintering can accelerate the phase transformation, leading to a significant decrease in ionic conductivity within 100 hours. Despite these findings, there is a lack of comprehensive quantitative analysis on the kinetics and morphological impacts of the c-t phase transformation at the Ni-YSZ interface. Moreover, the existing models often oversimplify an interplay between phase transformation and mechanical stress, restricting their ability to predict a long-term degradation behavior accurately. In this work, the phase field method was applied to simulate the c-t phase transformation in YSZ at the Ni-YSZ interface for elucidating its effect on the interface morphology and gaining insights into the SOFC degradation mechanisms.

Methods

A 3D microelastic phase field model was proposed to investigate the c-t phase transformation in YSZ. The evolution of both conserved and nonconserved order parameters were governed by the Allen-Cahn and Cahn-Hilliard equations, respectively. The total free energy was constructed under multiphysics coupling. The elastic constants for cubic YSZ and tetragonal YSZ were defined based on their respective crystal structures, with the values assigned to reflect their distinct mechanical properties. The misfit strain was determined using the lattice parameters for the cubic and tetragonal phases. The simulations were performed on a 64×64×64 voxel grid representing a single-crystal YSZ. The governing equations could be solved by the Fourier spectral method and a semi-implicit time-stepping scheme. The shape of the nucleus distribution was restricted to the YSZ surface to ensure a realism under SOFC operating conditions.

In this work, a dense YSZ disk with a diameter of 24 mm and a thickness of 0.5 mm was served as a substrate. After mechanical polishing and ultrasonic cleaning in acetone and ethanol, a stainless-steel shadow mask was employed to perform two perpendicular magnetron-sputtering depositions, producing a continuous Ni patterned film anode with a thickness of approximately 1 μm. For the cathode, a Pt paste was screen-printed onto the opposite side of the YSZ substrate. The discharge tests of the SOFC were carried out in a custom-designed alumina reaction chamber. Humidified 97% H₂ was supplied to the anode and pure O₂ to the cathode, both at a fixed flow rate of 50 mL_·min^-1. The cell was polarized at 0.7 V for 100 h. After testing, the Ni film was mechanically peeled off, regions where Ni detached were selected and the underlying YSZ surface morphology was examined by scanning electron microscopy (SEM).

Results and discussion

The phase field simulations reveal that tetragonal variants preferentially nucleate at near the Ni-YSZ interface. High stresses appear at around the Ni-YSZ interfaces with the development of the t-variants. At certain iterations, high local stresses of exceeding 2 GPa can occur, leading to a radial growth of the variants into the bulk. A significant morphological evolution occurs, forming ridge structures along the YSZ surface and submicron pores in the interior at near the Ni-YSZ interface. The agreement between simulation data and experimental results confirms that the c-t phase transformation causes these morphological evolutions. The formation of pores is particularly detrimental as it disrupts the ion conduction pathways and weakens mechanical stability of YSZ, which can be correlated to the degradation of the SOFC performance.

Conclusions

In this work, a phase field modeling was used to simulate the c-t phase transformation in YSZ at the Ni-YSZ interface in SOFC anodes. The simulations accurately reproduced the morphological evolutions during long-term operation of SOFC, i.e., the formation of surface ridges and internal submicron pores, providing important insights into the SOFC degradation mechanisms. The preferential nucleation of tetragonal variants at the surface and stress-driven void formation could highlight some key factors affecting interface stability. Although the simulations were conducted based on a single crystal, this work could lay a foundation for future research on polycrystalline and porous YSZ substrates. This work could indicate strategies for improving SOFC durability, such as optimizing material compositions or designing stress-relieving structures via identifying stress generation as a primary driver of degradation. The findings of this work could contribute to advancing SOFC technology towards sustainable and efficient energy solutions.

Key words

phase field method / solid oxide fuel cell / nickel-yttria-stabilized zirconia interface / cubic-tetragonal phase transformation / morphology evolution

引用本文

程凯, 焦震钧. 相场法模拟固体氧化物燃料电池阳极内Ni-YSZ界面处YSZ相变引发界面形貌演化. 硅酸盐学报, 2025 , 53 (10) : 2995 -3002 . DOI: 10.14062/j.issn.0454-5648.20250194

Kai CHENG, Zhenjun JIAO. Simulation of Evolution of Interface Morphology Induced by Yttria-Stabilized Zirconia (YSZ) Phase Transformation at Ni-YSZ Interface in Solid Oxide Fuel Cells Anode Using Phase Field Method[J]. Journal of the Chinese Ceramic Society, 2025 , 53 (10) : 2995 -3002 . DOI: 10.14062/j.issn.0454-5648.20250194

基金

收起

国家自然科学基金(11274218)
广东省高等院校创新团队计划(2021KCXTD006)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

2025年第53卷第10期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.14062/j.issn.0454-5648.20250194

接收时间：2025-03-21
首发时间：2026-10-21
出版时间：2025-10-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2025-03-21
修回日期：2025-04-16

基金

国家自然科学基金(11274218)

广东省高等院校创新团队计划(2021KCXTD006)

作者信息

哈尔滨工业大学（深圳）理学院，广东深圳 518055

通讯作者:

焦震钧（1982—），男，博士，教授。

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/gsyxb/CN/10.14062/j.issn.0454-5648.20250194

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT