Thermal conductivity calculation of magnetite using molecular dynamics simulation

Document Type : Research Paper


1 Faculty of Mechanical Engineering, Sahand University of Technology, Tabriz, Iran

2 Department of Mechanical Engineering, Azarbaijan Shahid Madani University, Tabriz, Iran


In the current research, thermal conductivity of magnetite (Fe3O4) has been calculated using molecular dynamic simulation. The rNEMD Molecular Dynamics Method provided in the LMMPS package is used for the simulation of the thermal conductivity. The effects of magnetite layer size and temperature on the thermal conductivity have been investigated. The numerical results have been validated by experimental data. Results show that the thermal conductivity of magnetite can be predicted appropriately using Buckingham potential function. Moreover, Thermal conductivity of magnetite is shown to be a decreasing function of temperature. The obtained results provide a benchmark for magnetite ferrofluid heat transfer simulations, which has been extensively increased in recent years.


[1] Masaaki Motozawa, Jia Chang, Tatsou Sawada, Yasuo Kawaguchi, “Effect of Magnetite Field on Heat Transfer in Rectangular Duct Flow of a Magnetite Fluid, ” 12th International Conference on Magnetic Fluids, Physics Procedia 9 (2010) 190-193.
[2] Mohammad Hassan Omidi, Mahboobeh Alibeygi, Farideh Piri, Mohammad Masoudifard, “Polystyrene/magnetite nanocomposite synthesis and characterization: investigation of magnetite and electrical properties for using as microelectromechanical systems (MEMS), ” Materials Science-Poland, 35(1), 2017, pp, 105-110
[3] J. Molgaard and W. W. Smeltzer, “Thermal conductivity of magnetite and hematite,” J. Appl. Phys., vol. 42, no. 9, pp. 3644–3647, 1971.
[4] SYDNEY P. Clark, “Handbook of physical constants,” The Geological Society of America Memoir, PP. 461-482, 1996.
[5] N.-W. Park et al., “Reduced temperature-dependent thermal conductivity of magnetite thin films by controlling film thickness.,” Nanoscale Res. Lett., vol. 9, no. 1, p. 96(8pp), 2014.
[6] H. Aminfar, N. Razmara, and M. Mohammadpourfard, “Molecular dynamics study of ferrofluid flow inside nanochannels under magnetic fields,” J. Comput. Theor. Nanosci., vol. 12, no. 9, pp. 2339–2347, Sep. 2015.
[7] A. R. bin Saleman, H. K. Chilukoti, G. Kikugawa, M. Shibahara, and T. Ohara, “A molecular dynamics study on the thermal transport properties and the structure of the solid-liquid interfaces between face-centered cubic (FCC) crystal planes of gold in contact with linear alkane liquids,” Int. J. Heat Mass Transf., vol. 105, pp. 168–179, Feb. 2017.
[8] A. Vegiri and E. I. Kamitsos, “Molecular dynamics study of structural reorganization by electro-thermal poling in sodium diborate glass,” J. Non. Cryst. Solids, vol. 472, pp. 14–24, Sep. 2017.
[9] L. Momenzadeh, B. Moghtaderi, O. Buzzi, X. Liu, S. W. Sloan, and G. E. Murch, “The thermal conductivity decomposition of calcite calculated by molecular dynamics simulation,” Comput. Mater. Sci., vol. 141, pp. 170–179, 2018.
[10] M. E. Tuckerman, “Statistical Mechanics: Theory and Molecular Simulation,” p. 720, 2010.
[11] F. Müller-Plathe, “A simple nonequilibrium molecular dynamics method for calculating the thermal conductivity,” J. Chem. Phys., vol. 106, no. 14, pp. 6082–6085, 1997.
[12] J. Severin and P. Jund, “Thermal conductivity calculation in anisotropic crystals by molecular dynamics: Application to α−Fe2O3,” J. Chem. Phys., vol. 146, no. 5, 2017.
[13] D. Chicot et al., “Mechanical properties of magnetite (Fe3O4), hematite (α-Fe2O3) and goethite (α-FeO·OH) by instrumented indentation and molecular dynamics analysis,” Mater. Chem. Phys., vol. 129, no. 3, pp. 862–870, 2011.
[14] S. Hess, J. B. Hayter, and R. Pynn, “A comparison of molecular dynamics and analytic calculation of correlations in an aligned ferrofluid,” Mol. Phys., vol. 53, no. 6, pp. 1527–1533, 1984.
[15] J. Vaari, “Molecular dynamics simulations of vacancy diffusion in chromium(III) oxide, hematite, magnetite and chromite,” Solid State Ionics, vol. 270, pp. 10–17, 2015.
[16] R. M. Cornell and U. Schwertmann, Iron Oxides. ISBN: 3-527-30274-3 . 2003.
[17] L. Minervini, M. O. Zacate, and R. W. Grimes, “Defect cluster formation in M2O3-doped CeO2,” Solid State Ionics, vol. 116, no. 3, pp. 339–349, Jan. 1999.
[18] G. V Lewis and C. R. A. Catlow, “Potential models for ionic oxides,” J. Phys. C Solid State Phys., vol. 18, no. 6, pp. 1149–1161, Feb. 1985.
[19] S. Y. Hideo Kaburaki, Ju Li, “Thermal Conductivity of Solid Argon by Classical Molecular Dynamics,” Mat. Res. Soc. Symp. Proc. Mater. Res. Soc., vol. 538, no. 0, pp. 503–508, 1999.
[20] D. K. Christen and G. L. Pollack, “Thermal conductivity of solid argon,” Phys. Rev. B, vol. 12, no. 8, pp. 3380–3391, 1975.
[21] M. Alaghemandi, E. Algaer, M. C. Böhm, and F. Müller-Plathe, “The thermal conductivity and thermal rectification of carbon nanotubes studied using reverse non-equilibrium molecular dynamics simulations,” Nanotechnology, vol. 20, no. 11, 2009.
[22] A. Soleimani, H. Araghi, Z. Zabihi, and A. Alibakhshi, “A comparative study of molecular dynamics simulation methods for evaluation of the thermal conductivity and phonon transport in Si nanowires,” Comput. Mater. Sci., vol. 142, pp. 346–354, 2018.
[23] J. R. Lukes and H. Zhong, “Thermal Conductivity of Individual Single-Wall Carbon Nanotubes,” J. Heat Transfer, vol. 129, no. 6, p. 705, 2007.
[24] M. Khalkhali and F. Khoeini, “Impact of torsion and disorder on the thermal conductivity of Si nanowires: A nonequilibrium molecular dynamics study,” J. Phys. Chem. Solids, vol. 112, pp. 216–221, 2018.