全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

M_∞=1及其附近翼型绕流的计算结果

本文全速势方程的有限差分数值计算结果,给出了0°和2°攻角下NACA0012翼型当地马赫数分布冻结时的自由流马赫数范围,以及M_∞稍大于1变至1.30时前方脱体激波的变化位置。     

跨声速翼型有粘／无粘强干扰计算

本文介绍一种二元跨声速激波-边界层强干扰的计算方法。边界层计算采用湍流边界层积分反方法,它借助Whitfield和Swafford提供的既适合附着流,也适合分离流的速度剖面表达式。跨声速无粘流用全速势方程模拟。通过边界上排溢速度来考虑粘性的影响,用有粘/无粘迭代得到粘性流解。本方法计算的结果与其它方法以及实验的结果进行了比较,证明该方法可以在工程上推广使用。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

计算机翼跨音速绕流的加罚Galerkin有限元法

1.引言 全速位Galerkin有限元法是跨音速有限元法中重要方法之一。处理超临界流时,为了自动捕捉到激波,通常采用上风技术,但因此导致了质量守恒条件得不到完全满足。文献[1]从Bateman变分原理出发,应用加罚函数方法,找出了强制满足质量守恒条件在Galerkin有限元方程中所应加入的附加项。本文将文献[1]的方法应用到超临界跨音速绕机翼流动,克服了一般Galerkin有限元法中由于采用上风技术而出现的激波不够陡和位     

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。     

管道中复杂马赫反射流场的高精度数值模拟

本文分析了一类高精度迎风TVD欧拉方程差分格式的精度,并用240×80个网格点进行马赫反射的计算,获得了目前实验和计算不易得到的两个激波三叉点的复杂波系结构。本文还研究了激波管上壁面对马赫反射的影响,并计算了轴对称扩张管道中激波的运动及激波与半球形障碍物的相互作用。计算结果分析指出,本文所采用的高精度迎风TVD格式与计算方法可以用来模拟具有复杂波系的工程问题研究。     

联合使用激波捕获—拟合法求解激波间断的跨声速流场

吸取激波捕获和激波拟合两种方法各自的优点,发展了一种求解有激波的跨声速流场的方法,并对回转叶栅中定常跨声速流场做了计算。计算中,通过捕获法确定初始激波位置,然后经反复修正,拟合出确定的流场通道激波。计算表明,这种方法可自动获得清晰的激波,流场中气流各参数分布合理。这一方法可适用于复杂边界和不同进口M数,计算时间仅比势函数方法多一倍左右。     

翼-身组合体跨音速绕流全位势积分方程数值计算

 用积分方程方法求解 Prandtl- Glauert算子表示的全位势方程 ,并计算了翼 -身组合体跨音速绕流。用 Murman- Cole差分格式计算空间场源强度 ,以捕捉激波。计算结果与相应的实验结果符合良好     

修正的跨声速小扰动方程以及与其相容的激波方程

本文给出了与修正的跨声速小扰动(MSD)方程相容的激波方程,计算证明它的极曲线与精确的激波极曲线十分接近。     

绕升力机身跨音速流动的有限元素法数值分析

 <正> 本文应用全速位方程最小压强积分的有限元素法解绕升力机身的跨音速流动。机身头部可以是尖的或带有进气道的。用头部伸出无穷长圆柱来模拟头部进气。尖头机身在元素形状处理上具有一定的复杂性,从数值实验角度考虑我们分别按尖头和头部伸出无穷长细圆柱来近似尖头头部两种方法对尖头机身进行计算。采用人工密度法捕获激波,但由于绕机身流线形状比较复杂,直接采用机翼问题中的人工粘性公式,解往往不收敛。我们推导出一种较精确的,适用于复杂流动情况的人工粘性计算公式,为了加速收敛,应用网格逐次加密技术,可使收敛速度提高2～3倍。本文还给出绕升力机身跨音速流动的远场速位解析式。     

三维跨声速流粘性—无粘相互作用问题的分析和计算

本文对跨声速粘性-无粘相互作用问题进行了分析,论证了如何选择湍流边界层反模式的输入量。在FL027多重网格基础上作了局部加密。计算例子表明,本文建议的半逆方法适用于弱激波情形,可用于实际设计之中。     

运动激波绕射过程中的粘性效应

应用 TVD格式求解二维非定常全 N- S方程 ,以研究运动激波绕射现象。考虑粘性后 ,绕凹角流动时 ,压缩拐角附近表面压力分布与实验结果更为符合 ;绕凸起物时 ,出现反射激波因边界层干扰根部加宽 ,前、后角点附近有局部分离与再附 ,并引起局部压力上升 (与无粘流场模拟相比 )等现象。结果表明 ,全 N- S方程可较好地模拟实际粘性流中复杂激波障碍物干扰流场特性     

大扰动跨音速流积分方程的导出及亚临界流的数值解法

本文对一方向上大扰动的平面跨音速势流微分方程提出了新的近似简化方案,从数学上完成了相应积分方程的推导;其中关于边界及激波条件作了小扰动近似处理;最后提供了关于亚临界流的算例。