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

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

输气管道壁面涂料减阻机理的实验研究

用IFA-300热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中不同壁面涂料的管道湍流边界层不同法向位置流向速度分量的时间序列信号,利用湍流边界层近壁区域对数律平均速度剖面与壁面摩擦速度、流体黏性系数等内尺度物理量的关系和壁面摩擦速度与壁面摩擦切应力的关系,在准确测量湍流边界层近壁区域对数律平均速度剖面的基础上,间接测量湍流边界层的壁面摩擦阻力．对不同壁面涂料的壁湍流脉动速度信号用子波分析进行多尺度分解,用子波系数的瞬时强度因子和平坦因子检测管道湍流边界层中的多尺度相干结构,提取不同尺度相干结构的条件相位平均波形,对比研究输气管道壁面涂料的减阻机理．     

Richtmyer-Meshkov不稳定性扰动增长率冲击模型的微尺度效应

本文通过考虑流体介质输运性质对激波理论的影响,针对微型激波管中的RichtmyerMeshkov(RM)不稳定性,分析了微尺度效应对三种冲击模型扰动增长率的影响。一维气体动力学计算结果表明:流动尺度对于Richtmyer模型、Meyer-Blewett(M-B)模型和Vandenboomgaerde-Mügler-Gauthier(V-M-G)模型中的线性扰动增长率有显著影响。当流动尺度由宏观状态逐渐减小至微尺度时,三种模型的扰动增长率均会经历从少量增长到明显下降,然后迅速上升的过程。微尺度条件下,V-M-G模型的扰动增长率相比宏观尺度有显著提高。与其他两种模型相比,修正后的V-M-G模型更合理地描述了微尺度效应对于线性扰动增长率的影响。此外,对于修正后的V-M-G模型,当入射激波马赫数较低时,扰动增长率受微尺度效应的影响更为明显。     

微型矩形凹槽的长方体通道的减阻和传热特性

在长方体通道底面沿展向方向间隔设置了微型矩形凹槽，凹槽的深度与边界层尺度相当。采用大涡模拟方法对长方体通道内的流动及传热特性进行研究。数值计算结果表明：在长方体通道内设置的微矩形凹槽可以诱导“突出效应”及二次涡，二次涡的作用类似于微型空气滚动轴承，因而可减小流阻，并使传热性能略有提高。研究表明：微凹槽导致了速度滑移，从而有效降低了通道底部附近流体速度梯度；造成低速条纹变宽，使高低速流体的混合受到抑制。微凹槽内产生的二次涡增加了黏性底层的厚度，且二次涡与微凹槽上方流体之间的滚动摩擦代替了壁面与流体之间的滑动摩擦。与没有布置微型矩形凹槽的长方体通道相比，布置微凹槽的长方体通道可在不影响传热效果的前提下达到6%以上的减阻率。     

微通道壁面突起对片上电泳分离流动特性影响数值计算研究

针对制造过程引起的通道壁面突起对微流控芯片电泳分离使用的影响进行数值计算和理论分析研究.分析了壁面突起产生机制及其对电泳分离效果的影响;阐述了毛细管电泳分离的物理模型并进行了离散化;编制了电泳分离数值计算程序,采用有限体积法计算微通道内电参数分布、缓冲溶液流场分布和样品区带分布;给出微通道壁面突起高度和突起宽度对电泳分离过程影响的计算结果.从计算结果可知:壁面突起高度是影响电渗流流速的主要因素,当壁面突起高度与微通道宽度的比值从0增加至0.2,电渗流流速变化幅度约为20%.     

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

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

页岩气和致密砂岩气藏微裂缝气体传输特性

页岩气和致密砂岩气藏发育微裂缝,其开度多在纳米级和微米级尺度且变化大,因此微裂缝气体传输机理异常复杂.本文基于滑脱流动和努森扩散模型,分别以分子之间碰撞频率和分子与壁面碰撞频率占总碰撞频率的比值作为滑脱流动和努森扩散的权重系数,耦合这两种传输机理,建立了微裂缝气体传输模型. 该模型考虑微裂缝形状和尺度对气体传输的影响. 模型可靠性用分子模拟数据验证.结果表明:(1)模型能够合理描述微裂缝中所有气体传输机理,包括连续流动,滑脱流动和过渡流动;(2)模型能够描述不同开发阶段,微裂缝中各气体传输机理对传输贡献的逐渐变化过程;(3)微裂缝形状和尺度影响气体传输,相同开度且宽度越大的微裂缝,气体传输能力越强,且在高压和微裂缝大开度的情况下表现更明显.      

微扩散管内幂律流体在不同外电场驱动下非稳态流动数值仿真

电场驱动下的非牛顿流体在微米级扩散管道内非稳态电渗流动特性是MEMS管设计人员关注的焦点,大部分实际液体可近似为用幂律模型描述的纯粘性流体,所以论文针对幂律流体在有限长微扩散管道内在两种不同形式的外加电场驱动下的非稳态电渗流动情况进行数值仿真.基于Ostwald-De Wael幂律模型和连续介质假说,采用高精度紧致有限差分离散二维完全Poisson-Boltzmann电势方程和Cauchy动量方程,对恒定电场及满足Maxwell方程的电场进行数值仿真,讨论了微扩散管中幂律流体在两种不同外加电场驱动下的瞬时流场分布的差异.结果表明,初始时刻固定扩散角和无量纲壁面电势,无量纲电动宽度的变化对幂律流体电渗流速度分布影响较大;在微扩散管上游等截面处,由恒定电场驱动及Maxwell电场驱动电渗流速度分布差别极小,在扩散管中下游则出现了明显的差别;由恒定电场驱动下的电渗流动在扩散管不同截面下的速度峰值相近,但Maxwell电场诱导的电渗流速度峰值则随管道半径变化出现较大差别.对于外加电场驱动的电渗流动,不同形式的外电场可使流场产生较大差别,而不同性质的流体也会形成不同的流场分布.     

THERMOELASTIC VIBRATION OF FLUID-CONVEYING MICROTUBES EMBEDDED IN ELASTIC MEDIUMS

研究热环境中被弹性介质包围的微米输流管道的横向振动问题. 根据Hamilton 原理及非线性热弹性理论建立管道横向振动控制方程,并利用复模态法对其进行求解,得到了系统的固有频率和屈曲失稳临界流速,讨论了环境温度和一些重要系统参数对管道振动特性的影响. 研究结果表明:环境温度变化、管道和流体的微尺度效应、管道外径及弹性介质刚度对输流微管道固有频率和临界流速都有很大影响.     

基于Cosserat理论的四边简支自由振动微平板尺度效应研究

微观下材料内部结构将极大地影响材料的力学性能,对微纳米器件中典型的微平板结构尺度效应进行研究具有十分重要的意义.论文基于Cosserat理论推导出了微平板自由振动的微分方程,并根据四边简支边界条件假设振型函数,给出了固有频率的计算公式,对不同尺寸微平板固有频率的尺度效应进行了仿真分析.结果表明,考虑了尺度效应的微平板自由振动固有频率要高于经典理论中的固有频率.当特征长度与微平板厚度大小相当时,微平板固有频率表现出明显的尺度效应,并随着特征长度的增加而增大.同时,自由振动的尺度效应将随着微平板厚度的减小而逐渐增强,振动模态及长宽比不影响尺度效应.论文的研究将为微结构与系统的应用提供一定的理论基础.     

Parameters analysis of fluid diffusion

The influences of fluid density, diffusivity, viscosity, width of the flow channel, travel distance, and flow velocity on fluid diffusion are analyzed theoretically and numerically. Concentration boundary layer is taken as the quantitative index of fluid diffusion in this work. The results show that diffusion is a function of travel distance, diffusivity, fluid density, and flow velocity. Diffusion is independent of the width of the channel. Viscous effect determines the velocity gradient and does not affect diffusion directly. The usually used Péclet number uL/D cannot govern the full condition of fluid diffusion. For two-fluids co-flowing in a two-dimensional straight channel with relative low viscous effect, diffusion is proportional to the square root of travel distance and diffusivity, and is inversely proportional to the square root of flow velocity.     

Diffusion velocity for a three-dimensional vorticity field

A discussion is presented on the existence of a diffusion velocity for the vorticity vector that satisfies extensions of the Helmholtz vortex laws in a three-dimensional, incompressible, viscous fluid flow. A general form for the diffusion velocity is derived for a complex-lamellar vorticity field that satisfies the property that circulation is invariant about a region that is advected with the sum of the fluid velocity and the diffusion velocity. A consequence of this property is that vortex lines will be material lines with respect to this combined velocity field. The question of existence of diffusion velocity for a general three-dimensional vorticity field is shown to be equivalent to the question of existence of solutions of a certain Fredholm equation of the first kind. An example is given for which it is shown that a diffusion velocity satisfying this property does not, in general, exist. Properties of the simple expression for diffusion velocity for a complex-lamellar vorticity field are examined when applied to the more general case of an arbitrary three-dimensional flow. It is found that this form of diffusion velocity, while not satisfying the condition of circulation invariance, nevertheless has certain desirable properties for computation of viscous flows using Lagrangian vortex methods. The significance and structure of the noncomplex-lamellar part of the viscous diffusion term is examined for the special case of decaying homogeneous turbulence.     

Photobleached-fluorescence imaging of microflows

A photobleached-fluorescence imaging technique for visualizing microscale flow fields and obtaining molecular diffusion and advection information has been developed. The technique tracks fluorophores in the region of a photobleached line in a planar microdevice and yields quantitative diffusive and advective transport data. Visualizations of two- and weakly three-dimensional electroosmotically and pressure-driven fluid flow fields are demonstrated using the photobleaching of fluorescein and fluorescein-dextran conjugates. Photobleached-fluorescence imaging tracks undisturbed fluorophores, functions in polymer and glass microfluidic devices, can take advantage of fluorescent conjugates present in biochemical assays, and has a photobleached region that is flow independent.     

Problem of Capillary Displacement for One Model of Three-Phase Filtration

The class of capillary pressures corresponding to a triangular tensor of capillary diffusion in a three-phase fluid is studied. Filtration with such a tensor is described by a parabolic system of equations degenerating on solutions. This system is integrodifferential because the desired quantities are the total flow rate and the phase-saturation distribution under conditions of a specified pressure drop on the boundaries of the flow region in one of the phases. It is shown that in the problem of capillary displacement, the degenerate system can be studied using a special maximum principle.     

Numerical Investigation of Magnetic Nanoparticles Distribution Inside a Cylindrical Porous Tumor Considering the Influences of Interstitial Fluid Flow

A numerical simulation of interstitial fluid flow and blood flow and diffusion of magnetic nanoparticles (MNPs) are developed, based on the governing equations for the fluid flow, i.e., the continuity and momentum and mass diffusion equations, to a tissue containing two-dimensional cylindrical tumor. The tumor is assumed to be rigid porous media with a necrotic core, interstitial fluid and two capillaries with arterial pressure input and venous pressure output. Blood flow through the capillaries and interstitial fluid flow in tumor tissues are carried by extended Poiseuille’s law and Darcy’s law, respectively. Transvascular flows are also described using Starling’s law. MNPs diffuse by interstitial fluid flow in tumor. The finite difference method has been used to simulate interstitial fluid pressure and velocity, blood pressure and velocity and diffusion of MNPs injected inside a biological tissue during magnetic fluid hyperthermia (MFH). Results show that the interstitial pressure has a maximum value at the center of the tumor and decreases toward the first capillary. The reduction continues between two capillaries, and interstitial pressure finally decreases in direction of the tumor perimeter. This study also shows that decreasing in intercapillary distance may cause a decrease in interstitial pressure. Furthermore, multi-site injection of nanoparticles has better effect on MFH.     

Transient thermal entrance heat transfer in laminar pipe flows with step change in pumping pressure

Analysis is made for the transient heat transfer phenomena in the thermal entrance region of laminar pipe flows. The transient results from both the change in flow field, a step change in pressure gradient from zero to a fixed value, and the change in thermal field, a step change in the inlet temperature. An exponential scheme has been employed to solve the energy equation with the presence of axial heat conduction in the fluid. In order to demonstrate the results more clearly, a modified Nusselt number is introduced. The unsteady axial variations of conventional Nusselt number, modified Nusselt number, bulk fluid temperature and pipe wall temperature are presented for water and air over a wide range of outside heat transfer coefficients. It is observed that the outside heat transfer coefficient has a significant influences on the transient heat transfer processes. The results can be comprehensively interpreted by the interactions among the axial convection, axial diffusion, and radial diffusion.     

一种骨小管中液体流动产生的流量及切应力模型

骨组织受力变形后其内部液体就会流动,同时在其微观结构——骨单元壁中扩散,并进一步产生一系列与骨液流动相关的物理效应,如流体剪切应力、流动电位等,这些物理效应被细胞感知并做出破骨或成骨等反应,来使骨适应外部载荷环境.鉴于骨组织产生的内部液体流动很难实验测定,理论模拟是目前的主要研究手段.基于骨单元的多孔弹性性质建立了骨小管内部液体的流动模型,该模型将骨单元所受的外部载荷与骨小管内部液体的压力、流速、流量和切应力联系起来,并进一步可以研究其力传导与力电传导机制.骨小管模型的建立分别基于中空和考虑哈弗液体的骨单元模型,并考虑了骨单元外壁的弹性约束和刚性位移约束两种边界条件.最终得到骨单元在外部轴向载荷作用下,骨小管内部液体的流量及流体切应力的解析解.结果表明:骨小管中的液体流量与流体切应力都正比于应变载荷幅值和频率,并由载荷的应变率决定.因此应变率可以作为控制流量和流体切应力的一种生理载荷因素.流量随着骨小管半径的增大而非线性增大,而流体切应力则随着骨小管半径的增大而线性增大.此外,在相同的载荷下,含哈弗液体的骨单元的模型中,骨小管中液体的流量和切应力均大于中空骨单元模型.     

Mass and heat transfer during solid sphere dissolution in a non-uniform fluid flow

We studied a nonisothermal dissolution of a solvable solid spherical particle in an axisymmetric non-uniform fluid flow when the concentration level of the solute in the solvent is finite (finite dilution of solute approximation). It is shown that simultaneous heat and mass transfer during solid sphere dissolution in a uniform fluid flow, axisymmetric shear flow, shear-translational flow and flow with a parabolic velocity profile can be described by a system of generalized equations of convective diffusion and energy. Solutions of diffusion and energy equations are obtained in an exact analytical form. Using a general solution the asymptotic solutions for heat and mass transfer problem during spherical solid particle dissolution in a uniform fluid flow, axisymmetric shear flow, shear-translational flow and flow with parabolic velocity profile are derived. Theoretical results are in compliance with the available experimental data on falling urea particles dissolution in water and for solid sphere dissolution in a shear flow.     

Vermischung von zwei Fluiden durch Diffusion in einer Potentialströmung

The problem may be formulated as follows: “Subject to investigation is a plane flow in a long straight channel or duct of constant width. A fluid enters the channel at one of its ends. Another fluid is supplied by two sources located on opposit sides of the duct. The flow being a potential one, the process of mixing the fluids is mere diffusion. The direction of diffusion considered is perpendicular to the direction of flow. The mixture leaves the channel at its other end. The aim is to determine the steady state concentration field inside the duct. “ The solution given bases ona generalisation according to the principles of conformity. This involves the introduction of “confrom models.?