王貴全

(荷蘭特文特大學(xué) 流體物理課題組，恩斯赫德 7500AE，荷蘭)

0 引 言

慣性粒子在壁面湍流中的流動(dòng)可用于研究沙粒在河流中的沉積、雨雪的沉積速度與分布、沙塵在地球以及火星近地面層中的傳輸以及空氣污染物在城市上空的分布等很多與地球物理、工程領(lǐng)域以及環(huán)境科學(xué)相關(guān)的重要問題[1-4]。實(shí)驗(yàn)研究主要側(cè)重于流體/粒子的平均速度和雷諾應(yīng)力，以及粒子的分布，然而由于顆粒聚集與群聚效應(yīng)以及壁面的影響，導(dǎo)致實(shí)驗(yàn)研究的難度增加，限制了實(shí)驗(yàn)研究的發(fā)展[5]。隨著計(jì)算能力的不斷提高，直接數(shù)值模擬壁面湍流耦合慣性點(diǎn)粒子可以達(dá)到中等湍流雷諾數(shù)Reτ～O(103)（根據(jù)壁面黏性尺度），耦合顆粒數(shù)目可以達(dá)到O(108)以上，可以看到直接數(shù)值模擬已經(jīng)成為一個(gè)有效的技術(shù)手段被應(yīng)用于探究湍流粒子兩相流動(dòng)這個(gè)非常復(fù)雜的物理問題，其不僅可以捕捉到最小的流動(dòng)尺度，分析顆粒在湍流中的分布狀態(tài)，也可以用于研究粒子對(duì)湍流的調(diào)制。

對(duì)于慣性粒子，本文假設(shè)顆粒密度遠(yuǎn)大于流體相密度，顆粒尺度遠(yuǎn)小于湍流最小耗散尺度，在上述假設(shè)條件下，慣性粒子可被近似為點(diǎn)粒子，這樣可以極大地簡(jiǎn)化顆粒相的數(shù)值模型。而隨著顆粒相濃度的下降，粒子與粒子間的相互作用力以及粒子對(duì)流體間的反作用力也可以逐漸的被忽略，這樣同樣也會(huì)提高計(jì)算效率。Elghobashi[6]根據(jù)顆粒相濃度的不同，將粒子與流動(dòng)的耦合分為：四向耦合（同時(shí)考慮粒子對(duì)流體的反饋力以及粒子之間的相互作用力）、雙向耦合（不考慮粒子之間的相互作用力）以及單向耦合（不考慮粒子對(duì)流體的反饋力以及粒子之間的相互作用力）。本文應(yīng)用直接數(shù)值模擬和拉格朗日點(diǎn)粒子跟蹤方法對(duì)慣性粒子在壁面湍流這一復(fù)雜問題進(jìn)行研究，主要關(guān)注粒子對(duì)湍流的調(diào)節(jié)機(jī)制及粒子分布。對(duì)湍流的調(diào)制主要研究粒子慣性對(duì)湍流內(nèi)層大尺度結(jié)構(gòu)LSMs以及外層超大尺度結(jié)構(gòu)VLSMs的影響。對(duì)粒子分布的研究包含兩部分：一是粒子聚集，這里指粒子在壁面垂直方向由于湍流場(chǎng)導(dǎo)致的濃度不均勻；二是粒子群聚形態(tài)，這里指粒子在水平截面上由于湍流場(chǎng)導(dǎo)致的不均勻分布形態(tài)。本文是對(duì)作者以往開展的相關(guān)研究工作進(jìn)行小結(jié)[7-13]，結(jié)構(gòu)安排如下：第1節(jié)介紹此工作采用的數(shù)值方法及驗(yàn)證過程（單向、雙向及四向耦合），第2節(jié)討論粒子對(duì)湍流內(nèi)層LSMs與外層VLSMs的調(diào)制機(jī)理（雙向耦合），第3節(jié)分析粒子在LSMs與VLSMs中的分布（雙向耦合），第4節(jié)關(guān)注粒子在重力驅(qū)動(dòng)下的沉積速度（單向耦合）。

1 數(shù)值模擬方法及驗(yàn)證

1.1 數(shù)值方法

1.1.1 流體相

采用直接數(shù)值模擬對(duì)不可壓縮流體求解，流向（x）和展向（z）采用周期性邊界條件并應(yīng)用譜方法求解，壁面垂直方向（y）采用二階有限差分格式，時(shí)間推進(jìn)采用三階Runge-Kutta格式。不可壓縮條件通過求解壓力泊松方程來滿足。求解的連續(xù)方程和動(dòng)量方程分別為：

其中ui為流體速度，p為壓力，F(xiàn)i為粒子對(duì)流體相的反作用力，v是 流體動(dòng)力學(xué)黏度， ρf是流體的密度。

1.1.2 顆粒相

當(dāng)粒子的密度遠(yuǎn)大于流體相密度、粒子的尺度遠(yuǎn)小于最小的湍流耗散尺度，粒子被近似為點(diǎn)粒子。此時(shí)只考慮Schiller-Naumann曳力模型，粒子的速度和軌跡可以通過積分如下方程：

在研究壁面湍流耦合慣性粒子這一問題時(shí)，最重要的控制參數(shù)為粒子的斯托克斯數(shù)（Stokes number,St），定義為粒子的特征時(shí)間尺度（ τp）和流場(chǎng)特征時(shí)間尺度（ τf）的比值。而由于壁面湍流邊界層內(nèi)部區(qū)域與外部區(qū)域存在不同的流場(chǎng)特征時(shí)間尺度，為了研究粒子與不同湍流尺度的耦合作用，定義了如下三個(gè)粒子斯托克斯數(shù)：

1.2 模型驗(yàn)證[7]

在二維耦合模型中，點(diǎn)粒子方法是將點(diǎn)力耦合到最近的歐拉網(wǎng)格點(diǎn)上，忽略粒子體積對(duì)流體的擾動(dòng)，除了上述傳統(tǒng)的點(diǎn)粒子方法，Capecelatro[14]介紹了一種體積平均方法，在計(jì)算粒子力的時(shí)候，不僅考慮曳力模型，同時(shí)考慮粒子大小對(duì)本地壓力和黏性應(yīng)力的影響，通過高斯核函數(shù)耦合到歐拉網(wǎng)格點(diǎn)上。Gualtieri[15]介紹了一種正則化點(diǎn)粒子方法，在計(jì)算中通過多極展開獲得球的斯托克斯流動(dòng)的一階偶極子來近似點(diǎn)粒子對(duì)遠(yuǎn)場(chǎng)的擾動(dòng)。上述兩種或其他的改進(jìn)點(diǎn)粒子方法一般都是嘗試通過引入模型來考慮粒子的體積效應(yīng)，但弊端是大大增加了計(jì)算量。為了驗(yàn)證傳統(tǒng)的點(diǎn)粒子方法在模擬壁面湍流耦合粒子問題時(shí)的準(zhǔn)確性，我們對(duì)傳統(tǒng)點(diǎn)粒子方法與體積平均方法[14]以及實(shí)驗(yàn)結(jié)果[16]進(jìn)行了比較，驗(yàn)證過程及參數(shù)如下：

（1） 首先與已發(fā)表的采用傳統(tǒng)點(diǎn)粒子方法的結(jié)果[17]進(jìn)行比較，目的是驗(yàn)證代碼編寫的正確性。參數(shù)為Reτ= 180、St+= 30、雙向耦合。與文獻(xiàn)[17]比較了流體相與粒子相的平均速度及脈動(dòng)速度，得到了幾乎一致的結(jié)果；

（2） 進(jìn)一步比較傳統(tǒng)點(diǎn)粒子方法與體積平均方法[14]，來了解忽略顆粒體積效應(yīng)對(duì)模擬結(jié)果的影響。參數(shù)為Reτ= 630、St+= 2030以及Reτ= 227、St+= 58.6，四向耦合。與文獻(xiàn)[14]比較了流體相與粒子相的平均速度及脈動(dòng)速度，得到了幾乎一致的結(jié)果；

（3） 將傳統(tǒng)點(diǎn)粒子方法、體積平均方法[14]與實(shí)驗(yàn)結(jié)果[16]進(jìn)行一對(duì)一的比較。參數(shù)為Reτ= 227、St+=58.6，模擬中采用四向耦合。比較了粒子濃度曲線，粒子聚集與群聚等穩(wěn)態(tài)值（箱形計(jì)數(shù)法，馮洛諾伊圖，概率密度函數(shù)，徑向分布函數(shù)），發(fā)現(xiàn)兩種模擬得到的結(jié)果近乎一致。但與實(shí)驗(yàn)結(jié)果相比，尤其是在近壁面處，粒子的壁面垂直脈動(dòng)速度以及粒子的濃度有差別，可能的原因有—數(shù)值模擬中沒有正確考慮湍流湍泳以及近壁面升力的影響[18]，或者實(shí)驗(yàn)中靜電效應(yīng)以及近壁面由于粒子聚集導(dǎo)致測(cè)量不夠精確。對(duì)點(diǎn)粒子模型在近壁面處的修正是個(gè)一直沒有解決的問題，一種可能適用的方法是將小粒子看作有限大小粒子進(jìn)行捕捉界面研究，但是小粒子的尺度遠(yuǎn)小于湍流最小耗散尺度，這樣做會(huì)導(dǎo)致難以承受的計(jì)算網(wǎng)格數(shù)量，據(jù)作者所知，目前還缺乏相關(guān)的研究工作。

綜上，在粒子壁面湍流中，在St+= O(10～103)區(qū)間內(nèi)，傳統(tǒng)的點(diǎn)粒子方法和改進(jìn)的點(diǎn)粒子方法獲得的粒子的速度與濃度一階、二階穩(wěn)態(tài)量以及流體相的速度一階、二階和三階穩(wěn)態(tài)量幾乎一致。由于傳統(tǒng)點(diǎn)粒子方法可以節(jié)省大量的計(jì)算時(shí)間，使得在直接數(shù)值模擬中可以引入更多數(shù)量的粒子，因此在本文的研究中，將采用這種數(shù)值方法。

1.3 計(jì)算域尺度選擇[8]

湍流內(nèi)層的大尺度結(jié)構(gòu)（LSMs）和湍流外層的超大尺度結(jié)構(gòu)（VLSMs）的存在導(dǎo)致流體的一階與二階穩(wěn)態(tài)值依賴于計(jì)算域的尺寸[19]。圖1展示了瞬時(shí)流向脈動(dòng)速度的云圖在湍流內(nèi)層與外層、四種計(jì)算域尺寸下的模擬結(jié)果：超大尺寸Lx= 12πh，Lz= 3πh；大尺寸Lx= 6πh，Lz= 3πh；中尺寸Lx= 6h，Lz= 3h；小尺寸Lx= 2.5h，Lz= 1.5h。在湍流內(nèi)層（y+= 55處水平截面），可以看到，小計(jì)算域圖1（d）只能捕捉一對(duì)高低速條帶（高低速條帶是由于流向渦結(jié)構(gòu)將近壁面低速流體“上揚(yáng)”到遠(yuǎn)壁面的位置，或?qū)⑦h(yuǎn)壁面高速流體“下掃”到近壁面的位置），中等以上的計(jì)算域圖1（a-c）可以得到大于三對(duì)以上的高低速條帶。在湍流外層（y+= 275處水平截面）：小計(jì)算域得不到正確的高低速條帶；中等計(jì)算域在展向方向可以捕捉到一對(duì)高低速條帶，但是流向方向有截?cái)?；大尺度?jì)算域的高低速條帶幾乎與超大計(jì)算域的結(jié)果相同，在流向方向幾乎沒有截?cái)唷?/p>

在加入粒子后，粒子一階（平均速度/平均濃度）、二階穩(wěn)態(tài)（脈動(dòng)速度）也依賴于計(jì)算域的尺寸[20]，在已發(fā)表的研究成果中，所采用的計(jì)算域尺寸也不盡相同（文獻(xiàn)[8]的表1列出了文獻(xiàn)中常用的計(jì)算域尺寸），可以看到，研究計(jì)算域尺寸對(duì)粒子一階與二階穩(wěn)態(tài)值的影響是非常重要的。我們對(duì)Reτ= 550、Reτ= 950以及St+= 2.41～908進(jìn)行了模擬比較[8]。結(jié)果顯示在單相流中，中等尺寸計(jì)算域、大尺寸計(jì)算域和超大尺寸計(jì)算域的流體的一階、二階穩(wěn)態(tài)值幾乎相同，而小尺寸計(jì)算域會(huì)低估流體的一階、二階穩(wěn)態(tài)值。當(dāng)加入粒子以后，流體相二階穩(wěn)態(tài)量在小尺寸計(jì)算域與大尺寸計(jì)算域相差較大，顆粒相的近壁面濃度在中尺寸計(jì)算域會(huì)被低估（10%～20%），尤其是對(duì)中高慣性的粒子（St+= 24.2～182），但顆粒相的二階穩(wěn)態(tài)在中尺寸計(jì)算域和大尺寸計(jì)算域得到的結(jié)果幾乎一致。因此，想要精確得到粒子在湍流內(nèi)層的濃度，需要采用大尺寸計(jì)算域。而對(duì)于模擬粒子在湍流外層的濃度，中等尺度計(jì)算域能得到和大尺度計(jì)算域相似的結(jié)果。

圖1 兩個(gè)壁面垂直位置的橫截面的瞬時(shí)流向脈動(dòng)速度云圖，y+ = 55和y+ = 275分別代表湍流內(nèi)層與外層，（a-d）對(duì)應(yīng)于超大尺寸，大尺寸，中尺寸以及小尺寸[8]Fig. 1 Contours of instantaneous streamwise velocity fluctuation at two xz planes (y+ = 55 and y+ = 275) representing the inner and outer layers. Panels (a-d) refer to single-phase flow with four decreasing domain sizes, respectively [8]

2 慣性粒子對(duì)壁面湍流的調(diào)制機(jī)理

當(dāng)慣性粒子的濃度達(dá)到可以對(duì)流體相產(chǎn)生影響的時(shí)候，被調(diào)制的流體相將會(huì)反之影響粒子的分布與聚集[5]。在壁面湍流中，LSMs和VLSMs分別控制湍流內(nèi)層和外層的流體動(dòng)力學(xué)過程，因此我們將關(guān)注粒子對(duì)LSMs和VLSMs的調(diào)制機(jī)理[9-10]。為了突出粒子慣性的影響，在本節(jié)的計(jì)算中沒有考慮重力的影響。

2.1 對(duì)LSMs的調(diào)制機(jī)理[9]

一般地，湍流調(diào)制依賴于粒子大小和湍流尺度的比值，同時(shí)依賴于粒子的響應(yīng)時(shí)間尺度和流動(dòng)特征時(shí)間尺度的比值。Gore[21]提出了用粒子與湍流積分尺度比來判斷湍流抑制還是增強(qiáng)。Tanaka[22]提出了用流動(dòng)“含能尺度”和粒子斯托克斯數(shù)構(gòu)造的粒子動(dòng)量數(shù)來預(yù)測(cè)湍流調(diào)制。盡管沒有統(tǒng)一的機(jī)理來解釋粒子的湍流調(diào)制過程，但在湍流內(nèi)層中，小慣性粒子傾向于促進(jìn)湍流、大慣性粒子傾向于抑制湍流，這個(gè)推論被廣泛認(rèn)可。我們通過研究粒子對(duì)近壁面擬續(xù)結(jié)構(gòu)的調(diào)制作用，對(duì)其動(dòng)力學(xué)演化過程進(jìn)行細(xì)致分析，目標(biāo)是嘗試解釋粒子對(duì)湍流的調(diào)制機(jī)理[9]。

在湍流內(nèi)層，近壁面相干結(jié)構(gòu)的演化是一個(gè)復(fù)雜的非線性過程，對(duì)其機(jī)理的研究一直持續(xù)到現(xiàn)在。其中Hamilton[23]提出的湍流自維持過程被廣泛的用于研究LSMs的演化過程。Hamilton[23]應(yīng)用了Jiménez[24]提出的最小湍流單元概念到平板Couette流動(dòng)中，這樣做的好處是可以控制展向計(jì)算域來孤立出一組大尺度條紋結(jié)構(gòu)（large-scale streaks，LSS）和大尺度渦結(jié)構(gòu)（large-scale vortex，LSV），這樣就可以把注意力放在處于同一組自維持過程中的LSS和LSV的相互轉(zhuǎn)換過程上，而不受處于其他組自維持過程中的LSS和LSV的影響，同時(shí)方便做能譜分析。而采用平板Couette流動(dòng)來研究湍流自維持過程卻不采用Poiseuille流動(dòng)的好處是，在上下壁面之間只存在一組湍流自維持過程，即只存在一組LSS和LSV[28]。

圖2展示了湍流自維持過程發(fā)生的位置及三種結(jié)構(gòu)的演化過程。自維持過程發(fā)生于湍流內(nèi)層中，包括三種結(jié)構(gòu)：大尺度渦結(jié)構(gòu)（LSV），大尺度條帶結(jié)構(gòu)（LSS），以及流向相關(guān)條帶結(jié)構(gòu)（x-dependent streaks）。三種結(jié)構(gòu)之間的三個(gè)轉(zhuǎn)換過程分別為：LSS形成，LSS分解，LSV再生。假設(shè)從LSV（流向不相關(guān)）開始，流體由于LSV的作用將流體“上揚(yáng)”或“下掃”，形成了本地高速/低速的條紋結(jié)構(gòu)LSS（流向不相關(guān)），而LSS是不穩(wěn)定的，其失穩(wěn)產(chǎn)生流向相關(guān)的條帶，流向相關(guān)條帶引入流向速度梯度（?ux/?x），其與流向渦量（ωx）相互作用，對(duì)應(yīng)于渦量方程中的渦旋拉伸項(xiàng)，渦旋拉伸項(xiàng)的增加會(huì)促進(jìn)LSV的生成，進(jìn)而形成了一個(gè)完整的湍流自維持過程。應(yīng)該說明對(duì)湍流自維持過程的研究是在譜空間中進(jìn)行的，上述簡(jiǎn)單的物理描述相對(duì)應(yīng)的數(shù)學(xué)表達(dá)式請(qǐng)參考文獻(xiàn)[9,23]。

圖2 壁面湍流的內(nèi)外分區(qū)[25-26]以及自維持過程[23,27-28]. 在湍流Couette流動(dòng)的最小單元中，右上的云圖顯示了“上揚(yáng)”和“下掃”的LSS，向量圖顯示了一對(duì)旋轉(zhuǎn)方向相反的LSV[9]Fig. 2 Sketch of various wall regions and boundary layers in wall units[25-26] and the regeneration cycle sub-steps[23,27-28]. In the top-right is the flow field in turbulent pCf(plane Couette flow) in miniunit. Large-scale vortex (LSV) is shown by vector fields and colored by scalar value of the streamwise vorticity. Large-scale streak (LSS) is shown by the contour of the streamwise velocity magnitude[9]

近三十年的研究證明湍流自維持過程可以被用來闡述近壁面LSMs的演化機(jī)理，因此我們[9]對(duì)粒子慣性及濃度對(duì)湍流自維持過程的調(diào)制進(jìn)行了細(xì)致的研究，發(fā)現(xiàn)低慣性粒子（Stturb= 0.056，用LSV的周轉(zhuǎn)時(shí)間尺度對(duì)粒子響應(yīng)時(shí)間無量綱化）在LSV中的滯留時(shí)間更長(zhǎng)，高慣性粒子（Stturb= 0.625）會(huì)比較均勻的分布在湍流場(chǎng)中。相對(duì)應(yīng)的，低慣性粒子會(huì)促進(jìn)湍流并縮短湍流自維持過程的周期，高慣性粒子會(huì)抑制湍流并增長(zhǎng)湍流自維持過程的周期。由于低慣性（/高慣性）粒子的存在會(huì)增加（/減少）LSV的強(qiáng)度，導(dǎo)致近壁面湍流的增強(qiáng)(/減弱)。

以上，我們通過分析慣性粒子對(duì)湍流自維持過程的調(diào)制機(jī)理來研究粒子對(duì)近壁面大尺度結(jié)構(gòu)的影響。然而在實(shí)際應(yīng)用中，例如風(fēng)沙流流動(dòng)的雷諾數(shù)可以高達(dá)Reτ= O（106）[29]，導(dǎo)致湍流的內(nèi)層區(qū)間（y+＜100）只占據(jù)了整個(gè)湍流空間的很小一部分（～100/Reτ），因此除了研究粒子對(duì)近壁面LSMs的調(diào)制機(jī)理，有必要研究更高雷諾數(shù)下慣性粒子對(duì)壁面湍流外層的調(diào)制過程。

2.2 對(duì)VLSMs的調(diào)制機(jī)理[10]

隨著流動(dòng)雷諾數(shù)的增加，實(shí)驗(yàn)中發(fā)現(xiàn)了湍流中的另外一種重要的特征結(jié)構(gòu)，即超大尺寸結(jié)構(gòu)(VLSMs)[30]。近年來由于直接數(shù)值模擬可以達(dá)到更高的雷諾數(shù)，VLSMs也成為了數(shù)值研究湍流方向的熱點(diǎn)[31-33]。VLSMs的流向長(zhǎng)度可以達(dá)到20倍壁面高度而且含有很高的能量，可以攜帶40%～65%的湍動(dòng)能和30%～50%的雷諾應(yīng)力[34-35]。由于高雷諾數(shù)湍流中外層占據(jù)大部分的空間，因此研究慣性粒子對(duì)VLSMs的影響有重要的應(yīng)用價(jià)值，比如沙塵對(duì)大氣近壁面層湍流的影響[2]。然而到目前為止，慣性粒子的湍流調(diào)制數(shù)值模擬研究還局限在相對(duì)小的雷諾數(shù)(Reτ＜ 300)和對(duì)湍流內(nèi)層結(jié)構(gòu)的研究，因此我們?cè)谖墨I(xiàn)[10]中研究了中等雷諾數(shù)(Reτ= 550，Reτ= 950)湍流中粒子慣性對(duì)VLSMs的影響。

圖3展示了慣性粒子對(duì)VLSMs調(diào)制可能存在兩個(gè)路徑：1）直接路徑，慣性粒子通過影響湍動(dòng)能傳輸方程直接調(diào)制VLSMs；2）間接路徑，慣性粒子通過調(diào)制湍流內(nèi)層LSMs，進(jìn)而調(diào)制湍流外層VLSMs。在譜空間里，直接路徑代表慣性粒子直接調(diào)制與VLSMs同尺度的湍流結(jié)構(gòu)，可以通過譜方法來分析湍動(dòng)能傳輸方程進(jìn)行研究；間接路徑代表慣性粒子通過調(diào)制內(nèi)層LSMs來達(dá)到調(diào)制外層VLSMs的目的。然而對(duì)于間接路徑，目前LSMs和VLSMs的相互作用還尚有爭(zhēng)論。一個(gè)觀點(diǎn)是LSMs和VLSMs的形成是由不同的機(jī)理導(dǎo)致的，理由是Jiménez[25]發(fā)現(xiàn)了LSMs的自維持機(jī)理不需要VLSMs的參與，Hwang[31]通過大渦模擬的概念進(jìn)行濾波測(cè)試，發(fā)現(xiàn)VLSMs可以在沒有LSMs的存在時(shí)產(chǎn)生，Rawat[32]證明VLSMs也是一種自維持結(jié)構(gòu)并且不從LSMs吸收能量；另一種觀點(diǎn)是VLSMs和LSMs相互耦合存在，理由是文獻(xiàn)[30,36]認(rèn)為VLSMs并不是一種新的結(jié)構(gòu)，而是由LSMs發(fā)展到湍流外層疊和而成，Toh[37]發(fā)現(xiàn)LSMs和VLSMs相互作用并彼此促進(jìn)的周期過程，Lee[33]在湍流內(nèi)層發(fā)現(xiàn)了能量反串的過程，即LSMs向VLSMs傳遞能量。

圖3 慣性粒子調(diào)制VLSMs的兩個(gè)路徑，云圖顯示流向截面的瞬時(shí)速度，對(duì)應(yīng)Reτ = 550的明渠流動(dòng)[10]Fig. 3 Schematic of two routes of inertial particle effects on VLSMs through direct impact on turbulent kinetic energy transportation or via indirect upscale energy transfer from LSMs. On the right are streamwise velocity contours in the cross-stream plane of Reτ = 550 open channel flow and the flow region[10]

我們研究[10]發(fā)現(xiàn)VLSMs增加與粒子慣性之間的非單調(diào)變化關(guān)系，即小慣性粒子（St+= 2.42根據(jù)湍流內(nèi)層時(shí)間尺度）和大慣性粒子（Stout= 6.0，8.2，根據(jù)湍流外層時(shí)間尺度）會(huì)增強(qiáng)VLSMs的強(qiáng)度（VLSMs強(qiáng)度通過計(jì)算預(yù)混譜中的雙模態(tài)得到），然而中等慣性粒子對(duì)VLSMs的能譜幾乎沒有影響。為了驗(yàn)證上述提出的雙路徑調(diào)制VLSMs的假設(shè)，首先對(duì)小慣性粒子與大慣性粒子與流動(dòng)的耦合進(jìn)行人為控制，即：

（1）小慣性粒子與湍流內(nèi)層雙方向耦合，但與湍流外層單方向耦合；

（2）小慣性粒子與湍流內(nèi)層單方向耦合，但與湍流外層雙方向耦合；

（3）大慣性粒子與湍流內(nèi)層雙方向耦合，但與湍流外層單方向耦合；

（4）大慣性粒子與湍流內(nèi)層單方向耦合，但與湍流外層雙方向耦合。

測(cè)試結(jié)果為僅有（1）與（4）會(huì)增加VLSMs的強(qiáng)度。此測(cè)試說明，小慣性粒子促進(jìn)VLSMs是通過其在湍流內(nèi)層的反饋力，而大慣性粒子促進(jìn)VLSMs是通過其在湍流外層的反饋力。進(jìn)一步，用譜方法分析了雷諾應(yīng)力輸運(yùn)方程，發(fā)現(xiàn)高慣性粒子對(duì)流體的反饋力作用于雷諾剪切應(yīng)力輸運(yùn)方程，促進(jìn)與VLSMs相同尺度的雷諾剪切應(yīng)力的生成，雷諾剪切應(yīng)力作用于流向湍動(dòng)能輸運(yùn)方程而直接增加VLSMs的強(qiáng)度。物理解釋為：VLSMs與高慣性粒子相互作用，粒子的反饋力促進(jìn)與VLSMs尺度相同的“上揚(yáng)”和“下掃”流動(dòng)，隨后這些超大尺度的“上揚(yáng)”和“下掃”結(jié)構(gòu)從平均流動(dòng)中吸收能量，最終直接促進(jìn)同等尺度VLSMs的生成；低慣性粒子增強(qiáng)湍流自維持過程[9]并增加LSMs的強(qiáng)度，通過能量反向級(jí)串，最終促進(jìn)VLSMs的生成。這種間接調(diào)制機(jī)理，即小粒子通過調(diào)制LSMs并最終調(diào)制VLSMs，支持了上述的第二種觀點(diǎn)[30,33,36-37]。

以上討論了慣性粒子的存在對(duì)LSMs和VLSMs的調(diào)制機(jī)理，除了粒子對(duì)流場(chǎng)的影響以外，研究粒子在LSMs和VLSMs的聚集狀態(tài)也對(duì)大渦模擬粒子流的準(zhǔn)確性及局限性提供借鑒意義。由于計(jì)算能力的限制，目前直接數(shù)值模擬耦合大量的慣性粒子（O（108））只能達(dá)到中等的雷諾數(shù)（Reτ＜1 000），遠(yuǎn)遠(yuǎn)不能達(dá)到實(shí)際情況中的超高雷諾數(shù)Reτ～O（106）。在高雷諾數(shù)的情況下，由于大渦模擬過濾掉小于計(jì)算網(wǎng)格的流動(dòng)結(jié)構(gòu)，導(dǎo)致計(jì)算中得到的粒子本地處的流體速度并不精確（即滑移速度不準(zhǔn)確），Marchioli[38]發(fā)現(xiàn)即使是在低雷諾數(shù)Reτ= 180的情況下，大渦模擬也會(huì)低估粒子近壁面的聚集，同時(shí)粒子的群聚形態(tài)也不同于直接數(shù)值模擬。而在高雷諾數(shù)的情況下，采用大渦模擬往往只能捕捉到VLSMs而過濾掉LSMs，想要準(zhǔn)確地應(yīng)用大渦模擬到高雷諾數(shù)湍流粒子流中[39-40]，有必要衡量這兩種重要流動(dòng)特征結(jié)構(gòu)對(duì)不同慣性的粒子分布的影響。

3 粒子在LSMs及VLSMs中的聚集與

采用了如下的策略來應(yīng)用直接數(shù)值模擬分離LSMs與VLSMs對(duì)粒子的影響：

（1）應(yīng)用截?cái)嘤?jì)算空間的方法分離出LSMs，比較粒子分布在只有LSMs及在同時(shí)有LSMs和VLSMs時(shí)的不同，見圖1；

（2）在計(jì)算過程中過程中，在譜空間分離出LSMs和VLSMs，人為的控制粒子與這兩種結(jié)構(gòu)流場(chǎng)的耦合，比較粒子分布在只有VLSMs及同時(shí)有LSMs和VLSMs時(shí)的不同，見圖4。

圖4 Reτ = 550單相流，y+ = 100水平截面以及四個(gè)邊界截面的瞬時(shí)流向脈動(dòng)速度云圖，“+”代表依據(jù)黏性尺度進(jìn)行無量綱化[11]Fig. 4 Instantaneous contours of streamwise velocity fluctuation on a wall-parallel plane at y+ = 100(and domain boundary walls) in single-phase flow，normalized by uτ[11]

對(duì)于測(cè)試（1），在第1.3節(jié)中已經(jīng)介紹了四種計(jì)算域尺寸對(duì)粒子的影響，重點(diǎn)關(guān)注粒子的一階與二階穩(wěn)態(tài)[8]。本節(jié)用截?cái)嘤?jì)算域作為一種分離LSMs的方法，關(guān)注VLSMs的存在對(duì)粒子分布的影響。首先根據(jù)文獻(xiàn)[41]，選擇計(jì)算域流向尺寸足夠大到包含一個(gè)LSMs并足夠小到排除外層VLSMs對(duì)內(nèi)層LSMs的影響，在Reτ= 550時(shí)，我們選擇截?cái)嘤?jì)算域尺寸為L(zhǎng)x+= 1 375、Lz+= 825，而大尺寸計(jì)算域尺寸為L(zhǎng)x+=10 367、Lz+= 3 456，通過比較二維空間的預(yù)混譜，截?cái)嘤?jì)算域得到的湍流內(nèi)層能量譜與大尺寸計(jì)算域以及文獻(xiàn)[42]的結(jié)果幾乎吻合，但湍流外層的能量譜由于截?cái)嘤虻挠绊懀〕叽缬?jì)算域由于尺寸太小顯然捕捉不到能量譜的雙模態(tài)形式。在加入慣性粒子（低慣性St+= 24.2和高慣性St+= 182）以后：一方面，在湍流內(nèi)層尤其是接近壁面處，粒子的濃度小于大尺寸計(jì)算域接近20%，這歸因于外層VLSMs的存在增強(qiáng)了湍流湍泳（由于脈動(dòng)速度在壁面垂直方向梯度誘導(dǎo)的粒子運(yùn)動(dòng)，粒子趨向于向湍流度弱的方向移動(dòng)），由于粒子的平均濃度在兩種計(jì)算域中保持相同，因此截?cái)嘤虻耐牧魍鈱恿Ｗ訚舛葎t小于大計(jì)算域中的粒子濃度；另一方面，我們通過馮洛諾伊圖來計(jì)算粒子在壁面垂直方向上的群聚形態(tài)[43]，發(fā)現(xiàn)無論對(duì)低或高慣性粒子，截?cái)嘤蚩偸勤呄蛴诘凸懒Ｗ拥娜壕坌?yīng)，尤其是在湍流外層區(qū)域（y+＜ 100），這同樣歸因于外層VLSMs的影響。

對(duì)于測(cè)試（2），在計(jì)算過程中我們?nèi)藶榈目刂屏Ｗ优c流場(chǎng)特定結(jié)構(gòu)的耦合，具體的步驟為：將速度場(chǎng)變換到譜空間，將流場(chǎng)分為L(zhǎng)SMs和VLSMs兩部分，再從譜空間變換到速度場(chǎng)（見圖4），粒子只與LSMs或VLSMs對(duì)應(yīng)的速度場(chǎng)耦合，這種人為控制的粒子-流場(chǎng)耦合在計(jì)算的每一步都進(jìn)行一次，因此計(jì)算時(shí)間較長(zhǎng)。通過以上方式可以得到三種粒子分布，即全流場(chǎng)下、LSMs流場(chǎng)下、VLSMs流場(chǎng)下的粒子分布（見圖5），可以看到對(duì)于小慣性粒子，其粒子群聚主要發(fā)生在湍流內(nèi)層，直觀上LSMs流場(chǎng)下粒子湍流內(nèi)層的群聚看起來與全模擬相似（圖5（a）與（c）比較）。而在湍流外層，由于其響應(yīng)時(shí)間尺度遠(yuǎn)小于特征結(jié)構(gòu)VLSMs的時(shí)間尺度，因此小慣性粒子在湍流外層中趨向于無慣性粒子；對(duì)于高慣性粒子，由于其響應(yīng)時(shí)間尺度與VLSMs的特征尺度相當(dāng)，因此高慣性粒子的群聚發(fā)生于湍流外層，然而即使是直觀上，VLSMs流場(chǎng)下的粒子在湍流外層的群聚形態(tài)也和全模擬相差很大（圖5（b）與（d）比較）。隨后對(duì)粒子聚集（粒子濃度曲線）以及粒子群聚（馮洛諾伊圖計(jì)算）進(jìn)行了比較，發(fā)現(xiàn)在湍流內(nèi)層，LSMs流場(chǎng)會(huì)低估慣性粒子的濃度，在湍流外層，VLSMs流場(chǎng)可以得到和全流場(chǎng)相近的慣性粒子的濃度。對(duì)于粒子群聚現(xiàn)象，LSMs流場(chǎng)結(jié)果更接近于全流場(chǎng)結(jié)果，而VLSMs流場(chǎng)結(jié)果與全流場(chǎng)結(jié)果相差很大。

根據(jù)測(cè)試（1）與測(cè)試（2），發(fā)現(xiàn)LSMs影響粒子的聚集并控制粒子的群聚。在大渦模擬粒子流中，過濾掉LSMs可以估計(jì)湍流外層的粒子濃度，但會(huì)低估湍流內(nèi)層的粒子濃度。大渦模擬捕捉不到正確的粒子群聚，導(dǎo)致粒子與流場(chǎng)的雙向耦合作用到不合適的流場(chǎng)結(jié)構(gòu)中，因此不適宜采用大渦模擬和粒子的雙向耦合研究粒子的反饋?zhàn)饔脤?duì)流場(chǎng)的調(diào)制。

圖5 粒子在三種流場(chǎng)中的分布[11]Fig. 5 Instantaneous snapshots of particle locations (black dots)[11]

4 影響粒子沉積的機(jī)理研究[12-13]

在一些地球物理相關(guān)測(cè)量中，能準(zhǔn)確估算雪、雨、沙塵等的表面沉積速率對(duì)研究降水量、土地沙化、城市污染有重要的實(shí)際應(yīng)用價(jià)值[44-45]。這些粒子的沉積過程，可以看作慣性粒子在重力作用下穿越壁面湍流場(chǎng)到達(dá)壁面這個(gè)簡(jiǎn)化的物理過程。當(dāng)粒子的St+～0時(shí)，粒子可以近似看做具有沉降速度的被動(dòng)標(biāo)量，Rouse[46]推導(dǎo)出了用于預(yù)測(cè)粒子垂直濃度的冪定律，冪指數(shù)與慣性粒子的斯托克斯沉積速度成一定比例，被廣泛應(yīng)用于地球物理的相關(guān)研究中。然而對(duì)于慣性粒子，文獻(xiàn)[46]給出的冪定律則不再適用，Richter[47]嘗試加入粒子對(duì)流速度的一階展開（見文獻(xiàn)[48]）去修改冪定律，通過與直接數(shù)值模擬結(jié)果比較，發(fā)現(xiàn)這樣考慮粒子慣性效應(yīng)的一階修正只使用于St+＜＜1的情況，對(duì)于St+＞1的慣性粒子，目前存在的理論模型仍然很難準(zhǔn)確預(yù)測(cè)粒子的垂直濃度。對(duì)于慣性粒子，文獻(xiàn)[46]結(jié)論不適用的原因是其假設(shè)粒子在湍流場(chǎng)中均勻分布、粒子的沉積速度等于其斯托克斯沉積速度，這樣的假設(shè)對(duì)于St+＞1的慣性粒子不再適用，因?yàn)榧词乖诰鶆蛞恢峦牧髦?，由于慣性粒子趨向于遠(yuǎn)離渦的中心，粒子傾向于被渦旋轉(zhuǎn)向重力方向加速，導(dǎo)致慣性粒子的沉降速度大于其斯托克斯沉積速度[49-50]。而在壁面湍流中，由于多種湍流尺度的存在，不同慣性的粒子會(huì)偏向聚集于湍流場(chǎng)的不同結(jié)構(gòu)中，影響粒子加速的因素變得更加復(fù)雜。

在文獻(xiàn)[12-13]中，我們首先通過推導(dǎo)相空間下的概率密度分布的主方程（phase-space master PDF equation）得到了控制粒子穩(wěn)態(tài)的輸運(yùn)方程，進(jìn)而推導(dǎo)出控制粒子沉降速度的方程（數(shù)學(xué)推導(dǎo)過程請(qǐng)參考文獻(xiàn)[12,51]）：

其中vp(t) 為粒子的沉降速度，y為壁面垂直方向，ρ為邊〈緣概率密度函〉 數(shù)用來描述粒子的空間分布，S=(vp(t)-〈vp(t)〉)2為粒子脈動(dòng)速度的協(xié)方差，up(t)為粒子處的流體速度。通過對(duì)此控制方程的分析，發(fā)現(xiàn)粒子的沉降速度是由如下五個(gè)因素來控制，對(duì)應(yīng)方程右側(cè)由左至右依次為：粒子加速項(xiàng)，粒子擴(kuò)散項(xiàng)，湍流湍泳項(xiàng)，粒子本地流場(chǎng)速度項(xiàng)以及重力加速項(xiàng)。對(duì)上述五個(gè)因素的物理描述為：粒子加速項(xiàng)是由粒子在壁面垂直方向的速度梯度導(dǎo)致的，粒子擴(kuò)散項(xiàng)是由粒子在壁面垂直方向的濃度梯度導(dǎo)致的，粒子本地流場(chǎng)速度項(xiàng)是由粒子在湍流中的偏向聚集現(xiàn)象導(dǎo)致的，重力加速項(xiàng)是由外力導(dǎo)致的。隨后通過直接數(shù)值模擬計(jì)算了上述五個(gè)因素在零粒子通量以及有限大小粒子通量?jī)煞N情況下對(duì)粒子沉降速度的貢獻(xiàn)。進(jìn)一步的，我們討論了擴(kuò)展文獻(xiàn)[46]中的模型從St+～0粒子到St+＞ 1粒子的可能性，發(fā)現(xiàn)模型中的不封閉項(xiàng)包括由于粒子聚集導(dǎo)致的漂移系數(shù)以及粒子壁面垂直方向的湍動(dòng)能。我們目前仍然尚無法提出滿意的模型用來預(yù)測(cè)St+＞ 1慣性粒子的壁面垂直方向的濃度曲線，這將是未來要開展的工作之一。

5 結(jié) 論

本文應(yīng)用直接數(shù)值模擬耦合拉格朗日點(diǎn)粒子這種相對(duì)簡(jiǎn)單的數(shù)值模型，對(duì)慣性粒子和壁面湍流這個(gè)復(fù)雜且具有重要實(shí)際應(yīng)用的問題開展了一些相關(guān)研究。重點(diǎn)關(guān)注的物理問題是粒子對(duì)壁面湍流的調(diào)制，以及粒子在兩種特征湍流結(jié)構(gòu)中的分布以及粒子的沉積速度。在未來的研究中，將繼續(xù)嘗試建模預(yù)測(cè)粒子垂向濃度曲線，熱對(duì)流效應(yīng)對(duì)壁面湍流及粒子分布的影響，粒子在復(fù)雜流場(chǎng)結(jié)構(gòu)中的運(yùn)動(dòng)（例如旋轉(zhuǎn)湍流），帶電/磁粒子在壁面湍流中的運(yùn)動(dòng)，以及應(yīng)用大渦模擬技術(shù)研究復(fù)雜地形下粒子分布等相關(guān)研究工作。

致謝：本文涉及到的工作是作者在美國(guó)圣母大學(xué)的研究?jī)?nèi)容，全部已公開發(fā)表到經(jīng)過同行評(píng)審的專業(yè)期刊或尚未經(jīng)過同行評(píng)審的預(yù)印版。本工作得益于美國(guó)陸軍研究實(shí)驗(yàn)室以及海軍研究辦公室的項(xiàng)目支持，計(jì)算資源來自于高性能計(jì)算現(xiàn)代化項(xiàng)目（HPCMP）以及圣母大學(xué)研究計(jì)算中心。相關(guān)的研究工作得益于幾位合作者的貢獻(xiàn)：美國(guó)圣母大學(xué)的Hyungwon John Park博士以及David H. Richter副教授，美國(guó)杜克大學(xué)的Andrew Bragg助理教授，美國(guó)明尼蘇達(dá)大學(xué)的Kee Onn Fong博士以及Filippo Coletti副教授（現(xiàn)工作于蘇黎世聯(lián)邦理工大學(xué)），美國(guó)密歇根大學(xué)的Jesse Capecelatro助理教授。感謝幾位優(yōu)秀學(xué)者訪問本人所在的美國(guó)圣母大學(xué)課題組時(shí)所提出的寶貴建議：美國(guó)加州大學(xué)洛杉磯分校的Marcelo Chamecki教授以及美國(guó)華盛頓大學(xué)的James J.Riley教授。感謝已發(fā)表論文和本論文匿名審稿人提出的寶貴建議。最后感謝相關(guān)專家的推薦及《空氣動(dòng)力學(xué)學(xué)報(bào)》編輯部的邀請(qǐng)。