李 蒙，徐振山，陳永平

(河海大學(xué)港口海岸與近海工程學(xué)院，江蘇 南京 210098)

廢水(如達(dá)標(biāo)廢水、濃鹽水、溫排水)向河流和海洋的排放是環(huán)境水力學(xué)和排放工程關(guān)注的重點(diǎn)問(wèn)題。陸源廢水通常以射流的形式排放到受納水體中，排放過(guò)程可劃分為排放近區(qū)的初始稀釋過(guò)程和排放遠(yuǎn)區(qū)的后續(xù)稀釋過(guò)程。Zhao等[1]認(rèn)為廢水在排放近區(qū)的初始稀釋度將決定排放遠(yuǎn)區(qū)污染物的影響范圍。因此，研究廢水在排放近區(qū)初始稀釋度具有現(xiàn)實(shí)意義。

定常排放的非振蕩射流是目前廢水排放的主要形式。Kelso等[2]發(fā)現(xiàn)，與靜水環(huán)境中射流相比，橫流環(huán)境中的非振蕩射流軸線流速衰減明顯加快，拓展寬度大大增加。Sau等[3-4]學(xué)者總結(jié)出非振蕩射流4種典型的渦動(dòng)結(jié)構(gòu)為剪切層渦、馬蹄形渦、反向?qū)ΨQ(chēng)漩渦(CVP)和尾渦，其中CVP結(jié)構(gòu)是橫流環(huán)境中非振蕩射流最為重要的特征之一。Plesniak等[5-6]發(fā)現(xiàn)射流-橫流速度比是影響非振蕩射流運(yùn)動(dòng)和擴(kuò)散特性的主要因素。

在排放近區(qū)，除受納水體的動(dòng)力外，射流本身的參數(shù)如射流水體速度、密度、排放角度也會(huì)對(duì)初始稀釋過(guò)程產(chǎn)生影響。正弦振蕩射流的研究近幾年來(lái)有了一定的進(jìn)展。Hsu等[7-8]采用粒子圖像測(cè)速技術(shù)(PIV)對(duì)橫流環(huán)境中正弦振蕩射流的流場(chǎng)進(jìn)行測(cè)量，將一個(gè)周期內(nèi)射流的運(yùn)動(dòng)過(guò)程劃分為4個(gè)典型模式：下沖模式、恒定流主導(dǎo)模式、射流主導(dǎo)模式和過(guò)渡模式。Marcum等[9-10]通過(guò)室內(nèi)物理試驗(yàn)分析了射流振蕩幅度和頻率對(duì)橫流-射流相互作用的影響，認(rèn)為可通過(guò)調(diào)節(jié)射流振幅和頻率來(lái)增大射流與周?chē)黧w的摻混和卷吸。波流環(huán)境中的非振蕩射流與橫流環(huán)境中的振蕩射流有著類(lèi)似之處。徐振山等[11-12]建立了波流共同作用下非振蕩射流的三維數(shù)學(xué)模型，并系統(tǒng)分析了非振蕩射流的三維流場(chǎng)結(jié)構(gòu)，發(fā)現(xiàn)表征射流最主要特征的CVP結(jié)構(gòu)發(fā)生了較大改變，非振蕩射流與周?chē)w的摻混大大增加，其擴(kuò)展寬度也明顯增加。張玉玲等[13]與丁宏偉等[14]分別進(jìn)行了波流環(huán)境中多孔射流與雙孔浮射流的研究，其相關(guān)結(jié)果可以為振蕩射流的研究提供參考。

截至目前，關(guān)于橫流環(huán)境中正弦振蕩射流的研究還局限于二維層面，三維的運(yùn)動(dòng)特性和稀釋規(guī)律仍需進(jìn)一步研究。本文采用大渦數(shù)學(xué)模型對(duì)橫流環(huán)境下的正弦振蕩射流進(jìn)行模擬，研究振蕩射流在橫流環(huán)境下的運(yùn)動(dòng)特性和稀釋規(guī)律，探討振蕩參數(shù)對(duì)射流的影響。

1 數(shù)學(xué)模型簡(jiǎn)介

1.1 LES控制方程和數(shù)值算法

大渦數(shù)值模擬(large-eddy-simulation，LES)的控制方程為空間平均的N-S方程和標(biāo)量輸運(yùn)方程：

(1)

(2)

(3)

式中：xi(i=1，2，3) 為笛卡爾坐標(biāo)系統(tǒng)中的x、y、z；ui(i=1，2，3)為xi方向的瞬時(shí)速度分量u、v、w；ρ為密度；p為瞬時(shí)靜水壓強(qiáng)；v為運(yùn)動(dòng)黏滯系數(shù)；gi為重力加速度；τij,SGS為亞格子剪應(yīng)力；c為濃度標(biāo)量；Γ為分子擴(kuò)散系數(shù)；qi,SGS為標(biāo)量亞格子紊流通量；Pr為分子Prandtl數(shù)。

τij,SGS可分為各向異性的偏斜部分τij和各向同性的對(duì)角部分τkk,SGS，τkk,SGS可與壓力項(xiàng)合并求解。應(yīng)力張量τij采用Smagorinsky模型[15]模擬表示為

(4)

(5)

式中：Sij為對(duì)稱(chēng)變形速度梯度張量；νt為渦黏性系數(shù)；Cs為Smagorinsky常數(shù)，一般取值為0.12～0.28，根據(jù)張玉玲[16]模擬率定可得Cs=0.175時(shí)更符合實(shí)際情況；Δ為濾波尺度，表示為Δ=(Δx1Δx2Δx3)1/3，Δx1=Δx、Δx2=Δy、Δx3=Δz分別為x、y和z方向上的空間步長(zhǎng)。

采用有限差分法中的分裂算子法求解控制方程，計(jì)算過(guò)程主要分為3步：第一步求解動(dòng)量方程中的對(duì)流項(xiàng)，采用二次向后特征線法和 Lax-Wendroff格式平均法對(duì)對(duì)流項(xiàng)進(jìn)行離散求解以提高計(jì)算精度；第二步采用時(shí)間前差、空間中心差分格式離散求解擴(kuò)散項(xiàng)；第三步采用投影法對(duì)壓力傳播項(xiàng)中的壓力和速度進(jìn)行耦合求解，并更新自由表面。

1.2 模型邊界條件

數(shù)值計(jì)算采用合成渦法(synthetic-eddy-method，SEM)生成振蕩射流邊界條件，并將得到的振蕩射流植入數(shù)值水槽中進(jìn)行模擬。底部采用可滑移邊界預(yù)測(cè)第一個(gè)網(wǎng)格節(jié)點(diǎn)的速度梯度，采用壁面函數(shù)法計(jì)算底部切應(yīng)力。側(cè)向邊界采用不可入邊界條件。

2 模型驗(yàn)證與模擬工況

2.1 橫流環(huán)境中射流驗(yàn)證

采用一組橫流環(huán)境中非振蕩射流的數(shù)據(jù)對(duì)模型速度與濃度結(jié)果進(jìn)行驗(yàn)證，具體參數(shù)如下：試驗(yàn)水深h=0.5 m，非振蕩射流平均速度w0=1.22 m/s，射流雷諾數(shù)Re=12 200，水流流速u(mài)0=0.133 m/s。

圖1與圖2分別為不同下游斷面上射流時(shí)均垂向速度、濃度沿水深分布數(shù)值模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)對(duì)比。在射流主體區(qū)域數(shù)值模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)吻合較好，試驗(yàn)數(shù)據(jù)與數(shù)值模擬結(jié)果誤差在可接受范圍內(nèi)。針對(duì)該模型的更多驗(yàn)證參見(jiàn)文獻(xiàn)[17]。

圖1 不同下游斷面上射流時(shí)均垂向速度數(shù)值模擬結(jié)果與試驗(yàn)數(shù)據(jù)對(duì)比Fig.1 Comparison of numerical simulation results and experimental data of jet time-averaged vertical velocity on different downstream cross sections

圖2 不同下游斷面上射流時(shí)均濃度數(shù)值模擬結(jié)果與試驗(yàn)數(shù)據(jù)對(duì)比Fig.2 Comparison of numerical simulation results and experimental data of jet time-averaged concentration on different downstream cross sections

2.2 模擬工況設(shè)置

振蕩射流的影響因素包括受納水體參數(shù)，為橫流流速u(mài)0；射流本身的參數(shù)包括射流平均速度w0，振蕩振幅wA，以及周期T(或頻率f=1/T)。用斯特勞哈爾數(shù)St表示射流周期T的影響(式(6))。定義在一個(gè)周期內(nèi)，t=1/4T為最大流速時(shí)刻，t=1/2T為速度逐漸減小的平均流速時(shí)刻，t=3/4T為最小流速時(shí)刻，t=T為速度逐漸增大的平均流速時(shí)刻。影響橫流環(huán)境中振蕩射流的因素有射流-橫流速度比Rjc(式(7))、振幅-射流速度比Raj(式(8))以及斯特勞哈爾數(shù)St(式(6)))：

(6)

(7)

(8)

式中：d為射流口孔徑，為0.01 m。

各組次的橫流以及振蕩射流的計(jì)算參數(shù)如表1所示。

表1 橫流環(huán)境中振蕩射流計(jì)算組次及參數(shù)

3 模擬結(jié)果與討論

3.1 射流對(duì)稱(chēng)縱剖面特征

從相位平均流場(chǎng)與時(shí)均流場(chǎng)兩個(gè)方面描述射流的對(duì)稱(chēng)縱剖面特征。

圖3為NJ04組次4個(gè)典型相位下對(duì)稱(chēng)縱剖面(y/d=0)上射流相位平均流速場(chǎng)與濃度場(chǎng)。從圖3可看出，在橫流作用下，振蕩射流呈現(xiàn)出明顯的間歇性污染物云團(tuán)現(xiàn)象，與波流環(huán)境中的非振蕩射流[13]類(lèi)似。間歇性污染物云團(tuán)分別標(biāo)記為A1～A3、B1～B3、C1～C3、D1～D3。A1-B1-C1-D1-A2可以看作在一個(gè)振蕩周期內(nèi)間歇性污染物云團(tuán)的形成和發(fā)展過(guò)程。1/4T時(shí)刻沖擊高度最小，彎曲程度最大，在這一階段，射流柱隨時(shí)間的推移而向上傾斜，污染物云團(tuán)逐漸發(fā)展擴(kuò)大；1/2T時(shí)刻射流速度減小，使得射流柱向下傾斜，當(dāng)污染物云團(tuán)向下游發(fā)展時(shí)，由于剪切層渦與橫流之間發(fā)生卷吸效應(yīng)，迎流側(cè)剪切層渦擴(kuò)大并變得不連貫；3/4T時(shí)刻沖擊高度最大，彎曲程度最小，射流柱被橫流沖擊產(chǎn)生明顯偏轉(zhuǎn)；T時(shí)刻射流速度逐漸增大，射流柱逐漸向上傾斜。這一過(guò)程在每個(gè)循環(huán)中持續(xù)進(jìn)行，與Shi等[18]的研究一致，振蕩射流的存在使得污染物云團(tuán)周期性地向下游傳播。根據(jù)Hsu等[8]的研究，1/4T時(shí)刻為射流主導(dǎo)模式，1/2T時(shí)刻為過(guò)渡模式，3/4T時(shí)刻為下沖模式，T時(shí)刻為橫流主導(dǎo)模式。

圖3 NJ04組次4個(gè)典型相位對(duì)稱(chēng)縱剖面(y/d=0)上射流相位平均流速場(chǎng)和濃度場(chǎng)Fig.3 Phase-averaged flow velocity field and concentration field of jet on symmetrical longitudinal section(y/d=0) under four typical phases of group NJ04

圖4是NJ00、NJ01、NJ02、NJ03組次對(duì)稱(chēng)縱剖面上射流的時(shí)均流速場(chǎng)與濃度場(chǎng)對(duì)比。相較于橫流環(huán)境中的非振蕩射流，振蕩射流的間歇性污染物云團(tuán)沖擊高度更高，彎曲程度更大，表明振蕩射流的擴(kuò)散速度更快，分布寬度更大。隨著射流振蕩振幅和周期的增大，間歇性污染物云團(tuán)沖擊高度更高，彎曲程度減小。由于振蕩射流增強(qiáng)了流場(chǎng)的非定常性，當(dāng)射流-橫流速度比相同時(shí)，振蕩射流與橫流的相互作用使得橫流環(huán)境中振蕩射流的時(shí)均速度矢量較非振蕩射流向下彎曲，整個(gè)射流體更貼近底床，表明在射流周期性振蕩的影響下，射流前行過(guò)程中的有效動(dòng)量較無(wú)振蕩時(shí)衰減得更快。另外，振蕩作用(振幅與周期)的加強(qiáng)使射流軌跡擺動(dòng)幅度增大，流線發(fā)生更加明顯的彎曲，垂向平均流速衰減得更快。

圖4 4組次對(duì)稱(chēng)縱剖面(y/d=0)上射流時(shí)均流場(chǎng)和濃度場(chǎng)Fig.4 Time-averaged flow field and concentration field of jet on symmetrical longitudinal section(y/d=0) of four groups

圖5為橫流環(huán)境中非振蕩射流與振蕩射流的時(shí)均流速分布，兩者的垂向流速基本呈現(xiàn)高斯分布。由圖5可見(jiàn)，污染物在向下游傳播的過(guò)程中，流速分布曲線的高斯分布更加明顯，分布寬度變寬，最大流速值減小且其垂向位置逐漸變高。非振蕩射流的最大流速均大于振蕩射流的最大流速，而振蕩射流的分布寬度明顯大于非振蕩射流，可見(jiàn)振蕩作用引發(fā)的擴(kuò)散效應(yīng)顯著。

圖5 橫流環(huán)境中非振蕩射流與振蕩射流下游斷面時(shí)均流速分布Fig.5 Time-averaged velocity distribution on downstream sections of non-oscillating jet and oscillating jet in cross flow environment

3.2 射流下游斷面特征

圖6分別為橫流環(huán)境中非振蕩射流與振蕩射流(組次NJ00和組次NJ03)各下游斷面上時(shí)均濃度場(chǎng)。從圖6明顯可以看出，振蕩射流的存在使得射流在垂向上的寬度明顯增加。

圖6 橫流環(huán)境中非振蕩射流與振蕩射流各下游斷面上時(shí)均濃度場(chǎng)Fig.6 Time-averaged concentration fields of non-oscillating jet and oscillating jet at each downstream section in cross flow environment

橫流環(huán)境中，由于CVP結(jié)構(gòu)的存在，非振蕩射流下游斷面時(shí)均濃度有兩個(gè)極大值；振蕩射流除上述兩個(gè)濃度極大值外，還存在一個(gè)由間歇性污染物云團(tuán)引起的極大值。通過(guò)對(duì)比發(fā)現(xiàn)，振蕩作用不僅使斷面上時(shí)均濃度的最大值及其垂向位置發(fā)生變化，也使一定濃度包絡(luò)線的面積發(fā)生了改變。為了對(duì)比各組次的稀釋規(guī)律，采用斷面濃度最大值Cm(或斷面最小稀釋度Sc)、斷面濃度最大值Cm的垂向位置Zm、由0.25Cm等值線決定的射流可視范圍(即寬度Rh和高度Rv)3個(gè)指標(biāo)描述橫流環(huán)境中射流下游斷面上的濃度特征，如圖7所示。最小稀釋度公式為

(9)

式中C0為初始濃度。

圖7 描述橫流環(huán)境中非振蕩射流和振蕩射流下游斷面上濃度特征的指標(biāo)Fig.7 Index to describe the concentration characteristics of non-oscillating jet and oscillating jet for downstream sections in cross flow environment

圖8為NJ00、NJ01、NJ02、NJ03組次的稀釋參數(shù)對(duì)比。與橫流環(huán)境中非振蕩射流相比，振蕩射流的濃度最大值較小，即斷面最小稀釋度較大，最小稀釋度垂向位置更靠近底床，斷面可視范圍面積顯著增加。由此可見(jiàn)，振蕩作用使得射流體與周?chē)w的相互作用加強(qiáng)，增強(qiáng)了射流的稀釋能力；并且振蕩作用增加了射流沖擊高度，使得射流的影響范圍變大，有利于射流與周?chē)w的混合。以斷面最小稀釋度為例，振蕩射流(NJ03)與非振蕩射流(NJ00)斷面最小稀釋度為35的水平位置分別為x/d=21和x/d=29，可見(jiàn)振蕩作用使射流得到了更快的稀釋。NJ01與NJ03的對(duì)比組也顯示了與上述現(xiàn)象相似的規(guī)律，說(shuō)明隨著振幅的增大，射流的稀釋能力也逐漸增強(qiáng)。

圖8 各組次稀釋度特性對(duì)比Fig.8 Comparison of dilution degree characteristics of each group

根據(jù)NJ01與NJ02對(duì)比斯特勞哈爾數(shù)(周期)的影響，從圖8可以看到斯特勞哈爾數(shù)對(duì)最小稀釋度影響較小，在距射流出口水平距離較小的位置(x/d<26)，兩組射流的斷面最小稀釋度基本一樣，之后發(fā)生一定變化——斯特勞哈爾數(shù)越小(周期越大)，最小稀釋度越大。最小稀釋度垂向位置與斷面可視范圍面積隨著斯特勞哈爾數(shù)減小分別呈現(xiàn)出降低與增大的現(xiàn)象。總體來(lái)說(shuō)，斯特勞哈爾數(shù)越小(周期越大)，振蕩射流稀釋效果越好。

4 結(jié) 論

a.在橫流環(huán)境下振蕩射流比非振蕩射流下游的CVP結(jié)構(gòu)主流線位置更靠下，間歇性污染物云團(tuán)沖擊高度更高，水體間摻混更加充分。當(dāng)前組次下，振蕩射流將斷面最小稀釋度35的位置從x/d=29處前移到x/d=21處，稀釋能力大大增強(qiáng)。

b.當(dāng)前組次下，振蕩幅度增加1倍，x/d=25處斷面可視范圍面積增加了近10%，稀釋能力顯著增強(qiáng)。表明在橫流環(huán)境中，振蕩射流的振蕩作用越大，即振幅與周期越大，其沖擊高度越高，彎曲程度越小，最小稀釋度越大，斷面可視范圍面積越大。

總體而言，橫流環(huán)境中的振蕩射流有利于污染物稀釋?zhuān)罄m(xù)在排污擴(kuò)散器設(shè)計(jì)時(shí)可采用振蕩射流模式。