计算流体力学CFD二维方腔流数值仿真系统

资 源 简 介

本项目深入展示了计算流体力学（CFD）的基本思想，重点在于将偏微分方程组形式的流体动力学控制方程通过数值手段进行离散与求解。 项目内含一个完整的Matlab编写的二维方腔驱动流（Lid-Driven Cavity Flow）模拟例子，通过编写数值解法程序，实现了对流体流动行为的定量模拟。 其实现过程涵盖了计算区域的网格划分、控制方程（Navier-Stokes方程）的离散化处理、以及压力项与速度项的解耦计算。

详 情 说 明

计算流体力学（CFD）基本思想与Matlab数值仿真系统

本系统是一个基于Matlab开发的计算流体力学（CFD）教学与科研演示工具。其核心目标是展示如何通过数值方法求解不可压缩Navier-Stokes（N-S）方程组，并以经典的二维方腔驱动流（Lid-Driven Cavity Flow）为例，实现了流体运动过程的定量模拟与可视化。

项目介绍

本项目通过自顶向下的设计思路，将复杂的流体动力学控制方程离散化为计算机可处理的代数方程组。系统聚焦于投影法（Projection Method）的应用，实现了压力场与速度场的解耦计算，完整覆盖了前处理（网格划分）、求解器（时间步迭代、SOR算法求解压力泊松方程）及后处理（流函数计算、剖面分析及云图绘制）的流体力学仿真全过程。

功能特性

• 全数值离散求解：采用中心差分格式处理对流项与扩散项，并利用显式投影法进行时间推进。
• 高效压力解法器：内置逐次超松弛（SOR）迭代法求解压力泊松方程，支持自定义松弛因子（omega）以提高收敛速度。
• 物理仿真精度：支持自定义流体密度、动力粘性系数及驱动速度，适用于不同雷诺数（Re）下的流动模拟。
• 流函数后处理：基于旋度场求解流函数方程，生成平滑的流线分布图。
• 多维度可视化：提供速度矢量场、速率云图、压力分布、流函数曲线、中心线速度剖面及残差收敛曲线等六大分析图表。

系统要求

• 软件环境：Matlab R2016b 或更高版本。
• 硬件要求：建议 4GB RAM 以上，计算量取决于网格密度。
• 依赖库：无需额外第三方工具箱，仅需Matlab基础运行环境。

使用方法

1. 在Matlab编辑器中打开主程序文件。
2. 根据需要修改参数设置部分（如网格数 nx/ny、顶盖速度 U_lid 或粘性系数 mu）。
3. 点击运行脚本，系统将开始时间推进计算，并在命令行实时输出迭代步数与残差数值。
4. 计算完成后，系统将自动弹出包含六个子图的仿真结果分析界面。

核心实现逻辑

程序运行遵循标准的CFD数值求解流程：

1. 网格与参数设置：在 0.1m x 0.1m 的物理区域内建立等间距笛卡尔网格（默认 51x51），初始化速度场和压力场矩阵。
2. 中间速度预估：忽略压力梯度项，利用显式离散化的动量方程计算中间速度（u_star, v_star）。其中对流项采用二阶中心差分，扩散项利用拉普拉斯算符的二阶中心差分实现。
3. 压力泊松方程求解：构建压力的源语项（基于中间速度的散度），利用SOR迭代法求解压力泊松方程。在边界处应用Neumann边界条件（边界压力梯度为零），确保流场质量守恒。
4. 速度修正：利用求解出的压力场梯度对中间速度进行修正，得到满足连续性方程的真实速度场。
5. 边界条件应用：在每一步计算中强制执行物理边界条件（顶盖水平移动速度固定，其余三面墙壁为无滑移零速条件）。
6. 收敛判定：通过监控速度场的均方根残差（RMS Residual）来判断流场是否达到稳态，若满足容差或达到最大步数则停止迭代。

关键函数与算法分析

• 投影法（Projection Method）：该算法是处理不可压缩流的核心。它将动量方程分解为两步：第一步计算包含粘性和对流作用的临时速度，第二步通过压力梯度修正速度使之无散。
• SOR超松弛迭代：在求解压力泊松方程时，相比于标准的Jacobi或Gauss-Seidel迭代，SOR通过引入超松弛因子（如1.5）显著加快了代数方程组的收敛速度。
• 流函数方程（Stream Function Solver）：在得到最终速度场后，程序通过计算速度的旋度得到涡量场，并再次利用泊松迭代求解流函数，从而精准地描绘出流体的运动轨迹和旋涡位置。
• 后处理可视化逻辑：系统利用 contourf 绘制填充云图展示速度和压力分布；利用 quiver 函数叠加矢量场展示流动趋势；通过提取中心线上的速度数组展示边界层效应，这与经典CFD文献中的基准测试保持一致。