【摘要】：在科學(xué)研究、工程應(yīng)用和日常生活中,流固耦合現(xiàn)象隨處可見(jiàn)。在這一過(guò)程中,流態(tài)取決于結(jié)構(gòu)體的形狀及其運(yùn)動(dòng),反過(guò)來(lái)結(jié)構(gòu)體的運(yùn)動(dòng)和變形又受制于流體施加的力�？蒲泻凸こ填I(lǐng)域急需能從數(shù)學(xué)上準(zhǔn)確闡明流固耦合過(guò)程所有細(xì)節(jié)的方法。由于流固耦合過(guò)程內(nèi)蘊(yùn)的非線性和非定常性等原因,使得采用經(jīng)典方法分析存在諸多困難。與此相反,流固耦合的數(shù)值模型在近十年間取得了飛速發(fā)展,特別的,作為當(dāng)下研究的熱點(diǎn),浸沒(méi)邊界-格子Boltzmann(IB-LB)方法被證明是研究流固耦合現(xiàn)象行之有效的途徑。浸沒(méi)邊界-格子Boltzmann法的研究已經(jīng)持續(xù)了十多載,在提出該方法之初的許多弊端都已經(jīng)有效解決,比如速度無(wú)滑移邊界無(wú)法準(zhǔn)確模擬等。該方法便于實(shí)現(xiàn),編程簡(jiǎn)單,健壯性好等優(yōu)點(diǎn)使之成為時(shí)下的研究熱點(diǎn)。然而時(shí)至今日,在將IB-LB方法應(yīng)用于工程實(shí)際時(shí),暴露出如下三個(gè)問(wèn)題:1.計(jì)算大規(guī)模復(fù)雜流固耦合問(wèn)題效率低下,計(jì)算通常持續(xù)數(shù)周之久;2.模擬高剛度彈性結(jié)構(gòu)穩(wěn)定性較差;3.模擬傳熱流固耦合問(wèn)題時(shí),在浸沒(méi)邊界處存在溫度誤差,致使邊界模擬精度下降。上述問(wèn)題限制了 IB-LB方法成功應(yīng)用于工程實(shí)際。為了發(fā)展現(xiàn)有IB-LB方法,并將改進(jìn)后方法推向水電工程應(yīng)用,本文主要圍繞如下幾個(gè)方面展開(kāi)工作:1.發(fā)展現(xiàn)有IB-LB方法,解決其穩(wěn)定性差和邊界模擬精度低的缺點(diǎn)。采用傳統(tǒng)IB-LB方法模擬溫度場(chǎng)時(shí),在浸沒(méi)邊界處會(huì)出現(xiàn)溫度誤差。為了消除溫度誤差、提高溫度邊界的模擬精度,提出一種基于迭代修正熱源的IB-LB方法。本模型基于Cheng的源項(xiàng)引入方式,熱源被視為未知量而采用迭代修正迫使溫度邊界條件精確滿足。經(jīng)算例驗(yàn)證,本方法具有一階空間收斂性,且當(dāng)LB松弛時(shí)間位于特定范圍時(shí),溫度誤差能得到有效消除。結(jié)合二維平板間的熱傳導(dǎo)算例,通過(guò)理論分析,確定了溫度誤差與LB松弛時(shí)間的關(guān)系式,揭示了溫度誤差的來(lái)源,且發(fā)現(xiàn)當(dāng)松弛時(shí)間ι1.624時(shí),本方法精度高于直接熱源法。此外,顯式或半隱式離散的IB-LB方法穩(wěn)定性較差,特別在計(jì)算高剛度彈性結(jié)構(gòu)時(shí)問(wèn)題突出。為了提高IB-LB方法的穩(wěn)定性,本文針對(duì)IB-LB的插值方程提出一種全隱離散格式。其中外力項(xiàng)采用Cheng所提方法引入,且由罰方法來(lái)處理邊界質(zhì)量。典型算例顯示,相對(duì)于傳統(tǒng)顯格式,本方法能將最大時(shí)間步長(zhǎng)提高200倍,且彈性體的最大抗彎剛度能提高40-200 倍。2.提出了一種在GPU平臺(tái)上并行IB-LB方法的策略。在格子Boltzmann(LB)方法中,對(duì)于遷移步驟所產(chǎn)生的內(nèi)存訪問(wèn)的非對(duì)齊問(wèn)題,采用一維線程塊和共享存儲(chǔ)器解決。在浸沒(méi)邊界(IB)法中,所有子過(guò)程,包括速度插值、邊界力計(jì)算和作用力擴(kuò)散,采用不同的GPU核函數(shù)實(shí)現(xiàn)。另外為了充分利用數(shù)據(jù)并行性和保證數(shù)據(jù)訪問(wèn)的合并行為,使用一維數(shù)組存儲(chǔ)浸沒(méi)邊界信息。典型算例顯示,所有GPU核函數(shù)的帶寬利用率均達(dá)到60%以上,且相對(duì)于CPU串行程序取得百倍加速比。3.在前述GPU并行基礎(chǔ)上,提出一維-三維耦合的IB-LB河庫(kù)瞬變流模擬方法,克服了傳統(tǒng)斷波法和數(shù)值解一維/二維淺水方程難以考慮涌浪橫/垂向運(yùn)動(dòng)、機(jī)組進(jìn)水口三維流動(dòng)和實(shí)際三維地形的缺點(diǎn)。其中,狹長(zhǎng)河道由一維淺水LB模型模擬,而三維流動(dòng)顯著區(qū)域則用三維自由液面LB方法模擬。實(shí)際河床地形從GIS軟件提取,且地形邊界條件由IB方法實(shí)現(xiàn)。模擬了某水利樞紐大江機(jī)組甩負(fù)荷、二江機(jī)組等流量運(yùn)行時(shí)上水庫(kù)的涌浪過(guò)程。得到了水庫(kù)水面線、各部位水位的歷時(shí)過(guò)程,確定了最大、最小涌浪高度及其衰減規(guī)律。由頻譜分析得到的水位波動(dòng)各主頻成分與由物理規(guī)律的分析結(jié)果一致,說(shuō)明本模型計(jì)算實(shí)際河床式電站瞬變流效果較好。4.為了將經(jīng)GPU加速的IB-LB方法應(yīng)用于水力機(jī)械瞬變領(lǐng)域,提出了一種動(dòng)邊界處理方法,并實(shí)現(xiàn)了貫流式水輪機(jī)甩負(fù)荷過(guò)渡過(guò)程的模擬。其中,水輪機(jī)和活動(dòng)導(dǎo)葉由基于剛體歐拉角的IB方法處理,而蝸殼和尾水管則由LB高效的曲面邊界格式實(shí)現(xiàn)。計(jì)算得到了轉(zhuǎn)輪轉(zhuǎn)速、力矩、軸向水推力和管道內(nèi)壓力的歷時(shí)曲線,并由此確定了最大轉(zhuǎn)速上升率和最大轉(zhuǎn)速上升時(shí)間等參數(shù),證明了本方法在水力機(jī)械數(shù)值仿真領(lǐng)域有較好的應(yīng)用前景,為今后計(jì)算實(shí)際工況過(guò)渡過(guò)程奠定基礎(chǔ)。本文從穩(wěn)定性、邊界模擬精度和計(jì)算效率三個(gè)方面發(fā)展現(xiàn)有IB-LB方法,突破IB-LB方法的應(yīng)用瓶頸,并將改進(jìn)后方法首次應(yīng)用于水電工程流場(chǎng)的計(jì)算。在現(xiàn)有成果基礎(chǔ)上,今后的研究重點(diǎn)可以從高雷諾數(shù)湍流模型入手,解決采用浸沒(méi)邊界法模擬湍流問(wèn)題的難題,推進(jìn)IB-LB廣泛服務(wù)于實(shí)際工程問(wèn)題。
【學(xué)位授予單位】：武漢大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2018
【分類號(hào)】：TV734
【圖文】：

Ｈｏｍｏｇｅｎｅｏｕｓ邋Ｍｏｄｅｌ邋Ｈｅｔｅｒｏｇｅｎｅｏｕｓ邋Ｍｏｄｅｌ逡逑圖１．１串行和并行程序執(zhí)行結(jié)構(gòu)的區(qū)別逡逑以往絕大多數(shù)的代碼被寫(xiě)成順序執(zhí)行的形式［２７］，也就是說(shuō)代碼按照其出現(xiàn)的２逡逑

司生產(chǎn)的ＧｅＦｏｒｃｅ邋ＧＴＸ邋２８０顯卡，這一顯卡配備了邋２４０個(gè)計(jì)算核心，其中每個(gè)核心是逡逑多線程、順序單指令執(zhí)行的。自２００３年來(lái)，眾核路線的產(chǎn)品，特別是ＧＰＵ引領(lǐng)了浮逡逑點(diǎn)運(yùn)算性能的競(jìng)賽，而圖１．２很好地說(shuō)明了這一情況。即便在多核ＣＰＵ的性能提升逡逑逐漸放緩的今天，由ＧＰＵ帶領(lǐng)的眾核產(chǎn)品仍能取得指數(shù)倍的上升。到２０１７年，市面逡逑上最先進(jìn)的ＧＰＵ和ＣＰＵ的單精度峰值浮點(diǎn)運(yùn)算性能之比是１０。如此大的性能差距逡逑來(lái)自于兩種計(jì)算核心設(shè)計(jì)理念的差別。如圖１．３所示，ＣＰＵ的設(shè)計(jì)是以最優(yōu)化循序代逡逑碼的執(zhí)行為目標(biāo)，其上復(fù)雜的邏輯控制單元能夠讓線程內(nèi)部指令按照并行或亂序的方逡逑式執(zhí)行。更重要的是，大片芯片面積被用來(lái)布置緩存和邏輯控制單元，從而能夠顯著逡逑減少指令執(zhí)行和內(nèi)存訪問(wèn)所產(chǎn)生的延遲。需要注意的是，不管是復(fù)雜的邏輯控制單元逡逑還是緩存都不能提高峰值運(yùn)算速度。反之，將絕大多數(shù)芯片面積貢獻(xiàn)給算數(shù)邏輯單元逡逑的ＧＰＵ顯然能夠支持更高的計(jì)算速度。逡逑１２０００邋－Ｉ￣，一￣■￣￣￣￣￣￣￣￣￣￣￣￣逡逑＾邐丨Ｍ中yU哎逡逑２邋１００００＇邐一ａ－邋（；ｉ

【參考文獻(xiàn)】

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

1 夏林生;程永光;蔡芳;張曉曦;;水泵水輪機(jī)四象限工作區(qū)流動(dòng)特性數(shù)值分析[J];水利學(xué)報(bào);2015年07期

2 張克鵬;;基于PowerFlow的重型載貨汽車(chē)前擾流板氣動(dòng)特性分析[J];汽車(chē)技術(shù);2014年10期

3 夏林生;程永光;張曉曦;張春澤;;燈泡式水輪機(jī)飛逸過(guò)渡過(guò)程3維CFD模擬[J];四川大學(xué)學(xué)報(bào)(工程科學(xué)版);2014年05期

4 楊桀彬;楊建東;王超;;基于空間曲面的水泵水輪機(jī)機(jī)組數(shù)學(xué)模型及仿真[J];水力發(fā)電學(xué)報(bào);2013年05期

5 ;CUDA-Zero:a framework for porting shared memory GPU applications to multi-GPUs[J];Science China(Information Sciences);2012年03期

6 楊建東;趙琨;李玲;吳鵬;;淺析俄羅斯薩揚(yáng)-舒申斯克水電站7號(hào)和9號(hào)機(jī)組事故原因[J];水力發(fā)電學(xué)報(bào);2011年04期

7 王國(guó)秉;;關(guān)于俄羅斯薩揚(yáng)·舒申斯克水電站事故的思考[J];山西水利科技;2010年02期

8 周大慶;吳玉林;劉樹(shù)紅;;軸流式水輪機(jī)模型飛逸過(guò)程三維湍流數(shù)值模擬[J];水利學(xué)報(bào);2010年02期

9 李金偉;劉樹(shù)紅;周大慶;吳玉林;;混流式水輪機(jī)飛逸暫態(tài)過(guò)程的三維非定常湍流數(shù)值模擬[J];水力發(fā)電學(xué)報(bào);2009年01期

10 劉延澤;常近時(shí);;燈泡貫流式水輪機(jī)裝置甩負(fù)荷過(guò)渡過(guò)程基于內(nèi)特性解析理論的數(shù)值計(jì)算方法[J];中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào);2008年01期

本文編號(hào)：2741677