解双曲方程的wscip方法.pptVIP

数值实验-一维Euler方程的Riemann问题 求解过程：i) 非对流阶段使用二阶Runge-Kutta方法 求解： ii) 对流阶段使用WSCIP方法： 注：在计算过程中我们根据速度的变差符号对压强进行了如下修正： 数值实验-一维Euler方程的Riemann问题 1）Test 1 The solution of WSCIP and the exact solution for (left),u (right) at t=0.2513,△x=0.0033,CFL=0.3 The solution of WSCIP and the exact solution for energy (left),p (right) at t=0.2513,△x=0.0033,CFL=0.3 2）Test 2 The solution of WSCIP and the exact solution for (left),u (right) at t=0.1667,△x=0.0033,CFL=0.3 The solution of WSCIP and the exact solution for energy (left),p (right) at t=0.1667,△x=0.0033,CFL=0.3 WSCIP方法的总结 优势： 与传统CIP方法相比，略去了计算导数方程过程，从而降低了计算复杂度，对解函数的光滑性没有任何限制； 稳定性比传统CIP方法要好，在大梯度或间断处一般不会出现振荡; 在没设置限制器的情况下，对接触间断、激波及稀疏波均有较高的分辨能力。 不足之处： 逼近阶略低于传统CIP； 捕捉激波时有相位滞后现象。 谢 谢！  解双曲方程的C1三次加权CIP方法（ The Solution of Hyperbolic Equations by C1 Weighted Spline Interpolation of Three Degree） 冯仁忠 北京航空航天大学数学与系统科学学院 2012/12/08 湘潭大学 报告内容 CIP方法的发展历史与现状 CIP方法原理及弱点 加权样条 加权样条CIP方法 数值实验 总结 CIP方法的发展历史与现状 CIP（Cubic Interpolated Propagation)方法是日本学者T.Yabe等人于上世纪八十年代构造的一种求解双曲方程的显式紧致格式。该方法的主要思想是根据解当前值及导数值在每个单元上构造三次插值多项式，借助插值多项式和对流方程的Lagrange不变量，估计下一时刻的解值及导数值。 优点：计算简单、低耗散、稳定、高CFL条件 发展历程简介： 一维→高维（1985,1987，1991） 守恒形式（CIP-CSL2，CIP-CSL3,CIP-CSL4，2000，2001） 三角网格上的高阶守恒形式（CIP/MM-FVM，2005,2010） 应用：可压及不可压流、多相流、大气模型 CIP方法原理 线性对流方程 的特征 ，其解为 如何估计 ？ CIP方法原理 三次插值重构： 其中 CIP方法对 的估计： × CIP方法原理 CIP格式中导数的选取: 假设方程解的导数由下列主控制方程确定 CIP方法是一个具有空间精度三阶的非守恒紧致显格式，使用中不需要象DG方法那样进行相邻单元间信息的交换。 CIP方法的弱点 CIP方法的弱点： 在大梯度附近产生虚假振荡 要求解的导数连续，对于存在间断的解是不可能有此性质 需要解导数方程来获得导数值，间接地增加了计算量 在模拟间断解时会出现数值解相位滞后现象 在空间剖分网格 上构造如下C1三次插值样条 加权样条(weighted splines) 其中 加权样条(weighted splines) s(x)满足下列条件： s(x)中的参数 通过下列方式确定 此时s(x)中的参数 满足 加权样条(weighted splines) 若权函数w(x)在 上是分段常数，则 其中 Theorem Let Then the interpolating spline s with a weight function which is piecewise constant on the partition