CFD学习笔记
文章平均质量分 97
Patrick_PS
这个作者很懒,什么都没留下…
展开
-
流过平板的超声速流动的CFD计算(附完整代码)
入门CFD,主要参考书目《计算流体力学基础及其应用》(John D.Anderson著,吴颂平等译)实现了第 10章流过平板的超声速流动的代码。利用有限差分法求解二维Navier-Stokes方程,采用MacCormack显示方法。代码总体不难,按照作者提供的思路编写即可。最后得到的结果中,压力一项与作者给出的有出入,在平板后缘靠近壁面的地方有震荡,猜测与边界条件有关:书中给出的壁面边界条件是对压力进行插值计算,而代码中选择了水平速度,垂直速度,密度和温度作为独立变量进行计算,在壁面处...原创 2020-08-19 22:33:32 · 6422 阅读 · 2 评论 -
不可压库艾特流的压力修正法求解(附完整代码)
入门CFD,主要参考书目《计算流体力学基础及其应用》(John D.Anderson著,吴颂平等译)实现了第 9.4节另一种数值方法:压力修正法 的代码。该方法通过采用三套不同的网格,分别在不同的网格上计算压力,速度水平分量,速度垂直分量,从而避免出现所求物理量出现“棋盘式分布”。因为采用了交错网格,记录起来比较麻烦,所以采用了Python中的3阶ndarray来表示网格组,具体说明如下。1、网格的示意图如下所示,其中黑点代表“压力(p)网格”,红叉代表“垂直速度(v)网格”,蓝点代表“...原创 2020-08-15 22:06:20 · 2139 阅读 · 1 评论 -
不可压缩库艾特流的隐式求解(附完整代码)
入门CFD,主要参考书目《计算流体力学基础及其应用》(John D.Anderson著,吴颂平等译)实现了第 9.3节数值方法:隐式Crank-Nicolson 的代码,问题比较简单,主要是托马斯算法的实现,原书中没有提到对角占优的问题,该算法具体可参考10分钟理解托马斯算法(tridiagonal matrix algorithm,Thomas algorithm)。代码运行结果如下: 不足之处,欢...原创 2020-08-13 03:17:48 · 2273 阅读 · 1 评论 -
二维超声速流——普朗特-迈耶稀疏波的流场CFD解(附完整代码)
入门CFD,主要参考书目《计算流体力学基础及其应用》(John D.Anderson著,吴颂平等译)实现了第 8.3节普朗特-迈耶稀疏波流场的数值解 的代码,采用的是MacCormack 方法,对守恒型方程求解。注意:代码最后运行的结果与作者在书中所提供的并不完全相同:1、当采用作者在书中提供的思路,即在流场下边界壁面处均采用向前差分计算时,结果是在壁面处速度大于稀疏波后速度,稀疏波后速度约710m/s,近壁面速度可达720m/s以上,而作者在书中指出近壁面速度约为706-70...原创 2020-08-12 05:12:53 · 3966 阅读 · 11 评论 -
喷管流动的守恒型CFD解法及激波捕捉(附完整代码)
入门CFD,主要参考书目《计算流体力学基础及其应用》(John D.Anderson著,吴颂平等译)实现了第 7.6节激波捕捉 的代码,采用的是MacCormack 方法,对守恒型方程求解;关于非守恒型方程,可见亚声速-超声速等熵喷管拟一维流动的CFD解法(附完整代码)。代码中增加了求解析解的功能。对求解析解过程的理解和思考如下:1、激波的存在可看做间断点,将激波前后的流动分别看做是两个等熵流动。如下图所示(图片截取自参考书,下标带0的压力表示总压,亦即初始压力)。激波前的等熵流动...原创 2020-08-09 03:17:57 · 4095 阅读 · 0 评论 -
亚声速-超声速等熵喷管拟一维流动的CFD解法(附完整代码)
入门CFD,主要参考书目《计算流体力学基础及其应用》(John D.Anderson著,吴颂平等译)实现了第 7.3节亚声速-超声速等熵喷管流动的CFD解法的代码,采用的是MacCormack 方法。代码中增加了动态显示无量纲温度和无量纲密度的功能,参考的是python中plot实现即时数据动态显示方法,最终结果如下:图中红点代表喉口所在位置。 ...原创 2020-08-07 02:04:06 · 4041 阅读 · 0 评论