气体-表面相互作用的分子动力学模拟研究

采用分子动力学模拟方法研究了气体分子Ar在光滑和粗糙Pt表面上的散射规律.提出了一种速度抽样方法,计算了不同温度条件下气体分子对光滑和粗糙表面的切向动量适应系数和吸附概率.结果显示：光滑表面条件下,气体分子的切向动量系数和吸附概率都随着温度的升高而降低;粗糙度对气体分子切向动量与表面的适应具有极大的促进作用,当粗糙度足够大时,切向动量适应系数的大小趋近于1.0,对温度的敏感性也逐渐降低.采用粒子束方法对气体分子在光滑和粗糙表面上的散射规律进行了定量分析.总结了散射过程中气体分子的典型轨迹和动量变化规律,将气体分子在光滑表面的散射分为两种类型：单次碰撞后散射和多次碰撞后散射.单次碰撞后散射的气体分子平均切向动量有所减小,而经过多次碰撞后散射的气体分子则倾向于保持原有的平均切向动量.对于粗糙表面,粗糙度的存在使气体分子与表面间的动量和能量适应更加充分,导致气体分子在较粗糙表面上散射后的平均切向动量大幅减小并接近于0,且气体分子在表面上经历的碰撞次数越多,其散射后的能量损失越严重.     

双甘氨肽与高定向热解石墨的散射动力学研究

本文通过反应力场来研究中性双甘氨肽与高取向热解石墨的碰撞动力学过程，分别模拟初始入射角度为0o、20o、45o和70o以及入射能量为481.5 kJ/mol和表面温度为677 K的情况，并且确定和分析了散射产物的角度分布、平动能分布、内能分布和表面滞留时间分布. 双甘氨肽是一种多原子分子，具有较多的低频振动模式和较强的表面吸附相互作用，这些使得碰撞过程的表面滞留时间较长和容易造成能量损失，尤其是表面法线方向. 由于碰撞分子的动量沿表面法线方向上明显地发生损失，而沿表面平行方向上的动量则会大部分保留，因此，其散射角通常会呈现出超镜面反射分布，并且末态的平动能要远小于所谓硬立方模型的预测值. 本文加深了多肽分子与高取向热解石墨碰撞及其能量转移过程的认识和理解，有助于设计在稀薄大气中收集这类大分子的中性气体浓缩器.     

质子与羟基碰撞的含时密度泛函理论研究

采用将含时密度泛函理论和分子动力学非绝热耦合的方法,研究了不同入射速度下质子与羟基碰撞的反应动力学.计算了碰撞前后质子动能和羟基动能的变化及羟基电子和质子的运动.计算结果表明,质子沿垂直羟基分子轴方向入射时,质子与羟基碰撞后,质子被反弹且动能损失并俘获了羟基中氧的一部分电子,而丢失部分电子的羟基则获得动能以伸缩振动的形式向计算边界平动.随着入射质子的初动能增加,质子从羟基中俘获的电子增多,碰撞后羟基的键长变长,羟基振动变强而伸缩振动频率降低.此外,还发现质子的入射方向对碰撞过程的激发动力学有很大的影响.质子从不同的方向入射时,质子的入射初动能越大,其损失的动能越多且损失的动能与入射初动能呈线性关系,而入射方向对质子动能损失的影响很小.在质子入射初动能较低(小于25 eV)的情况下,羟基获得的动能与质子入射初动能呈线性关系且与入射方向无关;在质子入射初动能较高(大于25 eV)时,当质子沿羟基分子轴方向入射时,羟基动能的增量远大于质子沿垂直于羟基分子轴方向入射时羟基动能的增量.     

壁面效应对纳米尺度气体流动的影响规律研究

采用分子动力学方法研究了过渡区纳米通道内的壁面力场对气体剪切流动的影响规律.在纳米尺度下,壁面力场对流场的主导作用更加显著,流动物理量对于壁面条件和系统温度的变化也更加敏感.壁面原子的运动采用Einstein模型模拟,结果表明随着壁面刚度的增加,气体在近壁面区域的速度峰值减小,气体分子与壁面的动量适应性变差.壁面粗糙度通过金字塔形模型来研究,发现无论是主流区域还是近壁区域,壁面粗糙度对流动的影响都非常明显.当粗糙单元高度增大时,气体分子在壁面处的聚集现象明显,与壁面完全动量适应.本文还研究了系统温度对纳米通道流动的影响,结果表明温度的影响是全局性的,温度的升高导致整个通道内流速降低,近壁区域气体密度减小,气-固动量适应性变差.     

双甘氨肽与高定向热解石墨的散射动力学研究（英文）

本文通过反应力场来研究中性双甘氨肽与高取向热解石墨的碰撞动力学过程,分别模拟初始入射角度为0°、20°、45°和70°以及入射能量为481.5 kJ/mol和表面温度为677 K的情况,并且确定和分析了散射产物的角度分布、平动能分布、内能分布和表面滞留时间分布.双甘氨肽是一种多原子分子,具有较多的低频振动模式和较强的表面吸附相互作用,这些使得碰撞过程的表面滞留时间较长和容易造成能量损失,尤其是表面法线方向.由于碰撞分子的动量沿表面法线方向上明显地发生损失,而沿表面平行方向上的动量则会大部分保留,因此,其散射角通常会呈现出超镜面反射分布,并且末态的平动能要远小于所谓硬立方模型的预测值.本文加深了多肽分子与高取向热解石墨碰撞及其能量转移过程的认识和理解,有助于设计在稀薄大气中收集这类大分子的中性气体浓缩器.     

质子与乙烯分子碰撞的密度泛函理论研究(英文)

采用含时局域密度近似与分子动力学相结合的方法研究了不同入射速度的质子与乙烯分子碰撞的动力学。计算了质子的能量损失及碰撞后乙烯分子的电子和离子的运动状态, 研究了质子的入射方向及入射动能对整个系统的碰撞动力学的影响。计算结果表明, 当入射质子的动能较小（E_k0<250 eV）时, 在相同的入射速度下, 当质子垂直于分子平面入射时, 系统的电离最大, 质子俘获的电子多; 当质子的入射动能E_k0>250 eV时, 质子的能量损失与入射方向有密切的关系。In the framework of the time dependent local density approximation (TDLDA),which applied to valence electrons, coupled non adiabatically to molecular dynamics of ions, the microscopic mechanisms of collisions between energetic protons and ethylene are studied. Not only the amount of energy lost of the projectile, but also the electron and vibration excitations of the target are identified. In addition, the influences of the collision orientation on the energy loss of the proton and excitation dynamics of ethylene are discussed. It is found that the ionization is enhanced and more electrons are captured by the proton when the proton with the impact energy less than 250 eV moves perpendicularly to the molecular plane. A strong relation between the proton energy lost and the impact orientation is obtained when the impact energy is larger than 250 eV.     

微纳气体流动中动量协调系数的分子动力学模拟

协调系数反映气体与壁面动量和能量交换规律,获得准确的协调系数以及明晰影响协调系数的因素十分重要.使用分子动力学方法提出了一种协调系数的统计方法,在等温泊肃叶流动中,探讨了不同温度及流速下切向和法向动量协调系数的变化规律.结果表明,切向动量协调系数对温度敏感,随着温度的上升其值逐渐减小并趋向于稳定,而流速不同时变化不大;法向动量协调系数受温度和流速影响很小,在等温流动中趋向于1.     

纳米通道内气体剪切流动的分子动力学模拟

采用分子动力学模拟方法研究了表面力场对纳米通道内气体剪切流动的影响规律.结果显示通道内的气体流动分为两个区域:受壁面力场影响的近壁区域和不受壁面力场影响的主流区域.近壁区域内,气体流动特性和气体动力学理论预测差别很大,密度和速度急剧增大并出现峰值,正应力变化剧烈且各向异性,剪切应力在距壁面一个分子直径处出现突变.主流区域的气体流动特性与气体动力学理论预测相符合,该区域内的密度、正应力与剪切应力均为恒定值,速度分布亦符合应力-应变的线性响应关系.不同通道高度及密度下,近壁区域的归一化密度、速度及应力分布一致,表明近壁区域的气体流动特性仅由壁面力场所决定.随着壁面对气体分子势能作用的增强,气体分子在近壁区域的密度和速度随之增大,直至形成吸附层,导致速度滑移消失.通过剪切应力与切向动量适应系数(TMAC)的关系,得到不同壁面势能作用下的TMAC值,结果表明壁面对气体分子的势能作用越强,气体分子越容易在壁面发生漫反射.     

二维分子动力学计算中协调系数的表达式

协调系数是反映气体与壁面动量和能量交换的重要参数,可通过分子动力学方法进行统计,其二维模拟由于计算量较三维大大减小而逐渐得到应用.法向动量协调系数表达式中的p_(nw)和能量协调系数表达式中的E_w分别表示气体分子在壁面漫反射后的平均法向动量和平均能量,利用气体动理论推导发现,二维的p_(nw)与三维相同,而E_w则由于降维小于三维的结果,导致能量协调系数的表达式在二维计算中有别于三维.依据协调系数的表达式,使用分子动力学方法模拟三维和二维系统中的导热问题,结果表明,二维能量协调系数与三维相比较小,而法向动量协调系数差别不大.     

温度和分子间作用对切向动量协调系数的耦合效应

切向动量协调系数(TMAC)是描述气体滑移流动的重要边界条件。运用非平衡分子动力学方法,构建了能够反映流体粒子与壁面粒子相互作用关系的物理模型。结果显示:当壁面存在吸附层或壁面无气体吸附层时,同一势能强度下TMAC值都随着温度的升高而降低;而在吸附层能够解吸附的温度,TMAC值发生突跃。在本文的模拟条件下,气体粒子离开壁面吸附的能力和壁面粒子及吸附层粒子热运动产生的粗糙度决定了TMAC值的分布。     

气-固界面作用物理模型的确定性实现方法

针对最新发展的气-固界面作用物理模型,提出一种确定性计算方法。算法的核心是：针对离散速度空间分布函数,采用确定性算法计算出单轮气-固碰撞作用的散射核函数矩阵,并通过考虑吸引势阱作用和迭代累加多次气-固碰撞作用,物理地反映气体分子入射/反射速度分布函数在气-固界面上的变化。与现有基于Monte Carlo随机采样的实现方法对比,该确定性实现方法可以快速且准确地计算边界上气体分子入射/反射速度分布函数的映射关系。因此,该方法可为以离散速度空间为基础的稀薄气体流动模拟方法提供先进的边界条件。     

IR reflection-absorption spectroscopy of CO adsorbed on palladium

Kinetic theory is used to compute the flux and relative translational kinetic energy incident upon a surface oscillating in a rarefied gas. The flux incident upon the oscillatory surface is deficient in low-velocity molecules from the gas during the reentry half of a vibration cycle, in which the surface moves into the gas, because all molecules which are to strike the oscillatory surface from the gas must cross the plane of maximum surface extension during that cycle. The deficiency is largely compensated for by recapture, during the reentry half of a cycle, of low-velocity molecules emitted during the recession half of the cycle. The result of these two opposing effects is that the average energy of gas-surface collisions, and therefore the temperature rise of an oscillatory surface is greater than that of a constant speed plate of the same rms velocity in the same gas. For argon at 300 K and 10^?3 torr incident upon a surface with an rms velocity of 3.3 × 10³ cm sec^?1 the apparent average temperature of incident molecules is 302.91 ° for an oscillatory surface and 302.63 ° for a constant speed plate. Measurements of the temperature rise of an oscillatory surface offer a way to measure thermal accomodation coefficients.     

Gas molecule-molecule interaction and the gas-surface scattering effect on the rarefied gas flow through a slit into a vacuum

The effect of the gas molecule-molecule interaction and the gas-surface scattering on the gas flow through a slit into a vacuum are investigated in a wide range of the gas rarefaction using the direct simulation Monte Carlo method. To study the gas molecule-molecule interaction influence, we used the variable hard sphere and variable soft sphere models defined for an inverse-power-law potential and the generalized hard sphere model defined for the 12–6 Lennard-Jones potential. The Maxwell, Cercignani-Lampis, and Epstein models were used to simulate the gas-surface scattering. This study demonstrates that the gas molecule-molecule interaction can have a significant influence on the rarefied gas flow through a slit, while the influence of the gas-surface scattering is negligibly small. The presented numerical results are in agreement with the corresponding experimental ones. The article is published in the original.     

Asymptotic Theory for the Time-Dependent Behavior of a Slightly Rarefied Gas Over a Smooth Solid Boundary

A time-evolution of a slightly rarefied monoatomic gas, namely a gas for small Knudsen numbers, which is perturbed slowly and slightly from a reference uniform equilibrium state at rest is investigated on the basis of the linearized Boltzmann equation. By a systematic asymptotic analysis, a set of fluid-dynamic-type equations and its boundary conditions that describe the gas behavior up to the second order of the Knudsen number are derived. The developed theory covers a general intermolecular potential and a gas-surface interaction. It is shown that (i) the compressibility of the gas manifests itself from the leading order in the energy equation and from the first order in the continuity equation; (ii) although the momentum equation is the Stokes equation, it contains a double Laplacian of the leading order flow velocity as a source term at the second order; (iii) a double Laplacian source term also appears in the energy equation at the second order; (iv) the slip and jump conditions are the same as those in the time-independent case up to the first order, and the difference occurs at the second order in the jump conditions as the terms of the divergence of flow velocity and of the Laplacian of temperature. Numerical values of all the slip and jump coefficients are obtained for a hard-sphere gas by the use of a symmetric relation developed recently.     

The effect of surface roughness on rarefied gas flows by lattice Boltzmann method

This paper studies the roughness effect combining with effects of rarefaction and compressibility by a lattice Boltzmann model for rarefied gas flows at high Knudsen numbers. By discussing the effect of the tangential momentum accommodation coefficient on the rough boundary condition, the lattice Boltzmann simulations of nitrogen and helium flows are performed in a two-dimensional microchannel with rough boundaries. The surface roughness effects in the microchannel on the velocity field, the mass flow rate and the friction coefficient are studied and analysed. Numerical results for the two gases in micro scale show different characteristics from macroscopic flows and demonstrate the feasibility of the lattice Boltzmann model in rarefied gas dynamics.     

Molecular dynamics simulations of translational thermal accommodation coefficients for time-resolved LII

Time-resolved laser-induced incandescence demands precise knowledge of the thermal accommodation coefficient, but little is known about the gas-surface scattering physics that underlies this parameter. This paper presents a molecular dynamics simulation that shows how the thermal accommodation coefficient is influenced by the gas molecular mass and gas temperature. The molecular dynamics results also define scattering kernels that can be used as boundary conditions in Direct Simulation Monte Carlo simulations of heat and momentum transfer between soot aggregates and surrounding gas molecules.     

Discrete Boltzmann Method with Maxwell-Type Boundary Condition for Slip Flow

The rarefied effect of gas flow in microchannel is significant and cannot be well described by traditional hydrodynamic models. It has been known that discrete Boltzmann model(DBM) has the potential to investigate flows in a relatively wider range of Knudsen number because of its intrinsic kinetic nature inherited from Boltzmann equation.It is crucial to have a proper kinetic boundary condition for DBM to capture the velocity slip and the flow characteristics in the Knudsen layer. In this paper, we present a DBM combined with Maxwell-type boundary condition model for slip flow. The tangential momentum accommodation coefficient is introduced to implement a gas-surface interaction model.Both the velocity slip and the Knudsen layer under various Knudsen numbers and accommodation coefficients can be well described. Two kinds of slip flows, including Couette flow and Poiseuille flow, are simulated to verify the model.To dynamically compare results from different models, the relation between the definition of Knudsen number in hard sphere model and that in BGK model is clarified.