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

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

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

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

微通道周期流动电位势及电粘性效应

求解了双电层的Poisson-Boltzmann方程和流体运动的Navier-Stokes方程,得到在周期压差作用下,二维微通道的周期流动电位势,流动诱导电场和液体流动速度的解析解．量纲分析表明,流体电粘性力与以下3个参数有关:1) 电粘性数,它表示定常流动时,通道最大电粘性力与压力梯度的比;2) 形状函数,它表示电粘性力在通道横截面的分布形态; 3) 耦合系数,它表示电粘性力的振幅衰减特征和相位差．分析结果表明,微通道周期流动诱导电场、流动速度与频率Reynolds数有关．在频率Reynolds数小于1时,流动诱导电场随频率Reynolds数变化很慢．在频率Reynolds数大于1时,流动诱导电场随频率Reynolds数的增加快速衰减．在通道宽度与双电层厚度比值较小情况下,电粘性效应对周期流动速度和流动诱导电场有重要影响．     

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

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

等边三角形微通道内层流的流动特性和换热特性的研究

研究等边三角形截面微通道内充分发展层流的流动特性和换热特性,基于Navier-Stokes方程的基本理论,在等边三角形一边向流体加入定常热流密度时,给出了微通道内充分发展层流的速度分布和温度分布的近似解,以及微通道内充分发展对流传热的摩擦因子和Nusselt数;并通过商业软件Fluent对微通道内的流动和换热进行数值模拟,得到通道内温度和速度的数值解,进而计算得到充分发展对流传热的摩擦因子和Nusselt数;二者进行对比,结果吻合很好,验证了计算结果的正确性.     

滑移垂直壁面驻点附近的混合对流边界层流动

就粘性不可压缩流体,研究垂直壁面的滑移,对壁面驻点附近稳定混合对流边界层流动的影响.假定表面温度和外部流动速度与到驻点的距离呈线性变化.首先,将偏微分的控制方程,转变为常微分方程组,然后应用打靶法进行数值求解.对不同数值的控制参数,按分顺流和逆流两种情况,分析和讨论了流动特性和热传导特征.结果表明,逆流时,在浮力参数的某一范围内出现双解;顺流时,解是唯一的.一般而言,速度滑移导致壁面热传导率增大,而热滑移使之减小.     

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

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

微通道中电渗流及微混合的离子浓度效应

电渗流广泛应用于微流控芯片中的流体输运与混合.该文提出了一种离子浓度梯度对电渗流及微混合产生影响的变量模型,采用有限元分析方法对微通道中电渗流及微混合的离子浓度效应进行了数值模拟,分别讨论了zeta电势、介电常数等对微通道内流场和浓度场的影响规律,定量分析了微混合效率.结果表明,当zeta电势和介电常数随浓度变化时,微通道中流场分布不均匀,离子分布不对称.当溶液浓度趋近1 mol/L时,溶液基本无法进入微通道.微混合效率随溶液间浓度差的增大而减小,而且浓度差越大越能在较短距离内到达充分混合.     

矩形微通道散热器流道的数值模拟及尺寸优化

微通道散热器具有体积小、流速小、压降小、散热高等优点,随着工业微型化的发展,微型散热器的应用越来越广泛.已有的研究表明,微通道的散热性能主要决定于微通道的几何参数和流体的流动情况,相对于三角形和梯形结构,矩形微通道具有更好的散热性能.基于ANSYS Workbench有限元软件,对长度为40 mm,不同截面尺寸的单通道内流体流动及传热性能进行了数值模拟,给出具有较小压降、较大散热效率的微通道尺寸.对优化后的模型计算分析,在一定流体流速和温度的初始状态下,基底给一定热通量,经过计算,散热器可运输的热通量较高,压降较低,热传递效率较大,散热器具有良好的工作性能.     

周期壁面电势调制下平行板微管道中的电磁电渗流动

研究了平行板微管道中二维磁流体(MHD)电渗流(EOF)在zeta电势调制下的流动.流体的流动是由两个外加水平电场和垂直磁场所产生的Lorentz力和电场力的组合驱动的.在滑移边界条件下,得到了流函数以及速度分布的解析解.详细讨论了速度随Hartmann数Ha、滑移长度B、电动宽度K等相关的无量纲参数量级变化的变化规律.结果表明,调制的壁面电势会产生一个垂直速度分量,从而导致涡旋的形成.此外,可以观察到,速度的大小随着滑移长度B和电动宽度K的增大而增大.值得注意的是,速度的大小随着Ha值的增大而减小,这与一维流动中Ha值存在临界值的情况不同.     

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

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

Electroosmotic flow and heat transfer in a heterogeneous circular microchannel

Asymmetries in boundary condition are inevitable in practice in microfluidic channels, despite being rarely addressed from theoretical perspectives. Here, by arriving at closed form analytical solutions, we bring out a unique coupling between asymmetries in surface charge and heat transfer in electroosmotically driven microchannel flows. For illustration, we assume that the channel is laterally composed of two parts, each having specified values of the zeta potential and the wall heat flux. Considering low zeta potentials, we obtain analytical solutions in terms of infinite series for the dimensionless forms of the electric potential, the velocity, and the temperature distributions. We demonstrate that, by carefully adjusting the governing parameters, a variety of flow patterns may be achieved, a property that is crucial in applications such as liquid-phase transportation and mixing. Moreover, we show that the average velocity is a linear function of both the zeta potential ratio and the coverage factor. We further show that the average Nusselt number increases when part of the channel having the larger heat flux enlarges and the zeta potential of the part having the smaller surface charge increases. Hence, the maximum heat transfer rates are achieved when the boundary conditions are symmetrical.     

Magnetogasdynamic natural convection in a long vertical microchannel

Since the transport behavior of ionized gases at the microscale could be influenced by an applied magnetic field with ease, microscale magnetogasdynamics (MGD) promises to be particularly advantageous for magnetically controllable microfluidic devices. The purpose of this study is to investigate how magnetic force affects the MGD natural convection within a long asymmetrically heated vertical planar microchannel. The fully developed solutions of the thermal-flow fields and their characteristics are analytically derived on the basis of the first-order slip and jump boundary conditions and then presented for the thermophysical properties of ionized air at the standard reference state flowing through the microchannel with complete accommodation. The calculated results reveal that magnetic force plays a damping role in flow and results in decreases in flow rate, average flow drag, and average heat transfer rate. In addition, it is interesting that because the flow near the core is suppressed and the shear stress on the wall surface is reduced by the magnetic effects, a flatter velocity profile could be achieved by a greater magnetic force. These magnetic effects could be further magnified by increasing gas rarefaction or increasing cooler wall temperature.     

Effect of variable viscosity and viscous dissipation on the Hagen‐Poiseuille flow and entropy generation

In this article, the steady‐state flow of a Hagen‐Poiseuille modelin a circular pipe is considered and entropy generation due tofluid friction and heat transfer is examined. Because of variationin fluid viscosity, the entropy generation in the flow varies. Inhis model, Arrhenius law is applied for temperature equation‐dependent viscosity, and the influence of viscosity parameters on the entropy generation number and distribution of temperature and velocity is investigated. The governing momentum and energy equations, which are coupled due to the dissipative term in the energy equation, were solved by analytical techniques. The solutions of equations via perturbation method and homotopy perturbation method are obtained and then compared with those of numerical solutions. It is found that the fluid viscosity influences considerably the temperature distribution in the fluid close to the pipe wall, and increasing pipe wall temperature enhances the rate of entropy generation. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 529–540, 2011     

Influence of radiation on MHD free convection from a vertical flat plate embedded in porous media with thermophoretic deposition of particles

Magneto-hydrodynamics and thermal radiation effects on heat and mass transfer in steady laminar boundary layer flow of a Newtonian, viscous fluid over a vertical flat plate embedded in a fluid saturated porous media in the presence of the thermophoresis particle deposition effect is studied in this paper. The governing equations are transformed by special transformations. Brownian motion of particles and thermophoretic transport are considered in the flow equations. The magnetic field is considered to be applied. Rosseland approximation is used to describe the radiative heat flux in the energy equation. The resulting similarity equations are solved numerically by the fourth-order Runge–Kutta method with shooting technique. Many results are obtained and representative set is displayed graphically to illustrate the influence of the various parameters on the wall thermophoretic deposition velocity, concentration, temperature and velocity profiles.     

Heat and mass transfer in MHD non-Darcian flow of a micropolar fluid over a stretching sheet embedded in a porous media with non-uniform heat source and thermal radiation

A mathematical analysis has been carried out to study magnetohydrodynamic boundary layer flow, heat and mass transfer characteristic on steady two-dimensional flow of a micropolar fluid over a stretching sheet embedded in a non-Darcian porous medium with uniform magnetic field. Momentum boundary layer equation takes into account of transverse magnetic field whereas energy equation takes into account of Ohmic dissipation due to transverse magnetic field, thermal radiation and non-uniform source effects. An analysis has been performed for heating process namely the prescribed wall heat flux (PHF case). The governing system of partial differential equations is first transformed into a system of non-linear ordinary differential equations using similarity transformation. The transformed equations are non-linear coupled differential equations which are then linearized by quasi-linearization method and solved very efficiently by finite-difference method. Favorable comparisons with previously published work on various special cases of the problem are obtained. The effects of various physical parameters on velocity, temperature, concentration distributions are presented graphically and in tabular form.