FLUENT二维电弧数值模拟

   电弧的数值模拟面临挑战，因为它涉及电磁场、流体动力学、热传递等多物理场的耦合效应。而Fluent作为一个强大的CFD软件，提供了应对这些复杂耦合问题的解决方案。对于一些对称的电弧现象，二维模型能够在保证计算精度的前提下大幅降低计算量，并提供一个初步的仿真分析。同时，二维仿真可以帮助研究人员快速理解电弧的主要物理行为，为三维仿真提供参考。本文将基于二维仿真，介绍如何使用Fluent进行二维电弧仿真，从模型的建立到结果的后处理。

1、几何建模

  在二维仿真中，几何模型通常可以简化为一个长方形区域，代表电弧周围的流体区域。

图1 求解流程

在SpaceClaim中编辑几何结构，建立的二维结构如图2所示，对其边界与区域进行命名。（采用Ctrl+G建立群组来命名边界，或者在网格划分中进行边界与区域的命名）、

图2 几何模型

2、网格划分

几何模型建立后，网格划分是下一步的关键工作。网格的质量直接影响仿真的准确性和计算速度。在二维电弧仿真中，常使用结构化网格（例如四边形网格）来保证计算的稳定性。网格划分时需要特别注意电弧区域的网格密度，因为该区域内的物理变化较大，需要更细致的网格来捕捉温度场、电流密度等参数的变化。

将绘制好的几何模型导入mesh中进行网格划分，将模型改为流体域，如图3所示。

图3 网格划分

3、物理模型设置

3.1 General（通用）设置

在 Fluent界面左侧求解对话框中的通用设置中完成设置。在此项目中可以更改单位、调整比例尺、检查网格质量、显示网格、设置求解器类型（压力基或密度基）、速度格式（绝对或相对）、设置时间项（瞬态或稳态）、2D空间（平面、轴对称、带旋流的轴对称）、设置重力等。在本文中，我们设置为瞬态，观察电弧变化状态、设置2D空间为轴对称。具体设置如图4所示。

图4  General设置

3.2开启能量方程、加载用户自定义函数

电弧仿真涉及多个物理场的耦合，包括电磁场、热传递等。Fluent可以通过多种模块处理这些物理场的耦合，例如：

电磁场：通过用户自定义函数（UDF）或者电磁场模块定义电流密度、电场强度和磁场分布。

热传递：使用Fluent自带的能量方程，可以精确模拟电弧中的热量传递过程。

根据图示提示，设置相关参数，导入用户自定义函数UDF，编译、加载。

图5  能量设置

图6  UDF设置

3.3 材料属性设置

材料属性对电弧的数值模拟至关重要。对于电弧区域，气体（通常为氩气、氦气或空气）的热导率、比热容、电导率等属性随温度的变化较大。需要使用温度依赖的材料属性来描述这些参数。工作物体和电极通常使用导电材料（如铜或钢），其电导率和热导率也需要根据仿真要求进行定义。Fluent允许在材料库中选择常见的材料，也可以手动定义特定的材料属性。

本文中我们将定义好的材料属性，使用用户自定义数据库导入。如图7所示。

图7 材料设置

3.4 湍流模型

电弧仿真通常涉及高温高速度的流体运动，这些条件下流体容易进入湍流状态。常用的湍流模型包括标准k-ε模型和RNG k-ε模型。RNG k-ε模型对于高曲率流体流动具有更好的预测效果，因此在复杂电弧流动仿真中，建议优先选择该模型。

3.5 单元区域条件

在arc电弧区域内设置源项，选择相应的UDF，材料选择我们自定义的材料属性，如图8所示。

图8 单元区域条件设置

4、边界条件及初始条件设置

4.1 电极和工作物体的边界条件

在电弧仿真中，电极和工作物体的边界条件设置是关键。电极一般设置为电流输入端，用户可以直接定义电流值或者电压值。工作物体则通常作为电流的输出端，设置为接地条件。需要注意的是，在电极和工作物体的交界处，可能存在较大的电流密度梯度，因此需要在这些区域进行局部网格加密。

4.2 流体边界条件

对于流体的边界条件，入口可以设定为速度入口或压力入口，出口可以设定为压力出口。通常建议使用压力入口和出口，因为这样可以更稳定地模拟气体在电弧作用下的流动情况。

4.3 热边界条件

电弧仿真中，电极和工作物体的温度分布也是非常重要的。可以通过设定热流密度或者直接定义边界温度来模拟电弧对材料的加热过程。如果涉及到冷却流体的使用，也需要在入口和出口设置相应的温度条件。

本文中入口设定为速度入口，温度1000K，出口设定为压力出口。anode(阳极)设置温度3000K，在UDS中设定用户自定义标量边界值。cathode(阴极)设置温度1000K，wall设置温度1000K。

5、求解设置

5.1 离散格式选择

离散格式直接影响仿真的精度和收敛性。对于电弧仿真，通常建议使用二阶迎风格式，因为这种格式在捕捉高梯度的流体参数变化时具有更高的精度。

5.2 亚松驰因子设置

将能量改为0.5，其他参数不变。设置限值，将最大静温更改为500000K。如图9所示。

图9 亚松驰因子设置

5.3 初始化

将初始化温度设置为8000K。

5.4 时间步长和收敛准则

电弧仿真通常为瞬态问题，因此时间步长的选择十分重要。时间步长过大会导致仿真结果不准确，而时间步长过小则会增加计算量。建议根据电弧的变化速率选择合适的时间步长，通常为1e-5到1e-6秒。收敛准则可以通过监控残差值以及关键物理量（如电流密度、温度等）的变化来判断，确保仿真达到稳定状态。

详见  Fluent模拟计算不收敛“Floating Point Exception”常见原因及其解决办法-CSDN博客

6、结果后处理

6.1 电弧温度分布

温度场是电弧仿真中最重要的结果之一。可以通过后处理工具提取电弧区域的温度分布，分析电弧核心区域的高温状态以及温度随距离的衰减情况。这对于理解电弧的加热能力以及热影响区的大小具有重要意义。

6.2 速度场和压力场

通过流体速度场和压力场的分布，可以了解电弧对周围气体的推动作用以及产生的压力波动。对于一些高速电弧，如等离子弧，流体速度可能非常高，仿真结果可以用于优化喷嘴设计以及冷却系统的布局。

6.3 电流密度和热通量

电流密度分布直接决定了电弧的加热效率。通过分析电流密度的分布，可以进一步优化电极的设计，保证电弧在整个工作区域内的稳定性。此外，热通量的分析可以帮助研究人员了解热传递的路径和效率，从而为焊接或切割工艺提供数据支持。

6.4计算结果

7、常见问题和解决方案

7.1 收敛问题

在电弧仿真中，收敛问题是常见的难点。残差不收敛的原因可能包括网格质量不佳、时间步长选择不当、初始条件设置不合理等。建议通过逐步细化网格、调整时间步长以及采用更精确的初始条件来解决这些问题。

详见  Fluent模拟计算不收敛“Floating Point Exception”常见原因及其解决办法-CSDN博客

7.2 网格敏感性分析

网格的分辨率对仿真结果有显著影响，特别是在电弧的核心区域。通过网格敏感性分析，可以评估网格尺寸对结果的影响，从而选择最合适的网格划分方案。在实际应用中，过于细致的网格会增加计算时间，因此需要在精度和计算量之间找到平衡。