8.2 动态网格 运动处理的第二级:流域形状本身在变——活塞下行、阀瓣抬升、结构在气流里颤振——网格必须跟着边界实时变形。本节讲动网格的三件套技术(光顺、重构、层变)各自的修理对象、崩溃机理与配置参数,并给出六自由度耦合与 UDF 运动律的接口位置。这是全书物理专题的技术压轴,也是最容易"算到半夜负体积崩溃"的地方。 网格变形的三种修理术 光顺(Smoothing)。边界动一点,网格点跟着弹性或弹簧式地挪一挪,拓扑不变。适合小幅、光滑的边界运动;变形积累到网格严重畸变时它无能为力——网格还是那批网格,只是被拉歪了。 重构(Remeshing)。畸变越过阈值就把局部单元拆掉重划。它救光顺的命:大位移、大变形场景(阀全开全关、大幅摆动)离不开它。
运动处理的第二级:流域形状本身在变——活塞下行、阀瓣抬升、结构在气流里颤振——网格必须跟着边界实时变形。本节讲动网格的三件套技术(光顺、重构、层变)各自的修理对象、崩溃机理与配置参数,并给出六自由度耦合与 UDF 运动律的接口位置。这是全书物理专题的技术压轴,也是最容易"算到半夜负体积崩溃"的地方。
光顺(Smoothing)。边界动一点,网格点跟着弹性或弹簧式地挪一挪,拓扑不变。适合小幅、光滑的边界运动;变形积累到网格严重畸变时它无能为力——网格还是那批网格,只是被拉歪了。
重构(Remeshing)。畸变越过阈值就把局部单元拆掉重划。它救光顺的命:大位移、大变形场景(阀全开全关、大幅摆动)离不开它。重构引入两类开销:每步重构的 CPU 时间,以及重构带来的插值误差与噪声——重构频率过高说明尺寸参数设得不合理,过低则等待你的就是负体积。
层变(Layering)。边界平移专用:像手风琴一样增减边界层的层数(活塞往复活塞顶网格拉伸压缩,挤掉或新增一层)。它只对近似平面的边界运动有效,但精度干净、无插值污染,缸内问题是它的主场。
三件套的正确用法是组合拳:光顺吃小幅变形、重构处理拓扑变化、层变负责平面推进。配置前先做"运动预演"——把运动律的幅值频率与网格参数干跑几个周期,观察最小网格质量,是动网格项目省一夜算力的最佳投资。
负体积是动网格崩溃的统一症状:单元面积或体积变成负数,方程组在数学上失效,求解器报错退出。病因清单按频率排序:其一,时间步过大,边界一步跨过整个薄网格区——步长必须与最小网格尺寸和运动速度匹配(库朗数思维同样适用于运动);其二,重构参数(最大偏斜度、尺寸范围)设得太松,坏单元没被及时重构;其三,运动律本身有病——UDF 里速度写成了位移、正负号方向颠倒、曲线在拐点处不光滑;其四,预解压没有预演,几何在极限位置本来就留不出变形空间。诊断顺序照此排:先查运动律,再跑预演,最后调步长与重构参数——顺序反了会把模型问题错当数值问题调一宿。
背景。某先导阀需要评估快速关闭过程中的水锤压力峰值,阀瓣行程八毫米、动作时间二十毫秒,流域形状大幅变化——MRF 无能为力,动网格正选。
操作。网格策略:阀口密封区网格大幅加密(这里是压力峰值与变形的冲突区),运动用 UDF 给定阀瓣位移曲线(实测行程-时间数据拟合成光滑函数);动网格三件套全开:光顺用弹性系数法、重构限定阀腔区、层变负责阀杆平移段。时间步按位移速率与最小网格尺度初估后取微秒级,每步内动网格更新与求解交替进行。压力监测点布置在阀口下游二倍管径与管道远端两处。
结果。全行程跑通无崩溃,阀口最大流速与压降出现在行程中段的节流窗口,关闭末段下游监测点捕捉到水锤压力波并以声速沿管道传播,峰值与简化水锤理论公式的量级一致(精细差异来自流道三维效应与压缩性设置)。重构区在最大开度附近触发最频繁,该时段的网格质量直方图仍在合格线内。
解读。案例示范了动网格项目的三个纪律动作:运动预演先行(预演发现了密封区极限位置网格余量不足并提前加密)、压力监测点先布后算(水锤结论的可信度取决于监测点位置的事先论证)、结果用理论量级交叉对表(三维仿真的峰值与一维水锤公式互为校验,这就是第九章确认思想的预演)。
变式。结构变形不是预设而是流场耦合出来的(颤振、柔性阀片),叠加六自由度求解器或流固耦合接口,结构侧刚度阻尼参数进入敏感清单;多个刚体相互接触(气门落座)需要接触模型支持,Fluent 动网格对接触的成熟度有限,必要时考虑专用运动求解器配合网格重叠(overset)技术。
重构参数的阈值怎么定。最大偏斜度重构阈值取合格线上浮一点(比如零点九触发重构),尺寸范围取初始网格尺寸的零点五到一点五倍——阈值是"提前拆房"的触发线,设得太松等于等塌方。
光顺的弹性系数怎么调。边界附近刚度大、远处刚度小是默认思路:收得不顺时降低远场刚度系数让网格"软"下来,收敛慢则反向。光顺参数的调法是手感活,预演环节就是练手感的地方。
层变失效(层数不增减)什么原因。层变对边界平行度要求苛刻:分离面不平行、层高与分裂阈值不匹配都会静默失效。先在预演里可视化层变是否触发,再调分裂与塌缩阈值。
崩溃前兆速查表(按出现顺序):
| 阶段信号 | 含义 | 动作 |
|---|---|---|
| 最小正交质量持续下滑 | 变形积累临近极限 | 收紧重构阈值或减步长 |
| 偏斜度直方图尾部变胖 | 局部网格濒临畸变 | 预演定位运动极限区加密 |
| 迭代内残差突然抬升 | 网格更新扰动求解 | 减步长、查运动律光滑性 |
| 负体积报错退出 | 彻底崩溃 | 回运动律与预演,勿盲目续算 |
关于六自由度接口再多说两句配置逻辑:六自由度求解器在动网格更新时读取当前流场对刚体的合力与合力矩,积分牛顿-欧拉方程得到下一时刻的运动状态,再交给动网格三件套执行变形——它是"流场驱动运动"的闭环。配置上有三件事要一次定对:刚体的质量与转动惯量(惯性参数错,运动形态整体失真);初始自由度的释放(约束不该动的自由度,否则数值噪声会把刚体带漂);以及时间步与刚体固有频率的匹配(步长下探到周期的百分之一量级起步)。预演阶段把流场关掉、只跑运动积分,先确认刚体响应合理,再开流场耦合。