Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Inducing Crystallinity of Metal Thin Films with Weak Magnetic Fields without Thermal Annealing
Previous Article in Journal
Tracing Phase Transformation and Lattice Evolution in a TRIP Sheet Steel under High-Temperature Annealing by Real-Time In Situ Neutron Diffraction
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Dynamics Investigation of Graphene Nanoplate Diffusion Behavior in Poly-α-Olefin Lubricating Oil

1
MIIT Key Laboratory of Aerospace Bearing Technology and Equipment, Harbin Institute of Technology, Harbin 150001, China
2
School of Mechatronics Engineering, Harbin Institute of Technology, Harbin 150001, China
3
Department of Mechanical Engineering, Northwestern University, Evanston, IL 60208, USA
*
Author to whom correspondence should be addressed.
Crystals 2018, 8(9), 361; https://doi.org/10.3390/cryst8090361
Submission received: 27 June 2018 / Revised: 27 August 2018 / Accepted: 9 September 2018 / Published: 13 September 2018
Figure 1
<p>Snapshot of EMD simulation of (<b>a</b>) nanofluid cell with (<b>b</b>) graphene and (<b>c</b>) 1-decene trimer molecules.</p> ">
Figure 2
<p>Graphene plates with various functional groups: (<b>a</b>) pure graphene (<b>b</b>) GO (<b>c</b>) MG and (<b>d</b>) AG.</p> ">
Figure 3
<p>Equilibrium curve for energy of the graphene-PAO nanofluids.</p> ">
Figure 4
<p>(<b>a</b>) MSD curves and (<b>b</b>) diffusion coefficients of graphene with various sizes at 300 K and 0.1 MPa.</p> ">
Figure 5
<p>Diffusion coefficients of graphene with various functional groups calculated under NVT ensemble at 300 K and 0.1 MPa.</p> ">
Figure 6
<p>Photographs of blended lubricating oil containing (<b>a</b>) graphene obtained 3 h, 10 h, 100 h after sonication and (<b>b</b>) graphene oxide 0 week, 1 week, 2 weeks after sonication at room temperature.</p> ">
Figure 7
<p>Radial distribution function for graphene-PAO molecules at 300 K and 0.1 MPa from EMD simulations.</p> ">
Figure 8
<p>Diffusion coefficients of graphene containing various number of functional groups calculated at NVT ensemble at 300 K and 1 atm.</p> ">
Figure 9
<p>Diffusion coefficients of graphene containing two functional groups with respect to increasing temperature at 0.1 MPa.</p> ">
Figure 10
<p>Diffusion coefficients of graphene containing two functional groups with respect to increasing pressures at 300 K.</p> ">
Versions Notes

Abstract

:
Graphene as a type of novel additive significantly enhanced the tribological performance of blended lubricating oil. However, the dispersibility of graphene with long-term stability in lubricating oil is still a challenge. Chemical modification for graphene, rather than using surfactants, provided a better method to improve the dispersibility of graphene in lubricants. In this study, the equilibrium molecular dynamics (EMD) simulations were carried out to investigate the diffusion behavior of graphene nanoplates in poly-α-olefin (PAO) lubricating oil. The effects of graphene-size, edge-functionalization, temperature, and pressure on the diffusion coefficient were studied. In order to understand the influence of edge-functionalization, three different functional groups were grafted to the edge of graphene nanoplates: COOH, COON(CH3)2, CONH(CH2)8CH3 (termed GO, MG, and AG, respectively). The EMD simulations results demonstrated that the relationships between diffusion coefficient and graphene-size and number of functional groups were linear while the temperature and pressure had a nonlinear influence on the diffusion coefficient. It was found that the larger dimension and more functional groups provided the lower diffusion coefficient. AG with eight CONH(CH2)8CH3 groups exhibited the lowest diffusion coefficient. Furthermore, the experimental results and radial distribution function for graphene-PAO illustrated that the diffusion coefficient reflected the dispersibility of nanoparticles in nanofluids to some degree. To our best knowledge, this study is the first time the diffusion behavior of graphene in PAO lubricating oil was investigated using EMD simulations.

1. Introduction

Mechanical system failure can originate from friction and wear. Recently, research on the aerospace application of tribology have been given attention because more and more non-renewable energy resources have been consumed to overcome friction [1,2]. To prolong service life, lowering the friction energy consumption by increasing the anti-wearability has been an ongoing trend [3]. In this case, researchers believed that tribological performance of lubricating oils could be significantly enhanced by means of adding various nanoparticles, such as SiO2 [4], Al2O3 [5], MoS2 [6,7], WS2 [8], zinc dialkyl-dithio-phosphate (ZDDP) [9,10], cyclopropane-carboxylic acid (CPCa) [11], and carbon nanomaterials [12,13,14,15].
Over the last decades, graphene has attracted significant attention owing to its remarkable physical and chemical properties, including excellent thermal conductivity, high strength, and superior chemical inertness [16,17]. Moreover, as a novel lubricating oil additive with good environmental compatibility, graphene has proved to enhance dramatically the tribological properties of the lubricating oils [18,19]. The research studied by Eswaraish et al. [20] demonstrated that extreme pressure capacity, frictional coefficient and wear resistance were respectively improved by 40%, 80%, and 33% after adding graphene into the engine oil. Wu et al. [21] investigated the tribological performance of 4010 aviation lubricants containing graphene nanoplates. Compared to base oil, the friction coefficient and wear scar diameter were reduced by 27% and 43% with the addition of graphene. In addition, Mao et al. [22] carried out experiments to study the influence of the micromorphology of graphene on lubrication performance. The results illustrated that a graphene lubricant with regular edges has the best lubrication properties as friction reduction and anti-wear ability were improved by 27.9% and 14.1% of those of base oil.
However, as observed in our previous study [21], graphene was easy to form aggregates in lubricating oils, which was a challenge for the further applications of graphene as additives. In this case, Zheng et al. [23] added surfactant (span-80, C24H44O6) into PAO4 lubricating oil to improve the dispersion stability of graphene. In other ways, the methods of chemical modification of graphene have been widely used to obtain stable blended nanofluids [24,25]. Wu et al. [26] attached red phosphorus to graphene plates using a mechanical synthesis method and the mixed lubricating oil was uniformly dispersed for two months. Zhang et al. [27] developed a novel boehmite/graphene oxide nanohybrid and was observed to keep stable dispersion in lubricant base oil (VHVI8) at least 48 h.
Previous research on the dispersion stability of graphene focused on macroscopic experiments. However, limited efforts have been devoted to the investigations on the dispersion behavior of graphene from the microscopic perspective. In order to understand the microscopic behavior and obtain details that are inaccessible from experiments, molecular dynamics (MD) simulations were selected by many researchers [28,29,30,31,32,33,34]. Saathoff et al. [35] explored how three different edge terminations affected the dispersion stability of graphene in n-methyl-pyrrolidone (NMP) via MD simulations. Konatham et al. [36] reported that functional groups grafted on the graphene plates were effective for the stable dispersion of graphene in liquid organic alkanes. Chen et al. [37] investigated the aggregation rate of two parallel graphene in dodecane. The results illustrated that graphene aggregated when parts of two graphene are in contact with each other. However, only a few efforts have been devoted to investigating the dispersibility of graphene in lubricating oils. Its detailed process from a microscopic perspective is still unclear, which limited the wide application of graphene as additives in the tribology area. In this study, the dispersibility of graphene was evaluated using a diffusion coefficient, similar to the method provided by Cha et al. [38]. The EMD simulations were employed to calculate the diffusion coefficient of graphene-in-PAO4 systems to reveal the effects of graphene-size, functional groups, temperature, and pressure on dispersibility of graphene in PAO lubricant.

2. Materials and Simulation Method

In this work, all models were established using Materials Studio 8.0 (Accelrys Software Inc., San Diego, CA, USA). The consistent valence forcefield (CVFF) implemented in the software package Materials Studio was used for both PAO molecules and graphene. For the blended nanofluids, we used 1-decene trimer, the major component of PAO 4 [39], as the base lubricating oil. As shown in Figure 1, the graphene-in-PAO4 system consisted of one cubic cell of 6.12 nm in length (containing 20,300, 20,306, 20,322, and 20,362 atoms for using graphene (G), GO, MG and AG with side length of 1 nm, respectively) and contained a single graphene and 220 1-decene trimer molecules. Periodic boundary conditions were employed in the three directions. Four different graphene sizes with 1, 2, 3 and 4 nm were selected in the extraction, which resulted in the 0.57, 1.84, 3.15, and 7.54 wt.% concentrations of graphene in nanofluids. Graphene (G), along with its modified alternatives GO (G-COOH) [40], MG (G-COON(CH3)2) [41], and AG (G-CONH(CH2)8CH3) [42] were employed to explore their diffusion behavior in PAO 4 lubricating oil, as displayed in Figure 2. Graphenes with 4 nm side length were used for investigating the effect of the type and number of functional groups, temperature and pressure on diffusion coefficient. The graphene was gradually modified to study the effect of number of functional groups (2, 4, 6, and 8 corresponding to 2.9, 5.8, 8.7, and 11.6%, respectively).
All the EMD simulations were performed via Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [43] in this study. N (number of molecules), P (pressure), and T (temperature) ensemble for all temperatures were selected with a Nose-Hoover [44,45] style to retain the system temperature. The non-bonded interactions contained both electrostatic and van der Waals (vdW) interactions in the EMD simulations, where the time step length was 1 fs. The vdW interaction with a 10 Å cut off was explained by the standard 12/6 Lennard-Jones potential.
E r = 4 ε [ ( σ r ) 12 ( σ r ) 6 ]   r < r c ,
where i is the depth of potential well, r represents the distance between two atoms, and σ is the distance where the interaction between atoms is equal to zero. The short range electrostatic interactions up to a distance of 10 Å were accounted for by applying the Coulombic law
E q = C q i q j σ r   r < r c ,
Here C is the energy conversion constant, and q i , q j represent the charges of two atoms, respectively. The long range electrostatic interactions beyond 10 Å were processed using the particle-mesh Ewald (PME) method with a precision of 10−5. In order to investigate the influence of temperature on diffusion behavior of graphene in PAO 4 lubricant, the simulations were performed over a temperature range of 300–400 K with a degree interval of 20 K. Energy minimization was conducted before the EMD simulations began. The system was equilibrated for 10 ns at an N (number of molecules), P (pressure), and T (temperature) ensemble until its energy fluctuations become stable, as shown in Figure 3. This provided the initial point for all simulations. After equilibration, the production steps of 10 ns were conducted at an N (number of molecules), V (volume), and T (temperature) ensemble and the time step was 1 fs.

3. Results and Discussion

Graphene was found to easily aggregate and form sediment in lubricating oil. This property limited the application of graphene as lubricants additives, although it was proved to significantly enhance the tribological performance. In order to improve the dispersibility of graphene in PAO lubricating oil, it is necessary to quantitatively describe the influences of graphene-size, edge-functionalization, temperature and pressure on the dynamic properties of graphene in lubricating oil. The diffusion coefficient of graphene was calculated to describe its dynamic behavior in the nano range. The mean square displacements (MSDs) was another critical factor to describe the microscopic dynamics of nanoparticles, which reflects the rheological properties of nanofluids. Therefore, the motion of graphenes was analyzed by calculating the MSDs and diffusion coefficient ( D ) using Equations (3) and (4) [46,47].
M S D ( t ) = 1 N i = 1 N | r i ( t ) r i ( t 0 ) | 2 ,
D = 1 6 lim t ( d d t M S D ( t ) ) ,
Here N is the total number of atoms, r i ( t ) is the position of atom i at time t and r i ( t 0 ) represents the original location of atom. The appropriate size of additives corresponded to its excellent dispersion stability. In this study, different dimensions of graphene with 1, 2, 3, and 4 nm side length were used in the simulations to study the effect of graphene-size on the diffusion coefficient.
Figure 4 showed the MSDs and diffusion coefficients of single graphene in the PAO 4 lubricant. The fluctuations of MSDs were very small, and indicated that the equilibrium time of 10 ns was a suitable time period, as shown in Figure 4a. The results showed that the smaller graphene moves more than the larger one. Figure 4b exhibited that the diffusion coefficient decreased with increasing graphene-size. Compared to 1 nm graphene-size, the diffusion coefficient of 4 nm was reduced by 88.18%. Similar results were also found in another study, which illustrated that the larger nanoparticles were more stable than that the smaller ones [38].
Figure 5 exhibited the diffusion coefficients of four types of graphene containing two functional groups, which were shown in Figure 2, the original graphene (G), GO (G-COOH), MG (G-COON(CH3)2), and AG (G-CONH(CH2)8CH3) at 300 K and 0.1 MPa. It was clear that the diffusion coefficients decreased as the graphene-size became larger. When the graphene-size increased to 4 nm from 1 nm, the diffusion coefficients for GO, MG and AG were reduced by 83.75%, 84.91% and 84.04%, respectively. In addition, it was found that functional groups had a significant effect on the diffusion coefficient. Compared to G, the diffusion coefficients of GO, MG and AG were reduced by 48.18%, 70.91% and 91.45%, respectively when the graphene-size was 1 nm. This result demonstrated that AG, showing the smallest diffusion coefficient, was most stable in base oil. Similarly, Cha et al. [36] found that low speeds of graphene in liquids had less opportunity to aggregate. The diffusion coefficient of GO was smaller than that of G, which estimated that GO had the better dispersibility in lubricating oil and is consistent with the experimental result. As depicted in Figure 6, GO kept uniformly dispersed in lubricating oil for one week while the G exhibited complete sedimentation after setting for 100 h. In order to further understand the relationship between the diffusion coefficient and dispersibility of graphene in lubricating oil, the radial distribution functions (RDFs) for graphene-PAO molecules were computed at 300 K and 0.1 MPa.
In order to investigate the interaction between graphenes and PAO molecules, the RDFs between the centers of mass of graphenes and PAO molecules were calculated using graphenes with 4 nm side length. The RDFs correspond to the data obtained during 10 ns after equilibrium under NPT ensemble. As depicted in Figure 7, there are five sets of peaks when the distances between graphenes and PAO molecules are 1.325, 2.325, 2.575, 3.525 and 4.024 Å. The intensity becomes lower with the longer distance between graphene and PAO molecules, and the first set of peaks are the most important to evaluate the association between graphenes and PAO molecules. The higher peak at shorter separation was representative of the stronger association between graphene and PAO molecules; namely, stronger association represented better dispersibility [36]. It was clear that AG exhibited the best dispersibility while G showed the weakest. Four types of graphene, G, GO, MG and AG were ranked according to their dispersibility in lubricating oil as: G < GO < MG < AG.
The influence of numbers of functional groups on diffusion coefficient was investigated using graphene with 1 nm of side length. The numbers of 2, 4, 6 and 8 of functional groups were grafted to the edge of graphene and the original graphene provided a baseline. Due to limitation of our computation ability, graphenes with side length of 4 nm were not selected in this study, although larger ones could reflect the dynamic properties better. Figure 8 shows that the diffusion coefficients of GO, MG, and AG linearly decreased with respect to increasing number of functional groups. In comparison with the original graphene, the reduction of diffusion coefficients for various graphenes containing different number of functional groups are listed in Table 1. The diffusion coefficients of GO, MG and AG containing eight functional groups were reduced by 87.27%, 96.36%, and 96.73%, respectively. The results illustrated that more functional groups resulted in better dispersibility for graphenes in lubricating oil. The reason might be attributed to the fact that graphene with more functional groups provided more intensive association between graphenes and PAO molecules, which allowed graphene to stably dispersed in oil.
Different graphenes with side length of 4 nm were employed to investigate the effect of temperature and pressure on the diffusion coefficient of the four types of graphenes in oil. In our previous tribology tests [21], we found that the flash temperatures were higher than 300 but lower than 345 K when tribo-pairs were lubricated by oil containing graphene. In this case, temperatures ranging from 300 to 400 K were selected to study the diffusion coefficients of graphenes. Figure 9 showed that diffusion coefficients of graphenes increased with respect to increasing temperature, and AG exhibited the smallest diffusion coefficients at all simulated temperatures. In addition, the effect of temperature on diffusion coefficient was nonlinear, which was different from that of graphene-size and number of functional groups. The diffusion coefficients for G, GO, MG, and AG were improved by 310%, 326%, 383% and 422% when the temperature increased to 400 K from 300 K, which illustrated that AG was the most easily affected by temperature, and then MG, GO, and G. The results could be explained by functional groups being more easily affected by temperature than original graphene with stable structure.
The relationship between diffusion coefficient of graphenes and pressure was investigated under pressures ranging from 0.1 to 500 MPa, as shown in Figure 10. The effects of pressure on diffusion coefficients of graphenes were nonlinear, which was similar to temperature. However, the diffusion coefficients of graphenes decreased as the high pressure was applied in simulations. A similar result was obtained by Michalis et al. [48]. Moreover, the original graphene exhibited the largest diffusion coefficients under all simulated pressures while AG the smallest. When the pressure was raised up to 500 MPa the diffusion coefficients for G, GO, MG, and AG were reduced by 98.7%, 97.4%, 95.1%, and 93.6%, respectively.

4. Conclusions

The EMD simulations were conducted to investigate the diffusion coefficients of four types of graphenes in graphene-PAO systems, which offered a comprehensive microscopic insight into the dispersion behavior of graphene nanoplates in lubricating oil. The relationship between dispersibility and diffusion coefficient was established and validated by experimental results, radial distribution function for graphene-PAO and the results presented by Cha et al. [38] The simulation results demonstrated that the diffusion coefficient linearly decreased with the larger size and more functional groups. In addition, the diffusion coefficient nonlinearly increased with respect to increasing temperature while nonlinearly decreased as the pressure increased. For AG, it was improved by 1052% when the temperature increased to 400 K from 300 K. Furthermore, four types of graphene were ranked according to their dispersibility in lubricating oil as: G < GO < MG < AG. This work provided a new method to evaluate dispersibility of nanoparticles in nanofluids. To our best knowledge, this study was the first time the diffusion behavior of graphene in PAO lubricating oil was investigated using EMD simulations.

Author Contributions

L.G. design the simulation; L.W. performed the simulation and wrote the paper; L.W. and L.M.K. analyzed the data. L.M.K. and B.S. revised the paper.

Funding

This work was supported by the Natural Science Foundation of China (Grant No. 51805109) and the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (Grant No. 51521003). The first author thanks financial support from China Scholarship Council (CSC, No. 201706120143).

Acknowledgments

The authors would like to thank Q. Jane Wang, Mechanical Engineering, Northwestern University, for her helpful discussion and suggestions. The authors also want to thank reviewers for their helpful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jost, H.P. Tribology—Origin and future. Wear 1990, 136, 1–17. [Google Scholar] [CrossRef]
  2. Holmberg, K.; Erdemir, A. Influence of tribology on global energy consumption, costs and emissions. Friction 2017, 5, 263–284. [Google Scholar] [CrossRef] [Green Version]
  3. Jost, H.P.; Schofield, J. The practical aspects of friction, lubrication and wear. Proc. Inst. Mech. Eng. 1981, 195, 151–173. [Google Scholar] [CrossRef]
  4. Jiao, D.; Zheng, S.; Wang, Y.; Guan, R.; Cao, B. The tribology properties of alumina/silica composite nanoparticles as lubricant additives. Appl. Surf. Sci. 2011, 257, 5720–5725. [Google Scholar] [CrossRef]
  5. Etemadi, H.; Shojaei, A.; Jahanmard, P. Effect of alumina nanoparticle on tribological performance of automotive brake friction materials. J. Reinf. Plast. Compos. 2014, 33, 166–178. [Google Scholar] [CrossRef]
  6. Wu, H.X.; Qin, L.G.; Dong, G.N.; Hua, M.; Yang, S.C.; Zhang, J.F. An investigation on the lubrication mechanism of MoS2 nano sheet in point contact: The manner of particle entering the contact area. Tribol. Int. 2017, 107, 48–55. [Google Scholar] [CrossRef]
  7. Aralihalli, S.; Biswas, S.K. Grafting of dispersants on MoS2 in base oil lubrication of steel. Tribol. Lett. 2013, 49, 61–76. [Google Scholar] [CrossRef]
  8. Pottuz, L.J.; Dassenoy, F.; Belin, M.; Vacher, B.; Martin, J.; Fleischer, N. Ultralow-friction and wear properties of IF-WS 2 under boundary lubrication. Tribol. Lett. 2005, 18, 477–485. [Google Scholar] [CrossRef]
  9. Aldana, P.U.; Dassenoy, F.; Vacher, B.; Mogne, T.L.; Thiebaut, B. WS2 nanoparticles anti-wear and friction reducing properties on rough surfaces in the presence of ZDDP additive. Tribol. Int. 2016, 102, 213–221. [Google Scholar] [CrossRef]
  10. Sharma, V.; Timmons, R.; Erdemir, A.; Aswath, P.B. Plasma-functionalized polytetrafluoroethylene nanoparticles for improved wear in lubricated contact. ACS Appl. Mater. Interfaces 2017, 9, 25631–25641. [Google Scholar] [CrossRef] [PubMed]
  11. Johnson, B.; Wu, H.X.; Desanker, M.; Pickens, D.; Chung, Y.W.; Wang, Q.J. Direct formation of lubricious and wear-protective carbon films from phosphorus- and sulfur-free oil-soluble additives. Tribol. Lett. 2018, 66, 2. [Google Scholar] [CrossRef]
  12. Wang, Z.Q.; Ren, R.R.; Song, H.J.; Jia, X.H. Improved tribological properties of the synthesized copper/carbon nanotube nanocomposites for rapeseed oil-based additives. Appl. Surf. Sci. 2018, 428, 630–639. [Google Scholar] [CrossRef]
  13. Alazemi, A.A.; Etacheri, V.; Dysart, A.D.; Stacke, L.E.; Pol, V.G.; Sadeghi, F. ultrasmooth submicrometer carbon spheres as lubricant additives for friction and wear reduction. ACS Appl. Mater. Interfaces 2015, 7, 5514–5521. [Google Scholar] [CrossRef] [PubMed]
  14. Zhai, W.Z.; Srikanth, N.; Kong, L.B.; Zhou, K. Carbon nanomaterials in tribology. Carbon 2017, 119, 150–171. [Google Scholar] [CrossRef]
  15. Shahnazar, S.; Bagheri, S.; Hamid, S.B.A. Enhancing lubricant properties by nanoparticle Additives. Int. J. Hydrogen Energy 2016, 41, 3153–3170. [Google Scholar] [CrossRef]
  16. Yu, W.M.; David, M.; Routbort, J.L.; Choi, S. Review and comparison of nanofluid thermal conductivity and heat transfer enhancements. Heat. Transf. Eng. 2008, 29, 432–460. [Google Scholar] [CrossRef]
  17. Lee, C.; Wei, X.D.; Kysar, J.W.; Hone, J. Measurement of the elastic properties and intrinsic strength of monolayer graphene. Science 2008, 321, 385–388. [Google Scholar] [CrossRef] [PubMed]
  18. Berman, D.; Erdemir, A.; Sumant, A.V. Graphene: A new emerging lubricant. Mater. Today 2014, 17, 31–42. [Google Scholar] [CrossRef]
  19. Rasheed, A.K.; Khalid, M.; Rashmi, W.; Guupta, T.C.S.M.; Chan, A. Graphene based nanofluids and nanolubricants—Review of recent developments. Renew. Sustain. Energy Rev. 2016, 63, 346–362. [Google Scholar] [CrossRef]
  20. Eswaraiah, V.; Sankaranarayanan, V.; Ramaprabhu, S. Graphene-based engine oil nanofluids for tribological applications. ACS Appl. Mater. Interfaces 2011, 3, 4221–4227. [Google Scholar] [CrossRef] [PubMed]
  21. Wu, L.P.; Gu, L.; Xie, Z.J.; Zhang, C.W.; Song, B.Y. Improved tribological properties of Si3N4/GCr15 sliding pairs with few layer graphene as oil additives. Ceram. Int. 2017, 43, 14218–14224. [Google Scholar] [CrossRef]
  22. Mao, J.Y.; Zhao, J.; Wang, W.; He, Y.Y.; Luo, J.B. Influence of the micromorphology of reduced graphene oxide sheets on lubrication properties as a lubrication additive. Tribol. Int. 2018, 117, 614–621. [Google Scholar] [CrossRef]
  23. Zheng, D.; Cai, Z.B.; Shen, M.X.; Li, Z.Y.; Zhu, M.H. Investigation of the tribology behaviour of the graphene nanosheets as oil additives on textured alloy cast iron Surface. Appl. Surf. Sci. 2016, 387, 66–75. [Google Scholar] [CrossRef]
  24. Yang, J.; Xia, Y.F.; Song, H.J.; Chen, B.B.; Zhang, Z.Z. Synthesis of the liquid-like graphene with excellent tribological properties. Tribol. Int. 2017, 105, 118–124. [Google Scholar] [CrossRef]
  25. Lin, J.S.; Wang, L.W.; Chen, G.H. Modification of graphene platelets and their tribological properties as a lubricant additive. Tribol. Lett. 2011, 41, 209–215. [Google Scholar] [CrossRef]
  26. Wu, X.H.; Gong, K.L.; Zhao, G.Q.; Lou, W.J.; Wang, X.B.; Liu, W.M. Mechanical synthesis of chemically bonded phosphorus-graphene hybrid as high-temperature lubricating oil additive. RSC Adv. 2018, 8, 4595–4603. [Google Scholar] [CrossRef]
  27. Zhang, L.; He, Y.; Feng, S.W.; Zhang, L.; Zhang, L.; Jiao, Z.L.; Zhan, Y.Q.; Wang, Y.J. Preparation and tribological properties of novel boehmite/graphene oxide nano-hybrid. Ceram. Int. 2016, 42, 6178–6186. [Google Scholar] [CrossRef]
  28. Moghaddam, M.B.; Goharshadi, E.K.; Moosavi, F. Structural and transport properties and solubility parameter of graphene/glycerol nanofluids: A molecular dynamics simulation study. J. Mol. Liq. 2016, 222, 82–87. [Google Scholar] [CrossRef]
  29. Goyenola, C.; Gueorguiev, G.K.; Stafstrom, S.; Hultman, L. Fullerene-like CSx: A first-principles study of synthetic growth. Chem. Phys. Lett. 2011, 506, 86–91. [Google Scholar] [CrossRef] [Green Version]
  30. Chen, S.H.; Sun, S.Q.; Li, C.L.; Pittman, C.U.; Lacy, T.E.; Hu, S.Q.; Gwaltney, S.R. Behavior of protruding lateral plane graphene sheets in liquid dodecane: Molecular dynamics simulations. J. Nanopart. Res. 2016, 18, 317–329. [Google Scholar] [CrossRef]
  31. Loya, A.; Stair, J.L.; Jafri, A.R.; Yang, K.; Ren, G.G. A molecular dynamic investigation of viscosity and diffusion coefficient of nanoclusters in hydrocarbon fluids. Comp. Mater. Sci. 2015, 99, 242–246. [Google Scholar] [CrossRef]
  32. Tang, H.; Zhao, Y.; Yang, X.N.; Liu, D.M.; Shan, S.J.; Cui, F.Y. Understanding the roles of solution chemistries and functionalization on the aggregation of graphene-based nanomaterials using molecular dynamic simulations. J. Phys. Chem. C 2017, 121, 13888–13897. [Google Scholar] [CrossRef]
  33. Yang, J.; Yang, X.N.; Li, Y.P. Molecular simulation perspective of liquid-phase exfoliation, dispersion, and stabilization for graphene. Curr. Opin. Colloid Interface Sci. 2015, 20, 339–345. [Google Scholar] [CrossRef]
  34. Sangiovanni, D.G.; Gueorguiev, G.K.; Kakanakova-Georgieva, A. Ab initio molecular dynamics of atomic-scale surface reactions: Insights into metal organic chemical vapor deposition of AlN on graphene. PCCP 2015, 20, 17751–17761. [Google Scholar] [CrossRef] [PubMed]
  35. Saathoff, J.D.; Clancy, P. Effect of edge-functionalization on the ease of graphene nanoribbon aggregation in solvent. Carbon 2017, 115, 154–161. [Google Scholar] [CrossRef]
  36. Konatham, D.; Striolo, A. Molecular design of stable graphene nanosheets dispersions. Nano Lett. 2008, 8, 4630–4641. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, S.H.; Sun, S.Q.; Li, C.L.; Pittman, C.U.; Lacy, T.E.; Hu, S.Q.; Gwaltney, S.R. Molecular dynamics simulations of the graphene sheet aggregation in dodecane. Comput. Mater. Sci. 2017, 19, 195–206. [Google Scholar] [CrossRef]
  38. Cha, J.; Kyoung, W. Molecular dynamics simulation of the effects of affinity of functional groups and particle-size on the behavior of a graphene sheet in nanofluid. Comput. Mater. Sci. 2017, 139, 202–208. [Google Scholar] [CrossRef]
  39. Liu, P.Z.; Lu, J.; Yu, H.L.; Ren, N.; Lockwood, F.E.; Wang, Q.J. Lubricant shear thinning behavior correlated with variation of radius of gyration via molecular dynamics simulations. J. Chem. Phys. 2017, 147, 084904. [Google Scholar] [CrossRef] [PubMed]
  40. Xu, Y.X.; Shi, G.Q. Assembly of chemically modified graphene: Methods and applications. J. Mater. Chem. 2011, 21, 3311–3323. [Google Scholar] [CrossRef]
  41. Ai, K.L.; Liu, Y.L.; Lu, L.H.; Cheng, X.L.; Huo, L.H. A novel strategy for making soluble reduced graphene oxide sheets cheaply by adopting an endogenous reducing agent. J. Mater. Chem. 2011, 21, 3365–3370. [Google Scholar] [CrossRef]
  42. Choudhary, S.; Mungse, H.P.; Khatri, O.P. Dispersion of alkylated graphene in organic solvents and its potential for lubrication applications. J. Mater. Chem. 2012, 22, 21032–21039. [Google Scholar] [CrossRef]
  43. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef]
  44. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef]
  45. Hoover, W. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef]
  46. Allen, M.P.; Tildesley, D.J. Computer Simulation of Liquids; Oxford University Press: Oxford, UK, 1989. [Google Scholar]
  47. Khakan, H.; Yeganegi, S. Molecular dynamics simulations of amide functionalized imidazolium bis (trifluoromethanesulfonyl) imide dicationic ionic liquids. J. Phys. Chem. B 2017, 121, 7455–7463. [Google Scholar] [CrossRef] [PubMed]
  48. Michalis, V.K.; Moultos, O.A.; Tsimpanogianis, I.N.; Economou, I.G. Molecular dynamics simulations of the diffusion coefficients of light n-alkanes in water over a wide range of temperature and pressure. Fluid Phase Equilib. 2016, 407, 236–242. [Google Scholar] [CrossRef]
Figure 1. Snapshot of EMD simulation of (a) nanofluid cell with (b) graphene and (c) 1-decene trimer molecules.
Figure 1. Snapshot of EMD simulation of (a) nanofluid cell with (b) graphene and (c) 1-decene trimer molecules.
Crystals 08 00361 g001
Figure 2. Graphene plates with various functional groups: (a) pure graphene (b) GO (c) MG and (d) AG.
Figure 2. Graphene plates with various functional groups: (a) pure graphene (b) GO (c) MG and (d) AG.
Crystals 08 00361 g002
Figure 3. Equilibrium curve for energy of the graphene-PAO nanofluids.
Figure 3. Equilibrium curve for energy of the graphene-PAO nanofluids.
Crystals 08 00361 g003
Figure 4. (a) MSD curves and (b) diffusion coefficients of graphene with various sizes at 300 K and 0.1 MPa.
Figure 4. (a) MSD curves and (b) diffusion coefficients of graphene with various sizes at 300 K and 0.1 MPa.
Crystals 08 00361 g004
Figure 5. Diffusion coefficients of graphene with various functional groups calculated under NVT ensemble at 300 K and 0.1 MPa.
Figure 5. Diffusion coefficients of graphene with various functional groups calculated under NVT ensemble at 300 K and 0.1 MPa.
Crystals 08 00361 g005
Figure 6. Photographs of blended lubricating oil containing (a) graphene obtained 3 h, 10 h, 100 h after sonication and (b) graphene oxide 0 week, 1 week, 2 weeks after sonication at room temperature.
Figure 6. Photographs of blended lubricating oil containing (a) graphene obtained 3 h, 10 h, 100 h after sonication and (b) graphene oxide 0 week, 1 week, 2 weeks after sonication at room temperature.
Crystals 08 00361 g006
Figure 7. Radial distribution function for graphene-PAO molecules at 300 K and 0.1 MPa from EMD simulations.
Figure 7. Radial distribution function for graphene-PAO molecules at 300 K and 0.1 MPa from EMD simulations.
Crystals 08 00361 g007
Figure 8. Diffusion coefficients of graphene containing various number of functional groups calculated at NVT ensemble at 300 K and 1 atm.
Figure 8. Diffusion coefficients of graphene containing various number of functional groups calculated at NVT ensemble at 300 K and 1 atm.
Crystals 08 00361 g008
Figure 9. Diffusion coefficients of graphene containing two functional groups with respect to increasing temperature at 0.1 MPa.
Figure 9. Diffusion coefficients of graphene containing two functional groups with respect to increasing temperature at 0.1 MPa.
Crystals 08 00361 g009
Figure 10. Diffusion coefficients of graphene containing two functional groups with respect to increasing pressures at 300 K.
Figure 10. Diffusion coefficients of graphene containing two functional groups with respect to increasing pressures at 300 K.
Crystals 08 00361 g010
Table 1. The reductions in diffusion coefficients for various graphene containing different numbers of functional groups at 300 K and 0.1 MPa (%).
Table 1. The reductions in diffusion coefficients for various graphene containing different numbers of functional groups at 300 K and 0.1 MPa (%).
Number2468
GO48.1860.9168.1887.27
MG70.9179.0990.096.36
AG91.4592.6494.0996.73

Share and Cite

MDPI and ACS Style

Wu, L.; Song, B.; Keer, L.M.; Gu, L. Molecular Dynamics Investigation of Graphene Nanoplate Diffusion Behavior in Poly-α-Olefin Lubricating Oil. Crystals 2018, 8, 361. https://doi.org/10.3390/cryst8090361

AMA Style

Wu L, Song B, Keer LM, Gu L. Molecular Dynamics Investigation of Graphene Nanoplate Diffusion Behavior in Poly-α-Olefin Lubricating Oil. Crystals. 2018; 8(9):361. https://doi.org/10.3390/cryst8090361

Chicago/Turabian Style

Wu, Lupeng, Baoyu Song, Leon M. Keer, and Le Gu. 2018. "Molecular Dynamics Investigation of Graphene Nanoplate Diffusion Behavior in Poly-α-Olefin Lubricating Oil" Crystals 8, no. 9: 361. https://doi.org/10.3390/cryst8090361

APA Style

Wu, L., Song, B., Keer, L. M., & Gu, L. (2018). Molecular Dynamics Investigation of Graphene Nanoplate Diffusion Behavior in Poly-α-Olefin Lubricating Oil. Crystals, 8(9), 361. https://doi.org/10.3390/cryst8090361

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