4. Fluent求解器基础设置

在换热器仿真中,几何与网格只是“骨架”,真正决定物理场是否可信的,是求解器启动后的一系列基础设置。本节以一台典型的管壳式换热器(壳程走高温导热油、管程走冷却水)为例,说明从启动 Fluent 到完成材料、操作条件、求解器类型、单位与参考值设定的完整过程。示例中的物性参数与工况均为教学假设值,实际工程应以设计资料或物性数据库为准。

一、启动流程与求解器选择

Fluent 提供两种启动方式:Workbench 中的 Fluent 模块(推荐,便于与 Meshing、Geometry 参数联动)和独立启动的 Fluent Launcher。在 Launcher 中需要确认三项关键选择:Dimension(换热器通常为 3D)、Processing Options(Serial 或 Parallel,管壳式换热器网格量常在百万级,建议 4~8 进程并行)、Solver(Pressure-Based 或 Density-Based)。

求解器类型的选择依据是流体是否可压、马赫数高低以及是否存在强激波。换热器内液体流速一般 0.5~3 m/s,气体也多为中低速,马赫数远小于 0.3,因此压力基求解器(Pressure-Based)是默认且正确的选择。若壳程为高压蒸汽且流速接近声速,才考虑密度基。压力基下还有两种算法:分离式(Segregated)与耦合式(Coupled)。分离式内存占用小、适合大网格;耦合式收敛更快、对强耦合流动(如自然对流显著、浮升力主导)更稳健。本案例壳程存在明显温差引起的浮升力,建议采用 Coupled 算法,压力速度耦合选用 Coupled 格式。

时间类型上,稳态换热器设计工况用 Steady;若关注启停、流量突变或结垢瞬态,则选 Transient。启动后还需在 General 面板确认 Gravity 是否勾选——壳程自然对流不可忽略时应勾选,并按坐标方向设定重力分量(如 Y 方向 -9.81 m/s²)。

二、材料属性定义

换热器仿真至少涉及三种材料:管程水、壳程导热油、管壁金属(如不锈钢)。在 Materials 中可从 Fluent 数据库复制(Copy)后修改,也可新建(Create/Edit)。关键物性包括:密度、比热、导热系数、粘度。若温度变化范围大(如壳程进出口温差超过 40 K),应使用温度相关的分段线性(piecewise-linear)或多项式拟合,而非常数。

以壳程导热油为例,教学假设其物性随温度变化如下表。定义时在 Density 下拉选择 piecewise-linear,依次输入温度—密度点对;比热与导热系数同理。粘度对传热影响显著,建议用分段线性或 Sutherland 类关系描述。

温度 (K)密度 (kg/m³)比热 (J/(kg·K))导热系数 (W/(m·K))粘度 (Pa·s)
30086020000.1300.0032
35083021500.1250.0018
40080023000.1200.0011

管程水可用常数物性(密度 998.2 kg/m³、比热 4182 J/(kg·K)、导热系数 0.6 W/(m·K)、粘度 0.001003 Pa·s),若考虑温差较大也可启用分段线性。管壁不锈钢取密度 8030 kg/m³、比热 502 J/(kg·K)、导热系数 16.27 W/(m·K)。定义完成后,需在 Cell Zone Conditions 中把各流体域与固体域分别指派到对应材料,这一步常被遗漏,导致后续温度场异常。

三、操作条件设置

操作条件(Operating Conditions)主要设定操作压力与重力。对于不可压或弱可压流动,操作压力作为参考压力,实际压力 = 表压 + 操作压力。换热器壳程若为常压,操作压力取 101325 Pa;若为带压运行(如 0.6 MPa 表压),应把操作压力设为 0,而把边界条件中的压力设为绝对压力,避免参考压力与边界压力重复叠加造成错误。

重力设置需与坐标方向一致。若几何建模时 Y 轴为竖直方向,则重力在 Y 方向取 -9.81 m/s²;若 Z 轴竖直,则改为 Z 方向。操作密度(Operating Density)在封闭或半封闭自然对流问题中建议指定,取壳程平均温度下的密度,可加快收敛并减少浮力计算误差。对于本案例,可勾选 Gravity 并指定 Operating Density 为 830 kg/m³。

四、单位与参考值设定

Fluent 默认使用 SI 单位,但界面可切换。建议在 General → Units 中统一检查:长度 m、压力 Pa、温度 K、速度 m/s。若几何以 mm 建模,导入后需在 Scale 中换算,否则边界速度与雷诺数会差三个数量级。

参考值(Reference Values)用于计算无量纲系数,如升力系数、阻力系数、努塞尔数。换热器仿真中常关注压降与传热系数,参考值至少应设定:参考面积(取壳程流通截面积或管束迎流面积)、参考速度(壳程入口速度)、参考密度(壳程平均密度)、参考温度(入口温度或平均温度)。这些值可在 Reference Values 面板中由 Fluent 自动从某入口计算(Compute from),也可手动输入。参考温度对努塞尔数定义影响很大,全文与后处理必须保持一致。

建议在完成上述设置后,先初始化并做 10~20 步迭代,观察残差与进出口质量守恒是否合理,再正式计算。若残差震荡或质量不守恒,优先检查材料指派、操作压力与重力方向,而不是急于调松弛因子。
Fluent 换热器求解器基础设置关系图 启动 → 求解器类型 → 材料 → 操作条件 → 单位/参考值 → 初始化检查 ① 启动与求解器 Dimension:3D Processing:Parallel 4~8 Solver:Pressure-Based 算法:Coupled + Coupled 时间:Steady / Transient ② 材料属性 管程水:ρ=998.2 壳程油:分段线性 管壁:不锈钢 16.27 指派到 Cell Zone 温度相关物性优先 ③ 操作条件 操作压力:101325 Pa 带压时操作压力设 0 Gravity:Y = -9.81 Operating Density 830 方向与建模坐标一致 ④ 单位与参考值 长度 m / 压力 Pa 温度 K / 速度 m/s 参考面积:壳程截面积 参考速度:入口速度 参考温度:入口温度 ⑤ 初始化与首轮检查(10~20 步) • 残差是否持续下降,连续性残差是否低于 1e-3 • 进出口质量守恒误差是否小于 0.5% • 壳程温度场是否出现非物理过冲或局部异常 关键参数对换热器结果的影响链 物性随温度变化 → 传热系数偏差 操作压力设定 → 绝对/表压混淆 重力方向 → 自然对流强弱 参考温度 → Nu 定义一致性 收敛性 → 结果可信度

完成上述五类设置后,换热器仿真才具备物理意义上的可解性。实际工程中,建议把材料物性来源、操作压力基准、重力方向与参考值记录在仿真日志中,便于后续复现与对比。下一阶段即可进入边界条件与求解方法的细化设置。