Self-consistent kinetic modeling of electron distribution function in a low-pressure inductively coupled plasma

Based upon the spatially inhomogeneous Boltzmann equation in two-term approximation coupled with electromagnetic and fluid model analysis for the recently developed inductively coupled plasma sources, a self-consistent electron kinetic model is developed. The electron distribution function, spatial distributions of the electron density and ionization rate are calculated and discussed.

The kinetic-theory for dilute solid liquid 2-phase flow

On the basis of a brief review of the continuum theory for macroscopic descriptions and the kinetic theory for microscopic descriptions in solid/liquid two-phase flows, some suggestions are presented, i.e. the solid phase may be described by the Boltzmann equation and the liquid phase still be described by conservation laws in the continuum theory. Among them the action force on the particles by the liquid fluid is a coupling factor which connects the phases. For dilute steady solid/liquid two-phase flows, the particle velocity distribution function can be derived by analogy with the procedures in the kinetic theory of gas molecules for the equilibrium state instead of being assumed, as previous investigators did. This done, more detailed information, such as the velocity probability density distribution, mean velocity distribution and fluctuating intensity etc. can be obtained directly from the particle velocity distribution function or from its integration. Experiments have been performed for dilute solid/liquid two-phase flow in a 4 x 6 cm2 sized circulating square pipe system by means of laser Doppler anemometry so that the theories can be examined. The comparisons show that the theories agree very well with all the measured data.

Phase relaxation model for 2-phase flows

A phase relaxation model (PRM) for 2-phase flows is presented in this paper on the basis of three principal assumptions. The basic equations for PRM arc derived from the Boltzmann equations for gas-partlcle mixture, The general characteristics and solving process of the PRM's basic equations are also presented and discussed. Many terms in the PRM's basic equations contain a factor ε= ρgρp/ρg+ρp2 which is an intrinsic small parameter for 2-phase mixture, with ρg and ρp being respectively the densities of gas and particle phases.This makes it possible to simplify the computation of the PRM's basic equations. The model is applied to for example, studying file steady propagation of shock waves in gas-particle mixture. The analysis shows that with an increase of shock wave strength the relaxation process behind a gasdynamics shock front becomes a kind of dynamics relaxation instead of the standard exponential relaxation process. A method of determining experimentally the velocity and tem...更多perature relaxation rates (or times) of gas-particle flows is suggested and analyzed.

The concentration-jump coefficient in a rarefied binary gas-mixture

The problem of the concentration jump of a vapour in the vicinity of a plane wall, which consists of the condensed phase of the vapour, in a rarefied gas mixture of that vapour (A) and another 'inert' gas (B), is considered. The general formulation of the problem of determining the concentration-jump coefficient for dA is given. In the Knudsen layer the simplest model of Boley-Yip theory is used to simplify the Boltzmann equations for the binary gas mixture. The numerical calculation of the concentration jump coefficient for dA for various values of evaporation coefficient of A is illustrated for the case of the equilibrium concentration of B being much greater than that of A, for which experimental data are available.

认识稀薄气体动力学

以通俗易懂的方式介绍了空气动力学当气体间断分子效应显著时发展起来的特殊分支－－稀薄气体动力学。讨论了非平衡现象与稀薄气体动力学的。通过与8速度气体模型的间断Boltzmann方程的对比，解释了Boltzmann方程碰撞项的物理意义和数学困难，简要综述了其一般解法。讨论了分在物体表面的反射和问题的边界条件，着重介绍了直接模拟Monte Carlo (DSMC)方法和为克服低速稀薄流动（如MEMS中流动）中模拟困难的信息保存（IP）方法。

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Monte Carlo Modeling of Micro Scale Gas Flows

An information preservation (IP) method has been used to simulate many micro scale gas flows. It may efficiently reduce the statistical scatter inherent in conventional particle approaches such as the direct simulation Monte Carlo (DSMC) method. This paper reviews applications of IP to some benchmark problems. Comparison of the IP results with those given by experiment, DSMC, and the linearized Boltzmann equation, as well as the Navier-Stokes equations with a slip boundary condition, and the lattice Boltzmann equation, shows that the IP method is applicable to micro scale gas flows over the entire flow regime from continuum to free molecular.

Single curved fiber sedimentation under gravity

Dynamics of single curved fiber sedimentation under gravity are simulated by using the lattice Boltzmann method. The results of migration and rotation of the curved fiber at different Reynolds numbers are reported. The results show that the rotation and migration processes are sensitive to the curvature of the fiber. (c) 2007 Elsevier Ltd. All rights reserved.

模拟MEMS中气体流动的工具——信息保存法

MEMS中气体流动因特征尺度小而是稀薄气体的领域,本文首先介绍处理低速稀薄气体流动的一些方法:线化Boltzmann方程,Lattice Boltzmann方法,加滑移边界的Navier-Stoks方程,以及DSMC方法,并讨论它们模拟MEMS中过渡领域低速流动所遇到的困难。信息保存法克服了流速低使得信息噪声比小引起统计模拟的困难,已成功模拟了一些一维和二维问题。MEMS中流速低和大的长宽比的特点还引起出入口边界条件相互影响需要协调的问题,通过微槽道流动的算例,在模拟中采用守恒形式的质量守恒方程和超松弛法成功地解决了这一问题。处理有温度变化的MEMS流动问题和跨越领域的混合算法是重要的问题,本文用信息保存法也进行了有益的尝试。

钱学森技术科学思想与力学

本书收录会议论文60篇。

　　本书内容大体分为三类。一类着重介绍了钱学森先生的风范，他的爱国主义，他的勤奋好学，他的治学精神和治学态度，他在学术上的重大贡献以及他的科学思想。第二类是介绍在他的领导或影响下，五十年来我国力学研究取得的成就及进一步的进展。第三类涉及我国力学界正在开展的其他方面的学术研究进展。我们的共同目的是继承优良传统，开拓我国力学在新世纪的发展。

　　读者对象：力学、航空航天、机械等领域的工程技术人员，研究人员。

稀薄气体动力学

src="http://img14.360buyimg.com/n4/15638/2572fc6a-0288-44f9-97c1-643e21743417.jpg" border="0" alt="" hspace="8" width="100" height="100" align="left" />通常的气体动力学方法，当气体分子的平均自由程与流场特征长度相比不可忽略时，不再适用，要采用稀薄气体动力学的方法。这适用于航天飞行器在高空飞行时受的力和热，也适用于微机电系统和真空系统等离子体材料加工等21世纪技术前沿领域。本书系统、简明地阐述稀薄气体动力学方法，给出方法的基础并着重介绍直接模拟Monte Carlo(DSMC)方法以及与低速稀薄气体流动相关的前沿课题。全书共分7章。前两章是作为学科的基础引入的，第1章以空气为对象对于分子能态结构、能态分布以极小篇幅作了简要概括的叙述，以作为了解稀薄气流非平衡现象物理基础的初步。第2章对包括双体碰撞、Boltzmann方程以及气体的平衡态等分子动理论的基础做了必要的讨论，其中包括了对唯像论分子相互作用模型、变径硬球（VHS）、变径软球（VSS）和概括化硬球（GHS）等模型的介绍。第3章讨论了各种分子和表面的相互作用模型，包括反映细致平衡的互易原理和基于此原理的CLL模型的阐述。第4章讨论自由分子流。第5章讨论应用于滑流领域的各连续介质方程及滑流边界条件，一些简单解以及热泳问题。第6章则较全面、概括地介绍了求解过程领域中的各种解析和数值方法。第7章介绍了直接模拟Monte Carlo(DSMC)方法，讨论了非平衡流动及低速稀薄流动等前沿课题，包括处理内能松弛、化学反应的方法、用于复杂流场通用软件的方法、低速稀薄流动的信息保存（IP）方法等。　　本书适合高等学校力学一航空航天专业高年级学生、研究生及从事气动力学和航天研究的科研人员参考阅读。

通常的气体动力学方法，当气体分子的平均自由程与流场特征长度相比不可忽略时，不再适用，要采用稀薄气体动力学的方法。这适用于航天飞行器在高空飞行时受的力和热，也适用于微机电系统和真空系统等离子体材料加工等21世纪技术前沿领域。本书系统、简明地阐述稀薄气体动力学方法，给出方法的基础并着重介绍直接模拟Monte Carlo(DSMC)方法以及与低速稀薄气体流动相关的前沿课题。

符号表
绪论
第1节 稀薄气体动力学的提出
第2节 气体的分子模型
第3节 分子平均自由程
第4节 流动的领域划分
第5节 非平衡现象与稀薄气体动力学
第6节 相似准则
第1章 分子结构与能态
第1节 双原子分子
第2节 分子的能态分布
第3节 分子的内能、内自由度和内能分布函数
第2章 分子动理论基础
第1节 速度分布函数
第2节 宏观量的表达
第3节 分子的双体碰撞模型
第4节 碰撞截面与分子模型
第5节 Boltzmann方程
第6节 碰撞积分与气体分子的总碰撞数
第7节 碰撞积分的计算
第8节 Maxwell输运方程——矩方程
第9节 Maxwell分布
第10节 气体的平衡态
第11节 8速度气体模型
第12节 混合气体
第3章 分子表面相互作用
第1节 引言
第2节 镜面反射与漫反射，适应系数
第3节 互易性原理
第4节 CLL分子表面相互作用模型
第4章 自由分子流
第1节 气体中的分子数目通量和动量通量
第2节 作用于物体的气动力
第3节 表面元素的热传导
第4节 自由分子流出与热流逸
第5节 Couette流动与平板间的传热问题
第6节 无碰撞Boltzmann方程的通解，非定常流动
第5章 连续介质模型
第1节 引言
第2节 基本方程
第3节 滑流边界条件
第4节 一些简单问题的求解
第5节 热蠕动与热泳
第6章 过渡领域
第1节 概述
第2节 线化的BoltzmanN方程
第3节 矩方法
第4节 模型方程
第5节 有限差分法
第6节 间断纵坐标方法
第7节 积分方法
第8节 直接模拟方法
第7章 直接模拟Monte方法
第1节 引言
第2节 碰撞的取样
第3节 DSMC方法求解问题实例
第4节 内能的激发与松弛
第5节 化学反应的模拟
第6节 复杂流场的计算，位置元方法
第7节 微尺度低速气体流动，信息保存法
附录I 气体的性质和分子性质
附录II 分布函数求矩遇到的积分
附录III 具有给定分布的随机数的取样
附录IV Couette问题程序
参考文献
主题词索引

Size- And Temperature-Dependent Thermal Expansion Coefficient Of A Nanofilm

The thermal expansion coefficient (TEC) of an ideal crystal is derived by using a method of Boltzmann statistics. The Morse potential energy function is adopted to show the dependence of the TEC on the temperature. By taking the effects of the surface relaxation and the surface energy into consideration, the dimensionless TEC of a nanofilm is derived. It is shown that with decreasing thickness, the TEC can increase or decrease, depending on the surface relaxation of the nanofilm.

Solute transport in nanochannels with roughness-like structures

Newfound attention has been given to solute transport in nanochannels. Because the electric double layer (EDL) thickness is comparable to characteristic channel dimensions, nanochannels have been used to separate ionic species with a constant charge-to-size ratio (i.e., electrophoretic mobility) that otherwise cannot be separated in electroosmotic or pressure- driven flow along microchannels. In nanochannels, the electrical fields within the EDL cause transverse ion distributions and thus yield charge-dependent mean ion speeds in the flow. Surface roughness is usually inevitable during microfabrication of microchannels or nanochannels. Surface roughness is usually inevitable during the fabrication of nanochannels. In the present study, we develop a numerical model to investigate the transport of charged solutes in nanochannels with hundreds of roughness-like structures. The model is based on continuum theory that couples Navier-Stokes equations for flows, Poisson-Boltzmann equation for electrical fields, and Nernst-Planck equation for solute transports. Different operating conditions are considered and the solute transport patterns in rough channels are compared with those in smooth channels. Results indicate that solutes move slower in rough nanochannels than in smooth ones for both pressure- driven and electroosmotic flows. Moreover, solute separation can be significantly improved by surface roughness under certain circumstances.

Near dipole-dipole effects on the propagation of few-cycle pulse in a dense two-level medium

The propagation behaviors, which include the carrier-envelope phase, the area evolution and the solitary pulse number of few-cycle pulses in a dense two-level medium, are investigated based on full-wave Maxwell-Bloch equations by taking Lorentz local field correction (LFC) into account. Several novel features are found: the difference of the carrier-envelope phase between the cases with and without LFC can go up to pi at some location; although the area of ultrashort solitary pulses is lager than 2 pi, the area of the effective Rabi frequency, which equals to that the Rabi frequency pluses the product of the strength of the near dipole-dipole (NDD) interaction and the polarization, is consistent with the standard area theorem and keeps 2 pi; the large area pulse penetrating into the medium produces several solitary pulses as usual, but the number of solitary pulses changes at certain condition. (C) 2005 Optical Society of America.