小攻角下翼型分離流結構的數值研究

     計算域的邊界包括翼型固壁、速度入口、自由出流出口、Interior以及固壁邊界。計算域的外邊界長度分別是翼型表面弦長和最大厚度的10倍距離。所采用的邊界條件為:固壁表面采用無滑移條件、進口、出口和Interior面由特征相容條件所確定、計算邊界由插值確定。
     計算區域離散為非結構網格，計算區域整體為類O形網格。網格數為95202 ,翼型上、下表面均為196點，第一條網格線距翼型表面最近距離為弦長的10－4倍。計算過程中，采用了區域加密的方法，對翼型周圍流場網格進一步加密。不同的來流攻角采用不同網格，以保證來流方向不變。

                  圖1 NACA4415翼型 攻角時局部網格
    計算基于RANS方法，采用二維穩態分離解法的隱式解法，空間離散格式采用二階迎風格式，壓力－速度耦合采用SIMPLE解法。
    3.計算結果與分析
    3.1升力系數、阻力系數及分離流場
    對NACA4415翼型在小攻角-4°～20°范圍下的分離流動結構進行了數值模擬，取來流速度和翼型弦長為特征量,遠前方來流馬赫數和雷諾數分別為Ma=0.024和Re=2.86×104。圖2為NACA4415翼型攻角-4°～20°的之間升力和阻力系數計算值。

攻角 /(°)圖2 在小攻角范圍下NACA4415翼型升力系數和阻力系數 曲線