張耀屹，徐績青，許錫賓

(重慶交通大學(xué)河海學(xué)院，重慶400074)

波浪輻射應(yīng)力的概念是由Longuet-Higgins和Stewart提出的，他們在牛頓流體力學(xué)定理和速度沿水深均勻分布假定的基礎(chǔ)上，將波浪輻射應(yīng)力定義為由于波浪現(xiàn)象而產(chǎn)生的動量超通量［1］;可見輻射應(yīng)力只是一種通俗的稱謂，它是被平均在水深上的力，并不能完全體現(xiàn)流場中流體具體某一點的應(yīng)力。

雷諾輸運方程表明:一個物質(zhì)體系內(nèi)某種流體廣延量的增長率，等于體系在該時刻所占的那個空間區(qū)域中同一物理量的增長率加上單位時間內(nèi)由區(qū)域邊界流出的該物理量的總通量;進而雷諾輸運方程導(dǎo)出的流體動量方程探索出:在空間的一個固定區(qū)域中，單位時間內(nèi)流體總動量的增加，等于該時間內(nèi)通過區(qū)域邊界面進入?yún)^(qū)域的動量(包括由于流體進出該區(qū)域所攜帶的動量的凈增量和由于表面力作用傳遞到區(qū)域內(nèi)的動量)，加上在外場力作用下引起的區(qū)域內(nèi)流體動量的改變［2］。

通過對以上2組概念的對比，不難看出可以運用雷諾輸運方程求出由于波浪現(xiàn)象而產(chǎn)生的動量超通量，進而表示出波浪的輻射應(yīng)力張量。

1 輻射動量應(yīng)力的產(chǎn)生

1.1 傳統(tǒng)輻射應(yīng)力張量的表達式及應(yīng)用［3－5］

1.1.1 傳統(tǒng)輻射應(yīng)力張量的表達式

通常，輻射應(yīng)力的張量表達式為:

式中:sxx表示作用在垂直于x軸平面的x方向的剩余動量流，稱作輻射應(yīng)力在x方向上的主分量;syy表示作用在垂直于y軸平面上的y方向的剩余動量流，稱作輻射應(yīng)力在y方向上的主分量;sxy表示作用在垂直于x軸平面上的y方向動量流，稱作輻射應(yīng)力在垂直于x軸平面上的切向分量;syx表示作用在垂直于y軸平面上的x方向的動量流，稱作輻射應(yīng)力在垂直于y軸平面上的切向分量。以不同的波理論的水質(zhì)點速度及波壓力公式可以計算得到不同的輻射應(yīng)力表達式。據(jù)微幅波理論，波向角為α?xí)r輻射應(yīng)力的張量形式表達式為:

1.1.2 傳統(tǒng)輻射應(yīng)力張量推導(dǎo)的減水方程

用傳統(tǒng)的輻射應(yīng)力公式求得規(guī)則波的減水方程為:

式中:A為波浪振幅;k為波數(shù);h為水深。

1.2 “輻射動量應(yīng)力”張量的表示

1.2.1 理論探討

Longuet-Higgins和Stewart將作用于單位面積水柱體的總動量流的時均值減去無波浪作用時的靜水壓力定義為波浪剩余動量流，又稱作輻射應(yīng)力［1］。而實際上，由雷諾輸運方程出發(fā)推出動量流ρvivj才是真正的應(yīng)力張量。沿整個控制體的外表面做曲面積分，得到的是一個合力，這才是“輻射應(yīng)力”的本質(zhì)特征。由于:

1)動量流的周期時均值仍然是一個張量(其變換法則不受時間的影響);

根據(jù)傳統(tǒng)輻射應(yīng)力的概念，結(jié)合雷諾輸運方程，定義:輻射動量應(yīng)力=-動量流+反向應(yīng)力;由此可見:輻射動量應(yīng)力是一個真正的應(yīng)力張量。其表達式可表示為:

即:

式中:ρvivj為動量流的應(yīng)力張量分量為反向應(yīng)力的張量表達式。

1.2.2 公式證明［2，6-7］

靜水壓力的曲面積分(對坐標)已與控制體的重力抵消，對于二階斯托克斯波的作用，靜水壓力曲面積分等于控制體的重力，如圖1。

圖1 任意運動流體質(zhì)點受力Fig.1 Arbitrary exercise fluid particles stress

從動量流出發(fā)得到的該方程符合積分形式的動量方程，還可變?yōu)槲⒎中问降倪\動方程，從而得出“輻射動量應(yīng)力”張量表達式(2)，式(2)在流體中的確發(fā)揮的是輻射應(yīng)力的作用。

2 用“輻射動量應(yīng)力”求解減水方程

“輻射動量應(yīng)力”在海岸動力學(xué)中的應(yīng)用，即波浪增減水方程的求解。

波浪增減水方程可以看成沿岸流方程的特例，沿x方向(見圖2)的沿岸流方程:

其中:Cij為“輻射動量應(yīng)力”，對應(yīng)二階斯托克斯波(pij為反向應(yīng)力)。

這時，假定u=0，平直海岸情況下波高在沿岸方向均勻一致，即u與項均為0。沿岸流方程可簡化為:

因為，波浪沿x方向傳播，即Cxy=0，Cxz=0，那么此方程在本質(zhì)上仍然表現(xiàn)為波浪正向傳播時所產(chǎn)生的增減水作用，即

因為，二階斯托克斯波的時均水面下降η*，所以波動壓力的時均值取z'為坐標原點，即:

破波帶外:

整理得:

綜上所述:

由此可見，以上用“輻射動量應(yīng)力”求解的減水方程的結(jié)果式(4)忽略掉高階無窮小量后與用傳統(tǒng)輻射應(yīng)力張量算出的減水方程的結(jié)果式(1)形式上是一致的;另外，式(4)比式(1)多了一個高階無窮小的修正量值，理論上增強了運算結(jié)果的可靠性。

3 結(jié)論

1)用波浪輻射動量應(yīng)力分析波浪減水現(xiàn)象，推導(dǎo)出的減水方程與傳統(tǒng)輻射應(yīng)力張量計算出的結(jié)果基本吻合，具有可靠性。

2)相對于傳統(tǒng)概念的“輻射應(yīng)力”，“輻射動量應(yīng)力”具有明顯的優(yōu)越性。它體現(xiàn)的是真正的應(yīng)力張量;張量各分量沿水深的變化一目了然;作用力的計算比較簡單;概念非常清楚，體現(xiàn)了力學(xué)本源;拓展了輻射應(yīng)力的適應(yīng)范圍，傳統(tǒng)的輻射應(yīng)力概念在牛頓流體力學(xué)定理和速度沿水深均勻分布假定的基礎(chǔ)上建立的，而本表達式?jīng)]有相關(guān)的限制。

3)運用輻射動量應(yīng)力還可來求解波浪的增水方程和探尋沿岸流的生成機理。

［1］Longuet-higgns M S，Stewartr W.Radiation stresses in water waves:a physical discussion with application［J］.Deep Sea Resarch，1964，11(5):529-562.

［2］莊禮賢，尹協(xié)遠，馬暉揚.流體力學(xué)［M］.合肥:中國科學(xué)技術(shù)大學(xué)出版社，1991:83-89.

［3］馬錫超，卜令濤，王璇.輻射應(yīng)力理論在近岸區(qū)的應(yīng)用及其研究進展［J］.海岸工程，2008，27(1):73-77.MA Xi-chao，BU Lin-tao，WANG Xuan.Radiation to the application in the stress theory and its research progress inshore coastal engineering［J］.Coastal Engineering，2008，27(1):73-77.

［4］鄒志利.海岸動力學(xué)［M］.北京:人民交通出版社，2010:107-122.

［5］溫秀媛.波浪輻射應(yīng)力計算及其對低頻波浪的影響［D］.大連:大連理工大學(xué)，2009.

［6］鄭永紅，沈永明，邱大洪.結(jié)合拋物型緩坡方程計算波浪輻射應(yīng)力［J］.海洋學(xué)報，2000，22(6):110-116.ZHENG Yong-hong，SHEN Yong-ming，QIU Da-hong.Combined wave radiation stress calculation with Parabolic mild-slope equation［J］.Acta Oceanologica Sinica，2000，22(6):110-116.

［7］鄭金海，嚴以新.波浪輻射應(yīng)力張量的垂向變化［J］.水動力學(xué)研究與進展，2001，16(2):246-253.ZHENG Jin-hai，YAN Yi-xin.The Wave radiation stress tensor vertical change［J］.Journal of Hydrodynamics，2001，16(2):247-254.

［8］李孟國，王正林，蔣德才.近岸波浪傳播變形數(shù)學(xué)模型的研究與進展［J］.海洋工程，2002，20(4):43-57.LI Meng-guo，WANG Zheng-lin，JIANG De-cai.Offshore wave propagation of deformation mathematical model research andprogress［J］.Maritime Engineering，2002，20(4):43-57.

［9］林祥，尹寶樹，侯一筠，等.輻射應(yīng)力在黃河三角洲近岸波浪和潮汐風(fēng)暴潮相互作用中的影響［J］.海洋與湖沼，2002，33(6):615-621.LIN Xiang，YIN Bao-shu，HOU Yi-jun，et al.Radiation stress in the Yellow River delta offshore wave and tidal storm surge of interaction effects［J］.Oceanologia and Limnologia Sinica，2002，33(6):615-621.

［10］蔡曉禹，凌天請，唐伯明，等.波浪對散體岸坡沖刷破壞的機理［J］.重慶交通學(xué)院學(xué)報，2006，25(2):73-76.CAI Xiao-yu，LING Tian-qing，TANG Bo-ming，et al.The waves to scour the damage mechanism of Loose material accumulation［J］.Journal of Chongqing Jiaotong University，2006，25(2):73-76.