基于有限体积法的二维欧拉方程数值计算.docVIP

基于有限体积法的二维欧拉方程数值计算 《计算报告》 【摘要】采有中心格式有限体积法数值求解了二维跨音速欧拉方程。使用贴体网格代替物面周围的连续流动区域，空间离散后得到了每一个网格上的半离散方程；时间推进采用了五步二阶Runge-Kutta格式；采用了当地时间步长法、隐式残值平均和焓阻尼等技术使数值计算加速收敛到定常状态。最后给出了NACA0012翼型绕流的数值计算结果和相应的计算程序。 §1计算方法 控制方程 惯性坐标系下，理想可压流体二维非定常欧拉方程守恒形式 （1） 其中 （2） 上式中,分别代表压强、密度、速度方向分量、速度方向分量、单位质量总能和单位质量总焓；实际求解方程(1)时，四个方程有六个未知量，方程组不封闭。考虑到气体的热力学性质，补充两个方程，一般形式为： （3） （4） 对于完全气体流动，满足下列关系式 （5） （6） 为比热比,对于空气取，这样便构成了封闭的方程组。 有限体积法空间离散 将（1）式在任意的控制体单元（边界和单位外法向矢量分别记为和）上进行积分得 （7） 定义 （8） 则应用高斯定理可得 （9） 上式便为积分形式的欧拉方程，它是一个一般表达式，对任意的控制体都适用，如果是固定的、体积不变化的控制体，有 （10） 将积分形式欧拉方程应用到某一个网格单元得 （11） 应用积分中值定理 或 （12） 为在网格单元上的平均值，并定义在网格中心处，如图1；为的面积，它与时间无关 （13） 令，则可得到方程（11）的半离散方程 （14） 其中表示流出控制体的净流量 （15） 上式右端的各通量项的计算是类似的，例如 （16） 所有带分数下标的物理量均取相邻量的平均值，例如 （17） 下面就均匀直角网格讨论上式中心格式的离散精度。对于均匀网格设方向的步长分别为，且有（如图2）： （18） （19） 代入方程（14）得 （20） 显然，对于均匀直角网格，上述中心平均有限体积格式的空间离散具有二阶精度，对于一般的曲线网格，如果网格足够光滑，即网格的大小和形状没有较突然的变化，该格式具有近似二阶精度。 边界条件 欧拉方程进行空间离散时，另一个重要的方面是各类边界条件的处理，处理不当会引起计算的结果不准确或求解过程不稳定。对于二维流动，需要处理的边界条件有：物面边界条件、远场边界条件和周期性边界条件。 （1）物面边界条件 在进行绕流计算时，物面是引起流动扰动的根源，因而不同的绕流问题呈现出不同的流态的原因，主要在于物面边界的不同。另一方面，只有正确计算物面的压强，才能比较准确的得到浸没于流体中物体的受力特征。总之，物面边界条件的处理在流动计算中是十分重要的。 对于无粘流动，物面边界条件表现为无穿透条件，即在物面上法向速度为零： （21） 将上式代入靠近物面的第一层网格，发现对流项的积分只剩下压强项的积分。因此，在处理物面边界条件时，需要给出物面上的压强。物面压强可由法向动量关系式或采用插值公式得到。数值经验表明，在物面附近网格变化不大的情况下，二者的计算精度相当，而由于进行法向动量关系式计算时需要的内存量较大，故这里采用二次插值得到物面压强，如图3。 (22) （2）远场边界条件 我们在进行控制空间方程离散时，只能在有限的区域内划分网格，故而在外边界处不能直接赋为均匀场。非定常欧拉方程属于双曲型方程，方程外行的特征波应该传向无穷远处而不应该反射回来，因而需要在外部边界处按特征波的传播方向进行处理。 这里采用当地一维Riemann不变量处理远场边界。向内和向外传播的一维Riemann不变量分别为：