粗糙微通道内气体流动的分子动力学研究

为研究粗糙度对微通道内气体流动及其边界滑移性质的影响,采用分子动力学模拟方法研究了氩气在0.1 μm铂通道内的流动,通道表面粗糙度由三角粗糙元阵列构成。气体流动的边界条件决定于2个准则数: A/λ和Kn=λ/L(其中A为壁面粗糙度、λ为气体分子的平均自由程、L为流动系统的尺度)。Maxwell基于Kn的滑移模型当A/λ<0.25时适用;A/λ<1时,气体流动存在边界速度滑移,体现出稀薄效应,A/λ≈1时为无滑移,A/λ>1时为等效负滑移.     

微尺度流动的扩展二阶滑移边界条件

本文提出了求解微流动Navier-Stokes方程的扩展二阶滑移边界条件.该边界条件依赖一个有效的分子平均自由程,并考虑了气体分子与壁面碰撞的影响,能够刻画Knudsen层内的流动.基于该边界条件,对微通道内的气体流动进行了研究,计算结果与实验数据和其他理论结果相吻合.     

表面粗糙度的分形特征及其对微通道内层流流动的影响

基于分形几何学,研究了表面粗糙度的分形特征.采用Weierstrass- Mandelbrot函数对多尺度自仿射的表面粗糙度进行了描述;建立了微通道内层流流动的三维模型并对表面粗糙度的影响进行了数值模拟,分析了雷诺数、相对粗糙度和分形维数对流动阻力特性的影响.研究结果表明,与常规尺度通道不同,粗糙微通道的Poiseuille数不再是常数,而是随雷诺数近似线性增加;相对粗糙度越大,流动产生的回流和分离所导致的流动压降越明显.在相同的相对粗糙度下,粗糙表面的分形维数越大,表面轮廓变化就越频繁,这也将导致流动阻 关键词： 粗糙度 层流阻力系数 微通道 分形     

微通道内气体流动的格子Boltzmann方法研究

采用格子Boltzmann方法模拟了微通道在滑移区内不同Knudsen数下的微气体Poiseuille流,分析了微气体流动的速度分布以及流量与压降的关系,并给出了相对滑移长度和Poiseuille数随Knudsen数的变化特性。研究结果表明,微气体Poiseuille流的速度轮廓呈抛物线分布,但是边界速度大于0,出现速...     

粗糙表面形貌对微尺度流动换热的影响

基于分形几何定量描述了多尺度自仿射的粗糙表面形貌,建立了微通道内层流流动换热的理论模型并对表面形貌的影响进行了数值模拟.研究表明,自仿射分形维数直接反映了表面轮廓的不规则度,对于两个具有相同统计粗糙度的轮廓,可能存在不同的分形维数;与常规尺度通道不同,雷诺数、粗糙高度和粗糙表面分形维数都对微通道内层流流动换热有着重要影...     

微通道中气体流动的格子Boltzmann数值模拟

采用可压缩格子Boltzmann模型及非平衡外推边界条件,数值模拟微通道中的气体在滑移区域(Kn≤0.1)内的流动,计算结果包括出口速度剖面、通道中心压力分布以及质量流率等,与理论结果及其他实验结果符合得很好.还模拟了180°弯曲通道中的气体流动.结果表明,滑移速度的存在抑制了边界层的分离,因此在弯曲处不存在漩涡.计算结果还表明,弯道的存在显著影响了气体的质量流率.     

粗糙纳通道内流体流动与传热的分子动力学模拟研究

为研究粗糙表面对纳尺度流体流动和传热及其流固界面速度滑移与温度阶跃的影响,本文建立了粗糙纳通道内流体流动和传热耦合过程的分子动力学模型,模拟研究了粗糙通道内流体的微观结构、速度和温度分布、速度滑移和温度阶跃并与光滑通道进行了比较,并分析了固液相互作用强度和壁面刚度对界面处速度滑移和温度阶跃的影响规律. 研究结果表明,在外力作用下,纳通道主流区域的速度分布呈抛物线分布,由于流体流动导致的黏性耗散使得纳通道内的温度分布呈四次方分布. 并且,在固体壁面处存在速度滑移与温度阶跃. 表面粗糙度的存在使得流体剪切流动产生了额外的黏性耗散,使得粗糙纳通道内的流体速度水平小于光滑通道,温度水平高于光滑通道,并且粗糙表面的速度滑移与温度阶跃均小于光滑通道. 另外,固液相互作用强度的增大和壁面刚度的减小均可导致界面处速度滑移和温度阶跃程度降低. 关键词： 速度滑移 温度阶跃 流固界面 粗糙度     

微流动的同伦分析法解析求解

建立了滑移流区二阶滑移流动的控制方程,引入相似变换将微尺度二阶滑移流动的控制方程变换为非线性常微分方程.采用同伦分析法解析求解了该方程,获得了微尺度滑移流动的速度分布,发现在不同Knudsen数下,充分发展段无量纲速度分布有同一交点.通过与文献结果对比验证了本文求解结果的正确性.为深入研究微尺度滑移流动的物理本质提供了数学支持.     

粗糙度与气体稀薄性对微尺度流动特性的影响

对氮气和氦气在粗糙微通道以及光滑微通道内流动进行了阻力特性实验研究。实验结果表明,即使在较小的相对粗糙度高度下,由于微通道中的粗糙度分布密集,会极大地增加流动阻力,这是导致文献中微通道流动阻力系数实验值相互偏差的主要原因之一;而对于滑移区的气体流动,气体稀薄性使流动阻力明显减小而导致流量增加。     

压缩性和稀薄性对微通道气体流动换热的影响

本文对于微通道内稀薄气体二维可压缩滑移流动建立了数学模型,采用连续介质流动控制方程与壁面速度滑移和温度跳跃边界条件相组合描述该问题,并利用SIMPLE算法求解,所得结果与文献进行了对比,在相同条件下得到了较高的一致性.文中利用该计算模型的计算结果分析了气体的压缩性和稀薄性对微通道气体流动的影响,结果表明在所计算工况下稀薄性的影响更大。     

Non-Maxwell slippage induced by surface roughness for microscale gas flow: a molecular dynamics simulation

Rarefied gas flows in rough microchannels are investigated by non-equilibrium molecular dynamics simulations. The surface roughness is modelled by an array of triangular modules. The Maxwell slip model is found to break down due to the surface roughness for gas flows in microchannels with large surface roughness. Non-Maxwell slippage shows that the slip length is smaller than that predicted by the Maxwell model and is nonlinearly related to the mean free path. For larger surface roughness and smaller Knudsen number, the non-Maxwell effect becomes more pronounced. The boundary conditions, generally including velocity slip, no-slip and negative slip, depend not only on the Knudsen number but also on the surface roughness. Simulation results show that A/λ?≈?1 is a good criterion to validate the no-slip boundary condition and A/λ?>?0.3 can be a criterion to judge the occurrence of non-Maxwell slippage, where A is the surface roughness size and?λ?is the mean free path of gas molecules. The permeability enhanced by the surface roughness may be responsible for the roughness-induced non-Maxwell slippage.     

粗糙微通道内液体流动的多尺度耦合模拟

采用分子动力学与有限容积法多尺度耦合算法对粗糙微通道内的液体Poiseuille流动进行了模拟。分析了粗糙元高度、分布以及几何形状对通道内流动速度和边界滑移长度的影响。结果表明：随着粗糙元高度的增加,流动速度和粗糙元间隙底部壁面上滑移长度均减小;粗糙元分布越密或不同几何形状粗糙元所对应的固壁原子数越多,滑移长度越小,但...     

随机散射屏的原子力显微镜形貌分析及其光散射特性

用原子力显微镜对三种不同粗糙度的随机散射屏的表面形貌进行了测量分析,发现它们在短程范围内具有明显的分形特征.对于粗糙度较大和较小的散射屏,分形特征分别以无规则的 高度调制和无特征大小的小颗粒的形式存在.用自仿射分形表面模型对散射屏的统计特性进 行了描述和拟合.光散射测量发现,散射光强在远轴区域按负幂函数下降,理论分析证明这 源于表面的分形结构;在近轴区域有散射亮环存在,用自仿射分形表面模型尚不能给出理论解释. 关键词： 随机散射屏 光散射 自仿射分形表面     

Theoretical and numerical studies of non-equilibrium slip effects on a catalytic surface

A new concentration slip model to describe the rarefied gas effect on the species transport in microscale chemical reactors was derived from the approximate solution of the Boltzmann equation. The present model is more general and recovers the existing models in the limiting cases. The analytical results showed that the concentration slip is dominated by two different mechanisms, the reaction induced concentration slip (RIC) and the temperature slip induced concentration slip (TIC). The magnitude of RIC slip is proportional to the product of the Damköhler number and Knudsen number. The impact of the velocity, concentration and temperature slips on the coupling between the surface catalytic reactions and the homogeneous gas phase reactions was examined using the detailed chemistry of hydrogen and methane within a wide range of accommodation coefficients in a two-dimensional microscale chemical reactor. The results showed that the impact of reaction induced concentration slip (RIC) effects on catalytic reactions strongly depends on the Damköhler number, the Knudsen number and the surface accommodation coefficient. It was found that the TIC slip had a strong effect on the fuel oxidation rates and the RIC slip dramatically changed the mass fraction distributions of radicals, especially when the mass accommodation coefficients were far less than unity.     

纳米通道粗糙内壁对流体流动行为的影响

纳米流动系统具有高效、经济等优势,在众多领域具有广泛的应用前景.因该类系统具有极高的表面积体积比,致使界面滑移效应对流动具有显著影响.本文采用分子动力学方法以两无限大平行非对称壁面组成的Poiseuille流动为对象,分析了壁面粗糙度与润湿性变化对通道内流体流动的影响.对于不同结构类型的壁面,需要通过水动力位置来确定固液界面位置,准确计算固液界面位置有助于更好地分析界面滑移效应.研究结果表明,上下壁面不对称会引起通道内流场参数分布的不对称,壁面粗糙度及润湿性的变化会影响近壁面附近流体原子的流动特性,由于壁面凹槽的存在,粗糙壁面附近的数密度分布低于光滑壁面一侧.壁面粗糙度及润湿性的变化会影响固液界面位置,肋高变化及壁面润湿性对通道中速度分布影响较大,界面滑移速度及滑移长度随肋高和润湿性的增大而减小;肋间距变化对通道内流体流动影响较小,界面滑移速度和滑移长度基本保持恒定.     

Surface correlations of hydrodynamic drag for transitionally rough engineering surfaces

Rough surfaces are usually characterised by a single equivalent sand-grain roughness height scale that typically needs to be determined from laboratory experiments. Recently, this method has been complemented by a direct numerical simulation approach, whereby representative surfaces can be scanned and the roughness effects computed over a range of Reynolds number. This development raises the prospect over the coming years of having enough data for different types of rough surfaces to be able to relate surface characteristics to roughness effects, such as the roughness function that quantifies the downward displacement of the logarithmic law of the wall. In the present contribution, we use simulation data for 17 irregular surfaces at the same friction Reynolds number, for which they are in the transitionally rough regime. All surfaces are scaled to the same physical roughness height. Mean streamwise velocity profiles show a wide range of roughness function values, while the velocity defect profiles show a good collapse. Profile peaks of the turbulent kinetic energy also vary depending on the surface. We then consider which surface properties are important and how new properties can be incorporated into an empirical model, the accuracy of which can then be tested. Optimised models with several roughness parameters are systematically developed for the roughness function and profile peak turbulent kinetic energy. In determining the roughness function, besides the known parameters of solidity (or frontal area ratio) and skewness, it is shown that the streamwise correlation length and the root-mean-square roughness height are also significant. The peak turbulent kinetic energy is determined by the skewness and root-mean-square roughness height, along with the mean forward-facing surface angle and spanwise effective slope. The results suggest feasibility of relating rough-wall flow properties (throughout the range from hydrodynamically smooth to fully rough) to surface parameters.     

由于表面粗糙引起的激光声表面波色散的实验和理论研究

表面粗糙是材料制造过程中必有的副产物, 粗糙表面会引起其中传播的声表面波的速度发生变化. 在利用激光声表面波对材料性质进行评估时, 常用宽带的激光声表面波速度频散特性对材料性质进行反演. 为了研究表面粗糙度是否能作为反演的特征参数之一, 本文建立了激光在表面粗糙样品中激发声表面波、聚偏氟乙烯换能器宽带接收声表面波的实验装置来研究不同粗糙度表面对声表面波速度的影响; 理论上建立了激光在粗糙表面中激发声表面波的计算模型, 利用有限元法得到声表面波的时域特征, 并进一步得到声表面波的速度色散曲线, 理论结果和实验结果能很好地拟合. 这为利用激光声表面波对表面粗糙的评估提供理论和实验依据. 关键词： 表面粗糙 激光声表面波 速度色散 聚偏氟乙烯传感器 有限元法     

Monte Carlo methods for turbulent tracers with long range and fractal random velocity fields

Monte Carlo methods for computing various statistical aspects of turbulent diffusion with long range correlated and even fractal random velocity fields are described here. A simple explicit exactly solvable model with complex regimes of scaling behavior including trapping, subdiffusion, and superdiffusion is utilized to compare and contrast the capabilities of conventional Monte Carlo procedures such as the Fourier method and the moving average method; explicit numerical examples are presented which demonstrate the poor convergence of these conventional methods in various regimes with long range velocity correlations. A new method for computing fractal random fields involving wavelets and random plane waves developed recently by two of the authors [J. Comput. Phys. 117, 146 (1995)] is applied to compute pair dispersion over many decades for systematic families of anisotropic fractal velocity fields with the Kolmogorov spectrum. The important associated preconstant for pair dispersion in the Richardson law in these anisotropic settings is compared with the one obtained over many decades recently by two of the authors [Phys. Fluids 8, 1052 (1996)] for an isotropic fractal field with the Kolmogorov spectrum. (c) 1997 American Institute of Physics.     

半干旱草原下垫面动量和感热总体输送系数参数化关系研究

本文利用锡林郭勒草原2008年春季近地层涡旋相关系统和铁塔的风、 温平均梯度观测资料, 分析了总体输送系数随梯度Richardson数的变化特征, 建立了动量总体输送系数随大气稳定度、近地层风速以及感热总体输送系数随大气稳定度和近地层气温的关系. 中性条件下, 半干旱草原植被下垫面动量总体输送系数与近地层大气动力状态之间存在明显的相互作用, 总体输送系数与近地面层风速之间满足二次曲线拟合关系; 风速较小时, 大气动力特征对地表粗糙度长度的改变不是很明显, 动量总体输送系数随气流增强而增大; 而当风速较大时, 强风速会使植被高度发生改变, 动量总体输送系数随气流增强而减小. 另外, 感热总体输送系数与近地层气温之间也存在二次曲线关系. 动量总体输送系数与近地层风速之间的关系、感热总体输送系数与近地层气温之间关系的建立为总体输送系数参数化提供了重要途径, 同时该方案避免了对动力学粗糙度长度和热力学粗糙度长度的求解. 关键词： 总体输送系数 参数化 湍流通量 相似性函数