参数设置
Navier-Stokes方程
$$\frac{\partial \vec{u}}{\partial t}+ (\vec{u}\cdot\nabla)\vec{u}= -\frac{1}{\rho}\nabla p + \nu\nabla^2\vec{u}$$结果 1
—
结果 2
—
计算结果
用有限差分法对二维不可压缩Navier-Stokes方程进行数值分析。实时可视化方腔流动和后向台阶流动。
计算结果
本模拟器求解的是二维、不可压缩流体的Navier-Stokes方程，它由动量方程和连续性方程组成。动量方程描述了流体微团的运动规律：
$$\frac{\partial \vec{u}}{\partial t}+ (\vec{u}\cdot\nabla)\vec{u}= -\frac{1}{\rho}\nabla p + \nu\nabla^2\vec{u}$$其中，$\vec{u}$是速度矢量，$t$是时间，$p$是压力，$\rho$是密度（常数），$\nu$是运动粘度。等式左边第一项是局部加速度，第二项是对流项；右边第一项是压力梯度力，第二项是粘性扩散力。
为了确保流体不可压缩，即流体微团体积不变，必须同时满足连续性方程（质量守恒）：
$$\nabla \cdot \vec{u} = 0$$这个方程意味着在任何一个点上，流进去的速度和流出来的速度必须相等，没有“空隙”产生。在模拟中，压力$p$的主要作用就是调整速度场，使其时刻满足这个约束条件。
航空航天设计：用于模拟飞机机翼周围的空气流动，分析升力、阻力和可能发生的失速现象，是优化机翼形状、减少油耗的关键步骤。
汽车空气动力学：分析汽车高速行驶时的外部流场，用于设计车身造型以减少风阻和风噪，同时确保发动机冷却和刹车系统的散热气流畅通。
建筑与环境风工程：预测高楼、桥梁周围的风场，评估风压分布和风致振动，防止出现令人不适的强风区或导致结构破坏的涡激振动。
电子设备散热：模拟芯片、散热片周围的气流和热对流，优化风扇位置和散热器鳍片设计，确保电子设备在安全温度下稳定运行。
开始使用这个模拟器时，有几个容易误解的地方。首先，人们常认为“雷诺数越高，越接近真实流动”，但这是不对的。提高Re确实会使涡流变得更复杂，但这只是在模拟“非定常且接近湍流的行为”。然而，这个模拟是二维的，并且是基于“层流”假设进行计算的。它与真实世界的三维湍流有根本区别，因此，例如将Re=10000的结果直接用于泵设计是危险的。请记住，这只是一个用于观察现象“趋势”的工具。
其次，关于网格尺寸N的设置。N=32或64时计算速度快，但在高雷诺数下，涡流形状可能会显得粗糙，有时甚至不稳定并出现振荡。这是因为数值扩散效应过强，抵消了原本的物理粘性效应。在实际工程中，网格分辨率也是模拟的生命线，必须进行“网格依赖性研究”，即不断细化网格直到结果不再依赖于网格的粗糙度。在这个工具中，例如，你可以尝试在Re=1000时，将N分别设为64、128、256，并观察左下角二次涡流的呈现方式，从而亲身体验其重要性。
最后，要理解“未计算压力”的含义。涡量-流函数法计算稳定且效果出色，但它不能直接用于那些“需要了解压力本身的现象”，例如机翼的升阻力或管道内的压力损失。这时，就需要使用直接求解速度和压力的“原始变量法”（如MAC法、SIMPLE法等）。这个工具是理解流动“模式”和“结构”的第一步。
这个盖驱动方腔流动的计算，看似简单，实际上与许多工程领域的基础技术直接相关。首先要提到的是电子设备冷却设计。在优化机箱内部气流或散热片周围的传热时，关于“封闭区域内循环流动”的知识非常有用。可以想象，方腔盖板的运动对应着风扇，加热区域对应着发热元件。
另一个领域是喷涂和涂层工艺。当通过狭缝挤出涂料并涂覆到基板上时，涂料的粘度（即实质上的雷诺数）和模具形状决定了能否形成均匀的膜层。在这个模拟器中，观察高粘度（低Re）时流动迟缓、低粘度（高Re）时产生复杂涡流的现象，就能体会到工艺控制的难点。
此外，在生物工程领域也有应用。例如，在分析人工心脏内的血液流动或动脉瘤血管内的循环（淤滞）时，方腔流动的知识常被用作第一近似模型。特别是在角落产生的二次涡流，会形成“清洗不良区域”，这为预测易形成血栓的位置提供了线索。
熟悉这个模拟器后，下一步可以深入探究“数值求解算法”本身。这个工具内部的计算大致分为两个阶段。首先，使用时间显式方法（例如欧拉法）推进涡量输运方程。接着，为了从得到的涡量求解泊松方程 $\nabla^2 \psi = -\omega$，使用了迭代法（如高斯-赛德尔法或SOR法）。增大N会导致计算变慢，正是因为这种迭代计算的次数增加了。
在数学上，从纳维-斯托克斯方程中消去压力项推导出涡量方程的过程中，向量恒等式（例如 $\nabla \times (\nabla \phi) = 0$ ）是关键。此外，求解泊松方程在数学形式上与电磁学中的静电势计算或结构力学中的挠度计算完全相同。了解“不同领域出现相同方程”这一事实，会极大地拓宽你对工程学的视野。
作为下一步推荐的主题是“利达福方腔流动”。这是一个更通用的模型，不仅顶盖运动，侧面也可能运动，可以观察到更复杂的涡流演变。或者，以这个2D知识为基础，进一步学习3D方腔流动，或者物体绕流的基础——“圆柱绕流”，你将遇到卡门涡街等更丰富的现象。无论如何，通过这个工具所获得的“改变参数如何影响物理现象”的直观感受，在未来使用更高级的CFD软件时，必定会成为你的宝贵能力。