Experimental Study on Heat Transfer and Adsorption Cooling Performance of MIL-101/Few Layer Graphene Composite

: MIL-101 is a promising metal-organic frameworks (MOFs) material in adsorption chiller application due to its high adsorption capacity for water and excellent adsorption/desorption cyclic stability. Few layer graphene (FLG) as the thermal conductive additive was added into MIL-101 to improve inferior heat transfer of MIL-101 in the adsorption cooling process. The heat transfer characteristic of MIL-101/FLG adsorber and the adsorption cooling performance of the MIL-101/FLG-water working pair were studied. Results show that thermal conductivity of MIL-101/20%FLG composite is 5.79-6.54 times that of MIL-101. Adding FLG is conducive to the formation of heat transfer channels in MIL-101/FLG adsorber and the rapid removal of adsorption heat. The heating and cooling rate of MIL-101/FLG adsorber is ~2.2 times that of MIL-101 adsorber. Under typical adsorption water chiller conditions, the speciﬁc cooling power (SCP) and coefﬁcient of performance (COP) of the MIL-101/FLG-water working pair is 72.2–81.0 W kg − 1 and 0.187–0.202, respectively, at desorption temperatures of 70 ◦ C and 90 ◦ C, which is 1.43–1.56 times higher than the MIL-101-water working pair. The excellent structural and adsorption/desorption cyclic stability of MIL-101/FLG composite is veriﬁed after 50 consecutive cycles. It can provide a promising adsorbent candidate (MIL-101/FLG composite) in adsorption water chiller process.

Currently, some possible application directions of adsorption water chiller (16~20 • C) include proton exchange membrane fuel cell [27], automotive engine coolant [28], vehicle air-conditioning [29] and data center cooling [30,31]. The heat generated in fuel cell stacks (temperatures of 70-95 • C), internal combustion engines (100-130 • C) and server chipsets (temperatures of 70-90 • C) all could be used to drive the adsorption chiller. Although the In summary, the heat transfer property of silica gel and a 13X molecular sieve as well as the adsorption cooling performance of these composite materials/water working pairs could be enhanced by adding thermal conductive materials (metal, metal foam, multi-walled carbon nanotubes, graphite, etc.), which has been reported. Graphene is also a promising thermal conductive carbon material [48][49][50]. The thermal conductivity of graphene [51,52] and FLG [53,54] can reach as high as 2000-5000 W m −1 K −1 and 500-2000 W m −1 K −1 , respectively. Most of the literature on the adsorption cooling performance of MOFs has been based on simulative calculations of adsorption isotherm and adsorption kinetics, and there are few relevant reports on experimental studies. So far, adding few layer graphene (FLG) into MIL-101 (Cr) (herein MIL-101) for the enhancement of the heat transfer characteristics of an MIL-101/FLG adsorber and the adsorption cooling performance of an MIL-101/FLG composite-water working pair have not been reported. Moreover, there is also no report on an experimental adsorption cooling process of MOF composite-water working pairs in an adsorption water chiller application.
In this paper, FLG was added into MIL-101 for the enhancement of thermal conduction. The thermal conduction property of MIL-101/FLG composite adsorbents was determined. Heat transfer characteristics of MIL-101/FLG adsorbers during the adsorption cooling process were investigated. The heat transfer mechanism of MIL-101/FLG was analyzed. The adsorption cooling performance of MIL-101/FLG-water working pairs was investigated and compared to MIL-101-water working pairs under specific working conditions of the adsorption water chiller process. Moreover, the adsorption cooling cyclic stability of MIL-101/FLG-water working pairs was determined. The heat transfer characteristics and adsorption cooling performance of MIL-101/FLG-water working pairs were improved by incorporation of thermal conductive FLG, which could be conducive to increasing the energy conversion efficiency in adsorption water chiller applications. The results of this paper provide basic research for the utilization of an MIL-101/FLG composite in adsorption water chiller processes.

Preparation of MIL-101/FLG Composites
MIL-101 was prepared by a hydrothermal method [37,55], while the preparation method of MIL-101/FLG composites was described in our previous work [56].

Characterization Methods of MIL-101/FLG Composites
The surface micro morphology of FLG was conducted by scanning electron microscopy (SEM, Japan Hitachi, S4800). The crystal defects of FLG were measured by Raman spectra (Japan Horiba, Labram HR800, shift of 1000-3000 cm −1 ). The crystal structures of FLG, MIL-101/20%FLG before and after the adsorption cooling cycles were analyzed by Xray diffraction (XRD, Japan Rikagu, SmartLab) with the diffraction angles ranging from 2-30 • . N 2 adsorption/desorption isotherms of MIL-101/20%FLG before and after the adsorption cooling cycles at 77 K were measured on the BELSORP-max (Japan microtrac bayer) adsorption analyzer. The surface area was calculated by BET equation and the pore volume was calculated by t-plot methods.

Measurement of Thermophysical Properties of MIL-101/FLG Composites
As shown in Figure 1a   Thermophysical properties of MIL-101 and MIL-101/FLG composite were measured at 30-130 °C on a TPS-2500S type thermal property instrument (Hot-Disk Co., Ltd., Uppsala, Sweden). The thermal conductivity and specific heat capacity were measured through the "ISO22007-2-2015" standard [57,58]. The theoretical calculation equations of thermophysical properties refer to our previous work [56]. The standard deviations (error bars) of experimental results are shown in the results and discussion.

Determination of Adsorption and Desorption Performance of MIL-101/FLG Composite
The measuring methods of adsorption isotherms of water on MIL-101 and MIL-101/FLG composites and thermogravimetric-differential analysis (TG-DTA) curves of water on MIL-101 and MIL-101/FLG composites were described in our previous work [56]. The adsorption isotherms of water on MIL-101/FLG and MIL-101 were determined at 25 °C on a BELSORP-max (volumetric method) adsorption analyzer, while the adsorption rate curves of water on MIL-101/FLG and MIL-101 were measured at 25 °C by a high vacuum adsorption instrument (gravimetric method) and the adsorption pressure was 2500-2600 Pa.

Determination Method of Adsorption Cooling Performance of MIL-101/FLG-Water Working Pairs
The heat transfer characteristic of MIL-101/20%FLG adsorbers and adsorption cooling performance of MIL-101/20%FLG-water working pairs were investigated on a cyclic adsorption cooling device (shown in Figure 2). Thermophysical properties of MIL-101 and MIL-101/FLG composite were measured at 30-130 • C on a TPS-2500S type thermal property instrument (Hot-Disk Co., Ltd., Uppsala, Sweden). The thermal conductivity and specific heat capacity were measured through the "ISO22007-2-2015" standard [57,58]. The theoretical calculation equations of thermophysical properties refer to our previous work [56]. The standard deviations (error bars) of experimental results are shown in the results and discussion.

Determination of Adsorption and Desorption Performance of MIL-101/FLG Composite
The measuring methods of adsorption isotherms of water on MIL-101 and MIL-101/FLG composites and thermogravimetric-differential analysis (TG-DTA) curves of water on MIL-101 and MIL-101/FLG composites were described in our previous work [56]. The adsorption isotherms of water on MIL-101/FLG and MIL-101 were determined at 25 • C on a BELSORP-max (volumetric method) adsorption analyzer, while the adsorption rate curves of water on MIL-101/FLG and MIL-101 were measured at 25 • C by a high vacuum adsorption instrument (gravimetric method) and the adsorption pressure was 2500-2600 Pa.

Determination Method of Adsorption Cooling Performance of MIL-101/FLG-Water Working Pairs
The heat transfer characteristic of MIL-101/20%FLG adsorbers and adsorption cooling performance of MIL-101/20%FLG-water working pairs were investigated on a cyclic adsorption cooling device (shown in Figure 2). The main components in an adsorption cooling device are the adsorber, condenser, evaporator, refrigerator, heating furnace, vacuum pump and level gauge (0-220 mm, accuracy of 0.1 mm). The central and outer wall temperatures of the adsorber were measured by a Pt-100 type resistance (−30-200 °C, accuracy of ±0.1 °C). The center and outer wall temperatures of the adsorber were real-time recorded by a WP-MD807-71-23-HL type (accuracy of ±0.1 °C) intelligent multiway inspection display controller. The adsorber-heat exchanger was a glass cylinder (inner diameter of 32 mm and height of 400 mm) and total mass of 314.2 g. As seen in Figure 1c,d, the MIL-101/FLG and MIL-101 particles (particle size of about 380-830 μm) for loading into the adsorber were obtained by crushing and sieving the cylindrical specimens formed by the tablet press. The loading mass of MIL-101/FLG and MIL-101 in the adsorber was 25.0 g and 24.5 g, respectively; meanwhile, the loading volume of MIL-101/FLG and MIL-101 adsorption beds was 0.156 L and 0.142 L, respectively. Adsorption cooling performance was directly investigated on the MIL-101/FLG composite or MIL-101 fixed-bed under the same system configuration. Refer to the literature for the detailed experimental procedures [59].
where Qc is the cooling capacity per unit of cycle time, W; Δmad is the amount of cycled water, kg; ΔHe is the vaporization latent heat of water at the evaporation temperature, kJ kg −1 ; Δt is the cycle time, s; Qevap is the evaporative cooling capacity, kJ.
Q input =C a m a T 2 -T 1 +C lc m lc T 2 -T 1 +C m m m T 2 -T 1 +m lc h d The main components in an adsorption cooling device are the adsorber, condenser, evaporator, refrigerator, heating furnace, vacuum pump and level gauge (0-220 mm, accuracy of 0.1 mm). The central and outer wall temperatures of the adsorber were measured by a Pt-100 type resistance (−30-200 • C, accuracy of ±0.1 • C). The center and outer wall temperatures of the adsorber were real-time recorded by a WP-MD807-71-23-HL type (accuracy of ±0.1 • C) intelligent multiway inspection display controller. The adsorber-heat exchanger was a glass cylinder (inner diameter of 32 mm and height of 400 mm) and total mass of 314.2 g. As seen in Figure 1c,d, the MIL-101/FLG and MIL-101 particles (particle size of about 380-830 µm) for loading into the adsorber were obtained by crushing and sieving the cylindrical specimens formed by the tablet press. The loading mass of MIL-101/FLG and MIL-101 in the adsorber was 25.0 g and 24.5 g, respectively; meanwhile, the loading volume of MIL-101/FLG and MIL-101 adsorption beds was 0.156 L and 0.142 L, respectively. Adsorption cooling performance was directly investigated on the MIL-101/FLG composite or MIL-101 fixed-bed under the same system configuration. Refer to the literature for the detailed experimental procedures [59].
where Q c is the cooling capacity per unit of cycle time, W; ∆mad is the amount of cycled water, kg; ∆H e is the vaporization latent heat of water at the evaporation temperature, kJ kg −1 ; ∆t is the cycle time, s; Q evap is the evaporative cooling capacity, kJ.
where Q input is the heat absorbed during desorption process, kJ; C a , C lc and C m are the average specific heat capacity of the adsorbent, water and adsorber-heat exchanger, re- spectively, J kg −1 • C −1 ; m lc and m m are the mass of desorption of water and adsorber-heat exchanger, respectively, kg; h d is the desorption enthalpy of water, kJ kg −1 ; T 1 and T 2 are the initial and final temperatures of the adsorber during the heating desorption process, respectively, • C.
where SCP is the specific cooling power of the system per unit mass of adsorbent, W kg −1 ; m a is the adsorbent mass, kg; COP is the coefficient of performance of the system. The standard deviations (error bars) of repeat experimental results are shown in the results and discussion. The main uncertainties of Q c , SCP and COP are mainly from the level gauge ( Figure 2). The relative error of results was calculated by Equation (6) [47].
Herein, the relative error is 1.1% for Q c , SCP and COP.

Thermal Conductivity of MIL-101/FLG Composites
Thermal conductivities of MIL-101/FLG composites with FLG mass fraction of 5-20% were investigated between 30 • C and 130 • C and the results are shown in Figure 3. where Qinput is the heat absorbed during desorption process, kJ; Ca, Clc and Cm are the average specific heat capacity of the adsorbent, water and adsorber-heat exchanger, respectively, J kg −1 °C −1 ; mlc and mm are the mass of desorption of water and adsorber-heat exchanger, respectively, kg; hd is the desorption enthalpy of water, kJ kg −1 ; T1 and T2 are the initial and final temperatures of the adsorber during the heating desorption process, respectively, °C.
where SCP is the specific cooling power of the system per unit mass of adsorbent, W kg −1 ; ma is the adsorbent mass, kg; COP is the coefficient of performance of the system. The standard deviations (error bars) of repeat experimental results are shown in the results and discussion. The main uncertainties of Qc, SCP and COP are mainly from the level gauge ( Figure 2). The relative error of results was calculated by Equation 6 [47].
Herein, the relative error is 1.1% for Qc, SCP and COP.

Thermal Conductivity of MIL-101/FLG Composites
Thermal conductivities of MIL-101/FLG composites with FLG mass fraction of 5-20% were investigated between 30 °C and 130 °C and the results are shown in Figure 3.  Seen from Figure 3, the thermal conductivity of MIL-101 is 0.1272-0.1486 W m −1 K −1 at 30-130 • C. The thermal conductivity of MIL-101/20%FLG is 0.8322-0.8603 W m −1 K −1 at 30-130 • C, which is 5.79-6.54 times that of MIL-101, respectively. The thermal conductivity of the MIL-101/FLG composite is below 0.5 W m −1 K −1 at a low adding amount of FLG (<10%) owing to a poor heat conduction channel. When the content of FLG in MIL-101/FLG composites reaches 20% and is evenly dispersed in MIL-101, the effective thermal conductive pathways can be formed in the MIL-101/FLG composite. Hence, the heat could be quickly transferred through the thermal conductive pathways composed of FLG with minimum thermal resistance, leading to a sudden change of thermal conductivity. The suitable content of FLG added in the MIL-101/FLG composite is 20%.

Adsorption and Desorption Performance of MIL-101/FLG Composite
The adsorption isotherms and adsorption rate curve of water on MIL-101/20%FLG composites were measured and shown in Figures S1 and S2, while the TG-DTA curves of MIL-101/20%FLG composites were measured and displayed in Figure S3.
As shown in Figure S1, the MIL-101/20%FLG composite could be suitable in adsorption cooling/chiller applications with an evaporation temperature above 15 • C due to the steep uptake of water at P/P s of 0.40-0.50. Based on the TG-DTA curves of water on MIL-101/20%FLG ( Figure S3), the desorption temperature of the MIL-101/20%FLG-water working pair was selected as 70-130 • C according to specific application conditions.

Temperature Characteristic Curves of MIL-101/FLG Adsorber in Adsorption Cooling Process
The heat transfer characteristic of the MIL-101/FLG adsorber was evaluated by the temperature characteristic curves of MIL-101/20%FLG adsorbers during the adsorption, desorption and cooling process on the adsorption cooling device ( Figure 2).

Temperature Characteristic Curves during Adsorption Process
The temperature characteristic curves of MIL-101 and MIL-101/20%FLG adsorbers during the adsorption process were measured at an ambient temperature of 25-30 • C, adsorption temperature of 25-30 • C, evaporation temperature of 17 • C and adsorption time of 120 min, with the results shown in Figure 4.
conductive pathways can be formed in the MIL-101/FLG composite. Hence, the heat could be quickly transferred through the thermal conductive pathways composed of FLG with minimum thermal resistance, leading to a sudden change of thermal conductivity. The suitable content of FLG added in the MIL-101/FLG composite is 20%.

Adsorption and Desorption Performance of MIL-101/FLG Composite
The adsorption isotherms and adsorption rate curve of water on MIL-101/20%FLG composites were measured and shown in Figures S1 and S2, while the TG-DTA curves of MIL-101/20%FLG composites were measured and displayed in Figure S3.
As shown in Figure S1, the MIL-101/20%FLG composite could be suitable in adsorption cooling/chiller applications with an evaporation temperature above 15 °C due to the steep uptake of water at P/Ps of 0.40-0.50. Based on the TG-DTA curves of water on MIL-101/20%FLG ( Figure S3), the desorption temperature of the MIL-101/20%FLG-water working pair was selected as 70-130 °C according to specific application conditions.

Temperature Characteristic Curves of MIL-101/FLG Adsorber in Adsorption Cooling Process
The heat transfer characteristic of the MIL-101/FLG adsorber was evaluated by the temperature characteristic curves of MIL-101/20%FLG adsorbers during the adsorption, desorption and cooling process on the adsorption cooling device ( Figure 2).

Temperature Characteristic Curves During Adsorption Process
The temperature characteristic curves of MIL-101 and MIL-101/20%FLG adsorbers during the adsorption process were measured at an ambient temperature of 25-30 °C, adsorption temperature of 25-30 °C, evaporation temperature of 17 °C and adsorption time of 120 min, with the results shown in Figure 4.  As shown in Figure 4, within the initial stage of adsorption process (adsorption time of 9 min), the temperature of MIL-101 adsorbers rises from 29.0 °C to 41.7 °C (maximum temperature rise of 12.7 °C), while the temperature of MIL-101/20%FLG adsorbers rises from 29.1 °C to 32.6 °C and the maximum temperature rise is merely 3.5 °C, which is significantly lower than MIL-101 adsorbers. Since the equilibrium adsorption capacity of MIL-101 for water (1.14 g g −1 ads) is higher than MIL-101/20%FLG (0.90 g g −1 ads), a larger amount of adsorption heat is released by the adsorption process of water on MIL-101 and a higher temperature rise is caused, but the maximum temperature rise for MIL- As shown in Figure 4, within the initial stage of adsorption process (adsorption time of 9 min), the temperature of MIL-101 adsorbers rises from 29.0 • C to 41.7 • C (maximum temperature rise of 12.7 • C), while the temperature of MIL-101/20%FLG adsorbers rises from 29.1 • C to 32.6 • C and the maximum temperature rise is merely 3.5 • C, which is significantly lower than MIL-101 adsorbers. Since the equilibrium adsorption capacity of MIL-101 for water (1.14 g g −1 ads ) is higher than MIL-101/20%FLG (0.90 g g −1 ads ), a larger amount of adsorption heat is released by the adsorption process of water on MIL-101 and a higher temperature rise is caused, but the maximum temperature rise for MIL-101/20%FLG adsorbers per unit mass of MIL-101 is just 0.179 • C g −1 and still far lower than that of MIL-101 adsorbers (0.508 • C g −1 ). Hence, the lower maximum temperature rise of MIL-101/20%FLG adsorbers is the combination of adsorption capacity (low heat release rate) and better heat conduction property of MIL-101/20%FLG, so the generated adsorption heat could quickly be transferred to the external environment. As the adsorption process goes on, the adsorption rate decreases while the adsorption heat reduces; it takes 30 min (adsorption time of 6-36 min) for MIL-101/20%FLG adsorbers to drop from 32.6 • C to 28.0 • C with a cooling rate of 0.153 • C min −1 , while the temperature of MIL-101 adsorbers just falls from 32.6 • C to 29.9 • C in the same duration (adsorption time of 40-70 min) and the cooling rate is merely 0.09 • C min −1 . At the low temperature region, the cooling rate of MIL-101/20%FLG adsorbers is 1.7 times that of MIL-101 adsorbers. The thermal conductivity of MIL-101 adsorbers increases, and the rapid removal of adsorption heat is promoted by adding FLG; hence, the temperature rise of MIL-101/20%FLG adsorbers can be reduced during the adsorption process as the cooling rate is accelerated.

Heating Curves in the Desorption Process
The temperature characteristic curves of MIL-101 and MIL-101/20%FLG adsorbers during the heating desorption process were compared at an ambient temperature of 25-30 • C, while the initial temperature of MIL-101 and MIL-101/20%FLG adsorbers was 28.7 • C and 29.0 • C, and the desorption temperature was 130 • C with a desorption time of 40 min. The results are shown in Figure 5.
101/20%FLG adsorbers per unit mass of MIL-101 is just 0.179 °C g −1 and still far lower than that of MIL-101 adsorbers (0.508 °C g −1 ). Hence, the lower maximum temperature rise of MIL-101/20%FLG adsorbers is the combination of adsorption capacity (low heat release rate) and better heat conduction property of MIL-101/20%FLG, so the generated adsorption heat could quickly be transferred to the external environment. As the adsorption process goes on, the adsorption rate decreases while the adsorption heat reduces; it takes 30 min (adsorption time of 6-36 min) for MIL-101/20%FLG adsorbers to drop from 32.6 °C to 28.0 °C with a cooling rate of 0.153 °C min −1 , while the temperature of MIL-101 adsorbers just falls from 32.6 °C to 29.9 °C in the same duration (adsorption time of 40-70 min) and the cooling rate is merely 0.09 °C min −1 . At the low temperature region, the cooling rate of MIL-101/20%FLG adsorbers is 1.7 times that of MIL-101 adsorbers. The thermal conductivity of MIL-101 adsorbers increases, and the rapid removal of adsorption heat is promoted by adding FLG; hence, the temperature rise of MIL-101/20%FLG adsorbers can be reduced during the adsorption process as the cooling rate is accelerated.

Heating Curves in the Desorption Process
The temperature characteristic curves of MIL-101 and MIL-101/20%FLG adsorbers during the heating desorption process were compared at an ambient temperature of 25-30 °C, while the initial temperature of MIL-101 and MIL-101/20%FLG adsorbers was 28.7 °C and 29.0 °C, and the desorption temperature was 130 °C with a desorption time of 40 min. The results are shown in Figure 5. 2%. This is due to the fact that a thermal conductive channel can be formed within the MIL-101/FLG composite by adding FLG into MIL-101, which improves thermal conduction performance of MIL-101 adsorbers, and the thermal contact resistance between MIL-101 particles is reduced. Therefore, the heat transfer of MIL-101 adsorbers could be effectively improved 2%. This is due to the fact that a thermal conductive channel can be formed within the MIL-101/FLG composite by adding FLG into MIL-101, which improves thermal conduction performance of MIL-101 adsorbers, and the thermal contact resistance between MIL-101 particles is reduced. Therefore, the heat transfer of MIL-101 adsorbers could be effectively improved by adding FLG during the heating desorption process, which would be beneficial to reduce desorption time.

Cooling Rate
At an initial temperature of~130 • C and final temperature of~29 • C, the temperature drop curves of MIL-101 and MIL-101/20%FLG adsorbers were measured in the cooling process. The results are shown in Figure 6. by adding FLG during the heating desorption process, which would be beneficial to reduce desorption time.

Cooling Rate
At an initial temperature of ~130 °C and final temperature of ~29 °C, the temperature drop curves of MIL-101 and MIL-101/20%FLG adsorbers were measured in the cooling process. The results are shown in Figure 6. As seen in Figure 6, the cooling time of an MIL-101/20%FLG adsorber consists of only 27 min to drop from 130.6 °C to 29.0 °C with a cooling rate of 3.76 °C min −1 , while an MIL-101 adsorber needs 60 min to drop from 130.4 °C to 29.2 °C with a cooling rate of 1.69 °C min −1 . The cooling rate of an MIL-101/20%FLG adsorber is 2.22 times that of an MIL-101 adsorber; hence, the cooling time is shortened by 55%. In the continuous adsorption cooling cycle, compared to an MIL-101 adsorber, the MIL-101/20%FLG adsorber needs a shorter time to switch to the next cycle, which is conducive to improving the SCP of an adsorption cooling system.

Heat Transfer Mechanism
The phase structure, SEM image and Raman spectra of FLG were determined and shown in Figure S4. The results indicate that FLG is a layered structure without surface defects. In combination with the XRD, SEM and Raman characterization of FLG ( Figure  S4), a schematic diagram of the structure and heat transfer mechanism for an MIL-101/FLG composite were analyzed in Figure 7. As seen in Figure 6, the cooling time of an MIL-101/20%FLG adsorber consists of only 27 min to drop from 130.6 • C to 29.0 • C with a cooling rate of 3.76 • C min −1 , while an MIL-101 adsorber needs 60 min to drop from 130.4 • C to 29.2 • C with a cooling rate of 1.69 • C min −1 . The cooling rate of an MIL-101/20%FLG adsorber is 2.22 times that of an MIL-101 adsorber; hence, the cooling time is shortened by 55%. In the continuous adsorption cooling cycle, compared to an MIL-101 adsorber, the MIL-101/20%FLG adsorber needs a shorter time to switch to the next cycle, which is conducive to improving the SCP of an adsorption cooling system.

Heat Transfer Mechanism
The phase structure, SEM image and Raman spectra of FLG were determined and shown in Figure S4. The results indicate that FLG is a layered structure without surface defects. In combination with the XRD, SEM and Raman characterization of FLG ( Figure S4), a schematic diagram of the structure and heat transfer mechanism for an MIL-101/FLG composite were analyzed in Figure 7.
It can be seen in Figure 7a that heat in FLG is both conducted by phonon (lattice vibration wave) heat conduction and electronic heat conduction. As seen in Figure 7b, MIL-101 is a three-dimensional framework mesoporous material, and heat conduction in MIL-101 is mainly achieved by phonon diffusion. Due to its high porosity, severe phonon scattering occurs when heat is transferred in the MIL-101 skeleton, and there is large contact thermal resistance between MIL-101 particles, leading to a low thermal conductivity of only 0.1272-0.1486 W m −1 K −1 . As a consequence, the heating and cooling rate of the MIL-101 adsorber was affected. As shown in Figure 7c, since interfacial mutual interactions between FLG and MIL-101 can be promoted after a ball-milling process, a π-π stacking effect may be formed between the benzene ring of MIL-101 framework and hex-atomic ring in the layered FLG; namely, certain interactions exist between FLG and MIL-101. Moreover, fast migration of free electrons can be promoted due to the two-dimensional plane structure of FLG (with large surface area). Therefore, thermal resistance between MIL-101 particles is reduced and the thermal conduction of MIL-101/20%FLG composite is effectively improved. As a result, the heating and cooling rate of MIL-101/20%FLG adsorber increases while heating and cooling time is correspondingly shortened. It can be seen in Figure 7a that heat in FLG is both conducted by phonon (lattice vibration wave) heat conduction and electronic heat conduction. As seen in Figure 7b, MIL-101 is a three-dimensional framework mesoporous material, and heat conduction in MIL-101 is mainly achieved by phonon diffusion. Due to its high porosity, severe phonon scattering occurs when heat is transferred in the MIL-101 skeleton, and there is large contact thermal resistance between MIL-101 particles, leading to a low thermal conductivity of only 0.1272-0.1486 W m −1 K −1 . As a consequence, the heating and cooling rate of the MIL-101 adsorber was affected. As shown in Figure 7c, since interfacial mutual interactions between FLG and MIL-101 can be promoted after a ball-milling process, a π-π stacking effect may be formed between the benzene ring of MIL-101 framework and hex-atomic ring in the layered FLG; namely, certain interactions exist between FLG and MIL-101. Moreover, fast migration of free electrons can be promoted due to the two-dimensional plane structure of FLG (with large surface area). Therefore, thermal resistance between MIL-101 particles is reduced and the thermal conduction of MIL-101/20%FLG composite is effectively improved. As a result, the heating and cooling rate of MIL-101/20%FLG adsorber increases while heating and cooling time is correspondingly shortened.

Adsorption Cooling Performance of MIL-101/20%FLG-Water Working Pair
The effects of working condition parameters on the adsorption cooling performances (cooling capacity and cooling power) of an MIL-101/20%FLG-water working pair were investigated on the adsorption cooling device (Figure 2). The working condition parameters of the MIL-101-water working pair are shown in Figures S5 and S6.

Working Condition Parameters Evaporation Temperature
The cooling capacity and cooling power (SCP) of the MIL-101/20%FLG-water working pair at different evaporation temperatures (10-20 • C) were investigated at an adsorption time of 60 min and desorption time of 60 min. The results are shown in Figure 8.
As seen in Figure 8, the cooling capacity and SCP of the system is 2.14 W and 87.5 W kg −1 at an evaporation temperature of 15 • C, which is 3.06 times that of an evaporation temperature of 10 • C (0.70 W and 28.6 W kg −1 ). Raising the evaporation temperature to 17 • C and 20 • C, the cooling capacity and SCP of an MIL-101/20%FLG-water working pair just increases by 10.3% and 10%. This is due to the fact that the evaporation pressure of refrigerant (water) rises, and the mass transfer impetus of water on MIL-101/20%FLG increases with the rise of evaporation temperature. Results show that the MIL-101/20%FLG composite is more applicable to an adsorption water chiller application [12,35] with an evaporation temperature of 15-20 • C.

Working Condition Parameters
Evaporation Temperature The cooling capacity and cooling power (SCP) of the MIL-101/20%FLG-water working pair at different evaporation temperatures (10-20 °C) were investigated at an adsorption time of 60 min and desorption time of 60 min. The results are shown in Figure 8. As seen in Figure 8, the cooling capacity and SCP of the system is 2.14 W and 87.5 W kg −1 at an evaporation temperature of 15 °C, which is 3.06 times that of an evaporation temperature of 10 °C (0.70 W and 28.6 W kg −1 ). Raising the evaporation temperature to 17 °C and 20 °C, the cooling capacity and SCP of an MIL-101/20%FLG-water working pair just increases by 10.3% and 10%. This is due to the fact that the evaporation pressure of refrigerant (water) rises, and the mass transfer impetus of water on MIL-101/20%FLG increases with the rise of evaporation temperature. Results show that the MIL-101/20%FLG composite is more applicable to an adsorption water chiller application [12,35] with an evaporation temperature of 15-20 °C.

Adsorption Time
The cooling capacity and SCP of an MIL-101/20%FLG-water working pair at different adsorption times (30-300 min) were investigated at an evaporation temperature of 17 °C and desorption time of 60 min. The results are shown in Figure 9. Seen from Figure 9, the cooling capacity of the system continuously increases with a prolonged adsorption time and reaches up to 2.23-2.40 W at an adsorption time of 90-120 min. The SCP of the system firstly increases and then drops rapidly with the extension of adsorption time. SCP reaches 98.0 W kg −1 under an adsorption time between 90 and 120 min. This is due to the fact that after the initial adsorption stage, the adsorption process of water on MIL-101/20%FLG is mainly in the form of slow capillary condensation in mesoporous cages. The adsorption rate slows down, and the increase in cooling capacity is offset by the extension of adsorption time. Therefore, a reasonable adsorption time of an MIL-101/20%FLG-water working pair is 90-120 min, which is shortened by 25% compared to an MIL-101-water working pair whose optimal adsorption time is 120-180 min (as shown in Figure S5). This is because the removal rate of adsorption heat and cooling rate are quicker for MIL-101/20%FLG. The adsorption rate of MIL-101/20%FLG for water is relatively increased (as shown in Figure S2).
Compared to the literature [61] (adsorption time of 90-100 min), the adsorption time in this work could be further shortened by optimizing the adsorber configuration to strengthen mass transfer properties. Seen from Figure 9, the cooling capacity of the system continuously increases with a prolonged adsorption time and reaches up to 2.23-2.40 W at an adsorption time of 90-120 min. The SCP of the system firstly increases and then drops rapidly with the extension of adsorption time. SCP reaches 98.0 W kg −1 under an adsorption time between 90 and 120 min. This is due to the fact that after the initial adsorption stage, the adsorption process of water on MIL-101/20%FLG is mainly in the form of slow capillary condensation in mesoporous cages. The adsorption rate slows down, and the increase in cooling capacity is offset by the extension of adsorption time. Therefore, a reasonable adsorption time of an MIL-101/20%FLG-water working pair is 90-120 min, which is shortened by 25% compared to an MIL-101-water working pair whose optimal adsorption time is 120-180 min (as shown in Figure S5). This is because the removal rate of adsorption heat and cooling rate are quicker for MIL-101/20%FLG. The adsorption rate of MIL-101/20%FLG for water is relatively increased (as shown in Figure S2).
Compared to the literature [61] (adsorption time of 90-100 min), the adsorption time in this work could be further shortened by optimizing the adsorber configuration to strengthen mass transfer properties.

Desorption Time
The cooling capacity and SCP of MIL-101/20%FLG-water working pairs at different desorption times (20-60 min) were investigated at an evaporation temperature of 17 • C and adsorption time of 90 and 120 min. The results are shown in Figure 10. As shown in Figure 10, the cooling capacity (a) and SCP (b) rapidly increase below a desorption time of 30 min due to the fast desorption rate, while the cooling capacity of the system is 2.74 W and 2.48 W at adsorption times of 90 min and 120 min, respectively. Prolonging the desorption time to 40 min, the cooling capacity of the system is 2.78 W and 2.49 W at adsorption times of 90 min and 120 min, respectively. However, there is no obvious improvement of cooling capacity at a desorption time greater than 40 min. When desorption time is in the range of 30-40 min, the SCP of the MIL-101/20%FLG-water working pair reaches its maximum of 113.5 W kg −1 at an adsorption time of 90 min. Along with desorption time further increasing, the SCP significantly declines. Therefore, the suitable desorption time of an MIL-101/20%FLG-water working pair is 40 min, which is shortened by 33.3% compared to an MIL-101-water working pair with a desorption time of 60 min (as shown in Figure S6). This is attributed to the higher thermal conduction of MIL-101/20%FLG compared to MIL-101, resulting in a faster heating rate and shorter heating time; thus, the desorption rate is increased and desorption time is decreased.

Adsorption Cooling Performance Comparison of MIL-101/20%FLG-Water and MIL-101-Water Working Pairs
The MIL-101/FLG composite may be more applicable to adsorption water chiller applications with an evaporation temperature of 15-20 °C according to sections 3.  Figure 11, while the results of cooling capacity, SCP and COP are listed in Table 1. As shown in Figure 10, the cooling capacity (a) and SCP (b) rapidly increase below a desorption time of 30 min due to the fast desorption rate, while the cooling capacity of the system is 2.74 W and 2.48 W at adsorption times of 90 min and 120 min, respectively. Prolonging the desorption time to 40 min, the cooling capacity of the system is 2.78 W and 2.49 W at adsorption times of 90 min and 120 min, respectively. However, there is no obvious improvement of cooling capacity at a desorption time greater than 40 min. When desorption time is in the range of 30-40 min, the SCP of the MIL-101/20%FLG-water working pair reaches its maximum of 113.5 W kg −1 at an adsorption time of 90 min. Along with desorption time further increasing, the SCP significantly declines. Therefore, the suitable desorption time of an MIL-101/20%FLG-water working pair is 40 min, which is shortened by 33.3% compared to an MIL-101-water working pair with a desorption time of 60 min (as shown in Figure S6). This is attributed to the higher thermal conduction of MIL-101/20%FLG compared to MIL-101, resulting in a faster heating rate and shorter heating time; thus, the desorption rate is increased and desorption time is decreased.

Adsorption Cooling Performance Comparison of MIL-101/20%FLG-Water and MIL-101-Water Working Pairs
The MIL-101/FLG composite may be more applicable to adsorption water chiller applications with an evaporation temperature of 15-20 • C according to refrigerator are shown in Figure 11, while the results of cooling capacity, SCP and COP are listed in Table 1. It is shown in Figure 11 that during the cooling process, the temperature drop the of evaporator and refrigerator for MIL-101/20%FLG-water working pairs is 3.6 °C and 3.3 °C, respectively, which means chilled water of 16-20 °C might be produced. As seen in Table 1, at desorption temperatures of 70 °C and 90 °C, the cooling capacity, SCP and COP of MIL-101/20%FLG-water working pairs is 1.77-1.98 W, 72.2-81.0 W kg −1 and 0.187-0.202, respectively, which is 1.43-1.56 times higher than MIL-101-water working pairs. As mentioned above, this is mainly because the thermal conduction of MIL-101 is strengthened, and the heat transfer of an MIL-101 adsorber during the adsorption and desorption process is improved by addition of FLG, which is conducive to reducing the cycle time of water on MIL-101. The results indicate that the cooling performance of MIL-101-water working pairs could be significantly enhanced by adding FLG, and the energy conversion efficiency in adsorption water chiller applications might be improved. This provides a promising adsorbent candidate (MIL-101/FLG composite) in adsorption water chiller applications such as data center cooling, proton exchange membrane fuel cell, automotive engine coolant and vehicle air-conditioning; especially, a better cooling effect would be achieved in the vehicle air-conditioning application (evaporation temperature of 15-20 °C and feasible heat source temperature of 100-130 °C).

Cyclic Stability
The cyclic cooling capacity of an MIL-101/20%FLG-water working pair was investigated at an adsorption temperature of 30 °C, evaporation temperature of 20 °C, desorption temperature of 90 °C, condensation temperature of 30 °C, adsorption time of 90 min and desorption time of 30 min. The result is shown in Figure 12. It is shown in Figure 11 that during the cooling process, the temperature drop the of evaporator and refrigerator for MIL-101/20%FLG-water working pairs is 3.6 • C and 3.3 • C, respectively, which means chilled water of 16-20 • C might be produced.
As seen in Table 1, at desorption temperatures of 70 • C and 90 • C, the cooling capacity, SCP and COP of MIL-101/20%FLG-water working pairs is 1.77-1.98 W, 72.2-81.0 W kg −1 and 0.187-0.202, respectively, which is 1.43-1.56 times higher than MIL-101-water working pairs. As mentioned above, this is mainly because the thermal conduction of MIL-101 is strengthened, and the heat transfer of an MIL-101 adsorber during the adsorption and desorption process is improved by addition of FLG, which is conducive to reducing the cycle time of water on MIL-101. The results indicate that the cooling performance of MIL-101-water working pairs could be significantly enhanced by adding FLG, and the energy conversion efficiency in adsorption water chiller applications might be improved. This provides a promising adsorbent candidate (MIL-101/FLG composite) in adsorption water chiller applications such as data center cooling, proton exchange membrane fuel cell, automotive engine coolant and vehicle air-conditioning; especially, a better cooling effect would be achieved in the vehicle air-conditioning application (evaporation temperature of 15-20 • C and feasible heat source temperature of 100-130 • C).

Cyclic Stability
The cyclic cooling capacity of an MIL-101/20%FLG-water working pair was investigated at an adsorption temperature of 30 • C, evaporation temperature of 20 • C, desorption temperature of 90 • C, condensation temperature of 30 • C, adsorption time of 90 min and desorption time of 30 min. The result is shown in Figure 12. As seen in Figure 12, the cooling capacity of an MIL-101/20%FLG-water working pair is maintained at 1.96-2.01 W after 50 continuous adsorption/desorption cycles.

Characterization of Cycled MIL-101/20%FLG
The XRD patterns of fresh MIL-101/20%FLG and cycled MIL-101/20%FLG are shown in Figure 13, while N2 adsorption/desorption isotherms and pore structure parameters are shown in Figure 14 and Table 2, respectively.  Figure 13, the major peaks of MIL-101 components are retained in MIL-101/20%FLG, which indicates well the coexistence between FLG and MIL-101. After continuous adsorption cooling cycles, the crystallinity of cycled MIL-101/20%FLG is rarely changed and the phase structure of MIL-101/20%FLG remains stable. As seen in Figure 12, the cooling capacity of an MIL-101/20%FLG-water working pair is maintained at 1.96-2.01 W after 50 continuous adsorption/desorption cycles.

Characterization of Cycled MIL-101/20%FLG
The XRD patterns of fresh MIL-101/20%FLG and cycled MIL-101/20%FLG are shown in Figure 13, while N 2 adsorption/desorption isotherms and pore structure parameters are shown in Figure 14 and Table 2, respectively. As seen in Figure 12, the cooling capacity of an MIL-101/20%FLG-water working pair is maintained at 1.96-2.01 W after 50 continuous adsorption/desorption cycles.

Characterization of Cycled MIL-101/20%FLG
The XRD patterns of fresh MIL-101/20%FLG and cycled MIL-101/20%FLG are shown in Figure 13, while N2 adsorption/desorption isotherms and pore structure parameters are shown in Figure 14 and Table 2, respectively.  Figure 13, the major peaks of MIL-101 components are retained in MIL-101/20%FLG, which indicates well the coexistence between FLG and MIL-101. After continuous adsorption cooling cycles, the crystallinity of cycled MIL-101/20%FLG is rarely changed and the phase structure of MIL-101/20%FLG remains stable.  It is shown in Figure 14 that the N2 adsorption uptake of cycled MIL-101/20%FLG and fresh MIL-101/20%FLG is similar. Seen from Table 2, the BET surface area of fresh MIL-101/20%FLG is decreased by 23.9%, and the pore volume is reduced by 24.1%, compared to MIL-101 [37], which is consistent with the adding proportion of FLG. MIL-101 is not significantly destructed by the addition of FLG. After continuous adsorption and desorption cycles, the porous structure of cycled MIL-101/20%FLG is similar to fresh MIL-101/20%FLG.

Conclusions
In this work, a novel MIL-101/FLG composite was prepared for improving the thermal conduction performance of MIL-101. The heat transfer characteristic of MIL-101/FLG adsorber was evaluated during adsorption, desorption and cooling process. The adsorption cooling performance of an MIL-101/FLG-water working pair was investigated for typical adsorption water chiller working conditions. The thermal conductivity of the MIL-101/20%FLG composite is 0.8322-0.8603 W m −1 K −1 . The heating and cooling rate of the MIL-101/20%FLG adsorber is 2.2 times as much as that of an MIL-101 adsorber; meanwhile, the heating and cooling time is effectively shortened by ~55%. The enhanced heat transfer of an MIL-101/20%FLG adsorber is conducive to the improvement of the cooling performance of MIL-101-water working pairs. Under adsorption water chiller working conditions, the cooling capacity, SCP and COP of an MIL-101/20%FLG-water working pair are 1.77-1.98 W, 72.2-81.0 W kg −1 and 0.187-0.202, respectively. Moreover, the excellent adsorption/desorption cyclic stability of MIL-101/FLG composites is successfully verified. This work provides a reference for improving the heat transfer performance of MOFs and gives a promising material candidate for the MIL-101/FLG composite in research directions for adsorption water chillers (16-20 °C), including data center cooling,  As shown in Figure 13, the major peaks of MIL-101 components are retained in MIL-101/20%FLG, which indicates well the coexistence between FLG and MIL-101. After continuous adsorption cooling cycles, the crystallinity of cycled MIL-101/20%FLG is rarely changed and the phase structure of MIL-101/20%FLG remains stable.
It is shown in Figure 14 that the N 2 adsorption uptake of cycled MIL-101/20%FLG and fresh MIL-101/20%FLG is similar. Seen from Table 2, the BET surface area of fresh MIL-101/20%FLG is decreased by 23.9%, and the pore volume is reduced by 24.1%, compared to MIL-101 [37], which is consistent with the adding proportion of FLG. MIL-101 is not significantly destructed by the addition of FLG. After continuous adsorption and desorption cycles, the porous structure of cycled MIL-101/20%FLG is similar to fresh MIL-101/20%FLG.

Conclusions
In this work, a novel MIL-101/FLG composite was prepared for improving the thermal conduction performance of MIL-101. The heat transfer characteristic of MIL-101/FLG adsorber was evaluated during adsorption, desorption and cooling process. The adsorption cooling performance of an MIL-101/FLG-water working pair was investigated for typical adsorption water chiller working conditions. The thermal conductivity of the MIL-101/20%FLG composite is 0.8322-0.8603 W m −1 K −1 . The heating and cooling rate of the MIL-101/20%FLG adsorber is 2.2 times as much as that of an MIL-101 adsorber; meanwhile, the heating and cooling time is effectively shortened by~55%. The enhanced heat transfer of an MIL-101/20%FLG adsorber is conducive to the improvement of the cooling performance of MIL-101-water working pairs. Under adsorption water chiller working conditions, the cooling capacity, SCP and COP of an MIL-101/20%FLG-water working pair are 1.77-1.98 W, 72.2-81.0 W kg −1 and 0.187-0.202, respectively. Moreover, the excellent adsorption/desorption cyclic stability of MIL-101/FLG composites is successfully verified. This work provides a reference for improving the heat transfer performance of MOFs and gives a promising material candidate for the MIL-101/FLG composite in research directions for adsorption water chillers (16-20 • C), including data center cooling, liquid cooling of proton exchange membrane fuel cell, automotive engine coolant and especially vehicle air-conditioning. It is foreseeable that the overall heat transfer and cooling capacity of MIL-101/FLG-water working pairs can be further improved in a new adsorption cooling device configuration.

Conflicts of Interest:
The authors declare no conflict of interest.