非线性发展方程无网格比高精度有限元方法研究

来源 :郑州大学 | 被引量 : 0次 | 上传用户:liang6666
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本论文对于几类非线性的发展型方程(如非线性的抛物方程、非线性的Schr(o)dinger方程、非线性Sobolev方程、非线性Ginzburg-Landau方程、非线性双曲方程),从协调有限元方法、非协调有限元方法、混合有限元方法等不同角度,利用不同技巧深入系统地研究了其线性化的全离散格式的构造、无网格比约束下的超逼近和超收敛分析以及相应数值实验.  主要的创新点具体表现在解决了以下几个问题:  (1)超收敛结果对方程解的光滑性要求比较高,但构造时间离散辅助问题(即时间离散方程系统)时,在多边形区域(例如矩形)下,就无法保证其解较强模的有界性.因此我们利用了一些特殊的、不同以往的技巧,在其解空间较弱的条件下得到无网格比超收敛的结论;巧用Taylor展开式对非线性项进行处理,以保证对时间步长(τ)的阶不丢失.  (2)由于选择的全离散格式是线性化的形式,在利用数学归纳法分析第n层的结果时需要用到第n-1层的结论,我们用一个统一的系数来控制每一个时间层的结果,这也是其数学归纳法成立的关键所在.  (3)构造了非线性双曲方程新的二阶格式,以此得到无网格比超收敛结果.而以往对非线性双曲方程的无网格比研究甚至连收敛性也没有见到报道.  (4)对一些特殊的非线性发展方程,抛弃分裂误差思想,采用一些新的技巧也证明了其无网格比超收敛性。  首先,利用低阶非协调元、协调元对非线性抛物方程分别构造了线性化的CN(Crank-Nicolson)格式和线性化的BE(Backward-Euler)格式,分析了其无网格比超逼近的性质.通过在不同的时间层引入不同的时间离散方程,把误差分裂成时间误差和空间误差两部分,有技巧性地得到了其时间误差,并由此给出了时间离散方程解的H2-模有界性.避开了时间离散方程的解在矩形区域下达不到H3-模有界这一困难.利用数学归纳法,获得了每一层数值解某种模的有界性,并统一了数学归纳法每层结果中的系数.创新性地利用Taylor展开式解决了由非线性项带来的误差估计的困难,同时保证结果中时间步长(τ)的阶不丢失.在估计空间误差时,利用更为精细的估计,得到结果O(h(h+(τ)))(h为空间网格参数,(τ)为时间步长),而不是传统的意义下的O(h2).进一步地,限制非线性抛物方程的右端项仅为满足局部的Lipschitz连续时,需要把数值解Unh的有界性条件加强至L∞-模有界,在此基础上再讨论其无网格比超逼近性.  其次,对于非线性Schr(o)dinger方程,给出了其线性化的全离散格式(BE格式和CN格式).为了克服由方程本身的虚数单位i带来的分析困难,利用两个相邻时间层相减的技巧,达到了比以前参考文献更高阶的时间误差,也相对应的得到了时间离散方程解Un的更好的有界性.利用投影算子和插值算子相结合的优势得到了原始变量un在H1(Ω)模意义下的无网格比超逼近结果.一方面,投影算子的引入成功的降低了时间离散方程在矩形区域下对其解的较高光滑性要求,得到了每一层数值解的有界性,并保证了其每一层数值解的存在唯一性.另一方面,插值算子的引入让我们可以利用插值后处理技术,得到整体超收敛结果.以上两点也充分显示了插值算子和投影算子相结合在数值分析中的重要作用.  再次,利用H1-Galerkin有限元方法分别研究了非线性Sobolev方程和非线性Ginzburg-Landau方程的无网格比超逼近性质.一方面,尝试给出一个非协调有限元对(EQrot1和零阶Raviart-Thomas单元)的重要引理.注意到非线性Sobolev方程的特点,构造一个线性化的CN格式,区别于引入时间离散方程的分裂误差做法,使用新的技巧得无网格比超逼近结果.另一方面,对非线性Ginzburg-Landau方程,利用协调有限元对(双线性单元和零阶Raviart-Thomas单元),建立了它的一个线性化的CN格式,给出了无网格比的超逼近结果.  最后,针对非线性双曲方程,首次构造出一个新的线性化二阶格式,证明了其截断误差的二阶性质.通过引入相对应的时间离散方程系统,得到了其方程解的正则性,并由此证明了它对非协调单元的无网格比超逼近性质.  需要特别指出的是本论文针对上述的每一部分,都给出了相对应的数值算例,所得数值结果进一步说明了所采用的方法是高效可行的.
其他文献
粒子群优化算法(Particle Swarm Optimization Algorithm)是基于群智能的启发式算法,它是模仿鸟集体飞行觅食和鱼群的行为,通过集体之间的协作使得最终群体达到最优.虽然每个个
等周不等式是最古老最优美的几何不等式,很多数学家给出了不同的证明,等周不等式是连接几何与分析的一个桥梁.一方面,它与分析中著名的Sobolev不等式等价.随着等周不等式的发展,
提出了一种新的制备大面积超导薄膜的方法.他人的方法是:基片运动或靶体倾斜运动或激光束运动.该方法是:基片不动, 采用光学变换系统,让脉冲激光按设定的象行星运动一样的轨
建筑智能化所各种子系统的高科技产品几乎已全面进入数字化时代,并纷纷采用IP通信网络架构。面对这种转变,如何做到更加合理经济地设计智能化系统的网络承载平台,并最终过渡
为了科学、准确地评价披碱草属牧草苗期的耐盐性,试验设0(对照)、0.4%、0.8%、1.2%4个土壤盐浓度水平,测定不同盐浓度胁迫下4种披碱草属牧草(肥披碱草、披碱草、垂穗披碱草、
K.pneumoniae菌间歇发酵甘油生产1,3-丙二醇的过程是一个复杂的生物化学反应过程,反应进行中伴随着3-羟基丙醛、乙酸、乙醇等物质的生成,同时细胞内的酶及基因也指导蛋白的合成,
在中国酒吧行业,苏荷是一家做到了顶级分数的酒吧。十年前,它还偏居南宁一隅;十年之后,近百家苏荷酒吧星火燎原,席卷大江南北,影响了整整一代人。  十年前,舞台上光彩照人的明星们往往讳言自己曾经的酒吧驻唱经历;十年后,包括成青、林军、尼克、李祥祥、金池等,一个又一个通过“中国好声音”、“中国最强音”、“中国梦之声”等选秀节目脱颖而出的歌手,他们都在舞台上很自信地说出过往经历,因为,他们曾经驻唱的地方是
近年来,全球温室效应日趋严重,据国际能源署(IEA)2012年公布数据显示,从2006年至2010年这五年里,全球co2排放量每年以6亿吨递增,2011年cO:排放量创近十年新高.随着经济迅速的
A uniaxial load experiment on coal rocks at different stress rates was carried out, based on the characteristics of acoustic emission (AE) signals in cracking c
本文首先从两个连续的2×2矩阵谱问题,KdV谱问题,AKNS谱问题出发,总结它们通过dressing链方法构造出离散方程的过程.在此基础上构造了与三个连续的3×3矩阵谱问题相联系的半离