徐志明，楊蘇武，朱新龍，張一龍，劉坐東

?

三角翼渦流發(fā)生器納米氧化鎂顆粒污垢特性

徐志明1，楊蘇武1，朱新龍1，張一龍2，劉坐東2

（1東北電力大學(xué)能源與動(dòng)力工程學(xué)院，吉林 吉林132012；2華北電力大學(xué)，北京 102206）

為了探究三角翼渦流發(fā)生器的納米氧化鎂顆粒污垢特性，選用粒徑為50 nm的氧化鎂顆粒配制的膠體溶液為研究對(duì)象，研究了在不同水浴溫度、顆粒濃度、流速、三角翼間距以及不同布置等工況下三角翼渦流發(fā)生器的污垢特性。結(jié)果表明：三角翼渦流發(fā)生器具有抑垢特性，水浴溫度、濃度、流速對(duì)其抑垢能力及結(jié)垢速率均有影響，水浴溫度升高、顆粒濃度降低以及工質(zhì)流速升高都會(huì)導(dǎo)致三角翼的抑垢能力增強(qiáng)。三角翼不同列間距的抑垢能力隨著布置方式的改變而改變，采用均勻分布布置形式時(shí)40 mm列間距的抑垢能力最強(qiáng)，而采用同列數(shù)布置時(shí)80 mm列間距的抑垢能力最強(qiáng)。

三角翼；氧化鎂；結(jié)垢；污垢特性；間距；膠體；納米粒子

引 言

顆粒污垢是換熱設(shè)備中污垢的一種，它是指懸浮在流體中的固體顆粒在換熱面上的積聚[1]。鎂垢則是其中比較常見的一種。許多學(xué)者對(duì)于顆粒污垢的特性進(jìn)行了研究。顧業(yè)梅等[2]對(duì)弧線管以及光管中顆粒污垢誘導(dǎo)期的影響因素進(jìn)行了實(shí)驗(yàn)研究。張冠敏等[3]從實(shí)驗(yàn)與理論兩個(gè)方面研究了板式換熱器顆粒污垢特性。李紅霞等[4]對(duì)強(qiáng)化管內(nèi)氧化鋁和氧化鐵顆粒污垢的污垢特性進(jìn)行了實(shí)驗(yàn)與模擬。朱華等[5]研究了螺紋管中實(shí)際冷卻水污垢和顆粒污垢的特性。王景濤等[6]從不同顆粒直徑、顆粒濃度、工質(zhì)流速和溫度工況這幾個(gè)角度對(duì)納米氧化鎂顆粒污垢在交叉縮放橢圓管內(nèi)的污垢特性進(jìn)行了實(shí)驗(yàn)研究。

渦流發(fā)生器作為一種被動(dòng)式強(qiáng)化換熱裝置已得到廣泛的應(yīng)用，它能產(chǎn)生渦旋，減薄邊界層，對(duì)污垢的沉積能起到一定的抑制作用。三角翼渦流發(fā)生器作為渦流發(fā)生器中應(yīng)用得比較廣泛的一種，在傳熱、空氣動(dòng)力學(xué)以及抑垢方面都得到了廣泛研究。王令等[7]、田林等[8]分別對(duì)矩形通道內(nèi)安裝縱向渦流發(fā)生器強(qiáng)化換熱進(jìn)行了模擬與實(shí)驗(yàn)研究。Kwak等[9]、Joardar等[10]分別對(duì)不同換熱器內(nèi)三角翼渦流發(fā)生器的換熱與壓降特性進(jìn)行了實(shí)驗(yàn)研究。Zhou等[11]在文獻(xiàn)中對(duì)矩形翼、三角翼、梯形翼和柱面梯形翼的流動(dòng)特性以及換熱特性進(jìn)行了論述。Hasan等[12]研究了三角翼渦流發(fā)生器在交叉流動(dòng)條件下Na2SO4析晶污垢的污垢特性，提出需要從抑垢與壓降兩個(gè)方面考慮渦流發(fā)生器的設(shè)計(jì)。劉坐東等[13]和王宇朋等[14]分別從模擬與實(shí)驗(yàn)兩個(gè)方面研究了3種翼形渦流發(fā)生器的CaCO3污垢特性，研究表明攻角的改變對(duì)于三角翼的抑垢變化影響不明顯，列間距是一個(gè)重要的影響因素。目前國(guó)內(nèi)外對(duì)于三角翼在顆粒污垢方面的研究較少，而顆粒污垢是污垢中主要成分之一，國(guó)內(nèi)外許多學(xué)者都從各個(gè)方面對(duì)其進(jìn)行了大量研究，三角翼目前應(yīng)用得也十分廣泛，因此研究三角翼的顆粒污垢特性是十分必要的。本研究采用顆粒粒徑為50 nm的納米氧化鎂顆粒配制純膠體溶液來(lái)模擬換熱面顆粒污垢的生成，以達(dá)到探索三角翼顆粒污垢特性的目的。

1 實(shí)驗(yàn)系統(tǒng)簡(jiǎn)介

1.1 實(shí)驗(yàn)系統(tǒng)與裝置

實(shí)驗(yàn)系統(tǒng)如圖1所示，工質(zhì)由水泵從低位水箱中打入高位水箱，一部分流入實(shí)驗(yàn)段，實(shí)驗(yàn)段為由厚0.5 mm、寬100 mm、長(zhǎng)1000 mm的304不銹鋼板制成的換熱面以及換熱面之間的密封材料構(gòu)成的一個(gè)長(zhǎng)1000 mm、寬100 mm、高8.5 mm的矩形通道。渦流發(fā)生器按一定尺寸布置在換熱面上。工質(zhì)在實(shí)驗(yàn)段中與恒溫水浴進(jìn)行換熱后流回低位水

圖1 實(shí)驗(yàn)系統(tǒng)

圖2 渦流發(fā)生器

箱進(jìn)行冷卻，另外一部分經(jīng)過(guò)溢流板溢流后直接 流入低位水箱。實(shí)驗(yàn)的溫度由PT100T型熱電偶 采集，其精度為±0.15℃。流量采用型號(hào)為L(zhǎng)DE-15SM2F100的電磁流量計(jì)采集。

實(shí)驗(yàn)采用的渦流發(fā)生器為高5 mm，高寬比為1:5的三角翼渦流發(fā)生器，材質(zhì)為0.3 mm厚的304不銹鋼片。渦流發(fā)生器如圖2所示。

1.2 實(shí)驗(yàn)原理

實(shí)驗(yàn)臺(tái)能夠直接得到實(shí)驗(yàn)段進(jìn)出水口的溫度1、2，恒溫水浴溫度0，以及實(shí)驗(yàn)中的流量v。

由傳熱方程式(1)以及熱平衡方程式(2)

可以得到傳熱系數(shù)的表達(dá)式(3)

污垢熱阻可以由式(4)得到

(4)

式中，f為污垢熱阻，m2·K·W-1；0、分別為清潔狀態(tài)下和有污垢狀態(tài)下實(shí)驗(yàn)段的總傳熱系數(shù)，W·m2·K-1。

2 實(shí)驗(yàn)結(jié)果分析

實(shí)驗(yàn)從水浴溫度、流速、納米氧化鎂溶液濃度、三角翼列間距以及不同布置方式這5個(gè)方面探究了三角翼渦流發(fā)生器的納米氧化鎂污垢特性。渦流發(fā)生器的布置方式如圖3所示：入口穩(wěn)定段150 mm，為出口穩(wěn)定段，前沿間距20 mm，為間距，迎流攻角90°，迎流向布置2列12排，箭頭方向?yàn)楣べ|(zhì)流動(dòng)方向。

2.1 水浴溫度對(duì)三角翼污垢特性的影響

在其他工況相同、不同水浴溫度的條件下對(duì)光板與三角翼渦流發(fā)生器進(jìn)行了實(shí)驗(yàn)。實(shí)驗(yàn)的水浴溫度分別為45℃和50℃。

圖4為實(shí)驗(yàn)得到的污垢熱阻曲線。從圖中可以

圖3 渦流發(fā)生器布置示意圖

圖4 水浴溫度對(duì)污垢熱阻的影響

看出，相同工況下加裝了三角翼渦流發(fā)生器的污垢熱阻總是小于光板，而隨著水浴溫度的升高，光板和三角翼的污垢熱阻都有所加大且結(jié)垢速率也明顯加快，光板（曲線A和B)的污垢熱阻漸近值增大了約58%，而三角翼只增大了約33%，三角翼（曲線C和D）的抑垢能力隨水浴溫度升高而增強(qiáng)。

相同工況下，三角翼渦流發(fā)生器具有明顯的抑垢效果，這是由于流體流經(jīng)渦流發(fā)生器會(huì)產(chǎn)生縱向渦，破壞邊界層，對(duì)壁面存在磨蝕和剪切作用，從而降低污垢沉積率，增大污垢的剝蝕率，導(dǎo)致污垢熱阻減小。隨著水浴溫度的升高，污垢熱阻漸近值和結(jié)垢速率也升高，這是由于水浴溫度的升高會(huì)引起換熱表面溫度的升高，使邊界層內(nèi)粒子的布朗運(yùn)動(dòng)加劇，這會(huì)導(dǎo)致更多的顆粒輸運(yùn)到壁面沉積，而使污垢熱阻與結(jié)垢速率升高。而加裝三角翼渦流發(fā)生器后，渦流發(fā)生器產(chǎn)生的縱向渦對(duì)壁面沖刷擾動(dòng)，破壞了壁面的邊界層，使壁面附近的流體處于不穩(wěn)定的狀態(tài)，在一定程度上抑制了顆粒向壁面的輸運(yùn)沉積，從而導(dǎo)致三角翼抑垢能力隨水浴溫度升高而增強(qiáng)。

2.2 濃度對(duì)三角翼污垢特性的影響

在其他工況不變的條件下分別對(duì)光板與三角翼進(jìn)行了不同污垢顆粒濃度的實(shí)驗(yàn)研究，實(shí)驗(yàn)采用的濃度為200、400 mg·L-1，實(shí)驗(yàn)得出的污垢熱阻曲線見圖5。從圖中可以看出，隨著濃度的升高，污垢熱阻明顯升高，而且結(jié)垢速率也明顯加快，三角翼渦流發(fā)生器的抑垢能力隨濃度上升而有所 減弱，400 mg·L-1下三角翼（曲線B）比光板（曲線D）的污垢熱阻漸近值減小了約11%，而200 mg·L-1下三角翼（曲線A）比光板（曲線C）減小了約53%。

圖5 濃度對(duì)污垢熱阻的影響

污垢熱阻以及結(jié)垢速率都隨濃度上升而升高，這是由于顆粒濃度越大溶液和換熱面的質(zhì)量濃度梯度就越大，會(huì)驅(qū)使更多的顆粒向換熱表面輸運(yùn)，從而導(dǎo)致結(jié)垢速率加快，污垢熱阻漸進(jìn)值增大[15]。隨濃度上升三角翼渦流發(fā)生器的抑垢能力減弱，這是由于工質(zhì)遇到渦流發(fā)生器后在其上部多形成橫向渦，導(dǎo)致上部形成一個(gè)回流滯止區(qū)，此區(qū)域污垢沉積率較大，而工質(zhì)濃度升高導(dǎo)致溶液中含有的納米氧化鎂粒子增加，相當(dāng)于加大了回流滯止區(qū)的沉積率，因而導(dǎo)致抑垢能力的減弱。

2.3 流速對(duì)三角翼污垢特性的影響

在其他工況不變的條件下對(duì)流體不同流速工況下的污垢特性進(jìn)行了實(shí)驗(yàn)研究，實(shí)驗(yàn)的流速分別為0.1、0.15 m·s-1，圖6為實(shí)驗(yàn)的污垢熱阻。從圖中可以看出，污垢熱阻隨流速增大而減小，而且結(jié)垢速率隨流速升高而降低，三角翼的抑垢能力隨流速升高而升高，0.1 m·s-1的流速下三角翼（曲線A）與光板（曲線B）相比污垢熱阻漸近值降低了約11%，而0.15 m·s-1的流速下三角翼（曲線C）比光板（曲線D）污垢熱阻減小了約32%。

流速對(duì)污垢形成的影響十分復(fù)雜，在污垢的輸運(yùn)、附著和剝蝕這3個(gè)階段流速都有參與。在污垢剛開始形成的階段溶液中污垢顆粒濃度較大，流速的增大一方面會(huì)為顆粒間的有效碰撞提供能量，導(dǎo)致更多的顆粒有機(jī)會(huì)相互碰撞團(tuán)聚[16]（膠體小顆粒之間發(fā)生碰撞相互結(jié)合成直徑更大的顆粒的過(guò)程），膠體粒子受到的重力影響可以忽略[17]，但其團(tuán)聚后由于粒徑增大會(huì)受到重力影響而沉降，沉降的大顆粒不易輸運(yùn)到壁面上，因而導(dǎo)致結(jié)垢速率減慢，污垢熱阻降低。另一方面流速的增大也會(huì)增大壁面受到的剪切力，導(dǎo)致污垢的剝蝕率升高，從而

圖6 流速對(duì)污垢熱阻的影響

使污垢熱阻減小。三角翼渦流發(fā)生器產(chǎn)生的縱向渦會(huì)隨流速增大而增強(qiáng)，其對(duì)壁面的沖刷與剪切效果也增強(qiáng)，從而導(dǎo)致三角翼的抑垢能力隨流速升高而增強(qiáng)。

2.4 均勻分布布置不同間距對(duì)三角翼污垢特性的影響

其他工況不變，在進(jìn)出口穩(wěn)定段近似相等的條件下，采用均勻分布的方式分別將間距為40、60、80 mm的渦流發(fā)生器布置在換熱面上，它們的排數(shù)分別為18排、12排、9排，對(duì)這3種布置的渦流發(fā)生器進(jìn)行實(shí)驗(yàn)。圖7為實(shí)驗(yàn)的污垢熱阻。從圖中可知，均勻分布布置下60 mm（曲線C）間距的污垢熱阻最大，80 mm（曲線B）次之，40 mm（曲線A）最小。

圖7 均勻分布下間距對(duì)污垢熱阻的影響

這可以從兩個(gè)方面進(jìn)行分析。一方面是由于三角翼渦流發(fā)生器產(chǎn)生渦旋的傳播距離是有限的，當(dāng)間距為40 mm時(shí)，三角翼產(chǎn)生的渦旋能夠傳遞到下一個(gè)渦流發(fā)生器，從而使渦旋能不斷傳播下去，加劇了對(duì)流體的擾動(dòng)，加強(qiáng)了對(duì)壁面的沖刷能力，導(dǎo)致污垢熱阻減??；當(dāng)間距在60 mm時(shí)，上一個(gè)渦流發(fā)生器產(chǎn)生的縱向渦處于即將平復(fù)的狀態(tài)，而流體流經(jīng)渦流發(fā)生器在迎流方向通常形成橫向渦，產(chǎn)生的縱向渦不能傳遞到下一個(gè)渦流發(fā)生器，而是在之前與下一個(gè)渦流發(fā)生器產(chǎn)生的橫向渦相遇，互相抵消，加大了渦流發(fā)生器前的滯止區(qū)，從而導(dǎo)致污垢熱阻升高；而當(dāng)間距為80 mm時(shí)，上一個(gè)三角翼產(chǎn)生的渦旋已經(jīng)完全平復(fù)，間距成為次要影響因素。另一方面，若在一定尺寸的換熱面上均勻分布渦流發(fā)生器，間距的減小必然導(dǎo)致排數(shù)的增加，從而增加了縱向渦對(duì)壁面的影響面積，有利于減小污垢熱阻。

圖8 相同排數(shù)下間距對(duì)污垢熱阻的影響

2.5 同排數(shù)布置不同間距對(duì)三角翼污垢特性的影響

在其他工況不變且不考慮均勻分布即忽略出口穩(wěn)定段的條件下，采用相同排數(shù)（9排）對(duì)3種間距進(jìn)行了實(shí)驗(yàn)，結(jié)果如圖8所示。從圖中可以發(fā)現(xiàn)，相同列數(shù)的工況下，80 mm（曲線A）間距污垢熱阻最小，40 mm（曲線B）次之，60 mm（曲線C）最大。

相同排數(shù)布置工況下，40 mm間距的污垢熱阻略大于80 mm間距，這是由于雖然40 mm間距下產(chǎn)生的縱向渦抑垢能力強(qiáng)，但是在換熱面面積足夠大的情況下，采用相同排數(shù)布置，40 mm間距布置下的渦流發(fā)生器產(chǎn)生的縱向渦對(duì)整個(gè)換熱面的影響范圍整體來(lái)說(shuō)卻比80 mm間距小許多，所以導(dǎo)致其污垢熱阻略大于80 mm。而對(duì)于60 mm間距，在這種布置情況下，其產(chǎn)生縱向渦的抑垢效果差，對(duì)整個(gè)換熱面的影響范圍也小，所以其抑垢能力仍然是最差的。

3 結(jié) 論

（1）換熱表面加裝三角翼渦流發(fā)生器對(duì)納米氧化鎂污垢存在抑垢作用，水浴溫度的升高會(huì)導(dǎo)致污垢熱阻升高以及結(jié)垢速率加快，三角翼渦流發(fā)生器的抑垢能力隨之增強(qiáng)。

（2）納米氧化鎂濃度的升高會(huì)導(dǎo)致污垢熱阻升高，結(jié)垢速率加快，三角翼渦流發(fā)生器的抑垢能力隨之減弱。

（3）流速的升高會(huì)導(dǎo)致污垢熱阻降低，結(jié)垢速率也隨之降低，三角翼渦流發(fā)生器的抑垢能力隨之增強(qiáng)。

（4）在換熱面尺寸一定的條件下，采用均勻分布的布置形式，間距40 mm的排數(shù)最多，抑垢效果最好，間距80 mm次之，間距60 mm抑垢能力最差。

（5）當(dāng)換熱面足夠大時(shí)，采用相同列數(shù)的布置方式，間距80 mm的抑垢效果最好，間距40 mm次之，間距60 mm抑垢效果最差。

（6）三角翼渦流發(fā)生器設(shè)計(jì)尺寸可能對(duì)實(shí)驗(yàn)結(jié)果存在一定的影響，通過(guò)對(duì)其進(jìn)行優(yōu)化應(yīng)可以提高其抑垢性能。

符 號(hào) 說(shuō) 明

——渦流發(fā)生器間距，mm

——納米氧化鎂膠體溶液濃度，mg·L

，

——分別為清潔狀態(tài)下和有污垢狀態(tài)下實(shí)驗(yàn)段的總傳熱系數(shù)，W·m·K

——流量，L·min

——污垢熱阻，m·K·W

，，

——分別為實(shí)驗(yàn)段入口溫度、出口溫度以及恒溫水浴溫度，℃

——實(shí)驗(yàn)流速，m·s

References

[1] Yang Shanrang (楊善讓), Xu Zhiming (徐志明), Sun Lingfang (孫靈芳). Heat Exchange Equipment Fouling and the Countermeasures (換熱設(shè)備污垢與對(duì)策) [M]. 2nd ed. Beijing: Science Press, 2004: 15

[2] Gu Yemei (顧業(yè)梅), Li Yun (李云), Zhang Bingqiang (張兵強(qiáng)). Experimient research on effect factor for particulate fouling induction [J].(鍋爐制造), 2011 (2): 19-21

[3] Zhang Guanmin (張冠敏), Li Guanqiu (李冠球), Li Wei (李蔚), Huang Tao (黃濤), Ren Yucheng (任毓程). Experimental and theoretical investigations about particulate fouling in plate heat exchangers [J].(工程熱物理學(xué)報(bào)), 2013, 34 (9): 1715-1718

[4] Li Hongxia (李紅霞), Li Guanqiu (李冠球), Li Wei (李蔚). Analysis of in tubes particulate fouling characteristic [J].:(浙江大學(xué)學(xué)報(bào): 工學(xué)版), 2012, 46 (9): 1671-1677

[5] Zhu Hua (朱華), Zhuang Bo (莊博), Li Wei (李蔚), Zhong Changxiong (鐘昌雄). A study of practical cooling water and particulate fouling characteristics in spiral-grooved tubes [J].(熱能動(dòng)力工程), 2008, 23 (2): 165-169

[6] Xu Zhiming (徐志明), Wang Jingtao (王景濤), Wang Lei (王磊), Zhang Yilong (張一龍), Liu Zuodong (劉坐東), Jia Yuting (賈玉婷). Experimental analysis on the particulate fouling characteristics of the alternating elliptical axis tube [J].(化工進(jìn)展), 2014, 33(4): 831-836

[7] Wang Ling (王令), Chen Qiuyang (陳秋煬), Zhou Yangeng (周硯耕), Wang Qiuwang (王秋旺), Huang Yanping (黃彥平), Xiao Zejun (肖澤軍). Heat transfer enhancement in rectangular narrow channel with periodically mounted longitudinal vortex generators on one sidewall [J].(核動(dòng)力工程), 2005, 26 (4): 344-347

[8] Tian Lin (田林), Bai Wei (柏巍), Xue Shanhu (薛山虎), Huang Zipeng (黃自鵬), Wang Qiuwang (王秋旺). Numerieal study of influence of longitudinal vortex generator on flow and heat transfer in rectangular channel [J].(工程熱物理學(xué)報(bào)), 2013, 34 (2): 324-327

[9] Kwak K M, Torii K, Nishino K. Simultaneous heat transfer enhancement and pressure loss reduction for finned-tube bundles with the first or two transverse rows of built-in winglets [J]., 2005, 29 (25): 625-632

[10] Joardar A, Jacobi A M. Impact of leading edge delta-wing vortex generators on the thermal performance of a flat tube, louvered-fin compact heat exchanger [J]., 2005, 48: 1480-1493

[11] Zhou G, Ye Q. Experimental investigations of thermal and flow characteristics of curved trapezoidal winglet type vortex generators [J]., 2012, 37: 241-248

[12] Hasan B O, Nathan G J, Ashman P J,. The use of turbulence generators to mitigate crystallization fouling under cross flow conditions [J]., 2012, 288: 108-117

[13] Xu Zhiming (徐志明), Liu Zuodong (劉坐東), Zhang Yilong (張一龍), Wang Yupeng (王宇朋), Zhang Zhongbin (張仲彬). Numerical simulation about influence of several kinds of winglet vortex generator structures on CaCO3fouling characteristics [J].(化工學(xué)報(bào)), 2013, 64 (7): 2468-2473

[14] Xu Zhiming (徐志明), Wang Yupeng (王宇朋), Liu Zuodong (劉坐東), Zhang Yilong (張一龍). Experimental investigation on the characteristic of restraining fouling of wing-type vortex generators [J].(工程熱物理學(xué)報(bào)), 2012, 33 (11): 1936-1939

[15] Xu Zhiming (徐志明), Dong Bing (董兵), Du Xiangyun (杜祥云), Wang Binglin (王丙林). Experimental study on particle fouling in plate heat exchangers [J].(動(dòng)力工程學(xué)報(bào)), 2013, 33(7): 539-543

[16] Yuan Wenjun (袁文俊), Zhou Yongmin (周永敏). Reasons for aggregation of nanoparticles and solutions [J].(材料導(dǎo)報(bào)), 2008, 22 (12): 59-61

[17] Feng Xusheng (馮續(xù)勝), Liu Hongguo (劉洪國(guó)), Hao Jingcheng (郝京誠(chéng)). Colloid Chemistry (膠體化學(xué)) [M]. Beijing: Chemical Industry Press, 2005: 63-71

Characteristics of nanometer magnesia particulate fouling of delta wing vortex generator

XU Zhiming1, YANG Suwu1, ZHU Xinlong1, ZHANG Yilong2, LIU Zuodong2

(1Energy and Power Engineering Institute, Northeast Dianli University, Jilin 132012, Jilin, China;2North China Electric Power University, Beijing 102206, China)

In order to explore the characteristics of nanometer magnesia particle fouling in delta wing vortex generator, colloidal solutions made of magnesia particles of 50 nm were used in the experiments. The fouling characteristics of the delta wing vortex generator under different conditions were investigated, including water bath temperature, particle concentration, flow velocity, spacing and arrangement of delta wings. The results show that the delta wing vortex generator can inhibit the fouling. The longitudinal vortex and stagnation zone produced by delta wing vortex generator and the characteristics of magnesia colloidal solution have a great influence on the fouling characteristics. Higher water bath temperature, lower magnesia solution concentration, and higher flow velocity can reduce the fouling rate. The fouling characteristics also depend on the arrangement of delta wings, which results in different longitudinal vortexes. With the uniform layout of delta wings, 40 mm spacing is the best for inhibition of dirt. With the same number of columns, 80 mm spacing is the best.

delta wing; magnesia; fouling; fouling characteristics; spacing; colloid; nanoparticles

date: 2014-06-26.

YANG Suwu, 1012461@qq.com

10.11949/j.issn.0438-1157.20140962

TK 124

A

0438—1157（2015）01—0086—06

國(guó)家自然科學(xué)基金項(xiàng)目（51076205）。

2014-06-26收到初稿，2014-09-02收到修改稿。

聯(lián)系人：楊蘇武。第一作者：徐志明（1959—），男，教授。

supported by the National Natural Science Foundation of China (51076205).