基于自由体积理论分析面团膜中的水分扩散特性

赵学伟，杨 毅，冯志强，赵 琼，刘兴丽，李望铭，张 华,

（1.郑州轻工业大学食品与生物工程学院，河南 郑州 450001；2.河南省冷链食品质量安全控制重点实验室，河南 郑州 450001；3.三全食品股份有限公司，河南 郑州 450044；4.郑州科技学院食品科学与工程学院，河南 郑州 450064）

对扩散特性研究是进行水分传质相关过程优化和品质控制的基础，扩散系数是表征扩散特性的主要物理量。近年来，动态水分吸附（dynamic vapor sorption，DVS）分析仪在研究水分扩散特性方面得到广泛应用。该法具有快速、准确等优点，但也存在局限性：由于气流速度相对较小，在样品尺寸较大时需要很长时间才能达到吸附平衡。

通常采用Fick第二扩散定律描述水分扩散过程，其实非Fick扩散在食品中也广泛存在。可以采用Langmuir扩散、黏弹性扩散、对数边界条件或分数阶偏微分方程等方法对非Fick扩散过程进行分析[1]，其中水分依赖性扩散系数是一种常见的方法。食品中关于扩散系数随水分含量变化的研究很多，但多数仅是给出变化趋势，并未建立两者之间的定量关系式[2-5]。关于水分依赖性扩散系数的表示式，通常认为扩散系数是水分含量的指数函数[6]。这种经验关系式往往对水分含量随时间变化的拟合较差，尤其在水分含量变化幅度较大时。因此，如何基于一定的理论假设建立扩散系数与水分含量之间的关系，就显得十分重要。

自由体积理论（free volume theory，FVT）认为物料体系内存在的自由体积是溶剂（本实验指水分）得以扩散的基础，一个体系内的自由体积越多，水分扩散越快。该理论基于溶剂和固相的自身特性表达扩散系数与水分含量的关系，具有较好的理论基础[7]。自由体积扩散模型有多种形式，其中Vrentas &Duda自由体积理论（VD理论）[8]在食品研究中的应用最广。黄秀玲等[9]基于该理论综述了小分子溶剂在食品包装材料中的扩散及影响因素。Ribeiro等[10]采用该模型在水分质量分数为0%时的极限形式研究聚糊精内的水分扩散。该模型还被用于模拟面包皮的水分吸湿[11]、麦芽糊精及其与乳清蛋白或酪蛋白混合体系[12]、花菜[13]、蔗糖、海藻糖、麦芽糊精以及乳糖膜[14-15]的干燥。

FVT的应用方式可分为两类：正向应用，既所有参数已知，据此计算扩散系数[14]或进一步模拟水分传递过程[13]；反向应用求参数，即先求出不同水分含量时的扩散系数，然后拟合该值求出FVT参数的最优值[10-11]。对于后一种情况，如何正确确定扩散系数与水分含量之间的关系很重要。规则区域法可以达到这个目的，但仅适用于干燥过程，且计算复杂[12,15]。对于DVS测定结果，通常认为一个相对湿度（relative humidity，RH）步幅内扩散系数为定值，根据多步RH的DVS实验结果确定扩散系数与水分含量的关系[10-11]。实际上，这仅是一个近似处理，如此求得的扩散系数只是一个表观值。直接将测得的水分含量与其预测值比较确定FVT参数，才是一个更为严谨的方法。

本研究以面团膜为研究对象，主要基于以下考虑：面团易做成膜状样品，便于后续的DVS测定；采用膜状样品，可以将实际上的三维传质问题简化为一维问题，便于后续的模拟计算，尤其在反演模拟求参数时，可以大大降低计算量。

基于上述分析，本研究首先采用DVS测定面团膜在20、40 ℃时9 个RH水平下的水分吸湿过程，然后根据FVT所确定的扩散系数与水分含量、温度的关系，直接采用2 个温度下的水分含量随时间变化的数据进行拟合；针对模拟结果的不足，对FVT模型进行适当改进。然后，从FVT角度分析扩散系数随水分含量变化的规律和机理。最后对低水分阶段的扩散机理进行深入讨论。

1 材料与方法

1.1 材料

金苑面粉 市购。

1.2 仪器与设备

Intrinsic型DVS分析仪 英国Hiden公司；电显千分尺中国Syntek公司。

1.3 方法

1.3.1 扩散过程描述

如果认为水分的传质过程只有扩散，水分浓度梯度为传质动力，通常采用Fick第二扩散定律描述水分含量随位置和时间的变化。对于膜状样品，如果膜的长、宽尺寸比厚度大得多，则可以简化为一维问题，见式（1）[16]。

式中：c为水分质量浓度/（kg/m3）；t为时间/s；D为扩散系数/（m2/s）；x为膜厚度方向上距中心线的距离/m。

对于一个具体的传递过程，如果物体与周围流体介质之间发生对流传质，则边界条件表示为式（2）。

式中：l为膜的厚度/m；hm为对流传质系数/（m/s），本研究采用前人针对DVS测得的对流传质系数（8×10-3m/s）[17]；cv.air为周围空气中水蒸气的质量浓度/（kg/m3）；cv.surf为与食品表面水分含量处于平衡状态的水蒸气质量浓度/（kg/m3）。

假设水蒸气是理想气体，则有PVV＝nRT，其中，Pv为水蒸气分压/Pa，V为气体体积/m3，n为V体积内水蒸气的物质的量；R为气体常数（8.314 J/（mol·K））；T为温度/K；则有其中，M为水分子摩尔质量（18 g/mol）；aw为水分活度；Pv.sat为饱和蒸汽压/Pa。则式（2）转变为式（3），其中，RHair为周围空气的相对湿度/%；aw.surf为与cv.surf达到平衡时的水分活度。

1.3.2 自扩散系数与互扩散系数的关系

Fick扩散方程中的扩散系数实际上是互扩散系数，FVT给出的是自扩散系数，对于双组分体系，两扩散系数之间存在如下关系（式（4））[8]：

式中：D为组分1、2的互扩散系数；D1、D2分别为组分1（本研究指水）、2的自扩散系数；ω1、ω2分别为组分1、2的质量分数/%；μ1为组分1的化学势。

对于大分子（D2）-溶剂（D1）体系，由于D1比D2大得多，Vrentas等[8]指出在溶剂浓度接近0时，方程（4）可以近似表示为：

通常认为在溶剂浓度高时该式也适用。

水分的化学势μ1可用式（6）表示，其中为纯水的化学势，代入式（5），同时考虑到ω2=1-ω1，则得式（7）：

引入参数Q，如式（8）所示，则有D＝D1Q，可以看出，互扩散系数是对自扩散系数进行热力学校正的结果，Q为热力学因子。

如果D1与水分含量、温度的关系已知，aw与水分质量分数的关系可由等温吸附模型确定，则对于一个恒温过程，可以由式（7）计算出一定温度和水分含量下的扩散系数D。本实验如无特殊强调，所说的扩散系数通常指互扩散系数。

1.3.3 基于FVT的自扩散系数

与传热不同，传质涉及到物质实体在介质中的空间移动，所以介质中一定要有空隙才能发生传质，这是FVT的基本假设。FVT扩散模型也有多种形式，主要介绍VD理论[8]。不论是单组分还是混合体系，其体积包括3 部分：分子实际占有体积、距分子较近但未被占有的间隙自由体积和距分子较远的孔自由体积。间隙自由体积内由于存在较高的斥力，不能用于传质；孔自由体积由于能垒较低可用于溶剂分子的扩散。因此，如无特殊说明，本实验中的自由体积均指孔自由体积。要发生扩散，还必须满足：孔的大小要达到扩散分子的大小，且扩散分子要有足够的能量。基于上述思想，Vrentas等[8]就大分子体系内溶剂小分子的自扩散系数得出如下表达式：

式中：Tg1、Tg2、Tgm分别为溶剂、大分子、两者混合体系的玻璃化转变温度/K；K11/γ1、K21分别为溶剂的自由体积参数，K12/γ2、K22分别为大分子的自由体积参数。

对于一个恒温吸湿过程，水分含量逐渐增大，混合体系可能会经历从玻璃态到橡胶态的转变。处于橡胶态时，固相内的自由体积由式（13）计算。处于玻璃态时，由于玻璃态属于热力学非平衡状态，这就导致与同等条件下达到平衡态时相比多出一部分自由体积FH2g/γ2，这部分被冻结的自由体积也会用于扩散，这时大分子内的总自由体积用式（14）计算[7]。

式（7）、（9）～（15）中共有17 个参数，将K11/γ1和K12/γ2看作复合参数，则独立的参数有15 个。这些参数都是基于材料自身的物性。关于水的参数，其值相对恒定，不同文献的差别不大。对于食品大分子，由于组分复杂，相关实验测定数据又少，计算较为困难。在食品相关研究中，一般先求出不同水分含量时的扩散系数，然后反演拟合求出参数。

1.3.4 面团膜的制作

制作水分质量分数约40%的面团，放置约1 h，让面筋充分形成。取少量面团，用手碾成直径约4 mm的长条，放在光滑的塑料板上，用不锈钢棒沿长条的轴线方向往复辊压多次，使中间区域的面团在塑料板上形成一层很薄的膜。室温自然干燥，约1 h即可达到半干状态，将面团膜分割成约1 cm见方的小块，再自然干燥约3 h得到面团膜样品。

1.3.5 吸湿实验

采用Intrinsic型DVS分析仪进行测定。测定前用电显千分尺测定样品膜的长宽和厚度，本研究采用的面团膜样品厚度在96～99 μm之间。由于DVS样品盘不镂空，采用不锈钢丝自制一个三脚支架，放在DVS的样品盘上，将样品膜放在支架上，以保证样品两面的吸湿基本同步。在20、40 ℃两个温度下测定，RH设定为0%，使样品在DVS测定腔内充分干燥，然后RH从0%开始以10%的步幅渐次升高到90%，平衡标准：连续5 min内样品质量的变化速率低于0.001%/min。每个实验重复1 次，求平均值后进行后续的数据分析。

采用COMSOL5.4a（COMSOL Inc.,Burlington,MA）软件中的稀相传质模块，模拟水分扩散过程；采用其优化模块，反演确定模型参数，反演时采用蒙特卡洛优化算法，最优容差为0.001。优化的目标函数OF见式（17）：

式中：n为吸湿实验中c-t数据对的数量；cnum、cexp分别为c的模拟值、测定值。

2 结果与分析

2.1 吸湿等温模型

根据DVS数据计算不同时间点对应的水分含量，取最后10 个数据的平均值作为平衡水分。20 ℃时，在RH 0-10吸湿阶段没有达到真正平衡，其余各步基本达到平衡。采用Peleg提出的公式M（t）=外推求RH0-10（数据表示RH值范围0%～10%，下同）段的平衡水分，其中M0、M（t）、ME分别为吸湿开始时、t时、平衡时的水分含量；k为拟合参数。可以采用Chung-Pfrost等温吸湿模型[18]模拟水分活度与干基水分含量之间的关系，采用水分质量分数代替干基水分含量，见式（18）。此处理便于求Q值，式（18）变换可得式（19），结合式（8）即可求出Q，其中a、b均为Chung-Pfost模型拟合参数。面团膜的吸湿平衡水分及其模型拟合结果如图1所示。

图1 两个温度下面团膜的吸湿平衡水分及其模拟结果Fig.1 Water absorption isotherms of dough film at two temperatures and its modeling results

面团膜的等温吸湿属于典型的Langmuir II型吸附。采用饱和盐水法40 ℃下对杜伦小麦粉进行测定，aw为0.216时的平衡水分质量分数约为14%（干基），aw为0.836时的平衡水分质量分数约为17.5%[19]。与此相比，本研究测得的平衡水分质量分数在低RH时偏低，这可能归因于测定方法的差异。采用DVS测定时，尽管可以设定一个吸湿平衡标准，但吸湿并没有真正达到平衡；在低RH时，很低的水分扩散系数又会加重由此带来的影响。采用DVS法测得的杜伦小麦粗粉的平衡水分质量分数更低[4]，同样是采用DVS法，面包皮[11]、淀粉膜和面筋膜[3]的吸湿平衡水分质量分数与本研究的结果基本一致。Chung-Pfost模型的拟合参数为：20 ℃时，a=4.727、b=19.327；40 ℃时，a=4.933、b=21.558。

2.2 面团膜内的非Fick扩散

采用对流边界条件，认为每RH步幅吸湿过程中的扩散系数为定值，根据式（1）反演确定扩散系数。结果发现，在中等水分质量分数范围内的拟合度较高，但对RH0-10、RH80-90的拟合效果不佳，见图2。为了便于结果展示，对吸水量进行归一化处理，以水分比表示。可以看出，RH80-90吸湿过程中，对初始阶段的预测值偏高，对后段预测值偏低；而对于RH0-10的吸湿过程，初始段的预测值偏低而后段偏高。可能由于相同RH步幅，中等RH水平时水分差较小（图1），这时恒定扩散系数假设带来的偏差不明显。

图2 RH0-10及RH80-90的吸湿过程及其Fick模型拟合Fig.2 Water absorption process by dough film at different RH levels and its modeling resulting using the Fickian law

关于食品中水分的非Fick扩散，对于104 μm厚的淀粉膜在RH0-10内的吸湿，Fick第二定律对吸湿初始阶段的预测值偏高，而对后段的预测值偏低[3]。软麦粉在RH60-65、RH60-75、RH60-85范围内的吸湿[5]，以及面包皮在RH60-70范围内的吸湿[11]，都观察到同样的现象。随水分含量增高扩散系数加速增大，在淀粉膜和面筋膜中也观察到类似的变化趋势[11]，但本研究在高RH时D值更大，可能由于样品未熟化。

2.3 基于自由体积理论的水分扩散特性分析

2.3.1 FVT参数取值

将前人采用DSC测定杜伦麦粗粉[24]、动态力学热分析法[25]以及差示扫描量热法[26]测定小麦粉得到的玻璃化转变温度汇总于图3，以式（16）进行回归，得到Tg2=412.6 K，A=595.5 K，与文献[11]中Tg2=415 K较为接近，但A值相差较大。

图3 小麦粉玻璃化转变温度汇总Fig.3 Glass transition temperatures for wheat flour reported in the literature

可以采用WFL 常数C1、C2，由＝2.303C1C2，K22＝C2计算出K12/γ2和K22。WLF普适常数C1＝17 K、C2＝51 K可用于低分子糖，但对于食品大分子来说，其WLF常数与普适常数有很大差别[10]。因此，本研究中将K22和K12/γ2作为待估参数。ξ和α2g也作为待估参数。其余各参数的取值汇总于表1，该表同时给出了待估参数的优化结果。

表1 FVT参数的取值Table 1 Values of parameters of free volume theory

2.3.2 拟合效果描述

根据FVT在最优参数下对该过程的模拟结果如图4所示。可以看出，根据VD理论计算的扩散系数能够在大部分水分质量浓度范围内较好地模拟面团膜的吸湿过程，但在低水分区域（主要集中在RH0-10吸湿段），模型的预测值偏低，尤其是20 ℃时的吸湿。在高水分区域可能出现预测值偏高的现象。

图4 面团膜在40（A）、20℃（B）的吸湿及其基于FVT的模拟结果Fig.4 Water absorption by dough films at 40 (A) and 20℃ (B) and its modeling results based on the free volume theory

本研究针对面团膜求得K22=141 K，Meinders等[11]拟合扩散系数求得K22=149 K，与本研究的结果十分接近。这时K22-Tg2=-271.6 K。面筋、淀粉的WLF常数C2约为30 K[22]，由此计算出的K22-Tg2值（-395、-540 K）比本研究低得多。Meinders等[11]将这种不一致归因于体系中可能存在低分子糖，因为低分子糖的K22-Tg2=-270 K[20]。

本研究求得K12/γ2=2.207×10-4cm3/g，比前人根据淀粉和面筋的WLF常数[22]估算出的K12/γ2值1.4×10-3cm3/g[11]低。分析显示，根据该值计算出的自由体积明显太高，说明该参数的取值偏高。另外，关于食品体系的多项研究中，该值为10-4cm3/g数量级[11-14]。本研究拟合求得面团膜的热膨胀系数α2g=7.464×10-5K-1，约为橡胶态时的1/4，这个比值基本合理[28]，Meinders等[11]假设α2=1×10-4K-1，为α2g的10 倍。

图4显示，FVT对2 个温度下RH0-10吸湿段的预测值偏低。从FVT的角度看，在低水分段又出现了另外的自由体积。如何计算这部分自由体积，现有文献中没有相关论述。本研究提出用式（20）计算，其中ω*、β为待定系数。在式（13）中再加入该自由体积项，拟合RH0-10段吸湿数据，得到ω*、β的优化值分别为0.037、0.019。改进后的拟合情况见图4插图。

2.3.3 扩散系数随水分含量的变化

为进一步了解水分在面团膜内的扩散过程，以表1中的参数，根据式（9）、（8）以及式（7）分别计算水分质量分数0.01～0.2时水分的自扩散系数D1、热力学因子Q以及互扩散系数D的值，40、20 ℃下的结果如图5所示。图中同时给出了面团膜在水分变化过程中如果始终处于热力学平衡时的自扩散系数。Q为热力学校正因子，其曲线形状主要是由吸湿等温线的形状决定，是偏心的钟形曲线。一般在0.01＜ω＜0.1范围内Q值大于1。D是Q对D1校正的结果，因此D与D1的变化趋势类似，分为3 个阶段：在水分很低时扩散系数随水分含量变化较小；中间区域，扩散系数随水分含量快速增大；后期，扩散系数增大更快（图中纵坐标采用的是对数坐标）。对于20 ℃，橡胶态基本处于观察区域之外，仅显示出2 个阶段。由于玻璃体内挟裹有额外的自由体积，实际的自扩散系数比热力平衡的值大得多。温度高时扩散系数较大，根据式（9），温度升高从两个方面导致扩散系数增大，一方面是提高水分子热运动的动能，从而增大逃离现有空穴的几率；另一方面，提高温度可以增大体系内的自由体积。

图5 40℃（A）和20℃（B）时自扩散系数、互扩散系数及热力学因子随水分质量分数的变化Fig.5 Variations in self-diffusivity coeffciient,mutual diffusivity coeffciient and thermodynamic factors with water content at 40 (A) and 20℃ (B)

总体来说，水分在面团膜内的自扩散系数值在10-14～10-11数量级范围内。比较相关小麦胚乳（1.8×10-10～1.2×10-9m2/s，水分含量20%～36%）[29]、小麦淀粉（0.7×10-11～2.8×10-11m2/s，水分含量26%～37%）[30]的研究结果，根据图5中D1的变化趋势，考虑水分含量对D1的影响，大致可以判定这些结果与本研究的相差不大。

关于扩散系数随水分含量变化的趋势，综合分析大致可以分为3 类。第1类，单调增大，采用FVT的文献中多是如此[11-14]，其他如淀粉膜和面筋膜[4]，D值变幅一般较大。第2类，钟形曲线，如杜伦麦粗粉[5]、淀粉颗粒[2]。第3类，先维持基本不变然后下降，如小麦粉及其主要组分[3]。后2 种趋势所涉及的文献中均采用DVS测定水分吸附过程，且D值变幅更小。如果在根据水分吸附数据推断D值时就认为D是水分含量的函数，那么得到的趋势线大概率是单调增大。本研究结果与第1类接近，但在低水分区域略有不同。

水在面团膜内的扩散系数主要集中在10-13～10-12数量级范围内，与淀粉粒（30 ℃，在10-13m2/s数量级）[31]、面筋膜（25 ℃，0.6×10-13～7×10-13m2/s）[4]的结果基本相当，但高于小麦粉（在10-15m2/s数量级）和小麦淀粉（在10-16m2/s数量级）[3]，而低于小麦粗粉（在10-12m2/s 数量级）[5]、杜伦麦粗粉聚集体（在10-10m2/s数量级）[32]。

2.3.4 自由体积随水分含量的变化

根据式（10）计算40、20 ℃时面团膜内的自由体积，如图6所示，由于体系的自由体积来自水和大分子，以及相态变化所挟裹的自由体积。热力平衡态时的自由体积来自水分和大分子的自由体积之和，在水分质量分数大于0.025时该值主要来自水分。体系进入玻璃态时，挟裹的自由体积开始出现，相比于40、20 ℃时玻璃化转变对自由体积的贡献更大。极低水分含量时，除玻璃化转变外，出现另一相态转变，导致自由体积的进一步增大。

图6 40℃（A）和20℃（B）时面团膜内的总自由体积随水分质量分数的变化及其分解Fig.6 Changes in total free volume in dough films with water content and its decomposition at 40 (A) and 20℃ (B)

面团膜内的总自由体积在0.035～0.07 cm3/g之间，温度升高导致自由体积增大，但水分含量变化带来的影响更大。由图6可以看出，水分含量在很低水平时其变化对总自由体积的影响不大，在玻璃态时影响较大，在橡胶态时影响最大。根据文献[12]中的FVT参数计算发现，40 ℃时水分质量分数从0.01增大到0.20时，麦芽糊精及其与乳清蛋白或酪蛋白混合体系内的自由体积从0.0031 cm3/g增大到0.0627 cm3/g。与本研究的结果较为接近。根据文献[11]中的FVT参数，40 ℃时面包皮内的自由体积由0.0783 cm3/g增加到0.1261 cm3/g，该值明显较大。Benczedi[33]根据淀粉-水混合体系的平均场状态方程计算出了25 ℃时无定型态淀粉中水分浓度与自由体积的对应值，本实验结合Benczedi等[34]给出的该体系中水分浓度与质量分数的对应值，计算出不同质量分数时的自由体积，发现在水分质量分数小于0.043时自由体积随水分质量分数减小几乎线性增大到0.016 cm3/g（对应的水分质量分数为0.0042），在水分质量分数大于0.043时，自由体积基本不变，在0.001 cm3/g到0.002 cm3/g之间，FVT难以解释这种现象，但这个临界值0.043与本研究ω*=0.037比较接近。麦芽聚糖中自由孔的大小在水分质量分数约为0.07时出现最小值[23]。另有研究显示，孔的大小与麦芽聚糖的比容成正比[35]，则可以推论出自由体积在低水分区域内随水分含量的减小而增大，这与本研究的观察结果一致。

3 讨论

3.1 低水分含量时的水分扩散系数变化特点

结果显示，扩散系数随水分含量降低而降低，达到一定程度后略有升高。将前人关于食品体系中水分扩散系数变化的相关数据汇总，见图7，涉及淀粉粒[31]、面包皮[11]、淀粉膜[4]、直链淀粉[36]和麦芽糖[37-38]。图中同时给出本研究所得面团膜在40 ℃时的扩散系数。相关研究[11,37-38]显示，在低水分区域，随水分含量降低有一个D值减速降低甚至升高的阶段。由于计算策略的不同，可以在RH0-10范围内求出多个D值，从而观察到这种现象，Meinders等[11]测定了RH0-5内的D值，也观察到了这种现象。Marousis等[39]采用干燥法测定淀粉粉体的表观扩散系数，D值较大，变幅却较小，见图7插图。在水分质量分数低于0.05时，仍能观察到D值随水分含量减小而增大的趋势。另外，采用核磁共振成像技术测定直连淀粉中水的扩散系数，在水分质量分数大于约0.24时，扩散系数加速增大；而对于40 ℃的面团膜，在水分质量分数大于约0.17也有类似现象；对于淀粉膜也有类似现象，对于麦芽糖，此现象于更低的水分含量时出现，见图7。

图7 不同食品材料中的水分扩散系数随其质量分数的变化汇总Fig.7 Variation in water diffusivity with water content for different food materials reported in the literature

3.2 低水分时体系变化

对于大分子体系，温度或水分含量的变化除引起玻璃化转变外，还可能引起多种相态转变，包括β-转变、γ-转变等。根据图3可以判定在图7中面团膜在水分质量分数小于约0.04的区域不可能发生玻璃化转变，因为体系此时已经处于玻璃态。湿基水分质量分数5%的糊化蜡质玉米淀粉在20 ℃发生β-转变[40]，干基水分质量分数5%的豌豆子叶在5～35 ℃间发生β-松弛[41]。本研究是在20、40 ℃下进行，结合前面推断的水分质量分数临界值ω*=0.037，推断面团膜在低水分区域发生了β-转变。多糖中葡萄糖残基C-5位置上羟甲基的旋转[42]或单体单元围绕糖苷键的旋转[43]是产生β-松弛的原因。

3.3 从β-转变到自由体积、水分扩散分析

水分子在玻璃态葡萄糖内通过2 种方式实现扩散[44]：1）跳到邻近水分子的位置；2）葡萄糖分子发生旋转，腾出位置，水分子跳至该空位。有研究证实对于聚碳酸酯玻璃体，β-松弛可能与WLF自由体积有关[45]。前面提到β-松弛的分子根源是某些侧链或结构单元的旋转，这些“旋转”也确能产生自由体积[46]，这将利于扩散发生。实际上，β-转变也是一种相态转变，与玻璃化转变一样，也没有达到热力学平衡，同样会挟裹一定的自由体积。

Preda等[47]根据β-松弛时间与扩散系数线性相关，认为β-松弛决定水分在聚酰胺内的扩散。通过比较扩散特性时间（分子完成一次随机跳跃的用时）与松弛时间也可用于推断松弛与扩散的关系。关于低水分面粉或面团的β-松弛时间鲜见报道，采用水分质量分数5.4%淀粉的松弛时间τ约为5×10-8s[48]，假定扩散特征时间与松弛时间相等，结合本研究40 ℃低水分时的扩散系数约为8×10-14m2/s，计算水分子在面团膜内的跳跃步长约为0.15 nm。平移被冻结的葡萄糖玻璃体中水分子的跳跃步长约为0.35 nm[44]。与这些结果相比，本研究得到的扩散步长明显偏低。上面估算扩散特征时间时采用淀粉的β-松弛时间可能比面团膜体系的实际值低，因为面团膜体系中除淀粉外还含有蛋白质，含水量为6%的弹性蛋白在27 ℃时的β-松弛时间在10-6s数量级上[49]，比淀粉大。

4 结论

VD理论对面团膜在低水分区的模拟结果偏低，改进后的模型能够较好模拟面团膜的整个吸湿过程，反演求得的自由体积参数值基本合理。模拟分析显示，水分扩散系数随水分质量分数增大而增大，且变幅较大，但与以往的研究结果不同之处在于，这个变化过程分为3 个阶段：第1个阶段，扩散系数略有减小；在第2、3阶段，扩散系数增大，来自水分的自由体积增多是主要因素；对于第1、2阶段，玻璃化转变所挟裹的自由体积会部分补偿水分含量较低而损失的自由体积。对低水分段水分扩散机理的深入分析显示，面团膜在第1阶段还可能发生了β-转变，由此带来另外的自由体积及扩散系数的进一步增大。采用不同于以往的确定扩散系数与水分含量关系的计算策略，得以揭示出扩散系数随水分含量变化呈现出3 个阶段这一特征。