Mathematical Modeling Analysis and Optimization of Key Design Parameters of Proton-conductive Solid Oxide Fuel Cells

A proton-conductive solid oxide fuel cell (H-SOFC) has the advantage of operating at higher temperatures than a PEM fuel cell, but at lower temperatures than a SOFC. This study proposes a mathematical model for an H-SOFC in order to simulate the performance and optimize the flow channel designs. The model analyzes the average mass transfer and species' concentrations in flow channels, which allows the determination of an average concentration polarization in anode and cathode gas channels, the proton conductivity of electrolyte membranes, as well as the activation polarization. An electrical circuit for the current and proton conduction is applied to analyze the ohmic losses from an anode current collector to a cathode current collector. The model uses relatively less amount of computational time to find the V-I curve of the fuel cell, and thus it can be applied to compute a large amount of cases with different flow channel dimensions and operating parameters for optimization. The modeling simulation results agreed satisfactorily with the experimental results from literature. Simulation results showed that a relatively small total width of flow channel and rib, together with a small ratio of the rib's width versus the total width, are preferable for obtaining high power densities and thus high efficiency. Nomenclature: Mass transfer area of cross section (m 2) Species molar concentration (mol/m 3) OPEN ACCESS Effective diffusivity (m 2 /s) , Knudsen diffusivity (i—species) (m 2 /s) Ordinary diffusivity (i and j—species) (m 2 /s) Pore size (µm) Grain size (µm) E Electromotive force from the Nernst Equation (V) E act,a Activation energy level at the anode side (J/mol) E act,c Activation energy level at the cathode side (J/mol) Faraday's constant 96,485.3 (C/mol) ℎ Convective mass transfer coefficient m/s H-SOFC Proton conducting Solid Oxide Fuel Cell Current Density A/m 2 Species molar flux (mol/s) Exchange current density (A/m 3) _ Current density Molar mass kg/mol Molar consumption rate mol/s Partial pressure R Gas constant 8.314 (J/mol/K) Average pore radius of the electrodes (µm) T Temperature (K) Velocity (m/s) Mole fraction Y Ratio of grain contact neck to the grain size Greek Symbols Tortuosity Porosity Electrode thickness (m) α Transfer coefficient Subscripts and Superscripts int Interface (electrode-functional layer) conc Concentration polarization c Cathode a Anode af Anode functional layer cf Cathode functional layer H 2 Hydrogen O 2


( )
Effective diffusivity (m 2 /s) , Knudsen diffusivity (i-species) (m 2 /s) Ordinary diffusivity (i and j-species) (m 2 /s) Pore size (µm) Grain size (µm) E Electromotive force from the Nernst Equation (V) E act,a Activation energy level at the anode side (J/mol) E act,c Activation energy level at the cathode side (J/mol) Faraday's constant 96,485.Oxygen N 2 Nitrogen S Surface ∞ Bulk flow

Introduction
Having an electrolyte conductive to oxide ions, a solid oxide fuel cell (O-SOFC) works at relatively high operating temperatures, which helps to maintain a low activation polarization and eliminate the use of expensive catalysts in a fuel cell [1].However, high operating temperatures also result in disadvantages including potential thermal fatigue/failure of the cell material and gas sealing, as well as the thermal stress in the ceramic cell components [2].With the solid oxide electrolyte being conductive to protons, a proton-conductive solid oxide fuel cell (H-SOFC) works at a relatively lower temperature [3] than a regular SOFC.A lower operation temperature helps to alleviate the problems of thermal stress and thermal expansion mismatch related to high operating temperatures in a regular SOFC.An H-SOFC also allows more utilization of the fuel (H 2 ) and thus a better efficiency than a regular SOFC [4].
A number of research efforts have been devoted to experimental studies on proton conductive solid oxide fuel cells.Some of them focused on property of materials used for H-SOFC.For the key component material of an H-SOFC, Zhao et al. studied the performance of H-SOFC using BaCe 0.7 In 0.3−x Y x O 3−δ as the electrolyte material [5].Ling et al. studied the fuel cell performance using a stable La 2 Ce 2 O 7 as the electrolyte material [6].Guo et al. studied the performance of a carbon dioxide-tolerant proton-conducting solid oxide fuel cell with a dual-layer electrolyte [7].H-SOFC was also tested using BaZr 0.8 In 0.2 O 3−δ as the proton-conductive electrolyte [8].For electrode materials, Lin et al. evaluated the performance of H-SOFC using BaCo 0.7 Fe 0.2 Nb 0.1 O 3−δ (BCFN) as the cathode material [9].Deng et al. investigated fuel cell performance using a cathode made of the material of PrBa 0.5 Sr 0.5 Co 2 O 5−δ [10].Zhao et al. studied the performance of a cobalt-free proton-conductive oxide fuel cell performance using Ba 0.5 Sr 0.5 Fe 0.8 Cu 0.2 O 3−δ as cathode material [11].For fabrication related issues, Tsai et al. studied the tortuosity in electrodes materials in an anode-supported H-SOFC [12].
While the materials of the H-SOFC component are very important, a good design and management of the flow field is also very helpful to a better performance of the fuel cells.The current work presents studies and design optimization of a H-SOFC via modeling and simulation analysis.The modeling will give an easy-to-approach and comprehensive analysis to the mass transfer, activation and ohmic polarizations, and can predict the fuel cell performance and thus optimize the flow channel designs [13].The simulation results are to be compared with the experimental results obtained from literature to validate the model.The power density and voltage output due to the optimization of the dimensions of flow channels and ribs will be presented.

Basic Aspects of an H-SOFC
A proton conducting solid oxide fuel cell consists of a metal oxide electrolyte sandwiched between two electrodes.Fuel is supplied on the anode side which is oxidized into protons in an electrochemical reaction as: (1) This produces electrons which move towards cathode via an external circuit driven by the potential difference between the two electrodes.The fuel consumption leads to concentration gradient of hydrogen (if fuel is a mixture) in anode and therefore feed to the anode-electrolyte interface by permeation [14].In cathode, air is supplied, which offers oxygen to react with the receiving electrons and protons from anode to form water from the following electrochemical reaction: (2) The anode material and catalyst are required to be highly active to improve the H-SOFC performance.These days, mixtures of nickel oxide and electrolyte are used as anode support materials [15].Nickel behaves like a catalyst which increases the rate of chemical reaction and oxidizes the fuel at the interface.The anode reaction occurs at the interface of anode and electrolyte.Together with the reactant phase, the anode and electrolyte forms the so-called triple phase boundary.Mixtures of NiO and electrolyte increase the triple phase boundary.Electrolyte material behaves like a proton conductor which conducts protons and transfers it across electrolyte layer.Good proton conductivity is important to the increase of the rate of the overall reaction [16].The electrolyte must be as thin as possible and have low activation energy.The perovskite structured compound (ABO 3 specially AZrO 3 ) family are considered especially good electrolyte materials.An ideal cathode material must be chemically non-reactive with the electrolyte.A cathode material must be porous with high oxygen permeability.Due to these reasons novel perovskite oxide materials are chosen as cathode.Recently, Ba 0.5 Sr 0.5 Zn 0.2 Fe 0.8 O 3−δ (BSZF) has been developed as a novel cobalt free oxygen permeable membrane with high permeation behavior and good chemical stability at high temperatures [17].
Different from an oxide-ion conductive SOFC, H-SOFCs have the electromotive force expressed as: ( Due to its high operation temperature, the water in H-SOFC is in vapor state and is assumed as an ideal gas, the same as other gas species.

Numerical Modeling to H-SOFC
Figure 1 shows the schematic of the typical elements/components of an H-SOFC, which includes the electrolyte, electrode components, flow channels, as well as the flow channel walls also acting also as local current collectors.BaCe 0.9 Y 0.1 O 2.95 (BCY) and BaCe 0.5 Zr 0.3 Y 0.16 Zn 0.04 O 3−δ (BCZYZn) are used as electrolytes for the current mathematic model [15,17].The anode material is NiO which has a substrate of NiO-BCY and NiO-BCZYZn.The cathode materials are Ba 0.5 Sr 0.5 Zn 0.2 Fe 0.8 O 3−δ (BSZF) and Ba 0.5 Sr 0.5 Co 0.8 Fe 0.2 O 3−δ (BSCF).The electrode and electrolyte materials properties used for the simulation are presented in Tables 1 and 2. The contact resistances between current collector and electrodes vary with operation temperatures, which are given in Table 3. Fuel was supplied on the anode side, and cathode was fed with air for the needed oxygen.The utilization of fuel and oxygen are 85% and 50%, respectively, and a unit length of 1.0 m of the channel length (normal to the paper in Figure 1) is      To obtain the electromotive force, the partial pressures of species during an electrochemical reaction are needed.The partial pressure of a species is proportional to its molar fraction in a mixture.Therefore, the mass transfer processes of species from bulk flow to the reaction site (the electrode and electrolyte interface) have to be analyzed.

Anode Side Mass Transfer Analysis
The convective mass transfer flow rate of hydrogen from bulk flow to the surface of anode is expressed in Equation ( 4).The consumed hydrogen flow rate ṅ equals to the current divided by 2F, where F is the Faraday constant.(4) where h H2 is convective mass transfer coefficient, is the average molar concentration of hydrogen in the bulk flow, is the average molar concentration of hydrogen at the anode surface, and A is the mass transfer area of the electrode in exposure to the flow channels.
The mass diffusion for hydrogen and vapor in porous electrodes are analyzed by calculating the mass transfer flow rate as given in Equations ( 5) and (6).For all the analyses, the mass fluxes from surfaces toward the bulk flow are designated to be positive.
(5) (6) where υ is an anode diffusion velocity due to the mass diffusion of hydrogen, C H2 and C H2O are the average molar concentration of hydrogen and water in the porous layer, respectively.
In the fuel channel, the hydrogen is consumed and protons conduct through the electrolyte to form water with oxygen on the cathode side.Therefore, there is no water vapor flux on anode, which gives J H2O = 0.The effective diffusivity of hydrogen-water system comprises of Knudsen and ordinary diffusion.
(7) (8) where τ and ε represent the tortuosity and porosity of porous electrode, respectively.Knudsen diffusivities D H2,k and D H2O,k are calculated using , where r e is the average pore radius of porous electrode.D H2−H2O is the binary diffusivity, which is obtained through the following equation [18], where V H2 , V N2 , V O2 and V H2O are 7.07, 17.9, 16.6 and 12.7, respectively [18].The total wall flow velocity due to mass transfer is directly related to summation of mass flux of all species and the total density at wall.Therefore, there is.
This is a basic equation describing mass transfer boundary conditions [19].Substituting diffusion velocity υ a back to Equations ( 5) and ( 6) and introducing the new effective diffusivities for hydrogen and vapor as given by Equations ( 11) and ( 12), we could obtain compact expression for mass flux as given by Equations ( 13) and ( 14).(11) (12) (13) (14) Given the molar consumption fluxes of hydrogen and water to left-hand sides of the above equations and considering the hydrogen molar concentration difference from the bulk flow to the anode surface, the concentration difference of hydrogen through the porous layer is given in Equation (15).Correspondingly, the concentration difference of water is given in Equation ( 16). ( 15) (16) To consider the mass diffusion resistance in anode functional layer, we can simply add the mass transfer resistance term of the functional layer, , to the parenthesis on the left-hand side of Equations ( 15) and ( 16).Here δ af is the thickness of the anode functional layer.The equation for is in the same form as that of Equation (11).However, the parameters used for calculating are from anode functional layer.Equations for mass fluxes for water vapor in the functional layer are similar to those in the anode layer, which are not presented here.
From the above analysis, it is clear that with the given current density and the concentration of species in the bulk flow, one can obtain the species concentration at the anode/electrolyte interface.

Cathode Side Mass Transfer Analysis
On the cathode side, there are three species, O 2 , N 2 and H 2 O involved in the mass transfer.Similar to the anode side, the same analysis should be applied to the cathode side to find out the concentration of O 2 , N 2 and H 2 O at electrolyte/cathode interface.The molar consumption rate of O 2 is known as the total current divided by 4F, where F is the Faraday's constant.Analogous to mass transfer, the mass transfer flux of oxygen from bulk flow to the cathode surface is similar to that of hydrogen at the anode side.
The mass transfer of oxygen, nitrogen and water in porous cathode layer are described using the following equations.(17) (18) (19) where υ c is an overall cathode mass diffusion velocity due to all species' mass diffusion on cathode.The mass transfer fluxes of oxygen and water are related to the current densities of the fuel cell.Nitrogen is not involved in any reaction, which has a flux of zero.The total mass diffusion velocity υ c is in the form of: (20) The effective diffusivities of species in mixture are given in Equations ( 21)-( 23).The ordinary diffusivity of one species in a mixture of more than two species is given by Equation (24).(24) Here, cathode side effective diffusivities for three of the species are given in Equations ( 25)-( 27).These diffusivities are used for cathode species' mass transfer, given by Equations ( 28)-(30) after substituting cathode velocity to the mass flux equations. ( ) ) ) To consider the oxygen mass diffusion resistance in cathode functional layer, a new term, R cfl , given in Equation ( 31) is introduced as follows, which can be added to the left-hand side of Equation (28) for cathode functional layer. (31 where δ cf is the thickness of the cathode functional layer.The equation for D cf−O2(eff) is of the same form of Equation ( 25), however, parameters for D cf−O2(eff) must be from the cathode functional layer.Similar process could be applied to solve for nitrogen and water on cathode side.The mass concentration for these species at the interface between electrolyte and functional layer could be solved using similar approach with consideration of a resistivity term of R cfl_N2 and R cfl_H2O .The resistivity terms for these two species are listed as follows, by Equations ( 32) and (33).
(32) (33) From the above analysis, it is clear that at any given current density, the molar concentrations of all species at the electrode/electrolyte can be obtained.Finally, the partial pressures of all species are related to the molar concentration in the following forms:

Activation Polarization
The activation polarization is given by the Butler-Volmer equation [20], as given in Equation (37), where J _current is the current density, J 0 is the exchange current density, is the transfer coefficient which is typically set as 0.5, is the number of electrons passed through the external circuit for every mole of fuel oxidation.For H-SOFC, the value of is 2. Therefore, the activation polarization could be expressed in terms of current density J, as shown in Equation (38).
The exchange current density is related to material property and operation conditions (porosity, pore size, temperature and pressure).According to literature [20], exchange current density for O-SOFC is expressed in terms of the effect of both micro structural properties and operating condition as given in Equations ( 39) and (40).In Equations ( 39) and ( 40), E act,a (=1.0 × 10 5 J/mol) and E act,c (=1.2 × 10 5 J/mol) are activation energy levels at anode and cathode, respectively; Y is the ratio of grain contact neck to the grain size and n is the porosity of electrodes; D p and D s represent the pore size and grain size; k a and k c are adjustable coefficients, which are reported as k a of 6.634 × 10 −8 and k c of 7.534 × 10 −8 in reference [20].Here we assume that Equations ( 39) and ( 40) are applicable for H-SOFC as well as O-SFOC.

Ohmic Loss
The flow of electrons and ions in fuel cell components results in ohmic polarization due to ohmic resistance in all the layers as well as electrical contact resistance [21] between layers.In order to have a precise ohmic loss analysis, an equivalent electrical circuit for a flow channel and its two walls (local current collectors) was constructed as shown in Figure 2.There are five layers in the electrode assembly.Electrode layers for both anode and cathode are porous material for reactants to diffuse through.Another layer is functional layer where electrochemical reaction takes place.Having functional layer is the recent technology that can enhance the electrochemical reaction by creating more morphological contact between electrode material and electrolyte material.The temperature of functional layer is slight higher than other layers due to exothermic reaction.The current conduction route from a typical anode-side current collector to a cathode-side current collector can be discretized into multiple segments.With the symmetric discretization of the flow channel, computational time for the circuit can be significantly reduced.The electrical potentials at all the nodes are calculated using Kirchhoff's current law, which states that the summation of current flow into the node should be zero.A similar method has been applied to calculate the ohmic loss in PEM fuel cells and O-SOFC [22][23][24].
The electromotive force shown in the equivalent circuit should be obtained from Equation (3) which also subtracts the activation polarizations of both anode and cathode.The species partial pressures obtained from the mass transfer analysis are used in Equation (3)._ 0 (1 ) exp( ) exp( )

Computational Procedures
The numerical calculation follows a particular procedure, which is outlined as a flow chart in Figure 3. First, all the physical properties and dimensions of the fuel cell are defined.With the known amount of reactants as well as the prescribed current density, the values of concentrations of species at intermediate layers were calculated through iterations of the equations for the mass transfer, which consequently converges.Using these values of concentration, the partial pressures of all species can be calculated which thus considers the concentration polarization in the electromotive force as given in Equation ( 3).The electrical circuit was then analyzed for ohmic losses by discretizing the single channel fuel cell structure with multiple nodes.The electromotive force shown in the electrical circuit is the value from Equation (3) subtracted with the activation polarization.The activation polarization was considered through Volmer-Butler's equation.From the solution of electrical circuit, the fuel cell voltage is obtained with the given current density and other conditions.Multiple calculations are conducted for a specific temperature and a range of current densities for the V-i curve.The power output from the fuel cell is easily obtained as the product of cell voltage and current, and therefore, the power density versus current density is also obtained for the fuel cell.

Validation to the Modeling
The analytical modeling is validated by comparing the simulation results with experimental data in literature.The cited experimental results are from references [15].The studies in the reference reported operation temperatures of the H-SOFC up to 700 °C.
Material properties of the H-SOFC from reference [15] are listed in Tables 4 and 5, which were used in the modeling simulation.The tortuosities used for simulation are properly assumed for H-SOFC electrodes.The conductivity of electrode changes with temperature insignificantly compared to that of electrolyte, therefore, it is treated as a constant, as given in Table 5.Current simulation model has a flow channels with width of 1.2 mm and rib width 0.6 mm for both anode and cathode channels.The present simulation results of V-I curves are compared with the experimental results from reference [15], as given in Figure 4.The simulation results agree with the experimental data very well, which indicates the validity of the model as well as the related physical properties of the fuel cell components from the literature.The contact resistances at interface at different operation temperatures are shown in Table 3.The comparison of the modeling results with experimental data from references has evident that the currently developed H-SOFC model is effective and reliable.In the following section the model will be used to optimize the geometrical designs of the flow channels in H-SOFC.

Optimization of the Dimensions of Flow Channels and Walls of an H-SOFC
A quite number of researchers have reported studies to flow channel and flow field designs, which have been recognized to affect the fuel cell performance in both PEM fuel cells [25] and SOFCs.For PEM fuel cells, it is widely known that serpentine flow channels are known as the "industry standard" due to the better performance compared to that of many other different designs [26].Serpentine flow channels are commonly adopted in PEM fuel cells [22].In the present work, serpentine flow channels are considered in the modeling and simulation.The dimensions of flow channels that affect the H-SOFC performance will be optimized.The material property of electrode and electrolyte from Tables 1-3 are used for the simulation.
Two important parameters, the total width L o , which includes the width of the flow channel and the width of one of the channel walls, and the ratio of the rib's width versus the total width, are to be optimized for the serpentine flow channels.The rib width is designated as Lc and the channel width is L o −L c , as seen in Figure 2. The rib ratio R Lc is then given as: (41) Before the optimization to the flow channel designs, the effect of the operating temperature is studied, and an operating temperature is selected for the optimization of the dimensions of the channels.For this study, the total width is 3 mm and the rib ratio R Lc is 0.3.The operation temperature ranges from 500 °C to 700 °C.As seen in Figure 5 the output power density increases with the increase of temperature, which is the same phenomenon as was shown in the model validation simulations and experimental data from literature.High operation temperature makes it possible to have greater reaction kinetics and greater proton conductivity, and also a reduced contact resistance between electrode and current collector, which thus results in higher output voltage.The flow channel geometry optimization simulation is conducted under operation temperature of 700 °C which allows it to achieve high power density among the simulated temperatures.Three groups of dimensions of channels are selected for optimization design, which have total width L o in 2 mm, 4 mm, and 6 mm, respectively.In each group of dimensions, five rib ratios are selected.Figure 6 presents the power density versus current density of the three groups of cells.
Figure 6a-c show a common characteristics that the cell power output increases with the decrease of the rib ratio for the studied cases under a constant total channel width L o .The maximum output power density in Figure 6a is around 0.64 W/cm 2 at a rib ratio of R Lc = 0.05, which is 70% greater than the case with a rib ratio of R Lc = 0.5.Obviously, the smaller rib ratio means more opening area of electrochemical reactions, and therefore, the fuel cell power density increases.On the other hand, with more and more decrease of the rib ratio the increase of the fuel cell power density becomes less significant, which implies that too small rib ratio should not be selected.A small width of the rib also means high contact resistance.It can also be difficult to fabricate the flow channels with very small walls.
At larger total widths, L o = 4 mm, and 6 mm, the decrease of the rib ratio also results in the increase of the fuel cell power density.However, the maximum power densities from the three cases of different total widths are different.With smaller total width, the maximum power density is higher than that with larger total width.The maximum power density in Figure 6a-c is 0.64 W/cm 2 , 0.6 W/cm 2 and 0.58 W/cm 2 , respectively.
To have an overview of the optimization of flow channels, Figure 7 shows the effect of rib ratio to the maximum power density at different operation temperatures, ranging from 500 °C to 700 °C.The selected another three cases of channel total widths are 1.5 mm, 3.0 mm, and 4.5 mm.The rib ratio in each case varies from 0.01 to 0.3.
It is seen from Figure 7 that at all the studied temperatures, the decrease of the rib ratio is preferable, which results in the increase of the maximum power density.At different operation temperatures, the optimum rib ratio can be slightly different.For example, at temperature of 500 °C, the optimum rib ratio is around 0.1, while at temperature of 700 °C the optimum rib ratio is smaller.Nevertheless, nearly one optimum rib ratio can be largely applicable to all the different operation temperatures, for example, a rib ratio of 0.05, in all the studied cases.To show the overview of the effect of both the total width L o and the rib ratio R Lc , Figure 8 presents a contour of the maximum power density as a function of the total width and rib ratio.It is clear that the peak of the maximum power density locates at the coordinate of both small total width and small rib ratio.For a fixed total width (summation of channel width and rib width), the increase of R Lc means that channel rib area increases which will tend to improve the current collection efficiency.However, increase of R Lc results in the decrease of channel flow area, which does not benefit to improvement of the cell power density.This is because that decrease of channel area will decrease the exposure of reactants to electrode layer and reduce the diffusion and reaction, which will result in the drop of power density.It is also found that the maximum power density from fuel cells is more sensitive to the variation of the rib ratio rather than the total channel width.The channel and rib width affect the performance of proton conducting SOFCs and it also has effect on conventional SOFCs and other type of fuel cells [13].

Conclusions
An analytical model was developed for an H-SOFC in order to study the fuel cell performance and also to optimize the flow channel designs.This is a one-dimensional model that was validated by experimental data from the literature and was used to simulate and optimize the flow channel dimensions, mainly two important parameters-the total width of flow channel with channel rib/wall and the ratio of the channel rib/wall to the total width.
It was found from the optimization study that the fuel cell maximum power density can be higher if one chooses a smaller total width as well as a smaller ratio of the width of channel rib/wall against the total width.However, this ratio has an optimum point, too small a value can also result in decrease of the maximum power density.To show this feature, a contour of the maximum power density as a function of the total width and the rib ratio is provided in the paper.
The optimum ratio of channel rib/wall versus the total width varies slightly with the difference of the operation temperature.Nevertheless, nearly one optimum rib ratio can be largely applicable to all the different operation temperatures.

2
modeling.Setting the utilization percentage of fuel and oxygen being constant, the flow rates of fuel and air vary with the current density.

Figure 1 .
Figure 1.Schematic of H-SOFC showing only one pair of flow channels on two sides of the MEA.

Figure 2 .
Figure 2. Equivalent electrical circuit for an H-SOFC.

Figure 3 .
Figure 3. Computational procedures for simulation of proton-conductive solid oxide fuel cells.

Figure 4 .
Figure 4. Comparison of present modeling results with experimental data in the reference [15].(a) Operation temperature varies from 500 °C to 600 °C; (b) Operation temperature varies from 650 °C to 700 °C.

Figure 5 .
Figure 5. Fuel cell output power density versus current density at various operation temperatures.(Total channel width L o = 3 mm; rib ratio R Lc = 0.3).

Figure 6 .
Figure 6.Cell voltage and power density versus current density at cell temperature of 700 °C.(a) L o = 2 mm; (b) L o = 4 mm; (c) L o = 6 mm.

Figure 7 .Figure 8 .
Figure 7.The maximum power densities obtained at channel designs with different rib ratios and different operation temperatures.(a) L o = 1.5 mm; (b) L o = 3.0 mm; (c) L o = 4.5 mm.

Table 1 .
Physical properties of the electrodes materials.

Table 2 .
Proton conductivity of the electrolyte at different operating temperatures.

Table 4 .
[15]rial properties from reference[15]used for the present validation of the model.

Table 5 .
Conductivities of electrode and electrolyte at different temperatures.