| 217 | 0 | 416 |
| 下载次数 | 被引频次 | 阅读次数 |
固体氧化物燃料电池(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% H2 was supplied to the anode and pure O2 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 poly crystalline 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.
[1] PAN Z H, SHEN J, WANG J Y, et al. Thermodynamic analyses of a standalone diesel-fueled distributed power generation system based on solid oxide fuel cells[J]. Appl Energy, 2022, 308:118396.
[2] WANG Y L, DUNCAN K, WACHSMAN E D, et al. The effect of oxygen vacancy concentration on the elastic modulus of fluorite-structured oxides[J]. Solid State Ion, 2007, 178(1/2):53–58.
[3] BUTZ B, LEFARTH A, ST?RMER H, et al. Accelerated degradation of 8.5mol%Y2O3-doped zirconia by dissolved Ni[J]. Solid State Ion,2012, 214:37–44.
[4] KHAN M Z, SONG R H, HUSSAIN A, et al. Effect of applied current density on the degradation behavior of anode-supported flat-tubular solid oxide fuel cells[J]. J Eur Ceram Soc, 2020, 40(4):1407–1417.
[5] SHIN J S, SAQIB M, JO M, et al. Degradation mechanisms of solid oxide fuel cells under various thermal cycling conditions[J]. ACS Appl Mater Interfaces, 2021, 13(42):49868–49878.
[6] SHAUR A, REHMAN S U, KIM H S, et al. Hybrid electrochemical deposition route for the facile nanofabrication of a Cr-poisoningtolerant La(Ni, Fe)O3–δ cathode for solid oxide fuel cells[J]. ACS Appl Mater Interfaces, 2020, 12(5):5730–5738.
[7] YANG W Y, PAN Z H, JIAO Z J, et al. Advanced microstructure characterization and microstructural evolution of porous cermet electrodes in solid oxide cells:A comprehensive review[J]. Energy Rev,2025, 4(1):100104.
[8] PADINJARETHIL A K, BIANCHI F R, HAGEN A, et al. Steam and polarization effects on Ni–YSZ electrode due to degradation under electrolysis and fuel cell operation[J]. J Power Sources, 2025, 632:236296.
[9] HATTORI M, TAKEDA Y, LEE J H, et al. Effect of annealing on the electrical conductivity of the Y2O3–ZrO2 system[J]. J Power Sources,2004, 131(1–2):247–250.
[10] LINDEROTH S, BONANOS N, JENSEN K V, et al. Effect of NiO-to-Ni transformation on conductivity and structure of yttria-stabilized ZrO2[J]. J Am Ceram Soc, 2001, 84(11):2652–2656.
[11] MCPHAIL S J, FRANGINI S, LAURENCIN J, et al. Addressing planar solid oxide cell degradation mechanisms:A critical review of selected components[J]. Electrochem Sci Adv, 2022, 2(5):e2100024.
[12] TIETZ F, SEBOLD D, BRISSE A, et al. Degradation phenomena in a solid oxide electrolysis cell after 9000 h of operation[J]. J Power Sources, 2013, 223:129–135.
[13] SHIRATORI Y, TIETZ F, BUCHKREMER H P, et al. YSZ–Mg O composite electrolyte with adjusted thermal expansion coefficient to other SOFC components[J]. Solid State Ion, 2003, 164(1–2):27–33.
[14] KISHIMOTO H, YASHIRO K, SHIMONOSONO T, et al. In situ analysis on the electrical conductivity degradation of NiO doped yttria stabilized zirconia electrolyte by micro-Raman spectroscopy[J].Electrochim Acta, 2012, 82:263–267.
[15] DEVILLE S, CHEVALIER J. Martensitic relief observation by atomic force microscopy in yttria-stabilized zirconia[J]. J Am Ceram Soc, 2003,86(12):2225–2227.
[16] LAGUNA-BERCERO M A, CAMPANA R, LARREA A, et al.Electrolyte degradation in anode supported microtubular yttria stabilized zirconia-based solid oxide steam electrolysis cells at high voltages of operation[J]. J Power Sources, 2011, 196(21):8942–8947.
[17] MAMIVAND M, ASLE ZAEEM M, EL KADIRI H. Phase field modeling of stress-induced tetragonal-to-monoclinic transformation in zirconia and its effect on transformation toughening[J]. Acta Mater,2014, 64:208–219.
[18] LEI Y K, CHENG T L, WEN Y H. Phase field modeling of microstructure evolution and concomitant effective conductivity change in solid oxide fuel cell electrodes[J]. J Power Sources, 2017, 345:275–289.
[19] WILSON J R, KOBSIRIPHAT W, MENDOZA R, et al. Threedimensional reconstruction of a solid-oxide fuel-cell anode[J]. Nat Mater, 2006, 5:541–544.
[20] CHEN H Y, YU H C, SCOTT CRONIN J, et al. Simulation of coarsening in three-phase solid oxide fuel cell anodes[J]. J Power Sources, 2011, 196(3):1333–1337.
[21] JIAO Z J, SHIKAZONO N. Simulation of solid oxide fuel cell anode microstructure evolution using phase field method[J]. J Electrochem Soc, 2013, 160(6):F709–F715.
[22] JIAO Z J, SHIKAZONO N. Simulation of nickel morphological and crystal structures evolution in solid oxide fuel cell anode using phase field method[J]. J Electrochem Soc, 2014, 161(5):F577–F582.
[23] JENSEN K V, WALLENBERG R, CHORKENDORFF I, et al. Effect of impurities on structural and electrochemical properties of the Ni–YSZ interface[J]. Solid State Ion, 2003, 160(1–2):27–37.
[24] WANG K, ZHAO M, REN X R, et al. High temperature mechanical properties of zirconia metastable t'-phase degraded yttria stabilized zirconia[J]. Ceram Int, 2019, 45(14):17376–17381.
[25] DA Y L, XIAO Y, ZHONG Z, et al. Predictions on conductivity and mechanical property evolutions of yttria-stabilized zirconia in solid oxide fuel cells based on phase-field modeling of cubic–tetragonal phase transformation[J]. J Eur Ceram Soc, 2022, 42(8):3489–3499.
[26] ZHU J, CHEN L Q, SHEN J, et al. Coarsening kinetics from a variable-mobility cahn-Hilliard equation:Application of a semi-implicit Fourier spectral method[J]. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics, 1999, 60(4 Pt A):3564–3572.
[27] SEOL D J, HU S Y, LI Y L, et al. Computer simulation of spinodal decomposition in constrained films[J]. Acta Mater, 2003, 51(17):5173–5185.
[28] CHEN L Q, SHEN J. Applications of semi-implicit Fourier-spectral method to phase field equations[J]. Comput Phys Commun, 1998,108(2–3):147–158.
[29] JIAO Z J, SHIKAZONO N. In operando optical study of active three phase boundary of nickel-yttria stabilized zirconia solid-oxide fuel cell anode under polarization[J]. J Power Sources, 2018, 396:119–123.
[30] ARTEMEV A, JIN Y, KHACHATURYAN A G. Three-dimensional phase field model of proper martensitic transformation[J]. Acta Mater,2001, 49(7):1165–1177.
基本信息:
DOI：10.14062/j.issn.0454-5648.20250194
中图分类号:O646;TM911.4
引用信息:
[1]程凯,焦震钧.相场法模拟固体氧化物燃料电池阳极内Ni–YSZ界面处YSZ相变引发界面形貌演化[J].硅酸盐学报,2025,53(10):2995-3002.DOI:10.14062/j.issn.0454-5648.20250194.
基金信息:
国家自然科学基金（11274218）; 广东省高等院校创新团队计划（2021KCXTD006）