对流扩散方程仿真器 返回
Interactive Tool — CFD Numerics

对流扩散方程仿真器
(数值稳定性可视化)

将中心差分、迎风差分、QUICK法的稳态解与分析解对比。实时观察Peclet数增大时中心差分产生振荡的现象。直观理解CFD数值格式。

流动与扩散参数
对流速度 u
m/s
正值(左→右方向)
扩散系数 D
m²/s
网格间距 Δx
m
网格数 N = 1.0/Δx
稳定 — Pe = 1.00
数值格式
边界条件
φ(0) = 1.0
φ(L) = 0.0
稳态解对比 中心差分
Cd
数值特性摘要
计算结果
Peclet数 Pe = uΔx/D
数值扩散系数 Dnum
最大误差 |φ_num − φ_exact|

理论说明 — 对流扩散与数值稳定性

稳态1维对流扩散方程分析解(φ(0)=1, φ(L)=0):

$$\phi(x) = \frac{e^{Pe \cdot x/L}- e^{Pe}}{1 - e^{Pe}}, \quad Pe_{total}= \frac{uL}{D}$$

中心差分格式在单元Peclet数 $Pe_{cell}= u\Delta x / D > 2$ 时会产生非物理振荡。迎风差分格式无条件稳定,但会引入 $D_{num}= u\Delta x/2$ 的数值扩散。QUICK格式介于两者之间,具有3阶精度,但在某些情况下仍会有振荡残留。

什么是稳态一维对流扩散方程

🙋
“对流扩散方程”听起来好复杂,它到底是什么呀?
🎓
简单来说,它描述了一个量(比如污染物浓度或温度)在“被流体带着跑”(对流)和“自己向四周散开”(扩散)两种作用下的平衡状态。在实际工程中,比如预测烟囱排出的烟雾如何在下风向扩散,就靠这个方程。你试着在模拟器里把“对流速度”滑块调大,看看右侧的曲线怎么变化,就能直观感受到“带着跑”的效果了。
🙋
诶,真的吗?我调大了速度,发现“中心差分”那条线变得歪歪扭扭,像锯齿一样,这是为啥?
🎓
这就是数值计算里一个经典的有趣现象!当对流作用远强于扩散作用时(专业上叫Peclet数很大),中心差分格式会失去稳定性,产生你看到的这些非物理的“数值振荡”,就像计算出来的浓度在上下乱跳,这在实际中是不可能发生的。工程现场常见的是改用“迎风差分”,你切换到它看看,锯齿是不是消失了?
🙋
锯齿是没了,但迎风差分的线好像比“精确解”那条线更平缓一些,这是不是也不对?
🎓
你的观察非常敏锐!这就是“数值耗散”或叫“假扩散”。迎风格式用“上游”的信息来保证稳定,但代价是人为地引入了额外的“抹平”效果,就像给图像加了模糊滤镜。你可以再试试第三个选项“QUICK格式”,它是一种高阶格式,在大部分区域能更贴合精确解,但在某些陡峭变化的地方可能也有小问题。改变这些参数,对比三种线的差异,你就能深刻理解CFD(计算流体力学)中选择格式时的权衡了。

物理模型与关键公式

本模拟器求解的是一维稳态对流扩散方程,它描述了物理量φ(如浓度)在对流和扩散共同作用下达到的平衡状态。

$$u \frac{d\phi}{dx}= D \frac{d^2\phi}{dx^2}$$

其中,$u$ 是对流速度(m/s),$D$ 是扩散系数(m²/s),$\phi$ 是待求的标量(如浓度)。方程左边代表对流项,右边代表扩散项。

问题的解析解(精确解)可以用来检验数值方法的准确性,它依赖于一个关键的无量纲数——佩克莱数(Peclet number)。

$$\phi(x) = \frac{e^{Pe \cdot x/L}- e^{Pe}}{1 - e^{Pe}}, \quad Pe = \frac{uL}{D}$$

这里,$Pe$ 是佩克莱数,表征对流与扩散作用的相对强弱。$L$是域的长度,$x$是位置坐标。$Pe$越大,对流越占主导,解的变化越剧烈,对数值方法挑战越大。

现实世界中的应用

环境工程与污染物扩散:用于预测工厂排放的废气或废水中的污染物在河流、大气中的分布情况,是环境影响评估的核心工具。

汽车与航空航天工程:模拟发动机燃烧室内燃料与空气的混合过程,或者车身表面附面层内热量和动量的输运,对优化设计和提高效率至关重要。

电子设备散热设计:分析芯片发热时,热量在散热鳍片和周围空气中的对流与扩散,帮助设计更高效的冷却系统,防止设备过热。

化工过程与反应器设计:计算反应物在管道或反应器内的浓度分布,确保化学反应能充分、高效地进行,是化工流程模拟的基础。

常见误解与注意事项

模型假设:本模拟器所用数学模型基于线性、均质、各向同性等简化假设。在将计算结果直接用于设计决策之前,务必确认实际系统是否满足这些假设。

单位与量纲:许多计算错误源于单位换算错误或数量级判断失误。请时刻注意各参数输入框旁标注的单位。

结果验证:始终将模拟器输出结果与物理直觉或手算结果进行核对。若结果出乎意料,请检查输入参数或采用独立方法进行验证。

进阶学习指引

深化理论:在本工具的简化模型基础上,进一步研究非线性效应、三维行为和时间依赖现象。阅读专业教材和学术论文,掌握严格的数学推导,是提升工程解题能力的关键。

数值方法:系统学习有限元法(FEM)、有限差分法(FDM)和有限体积法(FVM),理解商业CAE求解器的内部运行机制,这将显著提升您设置有效仿真的能力。

实验验证:理论和仿真结果必须通过实验数据加以验证。养成将计算结果与测量值进行对比的习惯,这正是V&V(验证与确认)的精髓所在。

CAE工具:准备好后,可进一步探索Ansys、Abaqus、OpenFOAM、COMSOL等业界主流工具。通过本模拟器培养的物理直觉,将帮助您更有效地配置和使用这些工具。