在 COMSOL® 中构建磁流体动力学多物理场模型

软件中的模型都是从零开始构建的,软件支持多物理场,因此用户可以按照自己的意愿轻松地组合代表不同物理场现象的模型。有时这可以通过使用软件的内置功能来实现,但有些情况下,用户需要做一些额外的工作。我们以构建磁流体动力学(MHD)模型为例介绍一下这个工作流程。

磁流体动力学的多物理场建模

MHD 现象的建模本质上是一个多物理场问题;必须用数值方法求解流体流动、电流和磁场之间的耦合问题。这些不同的场都是由偏微分方程描述的,可以通过有限元方法求解。
在 COMSOL® 中构建磁流体动力学多物理场模型的图1
施加电流时两个磁体之间通道中导电流体的 MHD 问题。

我们看看如何在一个相对简单的问题背景下进行建模:如上所述,绝缘的矩形通道内为不可压缩导电流体,这个通道连接两个流体静压相等的无限大容器(未建模)。有两个电极穿过流动通道在两侧伸出,通过施加电势差驱动电流通过流体。此外,在上方和下方分别放置一个圆形磁铁。磁体产生静磁场在 COMSOL® 中构建磁流体动力学多物理场模型的图2,使得具有导电性 在 COMSOL® 中构建磁流体动力学多物理场模型的图3以一定速度 在 COMSOL® 中构建磁流体动力学多物理场模型的图4移动通过该磁场,从而产生感应电流。在 COMSOL® 中构建磁流体动力学多物理场模型的图5。除了这些感应电流之外,由于电势场的边界条件,还会产生电流 在 COMSOL® 中构建磁流体动力学多物理场模型的图6,因此流体中的总电流变为:

在 COMSOL® 中构建磁流体动力学多物理场模型的图7

流经磁场的电流将对流体产生体积力 在 COMSOL® 中构建磁流体动力学多物理场模型的图8,并将流体从一个容器泵送到另一个容器。我们假设系统在稳定状态下运行。

耦合电场、磁场和流场

对于这个问题,我们需要求解流体中的偏微分方程组来描述电场和磁场。方程式为:

在 COMSOL® 中构建磁流体动力学多物理场模型的图9

在 COMSOL® 中构建磁流体动力学多物理场模型的图10

这组方程通过磁场和电场 接口(AC/DC模块的一部),使用安培定律和电流守恒 特征以及单独的速度(洛伦兹项) 特征求解。

在移动流体周围的空间中,没有电流,所以我们只需求解单矢量方程:

在 COMSOL® 中构建磁流体动力学多物理场模型的图11

其中在 COMSOL® 中构建磁流体动力学多物理场模型的图12是剩余磁通密度,它仅在磁域中非零。当单独求解上述方程时,请使用磁场和电场 接口中的安培定律 特征。

我们假设通道壁的属性不影响场,因此在模型中忽略它们。使用一组材料属性和边界条件来给出说明性结果。任何位置的磁场边界条件都是磁绝缘 条件, xy 平面除外,该平面采用理想磁导体 条件来利用系统的对称性。表示电极的域必须一直延伸到建模域的边界,接触磁绝缘 边界,以提供电流返回路径电压 接地 终端 条件应用于这些外表面,而电绝缘 条件应用于所有其他适用的边界。

此外,我们还需要求解通道中的流场。我们假设流动是层流,从而在通道域中求解纳维-斯托克斯方程。如果流动是湍流,我们可以添加一个湍流模型开放边界 条件应用于通道的两端,表压为零。对称 条件应用于 xy 平面。计算域如下图所示。

软件中的模型都是从零开始构建的,软件支持多物理场,因此用户可以按照自己的意愿轻松地组合代表不同物理场现象的模型。有时这可以通过使用软件的内置功能来实现,但有些情况下,用户需要做一些额外的工作。我们以构建磁流体动力学(MHD)模型为例介绍一下这个工作流程。

磁流体动力学的多物理场建模

MHD 现象的建模本质上是一个多物理场问题;必须用数值方法求解流体流动、电流和磁场之间的耦合问题。这些不同的场都是由偏微分方程描述的,可以通过有限元方法求解。
在 COMSOL® 中构建磁流体动力学多物理场模型的图13
施加电流时两个磁体之间通道中导电流体的 MHD 问题。

我们看看如何在一个相对简单的问题背景下进行建模:如上所述,绝缘的矩形通道内为不可压缩导电流体,这个通道连接两个流体静压相等的无限大容器(未建模)。有两个电极穿过流动通道在两侧伸出,通过施加电势差驱动电流通过流体。此外,在上方和下方分别放置一个圆形磁铁。磁体产生静磁场在 COMSOL® 中构建磁流体动力学多物理场模型的图14,使得具有导电性 在 COMSOL® 中构建磁流体动力学多物理场模型的图15以一定速度 在 COMSOL® 中构建磁流体动力学多物理场模型的图16移动通过该磁场,从而产生感应电流。在 COMSOL® 中构建磁流体动力学多物理场模型的图17。除了这些感应电流之外,由于电势场的边界条件,还会产生电流 在 COMSOL® 中构建磁流体动力学多物理场模型的图18,因此流体中的总电流变为:

在 COMSOL® 中构建磁流体动力学多物理场模型的图19

流经磁场的电流将对流体产生体积力 在 COMSOL® 中构建磁流体动力学多物理场模型的图20,并将流体从一个容器泵送到另一个容器。我们假设系统在稳定状态下运行。

耦合电场、磁场和流场

对于这个问题,我们需要求解流体中的偏微分方程组来描述电场和磁场。方程式为:

在 COMSOL® 中构建磁流体动力学多物理场模型的图21

在 COMSOL® 中构建磁流体动力学多物理场模型的图22

这组方程通过磁场和电场 接口(AC/DC模块的一部),使用安培定律和电流守恒 特征以及单独的速度(洛伦兹项) 特征求解。

在移动流体周围的空间中,没有电流,所以我们只需求解单矢量方程:

在 COMSOL® 中构建磁流体动力学多物理场模型的图23

其中在 COMSOL® 中构建磁流体动力学多物理场模型的图24是剩余磁通密度,它仅在磁域中非零。当单独求解上述方程时,请使用磁场和电场 接口中的安培定律 特征。

我们假设通道壁的属性不影响场,因此在模型中忽略它们。使用一组材料属性和边界条件来给出说明性结果。任何位置的磁场边界条件都是磁绝缘 条件, xy 平面除外,该平面采用理想磁导体 条件来利用系统的对称性。表示电极的域必须一直延伸到建模域的边界,接触磁绝缘 边界,以提供电流返回路径电压 接地 终端 条件应用于这些外表面,而电绝缘 条件应用于所有其他适用的边界。

此外,我们还需要求解通道中的流场。我们假设流动是层流,从而在通道域中求解纳维-斯托克斯方程。如果流动是湍流,我们可以添加一个湍流模型开放边界 条件应用于通道的两端,表压为零。对称 条件应用于 xy 平面。计算域如下图所示。


在 COMSOL® 中构建磁流体动力学多物理场模型的图25在 COMSOL® 中构建磁流体动力学多物理场模型的图26


计算域和边界条件。

流动将由流体中电流和磁场的相互作用产生的体积力 在 COMSOL® 中构建磁流体动力学多物理场模型的图27。这个力的表达式没有内置到软件中,所以在这里我们需要做一些手工操作。我们需要找到电流和磁场分量的内置表达式,可以通过查看方程视图 并生成报告来实现,如知识库条目中关于实现用户定义的多物理场耦合的描述。这些内置表达式用于定义流体上的体积力,如下面的屏幕截图所示。
在 COMSOL® 中构建磁流体动力学多物理场模型的图28
显示计算力分量的变量的屏幕截图。

最后,要将计算出的速度场耦合回电磁问题,请使用磁场和电场 接口中的速度(洛伦兹项) 特征,如下面的屏幕截图所示。请注意,软件会自动将流体速度场识别为此特征的输入。非常简单!这两个物理场之间的耦合现在完全实现了。
在 COMSOL® 中构建磁流体动力学多物理场模型的图29
显示速度如何耦合到 磁场和电场接口的屏幕截图。

MHD 问题的网格划分和求解

说到单元网格划分和单元阶次,这里一个重要的问题是模型的计算量。求解流体和周围域中的磁场和电场是模型中计算量最大的部分,因此我们希望将整个模型中的网格单元总数保持最少。基于线性静态问题的一些经验法则,我们可以说至少具有二阶单元是一个很好的起点。因此,我们将流体流动的离散化转换为 P2 + P2 离散化,这意味着速度和压力都用二阶基函数来描述。磁场和电场都用二阶离散化来描述。由于所有场都被离散化为至少二阶,因此几何形状的阶次也将自动变为二阶。对可选网格阶次和网格大小的全面调查留给有积极性的读者作为练习。

求解时,软件将自动采用所谓的分离方法,在确定电磁场和速度场之间来回切换,并计算这些场的线性子系统,每个子系统都有自己的优化迭代求解器。由于这种多物理场问题本质上是非线性的,因此了解解决此类问题时可能出现的难题以及如何解决这些难题通常也很有帮助,正如此条知识库条目中关于提高非线性稳态模型的收敛性所述。

多物理场分析的结果如下图所示。我们观察到明显的泵送效应:施加的电压导致电流流过流体,当这些电荷在磁场中移动时,它们会受到一个力的作用,这个力被传递给流体。
在 COMSOL® 中构建磁流体动力学多物理场模型的图30
由于 MHD 多物理场耦合引起的流体泵送的结果图。

简化 MHD 模型

到目前为止,我们已经建立了一个包括磁场、电流和流体流动的模型,我们考虑了所有物理场方程之间的双向耦合。也就是说,每一种物理场现象都会影响其他物理场现象。但事实证明,对于这种特殊情况,我们不需要这样做。接下来我们看一下其中的原因,以及它如何让我们的模型更简单。

如果我们回过头来看一下之前的所有控制方程,我们可以看到只有两个方程引入了物理场现象之间的耦合。方程 在 COMSOL® 中构建磁流体动力学多物理场模型的图31由于电流和磁场而对流体施加一个力,还有一个方程 在 COMSOL® 中构建磁流体动力学多物理场模型的图32。后一个方程表明,电流是由于外加电压边界条件以及导电流体通过磁场的运动而产生的。但是,如果我们假设前一项远大于后一项(即在 COMSOL® 中构建磁流体动力学多物理场模型的图33),那么我们将当前的方程简化为:在 COMSOL® 中构建磁流体动力学多物理场模型的图34。这意味着流体流动问题不会影响电流,流动方程可以与电磁场方程完全分开求解。也就是说,我们可以首先求解电磁场,一旦知道了电磁场,就使用这些场作为流动问题的输入,从而使问题单向耦合。

我们还可以进行额外的简化。严格地说,磁场是由磁铁和电流引起的。然而,对于我们这里分析的边界条件和材料特性,由于电流产生的磁场远小于由磁体引起的磁场。因此,我们可以做出简化的假设,即磁场仅仅是由于磁铁而产生的;也就是说,电流不会产生明显的磁场。这样,我们可以在无电流假设下求解磁场,并分别使用磁场,无电流 电流 接口求解电流。这些物理场接口具有与前面讨论的类似的一组边界和域条件。

磁场,无电流 接口定义了方程 在 COMSOL® 中构建磁流体动力学多物理场模型的图35,该方程的计算量远低于磁场和电场 接口中定义的方程组。此外,这个方程可以独立于电流求解。
在 COMSOL® 中构建磁流体动力学多物理场模型的图36


显示简化模型设置的屏幕截图。

上面的屏幕截图显示了分析这些简化后新模型的设置。流体上体积力的表达式将使用不同的变量名,但除此之外,该模型与之前非常相似。请注意,三个不同的物理场接口在三个单独的研究步骤中求解。磁场,无电流 电流 接口方程可以分别求解,两者都必须在层流 接口方程之前求解。
在 COMSOL® 中构建磁流体动力学多物理场模型的图37
简化的 MHD 模型的结果。

与完全耦合的情况相比,求解这种简化模型时,求解时间将大大减少,这是因为,物理场方程是分开求解的,软件不需要在它们之间进行迭代。从上面显示的结果我们可以看出,这些解几乎与之前未简化的情况相同。当然,我们所做的这些假设和简化确实有其局限性,因此对照完整的模型进行检查是没有坏处的,但 COMSOL Multiphysics 平台具有强大的功能和灵活性,我们可以用它来轻松构建简化模型和完整模型,对它们进行比较,并根据需要进行修改

来自https://cn.comsol.com/blogs/building-a-magnetohydrodynamic-multiphysics-model-in-comsol/作者Walter Frei

COMSOL

在 COMSOL® 中构建磁流体动力学多物理场模型的评论2条

在 COMSOL® 中构建磁流体动力学多物理场模型的相关案例教程

开放群:566811107(资料多,不仅限交流) 群一:836281296 群二:594368389 群三:1080606488 群四: 678357196 我的qq: 209870384有兴趣的可以加我,交流模型。 更多相关分析,可以查看以下链接 基于Comsol的锂电池针刺、内短路和过充仿真 应用COMSOL APP分析锂电池热失控蔓延防控措施 基于comsol的电芯电化学充放电膨胀分析 基于
comsol电弧MHD模型,包括电磁场、传热场和流体场四个物理场之间的多物理场耦合。
积分是数学模型中最重要的功能之一,特别是对数值仿真而言。例如,偏微分方程组 (PDEs) 就是由积分平衡方程派生而来。当需要对偏微分方程进行数值求解时,积分也将发挥非常重要的作用。本文介绍了 COMSOL 软件中可用的积分方法以及如何使用。 积分的重要性 COMSOL 使用了有限元方法,它将控制 PDE 转化为积分方程,换言之,就是弱形式。如果仔细观察一下 COMSOL 软件,您可能会发现许多边界
:870196003 1、裂隙岩体热-流-固耦合数值建模 感兴趣的话点“我想要”和我私聊吧~ 2、COMSOL with MATLAB 代码 二维随机裂隙, 2维随机裂隙生成。功能:可以实现多组不同方向,不同分布规律的裂隙生成(任意组数都可以)。需要输入的参数有:每组裂隙的迹长范围、分布规律(正态分布or随机均匀分布)、倾角。 可以一次生成多组不同产状的裂隙!一次成型,无需CAD等中间工具。 3、
初学者对于comsol一开始怎么学习怎么入门,从哪方面入手等不是很清楚,自己盲目的学习有时候会浪费很多时间,可能效果一般,下面是comsol的仿真案例及软件的基本操作方向,初学的同学可以参考以下内容 COMSOL 仿真实践(RF 及波动光学模块案例 Step by step 详解): 1、光子晶体能带分析、能谱计算、光纤模态计算、微腔腔膜求解; 2、类比凝聚态领域魔角石墨烯的 moiré 光子晶体
CAE工程师
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    2 2