沿面放電等離子體滅活蘋(píng)果汁中耐高滲酵母的模型

孫嘉莉，王 英，敖 羽，張文樂(lè)，張紹君，包曉瑋

（新疆農(nóng)業(yè)大學(xué)食品科學(xué)與藥學(xué)學(xué)院，新疆 烏魯木齊 830052）

在食品工業(yè)中，確保食品安全和延長(zhǎng)其保質(zhì)期是至關(guān)重要的。盡管傳統(tǒng)的熱處理方法在滅活食品中的微生物方面效果顯著，但這些方法往往會(huì)損害食品的營(yíng)養(yǎng)價(jià)值和感官品質(zhì)。因此，開(kāi)發(fā)既能有效滅活微生物又能保持食品品質(zhì)的新技術(shù)，已成為食品科學(xué)領(lǐng)域的一個(gè)重要研究方向。蘋(píng)果果汁因其營(yíng)養(yǎng)價(jià)值及獨(dú)特的口味，在世界范圍內(nèi)的消費(fèi)日益增長(zhǎng)[1]。然而酵母很容易使蘋(píng)果汁發(fā)生腐敗，如果不使用適當(dāng)?shù)姆椒▽?duì)其進(jìn)行滅活，就會(huì)對(duì)人體健康造成嚴(yán)重危害[2]。其中，耐高滲酵母（Zygosaccharomyces rouxii）是食品工業(yè)中著名的腐敗酵母，也是造成蘋(píng)果汁腐敗的原因之一[3]。Z.rouxii因?yàn)榫哂心退?、耐高滲透壓、耐弱酸抑菌劑以及對(duì)高溫環(huán)境的適應(yīng)性等特點(diǎn)，可導(dǎo)致高糖類食品發(fā)生變質(zhì)，從而使蘋(píng)果汁發(fā)酵變質(zhì)而腐敗，使產(chǎn)品品質(zhì)下降[4-6]。因此，探尋一種對(duì)蘋(píng)果汁中Z.rouxii安全、高效的滅活方法勢(shì)在必行。

等離子體被視為物質(zhì)的第四種存在狀態(tài)[7]，是含有多種活性物質(zhì)的電離氣體，由于原子或者原子團(tuán)失去電子后電離形成[8]。根據(jù)激發(fā)出的電子溫度，等離子體可以劃分為高溫等離子體和低溫等離子體。在低溫條件下，氣體中的原子或分子被激發(fā)而形成的帶電粒子云為放電等離子體。而沿面放電等離子體是指在介質(zhì)表面上形成的放電等離子體。當(dāng)在介質(zhì)表面施加足夠高的電場(chǎng)或電壓時(shí)，氣體分子會(huì)被激發(fā)并形成帶電粒子云，從而形成沿面放電等離子體。大氣壓低溫等離子體技術(shù)是對(duì)食品工業(yè)發(fā)展有促進(jìn)作用的新型非熱滅菌技術(shù)，因此廣受全球研究者的關(guān)注[9]。研究表明，由于滲透深度小，等離子體工藝對(duì)食物處理存在限制性。但這些研究同時(shí)表明，低溫等離子體具有對(duì)食品表面微生物滅活的潛力[10-11]。低溫等離子體技術(shù)作為非熱處理技術(shù)之一[9]，在滅菌的過(guò)程中有很多的優(yōu)點(diǎn)，如安全、高效、不易產(chǎn)生殘留[12]。放電等離子體處理速度快，處理能力強(qiáng)，已經(jīng)被證實(shí)可成功滅活火腿表面的單核細(xì)胞增生李斯特菌[13]、牛肉干表面的金黃色葡萄球菌[14]、鹽水鴨表面的大腸桿菌[15]等，因此，應(yīng)用沿面放電等離子體技術(shù)滅活蘋(píng)果果汁中的Z.rouxii是可行的。

由于受到實(shí)驗(yàn)條件的制約，確定適宜的等離子體處理工藝既費(fèi)時(shí)又費(fèi)力，并且在技術(shù)上實(shí)施起來(lái)有一定的難度。計(jì)算機(jī)模擬技術(shù)因其不需要花費(fèi)大量實(shí)驗(yàn)時(shí)間，能夠從不同數(shù)學(xué)建模中得到準(zhǔn)確結(jié)果，可用作等離子體傳質(zhì)過(guò)程快速分析的一種手段，方便等離子體系統(tǒng)設(shè)計(jì)和擴(kuò)充[16-17]。因此，利用計(jì)算機(jī)模擬技術(shù)對(duì)等離子體處理過(guò)程進(jìn)行研究是目前等離子體技術(shù)研究的一個(gè)新動(dòng)向。

因此，本研究采用數(shù)值模擬軟件COMSOL Multiphysics建立模擬處理過(guò)程中微生物滅活作用的數(shù)值模型，并在不同放電電壓下驗(yàn)證模型，分析不同氣泡大小和不同處理量對(duì)失活效果的影響，以探討研究沿面放電等離子滅活液體中微生物的效果。然后，將模型模擬所得失活曲線與實(shí)驗(yàn)所得失活曲線進(jìn)行比較，進(jìn)一步研究沿面放電等離子體對(duì)液體微生物的滅活影響。最后，通過(guò)數(shù)據(jù)分析確認(rèn)沿面放電等離子體在蘋(píng)果汁工業(yè)中滅活Z.rouxii的方法是否可行。該模型有利于沿面放電等離子體滅活微生物的工藝開(kāi)發(fā)，同時(shí)為沿面放電等離子體在食品工業(yè)上的應(yīng)用建立理論基礎(chǔ)，促進(jìn)等離子體技術(shù)在食品工業(yè)中的發(fā)展。

1 材料與方法

1.1 材料與試劑

紅富士蘋(píng)果，采購(gòu)于陜西省楊凌市，采收時(shí)間為9月中下旬，果實(shí)大小均一，表面無(wú)機(jī)械損傷，無(wú)病害或腐爛等問(wèn)題，將蘋(píng)果去皮榨汁后冷藏備用。

本研究采用的Z.rouxiiLB（B-WHX-12-54）由西北農(nóng)林科技大學(xué)提供，GenBank登錄號(hào)為KC544459。將Z.rouxiiLB菌株培養(yǎng)在（28.0±0.1）℃條件下的酵母浸出粉胨葡萄糖（yeast peptone dextrose agar，YPD）液體培養(yǎng)基（20 g/L蛋白胨、20 g/L葡萄糖和10 g/L酵母抽取物）中，振蕩培養(yǎng)48 h（150 r/min）。后將培養(yǎng)液在9 000 r/min條件下離心8 min，去除上清液，用無(wú)菌蒸餾水清洗2 遍。將所得酵母重懸在無(wú)菌蒸餾水中，可獲得穩(wěn)定生長(zhǎng)階段的微生物，約為107～108CFU/mL，后冷藏保存?zhèn)溆?。將Z.rouxiiLB接種到每個(gè)蘋(píng)果汁樣品中，并使每一個(gè)樣品中Z.rouxiiLB的菌落數(shù)量約為1.5×106CFU/mL。

葡萄糖 國(guó)藥集團(tuán)化學(xué)試劑有限公司；氯化鈉（分析純）四川西隴化工有限公司；瓊脂粉 上海源葉生物科技有限公司；蛋白胨 上海中秦化學(xué)試劑有限公司。

1.2 儀器與設(shè)備

JX-RS485-20PL型臭氧傳感器 北京金訊昌通電子科技有限公司；i-SPEED數(shù)碼相機(jī) 奧林巴斯（中國(guó)）有限公司。

沿面放電等離子體反應(yīng)器系統(tǒng)示意圖見(jiàn)圖1。沿面放電等離子體系統(tǒng)由高壓交流電源（0～30 kV、50 Hz）、接地電極和有機(jī)玻璃容器組成。反應(yīng)器容器直徑40 mm、高350 mm，采用有機(jī)玻璃圓筒。高壓電極為直徑1 mm、放電長(zhǎng)度100 mm的不銹鋼彈簧，并附著在直徑為12 mm、高度為350 mm的圓柱形石英管內(nèi)壁上。反應(yīng)器中的細(xì)胞懸浮液為接地電極，沿著石英管的縱軸方向沿面放電[18]。當(dāng)沿面放電等離子體系統(tǒng)啟動(dòng)時(shí)，干燥空氣由氣泵泵入圓柱形石英管，并在石英管底部通過(guò)曝氣器以150 L/h的流速進(jìn)入反應(yīng)器中進(jìn)行處理，總處理量為300 mL。

圖1 沿面放電等離子體反應(yīng)器示意圖Fig.1 Schematic diagram of surface discharge plasma reactor

1.3 方法

1.3.1 物理模型的建立

利用數(shù)值仿真軟件COMSOL Multiphysics建立沿面放電等離子體模型，如圖2所示，研究在21、15、11 kV放電電壓，直徑為5.46 mm的大氣泡和直徑為3.21 mm的小氣泡，反應(yīng)器直徑為40、50 mm和80 mm時(shí)的處理量條件下沿面放電等離子體滅活蘋(píng)果果汁中Z.rouxii的效果，并對(duì)此模型進(jìn)行驗(yàn)證實(shí)驗(yàn)。

圖2 模擬中使用的三維幾何體Fig.2 Three-dimensional geometry used in the simulation

沿面放電等離子體是可以產(chǎn)生多種反應(yīng)物質(zhì)的化學(xué)反應(yīng)，如H2O2、O3、·OH和·O[19]。因?yàn)槌粞跏歉鞣N等離子體中公認(rèn)的成分，并且具有成熟的測(cè)量方法。因此，本項(xiàng)研究中沿面放電等離子的建模參數(shù)采用臭氧濃度。等離子體放電時(shí)，出口處采用精準(zhǔn)度為2 mg/m3的臭氧傳感器每隔30 s測(cè)定一次臭氧濃度。在放電過(guò)程中，·O是在氧分子被激發(fā)并通過(guò)高能電子攻擊解離后產(chǎn)生的，而臭氧則由·O與O2反應(yīng)生成，如式（1）所示[19]，臭氧的濃度與放電電壓有關(guān)。

1.3.2 控制方程

1.3.2.1 湍流中的氣泡運(yùn)輸

在酵母經(jīng)臭氧滅活時(shí)，假定氣泡一直以終端速度連續(xù)移動(dòng)?？勺魅缦录僭O(shè)：

氣體密度于液體密度而言，可忽略；通過(guò)黏性力和壓力平衡，確定氣泡在液體中的相對(duì)運(yùn)動(dòng)；這兩個(gè)因素具有相同的壓力場(chǎng)。

因此，液體RANSk-ε湍流方程是可解的，同時(shí)臭氧氣泡的速度可由滑移模型進(jìn)行引導(dǎo)。臭氧氣泡的體積分?jǐn)?shù)可由混合平均連續(xù)性方程計(jì)算出壓力分布并由有效氣體密度的傳輸方程解出，如式（2）所示[20]，其中下標(biāo)l表示與液體相關(guān)的量，下標(biāo)g表示和氣體有關(guān)的量。

連續(xù)性方程為式（3），臭氧在水中的氣體輸送過(guò)程為式（4）：

式中：ρg為臭氧的氣體密度/（kg/m3）；為臭氧體積分?jǐn)?shù)/（m3/m3）；ug為氣體速度矢量/（m/s）；mgl為從氣體到液體的傳質(zhì)速率/（kg/（m3·s））；下同。

1.3.2.2 稀釋物種運(yùn)輸

參考馬國(guó)瑜等[21]對(duì)于擴(kuò)散問(wèn)題的處理方法，通過(guò)使用稀釋物種運(yùn)輸模型，能夠獲取溶劑中臭氧和酵母的濃度分布。這種模型能夠處理溶劑中物質(zhì)的運(yùn)輸和反應(yīng)過(guò)程。其中，物質(zhì)在溶劑中的運(yùn)輸和反應(yīng)主要在這一界面上進(jìn)行。菲克定律描述的擴(kuò)散現(xiàn)象，以及由流場(chǎng)引起的對(duì)流現(xiàn)象，共同為微生物提供了運(yùn)輸?shù)膭?dòng)力，上述內(nèi)容由式（5）、（6）進(jìn)行擬合：

式中：Di為物種擴(kuò)散系數(shù)/（m2/s）；ci為物質(zhì)濃度/（mol/m3）；u為流體速度矢量/（m/s）；Ri為物種反應(yīng)速率/（mol/（m3·s））；Ni為通量矢量/（mol/（m2·s））；i指代第i類物質(zhì)。

在運(yùn)輸?shù)倪^(guò)程當(dāng)中，臭氧與酵母菌兩者互相反應(yīng)且反應(yīng)速率相等，定義為式（7）：

式中：RO3為臭氧的反應(yīng)速率/（mol/（m3·s））；Ry為酵母菌液的反應(yīng)速率/（mol/（m3·s））；cy為酵母菌液的濃度/（mol/m3）；cO3為臭氧的濃度/（mol/m3）；kreac為反應(yīng)速率系數(shù)/（m3/（mol·s））。kreac由微生物的常數(shù)所確定，不能被預(yù)測(cè)。參考Zhou Pei等[22]的方法，其數(shù)值可以通過(guò)實(shí)驗(yàn)法得出，用m3/（mol·s）表示。對(duì)于已經(jīng)確定的氣體流速與等離子設(shè)備，kreac與N0（初始濃度）存在相關(guān)性。

1.3.3 初始和邊界條件

曝氣器位于反應(yīng)器底層，氣體流速可以通過(guò)氣體流量計(jì)得到穩(wěn)定的控制。通過(guò)k-ε湍流方程（2）中的重力項(xiàng)考慮重力：將壓力點(diǎn)約束（p＝0）添加到出口邊界（靠近大氣）上的點(diǎn)，規(guī)定壓力等級(jí)，提供一個(gè)界值，用于COMSOL啟動(dòng)模擬。

氣體入口處數(shù)學(xué)描述如式（8）～（11）所示：

式中：k為湍流動(dòng)能/（（mol·s2）/m2）；n為入口邊界處的法向量；ul,tang為壁面切向速度/（m/s）。

幾何體的頂端是可以忽略表面運(yùn)動(dòng)的自由表面，以液體的滑動(dòng)條件對(duì)表面近似模擬。臭氧可以經(jīng)過(guò)此邊界條件自由離開(kāi)反應(yīng)器。

在模擬流速恒定且進(jìn)口物質(zhì)通量率穩(wěn)定的臭氧源時(shí)，由于酵母菌僅在水中存在，因此進(jìn)入口濃度定為零。

1.3.4 模型參數(shù)

等離子體滅活酵母的過(guò)程在室溫條件下進(jìn)行，在此過(guò)程的建模中，臭氧的物理性質(zhì)具有重要意義[23-25]。

臭氧的物理性質(zhì)被確定為：

臭氧密度：ρg＝2.14 kg/m3；臭氧擴(kuò)散系數(shù)：Dl＝1.74×10-9m2/s；溫度：T＝290.15 K；氣體速度矢量：Vg＝8.33×10-6m3/s。

氣泡大小對(duì)等離子體滅活微生物過(guò)程中的數(shù)值模擬有重要約束作用。為了研究氣泡大小對(duì)滅活效果的影響，采用兩種不同大小的曝氣器以產(chǎn)生不同尺寸的氣泡，如圖3所示，氣泡大小采用照相法進(jìn)行測(cè)量。為了避免受蘋(píng)果汁色澤的影響，拍攝過(guò)程在純凈水中進(jìn)行。利用數(shù)碼相機(jī)（幀速率：2 000 fps），在反應(yīng)器的中間位置獲取氣泡圖像。然后使用圖像處理系統(tǒng)處理微氣泡，選取大約200 個(gè)氣泡同時(shí)統(tǒng)計(jì)其直徑大小，計(jì)算出Sauter平均直徑[25]，如式（12）所示：

圖3 氣泡照片F(xiàn)ig.3 Photographs of bubbles

式中：dbs為氣泡Sauter平均直徑/m；di為氣泡直徑/m；ni為di的氣泡直徑范圍內(nèi)的氣泡數(shù)量/個(gè)；根據(jù)計(jì)算，大氣泡的直徑平均為5.46 mm，小氣泡的直徑平均為3.21 mm。

1.3.5 微生物計(jì)數(shù)

等離子處理后的微生物通過(guò)細(xì)胞培養(yǎng)法進(jìn)行確定，并使用YPD培養(yǎng)基在28 ℃下培養(yǎng)72 h計(jì)算菌落數(shù)，一式三份進(jìn)行分析。使用一級(jí)動(dòng)力學(xué)模型[26]，如式（13）所示，得到失活曲線。

式中：N0為微生物初始濃度/（CFU/mL）；N為經(jīng)等離子體處理t時(shí)間后微生物的濃度/（CFU/mL）；t為等離子體處理時(shí)間/min；D為殺滅其中90%微生物所需要的時(shí)間/min。

1.3.6 模型驗(yàn)證

實(shí)驗(yàn)利用沿面放電等離子體反應(yīng)系統(tǒng)驗(yàn)證計(jì)算機(jī)模型，放電電壓（21、15、11 kV）。實(shí)驗(yàn)中，在固定時(shí)間取反應(yīng)器中含有酵母的蘋(píng)果汁樣品，并采用細(xì)胞培養(yǎng)法測(cè)得酵母濃度。所有實(shí)驗(yàn)進(jìn)行3 次，利用COMSOL軟件對(duì)酵母菌的存活曲線進(jìn)行分析。

1.3.7 模型應(yīng)用

1.3.7.1 不同氣泡大小對(duì)滅活效果的影響

采用兩種不同大小的曝氣器以產(chǎn)生不同尺寸的氣泡，以探究不同氣泡大小對(duì)滅活效果的影響。根據(jù)計(jì)算，選擇直徑為5.46 mm的大氣泡和直徑為3.21 mm的小氣泡。將初始菌液濃度定為1.5×102、1.5×104CFU/mL和1.5×106CFU/mL，在模型中輸入氣泡大小參數(shù)得到存活曲線，并將其與實(shí)驗(yàn)結(jié)果作對(duì)比。

1.3.7.2 不同處理量對(duì)滅活效果的影響

不同處理量對(duì)氣液傳質(zhì)過(guò)程的影響在反應(yīng)器的設(shè)計(jì)和放大的過(guò)程中起著重要作用[20]。本實(shí)驗(yàn)在擴(kuò)大反應(yīng)容器的半徑時(shí)，加工的蘋(píng)果汁的量也相應(yīng)地改變。當(dāng)反應(yīng)容器的直徑為40、50 mm和80 mm時(shí)，相應(yīng)的加工蘋(píng)果汁的體積分別擴(kuò)大了1、1.56 倍和4 倍。將參數(shù)輸入到模型中，并將模擬值與實(shí)驗(yàn)值進(jìn)行對(duì)比，進(jìn)而定量分析處理量對(duì)滅活效果的影響。

1.4 數(shù)據(jù)處理與統(tǒng)計(jì)分析

所有實(shí)驗(yàn)重復(fù)進(jìn)行3 次。圖表繪制使用軟件Origin 2021，數(shù)據(jù)使用SPSS 26進(jìn)行方差分析（ANOVA），數(shù)據(jù)平均值間的顯著性分析采用Tukey多重比較法，水平為P＝0.05，P＜0.05時(shí)效果影響顯著，P≥0.05效果影響不顯著。

2 結(jié)果與分析

2.1 模型參數(shù)

圖4體現(xiàn)了在具有差異的初始菌液濃度下進(jìn)行滅活實(shí)驗(yàn)，所獲得相應(yīng)的反應(yīng)速率系數(shù)。實(shí)驗(yàn)中放電電壓為21 kV，氣泡大小為3.21 mm，處理量為300 mL。結(jié)果表明，隨著初始菌液濃度的升高，kreac值降低，當(dāng)初始菌液濃度為1.5×102CFU/mL時(shí)，獲得最大kreac值，為0.256，而當(dāng)初始菌液濃度為1.5×106CFU/mL時(shí)，獲得最小kreac值，為0.202。所得結(jié)果可能是因?yàn)槲⑸锏亩逊e結(jié)構(gòu)可以在等離子處理過(guò)程中作為保護(hù)屏障[27]。實(shí)驗(yàn)數(shù)據(jù)顯示，kreac值與初始菌液濃度值N0值之間的線性關(guān)系可以用kreac＝-0.014 38lgN0+0.288 05（R2＝0.98）描述，這些數(shù)據(jù)進(jìn)一步被用于計(jì)算機(jī)模擬模型中，并應(yīng)用COMSOL軟件中的線性插值函數(shù)進(jìn)行輸入。

圖4 反應(yīng)速率系數(shù)（kreac）與酵母濃度（N0）之間的相關(guān)性Fig.4 Correlation between reaction rate coefficient (kreac) and yeast concentration (N0)

2.2 模型驗(yàn)證

當(dāng)放電電壓設(shè)置為21、15、11 kV時(shí)，通過(guò)臭氧傳感器測(cè)量得到入口處臭氧濃度分別為21.43、17.14 mg/L和12.86 mg/L。放電電壓的變化會(huì)影響活性物質(zhì)的濃度，即較高的區(qū)域濃度值與能量密度的增加相關(guān)[28]。Mai-Prochnow等[29]在白色念珠菌和小孢子菌的滅活過(guò)程中也發(fā)現(xiàn)了類似的結(jié)果。通過(guò)實(shí)驗(yàn)結(jié)果可得，放電電壓的升高可在達(dá)到同樣的微生物滅活效果時(shí)，等離子體的處理時(shí)間縮短。

圖5顯示了當(dāng)初始菌液濃度為1.5×104CFU/mL、氣泡大小為3.21 mm、處理量為300 mL時(shí)，在3 種不同的放電電壓21、15 kV和11 kV下，得到的實(shí)驗(yàn)和模擬的酵母失活曲線。將3 條失活曲線作比較，模擬結(jié)果與實(shí)驗(yàn)結(jié)果具有很好的重合性。酵母被滅活所需要的時(shí)間隨著放電電壓的升高及臭氧濃度的增加而減少。在細(xì)胞計(jì)數(shù)采樣過(guò)程中，存儲(chǔ)時(shí)間會(huì)使模擬結(jié)果和實(shí)驗(yàn)結(jié)果產(chǎn)生差異。即在液體中產(chǎn)生的副產(chǎn)物經(jīng)過(guò)等離子體處理后仍可以繼續(xù)發(fā)揮殺菌作用，這會(huì)導(dǎo)致模擬結(jié)果較實(shí)驗(yàn)結(jié)果稍高。在等離子體處理蘋(píng)果汁中的弗氏檸檬酸桿菌[26]的滅活過(guò)程中也觀察到了類似的結(jié)果。表1將3 種放電電壓下等離子體對(duì)Z.rouxiiLB滅活過(guò)程的模擬D值和實(shí)驗(yàn)D值進(jìn)行對(duì)比。當(dāng)放電電壓為21、15 kV和11 kV時(shí)，D值誤差分別為4.652 5%、4.615 8%和3.306 3%。

表1 3 種不同放電電壓下等離子體對(duì)Z.rouxii LB滅活過(guò)程的實(shí)驗(yàn)值和模擬值Table 1 Experimental and simulated values for the inactivation process of Z.rouxii LB by plasma under three different discharge voltages

圖5 不同放電電壓對(duì)等離子體滅活Z.rouxii LB的影響Fig.5 Effect of different discharge voltages on SDP inactivation of Z.rouxii LB

2.3 不同氣泡大小的影響

圖6顯示了在不同氣泡大小3.21 mm和5.46 mm下沿面放電等離子體處理蘋(píng)果汁中Z.rouxiiLB的滅活曲線。沿面放電等離子體放電電壓為21 kV，Z.rouxiiLB初始濃度為1.5×102、1.5×104、1.5×106CFU/mL，處理量為300 mL。

圖6 不同氣泡大小對(duì)等離子體滅活Z.rouxii LB的影響Fig.6 Effect of different bubble sizes on SDP inactivation of Z.rouxii LB

結(jié)果表明，在21 kV的放電電壓下，初始菌液濃度為1.5×102CFU/mL時(shí)，3.21 mm和5.46 mm氣泡大小模擬的D值分別為2.493 2 min和2.612 8 min，兩者之間相差4.575 1%；初始菌液濃度為1.5×104CFU/mL時(shí)，模擬D值分別為2.560 6 min和2.6712 min，兩者之間相差4.1455%；初始菌液濃度為1.5×106CFU/mL時(shí)，模擬D值分別為3.644 7 min和3.936 1 min，兩者之間相差7.402%，在處理過(guò)程中，測(cè)得的實(shí)驗(yàn)值比模擬值偏低，尤其是在初始菌液濃度較高，為1.5×106CFU/mL時(shí)。這可能是由于脫離處理后，等離子體的殘留活性物質(zhì)依然可以發(fā)生滅活作用。此現(xiàn)象與Magureanu[30]和郭儉[31]等的研究一致，其中的機(jī)理有待進(jìn)一步研究。在不同氣泡大小情況下D值誤差均不超過(guò)8%，實(shí)驗(yàn)結(jié)果與模擬結(jié)果相吻合。由此可見(jiàn)，氣泡大小的改變對(duì)等離子體滅活作用不會(huì)產(chǎn)生顯著的影響（P≥0.5）。

2.4 不同處理量的影響

圖7顯示了在反應(yīng)器直徑為40、50 mm和80 mm時(shí)，沿面放電等離子體處理蘋(píng)果汁中Z.rouxiiLB的模型及實(shí)驗(yàn)與模擬D值的比較。此時(shí)，加工蘋(píng)果汁的體積分別擴(kuò)大了1、1.56 倍和4 倍。沿面放電等離子體放電電壓為21 kV，Z.rouxiiLB初始濃度為1.5×104CFU/mL，氣泡大小3.21 mm。結(jié)果表明，在21 kV的放電電壓下，初始菌液濃度為1.5×104CFU/mL時(shí)，當(dāng)反應(yīng)容器的直徑為40、50 mm和80 mm時(shí)，實(shí)驗(yàn)D值分別為2.735、4.002 min和4.759 min，模擬D值分別是2.561、3.855 min和4.388 min。實(shí)驗(yàn)D值與模擬值誤差不超過(guò)8%，匹配良好。研究證明，處理量的變化對(duì)等離子體滅活微生物效果影響顯著（P＜0.05）。當(dāng)蘋(píng)果汁的加工量增加時(shí)，達(dá)到所需滅活效果的加工時(shí)間更長(zhǎng)。

圖7 不同處理量對(duì)等離子體滅活Z.rouxii LB的影響Fig.7 Effect of different sample volumes on SDP inactivation of Z.rouxii LB

3 結(jié)論

在本研究中，利用數(shù)值仿真軟件COMSOL Multiphysics建立計(jì)算機(jī)模型，研究蘋(píng)果果汁中Z.rouxiiLB在沿面放電等離子體反應(yīng)器中的滅活過(guò)程。依據(jù)此項(xiàng)研究所建立的計(jì)算機(jī)模型，得到了不同放電電壓條件下酵母的存活曲線。通過(guò)驗(yàn)證后的模型，進(jìn)一步研究不同氣泡大小和處理量對(duì)等離子體滅活耐高滲酵母的影響，模擬結(jié)果和實(shí)驗(yàn)結(jié)果在各個(gè)條件下誤差不超過(guò)8%。實(shí)驗(yàn)結(jié)果與模擬結(jié)果同時(shí)表明，隨著初始菌液濃度的增加，反應(yīng)速率系數(shù)減小，隨著放電電壓的增加，滅活酵母所需時(shí)間也將縮短。不同氣泡尺寸對(duì)酵母菌的滅活作用不顯著，處理量增加會(huì)使滅活時(shí)間延長(zhǎng)。在模擬結(jié)果的基礎(chǔ)上，可以對(duì)等離子體滅活酵母的處理過(guò)程最佳化，同時(shí)為等離子體技術(shù)的工程化應(yīng)用奠定理論基礎(chǔ)。