傳熱與流體流動(dòng)的數(shù)值計(jì)算

1、傳熱與流體流動(dòng)的數(shù)值計(jì)算傳熱與流體流動(dòng)的數(shù)值計(jì)算 美 S.V. 帕坦卡 著 同濟(jì)大學(xué)機(jī)械工程學(xué)院 朱 彤 本課程學(xué)習(xí)內(nèi)容本課程學(xué)習(xí)內(nèi)容 物理現(xiàn)象的數(shù)學(xué)描述 離散化方法 擴(kuò)散項(xiàng)處理 對(duì)流與擴(kuò)散 流場(chǎng)的計(jì)算 湍流數(shù)學(xué)模型 Fluent基礎(chǔ)知識(shí)介紹 參參 考考 書書 目目 傳熱與流體流動(dòng)的數(shù)值計(jì)算美 S.V. 帕坦卡 湍流是勛剛 湍流計(jì)算模型陳義良 數(shù)值傳熱學(xué)陶文銓 第一章第一章 引引 論論 1-1 研究的范疇 傳熱與流體流動(dòng)的重要性：遍及我們生活中的各個(gè) 方面 認(rèn)識(shí)和估計(jì)這些過(guò)程的必要：預(yù)報(bào)、控制 預(yù)測(cè)的本質(zhì)：說(shuō)明其中每一個(gè)物理量如何隨著幾何條件、 流量以及流體物性等的變化而改變的 目的：盡可能設(shè)

2、計(jì)一種具有完全的通用性能的數(shù)值方法 實(shí)驗(yàn)研究 全比例實(shí)驗(yàn) 模型實(shí)驗(yàn) 結(jié)果外推 測(cè)量?jī)x表精度 理論計(jì)算 理論預(yù)測(cè)出自于數(shù)學(xué)模型的結(jié)果 數(shù)學(xué)模型主要由一組微分方程組成 解析解和數(shù)值解 1-2 1-2 預(yù)測(cè)的方法預(yù)測(cè)的方法 1-2 1-2 預(yù)測(cè)的方法預(yù)測(cè)的方法 理論計(jì)算的優(yōu)點(diǎn) 成本低、速度快、資料完備、具有模擬真實(shí)條件的 能力、具有模擬理想條件的能力 理論計(jì)算的缺點(diǎn) 實(shí)際問(wèn)題分為 A類：能夠用合適的數(shù)學(xué)模型描述 B類：與A類相反的問(wèn)題湍流、多相流、NOX生成、 非牛頓流體流動(dòng) 預(yù)測(cè)方法的選擇 2-1 控制微分方程 微分方程的意義 各個(gè)微分方程都代表著一定的守恒原理。 每一個(gè)方程以一定的物理量做為它的

3、因變量， 方程本身則代表著那些影響該因變量的各個(gè)因 素之間必定存在著的某種平衡。 通常以單位質(zhì)量為基礎(chǔ)來(lái)表示各因變量。 例如：質(zhì)量分量、速度（單位質(zhì)量的動(dòng)量）、 比焓等。 第二章第二章 物理現(xiàn)象的數(shù)學(xué)描述物理現(xiàn)象的數(shù)學(xué)描述 div y xz J JJ J xyz 單位容積流出的凈流量 假設(shè)J代表一個(gè)典型因變量的流量密度 d d d x x J Jx x x J 以單位容積為基礎(chǔ)來(lái)表達(dá)一項(xiàng)變化速 率 代表在單位容積內(nèi)所包含的相應(yīng)廣延性 質(zhì)的大小 表示單位容積內(nèi)有關(guān)性質(zhì)的變化率 () t () t 化學(xué)組分的守恒化學(xué)組分的守恒 令ml代表一種化學(xué)組分l的質(zhì)量分量。當(dāng)存在速度場(chǎng)u是， 守恒表示為：

4、()div() llll mumJR t 單位容積內(nèi)化學(xué)組 分l的質(zhì)量變化率 組分l的對(duì)流 流量密度 擴(kuò)散流 量密度 單位容積化學(xué) 組分l的生成率 如果用菲克擴(kuò)散定律表示Jl， 擴(kuò)散系數(shù) 得到： lgradll Jm l ()div()div(grad) llll mummR t 能量方程能量方程 對(duì)于可以忽略粘性耗散作用的穩(wěn)態(tài)的低速流： 其中h是比焓，k是導(dǎo)熱系數(shù)，T是溫度，Sh是容積發(fā)熱率 對(duì)理想氣體以及固體和液體，將 c gradT=gradh 代入，得 到 其中c是定壓比熱。假設(shè)c為常數(shù)，即h=cT。 若u=0，則得到穩(wěn)態(tài)熱傳導(dǎo)方程： div()div(grad ) h k uhhS

5、c div()div( grad ) h uhkTS div()div(grad ) h kS uTT cc div( grad )0 h kTS 動(dòng)量方程動(dòng)量方程 由于必須同時(shí)考慮切應(yīng)力和正應(yīng)力，加之流體流動(dòng)有關(guān)的 斯托克斯粘性定律比菲克定律或傅里葉定律復(fù)雜，動(dòng)量方 程要復(fù)雜得多。用u表示x方向速度，有： 其中是粘度，p是壓力，Bx是x方向的單位容積內(nèi)體積力，Vx代表 除去以div(gradu)所表示的粘性力項(xiàng)之外的其他所有粘性力項(xiàng)。 ()div( u )div( grad ) xx p uuuBV tx 紊流的時(shí)間平均方程紊流的時(shí)間平均方程 人們假設(shè)：紊流中存在有相對(duì)平均值的快速而隨機(jī)的脈

6、動(dòng)。 由Reynold時(shí)均運(yùn)算所產(chǎn)生的附加項(xiàng)是：雷諾應(yīng)力，紊流 熱流密度，紊流擴(kuò)散流量密度等。 許多紊流模型采用紊流粘度或紊流擴(kuò)散系數(shù)的概念來(lái)表示 紊流應(yīng)力以及流量密度。結(jié)果，紊流的時(shí)間平均方程就具 有了與層流流動(dòng)方程完全相同的形式。 諸如粘度、擴(kuò)散系數(shù)以及導(dǎo)熱系數(shù)這樣一些層流交換系數(shù) 需要用相應(yīng)的有效（即層流加紊流）交換系數(shù)取代。 相當(dāng)于具有一個(gè)相當(dāng)復(fù)雜的粘度表達(dá)式的層流流動(dòng)方程。 紊流的動(dòng)能方程紊流的動(dòng)能方程 紊流“雙方程模型”：把紊流脈動(dòng)動(dòng)能k的方程作為其中 的方程之一。 K的擴(kuò) 散系數(shù) 紊流能量 的生成率 動(dòng)能的 耗散率 ()div()div(grad ) k kukkG t 通用微分

7、方程通用微分方程 其中因變量可以代表各種不同的物理量 質(zhì)量守恒或連續(xù)性方程： 擴(kuò)散系數(shù) 源項(xiàng)擴(kuò)散項(xiàng)對(duì)流項(xiàng)不穩(wěn)態(tài)項(xiàng) 直角坐標(biāo)的張量表達(dá)形式： ()div()div( grad )uS t div()0u t ()()() ()0 j jjj j j uS txxx u tx 可以代表無(wú)因次的變量 熱、質(zhì)傳遞，流體流動(dòng)，紊流以及有關(guān)的一些現(xiàn) 象的所有有關(guān)微分方程都可以看成通用方程的一 個(gè)特殊情況；可以只編寫一個(gè)求解通用方程的程 序，對(duì)不同意義的 重復(fù)使用這個(gè)程序； 對(duì)不同的 需要對(duì)相應(yīng)的和S分別賦以各自合適 的表達(dá)式，同時(shí)給出合適的初始條件和邊界條件。 自變量 一般來(lái)說(shuō)，因變量是三個(gè)空間坐標(biāo)與時(shí)間

8、的函數(shù) =（x，y，z，t） 其中x，y，z以及t都是自變量。 當(dāng)有關(guān)的物理量只與一個(gè)空間坐標(biāo)有關(guān)時(shí)，所研究 的問(wèn)題是一維的；當(dāng)問(wèn)題與時(shí)間無(wú)關(guān)時(shí)，叫做穩(wěn)態(tài) 的，否則叫做非穩(wěn)態(tài)或與時(shí)間有關(guān)的問(wèn)題。 另一種寫法： z=z（x，y，T） z是因變量，代表在位置（x，y）相對(duì)于溫度T的等溫 面高度。 2-2 2-2 坐標(biāo)的性質(zhì)坐標(biāo)的性質(zhì) 恰當(dāng)明智地選擇坐標(biāo)系統(tǒng)有時(shí)可以減少所需要的自變量數(shù)。 并非只能使用直角坐標(biāo)系，任何一種描述空間位置的方式都 是可以采用的。 例子： 1. 在一個(gè)靜止的坐標(biāo)系上看以恒定速度飛行的飛機(jī) 周圍的流體流動(dòng)是非穩(wěn)態(tài)的；但是相對(duì)于固定在飛機(jī) 上的移動(dòng)坐標(biāo)系而言，流動(dòng)是穩(wěn)態(tài)的。 2

9、. 在一圓管內(nèi)的軸對(duì)稱流動(dòng)于直角坐標(biāo)系內(nèi)是三維 的，但在r，z的圓柱極坐標(biāo)系內(nèi)則是二維的。 3. 坐標(biāo)變換可能用來(lái)進(jìn)一步減少自變量數(shù)量。 4. 改變因變量可能導(dǎo)致自變量數(shù)目的減少。 坐標(biāo)的合適選擇坐標(biāo)的合適選擇 如果在一個(gè)坐標(biāo)上的一個(gè)給定位置處的條件，要受該位置兩側(cè)條 件變化的影響的話，那么這個(gè)坐標(biāo)就是一個(gè)雙向的坐標(biāo)。 如果在一個(gè)坐標(biāo)上的一個(gè)給定位置處的條件只受該位置一側(cè)條件 變化的影響，這樣的坐標(biāo)就是一個(gè)單向的坐標(biāo)。 空間坐標(biāo)是雙向坐標(biāo)，但如果在一個(gè)坐標(biāo)方向上有很強(qiáng)的單向流 動(dòng)，也可以近似作為單向坐標(biāo)。對(duì)流單向，擴(kuò)散雙向 時(shí)間坐標(biāo)是單向坐標(biāo)。 拋物型表示一種單向的狀態(tài)；橢圓型表示雙向的概念。

10、 非穩(wěn)態(tài)導(dǎo)熱問(wèn)題實(shí)際上是時(shí)間坐標(biāo)上的拋物型和空間坐標(biāo)上的橢 圓型問(wèn)題；穩(wěn)態(tài)導(dǎo)熱對(duì)所有的坐標(biāo)都是橢圓型的。 討論單向雙向坐標(biāo)的動(dòng)機(jī)：如果可以用一個(gè)單向的坐標(biāo)來(lái)規(guī)定一個(gè)給 定的狀態(tài)，就有可能大大節(jié)省計(jì)算機(jī)的存儲(chǔ)量和時(shí)間。 單向與雙向的坐標(biāo)單向與雙向的坐標(biāo) 3-1 數(shù)值方法的本質(zhì) 一個(gè)微分方程的數(shù)值解系由一組可以構(gòu)成因 變量的分布的數(shù)所組成 類似于在實(shí)驗(yàn)室中進(jìn)行實(shí)驗(yàn) 定義 數(shù)值方法就是把計(jì)算域內(nèi)有限數(shù)量位置（叫做網(wǎng)格 結(jié)點(diǎn)）上的因變量值當(dāng)作為基本的未知量來(lái)處理。 任務(wù) 任務(wù)是提供一組關(guān)于這些未知量的代數(shù)方程并規(guī)定 求解這組方程的算法。 第三章第三章 離散化方法離散化方法 離散化的概念離散化的概念 把

11、注意力集中在網(wǎng)格結(jié)點(diǎn)處的值，用離散的值取代包含在 微分方程精確解中的連續(xù)信息。 網(wǎng)格結(jié)點(diǎn)上未知值的代數(shù)方程（離散化方程）是由支配 的微分方程推導(dǎo)而得。推導(dǎo)過(guò)程中，必須對(duì)網(wǎng)格結(jié)點(diǎn)之 間如何變化做某種假設(shè)。采用分段分布：一定的段僅僅 用一個(gè)小區(qū)域的內(nèi)部及邊界上的網(wǎng)格結(jié)點(diǎn)上的值來(lái)描述 該區(qū)間內(nèi)的變化。 一般將計(jì)算域分成一定數(shù)量的子域或單元。每個(gè)子域可以 有一個(gè)獨(dú)立的分布假設(shè)。 這種對(duì)空間和因變量所作的系統(tǒng)的離散化使得我們有可能 用比較容易求解的簡(jiǎn)單的代數(shù)方程取代控制微分方程。 離散化方程的結(jié)構(gòu)離散化方程的結(jié)構(gòu) 一個(gè)離散化方程是連接一組網(wǎng)格結(jié)點(diǎn)處值的代數(shù)關(guān)系式。 由支配的微分方程推導(dǎo)而得，表示與該微分

12、方程相同的 物理信息。 一定的離散化方程只與少數(shù)的幾個(gè)網(wǎng)格結(jié)點(diǎn)有關(guān)。在一個(gè) 網(wǎng)格結(jié)點(diǎn)處的值只影響與其緊相鄰的一些點(diǎn)上的分布。 結(jié)點(diǎn)數(shù)目變化很大時(shí)離散化方程的解趨近于相應(yīng)微分方程 的精確解。 當(dāng)網(wǎng)格結(jié)點(diǎn)緊挨在一起時(shí)，在相鄰點(diǎn)之間的變化就變得 很小，有關(guān)分布假設(shè)的實(shí)際細(xì)節(jié)就不重要了。 對(duì)于一個(gè)已知的微分方程，可能的離散化方程不是唯一的。 不同形式起因于分布假設(shè)以及推導(dǎo)方法的不同。 有限差分法和有限元法之間的區(qū)別來(lái)自選擇分布和推導(dǎo)離 散化方程的方法不同。 泰勒級(jí)數(shù)公式 如圖中網(wǎng)格結(jié)點(diǎn)。結(jié)點(diǎn) 2為結(jié)點(diǎn)1，3的中點(diǎn)， 在2周圍展開泰勒級(jí)數(shù)： 3-2 3-2 推導(dǎo)離散化方程的方法推導(dǎo)離散化方程的方法 恰好

13、在第三項(xiàng)之后截?cái)嗉?jí)數(shù)，兩方程相加相減得到： 代入微分方程就推出有限差分方程。 假設(shè)：的 變化多少 有點(diǎn)像x的 一個(gè)多項(xiàng) 式，從而 高階導(dǎo)數(shù) 項(xiàng)不那么 重要。 2 2 12 2 2 2 2 2 32 2 2 2 d1d d2d d1d d2d xx xx xx xx 31 2 2 132 22 2 d d2 d2 d() xx xx - - 變分公式變分公式 變分法證明：求解某些微分方程的問(wèn)題等效于使一稱之為 泛函的相關(guān)量最小化。即變分原理。 如果相關(guān)于因變量的網(wǎng)格點(diǎn)值使泛函最小，那么所得到的 條件即給出所需要的離散化方程。 普遍用于應(yīng)力分析的有限元法。但其適用范圍有限。 - - 加權(quán)余數(shù)法加權(quán)

14、余數(shù)法 令微分方程由 L()=0 表示。假設(shè)一個(gè)包含有若干不確定 參數(shù)的近似解，如： 其中a都是參數(shù)。上式代入微分方程留下一個(gè)余數(shù)： 假設(shè)：W是加權(quán)函數(shù)，選擇不同種類的加權(quán)函 數(shù)就可以得到不同類型的方法。 最簡(jiǎn)單的加權(quán)函數(shù)取W=1。即，對(duì)每一個(gè)控制容積，余數(shù)的積分 必須為0 0WRdx ( )RL 2 012 m m aa xa xa x 控制容積公式控制容積公式 控制容積公式可以看成是加權(quán)余數(shù)法的一種特殊形式。 把計(jì)算域分成許多互不重疊的控制容積，并使每一個(gè)網(wǎng)格 結(jié)點(diǎn)都由一個(gè)控制容積所包圍。對(duì)每一個(gè)控制容積積分微 分方程。應(yīng)用表示網(wǎng)格結(jié)點(diǎn)之間變化的分段分布關(guān)系來(lái) 計(jì)算所要求的積分。得到了一個(gè)

15、包含有一組網(wǎng)格結(jié)點(diǎn) 處的值的離散化方程。 按照這個(gè)原則所得到的離散化方程表示關(guān)于有限控制容積 的的守恒原理，這就象微分方程表示關(guān)于無(wú)窮小控制容 積內(nèi)的的守恒原理一樣。 優(yōu)點(diǎn)：所得到的結(jié)果將意味著任何一組的控制容 積內(nèi)，也就是整個(gè)計(jì)算域內(nèi)，諸如質(zhì)量、動(dòng)量以 及能量這樣一些物理量的積分守恒都可以精確地 得到滿足。 無(wú)論網(wǎng)格數(shù)量多少，這種積分守恒都滿足。 求解離散化方程以得到在網(wǎng)格結(jié)點(diǎn)上的因變量值 時(shí)，可以用兩種不同方法來(lái)表示結(jié)果。 有限元法和大多數(shù)加權(quán)余數(shù)法中把假設(shè)的由網(wǎng)格結(jié)點(diǎn) 上的值以及網(wǎng)格結(jié)點(diǎn)之間的內(nèi)插函數(shù)（或分布）所構(gòu) 成的的變化取作近似解。 有限差分法只考慮把網(wǎng)格結(jié)點(diǎn)上的值構(gòu)成解，而不去

16、管在網(wǎng)格結(jié)點(diǎn)之間是怎樣變化。類似于實(shí)驗(yàn)。 對(duì)微分方程中的不同項(xiàng)可以采用不同的分布假設(shè)來(lái)進(jìn) 行積分。 其中k是導(dǎo)熱系數(shù)，T是溫度，S是單位容積的發(fā)熱率。使用 下圖的網(wǎng)格結(jié)點(diǎn)群： 3-3 3-3 一個(gè)說(shuō)明性的例子一個(gè)說(shuō)明性的例子 一維穩(wěn)態(tài)熱傳導(dǎo)問(wèn)題方程： dd ()0 dd T kS xx dd ()()d0(3.11) dd e ew w TT kkS x xx dd d0 dd e w T kSx xx 兩種分布假設(shè)： 階梯形分布斜率 在控制容積面上是不確 定的； 線性分布網(wǎng)格結(jié) 點(diǎn)之間采用線性的內(nèi)插 函數(shù)。 分布曲線的假設(shè)分布曲線的假設(shè) Tp Tw TE 用分段線性分布計(jì)算方程（3.11）中

17、的dT/dx，所得方程： 離散化方程離散化方程 縮寫為： 比較方便的可把方程（3.13）看作： 沒有必要對(duì) 所有的量都 采用同樣的 分布函數(shù)。 () () ppEEWW e E e w W w PEW a Ta Ta Tb k a x k a x aaa bS x ppnbnb a Ta Tb ()() 0 ()() eEPwPW ew k TTkTT S x xx 上述有關(guān)選擇分布函數(shù)的自由度，最終導(dǎo)致不同變型的離 散化方程形式。當(dāng)網(wǎng)格結(jié)點(diǎn)的數(shù)目增加時(shí)，所有這些不同 形式的方程都會(huì)給出相同的解。 附加要求（可大大減少可以接受的公式的數(shù)目）： 即便是采用很粗的網(wǎng)格，解也總應(yīng)該滿足： 物理上真實(shí)

18、的性狀 一個(gè)真實(shí)的變化應(yīng)當(dāng)具有與準(zhǔn)確 變化相同的定性傾向，如圖。 總的平衡 對(duì)整個(gè)計(jì)算域應(yīng)當(dāng)滿足積分守恒。 熱流密度、質(zhì)量流量以及動(dòng)量通 量必須準(zhǔn)確地同相應(yīng)的源、匯建 立平衡，這種平衡對(duì)于任何數(shù)目 的網(wǎng)格結(jié)點(diǎn)都應(yīng)得到滿足。 物理上的真實(shí)性和總的平衡兩個(gè)約束 條件將用來(lái)指導(dǎo)選擇分布假設(shè)以及所 采用的有關(guān)措施。 指導(dǎo)原則指導(dǎo)原則 一般來(lái)說(shuō)，源項(xiàng)是因變量T本身的函數(shù) 由于離散化方程需要用線性代數(shù)的技術(shù)來(lái)求解，只能考慮一種線性的 函數(shù)關(guān)系。 源項(xiàng)的處理源項(xiàng)的處理 SC是常數(shù)部分，Sp是Tp的系數(shù)。 假設(shè)：Tp值代表整個(gè)控制容積內(nèi)的值，即采用了階梯式分布。 新的方程組： () () ppEEWW e E e w W w PEWp c a Ta Ta Tb k a x k a x aaaSx bSx c