李?翔,裴毅強,秦?靜, 2,王同金,趙樂文,詹樟松,鄭建軍,陳堂明
?
壁溫和機油油膜對噴霧撞壁后發(fā)展的影響
李?翔1,裴毅強1,秦?靜1, 2,王同金1,趙樂文1,詹樟松3,鄭建軍3,陳堂明3
(1. 天津大學(xué)內(nèi)燃機燃燒學(xué)國家重點實驗室,天津 300072;2. 天津大學(xué)內(nèi)燃機研究所,天津 300072;3. 重慶長安汽車股份有限公司動力研究院,重慶 401120)
為研究壁面溫度和機油油膜對缸內(nèi)直噴汽油(GDI)發(fā)動機噴霧撞壁后發(fā)展的影響,運用陰影法成像系統(tǒng)和相位多普勒激光測試系統(tǒng)(PDA)對噴霧撞壁后反彈高度、擴散距離以及特征點的粒徑粒速做了試驗研究.結(jié)果表明:壁溫的升高或機油油膜的存在對噴霧撞壁后反彈和擴散起明顯的促進作用.但是,當(dāng)壁面達到萊頓弗羅斯特效應(yīng)溫度后,撞壁噴霧的反彈高度和擴散距離又開始減?。跍氐纳呋驒C油油膜的存在,都會使噴霧軸線測點處法向速度向上的液滴比例增加,且噴霧邊緣測點處液滴的切向平均速度明顯增大.壁溫的升高使噴霧撞壁后液滴的粒徑分布向更小粒徑的方向偏移,液滴粒徑的算術(shù)平均值減小,但100,mm以上“超大液滴”的出現(xiàn)使軸線測點處液滴的索特平均直徑增大.機油油膜的存在使噴霧撞壁后液滴粒徑的算術(shù)平均值減小,但索特平均直徑增大.
缸內(nèi)直噴汽油發(fā)動機;噴霧撞壁;壁面溫度;機油油膜
缸內(nèi)直噴汽油(gasoline direct injection,GDI)發(fā)動機將燃油直接噴入氣缸,燃料的霧化、蒸發(fā)吸熱作用可使進氣得到冷卻,從而使發(fā)動機功率提高,在動力性和經(jīng)濟性方面比進氣燃油式噴射(port fuel injection,PFI)發(fā)動機表現(xiàn)出更大的優(yōu)勢和潛力.目前,GDI發(fā)動機已經(jīng)成為汽油機的研究熱點[1-3].
為滿足更嚴(yán)格的顆粒物排放法規(guī),GDI發(fā)動機采用較高的燃油噴射壓力,以獲得更小粒徑的液滴,從而使缸內(nèi)混合氣更加均勻,減少顆粒物生成[4-5].但是,提高噴油壓力會增加噴霧貫穿距,受制于缸徑尺寸、噴油定時等因素,較長的噴霧貫穿距易造成噴霧撞擊缸套或活塞頂?shù)默F(xiàn)象,這會使缸內(nèi)局部混合氣過濃,導(dǎo)致發(fā)動機顆粒物和碳?xì)渑欧叛杆偕撸視觿「滋妆诿娴臋C油稀釋,使機油黏度下降,潤滑性降低[6].
國內(nèi)外對于GDI發(fā)動機噴霧撞壁進行了大量的試驗研究,但主要集中在噴射策略、撞壁距離及撞壁角度對噴霧撞壁的影響方面[7-9].
不同的壁面條件對于噴霧撞壁后的發(fā)展有較大影響,關(guān)于壁面溫度和機油油膜對噴霧撞壁影響的研究較少.尤其缺乏以壁溫或機油油膜作為單一變量,進行其對撞壁噴霧后結(jié)構(gòu)發(fā)展和液滴粒徑粒速分布的影響規(guī)律的詳細(xì)試驗研究.
本文運用陰影法成像測試及PDA粒速粒徑測試相結(jié)合的手段,對293,K、333,K、373,K和473,K壁溫的干壁面及存在機油油膜的293,K濕壁面進行了試驗研究,為充分認(rèn)識壁溫和機油油膜對GDI噴油器噴霧撞壁后的宏觀結(jié)構(gòu)發(fā)展、特征點粒速粒徑分布提供了有價值的信息和指導(dǎo),為解決缸內(nèi)燃油撞壁導(dǎo)致的排放問題也提供了一定的理論參考.
1.1?噴油器和試驗環(huán)境參數(shù)
為了在陰影法光學(xué)試驗中得到更清晰的單束噴霧輪廓圖像,試驗中用工業(yè)修補劑將一個6孔GDI噴油器的其中5孔密封,改造為單孔噴油器.
提高噴射壓力會導(dǎo)致貫穿距的增加,目前GDI發(fā)動機噴射壓力最大值通常約15,MPa.噴油器改造前后,噴孔總面積的減小會影響噴霧出口速度,進而影響噴霧貫穿距[10].通過對比,改造后噴油器在7,MPa噴射壓力條件下與改造前噴油器15,MPa噴射壓力下其中單一油束的貫穿距、噴霧錐角基本吻合.
因此,本文噴油壓力定為7,MPa,噴油脈寬定為2.5,ms,試驗在293,K、0.1,MPa的常溫常壓環(huán)境條件下進行.按照SAE J2715標(biāo)準(zhǔn)[11],試驗燃料為正庚烷,并將其溫度控制在(294±2)K.為接近缸壁或活塞頂?shù)牟馁|(zhì),試驗選取表面粗糙度為0.4的可加熱鋁合金平板作為噴霧撞壁壁面,壁溫控制精度為±2,K.
1.2?陰影法光學(xué)試驗系統(tǒng)
噴霧撞壁后宏觀結(jié)構(gòu)獲得以陰影法原理為基礎(chǔ),可得到清晰的噴霧液相信息.如圖1所示,試驗用高速相機是Photron Fastcam SA1.1型相機,拍攝頻率選取10,000幀/s,分辨率為768×768像素.用Matlab軟件對圖片進行后處理,得到反彈高度與擴散距離,如圖2所示.
圖1?陰影法光學(xué)試驗系統(tǒng)示意
圖2?反彈高度和擴散距離示意
1.3?PDA試驗系統(tǒng)
圖3為PDA粒徑粒速測試系統(tǒng)示意.試驗過程中,476.5,nm、488.0,nm及514.5,nm 3種波長的激光由發(fā)射器照向噴霧,噴霧液滴將激光折射,接收器吸收折射光信號并記錄液滴的粒速粒徑信息.對于每個試驗工況,每進行一次噴霧后,就對壁面條件進行復(fù)原,以消除噴霧附壁油膜對壁面的影響,保持每次噴霧時壁面條件的一致性.
圖3?PDA試驗系統(tǒng)示意
為減少試驗的誤差,每個工況采集20,000個以上的液滴.對采集數(shù)據(jù)結(jié)果進行后處理,得到該工況下測試點液滴的粒速粒徑分布、索特平均直徑(32)與粒徑的算術(shù)平均值(10).試驗中將100,μm以上的液滴定義為“超大液滴”,因其粒徑遠大于普通液滴,它的出現(xiàn)會造成該測點32的數(shù)值明顯提高.
1.4?噴霧測量特征點的選取
如圖4所示,壁面位于噴孔正下方33,mm,且與噴霧中軸線垂直.在壁面正上方5,mm(即噴孔正下方28,mm)所在的平面選取PDA試驗特征測點和.點位于該平面與噴霧軸線相交處,從點沿平面水平向右8,mm選取點.點可直觀反映垂直壁面的噴霧軸線上近壁面處的液滴粒徑、粒速分布情況.點可直觀反映靠近噴霧邊緣近壁面處的液滴粒徑和粒速分布情況.
圖4?噴霧測量特征點A、B示意
2.1?壁面溫度和機油油膜對噴霧撞壁后宏觀結(jié)構(gòu)發(fā)展的影響
為盡可能排除測量誤差的影響,對每個工況點重復(fù)進行20次拍攝,并求平均值,下文圖5~圖8中用“I”的上、下沿分別表示該點在20次重復(fù)拍攝中數(shù)據(jù)的最大值和最小值.
壁面溫度對噴霧撞壁后液滴的反彈高度及擴散距離有明顯的影響.如圖5所示,在壁溫293,K、333,K、373,K條件下,噴霧反彈高度均在3.2,ms ASOI(噴油后時刻)達到最大值,分別為15.4,mm、22.8,mm和26.1,mm.之后,隨著液滴的氣化蒸發(fā)、環(huán)境氣體的阻力,噴霧高度開始緩慢下降.如圖6所示,噴霧擴散距離最大值由壁溫293,K時的43.3,mm增加到373,K時的63.4,mm,增大了46.42%,.噴霧撞擊壁面后,一部分液滴會發(fā)生反彈、破碎.在低于正庚烷沸點371.5,K的溫度范圍內(nèi),壁溫升高使液滴與壁面之間的熱交換增強,液滴蒸發(fā)速率明顯加快[12],蒸發(fā)產(chǎn)生向上的動力使液滴反彈能量增加.另外,壁面附近空氣受熱后流動性的增加也會加劇液滴的破碎趨勢,更小質(zhì)量的液滴在相同的動能作用下更易反彈和擴散.綜合上述作用,隨壁溫從293,K升高到373,K,撞壁噴霧的反彈高度和擴散距離明顯增大.
燃料的萊頓弗羅斯特效應(yīng)溫度L(K)與該燃料的沸點溫度B(K)的關(guān)系[12]為
???TL=1.06(TB-273)+325(1)
將常壓下正庚烷的沸點371.5,K代入式(1),得到其萊頓弗羅斯特效應(yīng)溫度為429.41,K.
因此,在473,K壁溫條件下,噴霧前端鋒面與壁面首先接觸的液滴會發(fā)生萊頓弗羅斯特效應(yīng),從而迅速沸騰形成蒸汽層浮于壁面表層,使其余大部分噴霧液滴避免了與高溫壁面直接接觸,液滴與壁面的熱交換作用減弱,蒸發(fā)速率減緩,反彈及擴散的初始動能大幅降低;而且,壁溫的進一步升高使得撞壁液滴蒸發(fā)氣化的比例提高,削弱了液相噴霧的反彈擴散趨勢.綜上,壁溫升高到473,K后,噴霧撞壁后的反彈高度和擴散距離有所減小,其最大值比壁溫373,K時分別減小了24.07%,和18.89%,,且反彈高度最大值的出現(xiàn)時刻提前至2.8,ms ASOI.
圖5?壁面溫度對噴霧撞壁后反彈高度的影響
圖6?壁面溫度對噴霧撞壁后擴散距離的影響
發(fā)動機內(nèi)流體潤滑的典型機油油膜厚度通常約為1~100,μm[13].因此,試驗中,濕壁面的機油油膜厚度定為60,μm.通過用微量天平測出平板涂機油前后的質(zhì)量差值Δ,用Δ除以壁面表面積得出油膜厚度.同時考慮本試驗中壁面粗糙度a為0.4,遠小于Δ/,故實際油膜厚度為
???h=Δm/S+Ra≈Δm/S(2)
如圖7和圖8所示,機油油膜對噴霧撞壁后反彈和擴散起明顯的促進作用.噴霧反彈高度最大值從15.4,mm增至27.0,mm,增大了75.32%,,而擴散距離最大值從38.1,mm增至51.5,mm,增大了35.17%,.如圖9[14]所示,噴霧單個液滴撞擊機油油膜后,將攤開擴散形成“皇冠狀”的腔體,腔體邊緣破碎后形成小液滴,其中夾帶著一部分機油.這數(shù)量巨大的具有反彈和擴散方向初始動能的小液滴,使得噴霧反彈和擴散的趨勢明顯增強.
圖7?機油油膜對噴霧撞壁后反彈高度的影響
圖8?機油油膜對噴霧撞壁后擴散距離的影響
2.2?壁面溫度和機油油膜對噴霧特征點速度的影響
圖10~圖13展示了壁面溫度和機油油膜對點液滴法向速度分布及點液滴切向平均速度的影響.圖10和圖12中設(shè)垂直向下速度方向為正,向上為負(fù).圖11和圖13中設(shè)平行向右速度方向為正,向左為負(fù).
圖9?“皇冠狀”腔體及破碎飛濺
圖10壁面溫度對點A處液滴法向速度概率分布的影響 Fig.10Effect of wall temperature on the normal velocity distribution of point A droplets
圖11?壁面溫度對點B處液滴切向平均速度的影響
由圖10可知,在293,K、333,K、373,K、473,K壁面條件下,點液滴法向速度方向向上所占的比例分別為3.5%,、16.8%,、53.2%,和58.8%,,即隨壁溫的升高,速度方向向上的液滴比例明顯增加.
壁溫升高至373,K和473,K時,速度分布曲線呈現(xiàn)單峰狀態(tài),峰值粒速在-5~0,m/s之間.這是由于較高的壁溫加劇了液滴撞壁后的破碎反彈趨勢,產(chǎn)生更多向上反彈的小液滴.液滴反彈至壁面上方5,mm測點的過程中,伴隨著氣化揮發(fā),且需要克服空氣阻力和自身重力,到達測點時的動能已急劇減少,所以其向上速度接近0,m/s.
在293,K與333,K壁溫條件時,點液滴粒速呈雙峰分布,向下速度較大(大于30,m/s)的粒徑所占比例分別高達69.2%,和50.5%,.
由圖11可知,點的切向平均速度隨壁溫升高而增大.但是,從373,K升高到473,K,速度僅增加了0.5,m/s.液滴切向速度的小幅增加有助于噴霧擴散距離的增大,但難以抵消萊頓弗羅斯特效應(yīng)和液滴高溫氣化對液滴擴散距離的抑制作用.所以,相比于壁溫373,K條件,壁溫473,K時噴霧撞壁后擴散距離明顯減小,如圖6所示.
圖12?機油油膜對點A處液滴法向速度的影響
圖13?機油油膜對點B處液滴切向平均速度的影響
由圖12和圖13可知,293,K干壁面、濕壁面條件下,點液滴法向速度方向向上的比例分別為3.5%,、10.5%,.對于點,293,K濕壁面比干壁面的切向平均速度大0.5,m/s.這是因為噴霧撞擊濕壁面機油油膜后,形成“皇冠狀”腔體,腔體邊緣破碎形成數(shù)量巨大的小液滴,使經(jīng)過點、分別具有垂直向上、水平向右速度的液滴數(shù)量明顯增多.
2.3?壁面溫度和機油油膜對噴霧特征點粒徑的影響
由圖14~圖17可知,隨著壁溫的升高,點液滴的粒徑分布向更小粒徑方向整體偏移,5,μm以下粒徑液滴的比例大幅增加,且2,μm 左右粒徑的峰值明顯上升.這是由于壁溫升高加劇了液滴撞壁后的破碎反彈趨勢,產(chǎn)生更多反彈向上的小液滴.在發(fā)動機工作過程中,小粒徑液滴比例的增加有利于提升缸內(nèi)油氣混合的均勻性,減少顆粒物的生成.
圖14?在293,K壁溫下點A處液滴的粒徑分布
圖15?在333,K壁溫下點A處液滴的粒徑分布
圖16?在373,K壁溫下點A處液滴的粒徑分布
圖17?在473,K壁溫下點A處液滴的粒徑分布
然而,由圖14~圖17可知,100,μm以上“超大液滴”的比例隨壁溫上升出現(xiàn)了逐步增長的現(xiàn)象.“超大液滴”的增加一方面是因為隨壁溫上升,反彈小液滴數(shù)量增加使無規(guī)則運動造成彼此的碰撞合并現(xiàn)象增加,另一方面是因為軸線處的噴霧濃度較高,入射方向的液滴較密集,反彈小液滴的增多會增加其與入射方向液滴發(fā)生碰撞合并的概率.在發(fā)動機工作過程中,“超大液滴”的出現(xiàn)易造成局部混合氣濃度過高,會促進碳煙的生成,對降低發(fā)動機顆粒物排放是不利的.
因此,在圖18中,小液滴數(shù)量比例的增加使10隨壁溫升高而逐步減小,從7.0,μm下降至3.2,μm,減小了54.29%,.但是,隨壁溫升高,“超大液滴”的明顯增加導(dǎo)致32從33.6,μm上升至 43.9,μm,增大了30.95%,.
圖18?壁溫對點A處液滴的D32和D10的影響
由圖19~圖22可知,隨壁溫升高,液滴撞壁后的破碎反彈趨勢加劇,產(chǎn)生更多反彈小液滴,點液滴粒徑分布向更小粒徑的方向整體偏移,2,μm左右粒徑的峰值明顯上升.
然而,點與噴霧軸線點明顯不同的是,點的“超大液滴”比例隨壁溫上升而逐步減少.這是因為點位于靠近噴霧油束邊緣的近壁面位置,噴霧撞壁后液滴的初始動能小于噴霧軸線點處,液滴撞壁后的反彈破碎趨勢減弱.因此,在壁溫293,K和333,K條件下,噴霧撞壁后存在少量未充分破碎的較大反彈液滴.壁溫升高至373,K和473,K后,噴霧撞壁后的破碎現(xiàn)象加強,但噴霧邊緣處反彈液滴的數(shù)量濃度明顯少于噴霧軸線點處,液滴發(fā)生碰撞合并的概率較低,難以形成“超大液滴”.
圖19?在293,K壁溫下點B處液滴的粒徑分布
圖20?在333,K壁溫下點B處液滴的粒徑分布
圖21?在373,K壁溫下點B處液滴的粒徑分布
圖22?在473,K壁溫下點B處液滴的粒徑分布
因此,在圖23中,10和32都隨壁溫的升高而呈現(xiàn)逐步減小的趨勢.由圖24、圖25分別與圖14、圖19的對比可知,壁面存在機油油膜的條件下,點、的粒徑分布均向更小粒徑的方向偏移,2,μm左右小粒徑液滴的比例明顯增加.因此,圖26和圖27中顯示出濕壁面條件下,點、的10的數(shù)值比干壁面條件下分別減小了12.35%,和6.93%.
圖23?壁溫對點B處的D32和D10的影響
圖24?在293,K濕壁條件下點A處的粒徑分布
圖25?在293,K濕壁條件下點B處的粒徑分布
圖26?機油油膜對點A處的D32和D10的影響
圖27?機油油膜對點B處的D32和D10的影響
相比于293,K干壁面,噴霧撞擊機油油膜濕壁面后的反彈飛濺液滴中會夾帶著一部分粒徑較大的機油油滴,造成點、處的100,μm以上“超大液滴”的比例有所增加.這就導(dǎo)致了圖26和圖27中濕壁面條件下,點、的32數(shù)值比干壁面分別增大了18.71%,和4.37%,.
(1) 相比于常溫壁面,在熱壁面條件下,撞壁噴霧的反彈高度和擴散距離明顯增大.但是,當(dāng)壁溫達到473,K,撞壁噴霧的反彈高度和擴散距離又開始減小,這與壁面產(chǎn)生萊頓弗羅斯特效應(yīng)和液滴高溫氣化有關(guān).機油油膜的存在使噴霧撞壁后反彈高度和擴散距離明顯增加.
(2) 壁溫的升高或者機油油膜的存在,都會使噴霧軸線測點處法向速度向上的液滴比例增加,且噴霧邊緣測點處的液滴切向平均速度明顯增大.
(3) 壁溫的升高使噴霧撞壁后液滴的粒徑分布向更小粒徑的方向偏移,液滴粒徑的算術(shù)平均值減?。?,“超大液滴”的出現(xiàn)使噴霧撞壁后軸線測點處液滴的索特平均直徑增大.
(4) 機油油膜的存在使噴霧撞壁后液滴粒徑的算術(shù)平均值減小,但飛濺的機油使“超大液滴”比例增加,造成噴霧撞壁后液滴的索特平均直徑增大.
[1] Chen L F,Stone R,Richardson D. Effect of the valve timing and the coolant temperature on particulate emissions from a gasoline direct injection engine fuelled with gasoline and with a gasoline-ethanol blend[J].:,2012,226(10):1419-1430.
[2] Park C,Kim S,Kim H,et al. Stratified lean combustion characteristics of a spray-guided combustion system in a gasoline direct injection engine[J].,2012,41(1):401-407.
[3] 馬?驍,何?旭,王建昕,等. 用激光誘導(dǎo)熒光法測量GDI發(fā)動機缸內(nèi)混合氣體分布[J]. 內(nèi)燃機工程,2010,31(4):1-5.
Ma Xiao,He Xu,Wang Jianxin,et al. In-cylinder mixture distribution measurement in a GDI engine using laser-introduced fluorescence[J].,2010,31(4):1-5(in Chinese).
[4] Sabathil D,Koenigstein A,Schaffner P,et al. The influence of DISI engine operating parameters on particle number emissions[C]//. Detroit,USA,2011:2011-01-0143.
[5] He X,Ratcliff M A,Zigler B T. Effects of gasoline direct injection engine operating parameters on particle number emissions[J].,2012,26(4):2014-2027.
[6] Kapus P,Ogris M,Hollerer P. Measures to reduce particulate emissions from gasoline DI engines[C]//. Detroit,USA,2011:2011-01-1219.
[7] Gold M,Stokes J,Morgan R,et al. Air-fuel mixing in a homogeneous charge DI gasoline engin[C]//. Detroit,USA,2001:2001-01-0968.
[8] 范錢旺,高?雅,董戰(zhàn)力,等. 非對稱多孔噴油器撞壁噴霧特性的試驗與模擬[J]. 內(nèi)燃機學(xué)報,2012,30(1):35-41.
Fan Qianwang,Gao Ya,Dong Zhanli,et al. Experiment and simulation on impinging spray characteristics of an asymmetrical multi-hole injector[J].,2012,30(1):35-41(in Chinese).
[9] 李相超,張玉銀,許?敏,等. 直噴汽油機缸內(nèi)噴霧濕壁問題研究[J]. 內(nèi)燃機工程,2012,33(5):17-23.
Li Xiangchao,Zhang Yuyin,Xu Min,et al. Investigation on wall wetting for GDI engine[J].,2012,33(5):17-23(in Chinese).
[10] 高?劍,蔣德明,廖世勇,等. 缸內(nèi)直噴汽油機高壓渦旋噴油器的油束模型[J]. 西安交通大學(xué)學(xué)報,2003,37(9):898-902.
Gao Jian,Jiang Deming,Liao Shiyong,et al. Spray model of high pressure swirl injector in the gasoline direct injection engines[J].,2003,37(9):898-902(in Chinese).
[11] Hung D,Hamngtom D L,Gandhi A H,et al. Gasoline fuel injector spray measurement and characterization—A new SAE J2715 recommended practice[C]//. Detroit,USA,2008:2008-01-1068.
[12] Mills A A,F(xiàn)ry J D. Rate of evaporation of hydrocarbons from a hot surface:Nukiyama and Leidenfrost temperatures[J].,1982,3:152-154.
[13] 武?彬. 汽油機低摩擦系統(tǒng)優(yōu)化及其對節(jié)能影響研究[D]. 長春:吉林大學(xué)汽車工程學(xué)院,2013.
Wu Bin. Study on Low Friction System Optimization of Gasoline Engine and Its Effects on Energy-Saving[D]. Changchun:College of Auto Engineering,Jilin University,2013(in Chinese).
[14] Cossali G E,Coghe A,Marengo M. The impact of a single drop on a wetted solid surface[J].,1997,22:463-472.
(責(zé)任編輯:孫立華)
Effects of Wall Temperature and Oil Film on the Process of a Gasoline Spray After Wall Impingement
Li Xiang1,Pei Yiqiang1,Qin Jing1, 2,Wang Tongjin1,Zhao Lewen1,Zhan Zhangsong3,Zheng Jianjun3,Chen Tangming3
(1.State Key Laboratory of Engines,Tianjin University,Tianjin 300072,China;2.Internal Combustion Engine Research Institute,Tianjin University,Tianjin 300072,China;3.Powertrain Engineering R&D Center,Chongqing Chang’an Automobile Limited Company,Chongqing 401120,China)
The aim of this experimental research is to study the effects of wall temperature and oil film on the process of a gasoline spray after wall impingement.The research focuses on the spray rebound height,diffusion distance,droplet size and velocity distribution after wall impingement using shadow photographic method and phase Doppler analyzer(PDA).The results show that the wall temperature or oil film plays a significant role in raising the spray rebound height and diffusion distance.But when the wall temperature is above the Leidenfrost temperature,the rebound height and diffusion distance begin to decrease.The normal and tangential velocities of droplet increases with the rise of wall temperature or the existence of oil film.As the wall temperature increases,the drop size distribution is towards smaller diameters and the average diameter of droplets decreases.Nonetheless,due to “huge droplet” whose diameter is bigger than 100,μm,the normal direction droplets’ Sauter mean diameter increases.The oil film has an influence on the droplet distribution,with the droplets’ average diameter decreasing while their Sauter mean diameter increasing.
gasoline direct injection(GDI)engine;spray impingement;wall temperature;oil film
10.11784/tdxbz201512065
TK411
A
0493-2137(2016)11-1195-08
2015-12-18;
2016-04-06.
國家科技支撐計劃資助項目(2014BAG10B01).
李?翔(1988—??),男,博士研究生,qilongzhu007@hotmail.com.
裴毅強,peiyq@tju.edu.cn.