【有限體積法】(5)

首先再回顧一下之前的高斯迭代裙盾,就拿T22舉例绞愚,由于計(jì)算是按著從左向右朱巨,從下向上的方向????進(jìn)行的史翘,所以計(jì)算T22的時(shí)候,就可以拿剛剛計(jì)算出來(lái)的T12和T21來(lái)計(jì)算T22蔬崩,如果說(shuō)最后結(jié)果是收斂的恶座,那么計(jì)算出來(lái)的T22會(huì)更接近于最終值。

對(duì)比常規(guī)的Jacobi的思想沥阳,可以發(fā)現(xiàn)高斯迭代的方法有這么一個(gè)特點(diǎn):“在同一個(gè)迭代步內(nèi)跨琳,邊界條件向計(jì)算域傳播的更快⊥┖保”也正是這一個(gè)特點(diǎn)造成了高斯迭代的收斂快于Jacobi迭代脉让。

????通過(guò)公式了解一下這句話(huà),首先是Jacobi迭代:

這是Guass迭代:

首先討論Jacobi迭代下:在同一個(gè)時(shí)間步長(zhǎng)內(nèi)功炮,不同的迭代次數(shù)下溅潜,受影響的位置是“一圈一圈”向內(nèi)傳播的。

可以把這種影響想象成波薪伏!內(nèi)部節(jié)點(diǎn)因?yàn)橛谐跏紬l件的制約滚澜,開(kāi)始不受到外部邊界的影響,這種影響是從邊界向內(nèi)部傳播的嫁怀。

在分析Guass-Seidel迭代:

同時(shí)也只考慮一個(gè)時(shí)間步長(zhǎng)设捐,在一次迭代過(guò)程中計(jì)算場(chǎng)每一個(gè)點(diǎn)的溫度計(jì)算都受到了控制體左側(cè)和下側(cè)控制體的影響,也就是說(shuō)在一個(gè)迭代步內(nèi)塘淑,左側(cè)和下側(cè)的影響就傳遍了整個(gè)計(jì)算域萝招;但是上面和右面的控制體在計(jì)算域中第一次迭代的影響卻還是一步一步。

第一次迭代影響方向用粉色表示:第二次影響方向用藍(lán)色表示存捺;

第一迭代:

第二次迭代:

這個(gè)時(shí)候槐沼,我們?cè)賮?lái)考慮一下這句話(huà):“在同一個(gè)迭代步內(nèi),邊界條件向計(jì)算域傳播的更快捌治「诠常”這也就從邊界影響傳播速度的角度解釋了為什么Guass迭代要快于Jacobi迭代。

1.線(xiàn)迭代

參考之前的(1)-(2)式肖油,給出線(xiàn)迭代的公式:

從公式上開(kāi)看兼吓,這和Guass迭代差不多,影響都是從兩邊一次迭代就傳遍了整個(gè)計(jì)算域构韵,剩下兩邊只能一格一格的傳遞影響周蹭。即:左側(cè)和右側(cè)均使用新迭代步的值趋艘,上下側(cè)使用上一個(gè)迭代步的值。

因?yàn)槊恳恍猩鲜鞘褂昧诵碌降闹敌桌剩栽诿恳恍芯涂梢员硎境扇龑?duì)角矩陣瓷胧,有了三對(duì)角矩陣就意味著可以用TDMA算法,所以這個(gè)計(jì)算方法又是一個(gè)新的方法:ADI算法棚愤。

照貓畫(huà)虎搓萧,這是左右的影響一迭代步傳遍全場(chǎng),當(dāng)然一步上下全場(chǎng)也是可以的宛畦!

2.ADI迭代

這里瘸洛,在對(duì)(3)做進(jìn)一步修改:

這樣就得到了“隱式線(xiàn)迭代”。

想想也知道次和,(4)隱式線(xiàn)迭代可以做到一個(gè)迭代步內(nèi)x軸雙方向全場(chǎng)反肋,y方向單方向全場(chǎng),所以這個(gè)隱式線(xiàn)迭代也叫做“x方向隱式線(xiàn)迭代”踏施。這就有意思了石蔗,高斯因?yàn)橛袃蓚€(gè)面是一迭代步影響遍布全場(chǎng),就比Jacobi迭代快畅形,那么隱格式線(xiàn)迭代有三個(gè)面的影響單步遍全場(chǎng)养距,它的迭代絕對(duì)是快于高斯迭代的!

既然有“x方向隱式線(xiàn)迭代”日熬,那就必定有“y方向隱式線(xiàn)迭代”棍厌,它的表達(dá)式,可以寫(xiě)成:

(5)隱式線(xiàn)迭代可以做到一個(gè)迭代步內(nèi)y軸雙方向全場(chǎng)竖席,x方向單方向全場(chǎng)耘纱。

現(xiàn)在觀察(4)-(5)式子,可以發(fā)現(xiàn)不管是x方向隱式線(xiàn)迭代還是y方向隱式線(xiàn)迭代怕敬,

??????三面影響單步遍全場(chǎng)揣炕,總有那么一個(gè)方向的影響是一步一步傳遞的帘皿。譬如:x方向隱式線(xiàn)迭代东跪,上側(cè)的影響就是一步一步傳遞的;y方向隱式線(xiàn)迭代鹰溜,右側(cè)的影響就是一步一步傳遞的虽填。

現(xiàn)在這么想,如果我們?cè)谝粋€(gè)迭代步中曹动,先進(jìn)行一次x方向隱式線(xiàn)迭代(此時(shí)上側(cè)的影響還是一步一步)斋日,緊接著進(jìn)行一次y方向隱式線(xiàn)迭代(這樣就保證了上方向也可以單次迭代影響遍全場(chǎng))。把這種操作就叫做交替方向隱式線(xiàn)迭代墓陈,ADI迭代恶守。

這個(gè)ADI迭代第献,對(duì)于一個(gè)未知的問(wèn)題,也許利用x方向隱式線(xiàn)迭代需要10次兔港,利用y方向隱式線(xiàn)迭代需要20次庸毫,此時(shí),ADI迭代綜合了這兩種迭代衫樊,需要n次就可以達(dá)到收斂效果飒赃,這里n介于10和20。所以科侈,ADI迭代不能減少計(jì)算量载佳,但是對(duì)于未知問(wèn)題,它可以很好的平衡x方向隱式線(xiàn)迭代和y方向隱式線(xiàn)迭代的計(jì)算量臀栈,所以通用性較強(qiáng)蔫慧。

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市权薯,隨后出現(xiàn)的幾起案子藕漱,更是在濱河造成了極大的恐慌,老刑警劉巖崭闲,帶你破解...
    沈念sama閱讀 206,378評(píng)論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件肋联,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡刁俭,警方通過(guò)查閱死者的電腦和手機(jī)橄仍,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,356評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門(mén),熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)牍戚,“玉大人侮繁,你說(shuō)我怎么就攤上這事∪缧ⅲ” “怎么了宪哩?”我有些...
    開(kāi)封第一講書(shū)人閱讀 152,702評(píng)論 0 342
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)第晰。 經(jīng)常有香客問(wèn)我锁孟,道長(zhǎng),這世上最難降的妖魔是什么茁瘦? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 55,259評(píng)論 1 279
  • 正文 為了忘掉前任品抽,我火速辦了婚禮,結(jié)果婚禮上甜熔,老公的妹妹穿的比我還像新娘圆恤。我一直安慰自己,他們只是感情好腔稀,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,263評(píng)論 5 371
  • 文/花漫 我一把揭開(kāi)白布盆昙。 她就那樣靜靜地躺著羽历,像睡著了一般。 火紅的嫁衣襯著肌膚如雪淡喜。 梳的紋絲不亂的頭發(fā)上窄陡,一...
    開(kāi)封第一講書(shū)人閱讀 49,036評(píng)論 1 285
  • 那天,我揣著相機(jī)與錄音拆火,去河邊找鬼跳夭。 笑死,一個(gè)胖子當(dāng)著我的面吹牛们镜,可吹牛的內(nèi)容都是我干的币叹。 我是一名探鬼主播,決...
    沈念sama閱讀 38,349評(píng)論 3 400
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼模狭,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼颈抚!你這毒婦竟也來(lái)了?” 一聲冷哼從身側(cè)響起嚼鹉,我...
    開(kāi)封第一講書(shū)人閱讀 36,979評(píng)論 0 259
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤贩汉,失蹤者是張志新(化名)和其女友劉穎,沒(méi)想到半個(gè)月后锚赤,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體匹舞,經(jīng)...
    沈念sama閱讀 43,469評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 35,938評(píng)論 2 323
  • 正文 我和宋清朗相戀三年线脚,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了赐稽。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 38,059評(píng)論 1 333
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡浑侥,死狀恐怖姊舵,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情寓落,我是刑警寧澤括丁,帶...
    沈念sama閱讀 33,703評(píng)論 4 323
  • 正文 年R本政府宣布,位于F島的核電站伶选,受9級(jí)特大地震影響史飞,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜考蕾,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,257評(píng)論 3 307
  • 文/蒙蒙 一祸憋、第九天 我趴在偏房一處隱蔽的房頂上張望会宪。 院中可真熱鬧肖卧,春花似錦、人聲如沸掸鹅。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,262評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)。三九已至葵姥,卻和暖如春荷鼠,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背榔幸。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 31,485評(píng)論 1 262
  • 我被黑心中介騙來(lái)泰國(guó)打工允乐, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人削咆。 一個(gè)月前我還...
    沈念sama閱讀 45,501評(píng)論 2 354
  • 正文 我出身青樓牍疏,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親拨齐。 傳聞我的和親對(duì)象是個(gè)殘疾皇子鳞陨,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,792評(píng)論 2 345