CN112861313A - 一种粮食床通风过程中温度和水分变化的模拟分析方法 - Google Patents

Abstract

本发明公开了一种粮食床通风过程中温度和水分变化的模拟分析方法，所述方法包括：基于粮堆焓平衡方程和水分平衡方程，采用离散化方法将粮食床分成N个节段，采用粮食平衡水分方程计算第1个节段的粮层含水率，然后采用迎风差分格式将温度和水分分配到第2、3、4、······N‑1、N节点的过程。能够合理估计变化的粮食温度和含水率，通风空气的温度和湿度通常在通风过程连续变化，本发明的方法模拟通风粮堆的热量和水分迁移，给出热量和质量平衡方程的逐步推导，就一维系统假定均一流速的气流通过粮食床，热量也仅在气流方向流动，既可阐明解决问题的原理，又可保证软件程序中简单而整体的代数式。

Description

一种粮食床通风过程中温度和水分变化的模拟分析方法

技术领域

本发明涉及食品科学技术领域，具体涉及一种粮食床通风过程中温度和水分变化的模拟分析方法。

背景技术

目前，国家加大了对粮食储运领域应用基础研究，一些学者开始探索粮堆的温度和水分迁移过程，当前国内公开的试验结果主要是采用商业计算流体动力学(CFD)软件进行模拟。该软件包含了连续性方程、动量守恒方程、对流传热方程和水分迁移方程，其运行计算机语言环境是Fluent软件中的Fortran语言,该软件非常昂贵，一般科技工作者和粮库无法利用。

探索粮堆温度和水分迁移过程的算法，意义重大。粮堆生态系统中生物与非生物成分之间相互作用，是动态的，每个成分持续地影响其它成分。通风的作用就是运动合适质量的空气穿过粮堆以调节储粮存在的条件。粮食储藏作业的目标就是保持粮食的品质，减少昆虫、螨类、霉菌危害粮食的可能性。为了成功实现这个目标，就要调控储粮生态系统微生境，造成粮食籽粒内、粮堆(空间或籽粒之间)空隙内的物理和化学条件对有害生物是敌对的，而有益于粮食品质保持。在实践应用方面，国家粮食局2002年发布了《储粮机械通风技术规程》(LS/T 1202-2002)，提出了允许降温通风的温度条件和湿度条件。在这个储粮技术规程及我国粮食储藏学科的专著中，对通风过程缺乏描述。粮堆降温通风作业往往伴随降水，温度迁移前沿区带早于水分迁移前沿区带，国内对粮堆降温通风期间温度和水分迁移的前沿区带和尾巴区带缺乏研究。与国外粮食品种单一、储存时间短的储粮比较，考虑我国的粮仓仓型、储粮区域、粮食种类品种、储存时间、粮食含水率及杂质含量等因素到通风过程模拟中，有利于节约能耗，降低粮堆水分损失，达到优粮优储的目标。

粮堆物理特性的知识提供了粮堆通风作业的原理基础。粮堆的物理特性是内在联系的，并被其它的物理、生物因素所影响。粮堆颗粒间隙的条件包括空气的湿度和温度、有毒气体如熏蒸剂、接触式杀虫剂。颗粒间隙的条件强烈影响粮堆籽粒的条件。粮食颗粒之间的空气温度和湿度是最关键的非生物因子变量，通过通风作业过程调控。粮堆内湿度和干球温度对应于籽粒间隙空气的湿球温度，湿球温度与昆虫群体生长的速率直接相关。在粮堆中霉菌倾向于如何发展，也是粮堆温度和湿度的函数。

发明内容

为此，本发明提供一种粮食床通风过程中温度和水分变化的模拟分析方法。

为了实现上述目的，本发明提供如下技术方案：

本发明提供一种粮食床通风过程中温度和水分变化的模拟分析方法，所述方法包括：基于粮堆焓平衡方程和水分平衡方程，采用离散化方法将粮食床分成N个节段，采用粮食平衡水分方程计算第1个节段的粮层含水率，然后采用迎风差分格式将温度和水分分配到第2、3、4、······N-1、N节点的过程。

本发明的一个实施例中，所述离散化方法为：在通风气流方向将粮食床分为N-1层，粮食床的高度是L，每个层长度是Δx，则方程为

Δx＝L/(N-1)

所述第1个节段的粮层含水率的计算方程为：

其中，M是平衡含水率，小数干基表示，ERH是平衡相对湿度，小数表示；系数A、B、C与之对应的国产粮食种类的如表1所示。

本发明的一个实施例中，所述粮堆焓平衡方程为：

式中，M，粮食含水率，干基小数表示；T，粮堆温度，℃；w，粮粒间隙空气的湿度比率，kg水/kg空气；x,一维方向的距离，m；t，通风时间，hr；H_w，粮堆湿化积分热，J/kg；h_s，粮堆吸附等热，J/kg；h_v，粮堆游离水汽化潜热，J/kg；h_v＝h_s+H_w。

C_a，空气的比热，1.0J/(g·K)；C_σ，粮食的比热，1.3J/(g·K)；C₁，水的比热，4.18J/(g·K)；ρ_a，空气密度，取为1.2kg/m³；ρ_σ，粮食籽粒的密度，取为1300kg/m³；ρ₁，粮堆籽粒间隙水蒸气的密度，等于ρ_σ·(1-ε)，取为780kg/m³；ε是粮堆空隙率，取0.4；k_eff，通过粮食床的热传导率，取为0.15W/(m·℃)；f_a，通过单位床面积的通风气流速率，kg/(m²·s)。

本发明的一个实施例中，所述粮堆水分平衡方程为:

式中，μ_a是含有湿蒸气的干空气速率，m/s。

将f_a＝ρ_a·μ_a代入粮食水分平衡方程得：

式中，假定水蒸汽和空气具有相同的速率。

本发明的一个实施例中，所述粮食床水分平衡方程通过下述质量方程解释水分通过粮食床任何附加的运输可能来自水蒸气的扩散：

式中，D_eff，通过粮食床的水蒸汽扩散率，取为5×10^-6m²/h；。

本发明的一个实施例中，所述粮层水分平衡方程进一步处理为：

式中，

分别是第p+1，p时刻步态后第i节点的粮食含水率；

分别是第p时刻步态后第i+1，i，i-1节点的空气湿度比率；在零时刻步态后第i节点的值

(p＝0)是已知的，该水分平衡方程用于计算时刻Δt之后的

本发明的一个实施例中，所述迎风差分格式为：

式中，w_i，w_i-1分别是第i、i-1节点的粮食颗粒间隙空气的湿度比率。

本发明的一个实施例中，所述方法包括采用中心差分格式计算粮食床湿度比率、温度和水分速率变化的二阶导数。

本发明的一个实施例中，所述方法还包括通过算法中对高含水率玉米粮堆涉及储藏真菌引起的呼吸热。

本发明具有如下优点：

本发明的粮食床通风过程中温度和水分变化的模拟分析方法，能够合理估计变化的粮食温度和含水率，通风空气的温度和湿度通常在通风过程连续变化，本发明的方法模拟通风粮堆的热量和水分迁移，给出热量和质量平衡方程的逐步推导，就一维系统假定均一流速的气流通过粮食床，热量也仅在气流方向流动，既可阐明解决问题的原理，又可保证软件程序中简单而整体的代数式。

本发明的模拟分析方法试验证明：采用湿度比率0.013kg/kg和20℃空气降温通风35℃和含水率13.5％的稻谷粮堆，稻谷粮堆平均温度和含水率在通风150小时分别是23.57℃和12.42％，在通风200小时分别是21.09℃和12.25％，在通风300小时分别是20.03℃和12.49％。

采用湿度比率0.013kg/kg和20℃空气通风35℃和含水率12.0％的大豆粮堆，粮堆平均温度和含水率在通风150小时分别是23.22℃和11.30％，在通风200小时分别是20.89℃和11.20％，在通风300小时分别是20.02℃和11.48％。

附图说明

为了更清楚地说明本发明的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是示例性的，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图引伸获得其它的实施附图。

本说明书所绘示的结构、比例、大小等，均仅用以配合说明书所揭示的内容，以供熟悉此技术的人士了解与阅读，并非用以限定本发明可实施的限定条件，故不具技术上的实质意义，任何结构的修饰、比例关系的改变或大小的调整，在不影响本发明所能产生的功效及所能达成的目的下，均应仍落在本发明所揭示的技术内容得能涵盖的范围内。

图1为本发明粮食床的离散化方法的结构示意图；

图2为本发明的稻谷和大豆粮食床降温通风期间各节段温度变化；

图3为本发明的稻谷和大豆粮食床降温通风期间各节段含水率变化。

具体实施方式

以下由特定的具体实施例说明本发明的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例1、通风粮堆的焓平衡方程和水分平衡方程的建立

粮堆通风过程包含了水蒸气速率(μ₁)和含有湿蒸气的干空气速率(μ_a)。可认为这两个速率是一样的，由于干空气和水蒸气是混匀的；湿蒸气沿着浓度梯度流动是可能的，它较干空气流动快或慢，其中它是被运载的。认为这两个速率均等于干空气的速率，即μ₁＝μ_a。粮食颗粒空隙之间空气的湿度比率(w)被近似为：

式中ρ_a是空气的密度，ρ₁是粮堆籽粒间隙湿蒸气的密度。

通过一维粮堆具有横切面积1m²的干空气的质量流速f_a由以下表达式给出，f_a＝ρ_aμ_a (2)

式中，μ_a是干空气的速率，m/s。

就焓平衡方程，包含粮食温度(T)、粮食含水率(M)和粮粒间隙空气湿度比率(w)，表达为：

式中，M，粮食含水率，干基小数表示；T，粮堆温度，℃；w，粮粒间隙空气的湿度比率，kg水/kg空气；x,一维方向的距离，m；t，通风时间，h；H_w，粮堆湿化积分热，J/kg；h_s，粮堆吸附等热，J/kg；h_v，粮堆游离水汽化潜热，J/kg；h_v＝h_s+H_w。

C_a，空气的比热，1.0J/(g·K)；C_σ，粮食的比热，1.3J/(g·K)；C₁，水的比热，4.18J/(g·K)；ρ_a，空气密度，取为1.2kg/m³；ρ_σ，粮食籽粒的密度，取为1300kg/m³；ρ₁，粮堆籽粒间隙水蒸气的密度，等于ρ_σ·(1-ε)，取为780kg/m³；ε是粮堆空隙率，取0.4；k_eff，通过粮食床的热传导率，取为0.15W/(m·℃)；f_a，通过单位床面积的通风气流速率，kg/(m²·s)。

粮堆水分平衡方程为:

代入f_a＝ρ_a·μ_a到方程4得：

推导方程5，假定水蒸气和空气具有相同的速率。但是Thorpe在1995年指出，水分通过粮食床任何附加的运输可能来自水蒸汽的扩散，能够被加入一个额外的术语质量守恒方程来解释：

式中，D_eff，通过粮食床的水蒸汽扩散率，取为5×10^-6m²/h；在流动方向的水分扩散项是可忽略的，即使当f_a小到0.0001kg/(m²·s)。在没有气流的粮食床水分扩散项是很重要的。调控粮堆热量和质量迁移的偏微分方法，方程3和6是偶联的，一个偏微分方程的解将影响另一偏微分方程的解。而且，

项不是线性的，不可能获得密闭形式的解，因此，粮食温度和含水率直接表达为以沿着粮食床距离和通风开始时间的函数。

数学方法描述的难度是，通风空气的温度和湿度随着时间发生任意变化。获得通风粮堆温度和粮食含水率响应的方法是利用数值解分析，解主导的偏微分方程。我们采用Kreyszig在1998年提出的简化数值解方法，用于解方程3和6。

本实施例中，粮堆和粮食床为相同客体。本发明的方法在通风粮堆控制热和水分迁移的方程数值解，如图1所示，粮食床的离散化方法，在通风气流方向粮食床分为N-1粮层，即具有N个节段，床的高度是L(m)，因此每个粮层的长度是Δx(m)，则方程为：

节点被步长距离Δx分开，总计N个节段，在第1个节段，通风所用的气流的温度和湿度比率是已知的，它们对应进入粮堆的空气温度和湿度比率。假定空气和粮食之间质量和热平衡，方程8能够计算第1个节段的粮层含水率。

式中，A、B和C为系数。发明人发现，现有技术中系数A、B和C的值不能够覆盖国内的同一种类粮食品种多的问题，且所述A、B和C的值与所述粮堆粮食的种类、吸着类型(吸附、解吸或吸着平均)有关。为此，在本公开的一个优选示例中，发明人采用静态称重的平衡水分测定方法，选定国产的17个品种的稻谷、14个品种的小麦、16个品种的玉米和10个品种的大豆，在10～35℃的温度范围内和11％～96％平衡相对湿度范围内，经过大量的测试，测定其相应的平衡水分等温线；进一步采用方程8，通过非线性回归拟合获得系数A、B和C的值。并给出了国外部分粮食种类的A、B和C的值，以便用于进口粮食的粮堆水分迁移和温度迁移的分析。

表1粮食种类的方程(8)的系数

在通风开始点，粮食床内粮食含水率和温度的分布(它们在整个粮食床内是不均匀的)是已知的。因此，温度和含水率被分配到节点2、3、4、······N-1、N。下一步是计算每个节点粮食温度和粮食含水率随时间变化。为了完成这个目标，对方程3和6出现的一阶和二阶导数作简单近似处理。例如，

式中w_i、w_i-1分别是第i、i-1节点的粮食颗粒间隙空气的湿度比率。这种特别形式的求差，被Patankar在1980年称为迎风差分格式，由于它使得方程的解非常稳定。偏微分fw/fx已经被近似，利用了目标点和这个点的迎风点(即正对气流的方向)的w值。

在

的近似值中，在T-ΔT/2和T+ΔT/2的H_w值用于计算导数值。由于选择的点在T点之后与之前均等间隔排列，这个原理用于计算偏微分，被称为中心差分格式。支配水分扩散通过粮食床的二级导数采用中心差分格式能够得到可靠的近似值。在这个情况下，二级导数定义为速率变化的速率变化。当Δx→0时，

式中，

分别是节点i-1与i、i与i+1之间湿度比率的平均梯度。作出下列近似式：

把它们代入方程10，结果是：

通过依次类推，温度随迁移距离x变化的导数可以近似化为：

正在通风的粮食床内，粮食含水率和温度从点到点是连续变化的，可从方程3和方程6中可以捕捉到。按照上述处理，含水率和温度这些变量仅在节点计算。换句话说，正在通风粮食床的条件仅在不连续点计算，这个问题已经被离散化。当Δx变得越来越小，数值解通常成为连续方程3和6的更好近似值。

粮食的含水率和温度也随时间连续变化，在这个方法中，它们仅在一定的不连续的时间间隔被评价。当通风系统打开之后，如果粮食的含水率和温度在时间间隔Δt之后被计算，这些值在时刻Δt、2Δt、3Δt、···、(p-1)Δt、pΔt、(p+1)Δt···被更新。这里p表示第p时刻的更新。仅在一定时刻计算变量的这个思路，用于离散化方程3和6的时间导数，这样，

式中

是第p+1时刻步态后第i节点的粮食含水率；

是第p时刻步态后第i节点的粮食含水率。最初的粮食含水率定义为

温度随时间的微分以同样的方式处理，则为，

解方程3和6，通过把方程9、13、14、15代入它们，给出离散形式。

水分守恒方程6接着变换为：

解方程18最简单的方法，是明确地更新粮食含水率，这样，

在解程序的开始，值

(p＝0)是已知的，因此，方程19容易用于计算时刻Δt之后的

由于唯一的未知量

发生在方程的左边，方程右边的所有变量是已知的，未知数被清楚给出。

在内部节点的粮食温度在下一个时刻步态之后相似地计算，同时粮食含水率被更新，即，

式中，C_a，空气的比热，1.0J/(g·K)；C_σ，粮食的比热，1.3J/(g·K)；c₁是水的比热，4.18J/(g·K)；w＝ρ₁/ρ_a；ρ_σ是粮食的密度，取为1300kg/m³；ρ₁是粮堆籽粒间隙水蒸气的密度，等于ρ_σ(1-ε)，可取为780kg/m³；ρ_a是空气的密度，取为1.2kg/m³；ε是粮堆孔隙率，取为0.4。

这个解的方法编成Visual Basic程序给出。需要注意的是，编写的计算机程序和上述方程之间的主要差异是，程序中不包含

表达，一是由于它的值不能够准确计算，二是有学者在1972年指出它的值可以忽略。

实施例2、计算粮粒间隙空气的湿度比率H(i)及粮食含水率M(i)

(1)粮食水分汽化的潜热(h_v)

h_v＝2501.33-2.363T (22)

(2)粮食籽粒间隙空气的相对平衡湿度ERH、湿度比率H(i)及湿化积分热H_w：

P＝P_s·ERH (31)

式中M_db是含水率，小数干基表示。

(3)计算粮食含水率M(i)

实施例3、呼吸作用对通风粮食床的热和质量传递的影响

焓守恒方程为：

式中，S₁是呼吸粮食释放水分的速率，kg/(m³·s)；Q_r，粮食底物的氧化热释放的速率[kJ/(m³·s)]，定义为反应产物焓之和与反应剂之间的差异。ε_r是小数表示的粮堆孔隙率，取为0.4；ε_σ＝1-ε_r。在这个点上，引入有效热分散率K_eff代替了有效热传导率k_eff，这突出强调热能量也被水分分散所运输，这来自粮食籽粒之间空气的随机流动途径，这增大了粮食床的表观热传导率。Thorpe和Whitaker在1992年采用热分散率K_eff研究通风粮堆的热和质量传递。这个方程的一维形式是：

按照化学成分组成，S₄＝-1.66S₁ (38)

式中，S₄是氮气等非反应气体释放的速率，kg/(m³·s)。

水分守恒方程变成：

水分守恒方程的一维形式是：

式中，S₁是粮食呼吸释放的水分，kg/(m³·s)；v_a是通风气流的速率，m/s。

呼吸热通过对试验数据拟合数学方程获得真菌呼吸热的计算方法。呼吸作用被认为是碳水化合物完全燃烧形成CO₂和H₂O，同时释放热量。1kg粮食底物氧化释放15778kJ热，形成1.47kg二氧化碳和0.6kg水。已经定量分析了呼吸速率的影响因子。这里采用Thompson在1972年对玉米呼吸的研究，他测定了干物质丧失是时间依赖的，经过时刻t秒后，干物质分数丧失dm为：

dm＝8.83×10^-4{exp(1.667×10^-6t)-1}+2.833×10^-9t (41)

方程41用于15.5℃、含水率25％(湿基)、破碎率30％的脱粒玉米。由于粮食温度和含水率随粮仓内位点发生变化，需要将真实时刻t值作图变为生理时刻t_p。真实时刻与生理时刻的关系方程为：

式中，M_M和M_T是粮堆含水率和温度的修正因子。

(1)温度修正因子M_T

当T≤15.5℃和M≤19％时，M_T＝32.2exp(-0.1044T-1.856) (43)

当T＞15.5℃和19＜M＜28％时，

M_T＝32.2exp(-0.1044T-1.856)+[(M-19)/100]exp(0.0183T-0.2847)(44)

当T＞15.5℃和M＞28％时，

M_T＝32.2exp(-0.1044T-1.856)+0.09exp(0.0183T-0.2847) (45)

(2)水分修正因子M_M

式中M_DB是粮食干基含水率，％。

方程41到46用于计算生理时刻t_p后干物质丧失dm,容易确定干物质丧失的速率：

如图2所示，采用本发明的算法分析国产稻谷与大豆粮堆降温通风期间的温度和水分迁移变化，分析结果如图2和图3所示。

表2给出了国产稻谷和大豆粮堆降温通风期间，粮堆各节段降低到20℃的时数及对应含水率。如表3所示，在粮堆含水率保持不变的情况下降温通风，稻谷和大豆粮堆由35℃降低到20℃，分别需要123和126小时。

如图2和图3所示，在连续通风条件下，在降温通风前100小时，稻谷粮堆前4层已经降低到20℃，粮堆平均温度和含水率分别是27.08℃和12.84％。降温通风到150小时，稻谷粮堆前7层已经降低到20℃，粮堆平均温度和含水率分别是23.57℃和12.42％。降温通风到200小时，稻谷粮堆前11层已经降低到20℃，粮堆平均温度和含水率分别是21.09℃和12.25％。降温通风到250小时，稻谷粮堆前15层已经降低到20℃，粮堆平均温度和含水率分别是20.20℃和12.33％。降温通风到300小时，稻谷粮堆第21层已经降低到20℃，粮堆平均温度和含水率分别是20.03℃和12.49％。

表2国产粮食降温通风各层(节段)降低到20℃的时数及对应含水率

表3国产粮食降温通风期间风机运转时数预测

在连续通风条件下，在降温通风前100小时，大豆粮堆前4层已经降低到20℃，粮堆平均温度和含水率分别是26.81℃和11.64％。降温通风到150小时，大豆粮堆前8层已经降低到20℃，粮堆平均温度和含水率分别是23.22℃和11.30％。降温通风到200小时，大豆粮堆前11层已经降低到20℃，粮堆平均温度和含水率分别是20.89℃和11.20％。降温通风到250小时，大豆粮堆前15层已经降低到20℃，粮堆平均温度和含水率分别是20.14℃和11.31％。降温通风到300小时，大豆粮堆第21层已经降低到20℃，粮堆平均温度和含水率分别是20.02℃和11.48％。

本发明是基于CFD的粮食床通风过程中温度和水分变化的模拟分析方法，能够合理估计变化粮食温度和含水率，通风空气的温度和湿度通常在通风过程连续变化，本发明的方法模拟通风粮堆的热量和水分迁移，给出热量和质量平衡方程的逐步推导，就一维系统假定均一流速的气流通过粮食床，热量也仅在气流方向流动，既可阐明解决问题的原理，又可保证软件程序中简单而整体的代数式。

虽然，上文中已经用一般性说明及具体实施例对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

Claims (9)

1.一种粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，所述方法包括：基于粮堆焓平衡方程和水分平衡方程，采用离散化方法将粮食床分成N个节段，采用粮食平衡水分方程计算第1个节段的粮层含水率，然后采用迎风差分格式将温度和水分分配到第2、3、4、······N-1、N节点的过程。

2.如权利要求1所述的粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，

所述离散化方法为：在通风气流方向将粮食床分为N-1层，粮食床的高度是L，每个层长度是Δx，则方程为：

Δx＝L/(N-1)

所述第1个节段的粮层含水率的计算方程为：

式中，M是粮食平衡含水率，小数干基表示，ERH是平衡相对湿度，小数表示；T是粮堆温度，℃；系数A、B、C如表1所示；

表1

3.如权利要求1所述的粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，所述粮堆焓平衡方程为：

式中，M，粮食含水率，干基小数表示；T，粮堆温度，℃；w，粮粒间隙空气的湿度比率，kg水/kg空气；x,一维方向的距离，m；t，通风时间，h；H_w，粮堆湿化积分热，J/kg；h_s，粮堆吸附等热，J/kg；h_v，粮堆游离水汽化潜热，J/kg；h_v＝h_s+H_w；

C_a，空气的比热，1.0J/(g·K)；C_σ，粮食的比热，1.3J/(g·K)；C₁，水的比热，4.18J/(g·K)；ρ_a，空气密度，取为1.2kg/m³；ρ_σ，粮食籽粒的密度，取为1300kg/m³；ρ₁，粮堆籽粒间隙水蒸气的密度，等于ρ_σ(1-ε)，取为780kg/m³；ε是粮堆空隙率，取0.4；k_eff，通过粮食床的热传导率，取为0.15W/(m·℃)；f_a，通过单位床面积的通风气流速率，kg/(m²·s)。

4.如权利要求1所述的粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，

所述粮堆水分平衡方程为:

式中，μ_a是含有湿蒸气的干空气速率，m/s；

将f_a＝ρ_aμ_a代入粮食水分平衡方程得：

式中，假定水蒸气和空气具有相同的速率。

5.如权利要求4所述的粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，

所述粮堆水分平衡方程通过下述质量方程解释水分通过粮食床附加的运输来自水蒸气的扩散：

式中D_eff是通过粮食床的有效水分扩散系数，取为5×10^-6m²/h。

6.如权利要求5所述的粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，

所述粮堆水分平衡方程进一步处理为：

式中，

分别是第p+1，p时刻步态后第i节点的粮食含水率；

分别是第p时刻步态后第i+1，i，i-1节点的空气湿度比率；在零时刻步态后第i节点的值

(p＝0)是已知的，该水分平衡方程用于计算时刻Δt之后的

7.如权利要求1所述的粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，

所述迎风差分格式为：

式中，w_i、w_i-1分别是第i、i-1节点的粮食颗粒间隙空气的湿度比率。

8.如权利要求1所述的粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，

所述方法包括采用中心差分格式计算粮食床湿度比率、温度和水分速率变化的二阶导数。

9.如权利要求1所述的粮食床通风过程中温度和水分变化的模拟分析方法，其特征在于，

所述方法还包括通过算法中对高含水率玉米粮堆涉及储藏真菌代谢产生的呼吸热。

Images