液态流化床传热模拟.doc

液态流化床传热模拟 摘要 在这中,的热联合连续与离散模型。床的模拟来流动和传热行为。两个对流换热模型结果(阿加沃,1988,化学工程科学,43期, 2501页，及李和梅森,2000,粉末技术112期，273页) 已使用辛普森和帕特的实验经过了严格的评估(1977,化学工程科学,32期,67页) 使用丽和梅森的对流换热模型最大传热系数预测与辛普森和帕特实验良好吻合。报告还发现, 床高于0.7有正常粒子速度 ,正常粒子速度对应于最大的传热系数。模拟的系统运动与混合表明,他们是不会受到的热行为。 进入、传出和穿过浸在液体固体微粒是工业过程中经常遇到的一个。传热典型的系统包括在食品加工混合物管道运输储罐放射性废料安全处理和加工利用三相床作为生物反应器此外, 在换热器管固体颗粒的存在可以是有益的,粒子可以提高传热系数减少的形成。虽然许多以往的研究都考虑了性质,他们的(例如作用在微观粒子力、粒子之间传热和粒子间流体传热),因为调查实验系统的多尺度复杂性和测量技术的局限性这个困难通过计算机模拟所谓连续离散模型在牛顿第二运动定律解决了离散粒子的运动,方程连续流体。 2 热连续离散联合模型 2.1标准CCDM公式? CCDM结合了预测粒子运动的离散元法（DEM）和模拟连续流体行为的计算流体动力学（CFD）。离散元法模拟是以牛顿运动第二定律为基础，采用数值时间推进计划，来评估单个粒子的轨迹和旋转。在每个时间步长使用适当的接触定律来计算接触力，并分解成法向和切向分量。DEM中的关键假设是，假设用一个足够小的时间步长，干扰的传播不能超过其近邻粒子。 对于连续流体流动的计算，采用简易方法来确定流体的流速和压力，从而求解当地平均的连续性和Navier - Stokes方程。通过对流体与粒子间相互作用力详细运用牛顿第三运动定律，这种对流体的计算流体力学的计算与单个粒子的离散元法解决方案相结合。这确保了两套不同的尺度上求解的方程得以耦合。原连续离散联合模型模型公式的细节，正如气体粒子系统模拟的发展那样，在[7,22,23]中给出。 在这项工作中，我们考虑两个与标准连续离散联合模型公式的重大偏离。首先，我们在这个模型中包括了一些更复杂的水动力作用，该水动力作用在液体粒子系统中起作用，但通常由于低流体粘度和密度在气体粒子系统被忽视。其次，在模型中包括了用于模拟传热的机制。下面这两套发展进行讨论。 2.2 液体粒子系统水动力作用 CCDM技术液固系统的应用更复杂流体粒子相互作用力的考虑有必要，在液体介质中这可以包括除了稳定的阻力附加质量，巴塞历史，剪升，旋转，压力梯度力，在气固系统占主导地位。[25]中所讨论的，在这个中，我们忽略了巴塞历史和剪切升力。这导致下面的表达式： 其中ρf为流体密度，di是粒子直径，εi是周围粒子的空隙率，t为时间，g为重力加速度，up,i uf,i分别是粒子的速度和粒子流体速度。表1总结了用于方程（1）的系数公式。 在液体粒子系统其他复杂的是，由于粒子表面间的粒间的碰撞气体粒子系统不同，流体的密度和粘度。液，每个粒子的，eliquid，被视为一个粒子数，[26]所给关系式： 其中eair是颗粒在空气中的系数，S是所谓的临界影响数，低于该临界数反弹不会发生。在这项中，基Gondret等人在一个广泛的流体和颗粒性质测量设置S等于10。pp,i是粒子密度2.3 流体能量和动量方程 为了考虑了液固之间的，首先当地平均流体的能量方程解决的相的温度场。在这项中，能量方程流体温度。假设流体的密度ρf，热容量C和热导率κf，是恒定的，类似于流体动量方程当地平均流体能量方程如下形式： 其中Q是体积流体粒子热源项。 流体密度的假设是特别重要，因为它允许能量方程的解和连续性方程解的。但是，这个假设是一个过于简单化的东西在许多系统流体的密度温度的函数变化。因此当地流体温度的差异会流体密度的差异，导致流体浮力和浮力驱动对流。 为了说明浮力驱动的对流热CCDM的影响，使用Boussinesq近似。这是假设流体密度的变化足够小流体方程的密度可视为常数。唯一例外的是动量方程，假设流体密度随温度线性变化。这种近似被认为是足以使许多传热和流体动力学问题准确，温度差异不是特别大。Boussinesq近似应用当地平均流体动量方程结果： 其中p和F的分别是，流体压力和流体粒子相互作用的体积力，τ是流体粘性应力张量，ρ为在参考温度T的流体密度，β是热膨胀系数。流体连续性方程不变，计算公式如下： 2.4.1 相接触的固体粒子之间的导热qc ij 热传导是热能通过介质从温度较高的区域向温度较低区域的自发转移。有许多有关运用离散元法求解接触粒子之间的热传导的研究，以及采用理论分析求解粒子团导热系数的推导。但是，这些作品，只涉及到光滑粒子之间的接触。应当指出，有大量的有关粗糙表面接