• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于大渦模擬耦合模式的小尺度海氣相互作用研究*

    2020-11-24 01:20:50孫丹譯
    海洋與湖沼 2020年6期
    關(guān)鍵詞:收支湍流動(dòng)能

    孫丹譯 李 爽

    基于大渦模擬耦合模式的小尺度海氣相互作用研究*

    孫丹譯 李 爽①

    (浙江大學(xué)海洋學(xué)院 舟山 316021)

    大氣和海洋是影響地球氣候系統(tǒng)的兩個(gè)重要因素, 它們之間的相互作用是海洋和大氣研究的重要課題, 海氣耦合模式則是研究海氣相互作用的重要工具, 而海氣耦合模式重點(diǎn)考慮的參數(shù)是海氣通量。針對(duì)傳統(tǒng)的大尺度海氣耦合模式缺少湍流尺度分析的問(wèn)題, 本文使用并行大渦模擬海氣耦合模式(The Parallelized Large-Eddy Simulation Model, PALM), 在小尺度上探究風(fēng)速對(duì)海氣通量及湍流動(dòng)能收支(Turbulence Kinetic Energy Budget, TKE Budget)的影響, 設(shè)置了5、10和15m/s三種地轉(zhuǎn)風(fēng)速度對(duì)大氣邊界層(Atmospheric Boundary Layer, ABL)和海洋混合層(Oceanic Mixed Layer, OML)進(jìn)行海氣耦合模擬。研究表明: 海氣通量的分布與風(fēng)速大小密切相關(guān), 風(fēng)速越大, 凈熱通量和浮力通量相對(duì)越大, 由于溫度上升導(dǎo)致海水蒸發(fā)加劇, 使得大氣的淡水通量增大; 海洋湍流動(dòng)能收支各項(xiàng)在近海面處受風(fēng)速影響較大, 且隨著深度加深而逐漸減弱。本研究初步展示了小尺度海氣耦合模式在海氣通量研究中的應(yīng)用, 對(duì)進(jìn)行小尺度海氣相互作用研究具有一定的意義。

    大渦模擬; 海氣通量; 海氣耦合; 湍流動(dòng)能收支

    大氣和海洋是地球系統(tǒng)中的兩個(gè)重要組成部分, 它們之間的相互作用是地球氣候系統(tǒng)的關(guān)鍵調(diào)節(jié)因素之一, 而單一的大氣或海洋模式無(wú)法完整描述此過(guò)程, 所以海氣耦合模式被提出并得到快速發(fā)展, 其發(fā)展經(jīng)歷了從全球到區(qū)域的過(guò)程, 其中全球海氣耦合模式主要用于研究大尺度海氣相互作用, 然而全球模式的分辨率無(wú)法準(zhǔn)確描述中小尺度海氣相互作用, 因此區(qū)域海氣耦合模式憑借其較高的分辨率, 能夠?qū)庀嗷プ饔眠M(jìn)行更細(xì)致的研究, 其廣泛應(yīng)用于數(shù)值模擬、預(yù)報(bào)和區(qū)域氣候變化動(dòng)力降尺度研究(Peng, 2012)。海氣通量則是海氣相互作用中的一個(gè)重要研究課題(陳烈庭, 1991), 所以在全球或區(qū)域的海氣耦合模式中考慮海氣通量的變化對(duì)深入研究海氣相互作用具有重要意義。Zhou等(2000)將海氣淡水通量的交換引入耦合模式Global-Ocean-Atmosphere- Land-System-model(GOALS)成功地控制了氣候漂移現(xiàn)象(Neelin, 1992; 陳克明等, 1997)。海氣界面處的動(dòng)量通量也叫風(fēng)應(yīng)力, 是海氣相互作用研究的重點(diǎn)之一, 也是海氣耦合模式中的一個(gè)重要參數(shù)(馮興如等, 2018)。Ma等(2010)在海氣耦合模式FOAM中通過(guò)控制變量證明風(fēng)應(yīng)力對(duì)氣候的重要影響。Zhou等(2015)將AGCM模式和一個(gè)混合層海洋模式耦合用以探究海氣耦合對(duì)大西洋多年代際振蕩的影響。

    區(qū)域海氣耦合模式的關(guān)鍵技術(shù)是在系統(tǒng)質(zhì)量和能量守恒的前提下實(shí)現(xiàn)高效的海氣參數(shù)的交互傳遞, 現(xiàn)如今耦合模式普遍采用“耦合器”技術(shù)(Coupler), 而另一個(gè)技術(shù)難題是側(cè)邊界條件的處理, 對(duì)于側(cè)邊界的處理直接影響模擬結(jié)果(Peng, 2012)。上述這些海氣耦合模式對(duì)于處理海氣之間的質(zhì)量、動(dòng)量和能量的交換問(wèn)題通常是依靠經(jīng)典參數(shù)化方案, 而且模式的空間分辨率并不能準(zhǔn)確描述小尺度的海氣相互作用。2001年海氣相互作用耦合邊界層低風(fēng)速實(shí)驗(yàn)(The Coupled Boundary Layers and Air-Sea Interaction- Low Winds, CBLAST-Low)項(xiàng)目在沿海部署了一套觀測(cè)系統(tǒng)用以提供海氣數(shù)據(jù), 取得了很好的觀測(cè)結(jié)果(Austin, 2002; Vickers, 2006; Edson, 2007)。但是相對(duì)于觀測(cè)系統(tǒng)造價(jià)昂貴以及單點(diǎn)觀測(cè)所帶來(lái)結(jié)果特殊性的特點(diǎn), 數(shù)值模擬作為一種新興的研究方式近年來(lái)發(fā)展迅速。因?yàn)榇髿夂秃Q笾械耐牧鲿r(shí)刻影響著海氣通量的交換, 所以模式需要充分考慮湍流作用。以往的湍流模式大多是雷諾平均模式(Reynolds Average Navier-Stokes, RANS), Peng等(2012)中介紹的區(qū)域海氣耦合模式大都使用RANS方法, 而大渦模擬(Large Eddy Simulation, LES)與RANS不同, 其原理是基于尺度的分離: 在模擬中直接計(jì)算大尺度的脈動(dòng), 對(duì)于小尺度的湍流效應(yīng)則采用參數(shù)化的方案(Smagorinsky, 1963; Deardorff, 1970), 因?yàn)槠漭^為精確且節(jié)省計(jì)算成本的特性, 已經(jīng)成為大氣和海洋領(lǐng)域中研究湍流的一種成熟方法。Esau(2014)采用并行大渦模擬模式(The Parallelized Large-Eddy Simulation Model, PALM), 研究了大氣邊界層(Atmospheric Boundary Layer, ABL)與海洋混合層(Oceanic Mixed Layer, OML)之間在無(wú)剪切對(duì)流情況下的湍流相互作用。

    湍流作為大氣和海洋中顯著的動(dòng)力過(guò)程, 是影響海氣通量的重要因素。大渦模擬允許湍流場(chǎng)的三維實(shí)現(xiàn), 可以用湍流動(dòng)能平衡方程中的耗散項(xiàng)、剪切項(xiàng)和傳輸項(xiàng)等來(lái)刻畫(huà)湍流過(guò)程, 這是大尺度模式和大多數(shù)實(shí)驗(yàn)研究所不能實(shí)現(xiàn)的(Maronga, 2015)。所以在湍流尺度運(yùn)用大渦模擬來(lái)探究風(fēng)速對(duì)海氣通量的影響和對(duì)深入了解海氣相互作用過(guò)程具有重要的意義。為此本文在Esau(2014)的基礎(chǔ)上引入風(fēng)速, 用大渦模擬海氣耦合模式探究風(fēng)速對(duì)海氣通量和海洋湍流動(dòng)能收支各項(xiàng)分布的影響。

    1 模式簡(jiǎn)介

    1.1 PALM簡(jiǎn)介

    PALM是由德國(guó)漢諾威萊布尼茨大學(xué)開(kāi)發(fā)的用于大氣和海洋邊界層的一種湍流分解大渦模擬模式, 已應(yīng)用于各種大氣和海洋邊界層的模擬超過(guò)15年(Huq, 2019)。經(jīng)過(guò)優(yōu)化, PALM可用于大規(guī)模并行計(jì)算機(jī)架構(gòu), 這也是它能進(jìn)行模式耦合的技術(shù)基礎(chǔ), 關(guān)于PALM的詳細(xì)介紹可見(jiàn)https://palm.muk.uni- hannover.de/trac/wiki。

    根據(jù)公式(1)可得水平平均的SGS-TKE方程(Skyllingstad, 2000)為:

    其中:

    是壓強(qiáng),m是次網(wǎng)格渦粘系數(shù)。公式(6)中等號(hào)右邊各項(xiàng)依次是: 速度剪切項(xiàng)、浮力項(xiàng)、次網(wǎng)格耗散項(xiàng)、湍流動(dòng)能傳輸項(xiàng)、壓力項(xiàng)和耗散項(xiàng)。下文中湍流動(dòng)能收支的耗散項(xiàng)包含了次網(wǎng)格耗散, 而湍流動(dòng)能傳輸項(xiàng)包含了次網(wǎng)格湍流動(dòng)能傳輸。

    1.2 海氣耦合模式簡(jiǎn)介

    大氣和海洋的時(shí)間步長(zhǎng)是單獨(dú)設(shè)置的, 不要求相等, 然后以用戶(hù)指定的頻率執(zhí)行耦合, 目前, 該耦合要求大氣和海洋的水平域的大小相等, 而耦合則是通過(guò)海面數(shù)據(jù)的雙向雙線性插值實(shí)現(xiàn)的。此外, 在大氣和海洋中, 由于模式穩(wěn)定時(shí)間的不同, 即動(dòng)能輸出達(dá)到穩(wěn)定的時(shí)間不同, 可以為大氣和海洋先執(zhí)行非耦合的前驅(qū)運(yùn)行, 然后再進(jìn)行耦合的重新啟動(dòng)運(yùn)行, 這可以減少計(jì)算量, 也是本文所使用的方法。

    1.3 模式設(shè)置

    本文的模式設(shè)置中大氣部分的網(wǎng)格數(shù)是100× 100×60(,,), 網(wǎng)格大小是1m×1m×5m; 海洋部分的網(wǎng)格數(shù)是100×100×100(,,), 網(wǎng)格大小是1m×1m× 1m。本文中PALM海氣耦合分為兩個(gè)部分: 前驅(qū)運(yùn)行和耦合運(yùn)行, 兩者的設(shè)置大致相同, 時(shí)間步長(zhǎng)方案均采用三階Runge-Kutta方法, 側(cè)邊界條件均是周期性邊界條件, 大氣(海洋)模式中的下(上)邊界采用自由滑移邊界條件, 大氣模式的表面溫度是270.6K, 溫度梯度是0.015K/m, 海洋模式的溫度是272.65K, 鹽度是34.8, 緯度設(shè)置為45°N。本文不考慮海洋分層和海洋背景流; 前驅(qū)運(yùn)行中大氣模式依靠風(fēng)速輸入啟動(dòng), 海洋模式則是在運(yùn)行初始加入一個(gè)微小擾動(dòng)用于模式啟動(dòng), 具體表現(xiàn)是在海表面以熱通量的形式輸入模式, 該擾動(dòng)不會(huì)影響模擬結(jié)果。前驅(qū)運(yùn)行中大氣模式的運(yùn)行時(shí)間是7200s, 海洋模式則是40000s; 耦合運(yùn)行的時(shí)間均是86400s, 耦合步長(zhǎng)均是30s。這樣設(shè)置的原因是大氣和海洋的時(shí)間尺度不一致且海洋的響應(yīng)時(shí)間比大氣慢, 所以在模式前驅(qū)運(yùn)行中的穩(wěn)定時(shí)間海洋比大氣長(zhǎng)很多。本文通過(guò)設(shè)置三種地轉(zhuǎn)風(fēng)速度(5m/s,10m/s,15m/s)來(lái)探究其對(duì)海氣通量及海洋湍流動(dòng)能收支的影響。

    1.4 模式驗(yàn)證

    Esau(2014)運(yùn)用大渦模擬海氣耦合模式探究ABL與OML之間的湍流相互作用, 然而他并沒(méi)有考慮風(fēng)速對(duì)海氣湍流相互作用的影響, 為此本文在Easu(2014)的基礎(chǔ)上引入風(fēng)速, 使用耦合模式探究風(fēng)速對(duì)海氣相互作用的影響。在此之前, 本文先重復(fù)Easu(2014)的結(jié)果, 來(lái)證明該耦合模式在海氣相互作用研究中的可行性。

    耦合模式初始設(shè)置: 網(wǎng)格數(shù)128×128×64(,,), 網(wǎng)格大小25m×25m×10m(大氣), 25m×25m×5m(海洋)。將ABL中位溫的初始垂直梯度設(shè)為0.015K/m, 初始相對(duì)濕度為0%(干燥空氣)。OML的初始溫度和鹽度剖面來(lái)自Esau(2014): 海面溫度為275.65K, 鹽度為34.8, 大氣邊界層表面溫度為270.6K, 耦合模擬時(shí)間為22h。湍流是由ABL中40—90m層和OML中-40—90m層規(guī)定的速度場(chǎng)中的弱擾動(dòng)引起的, 其中湍流數(shù)據(jù)每60s采樣一次, 二維數(shù)據(jù)每5min采樣一次。

    圖1顯示了溫度、凈熱通量、動(dòng)量通量在最后一小時(shí)的時(shí)間平均剖面, 虛線是Esau(2014)的結(jié)果, 實(shí)線是耦合模式運(yùn)行結(jié)果。在圖1中, 兩者溫度值的分布相差不大, 凈熱通量和動(dòng)量通量的分布也大致相同, 本文將兩者的差異歸結(jié)為以下幾點(diǎn):

    (1) 模式直接取Esau (2014)文獻(xiàn)中的溫鹽初始剖面, 并沒(méi)有得到其原始數(shù)據(jù), 所以導(dǎo)致結(jié)果存在微小偏差。

    (2) 模式版本經(jīng)過(guò)更新?lián)Q代。

    (3) 關(guān)于耦合模式中前驅(qū)運(yùn)行的細(xì)節(jié)文獻(xiàn)中沒(méi)有詳細(xì)說(shuō)明, 可能模式輸入存在些許偏差。

    圖1 溫度、凈熱通量和動(dòng)量通量的垂向分布

    注: 實(shí)線: 模式運(yùn)行的結(jié)果; 虛線: Esau(2014)的結(jié)果

    2 模擬結(jié)果

    2.1 風(fēng)速和流速

    耦合模式由前驅(qū)運(yùn)行在大氣和海洋部分中產(chǎn)生穩(wěn)定的湍流, 之后再進(jìn)行耦合, 考慮到耦合模式有一個(gè)啟動(dòng)過(guò)程, 將垂向動(dòng)能進(jìn)行水平平均得到圖2, 橫坐標(biāo)是時(shí)間, 縱坐標(biāo)是高度, 在圖2中大氣部分從模式啟動(dòng)開(kāi)始就一直穩(wěn)定, 而海洋需要響應(yīng)時(shí)間, 所以從圖2中可以看到整體模式穩(wěn)定時(shí)間在耦合開(kāi)始1h后。

    地轉(zhuǎn)風(fēng)速為5m/s時(shí), 大氣的垂向動(dòng)能在0—150m的區(qū)域內(nèi)變化明顯, 其最值約為-0.05m2/s2, 相應(yīng)的海洋受風(fēng)速影響的深度約為20m, 垂向動(dòng)能的最值約為在-2×10-5m2/s2。地轉(zhuǎn)風(fēng)速為10m/s時(shí), 大氣受風(fēng)速影響的高度約為250m, 垂向動(dòng)能的最值約為-0.1m2/s2; 海洋受風(fēng)速影響的深度為25m, 垂向動(dòng)能的最值在-5×10-5m2/s2。地轉(zhuǎn)風(fēng)速為15m/s時(shí), 大氣受風(fēng)速影響的高度增大為300m, 垂向動(dòng)能的最值為-0.2m2/s2; 海洋受風(fēng)速影響的深度約為35m, 垂向動(dòng)能的最值為-1×10-4m2/s2。

    除了垂向動(dòng)能的時(shí)間變化外, 本文還關(guān)注了大氣海洋中速度分量在垂向上的變化, 將數(shù)據(jù)進(jìn)行時(shí)間平均后得到大氣風(fēng)速和海洋流速垂向剖面圖3。圖3中實(shí)線表示方向, 虛線表示方向, 紅線表示風(fēng)速為5 m/s, 藍(lán)線表示風(fēng)速為10 m/s, 黑線表示風(fēng)速為15m/s; 為了便于展示, 本文將海洋的流速數(shù)值放大100倍。

    圖3展示了耦合模式中大氣和海洋中速度分量的分布。由于模式中只在方向添加地轉(zhuǎn)風(fēng), 所以在圖3中海洋中方向流速?gòu)?開(kāi)始。隨著風(fēng)速的增加, 海洋流速逐漸增大且呈現(xiàn)向下延伸的趨勢(shì), 并在-50m左右迅速減小。

    圖2 垂向動(dòng)能時(shí)間序列

    注: a—c: 大氣的垂向動(dòng)能; d—f: 海洋的垂向動(dòng)能; a, d: 風(fēng)速5m/s; b, e: 風(fēng)速10m/s; c, f: 風(fēng)速15m/s

    圖3 速度剖面圖

    注: 實(shí)線表示, 虛線表示; 為了便于展示, 將實(shí)際海洋流速數(shù)值×102后繪制此圖

    2.2 海氣通量

    PALM中的海氣耦合是ABL和OML之間的耦合, 耦合是通過(guò)一個(gè)大氣-海洋界面來(lái)保存每個(gè)網(wǎng)格單元內(nèi)的質(zhì)量(濕度、鹽度)、動(dòng)量和熱量的湍流通量。ABL在OML的頂部邊界提供了通量, 而OML在ABL的底部邊界提供了速度和溫度。

    通過(guò)模式所得數(shù)據(jù)對(duì)其進(jìn)行時(shí)間平均得到圖4:

    由圖4可以直觀的看到不同風(fēng)速下海氣通量的垂向變化。在海氣界面處, 風(fēng)速越大, 大氣中的凈熱通量和淡水通量也相應(yīng)增大, 而浮力通量變化不大; 海洋中則是凈熱通量和浮力通量增大, 鹽度通量為負(fù)表示海面蒸發(fā)加劇。而且隨著風(fēng)速的增大, 各通量在海洋中向下傳遞的距離越大; 而在大氣中, 凈熱通量和淡水通量也相應(yīng)變大, 浮力通量在0—100m變化不大, 在100—300m變化明顯。

    本文取大氣10m處和水下2m處的通量值進(jìn)行對(duì)比得到圖5。圖5a中: 海氣的凈熱通量隨風(fēng)速增大而增大, 且隨著海氣之間相互傳遞熱通量, 兩者關(guān)系呈正相關(guān); 圖5b中: 隨著風(fēng)速增大, 海面受熱蒸發(fā)加劇, 導(dǎo)致大氣淡水通量增大, 而海洋鹽度通量(負(fù))增大。

    圖4 凈熱通量、淡水通量和浮力通量的垂向剖面圖

    注: a: 大氣和海洋的凈熱通量; b: 大氣的淡水通量, 海洋中用鹽度通量代替; c: 大氣和海洋的浮力通量, 圖中所示的海洋浮力通量值是實(shí)際值×103

    圖5 不同風(fēng)速下大氣海洋的凈熱通量, 淡水通量(鹽度通量)的關(guān)系圖

    注: 大氣的凈熱通量和淡水通量取10m處的值; 海洋的凈熱通量和鹽度通量取-2m處的值; 實(shí)線: 通量的趨勢(shì)

    2.3 湍流動(dòng)能收支

    圖6展示了海洋湍流動(dòng)能收支情況, 根據(jù)公式(6)本文用紅線表示速度剪切項(xiàng), 綠線表示湍流動(dòng)能傳輸項(xiàng), 藍(lán)線表示壓力項(xiàng), 粉線表示浮力項(xiàng), 黑線表示耗散項(xiàng), 為了便于展示將各項(xiàng)數(shù)值×106。左圖為三個(gè)風(fēng)速下的海洋湍流動(dòng)能收支情況, 右圖均為左圖-35—0m處的局部放大圖。

    從圖6a, 6c, 6e可以發(fā)現(xiàn)海洋湍流動(dòng)能收支各項(xiàng)分別在-15、-25、-35m處迅速減小, 說(shuō)明不同風(fēng)速所能影響海洋的深度隨風(fēng)速增加而加深。在圖6b中, 剪切項(xiàng)與湍流動(dòng)能傳輸項(xiàng)的值在近海面處差別不大, 耗散項(xiàng)的值最大可以到達(dá)-4.5×10-7m2/s3; 在圖6d中, 剪切項(xiàng)與耗散項(xiàng)在近海面處的值達(dá)到最大且兩項(xiàng)大致持平, 最值大約為2×10-6m2/s3; 在圖6f中, 剪切項(xiàng)、傳輸項(xiàng)和耗散項(xiàng)達(dá)到最大值, 最值約為5× 10-6m2/s3。隨著風(fēng)速的增加, 近海面處的剪切與傳輸必然增大, 導(dǎo)致向下傳遞的距離越遠(yuǎn), 所產(chǎn)生的耗散也相應(yīng)增大, 且隨著深度加深逐漸減弱。

    圖6 湍流動(dòng)能收支

    注: a: 5m/s時(shí)海洋湍流動(dòng)能收支, b是a的區(qū)域放大圖; c: 105m/s時(shí)海洋湍流動(dòng)能收支, d是c的區(qū)域放大圖; e: 155m/s時(shí)海洋湍流動(dòng)能收支, f是e的區(qū)域放大圖; TKE: Turbulence Kinetic Energy, 湍流動(dòng)能

    3 討論

    3.1 耦合模式與非耦合模式的差別

    前文分析了耦合模式的結(jié)果, 在耦合模式中海洋是依靠大氣風(fēng)速驅(qū)動(dòng), 而在非耦合模式中海洋是依靠固定的動(dòng)量通量輸入驅(qū)動(dòng)。圖7展示了地轉(zhuǎn)風(fēng)為105m/s情況下, 海氣耦合與非耦合模式的垂向流速在,兩個(gè)方向的時(shí)間序列。圖7a和c分別表示耦合模式下的垂向流速在方向和方向的時(shí)間序列; 圖7b和d分別表示非耦合模式下的垂向流速在方向和方向的時(shí)間序列。

    通過(guò)比較得知: 在整個(gè)模式運(yùn)行期間, 耦合模式海洋部分受風(fēng)速影響較深, 而非耦合模式流速變化巨大的區(qū)域大多集中在靠近海面的地方。

    圖8展示了經(jīng)過(guò)時(shí)間平均后不同深度(-5、-15、-30、-40m)的垂向流速水平截面。從中可以發(fā)現(xiàn)耦合模式與非耦合模式從水下5m到水下30m區(qū)域中的垂向流速分布大致相同, 而在水下40m處時(shí), 非耦合模式的垂向流速減弱明顯, 這也說(shuō)明了在耦合模式下, 受到大氣持續(xù)影響, 能量向下傳遞的距離更遠(yuǎn)。

    圖7 垂向流速時(shí)間序列

    注: a: 耦合模式下方向垂向流速時(shí)間序列; b: 非耦合模式下方向垂向流速時(shí)間序列; c: 耦合模式下方向垂向流速時(shí)間序列; d: 非耦合模式下方向垂向流速時(shí)間序列

    圖8 垂向流速水平截面

    注:: 水平方向坐標(biāo);: 水平方向坐標(biāo);: 垂直方向坐標(biāo)

    圖9展示了耦合模式與非耦合模式的海洋湍流動(dòng)能收支情況, 實(shí)線表示耦合模式, 虛線表示非耦合模式的結(jié)果, 從中可以發(fā)現(xiàn)二者趨勢(shì)大致相同, 但是非耦合模式的湍流動(dòng)能收支各項(xiàng)略小于耦合模式, 且相較于耦合模式在-25m處趨于0, 非耦合模式趨于0的深度較淺, 這也符合圖7和圖8顯示的耦合模式下海洋受風(fēng)速影響區(qū)域較深的結(jié)果。

    3.2 波浪破碎與Langmuir環(huán)流

    海洋混合層受風(fēng)速影響存在兩個(gè)重要湍流過(guò)程: 波浪破碎和Langmuir環(huán)流。波浪破碎在海面附近產(chǎn)生大量的小尺度湍流, 而由風(fēng)驅(qū)動(dòng)的表面流切變與Stokes漂移的相互作用產(chǎn)生Langmuir環(huán)流, 這兩者對(duì)于海洋混合層有較大影響。Noh等(2004)對(duì)海洋混合層進(jìn)行了波浪破碎和Langmuir環(huán)流的大渦模擬, 揭示了波流相互作用在海洋混合層垂向混合過(guò)程中的重要作用。Li等(2013)使用大渦模擬研究?jī)烧咴诤1磉吔鐚又械淖饔? 模式結(jié)果與從CBLAST站點(diǎn)湍流測(cè)量分析得出的結(jié)論一致。波浪破碎與Langmuir環(huán)流均發(fā)生在海氣邊界上, 影響海氣通量交換, 未來(lái)的工作重點(diǎn)便是在耦合模式中加入波浪破碎和Langmuir環(huán)流, 考慮波流相互作用對(duì)海氣通量交換的影響。

    圖9 湍流動(dòng)能收支對(duì)比

    注: a: 耦合模式與非耦合模式的海洋湍流動(dòng)能收支圖, b是a的區(qū)域放大圖, 其中實(shí)線表示耦合模式結(jié)果, 虛線表示非耦合模式結(jié)果

    4 結(jié)論

    本文為探究小尺度海氣相互作用, 使用大渦模擬海氣耦合模式, 設(shè)置了5、10和15m/s三種地轉(zhuǎn)風(fēng)速度, 展示了不同風(fēng)速下海氣通量和OML湍流動(dòng)能收支各項(xiàng)的分布結(jié)果, 并且對(duì)比分析了耦合模式與非耦合模式的差別。主要結(jié)論如下:

    (1) 隨著風(fēng)速的增加, 垂向動(dòng)能在大氣和海洋中的分布呈現(xiàn)縱向加深的趨勢(shì); 凈熱通量在海氣界面處進(jìn)行交換, 在海洋中隨著深度的增加而減小, 受到風(fēng)速的影響, 隨風(fēng)速的增大而增加。海洋受熱加劇海水蒸發(fā), 從而持續(xù)向大氣輸送水汽導(dǎo)致海洋鹽度通量的增大。大氣的浮力通量在100—300m隨風(fēng)速增大而減小, 而海洋的浮力通量隨風(fēng)速增大而增加。

    (2) 海洋混合層中湍流動(dòng)能收支也響應(yīng)風(fēng)速的變化, 其中剪切項(xiàng)、耗散項(xiàng)和傳輸項(xiàng)隨著風(fēng)速的增大而增加, 且在區(qū)域-35—0m變化明顯, 由于海洋模式受到大氣風(fēng)速的驅(qū)動(dòng), 所以在近海表面各項(xiàng)達(dá)到最大值, 之后隨著深度逐漸減弱。

    (3) 非耦合模式中海洋依靠定值的動(dòng)量通量輸入驅(qū)動(dòng), 從而導(dǎo)致非耦合模式中流速在靠近海面的區(qū)域變化明顯, 之后向下迅速減弱; 而耦合模式中海洋依靠大氣模式的風(fēng)速持續(xù)驅(qū)動(dòng), 所以其流速分布顯示在海表面處較大, 而隨著深度緩慢減弱。湍流動(dòng)能收支的結(jié)果表明耦合模式各項(xiàng)略大于非耦合模式。

    對(duì)于大尺度的海氣耦合模式不能解釋小尺度海氣相互作用的問(wèn)題, 本文通過(guò)大渦模擬在湍流尺度展示了海氣通量和海洋湍流動(dòng)能收支情況, 為完善海氣耦合模式提供新的方法。隨著現(xiàn)今計(jì)算機(jī)技術(shù)的發(fā)展, 耦合模式所用的計(jì)算資源與非耦合模式相差不大, 而且耦合模式相比非耦合模式, 充分考慮海洋和大氣的協(xié)同作用, 使得結(jié)果更趨近于真實(shí)情況。

    馮興如, 李水清, 尹寶樹(shù), 2018. 海氣動(dòng)量通量研究綜述. 海洋科學(xué), 42(10): 103—109, doi: 10.11759/hykx20180725001

    陳克明, 金向澤, 張學(xué)洪, 1997. 關(guān)于海氣耦合模式氣候漂移及敏感性的一點(diǎn)探討. 海洋學(xué)報(bào), 19(2): 38—51

    陳烈庭, 1991. 海氣相互作用. 地球科學(xué)進(jìn)展, 6(5): 76—77

    Austin T C, Edson J B, McGillis W R, 2002. A network-based telemetry architecture developed for the Martha's Vineyard Coastal Observatory. IEEE Journal of Oceanic Engineering, 27(2): 228—234, doi: 10.1109/ JOE.2002.1002477

    Deardorff J W, 1970. A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers. Journal of Fluid Mechanics, 41(2): 453—480, doi: 10.1017/S002211 2070000691

    Deardorff J W, 1980. Stratocumulus-capped mixed layers derived from a three-dimensional model. Boundary-Layer Meteorology, 18(4): 495—527, doi: 10.1007/BF00119502

    Edson J, Crawford T, Crescenti J, 2007. The coupled boundary layers and air–sea transfer experiment in low winds. Bulletin of the American Meteorological Society, 88(3): 341—356, doi: 10.1175/bams-88-3-341

    Esau I, 2014. Indirect air–sea interactions simulated with a coupled turbulence-resolving model. Ocean Dynamics, 64(5): 689—705, doi: 10.1007/s10236-014-0712-y

    Huq S, De Roo F, Raasch S, 2019. Vertically nested LES for high-resolution simulation of the surface layer in PALM (version 5.0). Geoscientific Model Development, 12(6): 2523—2538, doi: 10.5194/gmd-12-2523-2019

    Li S, Li M, Gerbi G P, 2013. Roles of breaking waves and Langmuir circulation in the surface boundary layer of a coastal ocean. Journal of Geophysical Research: Oceans, 118(10): 5173—5187, doi: 10.1002/jgrc.20387

    Ma H, Wu L X, Li C, 2010. The role of southern high latitude wind stress in global climate. Advances in Atmospheric Sciences, 27(2): 371—381

    Maronga B, Gryschka M, Heinze R, 2015. The parallelized large-eddy simulation model (PALM) version 4.0 for atmospheric and oceanic flows: model formulation, recent developments, and future perspectives. Geoscientific Model Development, 8(8): 2515—2551, doi: 10.5194/gmd-8-2515- 2015

    Neelin J D, Latif M, Allaart M A F, 1992. Tropical air-sea interaction in general circulation models. Climate Dynamics, 7(2):73—104, doi: 10.1007/BF00209610

    Noh Y, Min H S, Raasch S, 2004. Large eddy simulation of the ocean mixed layer: the effects of wave breaking and Langmuir circulation. Journal of Physical Oceanography, 34(4): 720—735, doi: 10.1175/1520-0485(2004)034<0720: LESOTO>2.0.CO; 2

    Peng S Q, Liu D L, Sun Z B, 2012. Recent advances in regional air-sea coupled models. Science China Earth Sciences, 55(9): 1391—1405, doi: 10.1007/s11430-012- 4386-3

    Skyllingstad E D, Smyth W D, Crawford G B, 2000. Resonant wind-driven mixing in the ocean boundary layer. Journal of Physical Oceanography, 30(8): 1866—1890, doi: 10.1175/ 1520-0485(2000)030<1866:RWDMIT>2.0.CO; 2

    Smagorinsky J, 1963. General circulation experiments with the primitive equations: I. The basic experiment. Monthly Weather Review, 91(3): 99—164, doi: 10.1175/1520-0493 (1963)091<0099:GCEWTP>2.3.CO; 2

    Steinhorn I, 1991. Salt flux and evaporation. Journal of Physical Oceanography, 21(11): 1681—1683, doi: 10.1175/1520- 0485(1991)021<1681:SFAE>2.0.CO;2

    Vickers D, Mahrt L, 2006. Evaluation of the air-sea bulk formula and sea-surface temperature variability from observations. Journal of Geophysical Research: Oceans, 111(C5): C05002, doi: 10.1029/2005JC003323

    Zhou T J, Zhang X H, Yu Y Q, 2000. Response of IAP/ LASG GOALS model to the coupling of air-sea fresh water exchange. Advances in Atmospheric Sciences, 17(3): 473—486, doi: 10.1007/s00376-000-0037-2

    Zhou X M, Li S L, Luo F F, 2015. Air-sea coupling enhances the east Asian winter climate response to the Atlantic Multidecadal Oscillation. Advances in Atmospheric Sciences, 32(12): 1647—1659, doi: 10.1007/s00376-015- 5030-x

    LARGE EDDY SIMULATION COUPLING MODEL FOR SMALL-SCALE AIR-SEA INTERACTION

    SUN Dan-Yi, LI Shuang

    (Ocean College, Zhejiang University, Zhoushan 316021, China)

    Atmosphere and ocean are two important factors that affect the Earth's climate system, and their interaction is an important topic for research. Aiming at the problem that the traditional large-scale air-sea coupled model lacks turbulence scale analysis, this study uses the Parallel Large Eddy Simulation Model (PALM) to explore the effects of wind speed on air-sea flux and turbulence kinetic energy budget on a small scale. Three types of geostrophic wind speeds(5, 10, and 15m/s) were set up to simulate the air-sea coupling of the atmospheric boundary layer (ABL) and ocean mixed layer (OML). Study shows that the distribution of air-sea flux is closely related to wind speed. The larger the wind speed, the larger the net heat flux and buoyancy flux. As the temperature rises, the evaporation of seawater intensifies and the freshwater flux of the atmosphere increases. Ocean turbulence kinetic energy budget is greatly affected by wind speed near the sea surface, and gradually weaken as the depth deepens. This study shows the application of the small-scale air-sea coupled model for the study of air-sea flux, which has certain significance for the study of small-scale air-sea interaction.

    large eddy simulation; air-sea flux; air-sea coupling; TKE budget

    * 國(guó)家自然科學(xué)基金項(xiàng)目, 41876003號(hào); 國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目, 2017YFA0604102號(hào); 國(guó)家自然科學(xué)基金項(xiàng)目, 41830533號(hào), 41911530149號(hào)。孫丹譯, 碩士研究生, E-mail: sundanyi@zju.edu.cn

    李 爽, 博士生導(dǎo)師, 副教授, E-mail: lshuang@zju.edu.cn

    2019-12-26,

    2020-01-28

    P731.26

    10.11693/hyhz20191200277

    猜你喜歡
    收支湍流動(dòng)能
    新動(dòng)能,源自創(chuàng)新力
    金橋(2022年4期)2022-05-05 06:09:50
    關(guān)于未納入海關(guān)統(tǒng)計(jì)的貨物貿(mào)易收支統(tǒng)計(jì)研究
    重氣瞬時(shí)泄漏擴(kuò)散的湍流模型驗(yàn)證
    為構(gòu)建開(kāi)放創(chuàng)新新高地增添動(dòng)能
    “金企對(duì)接”轉(zhuǎn)換舊動(dòng)能
    商周刊(2018年15期)2018-07-27 01:41:16
    澎湃新動(dòng)能
    商周刊(2018年11期)2018-06-13 03:41:52
    2016年浙江省一般公共預(yù)算收支決算總表
    “收支兩條線”生與死
    曾被寄予厚望的“收支兩條線”
    “青春期”湍流中的智慧引渡(三)
    国产成人av激情在线播放| √禁漫天堂资源中文www| 18禁美女被吸乳视频| 亚洲精品久久成人aⅴ小说| 亚洲午夜精品一区,二区,三区| 国产成人精品久久二区二区免费| 国产成人啪精品午夜网站| 久久这里只有精品19| 老鸭窝网址在线观看| 黄色视频,在线免费观看| 黄片小视频在线播放| 日韩 欧美 亚洲 中文字幕| 三上悠亚av全集在线观看| 欧美精品亚洲一区二区| 窝窝影院91人妻| 搡老乐熟女国产| 精品无人区乱码1区二区| 午夜精品在线福利| 久久精品成人免费网站| 久久中文字幕人妻熟女| 久久久水蜜桃国产精品网| 久久精品aⅴ一区二区三区四区| 国产精品.久久久| 国产男靠女视频免费网站| 国产无遮挡羞羞视频在线观看| 精品久久久精品久久久| 91老司机精品| 国产精品永久免费网站| 亚洲综合色网址| 极品人妻少妇av视频| 亚洲一卡2卡3卡4卡5卡精品中文| 最近最新免费中文字幕在线| 大型av网站在线播放| 日韩欧美免费精品| 国产精品.久久久| 黄色视频不卡| 亚洲成av片中文字幕在线观看| 两个人免费观看高清视频| 日韩精品免费视频一区二区三区| 巨乳人妻的诱惑在线观看| www.自偷自拍.com| 天天躁夜夜躁狠狠躁躁| 美女视频免费永久观看网站| 巨乳人妻的诱惑在线观看| 免费在线观看视频国产中文字幕亚洲| 久久人人爽av亚洲精品天堂| 国产片内射在线| 黑人猛操日本美女一级片| 熟女少妇亚洲综合色aaa.| 久久草成人影院| 熟女少妇亚洲综合色aaa.| 91麻豆av在线| 国产单亲对白刺激| e午夜精品久久久久久久| 久久这里只有精品19| 性少妇av在线| 精品无人区乱码1区二区| 午夜精品久久久久久毛片777| 超色免费av| 超色免费av| 男人的好看免费观看在线视频 | 精品欧美一区二区三区在线| 国产成人免费观看mmmm| 十分钟在线观看高清视频www| 美女 人体艺术 gogo| 在线观看日韩欧美| 一进一出抽搐动态| 午夜福利,免费看| 手机成人av网站| netflix在线观看网站| 精品卡一卡二卡四卡免费| 99久久99久久久精品蜜桃| 亚洲自偷自拍图片 自拍| 亚洲一区二区三区不卡视频| 老司机深夜福利视频在线观看| 一本综合久久免费| 一本综合久久免费| 高清欧美精品videossex| 欧美乱色亚洲激情| 国产精品国产高清国产av | 国产免费现黄频在线看| 久久久久久久午夜电影 | 岛国毛片在线播放| 久久久久久久午夜电影 | 露出奶头的视频| 99国产综合亚洲精品| 免费女性裸体啪啪无遮挡网站| 久久精品国产亚洲av高清一级| 亚洲aⅴ乱码一区二区在线播放 | 97人妻天天添夜夜摸| 亚洲专区国产一区二区| 精品少妇一区二区三区视频日本电影| 亚洲成人免费电影在线观看| 国产无遮挡羞羞视频在线观看| 中文字幕制服av| 丝瓜视频免费看黄片| 精品亚洲成国产av| 久久久国产欧美日韩av| 亚洲国产看品久久| 另类亚洲欧美激情| 高清视频免费观看一区二区| 国产午夜精品久久久久久| 国产成人精品在线电影| 狠狠婷婷综合久久久久久88av| 69av精品久久久久久| 久久精品成人免费网站| 波多野结衣一区麻豆| 视频区图区小说| 国产精品亚洲av一区麻豆| 国产片内射在线| 一级毛片女人18水好多| 丝袜美足系列| 又黄又粗又硬又大视频| 激情在线观看视频在线高清 | 日本a在线网址| 18在线观看网站| 制服诱惑二区| 国产97色在线日韩免费| 国产精品免费大片| 黑人欧美特级aaaaaa片| 丁香六月欧美| 99久久99久久久精品蜜桃| 手机成人av网站| 精品一区二区三区视频在线观看免费 | 女人高潮潮喷娇喘18禁视频| 在线观看一区二区三区激情| 久久亚洲精品不卡| 久久天堂一区二区三区四区| 十八禁网站免费在线| 午夜福利在线免费观看网站| 欧美激情高清一区二区三区| 久久久久久久久久久久大奶| 久久久国产成人免费| 成人国语在线视频| 国产精品一区二区在线观看99| 婷婷精品国产亚洲av在线 | 亚洲精华国产精华精| 国产欧美日韩一区二区精品| 欧美丝袜亚洲另类 | 丝袜人妻中文字幕| 波多野结衣一区麻豆| av片东京热男人的天堂| 久久草成人影院| 最新美女视频免费是黄的| 99国产精品99久久久久| 亚洲自偷自拍图片 自拍| 水蜜桃什么品种好| 丝袜在线中文字幕| 人人妻,人人澡人人爽秒播| 天天影视国产精品| 国产蜜桃级精品一区二区三区 | 在线播放国产精品三级| 50天的宝宝边吃奶边哭怎么回事| 啦啦啦视频在线资源免费观看| 亚洲va日本ⅴa欧美va伊人久久| 久久国产精品影院| 中文字幕制服av| avwww免费| 91精品三级在线观看| 天天添夜夜摸| 亚洲国产欧美一区二区综合| 在线观看免费午夜福利视频| 成人国语在线视频| 久久热在线av| 91在线观看av| 在线观看免费视频日本深夜| 亚洲专区中文字幕在线| 国产不卡一卡二| 亚洲第一青青草原| 久久人妻福利社区极品人妻图片| 久久精品亚洲av国产电影网| 黄色片一级片一级黄色片| av网站在线播放免费| xxxhd国产人妻xxx| 欧美性长视频在线观看| 日本黄色视频三级网站网址 | 亚洲熟女精品中文字幕| 少妇的丰满在线观看| 国产精品一区二区在线不卡| 国精品久久久久久国模美| 免费在线观看视频国产中文字幕亚洲| 国产亚洲精品久久久久久毛片 | 亚洲色图 男人天堂 中文字幕| 免费不卡黄色视频| 久久午夜综合久久蜜桃| 在线观看舔阴道视频| 国产精品亚洲一级av第二区| 十八禁高潮呻吟视频| a级片在线免费高清观看视频| 黑丝袜美女国产一区| 久久天躁狠狠躁夜夜2o2o| 欧美日韩一级在线毛片| 在线播放国产精品三级| 欧美日韩国产mv在线观看视频| 亚洲五月色婷婷综合| 亚洲伊人色综图| 久久久精品免费免费高清| 久久久国产精品麻豆| 国产成人影院久久av| 99久久精品国产亚洲精品| 国产aⅴ精品一区二区三区波| 精品国产国语对白av| 亚洲三区欧美一区| 国产成人系列免费观看| 国产欧美亚洲国产| 亚洲五月天丁香| 成人国语在线视频| 午夜免费鲁丝| 老司机在亚洲福利影院| 人人妻,人人澡人人爽秒播| 中文字幕人妻丝袜一区二区| 高清黄色对白视频在线免费看| 视频区图区小说| 国产aⅴ精品一区二区三区波| 欧美人与性动交α欧美软件| 国产精品久久久久久精品古装| 国产xxxxx性猛交| 90打野战视频偷拍视频| 亚洲人成电影观看| 人妻久久中文字幕网| 国产1区2区3区精品| 午夜两性在线视频| 新久久久久国产一级毛片| 日韩三级视频一区二区三区| 麻豆av在线久日| 国产成人av教育| 精品久久久久久久久久免费视频 | 亚洲一区高清亚洲精品| 国产精品一区二区在线不卡| 精品国产国语对白av| 身体一侧抽搐| 久久精品国产亚洲av香蕉五月 | 黄色怎么调成土黄色| 国产精品电影一区二区三区 | 久久精品91无色码中文字幕| 日韩有码中文字幕| xxx96com| 国产亚洲av高清不卡| 精品久久久久久,| 十八禁网站免费在线| 天天操日日干夜夜撸| aaaaa片日本免费| 国产在线精品亚洲第一网站| 老熟妇乱子伦视频在线观看| 国产精品久久久av美女十八| 叶爱在线成人免费视频播放| 久9热在线精品视频| 国产高清videossex| 中文字幕精品免费在线观看视频| 欧美 日韩 精品 国产| 国产精品香港三级国产av潘金莲| 天天添夜夜摸| 极品人妻少妇av视频| 欧美最黄视频在线播放免费 | 制服人妻中文乱码| 亚洲欧美精品综合一区二区三区| 亚洲国产毛片av蜜桃av| 成人亚洲精品一区在线观看| 亚洲人成电影免费在线| 国产日韩一区二区三区精品不卡| 欧美日韩av久久| 国产xxxxx性猛交| 久久中文字幕一级| 日本一区二区免费在线视频| 我的亚洲天堂| 黄色成人免费大全| 国产又爽黄色视频| 久久中文字幕人妻熟女| 国产欧美日韩综合在线一区二区| 精品国产乱子伦一区二区三区| 9热在线视频观看99| 国产真人三级小视频在线观看| 精品国产乱码久久久久久男人| 免费在线观看亚洲国产| 极品教师在线免费播放| 久久精品亚洲av国产电影网| 黄色丝袜av网址大全| 国产精品一区二区精品视频观看| 搡老岳熟女国产| 国产精品av久久久久免费| 国产极品粉嫩免费观看在线| 黄色毛片三级朝国网站| 飞空精品影院首页| 1024香蕉在线观看| 久久久久国内视频| 多毛熟女@视频| 久久久久国产精品人妻aⅴ院 | 男人操女人黄网站| 我的亚洲天堂| 精品国产超薄肉色丝袜足j| 亚洲视频免费观看视频| 久久精品成人免费网站| 欧美中文综合在线视频| 欧美日韩精品网址| 日本欧美视频一区| 91av网站免费观看| 可以免费在线观看a视频的电影网站| 午夜91福利影院| 亚洲专区中文字幕在线| 大香蕉久久网| 国产欧美日韩综合在线一区二区| 欧美最黄视频在线播放免费 | 怎么达到女性高潮| 女人被躁到高潮嗷嗷叫费观| 午夜福利,免费看| 欧美日韩成人在线一区二区| 成人18禁高潮啪啪吃奶动态图| 午夜激情av网站| 久久精品熟女亚洲av麻豆精品| 欧美在线一区亚洲| 免费久久久久久久精品成人欧美视频| 亚洲色图av天堂| 久久精品aⅴ一区二区三区四区| 人成视频在线观看免费观看| 午夜久久久在线观看| 欧美激情久久久久久爽电影 | 精品久久蜜臀av无| 又黄又粗又硬又大视频| 亚洲国产毛片av蜜桃av| 最近最新中文字幕大全免费视频| 欧美成人免费av一区二区三区 | 三上悠亚av全集在线观看| 免费高清在线观看日韩| 日本vs欧美在线观看视频| 欧美人与性动交α欧美软件| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人精品无人区| 美女扒开内裤让男人捅视频| 老司机在亚洲福利影院| 久久天堂一区二区三区四区| 99riav亚洲国产免费| 久久国产精品人妻蜜桃| 一区二区三区激情视频| 亚洲精品美女久久久久99蜜臀| 欧美色视频一区免费| 国产成人免费无遮挡视频| 国产精品影院久久| 午夜亚洲福利在线播放| 亚洲精品国产精品久久久不卡| 男人操女人黄网站| 成人手机av| 久久ye,这里只有精品| 激情视频va一区二区三区| 精品人妻熟女毛片av久久网站| 久久香蕉激情| 亚洲av成人一区二区三| 美女 人体艺术 gogo| 咕卡用的链子| 搡老岳熟女国产| 老鸭窝网址在线观看| 成人影院久久| 身体一侧抽搐| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲一区二区三区欧美精品| 亚洲avbb在线观看| 亚洲精品一卡2卡三卡4卡5卡| 少妇被粗大的猛进出69影院| 日韩大码丰满熟妇| 夜夜躁狠狠躁天天躁| 亚洲 国产 在线| 精品国产乱码久久久久久男人| 老司机福利观看| 麻豆成人av在线观看| 精品乱码久久久久久99久播| 久久香蕉激情| 欧美乱色亚洲激情| 午夜福利一区二区在线看| 久久精品亚洲av国产电影网| 18禁裸乳无遮挡免费网站照片 | 国产欧美日韩一区二区三| 欧美成狂野欧美在线观看| 看免费av毛片| 国产精品影院久久| 黄频高清免费视频| √禁漫天堂资源中文www| 女同久久另类99精品国产91| 久久香蕉国产精品| 婷婷丁香在线五月| 变态另类成人亚洲欧美熟女 | 久久热在线av| 亚洲专区中文字幕在线| 极品教师在线免费播放| 免费久久久久久久精品成人欧美视频| 黄色成人免费大全| 欧美日本中文国产一区发布| 国产成人免费观看mmmm| 手机成人av网站| 久久中文看片网| 久久久精品国产亚洲av高清涩受| 一本综合久久免费| 丰满的人妻完整版| 女警被强在线播放| 老司机深夜福利视频在线观看| 一区二区日韩欧美中文字幕| 大香蕉久久成人网| 啦啦啦在线免费观看视频4| 涩涩av久久男人的天堂| 国产精品久久久久久人妻精品电影| 久久午夜综合久久蜜桃| 每晚都被弄得嗷嗷叫到高潮| 叶爱在线成人免费视频播放| 高清在线国产一区| 欧美日韩亚洲高清精品| 久久中文字幕一级| 久久久久精品国产欧美久久久| 国产欧美日韩一区二区三区在线| 精品国产乱码久久久久久男人| 日韩欧美在线二视频 | 日韩大码丰满熟妇| 国产1区2区3区精品| 丁香欧美五月| 久久性视频一级片| 日韩免费高清中文字幕av| 另类亚洲欧美激情| 亚洲人成电影免费在线| www日本在线高清视频| 人妻一区二区av| 欧美黄色淫秽网站| 国产视频一区二区在线看| 久久香蕉国产精品| 欧美另类亚洲清纯唯美| 欧美日韩亚洲高清精品| 午夜亚洲福利在线播放| 国产aⅴ精品一区二区三区波| 国产深夜福利视频在线观看| 中亚洲国语对白在线视频| 午夜免费鲁丝| 波多野结衣av一区二区av| 欧美成狂野欧美在线观看| 精品国产国语对白av| 日韩欧美一区二区三区在线观看 | 国产成人啪精品午夜网站| 久久人人97超碰香蕉20202| 交换朋友夫妻互换小说| 窝窝影院91人妻| avwww免费| 99久久国产精品久久久| 精品福利观看| 欧美激情极品国产一区二区三区| 国产男女内射视频| 国产欧美日韩一区二区三区在线| 男女之事视频高清在线观看| 我的亚洲天堂| 午夜福利欧美成人| 欧美色视频一区免费| 亚洲一区中文字幕在线| 久久久国产一区二区| 国产精品一区二区精品视频观看| av有码第一页| 国产精品欧美亚洲77777| 纯流量卡能插随身wifi吗| 视频在线观看一区二区三区| 精品久久久精品久久久| 欧美激情 高清一区二区三区| 99热网站在线观看| 一区二区三区激情视频| 亚洲九九香蕉| 男女高潮啪啪啪动态图| 国产深夜福利视频在线观看| 高清毛片免费观看视频网站 | 亚洲欧美激情综合另类| 又黄又粗又硬又大视频| 国产欧美日韩一区二区三| 中出人妻视频一区二区| 久久久久久免费高清国产稀缺| 日本a在线网址| 免费一级毛片在线播放高清视频 | 精品午夜福利视频在线观看一区| 老司机靠b影院| 精品福利永久在线观看| 亚洲午夜理论影院| 欧美激情 高清一区二区三区| 国产成人免费观看mmmm| 亚洲精品在线美女| 亚洲五月天丁香| 色综合婷婷激情| 99热只有精品国产| 高清av免费在线| 精品人妻1区二区| 中文欧美无线码| a级毛片黄视频| 成人亚洲精品一区在线观看| 成人18禁高潮啪啪吃奶动态图| 日韩成人在线观看一区二区三区| 欧美日韩成人在线一区二区| 久久久久国产精品人妻aⅴ院 | 看免费av毛片| 国产高清激情床上av| 最新在线观看一区二区三区| 国产成人欧美| 亚洲精品自拍成人| 精品久久久精品久久久| 捣出白浆h1v1| 久久久久国产一级毛片高清牌| www.自偷自拍.com| 男人的好看免费观看在线视频 | 高清黄色对白视频在线免费看| 美女福利国产在线| 在线免费观看的www视频| 国产精品亚洲av一区麻豆| 纯流量卡能插随身wifi吗| 高清毛片免费观看视频网站 | 婷婷丁香在线五月| 日韩欧美国产一区二区入口| 91国产中文字幕| 黄色成人免费大全| 国产伦人伦偷精品视频| 黄色毛片三级朝国网站| 亚洲第一青青草原| 欧美成人免费av一区二区三区 | 99在线人妻在线中文字幕 | 美女午夜性视频免费| 日本撒尿小便嘘嘘汇集6| 久久久久久久国产电影| 十八禁人妻一区二区| 成年版毛片免费区| 精品国产乱码久久久久久男人| 国产激情欧美一区二区| 1024香蕉在线观看| 欧美国产精品一级二级三级| 美女午夜性视频免费| 成人国语在线视频| av天堂在线播放| 波多野结衣一区麻豆| 亚洲精品在线美女| 久久精品成人免费网站| 在线十欧美十亚洲十日本专区| 狠狠婷婷综合久久久久久88av| 国产1区2区3区精品| 亚洲精品美女久久av网站| 精品第一国产精品| 日韩欧美一区视频在线观看| 91麻豆av在线| 午夜影院日韩av| 国产麻豆69| 999久久久国产精品视频| 久久香蕉激情| 天堂俺去俺来也www色官网| 老鸭窝网址在线观看| 亚洲色图av天堂| 最近最新中文字幕大全免费视频| 韩国av一区二区三区四区| 岛国毛片在线播放| 亚洲色图综合在线观看| 老汉色av国产亚洲站长工具| 国产在线精品亚洲第一网站| 久久精品国产清高在天天线| 一级片'在线观看视频| 国产真人三级小视频在线观看| 欧美丝袜亚洲另类 | 精品午夜福利视频在线观看一区| 国产日韩欧美亚洲二区| 亚洲人成伊人成综合网2020| 国产欧美日韩一区二区三区在线| 如日韩欧美国产精品一区二区三区| 国产亚洲精品久久久久久毛片 | 亚洲国产欧美网| 成人免费观看视频高清| 看片在线看免费视频| 国产一区二区三区综合在线观看| 淫妇啪啪啪对白视频| 精品国产美女av久久久久小说| 国内久久婷婷六月综合欲色啪| 国产又色又爽无遮挡免费看| 自拍欧美九色日韩亚洲蝌蚪91| 国产淫语在线视频| 国产精品久久电影中文字幕 | 国产精品自产拍在线观看55亚洲 | 岛国毛片在线播放| 久久青草综合色| 一本一本久久a久久精品综合妖精| 久久久久精品国产欧美久久久| 在线观看免费视频网站a站| 欧美中文综合在线视频| 午夜福利乱码中文字幕| 午夜成年电影在线免费观看| 久久狼人影院| 老司机福利观看| avwww免费| 免费av中文字幕在线| av视频免费观看在线观看| 伦理电影免费视频| 国产精品免费一区二区三区在线 | 欧美精品啪啪一区二区三区| 男女下面插进去视频免费观看| 热99国产精品久久久久久7| 狠狠狠狠99中文字幕| 一边摸一边抽搐一进一出视频| 高清黄色对白视频在线免费看| 激情在线观看视频在线高清 | 亚洲专区字幕在线| 久久人妻福利社区极品人妻图片| 美国免费a级毛片| 天堂中文最新版在线下载| 超碰成人久久| 久久99一区二区三区| 亚洲国产精品合色在线| 夜夜爽天天搞| 在线观看免费视频日本深夜| 久久精品亚洲av国产电影网| 久久婷婷成人综合色麻豆| 欧美日韩一级在线毛片| 欧美+亚洲+日韩+国产| 国产麻豆69| 国产亚洲精品第一综合不卡| 18禁裸乳无遮挡动漫免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人免费av一区二区三区 | 欧美日韩乱码在线| 色婷婷久久久亚洲欧美| 日韩人妻精品一区2区三区| 亚洲中文字幕日韩| 欧美日韩福利视频一区二区| 99久久人妻综合| 国产精品亚洲av一区麻豆|