求解双曲守恒律的紧-WENO杂交格式及RBF-FD间断检测方法
发布时间:2021-12-23 06:48
设计准确、鲁棒和高效的间断检测方法,在任意给定的网格单元和任意时间上衡量解的光滑度,是杂交格式的关键问题.本文从有限差分形式的径向基函数(RBF-FD)方法出发设计间断检测方法,为解决其阈值参数的选取问题,采用数据分析中Tukey’s箱线图方法增强算法的鲁棒性,最终提出一种准确、鲁棒和高效的RBF-FD间断检测方法.数值结果分析表明,RBF-FD间断检测方法在间断捕捉的准确性、鲁棒性和效率等方面优于已有的RBF方法和单调多项式插值方法,并且与五阶WENO-Z格式相比,应用RBF-FD方法的杂交格式可以加速2倍左右.
【文章来源】:数值计算与计算机应用. 2020,41(03)
【文章页数】:14 页
【部分图文】:
图1五阶WBIfQ模板,萁中,U:?=?〇,?为包含H个点的子模板,#?=:?P〇U??i?U袅为整体模板??
240??数值计算与计算机应用??‘2020?年??其中,e?=?0.2,抑=-弋=?5.终止时间为t?=?5.两侧边界为固定边界,其值与初始边界相??同.取网格节点数况=8〇〇.??密度??RJPWD方法的方lag随时间的演化??°-5??图2求解Sltod^Etemity问题的杂交格式中间断检测算法捕捉的间断??三种间断检测方法在f?=?5时捕捉的间断如图3所被,其右侧子图展示了?HBP-FD方法的??间断标识F1够随时间的演化过程,其中黑色部分为Fkg?=?1的位置.BBF-FD方法和RBF方??法均准确地检测出间断的位置,包括主激波和随时间演化的小激波的位置.另一方面,两种方??法并没有像MI方法那样将连续、光滑的髙频波区域误判为间断.??5.4.?Blast-Waves?问题??考虑含有强间断结构的Blast-Waves相互作用问题I?1,初始条件为??f?(?1??1000??0.0?<x?<?0.1,??:(簡?p)=彳(1??,〇,??0.01??),??0.1?<?x?<?0.9,??1?(?1??,〇,??100??),??0.9?<?x?<?1.0,??其中,1?=?1.4,M?=?800,虽?终时刻为f?=?0.038.在立=0和:c?=?1处使用反射边界条件.??注意初始条件中P(%〇)?=?1是一个常数,W此,该问题中需要对密度P和压力P应用间断??检测算法,最终的间断位置为两者的并集.在最后时刻,由紧-WENO格式计算的密度和WENO??标识的时间演化如图3所忒与MI方法相比,采用RBF_断检测算法的紧-WENO杂交格式??能够更准确地捕捉压力和密度的间断位置,并且两种RBF方法检测出的间断
本文编号:3548002
【文章来源】:数值计算与计算机应用. 2020,41(03)
【文章页数】:14 页
【部分图文】:
图1五阶WBIfQ模板,萁中,U:?=?〇,?为包含H个点的子模板,#?=:?P〇U??i?U袅为整体模板??
240??数值计算与计算机应用??‘2020?年??其中,e?=?0.2,抑=-弋=?5.终止时间为t?=?5.两侧边界为固定边界,其值与初始边界相??同.取网格节点数况=8〇〇.??密度??RJPWD方法的方lag随时间的演化??°-5??图2求解Sltod^Etemity问题的杂交格式中间断检测算法捕捉的间断??三种间断检测方法在f?=?5时捕捉的间断如图3所被,其右侧子图展示了?HBP-FD方法的??间断标识F1够随时间的演化过程,其中黑色部分为Fkg?=?1的位置.BBF-FD方法和RBF方??法均准确地检测出间断的位置,包括主激波和随时间演化的小激波的位置.另一方面,两种方??法并没有像MI方法那样将连续、光滑的髙频波区域误判为间断.??5.4.?Blast-Waves?问题??考虑含有强间断结构的Blast-Waves相互作用问题I?1,初始条件为??f?(?1??1000??0.0?<x?<?0.1,??:(簡?p)=彳(1??,〇,??0.01??),??0.1?<?x?<?0.9,??1?(?1??,〇,??100??),??0.9?<?x?<?1.0,??其中,1?=?1.4,M?=?800,虽?终时刻为f?=?0.038.在立=0和:c?=?1处使用反射边界条件.??注意初始条件中P(%〇)?=?1是一个常数,W此,该问题中需要对密度P和压力P应用间断??检测算法,最终的间断位置为两者的并集.在最后时刻,由紧-WENO格式计算的密度和WENO??标识的时间演化如图3所忒与MI方法相比,采用RBF_断检测算法的紧-WENO杂交格式??能够更准确地捕捉压力和密度的间断位置,并且两种RBF方法检测出的间断
本文编号:3548002
本文链接:https://www.wllwen.com/kejilunwen/yysx/3548002.html