2D Navier-Stokes 流动模拟器 返回
流体分析

2D Navier-Stokes 流动模拟器

用涡度流函数法对盖驱动腔流进行实时数值分析。改变雷诺数,直观学习CFD基础。

参数设置
雷诺数 Re
网格大小 N
时间步长速度
统计
计算结果
400
Re
0
迭代次数
0.00
最大涡度
0.00
最小涡度
5
steps/frame
涡度等值线(红=正旋转、蓝=负旋转) + 流线
负旋转
正旋转
黑色线条:流线(流函数等值线)
理论·主要公式

涡度输运方程:

$$\frac{\partial \omega}{\partial t}+ u\frac{\partial \omega}{\partial x}+ v\frac{\partial \omega}{\partial y}= \nu \nabla^2 \omega$$

泊松方程:

$$\nabla^2 \psi = -\omega$$

速度—流函数关系:

$$u = \frac{\partial \psi}{\partial y},\quad v = -\frac{\partial \psi}{\partial x}$$

2D Navier-Stokes 流动模拟器概述

🙋
盖驱动腔流是什么?教科书里经常看到这个。
🎓
对,四周是墙的盒子,只有上边界(盖子)运动,这是CFD最基础的基准问题。简单来说,盖子推动盒内流体,产生复杂的涡,我们研究这个涡如何变化。用这个模拟器,你可以拖动"雷诺数"滑块,实时看到涡如何变化。
🙋
改变雷诺数涡会怎样变?
🎓
大不一样。Re=100时盒中心有一个大涡。但Re升到1000左右时,中心的涡会向左上方移动,左下右下角还会出现小涡,叫"二次涡"。试试用模拟器调节参数,观察一下。
🙋
为什么用"涡度"和"流函数"而不是速度和压力?
🎓
问得好。这叫"涡度-流函数法",是2D CFD的经典手法。最大优势是:连续方程自动满足,不用算压力,计算简洁稳定。所以教学工具特别适合用这个方法。试试改变格子数N看看计算细度的影响。

常见问题

雷诺数(Re)表示移流和扩散强度的比率。Re小时粘性主导,流动稳定层流;Re大时容易产生涡,流场复杂。腔内主涡的位置和强度都会变。试试用滑块改变Re,对比观察。
主要原因是时间步长太大。可以减小Δt(时间步长参数),或检查CFL条件是否满足。另外,Re过大(如>10000)数值不稳定,先用Re=100~1000试试。
本模拟器基于2D、非压缩、涡度-流函数法的层流假设,是盖驱动腔的标准基准问题。低中等Re(<5000)与实验和高精度计算吻合很好;高Re时忽略了3D和湍流效应,主要用于学习流动特征。
右边或上方的设置面板有显示选项复选框,选中/取消"速度矢量""流线""涡度等值线"等即可切换。也可调整色彩类型和透明度。如果显示重叠看不清,关闭不需要的项会清楚一些。

实际应用

CFD软件验证·基准测试:商业软件(ANSYS Fluent、西门子Star-CCM+)或开源软件(OpenFOAM)开发和导入时,第一步必须计算这个盖驱动腔流确认代码正确。Ghia等人(1982)的数据是全球标准。

混合·搅拌工艺基础研究:用来理解槽内搅拌叶片搅动流体的行为。把盖的运动看作搅拌叶,通过模拟预测什么样的涡能有效混合。

微流体芯片设计:生物芯片(Lab-on-a-chip)等微小流路内的流动常具有2D性质。药液混合、粒子分离设计时常用这类2D NS求解器做初步检验。

建筑·环境风工程简易评估:建筑周边风的流动或密闭空间内气流可用2D模型简化理解。比如确定机房内热气滞留的位置(二次涡发生的角部)。

常见误区和注意事项

用这个工具时容易产生几个误解。首先,"Re越高越接近真实流体"——错。Re升高确实导致涡复杂化,但那是"近似非定常乱流行为",而这个模拟是2D层流,和3D真实乱流本质不同。比如,Re=10000的结果直接用于泵设计很危险。这工具是看流动"趋势"的。

其次,网格大小N的设置。N=32或64计算快,但高Re时涡形粗糙,有时还不稳定振荡。这是数值耗散强,掩盖了真实物理粘性。实务中"网格相关性研究"是标准,即反复加密网格直到结果不依赖网格。在这个工具里,试试Re=1000时分别用N=64、128、256计算左下二次涡,就能体会网格的重要性。

最后,"没有压力计算"的含义要理解。涡度-流函数法计算稳定高效,但翼的升阻力、管道压力损失等"压力本身是关键"的现象无法直接求。那需要用"原始变量法"(MAC法、SIMPLE法等)。这工具是理解流动"图案"和"结构"的第一步。

使用指南

  1. 在1~10000范围内调节雷诺数(Re)滑块,设定流体粘性与惯性力比率
  2. 选择格子分割数(n)在16~256间,离散化计算区域并决定空间分辨率
  3. 指定时间步长(dt)在0.001~0.1秒范围内,保证CFD稳定性,点击运行按钮启动盖驱动腔流分析

具体计算示例

Re=1000、格子数n=128的计算需约500步收敛。盖速1.0m/s、腔体1m×1m时,中心涡度振荡范围0.8~1.2rad/s。Re=5000时隅部产生二次涡,最大涡度>3.5rad/s。设定dt=0.01s保持数值稳定。

实务注意事项