收藏本站
《中国科学技术大学》 2012年
收藏 | 手机打开
二维码
手机客户端打开本文

WENO方法的研究及其在对流扩散方程中的应用

刘元元  
【摘要】:本文主要研究加权本质无振荡数值方法(以下简称WENO方法)及其在对流扩散方程中的应用。 WENO方法是近年来广泛流行的一种高分辨率数值方法,用于解决以对流为主的对流扩散方程,特别是双曲守恒律问题。WENO方法具有在光滑区可以得到任意高阶精度,在间断区保持稳定的、本质无振荡性质和陡峭的间断过渡等优点。 WENO方法的主要思想是使用所有候选模板上的低阶逼近的凸组合来得到高阶逼近。WENO是一种插值或重构过程,可以方便地应用到有限差分格式或有限体积格式的框架中。在实际应用中,WENO过程有不同的表现形式:WENO插值、WENO重构和WENO积分。WENO构造过程中,格式的稳定性和无振荡性质主要依赖于线性权的正定性。我们讨论线性权在不同形式的WENO过程中的正定性,并给出线性权的显式表达和正定区间,以便为进一步设计各种形式的WENO格式包括WENO插值,WENO重构及WENO对高阶导数的逼近,和WENO积分打好基础。 对流扩散方程是粘性流体力学及其他实际问题中常见的一类基本运动方程。对于这类方程,除了极少数简单情形,大部分问题目前还无法求得精确解,只能利用数值方法进行数值模拟。因此,对流扩散方程数值方法的研究具有十分重要的理论和实际应用意义。对求解对流扩散方程,特别是以对流为主的对流扩散方程或者退化抛物方程,数值方法的设计过程中遇到的最大困难是间断的出现或者陡峭的过渡层。标准的WENO格式是用来求解双曲守恒律方程和以对流为主的对流扩散方程中的对流项,通常在求解对流扩散方程的WENO格式中,只对方程中的一阶对流项用WENO过程逼近,而对二阶导数扩散项用标准的中心差分逼近。这个方法在大多数的数值模拟中可以工作得很好,但是在求解有间断解的退化抛物方程中,会产生振荡的数值解甚至不稳定的格式。我们对二阶导数项构造守恒的有限差分WENO逼近,以多孔介质方程为例,考虑求解有间断解的非线性退化抛物方程,并应用推广到一维和二维问题。非线性退化抛物方程的解具有双曲守恒律方程解的类似性质,比如陡峭界面和紧支集界面以有限的速度传播。常用的方法是通过变量代换,将二阶方程改写成一阶方程组,然后再应用WENO重构过程求解该方程组,但该方法的相应的模板和计算量都很大。所以我们采用一个守恒的流通量差分,对二阶导数项直接进行逼近。我们首先讨论不同模板选取时的线性权,详细介绍六阶和八阶有限差分WENO格式的构造,理论分析格式的精度,并给出格式的CFL稳定性条件,最后数值试验验证格式的高分辨率和本质无振荡性质。 对流扩散方程解的一个重要性质是满足严格最大值原理,我们同样希望数值解也满足最大值原理,否则就失去了解的实际物理意义,比如放射性核素迁移的计算。多孔介质方程的有限差分WENO格式的数值解虽然是本质无振荡的,但是在陡峭界面处不能严格保证解是非负的。因此,我们将求解标量守恒律方程的满足最大值原理的格式推广,对非线性对流扩散方程,构造一个非传统的满足最大值原理的五阶有限体积WENO格式,并将格式应用推广到二维问题。主要思想是将单元积分平均算子应用两次,介绍双重积分平均值的概念,在此基础上构造五阶有限体积WENO逼近。我们理论分析证明格式的高精度,并在一定CFL稳定性条件下满足最大值原理,最后数值试验验证格式的高分辨率、本质无振荡和满足最大值原理的性质。
【学位授予单位】:中国科学技术大学
【学位级别】:博士
【学位授予年份】:2012
【分类号】:O241.8

手机知网App
【参考文献】
中国期刊全文数据库 前2条
1 孙澈 ,沈慧;THE FINITE DIFFERENCE STREAMLINE DIFFUSION METHODS FOR TIME-DEPENDENT CONVECTION-DIFFUSION EQUATIONS[J];Numerical Mathematics A Journal of Chinese Universities(English Series);1998年01期
2 张强,孙澈;非线性对流扩散问题的差分-流线扩散法[J];计算数学;1998年02期
【共引文献】
中国期刊全文数据库 前10条
1 张飞,汪继文;一维溃坝问题的间断Galerkin方法[J];安徽大学学报(自然科学版);2005年02期
2 姚兴军;张凤阳;王正东;章文俊;;基于VOF和CSF方法的变截面毛细流动数值建模[J];半导体技术;2011年02期
3 李物兰,张大凯;求解对流扩散方程的组合差商算法[J];北方工业大学学报;2005年01期
4 吴开腾;张莉;;一类三维基于通量分裂的高精度MmB差分格式[J];兵工学报;2010年S1期
5 赵建丛,康彤;非定常线性双曲方程的差分──间断Galerkin方法及其稳定性[J];北京广播学院学报(自然科学版);1999年01期
6 张斌;刘宇;王革;;激波与物质界面相互作用的数值模拟[J];北京航空航天大学学报;2009年10期
7 周应国;申长雨;陈静波;;半结晶性聚合物熔体冷却过程双尺度模拟[J];北京科技大学学报;2007年02期
8 方智勇;端木玉;朱仁庆;;基于Level-set法的液舱液体晃荡数值模拟[J];船舶力学;2007年01期
9 张舒娟;侯华;杨晶;毛红奎;;充型过程的数值模拟技术[J];材料导报;2007年03期
10 李晶;高桂景;王平义;;长江神仙桥码头工程对航道的影响研究[J];重庆交通大学学报(自然科学版);2007年04期
中国重要会议论文全文数据库 前9条
1 王成;郝莉;宁建国;;多物质界面追踪方法研究[A];第八届全国爆炸力学学术会议论文集[C];2007年
2 吴开腾;张莉;;一类三维基于通量分裂的高精度MmB差分格式[A];第九届全国冲击动力学学术会议论文集(下册)[C];2009年
3 谢昌坦;王革;;强激波与气-水界面相互作用的数值模拟[A];中国计算力学大会'2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集[C];2010年
4 强洪夫;陈福振;高巍然;;基于SPH方法的表面张力修正算法及其应用[A];中国计算力学大会'2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集[C];2010年
5 王博;李文成;邓子辰;;基于间断有限元方法的高精度数值计算[A];中国计算力学大会'2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集[C];2010年
6 缪国平;朱仁传;程建生;王景全;;海上风电场建设与海洋工程装备研发中若干水动力学关键技术问题[A];第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议暨两岸船舶与海洋工程水动力学研讨会文集[C];2008年
7 汪迎春;汪德爟;;立面溃坝水气两相流计算[A];第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议暨两岸船舶与海洋工程水动力学研讨会文集[C];2008年
8 王革;关奔;;激波作用下R_(22)气泡射流现象研究[A];第十五届全国激波与激波管学术会议论文集(下册)[C];2012年
9 汪文帅;张怀;李小凡;;间断的Galerkin方法在地震波场数值模拟中的应用概述[A];中国科学院地质与地球物理研究所2013年度(第13届)学术论文汇编——地球深部结构与过程研究室[C];2014年
中国博士学位论文全文数据库 前10条
1 付天光;逐孔起爆技术应用基础研究[D];辽宁工程技术大学;2010年
2 闫铂;反应扩散方程的格子Boltzmann方法及数值模拟[D];吉林大学;2011年
3 张建影;用于复Ginzburg-Landau方程的格子Boltzmann方法[D];吉林大学;2011年
4 李刚;基于多种指示子的杂交WENO格式[D];南京大学;2011年
5 曹文炅;铸造充型过程SPH方法建模及数值模拟[D];华南理工大学;2011年
6 李翠娟;超大跨径CFRP主缆悬索桥合理结构体系研究[D];西南交通大学;2011年
7 林毅;自由表面流动问题数值方法的理论研究及应用[D];天津大学;2010年
8 马千里;CFD非结构化网格流场体可视化方法研究[D];国防科学技术大学;2011年
9 王晓峰;不可压缩黏性流动与强化传热问题的高精度紧致差分方法[D];武汉大学;2011年
10 刘江川;泥石流数学模型构建及危险性评价研究[D];天津大学;2011年
中国硕士学位论文全文数据库 前10条
1 侯建军;高温矿井热环境数值模拟及热害控制技术研究[D];河南理工大学;2010年
2 侯鹏;基于双曲型方程RKDG法在流场及声场中的数值模拟研究[D];哈尔滨工程大学;2010年
3 高慧颖;基于CFD的大型船舶球艏线型设计比较研究[D];哈尔滨工程大学;2010年
4 郭涛;偏流板性能的分析与计算[D];哈尔滨工程大学;2010年
5 刘丽莉;电磁场驱动下的液体流动实验与模拟研究[D];辽宁工程技术大学;2009年
6 张凤阳;倒装芯片封装中下填充流动行为的数值模拟[D];华东理工大学;2011年
7 刘伟;基于混合网格的高阶间断Galerkin/有限体积混合算法的研究[D];中国空气动力研究与发展中心;2010年
8 朱培模;高效节能搅拌器的数值模拟及性能预测[D];昆明理工大学;2008年
9 林雄伟;恒流量控制阀的设计及系列化研究[D];浙江大学;2010年
10 唐丽娜;抛物型发展方程的间断体积元方法[D];山东师范大学;2011年
【相似文献】
中国期刊全文数据库 前10条
1 张池平,崔明根;再生核空间中求解定态对流扩散方程[J];应用数学和力学;1994年10期
2 陆金甫;对流扩散方程的一些单调性差分格式[J];计算物理;1991年02期
3 梁栋;对流扩散方程的一类迎风格式[J];计算数学;1991年02期
4 张志跃;变系数对流扩散方程的交替分段显-隐式方法[J];山东大学学报(理学版);2002年02期
5 张大凯;李洪林;;求解对流扩散方程的两层半显式差分格式[J];贵州科学;1991年01期
6 窦红;对流扩散方程的一种显式有限体积-有限元方法[J];应用数学与计算数学学报;2001年02期
7 由同顺;求解对流扩散方程的ICT-MMOCAA差分解法[J];高等学校计算数学学报;2002年02期
8 邹秀芬;对流扩散方程的格点模型[J];计算物理;1996年03期
9 崔明荣;动边界对流扩散问题的特征有限差分法[J];山东大学学报(自然科学版);1995年03期
10 武周虎;对流扩散方程FTCS差分格式稳定性分析研究[J];四川水力发电;1988年02期
中国重要会议论文全文数据库 前10条
1 崔霞;;抛物型对流扩散方程第三类边值问题的有限体积数值分析[A];中国工程物理研究院科技年报(2003)[C];2003年
2 葛永斌;田振夫;吴文权;;三维非定常对流扩散方程的高精度全隐式多重网格方法[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
3 王革;谢昌坦;;应用五阶WENO格式和MGFM计算多介质可压缩流[A];第十四届全国激波与激波管学术会议论文集(下册)[C];2010年
4 高智;;对流扩散方程的摄动有限体积(PFV)方法及讨论[A];计算流体力学研究进展——第十一届全国计算流体力学会议论文集[C];2002年
5 林建国;刘颖;;求解对流扩散方程的局部解析算法[A];第十六届全国水动力学研讨会文集[C];2002年
6 徐立;;可压缩流体二维R-T不稳定性的直接数值模拟[A];中国空气动力学学会物理气体动力学专业委员会第十一届学术交流会会议论文集[C];2003年
7 丁家平;;地下水污染问题的对流扩散方程在运动质点坐标系上的数值解法[A];水工渗流研究与应用进展——第五届全国水利工程渗流学术研讨会论文集[C];2006年
8 林建国;殷佩海;;污染浓度数值解的任意阶精度格式[A];中国航海学会船舶机电专业委员会2000年度学术报告会论文集[C];2000年
9 李占松;朱士江;;黄河合理整治宽度研究的理论思考[A];2007重大水利水电科技前沿院士论坛暨首届中国水利博士论坛论文集[C];2007年
10 林建国;谢志华;由晓丹;;污染物对流扩散方程的Hopscotch方法及几种改进格式[A];第七届全国水动力学学术会议暨第十九届全国水动力学研讨会文集(上册)[C];2005年
中国博士学位论文全文数据库 前10条
1 刘元元;WENO方法的研究及其在对流扩散方程中的应用[D];中国科学技术大学;2012年
2 李刚;基于多种指示子的杂交WENO格式[D];南京大学;2011年
3 王晓玲;对流扩散方程的特征有限元方法以及气泡行为模拟研究[D];郑州大学;2010年
4 王红梅;对流扩散方程的特征有限元方法[D];山东大学;2012年
5 祝鹏;Hamilton-Jacobi方程与对流扩散方程的新算法[D];湖南大学;2010年
6 王凯欣;多孔介质流体的数值方法及其分析与计算[D];山东大学;2010年
7 张瑞;有限体积WENO格式及其应用[D];中国科学技术大学;2010年
8 武从海;流体力学高精度高分辨率差分格式的研究[D];南京航空航天大学;2012年
9 朱嵩;基于贝叶斯推理的环境水力学反问题研究[D];浙江大学;2008年
10 廉西猛;海水入侵问题的间断有限元方法[D];山东大学;2010年
中国硕士学位论文全文数据库 前10条
1 孙玉晓;对流扩散方程的有限差分法[D];西南石油大学;2011年
2 冯安住;间断系数对流扩散方程初值及其区域分解算法[D];长沙理工大学;2012年
3 王宝艳;定常对流扩散方程的一种新型差分格式[D];东北大学;2008年
4 丁晓燕;求解非定常对流扩散方程的高精度指数型差分方法[D];宁夏大学;2013年
5 李芹芹;时—空分数阶对流扩散方程对流项系数识别[D];山东大学;2012年
6 曹广满;对流扩散方程的一种新型紧致差分方法[D];天津师范大学;2010年
7 魏海燕;对流扩散方程的高精度有限差分方法[D];辽宁工程技术大学;2011年
8 田文义;空间分数阶对流扩散方程的多项式谱配置方法[D];兰州大学;2012年
9 姜伟;一维对流扩散方程的同伦分析解[D];华北电力大学;2012年
10 盛志明;对流扩散方程的并行差分格式[D];中国工程物理研究院;2010年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026