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

        ?

        求解伯格斯方程的幾種算法

        2016-04-29 00:00:00馬龍馬紅權(quán)陳書葉婷婷
        新課程·下旬 2016年7期

        摘 要:伯格斯方程(Burgers equation)是一個具有重要物理意義的數(shù)學模型。結(jié)合算例比較了基于不同徑向基函數(shù)(Matern和MQ)的局部特別解方法和Local Kansa method,分析了它們的計算誤差和優(yōu)劣。

        關(guān)鍵詞:Burgers方程;徑向基函數(shù);局部近似特別解方法

        一、引言

        對很多物理問題來說,伯格斯方程(Burgers equation)是一個非常有用的數(shù)學模型,比如激波、淺水波問題和交通流動力學問題等。而且,由于伯格斯方程是比較少的可以得到精確解的一類非線性偏微分方程,它又常被用來檢驗數(shù)值方法的好壞優(yōu)劣,這也使得伯格斯方程在計算機時代具有了重要的應(yīng)用價值。近年來,無網(wǎng)格方法求解伯格斯方程逐漸受到重視,它既不需要進行網(wǎng)格劃分,又可以有效提高計算的精度。其中,基于徑向基函數(shù)(Radial Basis Function, RBFs)的無網(wǎng)格方法具有形式簡單和各向同性等諸多優(yōu)點,并且具有較強的比較能力,在數(shù)學上得到了大量研究和成功運用。本文結(jié)合算例比較了基于不同徑向基函數(shù)(Matern和MQ)的局部近似特別解(LMAPS)方法以及Local Kansa method,在求解伯格斯方程近似解的可行性。

        二、算例

        考慮到分別用基于matern徑向基函數(shù)的LMAPS和基于MQ函數(shù)的LMAPS方法來求解方程(1),表1表示節(jié)點在單位正方形的規(guī)則計算區(qū)域上均勻分布,如圖1所示。

        三、比較分析

        分別取總節(jié)點數(shù)為121和441在t=0.4,Re=1,局部點ns=5的情況下,LMAPS分別采用Matern徑向基和MQ徑向基函數(shù), Local Kansa method方法獲得的最大絕對誤差MAE,最大相對誤差MRE和均方根誤差列表RMSE。由于Matern徑向基函數(shù)不含有形參c,所以不用像MQ函數(shù)作徑向基函數(shù)那樣去考慮形參c的取值,由下表可以看出不論是用MQ作徑向基函數(shù),還是選取Matern徑向基函數(shù),都能達到很高的近似精度,取得令人滿意的效果,但是采用Matern RBFs時獲得的各種誤差相對來講是最大的,這說明求解均勻區(qū)域點均勻分布的偏微分方程并不像求解非均勻分布的情況那樣能取得較高的近似精度。

        當節(jié)點數(shù)n=121時,采用Matern徑向基函數(shù)時LMAPS方法的最大絕對誤差,最大相對誤差和均方根誤差都達到了10-5,而采用MQ函數(shù)作徑向基函數(shù)的LMAPS方法有更高的近似精度,當總節(jié)點數(shù)都增加到n=441時,基于兩種不同徑向基函數(shù)的LMAPS方法和Local Kansa method的近似誤差都有著不同程度的提高,可見節(jié)點分布越密LMAPS方法和Local Kansa method的計算精度越高,同時采用MQ徑向基函數(shù)的LMAPS方法比用Matern徑向基函數(shù)的LMAPS方法和Local Kansa method的計算精度更高,當LMAPS方法采用Matern RBFs時近似誤差則會比Local Kansa method較大一些。

        但是隨著節(jié)點數(shù)增加時,不論是采用MQ徑向基函數(shù)LMAPS方法,還是Local Kansa method都需要隨著嘗試改變c的取值,以便取得最好的近似精度。這里需要說明的是LMAPS方法和Local Kansa method雖然都可以通過增加節(jié)點數(shù)提高計算精度,但由于會增加計算的量,故運算也需要更長的時間。

        下圖2(a)、(b)、(c)、(d)、(e)、(f)都是在正方形區(qū)域內(nèi)的點均勻分布,總點數(shù)分別取n=121和n=441,局部區(qū)域點ns=5,邊界點為ns=40和ns=80,雷諾數(shù)Re=1在t=0.4時分別采用兩種不同徑向基函數(shù)的LMAPS方法和Local Kansa method獲得的絕對誤差圖像。

        四、結(jié)論

        在求解Burgers方程時,基于全局性質(zhì)的特別解方法得到的矩陣是滿陣或者是稠密矩陣,這些矩陣往往是奇異的,如果用來解決大規(guī)模問題甚至是病態(tài)的。為了規(guī)避這些問題,人們找到了局部近似特別解方法(LMAPS)。本文采用三種方法在規(guī)則區(qū)域內(nèi)求解點均勻分布的伯格斯方程,三種方法的最大絕對誤差,最大相對誤差和均方根誤差都達到了10-5以上,證明都是有效的,誤差和計算精度都是令人滿意的。

        參考文獻:

        [1]J. M. Burger. A Mathematical Model Illustrating the Theory of Turbulence in Adv. In Appl. Mech.I, Academic Press,New York, 1948:171-199.

        [2] 程玉民.科學和工程計算的新方法:無網(wǎng)格方法[J].計算機輔助工程,2009(1).

        [3]Li S, Liu W K. Meshfree and partical methods and their

        applications.Appl. Mech.Rev.,2002,55(1):1~34.

        [4]J.C. Li, Y.C. Hon, C.S. Chen.Numerical comparisons of two meshless meth-ods using radial basis functions,Eng Anal Bound Elem,2002(26):205-225.

        在线观看亚洲av每日更新影片 | 老太脱裤让老头玩ⅹxxxx| 日韩av一区二区毛片| 亚洲国产精品久久性色av| 91蜜桃精品一区二区三区毛片| 日本视频一区二区三区三州| 亚洲麻豆视频免费观看| 国产精品免费一区二区三区四区 | 欧美精品在线一区| 蜜桃视频在线免费观看完整版| 中文字幕影片免费人妻少妇| 永久黄网站色视频免费看| 亚洲av无码国产精品麻豆天美 | 插入日本少妇一区二区三区| 色噜噜狠狠狠综合曰曰曰| 天天狠狠综合精品视频一二三区| 日韩人妻有码中文字幕| 91国内偷拍精品对白| 国产内射视频在线免费观看| 久久久久亚洲av综合波多野结衣| 女邻居的大乳中文字幕| 久久久久亚洲女同一区二区| 超清无码AV丝袜片在线观看| 五月开心六月开心婷婷网| 极品少妇xxxx精品少妇偷拍| 99久久国产福利自产拍| 国产chinese在线视频| 精品自拍偷拍一区二区三区| 免费日本一区二区三区视频 | 北岛玲中文字幕人妻系列| 一区二区午夜视频在线观看| 亚洲精品国产第一综合色吧| 亚洲av日韩专区在线观看| 亚洲国产精品午夜电影| 日本免费看一区二区三区| 亚洲精品乱码久久久久久中文字幕 | 国产乱精品女同自线免费| 国产成人亚洲综合色婷婷| 国产精品亚洲一区二区无码| 网红极品女神精品视频在线| 国产一区二区三区亚洲|