Comparative study between continuum and atomistic approaches of liquid flow through a finite length cylindrical nanopore

Steady state pressure driven flow of liquid argon through a finite length cylindrical nanopore was investigated numerically by classical Navier-Stokes (NS) hydrodynamic models and nonequilibrium molecular dynamics (MD) simulations. In both approaches, the nanopore was nominally 2.2 nm in diameter and 6 nm long. For the MD simulations, the intermolecular properties of the walls were specified independently from the liquid. Comparisons between the approaches were made in terms of the gross feature of total flow rate through the nanopore, as well as the more refined considerations of the spatial distributions of pressure, density, and velocity. The results showed that for the NS equations to predict the same trends in total flow rate with increasing pressure difference as the MD simulation, submodels for variations in density and viscosity with pressure are needed to be included. The classical NS boundary conditions quantitatively agreed with the flow rate predictions from MD simulations only under the condition of having a neutral-like solid-liquid interaction. Under these conditions, the NS and MD models also agreed well in streamwise distributions of pressure, density, and velocity, but not in the radial direction.     

Statistical-mechanical theory of rheology: Lennard-Jones fluids

The generalized Boltzmann equation for simple dense fluids gives rise to the stress tensor evolution equation as a constitutive equation of generalized hydrodynamics for fluids far removed from equilibrium. It is possible to derive a formula for the non-Newtonian shear viscosity of the simple fluid from the stress tensor evolution equation in a suitable flow configuration. The non-Newtonian viscosity formula derived is applied to calculate the non-Newtonian viscosity as a function of the shear rate by means of statistical mechanics in the case of the Lennard-Jones fluid. For that purpose we have used the density-fluctuation theory for the Newtonian viscosity, the modified free volume theory for the self-diffusion coefficient, and the generic van der Waals equation of state to compute the mean free volume appearing in the modified free volume theory. Monte Carlo simulations are used to calculate the pair-correlation function appearing in the generic van der Waals equation of state and shear viscosity formula. To validate the Newtonian viscosity formula obtained we first have examined the density and temperature dependences of the shear viscosity in both subcritical and supercritical regions and compared them with molecular-dynamic simulation results. With the Newtonian shear viscosity and thermodynamic quantities so computed we then have calculated the shear rate dependence of the non-Newtonian shear viscosity and compared it with molecular-dynamics simulation results. The non-Newtonian viscosity formula is a universal function of the product of reduced shear rate (gamma*) times reduced relaxation time (taue*) that is independent of the material parameters, suggesting a possibility of the existence of rheological corresponding states of reduced density, temperature, and shear rate. When the simulation data are reduced appropriately and plotted against taue*gamma* they are found clustered around the reduced (universal) non-Newtonian viscosity formula. Thus we now have a molecular theory of non-Newtonian shear viscosity for the Lennard-Jones fluid, which can be implemented with a Monte Carlo simulation method for the pair-correlation function.     

Equilibrium calculations of viscosity and thermal conductivity across a solid-liquid interface using boundary fluctuations

We calculate viscosity and thermal conductivity in systems of Lennard-Jones particles consisting of coexisting solid and liquid with different interface wetting properties using the recently developed equilibrium boundary fluctuation theory. We compare the slip length and equivalent liquid length obtained from these calculations with those obtained from nonequilibrium molecular dynamics. The equilibrium and nonequilibrium calculations of the slip length and the sum of the thermal equivalent lengths are in good agreement. We conclude that for both interfacial properties, the nonequilibrium simulations were probing the linear response. The significant dependence of the intrinsic equivalence length on the interfacial temperature difference used to generate the thermal gradient is explained as a consequence of the different thermodynamic states of the two interfaces.     

Shear viscosity of liquid copper at experimentally accessible shear rates: application of the transient-time correlation function formalism

We report on nonequilibrium molecular dynamics (NEMD) simulations results on the shear viscosity of liquid copper, modeled by a many-body embedded-atoms model potential. Because conventional NEMD methods are restricted to very high shear rates (at least of the order of 10(10) s(-1), that is several orders of magnitude larger than those accessible by experiment), previous work only provided access to the response of the fluid in the shear-thinning regime. Using the transient-time correlation function formalism, we show how NEMD simulations can be extended to study the rheological properties of liquid copper subjected to low, experimentally accessible, shear rates. Our results provide a full picture of the rheology of the system, in the Newtonian regime as well as in the shear-thinning regime.     

含氰基离子液体的合成、表征及流变性质研究

摘要合成、 表征了一系列新的含氰基咪唑类离子液体. 测定了该离子液体的密度、 熔点及溶解性等物理性质, 研究了其在稳态、 瞬态和动态条件下的流变行为. 结果表明, 当剪切速率在0.1～50 s^-1范围内时, 其粘度不随剪切速率的变化而变化, 但随温度升高而降低, 粘流活化能随取代基长度变化呈现规律性变化. 对于1-丁基-3-氰乙基咪唑六氟磷酸盐离子液体, 维持剪切速率不变时, 其剪切应力和粘度均不随时间变化, 且随着温度的升高而降低; 在动态条件下, 在线性粘弹区, 复合粘度和损耗模量G″ 随温度升高而降低. 关键词     

Rheology of a thermotropic polyester in the nematic and isotropic states

The rheological properties of a thermotropic polyester were determined in the nematic and isotropic states. In the isotropic state, the viscosity is almost constant and the polymer is only slightly elastic. The nematic phase has a lower viscosity than the isotropic, except at low frequencies or shear rates, where the viscosity increases as though the polymer had a yield stress. There is a marked dependence of the rheology on shear history. The effects of shearing can be erased by returning the material first to the isotropic state and then back to the nematic state. The results are discussed with reference to analogous observations in small-molecule liquid crystals and in thermotropic aromatic co-polyesters.     

Nonlinear shear of entangled polymers from nonequilibrium molecular dynamics

This study aims to use molecular dynamics (MD) simulations of Kremer–Grest (KG) chains to inform future developments of models of entangled polymer dynamics. We perform nonequilibrium MD simulations, under shear flow, for well‐entangled KG chains. We study chains of 512 and 1000 KG beads, corresponding to 8 and 15 entanglements, respectively. We compute the linear rheological properties from equilibrium simulations of the stress autocorrelation and obtain from these data the tube model parameters. Under nonlinear shear flow, we compute the shear viscosity, the first and second normal stress differences, and chain contour length. For chains of 512 monomers, we obtain agreement with the results of Cao and Likhtman (ACS Macro. Lett. 2015, 4, 1376). We also compare our nonlinear results with the Graham, Likhtman and Milner‐McLeish (GLaMM) model. We identify some systematic disagreement that becomes larger for the longer chains. We made a comparison of the transient shear stress maximum from our simulations, two nonlinear models and experiments on a wide range of melts and solutions, including polystyrene (PS), polybutadiene, and styrene–butadiene rubber. This comparison establishes that the PS melt data show markedly different behavior to all other melts and solutions and KG simulations reproduce the PS data more closely than either the GLaMM or Xie and Schweizer models. We discuss the performance of these models against the data and simulations. Finally, by imposing a rapid reversing flow, we produce a method to extract the recoverable strain from MD simulations, valid for sufficiently entangled monodisperse polymers. We explore how the resulting data can probe the melt state just before the reversing flow. © 2019 The Authors. Journal of Polymer Science Part B: Polymer Physics published by Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 1692–1704     

Nonequilibrium molecular dynamics simulation to describe the rotation of rigid, low aspect ratio carbon nanotubes in simple shear flow

In this paper, the rotation of short carbon nanotubes in simple shear liquid argon flow was investigated by nonequilibrium molecular dynamics (MD) simulation. In their simulations, nanotubes were described as rigid cylinders of carbon atoms. Lennard-Jones potential was employed to represent both argon-argon and argon-carbon interactions. Results show that time period of a nanotube as calculated from MD simulations is longer than what would be calculated from Jeffery's equation based on the aspect ratio of the cylinder. The difference is much higher at low shear rates and for small aspect ratios. Results also reveal that adding caps to an open-ended nanotube speeds up its rotation.     

Molecular dynamics simulation of liquid sulfur dioxide

A previously proposed model for molecular dynamics (MD) simulation of liquid sulfur dioxide, SO(2), has been reviewed. Thermodynamic, structural, and dynamical properties were calculated for a large range of thermodynamic states. Predicted (P,V,T) of simulated system agrees with an elaborated equation of state recently proposed for liquid SO(2). Calculated heat capacity, expansion coefficient, and isothermal compressibility are also in good agreement with experimental data. Calculated equilibrium structure agrees with X-ray and neutron scattering measurements on liquid SO(2). The model also predicts the same (SO(2))(2) dimer structure as previously determined by ab initio calculations. Detailed analysis of equilibrium structure of liquid SO(2) is provided, indicating that, despite the rather large dipole moment of the SO(2) molecule, the structure is mainly determined by the Lennard-Jones interactions. Both single-particle and collective dynamics are investigated. Temperature dependency of dynamical properties is given. The MD results are compared with previous findings obtained from the analysis of inelastic neutron scattering spectra of liquid SO(2), including wave-vector dependent structural relaxation, tau(k), and viscosity, eta(k).     

The effects of some factors on the rheological properties of the lyotropic liquid crystals formed in Brij97/NaDC/IPM/water system: compositions,temperature and polyphenols

In this paper, interaction between Brij97 and NaDC in mixed micelles was studied by surface tension method. Phase behavior for Brij97/NaDC/IPM/H₂O system was investigated at the desire molar ratio of α_NaDC?=?0.42 and the human body temperature of 37?°C. Hexagonal and cubic liquid crystal were found in this system analyzed by small angle x-ray scattering (SAXS). some factors were changed to affect the rheological properties of liquid crystal. Rheological investigations showed: The higher content of water and lower experimental temperature favored the greater values of shear viscosity and viscoelastic modulus; The sample was closer to the plastic fluid and was more likely to recover from disturbed state to equilibrium state at lower ratio of NaDC; The curcumin and TP can influence the shear viscosity, viscoelastic properties and phase transition temperature, even change the structure of liquid crystal. This supplies a way to adjust rheological properties, structure and even phase transition by altering those factors.     

Assessment of phenomenological models for viscosity of liquids based on nonequilibrium atomistic simulations of copper

The shear viscosity of liquid copper is studied using nonequilibrium molecular-dynamics simulations under planar shear flow conditions. We examined variation of viscosity as function of shear rate at a range of pressures (ca. 0 - 40 GPa). We analyzed these results using eight different phenomenological models and find that the observed non-Newtonian behavior is best described by the Powell-Eyring (PE) model: eta(gamma) = (eta(0)-eta(infinity))sinh(-1)(taugamma)(taugamma) + eta(infinity), where gamma is the shear rate. Here eta(0) (the zero-shear-rate viscosity) extracted from the PE fit is in excellent agreement with available experimental data. The relaxation time tau from the PE fit describes the shear response to an applied stress. This provides the framework for interpreting the shear flow phenomena in complex systems, such as liquid metal and amorphous metal alloys.     

Structural relaxation and rheological response of a driven amorphous system

The interplay between the structural relaxation and the rheological response of a simple amorphous system {a 80:20 binary Lennard-Jones mixture [W. Kob and H. C. Andersen, Phys. Rev. Lett. 73, 1376 (1994)]} is studied via molecular dynamics simulations. In the quiescent state, the model is well known for its sluggish dynamics and a two step relaxation of correlation functions at low temperatures. An ideal glass transition temperature of Tc=0.435 has been identified in the previous studies via the analysis of the system's dynamics in the framework of the mode coupling theory of the glass transition [W. Kob and H. C. Andersen, Phys. Rev. E 51, 4626 (1995)]. In the present work, we focus on the question whether a signature of this ideal glass transition can also be found in the case where the system's dynamics is driven by a shear motion. Indeed, the following distinction in the structural relaxation is found: In the supercooled state, the structural relaxation is dominated by the shear at relatively high shear rates gamma, whereas at sufficiently low gamma the (shear-independent) equilibrium relaxation is recovered. In contrast to this, the structural relaxation of a glass is always driven by shear. This distinct behavior of the correlation functions is also reflected in the rheological response. In the supercooled state, the shear viscosity eta decreases with increasing shear rate (shear thinning) at high shear rates, but then converges toward a constant as the gamma is decreased below a (temperature-dependent) threshold value. Below Tc, on the other hand, the shear viscosity grows as eta proportional, etax 1/gamma, suggesting a divergence at gamma=0. Thus, within the accessible observation time window, a transition toward a nonergodic state seems to occur in the driven glass as the driving force approaches zero. As to the flow curves (stress versus shear rate), a plateau forms at low shear rates in the glassy phase. A consequence of this stress plateau for Poiseuille-type flows is demonstrated.     

Viscosity and structure configuration properties of equilibrium and supercooled liquid cobalt

The shear viscosity of liquid cobalt at the pressure p = 1.5 bar and at the temperatures corresponding to equilibrium liquid and supercooled liquid states is measured experimentally and evaluated by means of molecular dynamics simulations. Further, the shear viscosity is also calculated within the microscopic theoretical model. Comparison of our experimental, simulation and theoretical results with other available data allows one to examine the issue about the correct temperature dependence of the shear viscosity of liquid cobalt. A strong correlation between the viscosity and the configuration entropy of liquid cobalt over the considered temperature range is found, which can be taken into account by the Rosenfeld’s model.     

Pressure dependence of viscosity

We reanalyze the pressure dependence of viscosity of liquids of constant composition under isothermal conditions. Based exclusively on very general considerations concerning the relationship between viscosity and "free volume," we show that, at moderate values of pressure, viscosity increases, as a rule, with increasing pressure, provided the liquid is in stable or metastable (undercooled) equilibrium states. However, even if the behavior of the viscosity is governed by free volume effects, deviations from a positive pressure dependence are possible, when the liquid's thermal expansion coefficient is negative. We derive an equation that allows one to quantitatively determine the pressure dependence of viscosity, which requires, in the simplest case, only the knowledge of the temperature dependence of viscosity at constant pressure, the thermal expansion coefficient, and the isothermal compressibility of the liquid. As an example, the negative pressure dependence of water in the range of temperatures 0-4 degrees C and of several silicate liquids, such as albite, jadeite, dacite, basalts, etc., could be explained in such a way. Other glass-forming liquids initially (for moderate pressures) show a positive pressure dependence of viscosity that changes to a negative one when subjected to high (approximately GPa) isostatic pressure. A detailed analysis of water and already mentioned silicate melts at GPa pressures shows that, in addition to free volume effects, other pressure induced structural transformations may have to be accounted for in a variety of cases. By this reason, the theoretical analysis is extended (i) in order to describe the pressure dependence of viscosity for systems that are in frozen-in thermodynamic nonequilibrium states (glasses, i.e., undercooled liquids below the glass transition temperature Tg) and (ii) to systems which undergo, in addition to variations of the free volume, pressure induced changes of other structural parameters. In such cases a decrease of viscosity with increasing pressure may occur, in principle, even if the thermal expansion coefficient is positive. In this way, the present analysis grants a general tool to estimate the pressure dependence of viscosity and supposedly settles the controversy in the current literature.     

Investigating the temperature dependence of the viscosity of a non-Newtonian fluid within lithographically defined microchannels

We present a study of the rheological phenomenology of a non-Newtonian glass former within hybrid microchannels above the vitrification region. We determined the temperature behavior of the viscosity, which is well fitted by a Vogel-Fulcher-Tamman law for shear rates between 4 x 10(-2) and 9 x 10(-1) s(-1). The microflow viscosity was compared with previously reported conductivity data of the investigated molecular system. Our findings provide an insight into the coupling between the structural dynamics in the bulk and that within the microchannels, suggesting lithographically defined microfluidic systems as promising tools for the investigation of the rheological properties of complex liquids.     

Formation of a high coverage (3 x 3) NO phase on Pd(111) at elevated pressures: interplay between kinetic and thermodynamic accessibility

Using in situ polarization modulation infrared reflection absorption spectroscopy and density functional theory calculations, a new high-coverage monomeric NO adsorption state on Pd(111) was observed and proposed to have a (3 x 3)-7NO structure. Formation of this high coverage NO phase was found to take place only at elevated pressure and temperature conditions showing that some of the accessible thermodynamic equilibrium states at elevated temperatures and pressures are thermodynamically unfavorable or kinetically hindered at lower temperatures and pressures. Our results emphasize the danger of extrapolating results from traditional surface science experiments performed under ultrahigh vacuum to elevated temperature and pressure conditions encountered in heterogeneous catalysis.     

Thermodynamics,viscous flow and relaxation dynamics of bulk glass-forming Pd alloys

Experimental results concerning the thermodynamic excess properties and the equilibrium shear viscosity of the undercooled liquid state are presented, which were consistently obtained on Pd-rich bulk metallic glasses of different stability by applying the concept of the limiting fictive temperature. These results are evaluated with respect to the stability criteria that are controlling vitrification. In addition, modulated-temperature calorimetry has been performed to obtain the average relaxation time of the cooperative liquid modes near the glass transition directly in dependence of temperature. The results are compared to data obtained from static non-equilibrium measurements and discussed with respect to the relation between fragility and thermal stability of Pd-rich bulk glass-forming alloys.     

Solubility of carbon dioxide,methane, and ethane in 1-butanol and saturated liquid densities and viscosities

A designed pressure–volume–temperature (PVT) apparatus has been used to measure the (vapor + liquid) equilibrium properties of three binary mixtures (methane +, ethane +, and carbon dioxide + 1-butanol) at two temperatures (303 and 323) K and at the pressures up to 6 MPa. The solubility of the compressed gases in 1-butanol and the saturated liquid densities and viscosities were measured. In addition, the density and viscosity of pure 1-butanol were measured at two temperatures (303 and 323) K and at the pressures up to 10 MPa. The experimental results show that the solubility of the gases in 1-butanol increases with pressure and decreases with temperature. The dissolution of gases in 1-butanol causes a decline in the viscosity of liquid phase. The saturated liquid density follows a decreasing trend with the solubility of methane and ethane. However, the dissolution of carbon dioxide in 1-butanol leads to an increase in the density of liquid phase. The experimental data are well correlated with Soave–Redlich–Kwong (SRK) and Peng–Robinson (PR) equations of state (EOSs). SRK EOS was slightly superior for correlating the saturated liquid densities.