《5阶weno求解sod激波管》

-29)?将(2-25)代入(2-26),得(2-30)解(2-27),得(2-31)模板3因为(2-32)?设(2-33)?将(2-29)代入(2-30),得(2-34)?解(2-31),得(2-35)理想权重设(2-36)?将(2-20)、(2-24)、(2-28)、(2-32)代入(2-33),得(2-37)?解(2-34),得(2-38)光滑因子的计算定义(2-39)?其中: 使用模板k得到的插值函数,即(2-40)将(2-37)代入(2-36),得(2-41)模板内导数的近似第一个模板内因为(2-42)设(2-43) 将(2-39)代入(2-40),得(2-44) 解(2-41),得(2-45) 设(2-46) 将(2-39)代入(2-43),得(2-47) 解(2-44),得(2-48)第二个模板内因为(2-49) 设(2-50)将(2-46)代入(2-47),得(2-51) 解(2-48),得(2-52) 设(2-53) 将(2-46)代入(2-50),得(2-54) 解(2-51),得(2-55)第三个模板内因为(2-56)设(2-57) 将(2-53)代入(2-54),得(2-58) 解(2-55),得(2-59) 设(2-60) 将(2-53)代入(2-57),得(2-61)?解(2-58),得(2-62)三、作业:EquationSection(Next)EquationSection3(题5.1)针对如下Sod激波管问题(3-65) (3-66)用5阶WENO格式计算其数值解,画出t=0.14时刻密度、速度及压力的分布;并与精确解进行比较(要求数值解与精确解画在同一张图上,便于比较)。要求:空间网格数100,时间推进格式选用3阶Runge-Kutta,时间步长自选。解:数值算法描述空间离散格式空间离散使用5阶WENO离散格式,对于模型方程(3-67)