之前用 STAR-CCM+ 导入复杂 CAD 模型时，总是卡在几何清理和表面网格这一步：破面多、小特征杂乱、边界命名不清晰，导致体网格生成失败，反复调整非常浪费时间。后来把 HyperMesh 作为前处理工具，把几何清理和表面网格工作前置，整个 CFD 流程才顺畅起来。这篇文章就把这套基础流程完整记录下来，覆盖从 HyperMesh 网格划分、STAR-CCM+ 体网格生成、物理模型设置、求解到后处理的闭环操作，适合刚接触 CFD 仿真的同学，也适合想规范前处理流程的工程师参考。
1. CFD 仿真与工具链基础
1.1 HyperMesh 与 STAR-CCM+ 的定位
HyperMesh 是 Altair 公司推出的一款有限元前后处理软件，在结构分析和 CFD 前处理中应用非常广泛。它的强项是几何清理、网格划分和网格质量控制，尤其是面对复杂 CAD 模型时，可以快速处理小孔、小圆角、干涉面等常见几何问题。
STAR-CCM+ 则是 Siemens 旗下 Simcenter 平台中的 CFD 仿真软件，基于有限体积法完成流动和传热等物理场求解。它最大的特点是把网格生成、求解、后处理放在同一个环境中，内置了多面体网格、棱柱层网格和表面重构等功能，对复杂外流场和内部流动都有很强的适应性。
在完整 CFD 流程中，这两个工具的定位并不重叠。HyperMesh 负责“把几何变成高质量的表面网格”，STAR-CCM+ 负责“把表面网格变成体网格并完成物理场求解”。很多人习惯直接在 STAR-CCM+ 中导入原始 CAD 几何，然后通过包面（Surface Wrap）来处理模型缺陷，这种方式虽然可行，但对于包含大量细节特征的三维模型，前处理过程会比较依赖软件自动算法，控制力不如 HyperMesh 手动清理来得直接。
1.2 CFD 仿真为什么需要“前处理”
CFD 仿真的结果是否可信，很大程度上取决于网格质量。一套穿透、重叠、畸变的网格，即使在求解器设置上花费再多精力，最终也可能得到发散或者精度极低的结果。
前处理阶段主要完成三件事：
- 几何清理：删除多余的小特征、修复破面和缝隙，让计算域在拓扑上是干净的。
- 计算域设计：确定流体区域范围，比如外流场需要足够大的远场边界，内部流场需要合理设置入口和出口。
- 表面网格生成：用三角形或四边形单元离散几何表面，为后续体网格生成提供封闭的输入。
STAR-CCM+ 自带的网格工具虽然也能完成这些操作，但 HyperMesh 在几何修复和网格控制上更灵活。特别是当模型来自不同 CAD 软件、格式转换后出现破损时，HyperMesh 的几何编辑能力往往比直接进入 STAR-CCM+ 更高效。
1.3 典型 CFD 仿真流程
一个完整的 CFD 仿真流程，通常可以拆成下面几个阶段：
- 几何准备：导入 CAD 模型，简化结构，创建流体计算域。
- 表面网格：在 HyperMesh 中划分三角形面网格，并规划边界命名。
- 体网格：在 STAR-CCM+ 中导入面网格，生成多面体网格和棱柱层网格。
- 物理设置：选择物理模型，设置边界条件、初始条件和求解参数。
- 求解计算：迭代求解，监控残差和关键物理量。
- 后处理：提取压力场、速度场、流线、阻力系数等结果。
这个流程具备通用性。无论是无人机外流场、飞机外流场，还是汽车空气动力学分析，核心步骤都大同小异。初学者先跑通一个最小闭环，后续再针对具体工程问题扩展细节，会更扎实。
2. 环境准备与版本说明
2.1 软件与硬件环境
本文的演示环境以常见 CFD 工作站为例：
- 操作系统：Windows 10/11 或 Linux 均可，STAR-CCM+ 在两种系统下都能稳定运行。
- HyperMesh：来自 Altair HyperWorks 套件，版本不同界面布局略有差异，但核心面板名称基本一致。
- STAR-CCM+：Simcenter STAR-CCM+，版本不同不影响整体流程理解。
- 硬件建议：处理器核心数越多越好，建议 8 核以上；内存建议 16GB 起步，复杂三维外流场建议 32GB 以上；显卡要求不高，能流畅旋转模型即可。
需要注意的是，HyperMesh 和 STAR-CCM+ 都属于商业软件，需要有效的许可证才能运行。学校用户可以通过学校授权申请试用，企业用户可以使用公司购买的正式许可。本文不涉及任何安装破解相关内容，只讨论正常使用流程。
2.2 文件格式与单位约定
HyperMesh 和 STAR-CCM+ 之间的数据交换，常见有两种路线：
- 几何路线：在 HyperMesh 中完成几何清理后，导出 STEP（.stp）或 IGES（.igs）文件，再导入 STAR-CCM+ 重新划分网格。
- 网格路线：在 HyperMesh 中直接生成表面网格，导出为 Nastran（.nas/.bdf）、Abaqus（.inp）或 STL 格式，再导入 STAR-CCM+ 生成体网格。
本文重点介绍第二种路线，也就是“HyperMesh 出表面网格 → STAR-CCM+ 出体网格”的配合方式。这样能最大化发挥两个软件各自的优势。
单位是一个很容易出错的点。建议在项目开始前就确定统一单位系统。STAR-CCM+ 默认使用国际单位制（SI），长度单位是米，质量单位是千克，时间单位是秒。如果 HyperMesh 中模型以毫米为单位，那么导出前一定要记录好比例关系，导入 STAR-CCM+ 时再通过缩放或单位设置恢复正确尺寸。
2.3 许可证与并行计算
STAR-CCM+ 支持多核并行计算。在图形界面运行时，可以直接在 Run 面板中设置并行进程数。如果需要批处理计算，可以使用命令行方式提交任务：
starccm+ -batch -np 8 run_simulation.java case.sim其中-batch表示后台批处理模式，-np 8表示使用 8 个进程，run_simulation.java是录制好的宏文件，case.sim是仿真工程文件。这种方式在服务器集群上更常用，本地工作站也可以使用。
3. 几何准备与流体域设计
3.1 导入 CAD 几何
在 HyperMesh 中，第一步是导入原始 CAD 模型。以无人机外流场分析为例，CAD 模型通常以 STEP 或 IGES 格式提供。操作路径一般是：
File → Import → Geometry导入后，先在 HyperMesh 中检查模型是否存在破面、缝隙、重叠面等几何缺陷。可以通过 Hide/Show 功能逐层查看零件，确认模型的流场表面是否完整。如果模型由多个部件组成，建议在导入前先确认不需要的零部件是否已经删除，比如无人机内部的电池、电路板、线缆等，对外流场分析的影响通常很小，但会显著增加网格划分难度。
3.2 几何清理与简化
几何清理是 CFD 前处理中最费时间、也最影响后续结果的一步。常见清理对象包括：
- 小圆角和小倒角：直径小于基础网格尺寸的圆角可以直接删除，它们对流动影响很小。
- 小孔：直径远小于模型特征尺寸的孔，可以用几何编辑工具填充。
- 破碎短边：多条短线相交导致的拓扑混乱，需要合并或重生成边。
- 重叠面：两个面之间存在干涉时，需要用曲面编辑工具修正。
HyperMesh 中提供了 Defeature、Edit Geometry 等工具，可以批量删除小特征。清理的原则是“保留对流动有影响的几何特征，去除对网格尺寸不敏感的小特征”。以无人机为例，机臂的倒角可以简化，但螺旋桨旋转区域附近如果有防护罩，则需要保留完整的几何外形。
3.3 创建外流场计算域
外流场分析需要创建一个足够大的计算域，用来模拟无界空间中的流动。计算域尺寸太小，边界会对绕流结果产生明显干扰；尺寸太大，网格数量会急剧增加，浪费计算资源。
经验上，计算域的特征尺寸通常是物体特征长度的 10 到 20 倍。以轴距 500mm 的小型四旋翼无人机为例，外流场可以设置为一个长 10m、宽 6m、高 6m 的长方体区域。来流方向从无人机前方进入，因此入口到无人机前端的距离建议留 4m 以上，出口到无人机后端的距离可以更长，让尾流充分发展。
在 HyperMesh 中，计算域可以通过创建实体几何得到。实际操作时，可以在 CAD 软件中提前建好带计算域的装配体，也可以直接在 HyperMesh 中创建实体长方体，然后与无人机模型放在同一个模型中。这样后续划分网格时，只需要对计算域外表面和无人机表面分别生成面网格。
4. HyperMesh 网格划分实战
4.1 规划 Component 与边界命名
在 HyperMesh 中，所有几何和网格都归属于不同的 Component。为了让 STAR-CCM+ 导入后能直接识别边界，建议在划分网格之前就规划好 Component 结构。
以无人机外流场为例，可以创建以下 Component：
| Component 名称 | 含义 | 后续边界类型 |
|---|---|---|
| UAV_BODY | 无人机机身表面 | Wall |
| FARFIELD_INLET | 计算域入口 | Velocity Inlet |
| FARFIELD_OUTLET | 计算域出口 | Pressure Outlet |
| FARFIELD_SIDE | 计算域侧面和顶底面 | Symmetry |
命名规范尽量统一。推荐使用大写字母加下划线，例如UAV_BODY、FARFIELD_INLET。避免出现中文名、空格和特殊符号，因为这些字符在跨软件传递时可能出现兼容性问题。
在每个 Component 中，需要把对应的几何面移动进去。HyperMesh 中可以用Organize功能把几何移动到指定 Component。这个步骤看似繁琐，但关系到后续 STAR-CCM+ 中边界条件的设置效率，值得投入时间。
4.2 生成三角形面网格
CFD 表面网格通常使用三角形单元，因为三角形对复杂曲面的适应性更好。STAR-CCM+ 导入三角形面网格后，可以自动生成多面体体网格，因此 HyperMesh 中不需要生成四边形网格。
在 HyperMesh 中，进入 2D 网格划分面板，通常路径是：
2D → Automesh然后在空白区域选择要划分的几何，设置单元尺寸。单元尺寸需要根据模型实际尺寸和计算资源确定。以无人机模型为例，基础尺寸可以设为 5mm，机臂、桨叶等小特征区域可以加密到 2mm，远场边界可以放大到 20mm。
HyperMesh 还支持基于曲率的自动细化，在曲率较大的区域自动加密网格单元，这对于保留模型几何特征非常重要。如果只使用统一尺寸，曲面过渡区域容易出现较差的三角形。
4.3 网格质量检查与修复
网格划分完成后，一定要做质量检查。CFD 网格对质量的要求比结构网格更严格，因为畸变网格会直接影响求解收敛性。
在 HyperMesh 中，可以使用Check Elements功能检查三角形单元：
- 最小角度：建议大于 20 度，最好大于 30 度。
- 最大角度：建议小于 120 度。
- 长宽比（Aspect Ratio）：建议小于 5:1。
- 塌陷系数：越接近 0 越好。
如果发现质量较差的单元，可以使用Quality Index面板进行局部加密或重绘，也可以手动修改节点位置。对于 CFD 仿真，重点关注最小角度和长宽比这两个指标。出现大量畸形单元时，不要试图在 STAR-CCM+ 中补救，建议回到 HyperMesh 中重新调整网格参数。
4.4 导出面网格
面网格质量检查通过后，就可以导出网格文件。在 HyperMesh 中导出 Nastran 格式文件，路径通常为：
File → Export → Simulation → Nastran导出时需要注意模板选择。HyperMesh 针对不同求解器提供了不同导出模板，这里应该选择与 CFD 或者通用 Nastran 格式匹配的模板，确保 2D 单元信息能正常写入文件。
导出后的 .nas 文件会包含所有三角形单元的节点和单元信息，同时也会保留 Component 名称。这个文件就是接下来要导入 STAR-CCM+ 的文件。
如果后续发现 STAR-CCM+ 导入后边界名称没有正确识别，可以检查导出时是否选择了正确的单元类型和模板，尽量不要用 STL 格式替代，因为 STL 格式会丢失边界命名信息。
5. STAR-CCM+ 中导入与体网格生成
5.1 导入面网格与单位设置
打开 STAR-CCM+ 后，通过以下路径导入 HyperMesh 导出的文件：
File → Import → Surface Mesh选择之前导出的.nas文件。在导入向导中，STAR-CCM+ 会询问单位信息。如果模型在 HyperMesh 中以毫米为单位，这里需要选择毫米，或者手动输入缩放比例，让软件将模型尺寸换算为米。
导入完成后，模型会出现在图形窗口中。此时可以打开表面网格（Surface Mesh）的 Representation 节点，查看网格信息。
以下是一个简化版的 STAR-CCM+ 宏示例，用于导入表面网格并创建基本物理连续体。不同版本 API 会有细微差异，实际使用建议先通过宏录制功能生成基础代码：
// 文件：setup_macro.java import star.common.*; import star.base.neo.*; import star.flow.*; import star.turbulence.*; public class setup_macro extends StarMacro { public void execute() { Simulation simulation = getActiveSimulation(); // 导入表面网格 simulation.getImportManager().importFile("C:/cfd/uav_mesh.nas", 0, "user"); // 创建物理连续体 PhysicsContinuum physics = simulation.getContinuumManager() .createPhysicsContinuum(); physics.enable(ThreeDimensionalModel.class); physics.enable(SteadyModel.class); physics.enable(GasModel.class); physics.enable(SegregatedFlowModel.class); physics.enable(ConstantDensityModel.class); physics.enable(KEpsilonTurbulence.class); } }这段代码只展示了导入和连续体创建的思路，实际项目中可以直接在 STAR-CCM+ 图形界面完成，不一定需要编写宏。
5.2 检查表面网格与边界
导入表面网格后，需要检查网格是否满足生成体网格的条件。最重要的一点是表面网格必须封闭，不能有自由边。
在 STAR-CCM+ 中，可以通过右键 Representation 节点，选择 Diagnose 来检查自由边和非流形边。也可以通过视图工具直接高亮显示自由边。如果存在自由边，说明表面网格存在缝隙，需要在 HyperMesh 中修复后重新导出，或者在 STAR-CCM+ 中执行 Surface Repair。
同时检查边界列表。如果规划正确，导入后的 Region 下应该能看到UAV_BODY、FARFIELD_INLET、FARFIELD_OUTLET、FARFIELD_SIDE等 Boundary。然后根据需要将对应边界设置为 Wall、Velocity Inlet、Pressure Outlet 或 Symmetry Plane。
边界类型设置非常关键，它决定了求解器如何施加边界条件。入口和出口如果设置错误，会导致流动方向错误或者压力异常。
5.3 生成多面体网格与棱柱层
STAR-CCM+ 的核心网格能力是多面体网格（Polyhedral Mesher）。相比四面体网格，多面体网格在相同精度下单元数量更少，收敛速度更快，对复杂几何的适应性也更好。因此，在 STAR-CCM+ 中推荐使用多面体网格。
设置体网格时，展开 Operation，创建新的 Mesh Operation，添加以下网格器：
- Surface Remesher：对导入的三角形表面网格进行重新光顺，提高表面网格质量。
- Polyhedral Mesher：生成多面体体网格。
- Prism Layer Mesher：在壁面附近生成棱柱层网格，用于捕捉边界层流动。
Prism Layer 参数是 CFD 网格的关键。对于外流场分析，层数通常设置为 3 到 5 层，总厚度需要根据预估的边界层厚度确定。层数太少，近壁面速度梯度捕捉不足；层数太多，网格数量激增，计算开销变大。
Base Size 的设定也很重要。如果导入的表面网格尺寸在 5mm 到 20mm 之间，Base Size 可以设为 0.005m 或 0.01m，然后在部件表面设置局部加密尺寸。
5.4 体网格质量检查
体网格生成完成后，还需要做一次质量检查。在 STAR-CCM+ 中，可以通过网格统计信息查看：
- 负体积单元数量：必须为 0。
- 体积变化率：建议大于 0.1。
- 单元倾斜度（Skewness Angle）：建议小于 85 度。
如果存在负体积，需要回到网格参数中调整棱柱层厚度或者表面网格尺寸。负体积单元会导致求解发散，即使数量很少也不能忽略。
如果棱柱层生成失败，常见原因是表面网格局部质量太差，或者棱柱层总厚度超过局部几何特征尺寸。此时可以适当减小棱柱层厚度，或者在 HyperMesh 中加密对应区域的表面网格。
6. 物理模型与边界条件设置
6.1 物理连续体
在 STAR-CCM+ 中，物理连续体（Physics Continuum）定义了计算域内流体的物理模型。以低速不可压外流场为例，可以按照下面方式设置：
- 空间维度：三维（Three Dimensional）
- 时间属性：稳态（Steady）
- 材料：空气（Gas）
- 流动模型：分离流（Segregated Flow）
- 密度属性：常密度（Constant Density）
- 湍流模型：k-epsilon 或 k-omega SST
对于初学阶段，k-epsilon 模型参数稳定、容易收敛，适合快速验证流程。如果关注分离流动或者旋涡结构，可以改用 k-omega SST。对于高速可压缩流动，比如飞机外流场，则需要启用 Coupled Flow 模型和理想气体状态方程，计算难度会明显提升。
6.2 边界条件设置
以外流场阻力分析为例，来流速度通常取 10 m/s。设置方式如下：
| 边界名称 | 边界类型 | 参数设置 |
|---|---|---|
| FARFIELD_INLET | Velocity Inlet | 速度 = 10 m/s，方向沿 X 正向 |
| FARFIELD_OUTLET | Pressure Outlet | 表压 = 0 Pa |
| UAV_BODY | Wall | 无滑移壁面，粗糙度忽略 |
| FARFIELD_SIDE | Symmetry Plane | 无需参数 |
需要注意的是，入口速度参考系要与计算域方向一致。如果模型放置角度有偏转，可以在入口边界中设置速度分量，而不是旋转模型。这样调整起来更灵活。
Pressure Outlet 表压通常设为 0 Pa，表示出口通向环境大气。对于外流场，如果计算域足够大，入口和出口边界对物体附近的流动影响很小。
6.3 初始化与停止条件
求解之前需要设置初始条件和停止条件。初始条件合理的话，求解器可以更快收敛。对于 10 m/s 来流，初始条件可以设置为：
- 速度：10 m/s
- 压力：0 Pa
- 湍动能：一个较小值，比如 0.001 m²/s²
- 湍流耗散率：一个较小值，比如 0.001
停止条件一般有两种方式：按迭代步数或按残差阈值。推荐设置最大迭代步数为 1000，同时开启残差监控，当残差下降到 1e-4 以下并且不再继续下降时，可以手动停止计算。
7. 求解与后处理
7.1 求解器参数
在 STAR-CCM+ 的求解器中，默认参数通常可以支持常规模型稳定运行。如果残差下降缓慢，可以适当调整：
- 降低速度或压力松弛因子，例如从 0.7 降到 0.3。
- 如果使用 k-epsilon 模型，降低湍流相关松弛因子。
- 更换初始化方法，先使用一步势流初始化，再进入完整求解。
批处理运行时，可以在命令行中直接调用宏完成求解和结果保存：
starccm+ -batch -np 8 solve_and_save.java case.sim这种方式适合在夜间或者服务器上集中计算，避免长时间占用图形界面。
7.2 收敛判断
CFD 求解是否收敛，不能只盯着残差曲线。建议同时监控以下三个指标：
- 残差曲线：一般要求下降到 1e-4 到 1e-5，并且保持平稳。
- 力和力矩系数：无人机机身的阻力和升力系数是否基本稳定，不再出现周期性大幅波动。
- 全局质量守恒：入口和出口的质量流量是否基本相等。
如果残差已经很低，但阻力系数仍然持续变化，说明流场可能处于周期性非定常状态，比如卡门涡街。此时稳态求解器很难收敛到一个固定解，需要切换到瞬态计算。
一个便捷的判断方式是，创建 Report 和 Monitor，实时绘制阻力系数随迭代步数变化的曲线。曲线近似水平时，可以认为计算已经收敛。
7.3 结果提取
后处理常用的操作包括：
- 创建压力云图：选择机身表面，用 Scalar Scene 显示压力分布。
- 创建速度场剖面：在无人机后方创建平面，显示速度云图。
- 创建流线：从入口边界发射流线，观察无人机周围的流动结构。
- 导出阻力系数：通过 Report 功能计算阻力，并导出为 CSV 文件。
在图形窗口中，可以用鼠标旋转和缩放模型，检查流动细节。如果发现无人机尾部出现大范围分离区，说明阻力较大，后续可以通过优化机身外形来改善。
8. 常见问题与排查思路
下面这个表格总结了 HyperMesh 与 STAR-CCM+ 联合仿真中最常见的问题：
| 问题现象 | 常见原因 | 解决思路 |
|---|---|---|
| 导入网格后尺寸偏大或偏小 | 单位设置不一致 | 确认导入时单位，必要时手动缩放 |
| 体网格生成失败 | 表面网格存在自由边或缝隙 | 在 HyperMesh 中修复几何并重新导出 |
| 边界名称全部丢失 | 导出格式不支持边界命名 | 改用 Nastran 格式，避免 STL |
| 棱柱层生成失败 | 表面网格质量差或厚度过大 | 减小棱柱层厚度，加密表面网格 |
| 残差持续不下降 | 边界条件或湍流模型不合理 | 检查边界设置，尝试调整初始条件 |
| 求解发散出现 NaN | 负体积或初始化不当 | 检查网格负体积，降低松弛因子 |
导入网格后尺寸异常是最容易踩的坑。比如 HyperMesh 中以毫米建模，导出后直接导入 STAR-CCM+，并且选择单位是米，模型就会放大 1000 倍。遇到这种情况，可以先测量模型上一个已知特征的长度，快速判断单位是否正确。
边界名称丢失的问题比较隐蔽。STL 格式只有三角形网格信息，无法保存边界名称，所以如果直接导出 STL，导入 STAR-CCM+ 后所有边界会合并成一个。建议使用 Nastran 格式导出面网格，并且在 HyperMesh 中规划好 Component 名称。
残差不下降时，不要一开始就认为是网格问题。可以先检查边界条件是否与物理场景一致，例如入口速度方向是否正确、压力出口位置是否太靠近物体。很多情况下，调整初始条件后残差会明显改善。
9. 最佳实践与工程建议
结合实际项目经验，下面这些建议能帮助你减少返工，提升仿真效率。
提前规划边界命名。在 HyperMesh 中划分网格之前，就把每个边界的 Component 名称定义好。命名规范建议采用全大写加下划线，例如UAV_BODY、FARFIELD_INLET。这样导入 STAR-CCM+ 后可以直接识别和设置边界条件，避免手动逐个选择几何面。
全程保持统一单位。建议在项目初期就确定单位系统，推荐使用 SI 单位。如果 CAD 模型以毫米为单位，导出和导入时要记录清楚，并利用软件的单位缩放功能转换。
先粗网格后细网格。第一次试算时，不要追求非常细的网格。先用较大的网格尺寸跑通完整流程，确认边界条件和物理模型设置无误，再逐步细化网格做网格无关性验证。这样能节省大量调试时间。
重视棱柱层参数。外流场分析中，壁面附近的网格对阻力和升力的计算精度影响很大。棱柱层层数、总厚度和近壁面第一层高度都需要根据预估的 y+ 值调整。对于不使用壁面函数的低雷诺数模型，第一层高度尤其关键。
保存多个版本。网格修改、物理模型调整、边界条件变更之前，都建议另存一份副本。仿真分析是一个反复迭代的过程，版本管理可以让你随时回退到之前可用的状态。
验证结果合理性。仿真完成后，可以用简单估算或实验数据交叉验证。比如无人机阻力系数是否在合理范围，压力分布是否符合常规气动规律。如果结果明显违背物理常识，优先检查网格和边界条件，而不是继续细化后处理。
10. 总结与下一步
本文从 HyperMesh 的几何清理和表面网格划分开始，逐步完成了 STAR-CCM+ 中的体网格生成、物理模型设置、求解和后处理，走通了一个完整的 CFD 仿真基础流程。重点内容包括边界规划、三角形面网格质量检查、多面体网格与棱柱层设置、边界条件配置和收敛判断方法。
下一步可以围绕三个方面继续深入：第一，做网格无关性验证，研究网格尺寸对阻力系数的影响；第二，对比 k-epsilon 和 k-omega SST 两种湍流模型的计算差异；第三，尝试瞬态计算，观察无人机尾流中的旋涡脱落现象。完成这些练习后，你会对 CFD 仿真的数据和误差来源有更深刻的理解。
实际的工程问题不会像教程示例这样干净，复杂几何、多物理场耦合、计算资源限制都是常见挑战。但只要掌握了这套基础流程，后续无论面对汽车外流场、飞机外流场还是散热分析，都可以沿着同样的思路去拆解和推进。先完成一个最小闭环，再逐步扩展，这是学习 CFD 仿真最有效的方式。
