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

        ?

        二維拋物型奇異攝動問題的移動網(wǎng)格方法

        2021-04-26 08:19:58
        工程數(shù)學(xué)學(xué)報 2021年6期
        關(guān)鍵詞:拋物平面數(shù)值

        周 琴

        (湖南涉外經(jīng)濟學(xué)院信息與機電工程學(xué)院,長沙 410205)

        0 引言

        二維拋物型奇異攝動問題在自然科學(xué)、工程等領(lǐng)域有著廣泛的應(yīng)用背景,它可以描述物質(zhì)擴散時的濃度變化規(guī)律,土壤力學(xué)中的滲流問題等等。該類問題的解通常在局部范圍內(nèi)會產(chǎn)生劇烈的變化,在均勻網(wǎng)格上獲取數(shù)值解不能很好地逼近實際問題。目前求解奇異攝動問題的非均勻網(wǎng)格方法主要有自適應(yīng)移動網(wǎng)格方法[?,?,?,?]和層適應(yīng)網(wǎng)格方法[5–6]等。其中自適應(yīng)移動網(wǎng)格方法通過網(wǎng)格的迭代變化使解變化較大的區(qū)域集中較多的網(wǎng)格點,能較好地反映解的性質(zhì)。

        奇異攝動問題的自適應(yīng)移動網(wǎng)格方法研究中,關(guān)于二維問題的研究成果相對較少。文獻[?]針對一維奇異攝動對流擴散方程組問題,利用有限差分方法,提出了基于等弧長分布的自適應(yīng)移動網(wǎng)格方法,并給出了一階后驗誤差估計。文獻[2–3]分別對一維奇異攝動時變對流擴散問題和一階奇異攝動初值問題進行分析,建立了基于等分布弧長的移動網(wǎng)格方法。文獻[?]提出將適合問題的層適應(yīng)網(wǎng)格作為初始網(wǎng)格來進行網(wǎng)格移動,比較了兩種移動網(wǎng)格方法。文獻[7–8]針對一階的雙曲守恒律方程提出了基于等分布原理的移動網(wǎng)格方法。文獻[?,?,?,?,?]研究了自適應(yīng)移動網(wǎng)格方法在流體力學(xué)等問題中的應(yīng)用.

        我們考慮一類含源項二維拋物型奇異攝動問題其中0<ε ?1, Ω=[a1,b1]×[a2,b2]。函數(shù)f(x,y,t)為已知的源匯函數(shù),設(shè)它滿足一定的正則性條件。本文將提出計算問題(??)的自適應(yīng)移動網(wǎng)格算法,并給出兩個算例的數(shù)值實驗結(jié)果。

        1 數(shù)值格式

        本節(jié)構(gòu)造移動網(wǎng)格上求解問題(??)的數(shù)值格式。設(shè)某時刻物理平面上網(wǎng)格分布為{(xj,k,yj,k)}, j= 0,1,··· ,N, k= 0,1,··· ,M,對應(yīng)數(shù)值解為{uj,k}。由于自適應(yīng)移動網(wǎng)格方法中的網(wǎng)格是迭代變化的,網(wǎng)格移動后產(chǎn)生的新網(wǎng)格為不規(guī)則的四邊形網(wǎng)格,不適用常規(guī)的數(shù)值方法作數(shù)值計算。因此,在某個時間層迭代的最終網(wǎng)格形成后,需要將不規(guī)則的四邊形網(wǎng)格轉(zhuǎn)化為計算平面上均勻分布的網(wǎng)格,以便于后續(xù)計算。記計算平面上均勻分布的網(wǎng)格為{(ξj,ηk)}, j= 0,1,··· ,N, k= 0,1,··· ,M。在計算過程中,計算平面上的網(wǎng)格總是保持不變,用來輔助求解問題。記

        2 網(wǎng)格迭代方法和解的重新分布

        3 移動網(wǎng)格算法和數(shù)值實驗

        基于前文構(gòu)造的數(shù)值格式及網(wǎng)格迭代時解的重新分布方法,下面給出移動網(wǎng)格方法求解模型方程(??)的算法。

        第4 步 如果tm+1≤T,令m=m+1,轉(zhuǎn)第2 步;否則,算法終止,得到T時刻的數(shù)值解。

        下面根據(jù)上述算法對兩個算例進行數(shù)值實驗。

        圖1 t=0.5, ε=0.001 時,例1 的自適應(yīng)移動網(wǎng)格分布

        圖2 為例2 中t= 0.2, ε2= 0.0001 時,空間網(wǎng)格數(shù)分別為50×50 和100×100 節(jié)點的自適應(yīng)網(wǎng)格分布。該問題的解在邊界x=1, y=1 附近變化比較劇烈,由圖2 可見自適應(yīng)網(wǎng)格在該邊界層比較集中。這些結(jié)果充分說明了對于二維拋物型奇異攝動問題,通過移動網(wǎng)格方法求解能自適應(yīng)地進行局部加密,能較好地體現(xiàn)解在局部區(qū)域的特性。

        圖2 t=0.2, ε2 =0.0001 時,例2 的自適應(yīng)移動網(wǎng)格分布

        來度量。圖3(a)和圖3(b)分別為50×50 空間網(wǎng)格點的自適應(yīng)移動網(wǎng)格和均勻網(wǎng)格上求解例1的誤差圖像(取t=0.5, ε=0.001 時)。為了便于比較,圖3(a)選取的網(wǎng)格節(jié)點為計算平面上均勻分布的網(wǎng)格,該網(wǎng)格是由物理平面上的不規(guī)則網(wǎng)格經(jīng)過坐標變換映射產(chǎn)生的。比較兩個圖象可以看出,移動網(wǎng)格上求解的精度優(yōu)于均勻網(wǎng)格上求解的精度。

        圖3 t=0.5, ε=0.001 時,例1 中自適應(yīng)移動網(wǎng)格與均勻網(wǎng)格上求解的逐點誤差比較

        表1 至表4 分別列出了兩個算例中t= 0.1,0.2,0.3,0.4,0.5 時刻,網(wǎng)格剖分數(shù)分別取50×50、100×100、200×200 時,在自適應(yīng)移動網(wǎng)格上獲取的數(shù)值解與真解的L2誤差和最大模誤差。小參數(shù)ε的取值如表中所示。由表中數(shù)據(jù)可知該方法對于參數(shù)ε較小的情況是有效的,具有一定的收斂性。

        表1 移動網(wǎng)格上求解例1 的L2 誤差(ε=0.001)

        表2 移動網(wǎng)格上求解例1 的最大模誤差(ε=0.001)

        表3 移動網(wǎng)格上求解例2 的L2 誤差(ε2 =0.0001)

        表4 移動網(wǎng)格上求解例2 的最大模誤差(ε2 =0.0001)

        綜上,對于含源項二維拋物型奇異攝動問題(??),使用自適應(yīng)移動網(wǎng)格方法求解既能體現(xiàn)解在局部區(qū)域的特性,同時又能保證理想的求解精度。說明該方法對于二維拋物型奇異攝動問題的求解具有有效性和優(yōu)越性。

        猜你喜歡
        拋物平面數(shù)值
        用固定數(shù)值計算
        高空拋物罪的實踐擴張與目的限縮
        法律方法(2022年2期)2022-10-20 06:45:28
        數(shù)值大小比較“招招鮮”
        關(guān)于拋物-拋物Keller-Segel類模型的全局解和漸近性
        不要高空拋物!
        高空莫拋物
        參考答案
        關(guān)于有限域上的平面映射
        基于Fluent的GTAW數(shù)值模擬
        焊接(2016年2期)2016-02-27 13:01:02
        參考答案
        国产av专区一区二区三区| 亚洲av永久无码精品放毛片| 欧美成人片一区二区三区 | 99久久国产综合精品女乱人伦| 麻豆激情视频在线观看| 男女av一区二区三区| 亚洲精品一区久久久久久| 久久中文字幕无码一区二区| 蜜臀av一区二区三区精品| 亚洲av天堂免费在线观看| 久久国产精品精品国产色婷婷| 欧美国产日本精品一区二区三区| 国产麻豆一区二区三区在线播放| 懂色av一区二区三区尤物| 熟女体下毛毛黑森林| 中文亚洲爆乳av无码专区| 92自拍视频爽啪在线观看| 2021亚洲国产精品无码| 射死你天天日| 妞干网中文字幕| 国产人妖在线观看一区二区三区| 领导边摸边吃奶边做爽在线观看 | 丝袜美足在线视频国产在线看| 精品久久人妻av中文字幕| 熟妇人妻AV中文字幕老熟妇 | 97色在线视频| 午夜国产在线精彩自拍视频| 日韩av无码社区一区二区三区| 丰满少妇高潮惨叫正在播放| 乱色视频中文字幕在线看| 二区三区日本高清视频| 成人综合网站| 99精品热6080yy久久| 亚洲天堂av高清在线| 色综合av综合无码综合网站| 84pao强力打造免费视频34| 中文少妇一区二区三区| 四虎永久在线精品免费网址| 国产高清乱理伦片| 日韩肥熟妇无码一区二区三区 | 无码人妻精品一区二区三区夜夜嗨 |