RESUMO
Understanding the dynamics of electric-double-layer (EDL) charging in porous media is essential for advancements in next-generation energy storage devices. Due to the high computational demands of direct numerical simulations and a lack of interfacial boundary conditions for reduced-order models, the current understanding of EDL charging is limited to simple geometries. Here, we present a network model to predict EDL charging in arbitrary networks of long pores in the Debye-Hückel limit without restrictions on EDL thickness and pore radii. We demonstrate that electrolyte transport is described by Kirchhoff's laws in terms of the electrochemical potential of charge (the valence-weighted average of the ion electrochemical potentials) instead of the electric potential. By employing the equivalent circuit representation suggested by these modified Kirchhoff's laws, our methodology accurately captures the spatial and temporal dependencies of charge density and electric potential, matching results obtained from computationally intensive direct numerical simulations. Our network model provides results up to six orders of magnitude faster, enabling the efficient simulation of a triangular lattice of five thousand pores in 6 min. We employ the framework to study the impact of pore connectivity and polydispersity on electrode charging dynamics for pore networks and discuss how these factors affect the time scale, energy density, and power density of capacitive charging. The scalability and versatility of our methodology make it a rational tool for designing 3D-printed electrodes and for interpreting geometric effects on electrode impedance spectroscopy measurements.
RESUMO
Recent experimental studies have utilized AC electric fields and electrochemical reactions in multicomponent electrolyte solutions to control colloidal assembly. However, theoretical investigations have thus far been limited to binary electrolytes and have overlooked the impact of electrochemical reactions. In this study, we address these limitations by analyzing a system with multicomponent electrolytes, while also relaxing the assumption of ideally blocking electrodes to capture the effect of surface electrochemical reactions. Through a regular perturbation analysis in the low-applied-potential regime, we solve the Poisson-Nernst-Planck equations and obtain effective equations for electrical potential and ion concentrations. By employing a combination of numerical and analytical calculations, our analysis reveals a significant finding: electrochemical reactions alone can generate asymmetric rectified electric fields (AREFs), i.e., time-averaged, long-range electric fields, even when the diffusivities of the ionic species are equal. This finding expands our understanding beyond the conventional notion that AREFs arise solely from diffusivity contrast. Furthermore, we demonstrate that AREFs induced by electrochemical reactions can be stronger than those resulting from asymmetric diffusivities. Additionally, we report the emergence of asymmetric rectified concentration fields (ARCFs), i.e., time-averaged, long-range concentration fields, which supports the electrodiffusiophoresis mechanism of colloidal assembly observed in experiments. We also derive analytical expressions for AREFs and ARCFs, emphasizing the role of imbalances in ionic strength and charge density, respectively, as the driving forces behind their formation. The results presented in this article advance the field of colloidal assembly and also have implications for improved understanding of electrolyte transport in electrochemical devices.
RESUMO
Porous electrodes are found in energy storage devices such as supercapacitors and pseudocapacitors. However, the effect of electrode-pore-size distribution on their energy storage properties remains unclear. Here, we develop a model for the charging of electrical double layers inside a cylindrical pore for arbitrary pore size. We assume small applied potentials and perform a regular perturbation analysis to predict the evolution of electrical potential and ion concentrations in both the radial and axial directions. We validate our perturbation model with direct numerical simulations of the Poisson-Nernst-Planck equations, and obtain quantitative agreement between the two approaches for small and moderate potentials. Our analysis yields two main characteristic features of arbitrary pore size: (i) a monotonic decrease of the charging timescale with an increase in relative pore size (pore size relative to Debye length); (ii) large potential changes for overlapping double layers in a thin transition region, which we approximate mathematically by a jump discontinuity. We quantify the contributions of electromigration and charge diffusion fluxes, which provide mechanistic insights into the dependence of charging timescale and capacitance on pore size. We develop a modified transmission circuit model that captures the effect of arbitrary pore size and demonstrate that a time-dependent transition-region resistor needs to be included in the circuit. We also derive phenomenological expressions for average effective capacitance and charging timescale as a function of pore-size distribution. We show that the capacitance and charging timescale increase with smaller average pore sizes and with smaller polydispersity, resulting in a gain of energy density at a constant power density. Overall, our results advance the mechanistic understanding of electrical-double-layer charging.
