Next Article in Journal
On the Performance of Portable NiMH Batteries of General Use
Previous Article in Journal
S, Se-Codoped Dual Carbon Coating and Se Substitution in Co-Alkoxide-Derived CoS2 Through SeS2 Triggered Selenization for High-Performance Sodium-Ion Batteries
Previous Article in Special Issue
Recent Advancements in Artificial Intelligence in Battery Recycling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhancing Mass Transport in Organic Redox Flow Batteries Through Electrode Obstacle Design

by
Joseba Martínez-López
1,
Unai Fernández-Gamiz
1,*,
Eduardo Sánchez-Díez
2,
Aitor Beloki-Arrondo
3 and
Íñigo Ortega-Fernández
3
1
Nuclear Engineering and Fluid Mechanics Department, University of the Basque Country UPV/EHU, Nieves Cano 12, 01006 Vitoria-Gasteiz, Spain
2
Centre for Cooperative Research on Alternative Energies (CIC EnergiGUNE), Basque Research and Technology Alliance (BRTA), Alava Technology Park, Albert Einstein 48, 01510 Vitoria-Gasteiz, Spain
3
TECNALIA Research & Innovation, Basque Research and Technology Alliance (BRTA), Mikeletegi Pasealekua 2, 20009 Donostia-San Sebastian, Spain
*
Author to whom correspondence should be addressed.
Batteries 2025, 11(1), 29; https://doi.org/10.3390/batteries11010029
Submission received: 14 December 2024 / Revised: 9 January 2025 / Accepted: 11 January 2025 / Published: 16 January 2025
(This article belongs to the Special Issue Towards a Smarter Battery Management System: 2nd Edition)

Abstract

:
This study examines the impact of incorporating obstacles in the electrode structure of an organic redox flow battery with a flow-through configuration. Two configurations were compared: A control case without obstacles (Case 1) and a modified design with obstacles to enhance mass transport and uniformity (Case 2). While Case 1 exhibited marginally higher discharge voltages (average difference of 0.18%) due to reduced hydraulic resistance and lower Ohmic losses, Case 2 demonstrated significant improvements in concentration uniformity, particularly at low state-of-charge (SOC) levels. The obstacle design mitigated local depletion of active species, thereby enhancing limiting current density and improving minimum concentration values across the studied SOC range. However, the introduction of obstacles increased flow resistance and pressure drops, indicating a trade-off between electrochemical performance and pumping energy requirements. Notably, Case 2 performed better at lower flow rates, showcasing its potential to optimize efficiency under varying operating conditions. At higher flow rates, the advantages of Case 2 diminished but remained evident, with better concentration uniformity, higher minimum concentration values, and a 1% average increase in limiting current density. Future research should focus on optimizing obstacle geometry and positioning to further enhance performance.

1. Introduction

The energy scene is transforming into one moving toward a sustainable future. Renewable resources such as solar and wind are considered cleaner substitutes for fossil fuel, yet their intermittency and unpredictability introduce significant challenges toward grid stability [1,2]. While energy systems are evolving and incorporating increasing shares of variable renewable energy resources, energy storage solutions will be important to maintain the balance between supply and demand and ensure reliability to accomplish a successful energy transition [3,4]. Energy storage will play one of the leading roles in this transition, allowing energy to be captured and provided later at needed times. Among the different storage systems developed, redox flow batteries (RFBs) have emerged as one of the most promising technologies for wide application because of their modular design, scalability, and flexibility [5,6]. Unlike traditional batteries, RFBs allow independent scaling of storage capacity and power output due to the separation of the energy storage and power conversion components [7,8]. This feature makes them extremely suitable for grid-scale applications and integration with renewable energy sources.
Among the variants of RFBs, vanadium-based RFBs (VRFBs) are the most commercially developed technology [9]. These VRFBs have demonstrated promising performance in terms of long cycle life, chemical stability, and grid-scale deployments [10,11]. However, VRFBs also have their own challenges, such as the energy density being limited by the solubility constraints of vanadium species, and their reliance on vanadium, a critical metal with price and supply uncertainties [12,13]. Limitations in the existing systems has led to researchers exploring other chemistries such as organic redox flow batteries (ORFBs). This includes using organic active species derived from abundant elements in the Earth’s crust to reduce some drawbacks in VRFBs. These molecules provide a sustainable solution for the substitution of inorganic compounds with environmentally more benign materials [14,15]. These systems provide multiple advantages, like higher solubility, tunability, and the potential for cost reduction, and hence, may lead to better sustainability and efficiency in energy storage [16].
In the area of organic redox chemistry, significant progress has been achieved in alkaline systems. Lin et al. [17] presented the synergy between 2,6-dihydroxyanthraquinone (2,6 DHAQ) as a negative electrolyte and ferricyanide ([Fe(CN6)]4−/[Fe(CN6)]3−) as a positive counterpart. This configuration enabled the achievement of power densities of 0.45 W cm−2 at room temperature and 0.7 W cm−2 at 45 °C. Other researchers have confirmed the effectiveness of anthraquinones and ferricyanides as excellent alternatives [18,19,20]. Other chemistries have furthered the state of the art in alkaline RFBs. Zou et al. [21] reported an S/Fe RFB with a remarkable volumetric capacity of 40.74 Ah L−1 that can achieve 99% coulombic efficiency and maintain an extremely low capacity fade of 0.0166% per cycle. Recent work has shown that Fe/Mn-based alkaline batteries have promise for sustainable energy storage. Shen et al. [22] reported a theoretical cell voltage of 1.43 V with close to ideal coulombic efficiency over 400 cycles; however, there was a gradual drop in voltage efficiency from 75.3% to 61.4%.
While experimental approaches provide critical insights, their high resource demands have shifted attention toward computational methods. Within these, computational fluid dynamics (CFD) has proved capable of modeling electrochemical devices successfully [23]. In this framework, numerous studies have been conducted to enhance mass transport using numerical simulations [24,25]. Pan et al. [26] applied a gradient decrease from the inlet to the outlet in channel width. The uniformity of active species was improved, which translated to higher power density and reduced polarization losses. Martinez-Lopez et al. [27] investigated the impact of electrode compression in convection using a 2D numerical model. Their findings revealed that applying 50% compression enhanced velocity profiles by 12.7%. Xu et al. [28] utilized a different distributor within a detailed 3D model. They evaluated the batteries both with and without a flow field and with serpentine and parallel patterns. They studied the performance metrics, overpotentials, pressure drops, and uniformity along the distributor, including a sensitivity analysis on the flow rate. Xu et al. also gave an overall efficiency metric for each geometry and then concluded that the serpentine flow pattern had a better performance and efficiency. Chu et al. [29] introduced a 3D numerical model utilizing 2,6-dihydroanthraquinone and ferrocyanide electrolytes, focusing on the effects of various electrode geometries, such as rectangular, trapezoidal, and sector shapes, on performance metrics like voltage, overpotentials, and efficiency. The sector-shaped electrode configuration showed the best results in mass transfer and power-based efficiency. Aparicio-Mauricio et al. [30] explored the influence of manifold design on electrolyte distribution, showing that uneven flow from dividing manifolds led to increased shear stress, higher pressure drops, and performance losses, including elevated voltage and current density deviations. Akuzum et al. [31] investigated the effects of channel obstructions and ramps on mass transport. They found that a ramped channel design enhances electrochemical performance by increasing electrolyte penetration at the electrode–membrane interface. However, adding obstructions worsened performance due to electrolyte bypassing. Notably, both tapered and obstructed channels result in lower pressure drops compared to an unmodified flow field. Messaggi et al. [32] studied the effects of obstruction placement on the flow channels. Results showed that electrolyte penetration was enhanced when the obstructions were placed on the channel side, whereas those placed at the electrode interface direct electrolyte better toward the channel outlet. Biomimetic flow fields inspired by leaf veins were proposed for VRFBs to improve electrolyte distribution by Liu et al. [33]. Among the designs, circular obstacles in the main channel achieved the best performance, with enhanced voltage efficiency and concentration uniformity, offering a promising approach for practical applications. The use of static mixers along variable compression has also been studied to improve performance on flow-through configurations [34,35]. The reactant distribution is shown to be enhanced, translating to cell voltage improvement, reduction of concentration overpotential and extended capacity.
Previous research has predominantly concentrated on flow-by configuration VRFBs, and the limited research on organic redox flow batteries has focused on clarifying the influence of electrode construction and manifold configuration. This works aims to expand the knowledge on the effects of introducing obstacles on flow-through configuration batteries. To address this gap, a two-dimensional stationary isothermal model is developed to analyze obstacle placement compared to a baseline design under different working conditions. The study evaluates the impact on discharge voltage, limiting current density, concentration uniformity, and pressure drop across different states of charge (SOC) and flow rates.

2. Model Description

A schematic representation of the ORFB is shown in Figure 1. The electrolytes are kept in two separate reservoirs. The ([Fe(CN6)]4−/[Fe(CN6)]3−) redox pair is stored in the posolyte tank, while 2,6 dihydroxyanthraquinone (2,6 DHAQ)/2,6 reduced DHAQ is stored in the negolyte tank. The electrolyte is circulated to the porous electrodes, which act as active sites for the electrochemical reactions, via peristaltic pumps during battery operation. An ion-selective membrane is placed between the electrodes, which allows only selected ions to flow through while preventing cross-contamination. The following electrochemical processes are taking place on the cathode and anode surfaces:
P o s i t i v e : F e ( C N ) 6 3 + e charge discharge F e ( C N ) 6 4
N e g a t i v e : 2.6 DHAQ + 2 e discharge charge 2.6 reDHAQ
A two-dimensional cell model has been developed for this work, which includes three domains: the cathode, the ion exchange membrane, and the anode. The key assumptions to reduce the complexity of the multiphysics behavior of the battery taken into consideration are:
  • The model is stationary.
  • Every property of the electrode and membrane is isotropic.
  • It is assumed that the entire cell is isothermal.
  • The electrolyte is considered incompressible.
  • There is no modeling of parasitic reactions.
  • Infinite dilute approximation is considered.
  • The membrane only permits K+ ions to pass through. All other ion crossover is disregarded.
The concept of obstacle placement compared to a base case is shown in Figure 2. The obstacles are rectangular prisms.

2.1. Governing Equations

2.1.1. Mass Transport

The conservation of mass for each species i is expressed by Equation (3):
t ( ε c i ) + · N i = S i
where c denotes the concentration of species i, ε denotes the electrode porosity, and S stands for the source term of each substance i, provided in Table 1.
The Nernst–Planck equation expresses the flux of charged species by diffusion, migration, and convection, represented as Ni:
N i = D i e f f c i z i u i c i F φ l + u c i
where D e f f is the effective diffusivity, z stands for the charge of species, u indicates ionic mobility, F represents the Faraday constant, and φ l is the liquid phase potential. u indicates the electrolyte velocity. Table 2 provides a summary of electrolyte characteristics.
D i e f f is obtained by means of the Bruggemann correlation:
D i e f f = ε 3 / 2 D i
The Nernst–Einstein equation is used to calculate ionic mobility:
u i = D i e f f R T
where R refers to the universal gas constant and T is the temperature.
Darcy’s Law is used for the convection term in porous electrodes:
u = K μ p
where K stands for the porous electrode permeability, μ is the dynamic viscosity of electrolyte, and p refers to pressure. The following is the Kozeny–Carman equation for permeability:
K = d f 2 ε 3 16 k c k ( 1 ε ) 2
where d f stands for the electrode fiber diameter and k c k is the Kozeny–Carman constant, listed with additional electrode parameters in Table 3.
Since only K + ions are permitted to pass through the membrane, the ion flux is determined by:
N K + = σ m e m F Δ φ m e m
where the potential and membrane conductivity are denoted by φ m e m and σ m e m , respectively.

2.1.2. Reaction Kinetics

The electrochemical reactions at the electrode surface are modeled using the Butler–Volmer model:
i r , p o s = i 0 , p o s a e x p ( 1 α p o s ) F η + R T e x p ( α p o s ) F η + R T
i r , n e g = i 0 , n e g a e x p ( 1 α n e g ) F η R T e x p ( α n e g ) F η R T
The positive and negative reaction rates are denoted by the terms i r , p o s and i r , n e g , respectively. Specific surface area of the porous electrode is denoted by a, η + and η are the positive and negative overpotential, charge transfer coefficients for positive and negative reactions are denoted by α p o s and α n e g , and exchange current density i 0 , p o s and i 0 , n e g are determined by:
i 0 , p o s = F k p o s ( c [ F e ( C N ) 6 ] 3 ) 1 α p o s ( c [ F e ( C N ) 6 ] 4 ) α p o s
i 0 , n e g = F k n e g ( c D H A Q 2 ) 1 α n e g ( c D H A Q 4 ) α n e g
where the reaction rate constants for the positive and negative sides are denoted by k p o s and k n e g , respectively.
The following formulas can be used to determine the overpotential of the positive and negative electrodes:
η + = φ s φ l E e q , +
η = φ s φ l E e q ,
where φ s represents the solid phase potential, and E e q , + and E e q , denote open circuit potentials of positive and negative reactions, respectively. These potentials can be obtained using the Nernst equation:
E e q , + = E e q , + + R T F l n c [ F e ( C N ) 6 ] 3 c [ F e ( C N ) 6 ] 4
E e q , = E e q , + R T 2 F l n c D H A Q 2 c D H A Q 4
where E e q , + and E e q , indicate the standard equilibrium potentials for the positive and negative side reactions, outlined in Table 4 alongside the other kinetic parameters.

2.1.3. Charge Conservation

The principle of electroneutrality requires that the total charge in the electrolyte equals zero, ensuring electrical neutrality:
i z i c i = 0
To satisfy charge conservation, current flow in the solid and liquid phases is coupled with electrochemical reactions through the following equation:
· i l = · i s = i r
which implies that the current leaving the electrolyte, i l , matches the current entering the electrode, i s , and collectively equals the electrochemical reaction rate, i r . The solid and liquid phase currents are expressed as follows:
i s = σ s e f f φ s
i l = F i z i N i
where σ s e f f represents the effective conductivity of the electrode, determined using the bulk conductivity of the electrode material, σ s , and calculated as follows:
σ s e f f = ( 1 ε ) 3 / 2 σ s

2.1.4. Boundary Conditions

Boundary conditions are defined using the x and y coordinates shown in Figure 1. At the external boundary of the anode, located at x = x 0 , the solid potential is set to zero by grounding this boundary:
φ s = 0 a t x = x 0
At the cell’s lower boundary, y = 0 , the species inlet flux is defined based on the flow rate Q, electrode width w e , and electrode thickness t e :
n · u = Q ε w e t e a t x 0 < x < x 1 a n d x 2 < x < x 3 y = 0
The species concentration at the electrode inlet remains constant:
c i = c i i n a t x 0 < x < x 1 a n d x 2 < x < x 3 y = 0
At all other boundaries, the mass flux is zero, with the exception of K + at the membrane–electrode interface:
n · N i = 0 a t y = 0 < y < y = h f o r x = ( x 0 , x 1 , x 2 , x 3 )
Here, h represents the height of the cell. At this location, a pressure outlet is defined:
p = p o u t a t x 0 < x < x 1 a n d x 2 < x < x 3 y = h
and the diffusive flux of all species is neglected:
D i e f f c i · n = 0 a t x 0 < x < x 1 a n d x 2 < x < x 3 y = h
At the membrane–electrode interface, charge continuity is maintained by defining the continuity of the electrolyte current density to the current density in the membrane:
n · i l = n · i l , mem a t y = 0 < y < y = h f o r x = ( x 1 , x 2 )
The top and bottom edges of both the membrane and electrodes are electrically insulated:
n · i l = 0 n · i s = 0 x 0 < x < x 3 y = 0 a n d y = h
A constant current density, controlled by the user, is applied to the external boundary of the positive electrode:
n · i s = i a v g a t x = x 3
All geometric and operational parameters are provided in Table 5.

2.2. Battery Performance Parameters

The concentration of each species is determined based on the user-defined parameter SOC, with the total concentrations on the positive and negative sides denoted as c 1 and c 0 , respectively:
c r e D H A Q = c 0 · S O C
c D H A Q = c 0 · ( 1 S O C )
c F e ( C N ) 6 4 = c 1 · ( 1 S O C )
c F e ( C N ) 6 3 = c 1 · S O C
A crucial factor in assessing species distribution is the uniformity factor, which measures how evenly a species is spread throughout the system. It quantifies the concentration variation of species i and is calculated as follows:
U i = 1 1 c i , avg 1 V c i c i , avg 2 d V
where c i , a v g denotes the average concentration of species i, and V represents the volume of the electrode.

2.3. Numerical Model

The simulation model was created using COMSOL Multiphysics 5.5. For the positive and negative electrodes, the Tertiary Current Distribution module was utilized, while the Secondary Current Distribution module was applied to the membrane. The fluid dynamics equations within the porous electrodes were addressed using the Brinkman equations, solved via a finite volume method. The PARDISO solver was selected, with a convergence criterion of 1 × 10 6 relative error to ensure high precision in the results. To optimize computational efficiency while maintaining accuracy, a mesh independence study was carried out, demonstrating that grid refinements beyond 3100 elements did not significantly alter the outcomes, as shown in Figure 3.

3. Results

3.1. Experimental Validation

To validate the numerical model, the results were compared with experimental data from Chu et al. [29], using a rectangular electrode with dimensions of 20 mm × 20 mm × 3.5 mm. The flow rate of the circulation pump was set to 60 ml min−1, with a current density of 40 mA cm−2, and a Nafion 212 cation exchange membrane was used as the ion-exchange medium. In these tests, the concentration of the positive active material was 200 mol m−3, while the concentration of the negative active material was 100 mol m−3. Figure 4 presents the cell potential during charge and discharge modes for both experimental and numerical models. The maximum error, which occurs at the endpoints of the charge-discharge curves, remained below 3.6%, with an average error of 1.45%, confirming the accuracy and reliability of the numerical model for subsequent analysis.

3.2. Case Studies

Two configurations were investigated to assess the influence of obstacle placement. Case 1 serves as the baseline, with no obstacles incorporated, providing a reference for evaluating the effects of the modified configuration. In Case 2, obstacles are introduced within the electrode structure. Two obstacles are introduced in each electrode, both obstructing 50% of the electrode thickness. One is located at the electrode outer boundary, and the second is placed at the electrode membrane interface. All obstacles have a dimension of 1 mm in electrolyte flow direction. The obstacle at the outer boundary is located at one-third of the electrode length, and the other one at two-thirds of the electrode length.
For both configurations, the concentration of active material in the positive and negative compartments is set to 400 and 200 mol m−3, respectively. Additionally, the concentrations of K+ and OH ions are maintained at 1000 mol m−3 on both sides.

3.3. Effects of SOC

Looking at Figure 5, discharge curves for Case 1 and 2 can be analyzed. The discharge curves reveal that Case 1 consistently exhibits a slightly higher voltage, 0.18% on average, compared to Case 2 across the entire range of state of charge. This behavior indicates that the baseline configuration benefits from a more direct flow path for the electrolyte, minimizing the hydraulic resistance and associated pressure drop. In contrast, the inclusion of obstacles in Case 2 leads to additional flow resistance and localized pressure gradients, which may impede uniform electrolyte penetration. However, when looking at the lower end of the SOC range, the difference between Case 1 and 2 decreases to 0.07%, see the zoom in Figure 5. While the obstacle configuration aims to enhance mass transport and mitigate concentration gradients by promoting better mixing, these potential advantages appear offset by the increased Ohmic polarization effects.
Despite the slightly lower discharge voltage observed in Case 2, this configuration demonstrates significant improvements in terms of concentration uniformity and minimum species concentration, as depicted in Figure 6a,b. Across the entire SOC range, Case 2 demonstrates a higher concentration uniformity factor compared to the baseline Case 1, which indicates enhanced distribution of active species within the electrode, see Figure 6a. For a SOC value of 0.2, Case 1 achieves a value of 88.7%, while for Case 2, this rises to 90.1%. The obstacles actually help spread electrolyte more evenly through the electrode, which prevents troublesome zones with low concentration. Though the advantage in concentration uniformity narrows as SOC increases, the design in Case 2 still provides better distribution than in Case 1, reinforcing stable electrochemical reactions across varying SOC conditions. One potential explanation for this trend is that at lower SOCs, the availability of reactants throughout the electrolyte is significantly limited, leading to pronounced concentration gradients across the cell. Under these conditions, reactant depletion near the electrode surface becomes a critical issue, particularly in regions where the flow is less effective at replenishing the consumed species. This is the operating point where obstacles can force portions of the electrolyte flow into regions that might otherwise receive minimal reactant supply, promoting a more uniform utilization of the electrode area. For higher SOCs, this becomes less critical and obstacles seem to be less effective. Furthermore, the higher minimum concentration observed in Figure 6b highlights the improved utilization of active species in Case 2. This suggests that the presence of obstacles mitigates stagnant regions and optimizes mass transport.
Better electrolyte distribution translates into higher limiting current density, see Figure 7. This improvement can be attributed to effectively reducing concentration gradients and mitigating localized depletion of active species. The more uniform concentration profile in Case 2, as illustrated in Figure 6a, corroborates this observation, as a well-distributed reactant supply ensures that the electrode regions farther from the inlet remain reactive even at higher current densities, even at near discharged state.

3.4. Effects of Flow Rate

A key consideration in evaluating the effectiveness of electrode design modifications is the associated pressure drop, see Figure 8, as it directly impacts the pumping energy required for system operation. In the case of obstacle placement (Case 2), a higher pressure drop is observed compared to Case 1, by an average of 14%. The increased hydraulic resistance in Case 2 can be attributed to the obstacles interrupting the electrolyte flow path, creating regions of recirculation, while these flow dynamics improve electrolyte penetration and mixing, as seen in the higher concentration uniformity and elevated limiting current density, they simultaneously impose a penalty on pumping efficiency. This trade-off highlights a challenge on optimizing the electrode structure to enhance electrochemical performance while minimizing the energy costs associated with fluid transport.
Regardless of the increase in pressure drop for Case 2, improvements in other metrics are obtained. Figure 9a shows the relation between concentration uniformity factor and flow rate for both configurations. Both curves exhibit a rising trend as the flow rate increases, reflecting enhanced reactant distribution across the electrode. However, there are notable differences in their behavior. Case 2 demonstrates a 1.8% improvement in concentration uniformity for the lowest flow rate of 10 mL min−1, and it maintains higher values across the studied range. However, this difference shrinks as the flow rate increases, down to 1.0% for 60 mL min−1. Looking at Figure 9b, benefits in minimum concentration are observed. Case 2 obtains higher minimum concentration values across the studied range, averaging a 15% increase, which translates to reduced local reactant depletion. Thanks to reducing local reactant depletion, limiting current density is raised, see Figure 10, by an average of 1% over the whole range. It is worth noting that when zooming in on Figure 10, even at lower flow rates, Case 2 consistently outperforms Case 1.

4. Conclusions

In this work, the effect of placing obstacles in the electrode structure of an organic redox flow battery on reactant distribution and local depletion is investigated. Two configurations were studied: Case 1 was designated as the control without obstacles, while Case 2 featured a design with obstacles to improve mass transport and uniformity.
Analyzing the discharge profiles, it is noted that Case 1 maintained a marginally higher voltage compared to Case 2 throughout the complete state-of-charge range; the average difference was about 0.18%. Such a voltage difference can be explained due to the fact that Case 1 had a simpler flow channel leading to lower hydraulic resistance and smaller pressure drops. Conversely, the obstacles in Case 2 introduced more resistance to the flow, which consequently generated localized pressure gradients. Additionally, obstacle placement means reducing electrode area, which translates to lower electronic conductivity and increased Ohmic polarization. However, this disadvantage was compensated at lower SOC values, where the gap in performance between the two configurations narrowed considerably. Although the discharge voltage was marginally lower in Case 2, there were marked improvements in the concentration uniformity and minimum concentration of active species with this configuration. The distribution of reactants became much more uniform with the obstacles placed in Case 2, especially at low levels of SOC, thereby preventing areas of low concentration that could damage the cell, while the benefit of concentration uniformity decreased with higher values of SOC, Case 2 still showed better distribution throughout the whole range.
The higher pressure drop of Case 2 is mainly due to the fact that the obstacles block electrolyte flow. This suggests a trade-off between the electrochemical performance and energy consumption: while the obstacles improved the mixing quality and reduced local depletion, the obstacles also created additional flow resistance that required more pumping energy.
Looking at the performance of the design with obstructions at different flow rates, Case 2 did better than Case 1, especially with lower flow rates. Although the advantage became smaller with an increase in flow rate, Case 2 still had better concentration uniformity and higher minimum concentration values within the studied range. Furthermore, having obstacles in Case 2 increased the limiting current density by 1% on average. Even at the lower flow rates, Case 2 showed better performance, indicating that the obstacles may enhance electrochemical efficiency at different operating conditions.
In conclusion, though the introduction of obstacles in the electrode structure increased the pressure drop, the enhanced concentration uniformity, limiting current density, and active species utilization bring huge benefits. It is indicated that obstacle placement could be a very promising strategy to optimize the performance of redox flow batteries. The introduction of strategically placed obstacles presents a simple yet scalable design modification that can be adapted to various operating conditions. Furthermore, the insights into the relationship between flow dynamics and electrochemical performance establish a foundation for future optimization efforts, including refining obstacle geometry and placement to balance improvements in mass transport with the associated pumping energy costs.

Author Contributions

Conceptualization, J.M.-L. and U.F.-G.; methodology, A.B.-A.; software, E.S.-D.; validation, J.M.-L. and U.F.-G.; investigation, U.F.-G.; resources, E.S.-D.; data curation, A.B.-A. and Í.O.-F.; Writing—Review and Editing, J.M.-L. and Í.O.-F.; visualization, A.B.-A.; supervision, E.S.-D.; project administration, U.F.-G.; funding acquisition, U.F.-G. All authors have read and agreed to the published version of the manuscript.

Funding

Government of the Basque Country: Elkartek CICe2022; KK-2022/00043 and ELKARTEK24/15, KK-2024/00035; CIEMAT: Energía eólica offshore para el ensayo y el desarrollo energético de energías renovables e hidrógeno verde; Government of the Basque Country: Mobility Lab Foundation; Government of the Basque Country: ITSAS-REM (IT1514-22).

Data Availability Statement

All data generated in the current study are available upon reasonable request to the corresponding author.

Acknowledgments

The authors are grateful for the support provided by SGIker of UPV/EHU. This research was developed under the frame of the Joint Research Laboratory on Offshore Renewable Energy (JRL-ORE).

Conflicts of Interest

Authors Eduardo Sánchez-Díez, Aitor Beloki-Arrondo and Íñigo Ortega-Fernández were employed by Basque Research and Technology Alliance (BRTA). The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Nomenclature

List of symbols
aSpecific surface area
cConcentration
DDiffusion coefficient
d f    Fiber diameter
EEquilibrium potential
FFaraday constant
hHeight
i 0 Exchange current density
i r Electrochemical reaction rate
KPermeability
K c k Kozeny–Carman constant
kReaction rate constant
NFlux of charged species
pPressure
QVolumetric flow rate
RConstant of ideal gases
SSource term
TTemperature
tTime
uMobility
uVelocity
wWidth
zSpecies charge
Greek
α Charge transfer coefficient
ε Electrode porosity
η Overpotential
φ Potential
σ Conductivity
μ Dynamic viscosity
Superscripts and subscripts
+Positive side
Negative side
Standard
a v g Average
eElectrode
e f f Effective
iSpecies
lLiquid
m e m Membrane
o u t Outlet
sSolid

Abbreviations

The following abbreviations are used in this manuscript:
CFDComputational fluid dynamics
RFBRedox flow battery
VRFBVanadium redox flow battery
ORFBOrganic redox flow battery

References

  1. Nguyen, T.; Savinell, R.F. Flow Batteries. Electrochem. Soc. Interface 2010, 19, 54. [Google Scholar] [CrossRef]
  2. Asiaban, S.; Kayedpour, N.; Ebneali Samani, A.; Bozalakov, D.; De Kooning, J.; Crevecoeur, G.; Vandevelde, L. Wind and solar intermittency and the associated integration challenges: A comprehensive review including the status in the Belgian power system. Energies 2021, 14, 630. [Google Scholar] [CrossRef]
  3. Mutz, M.; Perovic, M.; Gümbel, P.; Königer, V.; Taranovskyy, A.; Li, Y.; Beran, L.; Käfer, T.; Dröder, K.; Knoblauch, V.; et al. Toward a Li-Ion Battery Ontology Covering Production and Material Structure. Energy Technol. 2022, 11, 2200681. [Google Scholar] [CrossRef]
  4. Braff, W.A.; Mueller, J.M.; Trancik, J.E. Value of storage technologies for wind and solar energy. Nat. Clim. Chang. 2016, 6, 964–969. [Google Scholar] [CrossRef]
  5. Yu, J.; Fan, L.; Yan, R.; Zhou, M.; Wang, Q. Redox Targeting-Based Aqueous Redox Flow Lithium Battery. ACS Energy Lett. 2018, 3, 2314–2320. [Google Scholar] [CrossRef]
  6. Reed, D.; Thomsen, E.; Li, B.; Wang, W.; Nie, Z.; Koeppel, B.; Kizewski, J.; Sprenkle, V. Stack Developments in a kW Class All Vanadium Mixed Acid Redox Flow Battery at the Pacific Northwest National Laboratory. J. Electrochem. Soc. 2015, 163, A5211. [Google Scholar] [CrossRef]
  7. Arevalo-Cid, P.; Dias, P.; Mendes, A.; Azevedo, J. Redox flow batteries: A new frontier on energy storage. Sustain. Energy Fuels 2021, 5, 5366–5419. [Google Scholar] [CrossRef]
  8. Lourenssen, K.; Williams, J.; Ahmadpour, F.; Clemmer, R.; Tasnim, S. Vanadium redox flow batteries: A comprehensive review. J. Energy Storage 2019, 25, 100844. [Google Scholar] [CrossRef]
  9. Liu, Y.; Xu, J.; Lu, S.; Xiang, Y. Titanium Nitride Nanorods Array-Decorated Graphite Felt as Highly Efficient Negative Electrode for Iron-Chromium Redox Flow Battery. Small 2023, 19, e2300943. [Google Scholar] [CrossRef]
  10. Cunha, A.; Martins, J.; Rodrigues, N.; Brito, F. Vanadium redox flow batteries: A technology review. Int. J. Energy Res. 2014, 39, 889–918. [Google Scholar] [CrossRef]
  11. Bhattarai, A.; Whitehead, A.; Schweiss, R.; Scherer, G.; Skyllas-Kazacos, M.; Wai, N.; Nguyen, T.; Ghimire, P.; Oo, M.; Hng, H. Anomalous Behavior of Anion Exchange Membrane during Operation of a Vanadium Redox Flow Battery. ACS Appl. Energy Mater. 2019, 2, 1712–1719. [Google Scholar] [CrossRef]
  12. Skyllas-Kazacos, M.; Chakrabarti, B.; Hajimolana, Y.S.; Mjalli, F.; Saleem, M. Progress in Flow Battery Research and Development. J. Electrochem. Soc. 2011, 158, R55–R79. [Google Scholar] [CrossRef]
  13. Minke, C.; Kunz, U.; Turek, T. Techno-economic assessment of novel vanadium redox flow batteries with large-area cells. J. Power Sources 2017, 361, 105–114. [Google Scholar] [CrossRef]
  14. Leung, P.; Shah, A.A.; Sanz, L.; Flox, C.; Morante, J.R.; Xu, Q.; Mohamed, M.R.; Ponce de León, C.; Walsh, F.C. Recent developments in organic redox flow batteries: A critical review. J. Power Sources 2017, 360, 243–283. [Google Scholar] [CrossRef]
  15. Yu, L.C.; Luo, Y.C.; Feng, W.; Zhang, S.; Zhang, X. Fluorinated TEMPO: A new redox-active catholyte material for aqueous Zn-anode hybrid flow batteries. J. Mater. Chem. A 2023, 11, 18911–18921. [Google Scholar] [CrossRef]
  16. Sánchez-Díez, E.; Ventosa, E.; Guarnieri, M.; Trovò, A.; Flox, C.; Marcilla, R.; Soavi, F.; Mazur, P.; Aranzabe, E.; Ferret, R. Redox flow batteries: Status and perspective towards sustainable stationary energy storage. J. Power Sources 2021, 481, 228804. [Google Scholar] [CrossRef]
  17. Lin, K.; Chen, Q.; Gerhardt, M.R.; Tong, L.; Kim, S.B.; Eisenach, L.; Valle, A.W.; Hardee, D.; Gordon, R.G.; Aziz, M.J.; et al. Alkaline quinone flow battery. Science 2015, 349, 1529–1532. [Google Scholar] [CrossRef] [PubMed]
  18. Huang, S.; Zhang, H.; Salla, M.; Zhuang, J.; Zhi, Y.; Wang, X.; Wang, Q. Molecular engineering of dihydroxyanthraquinone-based electrolytes for high-capacity aqueous organic redox flow batteries. Nat. Commun. 2022, 13, 4746. [Google Scholar] [CrossRef] [PubMed]
  19. Peng, K.; Li, Y.; Tang, G.; Liu, Y.; Yang, Z.; Xu, T. Solvation regulation to mitigate the decomposition of 2,6-dihydroxyanthraquinone in aqueous organic redox flow batteries. Energy Environ. Sci. 2023, 16, 430–437. [Google Scholar] [CrossRef]
  20. Jin, S.; Fell, E.M.; Vina-Lopez, L.; Jing, Y.; Michalak, P.W.; Gordon, R.G.; Aziz, M.J. Near Neutral pH Redox Flow Battery with Low Permeability and Long-Lifetime Phosphonated Viologen Active Species. Adv. Energy Mater. 2020, 10, 2000100. [Google Scholar] [CrossRef]
  21. Zou, H.; Xu, Z.; Xiong, L.; Wang, J.; Fu, H.; Cao, J.; Ding, M.; Wang, X.; Jia, C. An alkaline S/Fe redox flow battery endowed with high volumetric-capacity and long cycle-life. J. Power Sources 2024, 591, 233856. [Google Scholar] [CrossRef]
  22. Shen, X.; Kellamis, C.; Tam, V.; Sinclair, N.; Wainright, J.; Savinell, R. An All-Soluble Fe/Mn-Based Alkaline Redox Flow Battery System. ACS Appl. Mater. Interfaces 2024, 16, 18686–18692. [Google Scholar] [CrossRef]
  23. Aramendia, I.; Fernandez-Gamiz, U.; Martinez-San-Vicente, A.; Zulueta, E.; Lopez-Guede, J.M. Vanadium Redox Flow Batteries: A Review Oriented to Fluid-Dynamic Optimization. Energies 2021, 14, 176. [Google Scholar] [CrossRef]
  24. Esan, O.C.; Shi, X.; Pan, Z.; Huo, X.; An, L.; Zhao, T. Modeling and Simulation of Flow Batteries. Adv. Energy Mater. 2020, 10, 2000758. [Google Scholar] [CrossRef]
  25. Sun, J.; Guo, Z.; Pan, L.; Fan, X.; Wei, L.; Zhao, T. Redox flow batteries and their stack-scale flow fields. Carbon Neutrality 2023, 2, 30. [Google Scholar] [CrossRef]
  26. Pan, L.; Xie, J.; Guo, J.; Wei, D.; Qi, H.; Rao, H.; Leung, P.; Zeng, L.; Zhao, T.; Wei, L. In-plane gradient design of flow fields enables enhanced convections for redox flow batteries. Energy Adv. 2023, 2, 2006–2017. [Google Scholar] [CrossRef]
  27. Martinez Lopez, J.; Aramendia, I.; Fernandez-Gamiz, U.; Sanchez-Diez, E.; Beloki, A.; Kurt, E.; Lopez-Guede, J.M. Computational Modeling of a 2D Vanadium Redox Flow Battery Cell. JOM 2024, 76, 130–140. [Google Scholar] [CrossRef]
  28. Xu, Q.; Zhao, T.S.; Leung, P.K. Numerical investigations of flow field designs for vanadium redox flow batteries. Appl. Energy 2013, 105, 47–56. [Google Scholar] [CrossRef]
  29. Chu, F.; Su, M.; Xiao, G.; Tan, Z.; Yang, G. Analysis of Electrode Configuration Effects on Mass Transfer and Organic Redox Flow Battery Performance. Ind. Eng. Chem. Res. 2022, 61, 2915–2925. [Google Scholar] [CrossRef]
  30. Aparicio-Mauricio, G.; Jiménez-Lima, L.; Rivero, E.P.; Cruz-Díaz, M.R. 3D modeling and simulation of an alkaline flow battery considering the tertiary current distribution on the electrodes: Anthraquinone-ferro/ferricyanide as organic electroactive species. J. Power Sources 2023, 582, 233533. [Google Scholar] [CrossRef]
  31. Akuzum, B.; Alparslan, Y.C.; Robinson, N.C.; Agar, E.; Kumbur, E.C. Obstructed flow field designs for improved performance in vanadium redox flow batteries. J. Appl. Electrochem. 2019, 49, 551–561. [Google Scholar] [CrossRef]
  32. Messaggi, M.; Gambaro, C.; Casalegno, A.; Zago, M. Development of innovative flow fields in a vanadium redox flow battery: Design of channel obstructions with the aid of 3D computational fluid dynamic model and experimental validation through locally-resolved polarization curves. J. Power Sources 2022, 526, 231155. [Google Scholar] [CrossRef]
  33. Liu, Y.; Huang, Z.; Xie, X.; Liu, Y.; Wu, J.; Guo, Z.; Huang, Q. Redox flow battery:Flow field design based on bionic mechanism with different obstructions. Chem. Eng. J. 2024, 498, 155663. [Google Scholar] [CrossRef]
  34. Gurieff, N.; Keogh, D.F.; Timchenko, V.; Menictas, C. Enhanced Reactant Distribution in Redox Flow Cells. Molecules 2019, 24, 3877. [Google Scholar] [CrossRef] [PubMed]
  35. Gurieff, N.; Keogh, D.F.; Baldry, M.; Timchenko, V.; Green, D.; Koskinen, I.; Menictas, C. Mass Transport Optimization for Redox Flow Battery Design. Appl. Sci. 2020, 10, 2801. [Google Scholar] [CrossRef]
  36. Dhillon, A. Multiphysics Modelling of an Alkaline All-Iron, All-Soluble Aqueous Redox Flow Battery. Master’s Thesis, University of Waterloo, Waterloo, ON, Canada, 2021. [Google Scholar]
  37. Chen, Z.; Yu, W.; Liu, Y.; Zeng, Y.; He, Q.; Tan, P.; Ni, M. Mathematical modeling and numerical analysis of alkaline zinc-iron flow batteries for energy storage applications. Chem. Eng. J. 2021, 405, 126684. [Google Scholar] [CrossRef]
  38. Wang, Q.; Qu, Z.G.; Jiang, Z.Y.; Yang, W.W. Experimental study on the performance of a vanadium redox flow battery with non-uniformly compressed carbon felt electrode. Appl. Energy 2018, 213, 293–305. [Google Scholar] [CrossRef]
  39. Knehr, K.W.; Agar, E.; Dennison, C.R.; Kalidindi, A.R.; Kumbur, E.C. A Transient Vanadium Flow Battery Model Incorporating Vanadium Crossover and Water Transport through the Membrane. J. Electrochem. Soc. 2012, 159, A1446. [Google Scholar] [CrossRef]
  40. You, D.; Zhang, H.; Chen, J. A simple model for the vanadium redox battery. Electrochim. Acta 2009, 54, 6827–6836. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the numerical model.
Figure 1. Schematic representation of the numerical model.
Batteries 11 00029 g001
Figure 2. Schematic representation of the (a) typical no-obstacle flow-through configuration and (b) the proposed obstacle configuration.
Figure 2. Schematic representation of the (a) typical no-obstacle flow-through configuration and (b) the proposed obstacle configuration.
Batteries 11 00029 g002
Figure 3. Grid dependency study.
Figure 3. Grid dependency study.
Batteries 11 00029 g003
Figure 4. Results of experimental validation [29].
Figure 4. Results of experimental validation [29].
Batteries 11 00029 g004
Figure 5. Discharge curves obtained for a flow rate of 60 mL min−1 and 50 mA cm−2.
Figure 5. Discharge curves obtained for a flow rate of 60 mL min−1 and 50 mA cm−2.
Batteries 11 00029 g005
Figure 6. The (a) concentration uniformity factor and (b) minimum concentration of [Fe(CN6)]3− in the positive electrode for a flow rate of 60 mL min−1 and 50 mA cm−2.
Figure 6. The (a) concentration uniformity factor and (b) minimum concentration of [Fe(CN6)]3− in the positive electrode for a flow rate of 60 mL min−1 and 50 mA cm−2.
Batteries 11 00029 g006
Figure 7. Limiting current density for a flow rate of 60 mL min−1 and 50 mA cm−2.
Figure 7. Limiting current density for a flow rate of 60 mL min−1 and 50 mA cm−2.
Batteries 11 00029 g007
Figure 8. Pressure drop values for a State of Charge of 50% and 50 mA cm−2 in discharge mode.
Figure 8. Pressure drop values for a State of Charge of 50% and 50 mA cm−2 in discharge mode.
Batteries 11 00029 g008
Figure 9. The (a) concentration uniformity factor and (b) minimum concentration of [Fe(CN6)]3− in the positive electrode for a State of Charge of 50% and 50 mA cm−2.
Figure 9. The (a) concentration uniformity factor and (b) minimum concentration of [Fe(CN6)]3− in the positive electrode for a State of Charge of 50% and 50 mA cm−2.
Batteries 11 00029 g009
Figure 10. Limiting current density for a State of Charge of 50% and 50 mA cm−2.
Figure 10. Limiting current density for a State of Charge of 50% and 50 mA cm−2.
Batteries 11 00029 g010
Table 1. Source terms of species.
Table 1. Source terms of species.
SpeciesPositive ElectrodeNegative Electrode
[Fe(CN6)]4− i r / F -
[Fe(CN6)]3− i r / F -
DHAQ2−- 2 i r / F
DHAQ4−- 2 i r / F
Table 2. Properties of electrolyte.
Table 2. Properties of electrolyte.
QuantitySymbolValueReferences
Diffusivity of [Fe(CN6)]4− D [ F e ( C N ) 6 ] 4 8.2 × 10 10 m2 s−1[36]
Diffusivity of [Fe(CN6)]3− D [ F e ( C N ) 6 ] 3 8.6 × 10 10 m2 s−1[36]
Diffusivity of DHAQ2− D D H A Q 2 4.8 × 10 10 m2 s−1[17]
Diffusivity of DHAQ4− D D H A Q 4 4.8 × 10 10 m2 s−1[17]
Diffusivity of K+ D K + 1.96 × 10 9 m2 s−1[37]
Diffusivity of OH D O H 5.2 × 10 9 m2 s−1[37]
Table 3. Electrode parameters.
Table 3. Electrode parameters.
QuantitySymbolValueReferences
Non-compressed electrode porosity ε 0.895[38]
Non-compressed electrode specific surface areaa 3.5 × 10 4 m2 m−3Fitted
Non-compressed electrode conductivity σ s 66.7 S m−1[39]
Kozeny–Carman constant K c k 4.28[40]
Table 4. Kinetic parameters.
Table 4. Kinetic parameters.
QuantitySymbolValueReferences
Standard equilibrium potential for positive reaction E e q , + 0.33 V[29]
Standard equilibrium potential for negative reaction E e q , −0.71 V[29]
Cathodic transfer coefficient α p o s 0.5[17]
Anodic transfer coefficient α n e g 0.5[17]
Rate constant for positive reaction k p o s 6 × 10 7 m s−1[29]
Rate constant for negative reaction k n e g 7 × 10 7 m s−1[29]
Table 5. Geometrical and operation parameters.
Table 5. Geometrical and operation parameters.
QuantitySymbolValueReferences
Electrode lengthh0.04 m[38]
Non-compressed electrode thickness t e 0.006 m[38]
Electrode width w e 0.04 m[38]
Membrane thickness t m 0.000183 m[38]
Flow rateQ60 mL min−1-
TemperatureT298 K-
Outlet pressure p o u t 0 Pa-
State of Charge S O C 50%-
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Martínez-López, J.; Fernández-Gamiz, U.; Sánchez-Díez, E.; Beloki-Arrondo, A.; Ortega-Fernández, Í. Enhancing Mass Transport in Organic Redox Flow Batteries Through Electrode Obstacle Design. Batteries 2025, 11, 29. https://doi.org/10.3390/batteries11010029

AMA Style

Martínez-López J, Fernández-Gamiz U, Sánchez-Díez E, Beloki-Arrondo A, Ortega-Fernández Í. Enhancing Mass Transport in Organic Redox Flow Batteries Through Electrode Obstacle Design. Batteries. 2025; 11(1):29. https://doi.org/10.3390/batteries11010029

Chicago/Turabian Style

Martínez-López, Joseba, Unai Fernández-Gamiz, Eduardo Sánchez-Díez, Aitor Beloki-Arrondo, and Íñigo Ortega-Fernández. 2025. "Enhancing Mass Transport in Organic Redox Flow Batteries Through Electrode Obstacle Design" Batteries 11, no. 1: 29. https://doi.org/10.3390/batteries11010029

APA Style

Martínez-López, J., Fernández-Gamiz, U., Sánchez-Díez, E., Beloki-Arrondo, A., & Ortega-Fernández, Í. (2025). Enhancing Mass Transport in Organic Redox Flow Batteries Through Electrode Obstacle Design. Batteries, 11(1), 29. https://doi.org/10.3390/batteries11010029

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop