流体 & 传热 博客文章
通过 CFD 仿真分析帆船航行的物理原理
在所有类型的船只中,我发现帆船是最迷人的一种,尤其是当帆是唯一的推进方式时。动能从风传递到帆,产生能够平衡阻力并使帆船在水中移动的升力。
石墨烯革命:第四部分
石墨烯可以在高真空下利用热分解制造。为了设计和优化这些高真空系统,工程师可能转向仿真寻求帮助,不过目前还没有多少仿真工具能胜任这一工作。让我们来看一下真空系统如何与石墨烯的生产相关,为什么要模拟它们,以及如何模拟。
什么是分子流?
提到真空环境,大家可能很自然地会想到外太空的恶劣环境。在地球上建立人工真空环境是一项非常具有挑战性的任务。事实证明,我们甚至无法实现星际真空的低压,分子流就是发生在这种低压环境下。
非牛顿流体:倒番茄酱的困惑
如果你喜欢吃番茄酱,那么很有可能你已经体验过我们称之为“要么全有要么全无的番茄酱困境”。想象一下,你伸手去拿一个新的番茄酱玻璃瓶,取下瓶盖,把瓶子倒过来时,然后——“什么都没有发生”。直觉上,我们会摇晃或轻拍瓶子,然后——突然你的食物就完全被涂上了番茄酱(当然,除非你的反应真的很好)。在这篇博客文章中,我们将通过仿真模拟瓶子中的番茄酱的流量来演示倾倒番茄酱时的“全有或全无”效果。 我的番茄酱在哪里? 番茄酱除了是一种特别受欢迎的调味品以外,也是一种非牛顿流体。这是它会出现上述奇怪现象的关键。作为一种非牛顿流体,番茄酱的黏度——以及其流动能力——会随着流体的速度梯度或剪切速率变化。当第一次把番茄酱瓶子倒过来的时候,由于番茄酱的初始黏度很高,很可能你只能得到非常小的流量,如果有的话。好在番茄酱还有一个有趣的剪切稀化特性。当在重力之外施加外力时,增加的剪切力将导致番茄酱的黏度降低,更容易流动。这就是为什么我们有时必须轻拍或摇晃瓶子才能享受番茄酱加薯条。 模拟瓶子中的非牛顿流体流动 为了满足我们的好奇心,以及我们的胃口,并进一步完善我们的番茄酱倾倒技术,我们使用COMSOL Multiphysics 软件和其附加的计算流体力学模块对这个过程进行模拟。在模拟中,我们做了几个简化的假设,包括瓶子的无限延伸(即番茄酱一直都有)和单相层流,忽略了瓶内滞留的空气。所有速度均以米/秒(m/s)为单位。 我们规定了瓶子的几何形状和层流的边界,包括无滑移壁面、压力入口(代表未建模番茄酱部分的重量)和压力出口(外部大气),如下图所示。 用于模拟的二维几何图形。显示了应用的边界条件。 我们还添加了体积力特性,以包括重力和沿着 y 轴摇晃瓶子 0.2 秒所产生的力。 为了给求解器提供一致的初始条件,重力产生的体积力通过阶跃函数从 0 上升到其最大值(过渡区域为 0.05 秒)。摇动动作是以 0.2 秒宽的矩形函数模拟的。在此期间,体积力加倍。 倾倒番茄酱 我们基于Koocheki 等人的工作,我们先做了一个仿真为之后模拟番茄酱的流动做准备。为了模拟番茄酱的流动,我们使用了非牛顿幂律来表示动力黏度。番茄酱配方通过在幂律方程中指定稠度系数和流动特性指数来模拟。下面的图像显示了添加1%黄原胶的番茄酱配方中较高的初始黏度和非牛顿流动性。 左:以p1(见上面的模型设置)为单位测量的黏度是时间的函数,时间间隔 [0,0.1]秒。右图:以p1(见上面的模型设置)为单位测量的黏度是时间的函数,时间间隔为[0.1,2]秒。番茄酱的剪切稀化特性很明显,尤其是在瓶子摇晃的[1,1.2]秒间隔内。黏度降低,从而允许更多的番茄酱流出瓶子。 左:在三个不同时间沿线 1(见上面的模型设置)测量的黏度。右图:平均出口速度与时间的函数关系。可以看出,番茄酱在重力的作用下慢慢开始流动,然后在摇晃瓶子时容易增加。 左:速度幅度动画从 0.05 到 2 秒。右:1.1 秒时的速度大小。 “番茄酱困境” 等待流速的变化对你来说可能很熟悉。一旦流速开始改变,就会快速变化,导致番茄酱泛滥,覆盖在你的食物上。流速顺其自然地增加导致了剪切力增加,进而降低黏度,导致更高的流速。这就是为什么番茄酱最初似乎是静止的,但后来突然又像水一样流动的原因。 番茄酱从瓶子中流出这个特殊的问题比一开始看起来要复杂得多。有许多变量会影响我们体验一个“绝对的” 番茄酱惊喜。如果你有兴趣进一步研究这种效应,建议你使用 COMSOL Multiphysics 和计算流体力学模块进行模拟实验。例如,你可以研究诸如倾斜角度、瓶子里的空气、甚至瓶子的设计等变量是如何影响番茄酱的倾倒体验的。
石墨烯革命:第二部分
我在之前的一篇博客中曾经介绍过石墨烯的一些奇异特性。石墨烯仅包含单层原子,这意味着任何石墨烯基结构都将有极高的纵横比;而高纵横比的几何也为模拟带来了各种挑战。 石墨烯的传热模拟 COMSOL 提供了多种工具来帮助模拟具有极高纵横比的几何及特征。最近,人们利用 COMSOL 对石墨烯“被子”进行了传热模拟,《Nature Communications》杂志的”用于高功率 GaN 晶体管热管理的石墨烯被子“一文介绍了该研究。论文作者使用 COMSOL Multiphysics 证实,可通过引入由薄层石墨烯 (FLG) 制作的额外散热通道,即顶面导热片,来显著提升 AlGaN/GaN 异质结场效应晶体管 (HFET) 的局部热管理。 COMSOL Multiphysics 的传热接口支持您使用薄层特征特征模拟极高纵横比的组件。这一特征仅求解表面切面处的传热方程,因此免去了在高纵横比层中使用极端细化网络的需求。使用此方法能够极大减少计算时间和内存使用。 薄层特征设定窗口。 石墨烯的电气模拟 从 2006 年开始,人们就已经使用 COMSOL 来研究石墨烯的电气特征。在这篇论文中,研究人员使用 COMSOL 来推导石墨烯基复合材料的面内和横向电导率。我们可以很轻松地在 COMSOL Multiphysics 中输入电导率的张量物理量。您仅需提供电导率张量元,它可以是温度或其他任何量的函数。 可以轻松在电流接口电流守恒特征的设定窗口中指定各向异性电导率。 石墨烯的结构力学模拟 COMSOL 也可以模拟石墨烯的结构力学应用。在这篇论文中,研究人员计算了石墨烯膜在压力差作用下产生的挠度和应变。可通过电气检测到带结构中的变化,这说明它可用于制造超灵敏压力传感器。结构力学模块的壳接口主要用于薄壁结构中的结构力学分析,因此非常适合此类应用。壳接口使用 Mindlin-Reissner 类公式,即考虑了横向剪切形变。这意味着我们无须对极薄结构进行网格剖分,就可以获得高精度的结果。 壳接口中材料模型的设定窗口。 相关示例模型 我们现在已经分析了热学、电气以及力学的模拟概念,您可能在想有什么模型能同时用到所有这些概念。确实有一个,您可以查看案例库中电路板加热模型,如下所示。 电路板加热模型演示了热学、电气以及力学模拟概念。 这一多物理场示例模拟了加热电路器件时电热的产生、传热,以及机械应力和变形。模型用到了传热模块的固体传热接口,AC/DC 模块的电流,壳接口,以及结构力学模块的固体力学和壳接口。
模拟汽车制动盘的加热过程
汽车需要制动器(俗称刹车)的原因显而易见:你一定不希望制动器发生故障。造成制动器故障的原因有很多,其中一个就是制动盘过热。
模拟锂离子电池的散热
阅读全文,简要了解如何在 COMSOL Multiphysics® 软件中通过三个顺序研究来创建锂离子电池模型。在文章结尾,我们为您附上了教学模型的下载链接。