2006-2007 > Academic research: refereed > Publication in refereed journalVersion of RecordPublishe

Jiang, W

Li, YS

Wai, WH

Zhan, JM

PolyU Institutional Repository

　　A?? 21?? 5?　　　　　　 ? ? ? ? ? ? ? ? ? 　　　　　　Ser. A, Vo .l 21, No. 5　　　2006? 9?　　　　　　　　 JOURNAL OF HYDRODYNAM ICS　　　　　　　　Sep.t , 2006????:1000-4874(2006)05-0613-06???????????????????*?? 1 , 　　??? 2 , 　　???2 , 　　??? 3(1. ???????????? ,?? 510610;2. ??????????????;3. ???????????? ,?? 510275)?　?:　??????????????????? , ??????????????????? , ??? 1978? 7? 20? 23?? 1978? 7? 22? 23?????????????????? , ?????????????????? ,?? , ???????????????????????? , ??????????????????????。?　?　?:　????;????;???;???????:　P731, O35　　　?????:AA large eddy smi ulation turbulence model forestuary using spline correctionJIANGW ei1 , 　　WEI Yong-kang2 ,LI Yok-sheung2 , 　　ZHAN Jie-m in3(1. Guangdong Prov incia l Resea rch Institu te ofW a te rConservancyandHydropow e r, Guangzhou 510610, China;2. Dep.t of C ivil? Structural Eng inee ring, The Hong Kong Poly technic Unive rsity,Hong Kong, China;3. Dep.t o fAppliedM echanics and Eng ineering, Zhong shanUn iversity,Guangzhou 510275, China)Abstrac t:　A la rge eddy sim ulation turbu lencemode l fo rPearl R iver e stuary is established in this paper. G rad ients o fw ater leve lin horizonta l p lan are co rrec ted using cubic spline. Dynam ic process in L ingd ing Yang o f Pear l R iver e stua ry from 11:00 pm in Ju ly20, 1978 to11:00 pm in July 22, 1978 is simu la ted. The num erica l results show tha t the proposed m ode l can be used to p redict theproce ss co rrectly. A t the sam e tim e, because the cubic sp line can natu ra lly be used to describe the ze ro g radien t o fw a te r leve l in hori-* ????:　2006-02-28????:　???????? (10272118), ??????? (RGC--Po lyU 5058 /02E)???????????????????????(20020558013),?????????(2KM 061045)????:　??(1963 ～ ),? ,????? ,????? , ??。DOI牶牨牥牣牨牰牥牱牰牤j牣cnki牣cjhd牣牪牥牥牰牣牥牭牣牥牥牴zon tal plan near the bank, them ode l can be used to describe the loca l vo rtexes near the bank c lea rly.Key words:　 large eddy simu lation, cubic sp line;FDM;estuary1　引言????????? 、?????????????? 。?????????????? , ???????????????。?? ,?????????????? 。?????????????????????????????????。?? ,????????????????????????????????? Princeton (B lumberg andM ellor, 1987)?? 。??? , ???? (La rge EddyS imula tion)?????????????????????。???????????????????????? (Moeng, 1984;Mason, 1994),???? ,????????????????? ,???????????????????????(McW illiams e t a.l , 1993;Sky lling stad and Denbo,1995;Zhan and Li 1998, 2001). ??????? ,????? Smagorinsky-Lilly?? ,?????????????? ,???????? S ij???? 。????????????? ,?? Kolmogorov -5 /3??(Li lly, 1967), ????????????????????????????。2　大涡模拟和垂直方向的涡黏参数化????????? ,?????????? ,??:  u x + –v y +  w z =0 (1)  u t + uu x + vu y + wu z - f–v =-1ρ –p x (2) –v t + uv x + vv y + wv z +f u =-1ρ –p y (3)  w t + uw x + vw y + ww z =-1ρ –p z - g (4)???????????? ,???????????? ,?????????????????? , ???????????????????(Zhan? Li, 1998, 2001)。?????? ,?????????u iuj =u i u j +L(ui , u j) +u′iu′j (5)?? , i, j =1? 2?L(ui , u j) =Δ212γ ui x u j x +Δ222γ u i y u j y +Δ232γ u i z u j z?? γ????? , ?? 6. 0。Δi? i??????? 。???????????????? ,???(5)? i, j=1 , 2, 3??? ,????????? ,?????? 。?? ,??????? ,?????????。?? , (5)???? uu, uv? vv 。?? ,????? , ?? Sm ago rinsky-Lilly????????(5)??? u′iu′j????? ,?u′iu′j -23δi, j E =- 12μ(  ui xj +  uj xi) (6)??μ=(CsΔ)2 (  u x)2 +( –v x)2 +12(  u y + –v x)2 , E = 12u′iu′j (7)614 ?　?　?　?　?　?　?　?　?　　　　　　　　　2006?? 5?? δi, j? δ??;Δ??????????? Δ =(Δ1Δ2)1 /2 , Δ1 =Δx? Δ2 =Δy;Cs??????(???? 0. 2),????? Ko lmogorov -5 /3??(Lilly, 1967)??; E ???????????? 。2 E??????? ,????????? ,????????????: P =–p +23 E (8)?????? ,L(ui , u j) =Δ212γ ui x uj x +Δ212γ u i y u j y (9)- ρu′w′ z =τsx - τbxh +η (10)- ρv′w′ z =τsy - τbyh +η (11)?? , τs , τb??????????????? 。???????????? ,?? ,?τbx =ραSu , 　　τby =ραSv (12)?? , S = u2 +v2 , α=g /C2 , ???? C =h16 /n ,n????? ,??? n =0. 012。?????????????? ,????????? η??? 。?? ,???????????: η t =-  [ (h +η) u] x -  [ (h +η)–v] y (13)  u t + [ u  u +L( u  u)] x - [ u –v +L( u –v)] y - f–v =- g η x +  x(μ  u x) +12  x[ μ(  u y + –v x)] -τbxρ(h +η) (14) –v t + [–v  u +L(–v  u)] x - [–v –v +L(–v –v)] y - f u =- g η y +12  x[ μ(  u y + –v x)] +  y(μ –v y) -τbyρ(h +η) (15)? 1　????3　数值求解方法???????????????????????????????? ,????????????? 、???????? 。?? ,??????????????????? ,??????? , ????? (? 1)??????????? ,?????????????? , ??????????????????? 。?????????? (13)????????????? 。?????:ηn+1 - ηnΔt =- ( [ (h +η) u] n+1 /2 x + [ (h +η)–v]n+1 /2 y (16)??()n+1 /2 = 12[ ()] n+1 +()n ] +O(Δt2)??(16)??:(ηn+1)k+1 - ηnΔt =- 12 615???:???????????????????( [ (h +η) u] n+1 x + [ (h +η)–v]n+1 y )k -12( [ (h +η) u] n x + [ (h +η)–v]n y ) (17)?? ,? k =0? , [ ()n+1 ] 0 =()n。?????? η???????? 10-5。???????????? , (ηn+1)k+1?? ηn+1。?? ,???  un+1? –vn+1。???????? ,?????????? ,??????? 0,??????????????????????????????? ,?? (i, j)???  η x? η y 。???? ,????? ,??????????? ,??????????????????? ,? , η t +c η x =0 (18)????????????? ,?? ,  η t≈ 0 。???????????????????。4　模型在河口问题的适用性和优点??? (?? 2? )???????? ,????? 2110 km2。????????———?????????????。?????? , ??? 4 km。??????? (???),?????。?????????????????? 。??? 22°45′(???? )????? 22°20′(?????)。 5m ?????????????????????? ,??????????????? ,?????? ,??????。 20m?????? ,?????? 20 km ,??????? ,????? 、?????? ,???????。??????? ,?????? ,??????? 。??????????? ,??? 、??。?? ,???? 、???????? (?? 16 m;?? 16km)———????(?? 10m;?? 34 km)--?????(?????????? ,?????? ,????? )???。??????? 11. 45 m。???? ,????????-????-???? ,???????????????。???? 8. 6 m ,????????????????????。? 2　???????????? 3　?????????????????????????? , ??????????? ,????????? 2?? 3。???? ,????????:Δx =Δy =400m , ????? Δt =10 s。??? 1978? 7? 20? 23?? 1978? 7? 22? 23?????????????????? ,?????????????? 4?? 5??。???? ,???????? ,?????????????? 。?????????????? ,?????? ,???????????616 ?　?　?　?　?　?　?　?　?　　　　　　　　　2006?? 5?? 4　???????????????? 5　??????????????0. 6????????????????????? 6　?????????? 7　???????????? ,???????? ,?????????????????? 。?????????????????(???? ,?????),???????????????? 、?????????????????? ,????????????????????????? ,??????????????????? ,?????????。? 6 ～? 7?????????????????????? ,???? ,??????????????????。???? ,?????????617???:???????????????????????????????????? ,???????????? 。?????? ,???????????? ,???? ,???????? ,?????????????????? 。5　结论????????????????????????????????????????????????。???????????? (? )?? 0?????? ,????????????????????????? 。?????????????????????????? 。?　?　?　?:[ 1] 　BLUMBERG A F and MELLOR G L. A desc ription of athree-dimensional co asta l ocean circulation m ode l[ M ] . In:Three-Dim ensional Coastal Ocean M ode l, Vo .l 4, ed itedby Heaps, pp. 208, Am erican Geophy sica l Union, W ash-ing ton, D. C. , 1987.[ 2] 　 LILLY D K. The representa tion of sm a ll-sca le turbulencein nume rical simu la tion experim ents[ C] . P roceedings ofthe tenth IBM sc ientific com puting sym posium on env iron-m ental sciences. Thom as J. W a tson Research Center,York town He igh ts.[ 3] 　MASON P J. Large-eddy simu lation:A critical rev iew ofthe technique[ J] . Q. J. R. M e teo ro .l Soc. , 1994, 120:1-26.[ 4] 　MCW ILLIAMS J C, GALLACHER P C , MOENG C H andWYNGAARD. J C. Mode ling the oceanic plane ta ry bound-ary laye r[ M ] . In:La rge Eddy Simu la tion of Comp lex Eng i-neering and Geophysical F low s, ed ited by Bo ris Ga lperinand S teven A. O rzag, 1993, 441-454, Cam bridge Univ.P ress, New York.[ 5] 　MOENG CH. A la rge-eddy-simu lation m ode l for the studyof p lane tary bounda ry-laye r turbulence[ J] . Jou rna l o f theA tmosphe ric Sc ience s, 1984, 41(13):2052-2062.[ 6] 　SKYLLINGSTAD E D and DENBO D W. An ocean large-eddy simu lation of Langmu ir circu lations and convec tion inthe surface m ixed lay er[ J] . J. Geophy s. Res. , 1995,100:8501-8522.[ 7] 　ZHAN Jie-m ing and L I Y S. A hydrodynam ic larg e-eddysim ulation-model using d irec tiona l eddy visco sity me thod invertica l direction. Journal of Hydrodynam ics, Se r. B,1998, 10(1):36-43.[ 8] 　ZHAN Jie-m ing, L IY S. A Large Eddy Simu la tion Turbu-lenceM ode l for Coasta l Seas and Shallow W a ter P rob lem s[ J] . Jou rnal o f H ydrodynam ics, Se r. B, 2001, 13(4):24-27.618 ?　?　?　?　?　?　?　?　?　　　　　　　　　2006?? 5?

A large eddy simulation turbulence model for estuary using spline correction

Abstract

Similar works

Full text

Available Versions

PolyU Institutional Repository