Trans Pore的三维模型：模拟多孔介质的干燥在传热传质.doc

Trans Pore的三维模型：模拟多孔介质的干燥在传热传质 Patrick Perré，Ian W. Turnerb E.N.G.R.E.F；森林实验室，试验产品队伍，Girardet 14街,F-54042南希Cedex,法国；数学科学学院昆士兰科技大学邮政信箱 Q400布里斯班 澳大利亚东部城市 摘要 在运输方程支配干燥过程耦合和高度非线性的性质，强调了在这一领域的适用性数值模拟Trans Pore，这是个独特的三维干燥模型,能够处理多孔介质的传热传质。在目前的工作中使用的方程组，是为在宏观尺度模拟多孔介质干燥更全面的物理描述。 一些模拟结果表明这样的一个工具能够提供新的可能性，尤其是该面临的是数效应研究了各向同性介质（一种轻质混凝土立方体）和各向异性介质（木板各向异性率和几何因素均较大。这些模拟证明，只有三维空间都能够捕捉到了在高温干燥情况下中等宽度的效果。各向异性介质中的超压的效果而厚度和宽度热场之间的耦合和压力驱动流  1.引言 在这个实验中我们将讲述一个新的三维的干燥模型。同时我们将从三个宏观状态变量（一个是能量守恒方程，一个是水分保持，另一个是气相的保持）来展开说明这个干燥过程。在这个构想中一些重要的物理现象包括毛细管和对流的流体流动，约束和气体扩散以及对流气体迁移，所需的导热性能及热蒸发。最后，等温曲线，热量和质量交换条件，在充分的饱和条件下形成这个数值模型。 由此产生的代码很复杂、也很费时间，这与大型网络所遇到的三维空间网络中所消耗的CPU的时间不差上下。为了抵消这种结果，采用了最先进的数字技术，比如用来解决网络的非线性离散的准牛顿法 研发三维的干燥模型是根据以前一些专家有关二维模型的报告所得出的结论而引发的。首先，木材干燥应力在干燥过程中发展依靠强烈板的宽度，而纵向的方向是非常重要的代表传递现象全面的无线电频率或模拟微波干燥可进行绝对必要的高度耦合的非线性方程组的完整解决干燥模型似乎是一个新的文。油藏模拟和污染建模可以有关的结构化和非结构化网格实现3 - D模型参考多孔介质的物理行为通量和本构关系表达式一些新的3 - D模型具体计算细节了在多孔介质支配传热传质现象的宏观方程 能量守恒 气体守恒 根据达西定律可以得出液体与气体的相对速度 已知相潜力和χ是深度标边界的外部干燥的样品表面提出的条件被认为是如下形式 Jw 和Je分别记录了在边界上总水的通量和总焓在干燥表面的外部压力是固定在大气的值Dirichlet边界条件Pg=P∞形成一个适当的非线性方程为这个主要变量 在这期间每个外部边界在计算区内的控制卷的非线性迭代，而根据守恒定律，所以方程必须解决掉。 为了减少总体所需的计算时间，在这个模型中我们采用了对称平面。把这个多孔性样本考虑在内的话，只有1/8的样品被计算过时间，同时所有流动的液体、蒸汽、空气和热气的通量在这个平面上检测到的是0。起初这些多空介质是有一定量的水分和热的分布。在这个过程中大气压力值一直都保持不变。 2.3数值策略用于运输模式分辨率 为了说明离散过程，同时也为了弄清楚这一过程，已决定从解析下液体守恒定律入手。这整个过程是和其他两个守恒定律是相同的。整个过程是把平衡定律与一个典型的立体的控制卷结合一起的，如图1 所示 从文章中看在整合了计算机内的每一个控制卷后，一个非线性方程组的结果使本系统内每一个公式可以转化成离散的守恒方程，比如：（7） 在方程7中表示的是总液体守恒量：。代表的是该点周围邻近节点p设置[绝对磁导率；迁移率和扩散率]在流量表达式是系统功能状态的变量。由于方程式7它是有必要近似于每一个梯度术语。对于正交网格，一定量的φ梯度可以和这个公式一样简单。但是需要注意的是当网格不是正角正交的时候，这个公式需要修改下，在这种情况下一个渐变的切向分量的修正项必须添加 在这个公式中假设张量条件是对角线。如果该系统的各向异性放进笛卡尔坐标中的主惯性轴中，这个假设是合理的。这个复杂的配置已经在其他地方详细讨论过了，根据木板从树上切下了的位置可以来旋转材料轴。向量表示的是几何量，如例子： 一个有效计算每个控制体积面标准的方法是运用式子 （可以参考下图表1以及face ABCD)。 注意下像术语表示的是at the th CV face对部件液体的流动张量评估。 离散模拟(7)液体的守恒定律是一个复杂的非线性函数的状态变量在(Sw, Pw, T, , , )节点p和周围邻居节点在控制卷考虑的范围内。表达式中的系数是把运输系数的非线性和控制体积面的相关几何因数集合在一起的。其中运输系数的非线性涉及到节点P和邻居节点nib的状态函数。注意，平流基本运输系数和扩散条件要区别对待，否则的话，方案可能不是单一的。 对于每一个变异系数来说，系统状态变量是被认为有一个确切的定义。事实上，六个变