Effect of Morphological Characteristics of Aggregates on Thermal Properties of Molten Salt Nanofluids

: Molten salt-based nanofluid is a thermal storage and heat transfer medium for concentrated solar thermal power plants formed by adding nanoparticles to molten salt, which can enhance the thermal performance of molten salt. However, the nanoparticles tend to aggregate in nanofluids, causing changes in thermal properties. In this work, molecular dynamics simulations were used to study the effect of morphological characteristics of aggregates on the thermal conductivity and specific heat capacity of molten salt-based nanofluids. The results show that the aggregated nanoparticles cause a greater increase in thermal conductivity and specific heat capacity than dispersed nanoparticles. Additionally, the increase in fractal dimension leads to thermal conductivity reduction, while there is no clear correlation between the fractal dimension and specific heat capacity. New insights into the thermal properties of aggregated nanofluids are provided by analyzing the contribution of material components, heat flux fluctuation modes, and energy compositions. It is found that the thermal conductivity of aggregated nanofluids is mainly dominated by the base liquid and collision term. However, the specific heat is not related to the variation in the contribution of different energy compositions. Moreover, compared to the dispersed nanofluid, the increased specific heat capacity of aggregated nanofluids is attributed to the thicker semi-solid layer. This study provides guidance for the design and control of the thermal properties of molten salt-based nanofluids.


Introduction
Solar energy is abundant and emits very few pollutants, making it an important energy source to alleviate energy problems [1].Concentrated solar thermal power (CSP) is a promising solar energy technology.Thermal energy storage (TES) is a necessary part of a CSP plant, where stored heat can be used for continuous operation of the CSP plant during the night and on cloudy days [2] and, thus, help to improve power supply quality [3].Molten salts are widely used as the thermal storage material and heat transfer medium in CSP systems due to their high latent heat, low vapor pressure, wide operating temperature range, and low cost [4].However, the low thermal conductivity of molten salts reduces the charge-discharge rate and decreases the efficiency of TES systems.The low specific heat capacity leads to problems such as the larger scale of the heat storage system and the greater consumption of heat storage materials.These increase the cost and reduce the economy of CSP systems [5,6].
Numerous studies have shown that the addition of nanoparticles can enhance the thermal performance of molten salts [7][8][9][10][11][12][13][14][15].However, nanoparticles in molten salt-based nanofluids tend to aggregate and settle easily [16][17][18], affecting the thermal performance of nanofluids.Only a few studies have focused on the effect of nanoparticle aggregation on thermal conductivity.In addition, the effect of aggregation on specific heat capacity is often overlooked.It has been found that nanoadditive aggregation could promote heat conduction in water containing metal oxide nanoparticles and carbon nanotubes [19].Similarly enhanced thermal conductivity has been reported for Fe 3 O 4 -water nanofluids.When the aggregates are properly oriented, and the volume fraction of aggregates accounts for 65% of the total nanoparticle, the enhancement of thermal conductivity is the highest [20].The enhancement in thermal conductivity is attributed to the formation of the heat conduction pathways, which are formed by aggregates and have low thermal resistance [7,21].In contrast, the research of Al 2 O 3 -water nanofluids [11], ZnO-propylene glycol nanofluids [10], and CuO-water as well as CuO-ethylene glycol nanofluids [9] found that aggregation limited the thermal conductivity of nanofluids.This may be attributed to the reduction in the specific surface area of nanoparticles, limiting the effective area of thermal interaction between the base liquid and nanoparticles.To date, the impact of aggregation on the thermal conductivity of nanofluids remains a controversial issue.
On the other hand, different morphological characteristics of aggregates, such as fractal dimension, size, and shape, have an influence on the thermal properties of aggregated nanofluids [22].The nanoparticle agglomerates can be characterized by fractal dimension, which represents the tightness of nanoparticle aggregation.The effect of the fractal dimension of nanoparticles on the specific heat capacity has been rarely studied, and only a few articles have investigated the fractal dimension's influence on thermal conductivity.Wang et al. [23] studied Cu-Ar nanofluid through molecular dynamics simulation and found that the lower the fractal dimension of the nanoparticle, the greater the thermal conductivity.Li et al. [24] and Du et al. [25] also found the same conclusion in Cu-water nanofluid.The low fractal dimension is more likely to produce effective heat channels and, hence, promotes the heat transfer process inside the aggregated nanofluids [23,24].However, all of the previous research concerns molecular liquids, such as Cu-Ar and Cu-water nanofluids.Whether the same conclusion can be drawn about ionic liquids, such as solar salt-based nanofluids, remains to be explored.Additionally, the underlying mechanism of the effect of agglomeration morphology on thermophysical properties is still unclear and needs to be studied further.
Solar salt (NaNO 3 :KNO 3 = 60:40 wt.%) is the most-used working fluid in CSP systems due to its low melting point and, thus, a large liquid phase range, resulting in a large operating range and minimal energy demand to remain in liquid form [26]. Therefore, it was selected as the base fluid in this work.In addition, SiO 2 nanoparticles, which have the advantages of easy fabrication, low cost, and chemical stability in solar salt, were selected as the additive.To explore the impact of the fractal dimension of aggregates on the thermal properties of nanofluids and reveal its mechanism, the thermal conductivity and specific heat capacity of SiO 2 -solar salt were calculated using molecular dynamics (MD) simulations.The deeper heat conduction and thermal energy storage mechanisms were obtained from the semi-solid layer, the contributions of different material components and heat flux fluctuation modes to thermal conductivity, and the contributions of different energy compositions to specific heat capacity.

Simulation Model and Method
The model of solar salt is shown in Figure 1a.Following the density of solar salt (NaNO 3 :KNO 3 = 60:40 wt.%) of 1.9 g/cm 3 [27], Na + , K + , and NO 3 − ions were randomly filled in a cubic box with a side length of 55 Å, which is large enough to avoid the finite size effect [28].Due to the limitations of computer resources, four spherical SiO 2 nanoparticles with a diameter of 18 Å were created, as shown in Figure 1b.Nanoparticle diameters around 2 nm have been used in previous MD simulations [29,30].The nanoparticles were randomly filled in the solar salt box to obtain a dispersed nanofluid, as shown in Figure 1c.Then, nanoparticle aggregates with different fractal dimensions were created.Figure 1d-g shows aggregated nanofluids under 8 wt.% loading of nanoparticles.The fractal dimension (D f ) represents the tightness of nanoparticle aggregation, which is given by Equation (1) [23].
Here, N P is the number of primary nanoparticles of radius R, and R g is the radius of gyration.
Here, NP is the number of primary nanoparticles of radius R, and Rg is the radius of gyration.All MD simulations were done with the open-source LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator, 2020OCT29) package [31].The time step was set to 0.5 fs, and periodic boundary conditions were applied in all three dimensions.The Buckingham potential [32] was employed to describe the interactions between ions in nitrate salts, the interactions among atoms within SiO2 nanoparticles, and the interactions between SiO2 nanoparticle atoms and nitrate salt ions [33].The Buckingham potential is expressed by Equation (2) [32].
Here, i and j represent different atoms, respectively.rij represents the distance between atoms i and j. qi and qj represent the number of charges of atoms i and j, respectively.ρij is an ionic-pair-dependent length parameter, and the coefficients Aij and Cij are two potential parameters related to energy.The parameters of the various atoms in Equation ( 2) are shown in Table 1.All MD simulations were done with the open-source LAMMPS (Large-scale Atomic/ Molecular Massively Parallel Simulator, 2020OCT29) package [31].The time step was set to 0.5 fs, and periodic boundary conditions were applied in all three dimensions.The Buckingham potential [32] was employed to describe the interactions between ions in nitrate salts, the interactions among atoms within SiO 2 nanoparticles, and the interactions between SiO 2 nanoparticle atoms and nitrate salt ions [33].The Buckingham potential is expressed by Equation (2) [32].
Here, i and j represent different atoms, respectively.r ij represents the distance between atoms i and j. q i and q j represent the number of charges of atoms i and j, respectively.ρ ij is an ionic-pair-dependent length parameter, and the coefficients A ij and C ij are two potential parameters related to energy.The parameters of the various atoms in Equation ( 2) are shown in Table 1.
Table 1.Parameters of Buckingham potential function [29].For the Buckingham potential function, the parameters of the interactions between the different atoms were calculated by the Lorentz-Berthelot mixing rule [32] In nitrate, the bonded interactions were described using the Harmonic potential [34] as given in Equation ( 3) and Charmm potential [35] as given in Equation ( 4).
In Equation (3), r is the N-O bond length, r 0 is the equilibrium bond length, and K r is the force constant.In Equation ( 4), the first term on the right of equal sign is the simple harmonic term for O-N-O bond angle bending, the second term on the right of equal sign is the O-O non-bonding repulsion term, θ is the O-N-O bond angle, θ 0 is the equilibrium bond angle, K θ is the force constant, r is the O-O distance, r UB is the equilibrium distance, and K UB is the Urey-Bradley force constant.The parameters of Harmonic potential and Charmm potential are shown in Table 2.The calculation of thermal conductivity adopts the equilibrium molecular dynamics (EMD) method.For solar salts in the CSP plant, the working temperature range is 563-861 K [36] at ambient pressure of 1 bar [37].In order to obtain an equilibrium state, the systems were equilibrated in an NPT ensemble for 0.5 ns, with a pressure of 1 bar and temperature of 573.15 K, and then under the NVT ensemble for 1 ns, followed by a final 2 ns simulation in the NVE ensemble.After that, the thermal conductivity was obtained in the NVE ensemble for another 0.5 ns.The thermal conductivity was obtained through the heat current autocorrelation function (HCACF) using the Green-Kubo formula [38], as shown in Equation (5).
Here, k B , V, and T are Boltzmann constants, the system volume and temperature, respectively.<J(0)J(t)> is the ensemble average of HCACF.J(t) denotes the heat flux of the system at the time t [39], and its expression is shown in Equation (6).
Here, v i and E i represent the velocity vector and total energy of atom i, respectively.r ij and f ij represent the distance and interaction force between atoms i and j, respectively.N denotes the total number of atoms in the system.It is worth noting that in Equation ( 6), the third term on the right of the equal sign represents the partial enthalpy.For a singlecomponent system, the partial enthalpy is 0. However, for a multi-component system such as solar salt-based nanofluids, the partial enthalpy is non-zero.This portion of heat flux does not contribute to thermal conductivity and needs to be subtracted [39].The calculation of specific heat capacity adopts the energy fluctuation method.In order to obtain an equilibrium state, the systems were equilibrated in the NPT ensemble for 0.5 ns, with the pressure of 1 bar and temperature of 573.15 K, and then under the NVT ensemble for 1 ns.After that, the specific heat capacity C was obtained in the NVT ensemble for another 1 ns.The C is calculated as Equations ( 7) and ( 8) [40].
Here, σ 2 E denotes the rise and fall of the total energy of the system under the NVT ensemble.<> represents the average value of the system with the number of atoms as the weight, and m is the atom mass.

Thermal Conductivity and Specific Heat Capacity
Figure 2a shows the thermal conductivity of solar salt and nanofluids.The obtained thermal conductivity of solar salt at 573.15 K is 0.487 W•m −1 •K −1 , which is similar to the experimental results of 0.500 W•m −1 •K −1 [41].It validates the reliability of the potential energy functions and simulation method chosen in this study.It can also be observed from Figure 2a that the addition of nanoparticles leads to an increase in thermal conductivity.A similar result has been reported for binary nitrate salt with SiO 2 nanoparticles in experiments [42].Moreover, compared with the dispersed nanoparticles, the addition of aggregated nanoparticles results in a higher thermal conductivity.This result is different from that observed in previous experiments [43][44][45].Nasiri et al. [43] and Kazemi et al. [44] found that the thermal conductivity of nanofluids with small agglomeration and low settling is larger than that of nanofluids with large agglomeration and high settling.The contradictory conclusions between simulations and experiments are attributed to different stabilities, i.e., no settling and high settling.As the nanoparticle aggregates settle, the concentration of nanoparticles in the system decreases; that is, the mass fraction of the nanofluid decreases, resulting in a decrease in thermal conductivity [46].However, in MD simulations, there is no sedimentation in aggregated nanofluids.Figure 2a also shows that the smaller the fractal dimension of the aggregates, the higher the thermal conductivity of the nanofluid.As the fractal dimension increases from 1 to 2, the thermal conductivity decreases by 14.15%.Similar results have been reported in Cu/water nanofluid [24] and Cu/Ar nanofluid [23].
Here, σ 2 E denotes the rise and fall of the total energy of the system under the NVT ensemble.<> represents the average value of the system with the number of atoms as the weight, and m is the atom mass.

Thermal Conductivity and Specific Heat Capacity
Figure 2a shows the thermal conductivity of solar salt and nanofluids.The obtained thermal conductivity of solar salt at 573.15 K is 0.487 W•m −1 •K −1 , which is similar to the experimental results of 0.500 W•m −1 •K −1 [41].It validates the reliability of the potential energy functions and simulation method chosen in this study.It can also be observed from Figure 2a that the addition of nanoparticles leads to an increase in thermal conductivity.A similar result has been reported for binary nitrate salt with SiO2 nanoparticles in experiments [42].Moreover, compared with the dispersed nanoparticles, the addition of aggregated nanoparticles results in a higher thermal conductivity.This result is different from that observed in previous experiments [43][44][45].Nasiri et al. [43] and Kazemi et al. [44] found that the thermal conductivity of nanofluids with small agglomeration and low settling is larger than that of nanofluids with large agglomeration and high settling.The contradictory conclusions between simulations and experiments are attributed to different stabilities, i.e., no settling and high settling.As the nanoparticle aggregates settle, the concentration of nanoparticles in the system decreases; that is, the mass fraction of the nanofluid decreases, resulting in a decrease in thermal conductivity [46].However, in MD simulations, there is no sedimentation in aggregated nanofluids.Figure 2a also shows that the smaller the fractal dimension of the aggregates, the higher the thermal conductivity of the nanofluid.As the fractal dimension increases from 1 to 2, the thermal conductivity decreases by 14.15%.Similar results have been reported in Cu/water nanofluid [24] and Cu/Ar nanofluid [23].
Figure 2b shows the specific heat capacities of solar salt, dispersed nanofluid, and aggregated nanofluids with different fractal dimensions at 573.15 K.The specific heat capacity of solar salt at 573.15 K is determined to be 1.610 J•g −1 •K −1 , which is close to the experimental value of 1.648 J•g −1 •K −1 obtained in the literature [47].As shown in Figure 2b, the specific heat capacity of nanofluids is increased compared to solar salt, whether adding agglomerated or dispersed nanoparticles.However, the fractal dimension of nanoparticle aggregates is not directly related to the specific heat capacity, which is different from the monotonous dependence of thermal conductivity on the fractal dimension.Figure 2b shows the specific heat capacities of solar salt, dispersed nanofluid, and aggregated nanofluids with different fractal dimensions at 573.15 K.The specific heat capacity of solar salt at 573.15 K is determined to be 1.610 J•g −1 •K −1 , which is close to the experimental value of 1.648 J•g −1 •K −1 obtained in the literature [47].As shown in Figure 2b, the specific heat capacity of nanofluids is increased compared to solar salt, whether adding agglomerated or dispersed nanoparticles.However, the fractal dimension of nanoparticle aggregates is not directly related to the specific heat capacity, which is different from the monotonous dependence of thermal conductivity on the fractal dimension.

Contributions of Different Material Components and Heat Flux Fluctuation Modes to Thermal Conductivity
In order to study the heat conduction mechanism of aggregate nanofluids, the contributions of nanoparticles and base liquid to the thermal conductivity were examined, respectively.The heat flux of nanoparticle J N and base liquid J B is shown in Equations ( 9) and ( 10) [48].Hence, according to the Green-Kubo formula, the thermal conductivity can be calculated as a sum of four terms, including the self-correlations (i.e., λ BB and λ NN ) and crossed correlations (including λ BN and λ NB ).
It is found that the contributions of crossed correlations are small, which has also been reported in Cu-Ar nanofluid [49].Therefore, Figure 3a only shows the contributions of nanoparticles and base liquid to thermal conductivity due to self-correlations.It can be seen that the contribution of base fluid is about 70%, which is much larger than aggregated nanoparticles.In addition, the contribution of base fluid decreases with the increase of fractal dimension.As the fractal dimension increases from 1.116 to 2, the contribution of λ BB decreases by 6.95%.This tendency is consistent with the thermal conductivity shown in Figure 2a.It indicates that the base liquid dominates the heat conduction of aggregated nanofluids.

Contributions of Different Energy Compositions to Specific Heat Capacity
To reveal the physical insights of specific heat capacity, the effects of different energy components on the specific heat capacity of the nanofluid were explored.The contributions of kinetic and potential energy to the specific heat capacity were calculated first.The results are shown in Figure 4a.We can see that the contribution of potential energy is larger than that of kinetic energy.We further decomposed the contribution of potential energy to the van der Waals Cvdwl, short-ranged part of coulombic interaction Ccoul, longranged part of coulombic interaction Clong, bond Cbond, angle Cangle, and improper Cimpro potentials, as shown in Figure 4b.It can be seen that the ranking of contribution is Clong < Cimpro < Cbond < Cangle < Cvdwl < Ccoul.The contribution of Ccoul is the largest and negatively correlated with the fractal dimension.As the fractal dimension increases from 1.116 to 2, the contribution of Ccoul decreases by 30.17%.Differently, the specific heat capacity changes with fractal dimensions are irregular.The correlation between specific heat and the variation in the contribution of different energy compositions needs to be further explored.We hope that our findings can stimulate interest in studying the relationship between the specific heat and energy composition contribution.On the other hand, the heat flux vector can be decomposed into three terms (i.e., the kinetic (K), potential (P), and collision (C) terms), as shown in Equations ( 11) and ( 12) [50].Therefore, the thermal conductivity can be calculated as a sum of nine terms, including self-correlations (i.e., λ KK , λ PP , and λ CC ) and crossed correlations (including λ PC , λ CP , λ PK , λ KP , λ KC, and λ CK ).
Here, α and β represent solar salt and SiO 2 nanoparticles.k and l represent different components.m i and v i are the mass and velocity of atom i, and h is the mean partial enthalpy.A is an arbitrary constant, and r c is the cutoff radius.
The results show that the contributions of crossed correlations are negligible, which is consistent with that reported in Pt/Xe nanofluid [50].Hence, the contributions due to self-correlations are given in Figure 3b.As can be seen, the contribution of λ KK is small, which means Brownian motion has minimal effect on thermal conductivity.Similar results have been found in dispersed nanofluids [30].Moreover, the enhancement in thermal conductivity of aggregated nanofluids is primarily dominated by λ CC , which is larger than 60%.The fractal dimension of aggregates shows a negative correlation with λ CC .As the fractal dimension increases from 1.116 to 2, the contribution of λ CC decreases by 5.55%.This trend is the same with the thermal conductivity shown in Figure 2a.It suggests that the collision term may be responsible for the reduced thermal conductivity of nanofluids with larger fractal dimensions.

Contributions of Different Energy Compositions to Specific Heat Capacity
To reveal the physical insights of specific heat capacity, the effects of different energy components on the specific heat capacity of the nanofluid were explored.The contributions of kinetic and potential energy to the specific heat capacity were calculated first.The results are shown in Figure 4a.We can see that the contribution of potential energy is larger than that of kinetic energy.We further decomposed the contribution of potential energy to the van der Waals C vdwl , short-ranged part of coulombic interaction C coul , long-ranged part of coulombic interaction C long , bond C bond , angle C angle, and improper C impro potentials, as shown in Figure 4b.It can be seen that the ranking of contribution is The contribution of C coul is the largest and negatively correlated with the fractal dimension.As the fractal dimension increases from 1.116 to 2, the contribution of C coul decreases by 30.17%.Differently, the specific heat capacity changes with fractal dimensions are irregular.The correlation between specific heat and the variation in the contribution of different energy compositions needs to be further explored.We hope that our findings can stimulate interest in studying the relationship between the specific heat and energy composition contribution.

Contributions of Different Energy Compositions to Specific Heat Capacity
To reveal the physical insights of specific heat capacity, the effects of different energy components on the specific heat capacity of the nanofluid were explored.The contributions of kinetic and potential energy to the specific heat capacity were calculated first.The results are shown in Figure 4a.We can see that the contribution of potential energy is larger than that of kinetic energy.We further decomposed the contribution of potential energy to the van der Waals Cvdwl, short-ranged part of coulombic interaction Ccoul, longranged part of coulombic interaction Clong, bond Cbond, angle Cangle, and improper Cimpro potentials, as shown in Figure 4b.It can be seen that the ranking of contribution is Clong < Cimpro < Cbond < Cangle < Cvdwl < Ccoul.The contribution of Ccoul is the largest and negatively correlated with the fractal dimension.As the fractal dimension increases from 1.116 to 2, the contribution of Ccoul decreases by 30.17%.Differently, the specific heat capacity changes with fractal dimensions are irregular.The correlation between specific heat and the variation in the contribution of different energy compositions needs to be further explored.We hope that our findings can stimulate interest in studying the relationship between the specific heat and energy composition contribution.

Semi-Solid Layer
The density distribution of ions around the nanoparticle was further calculated and is shown in Figure 5a.As can be seen, the first peak of density is the greatest, then the other peaks weaken.This indicates that the addition of nanoparticles results in a strong

Semi-Solid Layer
The density distribution of ions around the nanoparticle was further calculated and is shown in Figure 5a.As can be seen, the first peak of density is the greatest, then the other peaks weaken.This indicates that the addition of nanoparticles results in a strong attraction between nanoparticles and the base liquid, leading to the formation of a semi-solid layer.Destroying the microstructure of the semi-solid layer requires energy, thus increasing the specific heat capacity of the nanofluid.The thickness of the semisolid layer is measured from the first peak to the third trough [51] of density distribution.Figure 5b shows the thickness of the semi-solid layer around dispersed nanoparticles and aggregated nanoparticles.We can find that the thickness of the semi-solid layer around agglomerated nanoparticles is greater than that around dispersed nanoparticles.In addition, the quantitative relationship of the semi-solid layer thickness shown in Figure 5b is consistent with the specific heat shown in Figure 2b.It verifies that the semi-solid layer is related to the specific heat.However, the quantitative relationship of semi-solid layer thickness does not correspond to the thermal conductivity.This indicates that the change in heat conduction of aggregated nanofluids cannot be explained by the semi-solid layer.Similar results have been found in dispersed nanofluids [29,30].
Energies 2024, 17, x FOR PEER REVIEW 8 of 11 attraction between nanoparticles and the base liquid, leading to the formation of a semisolid layer.Destroying the microstructure of the semi-solid layer requires energy, thus increasing the specific heat capacity of the nanofluid.The thickness of the semi-solid layer is measured from the first peak to the third trough [51] of density distribution.Figure 5b shows the thickness of the semi-solid layer around dispersed nanoparticles and aggregated nanoparticles.We can find that the thickness of the semi-solid layer around agglomerated nanoparticles is greater than that around dispersed nanoparticles.In addition, the quantitative relationship of the semi-solid layer thickness shown in Figure 5b is consistent with the specific heat shown in Figure 2b.It verifies that the semi-solid layer is related to the specific heat.However, the quantitative relationship of semi-solid layer thickness does not correspond to the thermal conductivity.This indicates that the change in heat conduction of aggregated nanofluids cannot be explained by the semi-solid layer.Similar results have been found in dispersed nanofluids [29,30].Apart from the addition of nanoparticles, the aggregate morphology can also provide a degree of freedom to promote the heat transfer and energy conversion efficiency of molten salt.Filling nanoparticles results in an increase in both specific heat capacity and thermal conductivity of molten salt, while the aggregate morphology has different effects on heat storage and heat transfer properties.Different from the irregular dependence of specific heat capacity on aggregate morphology, there is a negative correlation between aggregate morphology and thermal conductivity.Higher thermal conductivity can be obtained for a nanofluid with smaller aggregates by chemically treating the nanoparticles [52], adding surfactants [53], and changing pH values [54].

Conclusions
In this paper, the thermal conductivity and specific heat capacity of solar salt-based nanofluids were investigated by MD simulations.The results show that the addition of SiO2 nanoparticles increases the thermal conductivity and specific heat capacity of solar salt, and the aggregated nanofluids have greater thermal conductivity and specific heat capacity than the dispersed ones.It is interesting to find that the aggregate morphological characteristics have different effects on the thermal conductivity and specific heat capacity of nanofluids.The thermal conductivity decreases with increasing fractal dimension, while the fractal dimension has no correlation with specific heat capacity.Based on the analysis of the contributions of different material components and heat flux fluctuation modes to thermal conductivity and contributions of different energy compositions to specific heat capacity, new insights into heat conduction and thermal energy storage were provided.It was found that the base fluid and collision term have the dominant impact on the heat conduction of aggregated nanofluids.Moreover, the thickness of the semi- Apart from the addition of nanoparticles, the aggregate morphology can also provide a degree of freedom to promote the heat transfer and energy conversion efficiency of molten salt.Filling nanoparticles results in an increase in both specific heat capacity and thermal conductivity of molten salt, while the aggregate morphology has different effects on heat storage and heat transfer properties.Different from the irregular dependence of specific heat capacity on aggregate morphology, there is a negative correlation between aggregate morphology and thermal conductivity.Higher thermal conductivity can be obtained for a nanofluid with smaller aggregates by chemically treating the nanoparticles [52], adding surfactants [53], and changing pH values [54].

Conclusions
In this paper, the thermal conductivity and specific heat capacity of solar salt-based nanofluids were investigated by MD simulations.The results show that the addition of SiO 2 nanoparticles increases the thermal conductivity and specific heat capacity of solar salt, and the aggregated nanofluids have greater thermal conductivity and specific heat capacity than the dispersed ones.It is interesting to find that the aggregate morphological characteristics have different effects on the thermal conductivity and specific heat capacity of nanofluids.The thermal conductivity decreases with increasing fractal dimension, while the fractal dimension has no correlation with specific heat capacity.Based on the analysis of the contributions of different material components and heat flux fluctuation modes to thermal conductivity and contributions of different energy compositions to specific heat capacity, new insights into heat conduction and thermal energy storage were provided.It was found that the base fluid and collision term have the dominant impact on the heat conduction of aggregated nanofluids.Moreover, the thickness of the semi-solid layer is related to specific heat capacity.This MD study is focused on the new phenomena and physical insights into the thermal properties of aggregated nanofluids on the nanoscale.It is expected that the findings can inspire others to pursue the thermal properties of aggregated nanofluids at meso and macro scales.

Figure 2 .
Figure 2. (a) Thermal conductivity and (b) specific heat of solar salt and nanofluids at 573.15 K.Figure 2. (a) Thermal conductivity and (b) specific heat of solar salt and nanofluids at 573.15 K.

Figure 2 .
Figure 2. (a) Thermal conductivity and (b) specific heat of solar salt and nanofluids at 573.15 K.Figure 2. (a) Thermal conductivity and (b) specific heat of solar salt and nanofluids at 573.15 K.

Figure 3 .
Figure 3. (a) Contributions of base fluid and nanoparticle to thermal conductivity.(b) Contributions of kinetic, potential, and collision terms to thermal conductivity.

Figure 3 .
Figure 3. (a) Contributions of base fluid and nanoparticle to thermal conductivity.(b) Contributions of kinetic, potential, and collision terms to thermal conductivity.

Figure 3 .
Figure 3. (a) Contributions of base fluid and nanoparticle to thermal conductivity.(b) Contributions of kinetic, potential, and collision terms to thermal conductivity.

Figure 4 .
Figure 4. (a) Contributions of the kinetic and potential to specific heat capacity.(b) Contributions of the van der Waals, short-ranged part of coulombic interaction, long-ranged part of coulombic interaction, bond, angle, and improper potentials to specific heat capacity.

Figure 4 .
Figure 4. (a) Contributions of the kinetic and potential to specific heat capacity.(b) Contributions of the van der Waals, short-ranged part of coulombic interaction, long-ranged part of coulombic interaction, bond, angle, and improper potentials to specific heat capacity.

3 ]Figure 5 .
Figure 5. Density distribution of (a) base fluid around nanoparticles and (b) statistics of the thickness of semi-solid layer.

Figure 5 .
Figure 5. Density distribution of (a) base fluid around nanoparticles and (b) statistics of the thickness of semi-solid layer.