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

《用有限体积方法求解欧拉方程》由会员分享，可在线阅读，更多相关《用有限体积方法求解欧拉方程(12页珍藏版)》请在人人文库网上搜索。

1、有限体积法求解二维可压缩Euler方程计算流体力学课程大作业老师： 夏健、刘学强 学生： 徐锡虎 学号： SQ 日期： 2010年2月5日 目 录一、内容摘要(2)二、流动控制方程(2)三、有限体积法的空间离散(2)四、人工耗散(3)五、时间离散(4)六、边界条件(5)七、计算结果(8)八、结论与展望(11)参考文献(11)一、内容摘要本文通过运用JAMESON有限体积法求解了二维定常和非定常可压缩Euler方程。程序实现语言为C+。其中，使用的网格是三角形非结构网格。在时间推进上使用的是四步龙库塔推进格式。推进的时间步长取的是当地的时间步长。为了消除迭代误差、round-off等误差，本文采。

2、用了添加人工耗散项的办法。另外，本文计算了NACA0012翼型在跨音速下不同迎角的情况，并与fluent软件的计算结果进行了比较，来验证程序的准确性。二、流动控制方程守恒形式的Euler方程：(1)其中x和y代表笛卡儿坐标系。W是守恒变量。(2)F,G表示通量, (3),P , H和E表示密度，压强，单元总焓和单元总能量。U,V表示笛卡儿坐标系下的速度矢量。这些量由理想气体的单位体积的总能量和总焓相互联系。(4)(5)三、有限体积法的空间离散计算域被划分为互不重叠的单元。在每个单元运用守恒形式的Euler方程。由于每个单元相对于时间都是不变的，所以等式(1)可以写成：(8)其中和S是单元的体积。

3、和边界。W是单元的平均值。在对上述方程进行时间离散前，先对空间进行离散，则方程(6)可以写为：(9)其中表示第k个单元的体积，是第k个单元的守恒变量。表示第k个单元的通量。方程(7)的右边项可以写成：(10)其中 (11)