整粒蚕豆（Vicia faba L.）射频（RF）加热本征介电特性：多频—水分—温度域表征、混合模型校正与预测建模

《Journal of Stored Products Research》：Intrinsic dielectric properties of faba bean for radio-frequency heating applications: Characterization and modeling

　　

研究人员对整粒蚕豆在多种条件下进行了射频（RF）介电特性表征；采用混合模型由体相测量估算本征介电常数；并对经复折射率（CRI）估算的介电特性在频率、水分和温度维度进行建模。介电常数模型达到 R² = 0.9955 和交叉验证 R²（CV-R²） = 0.9943；扩展介电损耗模型达到 R² = 0.9834 和 CV-R² = 0.9803。

**论文解读**

Oladimeji Azeez、Suka Thangaraju和Oon-Doo Baik围绕整粒蚕豆（Vicia faba L.）射频（radio frequency, RF）加热所需的本征介电特性开展研究。介电加热包括RF（1–300 MHz）和微波（300 MHz–300 GHz），通过偶极极化和离子传导实现快速体积能量沉积，较传统大宗农产品加热具有更小热梯度；RF加热效能由复相对介电常数（ε* = ε′ ? jε″）控制，其中ε′为介电常数、ε″为介电损耗因子。水分含量影响最主导，温度促进分子迁移和离子传导，RF频段中离子传导主导损耗；颗粒体系中体积密度决定空气—固体体积分数，进而影响测得体相介电特性。已有豆类及农产品介电数据在实验域、样品形态和是否校正间隙空气方面不同；体相介电特性是籽粒—空气的体积平均响应，受堆积状态、测量配置和加工尺度影响，跨条件迁移性差。复折射率（complex refractive index, CRI）、Landau–Lifshitz–Looyenga（LLL）和Lichtenecker（LE）等混合模型可降低间隙空气和体积密度影响，但整粒蚕豆在多个工业RF频率、水分和温度下系统表征并比较混合校正的研究缺乏，这限制了RF采后处理设计、模拟与优化。为此，研究人员在10.34–20.14%（wet basis, w.b.）、20–100 °C、5–40.68 MHz范围内测定体相介电特性，应用并比较CRI、LLL、LE以估算密度校正的本征电特性，量化模型选择影响并传播测量不确定度，建立ε′和ε″随频率、水分和温度变化的半经验与经验模型，为RF加热模拟和工艺设计提供连续性质函数。

**关键技术方法（不超过250字）**：样品为Snowbird蚕豆，来自C-Merak（Tisdale, SK, Canada），初始10.34%（w.b.），蒸馏水调理至10.34–20.14%（w.b.），4 °C平衡72 h。水分用ASAE S352.2烘干法；真密度（true density, ρ_t）用氦比重瓶、体密度（bulk density, ρ_b）用量杯，孔隙率由密度计算。平行板夹具与LCR（电感—电容—电阻）仪在5–40.68 MHz六个工业频率、20–100 °C三重复测量。以CRI、LLL、LE混合模型校正，Monte Carlo传播不确定度；建立ε′对数频率模型、机理型和经验型ε″模型，采用分组交叉验证、bootstrap（自举法）等评价。

**研究结果**

**物理特性**：经调理蚕豆水分为10.34–20.14% w.b.，重复标准差0.07–0.27个百分点；ρ_b由0.801降至0.737 g cm^?3；ρ_t由1.438轻微降至1.384 g cm^?3，12.22% w.b.时先升至1.441 g cm^?3；孔隙率由0.443升至12.22% w.b.的0.463，14.12% w.b.后维持约0.467。间隙空气近半，说明需混合校正；孔隙变化小，校正值差异主要来自测得体相介电特性。

**体相介电特性与混合校正**：全范围内体相ε′为2.516–28.803、ε″为0.157–26.348。27.12 MHz、25 °C下，水分由10.34升至20.14% w.b.使ε′由2.785升至8.745、ε″由0.205升至1.200；20.14% w.b.升温至100 °C时ε′、ε″分别升至20.498和13.300；5→40.68 MHz使ε″由26.348降至11.345。模型估计总体由CRI向LLL、LE增大，但高损耗时顺序不完全一致；最高水分和温度下模型分歧增大，LE出现复数对数变换相关相位效应。CRI因文献基准和较低不确定性作为参考，LLL、LE用于展示模型选择敏感性；本征值为模型依赖估计，而非无空隙固相直接测量。

**Monte Carlo不确定度传播**：对体相ε′、ε″、ρ_b、ρ_t抽样并传播至各混合模型。ε′的95%相对半宽中位数为CRI 12.26%、LLL 13.53%、LE 21.83%；ε″为17.31%、17.47%、20.97%。LE最大ε′半宽140.58%，CRI和LLL分别约28.35%和31.51%。CRI与LLL传播相近而CRI略低；不确定性取决于模型公式。高温高湿下CRI估计值高且区间相对稳定，说明其高值不能仅归因于随机测量变异，也可能反映水合基质增强响应和校正放大。研究区分实验输入变异、混合模型依赖和未建模效应，指出窄区间不等于更真实，模型间差异也不等于体相测量不可靠。

**CRI估算的本征介电特性**：CRI校正后ε′为4.21–83.14、ε″为0.36–85.46；27.12 MHz下分别为4.55–57.43和0.49–42.27。最低频率和最严酷水—温组合下ε″可接近或略超ε′，损耗角正切约为1，符合电导主导损耗。

**频率依赖性**：ε′和ε″均随频率增加而降低，频率依赖性在高水分下增强；70 °C、20.14% w.b.时ε′由5 MHz的54.77降至40.68 MHz的37.99，ε″由39.85降至18.62。ε″逆频依赖更强，反映潮湿生物材料中电导主导损耗。

**水分和温度依赖性**：ε′和ε″总体随水分和温度增加并呈强耦合。27.12 MHz下，10.34% w.b.时ε′由25 °C的4.84升至100 °C的10.07，20.14% w.b.时由21.84升至57.43；ε″在16.36%和20.14% w.b.、60–70 °C以上增加更显著，5 MHz、25 °C下由10.34% w.b.的0.86升至20.14% w.b.的4.21。温度升高可能增强RF能量吸收，但温度发展还受场分布、几何与堆积、热湿传递及发生器—负载耦合影响，数据未证明自维持加热失稳。

**模型拟合与性能**

**介电常数模型**：采用ε′ = a + b log₁₀(f)形式，a、b为水分和温度的poly32响应面（含三阶项与交互项）。a、b面R²为0.992和0.974；全数据R² = 0.9955、RMSE = 0.8534、MAE = 0.6248；留一水—温条件交叉验证R² = 0.9943、RMSE = 0.9529、MAE = 0.6958；留一温度水平R² = 0.9948；留一水分水平降至R² = 0.9493、RMSE = 2.8513。系数需保留足够精度，四位小数会降低性能。模型适用于已表征域内插值，对未表征水分水平需谨慎。

**机理型介电损耗模型**：基于电导主导损耗，ε″ = K(f/f₀)ⁿ，K为水分和温度的poly32函数；优化得n = 0.3323（95% bootstrap区间0.292–0.367），残差平方和在n = 1方向上升最陡。全数据R² = 0.9730、RMSE = 1.885、MAE = 1.119；留一条件交叉验证R² = 0.9668、RMSE = 2.089。K(M,T)随水分和温度非线性上升，最高水平最陡；n偏离1表明理想逆频电导不能完全描述色散，异质水合基质中的偶极和Maxwell–Wagner型界面极化可能参与，但未量化。

**经验型介电损耗模型**：在log₁₀(ε″)域拟合。精简模型R² = 0.9704、RMSE = 1.973、交叉验证R² = 0.9619、RMSE = 2.239；扩展模型含高阶交互项，全数据R² = 0.9834、RMSE = 1.4795、MAE = 0.7216，留一条件交叉验证R² = 0.9803、RMSE = 1.6086，留一温度R² = 0.9756；赤池信息准则（AIC）和贝叶斯信息准则（BIC）最低。然其高阶对数频率或水分项在留一频率或水分时设计矩阵降秩，不能唯一识别，故不适合外推到未表征频率或水分。

**交叉验证与比较评价**：在同一分组方案下，扩展经验模型精度最高，但留一水分水平无法唯一估计；机理模型在留一水分水平下较稳健（R² = 0.7032、RMSE = 6.2496），精简经验模型为R² = ?0.0509、RMSE = 11.7605。留一频率时精简经验模型RMSE = 2.059优于机理模型2.237，但其他评价不优。故扩展模型仅推荐在六个实测频率上对水分和温度插值，机理模型推荐用于5–40.68 MHz RF加热模拟及需对未表征水分稳健的情形；所有模型不应外推到研究范围外而未独立验证。

**讨论总结**：CRI参考地位基于已发表基准及较低传播不确定性，但模型间在高湿高温条件下分歧显著，本征值是模型依赖估计；混合物校正可能提高跨堆积条件可迁移性，然数据未证明RF加热行为与堆积或尺度无关。高湿样品ε″随温度显著升高提示加热过程中RF能量吸收增强的潜在可能，但非热失控证据。文献比较受样品形态、温度、水分和堆积差异限制。局限包括：未用独立数据验证；室温密度一直用于100 °C；CRI未与蚕豆直接固相测量比对；Monte Carlo区间未含混合模型结构不确定性及水分和温度赋值不确定性，包括90–100 °C终点水分匹配。后续需独立验证、温度依赖密度测量和基于这些模型的电磁—传热耦合模拟。

**研究结论翻译**：据作者所知，该研究首次在多个工业RF频率（5–40.68 MHz）、水分含量（10.34–20.14% w.b.）和温度（20–100 °C）下，对堆积整粒蚕豆进行系统表征，并比较间隙空气的混合校正。体相和CRI估算的ε′、ε″均随频率升高而降低，总体随水分和温度升高而增加，二者同时较高时耦合效应最强；CRI估算ε′为4.21–83.14、ε″为0.36–85.46（27.12 MHz下分别为4.55–57.43和0.49–42.27）。CRI因已发表基准且传播不确定性最低（与LLL接近）而被采用为参考：ε′和ε″的中位相对95%半宽分别为12.3%和17.3%，LLL为13.5%和17.5%，LE为21.8%和21.0%。但在最湿、最热条件下模型估计仍有分歧，故本征值为模型依赖估计，而非无空隙籽粒的直接测量。对数频率ε′模型校准R² = 0.9955，留一条件交叉验证R² = 0.9943，留一水分水平验证降至0.9493，适用于已表征域内插值。对于ε″，扩展经验模型最准确（交叉验证R² = 0.9803），但在留一频率或水分水平时不能唯一识别；机理型幂律模型（n = 0.3323；95% bootstrap区间0.292–0.367）精度较低（交叉验证R² = 0.9668），但在留一水分水平时仍较稳健（R² = 0.7032，精简经验模型为?0.0509）。因此，扩展模型仅推荐在六个实测频率上对水分和温度插值，机理模型推荐用于5–40.68 MHz的RF加热模拟。最湿样品中ε″随温度明显升高，提示加热过程中RF能量吸收可能增强，但温度发展还取决于场分布、堆积及热湿传递。密度校正可能相对体相测量提高可迁移性，但这些数据未证明与堆积或尺度无关的RF加热行为。所有模型尚未用独立数据验证；室温密度应用于100 °C；CRI未与蚕豆直接固相测量比对；Monte Carlo区间未含混合模型结构不确定性及水分和温度赋值不确定性（包括90–100 °C终点水分匹配）。独立验证、温度依赖密度测量和基于这些模型的电磁—传热耦合模拟，是推进可靠蚕豆RF采后处理设计的下一步。