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

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

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

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

粘流与无粘流的相互作用计算

本文总结了粘流/无粘流的各种计算方法和结果。重点在于介绍定常流动中的弱相互作用。首先叙述了弱相互作用的数学模型。给出了不可压流动和跨音速流动中粘流/无粘流相互作用的某些正耦合的计算结果。讨论了在分离区附近边界层正方法失效的原因。然后介绍了边界层反方法和适用于带分离的流动中半反方法耦合的粘流/无粘流的相互作用方法。文中也简单地总结了三维情况的应用和强相互作用。     

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

跨音速非等熵全位势方法及粘流/无粘流相互作用计算

跨音速流动是现代航空军机和民机设计中必须计算的一种流动。目前采用求解N-S方程的方法尚不能提供很经济而又非常精确的解,也就不能成为工程设计使用的常规手段。在实际飞行中,高雷诺数的事实表明,粘流/无粘流相互作用计算的方法可以成为工程设计计算的有力工具。 欧拉方程是描述无粘流的准确方程,但与全位势方程相比其求解过程要复杂得多,计算费用也高。为取得最优外形工程设计,要作大量的流动模拟计算,故推广全位势方法并提高其准确度是很有意义及具吸引力的。全位势方法假定流动是等熵和无旋的,故只适用于无激波或仅有较弱激波的流场计算。为改进传统全位势方法,应对其作两方面的修正。分析表明,旋度修正比熵修正更为高阶,作为改     

Euler方程与边界层积分方程耦合求解跨音速翼型绕流

 应用Euler方程求解跨音速翼型特性时考虑了粘性影响,粘性影响是通过边界层动量和能量积分方程求解的,即粘流/无粘流迭代方法。其中Euler方程采用LU-ADI方法求解;边界层方程均由正解法过渡到反解法,以解决强激波干扰区出现小分离泡的计算问题。计算中使用了贴体C网格,通过一定变换使其保持基本正交。计算结果表明,压力分布、摩阻系数分布与实验结果符合较好。     

改进跨音速位势方法的讨论

 简要地回顾并讨论了传统位势方法之不足。对改进传统位势方法的两种必要的修正(熵修正和涡修正)的量级在小扰动条件下作了分析和估值。结果表明涡修正是更高阶小量,因而作为第一步修正可采用非等熵位势流来近似。采用熵修正激波算子来计及非等熵的激波突跃条件,并考虑了由熵修正带来的尾迹条件的变化。算例表明非等熵位势方法能得到预期的改进位势方法的效果。     

超音速钝体无粘绕流的计算

研究用三维特征线相容性关系的迎风差分格式,数值模拟超音速和高超音速钝体无粘绕流流场。用钝体反方法计算头激波及其波后的流场,用轴对称流特征线法确定初值面上的流动参数。算例的计算结果与实验数据比较吻合     

用粘流－无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进。改进的半反－局部联立耦合方法，考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好。算例表明，翼型高升力状态计算必须包括尾流的作用，也应当对势流计算压力与实际压力之间的差别进行修正。     

三维特体绕流数值模拟中薄层模型的应用

介绍用于近似计算绕三维尖无粘流激波薄层的二层模型。本法以数值确定激波的位置。在计算不同结构形状物体的空气动力特性时能确定气体动力学载荷的近似值。采用这种方法实际上排除了超高速理论中其他已知方法对尖体几何形状所加的约束。     

考虑粘性相互作用的跨音速机翼分析

本文将一些已有的三维全位势方程计算方法和一种三维积分边界层方法作了弱相互作用的耦合计算。计算了三种机翼并与已有的试验数据作了比较。     

一类轴流叶轮机转/静干扰问题的 3D Euler 数值研究方法

本文将作者发展的定常N-S解程序,改编成时间精确Euler解程序,用于研究转/静干扰问题。对一台高负荷风扇级无粘流场进行时间精确模拟的结果表明,该方法可被应用于轴流叶轮机内非定常流动机理研究。      

A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.     

跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

一种基于三维粘性流场的乘波体生成方法

为了便于乘波体的优化设计工作,本文实现了一种基于任意三维粘性流场,在激波面上截取乘波体前缘线来生成乘波体外型的设计方法,并与NASA Langley研究中心的CFL3D软件的算例结果进行了比较检验,误差满足工程计算要求.由于在流场计算中直接计入了粘性,而不像大多数的研究方法那样通过无粘计算再进行粘性修正,而且对于高超声速乘波体设计,可以在程序中考虑真实气体效应、化学反应等因素,因此本文的设计方法更接近工程实际,更具有推广价值.     

跨音速翼型分离绕流的粘性/无粘干扰迭代计算

本文用有粘/无粘干扰迭代的概念计算了跨音速任意翼型的绕流问题。位流的速位方程用AF2格式求解,而边界层微分方程用C-S盒式法求解,逆算法的引用可以克服边界层方程在分离点处的奇性问题,对分离区湍流代数模型的修正可以得到与实验更吻合的结果。计算结果表明有粘/无粘干扰迭代概念在小分离泡的情况下也是适用的。     

Viscosity effects on weak irregular reflection of shock waves in steady flow

The viscosity effects on weak shock wave reflection are investigated with the Navier–Stokes and DSMC flow solvers. It is shown that the viscosity plays a crucial role in the vicinity of three-shock intersection. Instead of a singular triple point, in viscous flow there is a smooth three shock transition zone, where one-dimensional shock jump relations cannot be applied. At the flow parameters corresponding to the von Neumann reflection, when no inviscid three-shock solution exists, the computations predict an irregular shock-wave configuration similar to that observed previously in experiments. The existence of a viscous zone in the region of shock-wave interaction allows a continuous transition from the parameters behind the Mach stem to the parameters behind the reflected shock, which is impossible in the inviscid three-shock theory. Results of numerical simulations agree qualitatively with the conclusions of the approximate viscous theory, which describes the flow in the vicinity of the three-shock intersection.