帶Robin邊界條件的熱方程反邊值問題
發(fā)布時間:2020-07-22 07:13
【摘要】:數(shù)學(xué)物理反問題已成為具有交叉性的計算數(shù)學(xué)、應(yīng)用數(shù)學(xué)和系統(tǒng)科學(xué)中的一個重要研究方向.它在生命科學(xué)、材料科學(xué)、信號處理以及工業(yè)控制等領(lǐng)域具有重要而廣泛的應(yīng)用背景.熱傳導(dǎo)方程反邊值問題是一類重要的數(shù)學(xué)物理反問題,在熱成像、擴散光學(xué)成像和熒光成像等領(lǐng)域有廣泛的應(yīng)用.帶Dirichet邊界條件或Neumann邊界條件的熱傳導(dǎo)方程反邊值問題已有很多的研究結(jié)果,但對帶Robin邊界條件的情形研究較少.本文中利用定解區(qū)域的已知邊界r2上的測量數(shù)據(jù)重建未知邊界Γ1的形狀及Robin系數(shù)的熱傳導(dǎo)方程反問題.假設(shè)熱方程在已知的邊界r2上滿足Neumann邊界條件,而在未知的邊界Γ1上滿足Robin邊界條件,且Robin系數(shù)λ僅與空間變量有關(guān).首先,將所研究的反問題歸結(jié)為一個非線性不適定算子方程.其次,利用Newton方法將其線性化,再用正則化的最小二乘法進行求解.本文重點研究上述反問題的數(shù)值算法及其數(shù)值實現(xiàn).全文結(jié)構(gòu)如下:第一章,敘述熱傳導(dǎo)反邊值問題的研究背景和已有的工作,并介紹本文的研究工作;第二章,利用邊界積分方程方法求解正問題,將熱方程的解用單層位勢形式表示出來,再由單層位勢及其導(dǎo)數(shù)在邊界上的跳躍關(guān)系得到一個邊界積分方程組,利用配置方法對積分方程組進行離散,求出密度函數(shù),得到正問題的數(shù)值解.從而得到反問題中所需要的外邊界上的測量數(shù)據(jù);第三章,建立求解反問題的數(shù)值方法,將反邊值問題歸結(jié)為一個非線性不適定算子方程,再用Newton方法進行求解,迭代過程中需要計算算子在Frechet導(dǎo)數(shù)意義下的全微分,因此本章先給出該算子的全微分的理論推導(dǎo),再給出數(shù)值計算方法;第四章,建立反問題的數(shù)值實現(xiàn)方案,給出數(shù)值算例說明算法的可行性和有效性;第五章,總結(jié)本文的主要工作和研究方法,并對未來的相關(guān)研究問題做出展望.
【學(xué)位授予單位】:東南大學(xué)
【學(xué)位級別】:碩士
【學(xué)位授予年份】:2018
【分類號】:O175.8
【圖文】:
逡逑其中1/是邊界的單位外法向量(如圖1.1所示),A邋=邋A0r)邋>邋0是Robin系數(shù),本文中假設(shè)它逡逑僅與空間變量有關(guān),T是固定的常數(shù).逡逑@逡逑圖1.1區(qū)域示意圖逡逑該模型對應(yīng)的正問題是己知邊界6上的熱流/,區(qū)域乃,以及Robin系數(shù)A,求溫度分逡逑布而反問題是己知1^2上的熱流和相應(yīng)的溫度分布逡逑u{x,t)邋:=邋g{x,t),邐(x,t)邋6邋r2邋X邋(0,T),邐(1.2)逡逑重建未知邊界h的幾何信息及Robin系數(shù)A.具體而言,本文研宄的反問題分為三部分:逡逑1.
第四章數(shù)值實現(xiàn)例1.在花生型邊界上重建阻尼系數(shù)逡逑A邋=尸.1邐,逡逑^/1邋-0.2sin(2a)逡逑是半徑為1.5,中心在原點的圓.逡逑初始猜測取=邋1,三角多項式的展開項為尸=2,左圖是精確測量數(shù)據(jù)反演阻尼A的結(jié)果,按照上述步驟開始迭代,迭代次數(shù)為L邋=邋10,右圖是加J邋=邋0.01噪音水平的數(shù)據(jù)反演阻尼系數(shù)A的結(jié)果,迭代次數(shù)為L邋=邋13.由圖中可看到在精確測量數(shù)據(jù)下可較好的反演結(jié)果,將測量值加上誤差之后也可得到阻尼系數(shù)的大致數(shù)值.逡逑
.媝s6逡逑06*邋1邋1邋1邋'邋'邋'邋1邋邐08>邋1邋1-邋1邋1邋1—^邋1—逡逑0邋5邋10邋15邋2025303540邐0邋1邋2邋3邋4邋5邋6逡逑圖4.1精確數(shù)據(jù)邐帶1%噪音水平的誤差數(shù)據(jù)逡逑例2.在花生型邊界上重建阻尼系數(shù)為逡逑A邋=邋1邋+邋sin3邋a,逡逑徏是半徑為2,中心在原點的圓.逡逑初始猜想取A0邋=邋1,三角多項式的展開項為P邋=邋4,左圖是用精確測量數(shù)據(jù)反演阻逡逑尼系數(shù)A的結(jié)果,按照上述步驟開始迭代,迭代次數(shù)為L邋=邋10,右圖是將精確測量數(shù)據(jù)加逡逑上<5邋=邋0.01噪音水平的誤差數(shù)據(jù)反演阻尼系數(shù)A的結(jié)果,迭代次數(shù)為L邋=邋13.由下面兩圖逡逑可以看到在精確數(shù)據(jù)下可得到較為精確的數(shù)值,測量數(shù)據(jù)加上誤差之后也可得到大致的逡逑數(shù)值?逡逑
本文編號:2765493
【學(xué)位授予單位】:東南大學(xué)
【學(xué)位級別】:碩士
【學(xué)位授予年份】:2018
【分類號】:O175.8
【圖文】:
逡逑其中1/是邊界的單位外法向量(如圖1.1所示),A邋=邋A0r)邋>邋0是Robin系數(shù),本文中假設(shè)它逡逑僅與空間變量有關(guān),T是固定的常數(shù).逡逑@逡逑圖1.1區(qū)域示意圖逡逑該模型對應(yīng)的正問題是己知邊界6上的熱流/,區(qū)域乃,以及Robin系數(shù)A,求溫度分逡逑布而反問題是己知1^2上的熱流和相應(yīng)的溫度分布逡逑u{x,t)邋:=邋g{x,t),邐(x,t)邋6邋r2邋X邋(0,T),邐(1.2)逡逑重建未知邊界h的幾何信息及Robin系數(shù)A.具體而言,本文研宄的反問題分為三部分:逡逑1.
第四章數(shù)值實現(xiàn)例1.在花生型邊界上重建阻尼系數(shù)逡逑A邋=尸.1邐,逡逑^/1邋-0.2sin(2a)逡逑是半徑為1.5,中心在原點的圓.逡逑初始猜測取=邋1,三角多項式的展開項為尸=2,左圖是精確測量數(shù)據(jù)反演阻尼A的結(jié)果,按照上述步驟開始迭代,迭代次數(shù)為L邋=邋10,右圖是加J邋=邋0.01噪音水平的數(shù)據(jù)反演阻尼系數(shù)A的結(jié)果,迭代次數(shù)為L邋=邋13.由圖中可看到在精確測量數(shù)據(jù)下可較好的反演結(jié)果,將測量值加上誤差之后也可得到阻尼系數(shù)的大致數(shù)值.逡逑
.媝s6逡逑06*邋1邋1邋1邋'邋'邋'邋1邋邐08>邋1邋1-邋1邋1邋1—^邋1—逡逑0邋5邋10邋15邋2025303540邐0邋1邋2邋3邋4邋5邋6逡逑圖4.1精確數(shù)據(jù)邐帶1%噪音水平的誤差數(shù)據(jù)逡逑例2.在花生型邊界上重建阻尼系數(shù)為逡逑A邋=邋1邋+邋sin3邋a,逡逑徏是半徑為2,中心在原點的圓.逡逑初始猜想取A0邋=邋1,三角多項式的展開項為P邋=邋4,左圖是用精確測量數(shù)據(jù)反演阻逡逑尼系數(shù)A的結(jié)果,按照上述步驟開始迭代,迭代次數(shù)為L邋=邋10,右圖是將精確測量數(shù)據(jù)加逡逑上<5邋=邋0.01噪音水平的誤差數(shù)據(jù)反演阻尼系數(shù)A的結(jié)果,迭代次數(shù)為L邋=邋13.由下面兩圖逡逑可以看到在精確數(shù)據(jù)下可得到較為精確的數(shù)值,測量數(shù)據(jù)加上誤差之后也可得到大致的逡逑數(shù)值?逡逑
【參考文獻】
相關(guān)期刊論文 前1條
1 Kyoungsun Kim;Gen Nakamura;;NUMERICAL IMPLEMENTATION FOR A 2-D THERMAL INHOMOGENEITY THROUGH THE DYNAMICAL PROBE METHOD[J];Journal of Computational Mathematics;2010年01期
相關(guān)博士學(xué)位論文 前1條
1 王玉嬋;幾類拋物型方程反邊值問題的數(shù)值求解[D];東南大學(xué);2017年
本文編號:2765493
本文鏈接:http://sikaile.net/kejilunwen/yysx/2765493.html
最近更新
教材專著