直流微装置的磁流体动力学模拟

第十四章直流微装置的磁流体动力学模拟JAIMEH.LOZANOPARADA,WILLIAMB.J.ZIMMERMANDepartmentofChemicalandProcessEngineering,UniversityofSheffield,NewcastleStreet,SheffieldS13JDUnitedKingdomE-mail:[e-mail]1．引言磁流体动力学理论（MHD）研究电磁场中导电流体的交互作用。它在很多领域，包括热核反应[1]、太阳和太空等离子体[2]、火箭引擎[3]中都有着非常重要的作用。目前对MHD的研究兴趣越来越集中在芯片实验中的微尺度流动控制应用上[4]-[6]。驱动MHD微尺度泵的Lorentz力，在方向和大小上取决于施加的磁场B和电场E矢量。这种泵的主要特性就是可以控制局部流体流动，不需要力学设备就可以精确控制流体在微尺度流道网络中按照预定路径流动[7]。这种借助Lorentz力的局部流体控制方法使得流体控制变得十分灵活，例如流体可以双向流动、累积、减速甚至回退。与电动泵使用高的轴线电压相比[8]，MHD微型泵使用低的横向电场。低的发热量使其可以用于驱动对高温和电压敏感的生物流动过程。简单的电子设备就可以顺序控制复杂微流动中的各个独立微型泵。流动速度通过电磁场的强度来控制。似乎到目前为止仍没有关于MHD微型泵模拟的发表文章。下面我们将给出一些基于Galerkin有限元法的微型泵模拟结果，模拟过程在商业软件COMSOLMultiphysics3.2中实现。数值求解采用压力修正算法——SIMPLE，它首先假设一个压力场，然后通过求解不可压缩流动的Navier-Stokes方程得到速度场。这些速度不需要满足Possion型连续方程，所以对压力场的修正也带来速度场的修正，最终满足质量守恒。求解速度场的同时计算电势场方程。这会得到Lorentz力，然后将其反馈回N-S方程并作为体积力处理。连续耦合Lorentz力和速度场直到Newton迭代收敛。2．理论为了研究电场和磁场对微通道反应器（MR）里中性导电流体的影响，我们对不可压缩各向同性牛顿流体建立一个MHD方程。在实验室尺度，MHD控制方程可以解耦合为三个主要问题[9]：一是基于质量守恒和动量守恒的连续方程流动问题；二是基于Maxwell方程的电动力学问题；三是基于能量守恒的热问题。但是流体动力学和电动力学的解耦合是存在问题的，因为定义了共同的压力边界条件。不过也很容易发现，当流动速度足够慢的时候，压力分布可以一阶近似为静流体情况，这样就可以定义一个纯磁流体动力学问题。这样就考虑了磁场对速度场的影响，但是忽略了流场对磁力线对流的影响（这就是著名的弱耦合），除非入口流速增大到磁场Reynolds数Rm比1大。当Rm≈1时，磁力线对流约占静磁场B0的1％，但是如果Rm≈10时，磁力线扰动带来的磁场近似等于没有扰动的磁场B0，也就是说这时磁力线的对流变得不可忽略。对于我们这种情况，如果产生磁场线对流需要入口平均流速u0≈104m/s，这是不可能的。所以这里忽略流体流动带来的磁力线对流。我们模型中用到的MHD直流微型泵方程和参数都取自参考文献10。文献中采用基于有限体积法的PHOENICS代码，通过外部FORTRAN子程序连接到N-S方程中来求解Lorentz体积力，我们这儿没采用这种方法，而是使用基于有限元法的COMSOLMultiphysics3.2软件来计算。这样做的好处是不用写外部子程序，源项和体积力已经预置进去了，计算方法已经考虑到这一点并进行了优化。这样我们就能够通过Newton迭代更快收敛。以下矢量方程以无量纲形式描述了该问题：0;200tBt;xxL;yyR;0uuu;0BBB,(1)200ppuBR;00JJuB;2220BRHa;00uBR这里3(m)是流体密度，(siemens/m)是电导率，0u(m/s)和0B(Tesla)分别是速度和电磁场强度。Ha是Hartmann数，2Ha可以看作磁场与普通粘性力的比值，(Pa)p是压力，R是通道宽度，(Kg/m·s)是动力学粘度。根据以上假设和无量纲数，可以构建稳态条件下的MHD方程。连续性方程为：0u(2)这里()um/s是满足无滑移边界条件的自由发散速度场。作用在带电颗粒上的静电磁场力为Lorentz力，在不可压缩Navier-Stokes方程中表现为体积力222Re1pHaHauuuJB(3)这里J(Ampere/m2)是根据欧姆定律算出的电流密度。JuB(4)这里(Volts)是由方程(5)给出的电势向量。方程(4)右侧第一项是施加的外部电场，第二项是由于磁场感应产生的电场。方程(3)中的叉积JB是Lorentz力。电势由以下方程得到：2***uB(5)外部电场需要满足绝缘边界条件0n。电场和磁场由Maxwell方程控制。方程(4)对z积分后得到平均电流密度/20000/21()hyzhdzEBuhjj(6)通过对方程(6)的仔细分析我们能够明白当平均电流密度为零、负数或正数时的含义。假设我们在微通道中横向插入一对电极（平行与B）并将其短接，然后感应电场E0＝0，当B0为正值时00yj。现在我们在外部电路中加入一些电阻，并逐渐增大其阻值，这些平均电流将会沿y的负方向流动，随着电路电阻的增大而减小。当电阻足够大时（无限大电阻），外部电路相当于开路，也就是说00j。当00j时平均Lorentz力（00jB）起着“制动”的作用。当电流最大时制动作用也最大，此时外部电阻为零。对于有限外部电阻情况，可以从流体中获得电功，此时系统相当于电池发电机。一定要注意到当j0＝0时平均流动速度由E0=u0B0精确确定，所以通过测量感应电场就可以确定流动通量。我们假设以电源来代替外部电路中的电阻。这时电流沿y轴正向流过流体，Lorentz力沿着流体流动的方向，所以外部电源加速了流动。这就是电池泵的原理。在COMSOLMultiphysics中求解方程(2)-(5)可能需要一些技巧，最基本的是要确认使用的版本中有这些模块，动量守恒我们选用不可压缩Navier-Stokes模块，电势选用直流导电模块。3．几何形状我们的MHD系统位于一个微型反应器的矩形横截面反应段中。宽度和高度分别为300和100微米，长度约为宽度的10倍。根据轴线局部电磁场的改变将该域分为三部分（左、中、右）。磁场在中部保持常数不变，在左侧和右侧成指数递减。在划分网格时注意到大的横纵比可能会给内部边界流量守恒计算带来问题，所以需要通过调整空间坐标使得长度满足网格划分要求。流动结构和电场、磁场如图1所示。在x-y平面放置一对磁铁（它可以作为直流微型泵的电磁场），产生沿z轴正方向的磁场，流体沿x轴正方向流动，当它穿过中间区域时流动情况被磁场改变。感应电流和电场沿y轴方向。图1给出了起泵或阀作用时的内部流型，这取决于电极电势的正负号。在来流区域流体是充分发展的层流（图2左侧）。为方便起见我们假设流道高度比宽度或长度小很多，所以磁场方向的电流密度比其它方向要小很多。所以该问题可以简化为二维模型。图1微通道几何结构图2(a)速度分布和电势等高线；(b)流线；(c)Lorentz力矢量。N=5×103；Ha＝4；Re＝3.2×10-33.1边界条件表1边界条件入口2*0(1,)1yuyuR*(1,)0vy壁面**(,)(,)0uxRvxR**(1,)(1,)0RR电极*0(,)2axRVx(7)绝缘壁面*(,)02axRxn轴对称*(,0)0uxy*(,0)0x出口*(1,)0py4．在COMSOLMultiphysics3.2中实现打开COMSOLMultiphysics3.2，进入模型导航栏。表2给出了详细的模型建立步骤。表2MHD微装置模型模型导航栏从下拉菜单中选择2-D选择FluidDynamics|IncompressibleNavier-Stokes点击Add，点击Multiphysics按钮选择Electromagnetics和ConductivemediaDC完成Draw菜单从(-1,0)到(-0.11,0.5)绘制矩形从(-0.11,0)到(0.11,0.5)绘制矩形从(0.11,0)到(1,0.5)绘制矩形从Options菜单中选择Axes/GridSettings设定xmin=-1.2;xmax=1.2;ymin=-0.5;ymax=1Options菜单|ConstantsRe=2e-3;Ha=10;u0=1;B0=1;a=1/9;N=500;V0=-10Options菜单|ScalarExpressions|GlobalBL=B0*exp(9*(x+a))BR=B0*exp(9*(a-x))Phyiscs菜单|Boundarysettingsns模式入口，设定边界条件：Inflow/Outflowvelocity设定u0=1-(y/0.5)^2和v0=0轴线边界条件：Slip/Symmetry顶部壁面边界条件：Noslip电极边界条件：Noslip出口边界条件：Outflow/PressureDC模式入口边界条件：V＝0轴线边界条件：V＝0顶部壁面边界条件：ElectricInsulation点击边界条件：ElectricPotential=V0Phyiscs菜单|SubdomiansettingsNs模式。子域11/;1/^2;*^2*;0xyyNHaFBLVBLuF子域21/;1/^2;0*0^2*;0xyyNHaFBVBuF子域31/;1/^2;*^2*;0xyyNHaFBRVBRuFDC模式。子域11.0;9.0****xyQvBLBLvBLu子域21.0;xyQvu子域31.0;9.0****xyQvBRBRvBRuSolve菜单点击工具栏上的求解按钮（＝）5．计算结果交互作用变量（Re/Ha2）和Hartmann数的变化源自Reynolds数和磁场强度的变化。对泵和阀的几种情况进行了建模。在所有情况下都考虑了磁铁区域外部的非导电壁面。通过在流道两侧放置磁铁，引入了互相垂直的电场和磁场，为了说明交叉的电磁场是如何起到泵和制动作用的，我们绘制了每种情况下的速度和压力分布。5.1MHD直流微型泵图2给出了MHD微流道的流型和电势等高线(a)，流线(b)和Lorentz矢量(c)。由于上游和下游区域没有电极，壁面绝缘，电动势和y方向的静电场方向相反。电动势决定了那里磁场比较强，该区域以外磁通密度迅速降低且静电力起主要作用。所以系统表现为制动作用。这从图2中可以看出，当从上游进入中间区域时，速度分布变平了，因为流动进入恒定磁场区域（中心）。速度分布呈现M形，因为Lorentz力在靠近电极的地方更强，使得这里比中心的轴线速度阻力更小。这使得流体进入电极区域以后速度分布发生明显的抛物线型扭曲。在这些区域流动加速形成凹型分布。下游的最大速度向中心移动，使得M型分布更为明显。由于流向下游时电动力降低，使得最大速度分布向中心移动，直到流动再次恢复抛物线分布。电动势的等高线从电极最大值处向外扩散，沿轴线为零。图2(b)给出了相应的流线。由于Lorentz力将流体向壁面处挤压，流线在这里比较集中，然后流动加速以保证连续性。图2(c)给出了由交叉电流密度J和磁场强度B产生的Lorentz力矢量。磁场是外部施加的，由永磁体或电磁体产生。在这种情况下磁雷诺数很小，所以采用弱耦合，认为磁场不受流体流动影响。电流密度矢量由两部分组成，一部分是由电极静电场产生，另一部分是由电动势场U×B产生。在上游和下游区域中静电场力占主导地位，电场矢量指向电极（图3(b)），所以上游轴线电流沿x轴正方向，下