二维管道流动和开口方腔流动有限体积法MATLAB程序
本案例为典型的稳态二维管道和开口方腔流道有限体积法MATLAB程序,编程思路以及相关算法依据《有限体积法基础(第二版)》或者英文版《An Introduction to Computational Fluid Dynamics》。压力速度耦合采用SIMPLE算法,迭代计算采用TDMA算法。每个算例包含一个主程序“SIMPLE_cavity”和“SIMPLE_channel”和TDMA函数“TDMA_pressure”和 “TDMA_velocity”。主程序定义了计算域、物性参数、收敛标准、初始化、有限体积法基本系数计算、后处理显示等内容。两个TDMA函数用于计算SIMPLE算法的过程量(速度、压力修正)。程序完全原创编写,可读性很强,每一个知识点都可以对照前文所述教材,有需要的读者可以付费获取。另外,本案例程序旨在帮助读者理解有限体积法和SIMPLE算法,读者也可以在该程序基础上进行优化。下面是管道流和开口方腔流的计算结果。图1 管道流迭代曲线图2 管道流速度压力云图图3 管道流数值计算和理论结果对比图4 开口方腔流迭代曲线图5 开口方腔流速度压力云图、速度矢量以下是程序代码和源文件索取路径,有偿提供。程序部分内容截图如下来源:仿真与工程