收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

偏微分方程数值计算及增量未知元方法研究

杨爱利  
【摘要】: 本文基于微分方程的有限差分技术以及一致网格增量未知元方法,分别对一维和二维具有时间依赖系数的热方程以及一类一般的三维对流扩散方程进行了不同的研究。由于一致网格增量未知元方法可以很好地降低矩阵条件数,所以该方法的优越性在我们的理论分析和数值实验中都很好地体现了出来。非一致网格作为一种更为灵活的形式,对于许多问题,特别是边界层问题的求解,有着一致网格所无法比拟的优势。相应地非一致网格上的增量未知元方法便自然地引起了我们的注意。所以在本文的后半部分,我们从理论方面,以Dirichlet问题为例,对一维和二维非一致网格增量未知元方法下的系数矩阵条件数进行了详细得分析,并用数值算例对我们的理论分析进行了验证。 对于具有时间依赖系数的一维热方程,我们提出了一类增量未知元型半隐θ格式,仔细分析了这一新格式的稳定性,并对其误差进行了估计。结果显示,当θ趋近于1/2时,格式的稳定性条件会明显改善.而当θ趋近于零时,我们得到了一个能帮助我们恢复初始格式误差的条件。对于二维情形,我们依然用有限差分进行离散,构造了一类交替方向增量未知元型半隐(ADIUSI)格式。并在傅立叶方法的帮助下,对格式的稳定性进行了详细的分析。数值试验验证了理论分析的正确性,其结果表明这一新的格式,在某些问题上,会比经典的交替方向格式更为有效。 对于一类一般形式的三维对流扩散方程,在有限差分和增量未知元方法下,可以得到一个增量未知元型正定但非对称的线性方程组。其系数矩阵条件数要远远优于不用增量未知元方法的情形。考虑到该方法的这一优点,我们在文中将其与几种经典的迭代方法相结合,来求解上述线性系统.并从理论上对该系统的增量未知元型系数矩阵条件数进行了估计,然后通过数值试验验证了这几种增量未知元型迭代方法的有效性。 我们注意到上述的差分离散和增量未知元方法都是在一致网格上进行的,然而对于许多问题,例如边界层问题、流体力学问题等,一致网格上的差分离散已不能满足它们求解精度的需要。很自然地,要考虑非一致网格上的差分离散以及相应地非一致网格增量未知元方法。随之而来的问题是,这种非一致网格增量未知元方法能否也像一致网格增量未知元方法那样,能有效降低系数矩阵条件数?我们在本文后半部分对于该方法的一维和二维情形都进行了详细分析。理论结果表明该方法依然可以很有效地降低矩阵条件数,数值试验结果与我们的理论分析完全吻合。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 牛胜利;罗旭东;王建国;乔登江;;带电粒子在辐射带中的运动轨迹模拟[J];核技术;2011年05期
2 武刚;许永忠;杨芝林;陈世仲;赵吉祥;;地质异常体有限差分数值模拟与精度分析[J];能源技术与管理;2011年03期
3 周刚;郑丙辉;胡德超;雷坤;乔飞;;正交曲线坐标系二维浅水方程ELADI有限差分方法[J];水科学进展;2011年04期
4 李少华;王彦宾;;松软覆盖层对隐伏断层带围陷波的影响研究[J];高原地震;2011年02期
5 房克照;邹志利;刘忠波;;沙坝海岸上裂流的数值模拟[J];水动力学研究与进展A辑;2011年04期
6 李军锋;李国明;李如山;;地震波传播数值模拟中的频散校正补偿[J];价值工程;2011年19期
7 张超;杨春和;;有限差分强度折减法求解边坡稳定性[J];土木工程与管理学报;2011年02期
8 狄振华;罗振东;谢正辉;王爱文;;二维非饱和土壤水流方程基于POD方法的有限元格式[J];北京交通大学学报;2011年03期
9 丁虎;张国策;陈立群;;超临界轴向运动梁横向非线性振动频域分析[J];固体力学学报;2011年04期
10 尹丽蓉;余爱晖;;奇异两点边值问题的四次样条解[J];杭州师范大学学报(自然科学版);2011年04期
11 ;[J];;年期
12 ;[J];;年期
13 ;[J];;年期
14 ;[J];;年期
15 ;[J];;年期
16 ;[J];;年期
17 ;[J];;年期
18 ;[J];;年期
19 ;[J];;年期
20 ;[J];;年期
中国重要会议论文全文数据库 前10条
1 张伟;陈晓非;;有限差分方法中存在复杂形状地表时的自由表面条件[A];中国地球物理.2003——中国地球物理学会第十九届年会论文集[C];2003年
2 马汉东;袁宏生;;非定常喷流有限差分方法数值试验研究[A];第十届全国计算流体力学会议论文集[C];2000年
3 杨顶辉;牟永光;;二维TI介质中PSV波方程的高分辨率近似解析离散方法[A];1997年中国地球物理学会第十三届学术年会论文集[C];1997年
4 刘培国;刘克成;何建国;尹家贤;张光甫;毛均杰;;一种超宽带平面天线的时域分析[A];1999年全国微波毫米波会议论文集(下册)[C];1999年
5 杜守继;孙钧;职洪涛;翁慧俐;;软岩隧道锚杆支护机理的数值解析[A];锚固与注浆新技术——第二届全国岩石锚固与注浆学术会议论文集[C];2002年
6 张晓;宋晏蓉;李辉辉;张鹏;田金荣;张新平;;掺镱双包层光纤放大器的理论计算[A];2010年西部光子学学术会议摘要集[C];2010年
7 谢小碧;姚振兴;;二维不均匀介质中点源P-SV波响应的有限差分近似算法[A];中国科学院地球物理研究所论文摘要集(1988)[C];1989年
8 叶芳;褚春雷;高金耀;;广义差分法及其在地震正演模拟中的应用[A];中国地球物理学会第二十届年会论文集[C];2004年
9 杜守继;孙钧;;软岩地区公路隧道复合支护结构的FLAC解析[A];盛世岁月——祝贺孙钧院士八秩华诞论文选集[C];2006年
10 杜守继;翁慧俐;孙钧;职洪涛;;公路隧道复合支护机理的数值解析[A];2001年全国公路隧道学术会议论文集[C];2001年
中国博士学位论文全文数据库 前10条
1 杨爱利;偏微分方程数值计算及增量未知元方法研究[D];兰州大学;2008年
2 黎洪;稠油油藏地质特征及注蒸汽开采的有限差分方法研究[D];中国矿业大学(北京);2002年
3 王晓秋;近场光学显微成像数值模拟研究[D];大连理工大学;2004年
4 梅正阳;相机权益评价的数值算法研究[D];华中科技大学;2006年
5 庄平辉;分数阶动力方程的数值方法及其理论分析[D];厦门大学;2008年
6 胡晓毅;工业数学中关于热传导方程应用的若干反问题及其数值解法[D];复旦大学;2009年
7 詹毅;复杂有耗色散地层中的FDTD方法以及在冲击探地雷达中的应用[D];西安电子科技大学;2000年
8 徐振礼;非线性人工边界条件[D];中国科学技术大学;2007年
9 刘少斌;等离子体覆盖目标的电磁特性及其在隐身技术中的应用[D];国防科学技术大学;2004年
10 田敏;时间依赖偏微分方程的区域分解并行算法[D];山东大学;2007年
中国硕士学位论文全文数据库 前10条
1 仝海霞;我国可转换债券的定价理论与实证分析[D];首都经济贸易大学;2005年
2 黄师化;基于复合数值方法的水流动画模拟[D];安徽大学;2005年
3 刘怡旻;最优股指期货套利策略[D];厦门大学;2008年
4 苏家煜;二维多孔材料声学性能研究[D];北京交通大学;2008年
5 孙惠;一维依赖于时间的非线性Schr(?)dinger方程的有限差分方法[D];湘潭大学;2008年
6 刘勇;基于两因素的可转换债券定价模型模拟与实证研究[D];山东大学;2007年
7 李明启;基于偏微分方程的曲面修补[D];郑州大学;2009年
8 柏琰;对称正则长波(SRLW)方程的有限差分方法[D];南京航空航天大学;2005年
9 巴塞;亚式期权定价的差分方法[D];东北大学;2005年
10 潘沙;预处理方法在全速度及化学非平衡流场数值模拟中的应用[D];国防科学技术大学;2004年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978