Effect and mechanism of doped graphene nanosheets on phase transition properties of sodium nitrate

Lü Hao-Xiang; Feng Dai-Li; Feng Yan-Hui; Zhang Xin-Xin

doi:10.7498/aps.71.20220354

Abstract

Authors and contacts

FullText HTML

1. 引　言
2. 实验制备及测试方法
2.1 样品制备
2.2 表征及蓄热性能测量方法
3. GNS@NaNO3复合相变材料表征及蓄热性能测量结果
3.1 GNS@NaNO3复合相变材料的形貌
3.2 GNS@NaNO3复合相变材料的蓄热性能
4. 石墨烯纳米片掺杂对硝酸钠熔盐相变特性的影响机制
4.1 模型构建及方法
4.2 复合相变材料的均方位移
4.3 石墨烯纳米片附近硝酸钠熔盐的分布
5. 结　论

Supplements

References

Cited By

Abstract

Molten salt is regarded as one of most promising candidates for solar energy storage due to possessing stable properties and large energy storage densities. However, the intrinsically low thermal conductivity of molten salt has become a bottleneck for rapid heat storage and transport. The addition of nanoparticles is generally considered to be a most effective way to improve the thermal conductivity of molten salt phase change materials (PCMs), while the phase change enthalpies of the nanocomposite phase change materials usually show two opposite trends of enhancement or decrement. Furthermore, the reason for the abnormal change of phase change enthalpy has not been clear in the literature so far, so the mechanism of change needs to be further explored. In this work, graphene nanosheets (GNS)@NaNO₃ (sodium nitrate) nanocomposite phase change materials are prepared by the hydration ultrasonic method. The materials are characterized by scanning electron microscope, and the phase change characteristics are measured using differential scanning calorimeter. Molecular dynamics simulation is carried out to explain the mechanism for the formation of the NaNO₃ dense layer and the non-collateral decrease of the enthalpy from the microscopic level. With the increase of GNS mass fraction, the melting point of the GNS@NaNO₃ composite phase change material decreases slightly while the phase change enthalpy decreases significantly with a non-colligative trend. A 13.81% decrease of the theoretical phase change enthalpy is observed with a GNS doping ratio of 1.5%. The NaNO₃ clusters observed on the surface of GNS are considered to have not melted, thereby resulting in a reduction in the phase change enthalpy. The mechanism is further investigated by molecular dynamics simulation, showing that the strong van der Waals attraction between GNS and NaNO₃ leads the 2–4 Å-thick NaNO₃ dense layer to form in the vicinity of GNS. With the increase of GNS mass fraction, the centroid equivalent distance between the dense layer and GNS gradually increases, which leads their mutual attraction to first increase and then weaken. When GNS mass fraction is 1.5%, the centroid equivalent distance reaches the position closest to the potential well, leading to a strongest mutual attraction. In other words, the phase change enthalpy decreases most obviously at this mass fraction. Thus, some conclusions can be drawn as follows. The type of interaction between molten salt and nano-enhancers and the position of the potential well are the fundamental reasons for the thickness of molten salt dense layer and the reduction of phase change enthalpy. The calculation of the interaction energy can be used to guide the selection of the mass fraction of the nano-enhancers, so as to avoid the loss of core material cluster and phase change enthalpy caused by the introduction of the nano-enhancers to a greatest extent. The preparation cost of the composite phase change material can also be reduced to a certain extent.

Keywords:

PACS:

88.30.mj	(Composite materials)
61.46.-w	(Structure of nanoscale materials)
65.80.-g	(Thermal properties of small particles, nanocrystals, nanotubes, and other related systems)
31.15.xv	(Molecular dynamics and other numerical methods)

Authors and contacts

Beijing Key Laboratory of Energy Saving and Emission Reduction for Metallurgical Industry, School of Energy and Environmental Engineering, University of Science and Technology Beijing, Beijing 100083, China

Corresponding author: Feng Dai-Li, dlfeng@ustb.edu.cn ; Feng Yan-Hui, yhfeng@me.ustb.edu.cn ;

Funds: Project supported by the National Key R & D Program of China (Grant No. 2018YFA0702300) and the National Natural Science Foundation of China (Grant Nos. 52176054, 51876007).

FullText HTML

1. 引　言

熔盐制备成本低, 具有良好的热性能和较低的蒸汽压, 作为传热流体和热储能介质(TES)被广泛应用于聚光太阳能(CSP)装置中^[1]. 其中熔融硝酸盐由于具有较好的热稳定性和相对更低的成本是最常用的盐类^[2]. 为了提高TES的热容量, 已经提出了一系列方法来改善热特性^[3], 如降低熔点和提升热导率^[4-6]. 自从Shin和Banerjee通过掺杂二氧化硅纳米颗粒使得四元氯化物共晶盐的比热增加15%, 二元碳酸共晶盐的比热增加24%以来^[7,8], 添加纳米颗粒改善熔盐相变热特性的方法就得到了广泛的推广和尝试^[2,9-13]. 相变焓作为熔盐相变热特性的重要指标, 在引入纳米颗粒掺杂后的熔盐体系中却表现出增强或减弱的两种相反趋势, 且目前文献中对于熔盐相变焓异常变化的原因尚不明确^[1,11,14,15]. 最近的研究指出, 纳米颗粒掺杂后熔盐中存在致密的压缩液体层, 压缩液体层内熔盐的性质与层外熔盐存在一定的差异^[2,10,13]. 然而, 目前的研究对于压缩液体层的形成机理及其对熔盐基复合相变材料相变热特性的影响规律并未进行深入探讨.

因此, 本文针对掺杂纳米颗粒的熔盐基复合相变材料相变热特性的变化展开研究, 并利用分子动力学模拟方法进一步讨论压缩层的形成机理及其对复合相变材料相变特性的影响规律. 从分子间相互作用类型及势阱深度的角度阐释了压缩层的形成机理以及相变热特性的非依数性变化, 旨在通过本文的工作对于未来纳米颗粒的种类及掺杂比的选择提供一定的指导. 为了更直接地探明纳米颗粒与熔盐之间的相互作用, 排除不同盐类之间相互作用的影响, 本文选取太阳盐组分中具有较高相变焓(172 J·g^–1)^[16]的硝酸钠作为单一的相变芯材. 具有超高热导率(约5000 W·m^–1·K^–1)和高比表面积(约2630 m²·g^–1)^[17]的石墨烯纳米片兼备优异的力学性能以及优良的化学和热学稳定性^[18], 被认为是最有前途的改善相变材料热性能的纳米颗粒. 考虑到单层石墨烯的引入在降低相变材料过冷度方面的表现优于多层石墨烯^[19], 本文选取单层石墨烯纳米片作为纳米增强剂开展工作.

2. 实验制备及测试方法

2.1 样品制备

以纯度为98.5%的硝酸钠作为复合相变材料的基材, 厚度为6—8 nm、平均直径为15 μm的片层状石墨烯纳米片作为纳米增强剂制备了纳米片质量分数为0%, 0.5%, 1%, 1.5%, 2%的石墨烯纳米片@硝酸钠(GNS@NaNO₃)复合相变材料(Pure NaNO₃, 0.5%GNS@NaNO₃, 1%GNS@NaNO₃, 1.5%GNS@NaNO₃, 2%GNS@NaNO₃). 制备流程参照Jo和Banerjee^[9]的方法, 并进行了一些调整, 具体如图1所示. 首先在蒸馏水中加入石墨烯纳米片与阿拉伯树胶, 而后将悬浮液置于电磁搅拌器中搅拌10 min以保证石墨烯纳米片的均匀分散(650 r·min^–1, 常温). 完成搅拌后, 将悬浮液置于超声发生器中超声处理(40 Hz) 90 min, 以使石墨烯均匀分散, 并降解阿拉伯树胶. 随后将硝酸钠加入经超声处理过的石墨烯纳米片悬浮液中, 并将混合后的硝酸钠-石墨烯纳米流体置入超声处理器中超声处理180 min, 以使硝酸钠与石墨烯纳米片充分混合形成复合纳米流体. 完成超声处理后, 将溶液置于培养皿中在热板上快速蒸干, 以减少完全蒸发的持续时间, 防止石墨烯纳米片的凝聚. 最后在培养皿中残余少量水分时将培养皿放入鼓风干燥箱中加速干燥, 以保证得到石墨烯纳米片均匀分散的复合相变材料.

$图 1 石墨烯纳米片@硝酸钠复合相变材料的制备流程\r\nFig. 1. Preparation process of GNS@NaNO3.$

图 1 石墨烯纳米片@硝酸钠复合相变材料的制备流程

Fig. 1. Preparation process of GNS@NaNO₃.

下载: 全尺寸图片幻灯片

2.2 表征及蓄热性能测量方法

采用FEI (型号Quanta 450 FEG)扫描电子显微镜(SEM)对GNS@NaNO₃复合相变材料的表面形貌进行研究. 为更好地观察石墨烯纳米片在熔盐中的分散情况, 选取×5 k的放大倍数.

GNS@NaNO₃复合相变材料的熔融和蓄热行为采用带冷却模块的梅特勒-托利多差示扫描量热仪(DSC)进行研究. 测试过程由室温升至350 ℃, 升温速率设定为5 ℃·min^–1. 取蒸干后的复合相变材料5—10 mg置于标准陶瓷坩埚中进行测量. 样品熔化后迅速将其冷却至室温以得到热流变化曲线.

3. GNS@NaNO₃复合相变材料表征及蓄热性能测量结果

3.1 GNS@NaNO₃复合相变材料的形貌

纯硝酸钠以及掺杂不同质量分数石墨烯纳米片的复合相变材料的表面形貌如图2所示. 从图2(b)—图2(e)可以看出, 石墨烯纳米片在硝酸钠中分散良好. 实际上, 当纳米片质量分数较小时, 纳米片更像是附着在硝酸钠表面. 但随着质量分数的增加, 纳米片附近的熔盐包裹在石墨烯纳米片表面形成了团簇状的致密结构.

$图 2 石墨烯纳米片@硝酸钠复合相变材料的扫描电子显微照片　(a)石墨烯纳米片; (b) 纯硝酸钠; (c) 纳米片质量分数为0.5%的石墨烯纳米片@硝酸钠; (d) 纳米片质量分数为1.0%的石墨烯纳米片@硝酸钠; (e) 纳米片质量分数为1.5%的石墨烯纳米片@硝酸钠; (f) 纳米片质量分数为2.0%的石墨烯纳米片@硝酸钠\r\nFig. 2. Scanning electron micrograph of GNS@NaNO3: (a) GNS; (b) pure NaNO3; (c) GNS@NaNO3 when GNS weight concentration is 0.5%; (d) GNS@NaNO3 when GNS weight concentration is 1.0%; (e) GNS@NaNO3 when GNS weight concentration is 1.5%; (f) GNS@NaNO3 when GNS weight concentration is 2.0%.$

图 2 石墨烯纳米片@硝酸钠复合相变材料的扫描电子显微照片　(a)石墨烯纳米片; (b) 纯硝酸钠; (c) 纳米片质量分数为0.5%的石墨烯纳米片@硝酸钠; (d) 纳米片质量分数为1.0%的石墨烯纳米片@硝酸钠; (e) 纳米片质量分数为1.5%的石墨烯纳米片@硝酸钠; (f) 纳米片质量分数为2.0%的石墨烯纳米片@硝酸钠

Fig. 2. Scanning electron micrograph of GNS@NaNO₃: (a) GNS; (b) pure NaNO₃; (c) GNS@NaNO₃ when GNS weight concentration is 0.5%; (d) GNS@NaNO₃ when GNS weight concentration is 1.0%; (e) GNS@NaNO₃ when GNS weight concentration is 1.5%; (f) GNS@NaNO₃ when GNS weight concentration is 2.0%.

下载: 全尺寸图片幻灯片

3.2 GNS@NaNO₃复合相变材料的蓄热性能

纯硝酸钠及复合相变材料的热流测量结果如图3所示. 由纯硝酸钠的吸收峰得出其熔点为308 ℃, 相应的相变焓约为174 J·g^–1[20]. 有趣的是, 掺杂石墨烯纳米片后的复合相变材料熔点略有下降, 但相变焓却出现了显著的下降. 同时可以观察到, 复合相变材料相变焓的下降与石墨烯纳米片的质量分数存在一定的关系但并未呈现出依数性, 类似的情况同样出现在Engelmann和Hentschke^[13]的研究中. 在石墨烯纳米片质量分数为1.5%时复合相变材料的相变焓下降最为明显, 由原本的174 J·g^–1降为147.72 J·g^–1.

$图 3 硝酸钠和石墨烯纳米片@硝酸钠复合相变材料的热流曲线\r\nFig. 3. Heat flow curves of pure NaNO3 and GNS@NaNO3.$

图 3 硝酸钠和石墨烯纳米片@硝酸钠复合相变材料的热流曲线

Fig. 3. Heat flow curves of pure NaNO₃ and GNS@NaNO₃.

下载: 全尺寸图片幻灯片

为更好地对比不同样品的相变焓与石墨烯纳米片质量分数之间的关系, 计算了表观蓄热效率, 定义为复合相变材料的实际相变焓与理论相变焓的比值:

$\theta = \frac{{{H_{{\text{S}},{\text{comp}}}}}}{{\left( {1 - w} \right) \times {H_{{\text{S,PCM}}}}}} \times 100{\text{%}},$

(1)

式中, θ表示表观蓄热效率, w表示石墨烯纳米片的质量分数, H_S,comp和H_S,PCM分别代表复合相变材料和硝酸钠的相变焓. 如表1所列, 复合相变材料的表观蓄热效率均低于100%, 且在石墨烯纳米片质量分数为1.5%时最低, 仅为86.19%. 结合复合相变材料的表面形貌不难得出结论, 石墨烯片附近团簇的硝酸钠熔盐并未参与熔化, 从而导致了复合相变材料相变焓的降低.

表 1 硝酸钠和石墨烯纳米片@硝酸钠复合相变材料的相变热性能

Table 1. Phase change characteristics of pure NaNO₃ and GNS@NaNO₃.

样品	T_m/℃	H_S,comp/(J·g^–1)	Theoretical H_S,comp/(J·g^–1)	θ/%
Pure NaNO₃	308.00±0.57	174.00±1.62	174	100
0.5%GNS@NaNO₃	306.87±0.40	153.12±1.48	173.13	88.44
1.0%GNS@NaNO₃	307.91±0.55	165.24±1.68	172.26	95.92
1.5%GNS@NaNO₃	307.33±0.42	147.72±1.15	171.39	86.19
2.0%GNS@NaNO₃	307.46±0.60	160.08±1.53	170.52	93.88

下载: 导出CSV

| 显示表格

4. 石墨烯纳米片掺杂对硝酸钠熔盐相变特性的影响机制

上述实验结果表明, 石墨烯纳米片的掺杂导致部分硝酸钠熔盐在纳米片附近形成了致密的团簇结构, 从而阻碍了熔盐的熔化, 体系的熔点与相变焓均有所降低. 但相变焓的降低并未呈现出对石墨烯纳米片质量分数的依数性关系, 即相变焓并非随纳米片质量分数的增加而线性降低. 为进一步探索复合材料相变焓非依数性降低的微观机理, 建立了复合相变材料的微观结构模型, 利用分子动力学模拟方法从体系内原子位置分布、分子间相互作用的类型及强度等角度进行阐释.

4.1 模型构建及方法

以石墨烯纳米片质量分数为1.0%的复合相变材料为例, 其结构模型见图4. 石墨的结构参照Materials Studio结构库中的石墨结构模型进行构建(a = b = 2.46 Å, c = 6.80 Å). 随后在石墨初始构型的基础上建立超胞, 并取单层石墨分子作为石墨烯纳米片的结构模型. 硝酸钠单胞取自晶体库, 随后对其超胞得到硝酸钠结构模型. 硝酸钠纳米结构模型的可靠性通过对比硝酸钠熔点的模拟结果与实验结果证实^[21]. 为了模拟石墨烯纳米片在硝酸钠中的均匀弥散, 复合材料模型的构建采取将边长为10 Å的石墨烯纳米片随机引入硝酸钠结构模型的方法. 在构建掺杂不同质量分数纳米片的复合材料模型时, 保持体系的大小为60 Å × 60 Å × 60 Å不变, 在模拟盒子中按照特定的质量百分比填充硝酸钠以及石墨烯纳米片.

$图 4 石墨烯纳米片质量分数为1.0%的石墨烯纳米片@硝酸钠复合相变材料结构模型\r\nFig. 4. Structural model of GNS@NaNO3 when GNS weight concentration is 1.0%.$

图 4 石墨烯纳米片质量分数为1.0%的石墨烯纳米片@硝酸钠复合相变材料结构模型

Fig. 4. Structural model of GNS@NaNO₃ when GNS weight concentration is 1.0%.

下载: 全尺寸图片幻灯片

本文所有模拟工作均通过LAMMPS平台开展^[22]. 采用分子动力学模拟探究纯硝酸钠及石墨烯纳米片质量分数为0.5%, 1.0%, 1.5%, 2.0%的复合相变材料相变热特性变化的微观机理. 选用周期性边界条件, 时间步长设置为0.1 fs, 截断半径设置为11 Å. 体系首先在热等压系综(NPT)中弛豫100 ps以达到初始平衡, 弛豫温度设定为600 K. 相互作用能、分子数密度分布等计算结果的获得通过在微正则系综(NVE)中处理200 ps实现.

硝酸钠分子间的相互作用以Buckingham势和Coulombic势修饰, 库仑相互作用的长程修正计算使用PPPM方法^[10], 表示为

$v = {A_{ij}}{{\text{e}}^{ - \frac{{{r_{ij}}}}{{{\rho _{ij}}}}}} - \frac{{{C_{ij}}}}{{r_{ij}^6}} + \frac{{{q_i}{q_j}}}{{{r_{ij}}}} ,$

(2)

式中, A, ρ, C是长度与能量相关的常数项; v是相互作用能; r是点i和j之间的距离; q_i和q_j分别是点i和j处的原子电荷量. 与此同时, 为保持硝酸盐的平面分布, 柔性硝酸根分子内部的相互作用(键、角和反常二面体势)以harmonic势^[13]来描述, 表示为以(3)式计算的键以及以(4)式计算的角和反常二面体相互作用:

${V_{\text{b}}} = {K_{\text{b}}}{(r - {r_0})^2} ,$

(3)

$V = {K_\theta }{(\theta - {\theta _0})^2} + {K_{{\text{UB}}}}{(r - {r_{{\text{UB}}}})^2} ,$

(4)

式中, r和θ分别代表键的长度和角度; K_b, K_θ和K_UB分别为键长、键角以及Urey Bradley项的力常数; r₀和θ₀分别为平衡态下的键长以及键角; r_UB代表Urey Bradley项的平衡距离. 采用AIREBO力场^[23]来描述石墨烯纳米片内部C—C之间的相互作用力, 硝酸钠与石墨烯之间的相互作用可以用Lennard-Jones (L-J)势函数^[24]体现, 公式如下:

${E_{ij}} = 4{\varepsilon _{ij}}\left[ {{{\left( {\frac{{{\sigma _{ij}}}}{{{r_{ij}}}}} \right)}^{12}} - {{\left( {\frac{{{\sigma _{ij}}}}{{{r_{ij}}}}} \right)}^6}} \right] .$

(5)

根据标准Lorenz-Bethelot混合原则, ε_ij采用几何平均值, σ_ij采用算术平均值, 表达式如下:

${\varepsilon _{ij}} = \sqrt {{\varepsilon _{ii}} \cdot {\varepsilon _{jj}}} ,$

(6)

${\sigma _{ij}} = \frac{1}{2}\left( {{\sigma _{ii}} + {\sigma _{jj}}} \right) ,$

(7)

式中, ε_ij代表原子i与j之间的相互作用强度, σ_ij是原子i和j之间的距离^[25,26]. 详细的力场参数列于补充材料(online).

4.2 复合相变材料的均方位移

通过计算复合材料均方位移(MSD)的变化对其熔点下降的现象给出相应机理解释. 均方位移可以在一定程度上反映纳米增强剂与相变材料之间的相互作用, 可通过下式计算得到:

${\text{MSD}} = |{{\boldsymbol{r}}_i}(t) - {{\boldsymbol{r}}_i}(0){|^2} ,$

(8)

式中, r_i(t)和r_i(0)分别表示分子i在时刻t和时刻0的位置向量.

图5所示为不同纳米片质量分数的复合材料在300 K下的均方位移随模拟时间的变化. 模拟过程在NPT系综下进行, 保持各体系的温度一致, 对体系内原子赋予初始速度后弛豫100 ns, 期间输出体系的MSD.

$图 5 石墨烯纳米片@硝酸钠复合相变材料300 K下的均方位移\r\nFig. 5. MSD of GNS@NaNO3 composites at 300 K.$

图 5 石墨烯纳米片@硝酸钠复合相变材料300 K下的均方位移

Fig. 5. MSD of GNS@NaNO₃ composites at 300 K.

下载: 全尺寸图片幻灯片

从图5可以看出, 随着模拟过程的进行, 各体系的均方位移逐渐增大. 且随着石墨烯纳米片质量分数的增加, 体系均方位移的增长速率显著提升. 自扩散系数D是随温度T发生变化的函数, 该系数的值越大, 说明体系内粒子扩散得越快即体系内粒子的流动性越强. 考虑到均方位移MSD计算过程中已对体系内的粒子数做了平均, 自扩散系数D可简化为MSD/6. 由此可知, 均方位移斜率越大意味着自扩散系数越大. 这说明纳米片的引入削弱了硝酸钠内部分子间的相互作用, 体系内硝酸钠分子无规则运动加剧, 平均自由程增大, 原子迁移率增高, 分子间碰撞频率增大, 从而降低了其熔点.

4.3 石墨烯纳米片附近硝酸钠熔盐的分布

首先通过600 K下体系内的原子位置分布来从微观层面描述硝酸钠熔盐在石墨烯纳米片附近的团簇情况. 以石墨烯纳米片作为坐标原点. 从图6可以看出, 有更多的硝酸钠分子出现在石墨烯纳米片附近, 且存在2—4 Å的致密层. 这与实验观察到的大量硝酸钠分子团簇在石墨烯纳米片表面的现象一致, 表明确实存在部分硝酸钠熔盐未发生熔化行为, 从而导致了复合相变材料的实际相变焓低于理论值的现象.

$图 6 石墨烯纳米片附近的原子分布(纳米片质量分数为1.5%的石墨烯纳米片@硝酸钠)\r\nFig. 6. NaNO3 atomic distribution close to GNS (GNS@NaNO3 when GNS weight concentration is 1.5%).$

图 6 石墨烯纳米片附近的原子分布(纳米片质量分数为1.5%的石墨烯纳米片@硝酸钠)

Fig. 6. NaNO₃ atomic distribution close to GNS (GNS@NaNO₃ when GNS weight concentration is 1.5%).

下载: 全尺寸图片幻灯片

同时可以观察到, 在不同质量分数下, 石墨烯纳米片附近熔盐致密层的厚度存在一定差异. 我们统计了不同纳米片质量分数下复合相变材料内部硝酸钠致密层到石墨烯纳米片的质心等效距离, 结果如表2所列.

表 2 各原子到石墨烯纳米片的质心等效距离

Table 2. Centroid equivalent distance of atoms in NaNO₃ dense layer to GNS.

原子种类	Centroid equivalent distance/Å
原子种类	0.5%	1.0%	1.5%	2.0%
Na	3.31	3.11	3.42	3.92
N	3.23	3.63	4.17	4.53
O	3.23	3.55	4.23	5.10

下载: 导出CSV

| 显示表格

考虑到石墨烯纳米片与硝酸钠之间仅存在以L-J势修饰的范德瓦耳斯相互作用, 二者之间的作用力强度与其距离密切相关. 硝酸钠熔盐各原子与石墨烯纳米片之间L-J势的平衡位置分别为r_0C-Na = 3.071 Å, r_0C-N = 3.656 Å, r_0C-O = 3.283 Å, 势阱位置r_e = 1.123r₀分别为r_eC-Na = 3.45 Å, r_eC-N = 4.11 Å, r_eC-O = 3.68 Å. 石墨烯纳米片与硝酸钠之间相互作用的强弱随着二者之间距离的变化而发生变化. 图7为石墨烯纳米片质量分数为1.0%和1.5%的复合相变材料内部原子间L-J相互作用示意图. 通过对比两种掺杂比石墨烯纳米片附近硝酸钠致密层的质心等效距离在L-J势函数曲线的位置得出结论: 当石墨烯纳米片质量分数为1.0%时, C-N间相互作用为斥力, C-Na和C-O间为极弱引力, 此时二者之间吸引力最小, 石墨烯纳米片对硝酸钠相变潜热的影响最弱; 当石墨烯纳米片质量分数为1.5%时, C-N, C-Na和C-O间距离均接近势阱位置, 存在较强引力, 此时二者间吸引力最大, 石墨烯纳米片对硝酸钠相变潜热的影响最强.

$图 7 石墨烯纳米片@硝酸钠原子间L-J相互作用示意　(a) 石墨烯纳米片质量分数为1.0% 的石墨烯纳米片@硝酸钠; (b) 石墨烯纳米片质量分数为1.5%的石墨烯纳米片@硝酸钠\r\nFig. 7. L-J interaction between atoms in GNS@NaNO3: (a) GNS@NaNO3 when GNS weight concentration is 1.0%; (b) GNS@NaNO3 GNS weight concentration is 1.5%.$

图 7 石墨烯纳米片@硝酸钠原子间L-J相互作用示意　(a) 石墨烯纳米片质量分数为1.0% 的石墨烯纳米片@硝酸钠; (b) 石墨烯纳米片质量分数为1.5%的石墨烯纳米片@硝酸钠

Fig. 7. L-J interaction between atoms in GNS@NaNO₃: (a) GNS@NaNO₃ when GNS weight concentration is 1.0%; (b) GNS@NaNO₃ GNS weight concentration is 1.5%.

下载: 全尺寸图片幻灯片

也就是说, 石墨烯纳米片与硝酸钠熔盐之间存在的相互吸引作用导致了熔盐致密层的形成. 而掺杂不同质量分数纳米片的复合材料体系内致密层的厚度不同, 这使得二者之间的相互作用强度存在一定的差异, 从而导致了复合材料相变焓不同程度的降低.

与此同时, 本文计算了复合体系内石墨烯纳米片与硝酸钠之间的相互作用能. 相互作用能常用于表征材料间的相互吸引或排斥的强度, 仅与材料的本征特性相关. 计算公式^[27]如下所示:

${E_{{\text{inter}}}} = {E_{{\text{total}}}} - {E_{{\text{GNS}}}} - {E_{{\text{NaN}}{{\text{O}}_3}}} \text{, }$

(9)

式中E_total代表模拟体系的总势能, E_GNS和 $E_{\rm NaNO_3}$ 分别代表石墨烯纳米片和硝酸钠的势能. 从图8可以看出, 石墨烯纳米片与硝酸钠之间存在较强的相互吸引作用且随着石墨烯纳米片质量分数的增加, 纳米片与硝酸钠之间相互作用能的变化趋势与相变焓恰恰相反. 在纳米片质量分数为1.0%时二者之间的相互吸引作用最弱, 质量分数为1.5%时相互吸引作用最强. 由此可以推断, 纳米片与硝酸钠熔盐之间相互作用能的存在及其强弱是复合相变材料相变焓降低的根本原因.

$图 8 石墨烯纳米片与硝酸钠之间的相互作用能\r\nFig. 8. Interaction energy between GNS and NaNO3.$

图 8 石墨烯纳米片与硝酸钠之间的相互作用能

Fig. 8. Interaction energy between GNS and NaNO₃.

下载: 全尺寸图片幻灯片

5. 结　论

本文制备了质量分数分别为0%, 0.5%, 1%, 1.5%, 2%的石墨烯纳米片@硝酸钠复合相变材料, 采用SEM对复合材料进行表征, 并通过DSC测试得到了复合材料的熔点及相变焓. 开展分子动力学模拟并分析了掺杂石墨烯纳米片对复合相变材料相变热特性的影响机理和作用规律, 主要结论如下:

1)石墨烯纳米片附近存在团簇状的硝酸钠熔盐, 且随着纳米片质量分数的增加, 石墨烯纳米片@硝酸钠复合相变材料的熔点略有下降, 相变焓呈现非依数性降低. 纳米片质量分数为1.5%时下降最为明显, 仅为理论相变焓的86.19%. 这说明团簇在纳米片附近的熔盐致密层未发生熔化.

2)石墨烯纳米片的引入削弱了硝酸钠内部分子间的相互作用, 体系内硝酸钠分子无规则运动加剧, 平均自由程增大, 原子迁移率增高, 分子间碰撞频率增大, 从而降低了其熔点. 纳米片与硝酸钠之间较强的范德瓦耳斯吸引作用是纳米片附近2—4 Å的熔盐致密层形成的根本原因. 且随着石墨烯纳米片质量分数的增加, 熔盐致密层与石墨烯纳米片之间的质心等效距离逐渐增大, 二者之间的相互吸引作用呈现先增强后减弱的趋势. 在纳米片质量分数为1.5%时二者间的质心等效距离最接近势阱位置, 相互吸引作用最强, 熔盐分子的运动受限最为严重, 因此复合材料的相变焓降低最为明显.

3)基于以上两点结论, 熔盐与纳米增强剂之间的相互作用类型及其势阱位置是决定复合后熔盐致密层厚度及相变焓降低程度的根本原因, 可借助相互作用能的计算指导纳米增强剂的掺杂比选择, 最大程度避免纳米增强剂的引入造成的芯材团簇和相变焓损失, 同时一定程度上减少复合相变材料的制备成本.

Supplements

158801-20220354补充材料.pdf

References

[1]	Lee D, Jo B 2021 Int. J. Energy Res. 45 3231 Google Scholar
[2]	Xiong Y X, Wang Z Y, Sun M Y, Wu Y T, Xu P, Qian X, Li C, Ding Y L, Ma C F 2021 Int. J. Energy Res. 45 5248 Google Scholar
[3]	Kumar N, Hirschey J, LaClair T J, Gluesenkamp K R, Graham S 2019 J. Energy Storage 24 100794 Google Scholar
[4]	Choi D H, Lee J, Hong H, Kang Y T 2014 Int J. Refrig. 42 112 Google Scholar
[5]	Wang M R, Kang Q J, Pan N 2009 Appl. Therm. Eng. 29 418 Google Scholar
[6]	Kim Y S, Kim D, Martin K J, Yu C, Grunlan J C 2010 Macromol. Mater. Eng. 295 431 Google Scholar
[7]	Shin D, Banerjee D 2010 J. Heat Transfer 133 024501 Google Scholar
[8]	Shin D, Banerjee D 2011 Int. J. Heat Mass Transfer 54 1064 Google Scholar
[9]	Jo B, Banerjee D 2014 Acta Mater. 75 80 Google Scholar
[10]	Qiao G, Alexiadis A, Ding Y L 2017 Powder Technol. 314 660 Google Scholar
[11]	Xiong Y X, Wang Z Y, Xu P, Chen H B, Wu Y T 2019 Energy Procedia 158 5551 Google Scholar
[12]	Li Z, Li B R, Du X Z, Wu H W 2020 Renewable Energy 146 816 Google Scholar
[13]	Engelmann S, Hentschke R 2019 Sci. Rep. 9 1 Google Scholar
[14]	Li J F, Lu W, Zeng Y B, Luo Z P 2014 Sol. Energy Mater. Sol. Cells 128 48 Google Scholar
[15]	Liu Y S, Yang Y Z 2017 Appl. Therm. Eng. 124 533 Google Scholar
[16]	Yu Q, Lu Y W, Zhang C C, Wu Y T, Sunden B 2019 Sol. Energy Mater. Sol. Cells 201 110055 Google Scholar
[17]	Park S, Ruoff R S 2009 Nat. Nanotechnol. 4 217 Google Scholar
[18]	Morelos-Gomez A, Terashima S, Yamanaka A, Cruz-Silva R, Ortiz-Medina J, Sánchez-Salas R, Fajardo-Díaz J L, Muñoz-Sandoval E, López-Urías F, Takeuchi K 2021 Carbon 181 118 Google Scholar
[19]	Yan X X, Zhao H B, Feng Y H, Qiu L, Lin L, Zhang X X, Ohara T 2022 Compos. B Eng. 228 109435 Google Scholar
[20]	D’Aguanno B, Karthik M, Grace A, Floris A 2018 Sci. Rep. 8 10485 Google Scholar
[21]	Jayaraman S, Thompson A P, von Lilienfeld O A, Maginn E J 2010 Ind. Eng. Chem. Res. 49 559 Google Scholar
[22]	Plimpton S 1995 J. Comput. Phys. 117 1 Google Scholar
[23]	Stuart S J, Tutein A B, Harrison J A 2000 J. Chem. Phys. 112 6472 Google Scholar
[24]	Jones J E 1924 Proc. R. Soc. Lond. A 106 463 Google Scholar
[25]	Xiong Y H, Wu H, Gao J S, Chen W, Zhang J C, Yue Y N 2019 Acta Phys. Chim. Sin. 35 1150 Google Scholar
[26]	Li Z, Cui L, Li B R, Du X Z 2020 Int. J. Heat Mass Transfer 153 119578 Google Scholar
[27]	Yang M J, Stipp S S, Harding J 2008 Cryst. Growth Des. 8 4066 Google Scholar

Cited By

期刊类型引用(1)

陈海生，李泓，徐玉杰，陈满，王亮，戴兴建，徐德厚，唐西胜，李先锋，胡勇胜，马衍伟，刘语，苏伟，王青松，陈军，卓萍，肖立业，周学志，冯自平，蒋凯，尉海军，唐永炳，陈人杰，刘亚涛，张宇鑫，林曦鹏，郭欢，张涵，张长昆，胡东旭，容晓晖，张熊，金凯强，姜丽华，彭煜民，刘世奇，朱轶林，王星，周鑫，欧学武，庞全全，俞振华，刘为，岳芬，李臻，宋振，王志峰，宋文吉，林海波，李杰才，易斌，李福军，潘新慧，李丽，马一鸣，李煌. 2022年中国储能技术研究进展. 储能科学与技术. 2023(05): 1516-1552 .

百度学术

其他类型引用(1)

图 1 石墨烯纳米片@硝酸钠复合相变材料的制备流程

Figure 1. Preparation process of GNS@NaNO₃.

DownLoad: Full-Size Img PowerPoint

Figure 2. Scanning electron micrograph of GNS@NaNO₃: (a) GNS; (b) pure NaNO₃; (c) GNS@NaNO₃ when GNS weight concentration is 0.5%; (d) GNS@NaNO₃ when GNS weight concentration is 1.0%; (e) GNS@NaNO₃ when GNS weight concentration is 1.5%; (f) GNS@NaNO₃ when GNS weight concentration is 2.0%.

DownLoad: Full-Size Img PowerPoint

图 3 硝酸钠和石墨烯纳米片@硝酸钠复合相变材料的热流曲线

Figure 3. Heat flow curves of pure NaNO₃ and GNS@NaNO₃.

DownLoad: Full-Size Img PowerPoint

图 4 石墨烯纳米片质量分数为1.0%的石墨烯纳米片@硝酸钠复合相变材料结构模型

Figure 4. Structural model of GNS@NaNO₃ when GNS weight concentration is 1.0%.

DownLoad: Full-Size Img PowerPoint

图 5 石墨烯纳米片@硝酸钠复合相变材料300 K下的均方位移

Figure 5. MSD of GNS@NaNO₃ composites at 300 K.

DownLoad: Full-Size Img PowerPoint

图 6 石墨烯纳米片附近的原子分布(纳米片质量分数为1.5%的石墨烯纳米片@硝酸钠)

Figure 6. NaNO₃ atomic distribution close to GNS (GNS@NaNO₃ when GNS weight concentration is 1.5%).

DownLoad: Full-Size Img PowerPoint

Figure 7. L-J interaction between atoms in GNS@NaNO₃: (a) GNS@NaNO₃ when GNS weight concentration is 1.0%; (b) GNS@NaNO₃ GNS weight concentration is 1.5%.

DownLoad: Full-Size Img PowerPoint

图 8 石墨烯纳米片与硝酸钠之间的相互作用能

Figure 8. Interaction energy between GNS and NaNO₃.

DownLoad: Full-Size Img PowerPoint

表 1 硝酸钠和石墨烯纳米片@硝酸钠复合相变材料的相变热性能

Table 1. Phase change characteristics of pure NaNO₃ and GNS@NaNO₃.

样品	T_m/℃	H_S,comp/(J·g^–1)	Theoretical H_S,comp/(J·g^–1)	θ/%
Pure NaNO₃	308.00±0.57	174.00±1.62	174	100
0.5%GNS@NaNO₃	306.87±0.40	153.12±1.48	173.13	88.44
1.0%GNS@NaNO₃	307.91±0.55	165.24±1.68	172.26	95.92
1.5%GNS@NaNO₃	307.33±0.42	147.72±1.15	171.39	86.19
2.0%GNS@NaNO₃	307.46±0.60	160.08±1.53	170.52	93.88

DownLoad: CSV

表 2 各原子到石墨烯纳米片的质心等效距离

Table 2. Centroid equivalent distance of atoms in NaNO₃ dense layer to GNS.

原子种类	Centroid equivalent distance/Å
原子种类	0.5%	1.0%	1.5%	2.0%
Na	3.31	3.11	3.42	3.92
N	3.23	3.63	4.17	4.53
O	3.23	3.55	4.23	5.10

DownLoad: CSV

[1]	Lee D, Jo B 2021 Int. J. Energy Res. 45 3231 Google Scholar
[2]	Xiong Y X, Wang Z Y, Sun M Y, Wu Y T, Xu P, Qian X, Li C, Ding Y L, Ma C F 2021 Int. J. Energy Res. 45 5248 Google Scholar
[3]	Kumar N, Hirschey J, LaClair T J, Gluesenkamp K R, Graham S 2019 J. Energy Storage 24 100794 Google Scholar
[4]	Choi D H, Lee J, Hong H, Kang Y T 2014 Int J. Refrig. 42 112 Google Scholar
[5]	Wang M R, Kang Q J, Pan N 2009 Appl. Therm. Eng. 29 418 Google Scholar
[6]	Kim Y S, Kim D, Martin K J, Yu C, Grunlan J C 2010 Macromol. Mater. Eng. 295 431 Google Scholar
[7]	Shin D, Banerjee D 2010 J. Heat Transfer 133 024501 Google Scholar
[8]	Shin D, Banerjee D 2011 Int. J. Heat Mass Transfer 54 1064 Google Scholar
[9]	Jo B, Banerjee D 2014 Acta Mater. 75 80 Google Scholar
[10]	Qiao G, Alexiadis A, Ding Y L 2017 Powder Technol. 314 660 Google Scholar
[11]	Xiong Y X, Wang Z Y, Xu P, Chen H B, Wu Y T 2019 Energy Procedia 158 5551 Google Scholar
[12]	Li Z, Li B R, Du X Z, Wu H W 2020 Renewable Energy 146 816 Google Scholar
[13]	Engelmann S, Hentschke R 2019 Sci. Rep. 9 1 Google Scholar
[14]	Li J F, Lu W, Zeng Y B, Luo Z P 2014 Sol. Energy Mater. Sol. Cells 128 48 Google Scholar
[15]	Liu Y S, Yang Y Z 2017 Appl. Therm. Eng. 124 533 Google Scholar
[16]	Yu Q, Lu Y W, Zhang C C, Wu Y T, Sunden B 2019 Sol. Energy Mater. Sol. Cells 201 110055 Google Scholar
[17]	Park S, Ruoff R S 2009 Nat. Nanotechnol. 4 217 Google Scholar
[18]	Morelos-Gomez A, Terashima S, Yamanaka A, Cruz-Silva R, Ortiz-Medina J, Sánchez-Salas R, Fajardo-Díaz J L, Muñoz-Sandoval E, López-Urías F, Takeuchi K 2021 Carbon 181 118 Google Scholar
[19]	Yan X X, Zhao H B, Feng Y H, Qiu L, Lin L, Zhang X X, Ohara T 2022 Compos. B Eng. 228 109435 Google Scholar
[20]	D’Aguanno B, Karthik M, Grace A, Floris A 2018 Sci. Rep. 8 10485 Google Scholar
[21]	Jayaraman S, Thompson A P, von Lilienfeld O A, Maginn E J 2010 Ind. Eng. Chem. Res. 49 559 Google Scholar
[22]	Plimpton S 1995 J. Comput. Phys. 117 1 Google Scholar
[23]	Stuart S J, Tutein A B, Harrison J A 2000 J. Chem. Phys. 112 6472 Google Scholar
[24]	Jones J E 1924 Proc. R. Soc. Lond. A 106 463 Google Scholar
[25]	Xiong Y H, Wu H, Gao J S, Chen W, Zhang J C, Yue Y N 2019 Acta Phys. Chim. Sin. 35 1150 Google Scholar
[26]	Li Z, Cui L, Li B R, Du X Z 2020 Int. J. Heat Mass Transfer 153 119578 Google Scholar
[27]	Yang M J, Stipp S S, Harding J 2008 Cryst. Growth Des. 8 4066 Google Scholar

[1]	Liu Wang-Wang, Zhang Ke-Xue, Wang Jun, Xia Guo-Dong. Simulation study of drag force characteristics of nanoparticles in transition regime. Acta Physica Sinica, 2024, 73(7): 075101. doi: 10.7498/aps.73.20231861
[2]	Liu Xiu-Cheng, Yang Zhi, Guo Hao, Chen Ying, Luo Xiang-Long, Chen Jian-Yong. Molecular dynamics simulation of thermal conductivity of diamond/epoxy resin composites. Acta Physica Sinica, 2023, 72(16): 168102. doi: 10.7498/aps.72.20222270
[3]	Ming Zhi-Fei, Song Hai-Yang, An Min-Rong. Mechanical behavior of graphene magnesium matrix composites based on molecular dynamics simulation. Acta Physica Sinica, 2022, 71(8): 086201. doi: 10.7498/aps.71.20211753
[4]	Yang Ya-Fan, Wang Jian-Zhou, Shang Xiang-Yu, Wang Tao, Sun Shu-Yu. Molecular dynamics simulation of swelling properties of Ca-montmorillonite at high temperatures. Acta Physica Sinica, 2022, 71(4): 043102. doi: 10.7498/aps.71.20211565
[5]	Liu Qing-Yang, Xu Qing-Song, Li Rui. Effect of N-doping on tribological properties of graphene by molecular dynamics simulation. Acta Physica Sinica, 2022, 71(14): 146801. doi: 10.7498/aps.71.20212309
[6]	Pan Ling, Zhang Hao, Lin Guo-Bin. Molecular dynamics simulation on dynamic behaviors of nanodroplets impinging on solid surfaces decorated with nanopillars. Acta Physica Sinica, 2021, 70(13): 134704. doi: 10.7498/aps.70.20210094
[7]	Shi Chao, Lin Chen-Sen, Chen Shuo, Zhu Jun. Molecular dynamics simulation of characteristic water molecular arrangement on graphene surface and wetting transparency of graphene. Acta Physica Sinica, 2019, 68(8): 086801. doi: 10.7498/aps.68.20182307
[8]	Yuan Si-Wei, Feng Yan-Hui, Wang Xin, Zhang Xin-Xin. Molecular dynamics simulation of thermal conductivity of mesoporous α-Al2O3. Acta Physica Sinica, 2014, 63(1): 014402. doi: 10.7498/aps.63.014402
[9]	Si Li-Na, Wang Xiao-Li. A molecular dynamics study on adhesive contact processes of surfaces with nanogrooves. Acta Physica Sinica, 2014, 63(23): 234601. doi: 10.7498/aps.63.234601
[10]	Qi Yu, Qu Chang-Rong, Wang Li, Fang Teng. Liquid-liquid phase segregation process of Fe50Cu50 melt by molecular dynamics simulation. Acta Physica Sinica, 2014, 63(4): 046401. doi: 10.7498/aps.63.46401
[11]	Jian Zeng-Yun, Gao A-Hong, Chang Fang-E, Tang Bo-Bo, Zhang Long, Li Na. Molecular dynamics simulation of the critical and subcritical nuclei during solidification of nickel melt. Acta Physica Sinica, 2013, 62(5): 056102. doi: 10.7498/aps.62.056102
[12]	Chen Qing, Sun Min-Hua. Molecular dynamics simulation of isothermal crystallization dynamics in Cu nanocluster. Acta Physica Sinica, 2013, 62(3): 036101. doi: 10.7498/aps.62.036101
[13]	Sun Wei-Feng, Wang Xuan. Molecular dynamics simulation study of polyimide/copper-nanoparticle composites. Acta Physica Sinica, 2013, 62(18): 186202. doi: 10.7498/aps.62.186202
[14]	Li Ming-Lin, Lin Fan, Chen Yue. Study on the mechanical properties of carbon nanocones using molecular dynamics simulation. Acta Physica Sinica, 2013, 62(1): 016102. doi: 10.7498/aps.62.016102
[15]	Xia Dong, Wang Xin-Qiang. Structures and melting behaviors of ultrathin platinum nanowires. Acta Physica Sinica, 2012, 61(13): 130510. doi: 10.7498/aps.61.130510
[16]	Quan Wei-Long, Li Hong-Xuan, Ji Li, Zhao Fei, Du Wen, Zhou Hui-Di, Chen Jian-Min. Molecular dynamical simulation on the mechanical property of hydrogenated diamond-like carbon films. Acta Physica Sinica, 2010, 59(8): 5687-5691. doi: 10.7498/aps.59.5687
[17]	Cao Li-Xia, Shang Jia-Xiang, Zhang Yue. Molecular dynamics simulation of stress-induced martensitic phase transformation in NiAl. Acta Physica Sinica, 2009, 58(10): 7307-7312. doi: 10.7498/aps.58.7307
[18]	Xie Fang, Zhu Ya-Bo, Zhang Zhao-Hui, Zhang Lin. Molecular dynamics simulation of multi-wall carbon nanotube oscillators. Acta Physica Sinica, 2008, 57(9): 5833-5837. doi: 10.7498/aps.57.5833
[19]	Meng Li-Jun, Zhang Kai-Wang, Zhong Jian-Xin. Molecular dynamics simulation of formation of silicon nanoparticles on surfaces of carbon nanotubes. Acta Physica Sinica, 2007, 56(2): 1009-1013. doi: 10.7498/aps.56.1009
[20]	Li Rui, Hu Yuan-Zhong, Wang Hui, Zhang Yu-Jun. Molecular dynamics simulation of motion of single-walled carbon nanotubes on graphite substrate. Acta Physica Sinica, 2006, 55(10): 5455-5459. doi: 10.7498/aps.55.5455

期刊类型引用(1)

百度学术

其他类型引用(1)

158801-20220354补充材料.pdf

Catalog

Vol. 71, Issue 15 - 2022-08-05

Metrics

Abstract views: 4900
PDF Downloads: 92
Cited By: 2

姓名
邮箱
手机号码
标题
留言内容
验证码

Search

Article

留言板