【摘要】：多孔介質(zhì)中滲流驅(qū)動問題是一類典型和復雜的物理現(xiàn)象,這個問題與生活實際中的環(huán)境污染和油藏開采等問題是密切相關(guān)的,這也是當今科學家非常感興趣的研究領(lǐng)域之一.宏觀尺度下對這類問題的模擬一般是通過質(zhì)量守恒方程以及Darcy定律的變分形式來實現(xiàn)的,所獲得的數(shù)學模型主要是由一組非線性偏微分方程耦合而形成.其中主要包含流體的輸運和流動兩個子問題:一個是關(guān)于流動的壓力方程,另一個是關(guān)于輸運的濃度方程.Ewing在[62]中對多孔介質(zhì)的滲流驅(qū)動問題的數(shù)學模型和數(shù)值模擬做了比較全面的概述.這種滲流驅(qū)動問題的計算規(guī)模非常龐大,我們希望能夠獲得更高效的數(shù)值方法用以快速求解這種問題,提高其計算效率.本文主要探討不可壓縮流體混溶驅(qū)動問題的高效數(shù)值方法,構(gòu)造幾類不可壓縮滲流驅(qū)動問題的混合有限元方法和標準有限元方法耦合的兩層網(wǎng)格方法,討論和分析兩層網(wǎng)格算法的收斂性質(zhì).本文主要考慮混溶驅(qū)動問題中的三類不可壓縮滲流驅(qū)動問題.首先考慮只有分子擴散現(xiàn)象的驅(qū)動問題,然后研究具有重力系數(shù)項r(c)的滲流驅(qū)動問題,最后探討即有分子擴散又有分子彌散現(xiàn)象的這種復雜的不可壓縮滲流驅(qū)動問題.我們針對上述三種不可壓縮滲流驅(qū)動問題分別構(gòu)造混合有限元法和標準有限元法耦合的兩層網(wǎng)格方法,并且獲得一系列的逼近性質(zhì)和誤差估計.通過研究發(fā)現(xiàn)這種求解近似問題的兩層網(wǎng)格算法能夠保持與原問題的變分問題相同的精度,但是可以很大程度上節(jié)省計算時間.本文主要工作分為如下三個部分.第一部分,首先給出只有分子擴散現(xiàn)象的不可壓縮混溶驅(qū)動模型問題的弱形式,接著將RTk混合有限元逼近壓力和流速方程與標準有限元逼近濃度方程相結(jié)合,從而建立混合有限元離散格式.根據(jù)橢圓投影性質(zhì)和許進超教授[112]的對偶論證技巧給出混合有限元解的Lq模誤差估計.然后,在細網(wǎng)格上利用牛頓迭代思想,構(gòu)造混合有限元離散格式的快速兩層網(wǎng)格算法.它的主要思想是將細空間上一個耦合的非線性問題轉(zhuǎn)化為粗空間上一個非線性問題和細空間上一個(或者兩個)線性問題.同時,我們也獲得了濃度的Lq誤差估計和壓力與流速的收斂性質(zhì).最后,給出數(shù)值算例驗證所得到的理論結(jié)果.第二部分,我們設計具有重力系數(shù)項r(c)的不可壓縮滲流驅(qū)動問題混合有限元和標準有限元耦合的兩層網(wǎng)格算法,首先建立不可壓縮滲流驅(qū)動問題的混合有限元離散格式,通過一些逼近性質(zhì)和對偶論證方法等技巧,獲得混合有限元離散格式的Lq模估計.然后,對混合有限元離散格式引入兩層網(wǎng)格技巧,設計混合有限元法的兩層網(wǎng)格算法,進一步給出相應的收斂性估計.通過收斂性分析和數(shù)值算例,可以發(fā)現(xiàn)所構(gòu)造的兩層網(wǎng)格算法能夠保持混合有限元法的相同精度,但是大大節(jié)約了數(shù)值計算的時間.第三部分,我們主要探討具有復雜的分子擴散和彌散現(xiàn)象的不可壓縮滲流驅(qū)動問題,構(gòu)造擴散系數(shù)為非線性項的滲流驅(qū)動問題混合有限元和標準有限元耦合的兩層網(wǎng)格算法.首先,我們給出問題的弱形式;然后對壓力和流速方程利用RTk的混合有限元逼近,對濃度方程利用標準有限元逼近,我們獲得混合有限元離散格式并對有限元解進行理論分析;接著,我們設計混合有限元和標準有限元耦合的兩層網(wǎng)格方法,通過對算法的分析和探討獲得關(guān)于濃度Lq誤差估計和壓力與流速的收斂性質(zhì).
【圖文】：

邐（４．４．２）逡逑本例中，我們選取區(qū)域Ｑ邋＝邋［０，１］邋ｘ邋［０，１］和常數(shù)Ｔ邋＞邋０．時間步長Ｔ邋＝邋ｌ．Ｏｅ邋—５．逡逑為了簡單方便起見，假定空間步長／／Ｐ邋＝邐＝圪心＝／ｉｐ邋＝邋／ｉ．圖４．１和圖４．２分逡逑別給出了在時刻ｆ邋＝邋１．０ｅ＿３處兩層網(wǎng)格解＾的圖形，觀察表４．１中，我們針對算逡逑法３給出了兩層網(wǎng)格方法關(guān)于濃度，壓力和速度相應的｜｜ｃｎ－Ｃ＾｜｜邋，逡逑和｜｜ｐ＂邋－邐誤差．從誤差分析中可以觀察到當選取粗網(wǎng)格步長和細網(wǎng)格步長為逡逑把＝＆時，，關(guān)于丨丨＃邋—邋［＾丨＾和丨｜ｐｎ邋一尸廠｜丨的誤差估計為一階的收斂性，濃度逡逑｜｜ｃ＂－Ｃ＾｜｜的收斂階大約２階，這恰好和我們兩層網(wǎng)格算法理論分析中定理３．３．１和逡逑定理３．３．２的分析式相吻合的．接下來，我們來看算法４的數(shù)值算例．表４．２中獲逡逑得兩層網(wǎng)格算法４的數(shù)值解｜｜ｃｎ邋—邐，｜｜ｕ＂邋—邋Ｃ／Ｍｖ和｜丨壙—的誤差估逡逑計，從計算所得誤差的數(shù)值結(jié)果可以發(fā)現(xiàn)，當保持粗網(wǎng)格不變，分別選取合適的細逡逑網(wǎng)格步長
【學位授予單位】：湘潭大學
【學位級別】：博士
【學位授予年份】：2017
【分類號】：O241.82;O357.3

【相似文獻】

相關(guān)期刊論文 前10條

1 穆姿秀;;巧設驅(qū)動問題 促進深度學習——以“鋁的重要化合物”教學為例[J];中學化學教學參考;2016年18期

2 吳芳蘭;;多向思維,綻放火花——以《圖形中的規(guī)律》教學為例[J];湖北教育(教育教學);2016年11期

3 魯統(tǒng)超;兩相滲流驅(qū)動問題的配置方法[J];山東電力高等�？茖W校學報;1998年01期

4 袁南兒;微機控制器的總線驅(qū)動問題[J];浙江工學院學報;1989年03期

5 隋曉紅,石磊;運算放大器容性負載驅(qū)動問題研究[J];煤礦機械;2003年03期

6 程愛杰;完全可壓縮滲流驅(qū)動問題的特征線有限元方法及其誤差分析[J];高校應用數(shù)學學報A輯(中文版);1994年04期

7 龍曉瀚;劉偉;;滲流驅(qū)動問題的多步特征有限元方法[J];魯東大學學報(自然科學版);2006年04期

8 楊繼明;;不可壓滲流驅(qū)動問題混合有限元/間斷有限元耦合格式的誤差分析[J];工程數(shù)學學報;2011年03期

9 程愛杰;一類可壓縮流體驅(qū)動問題的有限元方法[J];系統(tǒng)科學與數(shù)學;1992年03期

10 袁益讓;油、水二相滲流驅(qū)動問題的變網(wǎng)格有限元方法及其理論分析[J];中國科學(A輯 數(shù)學 物理學 天文學 技術(shù)科學);1986年02期

相關(guān)會議論文 前1條

1 梁仙紅;何長江;馮其京;劉軍;李征;;自適應兩步歐拉方法及其對炸藥驅(qū)動問題的模擬[A];中國力學學會學術(shù)大會'2009論文摘要集[C];2009年

相關(guān)重要報紙文章 前8條

1 甘肅 黃浩;SONY CPD-E100驅(qū)動問題解決[N];電腦報;2002年

2 河南 趙永福;金網(wǎng)霸Modem的驅(qū)動問題[N];中國電腦教育報;2001年

3 中國社科院經(jīng)濟研究所 劇錦文;政策貴在落地生根[N];人民日報;2015年

4 中國科學院大學教授 柳卸林;創(chuàng)新驅(qū)動需要系統(tǒng)謀劃[N];經(jīng)濟日報;2013年

5 復仇公子;隨時隨地方便“打”[N];電腦報;2004年

6 胡繼妹;認識新趨勢 確立新理念 關(guān)注新重點 實現(xiàn)新目標[N];湖州日報;2015年

7 鄭杰;發(fā)現(xiàn)學校問題[N];中國教師報;2012年

8 本報記者 胡迎新;百姓嫌藥價高怎么辦？[N];健康時報;2004年

相關(guān)博士學位論文 前5條

1 劉尚;不可壓縮滲流驅(qū)動問題混合有限元與標準有限元耦合的兩層網(wǎng)格法[D];湘潭大學;2017年

2 李志濤;多孔及變形介質(zhì)中混溶及非混溶驅(qū)動問題的數(shù)值方法[D];山東大學;2007年

3 楊繼明;滲流驅(qū)動問題間斷有限元高效數(shù)值方法研究[D];湘潭大學;2007年

4 陳鳳欣;油藏數(shù)值模擬方法研究及算法設計[D];山東師范大學;2012年

5 林福泳;正交有限元及其在工程中的應用[D];天津大學;2003年

相關(guān)碩士學位論文 前9條

1 薛方云;不可壓滲流驅(qū)動問題特征混合有限元的超收斂性分析[D];湘潭大學;2005年

2 陳華;可壓縮滲流驅(qū)動問題的混合有限元和間斷Galerkin方法[D];湘潭大學;2004年

3 梁慶樂;不可壓縮滲流驅(qū)動問題的穩(wěn)定化間斷Galerkin方法[D];四川大學;2006年

4 封常榮;一維不可壓縮兩相滲流驅(qū)動問題的一類數(shù)值方法[D];煙臺大學;2008年

5 馬效培;油水兩相混溶驅(qū)動問題的一類區(qū)域分解方法[D];山東大學;2012年

6 陳瑾華;油水二相滲流驅(qū)動問題的低階非協(xié)調(diào)特征有限元方法研究[D];鄭州大學;2017年

7 席開華;兩相多組分驅(qū)動問題的數(shù)值模擬和應用[D];山東大學;2010年

8 丁永明;可壓縮可混溶驅(qū)動問題的MMOCAA格式[D];山東大學;2009年

9 王蘋;多孔介質(zhì)中混溶驅(qū)動問題的間斷有限體積元方法[D];山東師范大學;2012年

本文編號：2535616