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

        ?

        一類非線性Burgers方程組差分格式的計算穩(wěn)定性*

        2018-12-15 01:50:44馬永柳曹艷華
        吉首大學學報(自然科學版) 2018年5期
        關鍵詞:啟發(fā)性初值微分

        馬永柳,曹艷華

        (華北電力大學數(shù)理學院,北京 102206)

        在利用數(shù)值方法求解非線性偏微分方程時會產(chǎn)生計算不穩(wěn)定性.計算不穩(wěn)定性不同于線性不穩(wěn)定性,它是一種突變型的指數(shù)增長,且不會隨著時間步長和空間步長的改變而改變,因此在求解非線性偏微分方程時首先要分析其計算穩(wěn)定性.非線性偏微分方程初值問題的解,一般情況下只在時間T的局部范圍內(nèi)存在,即使初值充分光滑,其解的整體存在性也無法確定.相應的解在有限時間內(nèi)會失去正則性并產(chǎn)生奇解,稱為解的破裂[1-2].1968年,C W Hirt[3]提出了分析非線性偏微分方程差分格式計算穩(wěn)定性的方法,即啟發(fā)性分析方法.楊曉忠等[4]將啟發(fā)性分析方法用于判定KdV方程和Burgers方程的差分格式的計算穩(wěn)定性,并討論了三時間層和二時間層差分格式的計算穩(wěn)定性問題,得到了計算穩(wěn)定性的判斷依據(jù).Burgers方程是模擬沖擊波傳播和反射的非線性數(shù)學物理方程,被廣泛應用在流體力學、非線性聲學和氣體動力學等領域[5-9].對于非線性耦合Burgers方程組

        (1)

        曲娜等[10]僅分析了其差分格式的穩(wěn)定性而沒有考慮計算穩(wěn)定性.筆者將利用啟發(fā)性分析方法來推導該非線性耦合Burgers方程組差分格式具有計算穩(wěn)定性的必要條件.

        1 差分格式的表示

        向前差分格式為

        (2)

        選適當內(nèi)節(jié)點(xj,tk+1),則該節(jié)點處的微分方程組為

        向后差分格式為

        (3)

        2 差分格式計算穩(wěn)定性的啟發(fā)性分析

        將差分格式進行Taylor展開和自循環(huán)消元,轉(zhuǎn)化為等價的變形方程組,即修正微分方程組(Modified PDE).判斷該差分格式具有計算穩(wěn)定性的依據(jù)是其修正微分方程組等號右端的二階導數(shù)項為正耗散項,此方法稱為MPDE識別法.筆者將MPDE識別法運用于非線性耦合Burgers方程組(1),以向前差分格式(2)為例進行分析.由Taylor展開,可得

        進一步得到如下展開:

        (4)

        將Taylor展開結果(4)代入向前差分格式(2),可得

        略去上下標后可得

        (5)

        由方程組(5)可得

        對(6)式關于t求微分,可得

        (8)

        對(6)式關于x求微分,可得

        (9)

        對(7)式關于x求微分,可得

        (10)

        對(9)式關于x求微分,可得

        (11)

        將(7),(10),(11)式代入(8)式,可得

        (12)

        同理

        (13)

        將(12),(13)式代入方程組(5)可得向前差分格式(2)的修正微分方程組

        同理,向后差分格式(3)的修正微分方程組為

        由啟發(fā)性分析方法可知,只有當差分格式的修正微分方程等號右端的二階耗散項系數(shù)為正時,差分格式才具有計算穩(wěn)定性;因此有如下結論:

        定理1非線性耦合Burgers方程組(1)的向前差分格式(2)具有計算穩(wěn)定性的必要條件為

        定理2非線性耦合Burgers方程組(1)的向后差分格式(3)具有計算穩(wěn)定性的必要條件為

        3 數(shù)值實驗

        由定理1和定理2可知,非線性耦合Burgers方程組(1)的差分格式的計算穩(wěn)定性僅與u(x,t)和v(x,t)的初值及α,β,Δt有關.現(xiàn)取4個初值進一步驗證非線性耦合Burgers方程組(1)的普遍性差分格式的計算穩(wěn)定性與初值、差分格式結構之間的關系.4個初值分別為:

        取α=-1,β=2,-5≤x≤5,0≤t≤1,則非線性耦合Burgers方程組(1)化為

        取Δx=0.1,Δt=0.001,左邊界值為0,右邊界值用Lagrange一次插值方法外插.差分格式的計算穩(wěn)定性結果列于表1.

        表1 差分格式的計算穩(wěn)定性結果

        由表1可知,當α,β,Δt固定時,差分格式的計算穩(wěn)定性僅與u(x,t)和v(x,t)的初值有關,而不同的初值在x處取不同值時,就得到不同的計算穩(wěn)定性結果.

        4 結語

        利用啟發(fā)性差分格式對一類非線性耦合Burgers方程組的差分格式進行計算穩(wěn)定性分析,當差分格式的修正微分方程組等號右端的二階耗散系數(shù)大于0時,差分格式是穩(wěn)定的.非線性發(fā)展方程差分格式是否具有計算穩(wěn)定性,與方程解的性質(zhì)緊密相關,而本研究中的非線性耦合Burgers方程組差分格式的計算穩(wěn)定性,依賴于其耗散系數(shù)α,β,Δt及初值u(x,0),v(x,0)的選取.

        雖然本研究得到的判斷依據(jù)只是判斷一類非線性耦合Burgers方程組差分格式計算穩(wěn)定性的必要條件,即差分格式具有計算穩(wěn)定性必須滿足該判斷依據(jù),而滿足該判斷依據(jù)的差分格式未必具有計算穩(wěn)定性;但是,由于在計算中利用該判斷依據(jù)可以刪除不具有計算穩(wěn)定性的差分格式,避免盲目計算,因此計算穩(wěn)定性的必要條件是非常有效且實用的.

        猜你喜歡
        啟發(fā)性初值微分
        觀點,因“啟發(fā)性”而深刻
        基于啟發(fā)性提示語的高中數(shù)學解題教學策略
        具非定常數(shù)初值的全變差方程解的漸近性
        談舞蹈教學中啟發(fā)性語言的運用
        樂府新聲(2021年3期)2021-11-04 09:22:06
        擬微分算子在Hp(ω)上的有界性
        一種適用于平動點周期軌道初值計算的簡化路徑搜索修正法
        上下解反向的脈沖微分包含解的存在性
        三維擬線性波方程的小初值光滑解
        開啟學科教學的多元功能
        借助微分探求連續(xù)函數(shù)的極值點
        a级毛片在线观看| 国产精品九九热| 成人区人妻精品一区二区不卡网站 | 三个黑人插一个女的视频| 国产一级一区二区三区在线播放| 手机在线观看日韩不卡av| 亚洲精品中文字幕视频色| 国内永久福利在线视频图片| 久热re这里精品视频在线6| 亚洲av永久无码一区二区三区| 中文字幕人妻在线中字| 漂亮人妻被中出中文字幕久久| 欧洲日本一线二线三线区本庄铃| 国产成人精品三级麻豆| 亚洲熟妇av日韩熟妇av| 精品专区一区二区三区| 国产精品久久国产精品久久| 日韩精品成人一区二区三区| 亚洲女人的天堂网av| 日韩一区二区三区精品视频| 国产午夜免费高清久久影院| 少妇高潮太爽了在线视频| 亚洲日韩国产一区二区三区在线 | 亚洲一区二区丝袜美腿| 午夜精品久久99蜜桃| 成年美女黄的视频网站| 国产精品夜间视频香蕉| 性生交大全免费看| 精品国产福利久久久| 按摩女内射少妇一二三区| 国产一区二区av在线免费观看| www夜片内射视频在观看视频| 亚洲综合色无码| 久久精品国产亚洲av麻| 亚洲欧美性另类春色| 一区二区三区视频在线免费观看| 中国亚洲av第一精品| 国产老熟妇精品观看| 东京无码熟妇人妻av在线网址| 久久永久免费视频| 加勒比一本大道大香蕉|