几个偏微分方程保持守恒性或者耗散性的算法研究

来源 :南京师范大学 | 被引量 : 0次 | 上传用户:suals
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着科学技术的迅猛发展,偏微分方程被用于描述越来越多的物理学、化学、地质、材料等实际问题。由于大多数偏微分方程的精确解很难用解析形式表示,诸多学者纷纷从数值算法的角度对偏微分方程进行研究,设计了各种各样的数值算法。众所周知,很多偏微分方程自身都具有一些原问题的本质特征。如何设计高效、稳定且保持方程本身特性的数值算法,成为数值算法的重要目标。1.能量守恒和质量守恒不仅是力学系统的重要守恒量,也是评价一个数值算法好坏的重要标准。因此,构造保能量守恒和质量守恒的算法引起了诸多科研工作者的兴趣。我们发现,相关的基于有限差分方法和谱配置方法的工作已有较多的积累,本文尝试基于Galerkin谱元算法讨论保结构算法。我们证明了基于线性元半离散耦合非线性薛定谔方程的Hamiltonian性质。随后,我们严格证明了 Galerkin谱元算法保离散能量守恒律和质量守恒律。我们利用插值理论、投影理论等,证明了 Galerkin谱元算法的收敛性。为了提高计算效率,我们采用矩阵对角化方法和快速算法,对Galerkin谱元算法进行显式化处理,减少计算规模。最后,数值实验表明新算法不仅能得到较好的数值解,而且能保质量守恒和能量守恒。通过与文献中已有算法比较,本文算法能更好地保质量守恒和能量守恒。2.辛守恒是哈密尔顿系统的最重要的几何性质。早在1984年,冯康先生提出了辛算法的概念。辛算法在解的稳定性、有效性以及长时间精确模拟解的行为方面有出色的表现。因此,文献中出现了很多哈密尔顿偏微分方程辛算法的研究成果。据我们所知,多数辛算法主要基于有限差分方法、小波方法、Fourier拟谱方法,而没有基于Galerkin有限元的辛算法。文献中关于二维哈密尔顿系统的辛算法讨论也很少。另外,我们发现辛算法主要是隐式的,这会很大程度影响计算效率。因此,本文我们尝试为二维非线性薛定谔方程构造Galerkin分裂辛算法,即在空间方向用Galerkin有限元方法离散得到一个哈密尔顿半离散系统,再利用Strang分裂方法将半离散系统分裂成两个子系统,一个是线性哈密尔顿子系统,一个是非线性子系统,其中利用辛中点在时间方向离散线性子系统,而对非线性子系统精确求解,从而得到一个保辛的全隐的数值算法。考虑计算规模比较大,我们引入FFT技术和矩阵对角化方法,将数值算法进行显式化处理,克服了隐式算法计算效率低的弱点。我们理论证明了 Galerkin分裂辛算法的稳定性及守恒性。最后,数值实验表明新的算法不仅能得到较好的数值解,而且能保系统的质量守恒。通过与文献中已有数值算法的比较,展现了 Galerkin分裂辛算法具有计算效率高等优势。3.目前基于Hamiltonian保结构算法已经很多的工作。我们尝试对带能量耗散的偏微分方程进行保结构算法的研究。带能量耗散的偏微分方程方程包含更多的动力系统,并且应用更为广泛,如:Birkhoffian动力系统,高震荡系统,共形Hamiltonian系统等。接下来,我们主要考虑带能量耗散的Cahn-Hilliard方程。不可压缩多相流问题是多材料混合流体力学系统中的热门问题。Cahn-Hilliard方程是处理不可压缩多相流自由界面问题的最常用的相场模型。能量耗散是相场模型系统的重要性质。目前,文献中已有大量关于Cahn-Hilliard方程的研究成果。我们注意到已有的成果主要基于全局能量耗散的性质。本文,我们发现Cahn-Hilliard方程具有一个局部能量耗散的性质。局部能量耗散性质不依赖于任何边界条件,能最大程度保留了原系统的信息。基于这一发现,我们利用复合构造方法和离散变分导方法,尝试为Cahn-Hilliard方程构造了三个局部保结构算法,严格证明了三个算法保局部能量耗散的性质。当给定合适的边界条件时,三个算法保全局能量耗散的性质也得到了证明。最后,数值实验表明新的数值算法不仅能得到较好的数值解,而且保系统能量耗散的性质和质量守恒。4.对Cahn-Hilliard方程已有的大量数值研究成果研究发现,显式的数值算法易于计算,但往往对步长有一定的限制,而隐式的数值算法一般会影响计算效率,但通常保全局能量耗散性质。因此,保能量耗散的线性隐的数值算法成为较好的选择。因此,学者们尝试为Cahn-Hilliard方程构造线性隐的保全局能量耗散的算法。目前,凸分裂方法、线性稳定化方法和能量二次化方法是构造线性隐式的全局保能量耗散算法最常用的方法,而基于保局部能量耗散性质构造线性隐式的数值算法仍是空白。因此,本文我们尝试为带变系数的Cahn-Hilliard方程构造线性隐的保局部能量耗散的算法。类似的,我们证明了带变系数的Cahn-Hilliard方程具有一个局部能量耗散性质。基于这一发现,本文我们利用能量二次化方法与复合构造方法,尝试为带变系数的Cahn-Hilliard方程构造了四个线性隐的保局部能量耗散的算法,并利用离散莱布尼兹法则,理论上严格证明了新算法保局部能量耗散的性质。当给定适当边界条件时,我们证明了新算法保全局能量耗散和质量守恒的性质。
其他文献
研究了SUMMA罐采样-全二维气相色谱-飞行时间质谱法测定环境空气中有机硫化物。采用SUMMA罐采集样品,自动进样器进样,三级冷阱预浓缩样品,全二维气相色谱-飞行时间质谱法(GC
选用杂交粳稻甬优8号,以国家粮食丰产工程兴化、姜堰实施基地1.0hm2连片与6.67hm2连片超高产攻关方为依托,研究了13.5thm-2超高产群体特征,并探讨了群体形成机制。结果表明,
本文运用文献计量学的方法。对2003年《郑州大学学报》(哲社版)的载文、作者及引文情况进行了统 计和分析。从而反映出该刊的一些特点和学术影响,为进一步办好刊物提供依据。
我国古代书院"博学"、"审问"、"慎思"、"明辨"和"笃行"的治学理念,对推进基础教育改革,培养创新型人才具有重要的现实意义。我国当前的基础教育改革可以借鉴古代书院广博深厚
采用气相色谱-质谱联用(GC-MS)技术对海红豆油脂肪酸组成进行分析和鉴定,并对海红豆油的理化性质进行分析研究。结果表明:从海红豆油中共分离出24个峰,鉴识出其中16种脂肪酸,
本文重点叙述了20世纪以来,天津旧书业的发展概况,尤其是雷梦辰、许杏林、王辛笛、张志公、孙犁、黄裳、姜德明、许政扬、来新夏等人笔下的旧书业风情,对于"藻玉堂"主人王雨
花丝活力长短在一定程度上影响着玉米制种产量和杂交种产量,对玉米生产具有重要意义。为阐明玉米花丝活力的遗传机理,本研究利用来自优良玉米杂交种农大108的一套203个重组近
目的研究神经内科护理中实施优质护理服务的临床效果。方法通过分组在常规护理的基础上实施优质护理措施,比较两组患者的住院时间以及对护理工作的满意度。结果 通过对两组
采用高效液相色谱-串联四极杆飞行时间质谱法研究大豆油甘油三酯组成特性。对大豆油中16种甘油三酯进行了定性和相对含量测定,分别为LLnLn(0.97%~4.72%)、LLLn(9.06%~14.22%)、
<正>江泽民总书记在1999年6月到陕西视察时指出,“加快开发西北地区应该是全面的,要把水资源的开发和有效利用放在突出位置”.“再造一个山川秀美的西北地区”.为西北地区的