雙噴嘴水平對置撞擊流混合器內湍流流動及混沌特性

張建偉，董鑫，馬紅越，馮穎

（沈陽化工大學能源與動力工程學院，遼寧 沈陽 110142）

撞擊流（impinging streams，簡稱IS）概念于1961 年由白俄羅斯的科學家Elperin 首先提出[1]，其特點是可顯著強化氣固、氣液、氣氣以及液液系統的熱、質傳遞過程。撞擊流技術具有促進微觀混合和強化相間傳遞的特點，并且還在制取超細粉體方面的性能優于普通攪拌槽混合器[2]，在燃燒、干燥、粉碎、研磨和吸收方面也產生良好的效果[3]。撞 擊流混合器是一種新型高效的反應、混合設備， 其內部流場的運動行為與上述現象有著緊密的 關系。

混合設備內強化動量傳遞過程往往與其流場特性密切相關，只有掌握了設備內的流動規律才能有效地分析流體微團之間的混合和傳遞特性。因此，撞擊流混合器內流場特性研究備受關注。考慮到撞擊流混合器具有優良傳質與傳熱行為的關鍵是兩股流體高速撞擊形成徑向流動強化流體微團湍動和混合，通過對撞擊流混合器內壓力波動信號和濃度值的瞬態特性進行研究，發現混合器內的瞬態流動具有明顯的非線性特征[4-7]。本文作者前期研究表 明[8-10]，基于混合器內各種波動信號的小波分析和分形理論認識到反應器中流體動力學行為的自相似分形特征和復雜性特征。由于撞擊流混合器內流域的非線性很大程度體現在混合器內瞬時速度場的波動，因此研究撞擊流混合器內湍流流動及混沌特性具有重要的實際意義。本工作根據激光多普勒測速儀測得的速度場，利用湍流理論計算得到湍流特性參數，并且利用重構相空間理論恢復速度信號的吸引子，進而得到混沌特征參數，研究了撞擊流混合器結構參數以及進口處雷諾數對湍流流動及混沌特性的影響。

圖1 實驗測量系統示意圖

1 實驗裝置與測量

1.1 流場測量系統

采用激光多普勒測速儀測量了有機玻璃筒和噴嘴組成的雙噴嘴對置撞擊流混合器內的流場，實驗系統流程如圖1所示，有機玻璃圓筒直徑為140mm，壁厚為5mm，高為500mm，噴嘴固定在有機玻璃筒兩側中間位置。實驗裝置是由泵、蓄水箱、轉子流量計、閥門、雙噴嘴水平對置撞擊流混合器以及其連接管線組成的介質循環系統，采用由丹麥DANTEC 公司生產的三維五光束多普勒測速儀對瞬時速度場進行測量，此測速儀則由2017 型水冷氬離子激光器、多色激光束分離器、三維五光束光學發射接收探頭、光電接收器、三維信號處理器、三維移動坐標系統和主控計算機組成。在混合器外部設置有方型夾套，在夾套內注入水以補償折射的影響，可以減小因為折射引起的測量誤差[11]。實驗所用介質選取蒸餾水，配有示蹤粒子，示蹤粒子為空心玻璃球，其密度為1.05～1.15g/cm3，直徑為10～15μm。

1.2 實驗裝置介紹

為便于實驗測量，坐標系選取直角坐標系。以兩對置噴嘴軸線中點為原點，兩對置噴嘴軸線方向為x 軸，兩對置噴嘴水平徑向為y 軸，有機玻璃筒體軸線方向為z 軸，如圖2 所示。通過撞擊流體入口b1,2不斷注入流體，液體排凈口c 處于關閉狀態，待流體液面位于對置噴嘴上方并保持浸沒狀態開始進行實驗。文中實驗的噴嘴間距設定為1d、3d 和5d，進口處雷諾數Re 分別為6050、7260 和8470。其中d、L 和分別表示噴嘴直徑和噴嘴間距。

圖2 實驗裝置示意圖（單位：mm）

2 撞擊流混合器內的流動特性分析

2.1 混合器內速度場的分析

實驗采用激光多普勒測速儀測量了進口處雷諾數Re 為6050 時噴嘴間距L 分別為1d、3d 和5d 的撞擊流混合器內瞬時速度場。比較坐標系下量綱為1 的徑向速度在撞擊面中心線上的分布，如圖3 所示。由于兩股流體碰撞后，流體的運動形式轉變成徑向射流，即沿著z 軸兩側方向分布，其速度分布趨勢關于直線z=0 近似反對稱。點z=0 附近區域速度近似為零，形成駐點，可視為撞擊中心區。由于噴嘴間距較小時，兩股流體撞擊不完全，能量傳遞不明顯，而噴嘴間距較大時，兩股流體撞擊前產生能量損失，碰撞時速度變化較大，但是產生的徑向速度不大。所以隨著噴嘴間距的增加，徑向速度先增加后減小。L=3d 為較適宜的噴嘴間距。

實驗測量了噴嘴間距L 為3d、雷諾數Re 分別為6050、7260 和8470 的撞擊流混合器內速度場。對比坐標系下量綱為1 的徑向速度在撞擊面中心線上的分布如圖4 所示。與圖3 的徑向分布相同，由于兩股流體碰撞后形成沿z 軸兩側方向徑向射流， 其速度分布趨勢關于直線z=0 近似反對稱。而且在點z=0 附近區域速度近似為零，隨著雷諾數的增加，量綱為1 的徑向速度呈遞增趨勢。這是因為雷諾數較大時，兩流體射流速度較大，在碰撞時產生的較大的能量傳遞，從而使徑向速度逐漸增加。

圖3 不同噴嘴間距下的徑向速度

圖4 不同雷諾數下的徑向速度

圖5、圖6 分別為量綱為1 的軸向速度在截面中心線上各點不同噴嘴間距和不同雷諾數下的對比結果。其分布規律基本相似，都近似關于直線z=0對稱，在點z=0 附近達到各自的極小值，此處為駐點。如圖5 所示，由于噴嘴間距較小時，兩股流體撞擊不完全，能量傳遞不明顯，而噴嘴間距較大時，兩股流體撞擊前產生能量損失，碰撞時速度變化較大，但是產生的軸向速度不大。所以隨著噴嘴間距的增加，軸向速度先增加后減小。L=3d 為較適宜的噴嘴間距。圖6 中，隨著雷諾數的增加，量綱為1軸向速度呈遞增趨勢。隨著雷諾數的增加，軸向速度逐漸增大。

從圖3～圖6 可以看出，截面內無因次化軸向速度和量綱為1 化徑向速度的分布形式有明顯的區別。而不同噴嘴間距和不同雷諾數下軸向速度的分布形式相同，徑向速度亦是如此。

圖5 不同噴嘴間距下的軸向速度

圖6 不同雷諾數下的軸向速度

2.2 混合器內流體的湍流流動特性

2.2.1 湍流特性描述參數

在湍流狀態下，流場中任意一點的速度是脈動的。描述湍流流場中某點速度時，通常用時均速度來表示這一點的運動速度，而描述流場湍流特性時，一般用與脈動速度相關的參數來表示，常用的有速度脈動均方根（RSM）、湍流強度、湍動能等。這些參數的表達式分別為式（1）～式（4）。

圖7 不同噴嘴間距下的湍流特性參數

2.2.2 不同噴嘴間距對湍流參數的影響

為了具體說明噴嘴間距的變化對撞擊流混合器內徑向速度場的影響，現對撞擊面上沿z 軸方向x=0處的湍流特性參數進行研究。徑向速度脈動均方根、湍流強度、湍動能隨噴嘴間距變化的趨勢圖如圖7所示（其中，噴嘴間距L 分別為1d、3d 和5d）。

圖7(a)是不同噴嘴間距下徑向速度脈動均方根沿z 軸的分布曲線，其分布規律基本相似；分布曲線都關于直線z=0 近似對稱，在點z=0 附近達到各自的極小值。速度脈動均方根隨著噴嘴間距的增加，呈先增加后減小的趨勢。這是因為噴嘴間距較大時，撞擊速度較小，撞擊較不劇烈，間距過小，撞擊不完全，可以得到利于混合的合適的噴嘴間距。

在圖7(b)中，各個噴嘴間距下的湍流強度分布規律相似，即關于z=0 反對稱，都在徑向z=-10～25mm 之間的湍流強度有明顯的變化。隨著噴嘴間距的增加，其值亦先增加后降低。

由圖7(c)可知，湍動能的分布規律與徑向速度脈動均方根的分布相近。其分布規律基本相似：分 布曲線都關于直線z=0 近似對稱，在點z=0 附近達到各自的極小值，并且湍動能隨著噴嘴間距的增加，其值亦先增加后降低。分析主要是噴嘴間距增加使徑向速度脈動被弱化所致。

2.2.3 不同雷諾數對湍流參數的影響

為了具體說明雷諾數的變化對撞擊流混合器內速度場的影響，現對撞擊面上沿z 軸方向x=0 處的湍流特性參數進行研究。徑向速度脈動均方根、湍流強度、湍動能隨雷諾數變化的趨勢圖如圖8 所示（其中，雷諾數Re 分別為6050、7260 和8470）。

圖8(a)是不同雷諾數下徑向速度脈動均方根沿z 軸的分布曲線，其分布規律基本相似；分布曲線都關于直線z=0 近似對稱，在點z=0 附近達到各自的極小值。而速度脈動均方根隨著雷諾數的增加，峰值分別為0.14m/s，0.15m/s 和0.22m/s，峰值分別為0.14m/s，0.15m/s 和0.22m/s，其速度脈動呈增大趨勢，尤其雷諾數Re=8470 對應的徑向速度脈動均方根較其他兩種情況下的脈動均方根有明顯變化。這是因為雷諾數較大時，撞擊速度較大，撞擊劇烈。

圖8 不同雷諾數下的湍流特性參數

在圖8(b)中，不同雷諾數下的湍流強度分布規律相似，即關于z=0 反對稱，都在徑向z=-10～30mm之間的湍流強度有明顯的變化，其余位置對應的湍流強度變化幅度不大。隨著雷諾數增大，湍流強度也隨之增大；表明此時雷諾數增加對湍流的強化作用增強，分析主要是由于雷諾數增加使速度脈動被大大強化所致。

由圖8(c)可知，湍動能的分布規律與徑向速度脈動均方根的分布相近。其分布規律基本相似；分布曲線都關于直線z=0 近似對稱，在點z=0 附近達到各自的極小值。湍動能隨著雷諾數的增加，其增加的幅度比脈動速度增加的幅度更大。分析主要是雷諾數增加使徑向速度脈動被強化所致。增加雷諾數可增加高湍動能分布區域，該區域更有利于物料混合。

2.3 混合器內流體的混沌特性

2.3.1 混沌理論

混沌理論在徑流預報[12]、經濟預測[13]、通信信號調制識別[14]等領域中得到廣泛應用，并顯示出獨特優點，并且近年來應用于內部流場的分析。分析混沌時間序列，首先需要進行混沌特性識別，其主要方法是通過計算混沌信號奇異吸引子的特征參數來辨別，關聯維是吸引子維數及其幾何結構復雜程度的表征，Kolmogorov 熵則是關于混沌系統的初始信息損失速率的度量，對動力學復雜性的度量描述，最大Lyapunov 指數刻畫混沌軌道的局域不穩定性和對初始條件敏感依賴性的定量判據。

首先采用由Packard 等[15]提出的相空間理論重構出瞬時速度信號的吸引子，利用適用于實驗觀測數據且計算方法簡潔的G-P 算法[16-17]對分數維的值進行計算。然后，選用Grassberger 和Procaccia[18]根據q 階Renyi 熵[19]提出的關聯熵K 的算法來估計K 熵。最后，選用Wolf 方法[20]對最大Lyapunov 指數進行估值計算。

2.3.2 重構相空間的參數選取

延遲時間τ 是平衡重構相空間各點間獨立性與相關性的重要參數，若延遲時間τ 過大，則各向量之間無相關性，而若延遲時間τ 過小，重構相空間中的各點不能作為獨立向量應用于計算。互信息法既評價信號的獨立性，也考慮了線性相關與非線性相關對信號的影響[21]，并且互信息法的有效性得到眾多相關研究的證實[22]，所以應用互信息法計算延遲時間。

為了令重構吸引子等價于原吸引子拓撲，應該確定最佳嵌入維數。若嵌入維數過小，吸引子則不能容納系統的吸引子，導致系統的動力特性不能全面展示；若嵌入維數太大，可能將噪聲引入到額外的維數上，從而增加額外的計算量和計算難度。因此需要確定恰當的嵌入維數。為確定最佳嵌入維數和線性區間，現選擇流場內（0，27，-38）位置處的瞬時速度，根據嵌入維數從1遞增到12進行計算，得到關聯積分分布曲線圖9(a)所示。

圖9 不同嵌入維下關聯積分曲線和關聯

圖9(a)中所示的曲線分別為嵌入維數m 由1 至12（從上至下）對應的關聯積分曲線，隨著關聯維數增加，圖中的關聯積分曲線均收斂于同一斜率，可初步確定此非線性系統的混沌吸引子的收斂性。在圖中ln r ∈[ - 2.5, - 1.8]區間內的大多數關聯積分曲線，均具有相對一致的斜率，則可以認定該區間為線性區間，并且其中存在關聯維。在線性區間內，擬合各點對應不同嵌入維數的關聯積分曲線的斜率，然后獲得不同嵌入維下的關聯維值。利用該方法即可計算撞擊面上不同位置處的5 個點的關聯維，從而得到各點對應各階嵌入維下的關聯維值如圖9(b)所示。嵌入維m≤8 下各點關聯維值隨m 的增加呈線性增大，顯示出一致性，而嵌入維m≥9時，由于嵌入維過大，導致關聯維值在撞擊區氣泡的干擾下出現不穩定的波動。因此，本文均按嵌入維m=8 計算混沌特征參數。

2.3.3 不同噴嘴間距對混沌特性的影響

進口處雷諾數Re 為6050、噴嘴直徑為12mm時，分別在兩噴嘴間距L 為1d、3d、5d 時，對撞擊面進行PDA 測量，計算所有測量點軸向瞬時速度的混沌特征參數：關聯維、Kolmogorov 熵和最大Lyapunov 指數，結果如圖10 所示。為確定不同的 噴嘴間距對流場混沌特征參數大小的影響，本文作者取不同噴嘴間距下同一直線處各點的混沌特征參數進行對比，如圖11 所示（測量點的軸向和垂直徑向坐標分別為x =0mm、z =-22mm）。

圖10 噴嘴間距L=1d 時撞擊面的混沌特征參數

在圖10(a)中，關聯維值在2.32～6.84，同樣可以得出此工況下流場具有自相似結構且存在分形。同時在撞擊面上從z =30mm 到z =-60mm 區域內的關聯維值普遍高于z <-60mm 的區域，說明此區域各點流體的運動行為比較復雜，因為當兩股流體相向高速流動撞擊時，反復做減幅振蕩運動，同時由于撞擊過程是在浸沒狀態下進行，因此兩股流體的震蕩運動會改變其周圍流體的流動狀態，從而導致此區域內大部分流體的運動行為比較復雜，關聯維值偏大。可以判斷該區域為撞擊流混合器的撞擊區。

圖10(b)中，在z =10mm 至z =-10mm、y=5mm至y=10mm 區域內（紅色環形區域）Kolmogorov熵最大，此區域流體的初始能量以較大的速率耗散，可認為該區域是兩股流體在浸沒狀態下的撞擊中心面。因為兩股高速湍動流體劇烈撞擊后，改變原有的運動趨勢，由軸向運動轉化為徑向運動并向四周擴散，因此原有的流動狀態減弱變為其他形式的流動狀態，促使兩股流體達到彼此干擾、混合的狀態。在撞擊面上從y=40mm 至y=45mm 區域內（黃色區域），其Kolmogorov 熵明顯略大于除撞擊中心面的其他區域，此處是由軸向流動轉為徑向流動的流體經擴散后碰撞混合器壁面，隨后做反向運動并與向壁面方向擴散的流體發生二次碰撞而消耗能量所引起的，此時也起到了二次混合的效果。因此，通過觀察整個撞擊面上Kolmogorov 值分布狀況，可以得到混合器內流場速度信息耗散率較大的區域，同時驗證了撞擊流混合器內流體的流動趨勢。

圖10(c)中，整個撞擊面上各點的Lyapunov 指數值都大于零，證明該撞擊流混合器的流場具有混沌特性。縱觀整個撞擊面上的Lyapunov 指數值，其數值分布比較均勻，說明整個撞擊面內測量點的速度對初始條件均比較敏感，依賴性較大。

圖11 結果顯示，混沌特征參數（關聯維、Kolmogorov 熵和最大Lyapunov 指數）的3 個數值均隨噴嘴間距的增大，呈現出相同的變化規律。當噴嘴間距L 由1d 增加到3d 時，3 個數值均明顯增大，這是因為隨著噴嘴間距的增大兩股流體撞擊的接觸面積增大，從而兩股流體的初始信息損失速率增大，這有利于微觀混合。當噴嘴間距再次增加到5d 時，3 個數值卻明顯減小，這是因為當噴嘴間距增加到一定距離，兩股流體相撞的劇烈程度減小，產生的震蕩程度削弱，不利于微觀混合的進行。因此，當噴嘴間距L 為3d 時，此撞擊流混合器的微觀混合效果最好。

2.3.4 不同雷諾數對混沌特性的影響

噴嘴間距L 為3d，噴嘴直徑為12mm 時，進口雷諾數Re 分別為6050、7260 和8470，對撞擊面進行PDA 測量，計算所有測量點軸向瞬時速度的混沌特征參數：關聯維、Kolmogorov 熵和最大Lyapunov 指數，結果如圖12 所示。為確定不同的噴嘴間距對流場混沌特征參數大小的影響，本文作者取不同雷諾數下同一直線處各點的混沌特征參數進行對比，如圖13 所示（直線上各測量點的軸向和垂直徑向坐標分別為x = 0mm、z =10mm）。

圖11 不同噴嘴間距下撞擊面的混沌特征參數

圖 12 中關聯維、Kolmogorov 熵和最大Lyapunov 指數三者的數值云圖與圖10 具有大致相同的分布規律，此工況下流場具有自相似結構且存在分形，證明該撞擊流混合器的流場具有混沌特性，因此可以得出不同操作條件下撞擊流混合器具有相同的流場特性和運動趨勢，驗證了改變單一操作條件進行對比實驗的合理性。

在圖13(a)關聯維隨噴嘴進口雷諾數變化的3條曲線中，位于最上方的曲線對應的操作工況為Re=8470，其次為Re=7260，最下方的曲線對應為Re=6050。

再觀察圖13(b)和圖13(c)可發現，Kolmogorov 熵和最大 Lyapunov 指數中3 條曲線的排列位置與關聯維中的排列位置基本一致，除了個別點的值有所跳躍，這是由于撞擊區流體的復雜性和不穩定性造成的，但并不影響3 條曲線的整體排列趨勢。

圖12 雷諾數Re=6050 時撞擊面的混沌特征參數

通過上述現象可以得出，在噴嘴直徑和噴嘴間距保持不變的工況下，3 個混沌特征參數（關聯維、Kolmogorov 熵和最大 Lyapunov 指數）隨著單一變量噴嘴進口處雷諾數的增加，其數值具有相同的變化規律，均逐漸增大，撞擊流混合器的微觀混合效果逐步增強。這是因為流體初始速度的增大提高了流場的復雜程度和湍流動能，導致兩股流體在進行撞擊混合時，初始信息損失速率增大，彼此達到更好的碰撞、擾動，進行信息互換至最終混合。

3 結 論

圖13 不同雷諾數下混沌特征參數的變化趨勢

（1）采用湍流理論和混沌理論對浸沒撞擊流混 合器內速度信號進行分析，得到流場的關聯維值均在2.32～7.52，流場具有自相似結構且存在分形。并且其Kolmogorov 熵和最大Lyapunov 指數均大于零，因此得出該雙噴對置撞擊流混合器速度場具有混沌特性。

（2）噴嘴間距較大時，撞擊速度較小，撞擊較不劇烈，而噴嘴間距較小時撞擊不夠完全，噴嘴進口處雷諾數和噴嘴直徑保持不變時，湍流特性參數和混沌特性參數均隨著的噴嘴間距增大現增加后減小，可以得到利于混合的噴嘴間距L=3d。

（3）由于進口處雷諾數的增加，流體初始速度的增大提高了流場的復雜程度和湍流動能，噴嘴間距和噴嘴直徑保持不變時，湍流特性參數和混沌特性參數均隨著噴嘴進口處雷諾數的增加而成遞增趨勢，雷諾準數有助于提高撞擊流混合器的微觀混合效果。

符 號 說 明

D—— 分形維數

不久，竹韻悠悠醒來了，臉色慘白地對龍斌笑了笑，站起來走到電腦前，在百度搜索鍵上打下了雜志文章的標題，一搜索又令她氣憤難當，不但這篇文章上了網，而且點擊率，跟貼率極高。許多不明真相的網友大罵竹韻是騷貨，是禍水，要將她拉出去斬立決，什么難聽的話都有。竹韻突然覺得自己掉進了一個污水池里，淹得快要窒息了……

d—— 噴嘴直徑，mm

Ia—— a 方向湍流強度

k—— 湍動能，m2/s2

L—— 噴嘴間距，mm

m—— 嵌入維

N—— 測量次數

Re—— 雷諾數

u—— 時均總速度，m/s

uai—— 第i 次的該點a 方向瞬時速度，m/s

λ—— 最大 Lyapunov 指數

τ—— 延遲時間，s

[1] Elprein I T. Heat and mass transfer in opposing currents[J]. Engng. Physics，1961，6：62-68.

[2] 伍沅，陳煜. 撞擊流反應制取“超細”白炭黑[J]. 化工學報，2003，54（10）：1381-1386.

[3] 伍沅. 撞擊流——原理.性質.應用[M]. 北京：化學工業出版社，2006.

[4] 張建偉，伍沅，舒安慶，等. 浸沒循環撞擊流反應器的壓力脈動特性[J]. 化工學報，2005，56（2）：266-269.

[5] 張建偉，焦麗. 撞擊流反應器壓力波動的多尺度多分形特征分析[J]. 化工學報，2006，57（7）：1558-1564.

[6] 林柯利，畢榮山，譚心舜. 利用面激光誘導熒光技術研究噴射器內液-液湍流混合特性[J]. 高校化學工程學報，2009，23（1）：28-33.

[7] 張建偉，殷婉君，孫熙同，等. 基于希爾伯特-黃變換的撞擊流混合器濃度場特性分析[J]. 高校化學工程學報，2014，28（5）：971-978.

[8] Zhang Jianwei，Wang Xiangwang，Wang Chunyu. Multi-scale and multi-fractal characteristics of instantaneous velocity signals in impinging stream mixer[J]. The Canadian Journal of Chemical Engineering，2010，88（2）：172-178.

[9] 張建偉，宋艷芳. 撞擊流反應器流場速度時間序列的結構復雜性分析[J]. 機械工程學報，2011，47（22）：174-180.

[10] 張建偉，馬彥東，王顯旺，等. 撞擊流混合器內速度脈動的間歇性分析[J]. 實驗流體力學，2012，26（2）：6-11.

[11] 毛德明，馮連芳，許國軍，等. 用LDA研究攪拌釜內的流場[J]. 高校化學工程學報，1996，10（3）：258-263.

[12] 趙永龍，丁晶，鄧育任. 混沌分析在水文預測中的應用與展望[J].水科學進展，1998，9（2）：181-186.

[13] 馬軍海，賈湖，盛昭瀚. 非線性經濟時序的預測方法及其應用研究[J]. 管理科學學報，2001（4）：49-54.

[14] 何建華，楊宗凱，王殊. 基于混沌和神經網絡的弱信號探測[J].電子學報，1998，26（10）：33-37.

[15] Packard N H ，Crutchfield J R，Farmer J D，et al. Geometry from a time series[J]. Rev. Letters，1980，45：712- 716.

[16] Grebogi C，Ott E，Yorke J A. Final state sensitivity：An obstruction to predictability[J]. Pys. Lett. A.，1983，99：415-419.

[17] Grassberger P，Procaccia I. Characterization of strange attracors[J]. Phys. Rev. Lett.，1983，50：346-349.

[18] Grassberger P，Procaccia I. Estimation of the Kolmogorov entropy from a chaotic signal[J]. Phys. Rev. A，1983，28：2591-2593.

[19] Renyi A. Probability Theory[M]. Appendix. Amsterdam ：North Holl and Press，1970.

[20] Wolf A，Swift J B，Swinney H L，et al. Determining Lyapunov exponents from a series[J]. Physical D，1985，16：285-317.

[21] Fraser A M，Swinney H L. Independent coordinates for strange attractors form mutual information[J]. Phys. Rev. A，1986，33（2）：1134-1140. 趙貴兵，石炎福，段文鋒，等. 從混沌時間序列同時計算關聯維和Kolmogorov熵[J].計算物理，1999，16（3）：309-315.