發(fā)展方程的高精度有限元方法研究
本文選題:發(fā)展方程 切入點:非協(xié)調元 出處:《鄭州大學》2016年博士論文
【摘要】:本文主要針對幾類發(fā)展方程(諸如非線性Schr¨odinger方程、Benjamin-BonaMahony(BBM)方程、非定常不可壓縮Navier-Stokes方程、Cahn-Hilliard(CH)方程以及對流占優(yōu)擴散方程),分別從非協(xié)調有限元方法、協(xié)調和非協(xié)調混合元方法出發(fā),對其收斂性、超逼近和超收斂等進行了深入系統(tǒng)的研究.首先,將一個非協(xié)調四邊形單元(改進類Wilson元)運用于求解非線性Schr-¨odinger方程.通過該單元相容誤差在能量模意義下可以達到O(h3)階,比其插值誤差高兩階這一特殊性質,對于半離散以及兩種全離散格式(Backward-Euler(B-E)和Crank-Nicolson(C-N)格式),在廣義矩形網格下導出了最優(yōu)誤差估計和超逼近性質.進一步地,通過插值后處理技術,在矩形網格下得到了整體超收斂結果.最后給出數值算例來驗證理論的正確性和方法的有效性.其次,研究了BBM方程的低階非協(xié)調有限元方法及混合元新格式.一方面,利用EQrot1元的兩個特殊性質(插值算子等價于Ritz投影算子,且相容誤差為O(h2)階,比插值誤差高一階),分別在半離散和兩種全離散格式(B-E和C-N格式)下導出其能量模意義下的超逼近和整體超收斂結果.另一方面,對該方程構造一種新的非協(xié)調混合元格式,借助于零階Raviart-Thomas(R-T)元的高精度特性,同樣在半離散以及全離散格式下得到了相關變量的超逼近和超收斂結果,并給出數值例子檢驗理論分析的正確性.再次,針對非定常不可壓縮Navier-Stokes方程提出了一個低階非協(xié)調混合元方法.運用帶約束的旋轉Q1(CQrot1)元以及分片常數(Q0)元分別逼近速度u和壓力p,并得到了能量模下u以及L2模下p的超逼近和超收斂結果.同時還進行了數值實驗,所得數值結果與理論分析相吻合.而后,對CH方程建立一個新的非協(xié)調混合元分析框架.利用一類非協(xié)調元的特殊性質(相容誤差為O(h2)階,比插值誤差高一階)和插值后處理技術,分別在半離散和B-E全離散格式下得到了原始變量u和輔助變量p在能量模意義下的超逼近和整體超收斂結果.并給出數值例子來驗證該方法的有效性.最后,針對對流占優(yōu)擴散方程,構造了一個新的特征混合元格式.并給出關于原始變量u以及輔助變量p的收斂性分析,最終結合數值例子來驗證此方法的有效性.
[Abstract]:In this paper, several kinds of evolution equations (such as the nonlinear Schr odinger equation Benjamin-Bona Mahony BBM) equation, the unsteady incompressible Navier-Stokes equation Cahn-Hilliardsch equation and the convection-dominated diffusion equation are studied respectively from the nonconforming finite element method, the concordant and the non-conforming mixed element method.The convergence, superapproximation and superconvergence are studied systematically.Firstly, a nonconforming quadrilateral element (modified Wilson element) is applied to solve the nonlinear odinger equation.According to the special property that the compatible error of the element can reach OH3) order in the sense of energy mode, the error is two orders higher than its interpolation error.For semi-discrete and two full discrete schemes, Backward-Eulerian B-E and Crank-Nicolsonian C-N schemes, the optimal error estimates and superapproximation properties are derived under generalized rectangular meshes.Furthermore, the global superconvergence results are obtained under rectangular meshes by interpolation postprocessing technique.Finally, numerical examples are given to verify the correctness of the theory and the validity of the method.Secondly, the low order nonconforming finite element method for BBM equation and the new mixed element scheme are studied.On the one hand, by using two special properties of EQrot1 element (interpolation operator is equivalent to Ritz projection operator, and the consistent error is OH2),The results of superapproximation and global superconvergence in the sense of energy modules are derived under semi-discrete and two kinds of full discrete schemes (B-E and C-N schemes respectively), which are higher than interpolation errors.On the other hand, a new nonconforming mixed element scheme is constructed for the equation. By virtue of the high accuracy property of the zero order Raviart-Thomas-R-T element, the superapproximation and superconvergence results of the related variables are also obtained in semi-discrete and fully discrete schemes.Numerical examples are given to verify the correctness of the theoretical analysis.Thirdly, a low order nonconforming mixed element method is proposed for unsteady incompressible Navier-Stokes equations.The velocity u and the pressure p are approximated by the rotating Q 1 C Q rot1 element and the piecewise constant Q 0) element respectively. The superapproximation and superconvergence results of u and L 2 modes under the energy mode are obtained.At the same time, numerical experiments are carried out, and the numerical results are in good agreement with the theoretical analysis.Then, a new nonconforming mixed element analysis framework is established for Ch equation.By using the special properties of a class of nonconforming elements (the compatible error is OH2, which is one order higher than the interpolation error) and the interpolation post-processing technique,The results of superapproximation and global superconvergence of the original variable u and the auxiliary variable p in the sense of energy module are obtained in semi-discrete and B-E full discrete schemes respectively.A numerical example is given to verify the effectiveness of the method.Finally, a new characteristic mixed element scheme is constructed for convection-dominated diffusion equations.The convergence analysis of the original variable u and the auxiliary variable p is given. Finally, numerical examples are used to verify the validity of the method.
【學位授予單位】:鄭州大學
【學位級別】:博士
【學位授予年份】:2016
【分類號】:O241.82
【相似文獻】
相關期刊論文 前10條
1 何文明,崔俊芝,張燕;關于矩形元的超收斂性[J];湘潭大學自然科學學報;2000年04期
2 張鐵;導數小片插值恢復技術與超收斂性[J];計算數學;2001年01期
3 冷向;不完全三次非協(xié)調三角膜元的超收斂分析[J];山東工程學院學報;2001年03期
4 陳傳淼;矩形奇妙族有限元的超收斂性[J];中國科學(A輯);2002年07期
5 趙慶華,朱起定;有限元的u-強超收斂點[J];計算數學;2004年03期
6 朱起定,趙慶華;有限元超收斂新論(英文)[J];數學進展;2004年04期
7 姜子文,王述香;粘彈性方程的有限元超收斂結果[J];科學技術與工程;2005年16期
8 許艷;王仁宏;許志強;;一類超收斂數值差商公式研究[J];計算數學;2007年01期
9 趙慶華;;一種積分形式的流量重構算法的超收斂性[J];湖南大學學報(自然科學版);2008年06期
10 劉經洪;孫海娜;朱起定;;三三次長方體有限元的超收斂[J];中國科學(A輯:數學);2009年05期
相關會議論文 前3條
1 王東東;;結構振動分析的超收斂等幾何方法[A];中國力學大會——2013論文摘要集[C];2013年
2 袁駟;王枚;林永靜;袁征;;有限元(線)法超收斂應力計算的新思路[A];第十一屆全國結構工程學術會議論文集第Ⅰ卷[C];2002年
3 肖嘉;葉康生;袁駟;;線法二階方程組有限元求解的EEP超收斂算法[A];第六屆全國土木工程研究生學術論壇論文集[C];2008年
相關博士學位論文 前10條
1 吳超;Maxwell方程棱元超收斂研究[D];湘潭大學;2015年
2 廖歆;發(fā)展方程的高精度有限元方法研究[D];鄭州大學;2016年
3 孟令雄;有限元超收斂后處理技術[D];湖南師范大學;2006年
4 許傳云;一維映射中的超收斂分岔與普適性研究[D];云南大學;2014年
5 魏華yN;界面問題和Laplace-Beltrami問題中的有限元超收斂及網格生成優(yōu)化研究[D];湘潭大學;2012年
6 熊之光;插值系數有限元法的超收斂性[D];湖南師范大學;2004年
7 唐義軍;基于改進位移模式的有限元超收斂算法研究[D];湖南大學;2013年
8 魏繼東;有限元最佳超收斂后處理技術[D];湖南師范大學;2009年
9 裴麗芳;非協(xié)調有限元方法新模式及超收斂研究[D];鄭州大學;2014年
10 秦衡峰;基于CVT網格的有限元超收斂及其在自適應有限元方法中的應用[D];湘潭大學;2008年
相關碩士學位論文 前10條
1 王斌;腔體電磁散射TM模型的梯度恢復和超收斂分析[D];鄭州大學;2015年
2 蔡惠;直邊含簡支的環(huán)扇形薄板彎曲問題的哈密頓方法與有限元法[D];溫州大學;2015年
3 曾強;桿件自由振動有限元新型超收斂算法研究[D];清華大學;2015年
4 張}液,
本文編號:1722642
本文鏈接:http://sikaile.net/shoufeilunwen/jckxbs/1722642.html