# Diffusion and Gas Conversion Analysis of Solid Oxide Fuel Cells at Loads via AC Impedance

Authors: Robert U. Payne; Ying Zhu; Wenhua H. Zhu; Mark S. Timper; S. Elangovan; Bruce J. Tatarchuk
Journal: International Journal of Electrochemistry
 (2011)

---

## Abstract

Impedance measurements were conducted under practical load conditions in solid oxide fuel cells of differing sizes. For a 2 cm2 button cell, impedance spectra data were separately measured for the anode, cathode, and total cell. Improved equivalent circuit models are proposed and applied to simulate each of measured impedance data. Circuit elements related to the chemical and physical processes have been added to the total-cell model to account for an extra relaxation process in the spectra not measured at either electrode. The processes to which elements are attributed have been deduced by varying cell temperature, load current, and hydrogen concentration. Spectra data were also obtained for a planar stack of five 61 cm2 cells and the individual cells therein, which were fitted to a simplified equivalent circuit model of the total button cell. Similar to the button cell, the planar cells and stack exhibit a pronounced low-frequency relaxation process, which has been attributed to concentration losses, that is, the combined effects of diffusion and gas conversion. The simplified total-cell model approximates well the dynamic behavior of the SOFC cells and the whole stack.


---

## Body

## 1. Introduction

Solid oxide fuel cells (SOFCs) have become a more attractive solution for power generation due to their higher energy efficiency, fuel flexibility, high-quality waste heat, tolerance to carbon monoxide, and lower/zero emissions    [1–3]. Effort has been made for fuel cell R&D and commercialization over the last several decades, with increasing concerns of fossil fuel rapid consumption and global warming from CO2 exhaust. For consideration of successful fuel cell commercialization, there are two major challenges, high cost and low reliability [4].  It is significantly important to have deep understanding of the electrode structures, catalytic mechanisms, reaction limitations, and failure mode diagnostics [5] in order to approach a technical breakthrough. AC impedance spectroscopy has been widely used for fuel cell characterization, reaction mechanism, and performance degradation analysis among those fuel cell diagnostic tools    [6–8].The dynamic load response of the device must be well understood in order to implement fuel cells into a power system. Impedance spectroscopy is a measurement technique capable of assessing the dynamic response characteristics of an electrochemical system [9] and is performed by superimposing an AC signal on the DC output of an electrochemical cell and measuring the impedance over a spectrum of frequencies.  These data are often fitted with an equivalent circuit model, which is a circuit of electrical elements producing a similar load response to the device under investigation [10]. These models are useful for designing the power conditioning system, through which the electrochemical device supplies power to a load.Equivalent circuits are also beneficial for stack diagnostics, because the different circuit elements theoretically represent different physical relaxation processes occurring in the electrochemical cell(s)   [11, 12]; hence, fit data can be used to distinguish which processes are limiting cell behavior and to what extent.  AC impedance is further useful for assessing system design and component fabrication, monitoring system aging and state of health, and performing quality control.Using electrochemical impedance spectroscopy, An et al. [13] measured the half-cell performance of the SOFCs and performed the durability analysis of the composite cathode after 1000 h of aging. Oestergard and Mogensen [14] conducted the kinetics study of the oxygen reduction reaction via AC impedance spectroscopy. Kato et al. [15] carried out impedance measurements under practical power generation conditions in a disk-type SOFC. Liu et al.    [16] fabricated a micro-SOFC stack, and it was operated under single-chamber conditions. Impedance spectra data of each cell were measured and analyzed under open-circuit conditions. Because of equipment limitations and the desire to maintain constant temperature and reactant concentration, most measurements of AC impedance data have been performed near open circuit voltage (OCV).  These investigations offer little insight into the behavior of a polarized cell.  The few studies that have tested a cell under load have overcome the power limitations of their equipment by testing electrodes much smaller than commercial ones.  An additional benefit of small cell studies is that cell conditions (reactant concentration, temperature, current, potential, etc.) are less distributed than in larger cells. This hinders the ability to scale up a small cell model to predict the behavior of a larger cell, but the small cell results are a beneficial reference for the larger SOFC work.  Lang et al.[17] applied electrochemical impedance spectroscopy at different operating gases and different electrical current loads to characterize solid oxide fuel cell short stacks for mobile applications. The nature of losses including ohmic and polarization resistances of the electrodes were examined and determined by fitting the impedance spectra to an equivalent circuit.When operating in a stack, individual cells are expected to behave differently due to manufacturing inconsistencies and flow, current, and temperature distribution.  Testing cells under the most realistic conditions possible is desirable for sound understanding of working systems.  This paper presents results and analysis of impedance spectra from solid oxide fuel cells at electrical current loads, the results of which are augmented by button cell studies.  The improved total-cell models are proposed for the SOFC cells and fitted to impedance spectra data obtained from the fuel cell tests.  Further analysis and discussion on the diffusion and gas conversion processes are conducted for the total SOFC cells and larger planar cells in the stack.


## 2. Experimental

Two different sizes of SOFC cells (button cells and larger planar cells, Ceramatec, Inc.) were operated and tested in the experimental work.  Both consisted of a zirconia-based electrolyte (ca. 170 μm thick), a nickel-based anode, and a manganite-based cathode.  The cells were fed at ca. 15 psig with air on the cathode side and mixtures of hydrogen and nitrogen bubbled through water on the anode side.  Cell temperature was maintained with furnaces controlled by a LabVIEW (National Instruments) program.  Furnaces were cooled down to 500°C during long intervals between tests (e.g., overnight) to minimize cell aging.The button cells consisted of an 8 cm2 disc of electrolyte onto either side of which 2 cm2 electrodes were brushed.  A small reference electrode of the same material as the corresponding electrode was brushed on both sides of the electrolyte near the electrode.  The SOFC cells were tested using a Solartron Analytical 1255 A frequency response analyzer in conjunction with a Schlumberger SI-1286 potentiostat/galvanostat.  The tests were controlled and data were recorded through a LabVIEW program.  Separate spectra were acquired galvanostatically for the anode, cathode, and total cell at a given current.  In each case, the system was perturbed with a sinusoidal wave of 5 mV amplitude over a spectrum of 1.0 MHz to 10 mHz.  The independent variables studied were cell current, cell temperature, and anode feed composition.  Gas mixtures were altered using rotameters upstream from the water bubbler maintaining a total flow of 39 SCCM.  Outlet gas compositions were measured with a molecular sieve column in an Agilent 3000 A gas chromatograph and analyzed with Agilent’s Cerity Network Data System software.  Inlet compositions were assumed to be equal to those of the outlet at open circuit.The planar stack consisted of five cells, each with an electrode area of 61 cm2.  The stack was fed with 1.5 SLPM air and 2.2 SLPM differing anode mixtures of hydrogen and nitrogen bubbled through water.  Anode inlet compositions were measured using the same Agilent system as before.  Stack temperature was regulated to be 800°C using a similar furnace and LabVIEW program as in the button cell tests.  Impedance spectra were measured with an FC350 potentiostat (Gamry Instruments) in conjunction with a TDI-Dynaload RBL488 programmable load.  Data were acquired at different stack currents for the cells on both ends of the stack (Cells 1 and 5) and the middle cell (Cell 3), as well as for the entire stack.  The sinusoidal waves were generated galvanostatically with currents producing a maximum potential amplitude of 5 mV (35–140 mA RMS depending on the operating DC current).  All spectra collected from either cell size were fitted to equivalent circuits with Gamry’s Echem Analyst software using the Levenberg-Marquardt algorithm.Impedance for the button cell was measured for the total cell as well as separately for anode and cathode with respect to a reference electrode at each test condition.  The equivalent circuit models from button cell results are properly combined, simplified, and then applied to those of the large cells and stack.  The fitted circuit elements are assigned to physical phenomena to determine which and to what extent physical processes limit stack power output.


## 3. Results and Discussion

Total-cell impedance data were analyzed in terms of equivalent circuits based on separately measured button cell anode and cathode spectra.  The association of different circuit elements with physical processes occurring within the cell was verified by varying cell conditions.  Total button cell model fits were compared to those of the individual electrodes, to identify differences and accordingly adjust the total-cell equivalent circuit configuration.  The button cell analyses were then conducted for reaction mechanisms and related physical processes.


### 3.1. Equivalent Circuit Models

SOFC cell models are based on button cell spectra, because their configuration is more conducive to impedance measurements.  The picayune nature of the button cells minimized temperature and current distribution over the active surface.  Furthermore, reference electrodes bonded to the button cell electrolyte allowed for the anode and cathode to be separately probed.  Individual electrode spectra were analyzed in terms of equivalent circuit models, which were then combined into a total-cell circuit.  The full button cell model was modified based on experimental data and simplified for large cell analysis.


#### 3.1.1. Half-Cell Measurements

Impedance spectra were measured separately for anode and cathode at different temperatures, DC load currents, and anode gas compositions.  Anode and cathode data were well fitted to a three- and a two-RC equivalent circuit model, respectively (Figures1 and 2).  The cathode fits were obtained using a constant phase element (CPE) in place of the Cd.  In-depth analysis of the fitted circuit elements is ill-advised, because studies [18, 19] have indicated that bonding a reference electrode to the same face of a solid electrolyte as the working electrode likely produces erroneous data.  Researchers also showed that the high- and low-frequency intercepts may be shifted and that an extraneous arc may exist at the frequency limits.  As such, specific values of the half-cell fits are still capable of indicating general trends in the simulation. This data simulation and analysis are helpful for providing insight into the findings of related anode/cathode reactions and transport processes in the whole SOFC cell.  Furthermore, only Rct and Cd in the circuits are assumed to be ascribed to real phenomenon (anode and cathode activation), whether the low-frequency elements in the models are caused by cell geometry or adsorbed intermediates.  Overall, the total-cell spectra data is reliable without the effect of shape, size, and location of the reference electrodes.Figure 1
Nyquist plot for a button cell anode operating at 800°C and 175 mA cm−2 and fed with 3% H2O and H2 diluted with varying amounts of N2. Data fits were obtained with the anode equivalent circuit model in the right-hand corner of the figure.Figure 2
Nyquist plot for a button cell cathode fed with air and operating at varying temperature and load current. The data fits were obtained with the cathode equivalent circuit model in the right-hand corner of the figure.  Note that a CPE is used to fit the data in place of theCd.


#### 3.1.2. Total-Cell Equivalent Circuit Model

The frequency-dependent impedance was measured between the anode and cathode under each of the same conditions as for the half-cell tests.  A full-cell equivalent circuit is typically that of the anode added in series with the cathode, wherein the ohmic contributions are combined into a single resistor.  Such an arrangement is theoretically well fitted to total-cell data assuming its spectrum was obtained at the same conditions as the individual electrode spectra; however, Figure3 depicts a large arc dominating the low-frequency data from the full-cell that is not present in the summed half cells.  The extraneous arc is not unique to the button cells and exists in the larger cell spectra (Figure 4) as well, albeit with an abscissa shifted to a higher frequency.Figure 3
Nyquist plot of the total cell compared with the sum of the individual anode and cathode spectra for a button cell operating at 750°C and 175 mA cm−2 and fed with 3% H2O and balance H2.  A large low-frequency arc is unique to the measured total-cell impedance.Figure 4
Nyquist plot of a planar fuel cell operating at 800°C and 52 mA cm−2 and fed with 3% H2O, 45% H2, and balance N2.  The cell has been fitted to two different equivalent circuits: one with three RCcircuits (simple fit) and the other with four (best fit).To aid the identification of the process responsible for the low-frequency impedance arc, the button cell data have been fitted to two different equivalent circuits.  Both models consist of the cathode and anode equivalent circuits used for the single-electrode fits in series and fixed to the previously obtained values.  A singleRΩ is in series with both circuits along with an L for lead inductance.  To fit the low-frequency loop, a parallel resistor and capacitor were added in series to one of the circuits (Figure 5(a)), while a finite diffusion element (FDE) was added to the other (Figure 5(b)).  The impedance of the FDE [20] is approximately described as (1)ZFDE=(1Y0jω)tanh⁡(Bjω),
where(2)Y0=(nF)2ARigT(∑1ciDi)-1,(3)B=δDi.
In theory, the location of the FDE in the equivalent circuit is related to the diffusing species: ionic diffusion is in series with the  RΩ, while H2 or O2 diffusion would be in series with the Faradaic impedance branch of the respective electrode.  Table 1(a) compares fitted values for the FDE in series with RΩ (as in Figure 5(b)) against those in series with  Rcta, which are markedly similar.  Fortunately, because of the insignificant change in fitted values with changing circuit topography, information can be inferred from the fits without knowing the low-frequency limiting process a priori.  In Section 3.3,the low-frequency impedance arc is further discussed in details.Table 1
Values for a parallel RC as shown in Figure5(a) are fitted to the low-frequency impedance loops of button cells.
(a)
Varying hydrogen feed concentration at two temperatures
Physical elementsRlfClfY0BY0*B
*NotexH(T= 750°C)(mΩ cm2)(F cm−2)(S cm−2 s1/2)(s1/2)(S cm−2s1/2)(s1/2)47.3%138.540.6025.664.49325.714.493i = 175 mA cm−263.7%137.441.4026.134.41226.424.43597.0%120.151.8531.254.70031.434.713xH(T = 800°C)33.1%172.358.4526.765.44127.535.532i = 175 mA cm−245.5%146.978.5535.566.09033.325.83197.0%123.688.3537.985.64340.645.079*The finite diffusion element (FDE), fitted in series with the Faradaic impedance as shown in Figure5(b), is also listed in this Table 1(a), showing little difference between the parallel RC and the series fits.(b)
Varied load current density at two hydrogen feed concentrations
Physical elementsRlfClfY0BNotei(mA cm−2, xH=33.1%)(mΩ cm2)(F cm−2)(S cm−2 s1/2)(s1/2)0.000638.816.997.2156.078T = 800°C, xw = 3.0%, Balance N2175.0172.358.4526.765.441300.0126.475.6535.285.313i(mA cm−2, xH=45.5%)0.000620.416.427.3105.849T = 800°C, xw = 3.0%, Balance N2175.0146.978.5535.566.090300.0100.0103.445.155.595Equivalent circuit model for a total button cell consisting of the button cell anode and cathode equivalent circuits in series with (a) a parallel RC (simple fit,Rlf and Clf) or (b) an FDE element (best fit).
(a)(b)Regardless of whether the extra semicircle is modeled with an FDE or a parallel RC, the total-cell circuits (Figures5(a) and 5(b)) must be simplified in order to fit the larger cells and stacks.  Fitting a circuit with so many elements of unknown value is rife with error, which is not an issue for the button cells, because the anode and cathode elements were set equal to the values obtained from the half-cell fits.  In simplifying the circuit, the minimum number of elements producing an RC circuit required to fit the data would be three: an FDE (or parallel low-frequency RC) and a parallel RC for the anode as well as the cathode.  Figure 4 shows that the three-RC circuit model is too simple needing a fourth RC circuit to fit the data more exactly.  The additional elements are ascribed to the anode (likely adsorption), because its effect on half-cell impedance is the most pronounced of the omitted elements (Figure 1); hence, it is in series with the Faradaic impedance branch of the anode in the large cell equivalent circuit model depicted in Figure 6.  The values of the fitted circuit elements are listed in Table 2.  In essence, each electrode model has been simplified for the total-cell fit by removing the pair of elements least influencing the shape of their respective half-cell curves, which is further supported by the dubious low-frequency behavior of the reference electrode configuration.  Now that the total-cell equivalent circuit model has been simplified, the validity of the model may be tested by examining the variation of the fitted values with respect to changing cell conditions.Table 2
Fitted equivalent circuit elements using the circuit model as shown in Figure6, varying anode fuel feeds for three individual cells in a five-cell stack, and operating at 800°C, 52.46 mA cm−2.  Feeds were 38.5 ~ 97.0% H2, 3.0% H2O, and balance N2.
Physical elementsxH (Cell 1)xH (Cell 3)xH (Cell 5)38.5%64.6%97.0%38.5%64.6%97.0%38.5%64.6%97.0%RΩ (mΩ cm2)131612871261673.4661.9657.6805.2792.4792.4Rctθ (mΩ cm2)300.2256.2244.9287.6235.0220.5226.5189.5178.2Cdc (mF cm−2)11.5015.1313.5611.8012.1614.0710.2810.8111.42Rcta (mΩ cm2)148.8174.8181.4115.8116.6141.981.6287.7897.23Cda (mF cm−2)5.6465.0894.6736.5646.4485.7896.3085.7305.815Rθa (mΩ cm2)†−118.0−69.17−39.12−78.32−74.42−60.45−68.44−74.60−61.67Cθa (F cm−2)†−1.984−3.739−5.705−1.481−2.356−2.726−1.277−1.455−2.449Y0 (S cm−2 s1/2)1.2041.5992.0511.2921.5541.8571.3511.5951.907B (s1/2)1.4451.5431.6681.3131.3941.4671.3371.3831.457χ2 (*10−6)‡63.744.326.917197.890.418212676.4†Negative values of Rθ and Cθ indicate that the inductive behavior of anode adsorbed intermediates is more significant than the capacitive portion.  The same fits can be obtained without negative values, when the parallel Rθ-Cθ is replaced with a series RL-L moved in parallel with the anode Rcta     [21].‡χ2 is the nonlinear least squares residual.Figure 6
Large cell equivalent circuit model.  The circuit has four RC circuits associated with the parallel RCs (RctaCda, RctcCdc, and RθaCθa) and the finite diffusion element (ZFDE).


### 3.2. High-Frequency Fits

Half-cell and total-cell impedance spectra are dominated by kinetics in the range of 10 kHz–1 Hz for the button cells, which are henceforth referred to as the high-frequency range.  The button cell spectra demonstrate the expected trends in the variance of high-frequency impedance.  Figure1 shows the Rcta increasing as the anode feed gas is diluted with inert gas (N2).  Cooling a cell from 800°C to 750°C also increases cell impedance as evident in Figure 2, which is partially why SOFCs must be run at such high temperatures.  The shape of loop changes to the normal Nyquist plot with increase of temperature.  Barbucci et al. [22] conducted the analysis of oxygen reduction on porous composite electrodes for better understanding of the mechanism of the SOFC cathode behavior. It gives similar Nyquist plots at high-frequency loop for the half cell with composite cathode at different temperatures.As for the variance ofRct with respect to current, the energy required to surmount the activation barrier is smaller at a larger overpotential; hence, Rct is expected to decrease as the overpotential (polarization) is increased.  Note that 97% hydrogen anode feeds are used to verify changes in cell polarization (Figure 7), because when inert gases are present equivalent inlet compositions are not comparable for varying cell current due to reactant consumption.  While button cell electrode configuration may cause errors in the half-cell data, the charge-transfer behavior trended as expected and is applied to test the veracity of the large cell fits.Figure 7
Nyquist plot for a total button cell fed with 97% H2 and balance H2O and operating at varying temperature and load current.Impedance data collected for the SOFC total cell are graphed in a Nyquist plot as shown in Figure4, and the experimental data have the same upward trend at high-frequency region as impedance data plotted by Lang et al. [17] for a one-cell short stack.  Data are well fitted to the total-cell model, but the mechanistic trends are not always correct in the charge-transfer dominated part of the SOFC cell impedance spectra.  One reason is that the impedance spectra for the large cells are dominated by the low-frequency curve (Figure 4).  Another source of fit error is in the relaxation times of the anode and cathode being similar enough to cause their respective curves to overlap.  Because the electrodes are indistinguishable, the variation in Rcta and Rctc does not necessarily match the theoretical trends confirmed in the half-cell measurements.  The Rct fits vary with changing H2 concentration and load current as shown in Table 1, respectively.  The value of Rct of the cathode (ca. 150 mΩ cm−2 at 800°C, 175 mA cm−2, and 97.0% H2 in Figure 2) is approximately 3 times more higher than that of the anode (ca. 45 mΩ cm−2 at 800°C, 175 mA cm−2, and 97.0% H2 in Figure 1).  The lower valued resistance (Rct at the anode) decreases with increasing xH, which is expected as shown in smaller loops in Figure 1.If the parallelRct-Cd is instead viewed as the combined kinetic fit, then the elements vary with changing conditions as expected.  In fact, a further simplified equivalent circuit (Figure 5(a)) wherein the electrode process elements are replaced by a single parallel Rct-CPE is well fitted to the large total-cell data (Figure 4).  As a whole, the Rct behavior nearly follows the trends as we expected in the simplified circuit.  The high-frequency fits of the total cell are useful for determining the overall kinetic performance of the cell.


### 3.3. Low-Frequency Fits

The low-frequency loop of the total-cell impedance spectra is related to reactant concentration at the interface, due to either mass transfer or gas conversion.  The FDE and parallel RC are used to model diffusion and gas conversion, respectively, and are fitted to the data equally well.  Because the process responsible for the low-frequency arc cannot be determined from the shape of the curve, the variation of the fitted values with respect to cell conditions must be investigated to determine whether diffusion or gas conversion is the cause.


#### 3.3.1. Mass Transfer

The impedance of mass transfer by diffusion is quantified by the fitted values of the FDE, which respond as expected to variations in temperature and reactant concentration.  Heating the cell imparts more thermal energy to the system facilitating mass transport and decreasing impedance (or increasing diffusivity).  Recalling the dependence of (1) and (2) on diffusivity verifies the change in the fitted Y0 and B with respect to temperature at the same cell current and reactant concentration (Table 1(a)).  Note that electrolytic diffusivity is more sensitive to temperature variations (Di∝eT [23]) than gas phase diffusivity (Di∝T1.75 [24]) thereby dominating the T in the denominator of  (2) and in gas concentration (ci∝T-1).  The relationship between Y0 and reactant concentration is more straightforward, increasing impedance as the feed is diluted with inert gas as explicitly shown in (1) and verified for both cell sizes in Tables 1 and 2.  The approximately constant value of B with respect to changing hydrogen concentration is also as expected, as it is a variable that should not appreciably alter Di or δ.The variation of the FDE with respect to current is as expected for the button cells, but not so for the larger ones.  Data fits in Table1 show that button cell Y0 increases with increasing current while B remains relatively constant, which is predicted by the FDE equations.  The only variables affected by a change in current are the concentration terms in the summation of (2) explaining why B remains the same.  The value of the summation is largely dominated by the reactant at lower concentration (normally water for the anode); hence, Y0 is large at high current, where cw is increased from its initial amount due to the reduction of hydrogen to water.  The size of the cell likely contributes to the difference between the button and larger cells, because reactant depletion at the cell outlet is more significant at higher currents and expected to have a more deleterious effect on larger cells.To explain the paucity of a similarly large low-frequency arc in half-cell impedance spectra, the potential of the reference electrodes must be considered.  If the concentrations of all reactants at the references are constant, then the potential difference between the two should be independent of the current through the cell.  The curves plotted in Figure8 indicate that concentration in fact varies at the reference surface.  The absence of the low-frequency arc from anode or cathode impedance spectra is then justified as existing at working and reference electrode to the same extent, thus canceling each other out in the single electrode measurements.Figure 8
Potential difference between anode and cathode reference electrodes plotted as a function of load current.Since reference concentrations have been shown to vary with respect to current, whether some or all diffusing species are responsible for the behavior should be considered.  The difficulty in gas phase diffusion occurring equally at the reference as at the active electrode is in the large low-frequency impedance at low current.  Intuitively, the high impedance measured in the low frequency of the total cell is not likely to exist to an equal extent at the reference as at the active electrode due to gas concentration, because the change in gas concentration due to consumption is negligible near OCV.  As for electrolytic mass transport, high mass transfer impedance at low current is also observed in lanthanum strontium manganese oxide (LSM) electrodes, which are known to be poor bulk ion conductors at low overpotentials.  Impedance spectra obtained for thin, dense LSM films on yttria stabilized zirconia    [17, 25] have abscissa in a similar frequency regime as depicted in the present study for the total cell as shown in Figure 4.  These arcs have been attributed to oxide (or vacancy) diffusion through the film, as confirmed by the shift of the abscissa frequency at differing film thicknesses.  Even though this Ceramatec cathode is not composed of LSM, the properties of other manganite-based materials appear to be similar in the equivalent circuit model.Electrolytic mass transfer is further supported by comparing the fittedB value between button and larger cells.  For the specific cells in this study, the thicknesses of the mixed conducting layer of the electrodes were significantly different between the button cells and the large cells.  The total thickness (including mixed conducting components) for the large cell anode and cathode is approximately 52–55 μm and 77–80 μm, respectively, compared to 307–358 μm and 347–398 μm for the button cell.  The effect that thinner electrodes impose on mass transfer is visually evident in the low-frequency semicircle being fully resolved for the large cell impedance compared to that of the button cells over the same frequency spectrum (Figures 3 and 4).  The shift of the abscissa frequency between the two from approximately 10 mHz to 200 mHz is quantified by the B parameter of the FDE, which changes due to the diffusion layer thickness.  The ionic/interstitial oxygen (or vacancy) diffusion layer would be the sum of the thicknesses of the electrolyte and electrodes, assuming the majority of reactions occur on the outer edges of the electrodes.  The ratio of the button cell B to that of the large cell is approximately three at OCV, which is quite close to the ratio of the cathode thicknesses (2.76–3.04).  Comparing higher current levels is likely not valid, because of temperature changes due to reaction heat generation in the larger cells.


#### 3.3.2. Gas Conversion Impedance

Primdahl and Mogensen [26] observed a low-frequency impedance arc related to gas conversion, when the reference electrode was located in a different atmosphere from the working electrode.  They proposed that the passage of AC current changes gas concentration at the active surface, most noticeably at low frequencies, resulting in a Nernstian potential loss.  Such is not observed when the working and reference electrodes are at the same concentration, which could be the cause of the previously discussed reference i-V behavior (Section 3.3.1) in Figure 8.  Their study modeled the gas conversion impedance (GCI) with a parallel Rlf and Clf in series with the electrode equivalent circuit (similar to the equivalent circuit described in Section 3.1.2, Figure 5(a)) and approximately replicated the fitted values with a constantly stirred tank reactor (CSTR) model yielding the following relationships:(4)Rlf=RigT(2F)2J(1xw+1xH),(5)Clf=4F2PV(RigT)2A(1xw+1xH)-1
wherein the resistance and capacitance are area specific.  Applying Primdahl’s model to the cells in question indicates that gas conversion may be partially responsible for the extra arc.While a 61 cm2 planar cell is expected to deviate from a single tank CSTR model due to the previously discussed distributed conditions, the general variation of impedance with each of the independent variables should be unchanged.  Increases in cell current or hydrogen concentration predict that Rlf should decrease and Clf should increase, a tendency observed in the button cell fits in Tables 1(a) and 1(b).  The same impedance effect should be observed for decreased temperature; however, Clf fits decrease by no less than 44% for a temperature decrement of 50°C.  Temperature was not explicitly controlled for the larger cells to test whether the same anomaly is present, but hydrogen concentration (Table 3(a)) and cell current (Table 3(b)) were varied with only the former behaving as in the button cell tests.  Specifically, Clf is only slightly affected by increasing stack current from 1.639 mA cm−2 to 52.46 mA cm−2, although the model predictions and button cell data indicate a large increase should be observed.  This effect cannot be justified by the expected increase in cell temperature induced by increasing current, because the GCI model indicates an increment of ≫100°C would be necessary for such a drastic shift to be observed.  Notwithstanding the fact that the button cell fitted Clf exhibit opposite behavior (Table 1(a)) with respect to the model predictions (5) at increased temperature.  The opposite trend is not clear and further intensive work is necessary to examine the effect of temperature on the gas conversion impedance, especially on the value of Clf capacitance.Table 3
Values for a parallel RC fitted to the low-frequency impedance loops of three cells in a five-cell planar stack operating at 800°C.
(a)
The SOFCs were fed with 38.5 ~ 97% H2, 3% H2O, and balance N2, and operated at 1.69 mA cm−2
Physical elementsxH (Cell 1)xH (Cell 3)xH (Cell 5)38.5%64.6%97.0%38.5%64.6%97.0%38.5%64.6%97.0%Rlf(Ω cm2)1.4301.5381.3731.3401.4791.2631.2991.4271.301Clf(F cm−2)1.1851.5861.6580.94771.3421.4650.98831.2381.447Y0 (S cm−2 s1/2)1.1351.3321.4411.0181.2111.3971.0591.2021.338B(s1/2)2.0252.4872.4021.7272.2102.1551.7432.0982.149(b)
Varying load currents when fed with 97.0% H2 and balance H2O
Physical elementsi (mA cm−2, Cell 1)i (mA cm−2, Cell 3)i (mA cm−2, Cell 5)1.63952.46147.51.63952.46147.51.63952.46147.5Rlf (Ω cm2)1.3730.65700.58021.2630.62830.48621.3010.60820.4373Clf (F cm−2)1.6591.7161.6481.4651.4171.3801.4471.4451.895Y0 (S cm−2 s1/2)1.4422.0511.9251.3971.8571.8211.3371.9071.731B(s1/2)2.4021.6681.5542.1561.4671.3762.1481.4571.259(c)
Varying the anode feed ratej when the SOFCs were operated at 52.46 mA cm−2 and composed of 97.0% H2 and balance H2O, where the total j is assumed to be evenly distributed between the five cells in the stack.
Physical elementsj (μmol s−1 cm−2, Cell 1)j (μmol s−1 cm−2, Cell 3)j (μmol s−1 cm−2, Cell 5)3.1684.9793.1684.9793.1684.979Rlf (Ω cm2)0.86320.65700.81310.62830.77590.6082Clf (F cm−2)1.9101.7161.6481.4171.7341.446Y0 (S cm−2 s1/2)1.9562.0511.9001.8571.9521.907B(s1/2)2.0491.6681.8631.4671.8431.457The best evidence that gas conversion is likely responsible for the low-frequency impedance loop is listed in Table3(c), wherein an increase in the total anode gas flow decreases Rlf and decreases Clf  The relationship between gas conversion impedance and anode feed rate j has been demonstrated experimentally by Momma et al. [27], who applied it to determine the flow distribution through a 46-cell stack.  Whether the low-frequency arc is due entirely to gas conversion is not conclusive enough to apply the same method to test the flow distribution in the Ceramatec stack.Neither gas conversion nor mass transfer is likely to exclusively cause the low-frequency impedance loop.  The two can scarcely be distinguished from the appearance of the arc, because a parallel RC and FDE are similarly shaped on a Nyquist plot [28] as evident in the χ2 values of the fits.  The difference is in the Warburg-like part of the FDE, which appears as a 45° line on a Nyquist plot to the higher frequency side of the FDE abscissa [20].  Anode adsorption occurs in the same frequency regime as the Warburg part of the FDE, which obscures the shape of the low-frequency loop.  Exacerbating the choice of elements is the functionality of  (2) and (4), which indicates that Rlf and Y0 often vary inversely; however, such a relationship is not always observed in the values listed in Table 3.  The correlation between fitted values and changes in hydrogen concentration or temperature is more congruent with mass transfer, while gas conversion better explains that of gas flow rate.  The effects of gas conversion and mass transfer are most likely present concurrently being superimposed as an arc in the low-frequency region of a Nyquist plot.The utility of applying equivalent circuit models for mechanistically discriminating problems with individual cell performance is evident in analyzing the different large cells.  In the stack under investigation, Cell 3 is in the middle of the stack flanked by Cells  1 and 5 on the outsides.  The two outside cells are not identical, with Cell 5 having a thermally insulating layer adjacent to the outer wall that is unique to it.  Because the temperature of Cell 1 is lower than the other cells in the stack, it lacks the thermal energy necessary to drive the reaction as fast, thus limiting performance in a stack of cells connected in series.  One need not know the geometry of the system to deduce that Cell 1 is at a lower temperature than the others; such can be implied from the model fits.  While the increasedRct for the cell in question may indicate a decreased temperature, the trend may also be induced by inconsistent electrode microstructure, catalyst loading, and/or reactant concentration.  Coupling the Rct behavior with the increase in RΩ and B for Cell 1 relative to the other cells proves the root of the problem, because decreased electrolyte conduction and diffusivity are also indicative of lower temperatures.  The analyses of the fitted large cell data show how impedance spectroscopy can be used to mechanistically discriminate the factors that limit stack performance.


### 3.4. Fuel Cell Stack

Data for the entire fuel cell stack were also obtained under the same test conditions as the single cell tests previously described in Sections3.2 and 3.3.  Figure 9 depicts some of the obtained curves, each of which is similar in shape to the individual cells; hence, the same equivalent circuit that was fitted to the large cell data has been similarly applied to the stack data.  The key to the success of applying an individual cell model to a stack of cells is the distribution of  RC circuits of each physical relaxation process from cell to cell.  Although not precisely equal between the five cells, the frequency ranges at which each physical process dominates cell impedance are similar enough for the arcs to superimpose.  Had the degree of distribution been greater, the semicircles of the Nyquist plot would have appeared to be depressed (CPE-like); instead, the individual cell model is well fitted to the stack impedance data, as exhibited in this figure.Figure 9
Nyquist plot of a planar fuel cell stack fed with 97% H2 and balance H2O operating at 800°C and varying load current.  The data have been fitted to a four-RC equivalent circuit.While the individual cell tests were designed to determine which processes limit cell output and to what degree, the stack test is primarily intended to obtain a simple, purely electronic model to replicate the dynamic response of the stack.  In order to represent the stack the purely physical electronic elements, the FDE of the equivalent circuit (Figure6) has been replaced with a parallel RC, which has already been shown to fit single cell data equally well.  Obtaining fits for the circuit elements over the expected range of operation allows for the stack to be dynamically modeled over said range.  The utility in eliminating all distributed elements from the model is to facilitate circuit simulation [12], which can be used to design process controls and to assess the optimum configuration of operating conditions, power conditioning systems, and combination of multiple stacks.


## 3.1. Equivalent Circuit Models

SOFC cell models are based on button cell spectra, because their configuration is more conducive to impedance measurements.  The picayune nature of the button cells minimized temperature and current distribution over the active surface.  Furthermore, reference electrodes bonded to the button cell electrolyte allowed for the anode and cathode to be separately probed.  Individual electrode spectra were analyzed in terms of equivalent circuit models, which were then combined into a total-cell circuit.  The full button cell model was modified based on experimental data and simplified for large cell analysis.


### 3.1.1. Half-Cell Measurements

Impedance spectra were measured separately for anode and cathode at different temperatures, DC load currents, and anode gas compositions.  Anode and cathode data were well fitted to a three- and a two-RC equivalent circuit model, respectively (Figures1 and 2).  The cathode fits were obtained using a constant phase element (CPE) in place of the Cd.  In-depth analysis of the fitted circuit elements is ill-advised, because studies [18, 19] have indicated that bonding a reference electrode to the same face of a solid electrolyte as the working electrode likely produces erroneous data.  Researchers also showed that the high- and low-frequency intercepts may be shifted and that an extraneous arc may exist at the frequency limits.  As such, specific values of the half-cell fits are still capable of indicating general trends in the simulation. This data simulation and analysis are helpful for providing insight into the findings of related anode/cathode reactions and transport processes in the whole SOFC cell.  Furthermore, only Rct and Cd in the circuits are assumed to be ascribed to real phenomenon (anode and cathode activation), whether the low-frequency elements in the models are caused by cell geometry or adsorbed intermediates.  Overall, the total-cell spectra data is reliable without the effect of shape, size, and location of the reference electrodes.Figure 1
Nyquist plot for a button cell anode operating at 800°C and 175 mA cm−2 and fed with 3% H2O and H2 diluted with varying amounts of N2. Data fits were obtained with the anode equivalent circuit model in the right-hand corner of the figure.Figure 2
Nyquist plot for a button cell cathode fed with air and operating at varying temperature and load current. The data fits were obtained with the cathode equivalent circuit model in the right-hand corner of the figure.  Note that a CPE is used to fit the data in place of theCd.


### 3.1.2. Total-Cell Equivalent Circuit Model

The frequency-dependent impedance was measured between the anode and cathode under each of the same conditions as for the half-cell tests.  A full-cell equivalent circuit is typically that of the anode added in series with the cathode, wherein the ohmic contributions are combined into a single resistor.  Such an arrangement is theoretically well fitted to total-cell data assuming its spectrum was obtained at the same conditions as the individual electrode spectra; however, Figure3 depicts a large arc dominating the low-frequency data from the full-cell that is not present in the summed half cells.  The extraneous arc is not unique to the button cells and exists in the larger cell spectra (Figure 4) as well, albeit with an abscissa shifted to a higher frequency.Figure 3
Nyquist plot of the total cell compared with the sum of the individual anode and cathode spectra for a button cell operating at 750°C and 175 mA cm−2 and fed with 3% H2O and balance H2.  A large low-frequency arc is unique to the measured total-cell impedance.Figure 4
Nyquist plot of a planar fuel cell operating at 800°C and 52 mA cm−2 and fed with 3% H2O, 45% H2, and balance N2.  The cell has been fitted to two different equivalent circuits: one with three RCcircuits (simple fit) and the other with four (best fit).To aid the identification of the process responsible for the low-frequency impedance arc, the button cell data have been fitted to two different equivalent circuits.  Both models consist of the cathode and anode equivalent circuits used for the single-electrode fits in series and fixed to the previously obtained values.  A singleRΩ is in series with both circuits along with an L for lead inductance.  To fit the low-frequency loop, a parallel resistor and capacitor were added in series to one of the circuits (Figure 5(a)), while a finite diffusion element (FDE) was added to the other (Figure 5(b)).  The impedance of the FDE [20] is approximately described as (1)ZFDE=(1Y0jω)tanh⁡(Bjω),
where(2)Y0=(nF)2ARigT(∑1ciDi)-1,(3)B=δDi.
In theory, the location of the FDE in the equivalent circuit is related to the diffusing species: ionic diffusion is in series with the  RΩ, while H2 or O2 diffusion would be in series with the Faradaic impedance branch of the respective electrode.  Table 1(a) compares fitted values for the FDE in series with RΩ (as in Figure 5(b)) against those in series with  Rcta, which are markedly similar.  Fortunately, because of the insignificant change in fitted values with changing circuit topography, information can be inferred from the fits without knowing the low-frequency limiting process a priori.  In Section 3.3,the low-frequency impedance arc is further discussed in details.Table 1
Values for a parallel RC as shown in Figure5(a) are fitted to the low-frequency impedance loops of button cells.
(a)
Varying hydrogen feed concentration at two temperatures
Physical elementsRlfClfY0BY0*B
*NotexH(T= 750°C)(mΩ cm2)(F cm−2)(S cm−2 s1/2)(s1/2)(S cm−2s1/2)(s1/2)47.3%138.540.6025.664.49325.714.493i = 175 mA cm−263.7%137.441.4026.134.41226.424.43597.0%120.151.8531.254.70031.434.713xH(T = 800°C)33.1%172.358.4526.765.44127.535.532i = 175 mA cm−245.5%146.978.5535.566.09033.325.83197.0%123.688.3537.985.64340.645.079*The finite diffusion element (FDE), fitted in series with the Faradaic impedance as shown in Figure5(b), is also listed in this Table 1(a), showing little difference between the parallel RC and the series fits.(b)
Varied load current density at two hydrogen feed concentrations
Physical elementsRlfClfY0BNotei(mA cm−2, xH=33.1%)(mΩ cm2)(F cm−2)(S cm−2 s1/2)(s1/2)0.000638.816.997.2156.078T = 800°C, xw = 3.0%, Balance N2175.0172.358.4526.765.441300.0126.475.6535.285.313i(mA cm−2, xH=45.5%)0.000620.416.427.3105.849T = 800°C, xw = 3.0%, Balance N2175.0146.978.5535.566.090300.0100.0103.445.155.595Equivalent circuit model for a total button cell consisting of the button cell anode and cathode equivalent circuits in series with (a) a parallel RC (simple fit,Rlf and Clf) or (b) an FDE element (best fit).
(a)(b)Regardless of whether the extra semicircle is modeled with an FDE or a parallel RC, the total-cell circuits (Figures5(a) and 5(b)) must be simplified in order to fit the larger cells and stacks.  Fitting a circuit with so many elements of unknown value is rife with error, which is not an issue for the button cells, because the anode and cathode elements were set equal to the values obtained from the half-cell fits.  In simplifying the circuit, the minimum number of elements producing an RC circuit required to fit the data would be three: an FDE (or parallel low-frequency RC) and a parallel RC for the anode as well as the cathode.  Figure 4 shows that the three-RC circuit model is too simple needing a fourth RC circuit to fit the data more exactly.  The additional elements are ascribed to the anode (likely adsorption), because its effect on half-cell impedance is the most pronounced of the omitted elements (Figure 1); hence, it is in series with the Faradaic impedance branch of the anode in the large cell equivalent circuit model depicted in Figure 6.  The values of the fitted circuit elements are listed in Table 2.  In essence, each electrode model has been simplified for the total-cell fit by removing the pair of elements least influencing the shape of their respective half-cell curves, which is further supported by the dubious low-frequency behavior of the reference electrode configuration.  Now that the total-cell equivalent circuit model has been simplified, the validity of the model may be tested by examining the variation of the fitted values with respect to changing cell conditions.Table 2
Fitted equivalent circuit elements using the circuit model as shown in Figure6, varying anode fuel feeds for three individual cells in a five-cell stack, and operating at 800°C, 52.46 mA cm−2.  Feeds were 38.5 ~ 97.0% H2, 3.0% H2O, and balance N2.
Physical elementsxH (Cell 1)xH (Cell 3)xH (Cell 5)38.5%64.6%97.0%38.5%64.6%97.0%38.5%64.6%97.0%RΩ (mΩ cm2)131612871261673.4661.9657.6805.2792.4792.4Rctθ (mΩ cm2)300.2256.2244.9287.6235.0220.5226.5189.5178.2Cdc (mF cm−2)11.5015.1313.5611.8012.1614.0710.2810.8111.42Rcta (mΩ cm2)148.8174.8181.4115.8116.6141.981.6287.7897.23Cda (mF cm−2)5.6465.0894.6736.5646.4485.7896.3085.7305.815Rθa (mΩ cm2)†−118.0−69.17−39.12−78.32−74.42−60.45−68.44−74.60−61.67Cθa (F cm−2)†−1.984−3.739−5.705−1.481−2.356−2.726−1.277−1.455−2.449Y0 (S cm−2 s1/2)1.2041.5992.0511.2921.5541.8571.3511.5951.907B (s1/2)1.4451.5431.6681.3131.3941.4671.3371.3831.457χ2 (*10−6)‡63.744.326.917197.890.418212676.4†Negative values of Rθ and Cθ indicate that the inductive behavior of anode adsorbed intermediates is more significant than the capacitive portion.  The same fits can be obtained without negative values, when the parallel Rθ-Cθ is replaced with a series RL-L moved in parallel with the anode Rcta     [21].‡χ2 is the nonlinear least squares residual.Figure 6
Large cell equivalent circuit model.  The circuit has four RC circuits associated with the parallel RCs (RctaCda, RctcCdc, and RθaCθa) and the finite diffusion element (ZFDE).


## 3.1.1. Half-Cell Measurements

Impedance spectra were measured separately for anode and cathode at different temperatures, DC load currents, and anode gas compositions.  Anode and cathode data were well fitted to a three- and a two-RC equivalent circuit model, respectively (Figures1 and 2).  The cathode fits were obtained using a constant phase element (CPE) in place of the Cd.  In-depth analysis of the fitted circuit elements is ill-advised, because studies [18, 19] have indicated that bonding a reference electrode to the same face of a solid electrolyte as the working electrode likely produces erroneous data.  Researchers also showed that the high- and low-frequency intercepts may be shifted and that an extraneous arc may exist at the frequency limits.  As such, specific values of the half-cell fits are still capable of indicating general trends in the simulation. This data simulation and analysis are helpful for providing insight into the findings of related anode/cathode reactions and transport processes in the whole SOFC cell.  Furthermore, only Rct and Cd in the circuits are assumed to be ascribed to real phenomenon (anode and cathode activation), whether the low-frequency elements in the models are caused by cell geometry or adsorbed intermediates.  Overall, the total-cell spectra data is reliable without the effect of shape, size, and location of the reference electrodes.Figure 1
Nyquist plot for a button cell anode operating at 800°C and 175 mA cm−2 and fed with 3% H2O and H2 diluted with varying amounts of N2. Data fits were obtained with the anode equivalent circuit model in the right-hand corner of the figure.Figure 2
Nyquist plot for a button cell cathode fed with air and operating at varying temperature and load current. The data fits were obtained with the cathode equivalent circuit model in the right-hand corner of the figure.  Note that a CPE is used to fit the data in place of theCd.


## 3.1.2. Total-Cell Equivalent Circuit Model

The frequency-dependent impedance was measured between the anode and cathode under each of the same conditions as for the half-cell tests.  A full-cell equivalent circuit is typically that of the anode added in series with the cathode, wherein the ohmic contributions are combined into a single resistor.  Such an arrangement is theoretically well fitted to total-cell data assuming its spectrum was obtained at the same conditions as the individual electrode spectra; however, Figure3 depicts a large arc dominating the low-frequency data from the full-cell that is not present in the summed half cells.  The extraneous arc is not unique to the button cells and exists in the larger cell spectra (Figure 4) as well, albeit with an abscissa shifted to a higher frequency.Figure 3
Nyquist plot of the total cell compared with the sum of the individual anode and cathode spectra for a button cell operating at 750°C and 175 mA cm−2 and fed with 3% H2O and balance H2.  A large low-frequency arc is unique to the measured total-cell impedance.Figure 4
Nyquist plot of a planar fuel cell operating at 800°C and 52 mA cm−2 and fed with 3% H2O, 45% H2, and balance N2.  The cell has been fitted to two different equivalent circuits: one with three RCcircuits (simple fit) and the other with four (best fit).To aid the identification of the process responsible for the low-frequency impedance arc, the button cell data have been fitted to two different equivalent circuits.  Both models consist of the cathode and anode equivalent circuits used for the single-electrode fits in series and fixed to the previously obtained values.  A singleRΩ is in series with both circuits along with an L for lead inductance.  To fit the low-frequency loop, a parallel resistor and capacitor were added in series to one of the circuits (Figure 5(a)), while a finite diffusion element (FDE) was added to the other (Figure 5(b)).  The impedance of the FDE [20] is approximately described as (1)ZFDE=(1Y0jω)tanh⁡(Bjω),
where(2)Y0=(nF)2ARigT(∑1ciDi)-1,(3)B=δDi.
In theory, the location of the FDE in the equivalent circuit is related to the diffusing species: ionic diffusion is in series with the  RΩ, while H2 or O2 diffusion would be in series with the Faradaic impedance branch of the respective electrode.  Table 1(a) compares fitted values for the FDE in series with RΩ (as in Figure 5(b)) against those in series with  Rcta, which are markedly similar.  Fortunately, because of the insignificant change in fitted values with changing circuit topography, information can be inferred from the fits without knowing the low-frequency limiting process a priori.  In Section 3.3,the low-frequency impedance arc is further discussed in details.Table 1
Values for a parallel RC as shown in Figure5(a) are fitted to the low-frequency impedance loops of button cells.
(a)
Varying hydrogen feed concentration at two temperatures
Physical elementsRlfClfY0BY0*B
*NotexH(T= 750°C)(mΩ cm2)(F cm−2)(S cm−2 s1/2)(s1/2)(S cm−2s1/2)(s1/2)47.3%138.540.6025.664.49325.714.493i = 175 mA cm−263.7%137.441.4026.134.41226.424.43597.0%120.151.8531.254.70031.434.713xH(T = 800°C)33.1%172.358.4526.765.44127.535.532i = 175 mA cm−245.5%146.978.5535.566.09033.325.83197.0%123.688.3537.985.64340.645.079*The finite diffusion element (FDE), fitted in series with the Faradaic impedance as shown in Figure5(b), is also listed in this Table 1(a), showing little difference between the parallel RC and the series fits.(b)
Varied load current density at two hydrogen feed concentrations
Physical elementsRlfClfY0BNotei(mA cm−2, xH=33.1%)(mΩ cm2)(F cm−2)(S cm−2 s1/2)(s1/2)0.000638.816.997.2156.078T = 800°C, xw = 3.0%, Balance N2175.0172.358.4526.765.441300.0126.475.6535.285.313i(mA cm−2, xH=45.5%)0.000620.416.427.3105.849T = 800°C, xw = 3.0%, Balance N2175.0146.978.5535.566.090300.0100.0103.445.155.595Equivalent circuit model for a total button cell consisting of the button cell anode and cathode equivalent circuits in series with (a) a parallel RC (simple fit,Rlf and Clf) or (b) an FDE element (best fit).
(a)(b)Regardless of whether the extra semicircle is modeled with an FDE or a parallel RC, the total-cell circuits (Figures5(a) and 5(b)) must be simplified in order to fit the larger cells and stacks.  Fitting a circuit with so many elements of unknown value is rife with error, which is not an issue for the button cells, because the anode and cathode elements were set equal to the values obtained from the half-cell fits.  In simplifying the circuit, the minimum number of elements producing an RC circuit required to fit the data would be three: an FDE (or parallel low-frequency RC) and a parallel RC for the anode as well as the cathode.  Figure 4 shows that the three-RC circuit model is too simple needing a fourth RC circuit to fit the data more exactly.  The additional elements are ascribed to the anode (likely adsorption), because its effect on half-cell impedance is the most pronounced of the omitted elements (Figure 1); hence, it is in series with the Faradaic impedance branch of the anode in the large cell equivalent circuit model depicted in Figure 6.  The values of the fitted circuit elements are listed in Table 2.  In essence, each electrode model has been simplified for the total-cell fit by removing the pair of elements least influencing the shape of their respective half-cell curves, which is further supported by the dubious low-frequency behavior of the reference electrode configuration.  Now that the total-cell equivalent circuit model has been simplified, the validity of the model may be tested by examining the variation of the fitted values with respect to changing cell conditions.Table 2
Fitted equivalent circuit elements using the circuit model as shown in Figure6, varying anode fuel feeds for three individual cells in a five-cell stack, and operating at 800°C, 52.46 mA cm−2.  Feeds were 38.5 ~ 97.0% H2, 3.0% H2O, and balance N2.
Physical elementsxH (Cell 1)xH (Cell 3)xH (Cell 5)38.5%64.6%97.0%38.5%64.6%97.0%38.5%64.6%97.0%RΩ (mΩ cm2)131612871261673.4661.9657.6805.2792.4792.4Rctθ (mΩ cm2)300.2256.2244.9287.6235.0220.5226.5189.5178.2Cdc (mF cm−2)11.5015.1313.5611.8012.1614.0710.2810.8111.42Rcta (mΩ cm2)148.8174.8181.4115.8116.6141.981.6287.7897.23Cda (mF cm−2)5.6465.0894.6736.5646.4485.7896.3085.7305.815Rθa (mΩ cm2)†−118.0−69.17−39.12−78.32−74.42−60.45−68.44−74.60−61.67Cθa (F cm−2)†−1.984−3.739−5.705−1.481−2.356−2.726−1.277−1.455−2.449Y0 (S cm−2 s1/2)1.2041.5992.0511.2921.5541.8571.3511.5951.907B (s1/2)1.4451.5431.6681.3131.3941.4671.3371.3831.457χ2 (*10−6)‡63.744.326.917197.890.418212676.4†Negative values of Rθ and Cθ indicate that the inductive behavior of anode adsorbed intermediates is more significant than the capacitive portion.  The same fits can be obtained without negative values, when the parallel Rθ-Cθ is replaced with a series RL-L moved in parallel with the anode Rcta     [21].‡χ2 is the nonlinear least squares residual.Figure 6
Large cell equivalent circuit model.  The circuit has four RC circuits associated with the parallel RCs (RctaCda, RctcCdc, and RθaCθa) and the finite diffusion element (ZFDE).


## 3.2. High-Frequency Fits

Half-cell and total-cell impedance spectra are dominated by kinetics in the range of 10 kHz–1 Hz for the button cells, which are henceforth referred to as the high-frequency range.  The button cell spectra demonstrate the expected trends in the variance of high-frequency impedance.  Figure1 shows the Rcta increasing as the anode feed gas is diluted with inert gas (N2).  Cooling a cell from 800°C to 750°C also increases cell impedance as evident in Figure 2, which is partially why SOFCs must be run at such high temperatures.  The shape of loop changes to the normal Nyquist plot with increase of temperature.  Barbucci et al. [22] conducted the analysis of oxygen reduction on porous composite electrodes for better understanding of the mechanism of the SOFC cathode behavior. It gives similar Nyquist plots at high-frequency loop for the half cell with composite cathode at different temperatures.As for the variance ofRct with respect to current, the energy required to surmount the activation barrier is smaller at a larger overpotential; hence, Rct is expected to decrease as the overpotential (polarization) is increased.  Note that 97% hydrogen anode feeds are used to verify changes in cell polarization (Figure 7), because when inert gases are present equivalent inlet compositions are not comparable for varying cell current due to reactant consumption.  While button cell electrode configuration may cause errors in the half-cell data, the charge-transfer behavior trended as expected and is applied to test the veracity of the large cell fits.Figure 7
Nyquist plot for a total button cell fed with 97% H2 and balance H2O and operating at varying temperature and load current.Impedance data collected for the SOFC total cell are graphed in a Nyquist plot as shown in Figure4, and the experimental data have the same upward trend at high-frequency region as impedance data plotted by Lang et al. [17] for a one-cell short stack.  Data are well fitted to the total-cell model, but the mechanistic trends are not always correct in the charge-transfer dominated part of the SOFC cell impedance spectra.  One reason is that the impedance spectra for the large cells are dominated by the low-frequency curve (Figure 4).  Another source of fit error is in the relaxation times of the anode and cathode being similar enough to cause their respective curves to overlap.  Because the electrodes are indistinguishable, the variation in Rcta and Rctc does not necessarily match the theoretical trends confirmed in the half-cell measurements.  The Rct fits vary with changing H2 concentration and load current as shown in Table 1, respectively.  The value of Rct of the cathode (ca. 150 mΩ cm−2 at 800°C, 175 mA cm−2, and 97.0% H2 in Figure 2) is approximately 3 times more higher than that of the anode (ca. 45 mΩ cm−2 at 800°C, 175 mA cm−2, and 97.0% H2 in Figure 1).  The lower valued resistance (Rct at the anode) decreases with increasing xH, which is expected as shown in smaller loops in Figure 1.If the parallelRct-Cd is instead viewed as the combined kinetic fit, then the elements vary with changing conditions as expected.  In fact, a further simplified equivalent circuit (Figure 5(a)) wherein the electrode process elements are replaced by a single parallel Rct-CPE is well fitted to the large total-cell data (Figure 4).  As a whole, the Rct behavior nearly follows the trends as we expected in the simplified circuit.  The high-frequency fits of the total cell are useful for determining the overall kinetic performance of the cell.


## 3.3. Low-Frequency Fits

The low-frequency loop of the total-cell impedance spectra is related to reactant concentration at the interface, due to either mass transfer or gas conversion.  The FDE and parallel RC are used to model diffusion and gas conversion, respectively, and are fitted to the data equally well.  Because the process responsible for the low-frequency arc cannot be determined from the shape of the curve, the variation of the fitted values with respect to cell conditions must be investigated to determine whether diffusion or gas conversion is the cause.


### 3.3.1. Mass Transfer

The impedance of mass transfer by diffusion is quantified by the fitted values of the FDE, which respond as expected to variations in temperature and reactant concentration.  Heating the cell imparts more thermal energy to the system facilitating mass transport and decreasing impedance (or increasing diffusivity).  Recalling the dependence of (1) and (2) on diffusivity verifies the change in the fitted Y0 and B with respect to temperature at the same cell current and reactant concentration (Table 1(a)).  Note that electrolytic diffusivity is more sensitive to temperature variations (Di∝eT [23]) than gas phase diffusivity (Di∝T1.75 [24]) thereby dominating the T in the denominator of  (2) and in gas concentration (ci∝T-1).  The relationship between Y0 and reactant concentration is more straightforward, increasing impedance as the feed is diluted with inert gas as explicitly shown in (1) and verified for both cell sizes in Tables 1 and 2.  The approximately constant value of B with respect to changing hydrogen concentration is also as expected, as it is a variable that should not appreciably alter Di or δ.The variation of the FDE with respect to current is as expected for the button cells, but not so for the larger ones.  Data fits in Table1 show that button cell Y0 increases with increasing current while B remains relatively constant, which is predicted by the FDE equations.  The only variables affected by a change in current are the concentration terms in the summation of (2) explaining why B remains the same.  The value of the summation is largely dominated by the reactant at lower concentration (normally water for the anode); hence, Y0 is large at high current, where cw is increased from its initial amount due to the reduction of hydrogen to water.  The size of the cell likely contributes to the difference between the button and larger cells, because reactant depletion at the cell outlet is more significant at higher currents and expected to have a more deleterious effect on larger cells.To explain the paucity of a similarly large low-frequency arc in half-cell impedance spectra, the potential of the reference electrodes must be considered.  If the concentrations of all reactants at the references are constant, then the potential difference between the two should be independent of the current through the cell.  The curves plotted in Figure8 indicate that concentration in fact varies at the reference surface.  The absence of the low-frequency arc from anode or cathode impedance spectra is then justified as existing at working and reference electrode to the same extent, thus canceling each other out in the single electrode measurements.Figure 8
Potential difference between anode and cathode reference electrodes plotted as a function of load current.Since reference concentrations have been shown to vary with respect to current, whether some or all diffusing species are responsible for the behavior should be considered.  The difficulty in gas phase diffusion occurring equally at the reference as at the active electrode is in the large low-frequency impedance at low current.  Intuitively, the high impedance measured in the low frequency of the total cell is not likely to exist to an equal extent at the reference as at the active electrode due to gas concentration, because the change in gas concentration due to consumption is negligible near OCV.  As for electrolytic mass transport, high mass transfer impedance at low current is also observed in lanthanum strontium manganese oxide (LSM) electrodes, which are known to be poor bulk ion conductors at low overpotentials.  Impedance spectra obtained for thin, dense LSM films on yttria stabilized zirconia    [17, 25] have abscissa in a similar frequency regime as depicted in the present study for the total cell as shown in Figure 4.  These arcs have been attributed to oxide (or vacancy) diffusion through the film, as confirmed by the shift of the abscissa frequency at differing film thicknesses.  Even though this Ceramatec cathode is not composed of LSM, the properties of other manganite-based materials appear to be similar in the equivalent circuit model.Electrolytic mass transfer is further supported by comparing the fittedB value between button and larger cells.  For the specific cells in this study, the thicknesses of the mixed conducting layer of the electrodes were significantly different between the button cells and the large cells.  The total thickness (including mixed conducting components) for the large cell anode and cathode is approximately 52–55 μm and 77–80 μm, respectively, compared to 307–358 μm and 347–398 μm for the button cell.  The effect that thinner electrodes impose on mass transfer is visually evident in the low-frequency semicircle being fully resolved for the large cell impedance compared to that of the button cells over the same frequency spectrum (Figures 3 and 4).  The shift of the abscissa frequency between the two from approximately 10 mHz to 200 mHz is quantified by the B parameter of the FDE, which changes due to the diffusion layer thickness.  The ionic/interstitial oxygen (or vacancy) diffusion layer would be the sum of the thicknesses of the electrolyte and electrodes, assuming the majority of reactions occur on the outer edges of the electrodes.  The ratio of the button cell B to that of the large cell is approximately three at OCV, which is quite close to the ratio of the cathode thicknesses (2.76–3.04).  Comparing higher current levels is likely not valid, because of temperature changes due to reaction heat generation in the larger cells.


### 3.3.2. Gas Conversion Impedance

Primdahl and Mogensen [26] observed a low-frequency impedance arc related to gas conversion, when the reference electrode was located in a different atmosphere from the working electrode.  They proposed that the passage of AC current changes gas concentration at the active surface, most noticeably at low frequencies, resulting in a Nernstian potential loss.  Such is not observed when the working and reference electrodes are at the same concentration, which could be the cause of the previously discussed reference i-V behavior (Section 3.3.1) in Figure 8.  Their study modeled the gas conversion impedance (GCI) with a parallel Rlf and Clf in series with the electrode equivalent circuit (similar to the equivalent circuit described in Section 3.1.2, Figure 5(a)) and approximately replicated the fitted values with a constantly stirred tank reactor (CSTR) model yielding the following relationships:(4)Rlf=RigT(2F)2J(1xw+1xH),(5)Clf=4F2PV(RigT)2A(1xw+1xH)-1
wherein the resistance and capacitance are area specific.  Applying Primdahl’s model to the cells in question indicates that gas conversion may be partially responsible for the extra arc.While a 61 cm2 planar cell is expected to deviate from a single tank CSTR model due to the previously discussed distributed conditions, the general variation of impedance with each of the independent variables should be unchanged.  Increases in cell current or hydrogen concentration predict that Rlf should decrease and Clf should increase, a tendency observed in the button cell fits in Tables 1(a) and 1(b).  The same impedance effect should be observed for decreased temperature; however, Clf fits decrease by no less than 44% for a temperature decrement of 50°C.  Temperature was not explicitly controlled for the larger cells to test whether the same anomaly is present, but hydrogen concentration (Table 3(a)) and cell current (Table 3(b)) were varied with only the former behaving as in the button cell tests.  Specifically, Clf is only slightly affected by increasing stack current from 1.639 mA cm−2 to 52.46 mA cm−2, although the model predictions and button cell data indicate a large increase should be observed.  This effect cannot be justified by the expected increase in cell temperature induced by increasing current, because the GCI model indicates an increment of ≫100°C would be necessary for such a drastic shift to be observed.  Notwithstanding the fact that the button cell fitted Clf exhibit opposite behavior (Table 1(a)) with respect to the model predictions (5) at increased temperature.  The opposite trend is not clear and further intensive work is necessary to examine the effect of temperature on the gas conversion impedance, especially on the value of Clf capacitance.Table 3
Values for a parallel RC fitted to the low-frequency impedance loops of three cells in a five-cell planar stack operating at 800°C.
(a)
The SOFCs were fed with 38.5 ~ 97% H2, 3% H2O, and balance N2, and operated at 1.69 mA cm−2
Physical elementsxH (Cell 1)xH (Cell 3)xH (Cell 5)38.5%64.6%97.0%38.5%64.6%97.0%38.5%64.6%97.0%Rlf(Ω cm2)1.4301.5381.3731.3401.4791.2631.2991.4271.301Clf(F cm−2)1.1851.5861.6580.94771.3421.4650.98831.2381.447Y0 (S cm−2 s1/2)1.1351.3321.4411.0181.2111.3971.0591.2021.338B(s1/2)2.0252.4872.4021.7272.2102.1551.7432.0982.149(b)
Varying load currents when fed with 97.0% H2 and balance H2O
Physical elementsi (mA cm−2, Cell 1)i (mA cm−2, Cell 3)i (mA cm−2, Cell 5)1.63952.46147.51.63952.46147.51.63952.46147.5Rlf (Ω cm2)1.3730.65700.58021.2630.62830.48621.3010.60820.4373Clf (F cm−2)1.6591.7161.6481.4651.4171.3801.4471.4451.895Y0 (S cm−2 s1/2)1.4422.0511.9251.3971.8571.8211.3371.9071.731B(s1/2)2.4021.6681.5542.1561.4671.3762.1481.4571.259(c)
Varying the anode feed ratej when the SOFCs were operated at 52.46 mA cm−2 and composed of 97.0% H2 and balance H2O, where the total j is assumed to be evenly distributed between the five cells in the stack.
Physical elementsj (μmol s−1 cm−2, Cell 1)j (μmol s−1 cm−2, Cell 3)j (μmol s−1 cm−2, Cell 5)3.1684.9793.1684.9793.1684.979Rlf (Ω cm2)0.86320.65700.81310.62830.77590.6082Clf (F cm−2)1.9101.7161.6481.4171.7341.446Y0 (S cm−2 s1/2)1.9562.0511.9001.8571.9521.907B(s1/2)2.0491.6681.8631.4671.8431.457The best evidence that gas conversion is likely responsible for the low-frequency impedance loop is listed in Table3(c), wherein an increase in the total anode gas flow decreases Rlf and decreases Clf  The relationship between gas conversion impedance and anode feed rate j has been demonstrated experimentally by Momma et al. [27], who applied it to determine the flow distribution through a 46-cell stack.  Whether the low-frequency arc is due entirely to gas conversion is not conclusive enough to apply the same method to test the flow distribution in the Ceramatec stack.Neither gas conversion nor mass transfer is likely to exclusively cause the low-frequency impedance loop.  The two can scarcely be distinguished from the appearance of the arc, because a parallel RC and FDE are similarly shaped on a Nyquist plot [28] as evident in the χ2 values of the fits.  The difference is in the Warburg-like part of the FDE, which appears as a 45° line on a Nyquist plot to the higher frequency side of the FDE abscissa [20].  Anode adsorption occurs in the same frequency regime as the Warburg part of the FDE, which obscures the shape of the low-frequency loop.  Exacerbating the choice of elements is the functionality of  (2) and (4), which indicates that Rlf and Y0 often vary inversely; however, such a relationship is not always observed in the values listed in Table 3.  The correlation between fitted values and changes in hydrogen concentration or temperature is more congruent with mass transfer, while gas conversion better explains that of gas flow rate.  The effects of gas conversion and mass transfer are most likely present concurrently being superimposed as an arc in the low-frequency region of a Nyquist plot.The utility of applying equivalent circuit models for mechanistically discriminating problems with individual cell performance is evident in analyzing the different large cells.  In the stack under investigation, Cell 3 is in the middle of the stack flanked by Cells  1 and 5 on the outsides.  The two outside cells are not identical, with Cell 5 having a thermally insulating layer adjacent to the outer wall that is unique to it.  Because the temperature of Cell 1 is lower than the other cells in the stack, it lacks the thermal energy necessary to drive the reaction as fast, thus limiting performance in a stack of cells connected in series.  One need not know the geometry of the system to deduce that Cell 1 is at a lower temperature than the others; such can be implied from the model fits.  While the increasedRct for the cell in question may indicate a decreased temperature, the trend may also be induced by inconsistent electrode microstructure, catalyst loading, and/or reactant concentration.  Coupling the Rct behavior with the increase in RΩ and B for Cell 1 relative to the other cells proves the root of the problem, because decreased electrolyte conduction and diffusivity are also indicative of lower temperatures.  The analyses of the fitted large cell data show how impedance spectroscopy can be used to mechanistically discriminate the factors that limit stack performance.


## 3.3.1. Mass Transfer

The impedance of mass transfer by diffusion is quantified by the fitted values of the FDE, which respond as expected to variations in temperature and reactant concentration.  Heating the cell imparts more thermal energy to the system facilitating mass transport and decreasing impedance (or increasing diffusivity).  Recalling the dependence of (1) and (2) on diffusivity verifies the change in the fitted Y0 and B with respect to temperature at the same cell current and reactant concentration (Table 1(a)).  Note that electrolytic diffusivity is more sensitive to temperature variations (Di∝eT [23]) than gas phase diffusivity (Di∝T1.75 [24]) thereby dominating the T in the denominator of  (2) and in gas concentration (ci∝T-1).  The relationship between Y0 and reactant concentration is more straightforward, increasing impedance as the feed is diluted with inert gas as explicitly shown in (1) and verified for both cell sizes in Tables 1 and 2.  The approximately constant value of B with respect to changing hydrogen concentration is also as expected, as it is a variable that should not appreciably alter Di or δ.The variation of the FDE with respect to current is as expected for the button cells, but not so for the larger ones.  Data fits in Table1 show that button cell Y0 increases with increasing current while B remains relatively constant, which is predicted by the FDE equations.  The only variables affected by a change in current are the concentration terms in the summation of (2) explaining why B remains the same.  The value of the summation is largely dominated by the reactant at lower concentration (normally water for the anode); hence, Y0 is large at high current, where cw is increased from its initial amount due to the reduction of hydrogen to water.  The size of the cell likely contributes to the difference between the button and larger cells, because reactant depletion at the cell outlet is more significant at higher currents and expected to have a more deleterious effect on larger cells.To explain the paucity of a similarly large low-frequency arc in half-cell impedance spectra, the potential of the reference electrodes must be considered.  If the concentrations of all reactants at the references are constant, then the potential difference between the two should be independent of the current through the cell.  The curves plotted in Figure8 indicate that concentration in fact varies at the reference surface.  The absence of the low-frequency arc from anode or cathode impedance spectra is then justified as existing at working and reference electrode to the same extent, thus canceling each other out in the single electrode measurements.Figure 8
Potential difference between anode and cathode reference electrodes plotted as a function of load current.Since reference concentrations have been shown to vary with respect to current, whether some or all diffusing species are responsible for the behavior should be considered.  The difficulty in gas phase diffusion occurring equally at the reference as at the active electrode is in the large low-frequency impedance at low current.  Intuitively, the high impedance measured in the low frequency of the total cell is not likely to exist to an equal extent at the reference as at the active electrode due to gas concentration, because the change in gas concentration due to consumption is negligible near OCV.  As for electrolytic mass transport, high mass transfer impedance at low current is also observed in lanthanum strontium manganese oxide (LSM) electrodes, which are known to be poor bulk ion conductors at low overpotentials.  Impedance spectra obtained for thin, dense LSM films on yttria stabilized zirconia    [17, 25] have abscissa in a similar frequency regime as depicted in the present study for the total cell as shown in Figure 4.  These arcs have been attributed to oxide (or vacancy) diffusion through the film, as confirmed by the shift of the abscissa frequency at differing film thicknesses.  Even though this Ceramatec cathode is not composed of LSM, the properties of other manganite-based materials appear to be similar in the equivalent circuit model.Electrolytic mass transfer is further supported by comparing the fittedB value between button and larger cells.  For the specific cells in this study, the thicknesses of the mixed conducting layer of the electrodes were significantly different between the button cells and the large cells.  The total thickness (including mixed conducting components) for the large cell anode and cathode is approximately 52–55 μm and 77–80 μm, respectively, compared to 307–358 μm and 347–398 μm for the button cell.  The effect that thinner electrodes impose on mass transfer is visually evident in the low-frequency semicircle being fully resolved for the large cell impedance compared to that of the button cells over the same frequency spectrum (Figures 3 and 4).  The shift of the abscissa frequency between the two from approximately 10 mHz to 200 mHz is quantified by the B parameter of the FDE, which changes due to the diffusion layer thickness.  The ionic/interstitial oxygen (or vacancy) diffusion layer would be the sum of the thicknesses of the electrolyte and electrodes, assuming the majority of reactions occur on the outer edges of the electrodes.  The ratio of the button cell B to that of the large cell is approximately three at OCV, which is quite close to the ratio of the cathode thicknesses (2.76–3.04).  Comparing higher current levels is likely not valid, because of temperature changes due to reaction heat generation in the larger cells.


## 3.3.2. Gas Conversion Impedance

Primdahl and Mogensen [26] observed a low-frequency impedance arc related to gas conversion, when the reference electrode was located in a different atmosphere from the working electrode.  They proposed that the passage of AC current changes gas concentration at the active surface, most noticeably at low frequencies, resulting in a Nernstian potential loss.  Such is not observed when the working and reference electrodes are at the same concentration, which could be the cause of the previously discussed reference i-V behavior (Section 3.3.1) in Figure 8.  Their study modeled the gas conversion impedance (GCI) with a parallel Rlf and Clf in series with the electrode equivalent circuit (similar to the equivalent circuit described in Section 3.1.2, Figure 5(a)) and approximately replicated the fitted values with a constantly stirred tank reactor (CSTR) model yielding the following relationships:(4)Rlf=RigT(2F)2J(1xw+1xH),(5)Clf=4F2PV(RigT)2A(1xw+1xH)-1
wherein the resistance and capacitance are area specific.  Applying Primdahl’s model to the cells in question indicates that gas conversion may be partially responsible for the extra arc.While a 61 cm2 planar cell is expected to deviate from a single tank CSTR model due to the previously discussed distributed conditions, the general variation of impedance with each of the independent variables should be unchanged.  Increases in cell current or hydrogen concentration predict that Rlf should decrease and Clf should increase, a tendency observed in the button cell fits in Tables 1(a) and 1(b).  The same impedance effect should be observed for decreased temperature; however, Clf fits decrease by no less than 44% for a temperature decrement of 50°C.  Temperature was not explicitly controlled for the larger cells to test whether the same anomaly is present, but hydrogen concentration (Table 3(a)) and cell current (Table 3(b)) were varied with only the former behaving as in the button cell tests.  Specifically, Clf is only slightly affected by increasing stack current from 1.639 mA cm−2 to 52.46 mA cm−2, although the model predictions and button cell data indicate a large increase should be observed.  This effect cannot be justified by the expected increase in cell temperature induced by increasing current, because the GCI model indicates an increment of ≫100°C would be necessary for such a drastic shift to be observed.  Notwithstanding the fact that the button cell fitted Clf exhibit opposite behavior (Table 1(a)) with respect to the model predictions (5) at increased temperature.  The opposite trend is not clear and further intensive work is necessary to examine the effect of temperature on the gas conversion impedance, especially on the value of Clf capacitance.Table 3
Values for a parallel RC fitted to the low-frequency impedance loops of three cells in a five-cell planar stack operating at 800°C.
(a)
The SOFCs were fed with 38.5 ~ 97% H2, 3% H2O, and balance N2, and operated at 1.69 mA cm−2
Physical elementsxH (Cell 1)xH (Cell 3)xH (Cell 5)38.5%64.6%97.0%38.5%64.6%97.0%38.5%64.6%97.0%Rlf(Ω cm2)1.4301.5381.3731.3401.4791.2631.2991.4271.301Clf(F cm−2)1.1851.5861.6580.94771.3421.4650.98831.2381.447Y0 (S cm−2 s1/2)1.1351.3321.4411.0181.2111.3971.0591.2021.338B(s1/2)2.0252.4872.4021.7272.2102.1551.7432.0982.149(b)
Varying load currents when fed with 97.0% H2 and balance H2O
Physical elementsi (mA cm−2, Cell 1)i (mA cm−2, Cell 3)i (mA cm−2, Cell 5)1.63952.46147.51.63952.46147.51.63952.46147.5Rlf (Ω cm2)1.3730.65700.58021.2630.62830.48621.3010.60820.4373Clf (F cm−2)1.6591.7161.6481.4651.4171.3801.4471.4451.895Y0 (S cm−2 s1/2)1.4422.0511.9251.3971.8571.8211.3371.9071.731B(s1/2)2.4021.6681.5542.1561.4671.3762.1481.4571.259(c)
Varying the anode feed ratej when the SOFCs were operated at 52.46 mA cm−2 and composed of 97.0% H2 and balance H2O, where the total j is assumed to be evenly distributed between the five cells in the stack.
Physical elementsj (μmol s−1 cm−2, Cell 1)j (μmol s−1 cm−2, Cell 3)j (μmol s−1 cm−2, Cell 5)3.1684.9793.1684.9793.1684.979Rlf (Ω cm2)0.86320.65700.81310.62830.77590.6082Clf (F cm−2)1.9101.7161.6481.4171.7341.446Y0 (S cm−2 s1/2)1.9562.0511.9001.8571.9521.907B(s1/2)2.0491.6681.8631.4671.8431.457The best evidence that gas conversion is likely responsible for the low-frequency impedance loop is listed in Table3(c), wherein an increase in the total anode gas flow decreases Rlf and decreases Clf  The relationship between gas conversion impedance and anode feed rate j has been demonstrated experimentally by Momma et al. [27], who applied it to determine the flow distribution through a 46-cell stack.  Whether the low-frequency arc is due entirely to gas conversion is not conclusive enough to apply the same method to test the flow distribution in the Ceramatec stack.Neither gas conversion nor mass transfer is likely to exclusively cause the low-frequency impedance loop.  The two can scarcely be distinguished from the appearance of the arc, because a parallel RC and FDE are similarly shaped on a Nyquist plot [28] as evident in the χ2 values of the fits.  The difference is in the Warburg-like part of the FDE, which appears as a 45° line on a Nyquist plot to the higher frequency side of the FDE abscissa [20].  Anode adsorption occurs in the same frequency regime as the Warburg part of the FDE, which obscures the shape of the low-frequency loop.  Exacerbating the choice of elements is the functionality of  (2) and (4), which indicates that Rlf and Y0 often vary inversely; however, such a relationship is not always observed in the values listed in Table 3.  The correlation between fitted values and changes in hydrogen concentration or temperature is more congruent with mass transfer, while gas conversion better explains that of gas flow rate.  The effects of gas conversion and mass transfer are most likely present concurrently being superimposed as an arc in the low-frequency region of a Nyquist plot.The utility of applying equivalent circuit models for mechanistically discriminating problems with individual cell performance is evident in analyzing the different large cells.  In the stack under investigation, Cell 3 is in the middle of the stack flanked by Cells  1 and 5 on the outsides.  The two outside cells are not identical, with Cell 5 having a thermally insulating layer adjacent to the outer wall that is unique to it.  Because the temperature of Cell 1 is lower than the other cells in the stack, it lacks the thermal energy necessary to drive the reaction as fast, thus limiting performance in a stack of cells connected in series.  One need not know the geometry of the system to deduce that Cell 1 is at a lower temperature than the others; such can be implied from the model fits.  While the increasedRct for the cell in question may indicate a decreased temperature, the trend may also be induced by inconsistent electrode microstructure, catalyst loading, and/or reactant concentration.  Coupling the Rct behavior with the increase in RΩ and B for Cell 1 relative to the other cells proves the root of the problem, because decreased electrolyte conduction and diffusivity are also indicative of lower temperatures.  The analyses of the fitted large cell data show how impedance spectroscopy can be used to mechanistically discriminate the factors that limit stack performance.


## 3.4. Fuel Cell Stack

Data for the entire fuel cell stack were also obtained under the same test conditions as the single cell tests previously described in Sections3.2 and 3.3.  Figure 9 depicts some of the obtained curves, each of which is similar in shape to the individual cells; hence, the same equivalent circuit that was fitted to the large cell data has been similarly applied to the stack data.  The key to the success of applying an individual cell model to a stack of cells is the distribution of  RC circuits of each physical relaxation process from cell to cell.  Although not precisely equal between the five cells, the frequency ranges at which each physical process dominates cell impedance are similar enough for the arcs to superimpose.  Had the degree of distribution been greater, the semicircles of the Nyquist plot would have appeared to be depressed (CPE-like); instead, the individual cell model is well fitted to the stack impedance data, as exhibited in this figure.Figure 9
Nyquist plot of a planar fuel cell stack fed with 97% H2 and balance H2O operating at 800°C and varying load current.  The data have been fitted to a four-RC equivalent circuit.While the individual cell tests were designed to determine which processes limit cell output and to what degree, the stack test is primarily intended to obtain a simple, purely electronic model to replicate the dynamic response of the stack.  In order to represent the stack the purely physical electronic elements, the FDE of the equivalent circuit (Figure6) has been replaced with a parallel RC, which has already been shown to fit single cell data equally well.  Obtaining fits for the circuit elements over the expected range of operation allows for the stack to be dynamically modeled over said range.  The utility in eliminating all distributed elements from the model is to facilitate circuit simulation [12], which can be used to design process controls and to assess the optimum configuration of operating conditions, power conditioning systems, and combination of multiple stacks.


## 4. Conclusions

In the work presented, impedance spectroscopy was successfully applied to solid oxide fuel cells and stacks at different operating conditions.  Improved equivalent circuit models were proposed and applied for fitting of button cells and larger cells in the stack.  A relaxation process unique to the total cell data was measured at low frequency and has been suggested to be due to the combined effects of mass transfer and gas conversion.  The total-cell circuit model is further simplified and adjusted to simulate the SOFC larger cells in a five-cell stack.  The cause for an underperforming cell in the stack was attributed to temperature, based on a comparison with the data from two other cells in the same stack.  The entire fuel cell stack was fitted to a model developed from the single cells, only with distributed elements replaced with physical electrical elements.  The stack circuit diagram and associated physical elements are essential to simulate its dynamic behavior to facilitate system design.  The work contained herein is useful for assessing in situ stack dynamics, stack aging, state of health, system controls, component fabrication, and quality control.



---

*Source: 465452-2011-10-01.xml*