莊大偉，楊藝菲，胡海濤，丁國良(上海交通大學(xué)制冷與低溫工程研究所，上海 200240)

?

豎直平板間液橋形狀的觀測(cè)與預(yù)測(cè)模型開發(fā)

莊大偉，楊藝菲，胡海濤，丁國良
(上海交通大學(xué)制冷與低溫工程研究所，上海 200240)

摘要：了解析濕工況下液橋?qū)艹崾綋Q熱器性能的影響，需要對(duì)液橋形狀進(jìn)行實(shí)驗(yàn)研究并建立描述液橋形狀的方法。通過搭建可視化實(shí)驗(yàn)臺(tái)觀察了液橋在豎直平板間的形狀，并測(cè)量了液橋的接觸線和接觸角。實(shí)驗(yàn)研究了不同體積、不同平板間距、不同平板材料間形成的液橋，并通過實(shí)驗(yàn)數(shù)據(jù)開發(fā)了能夠描述液橋形狀的關(guān)聯(lián)式，包括接觸線伸長比和接觸角分別與Bond數(shù)的關(guān)聯(lián)式。還根據(jù)關(guān)聯(lián)式建立了描述液橋形狀的方法并將結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，對(duì)比結(jié)果顯示接觸線的描述方法平均誤差為3.4%，接觸角的描述方法平均誤差為7.9%。

關(guān)鍵詞：液橋形狀；接觸角；預(yù)測(cè)；模型；實(shí)驗(yàn)驗(yàn)證

2015-11-11收到初稿，2016-02-16收到修改稿。

聯(lián)系人：丁國良。第一作者：莊大偉（1987—），男，博士后。

Received date: 2015-11-11.

Foundation item: supported by the National Natural Science Foundation of China (51576122) and the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (51521004).

引　言

空調(diào)蒸發(fā)器的空氣側(cè)運(yùn)行在濕工況下[1-3]，使得濕空氣中的水蒸氣會(huì)在翅片表面冷凝形成液滴[4-6]；液滴會(huì)隨著冷凝過程不斷長大[7-9]，并不斷與相鄰翅片上的液滴發(fā)生接觸并連接形成大量的液橋[10]。這些液橋若不能及時(shí)排出翅片，就會(huì)堵塞換熱器空氣側(cè)流道、占用翅片換熱面積，從而嚴(yán)重影響換熱器的性能[8, 11-13]。要提高濕工況下?lián)Q熱器性能，就需要使形成的液橋能夠快速排出換熱器。為了能夠使液橋快速排出，就需要對(duì)液橋的形狀進(jìn)行了解，并建立描述液橋形狀的方法。

液橋形狀包括液橋-翅片界面形狀和液橋-濕空氣界面形狀，可通過液橋三相接觸線和接觸角來描述[14]：① 液橋-翅片界面是由三相接觸線圍成的平面；② 液橋-濕空氣界面是以接觸角為斜率的曲線繞接觸線一周圍成的曲面。因此要實(shí)現(xiàn)液橋形狀的描述，就必須開發(fā)液橋三相接觸線和接觸角的定量預(yù)測(cè)模型。

液橋三相接觸線和接觸角模型的開發(fā)需要能夠反映液橋所受表面張力、重力這兩個(gè)因素的耦合影響。針對(duì)接觸線和接觸角的現(xiàn)有研究主要針對(duì)受到表面張力和軸向重力作用下的液橋，研究包括實(shí)驗(yàn)[14-17]和數(shù)值模型[18-21]兩類。

已有實(shí)驗(yàn)研究包括對(duì)液橋輪廓的研究[14-16]和對(duì)接觸線和接觸角的測(cè)量[17]。Verges等[14]通過實(shí)驗(yàn)研究了液橋接觸線的平均曲率，結(jié)果表明受到表面張力和軸向重力的液橋與平板間的接觸線可以由圓形來表示。Ferrera等[15]采用TIFA方法對(duì)表面張力進(jìn)行了計(jì)算并通過實(shí)驗(yàn)分析了液橋形狀與表面張力之間的關(guān)系。Chen等[17]通過實(shí)驗(yàn)研究了接觸角滯后性對(duì)液橋輪廓的影響。

已有數(shù)值模型的研究包括對(duì)液橋接觸線和接觸角的描述方法[18-21]。液橋在給定接觸角以及忽略重力的情況下可以簡(jiǎn)化為旋轉(zhuǎn)軸對(duì)稱模型[18]。Vogel[20]分析了軸向重力對(duì)于液橋形狀的影響，結(jié)果表明軸向重力會(huì)導(dǎo)致液橋形狀的變形。Chen等[21]建立了能夠反映接觸角滯后性對(duì)液橋形狀影響的模型，并發(fā)現(xiàn)在給定的平板間距下液橋在平衡狀態(tài)可能會(huì)出現(xiàn)兩種不同的形狀。

上述關(guān)于液橋的研究是針對(duì)受表面張力和軸向重力影響的情況。實(shí)際換熱器豎直翅片間附著的液橋所受的重力方向?yàn)檠匾簶驈较蚍较颉较蛑亓ψ饔貌坏珪?huì)使得液橋接觸線沿豎直方向發(fā)生伸長現(xiàn)象，還會(huì)導(dǎo)致接觸線上各點(diǎn)的接觸角隨著方位角不斷變化。已有研究中圓形接觸線和等接觸角模型難以描述實(shí)際換熱器中受表面張力和徑向重力共同作用的液橋的形狀。

本研究目的是基于實(shí)驗(yàn)，提出描述徑向重力-表面張力耦合作用下豎直翅片間液橋的形狀的方程，包括液橋三相接觸線方程、接觸角方程。

1　實(shí)驗(yàn)設(shè)計(jì)

1.1實(shí)驗(yàn)?zāi)康?/p>

換熱器豎直翅片間附著的液橋的形狀如圖1所示。由于受到徑向重力作用，液橋的形狀沿z軸方向非對(duì)稱，從而導(dǎo)致接觸線沿豎直方向伸長且接觸角不斷變化。本文的目的為通過實(shí)驗(yàn)觀測(cè)豎直翅片間的液橋形狀，并測(cè)量液橋接觸線和接觸角；在此基礎(chǔ)上開發(fā)接觸線和接觸角的定量預(yù)測(cè)模型，從而實(shí)現(xiàn)對(duì)液橋形狀的描述。

圖1　翅片間液橋示意圖Fig. 1　Schematic diagram of water bridge between fins

1.2實(shí)驗(yàn)裝置

為了研究豎直平板間液橋的接觸線和接觸角，實(shí)驗(yàn)裝置需要滿足以下要求：① 液橋可以在豎直平行平板間形成；② 平板間距可以調(diào)整；③ 可以對(duì)液橋形狀進(jìn)行觀察和記錄；④ 可以從不同角度對(duì)液橋的接觸角進(jìn)行觀察和測(cè)量。

根據(jù)上述要求，實(shí)驗(yàn)裝置包括液橋形成系統(tǒng)、平板固定系統(tǒng)、拍攝系統(tǒng)、旋轉(zhuǎn)系統(tǒng)，完整實(shí)驗(yàn)裝置如圖2所示。

圖2　實(shí)驗(yàn)裝置Fig. 2　Schematic diagram of experimental apparatus1—syringe；2—motor；3—platform；4—bracket；5—retaining nuts；6—shims；7—x-y-z movable platform；8—CCD camera；9—light source；10—microscopic lens；11—computer；12—rotating arm；13—jack；14—protractor

液體注射系統(tǒng)包括自動(dòng)注射器和步進(jìn)電機(jī)。由于人工注射形成的液橋可以通過人為控制使形成的液橋趨于對(duì)稱，比冷凝形成的液橋更有利于實(shí)驗(yàn)觀測(cè)，因此本文采用人工注射的方式在平板間形成液橋。其中注射器的精度為±0.02 μl，注射器針頭直徑為0.3 mm。步進(jìn)電機(jī)可以通過計(jì)算機(jī)控制開停，從而實(shí)現(xiàn)對(duì)自動(dòng)注射器的上升和下降的調(diào)節(jié)。

平板固定系統(tǒng)包括旋轉(zhuǎn)平臺(tái)、支架、螺絲、墊片和移動(dòng)平臺(tái)。其中旋轉(zhuǎn)平臺(tái)旋轉(zhuǎn)精度為±10′；移動(dòng)平臺(tái)的精度為±0.01 mm，通過移動(dòng)平臺(tái)可以實(shí)現(xiàn)對(duì)平板位置的調(diào)節(jié)；該系統(tǒng)中的支架可以實(shí)現(xiàn)旋轉(zhuǎn)系統(tǒng)和固定系統(tǒng)的連接；平板可以通過螺絲固定在旋轉(zhuǎn)平臺(tái)上；4種不同的墊片用來調(diào)節(jié)兩個(gè)平板之間的間距。平板固定系統(tǒng)各部件結(jié)構(gòu)與安裝位置見圖3。

圖3　平板固定系統(tǒng)照片F(xiàn)ig. 3　Photos of solid sample fixed system

拍攝系統(tǒng)包括CCD相機(jī)、LED冷光源、透鏡和電腦。其中LED冷光源可發(fā)出亮度可調(diào)節(jié)的單色光，用來觀察液橋輪廓；透鏡的放大倍數(shù)為0.7～4.5，光學(xué)變形率為0.01%；CCD相機(jī)分辨率為752×480，可以拍攝液橋輪廓并儲(chǔ)存在計(jì)算機(jī)中。

旋轉(zhuǎn)系統(tǒng)包括旋轉(zhuǎn)架、千斤頂和量角器。旋轉(zhuǎn)架為10 mm厚的不銹鋼板；通過千斤頂可以使整個(gè)旋轉(zhuǎn)臺(tái)及旋轉(zhuǎn)臺(tái)上的系統(tǒng)進(jìn)行旋轉(zhuǎn)運(yùn)動(dòng)；通過量角器控制旋轉(zhuǎn)角度。

1.3實(shí)驗(yàn)樣件

考慮到換熱器翅片上冷凝所形成的液態(tài)水可認(rèn)為是蒸餾水，因此實(shí)驗(yàn)液體試樣選擇為蒸餾水。

本文實(shí)驗(yàn)的固體樣件采用了5種潤濕性不同的材料制成的固體平板，包括鋁、銅、不銹鋼、光學(xué)玻璃、有機(jī)玻璃。其中鋁、銅和不銹鋼是換熱器中的常見材料，光學(xué)玻璃和有機(jī)玻璃是透明的，所以能夠方便地觀察到液橋的接觸線形狀。5種材料的潤濕性見表1。

表1　5種固體材料表面的潤濕性Table 1　Wettability of five types of solid plane surfaces

1.4實(shí)驗(yàn)步驟

詳細(xì)的實(shí)驗(yàn)操作步驟如下。

（1）將平板與旋轉(zhuǎn)平臺(tái)通過4個(gè)螺絲進(jìn)行連接，連接過程可采用不同厚度的墊片實(shí)現(xiàn)豎直平行平板間距的調(diào)節(jié)。

（2）通過電機(jī)以及移動(dòng)平臺(tái)使注射器處于合適的位置，注射蒸餾水在豎直平行平板間形成液橋，并通過控制注射器來控制形成液橋的體積。調(diào)節(jié)透鏡和光源使CCD相機(jī)可以獲得清晰的液橋輪廓圖像。

（3）通過調(diào)整千斤頂使得旋轉(zhuǎn)架與其他機(jī)構(gòu)一起旋轉(zhuǎn)，同時(shí)將旋轉(zhuǎn)平臺(tái)進(jìn)行旋轉(zhuǎn)，保證液橋處于豎直向下方向。當(dāng)量角器達(dá)到所需角度后停止旋轉(zhuǎn)并拍攝該角度下的液橋輪廓。

（4）測(cè)量不同觀察角度下的液橋最大接觸角和最小接觸角，以及分別測(cè)量液橋接觸線長軸和短軸的長度。

1.5實(shí)驗(yàn)工況與誤差分析

為了通過實(shí)驗(yàn)對(duì)液橋接觸線和接觸角進(jìn)行研究，實(shí)驗(yàn)工況包括不同液橋體積（3～21 μl）、不同平板間距（0.5～2.0 mm）和不同觀察角度（0°～90°）。

接觸線長軸和短軸的測(cè)量誤差可以由像素分析法得出為0.005 mm。接觸角的測(cè)量誤差由接觸角測(cè)量?jī)x廠商提供為0.1°。

2　液橋接觸線的實(shí)驗(yàn)結(jié)果及預(yù)測(cè)模型

2.1接觸線的實(shí)驗(yàn)觀察結(jié)果

圖4是通過CCD相機(jī)拍攝的豎直平板間液橋的輪廓。圖中液橋體積為6 μl，平板間距為1 mm，平板材料為有機(jī)玻璃。

圖4　豎直平板間液橋輪廓Fig. 4　Contour of water bridge between vertical plates

從圖4(a) 中可得，液橋與豎直平板間的接觸線近似為橢圓形。圖4(b) 中頂部彎液面與底部彎液面形狀不同，此時(shí)液橋長度為2L，即為接觸線橢圓的長軸長。圖4(c) 為液橋的俯視圖，即方位角為90°的情況下的液橋輪廓照片，在該視圖下液橋的兩個(gè)彎液面近似相同，液橋此時(shí)的長度為2W，即為接觸線橢圓的短軸長。

2.2接觸線預(yù)測(cè)模型

根據(jù)實(shí)驗(yàn)觀測(cè)結(jié)果，接觸線可以表達(dá)為橢圓方程，該方程在圓柱坐標(biāo)系下可以表示為

其中，r是橢圓弦長; φ是方位角，和觀測(cè)角相等; L是長軸長度; W是短軸長度; β是伸長比, 表示液橋在重力方向上的伸長程度。

由式（1）和式（2）可知，影響接觸線形狀的關(guān)鍵因素是伸長比β。圖5顯示了不同Bond數(shù)Bo下的伸長比，Bo定義如下

其中，ρ是蒸餾水的密度; g是重力加速度; d是平板間距; σ是氣液交界面的表面張力；L和W分別是橢圓接觸面的半長軸和半短軸。根據(jù)實(shí)驗(yàn)數(shù)據(jù)擬合得到伸長比β關(guān)于Bo的擬合方程

2.3接觸線預(yù)測(cè)模型的精度驗(yàn)證

如圖5所示，伸長比β隨著Bo增加而增加，且近似為線性關(guān)系。圖中共有21個(gè)數(shù)據(jù)點(diǎn)，所有數(shù)據(jù)點(diǎn)與接觸線的預(yù)測(cè)模型的誤差均在±5%的范圍內(nèi)，平均誤差為3.4%。

圖5　接觸線伸長比預(yù)測(cè)值與實(shí)驗(yàn)值的對(duì)比Fig. 5　Comparison between predicted results and experimental data of β

3　液橋接觸角的實(shí)驗(yàn)結(jié)果及預(yù)測(cè)模型

3.1接觸角的測(cè)量結(jié)果及預(yù)測(cè)模型

根據(jù)接觸角的測(cè)量結(jié)果可知，對(duì)于不同的平板材料上形成的液橋的接觸角與方位角的關(guān)系均呈三次方的函數(shù)關(guān)系

其中，θ(φ)是不同方位角的接觸角; φ是方位角; A、B、C、D是待定系數(shù)，可根據(jù)如下邊界條件得到

其中，θmax是在φ=0時(shí)的最大接觸角; θmin是在φ=π時(shí)的最小接觸角。

由式（5）和式（6）可以計(jì)算出待定系數(shù)A、B、C和D，如式（7）所示

由式（7）可得，最大接觸角和最小接觸角是計(jì)算不同方位角下接觸角的關(guān)鍵參數(shù)。

根據(jù)實(shí)驗(yàn)測(cè)試結(jié)果發(fā)現(xiàn)，液橋的最大接觸角近似等于水在平板材料上的前進(jìn)接觸角；最小接觸角的表達(dá)式如式(8)所示，該公式預(yù)測(cè)值與實(shí)驗(yàn)值平均偏差為5.9%。

將式(8)代入式(7)，可以得到預(yù)測(cè)液橋接觸角的模型，如式(9)所示

3.2液橋接觸角預(yù)測(cè)模型的精度驗(yàn)證

圖6顯示了根據(jù)式（9）對(duì)液橋接觸角的預(yù)測(cè)值與實(shí)驗(yàn)值的對(duì)比。圖中共包含78個(gè)數(shù)據(jù)點(diǎn)，其中有95%的點(diǎn)處于誤差±20%的范圍內(nèi)，平均誤差為7.9%。從圖中可以看出，模型在接觸角范圍為20°～90°的范圍內(nèi)可以較好地預(yù)測(cè)液橋的接觸角。

圖6　接觸角預(yù)測(cè)結(jié)果與實(shí)驗(yàn)值的對(duì)比Fig. 6　Comparison between predicted results and experimental data of θ(φ)

4　結(jié)　論

（1）通過可視化實(shí)驗(yàn)對(duì)描述豎直翅片間液橋形狀的關(guān)鍵參數(shù)包括接觸線和接觸角進(jìn)行觀測(cè)發(fā)現(xiàn)，液橋接觸線形狀近似為橢圓形，而接觸角與方位角呈三次多項(xiàng)式關(guān)系。

（2）開發(fā)了基于橢圓方程的液橋三相接觸線預(yù)測(cè)模型，模型對(duì)于橢圓伸長比的預(yù)測(cè)值與實(shí)驗(yàn)值的平均誤差為3.4%。

（3）開發(fā)了基于Bo和翅片表面潤濕性的液橋接觸角預(yù)測(cè)模型，模型的預(yù)測(cè)值與實(shí)驗(yàn)值的平均誤差為7.9%。

符號(hào)說明

Bo——Bond數(shù)

d——平板間距，mm

g——重力加速度，9.81 m·s?2

L——長軸長度，mm

p——邊界點(diǎn)距離，pixel

r——橢圓半徑，mm

V——體積，μl

W——短軸長度，mm

β——伸長比

θ——接觸角，(°)

ρ——密度，kg·m?3

σ——表面張力，N·m?1

φ——方位角，(°)

下角標(biāo)

A——前進(jìn)

max——最大值

min——最小值

R——后退

References

[1]WANG C C, LEE W S, SHEU W J, et al. Investigation of the airside performance of the slit fin-and-tube heat exchangers [J]. International Journal of Refrigeration, 1999, 22(8): 595-603.

[2]劉金平, 戴紹碧, 葉立平. 析濕工況下波紋翅片管換熱器空氣側(cè)換熱與壓降特性 [J]. 化工學(xué)報(bào), 2009, 60(10): 2438-2445. LIU J P, DAI S B, YE L P. Air-side heat transfer and friction characteristics of wavy fin-and-tube heat exchangers under dehumidifying conditions [J]. CIESC Journal, 2009, 60(10): 2438-2445.

[3]COMINI G, NONINO C, SAVINO S. Modeling of conjugate conduction and heat and mass convection in tube-fin exchangers [J]. International Journal of Numerical Methods for Heat and Fluid Flow, 2008, 18(7/8): 954-968.

[4]楊藝菲, 莊大偉, 胡海濤, 等. 濕工況下平翅片平面凝水形成及運(yùn)動(dòng)過程的數(shù)值模擬與實(shí)驗(yàn)驗(yàn)證 [J]. 化工學(xué)報(bào), 2014, 65(S2): 140-147. YANG Y F, ZHUANG D W, HU H T, et al. Numerical simulation and experimental validation of water condensing and moving on plain-fin surface under dehumidifying conditions [J]. CIESC Journal, 2014, 65(S2): 140-147.

[5]WANG C C，LIN Y T, LEE C J. An airside correlation for plain fin-and-tube heat exchangers in wet conditions [J]. International Journal of Heat and Mass Transfer, 2000, 43(10): 1869-1872.

[6]李?yuàn)? 李慧珍, 康海軍, 等. 析濕正弦波紋翅片管束的熱質(zhì)交換特性研究 [J]. 西安交通大學(xué)學(xué)報(bào), 1996, 30(8): 57-63. LI B, LI H Z, KANG H J, et al. Experimental study on heat/mass transfer and pressure drop for dehumidifying sinusoid fin-and-tube heat exchanger [J]. Journal of Xi’an Jiaotong University, 1996, 30(8):57-63.

[7]YUN R, KIM Y B, KIM Y C. Air side heat transfer characteristics of plate finned tube heat exchangers with slit fin configuration under wet conditions [J]. Applied Thermal Engineering, 2009, 29(14/15): 3014-3020.

[8]MA X K, DING G L, ZHANG Y M, et al. Airside heat transfer and friction characteristics for enhanced fin-and-tube heat exchanger with hydrophilic coating under wet conditions [J]. International Journal of Refrigeration, 2007, 30(7): 1153-1167.

[9]郭軼波, 杜志敏, 陳萍, 等. 入口相對(duì)濕度對(duì)開縫翅片管換熱器空氣側(cè)換熱與壓降特性的影響 [J]. 化工學(xué)報(bào), 2010, 61(S2): 49-53. GUO Y B, DU Z M, CHEN P, et al. Effect of inlet air relative humidity on air side characteristics of slit fin-tube heat exchangers [J]. CIESC Journal, 2010, 61(S2): 49-53.

[10]ZHUANG D W, HU H T, DING G L, et al. Numerical model for liquid droplet motion on vertical plain-fin surface [J]. HVAC and R Research, 2014, 20(3): 332-343.

[11]馬小魁, 丁國良, 張圓明. 析濕工況下帶親水層開縫翅片管換熱器空氣側(cè)傳熱傳質(zhì)特性 [J]. 化工學(xué)報(bào), 2007, 58(8): 1911-1916. MA X K, DING G L, ZHANG Y M. Airside heat and mass transfer characteristics of split fin-and-tube heat exchangers under dehumidifying conditions [J]. Journal of Chemical Industry and Engineering(China), 2007, 58(8): 1911-1916.

[12]MIN J C, WEBB R L. Condensate carryover phenomena in dehumidifying finned-tube heat exchangers [J]. Experimental Thermal and Fluid Science, 2000, 22(3/4): 175-182.

[13]WANG C C, HSIEH Y C, LIN Y T. Performance of plate finned tube heat exchangers under dehumidifying conditions [J]. Journal of Heat Transfer, 1997; 119(1):109-17.

[14]VERGES M A, LARSON M C, BACOU R. Force and shapes of liquid bridges between circular pads [J]. Experimental Mechanics, 2001, 41(4): 351-357.

[15]FERRERA C, CABEZAS M G, MONTANERO J M. An experimental analysis of the linear vibration of axisymmetric liquid bridges [J]. Physics of Fluids, 2006, 18(8): 082105.

[16]ZHANG F Y, YANG X G, WANG C Y. Liquid water removal from a polymer electrolyte fuel cell [J]. Journal of the Electrochemical Society, 2006, 153(2): A225-A232.

[17]CHEN H, TANG T, AMIRFAZLI A. Effect of contact angle hysteresis on breakage of a liquid bridge [J]. European Physical Journal: Special Topics, 2015, 224(2): 277-288.

[18]VOGEL T I. Stability of a liquid drop trapped between two parallel planes(Ⅱ): General contact angles [J]. SIAM Journal on Applied Mathematics, 1989, 49(4): 1009-1028.

[19]ACERO F J, FERRERA C, CABEZAS M G, et al. Liquid bridge equilibrium contours between non-circular supports [J]. Microgravity Science and Technology, 2005, 17(2): 18-30.

[20]VOGEL T I. Stability of a liquid drop trapped between two parallel planes [J]. SIAM Journal on Applied Mathematics, 1987, 47(3): 516-525.

[21]CHEN H, TANG T, AMIRFAZLI A. Modeling liquid bridge between surfaces with contact angle hysteresis [J]. Langmuir, 2013, 29(10): 3310-3319.

Visualization and prediction model on shape of liquid bridge

ZHUANG Dawei, YANG Yifei, HU Haitao, DING Guoliang
(Institute of Refrigeration & Cryogenics Engineering, Shanghai Jiao Tong University, Shanghai 200240, China)

Abstract:Water bridge can affect the performance of fin-and-tube heat exchanger working under dehumidifying. In order to know the influence of water bridge on heat exchanger, the shape of water bridge needs to be studied. In this research, a visual experiment was used to develop the prediction model of water bridge shape. Based on the visual experiment, the contact line and the contact angles of water bridge were measured, covering various water volumes, solid plate separations and five different solid materials. The prediction model of contact line via correlated aspect ratio was developed, and the prediction model of contact angles via the correlated maximum and minimum contact angles was proposed. The mean deviation of aspect ratio of contact line was 3.4% and that of contact angles was 7.9%.

Key words:water bridge shape; contact angle; prediction; model; experimental validation

中圖分類號(hào)：TK 124

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：0438—1157（2016）06—2224—06

DOI：10.11949/j.issn.0438-1157.20151707

基金項(xiàng)目：國家自然科學(xué)基金項(xiàng)目（51576122）；國家自然科學(xué)基金創(chuàng)新研究群體科學(xué)基金項(xiàng)目（51521004）。

Corresponding author:Prof. DING Guoliang, glding@sjtu.edu.cn