基于微通道构型的微流体流动控制研究

微通道不仅仅是作为流体流动的单元, 更是进行流体控制的工具,微通道自身特性和特征用在实现微流体的驱动、进样、混合、分离以及液滴的产生、控制等方面已经表现出了良好的效果.由于微通道中比表面积非常大, 表面效应极大影响流体流动,近年来多数研究集中在应用表面效应来实现微流体驱动与控制,而以利用微通道结构特征实现流体流动控制为目标的研究成果相对较少.为了提高对通道构型作用的认识,主要介绍了基于微通道构型的无可动部件的流体微阀和基于微通道构型微小液滴的产生及流动控制器两个方面的发展情况,表明微通道构型在微流体控制中同样可以发挥重要作用,甚至有望带来微流控技术的突破.      

微纳米结构壁面对流场及动压润滑的影响研究

固体边界具有的微纳米结构将影响流体在近壁面处的流动行为,进而由于尺度效应改变流体在整个微间隙的流动或润滑规律.将壁面可渗透微纳米结构等效为多孔介质薄膜,采用Brinkman方程来描述流体在近壁面边界渗透层内的流动,并将其与自由流动区域的不可压缩流体Navier-Stokes控制方程耦合,在界面处的连续边界条件下求解和分析了速度分布规律和压力变化规律.针对恒定法向承载力的油膜润滑条件,进一步讨论了静止表面或运动表面的微纳米结构对近壁面流动行为的影响;并揭示了考虑壁面微纳米结构的流体动压润滑的油膜厚度和摩擦系数的变化规律.论文结果为具有可渗透微结构表面的微间隙流动与润滑提供了理论参考.     

微纳米尺度效应作用下充流微通道系统波动特性研究

结合非局部弹性应力/应变梯度耦合本构关系和流体非局部应力关系式,基于Euler梁理论,建立了充流微通道流固耦合波传导模型;根据耦合固体非局部应力/应变梯度弹性效应以及流体非局部效应,分别模拟了微通道和管腔内流体的尺度效应,推导得出了充流微通道在微纳米尺度的波动控制方程和边界条件。通过对控制方程的求解,分析了不同类型尺度效应对微通道的波动和振动特性的影响。结果显示,各类尺度效应对系统的动力学特性影响不同。微通道非局部弹性效应对波动产生阻尼,特别是对波长较短的波传导;而应变梯度弹性效应对波传导有促进作用,且该效应对波动的影响与波长无关;非局部效应和应变梯度效应对微通道刚度产生不同影响,非局部效应降低刚度,应变梯度效应增加刚度。     

流动聚焦实验

介绍了一种制备微纳米量级颗粒的流动聚焦技术,它的最基本的特点是从毛细管流出的液体由高速运动气体驱动经小孔聚焦形成稳定的锥,锥顶端射出的微射流因不稳定性破碎成小颗粒.实验在自行设计的装置上完成.分析了流动聚焦技术中影响锥-射流以及颗粒形貌的因素,总结了过程中装置的结构参数以及气体压力降、液体流量和物性等参数对射流的直径的影响.结果显示该技术制备的颗粒单分散性好,直径达到微米和亚微米量级,在工业方面有重大应用前景.     

纳米流体液滴撞击壁面铺展动力学特性研究

纳米流体液滴撞击固体壁面的铺展动力学特性是基于液滴沉积实现高效传热传质过程的关键因素,然而由于纳米流体的非牛顿流变特性及液滴内微流动与纳米颗粒的耦合作用,目前对纳米流体液滴撞击固体壁面的铺展动力学行为缺乏足够的认识.本研究利用了两步法分别配制了分散有3种纳米颗粒的均匀稳定纳米流体(碳纳米管、石墨烯、纳米石墨粉),并对流体的流变特性进行了测量分析.利用显微高速数码摄像技术捕捉了液滴撞击固体壁面的动态过程,通过图像处理技术分析铺展过程中液滴的无量纲高度、铺展因子及动态接触角,探究了液滴在韦伯数约为200及800时撞击壁面后铺展沉积形态的演变规律.研究表明,3种不同纳米颗粒的加入均使基液表现出明显的剪切变稀特性,在液滴撞击壁面的铺展过程中,流体的剪切黏度起重要作用,液滴的无量纲高度和铺展因子的变化幅度随着纳米流体剪切黏度的增大而减小.纳米流体液滴撞击疏水表面时能更快的达到平衡状态,液滴的惯性力主导着液滴的初始铺展阶段,液滴的铺展范围和速度随撞击速度的增大而增大.开展该研究能够为基于液滴沉积的增益冷却技术以及微型高导热及导电材料的制造提供理论依据和技术指导.      

边界膜剪切弹性模量对于微接触性能的影响

建立接触模型,理论分析了微接触中边界膜剪切弹性模量对于接触性能的影响. 接触区由两平行平面形成,属一维接触. 上接触表面为粗糙表面,具有矩形微凸体. 下接触表面为光滑平面. 两接触表面均处理成刚性表面. 微接触区中充满流体. 它分成两个子区,在微接触的出口区由于极小的接触间隙充满边界膜,在微接触的入口区由于接触间隙较大充满流体膜. 边界膜和流体膜行为决定整个微接触性能. 当膜厚较大时,这里边界膜可看成纳米级薄膜. 由于上接触表面处有限的剪应力承受能力,边界膜可于上接触表面滑移. 设下接触表面处剪应力承受能力很大而边界膜在下接触表面不滑移. 由于边界膜-接触表面间相互作用,边界膜黏度、密度和剪切弹性模量均沿膜厚变化,在理论分析中使用它们的等效值,这些值与边界膜厚度有关. 流体膜在两个接触表面均不发生滑移,分析中不考虑流体膜剪切弹性模量. 流体膜采用传统分析法. 给出了理论分析和若干变工况参数下的计算结果.      

纳米流体液滴撞击壁面铺展动力学特性研究

纳米流体液滴撞击固体壁面的铺展动力学特性是基于液滴沉积实现高效传热传质过程的关键因素,然而由于纳米流体的非牛顿流变特性及液滴内微流动与纳米颗粒的耦合作用,目前对纳米流体液滴撞击固体壁面的铺展动力学行为缺乏足够的认识.本研究利用了两步法分别配制了分散有3种纳米颗粒的均匀稳定纳米流体(碳纳米管、石墨烯、纳米石墨粉),并对流体的流变特性进行了测量分析.利用显微高速数码摄像技术捕捉了液滴撞击固体壁面的动态过程,通过图像处理技术分析铺展过程中液滴的无量纲高度、铺展因子及动态接触角,探究了液滴在韦伯数约为200及800时撞击壁面后铺展沉积形态的演变规律.研究表明,3种不同纳米颗粒的加入均使基液表现出明显的剪切变稀特性,在液滴撞击壁面的铺展过程中,流体的剪切黏度起重要作用,液滴的无量纲高度和铺展因子的变化幅度随着纳米流体剪切黏度的增大而减小.纳米流体液滴撞击疏水表面时能更快的达到平衡状态,液滴的惯性力主导着液滴的初始铺展阶段,液滴的铺展范围和速度随撞击速度的增大而增大.开展该研究能够为基于液滴沉积的增益冷却技术以及微型高导热及导电材料的制造提供理论依据和技术指导.     

压力驱动微通道流中尺度与扩散关系的研究

对压力驱动微通道流中的扩散进行了数值模拟,说明微尺度效应使流体扩散的不均匀性趋于明显;分析了在各种W/D情况下,距离管道壁面各个层面上的扩散尺度变化,说明了微管道中的对流和壁面效应对于扩散尺度的影响。本文结论对于对微通道的优化设计有参考价值。     

光热激励下微悬臂梁在流体中的振动研究

微悬臂梁结构广泛应用于微纳电子机械系统. 在实际应用中,涂层和工作环境的变化对微悬臂梁结构动态工作模式有着不容忽视的影响. 运用流体中双层微悬臂梁的光热振动模型,研究了在激光光热驱动下,金涂层微悬臂梁在不同流体中的振动特性. 理论上得到了微悬臂梁的温度场,光热驱动力和振动变形场的解析表达式. 研究结果表明,流体环境对微悬臂梁的光热振动谱有显著的影响,主要表现在共振频率的偏移和品质因子的变化两个方面. 相比较于悬臂梁在真空中的响应,当悬臂梁在空气中振动时,共振频率向低频产生微小的漂移（0.7%）,共振峰未发生明显变化;然而,当悬臂梁在液体中振动的时候,共振频率向低频产生巨大的漂移（58%~80%）,而且品质因子发生量级上的减小,共振峰发生了畸变. 本研究对微纳探测以及原子力显微镜等仪器的设计优化,有着一定的理论指导意义.      

基于改进LBM的气液自发渗吸过程中动态润湿效应模拟

微通道内气液自发渗吸是广泛发生在自然界及诸多工业领域的物理现象,而动态接触角是影响整个渗吸过程的关键因素.针对该问题,本文使用改进的伪势多相流格子玻尔兹曼方法 (LBM),直接捕捉微通道内气液自发渗吸过程中的实时接触角,并分析接触角的动态变化特性及其对渗吸长度的影响.首先,本文在原始的伪势多相流LBM的基础上耦合Peng-Robinson (PR)状态方程,改进流体-流体作用力以及流-固作用力格式,并采用精确差分方法将外力添加至LBM框架中.然后,通过校准模型的热力学一致性,模拟测试界面张力,静态平衡接触角等界面现象验证了模型的准确性.最后,基于建立的模拟方法,在水平方向上模拟微通道内气液自发渗吸过程.结果表明:渗吸过程中的接触角呈现动态变化特征,在渗吸初期,因受到惯性力的影响存在较大波动;随着渗吸距离的增大,其逐渐减小并趋近于静态平衡接触角.渗吸过程中的接触角与微通道尺寸及静态接触角有关,随着微通道宽度增大,实时的动态接触角与静态接触角相差大;随着静态接触角增大,实时的动态接触角与静态接触角的相差增大.此外,忽略动态接触角的Lucas-Washburn (LW)方程所预测的弯液面位置...     

纳米尺度下烷烃润滑薄膜分层现象和速度滑移机制的分子模拟研究

采用非平衡分子动力学模拟方法,研究了纳米尺度下十六烷润滑薄膜的分层现象和速度滑移现象,重点考察了剪切速度对速度滑移现象的影响规律,并对其中的微观机理作出了解释.研究结果表明:在铁壁面的限制情况下润滑薄膜出现了分层现象,当润滑薄膜厚度超过50?时,其中间区域呈现出体相均质流体的特征.润滑薄膜层间滑移和界面滑移的临界剪切速度分别为5.5和7.5?/ps,随着剪切速度增加,界面滑移程度增强,而层间滑移程度减弱.润滑薄膜第1和2层十六烷分子层间桥接分子数目决定层间滑移程度,随着剪切速度增加,桥接分子数目也相应增多,层间滑移程度随之减弱.     

Dynamics of boundary layer formation in a volcano channel in a cavitating high-viscosity magma flow

Based on the full mathematical model of a viscous magma melt flow ascending in the gravity field behind a decompression wave front, an unsteady two-dimensional axisymmetric problem of the melt state dynamics at the initial stage of an explosive volcanic eruption and specific features of the flow in the vicinity of the channel wall for the cases of stationary and dynamically increasing viscosity are studied. The evolution of the boundary layer is numerically analyzed for a constant melt viscosity equal to μ = 10 ³ , 10 ⁵ , and 10 ⁷ Pa · sec. It is demonstrated that a boundary layer is formed on the wall of the channel with a radius of 100 m as the melt viscosity is changed in the range of 10 ³ –10 ⁵ Pa · sec, and the boundary layer thickness increases from 2 to 15 m. As the magma viscosity increases to 10 ⁷ Pa · sec, the boundary layer chokes the major part of the channel, thus, locking the flow in the vicinity of the axis of symmetry of the channel almost over the entire channel length. Substantial changes in the flow structure caused by dynamically increasing viscosity are demonstrated by an example of the melt in the channel with a radius of 10 m. By the time t = 1.1 sec, the boundary layer thickness in the channel cross section at a height of approximately 1000 m reaches almost 8 m, the boundary layer acquires the shape similar to a “diaphragm,” penetrates inward the channel by 200 m (with the mass velocity ranging from 0 to 15 m/sec), and locks the flow in a zone with a radius of approximately 2 m around the axis of symmetry of the channel.     

Microscale Visual Study of End Effects at Permeability Discontinuities

The physical effect of multiphase fluid distribution and flow at permeability boundaries has not been fully investigated, particularly at the pore scale (1–100 μm), although such behaviour can significantly affect the overall scaled-up reservoir trapping capacity and production performance. In this article, microscale physical models have been used to qualitatively study the pore scale flow events at permeability boundaries, both high to low and vice versa, to gain a better understanding of the role of these boundaries and water saturation on multiphase displacement behaviour at the pore scale. We have used etched glass models of stripes of large and small (a factor of two) pores with circular matrix. Capillary pressure, which is the controlling parameter is itself dependant on pore size and its spatial distribution, the magnitude of the interfacial tensions and the wettability between the fluids and the solid surface of the models. Sometimes, the only way the non-wetting fluid can penetrate the boundary is through a fortuitous leakage, whereby the presence of an initial saturation reduces the controlling capillary pressure. Examples are demonstrated including mechanisms of end-effects and how capillary boundary resistance (due to capillary forces) can be broken down and fluid movement across the boundary can develop. These micromodel experiments show vividly that connate water can assist in these processes, particularly oil trapping and leakage of water across a permeability boundary.     

Effect of Relaxation and Retardation Time on Peristaltic Transport of the Oldroydian Viscoelastic Fluid

The influence of relaxation and retardation time on peristaltic transport of an incompressible Oldroydian viscoelastic fluid by means of an infinite train of sinusoidal waves traveling along the walls of a two-dimensional flexible channel is investigated. A perturbation solution is obtained for the case in which the amplitude ratio (wave amplitude to channel half-width) is small. The results show that the values of the mean axial velocity of an Oldroydian viscoelastic fluid is smaller than that for a Newtonian fluid. The reflux phenomena are discussed. It is found that the critical reflux pressure gradient decreases with increasing retardation time and increases with increasing relaxation time. Numerical results are reported for different values of the physical parameters of interest. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 86–95, November–December, 2005.     

Unsteady two-fluid flow and heat transfer in a horizontal channel

The problem of unsteady oscillatory flow and heat transfer of two viscous immiscible fluids through a horizontal channel with isothermal permeable walls has been considered. The partial differential equations governing the flow and heat transfer are solved analytically using two-term harmonic and non-harmonic functions in both fluid regions of the channel. Effects of physical parameters such as viscosity ratio, conductivity ratio, Prandtl number and frequency parameter on the velocity and temperature fields are shown graphically. It is observed that the velocity and temperature decrease as the viscosity ratio increases, while they increase with increases in frequency parameter. The effect of increasing the thermal conductivity ratio also suppresses the temperature in both fluid regions. The effect of periodic frequency on the flow is depicted in tabular form. It is predicted that both the velocity and temperature profiles decrease as the periodic frequency increases.     

The Effects of Fluid Viscosity on the Kinematics and Material Properties of <Emphasis Type="Italic">C. elegans</Emphasis> Swimming at Low Reynolds Number

The effects of fluid viscosity on the kinematics of a small swimmer at low Reynolds numbers are investigated in both experiments and in a simple model. The swimmer is the nematode Caenorhabditis elegans, which is an undulating roundworm approximately 1 mm long. Experiments show that the nematode maintains a highly periodic swimming behavior as the fluid viscosity is varied from 1.0 to 12 mPa s. Surprisingly, the nematode’s swimming speed (~0.35 mm/s) is nearly insensitive to the range of fluid viscosities investigated here. However, the nematode’s beating frequency decreases to an asymptotic value (~1.7 Hz) with increasing fluid viscosity. A simple model is used to estimate the nematode’s Young’s modulus and tissue viscosity. Both material properties increase with increasing fluid viscosity. It is proposed that the increase in Young’s modulus may be associated with muscle contraction in response to larger mechanical loading while the increase in effective tissue viscosity may be associated with the energy necessary to overcome increased fluid drag forces.     

Compressible pulsating convection through regular and random porous media: the thermoacoustic case

The effects of material, geometry, length and position of the porous channels on energy transfer in air-filled enclosures carrying a compressible pulsating wave are investigated. The pulsating fluid motion is created by an acoustic driver in a resonant chamber. Three different porous materials (Corning Celcor, Reticulated Vitreous Carbon (RVC), and Mylar plastic), three different geometries (square, open foam, and circular cross-section), six different lengths, “L” (varying between 1 and 6.5 cm, L = 0.01–0.068 λ, where λ is the wavelength of the fundamental acoustic mode), and eight different positions (hot end of the channel, varying between 0.5 and 8 cm) of the channels from the pressure anti-node is experimentally measured. The surface temperature distribution on the channel wall and temperature difference generated across the channel walls are measured while energy flow along the channel walls is calculated analytically. The experimental results are compared with a 1-D numerical code and found excellent agreement. The material, geometry, length, and position of the porous channel strongly affect the energy interactions between the porous channel and the working fluid. The temperature difference generated across the porous RVC channel increases as the porosity increases form 20 to 80 PPI; but decreases if the porosity increases further. Corning Celcor shows improved temperature difference generated across the channel as the length of the channel increases; but then decreases if the length is further increased. The results of this study are applicable to the design of thermoacoustic devices.     

Instability of the interface in two-layer flows with large viscosity contrast at small Reynolds numbers

The Kelvin–Helmholtz instability is believed to be the dominant instability mechanism for free shear flows at large Reynolds numbers. At small Reynolds numbers, a new instability mode is identified when the temporal instability of parallel viscous two fluid mixing layers is extended to current-fluid mud systems by considering a composite error function velocity profile. The new mode is caused by the large viscosity difference between the two fluids. This interfacial mode exists when the fluid mud boundary layer is sufficiently thin. Its performance is different from that of the Kelvin–Helmholtz mode. This mode has not yet been reported for interface instability problems with large viscosity contrasts.These results are essential for further stability analysis of flows relevant to the breaking up of this type of interface.     

Effects of variable fluid properties on natural convection of MHD fluid from a heated vertical wavy surface

In the present paper, the influence of temperature-dependent fluid properties, density, viscosity and thermal conductivity on MHD natural convection flow from a heated vertical wavy surface is studied. It is assumed that, the fluid density and the thermal conductivity vary as exponential and linear functions of temperature, respectively. However, the fluid viscosity is assumed to vary as a reciprocal of a linear function of temperature. The model analysis used here is more relevant to liquid flow. Using the appropriate variables, the wavy surface are transformed into a flat one. The transformed boundary layer equations are solved numerically, using implicit-Chebyshev pseudospectral method, for several sets of values of the physical parameters, namely, the temperature dependent fluid properties parameters, the magnetic parameter, the amplitude-wavelength ratio parameter, and the Prandtl number. The numerical values obtained for the velocity, temperature, shearing stress, and the Nusselt number are presented through graphs and tables for several sets of values of the parameters. The effects of the physical parameters on the flow and heat transfer characteristics are discussed. The results were compared with numerical solutions of previous works. The present results are found to be in good agreement.     

The effect of material properties on the performance of a new geometry PEM fuel cell

In this paper a computational dynamics model for duct-shaped geometry proton exchange membrane (PEM) fuel cell was used to investigate the effect of changing gas diffusion layer and membrane properties on the performances, current density and gas concentration. The proposed model is a full cell model, which includes all the parts of the PEM fuel cell, flow channels, gas diffusion electrodes, catalyst layers and the membrane. Coupled transport and electrochemical kinetics equations are solved in a single domain; therefore no interfacial boundary condition is required at the internal boundaries between cell components. This computational fluid dynamics code is used as the direct problem solver, which is used to simulate the 2-dimensional mass, momentum and species transport phenomena as well as the electron- and proton-transfer process taking place in a PEMFC that cannot be investigated experimentally. The results show that by increasing the thickness and decreasing the porosity of GDL the performance of the cell enhances that it is different with planner PEM fuel cell. Also the results show that by increasing the thermal conductivity of the GDL and membrane, the overall cell performance increases.