3-D Modeling of a Vanadium Redox Flow Battery with Focus on Feeding Rate Effects
Xiao-Guang Yang,Qiang Ye,Ping Cheng
DOI: https://doi.org/10.1149/ma2014-03/2/5
2014-01-01
ECS Meeting Abstracts
Abstract:Renewable energy sources, like wind and solar, are being deployed worldwide in large numbers than ever before in efforts to build a low-carbon society. These renewable sources, however, are random and intermittent in nature, posing great threat to the stability of electric grids. Energy storage is a key technology to solve this issue. [1] Among the various existing energy storage technologies, the allvanadium redox flow battery (VRFB) has received enormous attention from both academia and industry. A great deal of research has been conducted in recent years with efforts to enhance the VRFB performance, most of which focuses on material research, such as modification of the existing materials or searching for alternatives. Apart from the material properties, as a matter of fact, the performance of a VRFB is also strongly affected by the operating parameters, such as the feeding rate (FR) of both anolyte and catholyte. Ma et al. [2] reported in an experimental study that both the charging and discharging capacities of a VRFB stack dropped significantly with the decrease of the FR, which can be attributed to the high concentration overpotential at the low FR. Though increasing the FR can alleviate this overpotential, it is at the expense of increased consumption of pump power. Therefore, some tradeoffs between pump and performance losses are required to determine the optimal FR. An effective parameter to evaluate the magnitude of the FR is the reactant stoichiometry (RS), defined as the ratio between the amount of reactant actually supplied to the cell and the amount needed theoretically to sustain the output current. In most of previous studies, the FR was maintained constant in the entire charge/discharge cycle, even though the reactant concentration changes dramatically during the cycle. This can lead to a colossal waste of energy. For instance, if the FR corresponds to the RS of 2.0 at the state of charge (SOC) of 90% in the charging process, the corresponding RS at 10% SOC is as high as 18 under this FR. Therefore, optimized FR strategy is needed to enhance the system efficiency. For this purpose, it is vital to gain a deep understanding on the effect of RS over the cell performance. This kind of research, however, is rather scarce in literature. Due to the high cost and lengthy time requirements of experimental study, numerical modeling becomes an indispensable tool in investigating the VRFB performance and several VRFB models have been proposed in the past few years, which provide valuable tools in predicting the cell performance. Nevertheless, most of these studies utilize a constant FR during the entire cycle and thus do not investigate the effect of the RS. Moreover, most of these models are 2-D, which cannot reflect the variation of species and current distributions along the in-plane direction. In this study, we present a 3-D VRFB model and highlight the RS effect on the cell performance. Special attention is paid to the species and current distributions along the in-plane direction which cannot be reflected in previous 2-D models. Figure 1 plots the variation of cell voltage with the RS at the SOC of 50% in the charging process. The operating current density is 40 mA cm, and the electrolyte is composed of 1.5 M vanadium and 3 M sulfuric acid. The cell voltage predicted by a 2-D model at the same operating condition is also given in this figure for comparison. As shown in this figure, the cell voltage varies little at high RS, but increases dramatically when the RS becomes lower than 5.0, indicating dramatic rise of the concentration overpotential at low flow velocities. It can also be noted from this figure that the cell voltage predicted by the present 3-D model is higher than that by the 2-D model, especially at low RS values, indicating that the 3-D model can account for the additional performance loss due to the variation of species concentration along the in-plane direction. Figure 2 reveals the uniformity of the distributions of V concentration and current density in the negative electrode along the in-plane direction with different RS values at the same operating condition as that in Figure 1. It can be seen from this figure that both the species and current distributions become highly non-uniform at lower RS values.