跨音速全速位积分方程数值计算.pdfVIP

。 一 义 1993年2月1日 气薪研究皇实验 第l期 跨音速全速位积分方程数值计算 董 军 (626所) 摘 要 本文培出一种利用格椿公式将跨音速全速位方程化为积分方程式，然 后进行数值计舞的j昆合方法 该方法将拣分方程方法与有限差分方法相结合用于二 维跨音逮流计算。在数值计算中采用上风技术 目口引入人工粘性．以此来 自动捕获 激渡 为减小计算误差 菩导敢量均在变换平面列求得 该方法空间计算区域小 计算收敛速度快一数值计算结果令人满意 主 题 词：墨型、堕童蕉 坌立蓬、格林套式、教照计算 O 引 言 跨音速毵流数值计算是航空 航天等高技术中的重要课题。。近些年来 采用积 分胄程方法计算跨音速流动有了一定的发展 但这些工作主要是针对跨音速小 扰动方程丽开展的。本文从跨音速全速位方程出发，萃嘻甩格林公式将其他为全速位 积分方程式 积分方程最终可分为面源顼和场源项之和，丽豫删采坩内奇牲画源法 离散。场辣顼徽为压缩牲修正 然后采用迭代方法求解。对于超临界流动，梳场中 出现局部超音速区，为了自动捕获到濑波采用上风技术，即引入九工粘性顼。为减 小计算误差 计算公式中各导数量均程变换平面内带有限差分方法求得 然后用雅 可比变换到物理平面 t 该方法以速度为基本未知量。跨音速全速位积分方程方法 与有限差分方法相比 其优点是计算重较小，计算收敛速度快，且远场边界条件 自 动得到满足，从而空间计算区域小 1 全速位方程 定常，无粘、无量纲形式的全速位方程 ( 审．)．(p中 )=0 其中 Y—l p=[1+二——札 (1一』V中 )] ’ 2 能量方程 l+ (1_lv ) 式中 p为流体密度，为当地音谴 中为总的逮位，Il-为自由流马赫数 Y 为比热比(Y瑕1．4)。 流场中速位中(x，y)由自由流U-的贡献和一个附加扰动谜位场 (x，y)所组成 中(x，y)=(x×C05a y×sind 1Il(x y) (2) 式中a为翼型攻角。 边界条件 定常、无粘气流流动的物面边界条件 V审 × 0 远场边界条件 中 ． ( ) 一 。。 (4) 式中 亩为翼型表面外法线． 为平行与自由流U_舶单位矢量 。 2 全速位积分方程 由(2)式，全速位方程(1]可写成扰动速位形式 中。． =9 其中 g 一 — p．× 中． p X 中．) (6) p 将方程(5)观作泊橐方程 利用格林公式和边界条件(4)可将其能为如下形式 的积分方程 ’