本文關(guān)鍵詞：守恒型分裂—區(qū)域分解格式的方法和理論及其在多孔介質(zhì)中多組分污染問題的應(yīng)用

  更多相關(guān)文章： 質(zhì)量守恒 非重疊塊狀區(qū)域分解 算子分裂 修正迎風(fēng) 無條件穩(wěn)定 誤差分析 并行性 多組分污染 多孔介質(zhì)

【摘要】：近年來,由于工業(yè)廢水,農(nóng)業(yè)灌溉和生活用水的大量排放,已經(jīng)造成地下水污染.因此,涌現(xiàn)出大量地下水污染問題的研究工作(見[11,16,31,35,52,74,76,79,80,84,102,103]等).在組分污染問題中,組分的運(yùn)移不僅僅受達(dá)西速度,分子擴(kuò)散和機(jī)械彌散的影響,還要受到組分間的化學(xué)反應(yīng)影響.一般而言,組分間的化學(xué)反應(yīng)分為動力控制溶解-沉淀反應(yīng)(見[12,19,38,39,44,49,56,63,64,70,78]等)和化學(xué)平衡反應(yīng)(見[14,45,51,65,96,97]等),例如：絡(luò)化反應(yīng),氧化-還原反應(yīng),離子交換及氣體的溶解-解析反應(yīng)等.模擬地下水中多組分運(yùn)移模型(見[18,64,71]等)需要求解耦合的非線性偏微分方程組(描述水頭的拋物型方程和描述多組分濃度的對流擴(kuò)散方程等).數(shù)值方法(有限元,有限差分等)(見[9,10,30,40,48,50,77,89,98,100,102,103]等)是一種非常有力的工具.由于偏微分方程組的耦合性和非線性,研究區(qū)域規(guī)模大和計(jì)算時(shí)間長等原因,開發(fā)求解地下水中耦合的水頭和多組分濃度方程組的質(zhì)量守恒型的區(qū)域分解的有效數(shù)值方法是非常重要的.區(qū)域分解方法(見[1,2,3,4,15,27,32,72,85,91]等)把一個(gè)大區(qū)域分成多個(gè)小子區(qū)域,然后在每個(gè)小子區(qū)域上同時(shí)求解偏微分方程組,從而有效地進(jìn)行并行計(jì)算.一般分為重疊區(qū)域分解(見[33,60,61,85,94]等)和非重疊區(qū)域分解(見[5,20,21,24,25,26,32,53,55,57,82,83,101,108,110]等).由于非重疊區(qū)域分解在每步時(shí)間間隔上具有計(jì)算量小和通信時(shí)間短的優(yōu)勢,求解拋物型方程的非迭代的顯隱非重疊區(qū)域分解方法已取得重要成果.[20,53]提出了混合格式,即首先通過顯格式計(jì)算子區(qū)域的內(nèi)邊界,然后通過隱格式計(jì)算子區(qū)域的內(nèi)解.為改善穩(wěn)定性條件,[26]進(jìn)一步提出了求解拋物型方程的顯隱區(qū)域分解(EIDD)方法.首先利用時(shí)間多步法或空間高階格式計(jì)算內(nèi)邊界,然后通過隱式計(jì)算子區(qū)域內(nèi)解.通過引入隱式校正,[83,110]得到求解二維拋物型方程的無條件穩(wěn)定顯隱區(qū)域分解算法.[83]提出了二維拋物型方程校正的顯隱區(qū)域分解算法(CEIDD),并給出了“Z”字型區(qū)域分解上的理論分析.[82]提出了求解二維拋物方程的時(shí)間三層格式的區(qū)域分解算法.即首先通過時(shí)間外推計(jì)算子區(qū)域的內(nèi)邊界點(diǎn),然后通過隱格式計(jì)算子區(qū)域內(nèi)解,最后通過隱格式校正內(nèi)邊界點(diǎn).近來,通過結(jié)合非重疊塊狀區(qū)域分解和算子分裂技巧,[57]提出了一種有效的塊狀區(qū)域分解上求解高維拋物型方程的分裂-區(qū)域分解方法(S-DDM)即通過時(shí)間多步顯式計(jì)算內(nèi)邊界點(diǎn),其次通過一維算子分裂隱格式計(jì)算子區(qū)域內(nèi)解.針對對流擴(kuò)散問題,[108]提出了求解二維變系數(shù)對流擴(kuò)散方程的時(shí)間三層的區(qū)域分解算法.[55]提出了求解帶狀區(qū)域的對流擴(kuò)散方程的區(qū)域分解算法.即通過顯式迎風(fēng)格式計(jì)算內(nèi)邊界點(diǎn),其次通過隱式修正迎風(fēng)格式計(jì)算內(nèi)邊界點(diǎn).通過結(jié)合非重疊塊狀區(qū)域分解,算子分裂和一階迎風(fēng)格式,[25]提出了一種時(shí)間多步的迎風(fēng)分裂-區(qū)域分解方法(S-DDM)求解對流擴(kuò)散方程.但是,以上的區(qū)域分解算法[20,24,25,26,53,55,57,82,83,108,110]都不滿足質(zhì)量守恒這一重要的物理定律.保持質(zhì)量守恒的區(qū)域分解算法尤其在求解空間規(guī)模大和計(jì)算時(shí)間長的問題中是非常重要的.[21]提出了求解一維和二維帶狀區(qū)域分解拋物型方程的守恒型區(qū)域分解算法.即首先通過加權(quán)算子計(jì)算內(nèi)邊界通量,然后通過隱格式計(jì)算子區(qū)域內(nèi)解.此格式是條件穩(wěn)定,且給出的二維帶狀區(qū)域分解上的誤差估計(jì)非最優(yōu).[109]提出了求解一維常系數(shù)拋物型方程的質(zhì)量守恒型區(qū)域分解算法.盡管此格式滿足質(zhì)量守恒和無條件穩(wěn)定的特點(diǎn),但是它不能求解二維塊狀區(qū)域分解的拋物型方程.格式[21,109]受到一維和二維帶狀區(qū)域分解的嚴(yán)格限制,因此,研究塊狀區(qū)域分解上的守恒型區(qū)域分解算法是一項(xiàng)非常具有挑戰(zhàn)和實(shí)際意義的工作.全文分為四章.結(jié)構(gòu)如下在第一章中,我們研究了帶有Neumann邊界條件的拋物型方程的質(zhì)量守恒型分裂-區(qū)域分解方法.拋物型方程常被用來模擬地下水的水頭,油藏中的壓力,熱傳導(dǎo)中的溫度等.由于研究和開發(fā)實(shí)際問題的規(guī)模大,計(jì)算時(shí)間長和計(jì)算存儲量大,質(zhì)量守恒的區(qū)域分解方法求解拋物型方程是非常需要的.在非重疊塊狀區(qū)域分解上,通過結(jié)合算子分裂技巧和解和通量在交錯網(wǎng)格上的耦合性技術(shù),我們提出和分析了求解二維塊狀區(qū)域分解的拋物型方程的質(zhì)量守恒型分裂-區(qū)域分解方法(S-DDM).分割計(jì)算區(qū)域?yàn)槎鄠�(gè)非重疊的塊狀子區(qū)域,在每個(gè)塊狀子區(qū)域上,我們采用兩步法計(jì)算內(nèi)解和通量.第一步,首先通過半隱交界面通量格式(顯格式)計(jì)算子區(qū)域x-方向內(nèi)邊界中間層通量,然后通過x-方向分裂解和通量的耦合隱格式計(jì)算子區(qū)域內(nèi)解和通量,再校正計(jì)算交界面的通量.第二步,通過半隱交界面通量格式(顯格式)計(jì)算子區(qū)域y-方向內(nèi)邊界中間層通量,其次通過y-方向分裂解和通量的耦合隱格式計(jì)算子區(qū)域內(nèi)解和通量,然后校正計(jì)算新層交界面的通量.所提出的質(zhì)量守恒的S-DDM格式不僅具有非重疊區(qū)域分解,算子分裂和無條件穩(wěn)定的優(yōu)勢,更重要的是格式滿足質(zhì)量守恒.由于采用算子分裂和非重疊塊狀區(qū)域分解,分析格式的穩(wěn)定性和收斂性具有挑戰(zhàn)性.通過推導(dǎo)一些關(guān)于內(nèi)邊界上數(shù)值通量的輔助性引理,定義中間層精確解和利用矩陣?yán)碚?我們證明了質(zhì)量守恒的S-DDM格式是無條件穩(wěn)定.證明了格式滿足質(zhì)量守恒,并給出格式的誤差估計(jì).數(shù)值算例驗(yàn)證理論分析.在第二章中,我們研究了二維對流擴(kuò)散方程的守恒區(qū)域分解算法.對流擴(kuò)散方程常被用來描述諸如熱量傳輸.組分質(zhì)量運(yùn)移,氣溶膠運(yùn)動等現(xiàn)象.眾所周知,在實(shí)際應(yīng)用中,由于對流項(xiàng)占優(yōu)導(dǎo)致偏微分方程雙曲特性,采用標(biāo)準(zhǔn)有限差分格式會引起數(shù)值震蕩.對于對流擴(kuò)散問題,雖已有顯隱區(qū)域分解方法的研究工作[25,55,108].但是這些區(qū)域分解算法都不滿足質(zhì)量守恒.至今未見到求解對流擴(kuò)散方程的質(zhì)量守恒區(qū)域分解方法的文獻(xiàn).我們結(jié)合二階修正迎風(fēng)技巧和解和通量的耦合格式,首次提出了求解二維塊狀區(qū)域分解的對流擴(kuò)散方程的質(zhì)量守恒的修正迎風(fēng)分裂-區(qū)域分解方法(S-DDM)在每個(gè)子塊狀區(qū)域上,我們采用兩步法計(jì)算最終的解和通量.在第一個(gè)時(shí)間半步?首先通過修正的半隱通量格式計(jì)算子區(qū)域x-方向中間層內(nèi)邊界通量,再通過x-方向分裂的解和通量的耦合修正迎風(fēng)隱格式計(jì)算子區(qū)域內(nèi)解和x-方向通量,然后校正計(jì)算中間層內(nèi)邊界通量.在第二個(gè)時(shí)間半步,通過修正的半隱通量格式計(jì)算子區(qū)域x-方向中間層內(nèi)邊界通量,其次通過y-方向分裂的解和通量的耦合修正迎風(fēng)隱格式計(jì)算子區(qū)域內(nèi)解和y-方向通量,最后校正計(jì)算新層內(nèi)邊界通量.我們證明了所提出的格式滿足質(zhì)量守恒,證明格式是無條件穩(wěn)定,進(jìn)一步分析得到誤差估計(jì).通過數(shù)值算例驗(yàn)證格式質(zhì)量守恒性,收斂性和并行性.我們提出的修正迎風(fēng)S-DDM格式不僅具有塊狀區(qū)域分解,算子分裂,二階精度和無條件穩(wěn)定.而且成功地克服了格式[25,55,108]不滿足質(zhì)量守恒的嚴(yán)重問題.在第三章中,我們分析了擴(kuò)散系數(shù)依賴于時(shí)間和空間的變系數(shù)拋物型方程的一類質(zhì)量守恒型區(qū)域分解方法.通過結(jié)合內(nèi)邊界加權(quán)平均差分格式和x-方向和y-方向的分裂技巧,我們提出和分析了帶有反應(yīng)項(xiàng)的變系數(shù)拋物型方程的非重疊塊狀區(qū)域分解上的一類質(zhì)量守恒型S-DDM格式.在每個(gè)時(shí)間間隔上,采用分?jǐn)?shù)步算子分裂格式計(jì)算塊狀子區(qū)域的內(nèi)解.首先通過多步加權(quán)平均差分格式計(jì)算x-方向內(nèi)邊界通量,然后通過x-方向分裂隱式計(jì)算子區(qū)域的中間層內(nèi)解和x-方向通量.同樣,沿y-方向計(jì)算新時(shí)刻的解和y-方向通量.我們證明所提出的內(nèi)邊界加權(quán)平均差分的分裂-區(qū)域分解格式滿足全局質(zhì)量守恒,并證明格式是弱條件穩(wěn)定.我們分析了格式的收斂性并給出離散L2范數(shù)下的最優(yōu)誤差估計(jì)O)(△t+hx2+hy2+Hx2/5Hy5/2)數(shù)值算例驗(yàn)證格式收斂性,質(zhì)量守恒性,穩(wěn)定性和并行有效性.我們提出的變系數(shù)拋物型方程的質(zhì)量守恒型S-DDM格式成功地克服了格式[21,109]僅在二維帶狀區(qū)域分解上的的限制條件.在第四章中,我們研究了質(zhì)量守恒的區(qū)域分解方法在多孔介質(zhì)中地下水多組分污染問題中的應(yīng)用.描述地下水多組分污染問題的數(shù)學(xué)模型是由一族耦合非線性偏微分方程組構(gòu)成(描述水頭的拋物型方程和描述多組分污染物濃度的對流擴(kuò)散方程).由于實(shí)際應(yīng)用中地質(zhì)結(jié)構(gòu)的復(fù)雜性,研究區(qū)域規(guī)模大和計(jì)算時(shí)間長等原因.開發(fā)一種有效求解地下水多組分污染問題的質(zhì)量守恒型的區(qū)域分解方法是非常重要的.結(jié)合非重疊塊狀區(qū)域分解,算子分裂技巧和解和通量在交錯網(wǎng)格上的耦合性,在每個(gè)時(shí)間間隔上采用兩步法依次提出求解非重疊塊狀區(qū)域分解上的水頭和污染物組分濃度的質(zhì)量守恒型S-DDM迭代算法.第一步,首先通過半隱通量格式(顯格式)計(jì)算水頭的內(nèi)交界面通量,其次通過一維算子分裂解和通量耦合的隱格式計(jì)算水頭內(nèi)解和通量.第二步,利用已得的水頭和線性化濃度法,計(jì)算x-方向和：-方向的達(dá)西速度.第三步.結(jié)合修正迎風(fēng)技巧和定義新的守恒通量,提出求解污染物組分濃度方程的質(zhì)量守恒的修正迎風(fēng)S-DDM格式.我們所提出的迭代逼近不但具有非重疊塊狀區(qū)域分解和算子分裂的優(yōu)勢,而且保持質(zhì)量守恒.數(shù)值算例驗(yàn)證格式質(zhì)量守恒,收斂階,穩(wěn)定性和并行有效性.我們模擬了在不同情形下污染物組分Ca2+,CO32-,Cl-和Na+間的沉淀反應(yīng)(CaCO3沉淀)和輸運(yùn)情況.模擬了污染物HCl流經(jīng)含有CaCO3鹽和MgCO3鹽的多孔介質(zhì)時(shí),與其發(fā)生溶解反應(yīng)產(chǎn)生新污染物Ca2+,Mg2+和弱酸H2C03的過程.實(shí)驗(yàn)表明所提出的守恒型區(qū)域分解方法很好的在塊狀區(qū)域分解上模擬了多孔介質(zhì)多組分污染問題.
【學(xué)位授予單位】：山東大學(xué)
【學(xué)位級別】：博士
【學(xué)位授予年份】：2016
【分類號】：O241.82

【相似文獻(xiàn)】

中國期刊全文數(shù)據(jù)庫 前10條

1 儲德林;求解一般橢圓偏微分方程的重迭型區(qū)域分解方法[J];計(jì)算物理;1991年04期

2 凌國燦,凌國平,王運(yùn)平;鈍體分離旋渦流動的區(qū)域分解、雜交數(shù)值模擬——Ⅰ.理論方法及其應(yīng)用[J];中國科學(xué)(A輯 數(shù)學(xué) 物理學(xué) 天文學(xué) 技術(shù)科學(xué));1992年01期

3 儲德林,胡顯承;求解橢圓型問題的一種基于區(qū)域分解的預(yù)處理器的構(gòu)造[J];高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào);1992年04期

4 馮庭桂;三溫近似輻射傳輸計(jì)算的區(qū)域分解方法[J];計(jì)算物理;1992年S1期

5 王壽城;關(guān)于并行迭代區(qū)域分解算法收斂性的注記[J];應(yīng)用數(shù)學(xué);2001年04期

6 陽述林,莫則堯,沈隆鈞;基于幾何區(qū)域分解的三維輸運(yùn)問題并行迭代算法[J];計(jì)算物理;2004年01期

7 陳星玎;曹顯兵;劉曼曼;;基于幾何非協(xié)調(diào)分解的區(qū)域分解方法誤差分析[J];數(shù)學(xué)的實(shí)踐與認(rèn)識;2011年14期

8 張輝;宋博;蔣耀林;;一種新的局部時(shí)間積分的區(qū)域分解波形松弛算法[J];中國科學(xué):數(shù)學(xué);2012年05期

9 胡顯承,儲德林;求解橢圓型問題的一種基于區(qū)域分解的預(yù)處理器[J];清華大學(xué)學(xué)報(bào)(自然科學(xué)版);1991年06期

10 胡顯承,儲德林;基于區(qū)域分解的橢圓型問題預(yù)處理器中區(qū)域邊界算子的譜分析[J];清華大學(xué)學(xué)報(bào)(自然科學(xué)版);1992年03期

中國重要會議論文全文數(shù)據(jù)庫 前7條

1 張磊;張國新;;基于重疊區(qū)域分解思想的有限元程序并行化研究[A];北京力學(xué)會第十六屆學(xué)術(shù)年會論文集[C];2010年

2 劉杰;陳靜;遲利華;胡慶豐;;結(jié)合能群和區(qū)域分解的粒子輸運(yùn)并行計(jì)算[A];全國計(jì)算物理學(xué)會第六屆年會和學(xué)術(shù)交流會論文摘要集[C];2007年

3 高偉;薛旦;張?jiān)?廖意;;基于區(qū)域分解技術(shù)的艙室攪拌器性能分析[A];2014第二屆中國指揮控制大會論文集（上）[C];2014年

4 龔瑩嵐;張宇;夏靈;;基于雙域方程求解的心臟電興奮傳導(dǎo)仿真[A];中國生物醫(yī)學(xué)工程進(jìn)展——2007中國生物醫(yī)學(xué)工程聯(lián)合學(xué)術(shù)年會論文集（下冊）[C];2007年

5 楊張衛(wèi);陳明;葉劍;陳如山;;區(qū)域分解結(jié)合多層快速多極子在電磁散射中的應(yīng)用[A];2011年全國微波毫米波會議論文集（下冊）[C];2011年

6 廉波;王正華;易曉山;;基于NACA0012翼型結(jié)構(gòu)網(wǎng)格的CFD并行模擬[A];2013年中國信息通信研究新進(jìn)展論文集[C];2014年

7 丁玉波;;并行有限元技術(shù)在復(fù)合材料厚板連接強(qiáng)度分析中的應(yīng)用[A];第九屆中國CAE工程分析技術(shù)年會專輯[C];2013年

中國重要報(bào)紙全文數(shù)據(jù)庫 前1條

1 愛芹;山東德棉實(shí)施高難品種質(zhì)量承包制[N];中國紡織報(bào);2011年

中國博士學(xué)位論文全文數(shù)據(jù)庫 前10條

1 鄭開來;基于電磁場積分方程的區(qū)域分解方法研究[D];東南大學(xué);2015年

2 鄒陽;電磁場數(shù)值計(jì)算的新型徑向基無單元算法研究[D];華中科技大學(xué);2015年

3 梁金剛;反應(yīng)堆蒙卡程序RMC大規(guī)模計(jì)算數(shù)據(jù)并行方法研究[D];清華大學(xué);2015年

4 周忠國;守恒型分裂—區(qū)域分解格式的方法和理論及其在多孔介質(zhì)中多組分污染問題的應(yīng)用[D];山東大學(xué);2016年

5 常洛;拋物方程的區(qū)域分解并行算法[D];山東大學(xué);2005年

6 段勇;偏微分方程的無網(wǎng)格區(qū)域分解方法[D];復(fù)旦大學(xué);2005年

7 田敏;時(shí)間依賴偏微分方程的區(qū)域分解并行算法[D];山東大學(xué);2007年

8 陳孝明;電磁計(jì)算問題中的區(qū)域分解無單元算法研究[D];華中科技大學(xué);2012年

9 王婷;拋物問題的顯—隱有限差分區(qū)域分解并行算法[D];山東大學(xué);2008年

10 張建松;發(fā)展方程的重疊型區(qū)域分解并行算法[D];山東大學(xué);2008年

中國碩士學(xué)位論文全文數(shù)據(jù)庫 前10條

1 李文全;非重疊型偽譜區(qū)域分解方法[D];蘇州大學(xué);2015年

2 佟佶t@;微帶貼片天線的非共形區(qū)域分解有限元分析[D];北京理工大學(xué);2016年

3 楊虹霞;粘彈性流體流動的一類特征積分平均非重疊型區(qū)域分解方法[D];長沙理工大學(xué);2013年

4 何曉明;曲邊界區(qū)域上二階偏微分方程的基于區(qū)域分解和雙二次等參有限元的分裂外推法[D];四川大學(xué);2005年

5 劉琳;偽譜區(qū)域分解方法及其應(yīng)用[D];蘇州大學(xué);2013年

6 丁玉波;區(qū)域分解技術(shù)在飛機(jī)結(jié)構(gòu)強(qiáng)度計(jì)算中的應(yīng)用研究[D];上海交通大學(xué);2011年

7 王樂楠;區(qū)域分解的有限差分算法[D];山東大學(xué);2006年

8 焦建英;求解微分方程的區(qū)域分解兩重網(wǎng)格算法[D];西安理工大學(xué);2007年

9 王麗;基于區(qū)域分解的混合有限元分裂外推[D];四川大學(xué);2006年

10 王慶超;積分微分方程的區(qū)域分裂并行算法[D];山東大學(xué);2008年

，

本文編號：1270512