求解欧拉方程的c语言,用有限体积方法求解欧拉方程.doc

PAGE

4 -

实用标准

文档

有限体积法求解二维可压缩Euler方程

——计算流体力学课程大作业

老师： 夏健、刘学强

学生： 徐锡虎

学号： SQ09018013018

日期： 2010年2月5日

目 录

一、内容摘要……………………………………………………(2)

二、流动控制方程………………………………………………(2)

三、有限体积法的空间离散……………………………………(2)

四、人工耗散……………………………………………………(3)

五、时间离散……………………………………………………(4)

六、边界条件……………………………………………………(5)

七、计算结果……………………………………………………(8)

八、结论与展望…………………………………………………(11)

参考文献…………………………………………………………(11)

一、内容摘要

本文通过运用JAMESON有限体积法求解了二维定常和非定常可压缩Euler方程。程序实现语言为C++。其中，使用的网格是三角形非结构网格。在时间推进上使用的是四步龙—库塔推进格式。推进的时间步长取的是当地的时间步长。为了消除迭代误差、round-off等误差，本文采用了添加人工耗散项的办法。另外，本文计算了NACA0012翼型在跨音速下不同迎角的情况，并与fluent软件的计算结果进行了比较，来验证程序的准确性。

二、流动控制方程

守恒形式的Euler方程：

(1)

其中x和y代表笛卡儿坐标系。W是守恒变量。

(2)