CFD 模拟(潘钰 2016,被引 20)+ 实仓试验(徐玉斌 2005)+ 原理分析(黄秋菊 2005),恰构成领域研究的方法论闭环,共同底座是 LS/T 1202-2002。
山东建筑大学王远成/杨开敏系(平房仓→浅圆仓径向通风+专著)、国家粮科院张忠杰/吴子丹/赵会义系(横向通风、180t 中试仓)、河南工业大学陈桂香/杨国峰系。
玉米横向粮层阻力≈竖向 1/2;竖向通风能耗约为横向 1/4;实仓调质增水 0.5~1.05 个百分点、能耗 0.407~2.04 kWh/(t·1%水分);装粮方式可使阻力差达一个数量级(国际实测)。
各向异性/变孔隙率未进 CFD 源项;均匀性评价指标不统一;调质"适度结露"合规争议与规程滞后(≤3 m、限加工前);国际大跨度平房仓风道 CFD 验证空白——中国途径比体系恰是优势。
基点论文 1 所引 Garg & Maier (2006) 经 OpenAlex/Crossref 多轮检索未检索到;"Smith 1996 TIPM"实为 de Ville & Smith 载于 Appl. Math. Modelling。新论文建议改引 Lawrence & Maier (2011/2012)。
三篇论文分别代表"数值模拟—实仓试验—原理分析"三种研究类型,合起来正好是储粮机械通风研究的完整方法论谱系。
方法:平房仓(60×21×10.5 m,粮堆 6 m,小麦 5800 t)U 型一机三道(途径比 1.21/1.56/1.83)与圭字型地上笼对比;粮堆作多孔介质,Ergun 动量源项;ICEM 约 310 万网格,SIMPLE 算法,单位通风量 5~25 m³/(t·h)。
结论:U 型阻力与途径比成正比、均匀性成反比;圭字型阻力介于 1.56~1.83、均匀性介于 1.21~1.56 之间(Weltens 指数,39 测点)。
方法:中储粮宁乡直属库 4550 t 晚籼稻,3 台湿膜保水机组(单台 20000 m³/h,饱和效率>90%),环流 243 h + 直流 30 h 交替。
结论:整仓增水约 1 个百分点,单位能耗 2.04 kWh/(t·1%水分),收益比 1:10.38;水分进入糙米而非稻壳;直流调质效率为环流的 3.76 倍,适当结露可提效——突破了当时规程"不允许结露、粮层<3 m"的限制。
框架:平衡水分原理——加湿条件为输入空气水蒸气分压>粮食平衡水分对应分压;增湿量按含湿量差计算;四因素:送风状态、风机选择、通风方式(循环 vs 直流,下进上回优先)、风道设计(途径比越大通风死角越大)。其"智能机械通风系统"展望已在 2010 年后成真。
方法主线:粮堆作多孔介质区域 + Ergun 型粘性/惯性阻力源项(ε=1 退化为 N-S 方程),风道与空腔区用标准 k-ε 湍流,粮堆内按层流处理;验证手段分实仓测点、中试模拟仓、试验台、工程能耗核算四类。
张来林七粮种回归 P=a·h·v^b(1993);杨进/杨国峰多元回归 P=a·h^b·v^c,深度指数 b>1——阻力随堆高超线性增长(1998/2001);Shedd 式、Hunter 二项式、规程公式四口径并存(王俊安 1999 对比)。已核实
国内粮堆温度场数值模拟早期代表作;粮科院自研 FUNLET 软件做平房仓压力场模拟。已核实快照
赵会义 180 t 模拟仓实测:玉米横向粮层阻力约为竖向的 1/2(中国粮油学报 2016(12))——粮堆阻力各向异性的直接证据,可用 Hunter 二项式或 Shedd 幂函数计算。已核实
横向通风无死角、装粮高度大时阻力优势明显;竖向通风能耗仅约横向的 1/4(俞晓静 2020,建筑节能);组合型风道气流最均匀;集风管与粮面平齐且上粗下细最优。团队出版专著《储粮生态系统多场耦合与数值模拟原理和方法》(科学出版社 2021)。已核实
仓储稻谷热湿耦合传递及黄变(品质)数值预测——南昌浅圆仓全年自然储藏+机械通风的温-湿-品质联合模拟,标志"温-湿-品质-虫害"多场耦合范式成型。李倩倩 2022:分区横向谷冷实仓能耗 0.1 kW·h/(t·℃)(国标上限的 20%)、失水率仅 0.1%。已核实
建立浅圆仓/筒仓非均匀孔隙率模型:温度预测与试验吻合优于匀质假设,并解释筒仓"冷皮热芯"成因。此前国内 CFD 普遍采用各向同性匀质 Ergun 源项。已核实
COMSOL 多物理场路线快速增多(温度年变化模拟最大误差 0.91 ℃);李杰:25 m 直径浅圆仓径向支风道 12.5 m 降温速率最大(0.051 ℃/h),地槽体系圭形综合最优;林俊 2024:分形通风系统 7.5 天最大温差 2.8→1.3 K、均匀性指数 0.919。已核实
在国内公开文献中未见应用报道——与国际方法存在代差;CFD 与 AI 结合仅见 YOLOv5s 通风特征识别(2023)、LSTM-MPC 干燥控制(2023)等萌芽。未检索到
| 主题 | 代表文献 | 关键量化结果 |
|---|---|---|
| 途径比 vs 阻力/均匀性 | 潘钰等 2016(基点论文) | U 型:阻力∝K,均匀性∝1/K;圭字型介于 K=1.56~1.83(阻力)、1.21~1.56(均匀性)之间 |
| 横向 vs 竖向阻力 | 赵会义等 2015/2016/2019 | 玉米横向粮层阻力约为竖向 1/2(180 t 模拟仓)——各向异性 |
| 通风方式能耗 | 俞晓静 2020;杨泰 2022;李倩倩 2022 | 竖向能耗≈横向 1/4;联合通风能耗最小;分区横向谷冷 0.1 kW·h/(t·℃)、失水 0.1% |
| 系统总阻力 | 马春云、曹阳 2017 | 27 m 跨度高大平房仓横向通风总阻力 <1400 Pa |
| 单位通风量 vs 降温 | 陈桂香 2014;张淑媛 2024 | 降温时间与单位通风量近似反比(非严格线性);3→8 m³/(t·h) 显著加快;60×24 m 仓 96 h 通风 36 h 后大部降至 15 ℃ |
| 非均匀孔隙率 | 柳阔 2023;张修霖 2021 | 温度预测精度优于匀质模型;解释"冷皮热芯" |
| 均匀性量化 | 林俊 2024 | 分形风道最大温差 2.8 K→1.3 K,均匀性指数 0.919 |
判据核心是"大气绝对湿度 vs 粮堆平衡绝对湿度":LS/T 1202-2002 规定调质通风湿度条件为 PS2.5<PS1(按水分增 2.5 个百分点后的平衡绝对湿度<即时大气绝对湿度)——2.5 个百分点裕度即"适度饱和增湿"的工程预留。李兴军等进一步提出"调质通风窗口"自控判据(大气与粮堆平衡相对湿度差 ≥25% 时开机)。
"粮食施药调质机"改造、离心风机压送高湿空气;徐玉斌 2005 湿膜保水机组实仓完成装备化验证。
湿膜(等焓加湿、无过饱和)、超声波/超声雾化(下行式内循环,蒋社才 2011);"内结露"路线首创(金梅/陶诚 2010,史钢强 2013/2014 发展):不外加湿源,利用东北粮堆"冷心"冷源使循环空气在粮堆内适度结露放湿。
史钢强 2014 绝对湿度判据程序化;"十三五"重点研发课题 2016YFD0401002 研发稻谷横向通风保水调控工艺与智能控制系统,通风水分损失控制在 0.3% 以内(2021 年通过绩效评价);山东建筑大学毛煜东团队承担数值模拟子课题(75 万元);刘家琦 2023 硕论做玉米调质热湿迁移预测。
| 品种/规模 | 设备与风路 | 增水幅度 | 要点 |
|---|---|---|---|
| 稻谷 4550 t(2005) | 湿膜保水机组,环流+直流交替 | ≈1.0 个百分点 | 直流调质效率为环流 3.76 倍;收益比 1:10.38 |
| 粳稻 2489 t(2017) | 2×0.85 kW 轴流吸出式,仅自然空气 | +0.5 | 333 h / 566 kWh;上下 4 层水分均匀;整精米率 71.7%(+3%) |
| 玉米(哈尔滨 2013) | 轴流风机+"冷心"内结露 | 13.5%→14.1% | 全仓增重 21.05 t;但"内结露严重、易局部水分偏高,且违反操作规程" |
| 玉米整仓(2006) | 湿膜+智能控制+离心风机下行循环 | 12.9%→13.6% | 多源题录一致 |
| 稻谷 1380 t(中山,出库前) | 智能超声雾化,下行式内循环 | 11.16%→11.74% | 轮换出库减损场景 |
| 加湿调质原理试验(齐河库 2010) | 加湿通风 | 水分自下而上推进 | 约 4 d 达层间均衡(张传洪、王远成等) |
国际文献经 OpenAlex / Crossref / Semantic Scholar / ASABE elibrary 学术 API 逐一核实,标注"摘要已验证"者均读取到含数值的摘要原文。
| 文献 | 关键量化结果 | 验证 |
|---|---|---|
| Molenda et al. 2005, Trans. ASAE 48(3) — Ergun 方程适用性 | 玉米/小麦/大豆实测:压降<500 Pa/m 时 RMSE<23 Pa/m,至 1800 Pa/m 时 RMSE=76 Pa/m;装粮方式致容重差异使阻力变化达一个数量级 | 摘要已验证 |
| Kumar & Muir 1986, Trans. ASAE — 各向异性 | 0.077 m/s 下垂直阻力比水平高 60%(小麦)~115%(大麦);分层装粮比分点装粮高 25~75% | 摘要已验证 |
| Ramaj et al. 2022, Appl. Sci. — 自压实 | 孔隙率不可逆演化 0.413→0.391(1 h)→0.370(236 h);Matthies–Petersen 模型 R²=0.990 | 摘要已验证 |
| Neethirajan et al. 2006, Biosyst. Eng. — X-ray CT 孔隙结构 | 水平气路面积/长度比垂直大 100%(麦/大麦/亚麻);气路数量多 92~187%;"孔隙面积/气路数之比"是方向性阻力最佳预测量 | 摘要已验证 |
| Lawrence & Maier 2012, Appl. Eng. Agric. — 粮面构型 | 40→20 ℃ 用时:挖核 102 h、平顶 114 h、尖顶 186 h(延迟 84 h≈55%,直接换算风机时数与电耗) | 摘要已验证 |
| Khankari, Morey & Patankar 1995 — 水分迁移/结露风险 | 秋季将粮温冷至 0 ℃ 可使全年水分稳定;小仓水分迁移更早发生 | 摘要已验证 |
| Morales-Quiros et al. 2019 — 冷冻通风商用对照(2×1350 t) | 28→17 ℃ 用 175 h;39→17 ℃ 用 245 h;环境通风夏季降不到 25 ℃、失水 0.5%;害虫繁殖率显著受抑 | 摘要已验证 |
| Plumier & Maier 2021, Agriculture — 测点布局优化 | 测温电缆均值 vs 全仓均值差最高 11.4 ℃;不同布局作控制输入的决策差 3.6 ℃ | 摘要已验证 |
| Khatchatourian 系(巴西,2006/2008/2009)+ Binelo 2019 | 大型粮堆气流-热状态耦合、三维气流、各向异性渗透率;以入口压力廓线为设计变量的气流优化 | 元数据验证 |
| Yue & Zhang 2017;Nwaizu & Zhang 2021 — 孔隙尺度模拟 | 绕开 Ergun 宏观经验式,基于真实孔隙重构直接模拟气流分布 | 元数据验证 |
| Kaliyan et al. 2006/2007;Bareil 2018 — 通风×生物学终点 | 累积致死指数 CLI=1 害虫判据(5×76.2 t 仓验证);冷却-害虫致死联合预测 | 摘要已验证 |
这是报告的落点:若研究延续平房仓"途径比—阻力—均匀性"或调质通风方向,以下缺口可直接支撑论文创新点陈述。
赵会义 180 t 实测(横向≈竖向 1/2)、Kumar & Muir(垂直高 60~115%)、Neethirajan CT 孔隙级证据、Ramaj 自压实演化已充分证明粮堆阻力各向异性且随时间变化,但国内多数论文仍用各向同性匀质 Ergun 源项,仅柳阔 2023、张修霖 2021 起步。把"变孔隙率+各向异性渗透率"引入风道布置—均匀性研究,创新性直接成立。
变异系数、Weltens 指数、最大温差、降温速率、标准差拟合并存,横向不可比。可提出统一评价框架,并与 GB 50320 的 K=1.2~1.5 与山西口径(降温 1.5~1.8 / 降水 1.2~1.5)的取值分歧做规范符合性对照。
"内结露"效率高但按现行规程属违规操作;规程"加工前+≤3 m"与实仓 4~6 m、出库调质实践脱节;局部水分偏高幅度—霉变概率—安全水分裕量的定量关系无人给出。"适度结露的风险定量边界 + 出库调质验收指标"是最有价值的立论缺口。
湿膜/超声雾化/高压微雾/内结露/蒸发冷却各路线无按统一方法测试的加湿量、饱和效率、吨水电耗对比数据;单位通风量对调质效果的影响在规程中阙如。
多数验证以温度为主、测点稀疏,水分场三维全场验证缺乏;多孔跳跃边界、局部非热平衡模型未见应用。国际已有测点布局优化与模型降阶控制,国内"CFD 代理模型+通风优化"空缺,AI 多停留在视频识别与 LSTM 预测。
全部经实际检索获得;快照级条目建议正式引用前下载原文核对数字。
官方/已核实层:LS/T 1202-2002 条文(河北省粮局官网 PDF)、GB 50320-2014 第 3.3 节条文(规范库全文)、三篇基点论文原文、山西粮局口径(官网)、长治/安阳案例(省级官媒/政府门户)、国际文献中标注"摘要已验证"的条目(OpenAlex/Crossref/Semantic Scholar/ASABE elibrary)。
题录/快照层:维普检索页题录(作者/年卷期/被引可靠,数值结论来自摘要);搜索快照与文库转载(题录可靠,具体数字未经全文复核——尤其段义三 0.455、德州 0.407 两组能耗,正式引用前请下载原文核对)。
未检索到(勿引用):Garg & Maier (2006);Smith 1996 年 TIPM 气流论文原文;GB 50320-2024 条文;LS/T 1202 允许通风条件中具体温差数值原文;国内文献中粮层阻力公式系数 a、b 具体数值(公开摘要均未给出);调质装备统一测试对比数据。
经验性判断:"能耗排序自然大气<内结露<外加湿"由 7 组实仓数据归纳,样本有限;"研究空白"判断基于本次检索覆盖面,不排除未收录的内部资料与学位论文。