寻求大密度比、压力比的激波管问题的边界值

onecfd

我采用下面边界值计算了一维Euler方程的Riemann问题，现有人建议我试算大密度比、压力比的激波管问题，但我找不到合适的初值，请大家帮忙提供一下。
我已算的边界值：
QL＝（1，0，2.5），QR＝（0.125，0，0.25）

xtuyyy

把密度比改成1;100
让初始速度为零，
压力连续。
有点象2D Rayleigh-Taylor Instability.哈哈！

abcc

经典算例：（\rho, u, p）=(1, 0, 1000) ,x<=0.1
                        =(1, 0, 0.01), 0.1<x<=0.9
                        =(1, 0, 100),  x>0.9
计算区域[0, 1], 反射边界条件。

abcc

[这个贴子最后由abcc在 2004/06/28 08:38pm 第 1 次编辑]

这是比较难的一个算例，算算看。有很多文献都会计算它，附图供你参考：）[DISABLELBCODE]

abcc

忘记说了。这是计算到0.038s的结果。