趙麗娜，徐國(guó)賓

（天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津　300072）

基于廣義流和廣義力的河流能耗率推導(dǎo)

趙麗娜，徐國(guó)賓

（天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津300072）

基于非平衡態(tài)熱力學(xué)理論中局域熵產(chǎn)生和能耗函數(shù)的關(guān)系，以及能耗函數(shù)可以寫作不可逆過程中每個(gè)廣義流與其相應(yīng)的廣義力的乘積之和，構(gòu)造了河流的廣義力和廣義流，推導(dǎo)出河流的能耗函數(shù)數(shù)學(xué)表達(dá)式.通過對(duì)能耗函數(shù)做體積分，進(jìn)而推導(dǎo)出河流的能耗率.從而避開使用眾多復(fù)雜的非平衡態(tài)熱力學(xué)和高等流體力學(xué)公式，過程簡(jiǎn)便且易于理解.

河流；非平衡態(tài)熱力學(xué)；廣義流；廣義力；能耗率

河流是個(gè)開放的、復(fù)雜的動(dòng)力系統(tǒng)，不僅具有自動(dòng)調(diào)整功能趨向于動(dòng)力平衡，而且河流內(nèi)部能量趨向于按照一定規(guī)律分配.國(guó)內(nèi)外許多學(xué)者曾從能耗率的角度去研究河流水流的運(yùn)動(dòng)規(guī)律，其中有美籍華裔楊志達(dá)和國(guó)內(nèi)徐國(guó)賓等學(xué)者.Yang等［1-5］利用河流的勢(shì)能與熱力學(xué)系統(tǒng)的熱能之間的類同，認(rèn)為河流中唯一有用的能量是它的勢(shì)能，并進(jìn)一步假定河流中的勢(shì)能和高程分別相當(dāng)于熱力學(xué)系統(tǒng)中的勢(shì)能和絕對(duì)溫度，根據(jù)這些概念，提出了單位水流功率，并基于單位水流功率概念，給出了最小能耗率原理數(shù)學(xué)表達(dá)式.徐國(guó)賓等［6-9］基于非平衡態(tài)熱力學(xué)的最小熵產(chǎn)生原理，利用流體力學(xué)中的連續(xù)方程、運(yùn)動(dòng)方程和能量方程以及熱力學(xué)的吉布斯公式，經(jīng)過嚴(yán)謹(jǐn)?shù)臄?shù)學(xué)推演得到了流體最小能耗率原理數(shù)學(xué)表達(dá)式，進(jìn)一步完善了最小能耗率理論，同時(shí)指出最小能耗率原理等價(jià)于最小熵產(chǎn)生原理.

本文將基于非平衡態(tài)熱力學(xué)理論，構(gòu)造河流的廣義流和廣義力，推導(dǎo)出河流的能耗函數(shù)數(shù)學(xué)表達(dá)式，進(jìn)而推導(dǎo)出河流的能耗率，從而避開使用眾多復(fù)雜的非平衡態(tài)熱力學(xué)和高等流體力學(xué)公式，使推導(dǎo)過程簡(jiǎn)便且易于理解.

1　熵流和熵產(chǎn)生的回顧

當(dāng)系統(tǒng)存在有外力場(chǎng)且其內(nèi)部存在黏滯性流動(dòng)時(shí)，單位體積中的熵平衡方程［10］為

式中：SV為熵密度，表示單位體積中的熵；t為時(shí)間；u為平均流速；jq為熱流；T為溫度；kμ為單位質(zhì)量的第k組分的化學(xué)勢(shì)；jk為單位質(zhì)量的第k組分的物質(zhì)流；Mk為單位質(zhì)量的第k組分的分子量；Fk為單位質(zhì)量的第k組分的質(zhì)量力；vk，λ為參與第λ個(gè)反應(yīng)的反應(yīng)物k的化學(xué)計(jì)量系數(shù).

一個(gè)系統(tǒng)的總熵隨時(shí)間的變化可以寫作

式中：S為系統(tǒng)的總熵；V為系統(tǒng)的體積；Ω為系統(tǒng)的表面積；JS代表通過單位面積的熵的交換速率，簡(jiǎn)稱熵流；σ代表單位體積中熵產(chǎn)生的速率，即局域熵產(chǎn)生.

利用Gauss定律，由式（2）可得

比較式（1）和式（3）得到熵流和熵產(chǎn)生的表達(dá)式

其中，定義

式（4）和式（5）中的每一項(xiàng)都有明確的物理意義.式（4）右邊第1項(xiàng)代表由對(duì)流過程引起的熵流，第2項(xiàng)為由熱傳導(dǎo)引起的熵流，第3項(xiàng)為由擴(kuò)散過程引起的熵流.

現(xiàn)在來專門研究一下式（5）.式（5）右邊的第1項(xiàng)與熱傳導(dǎo)有關(guān)，第2項(xiàng)和擴(kuò)散過程有關(guān)，第3項(xiàng)和黏滯性流動(dòng)有關(guān)，第4項(xiàng)和化學(xué)反應(yīng)有關(guān).值得注意的是式中每一項(xiàng)本身都是由兩個(gè)因子的乘積組成.其中的一個(gè)因子和不可逆過程的速率有關(guān)，例如熱流jq、擴(kuò)散流jk、黏性應(yīng)力張量Π/T和化學(xué)反應(yīng)速率ωλ，這些速率因子可以稱為不可逆過程的廣義流，簡(jiǎn)稱“流”；式（5）右邊各項(xiàng)中的另一個(gè)因子和引起相應(yīng)“流”的推動(dòng)力有關(guān)，例如溫度梯度?（1/T）引起熱流，化學(xué)勢(shì)梯度-?（μk/T）以及外力MkFk/T引起擴(kuò)散流，流速梯度-?u引起黏滯性流動(dòng)，Aλ/T是化學(xué)反應(yīng)的推動(dòng)力.這些和推動(dòng)力有關(guān)的因子可以廣義地稱之為不可逆過程的廣義力，簡(jiǎn)稱“力”.如果用Jk代表第k種不可逆過程的流，用Xk代表第k種不可逆過程的力，則式（5）可以寫作如下一般形式

這就是說，局域熵產(chǎn)生可以寫作不可逆過程的廣義流和相應(yīng)廣義力的乘積之和的形式.

2　基于廣義流和廣義力的河流能耗率推導(dǎo)

根據(jù)非平衡態(tài)熱力學(xué)，局域熵產(chǎn)生和能耗函數(shù)存在下列關(guān)系［10］

式中φ為能耗函數(shù).

在實(shí)際應(yīng)用中，可以認(rèn)為局域熵產(chǎn)生等價(jià)于能耗函數(shù)，則根據(jù)式（7），能耗函數(shù)也可以表示為

需要注意的是，構(gòu)造能耗函數(shù)的廣義流和廣義力的原則和構(gòu)造熵產(chǎn)生的廣義流和廣義力的原則一樣，只不過構(gòu)造能耗函數(shù)的廣義流和廣義力時(shí)，需要兩者的乘積具有能耗函數(shù)的量綱［ML-1，T-3］，其中M、L及T分別代表質(zhì)量、長(zhǎng)度及時(shí)間.

一旦合理地確定了廣義力和廣義流之后，便可以由式（9）定量地計(jì)算由不可逆過程引起的能耗函數(shù).

2.1河流的廣義力和廣義流的選擇

河流中存在兩種流，即物質(zhì)流和能量流［11］.在不考慮熱交換的情況下，能量流是通過動(dòng)量傳輸實(shí)現(xiàn)的，物質(zhì)流是通過擴(kuò)散完成的.

眾所周知，單位體積的流體所具有的動(dòng)量為uρ，將通過給定流體空間邊界面的流體動(dòng)量通量定義為動(dòng)量流，則動(dòng)量流等于流體的動(dòng)量uρ和質(zhì)量中心速度u的乘積［11］，即

式中ρ是水流密度.

流速梯度的存在將發(fā)生動(dòng)量擴(kuò)散，這類擴(kuò)散由高速區(qū)指向低速區(qū)，于是動(dòng)量流對(duì)應(yīng)的廣義力是沿著流向的梯度［10］，即

式中l(wèi)為流向的坐標(biāo)軸.

對(duì)于質(zhì)量流，定義為密度乘以速度，即單位時(shí)間內(nèi)擴(kuò)散的質(zhì)量［11］

在河流中對(duì)應(yīng)于水流流動(dòng)的驅(qū)動(dòng)因子是重力.重力在l方向的分量為質(zhì)量流對(duì)應(yīng)的廣義力［11］，即

式中：g為重力加速度；i表示水面比降［12］，，對(duì)于平原河流，i的值很小時(shí)（i≤6°）sini=i［12］，于是得到

2.2量綱分析

在國(guó)際單位制中，對(duì)于力學(xué)問題，采用的3個(gè)基本量綱分別是長(zhǎng)度、時(shí)間和質(zhì)量，即L、T、M制，廣義力和廣義流乘積的量綱都可以由長(zhǎng)度、時(shí)間和質(zhì)量的量綱導(dǎo)出［12］.河流的動(dòng)量流和動(dòng)量力的乘積為，其量綱為河流的質(zhì)量流和質(zhì)量力的乘積為JmXm=ρu· gi，其量綱為，可以看出本文構(gòu)造的河流的廣義力和廣義流的乘積具有能耗函數(shù)的量綱.

2.3河流的能耗函數(shù)

根據(jù)式（9）可知，河流的能耗函數(shù)可以表示為

河流的水力比降［12］，則式（16）變?yōu)?/p>

式中：γ為水容重；J為水力比降.

對(duì)能耗函數(shù)求體積分，得到流體的能耗率的表達(dá)式為

沿水流方向取單位長(zhǎng)度l=1，并且設(shè)在單位長(zhǎng)度內(nèi)J是常數(shù)，則式（18）簡(jiǎn)化為

式中：A為過流斷面；Q為流量.

設(shè)式（19）中，過水?dāng)嗝鏋槌?shù)，所以式（19）可以簡(jiǎn)化為Φ=uJ.為了實(shí)用起見，對(duì)于大多數(shù)的天然河流，可以用水面比降i取代水力比降J［13］，因此楊志達(dá)定義的單位水流功率ui是式（19）的簡(jiǎn)化和實(shí)用形式.

3　結(jié) 語

本文基于非平衡態(tài)熱力學(xué)理論，能耗函數(shù)等于每個(gè)廣義流與其相應(yīng)的廣義力的乘積之和，構(gòu)造了河流的廣義力和廣義流，進(jìn)而推導(dǎo)得出河流能耗率表達(dá)式，該表達(dá)式與楊志達(dá)和徐國(guó)賓推導(dǎo)得出的能耗率形式一樣，但推導(dǎo)過程簡(jiǎn)便且易于理解，從另一個(gè)角度再次證明了河流的能耗率公式的確是Φ=γQJ .

［1］Yang C T. Potential energy and stream morphology ［J］. Water Resources Research，1971，7（2）：311-322.

［2］Yang C T. Unit stream power and sediment transport ［J］. Journal of the Hydraulics Division，ASCE，1972，98：1805-1826.

［3］Yang C T. Minimum unit stream power and fluvial hydraulics ［J］. Journal of the Hydraulics Division，ASCE，1976，102（7）：919-934.

［4］Yang C T，Song C C S. Theory of minimum rate of energy dissipation ［J］. Journal of the Hydraulics Division，ASCE，1979，105：769-784.

［5］Yang C T，Song C C S. Theory of Minimum Energy andEnergy Dissipation Rate ［M］. Houston：Gulf Publishing Company，1986.

［6］徐國(guó)賓，楊志達(dá). 基于最小熵產(chǎn)生與耗散結(jié)構(gòu)和混沌理論的河床演變分析［J］. 水利學(xué)報(bào)，2012，43（8）：948-956.

Xu Guobin，Yang C T. Analysis of river bed changes based on the theories of minimum entropy production dissipative structure and chaos ［J］. Journal of Hydraulic Engineering，2012，43（8）：948-956（in Chinese）.

［7］徐國(guó)賓，練繼建. 流體最小熵產(chǎn)生原理與最小能耗率原理（Ⅰ）［J］. 水利學(xué)報(bào)，2003（5）：35-40.

Xu Guobin，Lian Jijian. Theories of the minimum rate of energy dissipation and the minimum entropy production of flow（Ⅰ）［J］. Journal of Hydraulic Engineering，2003（5）：35-40（in Chinese）.

［8］徐國(guó)賓，練繼建. 流體最小熵產(chǎn)生原理與最小能耗率原理（Ⅱ）［J］. 水利學(xué)報(bào)，2003（6）：43-47.

Xu Guobin，Lian Jijian. Theories of the minimum rate of energy dissipation and the minimum entropy production of flow（Ⅱ）［J］. Journal of Hydraulic Engineering，2003（6）：43-47（in Chinese）.

［9］徐國(guó)賓，練繼建. 河流調(diào)整中的熵、熵產(chǎn)生和能耗率的變化［J］. 水科學(xué)進(jìn)展，2004，15（1）：1-5.

Xu Guobin，Lian Jijian. Changes of the entropy production and the rate of energy dissipation in river adjustment［J］. Advances in Water Science，2004，15（1）：1-5（in Chinese）.

［10］李如生. 非平衡態(tài)熱力學(xué)和耗散結(jié)構(gòu)［M］. 北京：清華大學(xué)出版社，1986.

Li Rusheng. Non-Equilibrium Thermodynamics and Dissipative Structures［M］. Beijing：Tsinghua University Press，1986（in Chinese）.

［11］王 錚. 利用耗散結(jié)構(gòu)理論分析河道演變［J］. 地理科學(xué)，1989，9（2）：173-180.

Wang Zheng. Analysis of river channel change with theory of dissipative structure ［J］. Scientia Geographica Sinica，1989，9（2）：173-180（in Chinese）.

［12］李家星，趙振興. 水力學(xué)［M］. 南京：河海大學(xué)出版社，2005. Li Jiaxing，Zhao Zhenxing. Hydraulics ［M］. Nanjing：Hohai University Press，2005（in Chinese）.

［13］Yang C T. Formation of riffles and pools ［J］. Water Resources Research，1971，7（6）：1567-1574.

（責(zé)任編輯：樊素英）

Derivation of Energy Dissipation Rate of River Based on Generalized Flux and Generalized Force

Zhao Lina，Xu Guobin
（State Key Laboratory of Hydraulic Engineering Simulation and Safety，Tianjin University，Tianjin 300072，China）

Based on non-equilibrium thermodynamics，local entropy production is a function of energy dissipation function，and the latter can be written as the sum of the product of the generalized forces and the corresponding generalized flux of an irreversible process.In this way，once the generalized forces and the corresponding generalized flux of river are formed，the formula of energy dissipation function of river can be obtained.Energy dissipation rate of river system is obtained by volume integral for energy dissipation function.In the entire derivation，the use of many complex equations of non-equilibrium thermodynamics and advanced fluid mechanics can be avoided，making it simple and easy to understand.

river；non-equilibrium thermodynamics；generalized flux；generalized force；energy dissipation rate

TV147

A

0493-2137（2015）12-1126-04

10.11784/tdxbz201406020

2014-06-08；

2015-04-10.

國(guó)家自然科學(xué)基金創(chuàng)新研究群體科學(xué)基金資助項(xiàng)目（51321065）；國(guó)家自然科學(xué)基金資助項(xiàng)目（50679053）.

趙麗娜（1983—），女，博士研究生，neimengguzhaolina@163.com.

徐國(guó)賓，xuguob@sina.com.