trans pore的三维模型模拟多孔介质的干燥在传热传质.docVIP

Trans Pore的三维模型：模拟多孔介质的干燥在传热传质 Patrick Perré，Ian W. Turner b E.N.G.R.E.F；森林实验室，试验产品队伍，Girardet 14街,F-54042南希Cedex,法国；数学科学学院，昆士兰科技大学；邮政信箱 2434 Q400布里斯班 澳大利亚东部城市 摘要 在运输方程支配干燥过程耦合和高度非线性的性质，强调了在这一领域的适用性数值模拟。这项工作主要是介绍关于最新版本的Trans Pore，这是个独特的三维干燥模型,能够处理多孔介质的传热传质。在目前的工作中使用的方程组，是为在宏观尺度模拟多孔介质干燥更全面的物理描述。 一些模拟结果表明这样的一个工具能够提供新的可能性，尤其是该研究面临的是数效应研究了各向同性介质（一种轻质混凝土立方体）和各向异性介质（木板，其各向异性率和几何因素均较大）。这些模拟测试证明，只有三维空间都能够捕捉到了在高温干燥情况下中等宽度的效果。确实，在这种条件下，纵向方向需要突出一个强烈的各向异性介质中的超压的效果，而相应的厚度和宽度是热场之间的耦合和压力驱动流所不可缺的。 1.引言 在这个实验中我们将讲述一个新的三维的干燥模型。同时我们将从三个宏观状态变量（一个是能量守恒方程，一个是水分保持，另一个是气相的保持）来展开说明这个干燥过程。在这个构想中一些重要的物理现象包括毛细管和对流的流体流动，约束和气体扩散以及对流气体迁移，所需的导热性能及热蒸发。最后，等温曲线，热量和质量交换条件，在充分的饱和条件下形成这个数值模型。 由此产生的代码很复杂、也很费时间，这与大型网络所遇到的三维空间网络中所消耗的CPU的时间不差上下。为了抵消这种结果，采用了最先进的数字技术，比如用来解决网络的非线性离散的准牛顿法以及用于解决稀疏矩阵系统的线性化的预条件共轭梯度。 迄今为止，已有一些综合二维的模型可以用来研究干燥工艺。这些模型已被应用于木材烘干领域和建筑材料，比如砖。这些模型为干燥制度提供了一些基本的知识不仅使人们更了解这一程序同时，提出该更多具有创新的新的建议。 研发三维的干燥模型是根据以前一些专家有关二维模型的报告所得出的结论而引发的。首先，木材干燥应力在干燥过程中发展依靠强烈板的宽度，而纵向的方向是非常重要的代表传递现象。其次，在这种模式之前保证全面的无线电频率或模拟微波干燥可进行是绝对必要的。 三维空间的模型在多孔介质中的使用并不陌生，但是据我们所知，综合性很强的计算模型目前在干燥研究这块领域还不存在，不过存在一些简单的用于研究木梁的吸收和解吸水的性能以及分析木料木块的结构的传输模型。高度耦合的非线性方程组的完整解决干燥模型似乎是一个新的有价值的文献。在地下水流领域，油藏模拟和污染建模可以提供有关的结构化和非结构化网格实现3 - D模型参考资源。虽然在理论上，从二维到三维是很顺其自然一气呵成的，但从以上所有研究中我们可以看出要实现模拟时间的高效计算，线性代数中一些微妙的数据分析和创新性概念是必不可少的。 在这项研究中已经有两项研究来检验三维模型的性能。在实验1中描述的是在高温传递下轻质混凝土的烘干；实验2 针对高温与低温两种情况所产生的结果进行了分析。尤其需要强调的是在整个烘干过程中只有三维才能够正确地捕捉到多孔介质的物理行为。的确，高温干燥超压在一个重要的板材纵向的一个板盖上有高度各向异性的多孔介质。厚度和宽度的板必须解释热耦合机制、质量扩散和达西流。 2.数学公式 这个部分将着重介绍下干燥守恒方程以及相关通量和本构关系表达式。解析下离散化过程、简要概述下数值求解过程。还有讨论一些新的3 - D模型具体计算细节。 2.1传输模型 Whitaker详细研究了在多孔介质支配传热传质现象的宏观方程。这些公式后来也相应的被应用于软木干燥过程的模拟中。在那里关于它们有呈现几种不同干燥配置的能力也已得倒了证实。这些方程如下： 液体守恒 能量守恒 气体守恒 根据达西定律可以得出液体与气体的相对速度 已知相潜力和χ是深度标就可以得出所求的数。其他符号还是原来的含义。 2.2.边界和初始条件 有两种边界条件需要考虑，一个是外部条件，另一个是对称面。边界的外部干燥的样品表面提出的条件被认为是如下形式： Jw 和Je分别记录了在边界上总水的通量和总焓的通量。在干燥表面的外部压力是固定在大气的值。在计算时起重要作用的其中一个变量是大气的平均密度，Dirichlet边界条件Pg=P∞不得不修改了，以便形成一个适当的非线性方程为这个主要变量： 在这期间每个外部边界在计算区内的控制卷的非线性迭代，而根据守恒定律，所以方程必须解决掉。 为了减少总体所需的计算时间，在这个模型中我们采用了对称平面。把这个多孔性样本考虑在内的话，只有1/8的样品被计算过时间，同时所有流动