附录B一维可压NS方程B.docx附录 B 一维可压缩黏性流动问题的数值解法与计算程序
一维可压缩黏性流动是气体动力学中最经典的黏性流动问题, 对它采用迎风
型 TVD差分算法进行数值求解。同时,为了初学者入门和练****方便 , 这里给出了
由 C
glj
minmod
l
1
l
1
,
l
1
l
1
( )
j
2
j
2
j
2
j
2
glj
1
g lj
,
l
0
l
j
1
l
1
2
()
j
1
2
j
2
0,
l
1
0
j
2
l
R
1
l
u
l
1 ( )
1
1
j
j
j
2
2
2
z ,
z
Q z
1
z2
,
z
,
( )
2
1 Q
l
r
l
2
,
非定常
2
i
1
i
1
l
2
2
( )
i 1
1
2
l
,
定常
Q
1
2
2
i
流通量矢量 f 的非线性 Jiacobian系数矩阵 A = R ΛR-1
为:
0
1
0
A u
f
3
u2
3
u
1 ()
u
2
u a2
a2
2 u3
1
3
2
u2
u
2
1
非线性 Jiacobian系数矩阵 A 的特征值为:
1
0
0
Λ
0
2
0
()
0
0
3
1
u,
2
u
a,
3
u
a
系数矩阵
A
的右特征矢量 R、R-1 为:
非线性 Jiacobian
1
1
1
R
u
u
a
u
a ( )
附录B一维可压NS方程B 来自淘豆网www.taodocs.com转载请标明出处.