跨音速翼型反问题的赝势函数变域变分有限元法.pdfVIP

跨音速翼型反问题的赝势函数 变域变分有限元法 李孝伟 刘高联 (上海大学上海市应用数学和力学研究所，上海200072) 摘要本文运用赝势函数变域变分有限元方法进行了跨音速翼型设计。赝势函数的引 入既保留了势函数模型在求解方面的优越性，同时又是势函数对有旋流动的一个自然的、物理 上相容、数学上求解简便的推广；变域变分是一种处理未知界面的独特工具，极适合于气动外 形的设计求解。踌音速翼型设计的实例表明，本文方法是有旋流动翼型设计的有效工具。 关键词踌音速，变域交分，有限元，赝势函数 1 前言 文献E23突破了势函数模型的局限性，引入了一种新的通用函数——赝势函数，既保持了 势函数模型在求解方面的优越性，又能有效模拟含有激波的跨音速有旋流动。该腰势函数是 势函数对有旋流动的一个自然的、物理上相容、数学上求解简便的推广。进一步地，文献[2]还 建立了赝势函数的一系列变分原理，为跨音速有旋流动的有限元方法求解奠定了理论基础。 另外，文献[3]针对未知界面问题还发展了一种变域变分有限元方法，它将未知界面结合 在变分泛函中，使其与求解流场的控制方程结合起来，从而未知界面的求解和流场的求解完全 耦合进行，因而是一种处理未知界面的独特工具，极适合于气动外形的设计求解。 本文就是运用赝势函数和变域变分有限元两大工具对翼型跨音速流动杂交命题进行了数 值研究。设计结果证明了赝势函数变域变分有限元方法的有效性和优越性。 2二维定常、无粘、有旋流动的无量纲化基本方程 连续方程 —a=(／。轧一D：0 d2‘ Y向动量方程 ·173· 署一壶刀dS／d咿+PdH／d尘=o (2) 能量方程 南詈+舭譬=Hc奶 (3) 流函数定义 蓑+阳z=o (4) Gibbs恒等式 (5) 吣=dE+Pd(1／p筹I，=T，丽aEIs一声 赝势定义 一2“1’ 东 Ⅱ“ (6) 丑～，老=U2--axl 其中，地为来流马赫数，E为内能，P为赝势，Ⅱ为涡势，S为熵，丁为温度，H为滞止焓， H(奶已知。 3赝势函数变域变分原理及有限元解法 为了计算方便，且不失普遍性，可认为边界变分只是沿着Y向进行，参看图l，于是 图1物面法向变分方向与弧长变化关系图 ·174· ＆一曲．7 ＆·；；·ds·Y．；击=曲如 那么参照文献[23，可以知道反命题的变分泛函J如下： (7) ，(尹’A，=Fd．mdy--舅。c莉+，cB+西乃，，山(f，c苁，y‘山 其中，F=再丧一罚老{1+z-尹[1一(老)2一(老+Ⅱ)2]) 如果去掉(7)式右端的第三项，得到的便是正命题的变分泛函。 12里；；半按一般的方法进行求解，不再重述。下面只叙述12；{；筹立的求解。 由式(7)知：