• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    分離變量法和函數(shù)逼近法在偏微分方程的解法分析

    2021-05-31 03:10:02楊冬成
    保山學(xué)院學(xué)報(bào) 2021年2期
    關(guān)鍵詞:熱傳導(dǎo)常數(shù)導(dǎo)數(shù)

    楊冬成

    (江蘇聯(lián)合職業(yè)技術(shù)學(xué)院,江蘇 鹽城 224001)

    目前,在許多領(lǐng)域中數(shù)學(xué)模型都可以用偏微分方程來(lái)描述。同時(shí)許多重要物理和力學(xué)的基本方程本身就是偏微分方程,例如熱傳導(dǎo)方程[1,2]。隨著物理科學(xué)研究現(xiàn)象在廣度和深度上的擴(kuò)展,偏微分方程的應(yīng)用范圍更加廣泛[3?5]。從數(shù)學(xué)自身的角度來(lái)看,偏微分方程的求解促使數(shù)學(xué)在函數(shù)論、變分法、級(jí)數(shù)展開(kāi)、常微分方程、代數(shù)、微分幾何等各個(gè)方面進(jìn)行發(fā)展[6]。分?jǐn)?shù)階微積分是近十幾年的研究熱點(diǎn),例如一些由分?jǐn)?shù)階微分方程表示的問(wèn)題比起用經(jīng)典微分方程表示更深刻、更貼合實(shí)際[7]。分?jǐn)?shù)階微積分隱隱有取代經(jīng)典微積分之勢(shì)。而本研究主要使用不同的方法對(duì)一個(gè)分?jǐn)?shù)階偏微分方程予以求解,并對(duì)所用方法予以比較。

    1 熱傳導(dǎo)方程求解方法原理

    1.1 分離變量法原理

    分離變量法是將一個(gè)偏微分方程分解為兩個(gè)或多個(gè)常微分方程,其中該方程只含一個(gè)變量。將方程中含有各個(gè)變量的項(xiàng)分離開(kāi)來(lái),并將原方程拆分成兩個(gè)或多個(gè)更簡(jiǎn)單的只含有一個(gè)自變量的常微分方程。再運(yùn)用線性疊加原理,將非齊次方程拆分成多個(gè)齊次的或易于求解的方程。邊值問(wèn)題變量x的取值范圍總可以規(guī)范化為[0,1]。分離變量法只能對(duì)齊次邊值條件直接求解,非齊次邊值條件需要通過(guò)變換轉(zhuǎn)化為齊次邊值條件(可參考文獻(xiàn)[1])。此處以方程(1)(fi(t)=0,i=1,2,g=0)為例說(shuō)明該方法。

    令u(x,t)=X(x)T(t),即假設(shè)方程(1)的未知函數(shù)具有變量分離的特點(diǎn)。代入方程(1)中得到

    可以注意到左邊是t的函數(shù),右邊是x的函數(shù)(a是常數(shù))。因此,僅當(dāng)兩邊都為常數(shù)才能相等。故可以設(shè)該常數(shù)為λ,于是有

    但是,把零邊值條件代入,發(fā)現(xiàn)必須要求λ>0并且滿(mǎn)足:λ=n2π2,n=0,±1,…。

    記λn=n2π2,對(duì)于此λn代入(3),解出T(t)=Ane?λna2t,An是積分常數(shù)。這樣

    而線性方程具有可加性,因此,可以得到疊加之后的解:

    再把初值條件作傅里葉級(jí)數(shù)展開(kāi),求出系數(shù)An,就最終得到了方程的解。

    1.2 積分變換法原理

    積分變換是把某一類(lèi)的函數(shù)f(ξ)經(jīng)過(guò)積分F(s)=∫K(ξ,s)f(ξ)dξ變成另一類(lèi)的函數(shù)F(s)的運(yùn)算。而本研究主要使用積分變化法中的傅立葉變換(Fourier變換)和拉普拉斯變換(Laplace變換)。

    1)Fourier變換的性質(zhì):

    線性性質(zhì),設(shè)F1(ω)=F[f1(t)],F(xiàn)2(ω)=F[f2(t)],其中α和β是常數(shù),則F[αf1(t)+βf2(t)]=αF1(ω)+βF2(ω)。

    位移性質(zhì)F[f(t±t0)]=e±iωt0F[f(t)]。

    微分性質(zhì),若f(t)在(?∞,+∞)上連續(xù)或者只有有限個(gè)可去的間斷點(diǎn),并且當(dāng)|t|→ +∞時(shí),f(t)→ 0,則F[f′(t)]=iωF[f(t)]。

    2)Laplace變換的性質(zhì):

    1.3 函數(shù)逼近方法原理

    Adomian分解法原理如下,把分?jǐn)?shù)階微分方程改寫(xiě)為:

    Dαu=Lu+f(t)+N(u),n?1≤α<n

    其中L代表線性算子,N代表非線性算子,f(t)則是與u無(wú)關(guān)的項(xiàng)。兩邊作用逆算子Jα,成為,其中把初值條件代入,第一項(xiàng)是確定的。令u(t)=,代入改寫(xiě)后的方程得下式:

    把最后一項(xiàng)展開(kāi)為λ的冪級(jí)數(shù),記為。把上式中λ各冪次的系數(shù)收集起來(lái),可得出,。遞推方式就可以解出各個(gè)uk(t),從而得到方程近似解。

    2 生物傳熱微分方程的求解結(jié)果

    2.1 生物傳熱方程概述

    經(jīng)典的Pennes傳熱方程是:

    其中ρb,cb,wb分別是組織密度,組織比熱和血液灌注率,而Tb,S分別是動(dòng)脈血液溫度和環(huán)境作用的熱量,S一般可由代謝產(chǎn)熱及組織內(nèi)熱構(gòu)成[2]。生物體胳膊、腿部近似為圓柱形,軀干也可看作圓柱,因此,適用圓柱上的偏微分方程。可以使用如下方程:

    Feller[3]研究一般的擴(kuò)散問(wèn)題,把生物組織內(nèi)熱量傳導(dǎo)引起的變化由QT=kDαxu表示時(shí),就得到分?jǐn)?shù)階熱傳導(dǎo)方程:

    顯然,當(dāng)α=2就成為經(jīng)典Pennes方程。

    Povstenko[4]則利用了不同的Fourier定律,這里0≤α< 1,導(dǎo)數(shù)是Caputo導(dǎo)數(shù)。得到另一種分?jǐn)?shù)階熱傳導(dǎo)方程是:

    顯然,當(dāng)α=0也是經(jīng)典Pennes方程。但是要注意,(8)與(9)是完全不一樣的方程,它們的分?jǐn)?shù)階導(dǎo)數(shù)作用對(duì)象不同。本文考慮(9)即時(shí)間分?jǐn)?shù)階導(dǎo)數(shù)方程問(wèn)題,并假設(shè)S=0:

    借助文獻(xiàn)[5]的數(shù)據(jù),假設(shè)該方程的初邊值條件為:

    其中L是生物體中心與表皮的距離,坐標(biāo)原點(diǎn)在中心血管位置,Ta是環(huán)境溫度。作變換v=u?Tb,初邊值問(wèn)題化為:

    2.2 傅里葉-拉普拉斯變換求解方程

    2.3 分離變量法求解[6]

    2.4 函數(shù)逼近方法求解

    3 結(jié)論

    本研究以經(jīng)典的Pennes傳熱方程求解為例。采用了有限Fourier變換、Fourier?Laplace變換和分離變量法等方法求解方程。其次,運(yùn)用函數(shù)逼近方法中的Adomian分解法求出了方程的近似解。

    猜你喜歡
    熱傳導(dǎo)常數(shù)導(dǎo)數(shù)
    一類(lèi)三維逆時(shí)熱傳導(dǎo)問(wèn)題的數(shù)值求解
    解導(dǎo)數(shù)題的幾種構(gòu)造妙招
    關(guān)于Landau常數(shù)和Euler-Mascheroni常數(shù)的漸近展開(kāi)式以及Stirling級(jí)數(shù)的系數(shù)
    熱傳導(dǎo)方程解的部分Schauder估計(jì)
    一類(lèi)非線性反向熱傳導(dǎo)問(wèn)題的Fourier正則化方法
    關(guān)于導(dǎo)數(shù)解法
    幾個(gè)常數(shù)項(xiàng)級(jí)數(shù)的和
    萬(wàn)有引力常數(shù)的測(cè)量
    導(dǎo)數(shù)在圓錐曲線中的應(yīng)用
    函數(shù)與導(dǎo)數(shù)
    西和县| 即墨市| 康定县| 克什克腾旗| 文登市| 柏乡县| 敖汉旗| 江山市| 淮南市| 犍为县| 柯坪县| 阳城县| 霸州市| 朝阳市| 界首市| 福州市| 溧阳市| 北宁市| 通辽市| 尼勒克县| 贵南县| 青神县| 毕节市| 嘉鱼县| 湘潭市| 绍兴县| 克什克腾旗| 新巴尔虎左旗| 广宗县| 芒康县| 蕉岭县| 调兵山市| 临洮县| 壶关县| 神木县| 岳西县| 晋宁县| 嘉兴市| 和政县| 苏州市| 玉山县|