若干自由邊界問(wèn)題的高效算法研究.pdf

1、自由邊界問(wèn)題,是自由邊界需要和定解問(wèn)題的解一起確定的問(wèn)題,它在科學(xué)與工程領(lǐng)域有很廣泛的應(yīng)用.由于它是非線性問(wèn)題,難以獲得解析解,因此數(shù)值模擬是研究自由邊界問(wèn)題的重要手段.本文的核心是構(gòu)造和分析反Stefan問(wèn)題和超導(dǎo)納米線單光子探測(cè)器熱電耦合模型的高效數(shù)值算法,并對(duì)所提算法進(jìn)行數(shù)值模擬驗(yàn)證其計(jì)算效率.
　　首先,我們研究了求解拋物型方程初邊值問(wèn)題的一類有限元方法的穩(wěn)定性估計(jì).利用文獻(xiàn)[1]中的連續(xù)型時(shí)空有限元方法,對(duì)于帶Diric

2、hlet邊界條件和混合邊界條件的拋物型方程初邊值問(wèn)題獲得相應(yīng)的離散化方法.使用內(nèi)蘊(yùn)的推導(dǎo)技巧,獲得了有限元解的H1(QT)-模估計(jì).這個(gè)估計(jì)在拋物型方程反問(wèn)題研究中有重要應(yīng)用,它也是我們分析后文提出的求解反Stefan問(wèn)題數(shù)值方法解的存在性和收斂性的基礎(chǔ).
　　其次,我們探討了求解反Stefan問(wèn)題的基于虛擬Robin邊界條件的區(qū)域嵌入法.本文借鑒文獻(xiàn)[2]中的區(qū)域嵌入法思想,將邊界變動(dòng)區(qū)域?(t)延拓到固定區(qū)域?,在虛擬邊界上引

3、入Robin邊界條件,用新的函數(shù)q替代未知的自由邊界函數(shù)s(t)作為控制變量,將定義在邊界變動(dòng)區(qū)域?(t)上的反Stefan問(wèn)題描述為一個(gè)線性算子方程.考慮到該算子方程的不適定性,我們通過(guò)正則化方法將其轉(zhuǎn)化為優(yōu)化問(wèn)題,在此基礎(chǔ)上使用連續(xù)型時(shí)空有限元方法實(shí)現(xiàn)最優(yōu)化問(wèn)題的離散,從而獲得求解反Stefan問(wèn)題的一個(gè)高效數(shù)值算法.證明了該方法連續(xù)情形和有限元情形問(wèn)題解的存在唯一性,并證明了有限元解收斂到連續(xù)問(wèn)題解.應(yīng)該指出的是文[2]雖然提出用

4、區(qū)域嵌入法求解一維空間變量情形的反Stefan問(wèn)題,但沒(méi)有給出任何理論分析結(jié)果.對(duì)于一維反Stefan問(wèn)題,還構(gòu)造了一個(gè)時(shí)空有限元數(shù)值解法,該方法具有執(zhí)行簡(jiǎn)便的優(yōu)點(diǎn).通過(guò)一系列的數(shù)值實(shí)驗(yàn)驗(yàn)證了這些算法的有效性.
　　最后,我們給出了超導(dǎo)納米線單光子探測(cè)器電熱耦合模型數(shù)值算法.這是一個(gè)關(guān)于溫度和電流的非線性自由邊界問(wèn)題,由非線性拋物界面方程和二階常微分方程耦合而成.在這個(gè)模型中,超導(dǎo)納米線單光子探測(cè)器近似認(rèn)為是一維結(jié)構(gòu),熱響應(yīng)是一個(gè)