四阶抛物问题的弱Galerkin有限元法
本文关键词:四阶抛物问题的弱Galerkin有限元法
更多相关文章: 四阶抛物方程 弱Galerkin有限元 弱函数空间 离散弱Laplace算子
【摘要】:本文主要用弱Galerkin有限元法来研究一个四阶抛物方程初边值问题的数值计算方法.考虑如下四阶抛物方程的初边值问题:ut+△2u=f,x∈Ω, 0≤t≤t,(0.1) u=au/an=0,x∈aΩ,0≤t≤t(0.2) u(·,0)=Ψ,x∈Ω. (0.3)其中△是Laplace算子,Ω是Rd(d=2,3)中有界开区域且其边界aΩ是Lipschitz连续的.令H=L2(Ω)表示平方可积函数构成的空间,具有标准的内积(·,·)和范数‖·‖.我们也用Hm=Hm(Ω)表示标准的Sobolev空间且则方程(0.1)-(0.3)的变分形式为:求u∈L2(0,t;H02(Ω))使得u(0,·)=Ψ,且满足如下方程(ut,v)+(△u,△v)=(f,v),(?)v∈H02(Ω). (0.4)本文将用弱Galerkin有限元法(简记为WG)来求解方程(0.1)-(0.3).WG有限元法的基本思路是:构造弱函数空间W(Ω)近似H2空间,然后定义一个弱Laplace算子△w用来近似标准的Laplace算子△,再利用变分方程(0.4)和适当的稳定子s(.,.)建立求解方程(0.1)-(0.3)的数值计算方法.令Th表示区域Ω的一族正则三角剖分,T是其中任意的一个三角元,其直径为h7,令h=max{hT}.对任意给定的非负整数k≥2,用Pk(T)表示T上次数不超过k的多项式集合,用Rk(e)表示边界e(?)αT上次数不超过k的多项式集合,那么离散弱函数空间Wk(T)表示如下Wk(T)={{v0,vb,vg}:v0∈Pk(T),vb∈Pk(e),vgg∈[Pk-1(e)]d,e(?)aT}.从而得到弱有限元空间Vh如下Vh={{v0,vb,vg:{v0,vb,vg}|T∈Wk(T),(?)T∈Th}.用Vh0表Vh的子空间,其函数值在aΩ上为0,即Vh0={v={v0,vb,vg}∈Vh,vb|e=0,vg·n|e=0,e(?)aT ∩aΩ}.对Vh中任意的uh={u0,ub,ug}和v={v0,vb,vg},引进双线性形式如下其中·,·)αT表示在区域边界αT上的L2内积.离散弱Laplace算子记为△w,对(?)v∈Wk(T),定义△wv∈Pk(T)满足如下方程(△wv,φ)T=(v0,△φ)T-vb,%溅铡aT+vg·n,φaT,(?)φ∈Pk(T).同时记离散内积则方程(0.1)-(0.3)的半离散弱Galerkin有限元法为:寻找一个弱函数uh∈L∞(0,t;Vh)满足uh(0)=QhΨ使得(vh,t(t,·),v0)+(△wuh,Δwv)h+s(uh,v)=(f,v0),(?)v={v0,vb,vg}∈Vh0.(0.5)令K为时间步长,tn=nk,n=1,…,N,其中NK=t.用Un∈Vh表示u(tn)的近似.则求解方程(0.4)的向后欧拉格式的弱Galerkin法为:求Un∈Vh(n≥1),满足U0=QhΨ使得(aUn,v0)+(△wUn,△wv)+s(Un,v)=(f(tn),v0),(?)v∈Vh0. (0.6)K在本文中我们得到以下半离散和全离散弱Galerkin格式的误差估计.定理0.1.令uh={u0,ub,ug}为(0.5)满足初值条件uh(0)=QhΨ的弱Galerkin有限元的解.假设方程(0.1)-(0.3)的精确解满足u∈Hmax{k+1,4}(Ω)令eh=uh-Qhu为弱Galerkin近似解和真解“的L2投影之间的误差,则存在常数C使得||eh||2+(?)01|||eh(·,t|||2dτ≤||eh(·,0)||2+Ch2k-2((?)01(||u||k=12+h2δμ,2||u||42ds),(0.7)且4(?)01||eh,τ||2dτ+|||eh|||2 ≤2|||eh(0,·)|||2+Ch2k-2(||u||k+12+||u(0,·)||k+12+hδk,2(||u||42+||u(0,·)||42) (0.8) +(?)(""uτ||k+12+hδk,2||uτ||42)dτ+(?)0(||u||k+12+hδk,2||u42)dτ).定理0.2.令u∈Hmax{1+k,4}(Ω)和Un分别为方程(0.1)-(03)和(0.6)满足U0=Qhu(t0)的解ehn=Un-qhu(tn)为全离散的弱Galerkin近似解和真解u的L2投影之间的误差,则存在常数C,使得
【关键词】:四阶抛物方程 弱Galerkin有限元 弱函数空间 离散弱Laplace算子
【学位授予单位】:吉林大学
【学位级别】:硕士
【学位授予年份】:2015
【分类号】:O241.82
【目录】:
- 摘要4-8
- Abstract8-13
- 第一章 绪论13-20
- 1.1 边值问题的Galerkin法13-16
- 1.2 Galerkin方程解的存在唯一性和误差估计16-20
- 第二章 求抛物方程的弱Galerkin有限元法20-26
- 2.1 抛物方程的有限元法20-22
- 2.2 抛物方程的弱Galerkin有限元法22-26
- 第三章 四阶抛物问题的弱Galerkin有限元法26-39
- 3.1 引言26-27
- 3.2 弱Laplace算子和离散弱Laplace算子27-28
- 3.3 半离散和全离散弱Galerkin有限元法28-30
- 3.4 误差分析30-37
- 3.4.1 半离散弱Galerkin有限元法的误差分析30-35
- 3.4.2 全离散弱Galerkin有限元法的误差分析35-37
- 3.5 数值结果37-39
- 参考文献39-40
- 致谢40
【相似文献】
中国期刊全文数据库 前10条
1 刘绍球 ,李显霖;有限元法浅说[J];系统工程与电子技术;1981年02期
2 王月f;《水资源中有限元法》——第四届国际水资源应用有限元法会议会议录介绍[J];地质科技情报;1983年04期
3 鲍四元;;本科生《有限元法》课程教学的现状与若干改革研究[J];课程教育研究;2014年03期
4 石艳华;侯录;席俊明;;有限元法在口腔医学中的应用[J];今日科苑;2007年08期
5 李昱君;;从语境的视角看有限元法的早期发展[J];自然辩证法通讯;2008年03期
6 鲁建霞;苟惠芳;;有限元法的基本思想与发展过程[J];机械管理开发;2009年02期
7 葛藤;;《有限元法》本科教学改革探讨[J];科技信息;2010年13期
8 张泉乐;薛荔;;非线性问题的有限元法[J];武汉理工大学学报;2010年15期
9 王小荣;林顺洪;丁剑平;陈小亮;;有限元法课程的教学改革探讨[J];重庆科技学院学报(社会科学版);2011年06期
10 徐建全;陈铭年;林大同;;《有限元法》本科课程综合改革与实践[J];机电技术;2011年06期
中国重要会议论文全文数据库 前10条
1 唐旭海;郑超;张建海;;基于虚结点的多边形有限元法[A];中国力学学会学术大会'2009论文摘要集[C];2009年
2 赵俊峰;周慎杰;王炳雷;王锡平;;应变梯度微梁的有限元法[A];力学与工程应用(第十三卷)[C];2010年
3 Liu Bin;Xu Ran;Zhang Yihui;Fang Daining;;原子有限元法在多物理场研究中的新进展(英文)[A];中国力学大会——2013论文摘要集[C];2013年
4 丁胜勇;邵国建;;基于多边形有限元法的混凝土细观模拟[A];中国力学大会——2013论文摘要集[C];2013年
5 王耀;黑宝平;杨在林;;用于非均匀介质瞬态分析的动态非均匀有限元法[A];中国力学大会——2013论文摘要集[C];2013年
6 金峰;方修君;;扩展有限元法及与其他数值方法的联系[A];第17届全国结构工程学术会议论文集(第Ⅰ册)[C];2008年
7 邓娇;王宏志;傅向荣;邓小环;宋孟燕;葛丽娜;田歌;;单自由度时空有限元法[A];北京力学会第19届学术年会论文集[C];2013年
8 徐春晖;李明瑞;;可离散有限元法中应变度量的选择[A];第五届全国计算爆炸力学会议论文摘要[C];2012年
9 张谧;徐春辉;;由于连续体破坏问题研究的可离散有限元法[A];北京力学会第19届学术年会论文集[C];2013年
10 张俊峰;李兴国;娄国伟;;有限元法在单片电路耦合设计中的应用[A];2007年全国微波毫米波会议论文集(下册)[C];2007年
中国重要报纸全文数据库 前1条
1 实习记者 南秀渊;有限元法研究的峥嵘岁月[N];新清华;2014年
中国博士学位论文全文数据库 前10条
1 张双文;时域有限元法及其截断边界条件的研究[D];西南交通大学;2008年
2 付振东;非均质饱和多孔介质准静态行为分析的耦合多尺度及耦合升尺度有限元法[D];大连理工大学;2010年
3 林志良;比例边界有限元法及快速多极子边界元法的研究与应用[D];上海交通大学;2010年
4 卜令方;增强有限元法及其在岩土工程非连续变形分析中的应用[D];浙江大学;2014年
5 林峰;瞬变电磁场中时空有限元法的研究[D];沈阳工业大学;2005年
6 曹凤帅;比例边界有限元法在势流理论中的应用[D];大连理工大学;2009年
7 郑恩希;非多项式最小二乘有限元法在几种散射问题中的应用[D];吉林大学;2012年
8 王辉;基于径向基函数的无网格数值方法及杂交Trefftz有限元法[D];天津大学;2007年
9 李霄琳;非均质材料的光滑多尺度有限元法研究[D];吉林大学;2015年
10 何小祥;基于有限元法的混合技术及其在复杂电大腔体电磁分析中的应用[D];南京航空航天大学;2003年
中国硕士学位论文全文数据库 前10条
1 李昱君;语境视角下的有限元法发展史[D];山西大学;2008年
2 韩冬月;四阶抛物问题的弱Galerkin有限元法[D];吉林大学;2015年
3 樊现行;光滑有限元法理论及算法研究[D];山东大学;2013年
4 万德涛;二维光滑扩展有限元法研究[D];湖南大学;2014年
5 李超;无网格—有限元法耦合理论的研究及其应用[D];西安理工大学;2008年
6 胡姝芳;非线性哈密顿系统的有限元法研究及比较[D];湖南师范大学;2009年
7 谢海;扩展有限元法的研究[D];上海交通大学;2009年
8 姚冬梅;间断伽辽金有限元法及弹塑性固体和饱和多孔介质中动力和波传播问题分析[D];大连理工大学;2002年
9 张刚;基于光滑有限元法消除体积锁定和剪切锁定问题研究[D];山东大学;2012年
10 陆云;拓扑有限元法在电磁场问题中的应用[D];广西大学;2006年
,本文编号:1088223
本文链接:https://www.wllwen.com/kejilunwen/yysx/1088223.html