流通矢量分裂

注:本文为贴出各种流通矢量分裂方法的方法/格式,随时记录,以便需要用时随时查看

流通矢量分裂

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=ua,λ4=u+aμ1=v,μ2=v,μ3=va,μ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±+(ua)λ3±+(u+a)λ4±2(γ1)vλ1±+vλ3±+vλ4±(γ1)(u2+v2)λ1±+(hau)λ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±+(va)μ3±+(v+a)μ4±(γ1)(u2+v2)μ1±+(hav)μ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 )

  • 3
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值