附录B一维可压N-S方程(B).doc

word
word
1 / 20
word
附录B 一维可压缩黏性流动问题的数值解法与计算程序
一维可压缩黏性流动是气体动力学中最经典的黏性流动问题，对它采用迎风型差分算法进展数值求解。同时，为了初学者入门和练****方便,这里给出了由语言和语言编写的、计算一维可压缩黏性流动问题的计算程序，供大家学****参考。
B-1 利用二阶迎风型差分格式求解一维可压缩黏性流动问题
在两端开口管道中充满了可压缩黏性流体。当黏性流体以超声速从左向右运动时，一定会在管道中形成一道正激波，。和分别为激波波前和波后的参数。该问题可简化为一维可压缩黏性流动问题。当数值解达到稳定时，在管道中可求解得到一道稳定的激波。
、初始条件和边界条件
设流体是黏性流体。一维可压缩黏性流动问题，在数学上可以用一维可压缩黏性流动方程组来描述。量纲为一的一维方程组为：
〔〕
其中
〔〕
word
word
2 / 20
word
〔〕
其中和分别是量纲为一的密度、速度、压力和单位体积总能，为流体的黏性项。为普朗特数〔此处公式中是有量纲量〕，为雷诺数，为比定容热容，为比定压热容，是量纲为一的量，称为气体绝热指数，为当地声速。
求解区域为。取。
初始条件：在时刻，，其他物理量采用线性插值得到。
边界条件：左边界处： 〔〕
右边界处：
〔〕
差分格式
一维方程组中的对流项采用的二阶精度迎风型差分格式：
〔〕
〔〕
其中。向量在第个特征方向上分量为：
word
word
4 / 20
word
〔〕
〔〕
〔〕
〔〕
〔〕
〔〕
流通量矢量的非线性系数矩阵为：
〔〕
非线性系数矩阵的特征值为：
〔〕
word
word
5 / 20
word
非线性系数矩阵的右特征矢量为：
〔〕
〔〕
一维方程组中的黏性项采用二阶精度中心差分格式。
采用用语言和语言对一维可压缩黏性流动问题编制了计算程序，并对雷诺数的流动进展了计算， 。
图 语言程序得到的结果
图 语言程序得到的结果
、速度、压力和单位质量内能的分布。由上述计算结果中可以看出，采用二阶精度迎风型差分格式计算一维可压缩黏性流动问题得到的数值解和经典文献中的结果是完全一致的。计算结果明确，迎风型差分格式能够准确地捕捉激波连续，计算效果较好。由于本问题中黏性较大
word
word
5 / 20
word
，所以计算得到的激波比拟光滑，有一定的宽度。一维可压缩黏性流动问题的解是连续、光滑的。
B-2 一维可压缩黏性流动问题的数值计算源程序

//------------------------------------------------------------------------
// 二阶迎风型差分格式求解一维可压缩黏性流动问题〔语言版本〕
//-------------------------------------------------------------------------
#include ""
#include ""
#define im 201 //网格数
//全局变量：
double Q[3][im],Qold[3][im] ; //Q: [rou, rou*u, E]
double rou[im],u[im],p[im],T[im],E[im],a[im] ;
double Pr,Re,cv,cp,Ma,dx,dt ;
//-------------------------------------------------------------------
void initial()
{
double xl,xr,x;
double ul,Tl,ur,Tr;//进出口的u，T值
int i;
dx=/(im-1) ;
dt=-6 ;
Pr= ;
Re= ;
Ma=