二維雙曲型方程的分組并行格式

1、二維雙曲型方程的分組并行格式論文導(dǎo)讀:：本文構(gòu)造了求解二維雙曲型方程的初邊值問(wèn)題的一組分組并行算法(GE、GEL、GER),格式的局部截?cái)嗾`差階一般為,穩(wěn)定性條件為.數(shù)值例子驗(yàn)證了理論結(jié)果.論文關(guān)鍵詞：二維雙曲型方程，分組顯式格式:穩(wěn)定性:截?cái)嗾`差0 引言設(shè)數(shù)學(xué)模型為： 由于對(duì)此方程的計(jì)算具有極強(qiáng)的方向性且僅具有單邊邊界條件,故對(duì)二維雙曲型串行差分格式的并行化是一件很不容易的事情.從已有文獻(xiàn)看(見(jiàn)文【1】、【2】、【3】、【4】),尚未發(fā)現(xiàn)二維雙曲型方程的并行化格式.本文利用一個(gè)二維顯格式、兩個(gè)二維顯隱格式和一個(gè)二維隱格式構(gòu)造了一組分組顯式格式,格式的局部截?cái)嗾`差階一般為,穩(wěn)定性條件為.1 構(gòu)

2、造并行差分格式設(shè)問(wèn)題其中于是在第時(shí)間層上的4個(gè)函數(shù)值可由第時(shí)間層上的8個(gè)函數(shù)值顯式地進(jìn)行計(jì)算：即：其中：；(8)將式(3)、(4)、(5)、(6)分別在、處進(jìn)行Taylor級(jí)數(shù)展開(kāi)(見(jiàn)文【5】)得它們的局部截?cái)嗾`差分別為: (9) (10) (11) (12) 其中:.我們可將每一時(shí)間層上的節(jié)點(diǎn)除左邊界點(diǎn)和下邊界點(diǎn)外按、其中：四個(gè)相鄰點(diǎn)組成一組,共分為組.并對(duì)每組都使用GE格式(2).且可顯式表示為：(13)其中：；(14)1.2 GEL格式當(dāng)為奇數(shù)時(shí),在靠近下邊界的每?jī)蓚€(gè)內(nèi)點(diǎn)、組成一組,采用式(2)中的第一式和第二式;在靠近左邊界的每?jī)蓚€(gè)內(nèi)點(diǎn)、組成一組,采用式(2)中的第一式和第三式;在點(diǎn)采

3、用式(2)中的第一式;在其余個(gè)節(jié)點(diǎn)處反復(fù)使用GE格式(2),就得GEL格式,其矩陣形式為:(15)其中：;； 、 、同前所述. (16) 1.3 GER格式當(dāng)為奇數(shù)時(shí),在右邊界上的每?jī)蓚€(gè)內(nèi)點(diǎn)、其中：組成一組,采用式(2)中的第一式和第三式;在上邊界上的每?jī)蓚€(gè)內(nèi)點(diǎn)、(其中：)組成一組,采用式(2)中的第一式和第二式;在點(diǎn)采用式(2)中的第一式;在其余個(gè)節(jié)點(diǎn)處反復(fù)使用GE格式(2),就得GER格式,其矩陣形式為:(17)其中:;同前所述. (18)由(9)(12)式可得：定理1 當(dāng)且時(shí),式(13) 、式(15) 、式(17)的精度一般為.2 穩(wěn)定性分析我們首先對(duì)GE格式的穩(wěn)定性進(jìn)行分析，由式(13

4、)的增長(zhǎng)矩陣的特征方程可得：(19)那么：(20)從而：(21)其中：為特征值 (22)因是一個(gè)下三角行列式,故顯然有:且解得:又當(dāng)時(shí), ,從而:遞推可知：對(duì)有界.從而可得:定理2 GE格式(13)的穩(wěn)定性條件為類似可得:定理3 當(dāng)取時(shí),GEL格式(15)、GER格式(17)是穩(wěn)定的.3 數(shù)值例子考慮二維雙曲型方程初邊值問(wèn)題: (23) 它的精確解為:.以下數(shù)值例子取對(duì)GE格式取，對(duì)GEL、GER格式取.表1 并行格式GE、GEL、GER的誤差() 格式 誤差 (0.4,0.4) (0.4,0.8) (0.4,1.2) (0.8,0.4) (0.8,0.8) (0.8,1.2) (1.2,0.

5、4) (1.2,0.8) (1.2,1.2) GE 絕對(duì)誤差 2.935e-002 3.279e-002 2.345e-002 3.279e-002 3.727e-002 1.694e-002 2.345e-002 1.694e-002 1.601e-002 相對(duì)誤差 3.031e-002 3.818e-002 3.774e-002 3.818e-002 6.133e-002 6.062e-002 3.774e-002 6.062e-002 1.930e-002 GEL 絕對(duì)誤差 2.337e-002 2.788e-002 2.026e-002 2.788e-002 3.397e-002 1.

6、582e-002 2.026e-002 1.582e-002 1.493e-002 相對(duì)誤差 2.399e-002 3.228e-002 3.243e-002 3.228e-002 5.560e-002 5.638e-002 3.243e-002 5.638e-002 1.777e-001 GER 絕對(duì)誤差 1.648e-002 1.997e-002 1.695e-002 9.926e-003 3.196e-003 8.515e-004 4.054e-002 4.626e-002 5.072e-002 相對(duì)誤差 1.679e-002 2.291e-002 2.700e-002 1.101e-002 4.931e-003 2.865e-003 5.915e-002 1.350e-001 2.274e-001 以上數(shù)值例子驗(yàn)證了理論分析的正確性,說(shuō)明了本文的格式可行性與有效性.參考文獻(xiàn)：【1】EVANS D J,ABDULLAH A R B.Group explicitmethod for p