拋物型方程的并行有限差分方法
發(fā)布時間:2021-07-07 15:24
拋物型偏微分方程,簡稱拋物型方程是一類重要的偏微分方程,在自然科學諸多領域,許多現(xiàn)象都是利用拋物型方程(方程組)描述的,例如粒子和能量的擴散,物化反應,種群遷徙,物質(zhì)相互作用等等,而且拋物型方程在工程領域也有著廣泛的應用。目前對于拋物型方程發(fā)展了許多行之有效的數(shù)值方法,其中有限差分方法是最早為科技工作者運用且理論較為完善的方法,它已成為求解拋物型方程的一種重要方法。隨著大型計算機(并行機)發(fā)展,傳統(tǒng)的有限差分方法求解拋物型方程暴露出許多不足之處,例如顯式格式計算步長受到嚴格的限制,隱式格式需要求解聯(lián)立方程組,不便于直接在并行機上應用。因此構造具有并行性,良好的穩(wěn)定性和計算精度的新型有限差分方法是本文主要研究工作。本文研究內(nèi)容可以分為五個部分。第一章主要介紹本文的研究背景、研究現(xiàn)狀和文章的撰寫結構安排。第二章主要構造了求解熱傳導方程的并行有限差分方法,新的并行算法由兩個區(qū)域分解算法組成,當?shù)趎時間層解已知,利用兩區(qū)域分解算法分別計算第(n+1)時間層數(shù)值解,然后對所得到的兩個數(shù)值解進行求平均值,令這個平均值為第(n+1)時間層數(shù)值解。相比傳統(tǒng)的并行算法,新算法在保證并行本性和穩(wěn)定性的同時...
【文章來源】:武漢大學湖北省 211工程院校 985工程院校 教育部直屬院校
【文章頁數(shù)】:89 頁
【學位級別】:博士
【部分圖文】:
圖2.2:?i?=?0.4時刻精確解和數(shù)值解比較??
?y,t)?=?etsin(-Kx)sin(7ry).?(3.42)??圖3.1記錄了?Algorithm?4,?C-N格式和文獻[60]中高階ADI算法在網(wǎng)格比r?=?1.2時間為??t?=?0.5時,在取不同空間內(nèi)點數(shù)的CPU計算時間。結果表明,Algorithm?4具有良好的并行效??率。圖3.2給出了例3.1.2在t?=?0.4時刻下的精確解以及針對網(wǎng)格比r*?=?1分別取不同空間步長??/I,在t?=?0.4時刻下的數(shù)值解,我們可以看出Algorithm?4求解例3.1.2得到的數(shù)值解是令人滿??意的。表3.8給出了在擴散系數(shù)與對流系數(shù)�。�?=?6?=?1、空間步長h?=?1/25時Algorithm?4和??Zhu?[35]兩算法的范數(shù)誤差比較
?72??Number?of?space?grid?interior?points??圖3.1:三種算法的CT>/7計算時間??考慮二維對流占有擴散問題如下:??例?3.1.2??—-V???(aVu)?+?V???(bu)?=?f(x,?y,?t),?(x,?y)?e?n,?t?e?(0,1],??u[x,?y,?Q)?=?sin[nx)sin〔7ry),?(x,?y)?6?fl,??'?(3.41)??u(0,?y,?t)?=?u(l,?y,t)=0,?ye?[0,1],?f?G?(0,1],??u(x,o,?t)?=?u(x,?l,?t)?=?0,?x?e?[o,?i],?t?e?(o;?l],??其中求解區(qū)域n?=?K〇,?:〇丨2,并且右端項為??f(x,t)?-?e1?\2a-K2sin{i:x)sin{-Ky)?+?b-KC〇s{-nx)sin(ny)?+?b-irsin{'Kx)cos{/Ky)?+?sin(irx)sin(Try)].??算例3.1.2的精確解為??u(x,?y,t)?=?etsin(-Kx)sin(7ry).?(3.42)??圖3.1記錄了?Algorithm?4
【參考文獻】:
期刊論文
[1]An Explicit-Implicit Predictor-Corrector Domain Decomposition Method for Time Dependent Multi-Dimensional Convection Diffusion Equations[J]. Liyong Zhu~(1,2),Guangwei Yuan~1 and Qiang Du~(3,*) 1 Laboratory of Computational Physics,Institute of Applied Physics and Computational Mathematics,Beijing,China. 2 Division of Mathematical Sciences,School of Physical and Mathematical Sciences,Nanyang Technological University,Singapore,637371. 3 Department of Mathematics,Pennsylvania State University,University Park, PA 16802,USA.. Numerical Mathematics:Theory,Methods and Applications. 2009(03)
[2]求解Burgers方程的一種交替分段隱式算法[J]. 王廷春,張魯明. 計算物理. 2005(02)
[3]Burgers方程的交替分組顯式迭代方法[J]. 金承日,劉家琦. 計算物理. 1998(05)
[4]求解對流-擴散方程的交替分段顯-隱式方法[J]. 陸金甫,張寶琳,徐濤. 數(shù)值計算與計算機應用. 1998(03)
[5]非線性拋物組具并行本性的某些實用差分格式[J]. 周毓麟,沈隆鈞,袁光偉. 數(shù)值計算與計算機應用. 1997(01)
[6]非線性拋物組具并行本性的一般差分格式[J]. 周毓麟,袁光偉. 中國科學(A輯 數(shù)學 物理學 天文學 技術科學). 1997(02)
[7]帶非線性邊界條件的熱傳導方程的整體解與爆破問題[J]. 王術,王新明. 應用數(shù)學. 1997(01)
[8]拋物型方程有限差分并行解法[J]. 張寶琳,陳勁. 數(shù)值計算與計算機應用. 1995(03)
[9]拋物型方程的交替分段Crank—Nicholson格式[J]. 張寶琳,蘇秀敏. 計算物理. 1995(01)
[10]求解擴散方程的交替分段顯-隱式方法[J]. 張寶琳. 數(shù)值計算與計算機應用. 1991(04)
本文編號:3269879
【文章來源】:武漢大學湖北省 211工程院校 985工程院校 教育部直屬院校
【文章頁數(shù)】:89 頁
【學位級別】:博士
【部分圖文】:
圖2.2:?i?=?0.4時刻精確解和數(shù)值解比較??
?y,t)?=?etsin(-Kx)sin(7ry).?(3.42)??圖3.1記錄了?Algorithm?4,?C-N格式和文獻[60]中高階ADI算法在網(wǎng)格比r?=?1.2時間為??t?=?0.5時,在取不同空間內(nèi)點數(shù)的CPU計算時間。結果表明,Algorithm?4具有良好的并行效??率。圖3.2給出了例3.1.2在t?=?0.4時刻下的精確解以及針對網(wǎng)格比r*?=?1分別取不同空間步長??/I,在t?=?0.4時刻下的數(shù)值解,我們可以看出Algorithm?4求解例3.1.2得到的數(shù)值解是令人滿??意的。表3.8給出了在擴散系數(shù)與對流系數(shù)�。�?=?6?=?1、空間步長h?=?1/25時Algorithm?4和??Zhu?[35]兩算法的范數(shù)誤差比較
?72??Number?of?space?grid?interior?points??圖3.1:三種算法的CT>/7計算時間??考慮二維對流占有擴散問題如下:??例?3.1.2??—-V???(aVu)?+?V???(bu)?=?f(x,?y,?t),?(x,?y)?e?n,?t?e?(0,1],??u[x,?y,?Q)?=?sin[nx)sin〔7ry),?(x,?y)?6?fl,??'?(3.41)??u(0,?y,?t)?=?u(l,?y,t)=0,?ye?[0,1],?f?G?(0,1],??u(x,o,?t)?=?u(x,?l,?t)?=?0,?x?e?[o,?i],?t?e?(o;?l],??其中求解區(qū)域n?=?K〇,?:〇丨2,并且右端項為??f(x,t)?-?e1?\2a-K2sin{i:x)sin{-Ky)?+?b-KC〇s{-nx)sin(ny)?+?b-irsin{'Kx)cos{/Ky)?+?sin(irx)sin(Try)].??算例3.1.2的精確解為??u(x,?y,t)?=?etsin(-Kx)sin(7ry).?(3.42)??圖3.1記錄了?Algorithm?4
【參考文獻】:
期刊論文
[1]An Explicit-Implicit Predictor-Corrector Domain Decomposition Method for Time Dependent Multi-Dimensional Convection Diffusion Equations[J]. Liyong Zhu~(1,2),Guangwei Yuan~1 and Qiang Du~(3,*) 1 Laboratory of Computational Physics,Institute of Applied Physics and Computational Mathematics,Beijing,China. 2 Division of Mathematical Sciences,School of Physical and Mathematical Sciences,Nanyang Technological University,Singapore,637371. 3 Department of Mathematics,Pennsylvania State University,University Park, PA 16802,USA.. Numerical Mathematics:Theory,Methods and Applications. 2009(03)
[2]求解Burgers方程的一種交替分段隱式算法[J]. 王廷春,張魯明. 計算物理. 2005(02)
[3]Burgers方程的交替分組顯式迭代方法[J]. 金承日,劉家琦. 計算物理. 1998(05)
[4]求解對流-擴散方程的交替分段顯-隱式方法[J]. 陸金甫,張寶琳,徐濤. 數(shù)值計算與計算機應用. 1998(03)
[5]非線性拋物組具并行本性的某些實用差分格式[J]. 周毓麟,沈隆鈞,袁光偉. 數(shù)值計算與計算機應用. 1997(01)
[6]非線性拋物組具并行本性的一般差分格式[J]. 周毓麟,袁光偉. 中國科學(A輯 數(shù)學 物理學 天文學 技術科學). 1997(02)
[7]帶非線性邊界條件的熱傳導方程的整體解與爆破問題[J]. 王術,王新明. 應用數(shù)學. 1997(01)
[8]拋物型方程有限差分并行解法[J]. 張寶琳,陳勁. 數(shù)值計算與計算機應用. 1995(03)
[9]拋物型方程的交替分段Crank—Nicholson格式[J]. 張寶琳,蘇秀敏. 計算物理. 1995(01)
[10]求解擴散方程的交替分段顯-隱式方法[J]. 張寶琳. 數(shù)值計算與計算機應用. 1991(04)
本文編號:3269879
本文鏈接:http://sikaile.net/kejilunwen/yysx/3269879.html
最近更新
教材專著