第11章 电动力学流体 WBJ ZIMMERMAN and JM MACINNES Department .DOC

第11章 电动力学流体 W.B.J. ZIMMERMAN and J.M. MACINNES Department of Chemical and Process Engineering, University of Sheffield, Newcastle Street, Sheffield S1 3JD United Kingdom E-mail: w.zimmerman@shef.ac.uk 本章探讨微流道网络中电动力学流体的多物理场模拟。对于在Y-连接处合并的长流道网络的模拟，它把长流道（0-D）耦合边界条件的代数约束和2-D Y-连接结合起来。采用耦合变量，这些作为多物理场的扩展来实施。 1．简介 本章的目的是演示COMSOL Multiphysics建立电动力学流体模型的功能。电动力学流动最先进的应用是微流体，其中，微量的化学物（纳升）被瞬时输运，达到微流道网络中复杂的交换和排序，通过严格控制以获得化学反应和组分变化的高产率。由物理化学现象得到的运动流体特别重要，因为它涉及到快的响应时间且运动的机械部件不会被损坏。微流体和微机械器件（MEMs）之间存在很大的重叠。例如，与柔性致动器壁面和流道侧相邻的运动大分子被看成是微流体，但它们也是分子机械，但其尺度太小不能被看作传统的运动部件。 为了建立尽管是最简单的电动力学模型，多物理场都是很重要的，它耦合了电势场，化学输运和动量输运（Navier-Stokes）。然而，第一个简单的方法就引入了边界条件的耦合来逼近电化学边界层的运动。尽管这样的耦合是线性的，但COMSOL Multiphysics习惯于得到能接受的模拟，它要求弱的边界约束。然而，COMSOL Multiphysics处理应用节点时耦合了线性的Dirichlet边界条件和场变量，但不依赖于用户所认定的弱约束。在第二章已经讨论过，使用拉格朗日乘法器技术，COMSOL Multiphysics对于一般边界和辅助边界很容易处理。然而，Femlab的早期版本采用一个进入Model COMSOL Multiphysics程序中特定的弱边界节点的方法，然后要求用户对所应用的节点所耦合的任何边界进行确认。对已建立的模型所产生的m脚本文件进行仔细的检查发现，COMSOL Multiphysics对速度场和电场线性耦合的弱边界约束是否需要进行自动检测，并实施必要的拉格朗日乘法器。 当用户希望使用衍生2.1. 背景 电动力学流体由电场和液体中荷电组分（离子）相互作用产生。它存在两种明显的作用：与具有净电荷的壁面相邻的双层区域中液体所受的电场力，和流体体相中单个离子的运动（在双层区域之外），该体相中一般没有净电荷。对于所感兴趣的流道尺寸（比如说大于1(m），该双层可被看作为微元，它对流动的影响相当于（MacInnes et al.[3]，MacInnes [2]）边界条件对速度ui，电场和相关化学组分的质量分数Y的应用。 ， (1) 式中ni为垂直于壁面的单位矢量。 需要求解的方程组包括动量方程，质量守恒方程，电荷守恒方程和组分方程。最简单的例子可以由无量纲变量来表示。 动量输运和连续方程（Navier-Stokes）： ， （2） . 包括电脉的组分输运方程： . （3） 电荷守恒方程： . （4） 满足式（4）的电场一定也满足Gauss定律（参看第7章公式（1）），它是决定流体中电荷密度与位置的关系的方程。在电动力学流体的典型条件中，为了满足电荷守恒方程（4）和动量平衡方程（2），电荷密度可以忽略不计。电导率和zeta电位取决于组分Y的浓度，这里假设它们存在线形关系：和，式中下标r代表所考虑的流体中所涉及的两种纯溶液属性的比值。 作为流动入口边界条件的电势、压力和组分浓度必须给定，流动出口的电势和压力也必须确定。出口边界的组分浓度未知，组分扩散是唯一将区域中的组分浓度场和出口边界的组分分布联系起来的项，所考虑的组分扩散要给出一个近似值。一般认为，出口边界的组分扩散可以忽略不计，例如采用存在于流率项中的扩散项的Neumann边界条件。 电场被看作为准稳态，也就是说电场随着速度和浓度的变化及时进行调整。以上的方程代表了一个普适的问题，它提供了一种数值实现方法，当该方法一旦被确定，它可以用来计算任何特定电动流体的问题。对于例子的实施，假设系数值为1/Re=30，1/Pe=0.03，（壁中zeta电位不变），z=0（组分Y没有电