本文基于 Fluent 2025R2 官方用户手册,系统梳理动网格的完整设置流程、三种核心更新方法(平滑、动态层铺、重构)及其详细参数、内燃缸专用功能、六自由度(6 DOF)求解器、隐式更新、接触检测、事件定义以及稳态动网格应用。
无论是瞬态还是稳态问题,设置动网格的标准步骤如下。注意:此流程仅包含动网格模型本身必需的步骤,其他模型(湍流、传热、多相流等)按常规设置即可。
启用瞬态或稳态流动
Setup → General准备相邻运动区域的网格拓扑
设置单元区域条件和边界条件
Setup → Cell Zone Conditions Setup → Boundary Conditions启用动网格模型并设置参数
Setup → Dynamic Mesh → Dynamic Mesh创建动态区域
Setup → Dynamic Mesh → Create/Edit...(可选)显示预设运动
Setup → Dynamic Mesh → Display Zone Motion... 定义事件(仅瞬态)
Setup → Dynamic Mesh → Events...保存 Case 和 Data 文件
File → Write → Case & Data预览网格运动
Setup → Dynamic Mesh → Preview Mesh Motion... 指定压力-速度耦合方案
设置自动保存
Solution → Calculation Activities → Autosave (可选)设置图形动画
Calculation Activities Task Page 使用“动网格库朗数”指导求解
Moving Mesh Courant Number指示网格运动在一个时间步内扫过的单元数量。 Fluent 提供三类网格更新方法,可组合使用:
⚠️ 注意:平滑、层铺和重构与网格自适应不兼容。
平滑方法通过移动内部节点来“吸收”边界运动,节点数量和连接关系保持不变。启用方式:在 Dynamic Mesh Task Page 中勾选S moothing,然后点击Settings...选择具体方法。
控制方程为扩散方程:
其中 u̅ 为网格位移速度,γ 为扩散系数。边界条件来自用户指定的边界运动,在变形边界上网格运动与边界相切(法向速度分量为零)。
扩散系数的两种形式
基于边界距离(Boundary Distance)
基于单元体积(Cell Volume)
数值求解方法
/define/dynamic-mesh/controls/s moothing-parameters/diffusion-fvm? yes线性求解器选项(AMG Stabilization):
CG(共轭梯度法)
GMRES(广义最小残差法)
BCGSTAB(双共轭梯度稳定法)
适用性:支持所有单元类型,但对多面体不如 RBF 稳健;计算成本高于弹簧平滑,但网格质量更好;平移运动处理优于旋转运动(旋转运动推荐 RBF);与边界层平滑方法、面区域重构方法不兼容;不支持变形共形周期边界或固体移动壁面。
将网格边理想化为相互连接的弹簧网络。边界节点位移产生与位移成正比的力。
k_fac 为弹簧常数因子,范围 0~1。平衡时净力为零,得到迭代方程:
收敛判据:Δx_rms^{m} / Δx_rms^{1} < Convergence Tolerance。
弹簧常数因子的影响:默认值1有一定阻尼,内节点位移可能超出运动边界(产生畸变单元);值0无阻尼,内节点位移更好地跟随边界,保持原始分布。
适用性:适用于运动以平动为主、主要垂直于边界的场景;非四面体/三角形单元区域需谨慎,多面体单元尤其容易产生高畸变 → 推荐 RBF;不支持变形共形周期边界或固体移动壁面。
将网格视为线弹性固体,控制方程:∇·σ=0,σ=2με+λ tr(ε)I,ε=½(∇d̅+(∇d̅)ᵀ)。泊松比 ν = λ/(2(λ+μ)),范围 -1.0~0.5,仅支持常数材料属性。
边界条件:若指定了变形边界的几何(平面、圆柱等),则法向变形为零,切向应力为零;若几何为unspecified(仅此方法和 RBF 支持),则边界也可在法向变形,所有方向应力为零。
适用性:计算成本比扩散平滑更高,但对某些网格和运动能更好地保持质量;支持 2D 三角形/四边形,3D 四面体/六面体/楔形/金字塔;不支持多面体单元、悬挂节点;与边界层平滑、面区域重构不兼容;不支持变形共形周期边界或固体移动壁面。
RBF 求解一个位移场,该位移场在指定运动的边界上满足位移值,然后插值到所有其他节点。
主要优势:
唯一支持
变形边界几何选项:可指定几何(平面、圆柱等)让节点沿几何切向滑动;也可选择unspecified,节点可任意方向移动。对于非共形界面,若两侧界面均为 unspecified 的变形动态区域,Fluent 会保证两个界面保持连接(仅 RBF 支持,且要求界面两侧为同一单元区域)。
高级设置:Relative Convergence Tolerance(控制位移场求解精度);Verbosity=1 在控制台打印绝对容差;Local S moothing:仅对运动边界附近的单元进行平滑,远离边界的网格不受影响,可大幅减少计算量。若启用局部平滑,通常还需要同时启用重构,因为运动被吸收在更小的区域中。边界层单元会在局部平滑之前独立平滑,以保持形状。
注意:RBF 在分布式并行计算中的扩展性不如其他方法;与边界层平滑方法不兼容。
对于扩散、线弹性固体或 RBF 平滑,可以指定从固定的参考位置进行平滑。这对于周期性或准周期性运动(如涡轮机械叶片颤振)非常有用,可使每个循环的网格质量保持一致。启用方式:在 Mesh S moothing Parameters 对话框中勾选S moothing From Reference Position。
注意:若同时启用 Layering 或 Remeshing,或激活了自动自适应判据,则该选项不可用;不要以 legacy 格式保存 case 文件,否则重启后参考位置会不同。
最简单的方法:将每个节点移动到其相邻节点的几何中心。
β 为拉普拉斯节点松弛因子。Fluent 仅在改善网格质量(畸变率降低)时才更新节点位置。可用于变形边界(3D 三角形单元,2D 线性单元),也用于 2.5D 重构方法中。
用于使边界层网格随相邻运动/变形面区域一起变形。启用方式:在 Dynamic Mesh Zones 对话框中为面区域勾选Deform Adjacent Boundary Layer with Zone。工作原理:边界层中每个单元的节点坐标与底层面对应节点使用相同的位移向量。
适用性:支持 3D 楔形/六面体边界层,2D 四边形边界层;仅与弹簧平滑兼容,不可与扩散、线弹性固体或 RBF 平滑同时使用;主要用于平动,大幅旋转会导致边界层单元畸变。
适用于棱柱形网格(3D 多面体/六面体/楔形,2D 四边形)。根据运动边界相邻层的高度来添加或删除网格层。
层分裂条件(膨胀时):hmin> (1+αs)·hideal
层合并条件(压缩时):hmin< αc·hideal
其中 αs为分裂因子,αc为合并因子,hideal为理想单元高度(可常数或随时间和曲轴角变化)。
分裂选项:
高度法(Height Based)
比率法(Ratio Based)
适用条件:运动边界相邻单元必须是棱柱形(即使单元区域内混合其他形状也可以);网格层必须由单侧面区域完全包围(除非使用滑移界面);不能与自适应同时使用(几乎在所有情况下)。
当边界位移相对于局部单元尺寸过大时,仅靠平滑会导致网格质量恶化甚至负体积。Fluent 会聚合不满足畸变率或尺寸判据的单元,并局部重新划分网格。
启用 Unified Remeshing 后,算法结合多种重构方法的优点,并默认尝试保持初始网格尺寸分布,简化设置,提高稳健性(尤其适合并行计算)。
高级设置(点击 Advanced...):
Maximum Cell Skewness
Maximum Face Skewness
Retain Size Distribution
尺寸控制(Sizing Controls)和棱柱控制(Pris m Controls)可通过 TUI 命令创建,用于精细控制网格尺寸分布或边界层参数(第一层高度、增长率、层数等)。在 Watertight Geometry 网格工作流中创建的尺寸/棱柱控制会自动导入统一重构。
若选择 Methods‑Based Remeshing,可单独启用以下子方法:
局部单元重构(Local Cell Remeshing)
局部面重构(Local Face Remeshing)
单元区域重构(Cell Zone Remeshing)
面区域重构(Face Region Remeshing)
2.5D 表面重构(2.5D Surface Remeshing)
用于局部重构时,基于边界单元尺寸生成背景网格,从而控制网格尺寸的空间分布。
背景网格:Fluent 在每个流体区域周围画一个大约两倍于区域尺寸的包围盒,根据最短特征长度和用户指定的 Resolution 细分,创建背景网格。在每个背景网格顶点 I 计算尺寸函数值:
然后通过拉普拉斯平滑,并插值到任意单元形心 P:
参数:
α(Variation)
β(Rate)
d_b
Resolution
启用尺寸函数后,单元被标记为需要重构的条件是:size ∉ [4/5·γ·size_b, 5/4·γ·size_b]。
Fluent 根据启用的模型和单元形状自动决定每个区域使用哪种更新方法:
注意:重构可能导致并行计算中各分区网格略有不同,从而引起非常微小的解差异。若模型包含大量壳传导区域,每次网格更新会删除并重建所有壳,非常耗时。
二阶时间格式可与重构、层铺和平滑一起使用。层铺事件和网格平滑期间保持二阶精度;重构事件发生时,时间推进精度自动降为一阶。若重构非常频繁(每个时间步都发生),实际就是一阶格式。可通过 TUI 命令define/dynamic-mesh/transient-settings/verbosity查看当前使用的时间格式。
在 3D 中,重构和/或平滑可能导致角部被倒角。若某些角部(特定角度范围)对流动很重要,可启用特征检测来保留它们。在 Dynamic Mesh Zones 对话框的 Geometry Definition 选项卡中,勾选 Feature Detection 并设置 Feature Angle:180°(默认)不保留任何角部,0°保留所有角部。适用性:仅与面区域重构配合使用,特征边上的节点不允许平滑。
启用 In‑Cylinder 选项后,需要设置:曲轴转速(Crank Shaft Speed)、起始曲轴角(Starting Crank Angle)、曲轴周期(Crank Period,通常720°四冲程或360°二冲程)、曲轴角步长(Crank Angle Step Size,默认0.5°)。
内置活塞运动函数 piston‑full:活塞位置作为曲轴角的函数
其中 rc曲柄半径,L 连杆长度,x_offset 活塞销偏移(有正负),θc= θs+ t·Ωshaft。
限制参数:Piston Stroke Cutoff(活塞行程截止值,控制层铺何时开始);Minimum Valve Lift(最小气门升程,防止气门与阀座接触时单元退化为扁平)。
气门升程剖面:使用 Profile 文件定义,必须包含 angle 和 lift 字段。示例:
((ex-valve 5 point) (angle 0 180 270 360 720) (lift 0.05 0.05 1.8 0.05 0.05))
滑移接口事件:气门关闭时(升程<最小值)删除滑移接口,并将接口区域转换为壁面;气门打开时重新创建滑移接口。这些通过 Dynamic Mesh Events 实现。
涡流/滚流输出:可启用 Write In‑Cylinder Output,设置输出频率、涡流中心计算方法(重心/固定/UDF)、涡流轴、滚流轴、交叉滚流轴等。
计算作用在刚体上的气动力和力矩(通过压力与剪切应力积分),加上重力、弹簧力、推力等外部载荷,求解刚体运动。适用于外挂物投放、弹射分离等。
启用:Dynamic Mesh Task Page → Options → Six DOF
设置刚体属性:在 Six DOF Properties 对话框中为每个运动物体创建属性集,必须定义:Mass(质量)、Inertia Tensor(惯性张量,相对于重心和刚体方向)。可选用:One DOF Translation(单自由度平动,可定义弹簧刚度、预紧力、运动范围)、One DOF Rotation(单自由度转动,需定义转动惯量)。
运动历史输出:勾选 Write Motion History,为每个运动物体生成历史文件。
配合隐式更新:若六自由度运动依赖于流场,建议启用隐式网格更新。若仍不稳定,可在边界区域的 Solver Options 选项卡中启用 Solution Stabilization。
对于网格运动依赖于流场的应用(如六自由度、流固耦合),隐式更新允许在一个时间步内多次更新网格,基于收敛中的流场解逐步调整网格位置,实现更强的耦合,提高鲁棒性。
启用:Dynamic Mesh Task Page → Options → Implicit Update
设置(Options 对话框 → Implicit Update 选项卡):
Update Interval
Motion Relaxation ω
限制:不能与密度基显式求解器、稳态问题、内燃缸同时使用;使用 UDF 计算运动时,UDF 必须在每次调用时使用当前流场,不能存储之前的信息。
用于检测运动表面是否与其他表面发生接触(距离小于阈值)。接触信息可在 UDF 中用于约束运动,或控制接触区域的流动。
启用:Dynamic Mesh Task Page → Options → Contact Detection
设置:从 Face Zones 列表中选择参与检测的面(仅壁面和 Rigid Body/User‑Defined 动态区域);Proximity Threshold:触发接触的距离阈值;UDF:接触时调用的函数(可选)。
流控方法:
Contact Zones(接触区法)
Contact Marks(接触标记法)
显示接触标记:通过 TUIdefine/dynamic-mesh/controls/contact-parameters/render-contactcells?启用,然后在 Contours 中选择 Cell Info... → Contact Cell Mark 显示。
用于在瞬态模拟的特定时刻(时间或曲轴角)执行预定动作。内燃缸模拟中尤其重要(气门开闭、滑移接口创建/删除)。
设置:Dynamic Mesh Events 对话框。增加 Number of Events,勾选事件,输入名称,指定触发时间(At Time)或曲轴角(Act Crank Angle),点击 Define... 选择事件类型并设置参数。点击 Preview... 可播放事件以验证。
可用事件类型:Change Zone Type、Copy Zone BC、Activate/Deactivate Cell Zone、Create/Delete Sliding Interface、Change Motion Attribute、Change Time Step Size、Change Under‑Relaxation Factors、Insert/Remove Boundary Zone Layer、Insert/Remove Interior Zone Layer、Insert/Remove Cell Layer、Execute Command、Replace Mesh、Inert EGR Reset、Diesel Unsteady Flamelet Reset 等。
内燃缸事件循环:事件在满足 θ_event = θ_c ± n·θ_period 时重复执行。
动网格也可用于稳态问题,例如:在不同阀门位置达到稳态后进行比较,或流固耦合界面达到稳态后更新网格。
与瞬态的主要区别:
不可用
更新方式:GUI:Solution → Run Calculation → Update Dynamic Mesh...;TUI:solve → mesh-motion。
建议:六自由度稳态问题可能需要调整伪时间步长,确保在更新前稳态解已较好收敛。