用智能物面对非定常分离流进行最优自适应控制.pdfVIP

用智能物面对非定常分离流进行最优自适应控制 42 用智能物面对非定常分离流进行 最优自适应控制。 吴锤结1，2，3王亮2 (1中国科学院力学研究所非线性力学国家重点实验室，北京100080 2解放军理工大学理学院流体力学研究中心，南京211101； 3北京大学湍流与复杂系统国家重点实验室，北京100871) 摘要本文以二维大攻角翼型绕流为研究对象，将非定常动边界计算流体力学方法 与最优拉剞方法有机结合，研究可压和不可压非定常流智能物面最优自适应流动控 制的理论与算法，结果表明：在给定合适的最优控制目标函数下，智能物面可最优地 实时改变形状，得到能显著提高翼型性能的最优翼型。最优翼型在非设计工况下的 气动性能也比对照翼型有所提高。本研究为三维智能机翼的研究奠定了基础，提供 了初步的理论和算法。 ．一 关键词 非定常分离流 最优控制 非定常动边界计算流体力学 智能物面流动控 八 制 1 引言 对各种复杂流动的主动控制，如微尺度白适应流动控制(MAFC)、电流变液控制、利用记 忆合金等智能材料进行控制等，是近十年来国际上迅猛发展起来的一个崭新方向。对流体力 学而言，非定常控制技术已成为当前最主要的研究方向。 本文以翼型动态失速控制为例，研究可压和不可压非定常流智能物面最优自适应控制的 原理和方法，对位于大攻角下的NACA0012翼型进行非定常实时控制，得到非定常最优翼型 形状。这一研究不仅对动态失速控制有重要作用，还使我们进入智能翼型研究的行列。所得 结果不难推广到其他复杂非定常分离旋涡流动的控制，对理解生物运动中的流体力学机理(如 鱼类游动)也有直接帮助。 2．计算流体力学方法与最优控制算法 2．1求解不可压Navier．Stokes方程的数值方法 本文采用的求解不可压Navier-Stokes方程和连续性方程的算法参见文献[4]。 2．2求解可压缩Navier．Stokes方程的数值方法 设无量纲形式的Navier—Stokes方程为： o中国科学院力学研究所非螭性力攀目采重点妄鼍室和目寡自拣科学基金赉助项目( 2003空气动力学前沿研究论文集 f挚+v． pU)=0 1 dt J挚+可(pu圆H)+VP=V·(ys) (1) 【dt J掣+v T) L Ot ((卢E+P)H)=V·(ps“)+V·(ⅣV 其中，E为总能量；s为变形率张量；卢为黏性系数。各自的定义分别为：E=T+号，p2(y一 sut姒nd定律给出产讪(最门篝等)。 我们所用的算法[1]在时间离散时采用显示格式，在空间离散时采用迎风格式。利用有限 体积一Galerkin方法来迭代求解如下守恒形式的Navier．Stokes方程： 掣+可．(F(w)一N(w))：0 (2) 其中，w=(P，pul，pu2，庳)一，F和N分别代表对流和黏性算子。 如图1所示，计算区域n既可以用三角单元离散为哺=U，t，又可 以用多边形单元离散为Wh∈vh。因此，我们可以联想到每一个ZUh ∈U，其中K是三角网格上的连续函数。在多边形单元上的分片 常数函数Wh4为毗ICi2i—打J，um。反过来，如果知道了分片常 数孤，那么Wh(S：)=Wh4『G。寻找wj∈(Vh)6满足V九∈Ⅵ有 图1求解方程所