Search

The zero-density shear viscosity of different types of short Lennard-Jones chains, up to the hexa-decamer, has been evaluated using a non-equilibrium molecular dynamics scheme. Simulations have been performed on chains of variable rigidities going from the fully flexible to the fully rigid chains. Very interestingly, it is found that there exists a universal relation (a power law) between the zero-density viscosity of the Lennard-Jones chains and their radius of gyration whatever the rigidity of the chain and for all tested temperatures (ranging from 2.5 to 6 in reduced units). Furthermore, for the studied range of temperature, it is shown that the zero-density viscosity of both fully flexible chains and fully rigid chains models can be obtained with an accuracy of a few percents knowing only the dimer viscosity and the length of the chain. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Sorer effect in multicomponent system (more than 2 components) are being studied starting from binary solutions. Possible predictive theories are extensions of the thermodynamic theories, molecular dynamics, but also by the mechanical statistics extensions. The latter are illustrated at the light of the experiments obtained in microgravity conditions. For binary melts, a very good concordance between measurements on samples processed in microgravity and the statistical theory is obtained under the form of thermodiffusion factors as function of mass misfit. The same approach is used for hydrocarbon ternary and quaternary solution processed in microgravity in the SCCO experiment. [ABSTRACT FROM AUTHOR]

In recent years, extensive research has used molecular dynamics simulations to investigate gas separation through nanoporous graphene (NPG) membranes. However, most studies have considered graphene membranes as rigid, overlooking the impact of their inherent flexibility. This study systematically quantifies the effect of graphene flexibility on gas permeation by comparing the diffusion of various gases through flexible and rigid single-layer NPG models. The results demonstrate that flexibility notably increases permeance, particularly for gases with larger molecular diameters/pore size ratios, by allowing gas molecules greater mobility within the pore. Interestingly, the effect of flexibility boils down to the expansion of the average pore size, and the detail of the membrane's vibrational dynamics is of little importance in quantifying permeance. Our work shows that accounting for flexibility in molecular models improves the alignment of simulation results with experimental data, emphasizing the importance of considering membrane flexibility in predictive models of NPG membrane performance.

In this paper, we study the permeation of polyatomic gas molecules through 2D graphene membranes. Using equilibrium molecular dynamics simulations, we investigate the permeation of pure gas compounds (CH4, CO2, O2, N2, and H2) through nanoporous graphene membranes with varying pore sizes and geometries. Our simulations consider the recrossing mechanism, often neglected in previous studies, which has a significant effect on permeation for intermediate pore size to molecular diameter ratios. We find that the permeation process can be decoupled into two steps: the crossing process of gas molecules through the pore plane and the escaping process from the pore region to a neighboring adsorption site, which prevents recrossing. To account for these mechanisms, we use a permeance model expressed as the product of the permeance for the crossing process and the probability of molecule escape. This phenomenological model is extended to account for small polyatomic gas molecules and to describe permeation regimes ranging from molecular sieving to effusion. The proposed model captures the temperature dependence and provides insights into the key parameters of the gas/membrane interaction controlling the permeance of the system. This work lays the foundation for predicting gas permeance and exploring membrane separation factors in 2D materials such as graphene., (© 2023 Author(s). Published under an exclusive license by AIP Publishing.)

Two-dimensional (2D) membranes based on perforated graphene have great potential in the field of separation of chemical species for a variety of applications, including gas treatment. In addition to recent experimental studies, several works simulate the mechanisms of gas permeation through this type of membrane using molecular dynamics, but few combine different techniques to ensure that their method of choice captures all relevant mechanisms. In particular, the re-crossing mechanism leading a gas molecule that has crossed the plane of the membrane to rapidly re-cross it in the opposite direction has never been documented. In this work, we study gas permeation through a simplified 2D membrane model. We combine equilibrium and non-equilibrium molecular dynamics simulations to quantify the impact of these re-crossing mechanisms on the values of the computed transport coefficients. Using non-equilibrium simulations as reference, we show that the equilibrium simulation techniques commonly used can lead to a significant overestimation of the transport properties of the membrane. We propose a simple method to probe the re-crossing dynamics during equilibrium simulations, making it possible to compute correct values of the transport coefficient without the need for non-equilibrium simulations. Furthermore, by analyzing the phenomenology observed in the simulations, we derive an analytical formula for the permeance that takes the form of an Arrhenius law with a non-trivial temperature dependent prefactor. In excellent agreement with our simulation results, this model provides a simple theoretical framework that captures the main mechanisms involved in gas permeation through 2D membranes, including the effect of re-crossing.

The predictive capabilities of some existing theoretical models to quantify thermodiffusion have been investigated in this work. To do so, the tests have been performed on two model fluids, the hard-sphere and the Lennard-Jones (including spheres and dimers) ones, exploring different mixtures and thermodynamic conditions thanks to extensive molecular simulations. It has been confirmed that the thermal diffusion factor should be expressed as the sum of one term related to the isotope effect and one term related to the "chemical" effects and that a kinetic term is required to quantify thermodiffusion from the gas state to the liquid state. In addition, regarding the isotope effects, it has been obtained that none of the available theoretical models are able to yield a reasonable prediction relatively to the molecular simulations results and that the moment of inertia contribution is one order of magnitude smaller than the mass contribution in the liquid state. Finally, concerning the chemical effects, it has been shown the Shukla and Firoozabadi model, complemented with a kinetic term, is probably the most reasonable option to estimate the chemical contribution to the thermal diffusion factor, even if it fails in capturing the effect of the asymmetry in size and in shape between the species. Overall, this works confirms that there is still a lack of a generic model able to predict accurately thermal diffusion factors, or equivalently Soret coefficient, in simple binary mixtures from the gas state to the liquid state., (© 2022. The Author(s), under exclusive licence to EDP Sciences, SIF and Springer-Verlag GmbH Germany, part of Springer Nature.)

We propose a simple and generic definition of a demarcation reconciling structural and dynamic frameworks when combined with the entropy scaling framework. This crossover line between gas- and liquid-like behaviors is defined as the curve for which an individual property, the contribution to viscosity due to molecules' translation, is exactly equal to a collective property, the contribution to viscosity due to molecular interactions. Such a definition is shown to be consistent with the one based on the minima of the kinematic viscosity. For the hard sphere, this is shown to be an exact solution. For Lennard-Jones spheres and dimers and for some simple real fluids, this relation holds very well. This crossover line passes nearby the critical point, and for all studied fluids, it is well captured by the critical excess entropy curve for atomic fluids, emphasizing the link between transport properties and local structure.

In this work, we propose a generic and simple definition of a line separating gas-like and liquid-like fluid behaviors from the standpoint of shear viscosity. This definition is valid even for fluids such as the hard sphere and the inverse power law that exhibit a unique fluid phase. We argue that this line is defined by the location of the minimum of the macroscopically scaled viscosity when plotted as a function of the excess entropy, which differs from the popular Widom lines. For hard sphere, Lennard-Jones, and inverse-power-law fluids, such a line is located at an excess entropy approximately equal to -2/3 times Boltzmann's constant and corresponds to points in the thermodynamic phase diagram for which the kinetic contribution to viscosity is approximately half of the total viscosity. For flexible Lennard-Jones chains, the excess entropy at the minimum is a linear function of the chain length. This definition opens a straightforward route to classify the dynamical behavior of fluids from a single thermodynamic quantity obtainable from high-accuracy thermodynamic models.

Petroleum reservoirs are remarkable illustrations of the impact of a thermal gradient on fluid pressure and composition. This topic has been extensively studied during the last decades to build tools that are required by reservoir engineers to populate their models. However, one can get only a very limited number of representative samples from a given reservoir and assessing connectivity between all sampling points is often a key issue. In some extreme cases, the whole reservoir fluid properties must be derived from a single point to define the field development plan. To do so, available models are usually not satisfactory as they need too many parameters and so cannot be considered as predictive tools. We propose in this work a comprehensive approach based on the irreversible thermodynamics principles to derive the relationships between pressure, chemical potentials and thermal gradients in porous media. It appears that there is no need for additional assumptions, it is just a matter of a making the right choices along theoretical developments. One of the most important steps is to express the full pressure gradient. As a final result, we obtain the chemical potential gradients for all components of the mixtures that can be easily translated in term of compositions through Equation of State modelling. The most important features of the final expressions are: i) the species relative separation in a thermal field is sensitive to the relative diffusion coefficients at stationary state. In porous media, the separation is sensitive to the permeability when the overall mobility is similar to diffusive mobility; ii) the magnitude of the separation depends on the residual entropy of the species; iii) the separation is not simply balanced by the average residual entropy. The balance is modified by the relative diffusion mobility of the components; iv) in low permeability porous media, the thermal gradient induces a pressure gradient proportional to the fluid residual entropy. As a validation, the proposed approach has been applied on a reservoir fluid subjected to a geothermal gradient and compared with non-equilibrium molecular dynamics simulation results at the stationary state.

Noble gases, and the way they fractionate, is a promising approach to better constrain origin, migration and initial state distributions of fluids in gas and oil reservoirs. Thermodiffusion, is one of the phenomena that may lead to isotope and elemental fractionation of noble gases. However, this effect, assumed to be small, has not been quantified, nor measured, in oil and gas under reservoir conditions. Thus, in this work, molecular dynamics simulations have been performed to compute the thermal diffusion factors of noble gases, in a dense gas (methane) and in an oil (n-hexane) under high pressures. Interestingly, it has been found that thermal diffusion factors, associated to both isotopic ( 36 Ar, 40 Ar) and elemental fractionations of noble gases ( 4 He, 20 Ne, 40 Ar, 84 Kr and 131 Xe) in gas and oil, could be expressed as linear functions of the reduced masses. Regarding the amplitude of the phenomena, it has been found that, in a stationary 1D oil or gas fluid column, thermodiffusion due to a typical geothermal gradient has an impact on noble gas isotopic and elemental fractionation which is of the same order of magnitude than gravity segregation, but opposite in sign. In addition, the relative impact of thermodiffusion on isotopic and elemental fractionations depends on the fluid type which is another interesting feature. Thus, these first numerical results on isotopic and elemental fractionation of noble gases by thermodiffusion in simple pure gas and oil emphasize their interest as natural tracers that could be used to improve the pre-exploitation description of oil and gas reservoirs.

With the advent of graphene material, membranes based on single-layer nanoporous solids appear as promising devices for fluid separation, be it liquid or gaseous mixtures. The design of such architectured porous materials would greatly benefit from accurate models that can predict their transport and separation properties. More specifically, there is no universal understanding of how parameters such as temperature, fluid loading conditions, or the ratio of the pore size to the fluid molecular diameter influence the permeation process. In this study, we address the problem of pure supercritical fluids diffusing through simplified models of single-layer porous materials. Basically, we investigate a toy model that consists of a single-layer lattice of Lennard-Jones interaction sites with a slit gap of controllable width. We performed extensive equilibrium and biased molecular dynamics simulations to document the physical mechanisms involved at the molecular scale. We propose a general constitutive equation for the diffusional transport coefficient derived from classical statistical mechanics and kinetic theory, which can be further simplified in the ideal gas limit. This transport coefficient relates the molecular flux to the fluid density jump across the single-layer membrane. It is found to be proportional to the accessible surface porosity of the single-layer porous solid and to a thermodynamic factor accounting for the inhomogeneity of the fluid close to the pore entrance. Both quantities directly depend on the potential of mean force that results from molecular interactions between solid and fluid atoms. Comparisons with the simulations data show that the kinetic model captures how narrowing the pore size below the fluid molecular diameter lowers dramatically the value of the transport coefficient. Furthermore, we demonstrate that our general constitutive equation allows for a consistent interpretation of the intricate effects of temperature and fluid loading conditions on the permeation process.

The computation of diffusion coefficients in molecular systems ranks among the most useful applications of equilibrium molecular dynamics simulations. However, when dealing with the problem of fluid diffusion through vanishingly thin interfaces, classical techniques are not applicable. This is because the volume of space in which molecules diffuse is ill-defined. In such conditions, non-equilibrium techniques allow for the computation of transport coefficients per unit interface width, but their weak point lies in their inability to isolate the contribution of the different physical mechanisms prone to impact the flux of permeating molecules. In this work, we propose a simple and accurate method to compute the diffusional transport coefficient of a pure fluid through a planar interface from equilibrium molecular dynamics simulations, in the form of a diffusion coefficient per unit interface width. In order to demonstrate its validity and accuracy, we apply our method to the case study of a dilute gas diffusing through a smoothly repulsive single-layer porous solid. We believe this complementary technique can benefit to the interpretation of the results obtained on single-layer membranes by means of complex non-equilibrium methods.

Compositional grading within a mixture has a strong impact on the evaluation of the pre-exploitation distribution of hydrocarbons in underground layers and sediments. Thermodiffusion, which leads to a partial diffusive separation of species in a mixture due to the geothermal gradient, is thought to play an important role in determining the distribution of species in a reservoir. However, despite recent progress, thermodiffusion is still difficult to measure and model in multicomponent mixtures. In this work, we report on experimental investigations of the thermodiffusion of multicomponent n -alkane mixtures at pressure above 30 MPa. The experiments have been conducted in space onboard the Shi Jian 10 spacecraft so as to isolate the studied phenomena from convection. For the two exploitable cells, containing a ternary liquid mixture and a condensate gas, measurements have shown that the lightest and heaviest species had a tendency to migrate, relatively to the rest of the species, to the hot and cold region, respectively. These trends have been confirmed by molecular dynamics simulations. The measured condensate gas data have been used to quantify the influence of thermodiffusion on the initial fluid distribution of an idealised one dimension reservoir. The results obtained indicate that thermodiffusion tends to noticeably counteract the influence of gravitational segregation on the vertical distribution of species, which could result in an unstable fluid column. This confirms that, in oil and gas reservoirs, the availability of thermodiffusion data for multicomponent mixtures is crucial for a correct evaluation of the initial state fluid distribution.

In this work, we have evaluated the applicability of the so-called thermodynamic scaling and the isomorph frame to describe the shear viscosity of Mie n-6 fluids of varying repulsive exponents (n = 8, 12, 18, 24, and 36). Furthermore, the effectiveness of the thermodynamic scaling to deal with binary mixtures of Mie n-6 fluids has been explored as well. To generate the viscosity database of these fluids, extensive non-equilibrium molecular dynamics simulations have been performed for various thermodynamic conditions. Then, a systematic approach has been used to determine the gamma exponent value (γ) characteristic of the thermodynamic scaling approach for each system. In addition, the applicability of the isomorph theory with a density dependent gamma has been confirmed in pure fluids. In both pure fluids and mixtures, it has been found that the thermodynamic scaling with a constant gamma is sufficient to correlate the viscosity data on a large range of thermodynamic conditions covering liquid and supercritical states as long as the density is not too high. Interestingly, it has been obtained that, in pure fluids, the value of γ is directly proportional to the repulsive exponent of the Mie potential. Finally, it has been found that the value of γ in mixtures can be deduced from those of the pure component using a simple logarithmic mixing rule.

In this article, the coupling between swelling and shear in liquid saturated slit nanopores is studied using molecular dynamics simulations on Lennard-Jones systems. First, the consistency of the simulations using thermodynamics and direct routes is validated when dealing separately with swelling and shear. Then, the coupling between swelling and shear is explored by displacing the solid walls in one direction while letting them move freely on the other. Results indicate that shear can induce swelling and vice versa because of the confined fluid phase structure. This phenomenon, which is neglected in poromechanics modeling, may be non-negligible in highly structured microporous systems, such as clays. It implies that the response to a variation in the external load can be a combination of volumetric and shear deformations, because of the fluid. Finally, we explore the behavior induced by solid walls moving at a constant velocity. Interestingly, when the wall velocity exceeds the swelling velocity, the instantaneous states of the system are no longer at equilibrium and the averaged pore width slightly increases with increasing shear rate.

The multi-component diffusive mass transport is generally quantified by means of the Maxwell-Stefan diffusion coefficients when using molecular simulations. These coefficients can be related to the Fick diffusion coefficients using the thermodynamic correction factor matrix, which requires to run several simulations to estimate all the elements of the matrix. In a recent work, Schnell et al. ["Thermodynamics of small systems embedded in a reservoir: A detailed analysis of finite size effects," Mol. Phys. 110, 1069-1079 (2012)] developed an approach to determine the full matrix of thermodynamic factors from a single simulation in bulk. This approach relies on finite size effects of small systems on the density fluctuations. We present here an extension of their work for inhomogeneous Lennard Jones fluids confined in slit pores. We first verified this extension by cross validating the results obtained from this approach with the results obtained from the simulated adsorption isotherms, which allows to determine the thermodynamic factor in porous medium. We then studied the effects of the pore width (from 1 to 15 molecular sizes), of the solid-fluid interaction potential (Lennard Jones 9-3, hard wall potential) and of the reduced fluid density (from 0.1 to 0.7 at a reduced temperature T* = 2) on the thermodynamic factor. The deviation of the thermodynamic factor compared to its equivalent bulk value decreases when increasing the pore width and becomes insignificant for reduced pore width above 15. We also found that the thermodynamic factor is sensitive to the magnitude of the fluid-fluid and solid-fluid interactions, which softens or exacerbates the density fluctuations.

This work aims at providing a tractable approach to model the local shear viscosity of strongly inhomogeneous dense fluids composed of spherical molecules, in which the density variations occur on molecular distance. The proposed scheme, which relies on the local density average model, has been applied to the quasi-hard-sphere, the Week-Chandler-Andersen and the Lennard-Jones fluids. A weight function has been developed to deal with the hard-sphere fluid given the specificities of momentum exchange. To extend the approach to the smoothly repulsive potential, we have taken into account that the non-local contributions to the viscosity due to the interactions of particles separated by a given distance are temperature dependent. Then, using a simple perturbation scheme, the approach is extended to the Lennard-Jones fluids. It is shown that the viscosity profiles of inhomogeneous dense fluids deduced from this approach are consistent with those directly computed by non-equilibrium molecular dynamics simulations.

This work focuses on a possible influence of a nanoporous medium on the thermodiffusion of a fluid "isotopic" mixture. To do so, we performed molecular dynamics simulations of confined Lennard-Jones binary equimolar mixtures using grand-canonical like and non-equilibrium approaches in sub- and super-critical conditions. The study was conducted in atomistic slit pore of three adsorbent natures for various widths (from 5 to 35 times the size of a molecule). The simulation results indicate that for all thermodynamic conditions and whatever the pore characteristics, the confinement has a negligible effect on the thermal diffusion factor/Soret coefficient. However, when considered separately, the mass diffusion and thermodiffusion coefficients have been found to be largely influenced by the pore characteristics. These two coefficients decrease noticeably when adsorption is stronger and pore width smaller, a behavior that is consistent with a simple hydrodynamic explanation.

In this work, we have considered and compared two molecular dynamics schemes widely used when studying a thin fluid film confined between solid surfaces and undergoing boundary shear. In the first approach, the non-equilibrium simulations are performed on a confined fluid explicitly connected to bulk reservoirs. In the second one, non-equilibrium simulations are carried out on the confined fluid only, in which the average density is deduced from a prior simulation in the grand canonical ensemble. We have found that the apparent properties (average density and effective viscosity) of a strongly confined Lennard-Jones liquid are significantly different using one scheme or the other when the solid surfaces induce a strong structure in the whole fluid, i.e., for small separations between the solid surfaces. Furthermore, the shear velocity dependence of the friction force has been found to be as well very sensitive to the approach chosen and can be well understood in terms of the fluid structure, which can even lead to a visco-plastic behavior of the fluid in some cases. Finally, it is shown that the first scheme is the only one usable to explore the history-dependence of the friction force as observed in experiments.

The zero-density shear viscosity of different types of short Lennard-Jones chains, up to the hexa-decamer, has been evaluated using a non-equilibrium molecular dynamics scheme. Simulations have been performed on chains of variable rigidities going from the fully flexible to the fully rigid chains. Very interestingly, it is found that there exists a universal relation (a power law) between the zero-density viscosity of the Lennard-Jones chains and their radius of gyration whatever the rigidity of the chain and for all tested temperatures (ranging from 2.5 to 6 in reduced units). Furthermore, for the studied range of temperature, it is shown that the zero-density viscosity of both fully flexible chains and fully rigid chains models can be obtained with an accuracy of a few percents knowing only the dimer viscosity and the length of the chain.

In this paper, molecular dynamics simulations of a simple Lennard-Jones fluid confined in narrow slit pores and undergoing shear have been performed. The aim is to investigate the effects of density inhomogeneities at the fluid-solid interfaces on the shear viscosity profiles. It has been found that the local viscosity was varying strongly with the distance from the solid walls for both dilute and dense fluid states with oscillations correlated to the density ones. To describe the computed viscosity profiles, we propose a scheme that uses the local average density model, combined with an adequate weight function, for the configurational viscosity and a semiempirical model for the translational viscosity. It is shown that the proposed approach is able to provide viscosity profiles in good agreement with those coming from simulations for different pore widths and for different fluid states (dilute to dense).

In this work, we describe two grand canonical-like molecular dynamics approaches to investigate mass diffusion phenomenon of a simple Lennard-Jones fluid confined between solid surfaces and in direct contact with reservoirs. In the first method, the density is used as the control variable in the reservoir whereas it is the pressure in the second method. Both methods provide consistent results, however, the constant density approach is the most efficient with respect to the computational time and implementation. Then, employing the constant density approach, we have studied the transient behavior of the diffusion process associated with the migration of one fluid into another one confined between parallel solid walls. Results have shown that the evolution of molar fraction of the invading fluid follows roughly a 1D diffusion model when the solid phase is weakly or moderately adsorbent with a characteristic time increasing when the pore width decreases. However, when the adsorption is high and the pore width small (i.e., below ten molecular sizes), the apparent mass diffusion in the adsorbed layer is reduced compared to that in the center of the slit pore. Hence, this mass diffusion process becomes a two-dimension phenomenon that must take into account an effective mass diffusion coefficient varying locally.

Using molecular dynamics simulations on inhomogeneous fluids, we have studied the effects of strong density inhomogeneities of varying wavelengths on the shear viscosity computed locally. For dense fluids, the local average density model combined with an adequate weight function yields a good description of the viscosity profiles obtained by simulations. However, for low density inhomogeneous fluids, the local average density model is unable to describe correctly the viscosity profiles obtained by simulations. It is shown that this weakness can be overcome by taking into account the density inhomogeneity in the local translational contribution to the viscosity using a density gradient like approach.

In this work, we have tested the efficiency of two scaling approaches aiming at relating shear viscosity to a single thermodynamic quantity in dense fluids, namely the excess entropy and the thermodynamic scaling methods. Using accurate databases, we have applied these approaches first to a model fluid, the flexible Lennard-Jones chain fluid (from the monomer to the hexadecamer), then to real fluids, such as argon and normal alkanes. To enlarge noticeably the range of thermodynamics conditions for which these scaling methods are applicable, we have shown that the use of the residual viscosity instead of the total viscosity is preferable in the scaling procedures. It has been found that both approaches, using the adequate scaling, are suitable for the Lennard-Jones chain fluid model for a wide range of thermodynamic conditions whatever the chain length when scaling law exponents and prefactors are adjusted for each chain length. Furthermore, these results were found to be well respected by the corresponding real fluids.

Molecular dynamics simulations of surface tensions of short flexible Lennard-Jones chains, composed of 2, 3, 4, and 5 segments, have been performed in this work. Using the simulation results, it is shown that the reduced surface tension depends only on the chain length and the reduced temperature. As a consequence, simple three parameters corresponding states using the acentric factor is shown to yield an excellent estimation of the reduced surface tension of the flexible Lennard-Jones chain fluid model. In addition, it has been noticed that the reduced surface tension of this fluid model is a unique function of the coexisting liquid and vapor reduced densities (i.e., there exist a universal Parachor behavior) for all chain lengths tested. When applied to real fluids, this universal behavior holds rather well for a large class of real species which can be nonspherical, nonlinear, and even polar. Only the surface tension of hydrogen-bonding compounds seems to largely deviate from this universal Parachor behavior. These interesting features of the surface tension, written in appropriate scaled forms, can probably be used to improve molecular models, in particular, those on which modern molecular based equations of state rely on.

Using nonequilibrium molecular dynamics simulations on simple Lennard-Jones binary mixtures, we have studied the behavior of planar fluid-fluid interfaces undergoing shear flow. When the miscibility is low enough, a slip together with a partial depletion have been noticed at the interface between the two fluid phases. The slip length can reach a value equal to some molecular diameters and the corresponding interfacial viscosity can be two times smaller than the value in the bulk. It is shown how the omission of this slip may lead to flow-rate misevaluation when dealing with a multiphase flow in a nanoporous medium even for non polymer fluids. In addition, using the simulation results, a simple relation between interfacial tension and interfacial viscosity is proposed for the monoatomic systems studied in this work. Finally, it is shown that the interfacial viscosity cannot be fully accounted for by estimating the local viscosity deduced from the local thermodynamic properties of the interface.

Nonequilibrium molecular dynamics simulations have been performed to estimate, analyze, and correlate the thermal conductivity of a fluid composed of short Lennard-Jones chains (up to 16 segments) over a large range of thermodynamic conditions. It is shown that the dilute gas contribution to the thermal conductivity decreases when the chain length increases for a given temperature. In dense states, simulation results indicate that the residual thermal conductivity of the monomer increases strongly with density, but is weakly dependent on the temperature. Compared to the monomer value, it has been noted that the residual thermal conductivity of the chain was slightly decreasing with its length. Using these results, an empirical relation, including a contribution due to the critical enhancement, is proposed to provide an accurate estimation of the thermal conductivity of the Lennard-Jones chain fluid model (up to 16 segments) over the domain 0.8

In a first part, interfacial properties of a pure monoatomic fluid interacting through the Mie n-6 potential (n=8, 10, 12, and 20) have been studied using extensive molecular simulations. Monte Carlo and molecular dynamics simulations have been employed, using, respectively, the test area approach and the mechanic route. In order to yield reference values, simulations have been performed with a cutoff radius equal to 10sigma, which is shown to be sufficient to avoid long range corrections. It is shown that both approaches provide results consistent with each other. Using the molecular simulations results, it is demonstrated that a unique scaling law is able to provide an accurate estimation of the surface tension whatever the repulsive exponent n, even far from the critical point. Furthermore, it is shown that the surface tension of the Mie n-6 fluid is as well accurately described by a unique Parachor's law. Density profiles are shown to be well represented by the tanh mean field profile, with slight deviations for the lowest temperatures and the smallest n. In addition, the interfacial width is shown to increase when n decreases (for a given reduced temperature) and to follow the usual scaling behavior for not too low temperature. In a second part, interfacial properties of the Mie n-6 fluid computed by the gradient theory, coupled with an equation of state based on the Barker-Henderson perturbation theory, have been compared with those obtained by molecular simulations. It is demonstrated that, even far from the critical point, the gradient theory is efficient to compute surface tensions and density profiles of this model fluid, provided the equation of state accurately model the phase behavior of the fluid involved (which is not the case for n=8 in this study).

Extensive nonequilibrium molecular dynamics (NEMD) simulations have been carried out in order to estimate the Newtonian shear viscosity of a fluid composed of short chains (up to 16 segments) of jointed spheres [Lennard-Jones chain (LJC)] over a large range of thermodynamic conditions. Using the NEMD results, it is shown that the zero-density contribution decreases with the chain length for a given temperature and is simply proportional to N-1/2 , where N is the number of spheres composing the chain. In addition, it has been noticed that the residual shear viscosity is proportional to the chain length. Then, using these results, a relation is proposed to correlate the shear viscosity of the LJC fluid using the LJ fluid (the monomer) as a reference. It is shown that this correlation is able to provide an excellent estimation of the LJC fluid viscosity compared to NEMD results for N

In this work, using extensive molecular dynamics simulations of several thermophysical properties, it is proposed to analyze possible relationships (in the corresponding state sense) between monoatomic fluids for which the repulsive interactions are modeled by an inverse n-power form, the Lennard-Jones 12-6 (LJ), or by an exponential one, the exponential-6 (Exp-6). To compare results between them, two possible definitions of Exp-6 potentials "equivalent" to the LJ one are proposed. In pure fluids, for a large range of thermodynamic conditions, the properties computed are the surface tension, liquid/vapor equilibrium densities, one-phase potential energy, pressure, isometric heat capacity, thermal pressure coefficient, self-diffusion, shear viscosity, and thermal conductivity. Additionally, thermodiffusion (Soret effect) has been considered in "isotopic" equimolar mixtures. It is shown that despite similarities exhibited by alike radial distribution functions, differences exist between the thermodynamic properties values provided by the LJ fluid and the two equivalent Exp-6 fluids. Nevertheless, quite surprisingly, when temperature and density are used as inputs, all three direct transport properties are shown to be nearly independent of the choice of the potential tested. Unexpectedly, these similarities hold even for thermodiffusion which is a priori very sensitive to the nature of the interactions. These results indicate that the use of an Exp-6 potential form to describe nonbonded/nonpolar interaction in molecular simulation is an alternative (more physically acceptable) to the LJ potential when dealing simultaneously with thermodynamic and transport properties. However, when only transport properties are considered (including thermodiffusion), the Exp-6 potential form should not lead to any differences compared to the LJ one.

Molecular dynamics (MD) simulations of direct and derivative thermodynamic properties of the Mie n-6 fluid (n=8, 10, and 12) have been performed for liquid to supercritical states. Using the results, an in depth test of the monomer-monomer interaction estimation of a recently derived statistical associating fluid theory of variable range (SAFT-VR) equation of state [Lafitte et al., J. Chem. Phys., 124, 024509 (2006)] has been carried out based on the Mie n-6 potential. For pure fluids, using an appropriate scaling, MD simulations show that density and isometric heat capacity are nearly independent of n, whereas sound velocity and thermal pressure coefficient tend to increase with n. In addition, the results show that predictions provided by the equation of state are consistent with those coming from MD and catch correctly the trends of each property with n except for the heat capacity. The comparison is next extended to binary mixtures with components differing only in the value of the n parameter and which demonstrate the reliability of the scheme (MX1b) used by Lafitte et al. to deal with this parameter in the SAFT-VR equation of state. In addition, a new empirical one-fluid approximation of the n parameter is proposed thanks to MD simulations, which very favorably compare with the one-fluid model on n previously proposed in the literature. The consistency of this approximation is addressed by making use of it in combination with the SAFT-VR Mie equation of state. It is shown that using such an approach, which is easier to handle than the MX1b one, yields slightly improved results compared to those of the MX1b.