Understanding the transition process of phase change and dehydration reaction of salt hydrate for thermal energy storage

Understanding the transition process of phase change and dehydration reaction of salt hydrate for thermal energy storage

Journal Pre-proofs Understanding the transition process of phase change and dehydration reaction of salt hydrate for thermal energy storage T.S. Yan, ...

2MB Sizes 0 Downloads 18 Views

Journal Pre-proofs Understanding the transition process of phase change and dehydration reaction of salt hydrate for thermal energy storage T.S. Yan, T.X. Li, J.X. Xu, J.W. Chao PII: DOI: Reference:

S1359-4311(19)31854-X https://doi.org/10.1016/j.applthermaleng.2019.114655 ATE 114655

To appear in:

Applied Thermal Engineering

Received Date: Revised Date: Accepted Date:

21 March 2019 29 September 2019 6 November 2019

Please cite this article as: T.S. Yan, T.X. Li, J.X. Xu, J.W. Chao, Understanding the transition process of phase change and dehydration reaction of salt hydrate for thermal energy storage, Applied Thermal Engineering (2019), doi: https://doi.org/10.1016/j.applthermaleng.2019.114655

This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

© 2019 Published by Elsevier Ltd.

Understanding the transition process of phase change and dehydration reaction of salt hydrate for thermal energy storage T.S. Yan, T.X. Li, J.X. Xu, J.W. Chao Institute of Refrigeration and Cryogenics, Shanghai Jiao Tong University, Shanghai 200240, China

Abstract Salt hydrate is one kind of the most promising materials used for both latent heat and thermochemical heat storage systems. Currently, there is a lack of fundamental understanding of the correlation mechanism of salt hydrate as phase change material (PCM) and thermochemical material (TCM). The dehydrating and melting processes of salt hydrates (Na2SO4·10H2O, CH3COONa·3H2O, MgSO4·7H2O, and SrBr2·6H2O) were carefully analyzed at different heating rates. The experimental results show that the melt of salt hydrate occurs during its dehydrating process at high heating rates, and thus leads to the formation of molten hydrate. Moreover, the molten hydrate can dehydrate to lower hydrate directly with a slower dehydrating rate and higher activation energy. Furthermore, the temperature-pressure equilibrium of salt hydrate is higher than its dehydrating equilibrium line before phase change, but it is closer to the dissolving line. An equilibrium equation is developed for describing the transition process from molten hydrate salt to directly dehydrating into lower hydrate, and also used to explain the equilibrium curve after phase change. These findings of mutual effects between phase change and thermochemical dehydration can provide the understanding on the transition process of salt hydrate for thermal energy storage. Key words: salt hydrate; phase change; thermochemical; thermal energy storage



Corresponding author. Tel.: +86-21-34206335; fax: +86-21-34206335. E-mail address: [email protected] (T.X. Li). 1

Nomenclatures PCM

phase change material

TCM

thermochemical material

Dimensional variables Tm

melting temperature [K]

Tp

peak temperature [K]

P0

standard atmospheric pressure[Pa]

R

gas constant [J/(mol·K)]

A

pre-exponential factor

Ea

apparent activation energy [kJ/mol]

ΔHm

melting enthalpy [kJ/mol]

△Sm

melting entropy [J/(mol·K)]

ΔfH0

standard enthalpy of formation [kJ/mol]

ΔfS0

standard entropy of formation [J/(mol·K)]

ΔHr0

reaction enthalpy [kJ/mol]

ΔSr0

reaction entropy [J/(mol·K)]

ΔHsol0

enthalpy of the solution [kJ/mol]

φ

integral enthalpy of dilution of the solution [kJ/mol]

L

Greek symbols α

conversion rate

β

heating rate [K/min]

x

saturation solubility [wt%]

γ

dissolved ratio of salt hydrate [wt%]

2

1. Introduction Energy deficiency is increasingly severe. To reduce the energy consumption and increase the energy efficiency are more and more urgent. Among the total energy consumption, the building energy consumption can occupy a large proportion. More than 40%~60% building energy consumption should end up with heat demand, so it is of great value to increase the heat utilization efficiency. However, many renewable energies, like solar energy and low grade thermal source, are transient, fluctuant and limited by region [1]. Thermal energy storage is an effective technology to solve these problems, which can be mainly divided into three types: sensible heat storage, latent heat storage and thermochemical energy storage. Sensible heat is stored through the temperature difference of storage medium, which is simple but faces problems like low thermal storage density and large thermal loss. Latent heat and thermochemical heat storage show high thermal storage density which have gained more and more attention. Salt hydrate (MpXq·nH2O) is recognized as promising thermal storage material and widely used as phase change material (PCM) [2, 3] and thermochemical material (TCM) [4,5] in recent years. Figure 1 shows the thermal energy storage density and working temperature of typical salt hydrate as PCM or TCM [6-8]. The latent heat energy storage using salt hydrate as PCM is based on melting/crystallizing process as phase change between solid state and liquid state, and the energy storage density (phase change enthalpy) usually varies between 100~300 kJ/kg [9]. Thermochemical energy storage using salt hydrate as TCM is based on breaking/reforming bonds of water molecule and salt in crystal structure [10] as de-/rehydration reaction between solid salt and vapor, and the energy density (reaction enthalpy) is usually as high as 1000~2000 kJ/kg [11].

3

2500 Thermochemical Storage

Salt hydrate using as TCM

Thermal Energy Density( kJ/kg)

Na2S9H2O/0.5H2O

2000 Na2SO410H2O

1500

KAl(SO4)212H2O LiClH2O

1000

CaCl2H2O MgSO47H2O/H2O MgCl2H2O/2H2O

SrBr2H2O/H2O

Na2S2O35H2O

Latent Storage

500

CH COONa3H2O Ba(OH)28H2O Na2SO410H2O 3 KAl(SO4)212H2O MgCl2H2O CaCl2H2O MgSO47H2O SrBr6H2O

0 0

50

Salt hydrate using as PCM

100 Temperature( ℃)

150

200

Figure 1. Thermal energy storage density and working temperature of salt hydrate as PCM or TCM (Temperature by 20mbar vapor pressure) [6-8]

Unlike pure salt or other ionic solids, the water within the hydrated salt crystal leads to the complex performances of melt and dehydration. The solid hydrated salt will melt into molten hydrated salt under phase change temperature. Specifically, the molten state is usually regarded as an aqueous solution. According to melting behaviors, salt hydrate can be classified into three categories: congruent melt in which the dehydrated salt can be completely soluble in the hydration water, incongruent melt where the salt is partially soluble in the hydration water, and semi-congruent melt when the liquid and solid are in equilibrium with different compositions for the hydrate transformation to a lower hydrated salt [6,12,13]. Incongruent melt always accompanies with the formation of lower hydrate salt, which is irreversible and leads to the loss of storage capacity. During the melting process of Mg(NO3)2·6H2O, the dihydrate will be formed and leads to instability. To solve this problem, Mg(NO3)2·6H2O and extra water were mixed to prepare nanoencapuslated salt hydrate, which can realize 100+ heat uptake/release cycles [14]. El-Sebaii conducted the melting test of MgCl2·6H2O in a sealed container under the extra water principle, more than 4

1000 melt/solidification cycles were conducted with stable phase change performance [15]. As far as hydrated salt dehydrating into lower hydrated salt and vapor, it is influenced by both temperature and vapor pressure. Many salt hydrates undergo multi-step dehydration. For example, the dehydration of MgCl2·6H2O has three intermediate hydrated salts, including tetrahydrate, dihydrate and monohydrate [16]. In addition, other adverse reactions could lead to the decomposition of salt hydrate. For instance, MgCl2·6H2O will decompose into Mg(OH)Cl and HCl if dehydrating temperature is too high [17]. MgSO4·7H2O will dehydrate to MgSO4·6H2O at 25-55 ºC, while an amorphous phase is formed until converts to crystalline MgSO4 at 300 ºC [18]. When the water diffusion is limited under MgSO4·7H2O grains, an aqueous solution was produced during dehydration [19]. According to the way water released, the intermediate phase and final product, Galwey [20] described the solid-state dehydration as six different ways of water evolution type. The crystal may undergo partially melting or completely melting during dehydration process, impervious outer layer will be formed, which limits the vapor diffusion. From above, there is plenty of reports on phase change and dehydration of salt hydrate independently, but the effect of melt on dehydration or effect of dehydration on melt is seldom concerned. In this paper, we want to study the melting process of salt as PCM and dehydrating process as TCM simultaneously, the correlation of salt hydrate melting and dehydrating is concerned. Transition processes of four salt hydrates (Na2SO4·10H2O, CH3COONa·3H2O, MgSO4·7H2O, and SrBr2·6H2O) were carefully analyzed at different heating rates by TGDSC and DSC, dehydrating and melting parameters were calculated. Additionally, we designed a vacuum reactor system to study the temperature and pressure evolution during phase change process. The latent heat was analyzed based on different transition routes. This study shows the link between phase change and dehydration of salt hydrate which

5

would contribute to deepen our understanding on transition process of salt hydrate and its application for thermal energy storage.

2. Materials and Methods 2.1 Materials and Characterizations Four typical salt hydrates, Na2SO4·10H2O (99%, J&K China Chemical Ltd.), CH3COONa·3H2O (99.5%, Energy Chemical Co.), MgSO4·7H2O (99.5%, Shanghai Macklin Biochemical Co., Ltd), SrBr2·6H2O (99%, Energy Chemical Co.), were selected to investigate transition processes of phase change and dehydration. Differential Scanning Calorimetry (DSC8000, Pekin-Elmer) was used to determine phase change enthalpy (ΔHm) and melting temperature (Tm) of salt hydrate. Sample mass with about 5~10 mg was sealed in aluminum dish. Different heating rates (0.5K/min, 1K/min, 2K/min, 5K/min, 10K/min) were tested. Simultaneous thermal analysis (Netzsch STA 449 F3) was used to determine dehydration characteristics of salt hydrate, including reaction temperature, mass loss, reaction enthalpy, etc. About 20 mg salt hydrate sample was placed inside an aluminum dish and heated from room temperature to 300 ºC at different heat rates (0.2K/min, 1K/min, 2K/min, 5K/min, 10K/min). Nitrogen was used as a purge gas with 20mL/min to remove the moisture from the crucible.

2.2 transition experimental setup for salt hydrate A vacuum experimental device system was designed to test temperature and pressure change during transition process as shown in figure 2. The system is mainly composed of a vacuum reactor, a liquid reservoir, thermostatic baths, a vacuum pump. The temperature

6

sensor, PT100 with an accuracy of 0.15 ºC and pressure sensor with an accuracy of 1 mbar were used to measure temperature and pressure evolution.

Figure 2. Experimental setup for testing the transition process of salt hydrate

Dynamic melting test, test procedures including 1) preparation step: salt hydrates with thickness about 1~2 cm and mass about 200~300g, were sealed at the vacuum reactor. Air within the reactor will degas to less than 1mbar by vacuum pump. Thermostatic bath (oil bath) is used to adjust the salt at initial temperature. 2) heating step: set the thermostatic bath at end temperature to heat the salt to melt. Typically, for Na2SO4·10H2O, temperature increases from 20ºC to 50ºC, MgSO4·7H2O from 30ºC to 70ºC, CH3COONa·3H2O from 30ºC to 80ºC, SrBr2·6H2O from 50ºC to 100ºC. Salt temperature and vapor pressure were recorded. Multi-step equilibrium test: test procedures are similar to dynamic melting test. However, the heating temperature increases step by step, salt temperature and vapor pressure will keep enough time to reach equilibrium at each steps.

7

3. Results and discussions 3.1. Dehydrating and melting process of salt hydrates To analysis the dehydrating process, salt hydrate was heated in an open environment to ensure fast vapor exchange. Dehydrating curves, including mass loss line(TG), mass loss rate line(DTG), heat flow line(DSC), are shown in figure 3, figure 4, and figures S1-S4. As for the melting process, salt hydrate was sealed in a close environment to avoid vapor evaporation. Melting process at different heating rates, shown in figures S5-S8, were conducted to compare with the dehydrating process and calculate the activation energy of dynamic transition process. From figure 3 and figure 4, it can be seen that melt can influence the dehydrating process as the heating rate increases. For sodium sulfate decahydrate, melting heat peak will appear at all tested heating rates (see figure 3a and figure S1). Only part of hydrate melt during dehydration as the melting heat ranges from 150 J/g to 200 J/g in STA, in comparison with 248J/g in DSC. The additional melt won’t influence the sequential dehydrating process seen from the comparison of mass loss rate and heat flow curves in figure S1. However, for sodium acetate trihydrate, dehydration is divided into two steps when the heating rate is larger than 1K/min. As mass loss rate (DTG) curve and heat flow curve shown in figure S2, melting heat will lead heat flow deviates from DTG curve. After molten hydrate formed, the subsequent dehydration shows lower dehydrating rate. With the heating rate increasing, the molten hydrate proportion would increase from 32% (mass loss ratio compare to total mass loss) at 2K/min to 92% at 10K/min.

8

(b)

(a) 100

i Na2SO410H2O(s)Na2SO4(s)

-9

Mass(%)

80 -10

60 40

2

ln(/Tp)

i

-50

Na2SO410H2O 0.2K/min

-13

Na2SO410H2O 10K/min

20

40

60

80

100

120

Na2SO410H2O(s)Na2SO410H2O(l)(STA) Na2SO410H2O(s)Na2SO410H2O(l)(DSC)

0.0026

140

Temp(C)

0.0028

0.0030

0.0032

1/Tp

(d) -9

i CH3COONa3H2O(s)CH3COONa(s)

100

ii CH3COONa3H2O(l)CH3COONa(s)

80

-10

60 i

20 40

ii

0

-12

-20 CH3COONa3H2O 0.2K/min

-40

CH3COONa3H2O(s)CH3COONa(s)

CH3COONa3H2O 1K/min

-60

CH3COONa3H2O 2K/min

-13

CH3COONa3H2O 5K/min

-80 -100

-11

2

Mass(%)

Na2SO410H2O(s)Na2SO4(s)

Na2SO410H2O 5K/min

Melting

(c)

Heat flow(J/g/K)

-12

Na2SO410H2O 2K/min

-150 -200

-11

Na2SO410H2O 1K/min

-100

ln(/Tp)

Heat (J/g/K)

0

CH3COONa3H2O 10K/min

20

40

60

80

100 120 140 160 180 200

0.0024

Temp(C)

CH3COONa3H2O(l)CH3COONa(s) CH3COONa3H2O(s)CH3COONa3H2O(l) (DSC)

0.0026

0.0028

0.0030

1/Tp

Figure 3. TG-DSC curves of (a) Na2SO4·10H2O and (c) CH3COONa·3H2O vary with the change of heating rate β(K/min). Kinetic results of thermal dehydration or melting steps for (b) Na2SO4·10H2O and (d) CH3COONa·3H2O

For magnesium sulfate heptahydrate shown in figure 4a and figure S3, the dehydrating process can be divided into three steps. The first step was dehydration from MgSO4·7H2O to MgSO4·6H2O by losing 1 mol water at temperature range of 30~50 °C, the second step was dehydration from MgSO4·6H2O to MgSO4·xH2O (x < 0.1) at temperature range of 50~260 °C, and the final step (260~280 ºC) completely lost all the water to become MgSO4. Results are in agreement with previous reports in literature [18]. When heating rate is larger than 1K/min, a melting peak appears at 45~65 °C which is similar to sodium sulfate decahydrate. With the heating rate increasing from 2 K/min to 10 K/min, the proportion of

9

0.0032

melt will increase from 2% to 70%. When the heating rate reaches 10 K/min, the multi-step dehydration process becomes indistinguishable.

(b)

(a)

-9 i MgSO47H2O(s)MgSO46H2O(s)

100

Mass(%)

ii MgSO46H2O(s)MgSO4(s)xH2O(s)

-10

80 60

-11

i

2

20 40

ln(/Tp)

ii

Heat(J/g/K)

0 -20

MgSO47H2O 0.2K/min

-40

MgSO47H2O 2K/min

MgSO47H2O 1K/min

-12

-13

MgSO47H2O 10K/min

50

100

150

200

250

300

Temp(C)

(c)

MgSO46H2OMgSO4(s) MgSO47H2O(s)MgSO47H2O(l)(STA)

MgSO47H2O 5K/min

-60 Melting

MgSO47H2O(s)MgSO46H2O(s)

-14 0.0026

MgSO47H2O(s)MgSO47H2O(l)(DSC)

0.0028

0.0030

0.0032

1/Tp

(d)

110 i SrBr26H2O(s)SrBr2H2O(s)

100

Mass(%)

-9.5

ii SrBr2H2O(s)SrBr2(s)

90 80

-10.0

60

2

ii

i

10 50

ln(/Tp)

70

Heat(J/g/K)

0 -10

-11.0

SrBr26H2O 1K/min

-20

SrBr26H2O(s)SrBr2H2O(s) 1~5K/min

SrBr26H2O 2K/min

-11.5

SrBr26H2O 5K/min

-30

SrBr26H2O 8K/min

-40 -50

-10.5

100

150

200

250

SrBr26H2O(s) SrBr2H2O(s) 2nd stage 8~12K/min SrBr26H2O(s) SrBr26H2O(l) (DSC)

SrBr26H2O 10K/min

50

SrBr26H2O(s) SrBr2H2O(s) 1st stage 8~12K/min

300

-12.0 0.0022

0.0024

0.0026

0.0028

1/Tp

Temp(C)

Figure 4. TG-DSC curves of (a) MgSO4·7H2O and (c) SrBr2·6H2O vary with the change of heating rate β(K/min). Kinetic results of thermal dehydration or melting steps for (b) MgSO4·7H2O and (b) SrBr2·6H2O

While for strontium bromide hexahydrate, two-step dehydration of SrBr2·6H2O occurs at low heating rate (1~5K/min). The first step is dehydrating five water from hexahydrate to monohydrate at 40~120 ºC. The second step is dehydrating from monohydrate to anhydrate at 150~250 ºC. When heating rate is higher than 8K/min, the melting effect can be observed. Seen from figure 4c and figure S4, melt leads to two mass loss peak occurring at hexahydrate dehydrating to monohydrate step. Compare mass loss 10

rate curve with heat flow curve, the first peak of heat flow is composed of dehydration and melt. The mass loss proportion of 1st peak and 2nd peak decreases with heating rate increasing, which means no stable intermediate product formed but the melt interrupts the dehydrating process. Transition parameters are summarized in table 1. The dehydrating process or melting process can be described by Arrhenius equation. The activation energy is used to evaluate the energy required to initiate reaction, which can be calculated by Kissinger equation [21, 22].  a d  A e RT f ( ) dt E

d(ln

(1)



) Tp 2 E  a 1 R d( ) Tp

(2)

Compare melting parameters in STA and DSC, for Na2SO4· 10H2O and MgSO4· 7H2O, only part of salt hydrate melt in STA, and activation energies are lower than melt in DSC. In addition, when melt happens, the activation energy of sequential dehydration is higher than the original dehydration. It is consistent with higher dehydrating temperature and lower dehydrating rate for molten hydrate directly dehydrating into lower hydrate shown in figure 3 and figure 4.

11

Table 1 Transition parameters of dehydrating and melting processes at different heating rates β(K/min). Heating rate β

Mass loss Δm

Dehydrating heat or melting heat

Peak temperature Tp

Activation energy

(K/min)

(%)

(kJ/kg)

(ºC)

(kJ/mol)

1. Na2SO4·10H2O(s)→Na2SO4(s)+10H2O(g)

0.2/1/2/5/10

53.9

1764.0

45.6/60.9/70.2/77.0/99.2

67.9

2. Na2SO4·10H2O(s)→Na2SO4·10H2O(l) by STA

0.2/1/2/5/10

-

160.3/202.8/199.0/153.7/150.9

35.9/39.8/42.7/46.0/49.9

224.3

3.Na2SO4·10H2O(s)→Na2SO4·10H2O(l) by DSC

1/2/5/10

-

248.0

34.3/35.4/37.4/38.9

386.4

1. CH3COONa·3H2O(s)→CH3COONa(s)+3H2O(g)

0.2/1/2/5/10

40.3

1334.0

44.8/67.0/71.2/81.5/80.1

75.3

2. CH3COONa·3H2O(l)→CH3COONa(s)+3H2O(g)-2nd

2/5/10

13.2/26.4/37.1

-

67.1/70.2/65.1

97.9

3. CH3COONa·3H2O(s)→CH3COONa·3H2O(l) by DSC

1/2/5/10

-

249.7

58.6/59.3/63.4/63.8

332.1

1.MgSO4·7H2O(s)→MgSO4·6H2O(s)+H2O(g)

1/2/5/10

5.9/4.5/1.5/-

229.0/182.5/133.6/99.1

45.0/47.9/49.4/51.6

303.4

2. MgSO4·7H2O(s)→MgSO4·7H2O(l) by STA

1/2/5/10

-

1.7/16.3/23.6/49.3

51.4/52.9/55.9/57.9

88.4

3.MgSO4·6H2O(s)→MgSO4·xH2O(s)+(6-x)H2O(g)

1/2/5/10

43.2/45.5/47.3/52.3

1638.3/1510.5/1323.4/1305.9

57.5/78.1/85.6/95.4/97.1

30.3

4. MgSO4·7H2O(s)→MgSO4·7H2O(l) by DSC

0.5/1/2/5/10

-

69.4

50.1/51.0/51.7/54.2/56.4

398.4

1.SrBr2·6H2O(s)→SrBr2·H2O(s)+5H2O(g)

1/2/5

25.2

791.0

80.5/91.1/108.1

59

1-1 SrBr2·6H2O(s)→SrBr2·H2O(s)+5H2O(g) 1st stage

8/10/12

20.95/18.62/15.86

123/126.1/128.1

97.5

154.4/156.9/158.8

134.1

89.4/89.1/89.2/91.3/93.6

502.2

Transition process Na2SO4·10H2O

CH3COONa·3H2O

MgSO4·7H2O

SrBr2·6H2O

667.0 1-2 SrBr2·6H2O(s)→SrBr2·H2O(s)+5H2O(g) 2nd stage

8/10/12

4.62/6.92/9.34

2. SrBr2·6H2O(s)→SrBr2·6H2O(l) by DSC

0.5/1/2/5/10

-

145.7

12

3.2 Transition route of salt hydrates From above results, when heating temperature is higher than phase change temperature, melt will happen and alter the dehydrating process. Molten salt hydrate directly dehydrating to lower hydrate could be described as follows: M p X q  nH 2 O(l )  M p X q  mH 2 O( s )  (n  m) H 2 O( g ) H r0', n  m   f H 0 ( M p X q  mH 2 O, s )   f H 0 ( M p X q  nH 2 O, l ) (n  m)   f H 0 ( H 2 O, g )

 f H 0 (M p X q  nH 2 O, l )   f H 0 (M p X q  nH 2 O, s)  H m (M p X q  nH 2O)

(3) (4)

S r0', n  m   f S 0 ( M p X q  mH 2 O, s )   f S 0 ( M p X q  nH 2 O, l ) (n  m)   f S 0 ( H 2 O, g ) (5)  f S 0 ( M p X q  nH 2 O, l )   f S 0 ( M p X q  nH 2 O, s ) 

 ln

H m ( M p X q  nH 2 O ) Tm

Gr' , mn H r0', mn Sr0', mn p    RT RT R p0

(6)

(7)

So, the dehydrating equilibrium line can be corrected based on thermal physical parameters [7,23-25] and DSC results (see in table 1) The temperature and pressure evolution during melting process was conducted to analysis its transition route. Figure 5 and figure S9 illustrate the temperature and pressure route of MgSO4·7H2O during melting process. The vapor pressure will increase with the salt temperature before melt, then a plateau (50.5ºC, 9.5kPa) was found when reaching phase change temperature. After completely melt, salt temperature and vapor pressure continually increases. It should be noted that the P-T curve of MgSO4·7H2O phase change process is above the solid-gas equilibrium line when TTm was close to a new proposed equilibrium line of MgSO4·7H2O (molten state) and MgSO4·6H2O (solid state), which can be seen in figure 5-6, while this behavior has not gained much attention by other researchers. 13

7(l)->(s)

water

Pressure(kPa)

10

O( H2 6

s) O(

H2

7

M

s)

SO 4

Mg

O4 gS

7(s)->(s)

1 20

30

40

50

Temp(



60

)

70

80

Figure 5. Melting process of MgSO4·7H2O under phase diagram, black dashed line stands for the dynamic phase change route, the black rectangle scatter is the equilibrium point gained from figure S9, 7(s)->6(s) means the solid state heptahydrate dehydrating into hexahydrate, 7(l) means the molten state of heptahydrate. 100 6(l)->1(s)

Pressure(kPa)

Water

10

Sr

6H Br 2

O(

s)

2

SrBr2H2O(s) 6(s)->1(s) 1 40

60

80

Temp(



100

120

)

Figure 6. Melting process of SrBr2·6H2O under phase diagram, black dashed line stands for the dynamic phase change route, the black rectangle scatter is the equilibrium point gained from figure S10, 6(s)->1(s) means the solid state hexahydrate dehydrating into monohydrate, 6(l) means the molten of state hexahydrate.

The dynamic melting process and multi-step equilibrium procedure of SrBr2·6H2O are 14

shown in figure S10, and the P-T line is drawn in figure 6. The salt temperature keeps stable at 89.6 ºC, while vapor pressure is 17.6 kPa, particularly, the vapor pressure starts to increase before salt temperature begins to increase. Similar to MgSO4·7H2O, the equilibrium P-T curve of SrBr2·6H2O before melt won’t follow the crystallization line, which is closer to the solubility line. After completely melt, the molten salt hydrate dehydrating into lower hydrate can explain the temperature and pressure relation. The same experimental phenomenon was also found for Na2SO4·10H2O and CH3COONa·3H2O shown in figure S11-S14.

3.3 Latent heat evaluation for salt hydrate As previously mentioned, the equilibrium state before melt of salt hydrate is above the crystalline line but close to the solubility line. Generally, salt hydrate melt is transition from crystalline phase to liquid phase which undergoes no chemical change. However, the dehydrating process accompanying with salt dissolving is suggested here to evaluate the melting process. Two-step melting equations are described as follow, M p X q  nH 2 O( s)  M p X q  mH 2 O( s )  (n  m) H 2 O(l ) M p X q  mH 2 O(s)  (n  m) H 2 O(l )  pM q  (aq)  qX p  (aq)  nH 2 O(l ) H r0, n  m   f H 0 ( M p X q  mH 2 O, s )   f H 0 ( M p X q  nH 2 O, s ) (n  m)+ f H 0 ( H 2 O, l )

(8)

0 H sol  { p f H 0 (M q  , aq)  q f H 0 ( X p  , aq)  m f H 0 ( H 2 O, l )}- f H 0 (M p X q  mH 2 O, s) (9)

H sol  H 0 sol   L

(10)

H m  H Sol  H r0, nm

(11)



M MP Xq

18x(n  m)  x(M M P X q  18m)

(12)

As the phase change point is the equilibrium point of both dehydrating and dissolving processes, so it is possibility to realize phase transition from solid phase to highly concentrated aqueous phase. According to equations (8) to (12), the phase change enthalpy 15

was calculated to evaluate the validity of the dehydrating-dissolving path as listed in Table 2. It can be seen that the calculated melting heat of Na2SO4·10H2O was only 3% compared to experimental results. It indicates the possibility of Na2SO4·10H2O undergoes dehydrating-dissolving to accomplish melting and results in an aqueous solution. In addition, due to high content and weak linked water molecule within sodium sulfate decahydrate crystal, molten hydrate state can be considered to concentrated salt aqueous, thus leading the unstable for phase change cycle. However, calculated melting heats of MgSO4·7H2O, CH3COONa·3H2O and SrBr2·6H2O differ widely, which indicates the dehydratingdissolving route is discrepant for them. Thus, the aqueous solution was not the final melting product.

Table 2. Thermal evaluation for the dehydrating-dissolving route [7,25-27] ΔHr0

ΔHsol

Dissolve

ΔHm,cal

ΔHm,exp

Error

kJ/mol

kJ/mol

ratio,γ

kJ/mol

kJ/mol

%

Na2SO4·10H2O

81.4

-8.795

0.56

76.5

79.27

3

MgSO4·7H2O

15.92

-13.19

0.13

14.18

16.45

14

CH3COONa·3H2O

37.02

-17.30

0.94

20.75

37.16

44

SrBr2·6H2O

69.24

-42.22

0.85

33.48

49.80

33

4. Conclusions Four salt hydrates (Na2SO4·10H2O, CH3COONa·3H2O, MgSO4·7H2O, and SrBr2·6H2O) were used as typical examples to investigate the relation of dehydrating process and phase change melting process. The influence of heating rate on the dehydrating of salt hydrate was carefully studied. High heating rate during dehydrating process leads to the formation of molten hydrate, and the molten hydrate shows a lower dehydrating rate and higher activation energy when it further directly dehydrates to lower hydrate. The equilibrium curve of molten hydrate salt for directly dehydrating into lower hydrate is 16

determined by its melting enthalpy and entropy. The temperature-pressure equilibrium of salt hydrate before melt is above its dehydrating equilibrium line but close to the dissolving line. Dehydrating-dissolving process is suggested to describe the phase change of Na2SO4·10H2O. The mutual effect of melt and dehydration has a crucial influence on the equilibrium state and the thermal decomposition kinetics of salt hydrate, which are essential for the design and operation of thermochemical storage system. However, complex composition and crystal structure of salt hydrate during phase change and dehydration need further research to fully understand the transition process of salt hydrate for thermal energy storage.

Acknowledgments This work was supported by the National Key R&D Program of China under the contract No.2018YFE0100300 and the National Natural Science Foundation of China under the contract No.51876117.

17

References [1]

Frazzica, A. and Freni, A. Adsorbent working pairs for solar thermal energy storage in buildings, Renew. Energy 2017, 110, 87-94.

[2]

Mohamed, S. A.; Al-Sulaiman, F. A.; Ibrahim, N. I.; Zahir, M. H.; Al-Ahmed, A.; Saidur, R.; Yılbaş, B. S.; Sahin, A. Z. A review on current status and challenges of inorganic phase change materials for thermal energy storage systems, Renew. Sustain. Energy Rev. 2017, 70, 1072-1089.

[3]

Pereira da Cunha, J.; Eames, P. Thermal energy storage for low and medium temperature applications using phase change materials – A review, Appl. Energy 2016, 177, 227-238.

[4]

Ma X. K.; Xu C.; Yu Z.; Ju X. A review of salt hydrate-based sorption technologies for long-term thermal energy storage, Chin. Sci. Bull. 2015, 60, 3569-3579.

[5]

Yu N.; Wang R. Z.; Wang L. W. Sorption thermal storage for solar energy, Prog. Energy Combust. Sci. 2013, 39, 489-514.

[6]

Sharma A.; Tyagi V. V.; Chen C. R.; Buddhi D. Review on thermal energy storage with phase change materials and applications, Renew. Sustain. Energy Rev. 2009, 13, 318-345.

[7]

Donkers P. A. J.; Sögütoglu L. C.; Huinink H. P.; Fischer H. R.; Adan O. C. G. A review of salt hydrates for seasonal heat storage in domestic applications, Appl. Energy 2017, 199, 45-68.

[8]

Glasser L. Thermodynamics of inorganic hydration and of humidity control, with an extensive database of salt hydrate pairs, J. Chem. Eng. Data 2014, 59, 526-530.

[9]

Kenisarin M.; Mahkamov K. Salt hydrates as latent heat storage materials: Thermophysical properties and costs, Sol. Energy Mater. Sol. Cells 2016, 145, 255286. 18

[10] N'Tsoukpoe K. E.; Rammelberg H. U.; Lele A. F.; Korhammer K.; Watts B. A.; Schmidt, T.; Ruck, W. K. L. A review on the use of calcium chloride in applied thermal engineering, Appl. Therm. Eng. 2015, 75, 513-531. [11] N’Tsoukpoe K. E.; Schmidt T.; Rammelberg H. U.; Watts B. A. and Ruck W. K. L. A systematic multi-step screening of numerous salt hydrates for low temperature thermochemical energy storage, Appl. Energy 2014, 124, 1-16. [12] Pielichowska K.; Pielichowski K. Phase change materials for thermal energy storage, Prog. Mater Sci. 2014, 65, 67-123. [13] Tatsidjodoung P.; Le Pierrès N.; Luo L. A review of potential materials for thermal energy storage in building applications, Renew. Sustain. Energy Rev. 2013, 18, 327349. [14] Graham M.; Shchukina E.; De Castro P. F. and Shchukin D. Nanocapsules containing salt hydrate phase change materials for thermal energy storage, J. Mater. Chem. A 2016, 4, 16906-16912 [15] El-Sebaii A. A.;Al-Heniti S.;Al-Agel F.;Al-Ghamdi A. A. and Al-Marzouki F. One thousand thermal cycles of magnesium chloride hexahydrate as a promising PCM for indoor solar cooking, Energy Convers. Manage. 2011, 52, 1771-1777 [16] Ferchaud C. J.; Zondag H. A.; J. Veldhuis B. J.; de Boer R. Study of the reversible water vapour sorption process of MgSO4·7H2O and MgCl2·6H2O under the conditions of seasonal solar heat storage, Journal of Physics: Conference Series 2012, 395, 012069 [17] Whiting G. T.; Grondin D.; Stosic D.; Bennici S.; Auroux A. Zeolite–MgCl2 composites as potential long-term heat storage materials: Influence of zeolite properties on heats of water sorption, Sol. Energy Mater. Sol. Cells, 2014, 128: 289295 19

[18] Van Essen, V. M.; Zondag, H. A.; Gores, J. C.; Bleijendaal, L. P. J.; Bakker, M.; Schuitema, R.; van Helden, W. G. J.; He, Z.; Rindt, C. C. M. Characterization of MgSO4 hydrate for thermochemical seasonal heat storage, J. Energy Eng. 2009, 131, 041014. [19] Donkers, P. A. J.; Beckert, S.; Pel, L.; Stallmach, F.; Steiger, M.; Adan, O. C. G. Water transport in MgSO4·7H2O during dehydration in view of thermal storage, J. Phys. Chem. C 2015, 119 (52): 28711-28720. [20] Galwey, A K. Structure and order in thermal dehydrations of crystalline solids, Thermochim. Acta 2000, 355, 181-238. [21] Kissinger, H. E. Reaction Kinetics in Differential Thermal Analysis, Anal. Chem. 1957, 29, 1702-1706. [22] Sun X.; Lee K. O.; Medina M. A. et al. Melting temperature and enthalpy variations of phase change materials (PCMs): a differential scanning calorimetry (DSC) analysis, Phase Transitions 2018, 91, 667-680. [23] Steiger, M.; Linnow, K.; Ehrhardt, D.; Rohde, M., Decomposition reactions of magnesium sulfate hydrates and phase equilibria in the MgSO4–H2O and Na+–Mg2+– Cl−–SO42−–H2O systems with implications for Mars, Geochim. Cosmochim. Ac. 2011, 75, 3600-3626. [24] Fopah-Lele A. and Tamba J. G. A review on the use of SrBr2·6H2O as a potential material for low temperature energy storage systems and building applications, Sol. Energy Mater. Sol. Cells 2017, 164, 175-187 [25] Steiger, M.; Asmussen, S. Crystallization of sodium sulfate phases in porous materials: The phase diagram Na2SO4–H2O and the generation of stress. Geochim. Cosmochim. Ac. 2008, 72, 4291-4306. [26] Lide D R. CRC handbook of chemistry and physics, Crc Press 2007, 257, 423. 20

[27] Trausel, F.; de Jong, A.J.; Cuypers, R. A Review on the Properties of Salt Hydrates for Thermochemical Storage, Energy Procedia 2014, 48, 447-452.

21

Supplementary material

100

60 i

Heat (J/g/K)

40 20

2

0

0

-20

-2

-40

-4

-60

-6

-80

Melting

-100

40

DTG(%/min)

Mass/%

i Na2SO410H2O(s)Na2SO4(s)

80

-8 60

80

100

Temp(C)

Figure S1 TG-DTG-DSC curves of Na2SO4·10H2O with heating rate at 2K/min. At 35~50 ºC, an obvious heat peak but restricted mass loss peak indicates melt happens.

i CH3COONa3H2O(s)CH3COONa(s)

100

60

Heat flow(J/g/K)

i

ii

20 40

1

0

0 -1

-20 -2

-40 -60

-3

DTG(%/min)

Mass(%)

ii CH3COONa3H2O(l)CH3COONa(s)

80

-4

40

60

80

100

120

140

Temp(C)

Figure S2 TG-DTG-DSC curves of CH3COONa·3H2O with heating rate at 2K/min. The dehydration process can be divided into two stage due to the melting process.

22

Mass(%)

100

i MgSO47H2O(s)MgSO46H2O(s) ii MgSO46H2O(s)MgSO4(s)

80 60 ii

i

20 40

Heat(J/g/K)

-1

-20 Melting

-40

-2

DTG(%/min)

0

0

-60 -3

40

60

80

100

120

140

Temp(C)

Figure S3 TG-DTG-DSC curves of MgSO4·7H2O with heating rate at 2K/min. At 50~60 ºC, an obvious heat peak but restricted mass loss peak indicates melt occurs 110 i SrBr26H2O(s)SrBr2H2O(s)

Mass(%)

100

ii SrBr2H2O(s)SrBr2(s)

90 80 70 i

60

ii

5 50

-2

-5 2nd peak

-10 -15

-4

1st peak Melting

50

100

DTG(%/min)

Heat(J/g/K)

0

0

-6

150

200

250

300

Temp(C)

Figure S4 TG-DTG-DSC curves of SrBr2·6H2O with heating rate at 10K/min. The mass loss curve can be divided into three peaks due to the melting process, the first heat peak is composed of melting heat and dehydrating heat.

23

Na2SO410H2O(s)Na2SO410H2O(l)

Heat flow(mJ/K)

10K/min 5K/min 2K/min

1K/min 0.5K/min

25

30

35

40

45

Temp(C)

Figure S5 DSC curves of Na2SO4·10H2O with different heating rates

Heat flow(mJ/K)

CH3COONa3H2O(s)CH3COONa3H2O(l)

10K/min 5K/min 2K/min 1K/min

45

50

55

60

65

70

Temp(C)

Figure S6 DSC curves of CH3COONa·3H2O with different heating rates

24

MgSO47H2O(s)MgSO47H2O(l)

10K/min

Heat flow(mJ/K)

5K/min

2K/min

1K/min

0.5K/min

40

45

50

55

60

65

70

Temp(C)

Figure S7 DSC curves of MgSO4·7H2O with different heating rates

SrBr26H2O(s)SrBr26H2O(l)

Heat flow(mJ/K)

10K/min

5K/min 2K/min

1K/min

0.5K/min

80

85

90

95

100

Temp(C)

Figure S8 DSC curves of SrBr2·6H2O with different heating rates

25

(b)

Tsalt(C)

70

80

25

70

20

60

20

50

15

40

10

5

30

5

0 60

20

60 15 50 10 40

30 0

20

40

Tsalt(C)

Tsalt Pvapor

30

Pvapor(kPa)

80

30 Tsalt Pvapor

25

Pvapor(kPa)

(a)

0 0

50

100

t (min)

150

200

250

300

t (min)

Figure S9 Salt temperature and vapor pressure evaluation of MgSO4·7H2O during (a) a dynamic melting process and (b) multi heating steps under equilibrium condition

(b) 30 Tsalt Pvapor

100

120

50 Tsalt Pvapor

110

25

40

100

15 10

60

90 Tsalt(C)

80

Pvapor(kPa)

Tsalt(C)

20

30

80 20

70 60

5 40 0

20

40

60

80

0 100

t (min)

10

50 40 0

100

200

300

400

t (min)

Figure S10 Salt temperature and vapor pressure evaluation of SrBr2·6H2O during (a) a dynamic melting process and (b) multi heating steps under equilibrium condition 26

500

0 600

Pvapor(kPa)

(a)

Na2SO4 (s)->Na2SO4(aq)

10

Pressure(kPa)

Water

Na2SO410H2O(s) N

O aS

4(s

)

2

10(s)->0(s) 1 10

20

30

Temp(



40

)

50

60

Figure S11 Melting process of Na2SO4·10H2O under phase diagram, black dashed line stands for the dynamic phase change route, the black rectangle scatter is the equilibrium point gained from figure S12, 10(s)->0(s) means the solid state decahydrate dehydrating into anhydrate.

(a)

(b) 20 Tsalt Pvapor

50

14

Tsalt Pvapor

50

12

15

30 5

20

Tsalt(C)

10

Pvapor(kPa)

Tsalt(C)

40 40

10 8

30

6 4

20

2 10

0 0

100

200

300

400

500

10 0

t (min)

20

40

60

t (min)

Figure S12 Salt temperature and vapor pressure evaluation of Na2SO4·10H2O during (a) a dynamic melting process and (b)multi heating steps under equilibrium condition

27

0 80

Pvapor(kPa)

60

3(l)->0(s) Water

Pressure(kPa)

10

a

CO

ON

O( H2 3

s)

a(s

)

N OO

C CH 3

CH 3 3(s)->0(s) 1 20

40

Temp(



)

60

80

Figure S13. Melting process of CH3COONa·3H2O, black dashed line stands for the dynamic phase change route, the black rectangle scatter is the equilibrium point gained from figure S14, 3(s)->0(s) means the solid state trihydrate dehydrating into anhydrate, 3(l) means the liquid state trihydrate.

80

(b)

20

70

25 80

Tsalt Pvapor

15

Tsalt Pvapor

20

70

50 40

5

Tsalt(C)

10

Pvapor(kPa)

Tsalt(C)

60 60

15

50

10

40

30

5

30 0 0

20

40

60

80

0

t (min)

100

200

300

400

500

t (min)

Figure S14 Salt temperature and vapor pressure evaluation of CH3COONa·10H2O during a dynamic melting process and multi heating steps under equilibrium condition

28

0 600

Pvapor(kPa)

(a)

Highlights

1. The correlation of phase change and dehydration reaction of salt hydrates are analyzed

2. Phase change has an influence on the dehydrating process at high heating rate.

3. The equilibrium route between temperature and pressure during melting is revealed.

4. Dehydrating-dissolving melting of salt hydrate is discussed.

29