一維黎曼問題數(shù)值解與計(jì)算程序

1、-A.1-一維 一維 問題數(shù)值解與計(jì)算程序 問題數(shù)值解與計(jì)算程序 Riemann一維 問題，即激波管問題，是一個(gè)典型的一維可壓縮無黏氣體動(dòng)力 Riemann學(xué)問題，并有 解析解。對(duì)它采用二階精度 兩步差分格式進(jìn)行數(shù)值求 MacCormack解。同時(shí)，為了初學(xué)者入門和練習(xí)方便,這里給出了用 語言和 編寫的 C Fortran77計(jì)算一維 問題的計(jì)算程序，供大家學(xué)習(xí)參考。 RiemannA-1 A-1 利用 利用 兩步差分格式求解一維 兩步

2、差分格式求解一維 問題 問題 MacCormack Riemann1.一維 問題 Riemann一維 問題實(shí)際上就是激波管問題。激波管是一根兩端封閉、內(nèi)部充 Riemann滿氣體的直管，如圖 A.1 所示。在直管中由一薄膜將激波管隔開，在薄膜兩側(cè)充有均勻理想氣體（可以是同一種氣體，也可以是不同種氣體） ，薄膜兩側(cè)氣體的壓力、密度不同。當(dāng) 時(shí)， 0 ? t氣體處于靜止?fàn)顟B(tài)。當(dāng) 時(shí)，薄膜 0 t ?瞬時(shí)突然破裂，氣體從高壓端沖向低壓端，同時(shí)

3、在管內(nèi)形成激波、稀疏波和接觸間斷等復(fù)雜波系。2.基本方程組、初始條件和邊界條件設(shè)氣體是理想氣體。一維 問題在數(shù)學(xué)上可以用一維可壓縮無黏氣體 Riemann方程組來描述。在直角坐標(biāo)系下量綱為一的一維 方程組為： Euler Euler(A.1) , 1 1 x t x? ? ? ? ? ? ? ? ? 0 u f其中 (A.2) 2 ,( )uu u pE E p u? ?? ?? ? ? ??

4、 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?u f這里 、 、 、 分別是流體的密度、速度、壓力和單位體積總能。理想氣 ? u p E體狀態(tài)方程：? ? ? ? ? ?2 2 1 1 1 2 p e E u v ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?（A.3）初始條件： ； 。 1 1 1 1, 0, 1 u p ? ? ? ? 2 2 2 0.125, 0, 0.1 u p ?

5、? ? ?圖 A.1 激波管問題示意圖-A.3-的前置人工濾波法能消除激波附近的非物理振蕩，計(jì)算效果很好。從圖 A.2 和圖 A.3 中可以看出通過激波后氣體的密度、壓力和速度都是增加的；在壓力分布中存在第二個(gè)臺(tái)階，表明在這里存在一個(gè)接觸間斷，在接觸間斷兩側(cè)壓力是有間斷的，而密度和速度是相等的。這個(gè)計(jì)算結(jié)果正確地反映了一維 問題的物理特性，并被激波管實(shí)驗(yàn)所驗(yàn)證。 RiemannA-2 A-2 一維 一維 問題數(shù)值計(jì)算源程序 問題數(shù)值計(jì)算

6、源程序 Riemann1. 語言源程序 C// MacCormack1D.cpp : 定義控制臺(tái)應(yīng)用程序的入口點(diǎn)。///*-----------------------------------------------------------------------------------------------------*利用 差分格式求解一維激波管問題（ 語言版本） MacCormack C* ---------

7、----------------------------------------------------------------------------------------------*///#include “stdafx.h“#include #include #include #define GAMA 1.4//氣體常數(shù)#define PI 3.141592654#define L 2.0//計(jì)算區(qū)域圖 A.2 采用 語