最新内容

模拟电渗流和双电层
微流控设备非常小,因此控制和混合设备内液体的微泵和微混合器不能包含任何移动部件。作为替代,它们必须利用电渗流工作。这篇博客,我将介绍电渗流和双电层的概念,以及如何在 COMSOL 中建立这些模型,并演示 2 个实例模型。 微流控设备需要无移动部件的泵和混合器 近年来,微流控片上实验室系统在将传统的实验室规模的化学和生物分析的尺寸缩小到毫米到几厘米大小的片上形式上发挥了重要作用。这些设备通常被称为微型全分析系统(μTAS),并在医疗诊断、药物测试和输送、法医分析、DNA 分析,甚至免疫分析和毒性监测等方面都有着深度应用。由于尺寸极小,这些设备具有如即时检测(POCT)和诊断等众多优势。微流控设备需要的液体体积更小,适用于没有大量样品或试剂昂贵的应用。这些设备还可以一次处理多个样本(称为并行处理),并且对功率的要求低。 在典型的片上实验室系统中,片上混合和流体控制非常重要,通常这些系统需要使用微泵来控制通道内的流体流动,以及微混合器来加速混合过程。这些芯片中的流体通道的尺寸通常在 1µm 到 500µm 之间变化。在这种长度尺度上,使用任何移动部件来建造泵和混合器都是不实际的。没有移动部件的设备也更加可靠。那么,如何在没有任何移动部件的情况下激活流动?答案是:电渗。 什么是“电渗”? 在微流控领域,流动通常由电场驱动。顾名思义,电渗 指由外加电势引起的液体在微通道上的运动。用电场驱动流动,可以制造出没有移动部件的 泵和混合器。 是什么在驱动电渗流动? 为了更好地理解电场是如何驱动流动的,我们首先要了解在微通道内非常靠近壁的地方,即流-固界面处发生了什么。大多数的片上实验室设备是由硅玻璃制成的。当与流体(可以是水或任何缓冲溶液)接触时,玻璃表面参与酸碱反应和离子交换——我们把这个复杂的过程称为表面化学。由于表面化学作用,玻璃表面获得负密度电荷。为了对这种表面化学进行连续描述,引入了双电层的概念。双电层反映了流-固界面处的不均匀的电荷(离子)分布,且在物体表面分为两层: 第一层,被称为表面电荷层,是由于化学反应而被表面吸收的离子(在这种情况下是负电荷)组成的。 第二层,被称为扩散层,由受电吸引力和热运动的影响而被吸引到表面的自由离子组成。第二层电荷屏蔽了表面电荷,其净电荷与表面电荷相等,但极性相反。 双电层。 上图总结了双电层 的结构,显示了离子的分布是到玻璃壁面距离的函数,以及双电层中的电位(上面的蓝线)与电中性体中的一个点的关系。如果我们仔细观察扩散层,可以注意到它可以进一步被分成由滑动平面隔开的两部分。这个平面将左边的不动的流体(附着在表面上)与在切向应力影响下自由运动的流体分开。然后,由于库仑力的作用,可以使用电场来诱导双电层中净电荷的运动。离壁更远的地方是第三层,即电中性体。 由于在微流体通道中很难在不破坏流动的情况下进行任何形式的测量,所以通常从计算的角度对这些芯片进行分析。那么,如何使用仿真软件对其进行建模呢? 双电层模拟 这个问题包含三个物理场: 静电 物理场包含方程、边界条件和空间电荷,用于求解电势问题。电场是由电势场(E=-\nabla V)的梯度恢复的。空间电荷是由阴阳离子的贡献加起来得到的。这些离子的浓度是通过稀物质传递 接口计算出来的。 稀物质传递 接口求解了混合物中稀释物质的质量输送问题,求解了物质浓度问题。化学物质是通过扩散(菲克定律)、迁移(当与电场耦合时。在这种情况下,电场由静电 接口计算)和对流(当与流体流动耦合时。在这里,流体流动是由层流接口 计算的)模拟的。 层流接口有方程、边界条件和体积力,使用的是纳维-斯托克斯方程边界条件模拟自由运动的流体,求解速度场和压力问题。体积力 \rho{e} E,其中 \rho{e} 是电荷密度,由静电 接口计算。 双电层的厚度一般在几纳米左右,离子浓度在壁附近呈指数级变化。由于双电层的厚度非常小,在这个区域使用近似值可能是有利的。COMSOL 包含一个电渗速度边界条件,它忽略了壁面和滑动平面之间的流场,并根据 Zeta 电位用 Helmholtz-Smoluchowski 关系分析计算壁面的速度: \textbf{u}=\frac{\epsilon{W}\zeta{0}}{\eta}\nabla_{T}V 由此创建的模型的计算要求将大大降低。对许多实际工程应用来说,这是一种非常有用的方法。 因此,我们建议您在使用 COMSOL Multiphysics 建立模型之前,首先计算德拜长度,即双电层的长度。如果这个长度比要模拟的几何长度尺度小得多,就使用电渗速度边界条件。如果不是,则使用传统的无滑移速度壁边界条件,并求解双电层中的流动问题。请记住,浓度随双电层的电位呈指数变化;如果不使用电渗速度边界条件,就必须使用细的边界层网格来求解双电层的突然变化。 在下一节中,我将向大家展示两个适用于此方法的例子。第一个例子是一个微型泵,它在整个几何中被解析,包括双电层。第二个例子,是一个使用电渗速度边界条件的微混合器。 应用实例 示例1:微泵 在这个例子中,微通道(长 60nm,高 10nm)的上下壁均带负电(-0.02 C/m^{2}),在入口(左边界,6mV)和出口(右边界,0V)使用电极驱动流动。由此产生的电势、空间电荷分布和速度场如下图所示: 用静电物理场接口计算的可视化电势图。 净空间电荷的绘图,即阴离子和阳离子的贡献之和。带正电的净电荷屏蔽了带负电的壁。 速度绘图显示了在库仑力的作用下,双电层中净电荷的运动。 为了介绍下一个例子,即微混合器,让我们先看看当上部和下部玻璃壁的部分(在下一幅图中以蓝色显示)带正电 (0.06 C/m^{2}) 而不是负电时会发生什么: 电势、空间电荷和速度场的模拟结果如下所示: 虽然空间电荷在带负电的壁周围带正电,但在带正电的壁面周围是负的。空间电荷的反转导致壁附近的电渗速度相反。这种相反的近壁速度导致在通道中引入漩涡(如流线图所示),这可以用来混合不同的化学物质。 例2:微混合器 在微观层面上,流动通常是高度有序的层流,湍流的缺乏使得扩散成为混合的主要机制。虽然小分子(以及快速扩散的物质)的扩散混合可以在几十微米的距离内发生,但较大的分子,如肽、蛋白质和高分子量的核酸的混合可能在可比的距离内需要几分钟到几小时的平衡时间。对于许多化学分析来说,这种延迟是不切实际的。这些问题导致对微流控系统的更有效的混合器的强烈需求。在下面的例子中,壁是带负电的,就像我们刚看过的微泵例子一样,在入口(左边边界)和出口(右边边界)使用电极来驱动流动。为了引入一些混合,四个额外的电极被放置在混合室的壁上。这四个电极在壁上引起波动的电渗速度: 说明混合器如何运行的示意图。在入口处使用了两种不同浓度的流体来研究混合过程。 当不施加电场时,流动为层流,扩散系数非常小,所以两种流体在出口处分离良好。 当施加交变电场时,由于流动中的交变漩涡,混合大大增加。 结束语以及后续操作 这篇博客简要介绍了电渗流和双电层的概念,并演示了如何在 COMSOL 中对这类问题进行建模。如果你想了解更多关于电渗微混合器模型的信息,可以从 COMSOL 案例库中下载该模型和模型文档。您可以在此处了解 […]

线性静态问题的网格剖分注意事项
本篇博客中,我们介绍了线性静态有限元问题的网格剖分注意事项。这是网格剖分技巧系列博客的第一篇,希望能帮您建立起对有限元模型剖分网格的信心。

求解线性稳态有限元模型
本篇博客是求解器系列的首篇博客,将介绍用于求解所有线性稳态有限元问题的算法。虽然我们在博客中基于一维有限元问题进行介绍,但所讲解的内容具有普适性,能帮助您理解博客系列中接下来将介绍的更加复杂的非线性多物理场的求解技巧。

管路连接处预应力螺栓的拉伸应力仿真
在计算拉伸应力之前模拟螺栓预拉伸是一项重要的分析,使用 COMSOL Multiphysics 可轻松进行建模。

使用全局方程模拟室内空气温度
前些天我遇到了一个有趣的问题,它促发了我想要写一篇博客的灵感,因为它使我想到了我们的 COMSOL Multiphysics 软件中一个非常强大但常常没有得到充分利用的功能:全局方程。

模拟涡轮静叶片的热应力
仿真软件能够帮助我们理解和优化组件设计。任何一个仿真都需要基于真实情景建立 模型,以足够详细地表征实际应用,从而获得特定应用或组件的相关信息。

为什么微波炉加热食物不均匀?
我们可能都经历过这样的场景:下班回到家,把昨晚的剩菜放在微波炉里,坐下来准备吃一顿简餐,结果却发现吃到的食物一口滚烫,一口冰冷。这样的经历不止一次促使我思考:为什么微波炉对食物的加热会这么不均匀?

使用 COMSOL Multiphysics® 模拟磁致伸缩效应
如果你曾经站在变压器旁边,可能听到过它发出的嗡嗡声,并怀疑附近是不是有蜜蜂。下次再听到这种声音时,你大可以放心,这不是蜜蜂,而是变压器铁芯的磁致伸缩发出嗡嗡声。 什么是磁致伸缩? 磁致伸缩是一种效应,它会使所有暴露在磁场中的磁性材料的形状发生变化。例如,磁致伸缩效应会使一块铁伸长 0.002%,使镍收缩 0.007%。这一现象曾经因为被用在第一次世界大战期间的声呐设计中而引起了广泛关注。进一步的研究,最终研制出了用于工程的磁致伸缩材料,例如 Terfenol-D,以及最近研制出的 Galfenol,它的伸长率高达 0.04% ~ 0.2%。 磁场引起的应变现象也称为正(磁致伸缩)效应。磁致伸缩效应可以追溯到原子级的相互作用,它是磁性材料中的磁能和机械能在受到磁场和机械应力时所发挥的平衡作用而产生的。下面的动画是对磁致伸缩材料内部情况的简单说明。 当对材料施加交变的磁场时,构成材料的微小椭圆体磁铁会随着磁场大小和方向的变化来回翻转。这些微小磁性体方向的改变表现为一种宏观应变。如果以典型的电力线频率(50Hz – 60Hz)交变磁场,材料中的交变应变会使它像扬声器一样工作,从而产生可以听见的声音。这就解释了变压器发出的嗡嗡声之谜。 这种双向磁机械耦合也会产生逆 效应,即作用在磁性材料上的应力可以通过调整这些微小磁体的方向来改变材料本身的磁性状态。正效应和逆向效应分别用于驱动类和传感类应用。 磁致伸缩材料的应用 从航空航天、石油生产到声学和 MEMS,磁致伸缩材料几乎可以应用在所有行业。下面列出了一些重要的商业应用: 声学设备 声呐 水听器 用于清洗、混合和乳化的超声波振动器 超声波摩擦焊接 驱动器 直线电机和旋转电机 尺蠖式驱动器 用于机床头部的位置控制器 燃油喷射系统 光学扫描系统 液压驱动器,例如伺服阀和泵 用于减小阻力的智能机翼中的主动后缘 传感器 位置传感器 非接触式扭矩传感器 磁场传感器 MEMS生物和化学传感器 振动控制 减振器 平台稳定器 图像稳定器 能量收集器 混合智能结构 带混合压电/磁致伸缩磁芯的 Tonpilz 换能器 混合压电/磁致伸缩复合驱动器和传感器 您还可以利用磁致伸缩效应把家里客厅的墙壁或窗户变成扬声器! 那么,如何在 COMSOL Multiphysics 中对这种有趣的现象进行建模呢? 在 COMSOL Multiphysics 中对磁致伸缩进行建模 对磁致伸缩型换能器进行建模的正确方法包括准确模拟磁和结构性能,并使用适当的材料模型模拟这些物理场之间的相互作用。COMSOL 中内置了预定义的物理场接口,可用于设置磁仿真和结构仿真。COMSOL 还支持灵活地设置用户自定义的本构关系,用数学的方式表示材料模型。 实验表明,正向和逆向磁致伸缩效应都是非线性的。当模拟那些在准静态条件下运行,但暴露在大范围机械力和磁场中设备时,建立完整的非线性响应模型可能很重要。在这类设备中,了解磁致伸缩磁芯在什么工作条件下饱和是有用的。这些信息可以为设计人员提供极限值,还可以解释实际的非线性行为,例如传感器灵敏度的变化或用户期望从磁致伸缩设备获得的驱动器最大力。 在某些已知频率和已知工作条件下工作的声学换能器中,可以使用线性本构定律简化材料模型。这些定律(或方程)是在假设换能器操作涉及围绕偏置点的小幅振荡的条件下推导的。而如果在建模方法中考虑这些实际因素,那么我们就能够轻松模拟磁致伸缩换能器在较宽工作频率范围内的响应。 在 COMSOL Multiphysics 中,可以同时设置非线性和线性本构方程模拟磁致伸缩器件。接下来,我想与大家分享一些我们对一个实验换能器进行模拟的结果。 模拟磁致伸缩换能器 典型的换能器有一个被驱动线圈包围的磁致伸缩磁芯。流过线圈的电流会产生磁场。传感器有一个钢制外壳,包围着驱动线圈和铁芯。磁芯连接到活塞上,活塞用于在启动器配置中将磁芯的位移传递到外部机械部件上,或在传感器配置中将负载从外部机械或声源传递到磁芯上。钢制外壳、活塞和铁芯形成了一个封闭的磁通路径。 对于非线性模型,我们使用了 Galfenol的典型材料表征曲线,并确定了重要设计参数的非线性,例如换能器的阻力。我们还能够探索驱动和传感行为的变化与各种磁场和作用在传感器上的拉伸和压缩载荷的函数关系。有关这个模型的更多信息,请查看 COMSOL 案例库中的非线性磁致伸缩换能器和传感器 教程案例。 非线性磁致伸缩换能器仿真中的位移幅度、驱动器和传感器曲线以及换能器阻力图。 对于线性模型,我们使用了 Terfenol-D 的典型材料参数,并生成了驱动器载荷线。我们还研究了换能器位移的幅度和相位,以及驱动线圈阻抗的频率响应。 线性磁致伸缩传感器仿真中的驱动器载荷线、线圈阻抗、位移幅度和位移相位图。 2013 年 COMSOL […]