发展方程的高精度有限元方法研究-计算数学专业论文.docx

原创性声明 本人郑重声明：所呈交的学位论文, 是本人在导师的指导下,独立进行研究所取 得的成果. 除文中已经注明引用的内容外, 本论文不包含任何其他个人或集体已经 发表或撰写过的科研成果. 对本文的研究作出重要贡献的个人和集体, 均已在文中 以明确方式标明. 本声明的法律责任由本人承担. 学位论文作者:  日期: 年 月 日 学位论文使用授权声明 本人在导师指导下完成的论文及相关的职务作品,知识产权归属郑州大学. 根据 郑州大学有关保留、使用学位论文的规定, 同意学校保留或向国家有关部门或机构 送交论文的复印件和电子版,允许论文被查阅和借阅；本人授权郑州大学可以将本 学位论文的全部或部分编入有关数据库进行检索, 可以采用影印、缩印或者其他复 制手段保存论文和汇编本学位论文. 本人离校后发表、使用学位论文或与该学位论 文直接相关的学术论文或成果时, 第一署名单位仍然为郑州大学. 保密论文在解密 后应遵守此规定. 学位论文作者:  日期: 年 月 日 摘 要 本文主要针对几类发展方程 (诸如非线性 Schro¨dinger 方程、 Benjamin-Bona- Mahony (BBM) 方程、非定常不可压缩 Navier-Stokes 方程、Cahn-Hilliard (CH) 方 程以及对流占优扩散方程), 分别从非协调有限元方法、协调和非协调混合元方法出 发, 对其收敛性、超逼近和超收敛等进行了深入系统的研究. 首先, 将一个非协调四边形单元 (改进类 Wilson元) 运用于求解非线性 Schr- o¨dinger 方程. 通过该单元相容误差在能量模意义下可以达到 O(h3) 阶, 比其插值 误差高两阶这一特殊性质, 对于半离散以及两种全离散格式 (Backward-Euler (B-E) 和 Crank-Nicolson (C-N) 格式), 在广义矩形网格下导出了最优误差估计和超逼近性 质. 进一步地, 通过插值后处理技术, 在矩形网格下得到了整体超收敛结果. 最后给 出数值算例来验证理论的正确性和方法的有效性. 其次, 研究了 BBM 方程的低阶非协调有限元方法及混合元新格式. 一方面, 利 用 EQrot 2 1 元的两个特殊性质 (插值算子等价于 Ritz 投影算子, 且相容误差为 O(h ) 阶, 比插值误差高一阶), 分别在半离散和两种全离散格式 (B-E 和 C-N 格式) 下导出 其能量模意义下的超逼近和整体超收敛结果. 另一方面, 对该方程构造一种新的非 协调混合元格式, 借助于零阶 Raviart-Thomas (R-T) 元的高精度特性, 同样在半离散 以及全离散格式下得到了相关变量的超逼近和超收敛结果, 并给出数值例子检验理 论分析的正确性. 1再次, 针对非定常不可压缩 Navier-Stokes 方程提出了一个低阶非协调混合元方 法. 运用带约束的旋转 Q1 (CQrot ) 元以及分片常数 (Q0) 元分别逼近速度 ?u 和压力 p, 1 并得到了能量模下 ?u 以及 L2 模下 p 的超逼近和超收敛结果. 同时还进行了数值实 验, 所得数值结果与理论分析相吻合. 而后, 对 CH 方程建立一个新的非协调混合元分析框架. 利用一类非协调元的 特殊性质 (相容误差为 O(h2) 阶, 比插值误差高一阶) 和插值后处理技术, 分别在半 离散和 B-E 全离散格式下得到了原始变量 u 和辅助变量 p 在能量模意义下的超逼 II 摘 要 近和整体超收敛结果. 并给出数值例子来验证该方法的有效性. 最后, 针对对流占优扩散方程, 构造了一个新的特征混合元格式. 并给出关于原 始变量 u 以及辅助变量 ?p 的收敛性分析, 最终结合数值例子来验证此方法的有效性. 关键词: 发展方程, 非协调元, 混合有限元方法, 超逼近, 超收敛. Abstract In this dissertation, we focus on several evolution equations, such as nonlinear Schr o¨dinger equation, Benjamin-Bona-Mahony (BBM) equation, nonstationary incompress- ible Navier-Stokes equation, Cahn-Hilliard (CH) equation and the convection-dominated diffusion equation, study the nonconforming finite element methods, conforming and non- conforming mixed fini