$\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方程。第一个方程是动量方程,描述了流体微团速度的变化;第二个是连续性方程,保证了流体的不可压缩性。
$$\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现象”。本工具旨在通过本质理解帮助学习流动现象。实际流动是三维的,例如盖驱动方腔流动在具有深度方向的三维空间中会产生复杂的涡结构(如泰勒-戈特勒涡等)。虽然二维模拟成本较低且有助于理解现象本质,但需谨记实际设计中三维分析不可或缺。