Data on flow cell optimization for membrane-based electrokinetic energy conversion.

This article elaborates on the design and optimization of a specialized flow cell for the measurement of direct conversion of pressure into electrical energy (Electrokinetic Energy Conversion, EKEC) which has been presented in Østedgaard-Munck et al. (2017) [1]. Two main flow cell parameters have been monitored and optimized: A) the hydraulic pressure profile on each side of the membrane introduced by pumps recirculating the electrolyte solution through the flow fields and B) the electrical resistance between the current collectors across the combined flow cell. The latter parameter has been measured using four-point Electrochemical Impedance spectroscopy (EIS) for different flow rates and concentrations. The total cell resistance consists of contributions from different components: the membrane [Formula: see text], anode charge transfer [Formula: see text], cathode charge transfer [Formula: see text], and ion diffusion in the porous electrodes [Formula: see text]. The intrinsic membrane properties of Nafion 117 has been investigated experimentally in LiI/I2 solutions with concentrations ranging between 0.06 and 0.96 M and used to identify the preferred LiI/I2 solution concentration. This was achieved by measuring the solution uptake, internal solution concentration and ion exchange capacity. The membrane properties were further used to calculate the transport coefficients and electrokinetic Figure of merit in terms of the Uniform potential and Space charge models. Special attention has been put on the streaming potential coefficient which is an intrinsic property.

Specifications Table

Value of the data

Flow cell optimization with respect to flow patterns, electrical resistance, and general cell design.

Data analysis method and results for Electrochemical Impedance Spectroscopy on operating flow cells.

Optimization of the solution concentration including measurements of:

Solution uptake

Internal solution concentration

Ion exchange capacity

Modelling of membrane transport properties using two classical pore models (Uniform potential and Space charge models).

Data

Flow cell pressure profile

The flow regime in the graphite plate channels (see Fig. 6) was evaluated considering the Reynolds number, Re, which is calculated as:

Fig. 6

Exploded view of the flow cell and cross sections of the inlet/outlet and pressure indicator ports.

Assuming: (i) 12 channels with equal fluid velocity (in m s-1) where is the circulation flow rate (in m3 s-1) and is the cross sectional area (in m2), (ii) solution density ~ 103 kg m-3 and viscosity ~ 10-3 Pa s (assumed equal to that of water at 25 °C). = 10-3 m is the hydraulic equivalent diameter of a single square channel (width mm and height mm) having wetted perimeter . For the conditions used in the present work Re = 16.7–267 for = 0.2–3.2 mL s-1 which ensures that all experiments were conducted in laminar flow conditions in the flow field channels.

The pressure drop in each half cell was measured with two pressure indicators positioned between the internal (inlet) and the outlet port. Streaming potential coefficient measurements (at zero current density) where performed both with co-flow and counter flow conditions (see Fig. 1a) and were used to validate the pressure differences measured in the system. The data can be seen in Fig. 1b scaled with the logarithmic mean pressure difference.

Fig. 1

a) Schematic of the counter-flow and co-flow conditions for the fluid flow in the interdigitated channels showing the profiles of the qualitative pressures (solid lines) and transmembrane pressure difference (dashed lines. b) Streaming potential coefficients versus the circulation flow rate in counter and co-flow scaled by the logarithmic mean pressure difference.

Flow cell electrochemical impedance spectroscopy

The Bode plots for varying concentrations and circulation flow rates can be seen in Fig. 2a and b, respectively. Corresponding Nyquist plots can be seen in ref [1].

Fig. 2

a) Bode plots for varying LiI/I2 concentrations (0.06–0.96 M) at fixed qcirc=2.9 mL s-1. b) Bode plot for varying circulation flow rates (0–3.2 mL s-1) at fixed cext=0.26 M LiI/I2.

Membrane properties

See Table 1.

Table 1

Experimental data for solution uptake (s), ion exchange capacity (iec), internal solution concentration of LiI saturated with I2 and immobile charge density (), and calculated internal solution concentrations (, ) for two different immobile charge densities ( and ).

c_ext(M)	s(g_solgdry pol−1)	iec(meqgpol−1)	X˜_t (mol m-3)	c˜_exp (M)	c˜_UP,sa(M)	c˜_UP,tb(M)
0.06	0.20	–	–	0.0028⁎	0.00084	0.0011
0.12	0.23	–	–	0.0051⁎	0.0034	0.0042
0.16	0.19	–	–	0.0062	0.0060	0.0075
0.21	0.24	–	–	0.012	0.010	0.013
0.26	0.24	0.80	3528	0.017	0.016	0.020
0.51	0.22	0.80	3445	0.045	0.060	0.075
0.74	0.18	0.83	3539	0.062	0.125	0.154
0.96	0.23	0.78	3189	0.144	0.206	0.252
Average	0.21±0.02	0.80±0.02	3425±163	–	–	–

Notes: ⁎ Conductivity below detection limit. a, b The immobile charge densities used for the Donnan equilibrium calculations were = 4270 and = 3400 mol m-3, respectively. The subscripts “s” and “t” indicates whether the ion exchange capacity used in the calculations comes from the manufacturer “specifics” (iec ~ 0.91 meq ) or measured with titration (iec ~ 0.80 meq ).

Electrical potential profiles and Donnan equilibrium

To have insights on the optimal solution concentration for EKEC processes, the internal solution concentration was determined theoretically considering two classical pore models (Uniform potential and Space charge models). Generally a nanocapillary of radius R with charged walls (, in C m-2) in equilibrium with an external electrolyte solution (with concentration ) develops cation and anion concentration profiles along the radial direction (r) which can be calculated from the Boltzmann distribution as:where (dimensionless) and are the valence of the ith ion and the thermal voltage ( (in V) with Rg, T and F being the gas constant (in J mol-1 K-1), absolute temperature (in K) and Faraday constant (in C mol-1). Here the ions are modelled as point charges and no effects related to the ion size and deviation on the dielectric constant of the medium are taken in consideration. The electrical potential (in V) can be calculated from the Poisson equation in cylindrical coordinates for a 1:1 electrolyte as:where is the permittivity of the medium (6.9×10-10 F m-1). Eq. (3) must be solved with the proper boundary conditions which, for the case of interest here, are constant charge density on the pore wall (r=R) and symmetry at r=0:

The solution of the Poisson equation (Eq. 3) is central in the “Space charge” model theory and the calculation of the internal solution concentration (as integral of on the radial direction) will be referred to as “SC” [2], [3], [4] In the limiting case of constant electrical potential profile inside the pore (which is strictly valid when the pore radius is smaller than the characteristic Debye length ) the internal solution concentration ( can be determined as:which coincide with the one calculated from Donnan equilibrium [5], [6]. In Eq. (5) is the total ion concentration inside the membrane. The approximation of constant concentration profile inside the pore (or overlapping electrical double layers) is the main assumption in the “Uniform potential” model and the internal concentration calculated with Eq. (5) will be referred to as “UP” [2], [3], [7], [8]

In Table 1 determined with and are given as and , respectively. The influence of on is shown in Fig. 3 together with derived from the SC and UP models (see ref. [4], [8]) with equal to 3400 and 4270 mol m-3 and pore radii of 2.3, 2.5 and 2.7 nm, respectively. The pore radii used in these calculations have been chosen to obtain a good fit with the experimental streaming potential data at infinite dilute conditions (see next section) and are similar to the one measured from SAXS experiments in Nafion [9].

Fig. 3

Internal LiI/I2 solution concentration () as function of the external solution concentration () measured and determined based on two different models; Space charge (SC) (green curves) and Uniform potential (UP) (yellow curves) with different immobile charge densities () and pore radii (). Experimentally determined internal solution concentrations are represented by solid symbols.

Systematically are lower than calculated from the Donnan equilibrium which means that the effective immobile charge density might be higher than the one estimated from titration (i.e. ). Alternatively small systematic errors on determining the solution uptake would reflect in relatively large changes in and . Fig. 3 clearly shows that the coion exclusion in the conditions adopted in the present work is close to the ideal case with the electrical double layers totally overlapping (UP model), while for thin electrical double layers compared to the pore radius (SC model) highly overestimate .

Streaming potential and figure-of-merit

Onsager proposed, based on the fundamental theorem of non-equilibrium thermodynamics, that, for sufficiently slow processes near equilibrium, the fluxes can be expressed as a linear combination of all the conjugated and non-conjugated driving forces [10], [11]. Using this framework the isothermal transport of an ideal electrolyte across a membrane can be described, in dimensionless formulation, as:in which: , and represent the (dimensionless) volumetric solution, ion and electrical current flux, respectively; while the (dimensionless) driving forces: , and are the total pressure, chemical potential and electrical potential gradient, respectively. Finally and represent the (dimensionless) flux of positive and negative ions, respectively. In general the cross-coefficients are concentration dependent and proven to be Onsager symmetric, hence . In the absence of a concentration difference, i.e. =0 and =const, the system can be greatly simplified and in the integral form it can be rewritten as:where Jv and I are the (dimensionless) volumetric solution flux and current density, respectively.

From Eq. (9) and Eq. (10) the most accessible experimental transport coefficients are the hydraulic permeability (), streaming potential coefficient () and membrane conductivity ():in which the (dimensional) cross-coefficients K can be calculated from the dimensionless counterparts as:where is the reference concentration (here 1 mol m-3 as in ref. [4]) and D is the average ion diffusion coefficient described as where and are the anion and cation diffusion coefficients, respectively.[7] For the calculation reported in this article the diffusion coefficients for the cation (Li+) and anion () were = 1.03·10-9 m2 s-1 [12] and = 1.10·10-9 m2 s-1, respectively.

determined based on the SC and UP models are shown in Fig. 4a and b, respectively. It can be seen that for the SC model at a fixed geometry and immobile charge density, decreases with increasing external solution concentration while for the UP model is more or less concentration independent. The experimentally determined streaming potential coefficients are similar to the calculated ones from the UP model with a =3400 mol m-3 and a pore radius of 2.5 nm.

Fig. 4

Streaming potential coefficient () calculated from the phenomenological transport coefficients for different immobile charge densities () and pore radii (). The phenomenological transport coefficients were determined using two different models; a) the Space charge (SC) and b) Uniform potential (UP), and plotted as function of the external LiI/I2 concentration . The experimentally determined streaming potential coefficients are represented by the solid symbols.

In the data reported in Fig. 4 it has been assumed that the equilibrium LiI+I2⇌LiI3 was forced completely towards the triiodide for all used LiI/I2 solutions and hence the triiodide diffusion coefficient was used in the calculations. A sensitivity analysis of the influence of the diffusion coefficient on the calculated was performed by changing the anion diffusion coefficient to the iodide = 2.00·10-9 m2 s-1 [13]. is concentration dependent [13] and the highest value (i.e. at infinite dilution) was chosen to obtain the largest deviation from. The percentage deviation between calculated with the two different anion diffusion coefficients were in most cases below 2% and increased at the most up to ~ 5% (for UP) and ~ 8% (for SC) for the most concentrated LiI/I2 solution.

The experimentally accessible transport properties described in Eqs. (11), (12), (13) are used to determine the Figure of merit () (Eq. (17)) according to the formulation given by Bentien et al. [14], which was derived using the framework proposed by Osterle and co-workers [2], [3], for the electrokinetic energy conversion to calculate maximum conversion efficiency (:

In Fig. 5 the calculated using both the SC and UP models are shown for different s and pore radii.

Fig. 5

Figure of merit () shown as function of the external LiI/I2 solution concentration for different immobile charge densities () and pore radii (). was calculated from the phenomenological transport coefficients derived considering the two different models; a) Space charge (SC) and b) Uniform potential (UP).

In the present work the UP model gives a better description of the experimentally determined and with respect to the SC model. It is important to recall that the quantitative estimation of calculated from the UP model will in most cases of interest (pore radii > 0.5 nm) highly overestimate [7].

Experimental design, materials and methods

Cell design

The flow cell described here was purposely designed and produced. The geometrical active area of the mounted membrane was 25 cm2. Fig. 6 depicts the exploded view of the flow cell and the graphite (Royal Elite New Energy Science & Technology, China) interdigitated flow field blocks where the ports for inlet, outlet and in-pattern pressure monitoring are shown. The inlet and outlet ports were 1/4 in. NPT threaded and coupled with straight 1/4 in. NPT to 1/4 in. compression fittings (Swagelok, PTFE). The ports for the pressure indicators were 1/8 in. threaded and 1/8 in. NPT to 1/8 in. compression fitting (Swagelok, PTFE) was used for the coupling with the instruments. The larger dimension for the fluid inlet and outlet was chosen to decrease the pressure drops in the external hydraulic circuit. Each cell (high pressure side and low pressure side) consisted of a stainless steel endplate, an insulator, a copper current collector, and a graphite block with machined flow field. The hydraulic sealing was ensured using an outer O-ring and an inner O-ring (M Seals, Viton) on the high and low pressure side, respectively. The membrane acted as an electrical insulator between the two electrodes. Two carbon paper sheets (FuelCellStore: Toray Carbon Paper, thickness: 190 µm) were placed in each graphite block between the flow field and the membrane. The flow field was lowered 0.3 mm with respect to the graphite block surface to ensure adequate space to host two carbon paper sheets (final compression of the carbon paper ~ 26% with an estimated porosity under compression ~ 74% [15], [16]). The interdigitated flow field, seen in Fig. 6, consisted of a single wall separating inlet and outlet streams forming 12 “dead-end” channels. This design forced the solution over the wall, through the porous carbon paper electrode, and close to the membrane where the electrokinetic phenomena took place.

Electrochemical impedance spectroscopy

For the Electrochemical Impedance Spectroscopy (EIS) analysis the used equivalent circuit and a representative spectrum can be seen in ref. [1] Figs. 1d and e, respectively. The equation of the equivalent circuit is:

Due to the imperfection of the capacitive elements in the flow cell, constant phase elements (CPE) have been used in the model. The exponent n measures how far the element (with C, A, and D representing the anode, cathode and diffusion limited element, respectively) is from a perfect capacitor (i.e. n = 1) or a resistor (n = 0). For this analysis n was allowed to vary between 0.7 and 1 [17]. In Table 2 all model parameters for the fits are given.

Table 2

Model parameters used to fit the experimental data regarding the flow cell resistance elements (R), capacitance elements (C) and constant phase element factor (0.7 ≤ n ≥ 1.0). These parameters describe an imperfect capacitor for membrane, anode, cathode, and diffusion limited parts of the flow cell electrical circuit with i = “mem”, “A”, “C”, and “D”, respectively. Parameters are shown for varying external LiI/I2 solution concentration (cext) and for varying circulation flow rate (qcirc).

Varying external concentration at fixedqcirc= 2.9 mL s-1
cext(M)	Rmem(Ω)	RA(Ω)	RC(Ω)	RD(Ω)	CA(F)	CC(F)	CD(F)	nA(-)	nC(-)	nD(-)
0.06	0.057	0.281	0.128	0.107	2.72	0.60	460	0.7	0.7	0.840
0.12	0.055	0.163	0.102	0.046	5.00	0.87	1509	0.7	0.7	0.925
0.16	0.054	0.121	0.092	0.036	6.54	0.93	2255	0.7	0.7	0.927
0.21	0.053	0.099	0.076	0.018	7.79	1.22	4400	0.7	0.7	1
0.26	0.050	0.086	0.063	0.022	8.00	1.28	6820	0.7	0.7	0.979
0.51	0.043	0.042	0.032	0.006	8.48	3.77	10,181	0.844	0.736	0.7
0.74	0.041	0.029	0.036	0	8.19	7.03	–	0.729	0.7	–
0.96	0.043	0.029	0.029	0	8.50	41.00	1	0.7	0.7	1
Varying externalqcircat fixedcext= 0.26 M LiI/I2
qcirc(mL s-1)	Rmem(Ω)	RA(Ω)	RC(Ω)	RD(Ω)	CA(F)	CC(F)	CD(F)	nA(-)	nC(-)	nD(-)
0	0.064	0.083	0.047	0.614	37.93	3.15	449	0.7	0.7	1
0.2	0.063	0.037	0.046	0.016	37.90	3.04	110	0.7	0.7	1
0.7	0.061	0.039	0.040	0.020	21.98	3.65	1662	0.7	0.7	1
1.3	0.059	0.035	0.031	0.013	12.81	4.25	2062	0.7	0.7	1
2.1	0.056	0.029	0.029	0	7.00	3.01	–	0.7	0.7	–
2.9	0.056	0.027	0.027	0	3.20	5.02	–	0.7	0.7	–
3.2	0.057	0.027	0.027	0	3.00	4.99	–	0.7	0.7	–

Solution uptake and internal solution concentration

Nafion 117 samples in the H+ form were dried at 50 °C and the dry masses were measured. Afterwards the samples were equilibrated overnight in LiI/I2 solutions with concentrations in the range 0.06 to 0.96 M. Before measuring the wet membrane masses the membranes were carefully blotted with tissue paper to remove excess solution from the surfaces. The solution uptake (s) was determined for all LiI/I2 concentrations as:where and are the wet and dry membrane masses, respectively. The average solution uptake for the LiI/I2 solutions in Nafion 117 was 21±2% (additional details can be seen in Table 1). No evident decrease in the solution uptake was observed due to membrane osmotic dehydration at elevated concentrations.

After the wet mass measurements the membranes were transferred to known volumes (Vsol=20 mL) of milli-Q water (~18 MΩ cm) for 20 min to elute excess amounts of LiI/I2 according to the method described in ref. [18]. The elution was repeated to ensure complete LiI/I2 removal. The conductivity of the two solutions was measured using a conductivity probe (eDAQ, platinum plate electrodes, cell constant k = 10 cm-1). The solution concentration was afterwards calculated from the measured conductivity by means of a calibration curve (range 0.012 to 1.6 mM). The conductivity of the second solution was for all measurements equal to that of pure milli-Q water and was therefore not included in the determination of the eluted solution concentration (). was used together with s to determine the concentration of sorbed LiI/I2 inside the membrane (). The internal solution concentration (henceforward the tilde symbol refers to properties inside the membrane) was calculated as:where represents the solution density inside the membrane which was approximated to the density of pure water (). The determined s and are shown in Table 1.

Ion exchange capacity – titration method

A titration method similar to the one described by Erbil and Baysal [19] was adopted for the determination of the effective ion exchange capacity (iec) at several LiI/I2 concentrations. The membrane samples were first equilibrated in the LiI/I2 solutions for 24 h. The membranes were then eluted in 2 × 20 mL milli-Q water for 20 min to remove excess amounts of LiI/I2. Afterwards the samples were transferred to 2 × 10 mL of 0.050 M HCl for one day per time. The solutions were combined (Vtot = 20 ml) and known sample volumes (Vsample) were titrated with 0.050 M NaOH (cNaOH) (Metrohm Autotitrator (916Ti-Touch)). The original 0.050 M HCl solution was titrated as well. The difference in titrant volume at the equivalence point between the pure HCl solution (Vtitr,HCl) and the sample solution (Vtitr,mem) represents the amount of protons exchanged with the membrane. The iec was determined for the most concentrated LiI/I2 solutions (0.26, 0.51, 0.74 and 0.96 M) as:

An average iec of 0.80±0.02 was found which is slightly lower than the value reported from the manufacturer ~ 0.91 meq [20], [21] (see Table 2). From iec and s the immobile charge density per unit volume of solution inside the membrane can be determined as:where also in this case . The iec from the titration experiments were used to determine the immobile charge densities () (see Table 1) with an average value of ~3400 mol m-3 while an immobile charge density of 4270 mol m-3 () was determined using iec=0.91 meq .

Subject area	Chemical engineering
More specific subject area	Membrane technology for energy conversion
Type of data	Tables, graphs, figures
How data was acquired	Electrochemical Impedance Spectroscopy (CH Instruments, CHI660E)Conductivity probe (eDAQ, platinum plate electrodes, cell constant k = 10 cm-1).Autotitrator (Metrohm Autotitrator (916Ti-Touch))
Data format	Raw, analysed
Experimentalfactors	Nafion 117 pretreatment: boiled in 3 wt% H2O2 for 1 h. Washed in boiling in milli-Q water for 10 min. Boiled in 0.05 M sulfuric acid for 30 min. Washed several times in boiling water.After pretreatment the membrane was stored in 1 M LiCl until use.
Experimental features	Construction of a specialized flow cell on which 4-point EIS was performed in flow using LiI/I2solutions.Solution uptake measured as the mass of: (wet membrane – dry membrane)/dry membrane.Internal membrane solution concentration determined by measuring the conductivity related to the amount of excess ions sorbed and later released from the membrane upon immersion in water.Ion exchange capacity done by washing a solution soaked piece of membrane in water. Then the membrane was transferred to a hydrochloric acid solution of known concentration. End concentration, determined by titration, then yields the ion exchange capacity
Data source location	Hangoevej 2, 8200 Aarhus N, Denmark
Data accessibility	Data is displayed within this article
Related research article	This Data in Brief article is submitted as a companion paper to:Østedgaard-Munck, D. N., Catalano, J., Kristensen, M. B., & Bentien, A. (2017). Membrane-based electrokinetic energy conversion. Materials Today Energy, 5, 118–125.