流體力學根本方程組總結流體力學根本方程組包括連續性方程、運動方程、組分質量守恒方程、能量方程、本構方程、狀態方程及通用形式守恒方程。雖各相關文獻都有介紹這些根本方程組，但多采用作者熟悉的方式表示，往往同一方程具有多種形式，而難于直觀比照。以下內容是對文獻報道的各種形式的總結和比照，并分析了它們之間的轉化關系，以期徹底理解〔切實掌握微分方程中每一項的物理意義〕流體力學根本方程組的數學物理意義，為離散計算該方程組打下根底。1連續性方程根據文獻ADDINEN.CITE<EndNote><Cite><Author>吳望一</Author><Year>1982</Year><RecNum>6</RecNum><record><rec-number>6</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">吳望一</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">流體力學〔上冊〕</style></title></titles><dates><year>1982</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">北京大學出版社</style></publisher><urls></urls></record></Cite></EndNote>[1]連續性方程可由四種方法得到，分別為拉格朗日法下對有限體積和體積元應用質量守恒定律、在歐拉法下對有限體積應用質量守恒定律及在直角坐標系中直接應用質量守恒定律。L法有限體積分析取體積為，質量為的一定流體質點團，那么有：〔1〕因為速度散度的物理意義是相對體積膨脹率及密度隨體導數，即：〔2〕〔3〕代入式〔1〕得〔4〕運用奧高定理〔5〕得〔6〕上式即是連續性方程的積分形式。假定被積函數連續，而且體積是任意選取的，由此可知被積函數必須等于零，即：〔7〕或〔8〕在直角坐標系中連續性方程為：〔9〕或〔10〕連續性方程〔10〕說明，密度變化〔隨時間和位置〕等于密度和體積變形的乘積ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陳敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流體力學〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化學工業出版社</style></publisher><urls></urls></record></Cite></EndNote>[2]。L法體積元分析考慮質量為的體積元，對其用拉格朗日觀點，根據質量守恒定律有：〔11〕〔12〕兩邊同除以，得〔13〕或寫成〔14〕上式說明要維持質量守恒定律，相對體積變化率必須等于負的相對密度變化率。1.3E法有限體積分析著眼坐標空間，取空間中以面為界的有限體積，那么稱面為控制面，為控制體。取外法線方向為法線的正方向，為外法線方向的單位矢量。考慮該體積內流體質量的變化，該變化主要以下兩方面原因引起。第一，通過外表有流體流出或流入，單位時間內流出流入變化的總和為：〔15〕第二，由于密度場的不定常性〔注意，歐拉觀點下空間點是固定的，密度的變化只由場的不定常性刻畫〕，單位時間內體積的質量將變化，變化量為：〔16〕上述兩者應相等，即〔17〕由于體積是任意的，且被積函數連續，那么〔18〕1.4E法直角坐標系分析單位時間內通過外表EFGH的通量為：通過外表ABCD的通量為：其他三對外表類似，另外，該控制體內質量的變化率為：那么〔19〕特殊情況下的連續性方程：定常態：不可壓縮流體：2動量方程任取一體積為的流體，它的邊界為。根據動量定理，體積中流體動量的變化率等于作用在該體積上的質量力和面力〔應力〕之和。單位面積上的面力，其中是二階對稱應力張量，所以不是通常指的在〔單位體積面元的法線方向〕方向的分量。單位質量上的質量力為。那么作用在該體積上的質量力和面力分別為〔20〕及〔21〕動量變化率為〔22〕上述動量變化率的表達式可有兩種處理方法ADDINEN.CITE<EndNote><Cite><Author>吳望一</Author><Year>1982</Year><RecNum>6</RecNum><record><rec-number>6</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">吳望一</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">流體力學〔上冊〕</style></title></titles><dates><year>1982</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">北京大學出版社</style></publisher><urls></urls></record></Cite></EndNote>[1]，如下〔1〕求解上式右邊第二項內對體積元的隨體導數，那么〔23〕〔2〕對動量變化率表達式右邊第二項應用質量守恒定律〔24〕由上可得兩種積分形式的動量方程，即〔25〕或〔26〕由上，動量方程的微分形式為：〔27〕或〔28〕微分方程中各項的物理意義為，表示單位體積上慣性力，為單位體積上的質量力，為單位體積上應力張量的散度，它是與面力等效的體力分布函數〔由奧高公式轉化而來〕。在直角坐標系下以應力表示的運動方程可采取以下形式〔29〕或〔30〕這兩種表達方式的等號左邊實際只差了一個連續性方程，由根本微分公式〔31〕得〔32〕由連續性方程知〔33〕所以有〔34〕上述運動方程是以應力表示的粘性流體的運動方程，它們對任何粘性流體，任何運動狀態都是適用的。但它沒有反映出不同屬性的流體受力后的不同表現。另外，方程數和未知量之數不等，運動方程有三個，加上連續性方程共四個，但未知量卻有九個〔六個應力張量分量〔九個張量分量因對稱關系減少為六個〕和三個速度分量〕，所以該方程組不封閉。為使該方程組可解，必須考慮應力張量和變形速度張量之間的關系〔將應力張量用速度分量表示出來〕，補足所需的方程ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陳敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流體力學〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化學工業出版社</style></publisher><urls></urls></record></Cite></EndNote>[2]。3本構方程本構方程是表征流體宏觀性質的一種微分方程，它是表達流體粘性定律的應力張量和變形速度張量之間的關系。最簡單的應力與應變之間的關系是牛頓流體作一維運動，即牛頓剪切定律：〔35〕要得到普遍意義上的廣義牛頓定律需作一定假設，而首先應理解流體速度分解定理和變形速度張量。文獻ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陳敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流體力學〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化學工業出版社</style></publisher><urls></urls></record></Cite></EndNote>[2]對速度分解定理雖作了較直觀的描述和推導但不嚴格，而文獻ADDINEN.CITE<EndNote><Cite><Author>吳望一</Author><Year>1982</Year><RecNum>6</RecNum><record><rec-number>6</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">吳望一</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">流體力學〔上冊〕</style></title></titles><dates><year>1982</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">北京大學出版社</style></publisher><urls></urls></record></Cite></EndNote>[1]對該局部內表達較詳細。3.1速度分解定理剛體運動包括平動和轉動兩局部，一般可表為〔36〕其中是剛體中選定一點上的平動速度，是剛體繞點轉動的瞬時角速度矢量，就是要確定速度那一點到點的矢徑。轉動角速度可用表示〔37〕故〔38〕流體運動除平動、轉動外還有變形運動。設微團內點的速度為，鄰域內任一點的速度為。將在點泰勒展開并略去二階無窮小項，得〔39〕顯然，是一個二階張量〔局部速度梯度張量〕，由張量分解定理可將該張量分解成對稱張量和反對稱張量之和，于是〔40〕所以〔41〕上式右邊第二、三項可具體表示為〔42〕及〔43〕其中〔44〕另外；所以〔45〕上式說明流體運動可分為平動、轉動和變形三種形式組成，稱為變形速度張量，該定理稱為亥姆霍茲〔Helmholtz〕速度分解定理。另外，流變學中常用應變速率張量來表示流體的變形和拉伸〔或壓縮〕，而用轉動張量表示轉動，它們與流體力學中的變形速度張量和轉動張量的關系是：3.2變形速度張量的物理意義寫出變形速度的表達式〔46〕經分析可得；其中，及是角變形速率，亦稱剪切應變速率〔稱拉伸應變速率〕。文獻ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陳敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流體力學〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化學工業出版社</style></publisher><urls></urls></record></Cite></EndNote>[2]中定義的剪切速率的值是這里的一半，這是有問題的，因為剪切速率本身的值應以這里為準，但變形速度張量內剪切變形的量值為該剪切速率的一半。由上可知，變形速度張量的對角線分量，，的物理意義分別是軸線上線段元的相對拉伸速度或相對壓縮速度。而非對角線分量的物理意義分別是與軸、與軸、與軸之間夾角的剪切速率的負值。3.3廣義牛頓定律及根本假設〔1〕運動流體的應力張量在運動停止后應趨于靜止流體的應力張量。據此將應力張量寫成各向同性局部和各向異性局部是方便的。〔47〕是除去后得到的張量，稱為偏應力張量。當運動消失時它趨于零。可見，偏應力張量和應力張量一樣也是對稱張量。〔2〕偏應力張量的各分量是局部速度梯度張量各分量的線性齊次函數。當速度在空間均勻分布時，偏應力張量為零；當速度偏離均勻分布時，在粘性流體中產生了偏應力，它力圖使速度回復到均勻分布情形。〔3〕流體是各向同性的，即流體性質不依賴于方向或坐標系的轉換。根據假設〔2〕，有〔48〕顯然是一四階張量，它是表征流體粘性的常數，共個。根據假設〔3〕，是各向同性張量且對，對稱，故〔49〕觀察上式可知，對也是對稱的，物性常數減少至只有2個即第二粘度和粘度，證明見下。將上式代入偏應力表達式〔反對稱項為零〕得〔50〕那么應力張量為〔51〕引進〔52〕那么〔53〕根據上式〔54a〕〔54b〕〔54c〕將上三式等號兩邊相加，得〔55〕對可壓縮流體，流體的體積在運動過程中發生膨脹或收縮，它將引起平均法應力〔由奧高公式可證某固定點處所有方向上法應力的平均值等于三個方向上法應力的平均值，這是一個不隨坐標系改變的不變量〕的值發生的改變，稱為第二粘性系數亦稱膨脹粘性系數。應用斯托克斯假定，即，那么本構方程為〔56〕〔57〕〔58〕一般處理的是不可壓縮流體，那么〔59〕〔60〕〔61〕在直角坐標系下有〔62〕這里，有的文獻中將應力張量用表示。將上述應力張量與變形速度張量的關系式代入運動方程，得即〔63〕寫成直角坐標系下的形式〔64〕在數值傳熱學中ADDINEN.CITE<EndNote><Cite><Author>陶文銓</Author><Year>2001</Year><RecNum>8</RecNum><record><rec-number>8</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">陶文銓</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">數值傳熱學〔第二版〕</style></title></titles><dates><year>2001</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">西安</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">西安交通大學出版社</style></publisher><urls></urls></record></Cite></EndNote>[3]或CFD計算ADDINEN.CITE<EndNote><Cite><Author>王福軍</Author><Year>2004</Year><RecNum>9</RecNum><record><rec-number>9</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">王福軍</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">計算流體動力學——</style><styleface="normal"font="Arial"size="100%">CFD</style><styleface="normal"font="default"charset="134"size="100%">軟件原理與應用</style></title></titles><dates><year>2004</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">清華大學出版社</style></publisher><urls></urls></record></Cite></EndNote>[4]中常把上式等號右邊表示分子粘性作用的三項做如下變化，以第一式為例其中據此，有及〔65〕其中廣義源項定義為當流體粘度不變且不可壓縮時〔牛頓流體〕，有所以運動方程簡化為〔66〕其中是運動粘度，亦是動量擴散系數，單位。本構方程和運動方程是緊密聯系在一起的，通過本構方程可將應力張量用變形速度張量表示出來，即應力可用應變速率表示，而應變速率實際由速度分量決定，故使運動方程和連續性方程原那么上封閉可解。需要注意的是，這里討論的本構方程僅局限于牛頓流體，符合廣義牛頓定律的流體稱為牛頓流體，否那么稱為非牛頓流體。非牛頓流體的本構方程不能用廣義牛頓定律描述，如對聚合物溶液等流體應該參考相關文獻。4能量方程由能量守恒知，體積內流體的動能和內能的變換率等于單位時間內質量力和外表力所作的功加上單位時間內給予體積的熱量。體積內流體的動能和內能的總和為：〔67〕其中是單位體積內的流體內能。質量力對體積內流體所作的功為〔單位時間內移動距離，點積求做功〕：〔68〕外表力對體積內流體所作的功為〔單位時間內移動距離，點積求做功〕：〔69〕單位時間內以熱傳導方式通過外表傳給體積的熱量為：〔70〕上式被積函數實際就是傅里葉熱傳導定律，即熱流密度矢量正比于傳熱面法向溫度梯度。單位時間內由于輻射或其他原因〔反響、蒸發等〕傳入內的總熱量為〔為單位時間內傳入單位質量的熱量分布函數〕：〔71〕將上述各式進行守恒計算，得〔72〕這是積分形式的能量守恒方程。求解體積分的隨體導數并運用奧高公式把面積分轉化為體積分可得微分形式的能量守恒方程，即因為質量守恒定律所以〔73〕另外〔74〕〔75〕那么能量方程的微分形式為：〔76〕或〔77〕或〔78〕上式各項的物理意義如下，左邊第一、二項代表內能和動能的隨體導數，右邊第一項為哪一項單位體積內的質量力做功，第二項是單位體積內面力所作的功，第三項是單位體積內熱傳導輸入的熱量，最后一項表示由于輻射或其他物理或化學原因的熱量奉獻。能量守恒方程的另一種形式為〔79〕此式的物理意義為：單位體積內由于流體變形面力所作的功加上熱傳導及輻射等其他原因傳入的熱量恰好等于單位體積內的內能在單位時間內增加。將該式進一步簡化，有設〔80〕為由于粘性作用機械能轉化為熱能的局部，稱為耗散函數〔dissipationfunction〕。另外，在考慮液體流體時，比焓與內能值可看作相等，即，壓力不作功。那么所以有〔81〕其中是單位體積內熱源或由于輻射或其他物理或化學原因的熱量奉獻。一般較小可以忽略。對液體及固體可以取，進一步取為常數，并把耗散函數納入到源項，于是〔82〕對不壓縮流體〔83〕對于可以忽略粘性耗散作用的穩態低速流，能量方程可以簡化為〔84〕及〔85〕取速度為零，那么可得到穩態的熱傳導方程〔對流項消失〕：〔86〕5狀態方程由連續性方程、運動方程、能量方程確定的未知量有六個，但方程數只有五個，為使方程組封閉需補充一個聯系的狀態方程：〔87〕6組分質量守恒方程在一個特定的系統中，可能存在質的交換，或者存在多種化學組分，每一種組分都需要遵守組分質量守恒定律，即系統內某種化學組分對時間的變化率，等于通過系統界面的凈擴散流量與由反響產生的生成率之和，可表為ADDINEN.CITE<EndNote><Cite><Author>戴干策</Author><Year>2005</Year><RecNum>7</RecNum><record><rec-number>7</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"charset="134"size="100%">戴干策</style></author><author><styleface="normal"font="default"charset="134"size="100%">陳敏恒</style></author></authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">化工流體力學〔第二版〕</style></title></titles><dates><year>2005</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">化學工業出版社</style></publisher><urls></urls></record></Cite><Cite><Author>帕坦卡</Author><Year>1980</Year><RecNum>10</RecNum><record><rec-number>10</rec-number><ref-typename="Book">6</ref-type><contributors><authors><author><styleface="normal"font="default"size="100%">S</style><styleface="normal"font="default"charset="134"size="100%">.</style><styleface="normal"font="default"size="100%">V</style><styleface="normal"font="default"charset="134"size="100%">帕坦卡</style></author></authors><subsidiary-authors><author><styleface="normal"font="default"charset="134"size="100%">張政</style></author></subsidiary-authors></contributors><titles><title><styleface="normal"font="default"charset="134"size="100%">傳熱與流體流動的數值計算</style></title></titles><dates><year>1980</year></dates><pub-location><styleface="normal"font="default"charset="134"size="100%">北京</style></pub-location><publisher><styleface="normal"font="default"charset="134"size="100%">科學出版社</style></publisher><urls></urls></record></Cite></EndNote>[2,5]〔88〕其中代表單位體積內組分的質量變化率，是組分的對流流量密度。代表擴散流量密度，它由Fick定律給出。是單位體積內組分的生成率。費克定律：〔89〕其中為擴散系數。將擴散定律代入守恒方程，得〔90〕7控制方程的通用形式前面在牛頓流體的根底上，即在采用牛頓流體本構方程的根底上推導分析了運動方程和能量守恒方程，獲得了較全面的流體力學方程組，同時也采用了張量不變性記法、張量分量記法及直角坐標記法三種不同方式來表示這些根本方程組，可以說各方程之間到達了初