Next Article in Journal
Rapid Ferric Transformation by Reductive Dissolution of Schwertmannite for Highly Efficient Catalytic Degradation of Rhodamine B
Previous Article in Journal
Flexible and Stretchable Bio-Integrated Electronics Based on Carbon Nanotube and Graphene
Previous Article in Special Issue
Energy Storage Analysis of a Mixed R161/MOF-5 Nanoparticle Nanofluid Based on Molecular Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Energy Storage of R1234yf, R1234ze, R134a and R32/MOF-74 Nanofluids: A Molecular Simulation Study

1
Key Laboratory of Low-grade Energy Utilization Technology & System, Ministry of Education, College of Power Engineering, Chongqing University, Chongqing 400044, China
2
School of Chemical Engineering, The University of Queensland, Brisbane, QLD 4072, Australia
3
Chongqing Key Laboratory of Heterogeneous Material Mechanics, College of Aerospace Engineering, Chongqing University, Chongqing 400044, China
*
Author to whom correspondence should be addressed.
Materials 2018, 11(7), 1164; https://doi.org/10.3390/ma11071164
Submission received: 16 June 2018 / Revised: 4 July 2018 / Accepted: 6 July 2018 / Published: 8 July 2018
(This article belongs to the Special Issue Nanotechnology in Renewable Energy)

Abstract

:
Thermal energy storage can be carried out by working fluid adsorbing and desorbing in porous materials. In this paper, the energy storage properties of four refrigerants, R1234yf, R1234ze, R134a and R32, with M-metal organic framework (MOF)-74 (M = Zn, Ni, Mg, Co) nanoparticles are investigated using molecular dynamics simulations and grand canonical Monte Carlo simulations. The results show that M-MOF-74 can adsorb more R32 and R134a than R1234yf and R1234ze, as the molecular structures of R32 and R134a are smaller than those of R1234yf and R1234ze. Mg-MOF-74 owns a higher adsorbability than the other MOFs. The energy storage properties of the studied refrigerants can be enhanced when the sum of thermodynamic energy change of MOF particles and the desorption heat of fluid in MOFs is larger than the enthalpy change of pure organic fluid. The R1234yf/M-MOF-74 (M = Co, Mg, Ni) nanofluid can store more energy than other refrigerants/M-MOF-74 (M = Co, Mg, Ni) nanofluid. The energy storage enhancement ratios of R1234yf, R1234ze and R134a with Mg-MOF-74 nanoparticles are higher than those of other M-MOF-74 (M = Co, Ni, Zn) materials.

1. Introduction

The energy crisis and environmental pollution are serious issues in human society, as a consequence of the devastating consumption of fossil fuels such as coal, oil and natural gas. In China, the waste heat generated in industries is equivalent to the combustion energy of 340 million tons of standard coal each year. Finding an efficient way to recover the low-grade energy, such as industrial waste heat, geothermal and solar energy, is, hence, essential to combat these issues [1].
The organic Rankin cycle (ORC) uses an organic refrigerant as a working fluid to recover the heat from low-grade energy resources and turns it into electricity, providing a potential way to use the low-grade energy [2]. Further, the efficiency of ORC in low-grade energy recovery could be remarkably enhanced by increasing the thermal properties of the working fluid. McGrail et al. [3] observed reversible adsorption of R123 molecules in metal organic framework (MOF) nanoparticles. They showed that adding MOF nanoparticles into refrigerant (also known as metal-organic heat carrier nanofluids (MOHCs)) could significantly improve the heat capacity of the working fluid and subsequently enhance the efficiency of ORC in recovering the low-grade energy. During the evaporation process in ORC, excess heat is absorbed by MOHCs from the heat resource for the desorption of adsorbate molecules from MOF particles. Further, excess heat is released by the MOHCs during the condensation process as a consequence of the exothermic adsorption of refrigerant molecules into MOF particles. Additionally, the thermal conductivity of the working fluid is enhanced by adding nanoparticles [4]. The reported research works [5,6] indicated that the enhancement of the thermal properties of nanofluid are mainly due to the effective medium theory. Besides, the Brownian motion, layering phenomenon, clustering, ballistic phonon motion, thermal boundary resistance and mass difference scattering will also impact the thermal properties of nanofluid [7,8].
MOFs are highly crystallized by metal ions or metal clusters and organic ligands through the coordination bond or intermolecular interactions, possessing good thermal stability, high specific surface area and strong adsorption affinity for adsorbates. The adsorptive performance of MOFs is found superior to the traditional nanoporous materials such as zeolites and activated carbons [9], making them a promising material for energy storage and gas separation [10,11]. Among a wide range of MOF structures, MOF-74 owns unique secondary building units that have a strong fault tolerance of different metal particles [12]. Nevertheless, its adsorptive performance for organic refrigerates is little understood [13]. Thus, M-MOF-74 (M = Zn, Ni, Mg, Co) is employed as the adsorbent added in the organic refrigerants in this work.
Hitherto, the development of refrigerants has evolved for four generations, with the third and fourth generation of refrigerants being most widely used in practice. In resemblance to the thermodynamic properties of R410a, R32 (difluoromethane, also known as HFC-32), a halogenated hydrocarbon, is treated as an ideal substitute for the second generation ozone-depleting substances [14]. In addition, R134a (tetrafluoroethane, also known as HFC-134a) with zero ODP (ozone depletion potential) is widely used in industry for its excellent cooling performance [15]. Unlike the third generation refrigerants, R32 and R134a, the fourth generation refrigerants are sustainable and environmental friendly, which are proposed to reduce the emission of fluorinated greenhouse gases [16]. Among the hydrofluoroolefins (HFOs) refrigerants, R1234ze (1,3,3,3-tetrafluoropropene, also known as HFO-1234ze) and R1234yf (2,3,3,3-tetrafluoropropene, also known as HFO-1234yf) have demonstrated to be the best working substances to replace the third generation refrigerant, R134a, for their zero ODP, low GWP (global warming potential) and short residence time in the atmosphere [17]. Meanwhile, it has been widely recognized that both R1234yf and R1234ze are promising working fluids in the refrigeration cycle, heat pump and ORC [18,19].
For comparative purposes, the adsorption and desorption heats of the third and fourth generation refrigerant molecules, R32, R134a, R1234yf and R1234ze, in MOF-74 were measured by using grand canonical Monte Carlo (GCMC) simulations. In this connection, the maximized heat capacity among the four different nanofluidic fluids comprising the organic refrigerant (R32, R1234yf, R1234ze and R134a) and the MOF-74 particles is unveiled with the aid of molecular dynamics (MD) simulations [20,21,22].

2. Materials and Methods

The thermal energy stored in MOHCs can be calculated following [3,23],
Δ h MOHCs = ( 1 x ) Δ h Fluid + x ( T 0 T 1 C p d T ) MOFs + x Δ h desorption
or be written as,
Δ h MOHCs = Δ h Fluid + x ( ( T 0 T 1 C p d T ) MOFs + Δ h desorption Δ h Fluid
where the enthalpy gain of MOHCs (ΔhMOHCs) consists of the enthalpy change of pure organic fluid (ΔhFluid), the energy change of the MOF nanoparticles (( T 0 T 1 C p d T )MOFs) and the enthalpy of desorption (Δhdesorption) of the refrigerant molecules for the temperature difference of T1T0. Here, T1 and T0 are the temperatures of the heating and cooling resources. Cp is the heat capacity of MOF particles, and x is the mass fraction of MOF particles in MOHCs. It can be concluded that the MOHCs are able to store more thermal energy than the pure fluid when the sum of thermodynamic energy change of MOF particles and desorption heat of fluid in MOFs are larger than the enthalpy change of pure organic fluid. In other words, as long as the second term in the right side of Equation (2) is positive, the energy storage in MOHCs is enhanced compared to that in the pure organic fluid.
Since the enthalpy of pure species, hFluid (R32, R1234yf, R1234ze and R134a), is available in NIST [24], the first term on the right side of Equation (2), ΔhFluid, can be determined as ΔhFluid = hFluid.T − hFluid.T0, where T0 and T are the reference and interested temperatures of the working fluid. Meanwhile, the heat capacity of MOF, Cp, was measured in our equilibrium molecular dynamic (EMD) simulations. The last part, Δhdesorption, was calculated as the isosteric heat of adsorption in our GCMC simulations. Subsequently, ΔhMOHCs can be readily calculated.
The EMD and GCMC simulations were performed by using Materials Studio [25]. The COMPASS force field [26] was adopted to describe the intra- and inter-molecular interactions. The Ewald summation method [27] was employed to correct the long-range Coulomb interactions.

2.1. MD Computational Details

M-MOF-74 is formed by divalent metal ions (Zn, Mg, Co, Ni, etc.) connected with ligand 2, 5-dihydroxy terephthalic acid. It forms a 2-dimensional six-party channel and a 3-dimensional space network structure like a honeycomb, and metal atoms are octahedrally coordinated (dominated by five oxygen atoms and a water molecule) [28]. The molecular structure of MOF-74 (M = Ni) employed in this work is depicted in Figure 1. The MOF-74 sample comprised of 1 × 1 × 4 unit cells, with 648 atoms (including 288 carbon atoms, 72 oxygen atoms, 216 hydrogen atoms and 72 Ni atoms) was used in our GCMC and EMD simulations. The lattice parameters for the MOF-74 structure are a = b = 25.7856 Å, c = 27.0804 Å, α = β = 90° and γ = 120°.
The molecular configurations of refrigerants R1234yf, R1234ze, R134a and R32 are presented in Figure 2. In addition, the distribution of the adsorption sites of R1234ze in Co-MOF-74 is shown in Figure 3, which was obtained at 353 K and 5 MPa.
In EMD simulations, the thermodynamic energies of M-MOF-74 nanoparticles at different temperatures (293 K, 313 K, 333 K, 353 K, 373 K and 393 K) were computed in canonical (NVT) ensemble via the Forcite module in Materials Studio. The time step was set as 1 fs, with a 200-ps run being used to equilibrate the MOF structures. A 500-ps run after equilibration was further conducted to obtain the thermodynamic energy. The assignment of charges in MOF structures is based on the COMPASS force field, keeping the structures electronically neutral. The Berendsen thermostat was employed to adjust the temperature of the system. The cut-off distance is 12.5 Å. Periodic boundary conditions were applied in the X, Y and Z directions in all of our GCMC and EMD simulations.

2.2. GCMC Computational Details

The GCMC simulations are able to calculate the adsorption isotherms (293 K, 313 K, 333 K, 353 K, 373 K and 393 K) of refrigerants (R1234yf, R1234ze, R134a and R32) in the M-MOF-74 nanoparticles via the Sorption module in Materials Studio. The pressure range of 1–10,000 kPa is investigated in this work. The fugacity is calculated by the Peng–Robinson equation. For each point of the adsorption isotherm, the equilibration time was 10,000 cycles with another 100,000 cycles for statistics. Here, the simulation runs 5,500,000 cycles to obtain the adsorption isotherm at each temperature (293 K, 313 K, 333 K, 353 K, 373 K and 393 K).

3. Results and Discussion

3.1. Thermodynamic Energy of MOF-74

Our previous work [29] indicated that the present method is an efficient approach to obtain the thermodynamic parameters of MOFs. The changes of the thermodynamic energies of M-MOF-74 structures at different temperatures relative to that at the reference temperature of 293 K are plotted in Figure 4. The thermodynamic energy of M-MOF-74 increases linearly with temperature. The increment of thermodynamic energy per unit temperature is the heat capacity of MOFs, Cp, i.e., ρCp-Co-MOF-74 = 3.34 (J/cm3·K), ρCp-Mg-MOF-74 = 3.74 (J/cm3·K), ρCp-Ni-MOF-74 = 2.14 (J/cm3·K), ρCp-Zn-MOF-74 = 2.62 (J/cm3·K). These values are larger than ρCp-MOF-5 = 1.24 (J/cm3·K), which was calculated in our previous work [29]. The difference of heat capacities between MOF-74 and other MOFs was expected to be a result of the difference in their structures and components [3,30].

3.2. Adsorption Isotherms

Figure 5 depicts the adsorption isotherms of refrigerants including R1234yf, R1234ze, R134a and R32 in M-MOF-74 at different temperatures, in which symbols represent the results measured in our GCMC simulations. Among the four refrigerants investigated, R32 was demonstrated to have the highest adsorption capacity in the M-MOF-74 within the temperature range considered, evident in Figure 5c,g,k,o, followed by R134a. This was due to the molecule structure of refrigerants. As is shown in Figure 2, the R32 molecule has a structure with a smaller effective diameter. R1234yf and R1234ze own structures with large effective diameters. At the high pressures, when the adsorption approached saturation, the adsorbate molecules experienced a strong repulsive force from the tightly-packed neighbouring molecules and the host structure [31,32]. This repulsive interaction should be further enhanced for the molecules having a larger size, leading to reduction of the adsorption volume available in MOF-74 for R1234yf, R1234ze and R134a compared to R32. Consequently, the highest adsorption capacity was achieved for R32 in the pores of MOF-74. On the other hand, the Mg-MOF-74 adsorbed more refrigerants than other MOF-74s. This is because the ionic radius of Mg2+ is smaller than those of Co2+, Ni2+ and Zn2+, which results in a larger adsorption volume in the Mg-MOF-74. More importantly, the interactions between Mg2+ and adsorbate molecules were stronger than the counterpart interactions for Co2+, Ni2+ and Zn2+ ions [33]. The enthalpy of desorption calculated by GCMC simulations is shown in Appendix A, confirming this statement.

3.3. Thermal Energy Storage Properties of MOHCs

As previously mentioned, the thermal energy storage in MOHCs can be calculated according to Equation (1) or Equation (2). The relationships between thermal energy storage and temperature difference (the temperature of the cold source is 293 K) of the adsorbates (R1234yf, R1234ze, R134a, R32)/MOF-74 mixture at 3500 and 5000 kPa are shown in Figure 6 and Figure 7, respectively. The thermal energy storage capacity of adsorbates (R1234yf, R1234ze, R134a, R32) increased with the mass ratio of MOF-74 nanoparticles in the working fluid. The percentage enhancement of the thermal energy was calculated as [ΔhMOHCs(ΔT) − ΔhFluid(ΔT)]/ΔhFluid(ΔT), where ΔT is the temperature difference relative to the reference temperature of 293 K.
In the M-MOF-74 (M = Co, Mg, Ni) structures, the energy storage enhancement ratio of R1234yf was higher than that of the other three adsorbates. However, in Zn-MOF-74, R134a achieved the highest energy storage enhancement ratio among the four refrigerants considered. The energy storage enhancement ratio of R32 in Zn-MOF-74 was higher than that of the other three M-MOF-74 (M = Co, Mg, Ni) structures, while the energy storage enhancement ratios of R1234yf, R1234ze and R134a in Mg-MOF-74 were higher than those of the other three M-MOF-74 (M = Co, Ni, Zn). It is worth noting that the energy storage enhancement ratio at 5000 kPa is higher than that of 3500 kPa. However, the negative enhancement of thermal energy storage properties is found in a certain temperature range in Figure 6c,g,i,j,k,l,o and Figure 7c,k,l,o. This is because the enthalpy difference of the refrigerant was larger than the sum of the thermodynamic energy change of M-MOF-74 and desorption heat at the corresponding temperature range.

4. Conclusions

In this work, molecular simulations including MD and GCMC simulations were employed to investigate the energy storage by the adsorption of adsorbates (R1234yf, R1234ze, R134a, R32) in MOF-74. From the results, the following conclusions could be drawn.
Among the four different refrigerants, R32 has the largest adsorption capacity in M-MOF-74 materials, followed by R134a. This is due to the molecular structures of R32 and R134a being smaller than those of R1234yf and R1234ze. Mg-MOF-74 adsorbs more refrigerant than the other three MOFs, as the adsorption interaction between Mg-MOF-74 and refrigerants is stronger than the other MOFs.
The energy storage enhancement ratio of R1234yf with M-MOF-74 (M = Co, Mg, Ni) nanoparticles was higher than other refrigerants with M-MOF-74 (M = Co, Mg, Ni) nanoparticles. The energy storage enhancement ratio of R134a with Zn-MOF-74 nanoparticles was enhanced compared to other refrigerants with Zn-MOF-74 nanoparticles; while the energy storage enhancement ratio of R32 in Zn-MOF-74 was higher than those of the other three M-MOF-74 (M = Co, Mg, Ni) structures. The energy storage enhancement ratios of R1234yf, R1234ze and R134a in Mg-MOF-74 were higher than those of the other three M-MOF-74 (M = Co, Ni, Zn). Besides, the negative enhancement of thermal energy storage properties is found under certain conditions of MOHCs.

Author Contributions

Conceptualization, C.L. and Q.L.; Methodology, J.H. and Q.L.; Software, J.H.; Validation, J.H., C.L., L.L. and Q.L.; Formal Analysis, J.H. and Q.L.; Investigation, J.H. and Q.L.; Resources, Q.L.; Data Curation, J.H., C.L., L.L. and Q.L.; Writing-Original Draft Preparation, J.H., C.L., L.L. and Q.L.; Writing-Review & Editing, J.H., C.L., L.L. and Q.L.; Visualization, J.H.; Supervision, Q.L.; Project Administration, Q.L.; Funding Acquisition, Q.L.

Funding

This work is supported by the National Natural Science Foundation of China (No. 51506013), the Fundamental Research Funds for the Central Universities of China (Nos. 2018CDXYDL0001 and 2018CDXYHK0016) and the Graduate Research and Innovation Foundation of Chongqing, China (No. CYS18040).

Acknowledgments

The authors wish to acknowledge Muralidharan Paramsothy (the Guest Editor) and the anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ΔhMOHCsThe enthalpy change of MOHCs (kJ/kg)
ΔhFluidThe enthalpy change of pure organic fluid (kJ/kg)
ΔhdesorptionThe enthalpy of desorption (kJ/mol)
TTemperature (K)
CpThe heat capacity of MOFs (kJ/(mol·K))
xThe mass fraction of MOFs particle in MOHCs
ρDensity (g/cm3)

Appendix A. Enthalpy of Desorption Calculated by GCMC Simulations

Table A1. The adsorption heat calculated by grand canonical Monte Carlo (GCMC) simulations (kJ/mol).
Table A1. The adsorption heat calculated by grand canonical Monte Carlo (GCMC) simulations (kJ/mol).
MOHCs293 KSim313 KSim333 KSim353 KSim373 KSim393 KSim
R1234ze/Co128.19127.90127.03127.26127.10126.84
R1234yf/Co117.61116.50116.47115.96115.47114.57
R32/Co75.6775.2073.8675.3973.6674.39
R134a/Co124.93123.61123.31123.14123.03122.40
R1234ze/Mg279.73279.42278.58278.25277.57276.74
R1234yf/Mg235.72234.14237.97233.39230.97230.39
R32/Mg177.06176.96174.94173.94173.61171.57
R134a/Mg257.38256.32258.36259.78258.09256.47
R1234ze/Ni127.00126.54126.26125.93125.02124.87
R1234yf/Ni116.32115.67115.37115.01114.14113.70
R32/Ni76.0075.8975.6475.5974.4171.17
R134a/Ni123.42123.42123.09122.68122.48122.20
R1234ze/Zn128.18127.67127.28126.82126.11125.84
R1234yf/Zn117.87117.37117.21116.61115.99115.39
R32/Zn116.89116.58116.27114.68114.06115.62
R134a/Zn155.11153.91153.90153.87153.45152.49

References

  1. Chan, C.W.; Ling-Chin, J.; Roskilly, A.P. A review of chemical heat pumps, thermodynamic cycles and thermal energy storage technologies for low grade heat utilisation. Appl. Therm. Eng. 2013, 50, 1257–1273. [Google Scholar] [CrossRef]
  2. Yamamoto, T.; Furuhata, T.; Arai, N.; Mori, K. Design and testing of the Organic Rankine Cycle. Energy 2001, 26, 239–251. [Google Scholar] [CrossRef]
  3. McGrail, B.P.; Thallapally, P.K.; Blanchard, J.; Nune, S.K.; Jenks, J.J.; Dang, L.X. Metal-organic heat carrier nanofluids. Nano Energy 2013, 2, 845–855. [Google Scholar] [CrossRef]
  4. Li, Q.; Yu, Y.; Liu, Y.; Liu, C.; Lin, L. Thermal properties of the mixed n-octadecane/Cu nanoparticle nanofluids during phase transition: A molecular dynamics study. Materials 2017, 10, 38. [Google Scholar] [CrossRef] [PubMed]
  5. Xuan, Y.; Li, Q. Heat transfer enhancement of nanofluids. Int. J. Heat Fluid Flow 2000, 21, 58–64. [Google Scholar] [CrossRef]
  6. Celen, A.; Cebi, A.; Aktas, M.; Mahian, O.; Dalkilic, A.S.; Wongwises, S. A review of nanorefrigerants: Flow characteristics and applications. Int. J. Refrig. 2014, 44, 125–140. [Google Scholar] [CrossRef]
  7. Milanese, M.; Iacobazzi, F.; Colangelo, G.; Risi, A.D. An investigation of layering phenomenon at the liquid-solid interface in Cu and CuO based nanofluids. Int. J. Heat Mass Transf. 2016, 103, 564–571. [Google Scholar] [CrossRef]
  8. Iacobazzi, F.; Milanese, M.; Colangelo, G.; Lomascolo, M.; Risi, A.D. An explanation of the Al2O3 nanofluid thermal conductivity based on the phonon theory of liquid. Energy 2016, 116, 786–794. [Google Scholar] [CrossRef]
  9. James, S.L. Metal-organic frameworks. Chem. Soc. Rev. 2003, 32, 276–288. [Google Scholar] [CrossRef] [PubMed]
  10. Henninger, S.K.; Habib, H.A.; Janiak, C. MOFs as adsorbents for low temperature heating and cooling applications. J. Am. Chem. Soc. 2009, 131, 2776–2777. [Google Scholar] [CrossRef] [PubMed]
  11. Murray, L.J.; Dincă, M.; Long, J.R. Hydrogen storage in metal–organic frameworks. Chem. Soc. Rev. 2009, 38, 1294–1314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Li, J.R.; Sculley, J.; Zhou, H.C. Metal-organic frameworks for separations. Chem. Rev. 2012, 112, 869–932. [Google Scholar] [CrossRef] [PubMed]
  13. Getman, R.B.; Bae, Y.S.; Wilmer, C.E.; Snurr, R.Q. Review and analysis of molecular simulations of methane, hydrogen, and acetylene storage in metal-organic frameworks. Chem. Rev. 2012, 112, 703–723. [Google Scholar] [CrossRef] [PubMed]
  14. Cavallini, A.; Censi, G.; Del Col, D.; Doretti, L.; Longo, G.A.; Rossetto, L. Experimental investigation on condensation heat transfer and pressure drop of new HFC refrigerants (R134a, R125, R32, R410A, R236ea) in a horizontal smooth tube. Int. J. Refrig. 2001, 24, 73–87. [Google Scholar] [CrossRef]
  15. Steven Brown, J.; Yana-Motta, S.F.; Domanski, P.A. Comparitive analysis of an automotive air conditioning systems operating with CO2 and R134a. Int. J. Refrig. 2002, 25, 19–32. [Google Scholar] [CrossRef]
  16. Fouad, W.A.; Vega, L.F. Next generation of low global warming potential refrigerants: Thermodynamic properties molecular modeling. AIChE J. 2018, 64, 250–262. [Google Scholar] [CrossRef]
  17. Molés, F.; Navarro-Esbrí, J.; Peris, B.; Mota-Babiloni, A.; Barragán-Cervera, Á. Theoretical energy performance evaluation of different single stage vapour compression refrigeration configurations using R1234yf and R1234ze(E) as working fluids. Int. J. Refrig. 2014, 44, 141–150. [Google Scholar] [CrossRef] [Green Version]
  18. Yataganbaba, A.; Kilicarslan, A.; Kurtbaş, İ. Exergy analysis of R1234yf and R1234ze as R134a replacements in a two evaporator vapour compression refrigeration system. Int. J. Refrig. 2015, 60, 26–37. [Google Scholar] [CrossRef]
  19. Janković, Z.; Sieres Atienza, J.; Martínez Suárez, J.A. Thermodynamic and heat transfer analyses for R1234yf and R1234ze(E) as drop-in replacements for R134a in a small power refrigerating system. Appl. Therm. Eng. 2015, 80, 42–54. [Google Scholar] [CrossRef]
  20. Deethayat, T.; Kiatsiriroat, T.; Thawonngamyingsakul, C. Performance analysis of an organic Rankine cycle with internal heat exchanger having zeotropic working fluid. Case Stud. Therm. Eng. 2015, 6, 155–161. [Google Scholar] [CrossRef]
  21. Barse, K.A.; Mann, M.D. Maximizing ORC performance with optimal match of working fluid with system design. Appl. Therm. Eng. 2016, 100, 11–19. [Google Scholar] [CrossRef] [Green Version]
  22. Ding, Y.; Liu, C.; Zhang, C.; Xu, X.; Li, Q.; Mao, L. Exergoenvironmental model of Organic Rankine Cycle system including the manufacture and leakage of working fluid. Energy 2018, 145, 52–64. [Google Scholar] [CrossRef]
  23. Li, Q.; Xiao, Y.; Shi, X.; Song, S. Rapid evaporation of water on graphene/graphene-oxide: A molecular dynamics study. Nanomaterials 2017, 7, 265. [Google Scholar] [CrossRef] [PubMed]
  24. National Institute of Standards and Technology (NIST). Thermophysical Properties of Fluid Systems. Available online: http://webbook.nist.gov/chemistry/fluid/ (accessed on 18 May 2018).
  25. BIOVIA. Materials Studio; Accelrys Software Inc.: San Diego, CA, USA, 2010. [Google Scholar]
  26. Sun, H. COMPASS: An ab initio force-field optimized for condensed-phase applications overview with details on alkane and benzene compounds. J. Phys. Chem. B 1998, 102, 7338–7364. [Google Scholar] [CrossRef]
  27. Ewald, P.P. Die Berechnung optischer und elektrostatischer Gitterpotentiale. Ann. Phys. 1921, 369, 253–287. [Google Scholar] [CrossRef]
  28. Britt, D.; Furukawa, H.; Wang, B.; Grant Glover, T.; Yaghi, O.M. Highly efficient separation of carbon dioxide by a metal-organic framework replete with open metal sites. Proc. Natl. Acad. Sci. USA 2009, 106, 20637–20640. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Hu, J.; Liu, C.; Li, Q.; Shi, X. Molecular simulation of thermal energy storage of mixed CO2/IRMOF-1 nanoparticle nanofluid. Int. J. Heat Mass Transf. 2018, 125, 1345–1348. [Google Scholar] [CrossRef]
  30. Mu, B.; Walton, K.S. Thermal analysis and heat capacity study of metal–organic frameworks. J. Phys. Chem. C 2011, 115, 22748–22754. [Google Scholar] [CrossRef]
  31. Liu, L.; Bhatia, S.K. Molecular simulation of CO2 adsorption in the presence of water in single-walled carbon nanotubes. J. Phys. Chem. C 2013, 117, 13479–13491. [Google Scholar] [CrossRef]
  32. Liu, L.; Nicholson, D.; Bhatia, S.K. Adsorption of CH4 and CH4/CO2 mixtures in carbon nanotubes and disordered carbons: A molecular simulation study. Chem. Eng. Sci. 2015, 121, 268–278. [Google Scholar] [CrossRef]
  33. Pham, T.; Forrest, K.A.; Banerjee, R.; Orcajo, G.; Eckert, J.; Space, B. Understanding the H2 sorption trends in the M-MOF-74 series (M = Mg, Ni, Co, Zn). J. Phys. Chem. C 2014, 119, 1078–1090. [Google Scholar] [CrossRef]
Figure 1. Schematic view of the 1 × 1 × 1 unit cell of metal organic framework MOF-74 (M = Ni).
Figure 1. Schematic view of the 1 × 1 × 1 unit cell of metal organic framework MOF-74 (M = Ni).
Materials 11 01164 g001
Figure 2. Atomic configurations of R1234yf (a); R1234ze (b); R134a (c) and R32 (d).
Figure 2. Atomic configurations of R1234yf (a); R1234ze (b); R134a (c) and R32 (d).
Materials 11 01164 g002
Figure 3. Adsorption sites of R1234ze in Co-MOF-74 at 353 K and 5 MPa.
Figure 3. Adsorption sites of R1234ze in Co-MOF-74 at 353 K and 5 MPa.
Materials 11 01164 g003
Figure 4. Variation of the relative change of the thermodynamic energy of M-MOF-74 with temperature, with M denoting Co, Mg, Ni and Zn.
Figure 4. Variation of the relative change of the thermodynamic energy of M-MOF-74 with temperature, with M denoting Co, Mg, Ni and Zn.
Materials 11 01164 g004
Figure 5. Adsorption isotherms of refrigerants in M-MOF-74 at different temperatures. (a) R1234ze/Co-MOF-74; (b) R1234yf/Co-MOF-74; (c) R32/Co-MOF-74; (d) R134a/Co-MOF-74; (e) R1234ze/Mg-MOF-74; (f) R1234yf/Mg-MOF-74; (g) R32/Mg-MOF-74; (h) R134a/Mg-MOF-74; (i) R1234ze/Ni-MOF-74; (j) R1234yf/Ni-MOF-74; (k) R32/Ni-MOF-74; (l) R134a/Ni-MOF-74; (m) R1234ze/Zn-MOF-74; (n) R1234yf/Zn-MOF-74; (o) R32/Zn-MOF-74; (p) R134a/Zn-MOF-74.
Figure 5. Adsorption isotherms of refrigerants in M-MOF-74 at different temperatures. (a) R1234ze/Co-MOF-74; (b) R1234yf/Co-MOF-74; (c) R32/Co-MOF-74; (d) R134a/Co-MOF-74; (e) R1234ze/Mg-MOF-74; (f) R1234yf/Mg-MOF-74; (g) R32/Mg-MOF-74; (h) R134a/Mg-MOF-74; (i) R1234ze/Ni-MOF-74; (j) R1234yf/Ni-MOF-74; (k) R32/Ni-MOF-74; (l) R134a/Ni-MOF-74; (m) R1234ze/Zn-MOF-74; (n) R1234yf/Zn-MOF-74; (o) R32/Zn-MOF-74; (p) R134a/Zn-MOF-74.
Materials 11 01164 g005aMaterials 11 01164 g005bMaterials 11 01164 g005c
Figure 6. The relationship between the thermal energy storage and temperature difference of the adsorbates (R1234yf, R1234ze, R134a, R32)/MOF-74 mixture at 3500 kPa. (a) R1234ze/Co-MOF-74; (b) R1234yf/Co-MOF-74; (c) R32/Co-MOF-74; (d) R134a/Co-MOF-74; (e) R1234ze/Mg-MOF-74; (f) R1234yf/Mg-MOF-74; (g) R32/Mg-MOF-74; (h) R134a/Mg-MOF-74; (i) R1234ze/Ni-MOF-74; (j) R1234yf/Ni-MOF-74; (k) R32/Ni-MOF-74; (l) R134a/Ni-MOF-74; (m) R1234ze/Zn-MOF-74; (n) R1234yf/Zn-MOF-74; (o) R32/Zn-MOF-74; (p) R134a/Zn-MOF-74.
Figure 6. The relationship between the thermal energy storage and temperature difference of the adsorbates (R1234yf, R1234ze, R134a, R32)/MOF-74 mixture at 3500 kPa. (a) R1234ze/Co-MOF-74; (b) R1234yf/Co-MOF-74; (c) R32/Co-MOF-74; (d) R134a/Co-MOF-74; (e) R1234ze/Mg-MOF-74; (f) R1234yf/Mg-MOF-74; (g) R32/Mg-MOF-74; (h) R134a/Mg-MOF-74; (i) R1234ze/Ni-MOF-74; (j) R1234yf/Ni-MOF-74; (k) R32/Ni-MOF-74; (l) R134a/Ni-MOF-74; (m) R1234ze/Zn-MOF-74; (n) R1234yf/Zn-MOF-74; (o) R32/Zn-MOF-74; (p) R134a/Zn-MOF-74.
Materials 11 01164 g006aMaterials 11 01164 g006bMaterials 11 01164 g006c
Figure 7. The relationship between the thermal energy storage and temperature difference of the adsorbates (R1234yf, R1234ze, R134a, R32)/MOF-74 mixture at 5000 kPa. (a) R1234ze/Co-MOF-74; (b) R1234yf/Co-MOF-74; (c) R32/Co-MOF-74; (d) R134a/Co-MOF-74; (e) R1234ze/Mg-MOF-74; (f) R1234yf/Mg-MOF-74; (g) R32/Mg-MOF-74; (h) R134a/Mg-MOF-74; (i) R1234ze/Ni-MOF-74; (j) R1234yf/Ni-MOF-74; (k) R32/Ni-MOF-74; (l) R134a/Ni-MOF-74; (m) R1234ze/Zn-MOF-74; (n) R1234yf/Zn-MOF-74; (o) R32/Zn-MOF-74; (p) R134a/Zn-MOF-74.
Figure 7. The relationship between the thermal energy storage and temperature difference of the adsorbates (R1234yf, R1234ze, R134a, R32)/MOF-74 mixture at 5000 kPa. (a) R1234ze/Co-MOF-74; (b) R1234yf/Co-MOF-74; (c) R32/Co-MOF-74; (d) R134a/Co-MOF-74; (e) R1234ze/Mg-MOF-74; (f) R1234yf/Mg-MOF-74; (g) R32/Mg-MOF-74; (h) R134a/Mg-MOF-74; (i) R1234ze/Ni-MOF-74; (j) R1234yf/Ni-MOF-74; (k) R32/Ni-MOF-74; (l) R134a/Ni-MOF-74; (m) R1234ze/Zn-MOF-74; (n) R1234yf/Zn-MOF-74; (o) R32/Zn-MOF-74; (p) R134a/Zn-MOF-74.
Materials 11 01164 g007aMaterials 11 01164 g007bMaterials 11 01164 g007c

Share and Cite

MDPI and ACS Style

Hu, J.; Liu, C.; Liu, L.; Li, Q. Thermal Energy Storage of R1234yf, R1234ze, R134a and R32/MOF-74 Nanofluids: A Molecular Simulation Study. Materials 2018, 11, 1164. https://doi.org/10.3390/ma11071164

AMA Style

Hu J, Liu C, Liu L, Li Q. Thermal Energy Storage of R1234yf, R1234ze, R134a and R32/MOF-74 Nanofluids: A Molecular Simulation Study. Materials. 2018; 11(7):1164. https://doi.org/10.3390/ma11071164

Chicago/Turabian Style

Hu, Jieyao, Chao Liu, Lang Liu, and Qibin Li. 2018. "Thermal Energy Storage of R1234yf, R1234ze, R134a and R32/MOF-74 Nanofluids: A Molecular Simulation Study" Materials 11, no. 7: 1164. https://doi.org/10.3390/ma11071164

APA Style

Hu, J., Liu, C., Liu, L., & Li, Q. (2018). Thermal Energy Storage of R1234yf, R1234ze, R134a and R32/MOF-74 Nanofluids: A Molecular Simulation Study. Materials, 11(7), 1164. https://doi.org/10.3390/ma11071164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop