Subscriber access provided by University of South Dakota
C: Surfaces, Interfaces, Porous Materials, and Catalysis
Efficiency of Carnot and Conventional Capacitive Deionization Cycles Daniel Moreno, and Marta C. Hatzell J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b05940 • Publication Date (Web): 07 Sep 2018 Downloaded from http://pubs.acs.org on September 7, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Efficiency of Carnot and Conventional Capacitive Deionization Cycles Daniel Moreno and Marta C. Hatzell∗ 1
George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Ga E-mail: *
[email protected] Phone: +1 404-385-4503
2
Abstract
3
Brackish water treatment characteristically suffers from low efficiency. Overcoming
4
this challenge requires a fundamental understanding regarding the energetics of both
5
the ion removal process and the trajectory of the complete deionization cycle. Evalu-
6
ating the processes and cycles in tandem can elucidate potential sources of inefficiency
7
(entropy production), and may promote new operational schemes leading to higher ef-
8
ficiency low saline separations. Here, using a Gouy-Chapman-Stern model we elucidate
9
the efficiency of two capacitive deionization cycles. In the first cycle, electroadsorption
10
occurs as the cell voltage increases, which resembles processes employed experimen-
11
tally. In the second cycle, a theoretical Carnot-like electroadsorption process occurs
12
while maintaining a constant number of ions (N=constant). We show that the pro-
13
jected thermodynamic efficiency of both cycles exhibited similar efficiency which range
14
from 2-4.5%, while operating with a low saline feed (Cfeed = 20 mM) stream. However,
15
at elevated feed concentrations (100 mM), the Carnot-like electroadsorption process
16
resulted in a significant improvement in thermodynamic efficiency (30%), when com-
17
pared with the conventional electroadsorption process (15%). While the Carnot-like
1
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
18
electroadsorption processes is not technically feasible, it may serve as a benchmark for
19
understanding thermodynamic limits associated with low saline separations achieved
20
through electrosorption.
2
ACS Paragon Plus Environment
Page 2 of 24
Page 3 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
21
The Journal of Physical Chemistry
Introduction
22
The inextricable connection between water and energy promotes investigations that re-
23
duce energy consumption during desalination processes. Moving toward the use of low saline
24
brackish water as a feed solution is one simple approach to achieve low specific energy con-
25
sumption (SEC) desalination. This is because the theoretical minimum energy required for
26
separating brackish waters (≤0.3 kWh/m3 ) is significantly less than that required for seawa-
27
ter separations (3-10 kWh/m3 ). 1,2 Brackish water separation is also increasingly favorable
28
in remote or inland regions which do not have reliable access to seawater. 3
29
A growing interest beyond attaining low SEC is achieving high desalination efficiency.
30
The efficiency by which desalination technologies are evaluated are based on thermodynamic
31
principles or exergy considerations. The thermodynamic energy efficiency (TEE) provides
32
information regarding how far a technology is from the theoretical minimum, whereas exer-
33
getic efficiencies (ExE) provide guidance on how far a technology is away from a practical
34
minimum, based on chemical exergy changes. The TEE has primarily been the guiding tool
35
to measure the efficiency of reverse osmosis (RO), with efficiencies improving from only 5%
36
in the 1970s to over 50% today. 4,5 It should be noted that achieving 100% is not practically
37
possible; therefore little gains in terms of energy efficiency can be made through solely im-
38
proving reverse osmosis for seawater desalination. This is one reason why there is greater
39
interest in mitigating energy consumption in alternative water treatment processes such as
40
pre or post treatment. 6
41
Despite the near ideal performance, RO does not maintain this efficiency when moving
42
toward low saline solutions (brackish water). 7 The SEC for brackish water RO decreases, 8–10
43
but the corresponding decrease in the Gibbs free energy of mixing is greater. This low
44
efficiency is often overlooked, as there is a greater emphasis on SEC. Furthermore, exergetic
45
efficiencies for brackish water RO are less than 10%, 11–13 with the highest reported value
46
of 16%. 14 This is why historically alternative technologies (electrodialysis and capacitive
47
deionzation) are discussed for brackish water treatment. 15–19 Most capacitive deionization 3
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
48
cycles yield efficiencies in the regime of 1-5% with feed solutions of on the order of 20
49
mM; 16,20–23 however, the aim of these investigations was not on optimizing TEE. In more
50
recent efforts which aimed to improve TEE, maximum values 9% were shown to be possible. 24
51
Here, we aim to analyze the expected thermodynamic and exergetic efficiency of two
52
different capacitive deionization cycles. In the first cycle, electroadsorption and desorption
53
occur as the cell voltage is incrementally increased. This cycle is practically feasible as one
54
can execute the cycle experimentally. The second cycle is purely theoretical, and is based on
55
prior definitions of a Carnot-like chemical cycle for capacitive mixing. 25 In the second cycle,
56
salt removal occurs while maintaining a fixed number of ions (N=constant). By examining
57
both cycles, we aim to discern if the electroadsorption processes within a given desalination
58
cycle (Carnot or conventional) can limit the maximum efficiency. Furthermore, we aim to
59
better describe thermodynamic limitations associated with dilute solution separations.
60
Theory
61
The desalination cycle employed in most capacitive deionization studies includes four
62
consecutive processes. However, just like with thermal energy systems, these cycles can
63
operate using different modes. Thus, achieving ion removal is possible by various different
64
processes. 26–28 Here, we investigate two thermodynamic deionization cycles.
65
In the first cycle, we employ four processes which can be executed in a conventional CDI
66
system. Process one consists of deionization up to a set voltage Vmax . Next in process two,
67
as the feed is switched to a brine the electrodes are continued to be charged at fixed voltage
68
(iso-V - switch). In process three, the electrodes are discharged while the concentration
69
decreases. Then in process four, the electrodes are brought back to the initial conditions
70
(Cfeed replaces Cbrine ) using a constant voltage process (iso-V - switch). Since this cycle
71
operates using conventional processes, we term this the conventional cycle (Figure 1a and
72
b). Evaluating energy consumption will take place through the use of V-σ curves (Figure
4
ACS Paragon Plus Environment
Page 4 of 24
Page 5 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
73
The Journal of Physical Chemistry
1e).
74
In the second cycle, we employ processes which are analogous to a Carnot refrigeration
75
cycle. The Carnot cycles for thermal systems represent the ideal (reversible) cycle. It is well
76
known that Carnot cycles cannot be executed in reality, but the theoretical insight provides
77
guidance on how far a real cycle is away from an ideal system. The chemical Carnot cycle
78
was previously described in prior work. 25 Briefly, process one is described by deionization
79
at constant number of ions (iso-N), process two is described by deionization at constant
80
chemical potential (iso-µ), process three is described by resalination at constant number of
81
ions (iso-N), and process four is described by another resalination step at constant chemical
82
potential (iso-µ) (Figure 1c and d). We will refer to this cycle as the Carnot cycle, and
83
energy consumption will be detailed on the µ-N diagram (Figure 1f).
84
For both cycles, we assume a quasi-equilibrium condition, which implies that the cycle
85
operates on an infinitesimally slow flow rate to minimize irreversibilities. This greatly sim-
86
plifies the model since it eliminates the need to specify time dependent parameters, such
87
as current and flow rate. Additionally, considerations regarding cell geometry are not con-
88
sidered, since all quantities here will be evaluated on a per area or per volume basis. The
89
model also does not consider Faradaic losses due to side reactions or resistive losses due to
90
components such as membranes or contacts.
91
In order to describe various iso-property (µ, N, V, C, T, σ) processes for a given capacitive
92
deionization cycle a description for the electric double layer is needed. Here we employed the
93
well documented Gouy-Chapman-Stern (GCS) model to describe the electric double layer
94
dynamics necessary for ion removal. 20 In each of the simulations the maximum voltage Vmax
95
was limited to 1.2 V. In addition, the feed (Cfeed ) remained constant while the deionized
96
concentration (Cmin =Cdeionized ) and water recovery ratio α=Vdeionized /Vfeed were varied. In
97
order to maintain an ion balance during the cycle, the resulting brine concentration Cbrine is
98
determined by:
5
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Cbrine =
Cf eed − αCmin 1−α
Page 6 of 24
(1)
99
Since the Carnot analog cycles is constrained to limits of only two concentrations, the high
100
concentration is the brine Cbrine and the low concentration is the diluate Cmin . The chemical
101
potential in the cell at any point is:
µ = RT ln(
C ) Cref
(2)
102
where T is temperature, R is the universal gas constant, and Cref represents a reference
103
concentration, taken to be 1 M. 25 The total number of ions in the cell (per electrode area)
104
N is:
N = Γ + CLe Nav
(3)
105
where Nav is Avogadro’s constant, Le represents the pore length between carbon particles
106
(pore volume/pore area), taken in this study to be 4 nm, and Γ is the excess surface charge.
107
Γ is determined as a function of both the electrode surface charge σ and additional crossover
108
charge σ * : √ Γ = σ 2 + σ ∗2 − σ ∗
σ∗ =
1 2πλB λD
(4)
(5)
109
where the Debye and Bjerrum lengths for the electrolyte, λB and λD determine the crossover
110
charge. 20 Crossover charge is defined as the charge at which attractive forces from counter
111
ions and repulsive forces from co-ions are balanced. 29 When crossover charge exceeds surface
112
charge, the increase in counter ions exceeds the decrease in concentration of co-ions, resulting
113
in a net increase in number of ions adsorbed. 30 If surface charge is smaller, then Γ is directly 6
ACS Paragon Plus Environment
Page 7 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
114
The Journal of Physical Chemistry
proportional to σ. The diffuse layer voltage difference ∆VD is given by: s
∆VD = 2VT sinh−1 σ
πkB 2Cf eed Nav
(6)
115
where kB is Boltzmanns constant. The thermal voltage VT is defined as VT = RT/F, where F
116
is Faraday’s constant. For the Stern layer voltage ∆Vst , Gauss law is used for a 1:1 electrolyte
117
(as in the case with NaCl):
∆Vst = σe/Cst
(7)
118
where Cst is the Stern capacity and e is the elementary charge. Previous studies using GCS
119
models have assumed a constant Stern capacity around 0.2 F/m2 . 31,32 Assuming a symmetric
120
cell configuration, the diffuse and Stern layer potentials can be summed and then doubled
121
to provide the total voltage ∆Vcell through the cell:
∆Vcell = 2VT (∆Vst + ∆VD )
(8)
123
For CDI cycles, thermodynamic energy efficiency (TEE) is the ratio of the Gibbs energy H of mixing ∆Gmix to the input cycle work Wcycle , 33,34 measured by V dσ. To equate units
124
with the Gibbs energy of mixing, σ is represented as a volumetric-based charge after dividing
125
by 2x the pore length Le . For any generic separations process, the Gibbs energy of mixing
126
is defined as: 35
122
∆Gmix = nRT
hC
0
α
ln
C − αC C C − αC i 0 min min 0 min − ln C0 (1 − α) α Cmin (1 − α)
ηT EE =
∆Gmix Wcycle
(9)
(10)
127
where n is the van ’ t Hoff Factor, which for 1 : 1 electrolytes assumes a value of 2. Thus,
128
taking relevant limits, mixing energy should maximize when Cmin is lowest and when α is 7
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
129
Page 8 of 24
highest.
130
In the case of Carnot analog cycle, the GCS theory imposes constraints on the allowable
131
values for α and Cmin (Figure S1). These limitations occur in the recharging stage, in order
132
to maintain the constant number of ions condition. The limit varies with the chosen feed
133
concentration, water recovery ratio, and maximum allowable voltage. As a result, recovery
134
ratios above the allowable range are not computed. Detail on how this limit was found is
135
included in Supporting Information.
136
Chemical exergy of a feed can be calculated using electrolytic solution models (which use
137
feed concentration references) or Szargut models (which use lithosphere as the reference). In
138
an ideal mixture, molar chemical exergy can be calculated as a function of concentration of
139
water as well as the salt solution:
eCh sol =
X
xi eCh i + RT
i
X
xi lnxi
(11)
i
140
where xi is the mole fraction of the species considered, and eCh is standard molar chemical
141
exergy of this species. Here, the Szargut model for exergy calculation reduces to only consider
142
the activities and molar concentrations of H2 O and NaCl: s Bch = nH2 O eH2 O + nN aCl eN aCl + RT (nH2 O ln(aH2 O ) + nN aCl ln(mN aCl γN aCl ))
(12)
143
where n represents the molar concentration of each species, mNaCl is the molality of NaCl
144
solution (in mol/kg of solvent) and activity coefficients are given as aH2 O for H2 O and γ NaCl
145
for NaCl. 36 Changes in chemical exergy will be accounted for during both the charging
146
and discharging stages. Exergy will be evaluated using the Szargut model, 37 while activity
147
coefficients will be evaluated using the Pitzer model. 38 This is discussed in further detail in
148
Supporting Information (also see Figure S2 for comparative plots of chemical exergy changes
149
at varying concentrations).
150
Exergy efficiency (here referred to as ExE) is determined as the ratio of the chemical
8
ACS Paragon Plus Environment
Page 9 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
151
exergy change during charging and discharging, to the input electrical exergy Bel , which we
152
take to be the equivalent of the input work (Wcycle ):
ηExE = ∆Bch /∆Bel = ∆Bch /Wcycle
153
Results and Discussion
154
Impact of capacitive desalination cycle
(13)
155
Differences in the energetics of a complete capacitive deionization cycle is indicated by
156
a change in the area of a V-σ or µ-N diagram. The larger the area represents a cycle with
157
greater energy consumption (Figure 1e and f). Experimentalists are able to easily evaluate V-
158
σ diagrams through extrapolating experimental voltage and current (charge) data, whereas
159
in practice µ-N diagrams cannot be measured as easily since the total number of ions in the
160
cell cannot be precisely determined. However, from the differential form of the Helmholtz
161
free energy (dF=µ dN+Vdσ), one can easily see that for a complete thermodynamic cycle
162
where a system returns to the initial state (no net change in energy dF=0), the two integrals
163
are equal in magnitude and opposite in sign (µ dN=-Vdσ). Therefore, both profiles are
164
effective for evaluating total energy consumed during a capacitive deionization cycle (Figure
165
S3).
166
By fixing the initial and final conditions (Cfeed =20 mM, Vmax =1.2 V, and Cmin =1 mM),
167
of a desalination cycle, the measured energy consumption can be normalized per the same
168
number of ions removed, independent of the cycle trajectory (thermodynamic cycle). Here,
169
the fixed conditions equate to a desalination test with 75% salt removal and 50% water
170
recovery (conditions commonly employed by experimentalists). For these conditions, the
171
energy or area of both the V-σ (Figure 2a) and µ-N (Figure 2b) graphs was slightly greater
172
for the Carnot analog than for the conventional operation. For this scenario, the energy
173
consumed was 0.96 kWh/m3 (93 kT/ion) for the Carnot cycle, whereas only 0.48 kWh/m3 9
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
174
(47 kT/ion) was consumed through the conventional operation.
175
The larger energy consumed during the complete cycle can be attributed to either an
176
increased energy consumption during charging, or to a reduction in energy recovery during
177
discharging. For the individual deionization process (not the entire cycle) the energy is
178
indicated by the area under the 1-2-3 curve and 1’ -2’ -3’ on Figure 2. This energy for ion
179
removal is 8.7 kWh/m3 for the conventional operation and 7.9 kWh/m3 for the Carnot analog.
180
Thus, the Carnot analog did consume less energy during the deionization process than when
181
the cycle had operated conventionally, indicating that adsorption is indeed more ideal, albeit
182
small. Thus the reduction in cycle energy consumption is due to the conventional operation
183
recovering a larger portion of that energy. The corresponding thermodynamic efficiency
184
for the conventional operation and Carnot analog was TEE=3.6% and TEE=1.8%. The
185
thermodynamic efficiency appears low, but this is in line with other reported values. 16,21–23 As
186
the cycles studied here operate without transient or geometric considerations, we are limited
187
in the full range of explorable design parameters. Employing these additional capabilities
188
while tending to optimal operating conditions can allow for an increase in efficiencies. 24
189
Impact of feed water concentration
190
To evaluate the effect of the various boundary conditions, the TEE was measured for both
191
the Carnot analog and the conventional operation as the water recovery (α) and percent of
192
salt removed (∆C/Cfeed ) was varied (Figure 3). In all tests, the starting concentration C0
193
remained fixed at either 20 mM or 100 mM, maximum cell voltage Vmax was constrained
194
to 1.2 V, and minimum cell voltage Vmin was 0.1 V. With the 20 mM feed, increasing the
195
percent of salt removed, and water recovery resulted in a greater TEE irregardless of cycle
196
(Figure 3a and b). The TEE again was on average two times greater for the conventional
197
operation than for the Carnot cycle, but never exceeded 4%. While the percent of salt
198
removal had a significant effect on the apparent maximum thermodynamic efficiency, the
199
water recovery did not. Increasing the water recovery (α) from 25% to 75% only increased 10
ACS Paragon Plus Environment
Page 10 of 24
Page 11 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
200
the TEE by 0.5% for the conventional operation, and 2% for the Carnot cycle. The increase
201
observed at higher salt removal is due to a greater Gibbs energy of mixing needed for the
202
cycle, which grows faster than the required cycle work.
203
When the feed increased from 20 mM to 100 mM, the TEE increased by nearly four times.
204
With the Carnot cycle, the constant number of ions criteria could not be maintained at high
205
percent salt removal or higher concentrations. To overcome this model limitation, additional
206
boundary conditions with respect to the minimum concentration and water recovery were
207
placed on each simulation, taking into account the input voltage (discussed in greater detail
208
in Supporting Information; see Figure S1 and S4). It should also be noted that most CDI
209
systems do not operate with 100 mM salt solutions, as the limited available surface area
210
would likely result in salt removal of only 1% to 50% depending on the size of the system.
211
However, the aim in investigating these conditions was to evaluate the thermodynamics of the
212
separations process, as there is a greater amount of thermodynamic separations data available
213
at higher concentrations. Furthermore, with the development of CDI with unlimited surface
214
area (e.g. flow electrodes), capacitive deionization technologies are beginning to evaluate
215
higher saline solutions. 39,40
216
With 100 mM feed solutions, the TEE increases substantially from only 5% at 20 mM
217
to 20% at 100 mM (Figure 3c and d). In addition to the increase in efficiency, the Carnot
218
analog achieves higher TEE at corresponding salt removal rates and water recovery. This is
219
because for this concentration range, the Carnot cycle operates at nearly optimal conditions.
220
The higher conditions allow for the minimum cell charge to approach zero during discharge
221
(with the constant number of ions criteria), which substantially reduces the total cycle work
222
consumed. At higher concentration ranges (≥ 100 mM), the theoretical efficiency even begins
223
to approach 100% (Figure 4a). The conventional operating cycle, on the other hand, achieves
224
maximum efficiencies around 5-10% (Figure 4b).
225
The effect of feed water concentration (Figure S5), shows that there is a sharp increase
226
in the TEE for the Carnot analog if high water recovery is possible, with TEE approaching
11
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
227
30%. Unlike the Carnot analog (Figure S5a), the conventional operation does not show a
228
dependence on water recovery (Figure S5b). The stronger dependence for the Carnot analog
229
is an indicator that the cycle is approaching its near-ideal limiting cases at higher concentra-
230
tion. Thus, moving to even higher salinity feed streams would increase the theoretical limit
231
to above 50% which is in line with observed TEE found in most state of the art membrane
232
based technologies. At higher concentration, it also becomes apparent that the Carnot ana-
233
log and conventional electrosorption process differ. With the same concentration limits, the
234
Carnot analog exceeds the traditional cycle’s efficiency of ≈15% by over 2x.
235
For the cycle limits established, the effect on charge over the various concentrations was
236
also examined. Charge for the surface, crossover, and total excess charge are nearly the same
237
magnitude in all cases (Figure S6). As a result, efficiencies are not affected by charge, but
238
rather by deviations within these values during cycling which change the work required to
239
operate (See Figure 1b).
240
Exergetic efficiency of capacitive deionization cycles
241
Similar to TEE, exergetic efficiency (ExE) for the conventional cycle exceeded those of
242
the Carnot analog using the baseline conditions (Cfeed =20 mM, Cmin =1 mM, Vcell =1.2 V,
243
α=50%). ExE also generally increased with salt removal at high and low water recoveries,
244
for both the conventional cycle and the Carnot analog (Figure S7). For this nominal test
245
case, the maximum value for ExE observed was 12%, at a nearly 80% salt removal for the
246
conventional cycle (Figure S7b). This is comparable to the maximum values observed in
247
prior studies on MCDI. 14,36 As this study was computational and assumed to run infinitely
248
slow, additional exergy losses that did not need to be considered for this study (i.e. pumps)
249
would need to be accounted for in experiments, and would lower exergetic performance.
250
However in both cycles, the ExE (Figure 4c and d) efficiency increased with an increase
251
in water recovery and generally with increased salt removed. Furthermore, theoretical ExE
252
values approached 100% for the Carnot cycle and 30% for the conventional operation mode. 12
ACS Paragon Plus Environment
Page 12 of 24
Page 13 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
253
The Journal of Physical Chemistry
Conclusions
254
Various system properties and processes contribute to the observed thermodynamic and
255
exergetic efficiency of a capacitive deionization cycle. Here, we evaluate two deionization
256
cycles. One mimics an executable approach whereby ion removal occurs during charging
257
with constant voltage switching, the other is a Carnot analog cycle where ion removal occurs
258
while maintaining a constant number of ions. Under low saline conditions the energetics
259
associated with voltage regulated electroadsorption were slightly more favorable, and resulted
260
in predicted thermodynamic efficiencies that approach 2-5%. However, when feed solutions
261
were increased to 100 mM the energetics of constant N electroadsorption was favored over
262
the conventional approach, and the maximum thermodynamic efficiencies approached 37%
263
(Figure 5). In all simulations, thermodynamic efficiencies improved as Cfeed , and ∆C (Cfeed -
264
Cmin ) increased. A survey of experimental capacitive deionization system performance will
265
reveal that actual desalination performance data is far below these values (≤ 1 %) for CDI
266
and (≤ 3%) for membrane CDI. However, this is due to the low effective surface area available
267
in lab scale systems, which limit the operation to low Cfeed , ∆C, and only 50% water recovery.
268
As experimentalists begin to further explore the CDI operational space, increased TEE may
269
be observed. Interestingly, this is similar to trends observed with thermal heat engines,
270
whereby increasing the T0 and maximizing ∆T results in higher thermodynamic efficiencies,
271
due to an increase in the theoretical Carnot limit. When evaluating the CDI TEE to a
272
Carnot analog, at low feed concentrations, the actual cycle perform similarly to that of a
273
Carnot analog. This indicates that while efficiencies are low, limitations may exist which
274
prevent substantial improvements. However, efficiency improvements are increased if the
275
feed stream concentration increases.
13
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
276
Acknowledgement
277
This material is based upon work supported by the National Science Foundation under
278
Grant No. (1821843) for MCH. This work is also supported by the ARCS graduate fellowship
279
to Daniel Moreno.
280
Supporting Information Available
281
The Supporting Information includes additional detail on constraints on the Carnot ana-
282
log cycle and calculation of exergy values. We also show additional modeling result detailing
283
the effect of water recovery, surface charge, and feed concentrations.
284
285
286
287
288
289
290
291
292
293
294
295
This material is available free of charge via the Internet at http://pubs.acs.org/.
References (1) Spiegler, K. S.; El-Sayed, Y. M. The energetics of desalination processes. Desalination 2001, 134, 109–128. (2) Oren, Y. Capacitive deionization (CDI) for desalination and water treatmentpast, present and future (a review). Desalination 2008, 228, 10–29. (3) Karagiannis, I. C.; Soldatos, P. G. Water desalination cost literature: review and assessment. Desalination 2008, 223, 448–456. (4) Kremen, S. S. Reverse osmosis makes high quality water now. Environ. Sci. Technol. 1975, 9, 314–318. (5) McGovern, R. K.; Lienhard, J. H. On the potential of forward osmosis to energetically outperform reverse osmosis desalination. J. Membr. Sci. 2014, 469, 245–250.
14
ACS Paragon Plus Environment
Page 14 of 24
Page 15 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
296
297
The Journal of Physical Chemistry
(6) Elimelech, M.; Phillip, W. A. The future of seawater desalination: energy, technology, and the environment. Science 2011, 333, 712–717.
298
(7) Walha, K.; Amar, R. B.; Firdaous, L.; Qu´em´eneur, F.; Jaouen, P. Brackish groundwater
299
treatment by nanofiltration, reverse osmosis and electrodialysis in Tunisia: performance
300
and cost comparison. Desalination 2007, 207, 95–106.
301
302
(8) Hrayshat, E. S. Brackish water desalination by a stand alone reverse osmosis desalination unit powered by photovoltaic solar energy. Renew. Energy 2008, 33, 1784–1790.
303
(9) Wilf, M.; Konstantin, M.; Chencinsky, A. Evaluation of an ion exchange system regen-
304
erated with seawater for the increase of product recovery of reverse osmosis brackish
305
water plant. Desalination 1980, 34, 189–197.
306
(10) Alghoul, M.; Poovanaesvaran, P.; Sopian, K.; Sulaiman, M. Review of brackish water
307
reverse osmosis (BWRO) system designs. Renew. Sust. Energ. Rev. 2009, 13, 2661–
308
2667.
309
(11) Cerci, Y. Improving the Thermodynamic and Economic Efficiencies of Desalination
310
Plants: Minimum Work Required for Desalination and Case Studies of Four Work-
311
ing Plants; Desalination and Water Purification Research and Development Program
312
Report; U.S. Department of the Interior, Bureau of Reclamation: Denver, USA, 2003.
313
(12) Cerci, Y. Exergy analysis of a reverse osmosis desalination plant in California. Desali-
314
315
316
317
318
nation 2002, 142, 257–266. (13) Aljundi, I. H. Second-law analysis of a reverse osmosis plant in Jordan. Desalination 2009, 239, 207–215. (14) Qureshi, B. A.; Zubair, S. M. Exergetic analysis of a brackish water reverse osmosis desalination unit with various energy recovery systems. Energy 2015, 93, 256–265.
15
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
319
320
(15) Lee, H.-J.; Sarfert, F.; Strathmann, H.; Moon, S.-H. Designing of an electrodialysis desalination plant. Desalination 2002, 142, 267–286.
321
(16) Zhao, R.; Porada, S.; Biesheuvel, P. M.; Van der Wal, A. Energy consumption in
322
membrane capacitive deionization for different water recoveries and flow rates, and
323
comparison with reverse osmosis. Desalination 2013, 330, 335 – 341.
324
(17) Suss, M. E.; Porada, S.; Sun, X.; Biesheuvel, P. M.; Yoon, J.; Presser, V. Water desali-
325
nation via capacitive deionization: what is it and what can we expect from it? Energy
326
Environ. Sci. 2015, 8, 2296–2319.
327
(18) Anderson, M. A.; Cudero, A. L.; Palma, J. Capacitive deionization as an electrochemical
328
means of saving energy and delivering clean water. Comparison to present desalination
329
practices: Will it compete? Electrochim. Acta 2010, 55, 3845–3856.
330
(19) Biesheuvel, P. M.; Bazant, M. Z.; Cusick, R. D.; Hatton, T. A.; Hatzell, K. B.;
331
Hatzell, M. C.; Liang, S., P.; Lin; Porada, S.; Santiago, J. G. et al. Ca-
332
pacitive Deionization–defining a class of desalination technologies. arXiv preprint
333
arXiv:1709.05925 2017.
334
335
(20) Biesheuvel, P. M. Thermodynamic cycle analysis for capacitive deionization. J. Colloid Interface Sci. 2009, 332, 258–264.
336
(21) Zhao, R.; Biesheuvel, P. M.; der Wal A, V. Energy consumption and constant current
337
operation in membrane capacitive deionization. Energy Environ. Sci. 2012, 5, 9520 –
338
9527.
339
(22) Dykstra, J. E.; Zhao, R.; Biesheuvel, P. M.; der Wal A, V. Resistance Identification
340
and rational process design in Capacitive Deionization. Water Res. 2018, 88, 358–370.
341
(23) Shang, X.; Cusick, R. D.; Smith, K. C. A Combined Modeling and Experimental Study
16
ACS Paragon Plus Environment
Page 16 of 24
Page 17 of 24 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
342
Assessing the Impact of Fluid Pulsation on Charge and Energy Efficiency in Capacitive
343
Deionization. J. Electrochem. Soc. 2016, 164, ES53–547.
344
(24) Hemmatifar, A.; Ramachandran, A.; Liu, K.; Oyarzun, D. I.; Bazant, M. Z.; Santi-
345
ago, J. G. Thermodynamics of Ion Separation by Electrosorption. Environ. Sci. Tech-
346
nol. 2018, In Press.
347
348
349
350
351
352
(25) Boon, N.; van Roij, R. ’Blue energy’ from ion adsorption and electrode charging in sea and river water. Mol. Phys. 2011, 109, 1229–1241. (26) Janssen, M.; H¨artel, A.; Van Roij, R. Boosting capacitive blue-energy and desalination devices with waste heat. Physical review letters 2014, 113, 268501. (27) Hatzell, M. C.; Hatzell, K. B. Blue Refrigeration: Capacitive De-ionization for Brackish Water Treatment. J. Electrochem. En. Conv. Stor. 2018, 15, 011009.
353
(28) Zhang, J.; Hatzell, K. B.; Hatzell, M. C. A Combined Heat-and Power-Driven Mem-
354
brane Capacitive Deionization System. Environ. Sci. Technol. Lett. 2017, 4, 470–474.
355
(29) Markovich, T.; Andelman, D.; Podgornik, R. Charged Membranes: Poisson-Boltzmann
356
theory, DLVO paradigm and beyond. arXiv preprint arXiv:1603.09451 2016.
357
(30) Rica, R. A.; Ziano, R.; Salerno, D.; Mantegazza, F.; van Roij, R.; Brogioli, D. Capacitive
358
mixing for harvesting the free energy of solutions at different concentrations. Entropy
359
2013, 15, 1388–1407.
360
361
362
363
(31) Biesheuvel, P. M.; van der Wal, A. Membrane capacitive deionization. J. Membr. Sci. 2010, 346, 256–262. (32) Biesheuvel, P. M.; van der Wal, A. Dynamic adsorption/desorption process model for capacitive deionization. J. Phys. Chem. C. 2009, 113, 5636 – 5640.
17
ACS Paragon Plus Environment
The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
364
(33) Wang, L.; Biesheuvel, P. M.; Lin, S. Reversible thermodynamic cycle analysis for ca-
365
pacitive deionization with modified Donnan model. J. Colloid Interface Sci. 2018, 52,
366
522–528.
367
(34) Wang, L.; Lin, S. Membrane Capacitive Deionization with Constant Current vs. Con-
368
stant Voltage Charging: Which is Better? Environ. Sci. Technol. 2018, 13, 4051–4060.
369
(35) Sandler, S. I. Chemical, Biochemical, and Engineering Thermodynamics; John Wiley
370
371
372
373
374
375
376
& Son, 2006. (36) Fritz, P. A.; Zisopoulos, F.; Verheggen, S.; Schro¨en, K.; Boom, R. Exergy analysis of membrane capacitive deionization (MCDI). Desalination 2018, 444, 162–168. (37) Szargut, J.; Morris, D. R.; Steward, F. R. Exergy analysis of thermal, chemical, and metallurgical processes; Hemisphere Publishing, New York, NY, 1987. (38) Pitzer, K. S. Thermodynamics of electrolytes. I. Theoretical basis and general equations. J. Phys. Chem. 1973, 77, 268–277.
377
(39) Hatzell, K. B.; Eller, J.; Morelly, S. L.; Tang, M. H.; Alvarez, N. J.; Gogotsi, Y.
378
Direct observation of active material interactions in flowable electrodes using X-ray
379
tomography. Faraday Discuss. 2017, 199, 511–524.
380
(40) Moreno, D.; Hatzell, M. The influence of feed-electrode concentration differences in
381
flow-electrode systems for capacitive deionization. Industrial & Engineering Chemistry
382
Research 2018, 57, 8802–8809.
18
ACS Paragon Plus Environment
Page 18 of 24
Conventional (b)
Win
++++++++++++++++++
Cfeed
_
+ _
+ _
_
_ +
_
_ _
Cdeionize
+ + + + + + + - - -- -- -- -- -- -- -- -- -- -- -- -- - -
I
+ _
_
+ _
N
+ _
+
_
__ _
_
_
Carnot-analog (d) _
_ +
_
_ _
+ + + + + + + - - - - - - - - - - - - - - - -
Process 1’-2’
Cdeionize
µ
Constant V
Cdeionize
+ + + + + -- -- -- -- -- -- -- -- -- -- -- -- -- --
V
++++++++++++++++++
Cfeed
_
__ ++
Process
Win
(c)
Conventional
++++++++++++++++++
Cfeed
Process
(e)
Win
Win
++++++++++++++++++
Cfeed
_
+ _
__ ++
_
+
__ _
_
_
Cdeionize
+ + + + + -- -- -- -- -- -- -- -- -- -- -- -- -- --
1
Process 2’-3’
(f) 2 3
Carnot-analog
2’
4
Charge
4’
1’
Chemical Potential
(a)
Voltage
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
Constant N
Page 19 of 24
Constant µ
3’
Number of Ions
Figure 1: (a) Ion removal in the IV cycle occurs through a constant current (iso-I) process followed by a (b) constant voltage (iso-V) process. (c) Ion removal in the Nµ (Carnot-like) cycle occurs through a constant number of ions (iso-N) process followed by (d) a constant chemical potential (iso-µ) process. (e) The trajectory of the cycles is evaluated on a voltage versus charge (σ), or (f) chemical potential (µ) versus number of ions (N) diagram.
19
ACS Paragon Plus Environment
The Journal of Physical Chemistry
(a)1.4
Carnot cycle Conventional Cycle
Voltage (V)
1.2
3’ 2 3
1.0
4’
0.8 0.6 0.4 0.2
2’ 1’ 14
0.0 0.0
0.1 0.2 0.3 0.4 0.5 Surface Charge (#/nm2)
(b)-2.5 Chemical Potential
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 20 of 24
-3.0 -3.5 -4.0 1 -4.5 -5.0 -5.5 -6.0 -6.5 -7.0 -7.5 0.0
4 1’
0.6
4’ 3
2’ 0.1
3’
2
0.2 0.3 0.4 0.5 Number of ions (#/nm2)
0.6
Figure 2: (a and b) N-µ and I-V capacitive deionization desalination cycles represented by a V-σ, and µ-N diagrams. Numbers represent each process within the desalination cycle. For both cycles the Cfeed used is 20 mM, with a Cmin of 1 mM and α = 0.5.
20
ACS Paragon Plus Environment
Page 21 of 24
(b )
(a ) 5
5
α= 0 .2 5 α= 0 .5 α= 0 .7 5 4
4 3
T E E (% )
T E E (% )
3 2
2
1
1 0
0 .0
0 .2
0 .4 ∆C / C
(c )
0 .6
0 .8
0
1 .0
fe e d
(d )
2 0
0 .0
0 .2
0 .4
0 .6
∆C / C
2 0
0 .8
1 .0
0 .8
1 .0
fe e d
1 5 T E E (% )
1 5 T E E (% )
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
1 0
1 0
5 5
0
0 .0
0 .2
0 .4 ∆C / C
0 .6
0 .8
0
1 .0
0 .0
0 .2
0 .4 ∆C / C
fe e d
0 .6 fe e d
Figure 3: Evaluation of thermodynamic energy efficiency with Cfeed = 20 mM feed (a) Carnot analog, and (b) conventional mode. Efficiency measured with Cfeed = 100 mM feed (c) Carnot analog and (d) conventional mode cycle.
21
ACS Paragon Plus Environment
The Journal of Physical Chemistry
(b )
1 .0
0 .8
0 .8
0 .6
( b ) 0 .6 fe e d
1 .0
∆C / C
∆C / C
fe e d
(a )
0 .4
1 .0
0 .4
0 .8 0 .2 fe e d
0 .2
0 .0
0 .2
0 .4
0 .6
0 .8
1 .0
0 .4 0 .0
1 .0
α
(c )
0 .6 0 .0
∆C / C
0 .0
(d )
0 .2
0 .4
0 .6
0 .8
α
1 .0
0 .2 0 .8
0 .8
fe e d
0 .6
∆C / C
fe e d
0 .0
∆C / C
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 22 of 24
0 .4 0 .2
0 .6
0 .0
0 .2
0 .4
0 .6
0 .8
T E E /E x E 1 6 4 2 1 1 7 4 3 1 .0 2 1 0 0 0 0 0 0 1 .0
(% ) 0 0 5 2 7 8 2 .5 .9 .2 .1 .3 .8 7 .5 6 .3 7 .2 4 .1 5 .1 0
α
0 .4 0 .2
0 .0
0 .0 0 .0
0 .2
0 .4
0 .6
0 .8
1 .0
0 .0
0 .2
α
0 .4
0 .6
0 .8
1 .0
α
Figure 4: Contour plots of thermodynamic energy efficiency (TEE) and exergetic efficiency (ExE) as a function of water recovery ratio α and salt removal ratio ∆C/Cfeed . In all cases shown here Cfeed = 20 mM. (a) TEE, Carnot analog cycle, (b) TEE, conventional mode cycle, (c) ExE, Carnot analog cycle, (d) ExE, conventional mode cycle.
22
ACS Paragon Plus Environment
Page 23 of 24
50 Carnot Analog Conventional Operation
40
TEE (%)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
The Journal of Physical Chemistry
30 20 10 0
0
20
40
60
80
100
Cfeed (mM)
Figure 5: Evaluation of thermodynamic energy efficiency (TEE) at varying initial feed concentrations for a diluate value of Cd = 5 mM and water recovery α = 0.5.
23
ACS Paragon Plus Environment
The Journal of Physical Chemistry
Graphical TOC Entry Efficiency of Desalination Cycles
High Salt
++++++++++++++++++
_ _ _ __ _ __ _ + ++ + + + + + -- -- -- -- -- -- -- -- -- -- -- -- -- --
Low Salt
Voltage
384
Carnot-analog 1’ 4’
Conventional 3 2 Chemical Potential
383
+ _
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 24 of 24
Win
2’
1 4 Charge
24
3’
Number of Ions
ACS Paragon Plus Environment