叶栅流动特性数值计算

叶栅流动特性数值计算》实验指导航空航天学部动力系2012年12月动量方程+div(pvu)=--P+div(Ugradv)+SdtdyMy能量方程+div(pwu)=-dtCzMzPdti实验一亚声速叶栅流动特性数值计算实验NumericalExperimentofFlowCharacteristicofSubsonicCascade一、实验目的1、初步学习计算流体力学理论知识，加深对叶栅流动理论的认识；2、初步掌握计算流体力学相关软件；3、掌握亚声速叶栅流动的数值模拟方法；4、掌握对数值计算结果的分析与处理方法。二、实验内容本实验内容包括以下三个部分：1、初步学习计算流体力学理论知识，包括N-S方程、离散方法和湍流模型。2、学习相关流体力学软件：GAMBIT和FLUENT；3、掌握亚声速叶栅流动的数值模拟方法，其中包括以下部分：1）计算域的建立与网格划分；2）亚声速叶栅流动的数值求解，；3）数值计算结果的处理与分析。三、计算流体力学理论知识1、N-S方程Navier-Stocks方程是描述粘性流体方程，简称N-S方程。其中，p是流体密度；p是压力；u，v，w是流体在t时刻，在点（x,y，z）处的速度分量。常数卩依赖于流体的性质，叫做粘性系数。S为源项。连续方程+div（Pu）=0状态方程P二P（P,T）i=i（p,T）2、离散方法1）有限差分法有限差分法是数值解法中最经典的方法。它是将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程（控制方程）的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。求差分方程组的解，就是微分方程定解问题的数值近似解，这是一种直接将微分问题变为代数问题的近似数值解法。2）有限元法有限元法与有限差分法都是广泛应用的流体动力学数值计算方法。有限元法是将一个连续的求解域任意分成适当形状的许多微小单元，并于各小单元分片构造插值函数，然后根据极值原理，将问题的控制方程转化为所有单元上的有限元方程，把总体的极值作为各单元极值总和，即将局部单元总体合成，形成嵌入了指定边界条件的代数方程组，求解该方程组就得到各节点上待求的函数值。3）有限体积法有限体积法又称为控制体积法。其基本思路是：将计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积，将待解微分方程对每一个控制体积积分，从而得出一组离散方程。其中的未知数是网格点上的因变量e。为了求出控制体积的积分，必须假定e值在网格点之间的变化规律。从积分区域的选取方法看来，有限体积法属于加权余量法中的子域法，从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。简言之，子域法加离散，就是有限体积法的基本方法。3、湍流模型1）Spalart—Allmars单方程模型因为湍流现象是高度复杂的，所以至今还没有一种方法能够全面、准确地对所有流动问题中的湍流现象进行模拟。在涉及湍流的计算中，都要对湍流模型的模拟能力及计算所需的系统资源进行综合考虑,选择合适的湍流模型进行模拟。Spalart—Allmars模型是一方程模型里面最成功的一个模型，最早被用于由壁面限制情况的流动计算中，特别是存在逆压梯度的流动区域内，对边界层的计算效果较好，因此经常被用于流动分离区附近的计算，Spalart-Allmars单方程模型是用于求解动力涡粘输运方程的相对简单的一种模型，它包含一组最新发展的单方程模型，在这些方程里不必要去计算和局部剪切层厚度相关的长度尺度。Spalart—Allmars模型是专门用于求解航空领域的壁面限制流动，对于受压力梯度作用的边界层流动，已取得很好的效果，在透平机械中的应用也越来越普遍。FLUENT对Spalart—Allmars进行了改进，主要改进是可以在网格精度不高时使用壁面函数。在湍流对流场影响不大，同时网格粗糙时可以选用这个模型。Spalart—Allmars模型的变量中v是湍流动粘滞率除了近壁区域，方程是：（pv）+丁（pvu）二G+—[丁｛（卩+pv）「｝+Cp（「）2]—Y+Sdtoxivooxoxb2oxvvivjjj这里G是湍流粘度生成的，Y是被湍流粘度消去，发生在近壁区域S是用户定义的。注意到湍流动vv能在Spalart-Allmaras没有被计算，但估计雷诺压力时没有被考虑。2）k一匕模型Standardk—8model是最简单的完整的二方程模型，在流场和热交换模拟中，由于其实用性、经济性和可靠的精确性而得到广泛的应用。它是半经验模型，其导出项取决于经验的积累和对现象的考虑。随后为了改进其性能，出现了RNGk一8model和realizablek一8model。RNGk一8model是用严格的统计技术推导出来的（称为重正规化理论），但有了以下一些细化：RNGk—8model中多了附加项，提高了模拟速度较高流动的精确度；旋涡对湍流的影响被包括在内，增强模拟旋涡流动的精确度；相对于Standardk一8model这一高雷诺数模型来说，这些特点使RNGk一8model在模拟大多数流动现象时比Standardk一8model精确可靠得多。3）雷诺应力模型（RSM）雷诺应力模型（RSM）是最精细的一种湍流模型，抛弃了等方涡粘假设。通过求解雷诺应力输运方程和扩散率方程，RSM使雷诺平均N-S方程封闭。这意味着对二维流动需要求解四个附加方程，对三维流动则需要求解七个附加方程。由于RSM与一方程模型和二方程模型相比，考虑了流线弯曲、旋涡、旋转和应力骤然变化的影响，因此对复杂流动更可能给出精确的预测。但在雷诺应力输运方程中，RSM预测的灵活性仍然受到模拟不同项的封闭假设所限制。四、叶栅理论知识1、叶型几何参数叶型是沿主气流的流动方向上叶片横截面的几何形状。在轴流风机和压气机中，叶片的叶型有机翼叶型以及圆弧板叶型等。叶型中线——叶型轮廓线的各内切圆的中心连线，称为叶型中线，如图1中的点划线所示。叶弦——连接叶型中线两端点的直线称为叶弦。叶型的前、后端点分别为叶型的前缘点和后缘点。b——弦长，叶型的长度。c――叶型厚度，叶型中线的法线与叶型轮廓线交点的距离，即叶型内切圆的直径，叶型厚度的最大值称为最大叶型厚度C。maxC.b=――叶型的相对厚度。max叶型弯度，叶型中线到叶弦间的垂直距离，叶型弯度的最大值称为最大叶型弯fmaxfb=――叶型的最大相对弯度位置。maxXf――最大弯度位置，叶型前缘点沿叶弦方向至叶型最大弯度处的距离。X――最大厚度位置，叶型前缘点沿叶弦方向至叶型最大厚度处的距离。CX二X.b=――最大厚度处的相对距离。cC'X叶型前缘方向角，叶型前缘点处叶型中线的切线与叶弦间的夹角。1X――叶型后缘方向角，叶型后缘点处叶型中线的切线与叶弦间的夹角。2e二x+X——叶型的弯曲角。12r,r――分别为叶型前缘和后缘半径。122、叶型气动参数叶栅前额线一一叶栅中各叶型前缘点的连线。叶栅后额线一一叶栅中各叶型后缘点的连线。t――栅距，二维平面叶栅中相邻两叶型上对应点的距离。工=b：t――叶栅稠度，叶栅稠度为叶栅的重要参数之一，在轴流风机和压气机中，t=0.2〜2，甚0――叶型安装角，为叶弦与叶栅额线间的夹角。i0――叶片安装角，为叶型外弦与叶栅额线间的夹角。0――叶型进口几何角，为通过叶型前缘点向叶型中线所做的切线与前额线间的夹角。1A0――叶型出口几何角，为通过叶型后缘点向叶型中线所做的切线与后额线间2Aw1——气流进口相对速度，即叶栅进口前未受叶栅干扰时的气流相对速度。w——气流出口相对速度，即叶栅出口外不受叶栅影响时的气流相对速度。2w=(w+w)2—―平均相对速度矢量。m120进口气流角，为气流进口相对速度与叶栅额线间的夹角。10出口气流角，为气流出口相对速度与叶栅额线间的夹角。20――平均相对速度矢量与叶栅额线间的夹角。mi攻角，为平均相对矢量与叶弦间的夹角。i=0-0——气流进口攻角。1A1§=卩一0气流出口落后角。2A2A0=0-0――气流转折角，表示气流流过叶栅流动方向发生的改变。21A0=0—0=(0—8)—(0—i)=0—0+i—8=0+i—8212k1k2k1k五、计算域的建立与网格划分1、计算域的建立首先建立S1流面计算域：根据平面叶栅流场特点，流场应为周期性流场，可以将数值模型简化为只分析两个叶片之间的流道之间的流场，叶栅叶型数据详见表1。本例将进口边界设置为1倍弦长，出口边界约为2倍弦长。这主要是由于叶栅流场对上游影响较小，叶栅尾迹对叶栅下游流场影响较大。采用UG对曲线进行数据拟合，再将拟合后的曲线导入AutoCAD进行调整，建立叶栅流场数值模型。也可以直接采用AutoCAD软件，将叶型的吸力面、压力面数据的X、Y坐标通过Excel合并之后，再输入AutoCAD中进行曲线拟合，并设定适当的进口和出口，建立叶栅流场模型。图2计算域2、网格划分将叶栅模型导入Gambit中，为了划分出较高质量的网格，在划分网格之前可以对计算域进行分区，分为叶栅上游，叶栅流道，叶栅下游三个部分，先画线再画面网格的顺序。（1）在叶栅流场参数梯度比较大的地方采用局部加密网格，叶片表面，叶栅前缘和尾缘处均采用加密网格，在加密网格区域尽量避免网格密度的突变。（2）在其它区域采用较疏的网格，以提高计算速率。要求网格原则上，最小正交性角度越接近90。越好，最大网格长度比越接近1越好，最大网格延展比尽量接近1。但由于边界层、尾迹等的存在，很难达到要求兼顾各个方面，所以按：最小网格正交性角度大于10。，最大网格长宽比小于5000，最大网格延展比小于10。。网格不能出现负网格，否则将无法进行计算。六、亚声速叶栅流动的数值求解1、湍流模型选择：选择标准k-8模型和雷诺应力模型2、边界条件给定（1）进口边界条件：工况下进口气流角47.82°，进口气流马赫数30m/s，攻角为-5°、0°、5°和10°，气流边界条件采用速度进口（Velocity_inlet）。（2）出口边界条件：气流出口边界条件采用压力出口（Pressure_outlet）.（3）固体壁面边界条件：叶盆、叶背均采用绝热的无滑移固体壁边界条件（Wall）。（4）周期性边界条件：叶栅前后延伸部分均采用周期性边界条件（Periodic）。3、控制参数的选择图3控制参数的设置4、收敛准则所监视的参数包括速度、能量、连续性残差、进口流量、出口流量。所设定的收敛标准为速度、能量连续性残差下降三个数量级以上，nut函数收敛标准定为106；进出口流相对误差小于）.5%，且流量不再发生变化。对于有大分离涡的流动（尤其在进出口处），流量收敛曲线会发生振荡，此时对于迭代中分离涡的位置和强度都会发生不同程度的变化，呈现非定常特性，因此流量也会随之发生变化。在这种情况下，也可认为计算收敛。对于有大分离的情况，总压恢复系数也会呈现周期性变化，这种情况下也可认为计算收敛。七、数值计算结果的处理与分析1、数值计算结果的显示2、数值计算结果的处理3、数值计算结果的分析八、数值计算准备及预习要求预习要求：详细阅读本实验指导书，达到1）熟悉本数值计算目的及基本原理2）清楚数值计算的基本操作步骤3）清楚本数值计算结果处理方法九、数值计算内容及格式1、数值计算目的2、数值计算内容3、数值计算原理4、数值计算具体步骤5、数值计算结果分析（包括压力、速度等参数的线图、云图、不同湍流模型结果差异分析等）十、相关参数1、叶型参数表1NACA65叶型参数吸力面压力面X(mm)Y(mm)X(mm)Y(mm)0.23920.6448800.007521.00480.4940878.1761.21.24880.5242473.9282.98081.68080.6043269.7044.3922.62560.7884865.485.69524.28960.883261.2646.86966.29440.946457.047.90168.281.04452.8088.77612.2241.283248.5769.48816.1441.536844.34410.00820.0561.780840.1210.34423.9682.00435.89610.45627.8882.231.69610.3631.8242.367227.52810.07235.7762.514423.3929.639.762.651219.2968.94443.7922.785615.2488.10447.8642.889611.2567.06851.9922.94247.33125.813656.1842.92485.4045.087260.4322.81923.50564.276864.7442.60321.38482.9869.122.25440.528242.077673.5441.74240.25921.6478.0160.720.156721.392800.007520.23920.64482、叶片几何参数叶片弦长(b)(mm)表2NACA65叶栅参数3080叶栅进口速度G/s)弦长/栅距(b/t)/\1.6展弦比(h/b)/、1.2叶片安装角丫^)60叶片几何出气角a2'丿2P-11.52叶片几何进气角a.J1p47.82进口附面层厚度(mm)123、网格划分要求4、边界条件5、湍流模型与数值计算监控