Thermal Conductivity and Stability of Novel Aqueous Graphene Oxide–Al 2 O 3 Hybrid Nanoﬂuids for Cold Energy Storage

Featured Application: The paper not only provides method and theoretical guidance for the preparation of long-term stable high-performance hybrid nanoﬂuids, but also gives support for their e ﬃ cient applications in phase change cold energy storage. Abstract: Thermal ice storage has gained a lot of interest due to its ability as cold energy storage. However, low thermal conductivity and high supercooling degree have become major issues during thermal cycling. For reducing the cost and making full use of the advantages of the graphene oxide–Al 2 O 3 , this study proposes heat transfer enhancement of thermal ice storage using novel hybrid nanoﬂuids of aqueous graphene oxide–Al 2 O 3 . Thermal conductivity of aqueous graphene oxide–Al 2 O 3 nanoﬂuid was measured experimentally over a range of temperatures (0–70 ◦ C) and concentrations. Thermal conductivity of ice mixing with the hybrid nanoparticles was tested. The inﬂuences of pH, dispersant, ultrasonic power and ultrasonic time on the stability of the hybrid nanoﬂuids were examined. A new model for the e ﬀ ective thermal conductivity of the hybrid nanoﬂuids considering the structure and Brownian motion was proposed. The results showed that pH, dispersant, ultrasonic power level and ultrasonication duration are important factors a ﬀ ecting the stability of the hybrid nanoﬂuids tested. The optimum conditions for stability are pH = 11, 1% SDS, 375 W ultrasonic power level and 120 min ultrasonic application time. The thermal conductivity of hybrid nanoﬂuids increases with the increase of temperature and mass fraction of nanoparticles. A newly proposed thermal conductivity model considering the nanoﬂuid structure and Brownian motion can predict the thermal conductivity of hybrid nanoﬂuids reasonably well. ◦ (001) the ﬀ raction at 2 θ = ◦ , 35.15 , 37.77 , 43.35 , 52.55 , , 68.21 and 77.2 is the crystal phase of α - Al 2 O 3 (012), (104), (110), (113), (024), (116), (300) and (119), respectively, which is consistent with the standard α -Al 2 O 3 spectrum [63]. hybrid hybrid nanoparticles hybrid nanofluid conductivity of hybrid nanofluids. new model for thermal conductivity of hybrid nanofluid of flat rate interface thermal resistance of the graphene oxide, radius of Al2O3, fluid viscosity, Brownian motion. Results presented in this of the influence of factors on the enhancement of thermal conductivity of hybrid nanofluids. (4) Comparison of thermal conductivity of the GO nanoﬂuid with the hybrid nanoﬂuid displays similar thermal conductivity. Using less GO in the hybrid nanoﬂuid could reduce the cost of the hybrid nanoﬂuid. Moreover, a low mass fraction of hybrid nanoparticles could enhance the thermal conductivity of ice. (5) The correlation generated for thermal conductivity of the hybrid nanoﬂuid provided an accurate prediction of the thermal conductivity of water-based GO–Al 2 O 3 hybrid nanoﬂuids. The new model for thermal conductivity of hybrid nanoﬂuid shows the inﬂuence of thickness, length, ﬂat rate and the interface thermal resistance of the graphene oxide, particle radius of Al 2 O 3 , ﬂuid viscosity, temperature and Brownian motion. Results presented in this study gave a better understanding of the inﬂuence of various important factors on the enhancement of thermal conductivity of hybrid nanoﬂuids.


Introduction
Thermal ice storage is a very important kind of phase change cold energy storage, which can be used in the solar and wind energy system to reduce the fluctuations in the energy flow. The ice storage system not only can run fluently under a decoupled source-load condition, but can reduce maximum output power, average power and energy consumption by 20-40% in comparison with conventional systems in centralized and decentralized environments [1]. The thermal ice storage is also employed in heating, ventilation and air conditioning (HVAC) cooling, food processing, chemical reactions or pharmaceutical processing, inlet air cooling of turbine and district cooling plant [2].
As one type of phase change material (PCM), ice storage, however, has higher supercooling degree and lower thermal conductivity that reduce heat transfer performance. A high supercooling degree leads to reduced performance of the PCM thermal energy storage system because of the requirement of a large operating temperature range [3]. Fortunately, many researchers found that addition of nanoparticles not only enhance the thermal conductivity, but also reduce the supercooling degree, and improve the heat transfer performance greatly [4]. Nonetheless, the stability of the nanoparticles suspended in base fluids during the thermal cycling is becoming a key problem.
Colloidal suspensions containing dispersed metal or nonmetal nanoparticles are called nanofluids, which are widely used in several heat transfer or performance improvement areas such as solar energy collectors [5][6][7][8][9][10][11][12][13], phase change energy storage [14,15], micro-and power electronics [16], magnetic refrigeration [17,18], diesel or biodiesel fuels [19,20], ventilation and air conditioning systems [21], reactors in power generating plants, cancer treatment therapy, metallurgical operations and combustion engines [22][23][24]. Nanofluids can be divided into two categories according to the types of nanoparticles. One is a single nanofluid formed by a single nanoparticle and base liquids, and the other is a hybrid nanofluid formed by mixing multiple nanoparticles and base liquids. Hybrid nanofluid is a solid-liquid mixture formed by dispersing two or more different nanoparticles suspended in a liquid, which displays enhanced thermophysical properties and rheological characteristics compared to mono nanoparticles nanofluid due to a synergistic effect [25]. Properly selected hybrid nanofluids can not only have higher thermal conductivity than single-phase nanofluids, but also have better properties such as less mechanical resistance, better stability and lower cost. These superior characteristics make them particularly attractive for applications in many heat transfer areas [22,23]. However, there are still several challenges that have impeded the development of this field as concluded in many review articles [23][24][25][26][27][28][29][30][31][32][33][34], such as disagreement in results of various researchers on the same types of nanofluids; preparation of long-term, stable and homogeneous nanofluids; a lack of theoretical understanding of the mechanisms accountable for variety in properties and high cost of the hybrid nanoparticles. In particular, stability of nanofluids is an essential characteristic that impacts greatly the thermophysical properties and rheological characteristics as well as their field applications. Surfactants have been used to keep good dispersibility of nanofluids. However, surfactants may impact the thermal conductivity and viscosity of the hybrid nanofluid as well [35]. Among various nanomaterials used in nanofluids, graphene oxide (GO) is hydrophilic, has a superior dispersion stability and high thermal conductivity and can significantly reduce the degree of supercooling of water in phase change thermal storage [36][37][38][39][40][41], but the preparation technology is complex and expensive. Nano-Al 2 O 3 is a common metal oxide material and is widely used in various nanofluids because the Al 2 O 3 nanofluid has lower viscosity and better rheological behavior as compared to other particles such as copper oxide [42][43][44]. Moreover, nano-Al 2 O 3 has obvious advantages such as high thermal conductivity, cheap and safe [45][46][47][48]. For example, nanofluids of aluminum oxide can enhance thermal conductivity from 0.3% to 38% with a particle size range from 5 to 80 nm [49][50][51]. Moreover, the supercooling degree of water is reduced under the combined effect of ultrasound and aluminum oxide nanoparticles. A 63.7% reduction in the supercooling degree of water can be obtained at Al 2 O 3 nanoparticle concentration of 0.2 wt % with the proper ultrasonic intensity [52]. However, long-term stability of aqueous Al 2 O 3 is not as good as that of graphene oxide.
To reduce the cost and make full use of the advantages of the two materials, our study prepares a novel aqueous Al 2 O 3 and GO nanosheet hybrid nanofluid for thermal ice storage as we have done in our previous work about its specific heat [53]. The thermal conductivity and stability under thermal cycles of the hybrid nanofluid are investigated, the thermal conductivity of the ice mixed with the hybrid nanoparticles is tested, and a new thermal conductivity model of the hybrid nanofluids considering the influence of Brownian and construction is proposed, which presents the theoretical understanding of the mechanisms responsible for the enhancement of thermal conductivity.

Experimental Apparatus and Procedure
The GO-Al 2 O 3 hybrid nanofluids were created by the two-step method [54] using sodium dodecyl sulfonate (SDS) as a surfactant. Different masses of the graphene oxide nanoplatelet (purity ≥99 wt %, number of layers <10, from Shenzhen Turing Evolution Technology Co., Ltd.) and nano-Al 2 O 3 (content ≥ 99.99%, average particle size (30 ± 5 nm)) were fully dispersed in deionized water (DW). High-powered ultrasonication was used to disperse nanoparticles. Sodium hydroxide solution was used to adjust the pH.
The hybrid nanomaterials were analyzed and characterized by an X-ray diffractometer (Model: XRD Ultima IV). WGZ-2000 turbidity meter (indication error of ±3% F.S and repeatability ≤0.5%, from Shanghai Xinrui Instrument Co., Ltd. China) was used to test the absorbance of nanofluids with time. It can measure the degree of scattering or attenuation of light generated by insoluble particulate matter suspended in water or a transparent liquid, and quantitatively characterize the content of these suspended particulate matter with NTU (Nephelometric Turbidity Unit) as the turbidity unit.
The transient hot wire method has been regarded as the most accurate method as well as the primary method for determination of thermal conductivity of a fluid and solid. The method has become a standard method in many countries [55][56][57]. Thermal conductivity is measured by a TC3000E transient hot wire instrument (Xiaxi Electronic Technology Co., Ltd. Xi'an, China), which has been employed widely by many researchers [58][59][60]. The instrument works according to transient hot-wire method that complies with standards of ASTM C1113, ASTM D5930, GB/T 10297, and GB/T 11205. A water bath with an accuracy of ± 0.05 • C was employed to control temperature.
X-Ray Diffraction (XRD) is very useful to show the crystal structure [61], the hybrid nanomaterials are characterized by XRD as indicated in Figure 1. It can be seen that the diffraction peak at 2θ = 11.06 • was the characteristic peak of graphene oxide (001) [ theoretical understanding of the mechanisms responsible for the enhancement of thermal conductivity.

Experimental Apparatus and Procedure
The GO-Al2O3 hybrid nanofluids were created by the two-step method [54] using sodium dodecyl sulfonate (SDS) as a surfactant. Different masses of the graphene oxide nanoplatelet (purity ≥99 wt %, number of layers <10, from Shenzhen Turing Evolution Technology Co., Ltd.) and nano-Al2O3 (content ≥ 99.99%, average particle size (30 ± 5 nm)) were fully dispersed in deionized water (DW). High-powered ultrasonication was used to disperse nanoparticles. Sodium hydroxide solution was used to adjust the pH.
The hybrid nanomaterials were analyzed and characterized by an X-ray diffractometer (Model: XRD Ultima IV). WGZ-2000 turbidity meter (indication error of ±3% F.S and repeatability ≤0.5%, from Shanghai Xinrui Instrument Co., Ltd. China) was used to test the absorbance of nanofluids with time. It can measure the degree of scattering or attenuation of light generated by insoluble particulate matter suspended in water or a transparent liquid, and quantitatively characterize the content of these suspended particulate matter with NTU (Nephelometric Turbidity Unit) as the turbidity unit.
The transient hot wire method has been regarded as the most accurate method as well as the primary method for determination of thermal conductivity of a fluid and solid. The method has become a standard method in many countries [55][56][57]. Thermal conductivity is measured by a TC3000E transient hot wire instrument (Xiaxi Electronic Technology Co., Ltd. Xi'an, China), which has been employed widely by many researchers [58][59][60]. The instrument works according to transient hot-wire method that complies with standards of ASTM C1113, ASTM D5930, GB/T 10297, and GB/T 11205. A water bath with an accuracy of ± 0.05 °C was employed to control temperature.

Effect of Different Factors on the Stability of the Hybrid Nanofluids
According to the theory of electrostatic stabilization and steric stabilization [64], nanofluids stability is affected greatly by many factors such as the pH, surfactants and ultrasonic power and time. In the following parts, the aforesaid factors are investigated.

Influence of pH
The pH of nanofluids is usually adjusted to increase the charges on the surface of the nanoparticles, which keeps the surface potential of the nanoparticles away from isoelectric point so as to enhance the electrostatic repulsion and stability of the nanofluids. The pH adjustment is suitable for electrostatic stabilization mechanism of colloidal particles that can be described by the DLVO

Effect of Different Factors on the Stability of the Hybrid Nanofluids
According to the theory of electrostatic stabilization and steric stabilization [64], nanofluids stability is affected greatly by many factors such as the pH, surfactants and ultrasonic power and time. In the following parts, the aforesaid factors are investigated.

Influence of pH
The pH of nanofluids is usually adjusted to increase the charges on the surface of the nanoparticles, which keeps the surface potential of the nanoparticles away from isoelectric point so as to enhance the electrostatic repulsion and stability of the nanofluids. The pH adjustment is suitable for electrostatic stabilization mechanism of colloidal particles that can be described by the DLVO

Effect of Different Factors on the Stability of the Hybrid Nanofluids
According to the theory of electrostatic stabilization and steric stabilization [64], nanofluids stability is affected greatly by many factors such as the pH, surfactants and ultrasonic power and time. In the following parts, the aforesaid factors are investigated. The pH of nanofluids is usually adjusted to increase the charges on the surface of the nanoparticles, which keeps the surface potential of the nanoparticles away from isoelectric point so as to enhance the electrostatic repulsion and stability of the nanofluids. The pH adjustment is suitable for electrostatic stabilization mechanism of colloidal particles that can be described by the DLVO (Derjaguin, Landau, Verwey and Overbeek) theory, i.e., the total interaction between electrostatic stabilized particles is the combination of van der Waals attraction and electrostatic repulsion [64].
Hybrid nanofluids of 0.1% mass fraction with pH = 3, 7, and 11 were prepared and ultrasonically dispersed for 60 min. Figures 4 and 5 show, respectively, the condition of the nanofluids after one and seven days. The hybrid nanofluid with pH = 3 has completely precipitated after 7 days, and the hybrid nanofluid with pH = 7 had a large amount of precipitation, while the hybrid nanofluid with pH = 11 had no obvious change and showed very good stability. The same trend is observed for bottom precipitation, as shown in Figure 5. The pH had a great influence on the stability of the hybrid nanofluid, and the hybrid nanofluid of pH = 11 was shown to have the best stability. This is mainly due to the fact that Al 2 O 3 nanoparticles hydrate to produce hydroxyl groups, and the hydroxyl groups dissociate to charge the surface of the nanoparticles [65]. There are many hydrophilic acidic functional groups on the surface of the graphene oxide nanoparticles, which have strong ability to adsorb anions. Under alkaline conditions, a large amount of OH − ions were continuously adsorbed on the surfaces of Al 2 O 3 and GO nanoparticles to form double electric layers around the nanoparticles, which produced repulsive forces between the particles and prevent agglomeration and precipitation caused by collision between particles. Moreover, larger electrostatic repulsion also increased the distance between particles, which decreased van der Waals attraction and, thus, reduced the probability of particles agglomerating and sedimentation, and further improved the stability of the nanofluids [66].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 23 (Derjaguin, Landau, Verwey and Overbeek) theory, i.e., the total interaction between electrostatic stabilized particles is the combination of van der Waals attraction and electrostatic repulsion [64]. Hybrid nanofluids of 0.1% mass fraction with pH = 3, 7, and 11 were prepared and ultrasonically dispersed for 60 min. Figures 4 and 5 show, respectively, the condition of the nanofluids after one and seven days. The hybrid nanofluid with pH = 3 has completely precipitated after 7 days, and the hybrid nanofluid with pH = 7 had a large amount of precipitation, while the hybrid nanofluid with pH = 11 had no obvious change and showed very good stability. The same trend is observed for bottom precipitation, as shown in Figure 5. The pH had a great influence on the stability of the hybrid nanofluid, and the hybrid nanofluid of pH = 11 was shown to have the best stability. This is mainly due to the fact that Al2O3 nanoparticles hydrate to produce hydroxyl groups, and the hydroxyl groups dissociate to charge the surface of the nanoparticles [65]. There are many hydrophilic acidic functional groups on the surface of the graphene oxide nanoparticles, which have strong ability to adsorb anions. Under alkaline conditions, a large amount of ions were continuously adsorbed on the surfaces of Al2O3 and GO nanoparticles to form double electric layers around the nanoparticles, which produced repulsive forces between the particles and prevent agglomeration and precipitation caused by collision between particles. Moreover, larger electrostatic repulsion also increased the distance between particles, which decreased van der Waals attraction and, thus, reduced the probability of particles agglomerating and sedimentation, and further improved the stability of the nanofluids [66].

Influence of Dispersant
One of the dispersing mechanisms for nanofluid stabilization is steric stabilization, which can be attained by the addition of a dispersant [64]. In order to enhance the stability of the nanofluid, polymer chains are attached to charged particle surfaces or polyelectrolyte chains are attached onto uncharged particle surfaces to form a steric barrier between nanoparticles to screen the interparticle van der Waals attraction [67,68]. Compared to electrostatic stabilization, steric stabilization has many advantages [64]: (1) It is widely used in various dispersion systems as a thermodynamic stabilization method; (2) It is applicable to high concentration of nanoparticles and it redisperses the temporarily aggregated particles; (3) The dispersion is not sensitive to categories of electrolyte; (4) It shows good compatibility for simultaneously dispersing hybrid particles within the same base fluid. SDS, a hydrophilic anionic surfactant with good penetration and dispersion, which has superior compatibility with anions and no-ions, was used in this study. The hybrid nanofluid of 0.1 wt % (mass ratio of GO and Al2O3 is 1:1), pH = 11, was ultrasonically dispersed for 90 min, and a different fraction of SDS was added. Figure 6 shows the state of the hybrid nanofluid.
It shows good compatibility for simultaneously dispersing hybrid particles within the same base fluid. SDS, a hydrophilic anionic surfactant with good penetration and dispersion, which has superior compatibility with anions and no-ions, was used in this study. The hybrid nanofluid of 0.1 wt % (mass ratio of GO and Al2O3 was 1:1), pH = 11, was ultrasonically dispersed for 90 min, and a different fraction of SDS was added. Figure 6 shows the state of the hybrid nanofluid samples with (left side) and without (right side) SDS after 3 and 30 days. It can be seen that after 30 days, the hybrid nanofluid without dispersant had more sediments and a faded color. In comparison, the hybrid nanofluid with dispersant SDS had neither obvious agglomeration nor color change. It is clearly shown that the SDS can keep the nanofluid stable for a relatively long time.
For quantitative analysis of the stability, the absorbance of the hybrid nanofluid was measured with time using the WGZ-2000 turbidity meter. At some intervals after preparation of the nanofluid, 4 mL of the upper layer of the nanofluid was taken and diluted three times with deionized water, and the absorbance was sequentially measured. The absorbance reduction ratio (ARR) was defined as: A0 and An are the initial absorbance and the nth hour absorbance of each nanofluid, respectively. Figure 7 indicates the influence of different concentration of SDS on ARR. It was found that the nanofluid's absorbance decreased rapidly (ARR increases) in the first four days, then decreased slowly after (slower ARR increase), as can be seen in Figure 7a. The stability of the hybrid nanofluid was improved greatly with the addition of proper amount of SDS. Moreover, the absorbance increased (ARR decreased) at first, then decreased (ARR increased) with the addition of the SDS

Influence of Dispersant
One of the dispersing mechanisms for nanofluid stabilization is steric stabilization, which can be attained by the addition of a dispersant [64]. In order to enhance the stability of the nanofluid, polymer chains are attached to charged particle surfaces or polyelectrolyte chains are attached onto uncharged particle surfaces to form a steric barrier between nanoparticles to screen the interparticle van der Waals attraction [67,68]. Compared to electrostatic stabilization, steric stabilization has many advantages [64]: (1) It is widely used in various dispersion systems as a thermodynamic stabilization method; (2) It is applicable to high concentration of nanoparticles and it redisperses the temporarily aggregated particles; (3) The dispersion is not sensitive to categories of electrolyte; (4) It shows good compatibility for simultaneously dispersing hybrid particles within the same base fluid. SDS, a hydrophilic anionic surfactant with good penetration and dispersion, which has superior compatibility with anions and no-ions, was used in this study. The hybrid nanofluid of 0.1 wt % (mass ratio of GO and Al 2 O 3 is 1:1), pH = 11, was ultrasonically dispersed for 90 min, and a different fraction of SDS was added. Figure 6 shows the state of the hybrid nanofluid.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 23 fraction of above 1% as shown in Figure 7b. This is due to SDS being attached on the surface of nanoparticles to form double electric layers that increase the distance between nanoparticles, which in turn decrease the van der Waals attraction potential. The particles can be completely shielded by anti-ions with a suitable increase of SDS, and the two electric layers are in electrostatic equilibrium by repulsive force, which keep good stability of the nanofluid [68,69]. However, if the SDS amount surpasses the critical micelles concentration (CMC), surfactant monomers can join up together to form micelles via the hydrophilic group instead of attaching on the surface of nanoparticle. These strongly charged micelles could repel the charged nanoparticles and prompt nanoparticles to agglomerate [70][71][72].
(a) (b)   It shows good compatibility for simultaneously dispersing hybrid particles within the same base fluid. SDS, a hydrophilic anionic surfactant with good penetration and dispersion, which has superior compatibility with anions and no-ions, was used in this study. The hybrid nanofluid of 0.1 wt % (mass Appl. Sci. 2020, 10, 5768 7 of 23 ratio of GO and Al 2 O 3 was 1:1), pH = 11, was ultrasonically dispersed for 90 min, and a different fraction of SDS was added. Figure 6 shows the state of the hybrid nanofluid samples with (left side) and without (right side) SDS after 3 and 30 days. It can be seen that after 30 days, the hybrid nanofluid without dispersant had more sediments and a faded color. In comparison, the hybrid nanofluid with dispersant SDS had neither obvious agglomeration nor color change. It is clearly shown that the SDS can keep the nanofluid stable for a relatively long time.
For quantitative analysis of the stability, the absorbance of the hybrid nanofluid was measured with time using the WGZ-2000 turbidity meter. At some intervals after preparation of the nanofluid, 4 mL of the upper layer of the nanofluid was taken and diluted three times with deionized water, and the absorbance was sequentially measured. The absorbance reduction ratio (ARR) was defined as: A 0 and A n are the initial absorbance and the nth hour absorbance of each nanofluid, respectively. Figure 7 indicates the influence of different concentration of SDS on ARR. It was found that the nanofluid's absorbance decreased rapidly (ARR increases) in the first four days, then decreased slowly after (slower ARR increase), as can be seen in Figure 7a. The stability of the hybrid nanofluid was improved greatly with the addition of proper amount of SDS. Moreover, the absorbance increased (ARR decreased) at first, then decreased (ARR increased) with the addition of the SDS fraction of above 1% as shown in Figure 7b. This is due to SDS being attached on the surface of nanoparticles to form double electric layers that increase the distance between nanoparticles, which in turn decrease the van der Waals attraction potential. The particles can be completely shielded by anti-ions with a suitable increase of SDS, and the two electric layers are in electrostatic equilibrium by repulsive force, which keep good stability of the nanofluid [68,69]. However, if the SDS amount surpasses the critical micelles concentration (CMC), surfactant monomers can join up together to form micelles via the hydrophilic group instead of attaching on the surface of nanoparticle. These strongly charged micelles could repel the charged nanoparticles and prompt nanoparticles to agglomerate [70][71][72].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 23 fraction of above 1% as shown in Figure 7b. This is due to SDS being attached on the surface of nanoparticles to form double electric layers that increase the distance between nanoparticles, which in turn decrease the van der Waals attraction potential. The particles can be completely shielded by anti-ions with a suitable increase of SDS, and the two electric layers are in electrostatic equilibrium by repulsive force, which keep good stability of the nanofluid [68,69]. However, if the SDS amount surpasses the critical micelles concentration (CMC), surfactant monomers can join up together to form micelles via the hydrophilic group instead of attaching on the surface of nanoparticle. These strongly charged micelles could repel the charged nanoparticles and prompt nanoparticles to agglomerate [70][71][72].

Influence of Ultrasonic Power
The absorbance was measured at different ultrasonic power for the hybrid nanofluid sample and graphene oxide aqueous nanofluid, at pH = 11, 0.1% nanoparticle concentration, 1% SDS and ultrasonic dispersion of 120 min. The relation between ARR and ultrasonic power is shown in Figures  8 and 9.

Influence of Ultrasonic Power
The absorbance was measured at different ultrasonic power for the hybrid nanofluid sample and graphene oxide aqueous nanofluid, at pH = 11, 0.1% nanoparticle concentration, 1% SDS and ultrasonic dispersion of 120 min. The relation between ARR and ultrasonic power is shown in Figures 8 and 9  Figure 8a shows rapid decrease in hybrid nanofluid absorbance in the first four days, indicated by the rapid increase of ARR, and relatively stable absorbance afterwards. Moreover, the absorbance firstly increased, then decreased as ultrasonic power increased, as can be inferred from Figure 8b.
The results indicate that the stability of the hybrid nanofluid reached optimum at 375 W of ultrasonic power, where the ARR was 17.16% lower than that at 450 W after 7 days. This is because the ultrasonic cavitation became stronger as ultrasonic power increased, which reduced the size of the nanoparticles and increased the spacing distance between the nanoparticles by the electrostatic and steric mechanism. The van der Waals attraction between nanoparticles is thus reduced, which hinders agglomeration and enhances the stability of the hybrid nanofluid. However, when the ultrasonic power is larger than the optimum value, it leads surfactant monomers to form micelles rather than attaching on the surface of nanoparticle. The intensity and frequency of nanoparticle collision increased greatly, and colliding nanoparticles generated larger agglomerates and precipitates.
For the GO nanofluid shown in Figure 9, the influence of ultrasonic power on the stability was consistent with the trend of the GO-Al2O3 hybrid nanofluids. Nevertheless, the ARR of the GO nanofluid was lower than that of the hybrid nanofluid, which means higher absorbance and better stability of GO nanofluid as compared to the stability of the hybrid nanofluid, which was slightly reduced with the addition of the alumina nanoparticles.    Figure 8a shows rapid decrease in hybrid nanofluid absorbance in the first four days, indicated by the rapid increase of ARR, and relatively stable absorbance afterwards. Moreover, the absorbance firstly increased, then decreased as ultrasonic power increased, as can be inferred from Figure 8b. The results indicate that the stability of the hybrid nanofluid reached optimum at 375 W of ultrasonic power, where the ARR was 17.16% lower than that at 450 W after 7 days. This is because the ultrasonic cavitation became stronger as ultrasonic power increased, which reduced the size of the nanoparticles and increased the spacing distance between the nanoparticles by the electrostatic and steric mechanism. The van der Waals attraction between nanoparticles is thus reduced, which hinders agglomeration and enhances the stability of the hybrid nanofluid. However, when the ultrasonic power is larger than the optimum value, it leads surfactant monomers to form micelles rather than attaching on the surface of nanoparticle. The intensity and frequency of nanoparticle collision increased greatly, and colliding nanoparticles generated larger agglomerates and precipitates.

Influence of Ultrasonic Time
For the GO nanofluid shown in Figure 9, the influence of ultrasonic power on the stability was consistent with the trend of the GO-Al2O3 hybrid nanofluids. Nevertheless, the ARR of the GO nanofluid was lower than that of the hybrid nanofluid, which means higher absorbance and better stability of GO nanofluid as compared to the stability of the hybrid nanofluid, which was slightly reduced with the addition of the alumina nanoparticles.  Figure 8a shows rapid decrease in hybrid nanofluid absorbance in the first four days, indicated by the rapid increase of ARR, and relatively stable absorbance afterwards. Moreover, the absorbance firstly increased, then decreased as ultrasonic power increased, as can be inferred from Figure 8b. The results indicate that the stability of the hybrid nanofluid reached optimum at 375 W of ultrasonic power, where the ARR was 17.16% lower than that at 450 W after 7 days. This is because the ultrasonic cavitation became stronger as ultrasonic power increased, which reduced the size of the nanoparticles and increased the spacing distance between the nanoparticles by the electrostatic and steric mechanism. The van der Waals attraction between nanoparticles is thus reduced, which hinders agglomeration and enhances the stability of the hybrid nanofluid. However, when the ultrasonic power is larger than the optimum value, it leads surfactant monomers to form micelles rather than attaching on the surface of nanoparticle. The intensity and frequency of nanoparticle collision increased greatly, and colliding nanoparticles generated larger agglomerates and precipitates.

Influence of Ultrasonic Time
For the GO nanofluid shown in Figure 9, the influence of ultrasonic power on the stability was consistent with the trend of the GO-Al 2 O 3 hybrid nanofluids. Nevertheless, the ARR of the GO nanofluid was lower than that of the hybrid nanofluid, which means higher absorbance and better stability of GO nanofluid as compared to the stability of the hybrid nanofluid, which was slightly reduced with the addition of the alumina nanoparticles.

Influence of Ultrasonic Time
The absorbance was studied at different ultrasonic time for the hybrid nanofluid sample and graphene oxide aqueous nanofluid, at pH = 11, 0.1% nanoparticle concentration, 1% SDS and ultrasonic power of 375 W. The relationship between ARR and ultrasonic time is shown in Figures 10 and 11. agglomerates when ultrasonic time is short. As the ultrasonic time increases, there is enough time to smash agglomerates into smaller nanoparticles, and the nanoparticles can maintain suitable spacing under the coaction of dispersants and ions in the nanofluid, so that the van der Waals force keeps balance with electrostatic repulsion. However, when the ultrasonic time is longer than the optimal value, similar to the influence of ultrasonic power, it produces surfactant micelles, and simultaneously increases the frequency of nanoparticle collision greatly, leading to agglomeration and precipitation. In comparison with ARR of the GO nanofluid shown in Figure 11, the ultrasonic time had a different influence on the graphene oxide nanofluid; the longer ultrasonic time the better the stability of the nanofluid, because the GO nanosheet was hydrophilic and had superior dispersion stability in water, so even longer ultrasonic time could not form large agglomerates and sediment.  smash agglomerates into smaller nanoparticles, and the nanoparticles can maintain suitable spacing under the coaction of dispersants and ions in the nanofluid, so that the van der Waals force keeps balance with electrostatic repulsion. However, when the ultrasonic time is longer than the optimal value, similar to the influence of ultrasonic power, it produces surfactant micelles, and simultaneously increases the frequency of nanoparticle collision greatly, leading to agglomeration and precipitation. In comparison with ARR of the GO nanofluid shown in Figure 11, the ultrasonic time had a different influence on the graphene oxide nanofluid; the longer ultrasonic time the better the stability of the nanofluid, because the GO nanosheet was hydrophilic and had superior dispersion stability in water, so even longer ultrasonic time could not form large agglomerates and sediment.    It can be seen from Figure 10 that the trend of change in ARR with holding time and ultrasonic time was similar to that of ultrasonic power; a decrease in absorbance (ARR increases) in the first 4 days, and stable afterwards. From Figure 10b, it can be inferred that the stability of the hybrid nanofluid was optimal at 120 min of ultrasonic time, where the ARR was 29.65% lower than that at 30 min after 7 days. This is due to the fact that the ultrasonic cannot completely break the agglomerates when ultrasonic time is short. As the ultrasonic time increases, there is enough time to smash agglomerates into smaller nanoparticles, and the nanoparticles can maintain suitable spacing under the coaction of dispersants and ions in the nanofluid, so that the van der Waals force keeps balance with electrostatic repulsion. However, when the ultrasonic time is longer than the optimal value, similar to the influence of ultrasonic power, it produces surfactant micelles, and simultaneously increases the frequency of nanoparticle collision greatly, leading to agglomeration and precipitation. In comparison with ARR of the GO nanofluid shown in Figure 11, the ultrasonic time had a different influence on the graphene oxide nanofluid; the longer ultrasonic time the better the stability of the nanofluid, because the GO nanosheet was hydrophilic and had superior dispersion stability in water, so even longer ultrasonic time could not form large agglomerates and sediment.

Influence of Thermal Cycling
The absorbance was measured, as shown in Figure 12, with different numbers of thermal cycle of solidification and melting for the hybrid nanofluid sample at pH = 11, 0.1% nanoparticle concentration, 1% SDS, ultrasonic power 300 W and ultrasonic time of 120 min. Figure 11. Relationship between the absorbance reduction rate of the GO nanofluid and time.

Influence of Thermal Cycling
The absorbance was measured, as shown in Figure 12, with different numbers of thermal cycle of solidification and melting for the hybrid nanofluid sample at pH = 11, 0.1% nanoparticle concentration, 1% SDS, ultrasonic power 300 W and ultrasonic time of 120 min.
The absorbance of hybrid nanofluids decreased slightly before the 6 th cycle, and stayed stable afterwards, with an overall decrease of 6.2% at the 10 th cycle as compared with the initial absorbance. This is mainly because some aggregation and deposition occurred as the balance between the van der Waals attraction and electrostatic repulsion was broken by crystallization and melting. Nevertheless, the decrease in stability was slight, owing to the screen function of dispersants surrounding the nanoparticles.

Reliability Verification
In order to verify the accuracy of the experimental system, the thermal conductivity of deionized water was measured from 0 to 70℃. Experimental results shown in Figure 13 are compared with reference values [73]. The maximum error was 1.12%, which was used as experimental uncertainty shown in the corresponding figures. The enhancement in thermal conductivity and the uncertainty of thermal conductivity enhancement of the hybrid nanofluid are defined as in our previous work [54] by Equations (2) and (3) respectively: where knf and kbf are the thermal conductivity of the nanofluid and base fluid, respectively. The uncertainty of the thermal conductivity enhancement of the nanofluids was calculated by where and are uncertainty of knf and kbf, respectively. The absorbance of hybrid nanofluids decreased slightly before the 6th cycle, and stayed stable afterwards, with an overall decrease of 6.2% at the 10th cycle as compared with the initial absorbance. This is mainly because some aggregation and deposition occurred as the balance between the van der Waals attraction and electrostatic repulsion was broken by crystallization and melting. Nevertheless, the decrease in stability was slight, owing to the screen function of dispersants surrounding the nanoparticles.

Reliability Verification
In order to verify the accuracy of the experimental system, the thermal conductivity of deionized water was measured from 0 to 70 • C. Experimental results shown in Figure 13 are compared with reference values [73]. The maximum error was 1.12%, which was used as experimental uncertainty shown in the corresponding figures. The enhancement in thermal conductivity and the uncertainty of thermal conductivity enhancement of the hybrid nanofluid are defined as in our previous work [54] by Equations (2) and (3) respectively: where k nf and k bf are the thermal conductivity of the nanofluid and base fluid, respectively. The uncertainty of the thermal conductivity enhancement of the nanofluids was calculated by where σ n f and σ b f are uncertainty of k nf and k bf , respectively. Appl. Sci. 2020, 10, x FOR PEER REVIEW 11 of 23  Figure 14 shows the influence of mass fraction of nanoparticle on thermal conductivity of the hybrid nanofluid. It indicates that the thermal conductivity of the hybrid nanofluid increased with the mass fraction of the nanoparticles in the range of 0-0.2%. Moreover, the influence of the mass fraction on thermal conductivity enhancement at a lower temperature was more significant than at a high temperature. At 0 °C, the thermal conductivity of 0.2 wt % hybrid nanofluids increased by 7.2% as compared to the base fluid, and the thermal conductivity enhancement of 0.2 wt % was about 2.1 times larger than that of 0.025 wt %. At a higher temperature, 70 °C, the thermal conductivity of 0.2 wt % hybrid nanofluid increased by 1.76% as compared to the base fluid, and the thermal conductivity enhancement was 14.7 times that of the 0.025 wt % hybrid nanofluid. It shows that the nanoparticle mass fraction had a significant influence on the thermal conductivity of the hybrid nanofluid, owing to the fact that thermal conductivities of GO and Al2O3 nanoparticles were larger than that of deionized water. The liquid molecules on the surface of the nanoparticles form ordered structure, which acts as a thermal bridge between the fluid and the solid particles, resulted in a higher thermal conductivity and enhanced energy transfer within the system. Meanwhile, Brownian motion of nanoparticles caused microconvection in the liquid, which enhanced the energy transfer.

Effect of Temperature
Thermal conductivity of the hybrid nanofluid increased with temperature, as displayed in Figure 15, which was consistent with the trend of thermal conductivity of deionized water with temperature. Moreover, the thermal conductivity enhancement of the hybrid nanofluid decreased with increasing temperature. The thermal conductivity enhancement at 0 °C was 28.5 times that at 70  Figure 14 shows the influence of mass fraction of nanoparticle on thermal conductivity of the hybrid nanofluid. It indicates that the thermal conductivity of the hybrid nanofluid increased with the mass fraction of the nanoparticles in the range of 0-0.2%. Moreover, the influence of the mass fraction on thermal conductivity enhancement at a lower temperature was more significant than at a high temperature. At 0 • C, the thermal conductivity of 0.2 wt % hybrid nanofluids increased by 7.2% as compared to the base fluid, and the thermal conductivity enhancement of 0.2 wt % was about 2.1 times larger than that of 0.025 wt %. At a higher temperature, 70 • C, the thermal conductivity of 0.2 wt % hybrid nanofluid increased by 1.76% as compared to the base fluid, and the thermal conductivity enhancement was 14.7 times that of the 0.025 wt % hybrid nanofluid. It shows that the nanoparticle mass fraction had a significant influence on the thermal conductivity of the hybrid nanofluid, owing to the fact that thermal conductivities of GO and Al 2 O 3 nanoparticles were larger than that of deionized water. The liquid molecules on the surface of the nanoparticles form ordered structure, which acts as a thermal bridge between the fluid and the solid particles, resulted in a higher thermal conductivity and enhanced energy transfer within the system. Meanwhile, Brownian motion of nanoparticles caused microconvection in the liquid, which enhanced the energy transfer.

Effect of Temperature
Thermal conductivity of the hybrid nanofluid increased with temperature, as displayed in Figure 15, which was consistent with the trend of thermal conductivity of deionized water with temperature. Moreover, the thermal conductivity enhancement of the hybrid nanofluid decreased with increasing temperature. The thermal conductivity enhancement at 0 • C was 28.5 times that at 70 • C for 0.025 wt %, while the thermal conductivity enhancement at 0 • C was 4 times that at 70 • C for 0.2 wt %. This means that the influence of temperature on the thermal conductivity enhancement was more significant at a lower concentration, even though a higher concentration led to a higher thermal conductivity enhancement. This can be attributed to the fact that Brownian motion and thermal diffusion of nanoparticles increased as the temperature increased.
nanoparticle mass fraction had a significant influence on the thermal conductivity of the hybrid nanofluid, owing to the fact that thermal conductivities of GO and Al2O3 nanoparticles were larger than that of deionized water. The liquid molecules on the surface of the nanoparticles form ordered structure, which acts as a thermal bridge between the fluid and the solid particles, resulted in a higher thermal conductivity and enhanced energy transfer within the system. Meanwhile, Brownian motion of nanoparticles caused microconvection in the liquid, which enhanced the energy transfer.

Effect of Temperature
Thermal conductivity of the hybrid nanofluid increased with temperature, as displayed in Figure 15, which was consistent with the trend of thermal conductivity of deionized water with temperature. Moreover, the thermal conductivity enhancement of the hybrid nanofluid decreased with increasing temperature. The thermal conductivity enhancement at 0 °C was 28.5 times that at 70 °C for 0.025 wt %, while the thermal conductivity enhancement at 0 °C was 4 times that at 70 °C for 0.2 wt %. This means that the influence of temperature on the thermal conductivity enhancement was more significant at a lower concentration, even though a higher concentration led to a higher thermal conductivity enhancement. This can be attributed to the fact that Brownian motion and thermal diffusion of nanoparticles increased as the temperature increased.

Comparison of Thermal Conductivity of GO Nanofluid with the Hybrid Nanofluid
The thermal conductivity of aqueous GO nanofluid is measured at the same mass fraction and temperature as the GO-Al2O3 hybrid nanofluid, as shown in Figure 16. It displays a similar trend and thermal conductivity values, suggesting the potential benefit for using less GO and reducing the cost of the hybrid nanofluid. There are many correlations for the thermal conductivity of hybrid nanofluids considering the influence of temperature and volume fraction of nanoparticles [23].
Mohammad et al. [74] proposed a formula as Equation (4) for thermal conductivity of DWCNT- Figure 15. Influence of temperature on thermal conductivity.

Comparison of Thermal Conductivity of GO Nanofluid with the Hybrid Nanofluid
The thermal conductivity of aqueous GO nanofluid is measured at the same mass fraction and temperature as the GO-Al 2 O 3 hybrid nanofluid, as shown in Figure 16. It displays a similar trend and thermal conductivity values, suggesting the potential benefit for using less GO and reducing the cost of the hybrid nanofluid.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 12 of 23 °C for 0.025 wt %, while the thermal conductivity enhancement at 0 °C was 4 times that at 70 °C for 0.2 wt %. This means that the influence of temperature on the thermal conductivity enhancement was more significant at a lower concentration, even though a higher concentration led to a higher thermal conductivity enhancement. This can be attributed to the fact that Brownian motion and thermal diffusion of nanoparticles increased as the temperature increased.

Comparison of Thermal Conductivity of GO Nanofluid with the Hybrid Nanofluid
The thermal conductivity of aqueous GO nanofluid is measured at the same mass fraction and temperature as the GO-Al2O3 hybrid nanofluid, as shown in Figure 16. It displays a similar trend and thermal conductivity values, suggesting the potential benefit for using less GO and reducing the cost of the hybrid nanofluid.

Fitting of Thermal Conductivity of the Hybrid Nanofluid
There are many correlations for the thermal conductivity of hybrid nanofluids considering the influence of temperature and volume fraction of nanoparticles [23].

Fitting of Thermal Conductivity of the Hybrid Nanofluid
There are many correlations for the thermal conductivity of hybrid nanofluids considering the influence of temperature and volume fraction of nanoparticles [23].
Equation (9) is the Maxwell model [79], which was developed based on the thermal conductivity of solid suspensions in liquids. It is a base correlation to predict the thermal conductivity of nanofluids [5].
The Takebi model [80] provided correlations for the thermal conductivity of hybrid nanofluid shown in Equation (11), which should have extensive applicability because it is based on the principle of Maxwell model and mixture rule that Then The comparison between experimental data and calculated values using the Takebi model is shown in Figure 17. The values of properties of the nanomaterials used in Equation (11) are listed in Table 1. The intrinsic thermal conductivity of GO was different for various fabrication processes and measurement methods [81], but it had very little influence on the fitting results even for 1000-3000 W/(m·K) using Equation (11) due to the low concentration of nanoparticles. Figure 17 shows that the relative error was larger at a lower temperature, with a maximum relative error of 6.5%.
Appl. Sci. 2020, 10 Regression Equation (12) was obtained considering the influence of temperature and mass fraction.
where ∅ is the mass fraction, from 0.025% to 0.2%; T is the temperature in ℃, ranging from 0 to 70 ℃ and a, b, c, d, e and g are the parameters fitted by experiment data, as shown in Table 2. The relationship between volume fraction and mass fraction of the particles can be obtained according to the following equation: where and are the density of the base fluid and density of the nanoparticle, respectively.  Figure 18 displays the fitted values calculated by the regression equation (12) and experimental data. The maximum relative error was 0.5% and R-square was 0.996. Compared with the Takebi model [80] for thermal conductivity of hybrid nanofluids, equation (12) can provide a more accurate prediction for the thermal conductivity of water-based GO-Al2O3 hybrid nanofluids, but it cannot reflect any mechanisms that explain the reason behind thermal conductivity enhancement.    [82] 3000 [81] 790 [53] Regression Equation (12) was obtained considering the influence of temperature and mass fraction.
where ∅ m is the mass fraction, from 0.025% to 0.2%; T is the temperature in • C, ranging from 0 to 70 • C and a, b, c, d, e and g are the parameters fitted by experiment data, as shown in Table 2.
The relationship between volume fraction and mass fraction of the particles can be obtained according to the following equation: where ρ b f and ρ np are the density of the base fluid and density of the nanoparticle, respectively.  Figure 18 displays the fitted values calculated by the regression Equation (12) and experimental data. The maximum relative error was 0.5% and R-square was 0.996. Compared with the Takebi model [80] for thermal conductivity of hybrid nanofluids, Equation (12) can provide a more accurate prediction for the thermal conductivity of water-based GO-Al 2 O 3 hybrid nanofluids, but it cannot reflect any mechanisms that explain the reason behind thermal conductivity enhancement. 10 − Figure 18 displays the fitted values calculated by the regression equation (12) and experimental data. The maximum relative error was 0.5% and R-square was 0.996. Compared with the Takebi model [80] for thermal conductivity of hybrid nanofluids, equation (12) can provide a more accurate prediction for the thermal conductivity of water-based GO-Al2O3 hybrid nanofluids, but it cannot reflect any mechanisms that explain the reason behind thermal conductivity enhancement. Mass fraction/% Thermal conductivity/W/(m.K) Figure 18. Comparison of fitted values and experimental data. Figure 18. Comparison of fitted values and experimental data.

Mechanism Model of Thermal Conductivity of the Hybrid Nanofluid
Many researchers found [54,83] that suspended nanoparticles enhance the energy transfer of nanofluids in two ways: the first is static effective thermal conductivity, which can be modeled according to the effective medium theory [84]; the second is the stochastic motions of the nanoparticles and the interfacial interactions between particles and the liquid molecules. Assuming that the two mechanisms can be separated, the total effective heat transfer in the nanofluidic suspension can be written as follows [85]: Al 2 O 3 is a spherical particle, thus, its static thermal conductivity can be established by the H-C model [86], which introduces the influence of particle shape on thermal conductivity. The model is as follows: k e f f ,EMT,1 where n is an empirical shape factor, and n is 3 for spherical particles. K e f f ,EMT,1 is the effective thermal conductivity of the Al 2 O 3 nanofluid. Considering the effect of thickness t, length L, interfacial thermal resistance R k and average flatness ratio η, the effective thermal conductivity of the GO nanofluid, K e f f ,EMT,1 , was calculated as Equation (16) [54,87]: For the enhancement of thermal conductivity by random motion of nanoparticles, especially in lower viscosity liquid, Equation (17) was employed by Xuan et al. [88] considering the effects of viscosity, temperature and properties of nanoparticles on thermal conductivity.
ρ np and C p,np are the density and specific thermal capacity of nanoparticles, respectively, µ is the viscosity of the base fluid, k B is the Boltzmann constant, T is the Kelvin temperature and r c is the mean radius of gyration of the cluster. For spherical Al 2 O 3 nanoparticles, r c is the radius. For the GO nanosheet, as in our previous study [54], r c is approximated to half of the equivalent diameter d p,eq of the nanosheet.
Non-spherical particles volume v non-sph can be deemed as theoretical spherical particles that have the same volumes as the non-spherical particles [89].
Combining the influence of static thermal conductivity and Brownian motion of the hybrid nanofluid, a new thermal conductivity model of the Al 2 O 3 -GO hybrid nanofluid, k eff,hnf , was developed as Equation (19): where ϕ is the volume fraction of Al 2 O 3 nanoparticle in the total hybrid nanoparticles. Subscript 1 and 2 are for Al 2 O 3 and GO, respectively. The exact values for L, t, η and R k are fitted based on the experimental data. According to our previous research work [54], the range of L, t, η and R k can be determined as: L ∈ (8 × 10 −6 , 6 × 10 −5 ), t ∈ (1.6 × 10 −9 , 6 × 10 −8 ), η ∈ (0.2, 0.8) and R k ∈ (1 × 10 −10 , 1 × 10 −6 ). The fitting results are shown in Table 3 when R-square was at the maximum with 95% confidence bounds. Moreover, L was approximately 20 µm according to SEM, t was about 5 nm as reported by the supplier. The values of the four parameters could be determined as L = 2.167 × 10 −5 m, t = 6.312 × 10 −9 m, η = 0.3211 and R k = 3.113 × 10 −9 m 2 Kw −1 . It is demonstrated that R k of the GO-Al 2 O 3 aqueous hybrid nanofluid was less than 1.5 × 10 −8 m 2 kw −1 of graphene and water interfacial thermal resistance [54]. This is because GO had better hydrophilicity and better wettability at the solid-liquid interface, which led to a lower interfacial thermal resistance. On the other hand, the interaction between the GO nanosheet and Al 2 O 3 nanoparticles might change the properties of the interface and reduce the interfacial thermal resistance. Figure 19 compares the calculated results using the new model with the experimental data; the average error was 0.18%, and the maximum deviation was 3.8% at 0 • C. The main reason is that the model underestimates the influence of Brownian motion at a lower temperature, as many researchers [83] found the temperature dependence of this model (∝ √ T) is still small and presents little agreement with several experimental studies, which means that the impact of Brownian motion is greater than that obtained in this work. It is demonstrated that Rk of the GO-Al2O3 aqueous hybrid nanofluid was less than 1.5 × 10 −8 m 2 kw −1 of graphene and water interfacial thermal resistance [54]. This is because GO had better hydrophilicity and better wettability at the solid-liquid interface, which led to a lower interfacial thermal resistance. On the other hand, the interaction between the GO nanosheet and Al2O3 nanoparticles might change the properties of the interface and reduce the interfacial thermal resistance. Figure 19 compares the calculated results using the new model with the experimental data; the average error was 0.18%, and the maximum deviation was 3.8% at 0 ℃. The main reason is that the model underestimates the influence of Brownian motion at a lower temperature, as many researchers [83] found the temperature dependence of this model (∝ √ ) is still small and presents little agreement with several experimental studies, which means that the impact of Brownian motion is greater than that obtained in this work.

Thermal Conductivity of Ice Containing the Hybrid Nanoparticles
The thermal conductivity of ice mixed with the hybrid nanoparticles may be different from that of pure ice. The thermal conductivity enhancement of the ice as functions of temperature and the mass fraction of hybrid nanoparticles is shown in Figure 20, and the standard thermal conductivity

Thermal Conductivity of Ice Containing the Hybrid Nanoparticles
The thermal conductivity of ice mixed with the hybrid nanoparticles may be different from that of pure ice. The thermal conductivity enhancement of the ice as functions of temperature and the mass fraction of hybrid nanoparticles is shown in Figure 20, and the standard thermal conductivity of ice refers to [90]. The thermal conductivity enhancement of ice follows a different behavior in comparison with that of the hybrid nanofluids. It seems that the low mass fraction of hybrid nanoparticles can enhance the thermal conductivity with an enhancement increase of about 8% at 0.025%, −10 • C. This is probably because the addition of hybrid nanoparticles changes the collective vibrations of atoms in the crystal lattice, and affects the transfer of phonons. However, only three plots were achieved at each temperature because of the influence of coronavirus, further research should be carried out to give more definite relations between mass fraction and the thermal conductivity enhancement.

Conclusions
Aqueous GO-Al2O3 hybrid nanofluids used in thermal ice storage were prepared and characterized. The influences of pH, dispersant, ultrasonic power and ultrasonic time on the stability of hybrid nanofluids were quantitatively analyzed. The thermal conductivities of the hybrid nanofluids were measured and correlated with temperature and concentration. Thermal conductivity of ice containing the hybrid nanoparticles was tested. A new model for the effective thermal conductivity of the hybrid nanofluids considering the structure and Brownian motion was proposed.
The key results of this study were as follows: (1) The optimal conditions for stability of the aqueous GO-Al2O3 hybrid nanofluids (mass fraction of GO: Al2O3 = 1:1) were: pH = 11, 1% SDS, ultrasonic power of 375 W and ultrasonication time of 120 min.
(2) With the increase in temperature and the mass fraction of the nanoparticles, thermal conductivity of the hybrid nanofluid increased. The influence of mass fraction on the thermal conductivity enhancement at a lower temperature was greater than that at a high temperature. At 0 °C, the thermal conductivity of 0.2 wt % hybrid nanofluids increased by 7.1% as compared to that of the base fluid, which means that this hybrid nanofluid had a potential advantage to enhance heat transfer performance.
(3) The nanoparticle mass fraction had a significant influence on the thermal conductivity of the hybrid nanofluid. The influence of temperature on thermal conductivity enhancement at a lower concentration was larger, even though at a higher concentration there was greater thermal conductivity enhancement.

Conclusions
Aqueous GO-Al 2 O 3 hybrid nanofluids used in thermal ice storage were prepared and characterized. The influences of pH, dispersant, ultrasonic power and ultrasonic time on the stability of hybrid nanofluids were quantitatively analyzed. The thermal conductivities of the hybrid nanofluids were measured and correlated with temperature and concentration. Thermal conductivity of ice containing the hybrid nanoparticles was tested. A new model for the effective thermal conductivity of the hybrid nanofluids considering the structure and Brownian motion was proposed. The key results of this study were as follows: (1) The optimal conditions for stability of the aqueous GO-Al 2 O 3 hybrid nanofluids (mass fraction of GO: Al 2 O 3 = 1:1) were: pH = 11, 1% SDS, ultrasonic power of 375 W and ultrasonication time of 120 min.
(2) With the increase in temperature and the mass fraction of the nanoparticles, thermal conductivity of the hybrid nanofluid increased. The influence of mass fraction on the thermal conductivity enhancement at a lower temperature was greater than that at a high temperature. At 0 • C, the thermal conductivity of 0.2 wt % hybrid nanofluids increased by 7.1% as compared to that of the base fluid, which means that this hybrid nanofluid had a potential advantage to enhance heat transfer performance.
(3) The nanoparticle mass fraction had a significant influence on the thermal conductivity of the hybrid nanofluid. The influence of temperature on thermal conductivity enhancement at a lower concentration was larger, even though at a higher concentration there was greater thermal conductivity enhancement.
(4) Comparison of thermal conductivity of the GO nanofluid with the hybrid nanofluid displays similar thermal conductivity. Using less GO in the hybrid nanofluid could reduce the cost of the hybrid nanofluid. Moreover, a low mass fraction of hybrid nanoparticles could enhance the thermal conductivity of ice.
(5) The correlation generated for thermal conductivity of the hybrid nanofluid provided an accurate prediction of the thermal conductivity of water-based GO-Al 2 O 3 hybrid nanofluids. The new model for thermal conductivity of hybrid nanofluid shows the influence of thickness, length, flat rate and the interface thermal resistance of the graphene oxide, particle radius of Al 2 O 3 , fluid viscosity, temperature and Brownian motion. Results presented in this study gave a better understanding of the influence of various important factors on the enhancement of thermal conductivity of hybrid nanofluids.