Navier-Stokes方程流场可视化 返回
流体力学 / CFD

Navier-Stokes方程流场可视化

在浏览器中实时模拟二维不可压缩流动。改变雷诺数和流动类型,通过颜色映射直观可视化速度场和涡量场。

Preset
参数
統計
计算结果
最大速度
0.00
U_lid
最大渦度
0.0
步数
0
Re (実効)
100
Sim
初期化中...
理论与主要公式

$\frac{\partial \mathbf{u}}{\partial t}+ (\mathbf{u}\cdot\nabla)\mathbf{u}= -\nabla p + \frac{1}{Re}\nabla^2\mathbf{u}$
$\nabla \cdot \mathbf{u} = 0$
渦度: $\omega = \partial v/\partial x - \partial u/\partial y$

什么是Navier-Stokes方程与盖驱动方腔流

🙋
“盖驱动方腔流”是什么?听起来好复杂。
🎓
简单来说,就是在一个方盒子里研究水流。盒子顶部像盖子一样滑动,带动里面的水流动。这是流体力学和CFD(计算流体力学)里最经典的测试案例之一。你可以在模拟器里选择“Lid-Driven Cavity”模式,直接看到这个方盒子里的流动动画。
🙋
诶,真的吗?那旁边的“雷诺数”滑块是干嘛的?我拖动它会发生什么?
🎓
雷诺数(Re)是惯性力和粘性力的比值。你可以把它想象成“流动的混乱程度”。试着把滑块从100慢慢拖到1000,你会看到流线图发生巨大变化。在Re=100时,流体粘性主导,盒子中心只有一个大而稳定的漩涡;但当Re=1000时,惯性力变强,盒子底部角落会生成新的、更小更复杂的漩涡结构。
🙋
好神奇!那另一个“障碍物绕流”模式呢?和这个有什么不同?
🎓
在实际工程中,比如汽车或桥墩的设计,我们更关心流体绕过物体时的状态。切换到“Flow Past Obstacle”模式,你会看到流体从左向右流动,遇到中间一个方形障碍物。改变雷诺数,观察障碍物后方“尾流”的变化:Re低时,流线平滑地分开再汇合;Re高时,障碍物后方会产生交替脱落的漩涡,也就是著名的“卡门涡街”,这在工程现场是导致振动和噪音的常见原因。

物理模型与关键公式

本模拟器求解的是描述粘性牛顿流体运动的控制方程——非定常、不可压缩Navier-Stokes方程。第一个方程是动量方程,描述了流体微团速度的变化;第二个是连续性方程,保证了流体的不可压缩性。

$$\frac{\partial \mathbf{u}}{\partial t}+ (\mathbf{u}\cdot\nabla)\mathbf{u}= -\nabla p + \frac{1}{Re}\nabla^2\mathbf{u}$$ $$\nabla \cdot \mathbf{u}= 0$$

其中,$\mathbf{u}$是速度矢量,$p$是压力,$t$是时间,$Re$是雷诺数。$(\mathbf{u}\cdot\nabla)\mathbf{u}$是非线性对流项,是流动复杂性的来源;$\frac{1}{Re}\nabla^2\mathbf{u}$是粘性扩散项。

在2D模拟中,我们经常使用“涡量”来直观观察漩涡结构。涡量是流体局部旋转强度的度量,其定义如下:

$$\omega = \frac{\partial v}{\partial x}- \frac{\partial u}{\partial y}$$

这里,$u$和$v$分别是速度在x和y方向的分量。涡量$\omega$为正表示逆时针旋转,为负表示顺时针旋转。模拟器中的彩色云图(例如从蓝色到红色)通常就表示涡量的大小和方向,让你一眼就能看出哪里是漩涡的核心。

现实世界中的应用

汽车空气动力学:在汽车设计初期,工程师使用CFD模拟气流绕过车身的状况,分析阻力系数、升力以及后视镜、A柱等部位产生的涡流,从而优化外形以降低油耗和风噪。模拟器中的“障碍物绕流”就是对此类问题的极度简化。

建筑与结构风工程:设计高层建筑、大桥或冷却塔时,必须研究风绕过结构时产生的涡旋脱落频率,避免与结构固有频率发生共振(即“颤振”),否则可能导致像1940年塔科马海峡大桥倒塌那样的灾难性后果。

电子设备散热:手机、电脑芯片的散热风扇和风道设计,本质上也是在一个密闭腔体内引导空气流动。这与“盖驱动方腔流”的物理模型高度相关,工程师通过模拟优化气流路径,防止局部过热。

化工与混合过程:在大型搅拌罐或反应釜中,通过旋转的桨叶(类似于移动的“盖”)来混合流体。模拟不同雷诺数(对应不同转速或流体粘度)下的流场,对于确保反应物均匀混合、提高生产效率至关重要。

常见误解与注意事项

开始使用本模拟器时,有几个需要注意的要点。首先,请勿误认为“色度图数值直接对应物理压力”。此处显示的是“涡量”,代表流体的“旋转强度”。压力是另一个物理量,本工具并未将其可视化。例如,障碍物后方蓝红相邻的区域,表示方向相反的涡流正在相互抗衡。

其次,请勿过度极端地调整参数。特别是若将雷诺数(Re)一次性设为10000等极大值,可能导致计算失稳、出现网格状噪声(数值振荡)甚至发散。在实际CFD应用中,也不会直接使用最终条件计算,而是从较低Re值逐步提升以确认解的收敛性。例如,若想观察盖驱动方腔在Re=10000时的流动特性,建议按1000→3000→5000…的顺序逐步提升参数。

最后,切勿认为“2D结果能完全代表3D现象”。本工具旨在通过本质理解帮助学习流动现象。实际流动是三维的,例如盖驱动方腔流动在具有深度方向的三维空间中会产生复杂的涡结构(如泰勒-戈特勒涡等)。虽然二维模拟成本较低且有助于理解现象本质,但需谨记实际设计中三维分析不可或缺。