基于混合笛卡尔网格流动偏转器数值研究.docVIP

基于混合笛卡尔网格流动偏转器数值研究 　　摘要：采用混合笛卡尔网格（Hybrid Cartesian Grid，HCG）方法对流动偏转器进行了数值模拟和参数优化。混合笛卡尔网格由背景笛卡尔网格和贴体结构网格结合而成。通过计算加装流动偏转器的二维NACA0012翼型流场，并与干净翼型流场对比，验证了流动偏转器能有效控制流动分离、抑制失速的产生。最后采用基因算法对流动偏转器参数进行优化，较大的提高了翼型升力系数；同时充分体现了混合笛卡尔网格方法在处理多物体外形结构上方便、高效的特点。 　　关键词：流动偏转器 流动分离控制 混合笛卡尔网格 基因算法 　　中图分类号：V211 文献标识码：A 文章编号：1007-9416（2015）02-0000-00 　　流动分离现象在实际的工程应用中普遍存在。由于逆压梯度等原因，当流体绕流物体时，从某个位置开始脱离物面，并在物面附近出现与主流方法相反的回流，即为流动分离。分离产生后，压差阻力增大、升力减小，导致失速产生，严重影响飞机的气动性能和操纵性能。因此，控制流动分离的研究在航空工程中有着重要意义。 　　流动分离控制方法多种多样，目前较为常见的方法主要有等离子体激励器、合成射流等。南京航空航天大学的明晓教授提出了一种新的流动控制技术――流动偏转器，用实验研究的方法完成了流动偏转器的有效性验证以及部分参数设计 [1，2]；本文使用二维雷诺平均Navier-Stokes（RANS）方程求解器，对其进行了进一步数值研究。 　　流动偏转器尺寸很小，对于这种外形较为复杂的多物体粘性湍流流动问题，单纯使用结构网格则需分块处理，较为复杂；非结构网格则难以模拟边界层的流动；而笛卡尔网格方法对于固壁面的处理存在一定困难。因此本文采用混合笛卡尔网格方法，即边界层布置贴体结构网格、在流场其他区域布置笛卡尔网格。 　　1 流动偏转器 　　流动偏转器由安装在翼型上翼面前缘附近的一组平行小薄片构成，结构如图1所示，其中，n为薄片个数，w、t分别为薄片的宽度和厚度，d为薄片间距，x为偏转器转轴中心到翼型前缘的距离，L为最下层薄片到轴心的距离，θ为薄片与偏转器轴线的夹角，φ为轴线与翼型弦线的夹角。 　　流动偏转器的工作原理是改变来流流向，使其向翼型上翼面偏转，流体能够继续附着翼型上表面流动；同时使翼面上方动量较高的流体与翼面附近动量较低的流体相混合，增加边界层的动量，用以抵消逆压梯度引起的损耗，从而延缓分离的产生。 　　2 数值方法 　　2.1 控制方程 　　其中，W为守恒变量，Fc为对流通量，Fv为粘性通量。采用有限体积方法求解方程（1），对流项的离散使用低数值耗散的AUSM+格式[3]，粘性项使用二阶中心格式，时间离散则使用LU-SGS隐式迭代格式。湍流模型采用SST-2003 k-ω模型。 　　2.2 混合笛卡尔网格 　　混合笛卡尔网格[4，5]是由围绕物体几何外形的贴体结构网格和填充流场其他区域的背景笛卡尔网格构成的。结构网格和笛卡尔网格之间的信息传递通过网格交接面完成。 　　混合笛卡尔网格的生成主要包含以下步骤： 　　（1） 围绕每一个几何外形生成贴体结构网格，这些网格之间可能存在重叠部分； 　　（2） 根据各个贴体结构网格的位置，判断网格重叠区域并将重叠网格删去[6]； 　　（3） 在流场计算域内的其他区域填充笛卡尔网格，网格的几何加密次数根据贴体网格最外层单元的尺寸决定，笛卡尔网格和贴体网格之间会有重叠区域； 　　（4） 对笛卡尔网格单元的类型进行判断，删去在贴体网格内部的单元，最终形成混合笛卡尔网格。 　　计算域的混合笛卡尔网格如图2所示。 　　网格交接面处的数据传递直接影响计算的收敛性和精度，这里采用“贡献单元”方法[6]。 　　3 基因算法 　　基因算法是模仿自然界生物进化机制的随机全局搜索和优化方法。它能在搜索过程中自动获取和积累搜索空间的信息，采用适者生存的原则，在每一代种群中借鉴自然遗传的方法选择个体，产生新的种群，使得种群不断进化。每一代种群代表问题的解。它是一种高效、鲁棒的全局优化方法[7]。 　　本文中，以升力系数为优化目标，以流动偏转器的参数为决策变量，采用简单基因算法（SGA），其流程如下： 　　（1）确定解空间：即流动偏转器的优化参数以及变化范围。 　　（2）编码：采用二进制编码，将各个参数变量单独编码成二进制数并连接在一起，形成一个新的二进制数，即为个体（也称染色体）。 　　（3）产生初始种群：在解空间内随机生成一定数量的个体，组成初始种群。个体数称为种群规模。 　　（4）计算个体的适应度值：以升力系数CL为适应度函数，计算种群中每一个个体的适应度值，并按从高到低对个体排序。 　　（5）产生下一代种群：通过选择、交叉（cros