五类黎曼间断问题数值模拟题目.docVIP

五类黎曼间断问题数值模拟 题目： 1．对于一维流Euler方程： 试用通量的矩阵分裂格式的思想，构造此方程的数值解法。左右边界不考虑均使用Neumann边界条件： 2．编写MatLab程序计算本章第六节中所述的5类黎曼间断问题中的任意三类问题。 3．模拟初场如下的非定常流，并尝试分析常数k取何数值可使流场最大马赫数大于或小于0.3： 计算过程如下： 1.通量的矩阵分裂格式： 将方程组： 离散后得： 方程空间导数项的差分计算转化为网格单元界面处的通量计算。 而双曲型问题的主要特点是存在特征线，并且扰动以特征速度沿特征线传播。不同的特征分量，传播方向不同。为体现迎风的思想，需要考虑这种方向性。以一维情况为例，由于影响来自上游，对于向右传播的分量，应该采用左边的值计算，对于向左传播的分量，采用右边的值计算。按照这种思想，需要进行特征分裂。 对于方程组，当地线性化后为：。可将写成的形式，其中。方程组左乘，有： 定义，使得：，则有： 将分解成正部和负部，即，其中： ， 对这两种方向不同的特征波，考虑其影响域，可以导出界面处通量的表达式。由于： 也可以写成： 对于特征形式的通量，同样有： 也即有： 界面处的通量为： 也可以直接由左侧或右侧得到： 上式左乘，得： 对从左右两侧得到的通量进行平均，则得： “—”表示代入Roe平均值进行计算。 Roe平均值： ，， ， 前次作业已经导出矩阵的形式，并求出其特征值：，且相应的矩阵 及向量的形式为： 因此，通量的矩阵分裂格式： 2.黎曼间断问题： 一维流动，初始条件：，即所有物理量在t=0时刻的值在x=0处有间断，左右两侧常数分布。这种初始条件的间断一般是不满足间断关系的，因而是不稳定的。在后，这种初始间断立即分解为若干满足间断关系的间断及中心稀疏波或激波。 由于初始情况不同，一般该问题可以分解为下图五种类型中的一种。（虚线表示接触间断，单条线表示激波，一束线表示稀疏波） Riemann问题解的类型的判别规则： 因此对于该问题，用roe格式离散，初始条件按上面的方式给定，边界条件采用Neumann边界条件： 即可编写程序进行计算，现计算出三类结果，罗列如下： （1）左、右都是激波， （2）左稀疏波，右激波， （3）中间间断面，两边稀疏波 ， 上面三幅图中，速度、压力、密度同时出现突变，对应于激波间断；速度、压力、密度同时出现连续变化，对应于稀疏波；而速度、压力不变，密度突变，对应接触间断。 3.带压力脉冲的初场计算： 按条件： 在程序中设定初场，即：，，，在每个时间步计算各节点的马赫数分布，实际计算中发现，在时间推进的过程中，马赫数分布逐渐在两侧出现尖峰，之后尖峰增大，直到达到某个最大值，然后再减小，并且尖峰的位置向两侧传播。因此在每个时间步求出各节点马赫数的极大值，看其是否超过0.3。计算结果如下： CFL=0.5，N=100时，，最大马赫数大于0.3。时，某一时刻的马赫数分布如下，最大马赫数0.30039 经过计算发现，不同的CFL数和节点个数N的设置对计算结果有影响。节点个数增多，局部马赫数大于0.3需要的k值减小；CFL数减小，需要的k值则增大。具体如下： （1）CFL=0.5，N=200时，时即出现了马赫数大于0.3 （2）CFL=0.3，N=100时，时出现了马赫数大于0.3