異重流潛入運(yùn)動(dòng)的剖面二維數(shù)值模擬(1)

1、    異重流潛入運(yùn)動(dòng)的剖面二維數(shù)值模擬(1)    本文采用變密度流基本方程和混合有限分析法，求解了突擴(kuò)邊界下的異重流潛入運(yùn)動(dòng)。通過(guò)對(duì)數(shù)值計(jì)算結(jié)果分析，討論了進(jìn)口密度弗汝德數(shù)Fre對(duì)潛入運(yùn)動(dòng)的影響，加深了對(duì)異重流潛入規(guī)律的了解。 關(guān)鍵詞：異重流 潛入運(yùn)動(dòng) 數(shù)值計(jì)算 混合有限分析法 1 前言異重流是自然界中常見(jiàn)的一種流動(dòng)現(xiàn)象，如水庫(kù)中的渾水異重流，河口地區(qū)出現(xiàn)的鹽水楔，工業(yè)中遇到的溫差異重流，底部突然下沉的小型沉沙池中的泥沙異重流等，在水庫(kù)、船閘引航道、河港、河口、沉沙池等實(shí)際工程中，渾水異重流的淤積問(wèn)題往往是必須

2、考慮的問(wèn)題1。由于異重流潛入發(fā)生處的邊界條件復(fù)雜，具有與一般明渠流不同的水流泥沙運(yùn)動(dòng)規(guī)律，且潛入處的水力泥沙因素與異重流的發(fā)生、發(fā)展及運(yùn)行狀況密切相關(guān)2，歷來(lái)廣為研究者們所關(guān)注。用一維的方法，再加上一些經(jīng)驗(yàn)關(guān)系解異重流運(yùn)動(dòng)，無(wú)疑具有其實(shí)際意義，但異重流運(yùn)動(dòng)具有較強(qiáng)的三維性，不是一維流動(dòng)理論基礎(chǔ)所能較好解決的，并且在用一維方法求解異重流運(yùn)動(dòng)時(shí)，常常碰到諸如上下水層交界面位置及交界面上水體之間的動(dòng)量和質(zhì)量交換系數(shù)等難以確定的困難。近年來(lái)，逐步采用數(shù)值模擬技術(shù)對(duì)異重流潛入運(yùn)動(dòng)進(jìn)行研究，如Yoon T.H。采用流函數(shù)渦量法對(duì)異重流的潛入運(yùn)動(dòng)進(jìn)行了模擬和研究3；方春明、韓其為用立面二維對(duì)異重流潛入流動(dòng)進(jìn)

3、行數(shù)值模擬，探討了影響異重流潛入條件的因素4。但由于問(wèn)題本身的復(fù)雜性，許多問(wèn)題仍有待進(jìn)一步研究。    考慮到異重流潛入運(yùn)動(dòng)的過(guò)程中，不僅具有初始動(dòng)量，而且還受到因密度差引起的浮力作用，與靜止環(huán)境中的浮力射流運(yùn)動(dòng)有相似之處，不同之處在于異重流潛入發(fā)生處的邊界條件復(fù)雜，因此，本文將異重流的潛入運(yùn)動(dòng)視為一種具有復(fù)雜邊界的浮力射流運(yùn)動(dòng)，從Navier Stokes方程出發(fā)，引入滿足布辛涅斯克假定的一組變密度方程，并采用混合有限分析法對(duì)該運(yùn)動(dòng)進(jìn)行一系列的數(shù)值模擬，通過(guò)對(duì)數(shù)值計(jì)算結(jié)果進(jìn)行討論和分析，進(jìn)一步加深對(duì)該流動(dòng)的認(rèn)識(shí)。 2 數(shù)學(xué)模型  &

4、#160;      圖1 二維異重流潛入運(yùn)動(dòng)的示意圖Schematic diagram of 2 Ddensity current plunging flow圖1所示：密度為(0 )的渾水從水平河道流入密度為0的清水矩形水庫(kù)中，由于密度分布不均勻，流場(chǎng)中就出現(xiàn)了浮力，在幾何形狀突然擴(kuò)大的流道中，形成了有浮力回流的流動(dòng)。圖中x、z分別表示水平和豎直方向，進(jìn)口渾水流速u0滿足拋物線分布。2.1 控制方程該模型遵循以下假定：忽略泥沙與水流之間的速度差異；根據(jù)布辛涅斯克近似，僅重力項(xiàng)計(jì)及密度變化；認(rèn)為泥沙在水平和豎直方向上擴(kuò)散系數(shù)相同；僅對(duì)異重流潛入

5、段進(jìn)行模擬，此時(shí)浮力作用較大，且計(jì)算區(qū)間不太長(zhǎng)，因此在這里只考慮異重流的對(duì)流和擴(kuò)散，不考慮其沿程沉降4。為使計(jì)算和討論具有一般性，以水庫(kù)初始水深H、進(jìn)口平均流速u0為特征尺度，引入無(wú)量綱量X=x/H,Z=z/H,U=u/u0,W=w/u0,T=tu0/H,P*=p/u20則控制方程的無(wú)量綱形式為    連續(xù)方程            (1)    運(yùn)動(dòng)方程   

6、60;        (2)        (3)    濃度方程            (4)    式中 =-0/0為渾水的有效密度，Re=U0H/v為雷諾數(shù)，F(xiàn)r=U0/為弗汝德數(shù)，Sc=v/D為斯密特?cái)?shù)，v、D分別為流體的運(yùn)動(dòng)粘性系

7、數(shù)和泥沙擴(kuò)散系數(shù)。2.2 初始和邊界條件2.2.1 初始條件初始時(shí)刻，水庫(kù)水體靜止，進(jìn)口斷面的渾水流速滿足拋物線分布，且認(rèn)為此斷面的密度分布均勻；初始?jí)毫Π挫o水壓力給出：P(X,Z)=1-Z。2.2.2 邊界條件研究的流動(dòng)區(qū)域共有進(jìn)口、出口、自由表面、上游固壁及底部五個(gè)不同的邊界，如圖1所示。進(jìn)口：渾水流速滿足拋物線分布，斷面密度均勻分布，即U=7.5(6-5Z)(Z-H-H0/H),W=0,=/0;    出口：認(rèn)為流動(dòng)已充分發(fā)展，則;壁面：給定無(wú)滑移邊界條件，即U=0,W=0,(上游)或(底部)；自由表面：對(duì)自由面作剛蓋假定，并采用對(duì)稱平面條件，即。

8、            摘本文采用變密度流基本方程和混合有限分析法，求解了突擴(kuò)邊界下的異重流潛入運(yùn)動(dòng)。通過(guò)對(duì)數(shù)值計(jì)算結(jié)果分         本篇論文是由3COME文檔頻道的網(wǎng)友為您在網(wǎng)絡(luò)上收集整理餅投稿至本站的，論文版權(quán)屬原作者，請(qǐng)不用于商業(yè)用途或者抄襲，僅供參考學(xué)習(xí)之用，否者后果自負(fù)，如果此文侵犯您的合法權(quán)益，請(qǐng)聯(lián)系我們。    3 計(jì)算方法3.1 求解區(qū)

9、域及網(wǎng)格劃分該流動(dòng)問(wèn)題是開(kāi)邊界問(wèn)題中的一種，當(dāng)出口采用充分發(fā)展假定時(shí)，求解區(qū)域應(yīng)保持足夠的長(zhǎng)度，以便在我們感興趣的區(qū)域準(zhǔn)確地模擬出流體的運(yùn)動(dòng)規(guī)律。依文獻(xiàn)5，本文計(jì)算區(qū)域長(zhǎng)度取為L(zhǎng)=20。因水庫(kù)進(jìn)口附近的流動(dòng)是主的，遠(yuǎn)離進(jìn)口的流動(dòng)是漸變的，故在數(shù)值計(jì)算中，X方向上采用非均勻網(wǎng)格，進(jìn)口附近的網(wǎng)格密，越遠(yuǎn)離進(jìn)口，網(wǎng)格越疏；在Z方向上，除了0、1附近的網(wǎng)格密些（與邊界條件有關(guān)），其余部分采用均勻網(wǎng)格。依控制方程及邊界區(qū)域等特點(diǎn)，采用原始變量法求解為便，故采用交錯(cuò)網(wǎng)格來(lái)離散并求解整個(gè)流場(chǎng)和密度場(chǎng)6。計(jì)算時(shí)，動(dòng)量方程和濃度方程的控制精度為10-7，連續(xù)方程的控制精度為10-4。3.2 計(jì)算格式采用混合有限

10、分析法的五點(diǎn)格式進(jìn)行離散。在局部單元格上，將非線性方程（2）（4）線性化，寫(xiě)成下列一般形式        (5)則各方程的系數(shù)可歸納為表1的形式。方程(5)的混合有限分析解為        (6)表1 控制方程系數(shù)Coefficients of governing equations              

11、;   E    2A    2B    S            U    Re    ReU    ReW        &#

12、160;       W    Re    ReU    ReW    -Re1/Fr2                ScRe    ScReU  

13、  ScReW    0            式中 Ci-1,j、Ci 1,j、Ci,j-1、Ci,j 1和Cp分別為混合有限分析法系數(shù)，具體表達(dá)見(jiàn)文獻(xiàn)7。其中，(S-Et)n為第n次的迭代值，Sn用中心差分計(jì)算，nt用向前差分表示，即()np=ij n 1-nij /t，將(i-1,j)、(i 1,j)、(i,j-1)、(i,j 1)依此用1、2、3、4表示，則式（6）可寫(xiě)為   &

14、#160;    (7)對(duì)整個(gè)計(jì)算區(qū)域，由式（2）（4）形成代數(shù)方程組，它們都是對(duì)角占優(yōu)的，可用追趕法求解。4 計(jì)算結(jié)果和分析4.1 Fre對(duì)潛入形態(tài)的影響異重流在潛入水庫(kù)的過(guò)程中，不僅具有初始動(dòng)量，而且還受到因密度差而引起的浮力作用，流動(dòng)是非定常和變密度的，此時(shí)影響流動(dòng)的一個(gè)重的特征數(shù)是進(jìn)口密度弗汝德數(shù)Fre，具體表達(dá)形式為        (8)式中u0為進(jìn)口處渾水的特征流速，L為特征長(zhǎng)度，在本文中，u0等于進(jìn)口處渾水的平均流速，L等于進(jìn)口水深H0，為進(jìn)口處渾水的密度，&#

15、160;           摘本文采用變密度流基本方程和混合有限分析法，求解了突擴(kuò)邊界下的異重流潛入運(yùn)動(dòng)。通過(guò)對(duì)數(shù)值計(jì)算結(jié)果分         本篇論文是由3COME文檔頻道的網(wǎng)友為您在網(wǎng)絡(luò)上收集整理餅投稿至本站的，論文版權(quán)屬原作者，請(qǐng)不用于商業(yè)用途或者抄襲，僅供參考學(xué)習(xí)之用，否者后果自負(fù)，如果此文侵犯您的合法權(quán)益，請(qǐng)聯(lián)系我們。0為水庫(kù)的清水密度。Fre反映了進(jìn)口水流慣性力作用與浮力作用之比：當(dāng)其它條件不

16、變時(shí)，F(xiàn)re越小，說(shuō)明渾水慣性力越小，浮力作用越大；反之亦然。為研究異重流潛入形態(tài)隨Fre的變化規(guī)律，本文在進(jìn)口單寬流量一定的情況下，通過(guò)改變進(jìn)口處渾水的密度，分別計(jì)算了Fre=0.25，0.501，1.00，1.58時(shí)異重流各特征量的變化情況，并繪出了進(jìn)口附近區(qū)域在T=0.8（T為無(wú)量綱化時(shí)間，T=tu0/h）時(shí)刻相應(yīng)的有效密度等值線圖（/0）和速度矢量圖（圖2和圖3），由圖可看出，F(xiàn)re 對(duì)異重流潛入流動(dòng)特性的影響是較大的。    圖2是不同F(xiàn)re的相對(duì)有效密度等值線圖：當(dāng)Fre=0.25時(shí)，負(fù)浮力占主導(dǎo)地位，渾水入庫(kù)后不久便下沉，在進(jìn)口附近形成一股

17、潛流；當(dāng)Fre=0.501，1.00，1.58時(shí)，入庫(kù)渾水的初始動(dòng)量作用逐漸增大，負(fù)浮力失去其主導(dǎo)作用，渾水繼續(xù)向前水平運(yùn)動(dòng)，直到水庫(kù)清水對(duì)其阻力與維持其運(yùn)動(dòng)的力平衡為止，此時(shí)因負(fù)浮力作用渾水開(kāi)始潛入，該位置通常也稱為潛入點(diǎn)。從該圖還可以看出，F(xiàn)re越大，渾水與水庫(kù)水體的混合越充分，在水庫(kù)中的擴(kuò)散程度越大，潛入點(diǎn)距進(jìn)口的距離越長(zhǎng)，越不容易形成異重流；反之，則越容易形成異重流。圖3是相應(yīng)的速度矢量隨Fre的變化圖。該圖反映了在進(jìn)口下游一定區(qū)域內(nèi)會(huì)形成較大的回流旋渦區(qū)，隨著Fre的改變，回流旋渦區(qū)的形狀、大小和中心位置相應(yīng)發(fā)生改變。      &#

18、160;  圖2 異重流潛入時(shí)的密度場(chǎng)(/n)隨Fre的變化(T=0.8)Density field (/n) with the variation of Fre(T=0。8)in the case of density current plunging當(dāng)Fre較小時(shí)，渾水在進(jìn)口附近區(qū)域很快吸附近岸，在水庫(kù)中形成一個(gè)逆時(shí)針?lè)较虻男郎u，且該回流旋渦區(qū)范圍較大，密度場(chǎng)的范圍相應(yīng)地也較大；隨著Fre的增大，在進(jìn)口附近將出現(xiàn)一順時(shí)針?lè)较虻男郎u，同時(shí)逆時(shí)針?lè)较虻幕亓餍郎u區(qū)范圍將隨之減小，密度場(chǎng)的范圍相應(yīng)地也減??；當(dāng)Fre增大到一定程度后，進(jìn)口附近區(qū)域僅出現(xiàn)順時(shí)針?lè)较虻幕亓餍郎u，進(jìn)口附

19、近水流的摻混作用劇烈。            圖3 異重流潛入時(shí)的流場(chǎng)隨Fre的變化(T=0.8)Velocity field with the variation of Fre(T=0.8)in the case of density current plunging    圖4 Xp/H與Fre之間的關(guān)系Releation between Xp and Fre為進(jìn)一步研究異重流潛入形態(tài)與Fre之間的關(guān)系，我們還利用數(shù)值計(jì)算的結(jié)果，在

20、圖4中點(diǎn)繪了潛入點(diǎn)與進(jìn)口之間的無(wú)量綱距離Xp/H與Fre之間的關(guān)系。該圖清楚地反映了兩者之間的關(guān)系：Fre<0.6時(shí)，Xp/H隨Fre變化不大，進(jìn)口附近的流動(dòng)主受負(fù)浮力作用的影響，此時(shí)的流動(dòng)為卷流型；Fre>1.5時(shí)，隨著Fre增大，Xp/H迅速增大，進(jìn)口附近的流動(dòng)主受渾水在進(jìn)口處的初始動(dòng)量所控制，此時(shí)流動(dòng)為射流型。在卷流型和射流型之間存在一個(gè)過(guò)渡區(qū)，過(guò)渡區(qū)的流動(dòng)不僅具有一定的初始動(dòng)量，而且還受到負(fù)浮力作用，此時(shí)流動(dòng)為浮力射流型。若不考慮過(guò)渡區(qū)，則Fre<1時(shí)，流動(dòng)為卷流型；Fre>1時(shí)，流動(dòng)為射流型，與文獻(xiàn)3的結(jié)論相同。    

21、4.2 異重流潛入的發(fā)展過(guò)程異重流潛入水庫(kù)的過(guò)程中，流動(dòng)是非定常的，流場(chǎng)和密度場(chǎng)均隨時(shí)間變化。為研究異重流潛入隨時(shí)間的變化規(guī)律，本文根據(jù)數(shù)值計(jì)算結(jié)果，繪制了Fre=1.58的情況下，進(jìn)口附近區(qū)域分別在T=0.8，1.6，2.4，3.2時(shí)刻的速度場(chǎng)和密度場(chǎng)，如圖5和圖6所示。    圖5是不同時(shí)刻下的速度矢量圖：T=0.8時(shí)，因初始動(dòng)量的作用，在進(jìn)口附近將產(chǎn)生一順時(shí)針?lè)较虻幕亓餍郎u，隨著時(shí)間的增加，該回流區(qū)的范圍逐漸增大；當(dāng)水庫(kù)中的水體產(chǎn)生的阻力阻礙其向前進(jìn)一步運(yùn)動(dòng)時(shí)，渾水開(kāi)始下沉（T=1.6時(shí)刻），在該回流右側(cè)的自由表面附近的區(qū)域內(nèi)，將出現(xiàn)一逆時(shí)針?lè)较虻?/p>

22、環(huán)流；從T=2.4和T=3.2的速度矢量圖我們還可看出，隨著時(shí)間的增加，順時(shí)針回流因受左側(cè)邊界的限制，發(fā)展逐漸減慢，而右側(cè)環(huán)流因異重流沿底部繼續(xù)地向前運(yùn)動(dòng)，異重流與水庫(kù)水體之間剪切力逐漸發(fā)展，范圍將逐漸增大，水庫(kù)水體表面將有倒流現(xiàn)象產(chǎn)生。圖6是不同時(shí)刻下的相對(duì)有效密度等值線圖（/0）：當(dāng)T=0.8時(shí)，進(jìn)口附近區(qū)域基本處于初始動(dòng)量作用的控制，相對(duì)有效密度值沿水平方向逐漸延伸，隨著時(shí)間的增加，進(jìn)口動(dòng)量作用越來(lái)越小，因密度差而產(chǎn)生的負(fù)浮力作用卻逐漸增強(qiáng)，/0受負(fù)浮力和回流的影響而又逐漸降落，當(dāng)T=3.2時(shí)，/0貼到底部，在底部形成一股底流，直此異重流形成了。   

23、60;     圖5 異重流潛入過(guò)程中速度場(chǎng)隨時(shí)間的變化(Fre=1.58)Velocity field with the variation of time(Fre            摘本文采用變密度流基本方程和混合有限分析法，求解了突擴(kuò)邊界下的異重流潛入運(yùn)動(dòng)。通過(guò)對(duì)數(shù)值計(jì)算結(jié)果分         本篇論文是由3COME文檔頻道的網(wǎng)友為您在網(wǎng)絡(luò)上收集整

24、理餅投稿至本站的，論文版權(quán)屬原作者，請(qǐng)不用于商業(yè)用途或者抄襲，僅供參考學(xué)習(xí)之用，否者后果自負(fù)，如果此文侵犯您的合法權(quán)益，請(qǐng)聯(lián)系我們。=1.58)in the process of density current plunging            圖6 異重流潛入過(guò)程中密度場(chǎng)隨時(shí)間的變化(Fre=1.58)Density field with the variation of time(Fre=1.58)in the process of density current

25、plunging    圖7 T=3.2時(shí)刻的水平速度分量沿程分布圖(Fre=1.58)Longitudinal distribution of horizantal velocity component(Fre=1.58)as T=3.24.3 潛入前后渾水流速的沿程變化圖7是本文依據(jù)數(shù)值計(jì)算結(jié)果，點(diǎn)繪的Fre=1.58的異重流潛入時(shí)，在T=3.2時(shí)刻的水平速度分量的沿程分布，從該圖可以看出，異重流在潛入的過(guò)程中，水平速度分布沿程發(fā)生了根本的改變：明流最大流速點(diǎn)一般在水面處，潛入后，流速分布型式變化，且最大流速點(diǎn)的相對(duì)位置下移，這與文獻(xiàn)8的分析相同。5 結(jié)論本文采用混合有限分析法求解了二維Navier Stokes方程和濃度方程，并根據(jù)數(shù)值計(jì)算的結(jié)果，對(duì)異重流的潛入運(yùn)動(dòng)規(guī)律進(jìn)行了分析，結(jié)