The present work uses the concept of a scaled particle along with the perturbation and variation approach, to develop an equation of state (EOS) for a mixture of hard sphere (HS), Lennar-Jones (L J) fluids. A su...The present work uses the concept of a scaled particle along with the perturbation and variation approach, to develop an equation of state (EOS) for a mixture of hard sphere (HS), Lennar-Jones (L J) fluids. A suitable flexible functional form for the radial distribution function G(R) is assumed for the mixture, with R as a variable. The function G(R) has an arbitrary parameter m and a different equation of state can be obtained with a suitable choice of m. For m = 0.75 and m = 0.83 results are close to molecular dynamics (MD) result for pure HS and LJ fluid respectively.展开更多
Using nonequilibrium molecular dynamics simulations, we study the non-Newtonian rhe-ological behaviors of a monoatomic fluid governed by the Lennard-Jones potential. Both steady Couette and oscillatory shear flows are...Using nonequilibrium molecular dynamics simulations, we study the non-Newtonian rhe-ological behaviors of a monoatomic fluid governed by the Lennard-Jones potential. Both steady Couette and oscillatory shear flows are investigated. Shear thinning and normal stress effects are observed in the steady Couette flow simulations. The radial distribution function is calculated at different shear rates to exhibit the change of the microscopic struc- ture of molecules due to shear. We observe that for a larger shear rate the repulsion between molecules is more powerful while the attraction is weaker, and the above phenomena can also be confirmed by the analyses of the potential energy. By applying an oscillatory shear to the system, several findings are worth mentioning here: First, the phase difference between the shear stress and shear rate increases with the frequency. Second, the real part of complex viscosity first increases and then decreases while the imaginary part tends to increase mono- tonically, which results in the increase of the proportion of the imaginary part to the real part with the increasing frequency. Third, the ratio of the elastic modulus to the viscous modulus also increases with the frequency. These phenomena all indicate the appearance of viscoelasticity and the domination of elasticity over viscosity at high oscillation frequency for Lennard-Jones fluids.展开更多
Both a free volume approach for Helmholtz free energy and a theoretically-based fitted formula for radial distribution function (rdf) of hard sphere solid are employed to describe the Helmholtz free energy of Lennard-...Both a free volume approach for Helmholtz free energy and a theoretically-based fitted formula for radial distribution function (rdf) of hard sphere solid are employed to describe the Helmholtz free energy of Lennard-Jones solid in the framework of the first order thermodynamic perturbation theory, which also is employed for the uniform LennardJones fluid. The dividing of the Lennard-Jones potential follows from the WCA prescription, but the specification of the equivalent hard sphere diameter is determined by a simple iteration procedure devised originally for liquid state, but extended to solid state in the present study. Two hundred sheiks are used in the rdf to get an accurate perturbation term.The present approach is very accurate for the description of excess Helmholtz free energy of LJ solid, but shows some deviation from the simulation for excess Helmholtz free energy of uniform LJ fluid when the reduced temperature kT/ε is higher then 5. The present approach is satisfactory for description of solid-liquid phase transition of the Lennard-Jones model.展开更多
文摘The present work uses the concept of a scaled particle along with the perturbation and variation approach, to develop an equation of state (EOS) for a mixture of hard sphere (HS), Lennar-Jones (L J) fluids. A suitable flexible functional form for the radial distribution function G(R) is assumed for the mixture, with R as a variable. The function G(R) has an arbitrary parameter m and a different equation of state can be obtained with a suitable choice of m. For m = 0.75 and m = 0.83 results are close to molecular dynamics (MD) result for pure HS and LJ fluid respectively.
文摘Using nonequilibrium molecular dynamics simulations, we study the non-Newtonian rhe-ological behaviors of a monoatomic fluid governed by the Lennard-Jones potential. Both steady Couette and oscillatory shear flows are investigated. Shear thinning and normal stress effects are observed in the steady Couette flow simulations. The radial distribution function is calculated at different shear rates to exhibit the change of the microscopic struc- ture of molecules due to shear. We observe that for a larger shear rate the repulsion between molecules is more powerful while the attraction is weaker, and the above phenomena can also be confirmed by the analyses of the potential energy. By applying an oscillatory shear to the system, several findings are worth mentioning here: First, the phase difference between the shear stress and shear rate increases with the frequency. Second, the real part of complex viscosity first increases and then decreases while the imaginary part tends to increase mono- tonically, which results in the increase of the proportion of the imaginary part to the real part with the increasing frequency. Third, the ratio of the elastic modulus to the viscous modulus also increases with the frequency. These phenomena all indicate the appearance of viscoelasticity and the domination of elasticity over viscosity at high oscillation frequency for Lennard-Jones fluids.
文摘Both a free volume approach for Helmholtz free energy and a theoretically-based fitted formula for radial distribution function (rdf) of hard sphere solid are employed to describe the Helmholtz free energy of Lennard-Jones solid in the framework of the first order thermodynamic perturbation theory, which also is employed for the uniform LennardJones fluid. The dividing of the Lennard-Jones potential follows from the WCA prescription, but the specification of the equivalent hard sphere diameter is determined by a simple iteration procedure devised originally for liquid state, but extended to solid state in the present study. Two hundred sheiks are used in the rdf to get an accurate perturbation term.The present approach is very accurate for the description of excess Helmholtz free energy of LJ solid, but shows some deviation from the simulation for excess Helmholtz free energy of uniform LJ fluid when the reduced temperature kT/ε is higher then 5. The present approach is satisfactory for description of solid-liquid phase transition of the Lennard-Jones model.