Next Article in Journal
Self-Discharge of a Proton Exchange Membrane Electrolyzer: Investigation for Modeling Purposes
Previous Article in Journal
Study of the Ecological Footprint and Carbon Footprint in a Reverse Osmosis Sea Water Desalination Plant
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Electro-Chemo-Mechanical Behaviors within the Dense BaZr0.8Y0.2O3−δ Protonic-Ceramic Membrane in a Long Tubular Electrochemical Cell

1
Department of Mechanical Engineering, Colorado School of Mines, Golden, CO 80401, USA
2
HyET Hydrogen USA LLC, Sacramento, CA 95826, USA
*
Author to whom correspondence should be addressed.
Membranes 2021, 11(6), 378; https://doi.org/10.3390/membranes11060378
Submission received: 23 April 2021 / Revised: 11 May 2021 / Accepted: 17 May 2021 / Published: 22 May 2021
(This article belongs to the Section Membrane Applications)

Abstract

:
This paper reports an extended Nernst–Planck computational model that couples charged-defect transport and stress in tubular electrochemical cell with a ceramic proton-conducting membrane. The model is particularly concerned with coupled chemo-mechanical behaviors, including how electrochemical phenomena affect internal stresses and vice versa. The computational model predicts transient and steady-state defect concentrations, fluxes, stresses within a thin BaZr 0.8 Y 0.2 O 3 δ (BZY20) membrane. Depending on the polarization (i.e., imposed current density), the model predicts performance as a fuel cell or an electrolyzer. A sensitivity analysis reveals the importance of thermodynamic and transport properties, which are often not readily available.

Graphical Abstract

1. Introduction

Materials such as doped barium zirconates (e.g., BaZr 1 x Y x O 3 δ , BZY) are perovskites that have good proton conductivity at intermediate temperatures ( T 500 °C) [1] and chemical stability in environments containing H 2 O and CO 2 [1,2,3,4]. Thus, these materials are suitable as proton-conducting membranes for applications that depend on hydrogen separations (e.g., fuel cells) [5,6]. Protonic ceramic membranes have also been employed to compress hydrogen up to 50 bar. Methane is reformed in the presence of steam on the Ni-BZY based ceramic-metal support, and the produced hydrogen is galvanically driven through the membrane and compressed on the other side of the membrane [7]. Although BZY materials are predominantly proton conductors in moist reducing atmospheres, they are mixed ionic-electronic conductors (MIEC) in oxidizing atmospheres with the presence of O-site polarons. As the temperature increases above 600 °C, the conduction through oxygen vacancies becomes not negligible. The present paper focuses on the chemo-mechanical behavior of BaZr 0.8 Y 0.2 O 3 δ , called BZY20.
The concentration variations associated with charged-defect transport within the membrane cause lattice-scale strain and volume deformation, often referred to as chemical expansion [8,9]. Because ceramic materials are brittle, defect-induced stresses can initiate distortions, cracks, and fracture. Moreover, as discussed in early works by Lee [10] and Larché and Cahn [11,12], local stress gradients also contribute to defect fluxes [13].
Atkinson et al. [14,15] analyzed the effects of defect-induced stress in solid-oxide fuel cells (SOFCs) and oxygen-ion-conducting membranes with alternative geometries and mechanical constraints. Euser et al. developed an extended Nernst–Planck–Poisson (NPP) model to evaluate stress in planar and radial oxygen-separation membranes [16,17,18]. Andersson et al. [19] and Han et al. [20] used X-ray diffraction to measure the chemical expansion of BZY for different dopant levels in moist and dry environments. Marrocchelli et al. [21] and Bishop et al. [22] reported detailed descriptions of chemical expansion in the perovskite materials. Dubois et al. [23] developed a coupled chemo-thermo-mechanical model using an extended Nernst-Planck formulation and evaluated the stresses associated with a protonic-ceramic fuel cell (PCFC) in a button-cell configuration. Ricote et al. [24] used high-temperature X-ray diffraction (HT-XRD) and thermogravimetric analysis (TGA) to study the effects of oxygen partial pressure on lattice parameters and oxygen non-stoichiometry in BaZr 0.9 Dy 0.1 O 3 δ (BZDy10).
The current study develops an electro-chemo-mechanical computational model based on an extended Nernst–Planck (NP) formulation. The one-dimensional radial model considers a tubular configuration with a thin dense BZY20 membrane supported by a relatively thick porous composite electrode (Figure 1). The model predicts transient profiles of defect concentrations and stress states as functions of applied electrical current. As with all models, results depend on physical parameters, including thermodynamic and transport properties. A sensitivity analysis seeks to identify the most important parameters and properties.

2. Defect Chemistry

The structures and behaviors of proton-conducting ceramics have been documented in numerous reviews [1,2,19,25,26,27,28,29,30,31]. Oxides with ABO 3 perovskite structures, such as BaZrO 3 , are among the most promising proton conductors. Effective proton conductivity is established by B-site doping with a trivalent cation such as yttrium (Y 3 + ) to replace a fraction of the Zr 4 + , thus introducing oxygen vacancies to preserve local charge neutrality. The present study considers 20% Y doping, BaZr 0.8 Y 0.2 O 3 δ . The needed thermodynamic, transport, and mechanical properties are derived from prior publications [1,2,4,19].
Figure 1 illustrates aspects of the tubular geometry and the transport processes. The profiles and fluxes of three mobile charged defects (protons OH O , oxygen vacancies V O and polarons O O ) depend on the current density via an external circuit. Depending on the polarization (i.e., the direction of the electrical current), the cell may operate as a fuel cell or as an electrolyzer [32]. As illustrated in Figure 1, the present study considers gas-phase mixtures of 20% H2O and 80% O2 on the exterior electrode (positrode) and 97% H2 and 3% H2O on the interior electrode (negatrode).
Defect-incorporation reactions at the membrane boundaries may be stated in Kröger–Vink notation as
1 2 H 2 + O O OH O ,
1 2 O 2 + O O × + V O 2 O O ,
H 2 O + V O + O O × 2 OH O .
Corresponding equilibrium constants for the incorporation reactions are
K p , H 2 = OH O L p H 2 1 / 2 O O L ,
K p , O 2 = O O L 2 p O 2 1 / 2 O O X L V O L ,
K p , H 2 O = OH O L 2 p H 2 O O O × L V O L ,
where p k are the gas phase partial pressures and the defect concentrations are represented in lattice units. The defect molar concentrations may be evaluated using the molar volume (for BZY20, V m = 4.57 × 10 5 m3 mol−1) as
X k L = X k V m .
The present model includes polaron traps [4]. In other words, some fraction of the otherwise highly mobile small polarons can be trapped, thus immobilized in the proximity of the yttrium dopant. Stated as a reaction,
Y Zr + O O Y Zr O O .
The associated equilibrium constant follows as
K p , Trap = Y Zr O O L Y Zr L O O L .
The present model also assumes that the gas phases are equilibrated on both sides of the cell,
H 2 + 1 2 O 2 H 2 O .
Because the gas phase is equilibrated, the defect-incorporation equilibrium constants are not all independent; they are constrained as
K p , H 2 2 K p , O 2 = K p , H 2 O K p , G .
The equilibrium constants for each reaction can be evaluated via thermodynamic properties as
K p = exp Δ G R T = exp Δ H R T exp Δ S R ,
where T is temperature, R is the gas constant, Δ G is the change in Gibbs free energy, Δ S is the change in defect entropy, and Δ H represents change in enthalpy. Table 1 lists the thermodynamic properties used in the present study [4].
Assuming the defect reactions at the membrane surfaces are at equilibrium, defect concentrations at surfaces are constants that depend on the gas phase composition. Evaluating the boundary concentrations must consider site and electroneutrality balances. Charge neutrality requires that
2 V O L + OH O L + O O L Y Zr L = 0 .
The ABO3 pervoskite structure further constrains the oxygen site balance as
V O L + OH O L + O O L + O O x L + Y Zr O O L = 3 .
The yttrium doping is fixed as Y Zr L = 0.2 , which also constrains the trap balance as
Y Zr L + Y Zr O O L = Y Zr L .
The defect concentrations at the boundaries are calculated by solving Equations (5), (6), (9), (13)–(15) simultaneously. Algorithmic details about the calculation may be found in Zhu et al. [4].

3. Extended NP Membrane Model

The defect concentrations are governed by a conservation principle, represented mathematically as
[ X k ] t + J k = ω ˙ k ,
where t is time, J k are defect transport fluxes, and ω ˙ k are the net molar production rates of the polaron traps within the membrane [4]. In the current work, standard Nernst–Planck fluxes are extended to incorporate the contribution from hydrostatic-stress gradients. The defect-transport fluxes include the effects of diffusion J k D , migration, J k M , and hydrostatic stress J k S ,
J k = J k D + J k M + J k S .
The flux may be represented as
J k = D k [ X k ] + z k F R T [ X k ] Φ e 3 β k R T [ X k ] σ h ,
where D k are defect diffusion coefficients, z k are the defect’s charge, Φ e is the electrostatic potential, β k is the coefficient of chemical expansion, and σ h is the hydrostatic stress. The defect charges are z V O = + 2 , z OH O = + 1 , and z O O = + 1 . Hydrostatic stress within the membrane is defined as
σ h = σ r r + σ θ θ + σ z z 3 .
The polaron traps are produced and consumed according the reaction represented by Equation (8). The rate of progress of the trap reaction may be evaluated via mass-action kinetics as
q ˙ Trap = k f [ Y Zr ] [ O O ] k b Y Zr O O ,
where k f and k b are forward and backward rate constants, which are related through the equilibrium constant. The present model assumes k f = 10 6 . The species production rate due to the trap reaction can be expressed as ω ˙ Y Zr = q ˙ Trap , ω ˙ O O = q ˙ Trap , and ω ˙ Y Zr O O = q ˙ Trap .

3.1. Chemical-Expansion Coefficient

The chemical-expansion coefficient β k is defined to be the change in lattice volume with respect to the changes in defect concentration. As reported by Andersson et al. [19], the present model uses β OH O , L = 0.033 (in lattice units). In molar units, β k = β OH O , L V m = 1.5 × 10 6 m3 mol−1. Ricote et al. [24] showed that the O-site polarons produce negligible chemical expansion. The present model attributes all chemical expansion to the protons, neglecting lattice distortion due to polarons and oxygen vacancies.

3.2. Electrostatic Potential Φ e

The electrostatic potential may be calculated by solving the Gauss equation,
E t = 1 ε 0 ε r i k z k F J k ,
where the electrostatic-potential field is E = Φ e , ε 0 and ε r are vacuum and relative permittivities, respectively, and i is the external current density. However, the present model implements an approximately equivalent method by enforcing strict local electroneutrality. In this limit, the local charge flux must vanish as
k z k F X k = 0 .
The electrostatic potential on one of the boundaries is set to a reference Φ e = 0 , while the other boundary is set to impose the net current density i through the membrane. In other words, the charge flux through the membrane is exactly balanced by the current density i through the external circuit,
k z k F J k = i .
This algebraic equation, together with Equation (16), is sufficient to determine the local electrostatic-potential gradient Φ e .

3.3. Defect Diffusion Coefficients

The charged-defect diffusion coefficients are represented in an Arrhenius equation as
D k = D k exp E k R T ,
where D k and E k are pre-exponential factors and activation energies, respectively. As reported by Zhu et al. [4], Table 2 lists the values of D k and E k for three mobile charged defect.

3.4. Stress, Strain, and Displacement

The isotropic stress–strain relationship including chemically induced stress may be written as
σ = λ tr ( ϵ ) I + 2 G ϵ ( 3 λ + 2 G ) k β k Δ [ X k ] I .
where σ is the stress tensor, ϵ is the strain tensor, I is the unit tensor, and Δ [ X k ] represents the change in concentration with respect to a reference zero-strain state. The current study does not incorporate the residual strain related to the fabrication process [5,23]. The Lamé constant λ and shear modulus G are defined as
λ = E m ν m ( 1 + ν m ) ( 1 2 ν m ) ,
G = E m 2 ( 1 + ν m ) .
In these expressions, E m is the Young’s modulus, and ν m is Poisson’s ratio of the membrane. The tube is assumed to be infinitely long, resulting in a plane-strain state. Consequently, the strain components in the axial z direction vanish,
ϵ z z = ϵ z r = ϵ z θ = 0 .
The strain–displacement relationships with respect to radial displacement u are expressed as
ϵ r r = u r , ϵ θ θ = u r , ϵ r θ = 0 .
Assuming a quasi-static stress field and the absence of body forces, the divergence of the stress must vanish,
σ = 0
In one-dimensional radial coordinates, the stresses can be represented as
σ r r r + 1 r σ r θ θ + σ r z z + σ r r σ θ θ r = 0 .
Considering the tube to be long and axisymmetric, variations in axial z and circumferential θ stresses can be neglected. The radial component of the stress equilibrium simplifies to
σ r r r + σ r r σ θ θ r = 0 .
In practice, protonic-ceramic membranes should be as thin as possible to promote high proton fluxes (usually on the order of tens of microns). Such very thin membranes must be supported with relatively thick porous support structures that also serve as electrodes (cf., Figure 1). In reducing environments (e.g., anodes of protonic-ceramic fuel cells) the most widely used supports are composites of Ni-BZY. Such structures provide the mechanical strength and stability. The percolating Ni phase serves dual roles as the hydrogen-reduction catalyst and the electrical conductor. Table 3 lists the dimensions of tube used for the present study. The membrane is 10 μm thick and the porous support is 1 mm thick. The outer electrode is specified to be thin, but not explicitly included in the present model.
The current study, which specifically focuses on the electro-chemo-mechanical behavior within the membrane, neglects defect transport and chemical expansion within the porous support. Table 3 also includes mechanical properties of BZY20 and Ni-BZY20. These properties were determined experimentally using the ultrasonic measurements [23].
The mechanical aspects of the present model follow the approach reported by Euser et al. [17]. By substituting the stress–strain relationships (Equation (25)) and strain–displacement relationships (Equation (29)), the equilibrium equation (Equation (32)) can be solved analytically to produce the local displacements u, assuming a perfect bond between support and membrane properly satisfies continuity of displacement and traction at the membrane boundary. The constants of integration can be evaluated assuming inner and outer surfaces of the tube are traction-free. Using these displacements, the local stress components within the membrane and porous support can be evaluated.

3.5. Computational Implementation

The computational model is formulated using the method-of-lines and a one-dimensional radial finite-volume mesh network. The model incorporates two-way coupling between the defect transport and stress in the membrane. Transport induced stress is evaluated using the analytical solution. Euser et al. [17] reported the radial displacement and stress for the tubular membrane, respectively, as
u ( r ) = 1 + v m 1 v m I ( r ) r + C 1 r + C 2 r ,
σ r r ( r ) = E m 1 v m I ( r ) r 2 + E m 1 + v m 1 2 v m C 1 E m 1 + v m C 2 r 2 ,
where C 1 and C 2 are constants of integration, which can be calculated using the mechanical boundary conditions, as discussed earlier. In addition,
I ( r ) = r i r r k β k Δ X k d r .
where r i is the inner radius of membrane. The transport problem is solved numerically. Equation (16) is effectively a parabolic partial differential equation that is solved in Matlab using the ode15i function. The second-order flux derivatives are approximated using conservative central differences. The migration J k M and stress J k M fluxes are discretized using an up-winding scheme, where the coefficient of [ X k ] behaves as an artificial velocity [33,34]. The simulations run efficiently and quickly on a typical personal computer.

4. Results and Discussion

4.1. Concentration Profiles

Figure 2 shows the predicted steady-state lattice-scale concentration profiles [ X ] L of charged defects for current densities ranging between i = 2 A cm 2 (electrolysis) and i = 2 A cm 2 (fuel cell). The proton concentrations are larger than, but comparable to, the vacancy and the trap concentrations. The O-site polaron concentrations are much lower. Even at open circuit, there is significant curvature in the defect-concentration profiles. The sense of the curvatures changes depending on the polarization (i.e., direction of the charge flux). Under fuel-cell polarization ( i > 0 ), the profiles reveal boundary-layer behavior near the O 2 /H 2 O side of the membrane. Interestingly, under very strong fuel-cell polarization, a local minimum appears in the proton profile.
Figure 3 illustrates transient behavior of the concentration profiles, beginning from open-circuit conditions as a current density of i = 1 A cm 2 is suddenly imposed. The characteristic time constant for the transient is on the order of one second. The vacancies, traps, and polarons behave as may be anticipated, transitioning smoothly between the two polarization states. The proton concentrations reveal somewhat more unusual behavior. The proton concentrations initially increase, before later approaching the lower concentrations at the i = 1 A cm 2 steady-state. Although the magnitudes of the concentration variations are not large, there could be implications on the local transient stresses.

4.2. Proton Flux Profiles

Proton flux is an important measure of membrane performance. Figure 4 shows the predicted proton flux as a function of current density for selected oxygen concentrations at the tube’s exterior electrode. Interestingly, under fuel-cell polarization ( i > 0 ), the flux is only weakly affected by the O 2 concentration. The effect is much more significant under electrolysis polarization. As may be anticipated, high steam concentrations (i.e., low O 2 concentration) lead to higher proton fluxes under electrolysis conditions.
At steady state, the defect-flux profiles must be essentially flat (i.e., very little spatial dependence). However, because of the radial coordinates, there will be very slight radial variations in the steady-state fluxes. The net proton flux comprises diffusion, migration, and stress contributions (Equation (17)). These flux contributions do vary spatially, even at steady state. Figure 5 plots separately the individual steady-state flux-contribution profiles under a range of imposed current densities. These results show that the migration contributions generally dominate. Especially under fuel-cell polarization, the diffusion and stress contributions have strong gradients near the outer boundary. The behaviors are consistent with those reported by Euser et al. [16,17] for La 0.6 Sr 0.4 Co 0.8 Fe 0.2 O 3 δ (LSCF) oxygen-transport membranes.
Figure 6 shows proton-flux profiles during a transient between open circuit and a suddenly imposed current density of i = 1 A cm 2 . Even under open-circuit conditions, there is a small proton flux. Upon polarization, the flux profiles have spatial gradients, with the highest fluxes near the hydrogen fuel boundary. However, within a few seconds, the proton-flux profile becomes spatially uniform again at a significantly higher value, driving protons from the fuel side toward the oxygen side.

4.3. Stress Profiles

Figure 7 shows hydrostatic stress profiles within the membrane for different imposed current densities. The stresses in the fuel cell mode ( i > 0 ) are considerably higher than in electrolysis mode ( i < 0 ). The fuel-cell stresses have local maxima within the membrane, slightly exceeding the high fuel-side boundary stresses. Moreover, the fuel-cell stress gradients are high, especially near the oxygen boundary. The curvature of the hydrostatic-stress profiles is positive for electrolysis, while it is negative at open circuit and in fuel-cell operation. Under strong electrolysis polarization, there can be a slight local minimum in the steady-state hydrostatic-stress profile. The present model attributes chemical strain to the protons. Comparing Figure 2a and Figure 7 shows that the correlation between the proton concentration profile and the hydrostatic stress profiles is evident.
Figure 8 shows profiles of the normal stress components that contribute to the hydrostatic stress. The radial stresses are compressive while hoop and axial stresses are tensile and nearly identical. The hydrostatic stress distribution is similar to hoop and axial stresses, although with lower magnitude since radial stresses are compressive (Figure 7). The stresses in the porous support are negligible, as chemical expansion in the support is omitted.
Figure 9 shows the transient response of the hydrostatic stress as a current density of i = 1 A cm 2 is suddenly imposed from an initial open-circuit steady-state condition. The local stresses decrease at early times, with local upward curvature at a time of around one second. By about five seconds, the steady-state hydrostatic-stress profile is achieved. The transient stress profiles behave similar to the transient proton concentrations profiles (Figure 3a).
In a practically operating cell, the gas compositions and temperature can vary owing to reagent depletion, dilution, etc. Cell electrochemical performance generally depends on gas-phase concentrations. Local stresses are also affected by gas-phase composition. Figure 10 shows the maximum hydrostatic stress as a function of oxygen concentration for the imposed current density of i = 1 A cm 2 . As the oxygen partial pressure at the cathode boundary increases, the maximum hydrostatic stress decreases considerably. This curve does not depend significantly on the current density. Consequently, maintaining high oxygen concentration could enhance mechanical integrity.

4.4. Sensitivity Analysis

There is considerable uncertainty in some of the physical properties for BZY20. Thus, it is interesting to investigate the sensitivity of predicted results to the model’s properties and parameters. Figure 11 is a tornado plot that shows the effects of varying individual properties by ± 10 % on the predicted proton flux under fuel-cell polarization. The nominal model-predicted proton flux is J OH O = 0.1036 mol s 1 m 2 . Under these conditions, the highest sensitivities were found to be the activation energy of the proton diffusion coefficient E OH O , the enthalpy of the steam-incorporation reaction Δ H H 2 O , the activation energy of the vacancy diffusion coefficient E V O , and the entropy of the steam-incorporation reaction Δ S H 2 O .
Figure 12 is a tornado plot showing sensitivities of model-predicted maximum hydrostatic stress to model parameters. The cell is operating at 600 C with an imposed fuel-cell current density of i = 1 A cm 2 . Under these conditions, the nominal maximum hydrostatic stress is σ h = 128.7 MPa. The highest sensitivities are to the thermodynamic properties (enthalpy Δ H and entropy Δ S ) of the defect-incorporation reactions. Interestingly, mechanical properties (e.g., Young’s modulus, E m ) show relatively small sensitivity.
The present model uses best-available properties [4]. However, results of the sensitivity analysis point to needs for measuring, improving, and validating thermodynamic and transport properties. Such research is ongoing, but is a long-term process.

5. Summary and Conclusions

The performance of BZY20 as a proton-conducting fuel-cell or electrolysis membrane is simulated using an extended Nernst–Planck model. The model considers the one-dimensional radial behavior in a long tubular cell. The Nernst–Planck fluxes depend on local gradients of defect concentration, electrostatic potential, and hydrostatic stress. The model predicts transient and steady state profiles of defect concentrations and stress within a thin (10 μm) BZY20 membrane that is supported on a relatively thick porous Ni/BZY20 composite electrode. The present model concentrates on defect transport within the dense membrane, assuming ideal electrodes. In other words, charge-transfer polarization at the electrodes is neglected. The three charged defects (protons, oxygen vacancies, and small polarons) are incorporated into the membrane via equilibrated defect-incorporation reactions.
An important measure of membrane performance is the proton flux. However, as a practical matter, internal stresses are also important because of potential damage mechanisms with brittle ceramic materials. The present model couples the electrochemical performance and chemo-mechanical performance. A sensitivity analysis reveals an ongoing need to validate material-specific thermodynamic, transport, and mechanical properties, thus improving predictive capabilities of the models.

Author Contributions

Conceptualization, K.T., A.D., J.R.B., S.R., H.Z. and R.J.K.; methodology, K.T., A.D., J.R.B. and R.J.K.; software, K.T. and A.D.; validation, K.T., A.D., J.R.B., S.R., H.Z. and R.J.K.; formal analysis, K.T., A.D., J.R.B., S.R., H.Z. and R.J.K.; writing—original draft preparation, K.T.; writing—review and editing, K.T., A.D., J.R.B., S.R., H.Z. and R.J.K.; visualization, K.T. and R.J.K.; supervision, J.R.B., S.R. and R.J.K.; project administration, J.R.B. and R.J.K.; funding acquisition, J.R.B. and R.J.K. All authors have read and agreed to the published version of the manuscript.

Funding

We gratefully acknowledge support from the Colorado School of Mines Foundation via the Angel Research Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kreuer, K. Proton-conducting oxides. Annu. Rev. Mater. 2003, 33, 333–359. [Google Scholar] [CrossRef] [Green Version]
  2. Fabbri, E.; Pergolesi, D.; Licoccia, S.; Traversa, E. Does the increase in Y-dopant concentration improve the proton conductivity of BaZr1−xYxO3−δ fuel cell electrolytes? Solid State Ion. 2010, 181, 1043. [Google Scholar] [CrossRef]
  3. Imashuku, S.; Uda, T.; Nose, Y.; Taniguchi, G.; Ito, Y.; Awakura, Y. Dependence of dopant cations on microstructure and proton conductivity of barium zirconate. J. Electrochem. Soc. 2009, 156, BI. [Google Scholar] [CrossRef]
  4. Zhu, H.; Ricote, S.; Duan, C.; O’Hayre, R.; Tsvetkov, D.; Kee, R. Defect Incorporation and Transport within Dense BaZr0.8Y0.2O3−δ(BZY20) Proton-Conducting Membranes. J. Electrochem. Soc. 2018, 165, F581–F588. [Google Scholar] [CrossRef]
  5. Duan, C.; Tong, J.; Shang, M.; Nikodemski, S.; Sanders, M.; Ricote, S.; Almansoori, A.; O’Hayre, R. Readily processed protonic ceramic fuel cells with high performance at low temperature. Science 2015, 349, 1321–1326. [Google Scholar] [CrossRef] [PubMed]
  6. Duan, C.; Kee, R.; Zhu, H.; Karakaya, C.; Chen, Y.; Ricote, S.; Jarry, A.; Crumlin, E.; Hook, D.; Braun, R.; et al. Highly durable, coking and sulfur tolerant, fuel-flexible protonic ceramic fuel cell. Nature 2018, 557, 65–77. [Google Scholar] [CrossRef] [PubMed]
  7. Malerød-Fjeld, H.; Clark, D.; Yuste-Tirados, I.; Zanón, R.; Catalán-Martinez, D.; Beeaff, D.; Morejudo, S.; Vestre, P.; Norby, T.; Haugsrud, R.; et al. Thermo-electrochemical production of compressed hydrogen from methane with near-zero energy loss. Nature Energy 2017, 2, 923–931. [Google Scholar] [CrossRef]
  8. Adler, S. Chemical expansivity of electrochemical ceramics. J. Am. Ceram. Soc. 2001, 84, 2117–2119. [Google Scholar] [CrossRef]
  9. Chen, X.; Yu, J.; Adler, S. Thermal and chemical expansion of Sr-doped lanthanum cobalt oxide (La1−xSrxCoO3−δ). Chem. Mater. 2005, 17, 4537–4546. [Google Scholar] [CrossRef]
  10. Li, J.C.M. Physical chemistry of some microstructural phenomena. Metall. Mater. Trans. A 1978, 9, 1353–1380. [Google Scholar]
  11. Larché, F.; Cahn, J. The effect of self-stress on diffusion in solids. Acta Metall. 1982, 30, 1835–1846. [Google Scholar] [CrossRef]
  12. Larché, F.C.; Cahn, J.W. The Interactions of composition and stress in crystalline solids. Acta Metall. 1985, 33, 331–357. [Google Scholar] [CrossRef]
  13. Yang, F. Interaction between diffusion and chemical stresses. Mater. Sci. Eng. 2005, 409, 153–159. [Google Scholar] [CrossRef]
  14. Atkinson, A. Chemically-induced stresses in gadolinium-doped ceria solid oxide fuel cell electrolytes. Solid State Ion. 1997, 95, 249. [Google Scholar] [CrossRef]
  15. Atkinson, A.; Ramos, T. Chemically-induced stresses in ceramic oxygen ion-conducting membranes. Solid State Ion. 2000, 129, 259. [Google Scholar] [CrossRef]
  16. Euser, B.; Berger, J.; Zhu, H.; Kee, R. Defect-Transport-Induced Stress in Mixed Ionic-Electronic Conducting (MIEC) Ceramic Membranes. J. Electrochem. Soc. 2016, 163, F264–F271. [Google Scholar] [CrossRef]
  17. Euser, B.; Berger, J.; Zhu, H.; Kee, R. Chemically Induced Stress in Tubular Mixed Ionic-Electronic Conducting (MIEC) Ceramic Membranes. J. Electrochem. Soc. 2016, 163, F1294–F1301. [Google Scholar] [CrossRef]
  18. Euser, B.; Zhu, H.; Berger, J.; Lewinsohn, C.; Kee, R. Electrochemical-mechanical coupling in composite planar structures that integrate flow channels and ion-conducting membranes. J. Electrochem. Soc. 2017, 164, F732. [Google Scholar] [CrossRef]
  19. Andersson, A.; Selbach, S.; Knee, C.; Grande, T. Chemical Expansion Due to Hydration of Proton-Conducting Perovskite Oxide Ceramics. J Am. Ceram. Soc. 2014, 97, 2654–2661. [Google Scholar] [CrossRef] [Green Version]
  20. Han, D.; Hatada, N.; Uda, T. Chemical Expansion of Yttrium-Doped Barium Zirconate and Correlation with Proton Concentration and Conductivity. J. Am. Ceram. Soc. 2016, 99, 3745–3753. [Google Scholar] [CrossRef]
  21. Marrocchelli, D.; Perry, N.; Bishop, S. Understanding chemical expansion in perovskite-structured oxides. Phys. Chem. Chem. Phys. 2015, 17, 10028–10039. [Google Scholar] [CrossRef] [PubMed]
  22. Bishop, S.; Marrocchelli, D.; Chatzichristodoulou, C.; Perry, N.; Mogensen, M.; Tuller, H.; Wachsman, E. Chemical Expansion: Implications for Electrochemical Energy Storage and Conversion Devices. Annu. Rev. Mater. Res. 2014, 44, 205–239. [Google Scholar] [CrossRef]
  23. Dubois, A.; Taghikhani, K.; Berger, J.; Zhu, H.; O’Hayre, R.; Braun, R.; Kee, R.; Ricote, S. Chemo-thermo-mechanical coupling in protonic ceramic fuel cells from fabrication to operation. J. Electrochem. Soc. 2019, 166, F1007–F1015. [Google Scholar] [CrossRef]
  24. Ricote, S.; Hudish, G.; O’Brien, J.; Perry, N. Non stoichiometry and lattice expansion of BaZr0.9Dy0.1O3−δ in oxidizing atmospheres. Solid State Ion. 2019, 330, 33–39. [Google Scholar] [CrossRef]
  25. Yamanaka, S.; Fujikane, M.; Hamaguchi, T.; Muta, H.; Oyama, T.; Matsuda, T.; Kobayashi, S.I.; Kurosaki, K. Thermophysical properties of BaZrO3 and BaCeO3. J. Alloy. Compd. 2003, 359, 109–113. [Google Scholar] [CrossRef]
  26. Schober, T.; Krug, F.; Schilling, W. Criteria for the application of high temperature proton conductors in SOFCs. Solid State Ion. 1997, 97, 369–373. [Google Scholar] [CrossRef]
  27. Hashim, S.; Somalu, M.; Loh, K.; Liu, S.; Zhou, W.; Sunarso, J. Perovskite-based proton conducting membranes for hydrogen separation: A review. Int. J. Hydrog. 2018, 43, 15281–15305. [Google Scholar] [CrossRef]
  28. Løken, A.; Ricote, S.; Wachowski, S. Thermal and chemical expansion in proton ceramic electrolytes and compatible electrodes. Crystals 2018, 8, 365. [Google Scholar] [CrossRef] [Green Version]
  29. Meng, Y.; Gao, J.; Zhao, Z.; Amoroso, J.; Tong, J.; Brinkman, K. recent progress in low-temperature proton-conducting ceramics. J. Mater. Sci. 2019, 54, 9291–9312. [Google Scholar] [CrossRef] [Green Version]
  30. Lei, L.; Zhang, J.; Yuan, Z.; Liu, J.; Ni, M.; Chen, F. Progress Report on Proton Conducting Solid Oxide Electrolysis Cells. Adv. Funct. Mater. 2019, 29, 1903805. [Google Scholar] [CrossRef]
  31. Medvedev, D. Trends in research and development of protonic ceramic electrolysis cells. Int. J. Hydrog. Energy 2019, 44, 26711–26740. [Google Scholar] [CrossRef]
  32. Zhu, H.; Kee, R. Membrane polarization in mixed-conducting ceramic fuel cells and electrolyzers. Int. J. Hydrog. 2016, 41, 2931–2943. [Google Scholar] [CrossRef] [Green Version]
  33. Kee, R.; Zhu, H.; Hildenbrand, B.; Vøllestad, E.; Sanders, M.; O’Hayre, R. Modeling the steady-state and transient response of polarized and non-polarized proton-conducting doped-perovskite membranes. J. Electrochem. Soc. 2013, 160, F290–F300. [Google Scholar] [CrossRef]
  34. Vøllestad, E.; Zhu, H.; Kee, R. Interpretation of defect and gas-phase fluxes through mixed-conducting ceramics using Nernst–Planck–Poisson and integral formulations. J. Electrochem. Soc. 2014, 161, F114–F124. [Google Scholar] [CrossRef]
Figure 1. Proton-conducting tubular cell consisting of a dense BZY20 electrolyte supported by porous Ni-BZY20 operating as a fuel cell.
Figure 1. Proton-conducting tubular cell consisting of a dense BZY20 electrolyte supported by porous Ni-BZY20 operating as a fuel cell.
Membranes 11 00378 g001
Figure 2. Charged-defect concentration profiles (lattice units) across the membrane for imposed current densities ( i = 2 , 1 , 0 , 1 , 2 A cm 2 ) at 600 C. (a) protons, (b) oxygen vacancies, (c) trapped polarons, and (d) O-site polarons.
Figure 2. Charged-defect concentration profiles (lattice units) across the membrane for imposed current densities ( i = 2 , 1 , 0 , 1 , 2 A cm 2 ) at 600 C. (a) protons, (b) oxygen vacancies, (c) trapped polarons, and (d) O-site polarons.
Membranes 11 00378 g002
Figure 3. Charged-defect concentration profiles (lattice units) during a transient between open circuit and a suddenly imposed current density of i = 1 A cm 2 . The cell temperatures is fixed at 600 C. (a) protons, (b) oxygen vacancies, (c) trapped polarons, and (d) O-site polarons.
Figure 3. Charged-defect concentration profiles (lattice units) during a transient between open circuit and a suddenly imposed current density of i = 1 A cm 2 . The cell temperatures is fixed at 600 C. (a) protons, (b) oxygen vacancies, (c) trapped polarons, and (d) O-site polarons.
Membranes 11 00378 g003
Figure 4. Protonic flux change as a function of current density for selected air side oxygen compositions ranging from 1% to 80% at 600 C. The oxygen concentration is balanced with steam at the boundary.
Figure 4. Protonic flux change as a function of current density for selected air side oxygen compositions ranging from 1% to 80% at 600 C. The oxygen concentration is balanced with steam at the boundary.
Membranes 11 00378 g004
Figure 5. Steady-state proton-flux contributions for selected current densities ( i = 2 , 1 , 0 , 1 , 2 A cm 2 ) at 600 C. (a) diffusion flux, (b) migration flux, and (c) stress-induced flux.
Figure 5. Steady-state proton-flux contributions for selected current densities ( i = 2 , 1 , 0 , 1 , 2 A cm 2 ) at 600 C. (a) diffusion flux, (b) migration flux, and (c) stress-induced flux.
Membranes 11 00378 g005
Figure 6. Transient proton-flux profiles within the membrane, beginning at open circuit mode and suddenly imposing a fuel-cell current density i = 1 A cm 2 . The cell temperature is fixed at 600 C.
Figure 6. Transient proton-flux profiles within the membrane, beginning at open circuit mode and suddenly imposing a fuel-cell current density i = 1 A cm 2 . The cell temperature is fixed at 600 C.
Membranes 11 00378 g006
Figure 7. Steady-state hydrostatic stress profiles for selected current densities with the cell operating at 600 C.
Figure 7. Steady-state hydrostatic stress profiles for selected current densities with the cell operating at 600 C.
Membranes 11 00378 g007
Figure 8. Steady-state stress-contribution profiles for selected current densities ( i = 2 , 1 , 0 , 1 , 2 A cm 2 ) with the cell operating at 600 C. (a) radial, (b) hoop, and (c) axial stress components.
Figure 8. Steady-state stress-contribution profiles for selected current densities ( i = 2 , 1 , 0 , 1 , 2 A cm 2 ) with the cell operating at 600 C. (a) radial, (b) hoop, and (c) axial stress components.
Membranes 11 00378 g008
Figure 9. Transient hydrostatic stress profiles within the membrane, beginning at open circuit and suddenly imposing a fuel-cell current density of i = 1 A cm 2 .
Figure 9. Transient hydrostatic stress profiles within the membrane, beginning at open circuit and suddenly imposing a fuel-cell current density of i = 1 A cm 2 .
Membranes 11 00378 g009
Figure 10. Maximum hydrostatic stress as a function of oxygen fraction at the outer boundary under i = 1 A cm 2 imposed current density. The oxygen fraction is balanced with steam. The cell is fixed at a temperature of 600 C.
Figure 10. Maximum hydrostatic stress as a function of oxygen fraction at the outer boundary under i = 1 A cm 2 imposed current density. The oxygen fraction is balanced with steam. The cell is fixed at a temperature of 600 C.
Membranes 11 00378 g010
Figure 11. Sensitivity of steady-state fuel-cell proton flux to 10 % changes in electro-chemo-mechanical properties. The imposed current density is i = 1 A cm 2 .
Figure 11. Sensitivity of steady-state fuel-cell proton flux to 10 % changes in electro-chemo-mechanical properties. The imposed current density is i = 1 A cm 2 .
Membranes 11 00378 g011
Figure 12. Sensitivity of maximum hydrostatic stress to 10 % changes in electro-chemo-mechanical properties. The imposed current density is i = 1 A cm 2 .
Figure 12. Sensitivity of maximum hydrostatic stress to 10 % changes in electro-chemo-mechanical properties. The imposed current density is i = 1 A cm 2 .
Membranes 11 00378 g012
Table 1. Enthalpy and entropy of defect incorporation reactions in BZY20. Reproduced with permission from Zhu et al. [4]. Copyright The Electrochemical Society, 2018.
Table 1. Enthalpy and entropy of defect incorporation reactions in BZY20. Reproduced with permission from Zhu et al. [4]. Copyright The Electrochemical Society, 2018.
Reaction Δ H
( kJ mol 1 )
Δ S
( J mol 1 K 1 )
1 2 H 2 + O O OH O −228.36−54.80
1 2 O 2 + O O × + V O 2 O O 115.31−45.89
H 2 O + V O + O O × 2 OH O −93.30−100.00
Y Zr + O O Y Zr O O −90.30−6.71
H 2 + 1 2 O 2 H 2 O −248.11−55.48
Table 2. Pre-exponential factors and activation energies for diffusivities of mobile charged defects in BZY20. Reproduced with permission from Zhu et al. [4]. Copyright The Electrochemical Society, 2018.
Table 2. Pre-exponential factors and activation energies for diffusivities of mobile charged defects in BZY20. Reproduced with permission from Zhu et al. [4]. Copyright The Electrochemical Society, 2018.
Charged Defects D k
(m 2 s 1 )
E k
(kJ mol 1 )
OH O 5.18 × 10 7 60.66
V O 2.03 × 10 7 85.19
O O 1.38 × 10 5 7.18
Table 3. Thermo-mechanical and geometrical properties of model.
Table 3. Thermo-mechanical and geometrical properties of model.
ParameterValue
Temperature, T 600 C
Molar volume, V m 4.57 × 10 5 m 3 mol 1 [4]
Outer support radius, r s 4 × 10 3 m
Inner membrane radius, r i 5 × 10 3 m
Outer membrane radius, r o 5.01 × 10 3 m
Membrane Young’s modulus, E m 205 GPa
BZY20
Support Young’s modulus, E s 108 GPa
Ni-BZY20 (Dry–25% porosity)
Membrane Poisson’s ratio, ν m 0.24
BZY20
Support Poisson’s ratio, ν s 0.27
Ni-BZY20 (Dry–25% porosity)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Taghikhani, K.; Dubois, A.; Berger, J.R.; Ricote, S.; Zhu, H.; Kee, R.J. Modeling Electro-Chemo-Mechanical Behaviors within the Dense BaZr0.8Y0.2O3−δ Protonic-Ceramic Membrane in a Long Tubular Electrochemical Cell. Membranes 2021, 11, 378. https://doi.org/10.3390/membranes11060378

AMA Style

Taghikhani K, Dubois A, Berger JR, Ricote S, Zhu H, Kee RJ. Modeling Electro-Chemo-Mechanical Behaviors within the Dense BaZr0.8Y0.2O3−δ Protonic-Ceramic Membrane in a Long Tubular Electrochemical Cell. Membranes. 2021; 11(6):378. https://doi.org/10.3390/membranes11060378

Chicago/Turabian Style

Taghikhani, Kasra, Alexis Dubois, John R. Berger, Sandrine Ricote, Huayang Zhu, and Robert J. Kee. 2021. "Modeling Electro-Chemo-Mechanical Behaviors within the Dense BaZr0.8Y0.2O3−δ Protonic-Ceramic Membrane in a Long Tubular Electrochemical Cell" Membranes 11, no. 6: 378. https://doi.org/10.3390/membranes11060378

APA Style

Taghikhani, K., Dubois, A., Berger, J. R., Ricote, S., Zhu, H., & Kee, R. J. (2021). Modeling Electro-Chemo-Mechanical Behaviors within the Dense BaZr0.8Y0.2O3−δ Protonic-Ceramic Membrane in a Long Tubular Electrochemical Cell. Membranes, 11(6), 378. https://doi.org/10.3390/membranes11060378

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