1. 模型导入与检查别让第一步就“翻车”咱们做仿真第一步就是把模型“请”进Fluent。这听起来简单但新手最容易在这里栽跟头。我见过太多人模型导进去一片空白或者坐标轴飞到了十万八千里外然后就开始怀疑人生觉得软件有问题。其实绝大多数问题都出在导入前的准备和导入后的检查上。首先说导入。原始文章里提到从Workbench的Mesh文件导入这是最常见也最推荐的方式因为它能最大程度地保持网格和几何的关联性。但这里有个细节你保存Mesh文件的路径最好全是英文不要有中文或特殊符号。Fluent对中文路径的支持有时候会“抽风”虽然不一定每次都出错但为了求稳养成用英文命名的习惯能省去很多莫名其妙的麻烦。当你点击“File” - “Read” - “Mesh…”后弹出的浏览窗口里如果发现模型没有立刻显示在图形窗口先别慌。这不一定代表导入失败。很多时候只是图形显示的比例尺不合适模型可能太大或者太小缩成了一个点或者视图方向不对。这时候你需要做的第一件事不是重新导入而是检查模型是否真的被读入了。你可以看看左侧的模型树在“Mesh”节点下应该能看到你导入的网格名称并且后面会跟着网格数量Cells, Faces, Nodes。如果有这些信息恭喜你模型已经成功“上车”了。接下来点击菜单栏的“Display” - “Graphics and Animations”然后在弹出的窗口里双击“Mesh”。在“Mesh Display”的设置面板里确保“Options”下的“Filled”和“Edges”都被勾选上然后点击“Display”。这样你的模型就应该以实体加边线的形式显示出来了。如果还是看不到试试点击视图工具栏的“Auto Scale”那个像放大镜里有个A的图标或者直接按键盘的“A”键让软件自动调整视图到适合模型的大小。模型显示出来只是第一步关键的检查才刚刚开始。你需要像验货一样仔细看看这个“货”对不对版。第一检查模型完整性。用鼠标旋转、缩放模型看看有没有明显的缺失面、破洞或者不应该存在的孤立的网格片。一个快速的方法是在“Mesh Display”面板里勾选“Options”下的“Outline”然后显示这样只会显示模型的外轮廓线更容易发现内部是否有异常的网格隔离区。第二也是原始文章里提到的重中之重检查坐标位置。尤其是做流固耦合、多物理场仿真时流体域和固体域的坐标原点必须对齐。如果两个模型是在不同坐标系下画的导入Fluent后位置对不上那后续的边界条件设置、数据传递全都会出错。你可以在“Display” - “Views…”里查看和设置模型的视图中心坐标。更稳妥的做法是在CAD建模或网格划分阶段就确保所有部件使用统一的全局坐标系原点。2. 仿真环境定调稳态还是瞬态这是个问题模型检查无误接下来就要给整个仿真“定调子”了。这个“调子”就是仿真的物理环境设置它决定了你是在模拟一个稳定不变的状态还是一个随时间变化的过程。这一步选错了后面所有计算都是白费功夫。稳态Steady vs. 瞬态Transient这是第一个重大抉择。很多新手会懵不知道该选哪个。我用一个简单的比喻你想知道一条河在固定水闸开度下长期稳定的水流样子那就用稳态。它计算的是物理量不再随时间变化的最终状态计算速度快资源消耗少。但如果你想研究水闸突然打开时洪水波是如何传播的或者心脏瓣膜在一个心动周期中的开合过程那就必须用瞬态。瞬态模拟时间演变的过程计算量巨大但能捕捉到动态细节。原始文章里提到了设置重力加速度这里我展开说一下重力加速度的设置与否以及其方向直接决定了你的模型是否考虑体积力如自然对流、沉降等。如果你模拟的是管道内的强制对流且管道水平放置重力影响可以忽略那关掉重力即不勾选“Gravity”完全可以还能简化计算。但如果你模拟的是建筑室内的暖气片散热自然对流为主或者水池中的颗粒沉降那重力必须打开并且方向要设对通常是Y轴负方向即-9.81 m/s²。在Fluent中你可以在“General”任务页面找到“Gravity”的复选框和设置栏。设置完这些我强烈建议你养成一个习惯点击“General”任务页面顶部的“Check”按钮。这个操作经常被忽略但它无比重要。Fluent会快速检查你的网格质量、边界条件类型是否合理并在控制台Console输出报告。报告里会告诉你网格的最小体积、最大网格纵横比等关键信息。如果报告里出现“Warning: Minimum volume is negative!”最小体积为负这样的错误说明你的网格存在严重问题比如有负体积网格必须返回Mesh模块修复否则计算根本无法开始。这个“Check”就像开车前的仪表盘自检花5秒钟可能避免后面5个小时的崩溃。3. 湍流模型选择为什么是标准k-ε物理环境设好了现在要告诉Fluent你用什么样的“数学工具”来描述流动。对于绝大多数工程流动问题尤其是涉及机械、暖通、航空航天等领域流动都是湍流。因此选择一个合适的湍流模型至关重要。原始文章直接选择了标准k-ε模型这是一个非常经典和通用的选择。但为什么是它它适合所有情况吗标准k-ε模型可以看作是流体仿真里的“瑞士军刀”它适用范围广计算稳定性好对计算机资源要求相对不高对于很多没有强旋流、强压力梯度、大曲率流动的工程问题比如管道流、外部绕流如汽车、飞机在大部分区域的流动它能给出相当不错的结果。它的核心思想是求解湍动能k和湍流耗散率ε两个输运方程来模拟湍流效应。在Fluent中你可以在“Models”任务页面下点击“Viscous”按钮在弹出的模型列表中选择“k-epsilon (2 eqn)”然后在下面的“k-epsilon Model”中保持“Standard”的选项即可。但是这把“瑞士军刀”也有它的局限。它不适合模拟强旋流如旋风分离器、强压力梯度导致的分离流、以及低雷诺数流动如近壁面处的精细流动。如果你模拟的是飞机机翼失速后的大范围分离或者燃烧室内的剧烈旋流标准k-ε模型可能会严重失真。这时候你可能需要考虑SST k-ω模型或者更高级的雷诺应力模型RSM。对于新手我的建议是如果你不确定先从标准k-ε模型开始。它是个安全的起点。完成一次计算后你可以通过观察残差曲线是否平稳收敛、流场结果是否物理合理比如速度云图没有诡异的彩色斑点来初步判断模型是否适用。如果发现收敛困难或者结果明显不对再回过头来考虑更换更复杂的模型。4. 材料与边界条件给仿真世界注入灵魂模型和环境是骨架湍流模型是肌肉而材料和边界条件就是注入这个仿真世界的灵魂和血液。这一步如果设错了仿真结果会和实际情况南辕北辙。首先是材料。Fluent的默认流体是空气默认固体是铝。对于很多模拟来说直接用空气属性没问题。但如果你模拟的是水、油、或者其他特殊气体怎么办你需要自定义材料。在“Materials”任务页面从“Fluid”或“Solid”的材料库里选择你需要的材料比如“water-liquid (h2o)”然后点击“Copy”到当前工程中。这里有个小坑注意材料的单位制。Fluent默认使用的是国际单位制SI但如果你从库中拷贝的材料属性单位不一致可能会导致计算错误。通常库里的材料属性都是SI制的但检查一下总没坏处特别是密度、粘度这些参数。接下来是重头戏边界条件Boundary Conditions。这是把物理世界约束条件告诉软件的关键步骤。原始文章以速度入口和压力出口为例这是最常见的组合。设置时你需要双击模型树里对应的边界名称如inlet, outlet。对于速度入口除了设置速度大小和方向我特别提醒注意两点第一“湍流”参数的设置方式。对于入口边界湍流强度Turbulent Intensity和水力直径Hydraulic Diameter是更工程化的输入方式。湍流强度可以根据经验公式估算比如对于管道流入口湍流强度I ≈ 0.16 * (Re_DH)^(-1/8)水力直径就是入口截面的特征长度。第二方向矢量。确保速度方向垂直于入口面并且指向计算域内部。如果方向设反了流体就“倒灌”了。对于压力出口通常保持默认的“Gauge Pressure”为0 Pascal即可表示出口静压为环境大气压。但如果你模拟的是有背压的系统比如发动机排气管就需要设置一个非零的出口压力。设置边界条件时脑子里要时刻想着你模拟的物理场景到底是什么样的每一个数字都要有物理依据不要随便填。5. 求解器设置与初始化让计算稳如泰山所有条件都设好了是不是马上就可以点“Calculate”开始算了别急最后这一系列求解器设置决定了你的计算是“一路顺风”还是“半路抛锚”。首先是求解方法Solution Methods。在“Solution Methods”任务页面你会看到压力-速度耦合方式、空间离散格式等选项。对于新手我建议保持默认的“SIMPLE”算法和“Second Order Upwind”离散格式。SIMPLE算法是稳态计算最常用的稳定性好。“Second Order Upwind”二阶迎风格式比一阶格式精度高能减少数值扩散让结果更准确。虽然计算量稍大一点但对于现在的电脑来说这点开销完全值得。接下来是监控器Monitors。这里你一定要设置残差监视器。残差曲线是你观察计算是否收敛的“心电图”。在“Monitors”下双击“Residuals”可以设置收敛判据。默认的各方程残差收敛标准是1e-3对于大多数工程问题够用了。你可以把能量方程的收敛标准设得更严格些比如1e-6如果能量交换是你的关注重点的话。除了残差强烈建议你创建力、流量或某点物理量的监视器。比如你模拟一个机翼可以创建一个升力、阻力的监视器模拟管道流可以创建一个出口质量流量的监视器。有时候残差曲线看起来波动很大但你看力的监视曲线可能已经稳定不变了这说明计算其实在物理层面上已经收敛了可以手动停止。然后是初始化Initialization。原始文章提到了“标准初始化”。这是最常用的方法它会根据你设置的边界条件为整个流场提供一个合理的初始猜测值。点击“Initialize”后建议先不要直接开始计算。进行一次“Patch”操作比如如果你知道入口附近流速高你可以在那个区域局部初始化一个更高的速度值这能显著加快计算收敛速度。在“Solution Initialization”面板里选择“Patch…”然后选择你要初始化的区域和变量值即可。6. 计算运行与自动保存设置好就让它跑起来万事俱备只欠“计算”。在“Run Calculation”任务页面你需要设置两个关键参数迭代步数Number of Iterations和报告间隔Reporting Interval。对于稳态计算你设置一个足够大的迭代步数比如1000步然后点击“Calculate”让它跑。你不需要设置时间步长。计算过程中要紧盯控制台窗口和残差曲线图。如果残差曲线平稳下降并最终低于你设置的收敛标准同时你设置的力、流量等监视器也基本不变了那就可以手动点击“Stop”停止计算。如果跑了设定的步数比如1000步还没收敛你可以继续增加步数再算。对于瞬态计算设置就复杂一些了。你需要设置时间步长Time Step Size和总步数Number of Time Steps。时间步长的选择非常有讲究步长太大计算会不稳定甚至发散步长太小计算时间会长得让你绝望。一个经验法则是时间步长应该小于流体穿过最小网格单元所需时间的1/5。你可以先设一个较小的步长试算几步观察残差如果稳定再尝试调大。总步数 总模拟时间 / 时间步长。最后也是保命的一招设置自动保存Autosave。尤其是对于瞬态计算或者需要跑很久的稳态计算一定要设你永远不知道电脑会不会突然死机、断电。你可以在“Calculation Activities”任务页的“Autosave”中设置每隔多少步或每隔多少物理时间自动保存一次案例和数据文件。这样即使中断你也可以从最近的一个保存点重启计算避免前功尽弃。根据我的经验稳态计算可以设置每100-200步保存一次瞬态计算可以根据你的数据输出频率来设置。