搖擺工況下窄矩形通道內(nèi)兩相沸騰摩擦壓降特性

陳沖，高璞珍，余志庭，陳先兵

（哈爾濱工程大學(xué)核安全與仿真技術(shù)國防重點學(xué)科實驗室，黑龍江 哈爾濱 150001）

搖擺工況下窄矩形通道內(nèi)兩相沸騰摩擦壓降特性

陳沖，高璞珍，余志庭，陳先兵

（哈爾濱工程大學(xué)核安全與仿真技術(shù)國防重點學(xué)科實驗室，黑龍江 哈爾濱 150001）

為了研究搖擺工況下窄矩形通道內(nèi)的兩相摩擦壓降特性，進行了一系列的熱工水力實驗和理論分析。結(jié)果表明，搖擺工況下流體會受到附加慣性力的作用且實驗回路的空間位置也會出現(xiàn)周期性的變化，兩相摩擦壓降梯度的波動振幅隨著搖擺角度和搖擺周期的增加而增加；隨著通道熱通量的增加或者系統(tǒng)壓強的減小，兩相摩擦壓降梯度的波動振幅和時均值逐漸增加。窄矩形通道內(nèi)的質(zhì)量流速隨著兩相摩擦壓降梯度的波動而波動，且具有相同的波動周期，由于流體加速和壓力傳播的速度不同，流量波動和摩擦壓降波動存在約1/4周期的相位差。

搖擺運動；窄矩形通道；相變；兩相流；傳熱

引 言

最近幾年，隨著船舶領(lǐng)域的快速發(fā)展，海洋條件對船舶動力裝置的影響越來越受到重視。船舶在海洋上航行時，海浪的作用將會使船舶發(fā)生搖擺、傾斜、俯仰和起伏運動，這些附加運動都會對船舶動力裝置中的流動傳熱特性產(chǎn)生明顯的影響。

過去幾年，很多學(xué)者從基礎(chǔ)實驗和理論分析方面研究了海洋工況對單相流體的影響特性。Murata等[1]、Ishida等[2]、Tan等[3-4]在搖擺條件下對單相自然循環(huán)的流動傳熱特性進行了研究，結(jié)果表明，由搖擺運動而產(chǎn)生的附加壓降將會導(dǎo)致自然循環(huán)回路的流量波動，并且流量波動的幅度隨著搖擺振幅和搖擺頻率的增加而增加，但是自然循環(huán)回路的平均流量隨著搖擺振幅和搖擺頻率的增加而減小。他們的結(jié)果還表明在搖擺工況下自然循環(huán)流量下降的主要原因是因為搖擺導(dǎo)致循環(huán)回路的阻力增加。Xing等[5-6]和Wang等[7-8]在強迫循環(huán)工況下研究了搖擺振幅、搖擺周期以及泵的驅(qū)動壓頭對流量波動振幅的影響，結(jié)果表明，在泵的驅(qū)動壓頭不變的工況下，流量的波動振幅隨著搖擺振幅和搖擺頻率的增加而增加，同時泵的驅(qū)動壓頭對流量以及摩擦阻力的波動有很大的影響。隨著循環(huán)回路驅(qū)動壓頭的增加，流量和摩擦壓降的波動振幅迅速減小，當(dāng)回路驅(qū)動壓頭大于10倍的附加壓降時，回路的流量將不會出現(xiàn)明顯的波動。高高高等[9]提出并建立了核動力裝置一回路冷卻劑受海洋條件影響的數(shù)學(xué)模型，此模型對流體在海洋條件下的受力進行了詳細的分析，并給出了附加壓降的積分關(guān)系式，為以后的海洋條件下理論模型的建立奠定了基礎(chǔ)。Yan等[10-13]數(shù)值模擬了搖擺工況下單相流體的流動傳熱特性，模擬結(jié)果表明，搖擺工況下流量波動主要受到附加慣性切向力的影響，而附加慣性向心力和科氏力的影響可以忽略不計；同時通道尺寸和Reynolds數(shù)對流體的流動傳熱特性有很大的影響，隨著通道尺寸的減小和Reynolds數(shù)的增加，搖擺對流體流動傳熱特性的影響效果逐漸減小。

Xing等[14]在搖擺工況下對空氣-水的流動特性進行了實驗研究，結(jié)果表明在層流和過渡區(qū)域瞬時摩擦壓降隨著搖擺運動而產(chǎn)生波動，同時隨著液相和氣相流速的增加摩擦壓降波動的振幅逐漸減小，在湍流區(qū)域瞬時摩擦壓降并沒有出現(xiàn)明顯的周期性波動。Jin等[15-17]研究了搖擺工況下窄矩形通道內(nèi)空氣-水的流動特性，結(jié)果表明隨著窄矩形通道窄縫高度的減小，附加壓降對兩相壓降的影響逐漸減小。通道的空泡份額隨著搖擺角度和搖擺頻率的增加而增加，同時搖擺工況下的時均空泡份額小于穩(wěn)態(tài)工況下的空泡份額。

綜上所述，搖擺運動對單相流體的流動特性、傳熱特性都有很大的影響，但是搖擺對兩相流體影響特性的相關(guān)研究比較少，大部分都集中在絕熱工況下空氣-水兩相混合物的實驗研究，而加熱工況下?lián)u擺對兩相沸騰摩擦壓降的影響特性還不是很清楚。緊湊式換熱器在船舶動力裝置中被廣泛使用，主要是因為緊湊式換熱器在單位體積內(nèi)具有較大的換熱面積和換熱效率。一般來說，船用緊湊式換熱器都是由一系列窄矩形通道組成，窄縫的高度為0.5～2 mm[18]。很多專家學(xué)者針對窄矩形通道和微通道進行了一系列的沸騰實驗研究[19-25]，但是這些實驗研究都是在穩(wěn)態(tài)工況下進行的，而海洋條件下窄矩形通道和微通道的熱工水力特性還沒有得到詳細的研究。為了更好地了解海洋條件下窄矩形通道內(nèi)的熱工水力特性，本文針對2 mm×40 mm的窄矩形通道進行了一系列的熱工水力實驗研究。

1 實驗裝置

機械搖擺熱工水力實驗裝置如圖1所示，主要由實驗回路、數(shù)據(jù)采集系統(tǒng)和機械搖擺臺組成，實驗回路的最大質(zhì)量流速為5000 kg·m?2·s?1，最大系統(tǒng)壓強為3 MPa，機械搖擺臺的最大載重為2 t。實驗回路主要由實驗主回路和冷卻回路組成，冷卻回路主要包括冷凝器、冷卻水泵、冷卻塔和水柜。實驗主回路主要包括預(yù)熱器、過濾器、窄矩形實驗段、冷凝器、主泵、穩(wěn)壓器、電磁流量計和高頻直流電源，高頻直流電源用來加熱窄矩形實驗段，其最大的加熱功率為100 kW，窄矩形通道的內(nèi)尺寸為2 mm×40 mm，長度為1100 mm。預(yù)熱器的主要作用是加熱冷卻水并維持實驗段的入口水溫，其最大的加熱功率為45 kW。整個實驗回路的系統(tǒng)壓強由穩(wěn)壓器通過調(diào)節(jié)氮氣的進氣量控制。實驗回路的流量是由電磁流量計進行測量，其測量誤差在±0.3%以內(nèi)；實驗段的出入口水溫由兩個N型熱電偶進行測量，測量誤差在±0.2℃以內(nèi)；實驗段的壓降用兩個不同的壓差傳感器進行測量，量程分別為50 kPa和60 kPa，測量誤差在±0.2%以內(nèi)。

機械搖擺臺由三相異步電機和減速箱提供動力并由曲柄搖桿機構(gòu)進行控制，搖擺角度隨時間的變化為

搖擺角速度和搖擺角加速度分別為

式中，θt、ωt、βt、θmax和T分別為瞬時搖擺角度、搖擺角速度、搖擺角加速度、最大搖擺角度和搖擺周期。

圖1 機械搖擺熱工水力實驗裝置Fig.1 Mechanical rolling thermal-hydraulic experimental facility

2 數(shù)據(jù)處理

在搖擺工況下實驗通道內(nèi)流體的受力分析如圖2所示，與穩(wěn)態(tài)工況相比搖擺工況下流體會額外受到附加慣性切向力、附加慣性法向力和科氏力的作用，這些附加慣性力會對流體產(chǎn)生相應(yīng)的加速度和附加壓降。

式中，ace為附加慣性法向加速度，m·s?2；ata為附加慣性切向加速度，m·s?2；aco為科氏加速度，m·s?2；u為流體的速度，m·s?2；r為搖擺軸到流體質(zhì)點的距離。

如圖2所示，實驗通道的附加壓降為

將式(4)～式(6)代入式(7)中得

由于整個測壓裝置都安裝在搖擺實驗臺上，而且連接窄矩形實驗段和壓差傳感器的引壓管中充滿了去離子水，所以在實驗臺搖擺過程中，引壓管中去離子水也會產(chǎn)生相應(yīng)的附加壓降，而且引壓管中的水溫和加熱實驗通道中的水溫存在差異，這就導(dǎo)致了實驗測量的壓差和實驗段中的實際壓差并不一致。在計算窄矩形通道的實際壓降時需要考慮附加壓降和密度差帶來的影響。

式中，Δpmea為壓差傳感器測量的壓差，Δptp,f為兩相摩擦壓降，Δptp,acc為兩相加速壓降，Δptp,g為兩相重位壓降，Δpadd,tot為整個測壓回路的附加壓降，ρc為引壓管中冷卻水的密度，Δpadd,21、Δpadd,42、Δpadd,34和Δpadd,13分別為實驗段21、42、34和13的附加壓降，如圖2所示。

圖2 搖擺工況下流體的受力分析Fig.2 Force analysis of fluid element under rolling motion

實驗通道內(nèi)的重位壓降和加速壓降由均相流模型得

式中，g為重力加速度，xe為平衡態(tài)含氣率，G為質(zhì)量流速，ρg和ρf分別為飽和汽、飽和水的密度。

3 實驗結(jié)果與討論

3.1 搖擺對兩相摩擦壓降梯度的影響

搖擺角度和周期對窄矩形實驗通道內(nèi)摩擦壓降梯度的影響如圖3和圖4所示，A15T10意味著搖擺角度為15°、搖擺周期為10 s。通道內(nèi)的摩擦壓降梯度隨著實驗臺的搖擺運動而產(chǎn)生明顯的波動，摩擦壓降梯度的波動周期與搖擺周期一致，同時摩擦壓降梯度的波動振幅隨著搖擺角度的增加而逐漸增加。這主要是因為隨著搖擺臺的運動，實驗主回路會出現(xiàn)明顯的空間變化，來自于加熱實驗通道的汽水混合物會在冷凝入口區(qū)域出現(xiàn)周期性的擁堵和疏通，如圖5所示。當(dāng)搖擺臺向負最大角度（?θmax）運動過程中，由于浮力的作用氣泡會聚集在冷凝器入口腔室的最高處，同時氣泡會聚集形成汽團，阻礙汽水混合物的流動，導(dǎo)致了實驗段出口至冷凝器入口的空間靜壓迅速升高，從而引起實驗段通道內(nèi)的飽和水溫升高，在相同的熱通量工況下，通道內(nèi)的沸騰程度會逐漸減小，質(zhì)量含氣率減小，從而導(dǎo)致摩擦壓降梯度減??；當(dāng)搖擺臺向正最大角度（+θmax）運動過程中，由于實驗回路的空間位置有利于擁堵在冷凝器入口的汽團迅速流入冷凝器冷卻，實驗段出口至冷凝器入口的空間靜壓減小，實驗通道內(nèi)的飽和水溫降低，沸騰程度增加，質(zhì)量含氣率增加，從而導(dǎo)致摩擦壓降梯度增加。隨著搖擺角度的增加，實驗回路空間位置的變化量增加，汽水混合物在冷凝器入口擁堵和疏通的程度增加，實驗通道內(nèi)質(zhì)量含氣率的波動振幅增加，導(dǎo)致了摩擦壓降梯度的波動振幅增加。

圖3 搖擺角度對兩相摩擦壓降梯度的影響Fig.3 Effect of rolling amplitude on two-phase frictional pressure gradient

圖4 搖擺周期對兩相摩擦壓降梯度的影響Fig.4 Effect of roll period on two-phase frictional pressure gradient

搖擺周期對兩相摩擦壓降梯度的影響如圖4所示，隨著搖擺周期的增加，兩相摩擦壓降梯度的波動振幅略微增加，這主要是因為隨著搖擺周期的增加，來自加熱通道的汽水混合物擁堵在冷凝器入口腔室的時間和程度都會增加，從而導(dǎo)致摩擦壓降梯度的波動振幅增加。隨著搖擺周期的減小，兩相摩擦壓降梯度的波動振幅不是很均勻，主要是因為隨著搖擺周期的減小，汽水混合物在冷凝器入口腔室的聚合和破裂比較劇烈，壓強的波動幅度會出現(xiàn)不均勻的現(xiàn)象。

圖5 實驗回路隨搖擺運動的軌跡Fig.5 Trajectory of experimental loop under rolling motion

3.2 熱通量對兩相摩擦壓降梯度的影響

搖擺工況下熱通量對兩相摩擦壓降梯度的影響如圖6所示，系統(tǒng)壓強為1.20 MPa，時均質(zhì)量流速為250 kg·m?2·s?1。兩相摩擦壓降梯度的時均值隨著熱通量的增加而增加，這主要是因為隨著熱通量的增加，氣泡的產(chǎn)生和脫離頻率增加，流體的湍流程度增加，沸騰強度增加，通道的質(zhì)量含氣率增加，兩相摩擦壓降梯度的時均值增加。兩相摩擦壓降的波動振幅隨著熱通量的增加而迅速增加，其主要原因是隨著熱通量的增加，窄矩形通道出口的質(zhì)量含氣率增加，汽水混合物擁堵在冷凝器入口腔室的程度增加，摩擦壓降梯度的波動振幅增加。

圖6 熱通量對兩相摩擦壓降梯度的影響Fig.6 Effect of heat flux on two-phase frictional pressure gradient

3.3 系統(tǒng)壓強對兩相摩擦壓降梯度的影響

搖擺工況下系統(tǒng)壓強對兩相摩擦壓降梯度的影響如圖7所示，通道的熱通量為341 kW·m?2，時均質(zhì)量流速為527 kg·m?2·s?1。兩相摩擦壓降梯度的時均值隨著系統(tǒng)壓強的增加而逐漸減小，其主要原因分為兩個方面：第一是隨著系統(tǒng)壓強的增加，通道內(nèi)工質(zhì)的飽和溫度升高，在相同的熱通量工況下，氣泡數(shù)量以及脫離頻率減小；第二是隨著系統(tǒng)壓強的增加，氣泡的平均直徑迅速減小，也就是說隨著系統(tǒng)壓強的增加，氣泡的數(shù)量和直徑都減小，通道內(nèi)的沸騰程度減小，質(zhì)量含氣率減小，兩相摩擦壓降梯度減小。隨著系統(tǒng)壓強的增加兩相摩擦壓降梯度波動的振幅減小，其主要原因是隨著系統(tǒng)壓強的增加，實驗通道出口質(zhì)量含氣率減小，汽水混合物在冷凝器入口腔室的擁堵程度減小。

圖7 系統(tǒng)壓降對兩相摩擦壓降梯度的影響Fig.7 Effect of system pressure on two-phase frictional pressure gradient

3.4 摩擦壓降梯度對質(zhì)量流速的影響

兩相摩擦壓降梯度對質(zhì)量流速的影響如圖8所示，質(zhì)量流速隨著兩相摩擦壓降梯度的波動而波動，它們具有相同的波動周期，而且其相位差約為1/4周期。主要是因為隨著摩擦壓降梯度的周期性變化，通道內(nèi)流體的加速度將會出現(xiàn)周期性的變化，同時流體的速度為流體加速度的積分，從而導(dǎo)致兩相摩擦壓降梯度和質(zhì)量流速之間的相位差。隨著摩擦壓降梯度的增加，通道內(nèi)的阻力增加，質(zhì)量流速將會減小，反之隨著摩擦壓降梯度的減小，通道內(nèi)的阻力減小，在主泵驅(qū)動壓頭的作用下，質(zhì)量流速將會增加。

圖8 質(zhì)量流速隨兩相摩擦壓降梯度的變化Fig.8 Mass flux changes with two-phase frictional pressure gradient

4 結(jié) 論

（1）窄矩形通道內(nèi)兩相摩擦壓降梯度的波動周期與搖擺周期一致，且波動振幅隨著搖擺角度和搖擺周期的增加而增加，影響兩相摩擦壓降梯度波動的主要原因是搖擺工況周期性地收變了實驗回路的空間位置，汽水混合物在冷凝器入口處出現(xiàn)周期性的擁堵和疏通，導(dǎo)致實驗通道內(nèi)的靜壓出現(xiàn)周期性的增加和減小，通道內(nèi)的質(zhì)量含氣率出現(xiàn)周期性的變化。

（2）兩相摩擦壓降梯度的時均值和波動振幅隨著熱通量的增加而增加，隨著系統(tǒng)壓強的增加而減小，其主要原因是窄矩形通道內(nèi)的沸騰程度和出口含氣率隨著熱通量的增加和系統(tǒng)壓強的減小而增加。

（3）窄矩形通道內(nèi)的質(zhì)量流速和兩相摩擦壓降具有相同的波動周期，且存在大約為1/4周期的相位差。

符 號 說 明

A——搖擺角度，（°）

F——附加慣性力，N

G——質(zhì)量流速，kg·m2·s?1

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

l——通道長度，m

p——系統(tǒng)壓強，MPa

Δptp——兩相壓降，kPa

q——熱通量，kW·m?2

T——搖擺周期，s

t——時間，s

xe——平衡態(tài)含氣率

β——角加速度，rad·s?2

θmax——最大搖擺角度，（°）

ρ——密度，kg·m?3

ω——角速度，rad·s?1

下角標(biāo)

acc——加速壓降

add——附加壓降

ave——時均值

c——冷卻水

f——摩擦壓降

g——重位壓降

max——最大值

mea——測量值

min——最小值

tp——兩相工況

[1] Murata H, Sawada K, Kobayashi M. Experimental investigation ofnatural convection in a core of a marine reactor in rolling motion [J].Journal of Nuclear Science and Technology, 2000, 37(6): 209-517.

[2] Ishida T, Yoritsune T. Effects of ship motions on natural circulation of deep sea research reactor DRX [J].Nuclear Engineering and Design, 2002, 215(1/2): 51-67.

[3] Tan S, Su G, Gao P. Experimental and theoretical study on single-phase natural circulation flow and heat transfer under rolling motion condition [J].Applied Thermal Engineering, 2009, 29(14/15): 3160-3168.

[4] Tan S, Su G, Gao P. Heat transfer model of single-phase natural circulation flow under a rolling motion condition [J].Nuclear Engineering and Design, 2009, 239(10): 2212-2216.

[5] Xing D, Yan C, Sun L, Wang C. Effect of rolling motion on single-phase laminar flow resistance of forced circulation with different pump head [J].Annals of Nuclear Energy, 2013, 54: 141-148.

[6] Xing D, Yan C, Sun L, Xu C. Effects of rolling on characteristics of single-phase water flow in narrow rectangular ducts [J].Nuclear Engineering and Design, 2012, 247: 221-229.

[7] Wang C, Li X, Wang H, Gao P. Experimental study on friction and heat transfer characteristics of pulsating flow in rectangular channel under rolling motion [J].Progress in Nuclear Energy, 2014, 71: 73-81.

[8] Wang C, Wang S, Wang H, Gao P. Investigation of flow pulsation characteristic in single-phase forced circulation under rolling motion [J].Annals of Nuclear Energy, 2014, 64: 50-56.

[9] Gao Puzhen(高高高), Pang Fengge(龐鳳閣), Wang Zhaoxiang(王兆祥). Mathematical model of primary coolant in nuclear power plant influenced by ocean conditions [J].Journal of Harbin Engineering University(哈爾濱工程大學(xué)學(xué)報), 1997, 18(1): 26-29.

[10] Yan B H, Gu H Y. CFD analysis of flow and heat transfer of turbulent pulsating flow in a tube in rolling motion [J].Annals of Nuclear Energy, 2011, 38(9): 1833-1841.

[11] Yan B H, Gu H Y, Yu L. Effects of rolling motion on the flow and heat transfer of turbulent pulsating flow in channels [J].Progress in Nuclear Energy, 2012, 56: 24-36.

[12] Yan B H, Yu L, Yang Y H. Effects of ship motions on laminar flow in tubes [J].Annals of Nuclear Energy, 2010, 37(1): 52-57.

[13] Yan B H, Yu L, Yang Y H. Heat transfer with laminar pulsating flow in a channel or tube in rolling motion [J].International Journal of Thermal Sciences, 2010, 49(6): 1003-1009.

[14] Xing D, Yan C, Sun L, Jin G, Tan S. Frictional resistance of adiabatic two-phase flow in narrow rectangular duct under rolling conditions [J].Annals of Nuclear Energy, 2013, 53: 109-119.

[15] Jin G, Yan C, Sun L, Xing D. Effect of rolling motion on transient flow resistance of two-phase flow in a narrow rectangular duct [J].Annals of Nuclear Energy, 2014, 64: 135-143.

[16] Jin G, Yan C, Sun L, Wang Y, Yan C. Research on frictional resistance of bubbly flow in rolling rectangular ducts [J].Nuclear Engineering and Design, 2014, 278: 108-116.

[17] Jin G, Yan C, Sun L, Xing D, Zhou B. Void fraction of dispersed bubbly flow in a narrow rectangular channel under rolling conditions [J].Progress in Nuclear Energy, 2014, 70: 256-265.

[18] Lee H J, Lee S Y. Heat transfer correlation for boiling flows in small rectangular horizontal channels with low aspect ratios [J].International Journal of Multiphase Flow, 2001, 27: 2043-2062.

[19] Kim S, Mudawar I. Flow condensation in parallel micro-channels(Ⅱ): Heat transfer results and correlation technique [J].International Journal of Heat and Mass Transfer, 2012, 55: 984-994.

[20] Sun L, Mishima K. An evaluation of prediction methods for saturated flow boiling heat transfer in mini-channels [J].International Journal of Heat and Mass Transfer, 2009, 52: 5323-5329.

[21] Jaeseon L, Issam M. Two-phase flow in high-heat-flux micro-channel heat sink for refrigeration cooling applications(Ⅱ): Heat transfer characteristics [J].International Journal of Heat and Mass Transfer, 2005, 48: 941-955.

[22] Wei Jinghua(魏敬華), Pan Liangming(潘良明), Xu Jianjun(徐建軍), Huang Yanping(黃彥平). Numerical simulation of subcooled flow boiling in vertical rectangular channel under additional inertial forces [J].CIESC Journal(化工學(xué)報), 2011, 62(5): 1239-1245.

[23] Yuan Dewen(袁德文), Pan Liangming(潘良明), Chen Deqi(陳德奇). Model for single bubble growth of subcooled flow boiling in vertical narrow rectangular channel [J].CIESC Journal(化工學(xué)報), 2009, 60(11): 2723-2728.

[24] Pan Liangming(潘良明), He Chuan(何川), Xin Daoming(辛明道), Wu Xiaohang(吳小航). Comparison of bubble behavior for double-side and single-side heating of subcooled flow boiling in narrow channels [J].Journal of Chemical Industry and Engineering(China)(化工學(xué)報), 2004, 55(9): 1519-1522.

[25] Chen C, Gao P, Tan S, Huang D. Effects of rolling motion on thermal-hydraulic characteristics of boiling flow in rectangular narrow channel [J].Annals of Nuclear Energy, 2015, 76: 504-513.

Two-phase frictional pressure drop characteristics of boiling flow in rectangular narrow channel under rolling motion

CHEN Chong, GAO Puzhen, YU Zhiting, CHEN Xianbing
(Fundamental Science on Nuclear Safety and Simulation Technology Laboratory,Harbin Engineering University,Harbin150001,Heilongjiang,China)

In order to investigate the two-phase frictional pressure drop characteristics of boiling flows in a rectangular narrow channel under rolling motion, a series of thermal hydraulic experiments and theoretical analysis are performed. The results demonstrate that the additional inertial force is imposed on the fluid and the space of experimental loop will vary periodically under rolling motion. The fluctuation amplitude of the two-phase frictional gradient increases with increasing rolling angle and rolling period. The fluctuation amplitude and time average value of the two-phase frictional pressure gradient increase with increasing heat flux, while it decreases with the increase of system pressure. The mass flux varies with the fluctuation of frictional pressure gradient at the same period. The phase change between the fluctuation of mass flux and frictional pressure gradient is approximately equal to 1/4 rolling period due to the velocity difference of the pressure propagation and mass flux increases.

rolling motion; rectangular narrow channel; phase change; two-phase flow; heat transfer

GAO Puzhen, gaopuzhen@sina.com

10.11949/j.issn.0438-1157.20150126

TL 334

：A

：0438—1157（2015）10—3874—07

2015-01-26收到收收，2015-05-21收到收收收。

聯(lián)系人：高高高。

：陳陳(1989—)，男，博博研究生。

國家自然科學(xué)基金項目(50806014)。

Received date: 2015-01-26.

Foundation item: supported by the National Natural Science Foundation of China(50806014).