|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
根据第十章自己写的Fortran源代码
主程序名:main.f
子程序名:bc.f---------------------边界条件
conver.f-----------------迭代收敛判断
dynvis.f-----------------计算粘性项
mac.f--------------------MacCormack求解器
mdot.f-------------------质量守恒判断
output.f-----------------结果输出
qcx.f,qcy.f--------------计算导热项
tauxx.f,tauyy.f,tauxy.f--计算粘性项
thermc.f-----------------计算导热系数
tstep.f------------------计算时间步长
解决问题:平面上的超声速流(有粘性边界层和激波)
计算方法:守恒型控制方程,显式MacCormack技术的Time Marching算法。
评述:为该书介绍的最后一个也是最复杂的一个流动类型,程序求解完全的二维N-S方程,考虑了粘性和导热。当来流Mach数为4时,计算结果与书上的结果一致,然而,我的程序不能计算Mach数为25的情况(即书中提供的另一种工况的比较结果),经过仔细检查,我认为自己的程序没毛病,可能是书中求解步长的公式给错了(eq.(10.16) in p.457),还望大家帮助检查一下。最后还要说明的是,该章在指导你写程序时,书上有些关键地方有印刷错误,会误导你,一定要小心。(dgvehicle@yahoo.com)
参考文献:Chapter 10, Computational Fluid Dynamics, The Basics with Applications, John D. Anderson, JR. McGraw-Hill, 2002, 4
最后的评述:
到此为止,我把根据该书写的所有程序都奉献给了大家,发布的一个重要目的是与大家交流CFD和编程技术,因为我的程序还会有许多bug,请不吝指出。如有任何问题或交流意向,请发信给我(dgvehicle@yahoo.com)。
发布这些程序的第二个目的是给大家一个建议,即,最好亲自阅读这本书,并在它的指导下,自己写这些程序(尽管可能你对这些流动类型或Fortran编程比较熟练,最好还是花时间做一做),当然,其中一些小问题可能会“憋”你好几天,希望我的这些程序在你碰到麻烦的时候会有所帮助。
再次强调,自己动手编程非常重要,许多人对SIMPLE算法很熟,试问,又有几个人自己一行行的写过这个程序?所以要自己动手,这样,你会发现收获更大。尽管会走些弯路,但最终你会笑傲江湖的。
CFD既是科学,也是艺术,尽量去享受它吧。 |
|