Thermodiffusion in model nanofluids by molecular dynamics simulations

In this work, a new algorithm is proposed to compute single particle (infinite dilution) thermodiffusion using nonequilibrium molecular dynamics simulations through the estimation of the thermophoretic force that applies on a solute particle. This scheme is shown to provide consistent results for model nanofluids in the liquid state (spherical nonmetallic nanoparticles+Lennard-Jones fluid) where it appears that thermodiffusion amplitude, as well as thermal conductivity, decreases with nanoparticle concentration. Then, by changing the nature of the nanoparticle (size, mass, and internal stiffness) and that of the solvent (quality and viscosity), various trends are exhibited. In all cases, the single particle thermodiffusion is positive, i.e., the nanoparticle tends to migrate toward the cold area. The single particle thermal diffusion coefficient is shown to be independent of the size of the nanoparticle (diameter of 0.8-4 nm), whereas it increases with the quality of the solvent and is inversely proportional to the viscosity of the fluid. In addition, this coefficient is shown to be independent of the mass of the nanoparticle and to increase with the stiffness of the nanoparticle internal bonds. Besides, for these configurations, the mass diffusion coefficient behavior appears to be consistent with a Stokes-Einstein-like law.