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

        ?

        雙層位勢(shì)問題新的基本解法公式

        2016-02-09 09:28:32張耀明
        關(guān)鍵詞:內(nèi)點(diǎn)位勢(shì)邊值問題

        郭 璇,張耀明

        (山東理工大學(xué)理學(xué)院,山東淄博255049)

        雙層位勢(shì)問題新的基本解法公式

        郭 璇,張耀明

        (山東理工大學(xué)理學(xué)院,山東淄博255049)

        基于單層位勢(shì)和疊加原理的傳統(tǒng)的基本解法,在求解某些有限域問題時(shí),虛邊界位置的選擇會(huì)受到一定的限制,在求解某些無限域問題時(shí)可能會(huì)無解。為此提出了基于雙層位勢(shì)和疊加原理的改進(jìn)的基本解法,避免了傳統(tǒng)基本解法的不足。該方法適合求解任何邊值問題,其特點(diǎn)是有限域問題和無限域問題的基本解法公式具有不同的形式。

        基本解法;雙層位勢(shì);位勢(shì)問題

        基本解法(MFS)是由Kupradze和Aleksidze[1]在20世紀(jì)60年代提出的一種著名的邊界型數(shù)值方法?;窘夥ǖ膬?yōu)點(diǎn)是不需要?jiǎng)澐志W(wǎng)格或者單元,也不需要邊界積分,同時(shí)避免了邊界元法中的奇異性和幾乎奇異性,具有程序設(shè)計(jì)簡單、計(jì)算速度快、計(jì)算精度高等優(yōu)點(diǎn),至今已被廣泛應(yīng)用于各個(gè)領(lǐng)域,取得了很好的效果[2-8]。然而,基于單層位勢(shì)和疊加原理的傳統(tǒng)的基本解法具有一定的局限性,如在求解某些有限域問題時(shí),虛邊界位置的選擇會(huì)受到一定的限制,在求解某些無限域問題時(shí)可能無解。此外,基本解法的數(shù)值解的精度與虛擬邊界的選取位置密切相關(guān),距離太近無法避免基本解的奇異性,距離太大有時(shí)可能導(dǎo)致線性系統(tǒng)矩陣高度病態(tài),不能求得問題的有效解。虛、實(shí)邊界之間的“距離選擇”一直是一個(gè)未能很好解決的問題,一般靠研究者的經(jīng)驗(yàn)和誤差實(shí)驗(yàn)進(jìn)行判斷。針對(duì)這些情況,文獻(xiàn)[9-10]分析了距離選擇與問題解的性態(tài)之間的關(guān)系,結(jié)果表明:當(dāng)問題的解可以解析開拓時(shí),虛、實(shí)邊界之間的“最大距離”不受限制。文獻(xiàn)[11]借鑒等價(jià)間接邊界積分方程的形式改進(jìn)了基于單層位勢(shì)基本解法公式,并對(duì)基本解法線性系統(tǒng)矩陣進(jìn)行了規(guī)則化,應(yīng)用于薄體結(jié)構(gòu)問題的研究,但沒有從理論上分析方法的合理性。本文從理論和數(shù)值實(shí)驗(yàn)兩方面對(duì)傳統(tǒng)基本解法的局限性進(jìn)行了分析,并在此基礎(chǔ)上提出了基于雙層位勢(shì)的基本解法。

        1 邊值問題

        本文假定Ω是R2中的一個(gè)有界域,Ωc是其開補(bǔ),Γ=?Ω是它們共同的邊界。有界域Ω上的邊值問題是[9]

        無限域Ωc上的邊值問題是

        式(1)~(2)中:▽2為拉普拉斯算子;n為邊界外法線;Γu,Γq分別是u和?u/?n已知的邊界。二維位勢(shì)問題控制方程的基本解為

        式中x=(x1,x2),y=(y1,y2)分別為場(chǎng)點(diǎn)和源點(diǎn),為兩點(diǎn)之間的距離。

        2 傳統(tǒng)的基本解法的局限性

        傳統(tǒng)基本解法的公式是

        這里l=l(y)是在點(diǎn)y處的任一方向。

        筆者發(fā)現(xiàn):傳統(tǒng)的基本解法式(4)不具有普遍適用性。對(duì)某些內(nèi)邊值問題,虛邊界(虛源位置)的選擇受到一定的限制,對(duì)于某些外邊值問題,基本解法失效。

        考慮一個(gè)內(nèi)邊值問題:

        這里Br(0)是以坐標(biāo)原點(diǎn)為中心、半徑為r<1的圓域,?Br(0)是Br(0)的邊界。顯然問題(5)的解析解是u(x)=x1+x2+1,于是u(0)=1。現(xiàn)在用基本解式(4)求解問題(5)。假設(shè)所選的虛源節(jié)點(diǎn)xj(j=1,…,n)位于單位圓周?B1上,那么對(duì)于基本解法式(4),無論{aj}取何值,總有uMFS(0)=0,因此對(duì)于任意的r<1,相應(yīng)的邊值問題(5)的解總不能由基本解式(4)求得,表明虛邊界的選擇受到一定的限制。

        數(shù)值結(jié)果也表明了這一點(diǎn)。在應(yīng)用傳統(tǒng)基本解式(4)求解問題(5)時(shí),取r=1/2,選擇44個(gè)均布的虛源點(diǎn)位于單位圓周上。表1給出了部分邊界節(jié)點(diǎn)上邊界通量?u/?n的數(shù)值結(jié)果,表2給出了內(nèi)點(diǎn)溫度的數(shù)值結(jié)果。從表1和表2可看出:無論是內(nèi)點(diǎn)計(jì)算,還是邊界點(diǎn)計(jì)算,結(jié)果都不正確,驗(yàn)證了前面的理論分析。

        表1 邊界法向通量的數(shù)值結(jié)果

        表2 內(nèi)點(diǎn)溫度的數(shù)值結(jié)果

        現(xiàn)在考慮如下外邊值問題

        例如,考慮如下邊值問題

        這里B1(0)是中心在原點(diǎn)、半徑為1的圓域,是其邊界

        在應(yīng)用傳統(tǒng)基本解式(4)求解問題(7)時(shí),44個(gè)虛源點(diǎn)均布在半徑為0.5的圓周上。表3給出了部分邊界節(jié)點(diǎn)上邊界通量?u/?n的數(shù)值結(jié)果,表4給出了內(nèi)點(diǎn)溫度的數(shù)值結(jié)果。從表3和表4可看出:內(nèi)點(diǎn)和邊界點(diǎn)的數(shù)值結(jié)果都不正確,驗(yàn)證了前面的理論分析。

        表3 內(nèi)點(diǎn)溫度的數(shù)值結(jié)果

        表4 邊界法向通量的數(shù)值結(jié)果

        3 基于雙層位勢(shì)的基本解法

        為了避免傳統(tǒng)基本解法的不足,本文提出基于雙層位勢(shì)的基本解法。對(duì)于內(nèi)邊值問題(1),基本解法公式是

        現(xiàn)在應(yīng)用基于雙層位勢(shì)的基本解法式(8)求解問題(5),仍取r=1/2。同樣地,選擇44個(gè)虛源點(diǎn)均布于單位圓周?B1(0)上。表5給出了內(nèi)點(diǎn)溫度的數(shù)值結(jié)果,表6給出了部分邊界點(diǎn)上的法向通量?u/?n的數(shù)值結(jié)果。從表5和表6可以看出:無論是內(nèi)點(diǎn)溫度的數(shù)值結(jié)果還是邊界量的數(shù)值結(jié)果都達(dá)到了很高的精度。這表明:基本解式(8)完全避免了傳統(tǒng)的基本解式(4)的弊病。

        對(duì)于外邊值問題(2),基于雙層位勢(shì)的基本解法公式是

        現(xiàn)在應(yīng)用基本解法式(8)來求解問題(3)。同樣地,選擇44個(gè)虛源點(diǎn)均布于半徑為0.5的圓周上。表7給出了內(nèi)點(diǎn)溫度的數(shù)值結(jié)果,表8給出了部分邊界點(diǎn)上的法向通量?u/?n的數(shù)值結(jié)果。從表7和表8可以看出:無論是內(nèi)點(diǎn)溫度的數(shù)值結(jié)果還是邊界量的數(shù)值結(jié)果都達(dá)到了很高的精度。這表明:基本解式(8)完全地避免了傳統(tǒng)的基本解式(4)的弊病。

        表5 內(nèi)點(diǎn)溫度的數(shù)值結(jié)果

        表6 邊界法向通量的數(shù)值結(jié)果

        表7 內(nèi)點(diǎn)溫度的數(shù)值結(jié)果

        表8 邊界法向通量的數(shù)值結(jié)果

        4 結(jié)束語

        傳統(tǒng)的基本解法基于單層位勢(shì)和疊加原理,具有一定局限性。本文從理論和數(shù)值計(jì)算兩方面進(jìn)行了分析,提出了基于雙層位勢(shì)和疊加原理的基本解法,可適合于求解任何邊值問題。

        [1] KUPRADZE V D,ALEKSIDZE M A.The method of functiona1 equations for the approximate solution of certain boundary value prob1ems[J].USSR Comput.Math.Math.Phys,1964(4):82-126.

        [2] FAIRWEATHERG,KARAGEORGHISA.The method of fundamenal solutions for elliptic boundary value problems[J].Adv Comput Math 1998(9):69-95.

        [3] POULLIKKASA,KARAGEORGHIS A,Georgiou G.Methods of fundamental solutions for harmonic and biharmonic boundary value problems[J].Comput.Mech,1998,21:416-423.

        [4] GOLBERGM A,CHEN CS.The method of fundamental solutions for potential,Helmholtz and diffusion problems.Boundary integral methods:numerical and mathematical aspects[M].Boston:WIT Press and Computational Mechanics Publications,1999.

        [5] KITAGAWA T.On the numerical stability of the method of fundamental solution applied to the Dirichlet problem[J].Jpn JAppl Math,1988(5):123-133.

        [6] GEORGE S A F,YOUSSEF F R.The method of fundamental solutions applied to3D elasticity problems using a continuous collocation scheme[J].Engineering Analysis with boundary elements,2009,33:330-341.

        [7] BOSELLIF D,OBRIST L,KLEISER.A multilayer method of fundamental solutions for Stokes flow problems[J].Journal of Computational Physics,2012,231:6139-6158.

        [8] GUIMARAES S,TELLES JC F.The method of fundamental solutions for fracture mechanics-Reissner’s plate application[J].Engineering Analysis with boundary elements,2009,33:1152-1160.

        [9] 孫煥純.無奇異邊界元法[M].大連:大連理工大學(xué)出版社,1999.

        [10]張耀明,孫煥純,楊家新.虛邊界元法的理論分析[J].計(jì)算力學(xué)學(xué)報(bào),2000,17(1):56-62.

        [11]王發(fā)杰,張耀明,公顏鵬.改進(jìn)的基本解法在薄體各向異性位勢(shì)Cauchy問題中的應(yīng)用[J].工程力學(xué),2016,32(2):18-24.

        (責(zé)任編輯陳 艷)

        A New Method of Fundamental Solution Formulation to Double Layer Potential Problems

        GUO Xuan,ZHANG Yao-ming

        (School of Mathematics,Shandong University of Technology,Zibo 255049,China)

        The conventional method of fundamental solution(MFS),based on the single layer potential and the superposition principle,has many shortcomings,for instance,for special interior problems its fictitious boundary location may be limited and for certain exterior problems it may fail.In this paper,a new method of fundamental solution formulation,which is built on the double layer potential and the superposition principle,is developed.The method avoids some disadvantages of the traditional method of fundamental solution and thus is suitable for solving any boundary value problems.

        method of fundamental solution(MFS);double layer potential;potential problems

        O241.8

        A

        1674-8425(2016)12-0165-06

        10.3969/j.issn.1674-8425(z).2016.12.026

        2016-05-27

        山東省自然科學(xué)基金重點(diǎn)資助項(xiàng)目(ZR2010AZ003)

        郭璇(1992—),女,碩士研究生,主要從事科學(xué)計(jì)算與力學(xué)研究,E-mail:m18369972890@163.com。

        郭璇,張耀明.雙層位勢(shì)問題新的基本解法公式[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué)),2016(12):165-170.< class="emphasis_bold">Citationformat:

        format:GUO Xuan,ZHANG Yao-ming.A New Method of Fundamental Solution Formulation to Double Layer Potential Problems[J].Journal of Chongqing University of Technology(Natural Science),2016(12):165-170.

        猜你喜歡
        內(nèi)點(diǎn)位勢(shì)邊值問題
        非線性n 階m 點(diǎn)邊值問題正解的存在性
        含Hardy位勢(shì)的非線性Schr?dinger-Poisson方程正規(guī)化解的多重性
        一類帶強(qiáng)制位勢(shì)的p-Laplace特征值問題
        帶有積分邊界條件的奇異攝動(dòng)邊值問題的漸近解
        含變號(hào)位勢(shì)的ρ-Kirchhoff型方程組無窮多個(gè)高能量解的存在性
        含位勢(shì)的非線性雙調(diào)和方程解的存在性
        基于罰函數(shù)內(nèi)點(diǎn)法的泄露積分型回聲狀態(tài)網(wǎng)的參數(shù)優(yōu)化
        基于內(nèi)點(diǎn)方法的DSD算法與列生成算法
        非線性m點(diǎn)邊值問題的多重正解
        一個(gè)新的求解半正定規(guī)劃問題的原始對(duì)偶內(nèi)點(diǎn)算法
        国产午夜激情视频在线看| 无遮高潮国产免费观看| 在线播放国产女同闺蜜| 国产香蕉一区二区三区在线视频| 精品九九视频| 日韩亚洲精选一区二区三区 | 久久精品国产清自在天天线| 亚洲av无码1区2区久久| 国产亚洲精品综合一区| 免费视频一区二区三区美女| 少妇真实被内射视频三四区| 亚洲国产精品成人无码区| AV教师一区高清| 日韩一区二区三区人妻中文字幕| 久久久久久人妻无码| 国产男女猛烈视频在线观看| 国产日韩精品一区二区在线观看播放 | 国产成人一区二区三区视频免费蜜 | 日韩精品无码一本二本三本色| 日韩精品无码中文字幕电影| 色诱久久av| 日本精品一区二区三区在线播放| 日本男人精品一区二区| 在办公室被c到呻吟的动态图| 天堂国精产品2023年| 中字无码av电影在线观看网站| 午夜亚洲精品一区二区| 完整版免费av片| 丰满老熟妇好大bbbbb| 99热这里只有精品69| 在线观看人成网站深夜免费| 曰日本一级二级三级人人| 99久久无码一区人妻| a级毛片内射免费视频| 免费一级a毛片在线播出| 男女上床免费视频网站| 国产又粗又猛又黄又爽无遮挡| 911精品国产91久久久久| 色噜噜精品一区二区三区| 日韩人妻免费视频一专区| 蜜桃av抽搐高潮一区二区|