计算双螺杆挤出机流体流动的共形区域变换-北京化工大学.PDF

第45卷 第6期 北京化工大学学报(自然科学版) Vol.45,No.6 2018年 Journal of Beijing University of Chemical Technology (Natural Science) 2018 计算双螺杆挤出机流体流动的共形区域变换 李奇强摇 黄晋阳* (北京化工大学 理学院,北京摇 100029) 摘摇 要:双螺杆挤出机机筒内部的流体流动区域截面是随时间变化的复杂三连通区域,采取共形映射方法将这样 的三连通区域映射成不随时间变化的圆界区域,进一步将流动区域变成一个不随时间变化的柱形区域,可大幅简 化双螺杆挤出机流体流动问题的计算。 选择以外边界为圆周、内边界为随时间旋转的两个椭圆周的平面区域上的 热传导问题为研究对象,用本文方法进行计算并比较了精确解和数值解,计算结果表明两解的绝对误差较小。 本 文方法可用于双螺杆挤出机流体流动问题的计算,并简化了计算过程。 关键词:双螺杆挤出机;三连通区域;热传导方程;共形映射 中图分类号:O241郾82;O29;TQ320郾5摇 摇 DOI:10.13543/ j.bhxbzr.2018.06.019 求解双螺杆挤出机流体流动问题的主要困难。 文献 引摇 言 [4 -6]研究了用区域变换方法求解单螺杆挤出机 双螺杆挤出机输送效率高而且稳定,对于物料 的流体流动问题,但其采用的方法只适用于双连通 的处理更均匀,是塑料等高聚物生产的重要设备,因 区域同心的情况,因此需要寻找更合适的共形映射 此双螺杆挤出机内部聚合物流体的流动性态和数值 方法来解决三连通区域的变换问题。 模拟是一个热门的研究课题。 多连通区域的共形映射的计算方法有多种[7-9] , [7-8] 双螺杆挤出机内部的流体流动满足的数学模型 通常是将区域映射到圆缝 。 本文采用文献[9]的 是三维的非牛顿流体动力学方程组,区域复杂且随时 方法将三连通区域映射到圆界区域,同时对其算法的 间变化,因此是求解难点,通常直接采用有限元[1-3] 约束条件进行修改,以保证同一杆长位置不同转角的 等方法或将数学模型降维的方法来进行流体流动的 截面均映射到同样的圆界区域,并更正了其中共轭调 数值求解,但缺点是难以精确有效。 在此流体区域上 和函数求解方程式的错误。 为了验证共形映射算法 求解双螺杆挤出机挤出问题的另一个重要途径是先 的可行性,本文以内部边界曲线是两个椭圆的双螺杆 将求解区域规则化,并在这个规则区域上进行流体计 挤出机流体流动区域为例求解该时变区域上的热传 算。 通常来说双螺杆挤出机流体流动区域的截面是 导方程的定解,并将数值解与精确解进行比较。 计算 一个三连通区域,有3个边界:外边界是一个类似于8 结果表明,两解的绝对误差较小,可以用于求解双螺 字形的外边界,2个内边界是两根螺杆的截面边界。 杆挤出机内部流体流动问题。 因此用合适的区域变换方法将螺槽区域变成规则化 1摇 共形映射的计算 的区域,将大大简化对双螺杆挤出机流体的计算。 将一个截面在一个转动周期中不同时刻的区域 1郾1摇 共形映射方法的理论公式 形状经过规则化变换变为相同的圆界区域,再将双 本文使用并改进文献[9]的方法。 首先考虑一 螺杆挤出机流体流动区域的截面转化成圆界区域 个多边相连的有界平面区域赘 ,其边界分量是光滑a