徐巖, 劉香全, 3, 4, 宋仁剛, 岑顯榮, 郭雙喜, 周生啟

?

擴散對流溫鹽臺階結構的數(shù)值模擬

徐巖1, 劉香全1, 3, 4, 宋仁剛1, 岑顯榮2, 郭雙喜2, 周生啟2

1. 山東科技大學 電子通信與物理學院, 山東 青島 266590; 2. 熱帶海洋環(huán)境國家重點實驗室(中國科學院南海海洋研究所), 廣東 廣州 510301; 3. 中國科學院半導體研究所 集成光電子學國家重點實驗室, 北京 100083; 4. 中國科學院大學 材料科學與光電技術學院, 北京 100049

熱鹽驅動下的擴散對流現(xiàn)象是海洋中高緯度海域普遍存在的一種現(xiàn)象, 對其進行數(shù)值模擬可更細致地研究海洋小尺度動力過程。文章分析了擴散對流的形成機制, 建立了二維方腔模型, 通過有限體積法求解控制方程, 對其分層現(xiàn)象進行了數(shù)值模擬。研究給出了流場的溫度及鹽度隨時間演化的關系, 展現(xiàn)了流場中速度的渦旋結構, 分析了溫鹽臺階結構的生成、合并的演化過程, 并對其物質和能量的輸運進行了初步的理論解釋。另外, 對不同熱流密度情況下的擴散對流現(xiàn)象進行了對比研究, 發(fā)現(xiàn)隨著熱流密度的增加, 臺階結構的演變速率變快, 而且上邊界冷卻對其演化速率具有促進作用, 但熱流密度的改變并沒有對臺階結構的演變趨勢產(chǎn)生明顯的影響。

擴散對流; 溫鹽臺階; 鹽度; 溫度

在自然界及工業(yè)生產(chǎn)中存在著諸多的雙擴散現(xiàn)象, 如海洋中受溫鹽共同作用的擴散對流及鹽指現(xiàn)象, 地殼中軟流層的運動, 大氣中熱和污染物的擴散, 工業(yè)生產(chǎn)中晶體的成長、金屬的固化等(鄭珺婷等, 2002)。當流體中的兩種(或多種)組分具有不同的擴散率時, 容易產(chǎn)生雙擴散現(xiàn)象。在海洋環(huán)境中, 海水的熱擴散率(≈1.4×10-7m2·s-1)比鹽擴散率(≈1.4×10-9m2·s-1)大兩個量級, 且當兩者對海水的密度梯度的貢獻作用相反時, 在熱驅動和鹽驅動下形成雙擴散對流, 從而促進海洋內部混合, 并進一步影響著海洋中物質和能量的輸運(Schmitt, 1994)。

在雙擴散的室內實驗中, 線性層結的鹽水系統(tǒng)在底部加熱、頂部冷卻時, 會產(chǎn)生明顯的溫鹽臺階結構, 形成“擴散成層”現(xiàn)象(Turner, 1968; Huppert et al, 1979; Lewis et al, 1982; 萬偉等, 2014; Guo et al, 2016)。在海洋中, 當?shù)蜏氐望}的水位于高溫高鹽的水之上時, 就會產(chǎn)生擴散對流現(xiàn)象, 在溫鹽廓線中呈現(xiàn)臺階結構(鮑征宇, 1989; Guthrie et al, 2015; Kimura et al, 2015; Xie et al, 2015; Zhou et al, 2016; Bebieva et al, 2016; Shibley et al, 2017; Guthrie et al, 2017), 如圖1a所示。具體形成機制如下: 如圖1b所示, 在界面之下的高溫高鹽海水受到擾動, 從暖而咸的深層環(huán)境進入到界面之上的冷且淡的淺層環(huán)境, 由于熱擴散速率較大, 該海水微團就會迅速地向周圍環(huán)境釋放熱量, 直至溫度降低到與周圍海水微團溫度相同。但是此時由于鹽擴散速率較小, 使得該海水微團基本保持原來的鹽度不變, 深層的海水微團上升一段距離后, 與周圍海水微團密度一致; 若上升的距離較大, 使該海水微團的密度大于周圍海水微團, 則海水微團在重力的作用下向下運動, 最終達到比較穩(wěn)定的狀態(tài), 形成有限高度且相對均勻的對流層。類似地, 界面上的低溫低鹽的海水微團受到擾動, 從冷而淡的淺層環(huán)境進入界面之下的暖而咸的深層環(huán)境, 由于受熱而密度不斷減小, 在下降一段距離后, 密度與周圍海水一致, 之后繼續(xù)吸收熱量而密度重新變小, 微團又開始向上運動。如此反復, 在界面的上下形成了對流混合層?；旌蠈觾鹊臏囟取Ⅺ}度和密度均勻, 以對流運動為主; 界面處則具有較大的梯度, 以分子擴散運動為主, 最終形成臺階結構。

圖1 擴散對流

a. 北冰洋的溫度臺階結構(Neal et al, 1969)。圖中的矩形陰影表示需要放大后才能觀察到溫度臺階結構的剖面, 圓環(huán)表示溫度臺階界面的雙層結構。圖的左下橫軸(雙箭頭線)的長度表示溫度的漲落幅度; b. 擴散對流形成機制示意圖

Fig. 1 Thermal staircases in the Arctic Ocean (a), and illustration of formation mechanism of diffusive convection (b)

在對溫鹽臺階問題的研究中, 人們采用了多種方法, 如流體力學室內實驗、海洋現(xiàn)場觀測以及數(shù)值模擬等。Turner等(1964)通過實驗演示了溫鹽臺階結構的生成過程, 并對其進行了定性的描述。1968年, Turner (1968)通過實驗, 發(fā)現(xiàn)了第二個混合層出現(xiàn)之前, 第一個混合層的厚度和第一個擴散界面的溫差、鹽度差隨時間的變化與1/2的標度律相符合, 且對第二個混合層出現(xiàn)的臨界條件進行了研究。1972年, Huppert等(1972)的研究表明, Turner在室內實驗進行的溫鹽臺階結構研究得到的熱通量經(jīng)驗公式可應用于海洋中大尺度流動的研究。這些研究確立了擴散對流溫鹽臺階熱通量的基本規(guī)律, 即和溫差的4/3標度律。Marmorino等(1976)以及Kelley (1990)根據(jù)室內實驗數(shù)據(jù), 指出每個臺階都可以理解為一個包含對流環(huán)的對流層, 并確定了熱通量和密度比的關系, 進而完善了熱通量參數(shù)化方案。Carpenter等(2012)利用直接數(shù)值模擬方法研究了擴散對流界面的動力過程, 發(fā)現(xiàn)界面的熱輸運以分子熱擴散為主。Guthrie等(2015)利用北冰洋擴散對流的觀測數(shù)據(jù)對Kelley(1990)的參數(shù)化方案進行了驗證。以上的這些工作主要關注擴散對流的界面穩(wěn)定性及熱通量問題, 而對于臺階演化的完整過程還有待做進一步的細致研究。

本文在前人的基礎上, 通過數(shù)值方法模擬了擴散對流臺階結構的生成、演化、合并以及最后消失的整個過程。通過設置不同邊界條件, 來討論熱流密度對擴散對流臺階結構演變速率的影響規(guī)律。

1 模型建立

1.1 物理模型

用數(shù)值計算的手段對擴散對流進行模擬具有多種優(yōu)勢, 能夠精確地描述流場及標量場, 且嚴格控制邊界條件(邢元明等, 2013)。Turner(1965)指出, 雙擴散的穩(wěn)定性可以用密度比進行衡量, 定義如下:

式(1)中,為溫度,為鹽度,為流體密度,=-(1/)?/?和=-(1/)?/?分別為熱膨脹系數(shù)和鹽收縮系數(shù),S=?/?和T=?/?分別為鹽度和溫度的垂向梯度。定義為流體鹽擴散率與熱擴散率之比, 即=, 普朗特數(shù)=/為流體黏性與熱擴散率之比。線性穩(wěn)定理論表明, 1<R<(+1)/(+)是雙擴散不穩(wěn)定性發(fā)生的條件。在海洋中, Kelley等(2003)認為當密度比R在1~10范圍內時, 溫鹽臺階結構更容易產(chǎn)生, 并且當它越接近于1時臺階結構越明顯(Kelley et al, 2003; 屈玲, 2014)。

到目前為止, 對雙擴散的理論研究, Boussinesq模型得到了廣泛的認可。Boussinesq假設認為密度的變化并不顯著改變流體的性質(賴錫軍, 2004); 同時，在動量方程中, 密度的變化對慣性力項、壓力差項和黏性力項的影響可忽略不計, 而僅考慮其對質量力項的影響。假設流體為不可壓縮流體, 密度的線性狀態(tài)方程為:

表達式(2)中,0、0、0分別代表參考溫度、參考鹽度及參考密度。

在Boussinesq假設下, 雙擴散系統(tǒng)的控制方程為:

1.2 數(shù)值方法

建立擴散對流的二維方腔模型, 如圖2所示。設定其長和寬均為10cm, 底部加熱, 頂部冷卻(即底部和頂部施加固定熱流密度邊界條件), 左右壁面為絕熱壁面, 滿足?/?=0, 各邊界為無滑移邊界, 且滿足?/?=0 (詹杰民等, 2002; 盧玉華等, 2006)?？紤]重力的作用, 其方向為沿軸的負方向。設定初始鹽度隨高度線性遞減, 初始溫度則均勻分布, 由此形成線性層結的初始密度分布。

圖2 擴散對流數(shù)值模型

控制方程的離散采用有限體積法, 為了在保證計算穩(wěn)定性的前提下盡量提高計算精度, 本文在計算時對于擴散項采用一階迎風格式, 而對流項則采用中心差分格式; 在時間離散方面, 傳統(tǒng)的三時間層蛙跳格式在氣象和海洋數(shù)值預報中使用較廣泛, 但存在計算解及時間濾波的精度降階問題, 為彌補此不足, 可采用回溯時間差分格式(Feng et al, 2001)或Adams-Bashforth格式(Durran, 1991)等替代方案。本文采用Adams-Bashforth格式, 相比蛙跳格式它能更大程度地保持方程組的計算準確性。計算在方形網(wǎng)格上進行, 網(wǎng)格數(shù)為100×100, 時間步長設為0.01s。在擴散對流模擬中, 設二維方腔中的初始溫度為參考溫度(300K), 頂部的初始鹽度為0, 底部的初始鹽度為34‰。

用Paraview軟件對結果進行后處理時, 鹽度在0~34‰范圍內等間距的取50條等值線, 繪制中截線(=5cm)上的鹽度變化曲線。

2 結果與討論

2.1 臺階結構的演化過程

當加熱和冷卻的熱流密度均為2200W·m-2時, 觀察二維方腔系統(tǒng)鹽度的變化, 情況如下(圖3):

當時間為0.169h時, 雙擴散不穩(wěn)定性開始發(fā)生, 在方腔的頂部和底部出現(xiàn)了羽流結構; 時間為0.433h時, 頂部和底部已出現(xiàn)了明顯的臺階結構, 在高度約為1.1cm、1.9cm、8.2cm、8.9cm處出現(xiàn)4個較為明顯的擴散界面, 1.9~8.2cm的中間部分鹽度仍呈線性分布。這些臺階結構基本是關于=5cm這條直線對稱的, 這是由于模型的邊界設置的對稱性所致。

當時間約為0.839h時, 在高度約為1.7cm、2.3cm、3.2cm、6.9cm、7.7cm、8.4cm位置處出現(xiàn)了6個擴散界面, 較之前相比, 兩個主穩(wěn)定層(最底和最高的均勻層)的厚度增加了, 而兩個次穩(wěn)定層(頂端主穩(wěn)定層的下端, 底端主穩(wěn)定層的上端)的厚度減小了。

當時間為0.989h時, 在高度約為1.8cm、2.4cm、3.2cm、3.9cm、6.2cm、6.7cm、7.5cm、8.1cm位置處出現(xiàn)了8個擴散界面, 3.9~6.2cm的中間部分鹽度仍呈線性分布, 但是臨近界面處有了波動。

圖3 擴散對流模擬不同特征時刻的鹽度分布圖

當時間為1.081h時, 在高度約為2.0cm、3.2cm、4.1cm、4.5cm、6.0cm、6.6cm、8.0cm位置處出現(xiàn)了7個擴散界面; 當時間為1.392h時, 在高度約為2.2cm、3.4cm、4.3cm、5.4cm、6.5cm、7.8cm位置處出現(xiàn)了6個擴散界面, 存在著7個鹽度相對均勻的階層, 中間均勻層的厚度基本相同, 如圖4所示。在接下來約0.5h內, 擴散對流保持著相對穩(wěn)定性, 臺階的數(shù)目和厚度基本不變。隨著主穩(wěn)定層的不斷增厚, 兩個次穩(wěn)定層的厚度不斷減小。

當時間約為1.989h時, 上下端的次穩(wěn)定層消失, 只剩下5個相對穩(wěn)定的階層。當時間約為2.747h時, 下端新的次穩(wěn)定層消失, 剩下4個相對穩(wěn)定的階層, 如圖5所示。當時間為3.858h時, 上下兩端的主穩(wěn)定層向中間擴展, 變成3個相對穩(wěn)定的階層。中間穩(wěn)定的階層不斷縮小。當時間約為4.811h時, 中間均勻層消失, 變成2個相對穩(wěn)定的階層, 如圖6所示。然后擴散界面的位置基本保持不變, 上下兩層的鹽度差逐漸減小。當鹽度差減小到一定程度時, 交界面向上移動, 鹽度差依然在減小, 當時間約為17.106h時, 只存在一個穩(wěn)定的均勻層, 整個體系的鹽度為17‰。Carpenter等(2012)在兩層“耗盡” (“run-down”)模擬實驗中也關注了擴散對流的演化過程: 在“耗盡”實驗中, 上下層臺階是自由發(fā)展的,期間擴散界面的溫差和鹽度差均逐漸減小, 密度比不斷增大。相比之下, 本文模擬的擴散對流的演化經(jīng)歷了主穩(wěn)定層的不斷擴展和次穩(wěn)定層的不斷消失, 期間擴散界面的溫差基本維持不變而鹽度差逐漸減小, 密度比不斷減小, 與“耗盡”實驗的演化機制有所區(qū)別。

圖4 擴散對流模擬1.392h后鹽度狀態(tài)圖

a. 等值線圖; b. 曲線圖

Fig. 4 Salinity distribution after 1.392 h in diffusive convection simulation: (a) contour plot, and (b) line plot

圖5 擴散對流模擬2.747h后鹽度狀態(tài)圖

a. 等值線圖; b. 曲線圖

Fig. 5 Same as Fig. 4, except after 2.747 h

圖6 擴散對流模擬4.811h后鹽度狀態(tài)圖

a. 等值線圖; b. 曲線圖

Fig. 6 Same as Fig. 4, except after 4.811 h

在模擬過程中, 觀察二維方腔系統(tǒng)溫度的變化, 發(fā)現(xiàn)溫度曲線的變化趨勢與鹽度曲線的趨勢基本是一致的。萬偉等(2014)在擴散對流的室內實驗中也觀察到了類似的溫度臺階演化過程, 但由于缺少速度測量, 無法展現(xiàn)擴散對流的渦旋結構, 而本文的模擬工作可以彌補這些不足。觀察二維方腔系統(tǒng)流場速度的變化, 流速場中不斷伴隨著渦旋的生成與消亡。選取4.811h時的速度分布矢量圖, 如圖7a所示, 在擴散交界面處流動較為微弱, 說明在擴散界面內部主要以分子擴散的形式輸運物質, 以分子熱傳導的形式輸運熱量; 而在均勻的混合層內部渦旋較為明顯, 說明主要以對流的形式輸運物質和能量。選取17.106h時的速度矢量圖, 如圖7b所示, 此時整個系統(tǒng)的溫度和鹽度達到均勻穩(wěn)定狀態(tài), 形成了一個大尺度的環(huán)流。

圖7 流場速度分布矢量圖

a. 4.811 h時; b. 17.106 h時

Fig. 7 Vector plot of simulated velocity field: (a) 4.811 h, and (b) 17.106 h

2.2 不同熱流密度下對比討論

本模擬在對二維方腔上下邊界的加熱和冷卻設置時, 選取了不同的熱流密度。具體設置情況如下。

情況①: 上下熱流密度均為1000W·m–2;

情況②: 上下熱流密度均為2200W·m–2;

情況③: 上下熱流密度均為3000W·m–2;

情況④: 上冷卻熱流密度為1000W·m–2, 下加熱熱流密度為2000W·m–2。

我們關注不同工況下溫鹽臺階的數(shù)目隨時間的變化規(guī)律, 模擬結果如圖8所示。

圖8 不同熱流密度下擴散對流溫鹽臺階結構模擬結果對比

情況①—③中的上下邊界熱流密度在數(shù)值上具有對稱性。在圖8中可以看出, 與情況②的模擬結果相比, 情況①出現(xiàn)相同擴散界面的時間明顯增加, 說明臺階結構的生成及合并速率明顯減慢; 情況③的臺階結構的生成及合并速率明顯加快, 演變時間明顯變短。由擴散界面的變化趨勢可以看出, 在這3種情況下, 臺階結構的演變過程是類似的, 前期是臺階結構的增加, 后期是臺階結構的合并, 但臺階的數(shù)目及速率存在差異。

情況④中, 由于下邊界的熱流密度值比上邊界更大, 擴散對流溫鹽臺階結構的對稱性被打破。情況④與情況②相比, 在下邊界加熱熱流密度接近, 上邊界冷卻熱流密度差一倍的情況下, 情況②的曲線位于情況④曲線的左側, 如圖8所示, 表明情況②相同擴散界面出現(xiàn)的時間短, 演變速率較快。該現(xiàn)象說明了情況②中上邊界較大的冷卻熱流密度對溫鹽臺階結構的演化具有一定的促進作用。但是這兩種情況擴散界面數(shù)隨時間的演化趨勢是類似的, 并沒有明顯的不同。另外, 情況④在臺階演變過程中, 下方主穩(wěn)定層的厚度一直大于上方主穩(wěn)定層的厚度, 臺階結構不再是上下對稱的形式, 如圖9所示, 該現(xiàn)象是由于上下邊界熱流密度不同所致, 下邊界較大的熱流密度使得下方主穩(wěn)定層具有較快的演變速率。

圖9 情況④擴散對流模擬2.764h后鹽度狀態(tài)圖

a. 等值線圖; b. 曲線圖

Fig. 9 Salinity distribution after 2.764 h in diffusive convection simulation of case 4: (a) contour plot, and (b) line plot

3 結論

本文對擴散對流現(xiàn)象進行了數(shù)值模擬, 模擬時采用線性層結的鹽水系統(tǒng), 底部加熱, 頂部冷卻。在擴散對流溫鹽臺階的演化過程中, 溫度和鹽度的變化趨勢基本一致, 存在著溫鹽臺階的生成與合并。約在模擬開始的1h之前, 臺階的生成速率要大于臺階的合并速率, 表現(xiàn)為臺階數(shù)目的增多; 之后達到一個相對穩(wěn)定的階段, 接下來臺階結構不再生成, 只存在合并過程。約在1.5h之后, 次穩(wěn)定層厚度逐漸減小, 直到消失, 在合并的過程中, 歷時較長, 最終整個系統(tǒng)達到一個溫度和鹽度相對均勻穩(wěn)定的狀態(tài), 形成一個大尺度的環(huán)流。

在整個溫鹽臺階結構的演變過程中, 溫度和鹽度曲線的變化趨勢基本是一致的。對于流場中物質和能量的輸運而言, 在擴散界面內部, 主要以分子擴散的形式輸運物質, 以分子熱傳導的形式輸運熱量; 而在均勻的混合層內部主要以對流的形式輸運物質和能量。

另外, 在不同熱流密度情況下, 對擴散對流溫鹽臺階結構的演變過程進行了對比討論, 分析了熱流密度對溫鹽臺階結構形成的影響。模擬結果表明: 熱流密度的變化主要是改變了溫鹽臺階結構的演變速率, 并不影響演變的趨勢, 底部的加熱和頂部的冷卻都可以有效的加快其演變速率。

鮑征宇, 1989. 雙擴散對流理論簡介及國外研究現(xiàn)狀[J]. 地質科技情報, 8(4): 35–42. BAO ZHENGYU, 1989. Brief double-diffusive convection theory and its recent development abroad[J]. Geological Science and Technology Information, 8(4): 35–42 (in Chinese with English abstract).

賴錫軍, 2004. 溫鹽分層流的非結構網(wǎng)格數(shù)值模擬[D]. 南京: 河海大學. LAI XIJUN, 2004. Numerical simulation of stratified flows with combined temperature and salinity on unstructured grids[D]. Nanjing: Hohai University (in Chinese with English abstract).

盧玉華, 詹杰民, 2006. 三維方腔溫鹽雙擴散的格子Boltzmann方法數(shù)值模擬[J]. 物理學報, 55(9): 4774–4782. LU YUHUA, ZHAN JIEMIN, 2006. Three-dimensional numerical simulation of thermosolutal convection in enclosures using lattice Boltzmann method[J]. Acta Physica Sinica, 55(9): 4774–4782 (in Chinese with English abstract).

屈玲, 2014. 海洋擴散對流溫鹽臺階機理及演化的研究[D]. 北京: 中國科學院大學. QU LING, 2014. The fundamental study of diffusive convection staircase and its evolution in the ocean[D]. Beijing: University of Chinese Academy of Sciences (in Chinese with English abstract).

萬偉, 屈玲, 周生啟, 2014. 雙擴散對流中臺階結構的實驗研究[J]. 力學學報, 46(2): 217–223. WAN WEI, QU LING, ZHOU SHENGQI, 2014. Laboratory studies on the staircase structure of double-diffusive convection[J]. Chinese Journal of Theoretical and Applied Mechanics, 46(2): 217–223 (in Chinese with English abstract).

邢元明, 楊磊, 管玉平, 2013. 海表擴散層中氣體行為的分子動力學模擬[J]. 熱帶海洋學報, 32(2): 82–87. XING YUANMING, YANG LEI, GUAN YUPING, 2013. Molecular dynamics simulation of gas behaviors in seawater diffusive layer beneath the air-sea interface[J]. Journal of Tropical Oceanography, 32(2): 82–87 (in Chinese with English abstract).

詹杰民, 李毓湘, 2002. 溫鹽雙擴散均衡場中的振蕩現(xiàn)象[J]. 物理學報, 51(4): 828–834. ZHAN JIEMIN, LI YUXIANG, 2002. Oscillation phenomena in a thermohaline double-diffusive convection system[J]. Acta Physica Sinica, 51(4): 828–834 (in Chinese with English abstract).

鄭珺婷, 詹杰民, 2002. 雙擴散對流系統(tǒng)問題的研究進展[J]. 力學進展, 32(3): 415–424. ZHENG JUNTING, ZHAN JIEMIN, 2002. A review on double-diffusive convection[J]. Advances in Mechanics, 32(3): 415–424 (in Chinese with English abstract).

BEBIEVA Y, TIMMERMANS M L, 2016. An examination of double‐diffusive processes in a mesoscale eddy in the Arctic Ocean[J]. Journal of Geophysical Research: Oceans, 121(1): 457–475.

CARPENTER J R, SOMMER T, WüEST A, 2012. Simulations of a double-diffusive interface in the diffusive convection regime[J]. Journal of Fluid Mechanics, 711: 411–436.

DURRAN D R, 1991. The third-order Adams-Bashforth method: An attractive alternative to leapfrog time differencing[J]. Monthly Weather Review, 119(3): 702–720.

FENG GUOLIN, CAO HONGXING, DONG WENJIE, et al, 2001. A new difference scheme with multi-time levels[J]. Chinese Physics, 10(11): 1004–1010.

GUO SHUANGXI, ZHOU SHENGQI, QU LING, et al, 2016. Laboratory experiments on diffusive convection layer thickness and its oceanographic implications[J]. Journal of Geophysical Research: Oceans, 121(10): 7517–7529.

GUTHRIE J D, FER I, MORISON J, 2015. Observational validation of the diffusive convection flux laws in the Amundsen Basin, Arctic Ocean[J]. Journal of Geophysical Research: Oceans, 120(12): 7880–7896.

GUTHRIE J D, FER I, MORISON J H, 2017. Thermohaline staircases in the Amundsen Basin: Possible disruption by shear and mixing[J]. Journal of Geophysical Research: Oceans, 122(10): 7767–7782.

HUPPERT H E, TURNER J S, 1972. Double-diffusive convection and its implications for the temperature and salinity structure of the ocean and Lake Vanda[J]. Journal of Physical Oceanography, 2(4): 456–461.

HUPPERT H E, LINDEN P F, 1979. On heating a stable salinity gradient from below[J]. Journal of Fluid Mechanics, 95(3): 431–464.

KELLEY D E, 1990. Fluxes through diffusive staircases: A new formulation[J]. Journal of Geophysical Research: Oceans, 95(C3): 3365–3371.

KELLEY D E, FERNANDO H J S, GARGETT A E, et al, 2003. The diffusive regime of double-diffusive convection[J]. Progress in Oceanography, 56(3–4): 461–481.

KIMURA S, NICHOLLS K W, VENABLES E, 2015. Estimation of ice shelf melt rate in the presence of a thermohaline staircase[J]. Journal of Physical Oceanography, 45(1): 133–148.

LEWIS W T, INCROPERA F P, VISKANTA R, 1982. Interferometric study of stable salinity gradients heated from below or cooled from above[J]. Journal of Fluid Mechanics, 116: 411–430.

MARMORINO G O, CALDWELL D R, 1976. Heat and salt transport through a diffusive thermohaline interface[J]. Deep Sea Research and Oceanographic Abstracts, 23(1): 59–67.

NEAL V T, NESHYBA S, DENNER W, 1969. Thermal stratification in the Arctic Ocean[J]. Science, 166(3903): 373–374.

SCHMITT R W, 1994. Double diffusion in oceanography[J]. Annual Review of Fluid Mechanics, 26: 255–285.

SHIBLEY N C, TIMMERMANS M L, CARPENTER J R, et al, 2017. Spatial variability of the Arctic Ocean's double‐diffusive staircase[J]. Journal of Geophysical Research: Oceans, 122(2): 980–994.

TURNER J S, 1965. The coupled turbulent transports of salt and and heat across a sharp density interface[J]. International Journal of Heat and Mass Transfer, 8(5): 759–767.

TURNER J S, STOMMEL H, 1964. A new case of convection in the presence of combined vertical salinity and temperature gradients[J]. Proceedings of the National Academy of Sciences of the United States of America, 52(1): 49–53.

TURNER J S, 1968. The behaviour of a stable salinity gradient heated from below[J]. Journal of Fluid Mechanics, 33(1): 183–200.

XIE LINGLING, LI MINGMING, LI MIN, 2015. Observational analysis of the double-diffusive convection in the deep Canada Basin[J]. Acta Oceanologica Sinica, 34(11): 71–79.

ZHOU SHENGQI, LU YUANZHENG, SONG XUELONG, et al, 2016. New layer thickness parameterization of diffusive convection in the ocean[J]. Dynamics of Atmospheres and Oceans, 73: 87–97.

Numerical simulation of diffusive convection staircase

XU Yan1, LIU Xiangquan1, 3, 4, SONG Rengang1, CEN Xianrong2, GUO Shuangxi2, ZHOU Shengqi2

1. College of Electronic, Communication and Physics, Shandong University of Science & Technology, Qingdao 266590, China; 2. State Key Laboratory of Tropical Oceanography (South China Sea Institute of Oceanology, Chinese Academy of Sciences), Guangzhou 510301, China; 3. State Key Laboratory of Integrated Optoelectronics, Institute of Semiconductors, Chinese Academy of Sciences, Beijing 100083, China; 4. College of Materials Science and Optoelectronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China

Diffusive convection driven by temperature and salinity is a ubiquitous phenomenon in the mid and high latitudes of the ocean. Numerical simulation of diffusive convection can be used to study small-scale dynamic processes of the ocean. In this study, we analyze the formation mechanism of diffusive convection and establish a two-dimensional cavity model. The governing equation of the model is solved by the finite volume method, and the process of layering phenomenon is simulated numerically. The temporal evolution of temperature and salinity of the flow field is given, and the vortex structure of the velocity field is simulated. The processes of the formation and merging of diffusive convection staircases are analyzed, and a theoretical explanation for mass and energy transport processes is given. In addition, the diffusive convection under different heat flux density is discussed. It is found that the growth rate of diffusive convection staircases increases with the increase of heat flux density, and the cooling from the upper boundary plays an important role in accelerating its growth. The variation of heat flux density does not have a significant impact on the evolution of diffusive convection staircases.

diffusive convection; numerical simulation; salinity; temperature

2018-03-19;

2018-05-22. Editor: SUN Shujie

National Natural Science Foundation of China (11547037, 11604181, 41706029, 41776033); Natural Science Foundation of Guangdong Province (2016A030311042, 2016A030313155); Open Project Program of the State Key Laboratory of Tropical Oceanography (South China Sea Institute of Oceanology Chinese Academy of Sciences)(LTO1710)

P731.26

A

1009-5470(2019)01-0011-08

10.11978/2018028

2018-03-19；

2018-05-22。孫淑杰編輯

國家自然科學基金項目(11547037、11604181、41706029、41776033); 廣東省自然科學基金項目(2016A030311042、2016A030313155); 熱帶海洋環(huán)境國家重點實驗室(中國科學院南海海洋研究所)開放基金資助課題(LTO1710)

徐巖(1970—), 男, 山東省青島市人, 教授, 從事凝聚態(tài)物理中的非線性現(xiàn)象和統(tǒng)計物理研究。E-mail: x1y5@163.com

岑顯榮, 男, 博士, 主要從事海洋湍流混合的研究。E-mail: xrcen@scsio.ac.cn

CEN Xianrong. E-mail: xrcen@scsio.ac.cn