用积分方程方法决定胶粒之间的空耗势

用硬球模中性胶体粒子 ,数值求解双组分Ornstein Zernike积分方程 (当大的中性胶体粒子的浓度为零时 ) ,用来决定悬浮在溶剂 (用小的硬球模拟 )中两个胶粒之间的空耗势 .所预言的空耗势与文献的模拟数据和实验数据能很好地符合 .研究发现 ,基于空耗势的有效一组分Hansen Verlet一相相变标准完全不能预言双组分系统的液 固相变 .讨论了导致这种现象的原因 :Hansen Verlet一相相变标准不能自地处理有效一组分系统中固相与液相的体积能     

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

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

铋的固相及液相多相状态方程研究

铋在高温高压下存在一系列复杂相结构,相变伴随着密度、内能等物理性质的改变.采用半经验三项式Helmholtz自由能表达式,构建了铋的五个固相及液相的多相状态方程,其中离子热振动自由能计算基于经典平均场理论给出的平均场势函数模型开展.研究结果表明,状态方程计算给出的铋的相图、等温压缩线、液相的温度-密度曲线以及冲击Hugoniot线等均与实验测量符合较好,故可认为本文构建的铋的多相状态方程具有良好的参数合理性以及模型适用性.     

接枝聚合物对小分子的选择性吸附研究

文章应用密度泛函理论研究接枝于壁面的方阱链对二元小分子混合物的选择性吸附特性. 系统的Helmholtz剩余自由能泛函被表示为硬球排斥和方阱吸引两部分贡献之和,分别由硬球链流体状态方程和变阱宽方阱链流体状态方程的简单加权密度近似来进行计算. 用此理论方法,分别考察了接枝聚合物的结构性质,以及不同温度下接枝分子层对二元方阱流体的选择性吸附性能. 结果表明:分子刷厚度随接枝密度线性增长而随温度非线性增加,并且在高温下趋于饱和;在较低温度下,接枝聚合物刷能表现出很好的选择性吸附能力,当聚合物刷被加热到高于饱和温度时,该能力将大幅度地减弱. 关键词： 密度泛函理论 接枝聚合物 选择性吸附 方阱链     

CO2气、液、固三相饱和态下pVT的计算

本文基于分子聚集理论和硬球引力场模型应用分子热力学方法导出CS-vdW状态方程,应用上述方程计算了CO2气-液和气-固平衡相的饱和态下pVT性质.取得了很满意的结果.本方法的优点是使用简便,可供工程计算使用.     

应用胞腔理论计算纯物质的汽化热

1引言胞腔模型或称自由体积模型早在本世纪30年代就已经提出［’－‘］，但由于该模型的重要参数，即自由体积一直无法予以确定，以至于就无法用该模型来计算物质的一系列性质、作者曾于1985年利用胞腔型状态方程导出两个物质特性参数一内压因子和构形因子问，本文根据上述状态方程和两个特性参数导出了计算物质的自由体积公式。从此即可应用自由体积参数计算物质的热力学性质。本文是介绍应用自由体积参数计算纯物质的汽化热的方法。2自由体积模型的液体理论和自由体积的计算自由体积理论最早由Lennard－Jones和Devon－shire用来将液体的…     

突扩圆管内液-固两相流固体颗粒运动特性的DPM数值模拟

采用考虑颗粒碰撞的欧拉一拉格朗日数值模拟方法(DPM),对水平突扩圆管中液固两相流固体颗粒的碰撞过程进行了数值计算.在模型中,对液相采用欧拉法建立控制方程,对离散颗粒采用拉格朗日方法模拟.采用硬球模型描述颗粒间的碰撞作用.计算结果表明,该模型可以真实地模拟液固两相流中固体颗粒运动的动态变化过程以及颗粒的非均匀分布特征,从单颗粒层次上提供颗粒的运动信息,这有助于深入研究液固两相流中固体颗粒的运动规律.     

PETN炸药爆轰产物状态方程的理论研究

 采用van der Waals等效单组分流体模型和Ross硬球微扰理论软球修正模型,计算爆轰气相产物的状态方程;用石墨相、金刚石相、类石墨液相和类金刚石液相4种相态描述凝聚成分,由Gibbs自由能最小原理确定不同状态下的凝聚产物相态。对爆轰产物混合系统采用自由能最小原理,通过化学平衡方程组求解炸药爆轰产物系统的平衡组分。使用该理论计算PETN炸药Chapman-Jouguet（CJ）点的爆轰参数,其值与实验值符合得很好;同时计算了以CJ点为起始点的等熵卸载线,并与传统的Jones-Wilkins-Lee（JWL）状态方程的计算结果进行比较,发现计算的γ值是单调递减的,而JWL状态方程计算的γ值却出现了“双峰”现象。分析认为,传统的JWL状态方程给出的“双峰”变化,是由其函数形式自身决定的,并不对应实际物理过程。     

基于统计物理的爆轰产物物态方程研究

用van der Waals等效单组分流体模型和Ross硬球微扰理论软球修正模型，计算爆轰气相产物的状态方程；用石墨相、金刚石相、类石墨液相和类金刚石液相４种相态描述凝聚成分，由Gibbs自由能最小确定不同状态下的凝聚产物相态.对爆轰产物混合系统采用Gibbs自由能最小原理，通过化学平衡方程组求解炸药爆轰产物系统的平衡组分，计算结果与Becker-Kistiakowsky-Wilson (BKW)和Lennard-Jones-Devonshire结果相近.使用该理论对炸药的爆轰参数做了预言，与BKW，Jo 关键词： 爆轰产物 物态方程 化学平衡方程组     

基于Helmholtz自由能模型的聚乙烯的完全物态方程

基于Helmholtz自由能建立了聚乙烯的完全物态方程,通过该模型计算获得了聚乙烯的150 GPa压力范围内的冲击Hugoniot关系、冲击波温度-压力关系,计算结果与已有实验结果和分子动力学计算结果均符合较好,表明构建的物态方程对描述聚乙烯离解相变压力150 GPa内的热力学量具有很好的适用性.     

扩展桥密度泛函近似于非匀一非硬球流体

为了将非匀一硬球流体的桥密度泛函近似扩展到非匀一非硬球流体,提出了一个理论方案．所得的LJ流体的密度泛函近似计算简单,精确．特别是密度泛函近似仅仅需要共存体相流体的二阶直接相关函数作为输入,因而可以应用于超临界与亚临界温度．所提出的理论方案可以认为是热力学理论的非匀一对应物．     

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.     

Structural and thermodynamic properties of a multicomponent freely jointed hard sphere multi-Yukawa chain fluid

The product-reactant Ornstein-Zernike approach, represented by the polymer mean-spherical approximation (PMSA), is utilized to describe the structure and thermodynamic properties of the fluid of Yukawa hard sphere chain molecules. An analytical solution of the PMSA for the most general case of the multicomponent freely jointed hard sphere multi-Yukawa chain fluid is presented. As in the case of the regular MSA for the hard sphere Yukawa fluid, the problem is reduced to the solution of a set of nonlinear algebraic equations in the general case, and to a single equation in the case of the factorizable Yukawa potential coefficients. Closed form analytical expressions are presented for the contact values of the monomer-monomer radial distribution function, structure factors, internal energy, Helmholtz free energy, chemical potentials and pressure in terms of the quantities, which follows directly from the PMSA solution. By way of illustration, several different versions of the hard sphere Yukawa chain model are considered, represented by one-Yukawa chains of length m, where m = 2, 4, 8, 16. To validate the accuracy of the present theory, Monte Carlo simulations were carried out and the results are compared systematically with the theoretical results for the structure and thermodynamic properties of the system at hand. In general it is found that the theory performs very well, thus providing an analytical route to the equilibrium properties of a well defined model for chain fluids.     

Theoretical Investigation of Uniform and Non-uniform Penetrable Sphere Fluid

A bridge function approximation is proposed for a single-component fluid consisting of penetrable sphere interacting via a potential that remains finite and constant for center-center distance smaller than the particle diameter and is zero otherwise. The radial distribution function from the Ornstein-Zernike integral equation combined with the present bridge function approximation is in satisfactory agreement with the corresponding simulation data for all of the investigated state points. The presently calculated excess Helmholtz free energy respectively based on virial route and compressibility route is highly self-consistent, and is in very good agreement with simulational results for the case of low temperatures. The present bridge function approximation, combined with the bridge density functional approximation, can reproduce very accurately density profiles of the penetrable sphere fluid confined in a hard spherical cavity for all the cases where simulational results are available.     

Theoretical Investigation of Uniform and Non-uniform Penetrable Sphere Fluid

A bridge function approximation is proposed for a single-component fluid consisting of penetrable sphere interacting via a potential that remains finite and constant for center-center distance smaller than the particle diameter and is zero otherwise. The radial distribution function from the Ornstein-Zernike integral equation combined with the present bridge function approximation is in satisfactory agreement with the corresponding simulation data for all of the investigated state points. The presently calculated excess Helmholtz free energy respectively based on virial route and compressibility route is highly self-consistent, and is in very good agreement with simulational results for the case of low temperatures. The present bridge function approximation, combined with the bridge density functional approximation,can reproduce very accurately density profiles of the penetrable sphere fluid confined in a hard spherical cavity for all the cases where simulational results are available.     

TPT2 and SAFTD equations of state for mixtures of hard chain copolymers

We present the second-order thermodynamic perturbation theory (TPT2) and the dimer statistical associating fluid theory (SAFTD) equations of state for mixtures consisting of hetero-nuclear hard chain molecules based on extensions of Wertheim's theory for associating fluids. The second-order perturbation theory, TPT2, is based on the hard sphere mixture reference fluid. SAFTD is an extension of TPT1 (= SAFT) and is based on the non-spherical (hard disphere mixture) reference fluid. The TPT2 equation of state requires only the contact values of the hard sphere mixture site-site correlation functions, while the SAFTD equation of state requires the contact values of site-site correlation functions of both hard sphere and hard disphere mixtures. We test several approximations for site-site correlation functions of hard disphere mixtures and use these in the SAFTD equation of state to predict the compressibility factor of copolymers. Since simulation data are available only for a few pure copolymer systems, theoretical predictions are compared with molecular simulation results for the compressibility factor of pure hard chain copolymer systems. Our comparisons show a very good performance of TPT2, which is found to be more accurate than TPT1 (= SAFT). Using a modified Percus-Yevick site-site correlation function SAFTD is found to represent a significant improvement over SAFT and is slightly more accurate than TPT2. Comparison of SAFTD with generalized Flory dimer (GFD) theory shows that both are equivalent at intermediate to high densities for the compressibility factor of copolymer systems investigated here.     

The role of higher-order terms in perturbation approaches to the monomer and bonding contributions in a SAFT-type equation of state for square-well chain fluids

A perturbation theory for square-well chain fluids is developed within the scheme of the (generalised) Wertheim thermodynamic perturbation theory. The theory is based on the Pavlyukhin parametrisations [Y. T. Pavlyukhin, J. Struct. Chem. 53, 476 (2012)] of their simulation data for the first four perturbation terms in the high temperature expansion of the Helmholtz free energy of square-well monomer fluids combined with a second-order perturbation theory for the contact value of the radial distribution function of the square-well monomer fluid that enters into bonding contribution. To obtain the latter perturbation terms, we have performed computer simulations in the hard-sphere reference system. The importance of the perturbation terms beyond the second-order one for the monomer fluid and of the approximations of different orders in the bonding contribution for the chain fluids in the predicted equation of state, excess energy and liquid–vapour coexistence densities is analysed.     

Statistical Mechanics Approach for Uniform and Non-uniform Fluid with Hard Core and Interaction Tail

One recently proposed self-consistent hard sphere bridge functional was combined with an exponential function exp(-cr) and a re-normalized indirect correlation function to construct the bridge function for fluid with hard core and interaction tail. In the present approach, the adjustable parameter α was determined by the thermodynamic consistency realized on the compressibility modulus, the re-normalization of the indirect correlation function was realized by a modified Mayer function with the interaction potential replaced by the perturbative part of the interaction potential. As an example, the present bridge function was combined with the Ornstein-Zernike (OZ) equation to predict structure and thermodynamics properties in very good agreement with the simulation data available for Lennard-Jones (L J). Based on the universality principle of the free energy density functional and the test particle trick, the numerical solution of the OZ equation was employed to construct the first order direct correlation function of the non-uniform fluid as a functional of the density distribution by means of the indirect correlation function. In the framework of the density functional theory, the numerically obtained functional predicted the density distribution of LJ fluid confined in two planar hard walls that is in good agreement with the simulation data.