基于CFD技術的瓶裝啤酒隧道式巴氏殺菌過程研究

程　菊安家彥董文勇王　越

(1.大連工業大學生物工程學院，遼寧 大連　116034；2. 大連市第21中學，遼寧 大連　116021)

?

基于CFD技術的瓶裝啤酒隧道式巴氏殺菌過程研究

程菊1安家彥1董文勇2王越1

(1.大連工業大學生物工程學院，遼寧 大連116034；2. 大連市第21中學，遼寧 大連116021)

通過數值模擬瓶裝啤酒隧道式巴氏殺菌的過熱階段(巴氏殺菌加熱部分的最后一個階段)，確定瓶中的冷核(slowest heating zone，SHZ)，并運用L9(34)正交試驗，分析噴淋水溫度、噴口處噴淋水湍流強度以及瓶子的運行速度三因素對瓶中啤酒溫度分布均勻性的影響。結果表明：冷核起先是位于瓶底，后向上移動，但不超過灌裝高度的1/2；在加熱至7 min時，噴淋水溫度、噴口處噴淋水湍流強度影響顯著，瓶子的運行速度則不顯著，試驗的可能的最優條件：噴淋水溫度為65 ℃，噴口處噴淋水湍流強度為3%，瓶子運行的速度為3 mm/s。

啤酒；巴氏殺菌；數值模擬；冷核

隧道式巴氏殺菌是啤酒工業化生產過程中，殺死啤酒中的微生物以提高其生物穩定性的最常用的方法之一[1]。隧道式巴氏殺菌機包括加熱區(包括預熱區和過熱區)、保溫區和冷卻區三大溫區[2-3]。在加熱區，啤酒溫度大約由2 ℃上升至60 ℃；在保溫區，維持在60 ℃；經冷卻區處理后，產品的溫度由60 ℃降至21～27 ℃[4]。瓶裝啤酒在巴氏殺菌過程中，啤酒以酒瓶為界面進行傳熱：靠近瓶壁的啤酒溫度升高，密度減小，而中心部分溫度相對較低，密度相對較大，產生自然對流[5-6]。這就形成了冷核。冷核是在食物產品熱處理過程中，吸收熱量最少(或殺菌強度最弱)的區域[7-8]。因此，在啤酒巴氏殺菌的過程中，既要保證冷核達到滅菌要求，又要使瓶中除冷核以外的區域不會過巴氏殺菌，以免改變啤酒的感官，造成啤酒品質的降低。這就要求瓶中啤酒溫度分布要均勻。

目前研究瓶中啤酒的溫度分布及冷核的方法有兩種：① 在線檢測瓶中溫度和數值模擬啤酒巴氏殺菌過程。在線監測過程中，溫度探針會影響瓶中液體的流動和溫度分布[9-10]，同時只能檢測一個位置，而冷核的溫度和位置會隨時間和包裝容器結構變化而改變[11-12]，且在線監測的費用較高，操作復雜。② 借助計算流體力學(computational fluid dynamics，CFD)技術，數值模擬啤酒巴氏殺菌過程，從而得到包裝容器內液體的傳熱規律，可有效控制食品的品質。而且CFD能縮短時間，減少花費[13-14]。目前，CFD模擬已被用于研究啤酒的巴氏殺菌過程，如：優化啤酒巴氏殺菌機的隧道結構[15]；啤酒瓶方位的變化對巴氏殺菌單位(pasteurization unit，PU)的影響[2]；瓶內對流流動對巴氏殺菌效果和老化效果的影響[3]；瓶裝和罐裝對包裝內啤酒溫度分布、速度分布和PU的影響[16]；啤酒巴氏殺菌結束后，啤酒的出口溫度和PU值是否達標[17]；隨著加熱的進行瓶內溫度分布和SHZ位置的變化以及啤酒中酵母的死滅程度(用PU表示)[8]等。

在巴氏殺菌過程中，瓶內啤酒的冷核和溫度分布受加熱時間、瓶子間距、瓶子移動的速度、噴淋水的溫度和噴口處噴淋水湍流強度的影響[18-19]。目前，已有關于瓶中啤酒冷核隨時間變化[8]及瓶中啤酒加熱速率隨噴淋水的噴淋強度變化的研究[19]，但未見瓶子移動的速度、噴淋水的溫度和噴口處噴淋水湍流強度對瓶中啤酒溫度分布的影響的研究。因此本研究擬考察隧道式巴氏殺菌的過熱階段，運用L9(34)正交試驗，分析噴淋水溫度、噴口處噴淋水湍流強度以及瓶子的運行速度三因素對瓶中啤酒溫度分布均勻性的影響，為瓶裝啤酒巴氏殺菌提供理論基礎。

1　材料與方法

1.1CFD模擬的驗證方法

選取高274 mm，直徑72 mm，容積為500 mL的啤酒瓶為研究對象，其模型見圖1。試驗所用的巴氏殺菌機型號為SBW40A(廣東輕工機械二廠有限公司)。壓好蓋的瓶裝生啤酒通過巴氏殺菌機的隧道時，用ZH5793 PU計(北京中慧天誠科技有限公司)分別測定隧道中不同位置上的瓶中啤酒的溫度(軸線上灌裝高度1/3處溫度)，計算過熱階段瓶中啤酒不同時間酒溫的平均值。

圖1　啤酒瓶模型

1.2幾何模型

因計算機的計算能力有限，動網格的計算復雜，在三維條件下，整個模型很難進行計算，而在二維條件下既能節約計算時間，又能節省人力和物力。所以整個模型是二維的。

巴氏殺菌機模型網格的建立使用前處理軟件Gambit 2.3.16，模型的網格數為209 319，流體計算使用Fluent 12.0。噴淋口與瓶口的距離為300 mm，啤酒瓶個數為142個(設置了邊界條件的為3個)，瓶子間距為16 mm。二維網格模型的部分示意圖見圖2。

1.3控制方程(governing equation)

運用Ansys Fluent 12(2009)軟件求解模型及相關邊界條件的連續性方程、動量方程、能量方程[20]。

圖2　二維網格模型的部分示意圖

(a) 連續性方程(continuity equation)或質量守恒方程(conversation of mass equation)：

(1)

式中：

ρ——密度，kg/m3；

t——時間，s；

x——軸向坐標；

u——軸向速度，m/s；

v——徑向速度，m/s；

r——徑向坐標；

Sm——源項，指加入到連續相的質量(例如液滴的蒸發)，也指其他自定義源項，kg。

(b) 動量守恒方程：

軸向的動量守恒方程：

(2)

徑向的動量守恒方程：

(3)

其中：

(4)

式中：

P——靜壓，Pa；

μ——黏度，kg/(m·s)；

Fx和Fr——外部體積力，N。

(c) 能量守恒方程：

(5)

其中：

(6)

(7)

式中：

τij——應力張量；

Sh——化學反應熱和其他體積熱，J；

keff——有效導熱系數，W/(m·℃)；

Jj’——組分j’擴散通量，mol/(m2·s)。

1.4模擬條件

啤酒、玻璃啤酒瓶、空氣和水的熱物理性質列于表1[8, 16, 21-23]。黏度模型采用k-epsilon中RNG模型，噴嘴及其附近的長方形區域為動網格，運動速度大小為5 mm/s。瓶內的初始溫度為48 ℃，瓶內壓力為1.55×105Pa。壓力和速度的耦合計算采用SIMPLEC算法。

1.5L9(34)正交試驗

瓶中啤酒溫度分布均勻性和冷核受噴淋水溫度、噴口處噴淋水湍流強度和瓶子運行速度的影響，因此確定四因素三水平的正交試驗，其試驗方案見表2。

表1　空氣、啤酒、玻璃和水的熱物理性質

表2　L9(34)正交試驗方案

2　結果與討論

2.1CFD模擬的驗證

為了確保CFD模擬的可靠性，先對其進行實驗驗證。啤酒在通過巴氏殺菌過熱區(工藝參數：噴淋水溫度為61 ℃，初始溫度48 ℃，加熱9 min)時，將啤酒瓶軸線上灌裝高度1/3位置處的溫度實驗值和模擬值進行對比，結果見圖3。由圖3可知，模擬值與實驗值的溫度曲線走勢相同，它們之間的最大標準差為0.21，所以可以認為它們是吻合的?？梢娔M所設置各項參數是合理的，可以進行接下來的數值模擬。

2.2冷核的確定

為了考察瓶內啤酒在過熱區的溫度分布情況，設過熱區噴淋水的溫度為61 ℃，噴口處噴淋水湍流強度為1%，瓶子運行速度為3 mm/s，按照方法1.2和1.3進行數值模擬，結果見圖4。

圖3　模擬溫度與實驗溫度曲線的比較

圖4　瓶中啤酒溫度分布云圖

由圖4可知，1，5，9 min時瓶中啤酒有明顯溫度分布梯度，具有冷核，且冷核位置是隨著時間在變化的。1 min時冷核位于瓶底，這與傳統的理論相符；但隨著時間的進行，冷核位置在不斷的變化，到5 min時，冷核的位置位于啤酒灌裝高度1/4～1/2。瓶內啤酒在加熱過程中由于局部溫度、密度等的差異，引起啤酒的自然對流，從而影響傳熱效果，溫度分布不均勻，因此有冷核存在?？拷勘诘钠【葡仁軣?，溫度較瓶中心的啤酒溫度更高，所以冷核一直位于啤酒瓶軸線附近的區域。而且隨著加熱的進行，瓶底也慢慢被加熱，靠近瓶底部的啤酒溫度也慢慢地上升，所以冷核從瓶底向瓶底以上部分移動。

2.3L9(34)正交試驗結果

為了考察瓶裝啤酒通過巴氏殺菌機的過熱區時，巴氏殺菌參數—噴淋水溫度、噴口處噴淋水湍流強度以及瓶子的運行速度，對瓶中啤酒溫度分布均一性的影響，利用L9(34)正交試驗，數值模擬啤酒巴氏殺菌的過熱階段，分析啤酒冷核面積占啤酒總面積的比例。由于隨著加熱時間的延長，冷核的溫度在不斷升高，根據9組試驗的不同時間的冷核溫度，將 1，3，5，7，9 min時冷核溫度的取值范圍分別設為小于：50.42，54.54，57.03，58.57，59.51 ℃。分別計算9組試驗在1，3，5，7，9 min時啤酒冷核面積占啤酒總面積的比例，結果列于表3。按照所設的冷核溫度范圍，對所得的試驗數據進行方差分析，發現1，3，5，9 min時噴淋水溫度、噴口處噴淋水湍流強度以及瓶子的運行速度對瓶中啤酒溫度分布沒有影響或影響不顯著。這是因為在加熱的1～5 min時，瓶中啤酒加熱時間短，溫度上升不顯著；加熱9 min時，瓶中啤酒溫度分布已經均勻。7 min時試驗數據的方差分析及極差分析結果分別見表4、5。

經過計算及方差分析，取檢驗水平α=0.05，由表4可知，噴淋水溫度和噴口處噴淋水湍流強度對瓶中啤酒溫度分布均一性的影響顯著，啤酒瓶運行速度則影響不顯著。由表5可知，三因素的影響主次順序為：B>A>C，即噴口處噴淋水湍流強度>噴淋水溫度>啤酒瓶運行速度，結果與方差分析的相符，而可能的最優條件組合為A3B3C1，即噴淋水溫度為65 ℃，噴口處噴淋水強度為3%，瓶子運行的速度為3 mm/s。

表3　　　　二維啤酒瓶模型中的啤酒冷核面積占啤酒總面積的比例

?*表示影響顯著。

表5　7 min時的極差分析

3　結論

本研究利用CFD技術模擬了500 mL瓶裝啤酒的隧道式巴氏殺菌過程的過熱階段，研究噴淋水溫度、噴淋口噴淋水湍流強度和瓶子運行速度對瓶中啤酒溫度分布的影響，結論如下：

(1) 運用商業軟件Ansys Fluent12模擬的巴氏殺菌過熱階段溫度曲線，得到了實驗溫度曲線的驗證，說明該模擬方法和模擬參數設置是合理的。

(2) 在加熱過程中，瓶中啤酒存在冷核，表明瓶中啤酒的溫度分布不均勻；冷核的位置隨著加熱的進行向上移動，但不超過灌裝高度的1/2。

(3) 利用L9(34)正交試驗研究啤酒巴氏殺菌的過熱階段影響因素結果表明，噴口處噴淋水湍流強度和噴淋水溫度對瓶中啤酒溫度分布均一性的影響顯著，瓶子的移動速度則影響不明顯。

[1] BUZRUL S. A suitable model of microbial survival curves for beer pasteurization[J]. LWT - Food Science and Technology, 2007, 40(8): 1 330-1 336.

[2] AUGUSTO P E D, PINHEIRO T F, CRISTIANINI M. Using computational fluid-dynamics (CFD) for the evaluation of beer pasteurization: effect of orientation of cans[J]. Food Science and Technology, 2010, 30(4): 980-986.

[3] HORN C S, FRANKE M, BLAKEMORE F B, et al. Modelling and simulation of pasteurization and staling effects during tunnel pasteurization of bottled beer[J]. Food and Bioproducts Processing, 1997, 75(1): 23-33.

[4] ENGELMAN M S, SANI R L. Finite-element simulation of an in-package pasteurization process[J]. Numerical Heat Transfer, 1983, 6(1): 41-54.

[5] TEIXEIRA A A, DATTA A K. Numerically predicted transient temperature and velocity profiles during natural convection heating of canned liquid foods[J]. Journal of Food Science, 1988, 53(1): 191-195.

[6] ZECHMAN L G. Location of the slowest heating zone for natural convection heating fluids in metal containers[J]. Journal of Food Science, 1989, 54(1): 205-209.

[7] KUMAR A, BHATTACHARYA M, BLAYLOCK J. Numerical simulation of natural convection heating of canned thick viscous liquid food products[J]. Journal of Food Engineering, 1990, 55(5): 1 403-1 411.

[8] BHUVANESWARI E, ANANDHARAMAKRISHNAN C. Heat transfer analysis of pasteurization of bottled beer in a tunnel pasteurizer using computational fluid dynamics[J]. Innovative Food Science & Emerging Technologies, 2014, 23: 156-163.

[9] KANNAN A, SANDAKA P C G. Heat transfer analysis of canned food sterilization in a still retort[J]. Journal of Food Engineering, 2008, 88(2): 213-228.

[10] FARID M, ABDUL GHANI A G. A new computational technique for the estimation of sterilization time in canned food[J]. Chemical Engineering and Processing: Process Intensification, 2004, 43(4): 523-531.

[11] GHANI A G A, FARID M M, CHEN X D. Theoretical and experimental investigation of the thermal inactivation of bacillus stearothermophilus in food pouches[J]. Journal of Food Engineering, 2002, 51(3): 221-228.

[12] VARMA M N, KANNAN A. CFD studies on natural convective heating of canned food in conical and cylindrical containers[J]. Journal of Food Engineering, 2006, 77(4): 1024-1036.

[13] ASLAM BHUTTA M M, HAYAT N, BASHIR M H, et al. CFD applications in various heat exchangers design: a review[J]. Applied Thermal Engineering, 2012, 32: 1-12.

[14] NORTON T, SUN Da-wen. Computational fluid dynamics (CFD)-an effective and efficient design and analysis tool for the food industry: a review[J]. Trends in Food Science & Technology, 2006, 17(11): 600-620.

[15] DILAY E, VARGAS J V C, AMICO S C, et al. Modeling, simulation and optimization of a beer pasteurization tunnel[J]. Journal of Food Engineering, 2006, 77(3): 500-513.

[16] 王亮, 于艷艷, 馬曉彬, 等.罐裝啤酒與瓶裝啤酒的巴氏殺菌過程數值模擬[J]. 食品與發酵工業, 2014, 40(4): 42-46.

[17] 洪曉敏, 許薔. 啤酒巴氏殺菌過程的CFD數值分析[J]. 食品與機械, 2016, 32(2): 160-164.

[18] LEWICKI P P, WALCZAK W, GOSS B. Heat transfer in a tunnel pasteuriser. part i. factors affecting the rate of heating of liquid in a bottle[J]. Journal of Food Engineering, 1983, 2(4): 309-322.

[19] 樊延敏, 王長偉. 淺談啤酒殺菌機PU值控制系統[J]. 啤酒科技, 2003(7): 36, 38.

[20] 張凱, 王瑞金, 王剛. Fluent技術基礎與應用實例[M]. 2版. 北京: 清華大學出版, 2010: 33-36.

[21] AUGUSTO P E D, CRISTIANINI M. Numerical simulation of packed liquid food thermal process using computational fluid dynamics(CFD)[J]. International Journal of Food Engineering, 2011, 7(4): 1-22.

[22] BAILEY R T, ELBAN W L. Thermal performance of aluminum and glass beer bottles[J]. Heat Transfer Engineering, 2008, 29(7): 643-650.

[23] ERDOGDU F, TUTAR M. A computational study for axial rotation effects on heat transfer in rotating cans containing liquid water, semi-fluid food system and headspace[J]. International Journal of Heat and Mass Transfer, 2012, 55(13/14): 3 774-3 788.

Study on tunnel pasteurization of bottled beer based on CFD technology

CHENG Ju1ANJia-yan1DONGWeng-yong2WANGYue1

(1.SchoolofBiologicalEngineering,DalianPolytechnicUniversity,Dalian,Liaoning116034,China;2.No.12MiddleSchoolofDalian,Dalian,Liaoning116021,China)

In the industrial processing of beer, pasteurization is necessary the to guarantee its biological stability. The superheat phase (the last part of heating phase) was numerically simulated, in which the slowest heating zone (SHZ) inside the bottle was determined, and then the influence of three parameters on the uniformity of temperature distribution inside the bottled beer during pasteurization, i.e. the temperature of spray water, the turbulent intensity of spray water at the spray nozzles and the running velocity of bottles were analyzed using L9(34) orthogonal experiments. The results indicated that SHZ lay at the bottom of the bottle at the beginning and then moved upward without exceeding 1/2 of filling height of beer. Besides, after heating for 7 min, the temperature and the turbulent intensity of spray water at the spray nozzles affected the uniformity of temperature distribution inside the bottled beer obviously, while the running velocity of bottle had no significant effect on it. Our results showed that the probable optimum condition was 65 ℃ temperature of spray water and 3% turbulent intensity of it at the spray nozzles as well as 3 mm/s running velocity of bottle.

beer; pasteurization; numerical simulation; slowest heating zones (SHZ)

程菊，女，大連工業大學在讀研究生。

王越(1972—)，女，大連工業大學副教授，博士。

E-mail：wydqs2004@163.com

2016—02—16

10.13652/j.issn.1003-5788.2016.08.024