亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        求解對流擴散反應(yīng)方程的一種高精度緊致差分方法

        2021-07-14 02:04:06楊苗苗葛永斌
        關(guān)鍵詞:四階對流差分

        楊苗苗, 葛永斌

        (寧夏大學 數(shù)學統(tǒng)計學院,寧夏 銀川750021)

        對流擴散反應(yīng)方程作為一種基本的運動方程,在很多科學和工程技術(shù)領(lǐng)域都有著非常廣泛地應(yīng)用,可用于描述流體的流動與傳熱、燃燒與爆炸、核工業(yè)中反應(yīng)堆的冷卻及工業(yè)生產(chǎn)中化學堆的反應(yīng)等現(xiàn)象.因此對于該類方程精確數(shù)值解的研究具有很高的理論價值和實際作用.

        求解該類方程的方法有很多種,最常用的是有限差分法.目前已經(jīng)發(fā)展了很多具有高精度并且緊致的差分格式,如針對對流擴散方程:Sun等[1]利用理查森外推法與算子插值法改進了四階緊致差分公式,使格式具有六階精度;田振夫[2]在迎風變換的基礎(chǔ)上,將對流擴散方程轉(zhuǎn)化成含源純擴散方程,并利用Hermite法推導出指數(shù)型的差分格式,具有四階精度;梁昌弘等[3]利用3個點上的未知函數(shù)值及四階Padé格式將一階導數(shù)代換的方法,借助顯式緊致格式和隱式緊致格式的思想,得到了一種性能更優(yōu)的四階混合型格式;Tian等[4]建立了一種可用于求解定常對流擴散方程的高精度指數(shù)型有限差分格式,適合大梯度問題或邊界層問題的求解;Chen等[5]在差分格式的構(gòu)造中將攝動方法與一階迎風格式相結(jié)合,得到了求解二維定常對流擴散方程的一種高精度格式;Gupta等[6]采用泰勒級數(shù)展開法得到具有四階精度的有限差分格式;Spotz[7]則利用截斷誤差余項修正的方法給出了一種多項式型的四階緊致差分格式.盡管文獻[6-7]格式的推導方法不一樣,但本質(zhì)上是同一格式;Zhao等[8]將該方法推廣到變系數(shù)的對流擴散反應(yīng)方程中,還給出了所得格式的收斂性分析;劉明會[9]基于原方程代入和四階緊致差分公式法,將二階導數(shù)代替后得到了一種新的格式.文獻[8-9]格式盡管推導過程及方法不一樣,但本質(zhì)上仍是同一格式.田振夫[10]用截斷誤差余項修正法,將余項中的三階導數(shù)和四階導數(shù)用一階和二階導數(shù)替換,代入原方程整理便得到四階精度的差分格式,將該方法推廣到二維后,則得文獻[11]中的格式;其他針對對流擴散反應(yīng)方程的數(shù)值求解方法還有,魏劍英[12]首先由常系數(shù)變易法得到3個點上的指數(shù)型差分格式,再利用源項在中心點做二階泰勒級數(shù)展開,得到求解常系數(shù)的四階差分格式;田芳等[13]借助常系數(shù)指數(shù)型高精度緊致差分格式,采用殘量修正法得到變系數(shù)對流擴散反應(yīng)方程的緊致差分格式;田芳等[14]還基于泰勒級數(shù)展開,給出了一、二階導數(shù)的表達式,構(gòu)造了非均勻網(wǎng)格上的對流擴散反應(yīng)方程的多項式型的差分格式;蘭斌等[15]采用坐標變換法將原方程由物理空間的非均勻網(wǎng)格變換為計算空間的均勻網(wǎng)格,再結(jié)合中心差分格式得到了求解對流擴散反應(yīng)方程的緊致差分格式.

        本文針對一維、二維定常對流擴散反應(yīng)方程,基于截斷誤差余項修正的方法推導建立了一種新的緊致差分格式,具有四階精度.盡管本文方法與文獻[7-11]一樣采用到了截斷誤差余項修正法,但本文僅用到了一階導數(shù)的表達式來修正未知函數(shù)中的三階和四階導數(shù)項,而沒有用到二階導數(shù)項,這樣做的目的是使格式具有更小的耗散誤差.

        1 一維問題差分格式的建立

        首先,考慮如下一維對流擴散反應(yīng)方程:

        其中a>0為擴散項系數(shù),一般是常數(shù);u(x)是未知函數(shù),p(x)是對流項系數(shù),s(x)是反應(yīng)項系數(shù),f(x)為源項,并且要求u(x)、p(x)、s(x)、f(x)具有充分的光滑性.當s(x)≡0時,模型方程(1)為對流擴散方程.

        將區(qū)間[c,d]分為N等份:x i=c+ih,i=0,1,…,N,定義網(wǎng)格步長h=(d-c)/N;中心差分算子:

        由泰勒級數(shù)展開可得:

        將(3)和(4)式代入(1)式得

        由(1)式得

        對(6)式兩邊關(guān)于x求導可得

        則(7)式為文獻[7-10]中處理三階導數(shù)項的方法,即三階導數(shù)由一階和二階導數(shù)來表示.本文為了減小耗散誤差,將(6)式代入(7)式并整理可得

        (8)式為本文處理三階導數(shù)項的方法,即u xxx i僅由一階導數(shù)來表示.這是本文格式區(qū)別于文獻[7-10]的地方.同理,對四階導數(shù)的處理如下:

        即u xxxx i也僅用到一階導數(shù)u x i,沒有如文獻[7-10]一樣用到二階導數(shù)u xx i.

        將(8)和(9)式代入(5)式并舍去 Ο(h4)得

        對上式中的對流項和反應(yīng)項系數(shù)的導數(shù)利用中心差分公式(2)離散可得由該過程可知,本文格式的截斷誤差為Ο(h4),即格式(10)在空間上可達四階精度.對于f(x)的導數(shù)項的計算,可直接代入精確解,也可用中心差分離散,并不會影響格式的整體精度.另外,注意到本文一維格式只用到3個網(wǎng)格點,形成三對角型的代數(shù)方程組,可采用追趕法進行求解.

        2 誤差分析

        下面分別對格式(10)的耗散誤差和色散誤差進行分析.

        為了簡化計算,將本文格式(10)化為常系數(shù)

        則有

        同理

        將(12)~(14)式代入(11)式得

        化簡得

        因此,本文差分格式的特征函數(shù)可表示為

        其中

        另外,為了進行對比分析,通過計算可得文獻[3]中MHOC格式的特征函數(shù)為

        其中

        文獻[8]中FOC格式的特征函數(shù)為

        其中

        圖1~4分別給出了當η=1,Pe取不同數(shù)值時,MHOC格式[3]、FOC格式[8]和本文格式的耗散性及色散性的誤差分析.不難發(fā)現(xiàn),對耗散誤差而言,文本格式比MHOC格式和FOC格式的耗散誤差均?。坏蜕⑿远?,本文格式比FOC格式要好,但不及MHOC格式.

        圖1 當Pe=10時耗散誤差Fig.1 The dissipation error when Pe=10

        圖2 當Pe=1 000時耗散誤差Fig.2 The dissipation error when Pe=1 000

        圖3 當Pe=10時色散誤差Fig.3 The dispersion error when Pe=10

        圖4 當Pe=1 000時色散誤差Fig.4 The dispersion error when Pe=1 000

        3 二維問題差分格式的建立

        考慮如下二維對流擴散反應(yīng)方程

        其中,a,b>0為擴散項系數(shù),一般是常數(shù);u(x,y)是未知函數(shù),p(x,y)和q(x,y)是對流項系數(shù),s(x,y)是反應(yīng)項系數(shù),f(x,y)為源項,并且要求u(x,y)、p(x,y)、q(x,y)、s(x,y)、f(x,y)具有充分的光滑性.

        將區(qū)間[c,d]分為N等份:x i=c+ih,y j=c+jh,i,j=0,1,…,N,定義網(wǎng)格步長h=(d-c)/N.

        將(15)式可轉(zhuǎn)化成如下2個一維形式的方程:

        顯然方程(16)和(17)從形式上看可認為是一維的.因此,可直接利用一維對流擴散反應(yīng)方程的四階精度差分方法對其進行離散,即可得在x方向有

        其中

        對(16)式中f1(x,y)兩邊關(guān)于x求偏導,可得:

        將(20)和(21)式代入(19)式整理得

        同理可得在y方向有

        其中

        于是有

        注意(15)式為一種輔助關(guān)系,由(18)和(23)式的離散整理可得

        其中

        定義中心差分算子:

        將(26)式代入(25)式進行離散:

        整理可得(15)式的離散格式,即為本文格式:

        其中

        其中,A0,A1,…,A8中的函數(shù)p、q、s及其導數(shù)均在(i,j)點處取值.關(guān)于其一階和二階導數(shù)項的計算,均采用中心差分公式進行計算.

        由推導過程可知,該格式的截斷誤差為O(h4),即格式(27)在空間上可達四階精度.另外,注意到本文二維格式只用到9個網(wǎng)格點,故所得格式為緊致格式.

        4 數(shù)值算例

        為了驗證本文一維和二維格式的精確性和有效性,現(xiàn)將以下有精確解的數(shù)值算例分別與MHOC格式[3]、Chen格式[5]、HOC格式[7]和FOC格式[8]的最大絕對誤差的數(shù)值結(jié)果進行比較.其中,將一維和二維問題的最大絕對誤差和收斂階定義如下:

        式中的uc和ue分別代表數(shù)值解和精確解,h1和h2代表不同空間步長,Error1和Error2分別對應(yīng)步長為h1和h2時的最大絕對誤差.

        問題1[3]

        該問題的精確解為

        問題2[8]

        該問題的精確解為

        表2 問題2的最大絕對誤差和收斂階Tab.2 The maximum absolute error and convergence rate for Problem 2

        問題3[7]

        該問題的精確解為

        表3 問題3的最大絕對誤差和收斂階Tab.3 The maximum absolute error and convergence rate for Problem 3

        問題4

        該問題的精確解為

        表4 問題4的最大絕對誤差和收斂階Tab.4 The maximum absolute error and convergence rate for Problem 4

        表1~4列出了針對一維問題和二維問題1~4的4種格式在不同h下的最大絕對誤差和收斂階.不難發(fā)現(xiàn),盡管以上4種格式的精度基本上都達到了四階,但是本文計算所得的最大絕對誤差要更小些,即本文格式的計算結(jié)果更接近于精確解,也具有更高的計算精度.需要說明的是,對于MHOC格式[3]、Chen格式[5]、HOC格式[7]和FOC格式[8]的計算結(jié)果,采用了Phoebe Solver軟件進行計算得到.Phoebe Solver軟件為Web版,網(wǎng)址為www.phoebesolver.com.對于問題4沒有采用Chen格式[5]和HOC格式[7]計算是因為該問題是對流擴散反應(yīng)方程,而這2篇文獻中的方程模型為對流擴散方程,不包括反應(yīng)項.因此,這2篇文獻所提格式不適用于該問題的求解.

        表1 問題1的最大絕對誤差和收斂階Tab.1 The maximum absolute error and convergence rate for Problem 1

        5 結(jié)論

        本文首先針對方程(1)提出了一種新的緊致差分格式,具有四階精度,即在空間上采取泰勒級數(shù)展開法和截斷誤差余項修正法,先利用一階導數(shù)項來表示由原方程得到的未知函數(shù)的二階導數(shù)項,再表示出誤差余項中的三階和四階導數(shù)項,最終得到了一個三點四階的緊致差分格式,格式對應(yīng)的三對角線型方程組可采用追趕法進行求解.接下來,對本文格式的耗散誤差和色散誤差進行了分析,與文獻中格式的對比結(jié)果表明本文格式具有更小的耗散誤差.然后將該方法推廣到二維,得到了求解二維對流擴散反應(yīng)方程的一種具有四階精度的緊致差分格式.最后給出了數(shù)值實驗,通過與文獻中格式數(shù)值結(jié)果的比較發(fā)現(xiàn)本文格式具有更小的計算誤差和更高的計算精度.

        致謝寧夏自治區(qū)重點研發(fā)項目(2018BEE03007)和寧夏大學研究生創(chuàng)新項目(GIP2019010)對本文給予了資助,謹致謝意.

        猜你喜歡
        四階對流差分
        四階p-廣義Benney-Luke方程的初值問題
        齊口裂腹魚集群行為對流態(tài)的響應(yīng)
        數(shù)列與差分
        基于ANSYS的自然對流換熱系數(shù)計算方法研究
        二元驅(qū)油水界面Marangoni對流啟動殘余油機理
        基于差分隱私的大數(shù)據(jù)隱私保護
        帶參數(shù)的四階邊值問題正解的存在性
        相對差分單項測距△DOR
        太空探索(2014年1期)2014-07-10 13:41:50
        差分放大器在生理學中的應(yīng)用
        基于對流項的不同非線性差分格式的穩(wěn)定性
        中文字幕人妻丝袜乱一区三区| 国产激情视频在线观看大全| 狠狠色噜噜狠狠狠8888米奇| 国产高清在线精品一区二区三区 | 亚洲一区二区视频免费看| 亚洲av午夜一区二区三| 风韵多水的老熟妇| 国产99re在线观看只有精品| 日韩乱码精品中文字幕不卡| 少妇激情一区二区三区99| 曰本无码人妻丰满熟妇啪啪| 精品人妻无码视频中文字幕一区二区三区| 国产免费三级三级三级| 国产一级黄色录像大片| 亚洲精品久久区二区三区蜜桃臀| 欧美精品一区二区性色a+v| 日韩极品视频在线观看免费| 一区二区三区免费观看日本| 国产精品永久免费| 日韩精品中文字幕无码专区| 一区二区久久精品66国产精品| av黄页网国产精品大全| 亚洲欧美一区二区三区在线| 91久久精品无码人妻系列 | 亚洲午夜久久久久中文字幕久| 久久九九精品国产不卡一区| 少妇爆乳无码专区| 精品无码国产污污污免费网站| 久久久精品国产亚洲麻色欲| 国产性虐视频在线观看| 亚洲国产成人精品无码一区二区| 精品 无码 国产观看| 亚洲av手机在线一区| 久久精品国产亚洲av无码偷窥| 久久丫精品国产亚洲av| 色婷婷亚洲十月十月色天| 中文字幕一区二区精品视频| 免费a级毛片永久免费| 久久中文字幕亚洲精品最新 | 日韩在线精品视频一区| 任你躁国产自任一区二区三区|