全蛋液双频超声真空干燥的干燥特性及数学模型分析

奥鹏网院作业

全蛋液双频超声真空干燥的干燥特性及数学模型分析全蛋液双频超声真空干燥的干燥特性及数学模型分析

白喜婷1,2，侯亚玲1，朱文学1,2,*，孙国峰1

（1.河南科技大学食品与生物工程学院，河南 洛阳 471023；2.农产品干燥技术与装备河南省工程技术研究中心，河南 洛阳 471023）

摘 要：以全蛋液为研究对象，探讨在双频超声28 kHz+28 kHz的条件下，干燥温度、超声波功率、真空度对全蛋液双频超声真空干燥特性的影响，通过研究碘的释放量来讨论超声波功率对超声空化的影响，以及建立双频超声真空干燥数学模型。结果表明：干燥速率随着干燥温度的升高而升高，但干燥温度不宜太高；提高超声波功率和真空度均有利于提高干燥速率。碘释放量结果表明提高超声波功率可有效提高超声空化效果。有效水分扩散系数与干燥温度、超声波功率和真空度均呈正相关，计算得双频超声真空干燥活化能为9.76 kJ/mol，干燥能耗较低。Page模型能较好反映全蛋液双频超声真空干燥过程。该实验结果可为双频超声真空干燥的研究提供理论依据。

关键词：全蛋液；双频超声；真空干燥；干燥特性；数学模型

全蛋粉是以新鲜蛋液为原料，经多道工序处理后干燥加工得到的粉末状物品[1]。超声真空干燥技术是一种将超声波处理与真空技术相结合的干燥新技术，超声波在液体内的作用主要包括超声波的热作用、机械作用和空化作用[2]，真空技术可为物料干燥提供低温低压的干燥环境。Tekin等[3]研究超声辅助真空干燥对红辣椒干燥速率的影响，发现使用超声真空干燥可有效提高红辣椒的干燥速率，且红辣椒并未发生显著的生物活性化合物降解。Su Ya等[4]使用超声波辅助真空油炸薯片，干燥速率显著提高，且油炸马铃薯片的品质得到极大改善，能耗也有所降低。和大奎等[5]研究超声真空干燥对地黄浸膏的影响，发现使用超声真空干燥地黄浸膏可提高干燥速率。

上述研究展现了超声真空干燥技术在干燥过程中广泛的应用前景。但这些研究以单频超声干燥应用为主，为解决单频超声场中存在的超声不均匀问题，可采用双频超声。有关双频超声的研究主要在提取方面。张凡等[6]采用双频复合超声强化糖液结晶成核发现，使用双频复合超声时糖液结晶成核速率显著提高。李凯等[7]研究双频超声强化合成蔗糖月桂酸单酯工艺，结果表明采用双频组合超声时单酯产率得到了有效提高。Hu Aijun等[8]对马铃薯和小米淀粉使用单、双频超声处理，发现双频超声比单频超声对颗粒改性更有效。

目前双频超声干燥方向的应用较少，因此本实验以全蛋液为研究对象，在双频28 kHz+28 kHz的条件下，研究了干燥温度、超声波功率、真空度对全蛋液双频超声真空干燥特性的影响，通过研究碘的释放量来讨论超声波功率对超声空化的影响，以及建立双频超声真空干燥数学模型，以期为双频超声真空干燥的实际应用提供理论参考。

1 材料与方法1.1 材料与试剂

新鲜鸡蛋购于河南省洛阳市丹尼斯超市，并贮藏于2～4 ℃的冰箱中。根据GB 5009.3—2016《食品安全国家标准 食品中水分的测定》，测得蛋液样品的初始干基水分含量为3.27 kg/kg[9]。

碘化钾 天津市科密欧化学试剂有限公司；实验用水均为超纯水。

1.2 仪器与设备

TD-50002型电子天平 余姚市金诺天平仪器有限公司；DZF-6050型真空干燥箱 上海一恒科学仪器有限公司；KMD-M1型超声波发生器 深圳市科美达超声波设备有限公司；T6新世纪紫外-可见分光光度计 北京普析通用仪器有限公司。

该实验采用的双频超声真空干燥装置示意图如图1所示。本实验所用干燥设备是以马怡童等[10]所用设备为基础进行改良所得。所用双频超声联合真空干燥装置主要由超声系统、反应系统和真空系统组成。超声系统主要包含超声波发生器和超声波换能器。每个超声波发生器最大输出功率为900 W（连续可调），超声波换能器选用两个相同频率的喇叭形陶瓷压电超声波振子，其谐振频率为（28.0±0.5）kHz。超声波换能器可将电信号转换为相应的机械运动，其通过环氧树脂专用胶黏剂安装在反应器上，超声波辐射反应器内的物料，加剧了内部分子振动，从而达到增强真空干燥的效果[11]。不锈钢反应器规格为：长200 mm、宽100 mm、高100 mm。真空系统由真空干燥箱提供。在该实验中采用每组槽式超声换能器的两个频率同时作用，即双频复合形式。

     

图1 双频超声真空干燥装置示意图
Fig. 1 Schematic of dual-frequency ultrasonic vacuum drying device

1.3 方法

1.3.1 样品制备与预处理

将购买的新鲜鸡蛋清洗、消毒后，在室温下将分离的蛋液搅拌均匀并过滤得到原料。将处理好的样品均匀撒在超声波接收装置内，每次实验用量为100 g。

1.3.2 全蛋液双频超声真空干燥

本研究选取干燥温度、超声波功率、真空度为实验因素，分别进行单因素试验。选择超声波频率28 kHz+28 kHz，超声波功率为100 W，真空度为-0.1 MPa，全程超声，干燥温度分别取30、40、50、60、70 ℃，研究温度对干燥特性的影响；选择超声波频率28 kHz+28 kHz，温度为50 ℃，真空度为-0.1 MPa，全程超声，超声波功率分别取40、60、80、100 W，研究超声波功率对干燥特性的影响；选择超声波频率28 kHz+28 kHz，超声波功率为100 W，温度为50 ℃，全程超声，真空度分别取-0.08、-0.09、-0.10 MPa，研究真空度对干燥特性的影响。在干燥过程中每隔20 min，快速取出样品称其质量，直至前后两次质量读数基本不变时，干燥结束，每组实验重复3 次。

1.3.3 水分比和干燥速率的测定

AGPS全称：辅助卫星定位系统。AGPS是在接收 GPS信号的基础上，通过手机定位服务器作为辅助服务器来协助GPS接收器完成测距和定位服务。AGPS定位的优点是定位精度高，在正常的GPS工作环境下可达10m，并且首次捕获GPS信号较短，只需几秒。

水分比（moisture ratio，MR）下降的速率可用于反映样品干燥的快慢，可用式（1）计算[12]。

长期以来，在经济犯罪案件的侦查中有一种错误的认识:经济犯罪案件的侦查一般不存在可供勘查的现场。这种错误认识，致使经济犯罪侦查在理论上游离于传统侦查学之外，在实践中造成一些本应勘验的经济犯罪现场没有及时勘验，使得本能成功侦破的案件未能顺利破获甚至陷入侦查僵局。实际上，包括经济犯罪侦查在内的每一起犯罪案件的侦查都要进行现场勘查，只不过有的更侧重于现场勘验，有的更侧重于现场访问。除此二者之外，现场勘查还有一个现场分析的过程。

   

式中：Mt表示干燥至t时刻的干基水分含量/（kg/kg）；M0表示初始干基水分含量/（kg/kg）；Me表示平衡干基水分含量/（kg/kg）。

从海岸潮汐流中漂来的漂流木还会被困在北太平洋的环流带中，顺着洋流漂向更远的西部。在阿拉斯加西南部的亚北极苔原地区的尤皮克人中流传的圣歌、民谣和传说中，都讲述着浮木对于当地人重要性的故事。用漂流木建造的房屋给他们提供了庇护之所，用木材燃起的篝火照亮了北极的长夜，他们还用这些漂来的木材做成了雕刻精美的萨满面具，丰富了他们的精神生活。在阿拉斯加大陆和西伯利亚之间没有树的阿留申群岛上，当地人用从太平洋西北地区漂来的黄色雪松做成盖世无双的“海豚皮船”，成为克莱默用来进行水上调研和如今体育比赛项目中所用的现代皮划艇的前身。

由于Me与Mt和M0相比可忽略，因此式（1）可简化为式（2）。

   

干燥速率用公式（3）计算。

   

式中：Md表示干基水分含量/（kg/kg）；Md,i表示i时刻的干基水分含量/（kg/kg）；Md,i+1表示i+1时刻的干基水分含量/（kg/kg）。

1.3.4 有效水分扩散系数和活化能的测定

有效水分扩散系数可根据Fick扩散定律（公式（4））计算[13]。

(2)江河漂流生态游。张家界的茅岩河发源于湖南四大水系之一澧水，该水系从湖北流出，湍急回旋地奔腾在山峡里，到永定区温塘镇后，再流入较为开阔的盆谷里，非常适合开展漂流活动。除此之外，慈利县江垭镇的娄江也有着“小三峡”之称，距离武陵源风景区约20 km，距张家界永定区约70 km，该地是典型的喀斯特岩溶地貌，由于地形奇异，悬崖峭立，形成了滩多浪急的水势，在江上漂流的市民众多，惊险刺激。产品的开发离不开宣传，当地政府可以联合开发商举办漂流赛事，邀请媒体发布新闻，吸引更多的游客参与到漂流运动中。

   

式中：Deff表示有效水分扩散系数/（m2/s）；b为物料厚度的一半/m；t表示干燥时间/s；n表示迭代数。

   

对于长时间干燥，式（4）可简化为式（5）。对式（5）两边取自然对数，可得式（6），对ln MR与t进行线性拟合，根据所得斜率求出有效水分扩散系数Deff。活化能表示一个干燥进程发生所需的最低能量，可用阿伦尼乌斯公式（式（7））进行计算。

1.2.3 呼吸道烧伤治疗 临床研究证明,呼吸道烧伤的治疗关键在于保持患者呼吸道畅通,并及时清除坏死、脱落的黏膜组织,避免发生肺部感染。

   
   

式中：Deff表示有效水分扩散系数/（m2/s）；D0表示阿伦尼乌斯公式的指数前因子/（m2/s）；Ea表示活化能/（kJ/mol）；R表示摩尔气体常数，值为8.314×10-3 kJ/（mol·K）；T表示干燥温度/K。

将式（7）两边取自然对数得公式（8），对ln Deff与1/T进行线性拟合，活化能可根据斜率求出。

   

1.3.5 超声空化效应的测定

(3)筒身每砌完一步架(1.25 m)，应检查1次筒身中心线、圆周、垂直度、坡度。测中心线时，用激光垂准仪向上投测中心点，在操作平台上做上标记，并以此作中心。用带水平的尺测定半径，用坡度托线板找坡度。按每砌3～10皮砖检查1次，控制圆周、垂直度和坡度。如发现偏差过大，应及时纠正。

超声波在液体中产生的空化作用是超声强化物料干燥过程的主要因素，因此可通过研究超声的空化产额来说明干燥过程中超声强化效果[6]。本实验选择碘释放法研究超声波功率对超声空化的影响。当超声波作用于含一定溶解空气的碘化钾溶液时，碘离子会转化为碘单质析出，此时可采用紫外分光光度计来测定碘在354 nm波长处的吸光度，通过测定反应释放的碘单质的量来考察超声空化效应的强弱[6]。吸光度越大，说明此时溶液中碘单质的量越多，则超声空化效应越强。参数设计如下：选择28 kHz+28 kHz双频复合超声，功率为40、60、80、100 W，温度25 ℃。将0.2 mol/L的碘化钾溶液100 mL加入到反应器中，两个频率同时开20 min后，取样测量溶液在354 nm波长处的吸光度，每组实验重复3 次。

1.3.6 干燥曲线的拟合

选择国内外常用的10 种食品薄层干燥数学模型对全蛋液双频超声真空干燥实验数据进行拟合[14-23]。具体数学模型见表1。

表1 食品薄层干燥数学模型
Table 1 Mathematical models for describing thin layer drying of foods

     

模型名称 模型方程Lewis MR＝exp（-kt）Page MR＝exp（-ktn）Simplified Fick’s MR＝kexp（-c（t/L2））Henderson and Pabis MR＝aexp（-kt）Logarithmic MR＝aexp（-kt）+c Two term MR＝aexp（-k0t）+bexp（-k1t）Modified two term-I MR＝aexp（-k0t）+（1-a）exp（-k1t）Two-term exponential MR＝aexp（-kt）+（1-a）exp（-kat）Midilli MR＝aexp（-ktn）+bt Wang and Singh MR＝1+at+bt2

根据表1中的10 种干燥模型对全蛋液双频超声真空干燥数据进行拟合，拟合程度用决定系数R2、误差平方和（sum of squares due to error，SSE）和均方根误差（root mean square error，RMSE）来表示，R2越大，SSE和RMSE越小，说明拟合程度越好。R2、SSE和RMSE的计算分别见公式（9）～（11）。

想开展事业单位的政治思想工作，促进其文明建设，必须建立政工队伍，开展有效的政工工作，以适应社会主义经济建设需要。我国特色社会主义的建设与发展，离不开政工工作的有效发展。政治思想工作为其他工作的有效展开提供了思想基础，充分开展事业单位的政治思想工作可以更好促进我国的社会经济发展。因此，事业单位政工工作的有序开展，在极大程度上保障了事业单位的稳定运行，并使得事业单位整体在职员工的创新能力、工作激情等各项能力全面协调发展，提升了事业单位全体职工的思想政治道德高度。

   

式中：N为实验数据个数；MRpred,i为预测水分比；MRexp,i为实测水分比；为实测水分比的算术平均值。

1.4 数据处理与分析

本实验采用Excel 2010软件进行数据整理；采用Origin 8.5软件作图；采用PASW Statistics 18.0统计分析软件的方差分析法对实验数据进行差异显著性分析，P＜0.05表示差异显著；采用Matlab 2014a数学软件进行模拟。

2 结果与分析2.1 全蛋液双频超声真空干燥的干燥特性

2.1.1 温度对全蛋液双频超声真空干燥特性的影响

     

图2 温度对全蛋液双频超声真空干燥特性的影响
Fig. 2 Effect of temperature on dual-frequency ultrasonic vacuum drying characteristics of whole egg liquid

通常温度是决定干燥速率的主要因素，不同温度条件下，干燥温度越高，水分比下降速率越快，干燥所需时间越短。由图2A、B可知，与干燥温度30 ℃的干燥时间180 min相比，干燥温度40、50、60、70 ℃的干燥时间分别降低至160、140、120、100 min，缩短了11.1%、22.2%、33.3%、44.4%。与干燥温度30 ℃的平均干燥速率0.016 47 g/（g·min）相比，干燥温度40、50、60、70 ℃的平均干燥速率分别提高至0.018 18、0.020 67、0.023 33、0.027 72 g/（g·min），分别提高了10.4%、25.5%、41.7%、68.3%。温度越高，干燥过程中所能达到的最高干燥速率越大[24]。这是因为随着干燥温度的升高，物料内部水分子的运动越来越剧烈，进而有利于水分在物料中的移动，且水分的蒸发量增加有利于干燥过程进行[25]；另一方面，干燥介质温度的升高使得干燥介质与样品之间的温度梯度增大，因此热流密度也增加，进而干燥过程加快[26]。

在干燥中期，不同温度处理组之间整体水分比差异显著，说明干燥温度对于缩短干燥时间和提高干燥速率是有效的，进而可提高能量效率。在干燥后期，不同温度处理组之间水分比差异不显著，说明此时干燥已基本结束。当温度为70 ℃时，样品明显焦黄，而温度50 ℃与60 ℃之间差异不明显，若温度过低，干燥物料所需时间则会变长。考虑到能耗问题，选择温度为50 ℃进行后续实验。

2.1.2 超声波功率对全蛋液超声真空干燥特性的影响

本文对该站1～13轴基坑土方开挖过程中围护结构及周边环境变形监测情况进行分析。监测点平面布置图见图2。

     

图3 超声波功率对全蛋液双频超声真空干燥特性的影响
Fig. 3 Effect of ultrasonic power on dual-frequency ultrasonic vacuum drying characteristics of whole egg liquid

由图3A可知，与超声波功率40 W的干燥时间340 min相比，超声波功率60、80、100 W的干燥时间分别降低至220、180、140 min，缩短了35.3%、47.1%、58.8%。与超声波功率40 W的平均干燥速率0.009 10 g/（g·min）相比，超声波功率60、80、100 W的平均干燥速率分别提高至0.013 61、0.016 02、0.020 67 g/（g·min），提高了49.6%、76.0%、127.1%。在超声波强化作用下，干燥速率随着超声波功率的增加而上升，超声功率越大，干燥时间越短。这是因为，一方面，随着超声波功率的增加，物料内部分子之间振动加剧，超声所产生的空化作用、机械效应和热效应等都增强，从而提高了全蛋液内部水分的流动性，降低了水分扩散阻力[27]；另一方面，超声波功率的增加有利于提高水分从样品内部向表面的移动速率，促进热渗透及传热速率，导致更快的水分蒸发速率，从而能够提高传质，进而显著缩短干燥时间[28-29]。

从干燥速率曲线（图3B）可以看出，干燥初期，干燥速率呈现增加趋势，这是因为开始时，物料含水率较高，其内部水分大部分属于自由水，自由水流动性较强，利于超声波的传播与运输，可进一步加快超声波能量转化为热能，提高水分迁移速率，进而提高干燥速率[30]。随着干燥的进行，干燥速率呈现下降趋势，这是因为随着大量自由水被脱除，干燥以脱除半结合水和结合水为主，半结合水和结合水与蛋白质等大分子物质结合较紧密，其流动性较弱，不易被脱除[31]。同时，随着水分含量的下降，物料内部水分流动性减小，超声波的衰减系数逐渐增大，向物料内部传递能量阻力增大，其机械效应和空化效应随之减弱，因此超声的强化作用减弱甚至消失，进而导致干燥速率降低。结果表明，随着水分含量的下降，超声对全蛋液干燥速率的影响减小，与文献[32]所得结论一致。

各个超声波功率处理组之间水分比差异显著，这说明增加超声功率可加快全蛋液内部水分扩散，提高干燥速率；在蛋液干燥后期，各个处理组的干燥速率之间的差距越来越小。选择超声波功率为100 W进行后续实验。

在20世纪八九十年代，我国一些城市开始建立垃圾集中处置的准卫生填埋场，由于这些填埋场设计施工时，普遍采用简易填埋方式和传统填埋作业工艺，对周边环境存在不同程度的污染，难以满足现代卫生填埋场技术标准要求［1］。目前第1批填埋场大多运行了10多a，已逐渐接近封场状态，面临着老旧简易垃圾填埋场封场及渗沥液阻隔控制的难题［2］，亟待在原有垃圾场的基础上利用柔性垂直防渗技术，结合封场，解决渗沥液渗入地下，填埋气外泄，并向四周扩散，污染周边环境的问题。

2.1.3 真空度对全蛋液超声真空干燥特性的影响

     

图4 真空度对全蛋液超声真空干燥特性的影响
Fig. 4 Effect of vacuum degree on dual-frequency ultrasonic vacuum drying characteristics of whole egg liquid

由图4可知，与真空度-0.08 MPa组的干燥时间280 min相比，真空度-0.09、-0.10 MPa的干燥时间分别降低至180、140 min，缩短了35.7%、50.0%。与真空度-0.08 MPa组的平均干燥速率0.010 89 g/（g·min）相比，真空度-0.09、-0.10 MPa的平均干燥速率分别提高至0.016 33、0.020 67 g/（g·min），提高了50.0%、89.8%。真空度较高时，全蛋液的水分比下降快，干燥速率高；真空度较低时，水分比下降缓慢，干燥速率低。这是因为在同一干燥温度下，由于全蛋液内部的绝对压力不变，真空度高时，全蛋液表面水分蒸气压差会变大，蒸发速率会加快，因此干燥速率会增加[33]。另一方面，随着真空度的增加，全蛋液所含水分的沸点也会相应降低，这会导致全蛋液表面水分的蒸发速率加快，进而提高干燥速率[5]。方差分析结果表明，各个真空度处理组之间水分比差异显著，这说明提高真空度有利于提高干燥速率。

2.2 有效水分扩散系数和活化能

表2 不同条件下全蛋液双频超声真空干燥的有效水分扩散系数
Table 2 Effective water diffusion coefficients for dual-frequency ultrasonic vacuum drying of whole egg liquid under different conditions

     

实验条件 Deff/（m2/s）温度/℃30 5.71×10-8 40 5.91×10-8 50 7.18×10-8 60 7.38×10-8 70 9.00×10-8images/BZ_170_1741_1030_1744_1033.png超声功率/W 40 2.21×10-8 60 3.35×10-8 80 4.92×10-8 100 7.18×10-8真空度/MPa-0.08 2.69×10-8-0.09 4.69×10-8-0.10 7.18×10-8

由表2可知，当干燥温度为变量时，与干燥温度30 ℃时的有效水分扩散系数5.71×10-8 m2/s相比，干燥温度40、50、60、70 ℃时的有效水分扩散系数分别提高至5.91×10-8、7.18×10-8、7.38×10-8、9.00×10-8 m2/s，相应提高了3.5%、25.7%、29.2%、57.6%；当超声功率为变量时，与超声功率40 W组的有效水分扩散系数2.21×10-8 m2/s相比，超声功率60、80、100 W组的有效水分扩散系数分别提高至3.35×10-8、4.92×10-8、7.18×10-8 m2/s，相应提高了51.6%和1.23、2.25 倍；当真空度为变量时，与真空度-0.08 MPa组的有效水分扩散系数2.69×10-8 m2/s相比，真空度-0.09、-0.1 MPa组的有效水分扩散系数分别提高至4.69×10-8、7.18×10-8 m2/s，相应提高了74.3%和1.67 倍。由此可以说明，干燥温度、超声波功率和真空度均能改变全蛋液双频超声真空干燥过程，且随着干燥温度、超声波功率和真空度的增加，全蛋液双频超声真空干燥的有效水分扩散系数均有所增加，该实验结果与干燥速率随干燥条件变化规律一致。

实验材料有芹菜、芦蒿、无机碱,实验仪器有电动搅拌机(7312-I型,上海标本模型厂制造)、恒温水浴锅(HH-S,)、旋转粘度计(NDJ-1型，上海昌吉地质仪器有限公司)、电子分析天平(AB204-N型，梅特勒-托利多仪器有限公司).

将ln Deff和1/T曲线进行线性拟合，通过拟合直线的斜率，计算出全蛋液双频超声真空干燥的活化能为9.76 kJ/mol（R2＝0.93）。全蛋液双频超声真空干燥的活化能与文献[10]相比，活化能较低，说明使用双频干燥超声真空干燥，能够减小能耗，达到节能的效果。

2.3 超声波功率对空化效应的影响

从图5可以看出，在其他条件相同的情况下，超声波功率60、80、100 W时释放的碘量分别是超声功率40 W时的1.8、5.2、8.8 倍。这是因为随着超声波功率的增加，溶液与空气的传质面积增大，溶液中的机械扰动增强，溶液内分子间的作用增强，同时更多的空气进入到溶液中，这些变化均使得碘离子更容易转化为碘单质，则表现为A354 nm增加。因此，随着超声波功率的增加，溶液中声空化产额也在增加，超声空化效应增强，更有利于强化物料干燥。该实验结果进一步证明了水分比与干燥时间关系随超声波功率变化规律。

     

图5 超声波功率对空化效果的影响
Fig. 5 Effect of ultrasonic power on cavitation efficiency

2.4 全蛋液双频超声真空干燥的数学模型拟合结果

2.4.1 干燥模型拟合参数

表3 不同超声波功率下各干燥模型的统计参数和系数
Table 3 Statistical parameters and coefficients of various drying models under different ultrasonic powers

     

模型名称 功率/W 模型系数 R2 SSE RMSE Lewis 40 k＝0.005 979 0.919 6 0.177 60 0.099 33 60 k＝0.009 018 0.906 9 0.13758 0.111 84 80 k＝0.012 090 0.916 4 0.108 77 0.109 93 100 k＝0.018 360 0.922 8 0.086 75 0.111 32 Page 40 k＝0.000 119 5 n＝1.753 0.997 0 0.005 55 0.019 23 60 k＝0.000 124 3 n＝1.901 0.998 7 0.001 45 0.013 44 80 k＝0.000 243 2 n＝1.869 0.997 9 0.002 38 0.017 24 100 k＝0.000 336 4 n＝1.972 0.999 5 0.000 47 0.008 87 Simplified Fick’s 40 k＝1.004 c＝0.000 324 L＝0.026 0 0.911 1 0.174 72 0.104 49 60 k＝1.008 c＝0.002 206 L＝0.492 8 0.909 8 0.133 33 0.121 70 80 k＝1.001 c＝0.002 183 L＝0.425 1 0.916 5 0.108 79 0.124 61 100 k＝1.1003 c＝0.002 835 L＝0.392 3 0.923 7 0.085 84 0.131 00 Henderson and Pabis 40 a＝1.143 k＝0.006 81 0.939 8 0.125 56 0.085 94 60 a＝1.137 k＝0.010 23 0.924 1 0.101 96 0.100 98 80 a＝1.128 k＝0.013 51 0.928 9 0.082 29 0.101 42 100 a＝1.1034 k＝0.019 94 0.925 0 0.072 24 0.109 73 Logarithmic 40 a＝1.869 k＝0.002 538 c＝-0.809 4 0.989 4 0.020 91 0.036 15 60 a＝2.380 k＝0.002 863 c＝-1.320 0 0.984 8 0.018 32 0.04512 80 a＝1.812 k＝0.005 372 c＝-0.745 6 0.979 8 0.020 42 0.054 01 100 a＝1.423 k＝0.011 250 c＝-0.356 4 0.958 2 0.033 53 0.081 89 Two term 40 a＝-0.248 1 k0＝1.069 0 b＝1.247 0 k1＝0.007 42 0.959 1 0.090 26 0.077 57 60 a＝-0.324 2 k0＝0.999 1 b＝1.323 0 k1＝0.011 86 0.945 8 0.058 22 0.085 31 80 a＝1.398 2 k0＝0.016 5 b＝-0.399 4 k1＝1.174 0.961 8 0.033 17 0.074 36 100 a＝1.384 0 k0＝0.022 2 b＝-0.385 2 k1＝1.135 0. 975 1 0.003 26 0.028 56 Modified two term-I 40 a＝1.247 0 k0＝0.007 419 k1＝1.097 0.954 0 0.090 26 0.075 11 60 a＝-0.323 1 k0＝1.046 k1＝0.011 86 0.951 8 0.058 22 0.080 43 80 a＝-0.398 4 k0＝1.733 k1＝0.016 51 0.967 2 0.033 18 0.068 84 100 a＝-0.571 0 k0＝1.509 k1＝0.026 99 0.977 0 0.018 50 0.060 82

续表3

     

模型名称 功率/W 模型系数 R2 SSE RMSE Two-term exponential 40 a＝2.107 k＝0.009 87 0.981 8 0.037 98 0.047 26 60 a＝2.177 k＝0.015 51 0.982 1 0.024 02 0.049 01 80 a＝2.203 k＝0.020 76 0.987 7 0.014 21 0.042 14 100 a＝2.275 k＝0.031 83 0.990 2 0.009 47 0.039 73 Midilli 40 a＝0.988 4 k＝0.000 193 2 n＝1.628 b＝-0.000 222 7 0.993 7 0.013 16 0.027 82 60 a＝0.991 1 k＝0.000 171 3 n＝1.808 b＝-0.000 261 4 0.996 9 0.004 16 0.020 39 80 a＝1.001 0 k＝0.000 343 1 n＝1.776 b＝-0.000 207 2 0.998 6 0.001 23 0.014 34 100 a＝0.993 5 k＝0.000 301 6 n＝1.997 b＝-0.000 048 2 0.999 3 0.000 42 0.010 30 0.987 2 0.026 72 0.039 64 60 a＝-0.005 710 b＝4.133×10-6 0.981 9 0.024 26 0.049 26 80 a＝-0.008 076 b＝1.238×10-5 0.978 5 0.024 91 0.055 80 100 a＝-0.013 160 b＝4.124×10-5 0.968 8 0.030 02 0.070 73 40 a＝-0.003 964 b＝2.833×10-6 Wang and Singh

使用Matlab软件，根据表1中所列的数学模型对干燥实验数据进行拟合分析，以干燥条件变量为超声波功率示例，各干燥模型的干燥常数和决定系数R2、SSE和RMSE值如表3所示。Page模型和Midilli模型的R2值均大于0.99，SSE分别为0.000 47～0.005 55和0.000 42～0.013 16，RMSE值分别为0.008 87～0.019 23和0.010 30～0.027 82，均较优。说明Page模型和Midilli模型的拟合效果均较好。对其他干燥条件的实验数据进行拟合分析，结果同上，对比两个模型的各个指标参数，考虑到拟合效果和实际条件，选择Page模型作为全蛋液双频超声真空干燥模型。

碳关税的本质是一种出口税，其征收体现在出口商品的价格变化中，本文以中国出口至美国商品征收碳关税为例，则出口方程可由式(6)

2.4.2 Page模型的求解

表4 不同干燥条件下Page模型的统计参数和系数
Table 4 Statistical parameters and coefficients of the Page model under different drying conditions

     

干燥温度/℃ 功率/W 真空度/MPa Page模型系数 R2 SSE RMSEk n超声波50 100 -0.08 0.000 553 9 1.579 0.999 3 0.001 03 0.008 89 50 100 -0.09 0.000 187 0 1.935 0.999 6 0.000 47 0.007 66 50 100 -0.10 0.000 336 4 1.972 0.999 5 0.000 47 0.008 87 50 40 -0.10 0.000 119 5 1.753 0.997 0 0.005 55 0.019 23 50 60 -0.10 0.000 124 3 1.901 0.998 7 0.001 45 0.013 44 50 80 -0.10 0.000 243 2 1.869 0.997 9 0.002 39 0.017 28 30 100 -0.10 0.000 123 1 2.086 0.999 2 0.000 93 0.011 54 40 100 -0.10 0.000 321 4 1.920 0.999 0 0.001 07 0.012 37 60 100 -0.10 0.000 650 4 1.830 0.998 4 0.001 24 0.015 72 70 100 -0.10 0.000 528 1 1.949 0.999 1 0.000 60 0.012 21

根据Page模型对不同干燥条件下的全蛋液双频超声真空干燥实验数据进行拟合，得到相应的干燥常数如表4所示。Page模型中的干燥常数k和n与干燥温度（T/℃）、超声波功率（P/W）和真空度（V/MPa）相关，为得到干燥温度、超声波功率和真空度对模型的影响，利用PASW Statistics 18.0软件对Page模型中的干燥常数进行二次多项式回归分析。k和n的计算见式（12）和式（13）。

总之，法国在艺术与文化教育方面出台了一系列政策，建立起一套行之有效的工作体制，积累了大量的经验，其将艺术和文化教育作为基本公民权利的思路、做法对我国进一步开展美育具有重要的参考价值。

   
   

采用PASW Statistics 18.0软件中提供的多元线性回归中的向后回归法，剔除不显著的影响因素（P＞0.05），得到Page模型中干燥常数k和n的回归方程，结果如式（14）、（15）所示。

   

2.4.3 Page模型的验证

     

图6 Page模型实验值与预测值比较
Fig. 6 Comparison between experimental values and predicted values from Page model

选择干燥温度40 ℃、超声波功率100 W、真空度-0.10 MPa，干燥温度50 ℃、超声波功率100 W、真空度-0.10 MPa以及干燥温度50 ℃、超声波功率100 W、真空度-0.08 MPa的干燥实验值和最终模型预测值进行验证比较，结果如图6所示。实验值与Page模型预测值的吻合程度较高，说明Page模型的拟合程度较好。因此表明实验所建立的Page模型准确可靠，能够用于预测全蛋液双频超声真空干燥过程中不同实验条件下的水分变化。

但我们从上面所举例句中可以看出，德都方言的宾格、从格、造格、方位格等格语法范畴在青海汉话里都有其相对的语法形式来表现，而与现代汉语相背离，只是这种语法形式还不像蒙古语那样完善，即表示语法意义不同的“格”，起码应该有相应的、比较固定的不同的语法形式，但青海汉话还没有完全区别开来，如相当于阿尔泰语言的“宾格”、“方位格”，青海话里只用[xa](哈)这一个语音形式而没有形式上的固定区别的符号，尽管如此，在青海汉话里“格”语法范畴虽然不够完善，但明确地说明了蒙古语对于青海汉话的影响是深刻的。此外，青海汉话还有一种现象，跟蒙古语的“界限副动词”相类似。例如：

3 结 论

全蛋液双频超声真空干燥特性研究结果表明，干燥速率均随着干燥温度、超声波功率、真空度的升高而升高。有效水分扩散系数与干燥温度、超声波功率和真空度均呈正相关，且全蛋液双频超声真空干燥的活化能较低，可达到节能的效果。碘释放量测试结果表明提高超声波功率可有效提高超声强化效果。对不同干燥数学模型进行拟合分析，结果显示，Page模型能很好地反映全蛋液的干燥过程，且Page模型能较好地预测全蛋液双频超声真空干燥过程中水分比的变化规律。

参考文献：

[1] 迟玉杰, 沈青, 赵英, 等. 提高全蛋粉速溶性的研究[J]. 中国家禽,2016, 38(12): 1-3. DOI:10.16372/j.issn.1004-6364.2016.12.001.

[2] ZHAO Y T, WANG W W, ZHENG B D, et al. Mathematical modeling and influence of ultrasonic pretreatment on microwave vacuum drying kinetics of lotus (Nelumbo nucifera Gaertn.) seeds[J].Drying Technology, 2017, 35(5): 553-563. DOI:10.1080/07373937.2016.1193512.

[3] TEKIN Z H, BASLAR M. The effect of ultrasound-assisted vacuum drying on the drying rate and quality of red peppers[J]. Journal of Thermal Analysis and Calorimetry, 2018. 132(2): 1131-1143.DOI:10.1007/s10973-018-6991-7.

[4] SU Ya, ZHANG Min, ZHANG Weiming, et al. Ultrasonic microwaveassisted vacuum frying technique as a novel frying method for potato chips at low frying temperature[J]. Food and Bioproducts Processing,2018, 108: 95-104. DOI:10.1016/j.fbp.2018.02.001.

[5] 和大奎, 朱文学, 于斌, 等. 地黄浸膏超声真空干燥特性和动力学研究[J]. 江苏农业科学, 2017, 45(13): 157-164. DOI:10.15889/j.issn.1002-1302.2017.13.045.

[6] 张凡, 杨日福, 单佳维, 等. 双频复合超声强化无水葡萄糖溶液结晶成核研究[J]. 声学技术, 2015, 34(6): 515-520. DOI:10.16300/j.cnki.1000-3630.2015.06.009.

[7] 李凯, 蒙丽丹, 苏佳廷, 等. 双频超声强化酯交换合成蔗糖月桂酸单酯工艺[J]. 应用化学, 2016, 45(7): 1224-1227; 1231.

[8] HU Aijun, LI Yang, ZHENG Jie. Dual-frequency ultrasonic effect on the structure and properties of starch with different size[J]. LWTFood Science and Technology, 2019, 106: 254-262. DOI:10.1016/j.lwt.2019.02.040.

[9] 国家卫生和计划生育委员会. 食品安全国家标准 食品中水分的测定: GB 5009.3—2016[S]. 北京: 中国标准出版社, 2016: 1-2.

[10] 马怡童, 朱文学, 白喜婷, 等. 超声强化真空干燥全蛋液的干燥特性与动力学模型[J]. 食品科学, 2018, 39(3): 142-149. DOI:10.7506/spkx1002-6630-201803022.

[11] 杨日福, 耿琳琳, 范晓丹. 超声协同静电场提取黄花菜中总黄酮的研究[J]. 声学技术, 2017, 36(1): 32-37. DOI:10.16300/j.cnki.1000-3630.2017.01.007.

[12] 魏彦君. 南美白对虾超声波辅助热泵干燥动力学及品质特性研究[D]. 淄博: 山东理工大学, 2014: 13-14.

[13] 李叶贝, 任广跃, 屈展平, 等. 马铃薯小麦复合面条热泵干燥特性及数学模型的研究[J]. 中国粮油学报, 2019, 34(10): 7-15.

[14] PIRONE B N, MICHELIS A D, SALVATORI D M. Pretreatments effect in drying behaviour and colour of mature and immature‘Napolitana’ sweet cherries[J]. Food and Bioprocess Technology,2014, 7(6): 1640-1655. DOI:10.1007/s11947-013-1238-x.

[15] ZHAO H B, YANG Z H, TAO Z H. Drying kinetics of continuous and intermittent heat pump drying of green soybean seeds[J]. International Journal of Food Engineering, 2017, 13(11): 1-15. DOI:10.1515/ijfe-2017-0182.

[16] YU H M, ZUO C C, XIE Q J. Drying characteristics and model of chinese hawthorn using microwave coupled with hot air[J].Mathematical Problems in Engineering, 2015, 2015(3): 480752.DOI:10.1155/2015/480752.

[17] BORAH A, HAZARIKA K. Simulation and validation of a suitable model for thin layer drying of ginger rhizomes in an induced draft dryer[J]. International Journal of Green Energy, 2017, 14(13): 1150-1155. DOI:10.1080/15435075.2017.1369418.

[18] KADAM D M, DHINGRA D. Mass transfer kinetics of banana slices during osmo-convective drying[J]. Journal of Food Process Engineering, 2011, 34(2): 511-532. DOI:10.1111/j.1745-4530.2009.00373.x.

[19] ERTEKIN C, FIRAT M Z. A comprehensive review of thin-layer drying models used in agricultural products[J]. Critical Reviews in Food Science and Nutrition, 2017, 57(4): 701-717. DOI:10.1080/1040 8398.2014.910493.

[20] DA SILVA B G, FILETI A M, TARANTO O P. Drying of brazilian pepper-tree fruits (Schinus terebinthifolius Raddi): development of classical models and artificial neural network approach[J]. Chemical Engineering Communications, 2015, 202(8): 1089-1097. DOI:10.1080/00986445.2014.901220.

[21] G Ó RNICKI K, KALETA A, BRYŚ A, et al. Thin-layer drying of sawdust mixture[J]. Polish Journal of Chemical Technology, 2016,18(4): 65-70. DOI:10.1515/pjct-2016-0072.

[22] ARGO B D, SANDRA S, UBAIDILLAH U. Mathematical modeling on the thin layer drying kinetics of cassava chips in a multipurpose convective-type tray dryer heated by a gas burner[J]. Journal of Mechanical Science and Technology, 2018, 32(7): 3427-3435.DOI:10.1007/s12206-018-0646-2.

[23] NADHARI W N A W, HASHIM R, SULAIMAN O, et al. Drying kinetics of oil palm trunk waste in control atmosphere and open air convection drying[J]. International Journal of Heat and Mass Transfer,2014, 68: 14-20. DOI:10.1016/j.ijheatmasstransfer.2013.09.009.

[24] 巨浩羽, 肖红伟, 郑霞, 等. 干燥介质相对湿度对胡萝卜片热风干燥特性的影响[J]. 农业工程学报, 2015, 31(16): 296-304.DOI:10.11975/j.issn.1002-6819.2015.16.040.

[25] 刘云宏, 苗帅, 孙悦, 等. 接触式超声强化热泵干燥苹果片的干燥特性[J]. 农业机械学报, 2016, 47(2): 228-236. DOI:10.6041/j.issn.1000-1298.2016.02.030.

[26] DO NASCIMENTO E M G C, MULET A, ASCHERI J L R, et al.Effects of high-intensity ultrasound on drying kinetics and antioxidant properties of passion fruit peel[J]. Journal of Food Engineering, 2016,170: 108-118. DOI:10.1016/j.jfoodeng.2015.09.015.

[27] 朱宗升. 超声波-热泵联合干燥热力学分析及青豆种子干燥特性研究[D]. 天津: 天津大学, 2015: 92-93.

[28] SABAREZ H T, GALLEGO-JUAREZ J A, RIERA E. Ultrasonicassisted convective drying of apple slices[J]. Drying Technology,2012, 30(9): 989-997. DOI:10.1080/07373937.2012.677083.

[29] 刘云宏, 孙畅莹, 曾雅. 直触式超声功率对梨片超声强化热风干燥水分迁移的影响[J]. 农业工程学报, 2018, 34(19): 284-292.DOI:10.11975/j.issn.1002-6819.2018.19.036.

[30] AZOUBEL P M, DO AMPARO MELO BAIRNA M, DA ROCHA AMORIM M, et al. Effect of ultrasound on banana cv Pacovan drying kinetics[J]. Journal of Food Engineering, 2010, 97(2): 194-198.DOI:10.1016/j.jfoodeng.2009.10.009.

[31] LIU Y H, SUN Y, YU H C, et al. Hot air drying of purple-fleshed sweet potato with contact ultrasound assistance[J]. Drying Technology,2017, 35(5): 564-576. DOI:10.1016/j.jfoodeng.2009.10.009.

[32] UDOMKUN P, ARGYROPOULOS D, NAGLE M, et al. Single layer drying kinetics of papaya amidst vertical and horizontal airflow[J].LWT-Food Science & Technology, 2015, 64(1): 67-73. DOI:10.1016/j.lwt.2015.05.022.

[33] 杨飞, 何正斌, 赵阳, 等. 超声波-真空协同干燥自由水迁移速率[J].东北林业大学学报, 2012, 40(12): 103-107. DOI:10.3969/j.issn.1000-5382.2012.12.024.

Drying Characteristics and Mathematical Modeling of Whole Egg Liquid Dried by Dual-Frequency Ultrasonic Vacuum Drying

BAI Xiting1,2, HOU Yaling1, ZHU Wenxue1,2,*, SUN Guofeng1
(1. College of Food and Bioengineering, Henan University of Science and Technology, Luoyang 471023, China;2. Agricultural Product Drying Technology and Equipment of Henan Engineering Technology Research Center, Luoyang 471023, China)

Abstract: Whole egg liquid was dried by double-frequency (28 kHz + 28 kHz) ultrasonic vacuum drying. The effects of drying temperature, ultrasonic power, and vacuum degree on the drying characteristics of whole egg liquid was investigated.Also, the effect of ultrasonic power on the cavitation yield was explored by measuring the amount of iodine released.A mathematical model for the drying process was established. The results showed that the drying rate increased with increasing drying temperature, but the drying temperature should not be too high. Increasing the ultrasonic power and the vacuum degree was conducive to increasing the drying rate. The results of iodine release showed that increasing the ultrasonic power effectively improved the ultrasonic enhancement effect. The effective moisture diffusion coefficient was positively correlated with the drying temperature, ultrasonic power and vacuum degree. The calculated activation energy of dual-frequency ultrasonic vacuum drying was 9.76 kJ/mol, and the energy consumption for drying was low. The Page model could better reflect the dual-frequency ultrasonic vacuum drying process of whole egg liquid. The results from this study provide a theoretical basis for the study of dual-frequency ultrasonic vacuum drying.

Keywords: whole egg liquid; dual-frequency ultrasound; vacuum drying; drying characteristics; mathematical modeling

收稿日期：2019-12-21

基金项目：河南省高校科技创新团队支持计划项目（17IRTSTHN016）

第一作者简介：白喜婷（1966—）（ORCID: 0000-0002-0578-1980），女，副教授，硕士，研究方向为农产品加工与技术。E-mail: bxting@126.com

*通信作者简介：朱文学（1967—）（ORCID: 0000-0002-6362-0966），男，教授，博士，研究方向为农产品加工工程技术、干燥技术与设备。E-mail: zhuwenxue67@126.com

DOI:10.7506/spkx1002-6630-20191221-249

中图分类号：TS253.43

文献标志码：A

文章编号：1002-6630（2020）11-0157-08

引文格式：白喜婷, 侯亚玲, 朱文学, 等. 全蛋液双频超声真空干燥的干燥特性及数学模型分析[J]. 食品科学, 2020, 41(11):157-164. DOI:10.7506/spkx1002-6630-20191221-249. http://www.spkx.net.cn

BAI Xiting, HOU Yaling, ZHU Wenxue, et al. Drying characteristics and mathematical modeling of whole egg liquid dried by dual-frequency ultrasonic vacuum drying[J]. Food Science, 2020, 41(11): 157-164. (in Chinese with English abstract)DOI:10.7506/spkx1002-6630-20191221-249. http://www.spkx.net.cn