注:本文为贴出各种流通矢量分裂方法的方法/格式,随时记录,以便需要用时随时查看
流通矢量分裂
1. Steger-Warming 分裂
x
,
y
x,y
x,y方向的特征值分别为
λ
1
=
u
,
λ
2
=
u
,
λ
3
=
u
−
a
,
λ
4
=
u
+
a
μ
1
=
v
,
μ
2
=
v
,
μ
3
=
v
−
a
,
μ
4
=
v
+
a
\lambda_1=u,\lambda_2 = u,\lambda_3 = u-a,\lambda_4=u+a \\ \mu_1=v,\mu_2 = v,\mu_3 = v-a,\mu_4=v+a
λ1=u,λ2=u,λ3=u−a,λ4=u+aμ1=v,μ2=v,μ3=v−a,μ4=v+a
使用Steger-Warming分裂
x
x
x方向分裂后的流通矢量为
f
±
(
u
)
=
ρ
2
γ
[
2
(
γ
−
1
)
λ
1
±
+
λ
3
±
+
λ
4
±
2
(
γ
−
1
)
u
λ
1
±
+
(
u
−
a
)
λ
3
±
+
(
u
+
a
)
λ
4
±
2
(
γ
−
1
)
v
λ
1
±
+
v
λ
3
±
+
v
λ
4
±
(
γ
−
1
)
(
u
2
+
v
2
)
λ
1
±
+
(
h
−
a
u
)
λ
3
±
+
(
h
+
a
u
)
λ
4
±
]
\pmb{f}^\plusmn(\pmb{u}) = \frac{\rho}{2\gamma}\begin{bmatrix} 2(\gamma-1)\lambda_1^\plusmn+\lambda_3^\plusmn+\lambda_4^\plusmn \\ 2(\gamma-1)u\lambda_1^\plusmn+(u-a)\lambda_3^\plusmn+(u+a)\lambda_4^\plusmn \\ 2(\gamma-1)v\lambda_1^\plusmn+v\lambda_3^\plusmn+v\lambda_4^\plusmn \\ (\gamma-1)(u^2+v^2)\lambda_1^\plusmn+(h-au)\lambda_3^\plusmn+(h+au)\lambda_4^\plusmn \end{bmatrix}
f±(u)=2γρ
2(γ−1)λ1±+λ3±+λ4±2(γ−1)uλ1±+(u−a)λ3±+(u+a)λ4±2(γ−1)vλ1±+vλ3±+vλ4±(γ−1)(u2+v2)λ1±+(h−au)λ3±+(h+au)λ4±
y
y
y方向分裂后的流通量矢量为’’
g
±
(
u
)
=
ρ
2
γ
[
2
(
γ
−
1
)
μ
1
±
+
μ
3
±
+
μ
4
±
2
(
γ
−
1
)
u
μ
1
±
+
u
μ
3
±
+
u
μ
4
±
2
(
γ
−
1
)
v
μ
1
±
+
(
v
−
a
)
μ
3
±
+
(
v
+
a
)
μ
4
±
(
γ
−
1
)
(
u
2
+
v
2
)
μ
1
±
+
(
h
−
a
v
)
μ
3
±
+
(
h
+
a
v
)
μ
4
±
]
\pmb{g}^\plusmn(\pmb{u}) = \frac{\rho}{2\gamma}\begin{bmatrix} 2(\gamma-1)\mu_1^\plusmn+\mu_3^\plusmn+\mu_4^\plusmn \\ 2(\gamma-1)u\mu_1^\plusmn+u\mu_3^\plusmn+u\mu_4^\plusmn \\ 2(\gamma-1)v\mu_1^\plusmn+(v-a)\mu_3^\plusmn+(v+a)\mu_4^\plusmn \\ (\gamma-1)(u^2+v^2)\mu_1^\plusmn+(h-av)\mu_3^\plusmn+(h+av)\mu_4^\plusmn \end{bmatrix}
g±(u)=2γρ
2(γ−1)μ1±+μ3±+μ4±2(γ−1)uμ1±+uμ3±+uμ4±2(γ−1)vμ1±+(v−a)μ3±+(v+a)μ4±(γ−1)(u2+v2)μ1±+(h−av)μ3±+(h+av)μ4±
其中两个方程中的
λ
i
±
,
μ
i
±
\lambda^\plusmn _i,\mu^\plusmn_i
λi±,μi±
λ
i
±
=
1
2
(
λ
i
±
λ
i
2
+
ϵ
2
)
μ
i
±
=
1
2
(
μ
i
±
μ
i
2
+
ϵ
2
)
\lambda^\plusmn_i = \frac{1}{2}(\lambda_i\plusmn\sqrt{\lambda_i^2+\epsilon^2})\\ \mu^\plusmn_i = \frac{1}{2}(\mu_i\plusmn\sqrt{\mu_i^2+\epsilon^2})
λi±=21(λi±λi2+ϵ2)μi±=21(μi±μi2+ϵ2)