不可压缩圆柱绕流_坐标变换_SOR方法选编.docxVIP

不可压缩圆柱绕流摘要本文对无粘不可压圆柱绕流问题进行数值模拟，并将数值结果和理论解进行对比。在给定圆柱固壁和远场速度的边界条件下，通过求解流函数的Laplace方程，得到域内的流函数，进而计算出速度场。数值模拟的主要步骤有网格划分，物理域和计算域坐标变换，对流函数方程进行差分离散，SOR迭代方法求解和后处理绘出速度场。本文将先建立本问题的数学提法和理论解，随后基于数学提法进行数值求解，最后将理论和数值进行对比。理论和数值得到的速度场表现出了一致性，并且均反映出无粘假设下的流动和实际流动间的巨大差异。关键词不可压缩；圆柱绕流；坐标变换；SOR方法1 问题提出考虑图1所示无粘不可压圆柱绕流，来流速度为。选取合适的计算域范围和合理的边界条件，通过求解流函数方程或者速度势方程求圆柱周围的流场并与理论解相比较。图1无粘不可压圆柱绕流示意图2理论求解在圆柱圆心处取笛卡尔坐标，对于不可压缩二维流动，存在流函数有无穷远处来流，为无旋流动，根据Kelvin-Lagrange定理可知，本问题中，流场始终无旋。由，代入可得主控方程定解的边界条件有两个，一个是圆柱的固壁条件，考虑到固壁的不可穿透性有上式中为圆柱半径，第二个边界条件是无穷远处速度分布在圆柱坐标系中求解上述问题，坐标定义如图2，可得到坐标变换关系和微分关系图2坐标关系说明根据进行坐标转换，可得到主控方程和边界条件如下通过分离变量法[1]，可