Perturbation and variational approach for the equation of state for hard-sphere and Lennard-Jones fluids

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), Lennard–Jones (LJ) 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.     

A comprehensive approach to an equation of state for hard spheres and Lennard-Jones fluids

We present a simple method of obtaining various equations of state for hard sphere fluid in a simple unifying way.We will guess equations of state by using suitable axiomatic functional forms (n=1,2,3,4,5) for surface tension S n m (r),r ≥ d/2 with intermolecular separation r as a variable,where m is an arbitrary real number (pole).Among the equations of state obtained in this way are Percus-Yevick,scaled particle theory and Carnahan-Starling equations of state.In addition,we have found a simple equation of state for the hard sphere fluid in the region that represents the simulation data accurately.It is found that for both hard sphere fluids as well as Lennard-Jones fluids,with m=3/4 the derived equation of state (EOS) gives results which are in good agreement with computer simulation results.Furthermore,this equation of state gives the Percus-Yevick (pressure) EOS for the m=0,the Carnahan-Starling EOS for m=4/5,while for the value of m=1 it corresponds to a scaled particle theory EOS.     

Equation of state for pure fluids at critical temperature

In this paper, we employ the concept of probability for creating a cavity with diameter d in fluid along with the perturbation and variation approach, and develop an equation of state (EOS) for a hard sphere (HS) and Lennard– Jones (LJ) fluids. A suitable axiomatic form for surface tension S(r) is assumed for the pure fluid, with r as a variable. The function S(r) has an arbitrary parameter m. S(r) = A + B(d/r)/[1 + m(d/r)]. We use the condition in terms of radial distribution function G(λd, η) containing the self-consistent parameter λ and the condition of continuity at r = d/2 to determine A and B. A different EOS can be obtained with a suitable choice of m and the EOS has a lower root-mean-square deviation than that of Barker–Henderson BH2 for LJ fluids.     

Critical anomaly and finite size scaling of the self-diffusion coefficient for Lennardben Jones fluids by non-equilibrium molecular dynamic simulation

We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard-Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with temperature decreases by increasing density. For density ρ^* = ρσ³ = 0.84 we observe a peak at the value of the self-diffusion coefficient and the critical temperature T^* = kT/ε = 1.25. The value of the self-diffusion coefficient strongly depends on system size. The data of the self-diffusion coefficient are fitted to a simple analytic relation based on hydrodynamic arguments. This correction scales as N^-α, where α is an adjustable parameter and N is the number of particles. It is observed that the values of α < 1 provide quite a good correction to the simulation data. The system size dependence is very strong for lower densities, but it is not as strong for higher densities. The self-diffusion coefficient calculated with non-equilibrium molecular dynamic simulations at different temperatures and densities is in good agreement with other calculations from the literature.     

Computation of high-order virial coefficients in high-dimensional hard-sphere fluids by Mayer sampling

The Mayer sampling method was used to compute the virial coefficients of high-dimensional hard-sphere fluids. The first 64 virial coefficients for dimensions 12 < D ? 100 were obtained to high precision, and several lower dimensional virial coefficients were computed. The radii of convergence of the virial series in 13, 15, 17 and 19 dimensions agreed well with the analytical results from the Percus–Yevick closure.     

Thermodynamic study of fluid in terms of equation of state containing physical parameters

We introduce a simple condition for one mole fluid by considering the thermodynamics of molecules pointing towards the effective potential for the cluster.Efforts are made to estimate new physical parameter f in liquid state using the equation of state containing only two physical parameters such as the hard sphere diameter and binding energy.The temperature dependence of the structural properties and the thermodynamic behavior of the clusters are studied.Computations based on f predict the variation of numbers of particles at the contact point of the molecular cavity(radial distribution function).From the thermodynamic profile of the fluid,the model results are discussed in terms of the cavity due to the closed surface along with suitable energy.The present calculation is based upon the sample thermodynamic data for n-hexanol,such as the ultrasonic wave,density,volume expansion coefficient,and ratio of specific heat in the liquid state,and it is consistent with the thermodynamic relations containing physical parameters such as size and energy.Since the data is restricted to n-hexanol,we avoid giving the physical meaning of f,which is the key parameter studied in the present work.     

一种简化维里型状态方程预测高温甲烷PVT关系

为满足描述爆轰环境下高温气体高温、中高压状态的需求,本文提出了一种基于Lennard-Jones(LJ)势能函数的简化维里型状态方程Han-Long(HL).应用HL状态方程计算了甲烷1000 K以上112组理论和实验数据,计算所得体积平均绝对偏差约为1%,最大误差为3.28%,远低于DMW状态方程和BS状态方程的计算偏差.采用HL状态方程计算了甲烷冲击试验的热力学数据,计算所得体积偏差均小于3%.结果表明,HL状态方程能够很好的描述高温甲烷的热力学状态.     

An equation of state contribution for dipolar and quadrupolar square-well fluids

A dipolar–quadrupolar contribution to the residual Helmholtz energy for a polar square well (a square well plus either a point dipole or a point quadrupole) fluid is developed based on the Padé approximation. Taking the square well system as reference, the contribution is formulated using an expansion for radial distribution function of the reference system. In addition to square well potential parameters the contribution depends only on dipole and quadrupole moments. This term is added as perturbation to a generalized equation of state for square well fluids. The results are then compared with the available simulation data in the literature. With the new equation obtained, it was possible to predict liquid–vapour equilibrium properties and critical properties of polar square well fluids more accurately than with available perturbation theories for multipolar square well systems. Application of the equation of state to a real dipolar (water) and a real quadrupolar (carbon dioxide) fluid indicated that the polar contribution greatly improved the predictions of saturation properties. Accurate prediction of critical properties for polar square well fluids remains as a challenge. This work can be useful in the development of better equations of state.     

混合物状态方程的计算

多介质流体动力学过程的数值模拟往往涉及混合物状态方程的计算. 做图法和Newton 法是混合物状态方程计算常采用的方法, 前者虽直观精度却差, 后者计算效率高却只具有局部收敛性, 当解与其初始猜测值相差较远时Newton法不一定能够获得收敛解. 为此, 本文给出一种具有大范围收敛性的嵌入算法(imbedding method)求解混合物状态方程, 其基本思想是通过引入嵌入参数, 将待解的混合物状态方程和易解的混合物状态方程线性组合, 构成嵌入方程组, 当嵌入参数从0连续地变化到1 时, 嵌入方程组的解由易解的混合物状态方程的解连续地变化为待解的混合物状态方程的解. 嵌入方程组可由Newton法迭代求解, 也可转化为以嵌入参数为自变量的常微分方程组, 从而易于由成熟的计算方法如梯形法等进行求解. 进一步利用热力学基本关系, Maxwell形式的微分方程描述了压力和温度随嵌入参数的演化速率与应变速率和组分质量分数演化速率的关系. 对铅锡混合物热力学量的计算表明了本文算法的有效性.     

一种计算氩等离子物态方程的简单模型

将Thomas-Feimi统计模型电离势的数值结果进行函数逼近，给出一个便于近似数值求解的解析表达式和计算电离度的近似计算方法，计算了Ar元素LTE情况下的电离度和物态方程，结果与Saha模型的计算结果和实验结果符合较好.所提出的简单模型也适用于计算混合物物态方程，可以在电磁发射技术领域中的强电离等离子体中有广阔的应用前景. 关键词： 等离子体 物态方程 电离势 电离度     

Chemical potentials and phase equilibria of Lennard-Jones chain fluids

Computer simulations (molecular dynamics) were performed for ensembles of flexible tangent Lennard-Jones chains consisting of n sites (n = 1, 2, 4, 8, and 16). From these simulations, the orthobaric liquid and vapour densities were calculated not only with the traditional method of simulating a liquid film in coexistence with vapour, but also using the rigorous thermodynamic condition of satisfying the chemical potential equality between the phases in equilibrium. The agreement with literature data, as far as such exist, is excellent.     

H-theorem for the hard-sphere gas

The approach to equilibrium of the hard-sphere gas is discussed from the master-equation point of view. AnH-theorem is established, which is valid for arbitrary initial conditions.     

Universal cubic equation of state and contact values of the radial distribution functions for multi-component additive hard-sphere mixtures

A universal cubic equation of state (UC EOS) is proposed based on a modification of the virial Percus-Yevick (PY) integral equation EOS for hard-sphere fluid. The UC EOS is extended to multi-component hard-sphere mixtures based on a modification of Lebowitz solution of PY equation for hard-sphere mixtures. And expressions of the radial distribution functions at contact (RDFC) are improved with the form as simple as the original one. The numerical results for the compressibility factor and RDFC are in good agreement with the simulation results. The average errors of the compressibility factor relative to MC data are 3.40%, 1.84% and 0.92% for CP3P, BMCSL equations and UC EOS, respectively. The UC EOS is a unique cubic one with satisfactory precision among many EOSs in the literature both for pure and mixture fluids of hard spheres.     

混合物物态方程的计算

在采用体积相加原理计算混合物物态方程的基础上,建立了一种物理模型确定混合物温度。根据混合物中各组分温度和压强平衡条件,采用压强-密度迭代方法计算给出混合物物态方程,编制了两种组分的混合物物态方程计算程序。为检验建立的温度模型的合理性及程序的有效性,分析了不同密度、温度状态的氢（H2）和钨（W）组成的混合物状态参量,计算了以下情形及其组合情形的混合物物态方程：H2和W以不同质量比混合;质量比固定,单组分状态不同;温度区间和密度区间不同。研究表明:实际应用中在建立的混合物温度模型基础上确定的混合物物态方程是合理的。     

The infinite-volume ground state of the Lennard-Jones potential

We consider a finite chain of particles in one dimension, interacting through the Lennard-Jones potential. We prove the ground state is unique, and approaches uniform spacing in the infinite-particle limit.Research supported in part by NSF grant MCS 78-01520.     

A simple and accurate equation of state for hard-disk fluids

A long-standing problem for hard-disk fluids is to find a simple and accurate equation of state. By observing the known updated virial coefficients from up to , we obtain an accurate empirical formula, with errors less than 0.5%. We further assume that the empirical formula is valid for all the virial coefficients. Using this assumption, we obtain a simple and accurate equation of state. In the density range 0<ρa²≤0.83, the predicted pressures are in good agreement with the simulation results. For 0<ρa²≤0.70, the errors are less than 0.007%. Even at ρa²=0.83, the error is only 0.5%.     

Determination of the structure factor for the rare gas fluids based on the ORPA theory and LIR equation of state

Linear isotherm regularity is applied to generate an expression for the direct correlation function (DCF), based on the optimized random-phase approximation theory. We use the Lennard-Jones potential in the modelling of real fluids, in which its molecular parameters are state-dependent. Based on the perturbation theory, we assume that the core contribution of the DCF is related to the geometric effect via a linear form, which arises from the excluded volume, and the tail contribution is related to the long-range intermolecular interactions of the system via a non-linear form. We use this model to predict the behaviour of the structure factor, S(k), in a wide range of k values for the xenon and krypton liquids. Finally, we compare our results with the experimental and theoretical data available in literature. The model is also successful in the presentation of the Ornstein–Zernike behaviour of S(k) at the low-k region.     

混合物物态方程的计算

 在采用体积相加原理计算混合物物态方程的基础上,建立了一种物理模型确定混合物温度。根据混合物中各组分温度和压强平衡条件,采用压强-密度迭代方法计算给出混合物物态方程,编制了两种组分的混合物物态方程计算程序。为检验建立的温度模型的合理性及程序的有效性,分析了不同密度、温度状态的氢（H2）和钨（W）组成的混合物状态参量,计算了以下情形及其组合情形的混合物物态方程：H₂和W以不同质量比混合;质量比固定,单组分状态不同;温度区间和密度区间不同。研究表明:实际应用中在建立的混合物温度模型基础上确定的混合物物态方程是合理的。     

统一状态方程描述同位素气体

通过对分子动力学模拟公式和同位素相互作用势特点的分析,提出用一个统一的状态方程描述同位素气体的P-V-T特性.进而选用Benedict-Webb-Rubin方程作为统一的状态方程,利用H2气体的135组实验值确定其八个参数.这个方程的计算结果与H2和D2气体的实验值符合良好.T2气体目前尚未见到成套的P-V-T实验数据报道,不能直接与实验结果对比,但与T2气体的分子动力学模拟计算结果一致.     

用多次冲击压缩方法研究稠密氢氦等摩尔混合气体的物态方程

采用液氮冷却以及气体增压技术，制备了～30MPa及～90K初始状态的氢、氦等摩尔混合气体样品.以二级轻气炮作为加载工具，用不同灵敏度设置的两套多通道瞬态高温计系统获得完整、清晰的稠密氢、氦混合气体多次冲击压缩过程的光谱辐射强度信号.并建立起相应的实验数据处理和分析技术，获得了5—140GPa范围内氢、氦混合气体一至五次冲击雨贡纽物态方程，以及一次、二次和四次冲击温度实验数据.流体变分理论和离解模型用来分析和解释所获得的测量结果. 关键词： 氢、氦混合气体 多次冲击压缩 光谱辐射强度历史 物态方程