矩形突扩管道流动的数值模拟-计算数学专业毕业论文.docxVIP

矩形突拶蛰道流动蛉数壤搂毂 矩形突拶蛰道流动蛉数壤搂毂 燃要 P$65452 中文提要 本文对突扩镑道的浚动进行了数值研究。采用有限容积法在交错嘲掺上对 流动控制方程离散，然后采用SIMPLE算法求解。在处理对流项的离散时，提 出了一个改进黪簿分辨零缀合据式COPMD。荣躅GMRES算法求魍褰教瑟褥 的大型非对称稀疏代数方程组。比较了二维矩形突扩管道流动分别用FUD、 COPMD格式诗薜结果瓣差冥。最岳数慧模毅了三维矩形突扩繁遘兹滚动，褥 到了与实验可比的结果。 关键词：突扩管道流动，有限容积法，SIMPLE算法，商分辨率格式， (jMRES鳞法 作 者：周翔 撵罨老瑟：凌舀平 矩形突势管避滚动熬数擅搂撵 矩形突势管避滚动熬数擅搂撵 A堍嘲l Abstract The flow in a pipe with sudden expansion is numerically studied in this paper．The finite volume method iS used to discrete the flow controling equations on staggered gfid and the SIMPLE algorithm is adopted to solve the equations．An improved hi曲一resolution composite scheme COPMD is proposed to discreting the convective term。The GMRES method is used to solve the large non—symmetric sparse algebraic equations．The difference between the results，using the FUD and COPMD schemes respectively for the two—dimensional case，are given．Finally，the three-dimensional pipe flow is simulated，and the computed results are qualitatively agreed with the experiments。 Key words：pipe flow with sudden expansion，finite volume method，SIMPLE algorithm，high-resolution numerical scheme，GMRES method Written by Zhou Xiang Supervised by Prof．Ling Ouo。ping Il 矩形突扩管道流动的数值模拟 矩形突扩管道流动的数值模拟 §1引言 §1．引言 流动分离现蒙一值是流体力学最重要的研究课题之一。突扩矩形管道的流 动困其几何形状简单，常被逸作实验斌数值研究该现象的铡子。擎鞠盼数值番拜 究都只对二维流动的情形。随着计算枫和计算技术豹发展，现在正逐渐转向三 维流动的研究，例如研究侧壁对流动结构的影响【11等。本文拟在已有二维突扩 管道滤动计算豹基础土，进一步提高数值计葵躲精度、稳定燃积效率，并开餍 对三嘏流动的数值研究。 本文采瘸鸯限襄羧法在交锩网搀土瓣流动控漤方程凑数，然嚣采愆 SIMPLE算法求解方程。离散流动控制方程中各项时，扩散项一般采用中心靛 分格式，两蹲流颈的离散格式籍存多荦争选择。为使繇逸差努格式稽凌较高置具 有有界性，本文采用高分辨率组合格式。判断格式是否具有有界性，有【2】中的 对流有界性准则(cBe)与【3】中的拓宽的对流有界性准鲻(EcBC)两种翔据。本文 通过数值试验，在ECBC撼围之外给出了一个改进的离分辨率缀合格式 COPMD。在采用SIMPLE算法计算过程中，需要多次求解不同的大型代数方 程组。在对滚颂戆离教穆式采翅高分辨率组合掺式爱，掰要墩解方稷组豹系数 矩阵将变为至少九对角阵(二维)或十三对角阵(三维)。本文采用适合求解大型日} 对称耩菠霞数方程组懿GMRES雾法枝替传统豹基于TDMA算法豹交警方囱迭 代法求解所得的代数方程组。 本文第二节奔绍了二维突扩矩形餐遂流动静数学模壅；第三节蕊荦介绍了 求解流体流动问题的有限容积法和SIMPLE冀法；第四节引入了对流项离散的 高分辨率组合格式并键出