若干四阶非线性偏微分方程的数值解法

发布时间：2018-01-26 05:22

  本文关键词： 四阶非线性抛物方程 差分法 拟谱方法 Fourier谱方法 长时间行为　出处：《吉林大学》2017年博士论文　论文类型：学位论文

【摘要】：高阶非线性抛物方程作为数学模型可描述物理、化学、地理、环境科学等很多领域出现的现象,应用广泛,是非线性学科的一个重要组成部分.本文研究了三类具有应用背景的四阶非线性偏微分方程的数值解法.理论上主要研究数值方法的先验估计和收敛性,特别是数值解的H~2模先验估计,这是四阶方程收敛性分析的基础.此外,文中也给出了数值实验以验证理论上的结果.首先,本文研究了描述晶体表面增长的方程的差分法.关于该方程数值解法的研究不多,笔者曾经和其他作者合作研究过该问题的有限元法,但由于三次有限元形成刚度矩阵的计算量过大,因而又进一步来研究方程的差分法.本文针对一维模型和二维模型各给出了一种差分格式,分别证明了差分解的存在唯一性、先验估计和收敛性,并给出了数值算例.该方程的非线性项是有理式形式,直接处理比较复杂,但通过一次分部积分可以大幅降低非线性项的复杂度.在有限元法的研究中是用连续内积进行分析,从而可以直接用分部积分处理非线性项.但在构造差分格式时,如果只是简单移植,那么非线性项在离散内积下不满足分部积分.在不影响计算复杂度的前提下,为了能采取与有限元法相似的处理方法,本文针对该方程的一维模型构造了一种半隐半显型线性化差分格式,使得非线性项在相应的离散内积下可以进行一次分部积分,从而降低非线性项的处理难度.在理论上证明了该数值格式在时间方向上为一阶收敛,在空间方向上为二阶收敛,数值算例验证了这一结果.针对二维模型本文构造了 Crank-Nicolson型的差分格式,并定义了与之匹配的离散内积,使得非线性项也可以进行类似的处理.所构造的Crank-Nicolson型的差分格式虽然比线性化差分格式的计算量大,但在时间上达到二阶收敛,综合效果较好.数值结果也与该理论结果相吻合.其次,本文研究了二维扩展的Fisher-Kolmogorov(EFK)方程u_t + γ△~2u-△u + f(u)= 0的拟谱方法.关于该方程的数值解法有较多的研究成果,方法涉及各类有限元法和差分法.本文先对空间进行离散,构造了恰当的求解空间和离散内积,使得求解空间SN以及离散内积的定义不仅与方程的初边值条件相匹配,同时也满足SN空间中基函数的正交性.在此基础上,构造了半离散拟谱近似格式,并证明了半离散解的存在性、先验估计及收敛性.进而又建立了Crank-Nicolson型的全离散拟谱格式,证明了全离散解的先验估计,并给出了按L2范数的最优阶误差估计.针对EFK方程的非线性项f(u)= u~3-u,在H~2半模先验估计的分析中引入能量泛函E_h(t)=γ/2||△u_h||~2+1/2||%絬_h||~2+(H(u_h),1)_h,u_h∈sN,其中H(u_h)=1/4(1-u_h~2)~2,并且H'(u_h)=f(u_h).利用能量泛函随时间递减的性质证得H~2半模的先验估计.全离散解的先验估计证明中也构造和借助了相应的离散能量泛函.给出了 Fourier拟谱方法的收敛性结果,并用数值算例验证了相应的理论结果.最后,本文研究了一类对流Cahn-Hilliard方程u_t+γ△~2u-△f(u)%絞(u)=0的Fourier谱近似解的长时间行为,其f(u)为关于u的四次多项式,g(u)为关于u的五次多项式.文中给出了谱近似解的H~2模长时间先验估计的证明,并对有限时间段分析了数值解的收敛性.数值解的先验估计是理论分析的基础,但由于该方程非线性项的特殊性,得到其数值解的H0模先验估计是难点.为此本文针对非线性项%絞(u)给出了一种反对称格式,并利用该格式构造了半离散谱近似格式,进而建立了全离散格式.由于全离散格式中g(u)利用了反对称格式,所以有效地排除了该非线性项的影响,得到了 H0模的长时间先验估计.数值解H~2半模的先验估计中通常会借助于能量泛函,但该方程和H~2半模对应的能量泛函难以找到.由于有了H0模的先验估计,所以本文在H~2半模长时间先验估计中,对非线性项的影响,采用充分地利用四阶导数项和H0模的估计,并多次使用Nirenberge不等式的分析方法,得到了较好的理论结果.
[Abstract]:In this paper , we study the finite element method for solving the nonlinear term . In this paper , we study the finite element method for solving the nonlinear term by using the finite element method . In this paper , we have given a priori estimate of the H ~ 2 half - norm by means of the energy functional function .

【学位授予单位】：吉林大学
【学位级别】：博士
【学位授予年份】：2017
【分类号】：O241.82

【相似文献】

相关期刊论文 前10条

1 林志瑗;场波问题的一些数值解法[J];电子器件;1980年03期

2 路新玲;张启敏;;固定资产投资系统数值解的收敛性[J];宁夏大学学报(自然科学版);2007年04期

3 何剑;;随机Kuramoto-Sivashinsky方程数值解的收敛性[J];太原师范学院学报(自然科学版);2009年04期

4 林鹏程;关于二维薄片振动问题的数值解[J];福州大学学报;1961年00期

5 耿兆铨;非恒定水流数值解空间失稳问题分析[J];水利学报;1985年07期

6 崔耀东 ,林兆荣;钣料压延数值解[J];南京航空航天大学学报;1986年01期

7 杨年甫;张质文;;光贴片实验全场应力的数值解法[J];实验力学;1991年01期

8 徐永权;微分方程的数值解法[J];天津理工学院学报;1996年02期

9 闵弘,章光,王水林,王元汉;一种提高数值解精度的方法探讨[J];岩石力学与工程学报;2002年01期

10 项林川;一个追线问题的数值解法及推广[J];大学物理;2004年08期

相关会议论文 前10条

1 顾本立;顾宁;;二维逆散射成象数值解[A];1987年全国微波会议论文集（上）[C];1987年

2 石文静;孟庆春;;含穿透分层双金属梁数值解与片条合成能量解法的比较[A];中国力学学会学术大会'2009论文摘要集[C];2009年

3 杨吉生;徐立军;孙天宇;叶声华;;非常规形体散射近场的数值解法[A];1995年全国微波会议论文集（上册）[C];1995年

4 吴强;;GSM(1，1)模型的结构及数值解[A];数学及其应用文集——中南模糊数学和系统分会第三届年会论文集（上卷）[C];1995年

5 张智勇;陈予恕;;VRKF方法在求解某类非线性常微分方程数值解中的应用[A];第九届全国动力学与控制学术会议会议手册[C];2012年

6 谭建国;王洪礼;;泊松白噪声激励的随机动力学系统的数值解法[A];第八届全国动力学与控制学术会议论文集[C];2008年

7 张天元;莫俊永;;核三体Faddeev方程数值解[A];第八届全国核物理会议文摘集（上册）[C];1991年

8 刘心铸;吴宣志;董英君;;直流稳定电场的数值解法及其应用[A];中国地质科学院矿床地质研究所文集（2）[C];1981年

9 翟俊辉;陈忠安;;平面弹性接触问题的一种数值解法[A];第14届全国结构工程学术会议论文集（第一册）[C];2005年

10 卢文全;;关于光线方程数值解的精度问题[A];面向21世纪的科技进步与社会经济发展（上册）[C];1999年

相关博士学位论文 前9条

1 刘凤楠;若干四阶非线性偏微分方程的数值解法[D];吉林大学;2017年

2 梁芬;一些反问题的数值解法研究[D];汕头大学;2010年

3 曾光;椭圆型方程的数值解法与稳定性分析[D];电子科技大学;2011年

4 邵和助;Schr（？）dinger方程的数值解法[D];上海大学;2010年

5 邹战勇;一类离散HJB方程的数值解法[D];湖南大学;2008年

6 赵晓朋;高阶非线性抛物方程解的性质及其数值解法[D];吉林大学;2013年

7 孙越;广义椭球函数解析解与数值解的研究[D];北京邮电大学;2015年

8 陈宇;若干反问题模型的数值求解[D];上海交通大学;2009年

9 白冬梅;非线性Schr(?)dinger方程及其相关耦合方程的数值解法[D];南京航空航天大学;2010年

相关硕士学位论文 前10条

1 左一凡;无限时区正倒向随机微分方程及永续债券定价的数值解法[D];复旦大学;2014年

2 谢梦玲;一类非线性多延迟偏微分方程的数值解法[D];华中科技大学;2014年

3 穆锦荣;改进的数值方法及简支梁问题的数值解[D];青海民族大学;2016年

4 吴晓婷;木材热压及干燥过程的两种数值方法[D];东北林业大学;2016年

5 李盘润;随机固定资产模型数值解[D];北方民族大学;2016年

6 李运利;基于Bernstein多项式和样条函数的高阶常微分方程数值解法研究[D];合肥工业大学;2010年

7 白雪;关于Gilson-Pickering方程数值解及分支研究[D];哈尔滨工业大学;2012年

8 王倩;年龄相关的单一种群动态模型数值解法[D];新疆大学;2005年

9 李希文;马尔可夫调制随机脉冲微分方程数值解的收敛性与稳定性[D];哈尔滨工业大学;2009年

10 周曼;一类热导率的反演数值解法[D];浙江大学;2007年

，

本文编号：1464793