胡超群,朱 亮,2,沙作良,王彥飛,楊立斌,趙楠楠
(1.天津科技大學(xué)化工與材料學(xué)院,天津300457;2.邁安德集團(tuán)有限公司)
雜質(zhì)離子對(duì)共晶冷凍結(jié)晶過程中冰晶結(jié)晶特性的影響*
胡超群1,朱 亮1,2,沙作良1,王彥飛1,楊立斌1,趙楠楠1
(1.天津科技大學(xué)化工與材料學(xué)院,天津300457;2.邁安德集團(tuán)有限公司)
以硝酸鉀水溶液為研究對(duì)象,考察了雜質(zhì)離子對(duì)共晶冷凍結(jié)晶過程中冰晶結(jié)晶特性的影響。首先,建立并驗(yàn)證了共晶冷凍結(jié)晶過程冰點(diǎn)的測(cè)定方法,并基于Extend-UNIQUAN模型對(duì)硝酸鉀體系的冰點(diǎn)進(jìn)行模擬計(jì)算。同時(shí)考察了不同種類雜質(zhì)離子對(duì)硝酸鉀水溶液冰點(diǎn)、結(jié)晶點(diǎn)、以及冷凍介穩(wěn)區(qū)寬度的影響,實(shí)驗(yàn)研究結(jié)果表明:硝酸鉀水溶液中離子雜質(zhì)濃度越高,冰點(diǎn)、結(jié)晶點(diǎn)的溫度越低,介穩(wěn)區(qū)寬度越寬;雜質(zhì)離子濃度相同的硝酸鉀水溶液其冰點(diǎn)溫度隨著雜質(zhì)離子半徑減小而減小,結(jié)晶點(diǎn)溫度和介穩(wěn)區(qū)寬度與陽離子雜質(zhì)價(jià)態(tài)成反比,與陰離子雜質(zhì)價(jià)態(tài)成正比。并利用基于經(jīng)典成核理論所推導(dǎo)的介穩(wěn)區(qū)寬度模型方程對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行模擬計(jì)算,模擬結(jié)果與實(shí)驗(yàn)值符合較好,計(jì)算得到硝酸鉀水溶液的冰晶成核級(jí)數(shù)為1.46。
冰晶;介穩(wěn)區(qū);硝酸鉀;共晶冷凍結(jié)晶
共晶冷凍結(jié)晶是利用冰與水溶液之間的固液相平衡原理,將水以冰晶的形式析出并將其從液相中分離出去,從而獲得純凈的冰晶和鹽的一種結(jié)晶技術(shù)。相較于其他的結(jié)晶手段,該技術(shù)可在低溫下操作,氣液界面小,此外還具有微生物繁殖、溶質(zhì)的劣化及揮發(fā)性芳香成分損失可控制在極低水平等優(yōu)點(diǎn)。特別適用于濃縮熱敏性液態(tài)食品、生物制藥、要求保留天然色香味的高檔飲品及中藥湯劑等。水凍結(jié)所需要熱量為355 kJ/kg,100℃和0℃時(shí)水蒸發(fā)所需的熱量為2 248 kJ/kg和2 495 kJ/kg,由此可見冷凍結(jié)晶消耗的能量約為蒸發(fā)法的0.15倍。因此,從理論上講冷凍結(jié)晶是一種節(jié)省能量的結(jié)晶手段[1-2]。
共晶冷凍結(jié)晶已成功應(yīng)用于廢水處理工藝并取得許多進(jìn)展,其中廢水中的雜質(zhì)離子對(duì)體系的冰點(diǎn)、介穩(wěn)區(qū)寬度均有很大影響。冰點(diǎn)是溶液的一個(gè)重要屬性,對(duì)處理廢水工藝具有指導(dǎo)意義,介穩(wěn)區(qū)是結(jié)晶過程控制的重要環(huán)節(jié)。在工業(yè)結(jié)晶的實(shí)際操作中,過飽和度通常被控制在介穩(wěn)區(qū)范圍內(nèi)進(jìn)行,這樣能使得到的產(chǎn)品晶形好且粒度均勻。體系介穩(wěn)區(qū)的寬度是選擇合適的過飽和度操作條件的重要依據(jù),還是結(jié)晶器的設(shè)計(jì)的關(guān)鍵要素[3]。因此,探究雜質(zhì)離子對(duì)冰點(diǎn)、介穩(wěn)區(qū)寬度的影響具有重要的意義。
筆者以硝酸鉀為原料,考察了共晶冷凍結(jié)晶過程冰點(diǎn)的測(cè)定方法及不同種類的雜質(zhì)離子對(duì)冰點(diǎn)、結(jié)晶點(diǎn)、介穩(wěn)區(qū)寬度的影響,從而確定共晶冷凍結(jié)晶過程中的一些熱力學(xué)及動(dòng)力學(xué)參數(shù),該參數(shù)可用于共晶冷凍結(jié)晶過程中的控制。在實(shí)際生產(chǎn)中,可用于改善操作條件、減低成本、提高產(chǎn)品質(zhì)量。
1.1 原料和儀器
原料與試劑:硝酸鉀、硝酸鎂、硝酸鋁、氯化鉀、硫酸鉀,均為分析純,博歐特(天津)化工貿(mào)易有限公司提供;實(shí)驗(yàn)用水:去離子水(天津盛騰達(dá)化學(xué)試劑公司)。
儀器:500 mL結(jié)晶器、78-1型數(shù)顯攪拌器、LAUDA-RP845型循環(huán)水浴。冷凍結(jié)晶的實(shí)驗(yàn)裝置見圖1。
圖1 共晶冷凍結(jié)晶的實(shí)驗(yàn)裝置
1.2 理論分析
1.2.1 硝酸鉀水溶液冰點(diǎn)測(cè)定原理
冰點(diǎn)是水和冰平衡共存時(shí)的溫度,在硝酸鉀水溶液的冷凍過程中,當(dāng)溫度達(dá)到或低于冰點(diǎn)溫度時(shí)溶液卻不結(jié)冰的現(xiàn)象稱之為過冷現(xiàn)象。處于過冷現(xiàn)象的溶液極不穩(wěn)定,任一擾動(dòng)便可觸發(fā)其結(jié)晶。過冷的硝酸鉀水溶液突然爆發(fā)成核生成冰晶,體系由于發(fā)生相變而釋放潛熱,使溶液的溫度回升?;厣聊骋稽c(diǎn)時(shí),溫度有短暫的平衡,此后溶液溫度將繼續(xù)下降,這段較平穩(wěn)的溫度即為溶液的冰點(diǎn)。
1.2.2 硝酸鉀水溶液冰點(diǎn)的計(jì)算
根據(jù)物理化學(xué)的基本原理,溶液的冰點(diǎn)與溶液中溶劑的活度[4]有關(guān):
在Extend-UNIQUAC模型中,組分i的活度系數(shù)γi由組合項(xiàng)、剩余項(xiàng)及德拜-休克爾項(xiàng)3部分組成,分別反映分子大小和形狀、分子間的交互作用及靜電作用對(duì)γ的貢獻(xiàn)。
其中,對(duì)于離子態(tài)的靜電長(zhǎng)程作用項(xiàng)[6]為:
運(yùn)用Gibbs-Duhem方程,得到溶劑的活度系數(shù):
式中,Ms為溶劑的相對(duì)分子質(zhì)量;Aφ為 Debye-Hückel常數(shù);ρ為方法參數(shù);是以物質(zhì)的量分?jǐn)?shù)為標(biāo)度的離子強(qiáng)度;Zi為溶質(zhì)離子的電荷數(shù)。其中Aφ、ρ可查文獻(xiàn)[7]。
組合項(xiàng)和剩余項(xiàng)分別為:
其中:
式中,θi和φi是物質(zhì)i的平均面積分?jǐn)?shù)和體積分?jǐn)?shù);ri和qi的數(shù)值根據(jù)分子的vanderWalls體積和表面積算出[8];Z為晶格配位數(shù),其值取10;uij是分子i-j的相互作用能,uij=uji。 由式(2)得出水的活度系數(shù),可由式αi=γixi計(jì)算出水的活度,再代入式(1)中即可計(jì)算出溶液的冰點(diǎn)。
1.3 實(shí)驗(yàn)步驟
1.3.1 硝酸鉀水溶液冰點(diǎn)測(cè)定過程
開啟恒溫裝置,調(diào)至預(yù)定溫度恒溫。向結(jié)晶器中加入200.00 g蒸餾水,再加入10.53 g的硝酸鉀,在400 r/min的速率下充分?jǐn)嚢枰源_保硝酸鉀完全溶解,隨后先快速將溫度降至0℃,再以0.05℃/min的降溫速率將溶液從0℃開始冷凍。當(dāng)溶液開始析出冰時(shí),記錄因冰析出釋放熱量而使溶液溫度回升至最高點(diǎn)時(shí)的溫度,此后用帶有濾嘴的注射器抽取上清液,用烘干法計(jì)算溶液回升至最高點(diǎn)時(shí)硝酸鉀的濃度,進(jìn)而得到在此濃度下硝酸鉀溶液的冰點(diǎn)。
1.3.2 硝酸鉀水溶液介穩(wěn)區(qū)測(cè)定過程
開啟恒溫裝置,調(diào)至預(yù)定溫度恒溫。向結(jié)晶器中加入200.00 g蒸餾水,再加入10.53 g的硝酸鉀,在400 r/min的速率下充分?jǐn)嚢枰源_保硝酸鉀完全溶解,隨后先快速將溫度降至0℃,再以0.05℃/min的降溫速率將溶液從0℃開始冷凍。當(dāng)溶液開始析出冰晶時(shí),記錄結(jié)晶點(diǎn)的溫度,重復(fù)實(shí)驗(yàn)3次,記錄實(shí)驗(yàn)數(shù)據(jù)。
2.1 硝酸鉀水溶液冰點(diǎn)的測(cè)定結(jié)果與分析
實(shí)驗(yàn)考察了硝酸鉀水溶液冰點(diǎn)的測(cè)定方法,將實(shí)驗(yàn)測(cè)定的冰點(diǎn)數(shù)據(jù)與文獻(xiàn)提出的數(shù)據(jù)通過Extend-UNIQUAN模型對(duì)比,以確定測(cè)定冰點(diǎn)的方法是否可行。
圖2為質(zhì)量分?jǐn)?shù)為5%的硝酸鉀溶液冷凍結(jié)晶過程中溫度隨時(shí)間變化的曲線。由圖2可見,當(dāng)溫度降至結(jié)晶點(diǎn)溫度-5.61℃時(shí),溶液中開始有冰晶析出,由于冰析出釋放熱量使溶液的溫度回升至-1.45℃,在此溫度下維持短暫的恒定后繼續(xù)下降,得到平衡點(diǎn)時(shí)硝酸鉀的質(zhì)量分?jǐn)?shù)為5.44%(表1);配制質(zhì)量分?jǐn)?shù)為5.44%的硝酸鉀溶液,將溫度恒定至-1.45℃時(shí),水和冰平衡共存,溶液濃度恒定不變,通過驗(yàn)證得出5.44%硝酸鉀溶液的冰點(diǎn)為-1.45℃,與文獻(xiàn)[9]中5.0%的硝酸鉀溶液的冰點(diǎn)-1.4℃幾乎一致(見表2),由此可認(rèn)為該冰點(diǎn)測(cè)定方法是可行的。
圖2 5.0%的硝酸鉀水溶液共晶冷凍結(jié)晶過程中溫度隨時(shí)間變化曲線
表1 KNO3-H2O體系相平衡和介穩(wěn)區(qū)數(shù)據(jù)
表2 KNO3-H2O體系冰點(diǎn)的實(shí)驗(yàn)值與文獻(xiàn)值的對(duì)比
運(yùn)用上述測(cè)定冰點(diǎn)的方法,得到不同質(zhì)量分?jǐn)?shù)(w)的硝酸鉀水溶液的冰點(diǎn)(T),并將其與Extend-UNIQUAN模型計(jì)算出硝酸鉀溶液的冰點(diǎn)對(duì)比,結(jié)果見圖3。由圖3可見,運(yùn)用Extend-UNIQUAN模型計(jì)算得到的冰點(diǎn)值與本研究提出的測(cè)定冰點(diǎn)方法得到的冰點(diǎn)值基本一致,進(jìn)一步驗(yàn)證了測(cè)定方法的可行性。將冰點(diǎn)溫度與濃度數(shù)據(jù)經(jīng)二次多項(xiàng)式擬合得到T=-209.376 4w2+0.808 0w-0.801 1,R2=0.991 5。
圖3 硝酸鉀水溶液冰點(diǎn)模型計(jì)算值與實(shí)驗(yàn)值對(duì)比
2.2 含不同雜質(zhì)離子的硝酸鉀水溶液冰點(diǎn)的測(cè)定結(jié)果與分析
為考察不同雜質(zhì)離子對(duì)硝酸鉀水溶液的冰點(diǎn)溫度的影響,先配制質(zhì)量分?jǐn)?shù)為5.44%的硝酸鉀溶液,向溶液中添加不同濃度的Na+、Mg2+、Al3+等作為陽離子雜質(zhì)和不同濃度的Cl-、SO42-等作為陰離子雜質(zhì),設(shè)定攪拌速率為400 r/min,確定溶液充分混合后,以0.05℃/min的降溫速率將溶液從0℃開始冷凍,在冷凍結(jié)晶過程中記錄冰點(diǎn)溫度,同時(shí)利用Extend-UNIQUAN模型對(duì)體系的冰點(diǎn)進(jìn)行計(jì)算預(yù)測(cè),結(jié)果見圖4。由圖4結(jié)果可知,運(yùn)用Extend-UNIQUAN模型計(jì)算的冰點(diǎn)值的平均相對(duì)偏差為5%左右,在低濃度時(shí)其計(jì)算值與實(shí)驗(yàn)值更符合,高濃度時(shí)偏差略大。從圖4還可見,同種離子的濃度越高,溶液的冰點(diǎn)溫度越低;離子半徑越小,冰點(diǎn)溫度越低。因?yàn)樗闲?yīng)對(duì)離子半徑較小、電荷較高的溶液體系影響更大,因離子半徑r(Al3+)<r(Mg2+)<r(Na+)<r(Cl-)<r(SO42-),所以在摻雜的離子雜質(zhì)中,摻加Al3+的硝酸鉀水溶液的冰點(diǎn)值最低。
圖4 含不同離子雜質(zhì)的5.44%硝酸鉀水溶液冰點(diǎn)模型計(jì)算值與實(shí)驗(yàn)值對(duì)比曲線
將硝酸鉀水溶液的冰點(diǎn)溫度分別與雜質(zhì)離子的濃度數(shù)據(jù)經(jīng)二次多項(xiàng)式擬合得到通式:
式中,參數(shù)a、b和c見表3。
表3 含不同離子雜質(zhì)的5.44%硝酸鉀溶液冰點(diǎn)方程參數(shù)表
2.3 含不同雜質(zhì)離子的硝酸鉀水溶液結(jié)晶點(diǎn)的測(cè)定結(jié)果與分析
為考察不同雜質(zhì)離子對(duì)硝酸鉀水溶液結(jié)晶點(diǎn)溫度的影響,先配制質(zhì)量分?jǐn)?shù)為5.44%的硝酸鉀溶液,向溶液中添加不同濃度的雜質(zhì)離子(Na+、Mg2+、Al3+、Cl-、SO42-),設(shè)定攪拌速率為400 r/min,待溶液充分混合后,以0.05℃/min的降溫速率將溶液從0℃開始冷凍,在冷凍結(jié)晶過程中記錄冰點(diǎn)溫度,結(jié)果見圖5。
將硝酸鉀水溶液的結(jié)晶點(diǎn)溫度分別與雜質(zhì)離子的濃度數(shù)據(jù)經(jīng)二次多項(xiàng)式擬合得到通式:
圖5 不同離子雜質(zhì)對(duì)5.44%硝酸鉀水溶液結(jié)晶點(diǎn)溫度的影響
式中,參數(shù)a、b和c見表4。
表4 含不同離子雜質(zhì)的5.44%硝酸鉀溶液結(jié)晶點(diǎn)方程參數(shù)表
從圖5和表4可以看出,同種離子濃度越高,溶液的結(jié)晶點(diǎn)溫度越低;陽離子雜質(zhì)價(jià)態(tài)越高,溶液的結(jié)晶點(diǎn)溫度越低,陰離子雜質(zhì)價(jià)態(tài)越低,溶液的結(jié)晶點(diǎn)溫度越低。
2.4 硝酸鉀水溶液的介穩(wěn)區(qū)測(cè)定結(jié)果與分析
2.4.1 降溫速率對(duì)硝酸鉀水溶液介穩(wěn)區(qū)寬度的影響
為考察不同降溫速率對(duì)硝酸鉀水溶液介穩(wěn)區(qū)寬度的影響,實(shí)驗(yàn)在測(cè)量介穩(wěn)區(qū)時(shí)保持一定的濃度、降溫速率(0.05℃/min)等條件,實(shí)驗(yàn)結(jié)果見圖6。
根據(jù)經(jīng)典成核理論和J.Nyvlt[10]的實(shí)驗(yàn)結(jié)果,當(dāng)攪拌速率一定時(shí):
式中:b為降溫速率,℃/h;m為成核級(jí)數(shù);k為成核速率常數(shù);ΔTmax為介穩(wěn)區(qū)寬度,℃。
將ln b與ln(ΔTmax)作圖可得一條直線,直線的斜率為成核級(jí)數(shù)m的倒數(shù)。采用式(7)對(duì)不同降溫速率下硝酸鉀水溶液的介穩(wěn)區(qū)寬度進(jìn)行關(guān)聯(lián),得到不同冰點(diǎn)下的模型方程,結(jié)果見表5。
表5 不同冰點(diǎn)的介穩(wěn)區(qū)寬度模型參數(shù)
由表5可見,成核級(jí)數(shù)m值變化不大,說明硝酸鉀水溶液結(jié)晶與冰點(diǎn)溫度基本無關(guān),取平均值m= 1.46。選擇上述3組模型方程的計(jì)算數(shù)據(jù)與實(shí)驗(yàn)值作對(duì)比,結(jié)果見圖6。從圖6可見,模型方程的計(jì)算值與實(shí)驗(yàn)值十分接近,符合經(jīng)典成核理論,即隨著降溫速率的增大,介穩(wěn)區(qū)寬度增大。
圖6 不同降溫速率下介穩(wěn)區(qū)寬度模型方程計(jì)算值與實(shí)驗(yàn)值對(duì)比曲線
2.4.2 不同雜質(zhì)離子對(duì)硝酸鉀水溶液介穩(wěn)區(qū)寬度的影響
介穩(wěn)區(qū)寬度與攪拌速率、降溫速率、溫度及測(cè)定方法等因素息息相關(guān)。實(shí)驗(yàn)在測(cè)量介穩(wěn)區(qū)時(shí),選擇質(zhì)量分?jǐn)?shù)為5.44%、攪拌速率為400 r/min、降溫速率為0.05℃/min。在此條件下,考察了含有不同雜質(zhì)離子的硝酸鉀溶液的介穩(wěn)區(qū)寬度(ΔT),為工藝實(shí)驗(yàn)提供了熱力學(xué)數(shù)據(jù),實(shí)驗(yàn)結(jié)果見圖7。
圖7 不同離子雜質(zhì)對(duì)5.44%硝酸鉀溶液介穩(wěn)區(qū)寬度的影響
從圖7可以看出,若硝酸鉀溶液中摻雜雜質(zhì)則溶液的介穩(wěn)區(qū)寬度會(huì)增大,同時(shí)從圖7中可以看出離子雜質(zhì)濃度越高,溶液的介穩(wěn)寬度越寬;陽離子價(jià)態(tài)越高,溶液介穩(wěn)區(qū)寬度越寬,陰離子雜質(zhì)價(jià)態(tài)越低,溶液介穩(wěn)區(qū)寬度越寬。
通過對(duì)硝酸鉀溶液冷凍結(jié)晶過程進(jìn)行較為系統(tǒng)的研究,得到結(jié)論:1)將通過實(shí)驗(yàn)建立的測(cè)定冰點(diǎn)的方法測(cè)得的冰點(diǎn)值,分別與文獻(xiàn)值和 Extend-UNIQUAC模型計(jì)算值作對(duì)比,結(jié)果表明測(cè)定冰點(diǎn)方法是可行的。2)溶液中同種離子雜質(zhì)的濃度越高,其冰點(diǎn)溫度越低,不同離子雜質(zhì)條件下,離子半徑越小,冰點(diǎn)溫度越低。3)溶液中陽離子雜質(zhì)的價(jià)態(tài)越高,其結(jié)冰點(diǎn)溫度越低,介穩(wěn)區(qū)寬度越寬;反之陰離子雜質(zhì)價(jià)態(tài)越低,冰點(diǎn)和結(jié)冰點(diǎn)溫度越低,介穩(wěn)區(qū)寬度越寬。4)用經(jīng)典成核理論推導(dǎo)了介穩(wěn)區(qū)寬度的模型方程,結(jié)果與實(shí)驗(yàn)值符合較好,成核級(jí)數(shù)m為1.46。
[1] 李亞,孫瀟,孫衛(wèi)東.冷凍濃縮技術(shù)的應(yīng)用及研究進(jìn)展[J].食品與生物,2008(3):9-10,16.
[2] 謝紹雷,張超,紀(jì)律,等.硫酸鎂亞型鹽田鹵水低溫冷凍結(jié)晶規(guī)律的研究[J]無機(jī)鹽工業(yè),2015,47(2):35-38.
[3] Reddy ST,LewisAE,Witkamp GJ,etal.RecoveryofNa2SO4·10H2O from a reverse osmosis retentate by eutectic freeze crystallisation technology[J].Chemical Engineering Research and Design,2010,88(9):1153-1157.
[4] Lee LS,Sun SL,Lin C L.Predictionsof thermodynamic properties ofaqueous single-electrolyte solutionswith the two-ionic-parameter activity coefficientmodel[J].Fluid Phase Equilibria,2008,264(1/ 2):45-54.
[5] AliHaghtalab,Kiana Peyvandi.Generalized electrolyte-UNIQUACNRFmodel for calculation ofsolubility and vapor pressureofmulticomponent electrolytes solutions[J].Journal of Molecular Liquids,2012,165:101-112.
[6] AliHaghtalab,Kiana Peyvandi.Electrolyte-UNIQUAC-NRFmodel for the correlation of the mean activity coefficient of electrolyte solutions[J].Fluid Phase Equilibria,2009,281(2):163-171.
[7] Zhao E S,Yu M,SauvéR E,et al.Extension of theWilsonmodel to electrolyte solutions[J].Fluid Phase Equilibria,2000,173(2):161-175.
[8] Garcia A V,Thomsen K,Stenby E H.Prediction ofmineral scale formation in geothermaland oilfield operationsusing the Extended UNIQUACmodel.PartⅡ.Carbonate-scalingminerals[J].Geothermics,2006,35(3):239-284.
[9] 鄧天龍,周桓,陳俠.水鹽體系相圖原理及其應(yīng)用[M].北京:化學(xué)工業(yè)出版社,2013.
[10] Nyvlt J.Kinetics of nucleation in solutions[J].Journal of Crystal Growth,1968,3/4:377-383.
Effectof impurity ionson crystallization propertiesof ice crystals in eutectic freeze crystallization process
Hu Chaoqun1,Zhu Liang1,2,Sha Zuoliang1,Wang Yanfei1,Yang Libin1,Zhao Nannan1
(1.College ofChemical Engineering and MaterialsScience,Tianjin University of Science&Technology,Tianjin 300457,China;2.Myande Group Limited)
The effectof impurity ions in potassium nitrate solution on the crystallization properties of ice crystals during the eutectic freeze crystallization was investigated.At first,the determination method of freezing point was developed and verified,and the freezing point predicted modelwas developed based on Extend-UNIQUANmodel.With the determination method of freezing point,the influence of impurities ionson the freezing point,crystallization point,andmetastablewidth was investigated.Results showed that the higher concentrations of impurities in potassium nitrate aqueous solution,the lower temperatureof freezingpointand crystallization pointand thewiderofmetastable zone.However,given the same concentrations of different impurity ions in potassium nitrate aqueous solutions,the temperature of freezing pointwas inversely proportional to the radius of ions.In addition,the temperature of crystallization point decreased with increase in valence of cations and decreased in valence of anions,aswell as the width ofmetastable zone.The simulationmodel ofmetastable zone width was derived by means of the classical nucleation theory.The calculated results were in good agreement with the experimental data,and the nucleation orderof ice crystal in potassium nitrate aqueous solutionswasestimated to be 1.46.
ice crystal;metastable zone;potassium nitrate;eutectic freeze crystallization
TQ131.13
A
1006-4990(2017)05-0025-05
2016-11-22
胡超群(1993— ),女,碩士研究生,主要研究方向?yàn)楣I(yè)結(jié)晶。
朱亮
天津市應(yīng)用基礎(chǔ)與前沿技術(shù)研究計(jì)劃(重點(diǎn)項(xiàng)目)(14JCZDJC40900)。
聯(lián)系方式:zhuliang@tust.edu.cn