微通道液体流动双电层阻力效应

采用数值方法求解双电层的Poisson-Boltzmann方程和液体运动的Navier-Stokes方程,研究微通道双电层对压强梯度液体流动的阻力效应． 量纲分析表明,双电层阻力大小可以用一个无量纲的电阻力数表示．它与液体的介电系数、固体表面的zeta电位平方成正比,与液体的动力粘性系数、电导率以及微通道的宽度平方成反比．在计算流动诱导的流动电位势和电阻力时,提出电流密度平衡条件,可以消除传统电流平衡条件导致的固壁附近产生局部回流的不合理物理现象．还给出不同电阻力数的微通道流量、流量损失率、速度剖面的数值结果,合理解释了双电层对微通道液体流动的阻力效应．     

异质材料有限长微通道电渗流热效应

采用数值方法,分析有限长PDMS/玻璃微通道电渗流热效应．数值求解双电层的Poisson-Boltzmann方程,液体流动的Navier-Stokes方程和流-固耦合的热输运方程,分析二维微通道电渗流的温度特性．考虑温度变化对流体特性（介电系数、粘度、热和电传导率）的反馈效应．数值结果表明,在通道进口附近有一段热发展长度,这里的流动速度、温度、压强和电场快速变化,然后趋向到一个稳定状态．在高电场和厚芯片的情况下,热发展长度可以占据相当一部分的微通道．电渗流稳定态温度随外加电场和芯片厚度的增加而升高．由于壁面材料的热特性差异,在稳定态时的PDMS壁面温度比玻璃壁面温度高．研究还发现在微通道的纵向和横向截面有温度变化．壁面温升降低双电层电荷密度．微通道纵向温度变化诱发流体压强梯度和改变微通道电场特性．微通道进流温度不改变热稳定态的温度和热发展长度．     

微通道内电渗压力混合驱动幂律流体流动模拟

为了研究微通道内电渗压力混合驱动幂律流体的流动特性,建立了微通道内电渗压力混合驱动幂律流体的计算模型,其双电层电势、流体的流场分布分别由Poisson-Boltzmann(P-B)方程和Navier-Stokes(N-S)方程描述.讨论了无量纲Debye(德拜)参数K、壁面ζ*电势和幂律指数n对流体流动特性和Poiseuille数的影响.结果表明,当压力梯度与外加电场方向一致(Γ0)时,剪切变稀流体的速度大于剪切变稠流体;压力梯度与外加电场方向相反(Γ0)时,结果相反.Poiseuille数是无量纲Debye常数K、壁面ζ*电势和幂律指数n的增函数.     

矩形微通道热沉内单相稳态层流流体的流动与传热分析

采用解析方法分析了矩形微通道热沉内单相稳态层流流体的流动与传热．基于y方向流速和导热不变的假设,建立流体在矩形微通道内流动的流速方程和传热的温度方程,进而推导出Nusselt数和Poiseuille数的理论表达式．通过计算结果可以看出,推导的Nusselt数和Poiseuille数的解析解与其他文献的结果吻合较好,而且当宽高比趋于无穷大时,Nusselt数和Poiseuille数分别趋近于8.235和96,这与其他文献结果完全相同．在Reynolds数相同时,摩擦因数随着宽高比的增加而增加,而在相同宽高比时,摩擦因数随Reynolds数的增加而减小．     

具有正弦粗糙度的环形微管道中脉冲流动

研究了具有正弦粗糙度的环形微管道中脉冲流动,其中壁面粗糙度用小振幅的正弦波表示,不可压缩粘性脉冲流动由周期振荡的压力梯度驱动,运用摄动展开法求解了柱坐标系下的动量方程,获得了环形微管道内脉冲流动的近似解析速度及其体积流率.在此基础上,研究了相关无量纲参数,如Reynolds(雷诺)数Re、压力梯度振幅A、正弦波状粗糙的小振幅ε、内外半径之比α、相位差β及其波数λ对速度u及平均体积流率Φ_m的影响.结果表明,剖面速度随A的增大而增大,随Re的增大而减小,相位滞后χ随振荡Reynolds数Re的增大而增大.     

两平行圆盘间或两圆球间存在二阶流体时的挤压流动

为了进行湿颗粒群的离散元模拟,研究两圆球颗粒间二阶流体在挤压流动时的法向粘性力．首先用小参数法对两平行圆盘间二阶流体挤压流动的速度场和正应力分布进行了近似分析,然后用类似的方法,分析任意两圆球间二阶流体的挤压流动,得到了压力分布和法向粘性力的解析解．     

微极流体向受热面的MHD驻点流动

分析了有均匀横向磁场作用时,导电微极流体垂直冲击受热面时形成的二维驻点流动问题．应用适当的相似转换,将连续、动量、角动量及热量的控制方程,及其相应的边界条件,简化为无量纲形式．然后,利用以有限差分离散化为基础的算法,求解简化了的自相似非线性方程．用Richardson外推法,进一步求精其结果．以图表形式表示磁场参数、微极性参数、Prandtl数对流动和温度场的影响,说明了其解的重要特性．研究表明,随着磁场参数的增大,速度和热边界层厚度变小了．与Newton流体相比较,微极流体的剪应力和传热率出现明显的减少,这对聚合物生产过程中流体的流动和热量控制是有益的．     

电渗流中传热传质过程与熵的分析

在壁面存在恒定热通量条件下,分析微通道内电渗流中传热传质过程与熵的生成.建立数值计算模型,分别采用Poisson-Boltzmann方程、Navier-Stokes方程、Nernst-Planck方程和能量方程来描述微通道内双电层电势、流场、离子浓度和温度的分布情况.引入熵产生,进一步研究不同流动参数对流体传热过程的作用,讨论不同流动参数下各热效应的变化规律,并具体分析热效应参数对流体总熵增加及各部分热效应对总熵比重的影响.结果表明,动电参数与Joule(焦耳)热系数的增大会使得传热性能减弱,动电参数对传热性能影响更为明显;流体的总熵为动电参数、传质系数和质量弥散系数的增函数.     

计及Hall和离子滑移电流的旋转流体在与时间相关初值条件下的MHDCouette流动

考虑Hall和离子滑移电流的影响,在旋转系统中研究导电流体非稳定的MHD Couette流动．在小数值磁场Reynolds数假定下,推导出基本的控制方程,使用著名的Laplace变换技术,数值地求解该基本方程．分两种情况:磁场相对于流体或者移动平板固定时,得到速度和表面摩擦力统一的闭式表达式．用图形讨论了问题的不同参数,对速度和表面摩擦力的影响．所得结果显示,主流速度u和次生速度v随着Hall电流而增大．离子滑移电流的增大,也会导致主流速度u的增大,但会使次生速度v减小．还给出了旋转、Hall和离子滑移参数的综合影响,确定了次生运动对流体流动的贡献．     

微极液体在两个多孔圆盘间的MHD流动及其热传导

两个平行的无限大多孔圆盘,圆盘表面有均匀注入时,数值地研究圆盘间不可压缩导电微极流体,在横向外加磁场作用下的轴对称稳定层流.运用von Krmn的相似变换,将非线性运动的控制方程转化为无量纲形式.使用基于有限差分格式的算法,在相应的边界条件下,求解简化后耦合的常微分方程组.讨论Reynolds数、磁场参数、微极参数和Prandtl数,对流动速度和温度分布的影响.在特殊情况下,所得结果与已有文献的工作有着很好的一致性.研究表明,圆盘表面的传热率随着Rynolds数、磁场参数和Prandtl数的增加而增加;剪切应力随着注入的增加而减少,但它随着外部磁场的加强而增加.和Newton流体相比较,微极流体的剪切应力因素较弱,有利于聚合体加工过程中流动和温度的控制.     

储层岩石周期性电渗流特性的微观机制

以双电层电位理论和电渗流动的动量方程为基础,结合储层岩石平行毛管束模型,推导出岩石孔隙内周期性电渗流的解析式,揭示了储层中电渗效应的微观机制,分析了非密闭储层岩石中宏观电渗Darcy速度及密闭储层中电渗压力系数频散特性的影响因素．数学模拟结果表明:储层岩石孔隙中,周期性电渗流速度剖面在频率较高时呈“波浪”状;孔隙度越大,电渗Darcy速度模值越大,其相位也越大,而电渗压力系数数值越小．储层岩石的溶液浓度越小或阳离子交换量越大,电渗Darcy速度模值和电渗压力系数数值越大,但对电渗Darcy速度的相位没有影响．     

柔性圆柱形微管道内的电动流动及传热研究

研究了在纯压力驱动下,流体通过壁面带有某种电荷的聚电解质层（PEL）的微管道，即柔性微管道的电动流动和热传输特性.基于先前得到的电势和速度的解析解以及流向势的数值解,在热充分发展的情况下, 假设壁面热流恒定,利用有限差分法求解了包括黏性耗散和Joule(焦耳)热影响下的能量方程，获得了无量纲温度数值解.通过数值计算，给出了相关的无量纲参数对速度、温度以及Nusselt(努赛尔)数的影响.研究表明,当其他参数固定时,无量纲速度和温度随着无量纲聚电解质层厚度d的增大而减小,随着聚电解质层中等效双电层厚度与双电层厚度之比K_λ的增大而增大；Nusselt数随着Joule热系数S的增大而减小,随无量纲聚电解质层厚度d的增大而减小,随着K_λ的增大而增大.     

Heat transfer enhancement for boundary layer flow over a wedge by the use of electric fields

The effect of an electric field is investigated for heat transfer properties in a laminar, incompressible, non-isothermal boundary layer gas flow over a wedge. The governing boundary layer equations are reduced to an ordinary differential equation system using similarity transformations. The reduced equations are solved numerically for different values of electric and flow field parameters characterizing the ratio of electric force to fluid inertia force, Joule heating and ion kinetic work. For specific electric field function forms, leading to similarity solutions, velocity boundary layers are observed to become thinner and heat transfer properties are shown to be enhanced near the wall. The level of enhancement is controlled by the electric body force with additional effects of Joule heating and ion kinetic work on the bulk flow. The effects of low and high Prandtl numbers are also demonstrated. Heat transfer enhancement is observed to increase with increasing Prandtl number.     

The applicability of continuum models in the transitional regime of hypersonic flow over blunt bodies

Hypersonic rarefied gas flow over blunt bodies in the transitional flow regime (from continuum to free-molecule) is investigated. Asymptotically correct boundary conditions on the body surface are derived for the full and thin viscous shock layer models. The effect of taking into account the slip velocity and the temperature jump in the boundary condition along the surface on the extension of the limits of applicability of continuum models to high free-stream Knudsen numbers is investigated. Analytic relations are obtained, by an asymptotic method, for the heat transfer coefficient, the skin friction coefficient and the pressure as functions of the free-stream parameters and the geometry of the body in the flow field at low Reynolds number; the values of these coefficients approach their values in free-molecule flow (for unit accommodation coefficient) as the Reynolds number approaches zero. Numerical solutions of the thin viscous shock layer and full viscous shock layer equations, both with the no-slip boundary conditions and with boundary conditions taking into account the effects slip on the surface are obtained by the implicit finite-difference marching method of high accuracy of approximation. The asymptotic and numerical solutions are compared with the results of calculations by the Direct Simulation Monte Carlo method for flow over bodies of different shape and for the free-stream conditions corresponding to altitudes of 75–150 km of the trajectory of the Space Shuttle, and also with the known solutions for the free-molecule flow regine. The areas of applicability of the thin and full viscous shock layer models for calculating the pressure, skin friction and heat transfer on blunt bodies, in the hypersonic gas flow are estimated for various free-stream Knudsen numbers.     

High-Reynolds-number effects on turbulent scalings in compressible channel flow

The effect of the Reynolds number in a supersonic isothermal channel flow is studied using a direct numerical simulation (DNS). The bulk Mach number based on the wall temperature is 1.5, and the bulk Reynolds number is increased up to Re_τ ≈︁ 1000. The use of van Driest velocity transformation in the presence of heated walls has been questioned due to the poor accuracy at low Reynolds number. For this reason alternative transformations of the velocity profile and turbulence statistics have been proposed, as, for instance, semi-local scalings. We show that the van Driest transformation recovers its accuracy as the Reynolds number is increased. The Reynolds stresses collapse on the incompressible ones, when properly scaled with density, and very good agreement with the incompressible stresses is found in the outer layer. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effect of an external force on the by-pass transition mechanism in internal flows electrical double layer effect in mini and micro channels

The effect of the electric double layer (EDL) on the bypass transition mechanism is explored through direct numerical simulations. The electrokinetic effects destabilize rapidly the flow when the local disturbance and/or the Reynolds number are respectively strong and large enough to overcome the transient growth regime. It is found that a weak perturbation quickly leads to the transition through bypass and non-linear interactions under the EDL effect, while an order of magnitude larger disturbance is incapable to destabilize the macro-scale flow. The EDL develops some new transitional wall structures during the bypass process. It is concluded that electrostatic effects can be efficiently used to enhance the mixing and heat transfer in microchannels, providing that the diffuse layer is large enough.     

Laminar magnetohydrodynamic duct flow in the presence of a magnetic dipole

We study magnetohydrodynamic flow of a liquid metal in a straight duct. The magnetic field is produced by an exterior magnetic dipole. This basic configuration is of fundamental interest for Lorentz force velocimetry (LFV), where the Lorentz force opposing the relative motion of conducting medium and magnetic field is measured to determine the flow velocity. The Lorentz force acts in equal strength but opposite direction on the flow as well as on the dipole. We are interested in the dependence of the velocity on the flow rate and on strength of the magnetic field as well as on geometric parameters such as distance and position of the dipole relative to the duct. To this end, we perform numerical simulations with an accurate finite-difference method in the limit of small magnetic Reynolds number, whereby the induced magnetic field is assumed to be small compared with the external applied field. The hydrodynamic Reynolds number is also assumed to be small so that the flow remains laminar. The simulations allow us to quantify the magnetic obstacle effect as a potential complication for local flow measurement with LFV. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Mixed Analytical/Numerical Method for Velocity and Heat Transfer of Laminar Power-Law Fluids

This paper presents a relatively simple numerical method to investigate the flow and heat transfer of laminar power-law fluids over a semi-infinite plate in the presence of viscous dissipation and anisotropy radiation. On one hand, unlike most classical works, the effects of power-law viscosity on velocity and temperature fields are taken into account when both the dynamic viscosity and the thermal diffusivity vary as a power-law function. On the other hand, boundary layer equations are derived by Taylor expansion, and a mixed analytical/numerical method (a pseudo-similarity method) is proposed to effectively solve the boundary layer equations. This method has been justified by comparing its results with those of the original governing equations obtained by a finite element method. These results agree very well especially when the Reynolds number is large. We also observe that the robustness and accuracy of the algorithm are better when thermal boundary layer is thinner than velocity boundary layer.