收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

基于特征思想的高分辨率格式的研究和应用

郭彦  
【摘要】: 本文为数值求解依赖时间的偏微分方程提出两类基于特征思想的高分辨率格式。从而我们主要考虑两部分内容。首先基于CIP方法和高阶紧致方法,我们提出一种新型特征插值高阶紧致差分方法,并将其运用于数值模拟非线性波动方程。然后我们对守恒型双曲方程提出一种基于特征思想的有限体积方法,并应用于数值求解Eulelr方程组和浅水方程组。 论文的第一大部分,我们对非线性波动方程提出一种非守恒Semi-Lagrangian方法。这种格式的主要构造思想如下:根据两相邻节点的节点值和导数值构造三次多项式,然后运用Semi-Lagrangian的方法求解节点沿特征线移动的位置,由此三次多项式和节点位置得到时间上的推进。不同于传统的CIP格式,我们利用最初由Lele提出的高阶紧致格式直接利用节点值求解一阶导数值。最后我们将这种格式运用于数值模拟Burgers方程和KdV方程验证格式的高分辨率的性质。通过一系列的数值试验验证格式的有效性。最后将这种格式推广到一维的Euler方程组。 论文的第二大部分,我们对双曲守恒方程组构造有限体积格式,由于Euler方程组在空气动力学中的的重要性,我们主要基于Euler方程组构造格式。这种格式的空间导数通过有限体积离散,时间上采用Simpson公式离散。其中有限体积格式中的通量值通过单元边界点的值得到,而该单元边界点的值是通过中心加权重构得到。求解此节点值的主要思想如下:首先运用三阶或四阶Runge-Kutta方法求解特征方程,从而得到节点沿特征线的位置。然后采用CWENO重构得到多项式,我们分别构造三阶和五阶CWENO格式构造多项式,最后利用节点位置和此多项式得到节点值。这种高分辨率的有限体积格式结合特征思想和中心加权格式的优点。然后我们将这种基于特征思想的有限体积格式应用于一维的标量和Euler守恒方程,通过经典的算例验证格式的高分率和收敛的性质。 最后将这种基于特征思想的高阶有限体积格式应用于一些浅水问题,从而验证格式,如一维溃坝问题,临界左稀疏波问题,两稀疏波中间几乎为干底问题,干溃坝问题,生成干溃坝问题等。另外,我们同时将这种格式按照分裂方法推广到二维守恒律情形。并通过对二维的溃坝问题的应用验证格式的有效性,通过与WAF格式以及其它高阶数值格式相比较可以看到该格式得到非常好的数值结果。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 傅凱新;双調和方程数值解中“松弛——消元法”的应用[J];中山大学学报(自然科学版);1960年03期
2 姜长生;控制系统中的矩阵方程的数值解[J];电光与控制;1985年02期
3 吕炯兴;关于矩阵方程AX-XB=C的解[J];南京航空航天大学学报;1986年04期
4 李新芳;一个静电场问题的数值解[J];内蒙古师范大学学报(自然科学汉文版);1986年02期
5 李亮斌;杨翔翔;;在辐射和对流条件下肋片传热的研究[J];华侨大学学报(自然科学版);1989年04期
6 孙秀贞;;两平行圆形电流间的磁场分布的数值解及其匀场区[J];福建师范大学学报(自然科学版);1989年02期
7 姜宗林;;绕流问题的区域收缩法及其数值解[J];空气动力学学报;1989年01期
8 胡纪阳,胡瑞安;Mandelbrot数集和非线性过程的可视性图象[J];计算机辅助设计与图形学学报;1991年04期
9 章自振;;用正交多项式求线性微分方程的数值解[J];河南科技大学学报(自然科学版);1991年03期
10 郗大光,王荣生,陈家骅;一种求解气动双曲方程的高精度差分格式DEM[J];空气动力学学报;1992年04期
11 金保侠;一类TVD格式的熵强迫函数及熵条件[J];计算数学;1993年04期
12 何宏舟;杨翔翔;;梯形肋片端部第三类边界条件的换热研究[J];华侨大学学报(自然科学版);1993年03期
13 周军;矩阵Lyapunov方程的直接数值解[J];兰州大学学报(自然科学版);1994年03期
14 曹建胜,傅少川;约束矩阵方程组的解[J];山东工业大学学报;1998年03期
15 钟坦谊,邓中兴;第一类Fredholm积分方程的Hermite数值解[J];工程数学学报;1999年02期
16 蒋虎,黄王成,胡小工;二体问题状态转移矩阵的分析解和数值解的比较[J];云南天文台台刊;1999年04期
17 罗兴钧;Fourier型问题的数值解[J];赣南师范学院学报;2001年06期
18 孙光辉;利用计算机模拟RC电路实验[J];青岛大学学报(自然科学版);2001年02期
19 柯璇,张祖全,张绪宽;精确聚焦非球面透镜数值解及其计算机模拟检验[J];应用光学;2001年06期
20 杨明辉;矩阵方程AX+XB=C的一种数值算法[J];南京林业大学学报(自然科学版);2002年06期
中国重要会议论文全文数据库 前10条
1 任永亮;李鸿晶;王通;;结构固有振动分析的广义微分求积法[A];第一届全国工程安全与防护学术会议论文集[C];2008年
2 齐朝晖;陈礼;;多体系统动力学数值分析的力学解释初探[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
3 张天元;莫俊永;;核三体Faddeev方程数值解[A];第八届全国核物理会议文摘集(上册)[C];1991年
4 常静;李永安;刘学亭;;采暖散热器房间室内温度场分布研究[A];2006年山东省制冷空调学术年会“格力杯”优秀论文集[C];2006年
5 傅向荣;赵东;蒋东英;王宏志;赵阳;田歌;张浩;;楔形板的解析解及其数值解答探讨[A];北京力学会第十六届学术年会论文集[C];2010年
6 罗正钱;黄朝红;蔡志平;许惠英;刘孙丽;;掺磷拉曼光纤激光器的解析解及数值模拟[A];全国第十二次光纤通信暨第十三届集成光学学术会议论文集[C];2005年
7 石文静;孟庆春;;含穿透分层双金属梁数值解与片条合成能量解法的比较[A];中国力学学会学术大会'2009论文摘要集[C];2009年
8 刘乐民;尚新春;;磁致伸缩作动器振动问题的数值分析[A];北京力学会第13届学术年会论文集[C];2007年
9 马家驹;;线接触副凸度设计的研究 (Ⅰ)计算方法[A];第五届全国摩擦学学术会议论文集(上册)[C];1992年
10 刘长学;;200000阶稀疏矩阵方程的求解问题[A];土木工程中计算机应用文集——中国土木工程学会计算机应用学会成立大会暨第一次学术交流会论文集[C];1981年
中国博士学位论文全文数据库 前10条
1 高峰;代数函数逼近及其在微分方程数值解中的应用[D];大连理工大学;2004年
2 谭建国;泊松白噪声激励下的随机时滞系统的数值解及应用研究[D];天津大学;2010年
3 纪晓东;半空间中圆形洞室对弹性波的散射[D];天津大学;2005年
4 陈娓;关于自然电位测井的若干问题[D];复旦大学;2008年
5 沈淑君;分数阶对流—扩散方程的基本解和数值方法[D];厦门大学;2008年
6 王志勇;随机泛函微分方程的稳态数值解研究[D];华中科技大学;2008年
7 冯晓莉;椭圆和抛物方程的几类不适定问题[D];兰州大学;2010年
8 赵桂华;几类带跳随机微分方程数值解的收敛性和稳定性[D];哈尔滨工业大学;2009年
9 梅欢;谱元方法求解不可压缩流体流动及流动线性稳定性分析[D];重庆大学;2012年
10 康洪朝;高振荡问题的高效数值算法研究及实现[D];中南大学;2012年
中国硕士学位论文全文数据库 前10条
1 石兰芳;奇摄动问题的数值解[D];安徽师范大学;2002年
2 齐红妍;带跳BSDE的数值解法及其在金融中的应用[D];北方工业大学;2014年
3 张希;有损均匀传输线数值解的研究[D];重庆大学;2003年
4 谢支钢;边坡稳定性分析的数值解研究[D];武汉理工大学;2003年
5 钟澎洪;一类Landau-Lifshitz和Ginzburg-Landau方程的精确解与数值解[D];云南民族大学;2009年
6 刘云举;两类流体方程的均质化和数值解[D];华中科技大学;2012年
7 傅敏红;一种用于土—结构相互作用体系非线性地震动力响应分析的基于数值解与解析解耦合的新方法[D];厦门大学;2014年
8 廖肇源;样条函数在微分方程数值解中的若干应用[D];大连理工大学;2006年
9 NKOSINATHI DLAKAVU;抽水蓄能电站叉管损失影响的研究[D];河海大学;2006年
10 岳华;一类三阶常微分方程排放和吸收流边值问题的数值近似[D];吉林大学;2007年
中国重要报纸全文数据库 前3条
1 陆彦文;举世瞩目的“有生于无”难题[N];科技日报;2002年
2 杨振宁;他永远脚踏实地[N];光明日报;2001年
3 郝敏;数学有什么用?[N];中华读书报;2003年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978