卡门涡街不收敛?圆柱绕流5步出图

做圆柱绕流,很多人一上来入口全速 1m/s、时间步拉满,结果初始瞬间整个流场震荡、求解器直接不收敛;好不容易跑起来,涡街要么不脱、要么乱成一锅粥。问题还是在瞬态层流的设置细节上。这篇把 COMSOL 里 Re=100 卡门涡街案例拆成 5 步。


1. 先判流态:Re=100 才出卡门涡街

概念上,本模型是流过一个长圆柱体的非定常不可压缩流动,圆柱横放在通道里、略微偏离中心,故意让对称流不稳定。雷诺数小于 100 时流动是定常的;Re 到 100 左右就能形成充分发展的卡门涡街。所以物理场选二维层流(spf)+ 瞬态研究,别选稳态——稳态根本看不到周期性脱涡。

卡门涡街不收敛?圆柱绕流5步出图的图1


2. 参数与入口:平滑阶跃别让速度瞬间拉满

参数给全:通道高 h=0.41m、长 L=2.2m、宽 1m、圆柱半径 r=0.05m、入口速度 v1=1m/s。入口用抛物线剖面:6·v1·y·(h−y)/h²,顶部设压力出口。最关键的收敛技巧:用一个平滑阶跃函数 step1(位置 0.05、过渡区间 0.1)去乘入口速度,让速度从 0 平滑爬升,避免"初始瞬间全速"造成的首步不收敛。

3. 网格:物理场控制选"更细",自动出边界层

网格别手动撒。用物理场控制网格,单元大小选"更细",COMSOL 会自动在圆柱壁面和边界层加密。这个案例自动生成约 24638 个三角形单元 + 1772 个四边形边界层单元,最小质量 0.218。手动网格要么太粗壁面边界层没解出来,要么网格量爆炸。

图片

卡门涡街不收敛?圆柱绕流5步出图的图3


4. 分时步:前期大步、后期小步

瞬态求解最讲究时间步:0~3.4s 步长 0.2、3.5~7s 步长 0.02——前期流场还没发展开大步省时间,后期要捕捉涡街脱泄切小步。求解器时间步选"中级/严格"相对容差,相对容差留 0.005。反例是激光加热这类场景要反过来:一开始变化剧烈用小步、后期平稳放粗。

图片

卡门涡街不收敛?圆柱绕流5步出图的图5


5. 后处理:升阻力系数和中心线监测

跑稳后做三件事:加粒子追踪看涡街脱落形态;定义升/阻力系数,$C_l = \frac{2F_y}{\rho v_1^2 \cdot 2r \cdot w$(w 为展向单位宽度),监测它的周期性振荡;再沿通道中心线画一条节点线监测速度随时间的变化。升阻力出现稳定的周期波形,才算涡街真正发展充分。

小结

Re=100 二维层流瞬态;入口抛物线乘平滑阶跃函数防初始震荡;网格用物理场控制"更细"自动出边界层;分时步前期 0.2、后期 0.02;后处理看升阻力周期和中心线速度。五步走稳,卡门涡街一次成型。


【配套资料领取】

扫码添加小助理,私信回复关键词【圆柱绕流】,领取本篇配套资料包:①Re=100 几何与速度参数表 ②平滑阶跃 step 函数设置方法 ③抛物线入口与压力出口边界设置 ④物理场控制"更细"网格统计(24638 三角形/1772 四边形)⑤分时步策略与求解器容差 ⑥升阻力系数公式与中心线监测步骤。

卡门涡街不收敛?圆柱绕流5步出图的图6
查看全文
我的头像
默认 最新
ansys结构交流群