基于流型结构构造的混合过程强化设计方法.docxVIP

基于流型结构构造的混合过程强化设计方法 设备优化设计是实现化工过程的有效手段之一。传统的化工设备优化设计往往取决于设计师的经验。在优化方面，大多数优化方案是基于现有的结构提出的，然后根据实验或数值计算方法将不同的方案进行比较，以选择最合适的结构。由于不同的经验背景，不同的设计师会提出不同的优化方案。因此，这种传统的设备优化概念往往缺乏理论基础，无法实现一般功能。因此，有必要在化工设备上建立优化的一般方法，这是目前化工设备结构尚未完全解决的一个根本性问题。 针对上述问题,许多研究者进行了有益的尝试. 其中,从热力学角度出发进行传递过程的强化已经取得了广泛的应用.Bejan 基于前人所做的工作,本课题组提出基于特定流型结构构造的传递过程强化以及相关设备的优化设计思路.作为该方法的初步探索,笔者主要针对气体混合过程的强化对相关设备的优化设计进行了研究. 1 熵产优化方法 热力学第二定律揭示由于过程的不可逆性会造成有用功的损失,表现为系统的熵产.那么,从热力学角度上讲,降低系统熵产可以降低有用功的损失, 这是一种行之有效的提高热力学效率的方法.系统的熵产是与系统内部的流动方式以及热量传递过程密切相关,通过优化设备结构改变流动以及传热方式,降低系统的熵产,从而提高换热效果. 本文尝试将熵产寻优方法应用于混合过程优化中,以熵产为目标函数,在一定约束条件下推导出熵产极值的必要条件,并通过添加体积力来改变流型结构,使熵产达到极值. 1.1 热质传递的相互作用 不可压缩牛顿流体单位体积的熵产 式中:i、j为速度方向;下标“α ”为组分.式(1)中右边第1项是由于流动阻力产生的熵产;第2项和第3项分别是由于质量扩散和热量扩散产生的熵产;第4项是由于热质传递的相互作用引起的熵产;第5项是由于化学反应产生的熵产. 本文不涉及到化学反应,假设对流传质过程为等温过程,可以忽略传热过程造成的熵产.因此,式(1) 简化为 理想气体混合物中任意组分的化学势为 式中:μ p 根据Fick定律,传质通量表示为 由于质量扩散引起的熵产表示为 根据伯努利方程,对于流体系统机械能的损失, 或者说主要的流体压降是由流体之间的剪应力造成的,而这也会导致系统熵产,即流动阻力熵产. 1.2 气体混合物的确定 (1) 混合过程在常压下进行,气体混合物视为不可压缩理想气体混合物; (2) 忽略密度差引起的对流过程,气体混合物视为均质混合物; (3) 不考虑温度的影响. 基于以上假设,氨气/空气二元混合物的熵产为 稳态层流下,混合过程可以由连续性方程、组分输运方程以及运动方程描述,即 熵产寻优目标需通过改变混合器内部流场来实现,本文通过体积力F来改变流场,进而实现使熵产取极值的目标. 1.3 以扩散为中心的传递机制. 通过改变速度场而提高系统的传热传质效率是一种行之有效的方法 压降的约束表示为 混合过程的强化可以描述为:在约束条件下(连续性方程、组分输运方程以及一定的黏性耗散值)构造特定的流型结构,调整系统熵产.由扩散而造成的熵产能反映扩散进行的程度,因此在不同的传递机制下可以通过调整熵产达到强化传递的效果.在以扩散为主导的传递过程中,增大熵产意味着扩散过程得到强化;而在以对流为主导的传递过程中,降低熵产意味着对流过程促进传质从而抑制由浓度梯度引起的扩散.因此基于熵产分析的混合过程强化是典型的变分问题. 构造相应的拉格朗日乘子方程为 方程(13)对速度求变分,得 根据文献[20,23,28]报道,乘子C 对比方程(11)和(16),可得虚拟体积力表达式为 至此,得到了优化流场的控制方程组,包含连续性方程(9)、组分输运方程(10)、运动方程(11)、(17) 以及乘子方程(13).通过求解该方程组,可改变流型结构,得到熵产极值时的流场. 乘子C 2 混合过程的优化 2.1 混合过程的优化 为了节省计算资源,将管式混合器的某一中心截面作为计算模型.图1为气体混合器的二维物理模型,长度和高度分别为6,cm和1,cm.进口分为3部分,中间为氨气进口,上下分别为含氨气质量分数为1%的空气进口.为了避免进口和出口影响混合过程的强化,将混合器的中间区域作为优化区域,优化区域距离进口1,cm,总长度为2,cm. 混合过程操作温度恒定为293,K,操作压力为常压.混合物的密度恒定为进料组成下的平均密度,黏度为1.79×10 2.2 边境条件 入口边界条件:速度入口 出口边界条件:压力出口 壁面边界条件:采用无滑脱边壁条件 2.3 网格数目对混合均匀度的影响 首先根据物理模型建立计算模型.采用Fluent6.3.26进行数值计算,速度和压力离散方程组采用SIMPLE算法求解,离散过程的对流项和扩散项采用QUICK差分格式. 网格数对数值模拟至关重要,合理的网格数目能使模拟计算快速而准确.本文划分了数