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

PAGE \* MERGEFORMAT11  MACROBUTTON MTEditEquationSection2 Equation Chapter 1 Section 1 SEQ MTEqn \r \h \* MERGEFORMAT  SEQ MTSec \r 1 \h \* MERGEFORMAT  SEQ MTChap \r 1 \h \* MERGEFORMAT 不可压缩圆柱绕流 摘 要 本文对无粘不可压圆柱绕流问题进行数值模拟，并将数值结果和理论解进行对比。在给定圆柱固壁和远场速度的边界条件下，通过求解流函数的Laplace方程，得到域内的流函数，进而计算出速度场。数值模拟的主要步骤有网格划分，物理域和计算域坐标变换，对流函数方程进行差分离散， SOR迭代方法求解和后处理绘出速度场。本文将先建立本问题的数学提法和理论解，随后基于数学提法进行数值求解，最后将理论和数值进行对比。理论和数值得到的速度场表现出了一致性，并且均反映出无粘假设下的流动和实际流动间的巨大差异。 关键词 不可压缩；圆柱绕流；坐标变换；SOR方法 1 问题提出 考虑 REF _Ref376806707 \h 图 1所示无粘不可压圆柱绕流，来流速度为u∞=20m/s,v∞=0。选取合适的计算域范围和合理的边界条件，通过求解流函数方程或者速度势方程求圆柱周围的流场并与理论解相比较。 图  SEQ 图 \* ARABIC 1 无粘不可压圆柱绕流示意图 2 理论求解 在圆柱圆心处取笛卡尔坐标，对于不可压缩二维流动，存在流函数有  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) 无穷远处来流u∞=20m/s，为无旋流动，根据Kelvin-Lagrange定理可知，本问题中，流场始终无旋。由，代入 GOTOBUTTON ZEqnNum783236 \* MERGEFORMAT  REF ZEqnNum783236 \* Charformat \! \* MERGEFORMAT (1)可得主控方程  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 2) 定解的边界条件有两个，一个是圆柱的固壁条件，考虑到固壁的不可穿透性有  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 3) 上式中为圆柱半径，??二个边界条件是无穷远处速度分布  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 4) 在圆柱坐标系中求解上述问题，坐标定义如 REF _Ref405467412 \h 图 2，可得到坐标变换关系和微分关系  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 5)  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 6) 图  SEQ 图 \* ARABIC 2 坐标关系说明 根据 GOTOBUTTON ZEqnNum360833 \* MERGEFORMAT  REF ZEqnNum360833 \* Charformat \! \* MERGEFORMAT (5) GOTOBUTTON ZEqnNum998870 \* MERGEFORMAT  REF ZEqnNum998870 \* Charformat \! \* MERGEFORMAT (6)进行坐标转换，可得到主控方程和边界条件如下  MACROBUTTON MTPlaceRef \* MERG