航空六三一所混合网格自适应计算技术研究进展.pdfVIP

中国第一届近代空气动力学与气动热力学会议论文集 CARS一200600074 航空六三一所混合网格自适应计算技术研究进展 邱征齐志敏白文 (中国航空工业第六三一研究所航空气动数值模拟重点实验室，西安710068) 摘要：实现了四种自适应探测器：基于速度梯度的尾流探测器、基于总压比的激渡撩测器、基于熵的涡探测 器、基于速度梯度张量特征值分析的涡探测器。利用尾流探测器，实现了P,AE2822翼型尾流区域附近的网 格自动加密。通过比较速度梯度自适应加密前后的网格，说明基于速度梯度的尾流探测器能够较好地捕 捉尾流区域。利用激渡探测器．实现了ONERAM6机翼激渡区域附近的网格自动加密。通过比较总压比 探测器自适应加密前后的网格。说明基于总压比的激渡探测器能够较好的探铡激渡区域，从而改进机鬟表 面的压力分布。利用涡探测器，实现了三角翼祸流区域附近的网格自动加密。通过对比不同的祸探利器 所产生的网格和无粘计算结果，说明涡探测器都能较好地捕捉涡区域，并且改进三角翼上的压力分布。其 中，基于速度梯度张量特征值分析的涡探测器更关心祸核，并能长距离探测尾涡。基于速度梯度张量特征 值分折的提探测嚣是局部探测器，阉值变化影响的主要是涡区域：而基于熵的涡探测器是全局操测器，荫 值变化影响的是整体网格。 关键词：自适应；探测器；加密；墒；速度梯度张量特征值分析 度张量特征值分析的涡探溅器‘2“1，从网格和无 1引言 粘计算结果两个角度对比分析了三角翼涡流模 自适应网格方法主要分为h一方法、P一方 拟的效果。并与NASA的无粘计算结果相比较。 法和r--方法三种类型，其中h一方法是通过细 2算法 划网格单元实现的，又可分为NodeInsertion、Re- H—Rcfmement方法可以通过不同的探测器 ￡nement和Remeshing三种方法。与h一方法相 关的自适应准则有物理量梯度准则和离散误差 实现。 估计准则u，。 2．1基于速度梯度的尾流探测器的实现 为了提高尾流区域网格分辨率，本文采用基 速度梯度可以作为尾流探测器。定义如下： 于速度梯度的自适应罔格加密方法。在尾流部 VSensor= 分区域，虽然马赫数变化不大，但是速度方向变 化却很大。因此梯度计算使用的是速度矢量的 各个分量，而不是速度矢量的模或马赫数。 其中第i点速度为吒，来流速度定义为‰。 在三维Euler方程混合网格自适应计算中，为了确保自适应探测器的阚值与算例的无关性， 激波探测器的选择可以是速度梯度、密度梯度、压 根据来流速度做归一化处理。 强梯度或者三者的综合，但在三维N—S方程混 2．2基于总压比的激波探测器的实现 合网格自适应计算中，如果选择速度梯度作为激 总压比可以作为激波探测器。该探测器定 波探测器，由于附面层黏性，将沿固壁附近的网格 义为‘’】 ～ 充分加密。而激渡附近网格却加密不足。为了提 老=，(一+(孚)讹2)寿 高激波区域的网格分辨率，准确地预测激渡位置， 本文采用基于总压比探测器的自适应网格加密方 总压比反映了压力经过激波位置后的损失。 法，同时采用了区域自适应进行补偿加密。 这个值越小，说明压力损失越大。 为了提高涡流区域的网格分辨率，准确地预 2。3基于熵的涡探测器的实现 测涡流，可以采用基于不同探测器的自适应网格 嫡经常放用来绘制涡等值面。同样，它也可 加密方法。本文分别采用基于熵和基于速度梯