“…Thus, for the sake of simplicity and primarily to capture the underlying physics behind pH gradient, we have neglected any nonlinear effects such as finite ion size, a dielectric decrease, and temperature-dependent parameters. 31–33 The conservation of ionic species follows the Nernst-Planck equation, which, in a steady-state can be expressed as
Here, c i , N i , z i , D i respectively, represent the concentration, flux, valence, and diffusivities of i th, i = 1, 2, 3, and 4 ionic species, respectively, and R , F , and T are the universal gas constant, Faraday constant, and absolute temperature, respectively. The electric potential ϕ is governed by the local space charge density of the electrolyte solution and it follows the Poisson equation:
where, j is a binary variable, value of which is either 1 or 0 depending upon the region within the PEL or outside the PEL, respectively.…”