气液态分子运动的Hurst指数与相图研究.pdfVIP

中国I程热物理年会 I群热力学与能源利埘 2∞4年学术会议 编号：04103} 气液态分子运动的Hurst指数及相图研究 万黎明刘朝杨智勇张新铭曾丹苓 萝庆大学动力r程率院重庆400044 Tel：(023 E-mail：liuchao@cta．cq．ca 摘要：本文运用分了动力学模拟方法，模拟j-小刷温度条件下饱和气体和液体氩分予的运动过程。对 小川状态下分子运动在速度空阃的演化耜尉进f于丁比较。根据分形理论，对实际分了的分数布朗运动 进行九寸论。应用本文提出的约化趣模法，对分了分数布朝运动的Hut-st指数和轨迹分维数进行r数 值统计计算。 关键词： 分了动力学模拟轨迹分维数Huru指数相幽 1引言 分形几何是研究自然界中常见的、不规则的形体表现及其内在规律的非线性科学。 其研究对象是具有“自相似性”、“自仿射性”的分形体，其分形度量为分维数。分维数 的整数部分代表其拓扑维数．小数部分代表分形体的复杂程度和其充满空间的能力111。 分子动力学模拟(MD)方法的基本原理是，在～定宏观约菜条件下通过求解系统中分 j二的运动，获得分子运动的轨迹、位置和速度随时问演化的信息规律，模拟系统随时问 演变的微观过程，并统计计算系统的平衡状态参数和输运性质吲。随着计算机运算能力 的提高，MD模拟在统计物理等方面得到了迅速发展和广泛应用。 2分子动力学模拟方法 本文模拟对象为氩．势能函数采用L-J型双体势能函数1w： u(r)=46k，ry2一(o，r户I (1) 模拟体系由1372个氩分子构成，在初始时刻气液相密度分别为模拟温度对应的气液相 饱和值．呈FCC面心立方晶格结构分布于模拟体系中间。模拟体系为NvT正则系综， 采用周期边界条件，运动方程用Verlet算法计算，模拟计算的无因次时阔步长0．02，无 因次截断半径3．5。 3不同温度下饱和态氯分子速度相图 速度相图反映的是分子速度在速度空间上随时间的演化过程。本文对温度为84、95、 105、115、125、135和145KF平衡的饱和气体及液体氧系统进行了模拟。 图1为T=$4K，一个液态氩分子在14万模拟步内的速度随时间演化相图。模拟发 40 观，不同温度F液态分子速度相幽基本扫同，与幽1所示类似。均儿中间密集的对称雪 球形。液态分子速度相图的特点是：密集性和连贯性。这是由于液体内部分子密集，相 互作用强而产生的结果。 闰2为T=145K时。一个气态氯分子在14万模拟步内的速度随时间演化相图。不同 温度卜_气态分子速度相图基本相同。呈中间密集的对称雪球形。与液态比较，较为稀疏。 在温度较低时，这种稀疏性更明显．如图3所示。气态分子速度相图的特点：稀疏性和 跷跃性，这两个特性在湿度较低对会更为明显。这是因为在气体分子内部。分子间的相 互作用较弱。 崮284K下的液态氨分子速度相I皇l 圈2 145K下的气态氩分子速度相图 瞬3 90K下的气杏董分子速度相圈 图4a84K下的气态氩分子速度相图 在温度接近与三相点时(氩的三相点温度为83．78K)，速度相图出现了不同的特点， 如图4a所示．在相图中存在奇异吸引子。速度相图上并非在整个时间度都有奇异吸引 的，仅出现在一定模拟时间段上。这表明在该段时间内，气体分子的运动具有与通常气 体分子不同的运动形态．但并不是每个分子的相圈都会出现奇异吸引子，如图4b所示。 某个分子在捆图上一定时间出现的吸引子，反映出此时分子运动处于与通常气体分 子的不同运动状态。而吸引子内速度在速度空间上随时间的演化方式与气体有很大的不 同。吸引子的演化方式更接近于液态分子的演化方式，是密集和连贯的。从而可以推断： 在微观上，在接近三相点附近的气体分子．并不是以气态的形式稳定存在的，它可以处 下液态的运动形态。或者说．在三相点附近的气体，会在其内部形成液态微团。当分子