Sod激波管 五阶WENO 流通矢量分裂 三阶RK

一、流通矢量分裂

1.使用密度 d 和压力 p 构造声速 c

2.求出特征值 λ1=u ，λ2=u+c ，λ3=u-c

3.矢量分裂，将特征值 λ1， λ2， λ3 ，epsl=0.000001

以λ+ =（（λ² + epsl²）^0 .5+λ ）/2  ，λ- = （（λ² - epsl²）^0.5 -λ ）/2 将三个特征值进行分裂，

4.使用特征值以及对应矩阵求出 F+=。。。 F-=。。。 （F为ρ*u，ρ*u*u+p，以及u*（E+p））

二、五阶WENO

1.构造模板 a=f j-2 ，b=f j-1 ，c= f j ，d= f j+1 ， e=f j+2

2.求光滑度 ISk ：

IS1 = 13/12*（a-2*b+c)*(a-2*b+c）+ 1/4*(a-4*b+3*c)*(a-4*b+3*c)

IS2 = 13/12*(a-2*b+c)*(a-2*b+c) + 1/4*(a-c)*(a-c)

IS3 = 13/12*(a-2*b+c)*(a-2*b+c) + 1/4*(3*a-4*b+c)*(3*a-4*b+c)

3.求权值ck

求h j+1/2 在（a，b，c，d，e）下的差分格式，联立三个三阶差分模板（a，b，c），（b，c，d），（c，d，e），求解出权值c1，c2，c3

4.求 ak

eps=0.000001

a1=c1/(eps+IS1(i, a, b, c)^2

a2=c2/(eps+IS2(i, b, c, d)^2

a3=c3/(eps+IS3(i, c, d, e)^2

4.求 wk

w1=a1/（a1+a2+a3）

其余同理

5.构造新通量

    q1=1/3*a - 7/6*b + 11/6*c
    q2=-1/6*b + 5/6*c + 1/3*d
    q3=1/3*c + 5/6*d - 1/6*e

h j+1/2 = w1*q1 + w2*q2 + w3*q3

6.求出差分

dF/dx=（h j+1/2  -  h j-1/2）/Δx

三、三阶RK推进

1.求出L（U）

原公式为 dU/dt + dF/dx = 0

变换得 L（U） = - Δt*（dF/dx）=-Δt*（h j+1/2  -  h j-1/2）/Δx

2.在边界处构造虚单元，虚单元参数与边界单元参数随时相同

3.代入第一次推进，得到ρ1，p1，u1，c1，再次构造特征值并分裂，而后再WENO差分，得到新L（U）后二次推进，后续步骤类似。