多相流是自然界普遍存在的現(xiàn)象,對(duì)多相流的模擬是科學(xué)和工農(nóng)業(yè)應(yīng)用中常見的問題。多相流是一種復(fù)雜的流體系統(tǒng),常常伴隨著組分?jǐn)U散、相變、界面產(chǎn)生等物理和化學(xué)變化過程,之前用傳統(tǒng)的計(jì)算流體力學(xué)(CFD)方法如VOF、Level Set等對(duì)其進(jìn)行數(shù)值模擬,雖可描述少數(shù)大的界面運(yùn)動(dòng),但難以跟蹤大量細(xì)小、分散的界面。晶格Boltzmann方法(Lattice Boltzmann method,簡(jiǎn)稱LBM)是近三十年來國際上發(fā)展起來的一種流體系統(tǒng)建模和模擬的新方法,也是一種兼具了微觀分子動(dòng)力學(xué)模型和宏觀連續(xù)模型優(yōu)點(diǎn)的介觀方法,與傳統(tǒng)的CFD方法相比有明顯的優(yōu)勢(shì),如:有清晰的物理背景、算法簡(jiǎn)單、計(jì)算效率高、有高度的并行性等。經(jīng)過近30年來的發(fā)展與改進(jìn),LBM已成為計(jì)算流體力學(xué)中不可替代的新方法,在研究多相流領(lǐng)域中占有重要地位,尤其是在包含界面動(dòng)力學(xué)、相變等問題的復(fù)雜系統(tǒng)中進(jìn)行建模。然而,由于多相流問題本身的復(fù)雜性,進(jìn)一步運(yùn)用LBM于實(shí)際中還存在許多挑戰(zhàn),如:系統(tǒng)平衡時(shí)虛速度大、密度比較小等。近年來,人們?cè)谌绾卧黾酉到y(tǒng)模擬的穩(wěn)定性和密度比,降低平衡時(shí)的虛速度,以及提高熱力學(xué)一致性和伽利略不變性等方面做了大量的工作。一些研究者從力項(xiàng)技術(shù)的角度進(jìn)行了研究,認(rèn)為宏觀力在離散化時(shí)出現(xiàn)了偏差,并從恢復(fù)宏觀方程的角度進(jìn)行了分析,從而提出了許多的力項(xiàng)技術(shù),其中以精確差分方法為優(yōu)。還有不少研究者從非理想相互作用力的形式方面進(jìn)行了深入研究,目前計(jì)算相互作用力的方法主要有兩種,一個(gè)是以Shan-Chen為始的偽勢(shì)模型,另一個(gè)是從熱力學(xué)自由能角度出發(fā)的非理想力模型。自由能非理想力模型具有堅(jiān)實(shí)的熱力學(xué)理論基礎(chǔ),由該作用力模擬的兩相共存密度曲線雖與Maxwell等面積的理論值吻合較好,但仍未獲得理想的結(jié)果。在偽勢(shì)模型中,研究者雖然不斷對(duì)其進(jìn)行改進(jìn)與完善,但數(shù)值穩(wěn)定性、虛速度等問題仍然未能很好地解決。近年來有研究者提出組合勢(shì)的偽勢(shì)模型,因其精度高而廣為流傳與運(yùn)用,然而,本文在進(jìn)一步探究后發(fā)現(xiàn)該模型存在些問題,仍有深入研究和改進(jìn)的必要。首先,我們分析幾乎所有的非理想相互作用力模型,都涉及到梯度的計(jì)算,而密度梯度的計(jì)算在非理想相互作用計(jì)算中尤為重要,如果計(jì)算不準(zhǔn),將會(huì)對(duì)多相流模擬結(jié)果產(chǎn)生重大影響。目前為止,關(guān)于非理想作用力中計(jì)算梯度的算法屬傳統(tǒng)的中心差分算法最流行,因而少有人深入考究傳統(tǒng)的中心差分算法是否精確,但數(shù)值微分本身是典型的不適定問題,計(jì)算中稍有不慎就可能會(huì)產(chǎn)生很大誤差,甚至?xí)a(chǎn)生數(shù)值不穩(wěn)定問題。接著,我們分析了傳統(tǒng)中心差分算法在多相流模擬中的影響。典型的多相流平衡密度分布近似地呈雙曲正切函數(shù)分布,通過分析數(shù)值導(dǎo)數(shù)的精度,發(fā)現(xiàn)中心差分法(CDM)計(jì)算的數(shù)值導(dǎo)數(shù)與其理論解相比存在較大的偏差。再者,既然中心差分與理論解析存在偏差,那么是否存在更好的數(shù)值計(jì)算方法?我們聯(lián)想到在有限差分中常被用于數(shù)值計(jì)算的高階差分,經(jīng)分析,高階差分法(HDM)計(jì)算的數(shù)值導(dǎo)數(shù)與理論解析可以很好地吻合。因此,我們引入HDM來計(jì)算粒子間的相互作用,替代了常用的CDM,并提出了基于高階差分的偽勢(shì)多相流晶格Boltzmann方法。本模型不僅保留了偽勢(shì)模型的優(yōu)點(diǎn),實(shí)現(xiàn)簡(jiǎn)單,效率高,并行性好等特點(diǎn),而且具有較高的精度。為了驗(yàn)證該模型的性能,進(jìn)行了多個(gè)方面多相流模擬,包括靜態(tài)和動(dòng)態(tài)兩種情況。經(jīng)理論和數(shù)值試驗(yàn),本文提出的基于高階差分的偽勢(shì)多相流晶格Boltzmann方法主要具有如下優(yōu)勢(shì)和特色:1.在大溫度范圍和大密度比的情況下,可以實(shí)現(xiàn)常用的狀態(tài)方程的熱力學(xué)一致性,且沒有任何組合相互作用和任何額外的可調(diào)參數(shù)。2.在文中分析后得出,Shan-Chen的有效密度與Zhang的勢(shì)函數(shù)在數(shù)學(xué)理論分析上是一致的,可以同等地描述所有選擇的狀態(tài)方程的流體相變,因此,相互作用勢(shì)可任選其一。3.前人提出了很多的力項(xiàng)技術(shù),其中Kupershtokh提出的精確差分因其簡(jiǎn)單、準(zhǔn)確且與弛豫時(shí)間無關(guān)等優(yōu)點(diǎn)受到廣大研究者的青睞,故本文采取的力項(xiàng)技術(shù)就是精確差分。4.創(chuàng)新引入新的計(jì)算相互作用的數(shù)值方法。本文通過一系列的分析發(fā)現(xiàn),傳統(tǒng)的中心差分并不是模擬多相流計(jì)算梯度時(shí)最好的數(shù)值方法。因此,引入了一種數(shù)值計(jì)算方法——高階差分,該方法無論從數(shù)學(xué)理論上還是雙曲正切函數(shù)的模擬,結(jié)果都比傳統(tǒng)的中心差分好�？梢�,運(yùn)用高階差分來模擬多相流要比傳統(tǒng)的中心差分應(yīng)該好很多。5.將新模型在兩相共存與真實(shí)多相流中進(jìn)行模擬,并初步應(yīng)用液滴飛濺的實(shí)例充分說明該模型更可行、實(shí)用。數(shù)值實(shí)驗(yàn)中可發(fā)現(xiàn),運(yùn)用高階差分模擬的結(jié)果很好。6.新模型穩(wěn)定準(zhǔn)確,它的最大密度比甚至達(dá)到了10~9,足以滿足幾乎所有實(shí)際應(yīng)用的要求,故用高階差分的方法應(yīng)用于多相流的研究具有足夠的優(yōu)勢(shì),并期望獲得一些更有趣的結(jié)果。通過上述分析,不僅說明新模型擁有堅(jiān)實(shí)的理論基礎(chǔ)和良好的數(shù)值性能,而且體現(xiàn)了高階差分提出的必要性,以期在多相流數(shù)值計(jì)算領(lǐng)域中得到廣泛地應(yīng)用。
【學(xué)位單位】：廣西師范大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位年份】：2018
【中圖分類】：O359
【文章目錄】：
摘要
ABSTRACT
第1章 緒論
    1.1 研究背景及意義
    1.2 晶格Boltzmann方法對(duì)多相流研究的現(xiàn)狀及其發(fā)展趨勢(shì)
    1.3 論文的主要工作及創(chuàng)新點(diǎn)
    1.4 論文的組織與結(jié)構(gòu)
第2章 偽勢(shì)晶格Boltzmann方法及多相流模型基本理論
    2.1 晶格Boltzmann方法基本要素
        2.1.1 晶格Boltzmann方法基本模型
        2.1.2 單馳豫晶格Boltzmann方法（LBGK）的多尺度展開
    2.2 邊界條件
        2.2.1 反彈邊界法
        2.2.2 周期邊界法
    2.3 單組份偽勢(shì)多相流晶格Boltzmann方法
        2.3.1 Shan-Chen單組份偽勢(shì)模型
        2.3.2 Zhang等人的作用力偽勢(shì)模型
        2.3.3 Yuan等人的大密度比偽勢(shì)模型
        2.3.4 Kupershtokh等人改進(jìn)的偽勢(shì)模型
        2.3.5 Maxwell等面積法則
    2.4 單組份偽勢(shì)晶格Boltzmann模型的改進(jìn)
        2.4.1 幾種常用的狀態(tài)方程
        2.4.2 使用常規(guī)非理想氣體狀態(tài)方程
        2.4.3 調(diào)節(jié)狀態(tài)方程
    2.5 小結(jié)
第3章 中心差分與高階差分的偽勢(shì)晶格Boltzmann方法的基本分析
    3.1 引言
    3.2 傳統(tǒng)中心差分?jǐn)?shù)值方法
    3.3 傳統(tǒng)的中心差分算法在多相流模擬中可能存在不可忽略的誤差
        3.3.1 van der Waals流體的基本理論值
        3.3.2 傳統(tǒng)中心差分在van der Waals流體平衡態(tài)下的應(yīng)用
    3.4 高階差分?jǐn)?shù)值方法
        3.4.1 高階差分的基本理論
        3.4.2 高階差分在van der Waals流體平衡態(tài)下的應(yīng)用
    3.5 高階差分與中心差分在雙曲正切函數(shù)模擬中的驗(yàn)證對(duì)比
    3.6 小結(jié)
第4章 基于高階差分的偽勢(shì)多相流晶格Boltzmann方法
    4.1 引言
    4.2 確定相互作用力形式
    4.3 CDM和HDM在實(shí)際多相流中的性能
    4.4 模擬兩相共存密度曲線
    4.5 小結(jié)
第5章 基于高階差分模型的檢驗(yàn)與初步應(yīng)用
    5.1 概述
    5.2 模型的驗(yàn)證
        5.2.1 穩(wěn)定性與效率驗(yàn)證
        5.2.2 表面張力與Laplace定律驗(yàn)證
        5.2.3 虛速度
    5.3 模型的初步應(yīng)用
    5.4 小結(jié)
第6章 總結(jié)與展望
    6.1 全文總結(jié)
    6.2 研究展望
參考文獻(xiàn)
攻讀碩士學(xué)位期間的科研成果和其他獲獎(jiǎng)情況
致謝

【相似文獻(xiàn)】

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

1 湯波;李俊峰;王天舒;;帶自由面流體運(yùn)動(dòng)的單相格子Boltzmann方法模擬[J];清華大學(xué)學(xué)報(bào)(自然科學(xué)版)網(wǎng)絡(luò).預(yù)覽;2008年11期

2 靳遵龍;王元?jiǎng)P;王永慶;;基于格子Boltzmann方法的固定方柱繞流研究[J];鄭州大學(xué)學(xué)報(bào)(理學(xué)版);2018年03期

3 劉強(qiáng);謝偉;張?jiān)俜?田斌斌;邱遼原;;熵格子Boltzmann方法的重整化群代數(shù)湍流模型[J];船舶力學(xué);2017年11期

4 田忠偉;趙志剛;周亞;趙春波;;基于格子Boltzmann方法的裂隙介質(zhì)滲流模擬研究[J];煤炭技術(shù);2017年11期

5 黃兵方;聞炳海;邱文;趙琬玲;陳燕雁;;基于晶格Boltzmann方法的接觸角實(shí)時(shí)測(cè)量研究[J];廣西師范大學(xué)學(xué)報(bào)(自然科學(xué)版);2018年01期

6 霍宇濤;饒中浩;;基于偽焓法的多孔介質(zhì)固液相變格子Boltzmann方法[J];工程熱物理學(xué)報(bào);2018年05期

相關(guān)博士學(xué)位論文 前2條

1 錢尼貴;基于格子Boltzmann方法的路面路基細(xì)觀滲流特性及試驗(yàn)研究[D];華南理工大學(xué);2018年

2 劉冰冰;格子Boltzmann方法氣泡泵提升管內(nèi)氣泡運(yùn)動(dòng)行為研究[D];大連海事大學(xué);2017年

相關(guān)碩士學(xué)位論文 前10條

1 趙琬玲;基于高階差分的偽勢(shì)多相流晶格Boltzmann方法的研究[D];廣西師范大學(xué);2018年

2 耿孟孟;基于晶格Boltzmann方法的液滴彈跳現(xiàn)象的研究[D];廣西師范大學(xué);2018年

3 朱洪宇;基于漲落-格子Boltzmann方法的微顆粒布朗運(yùn)動(dòng)和泳動(dòng)的模擬研究[D];西安建筑科技大學(xué);2018年

4 周旋;基于化學(xué)勢(shì)的多相流晶格Boltzmann方法的研究[D];廣西師范大學(xué);2017年

5 焦晗;基于改進(jìn)的浸入邊界—格子Boltzmann方法圓柱繞流研究[D];哈爾濱工程大學(xué);2016年

6 周喻偉;格子Boltzmann方法對(duì)圓形空間自然對(duì)流的數(shù)值模擬[D];上海理工大學(xué);2014年

7 蔡銀娟;格子Boltzmann方法在含植被河道水流數(shù)值模擬中的應(yīng)用研究[D];北方民族大學(xué);2017年

8 韓永賓;含間斷系數(shù)擴(kuò)散方程的格子Boltzmann方法[D];湘潭大學(xué);2017年

9 李昕;格子Boltzmann方法研究微裂縫形態(tài)特征對(duì)滲流的影響[D];電子科技大學(xué);2017年

10 李紅群;基于格子Boltzmann方法的圓柱渦激振動(dòng)及傳熱特性研究[D];哈爾濱工程大學(xué);2016年

本文編號(hào)：2888228