跨音速压气机叶栅流动有粘-无粘相互作用的迭代计算

 本文利用较为先进的附面层内法向速度分布等半经验相关式,发展了一种适用于大逆压梯度的附面层积分方法。将激波近似地用大逆压梯度代替,对湍流附面层进行了快速计算,再将求得的附面层参数与无粘主流计算通过建立适当的数学模型结合起来,反复迭代后将计算结果与实验数据及无粘计算的结果比较,可以看出有粘-无粘迭代对无粘计算结果有明显改进。     

跨音速翼型绕流有粘-无粘迭代计算

本文提出一种新的跨音速有粘-无粘干扰迭代的算法。 外流计算利用文献[1]提出的一种快速算法。粘流(附面层)用Cebeci-Smith盒式法,湍流用简单的代数旋涡粘性模型描述。跨音速有粘-无粘迭代是以位移厚度为迭代参数,不断迭代直至收敛。 本文对零攻角无分离NACA0012翼型进行了计算,结果与实验吻合得比较满意,且收敛快,数值稳定,可期在实中得到应用。     

平面跨音叶栅中粘流-无粘流的强相互干扰计算

为了设计高效率高负荷的跨音叶轮机械,气动计算中必须计入粘性的影响。当叶片表面流动不发生分离,或分离区很小时,流道内的流动可以分为叶片表面以及尾缘下游的粘性薄层——边界层和尾迹区以及除此以外的无粘流区。因而可以分别应用欧拉方程组及边界层方程组求解这两个区内的流动。无粘流-粘流的相互影响是通过无粘流计算及边界层计算的一系列迭代计算来实现的:由给定的气动参数及叶栅几何形状进行无粘计算,得到叶片及尾迹表面的流动参数分布,作为首次边界层及尾迹计算的输入数据,求得边界层及尾迹的排挤厚度δ_1,并以此对无粘流动的几何边界进行修正,作为下一次无粘计算的流场边界。这样的迭代重复进行,直到满足全场的总体收敛条件为止,其迭代框图如图所1示。     

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

高负荷对转压气机中尾迹传播机制研究

为了研究尾迹在高负荷对转压气机中的转播机制,采用频域的非线性谐波方法对两级高负荷对转压气机进行了非定常流计算。结果表明:在设计点,通过与定常流计算的结果对比,发现第二级转子在0%~25%的近叶根区域和67%~100%近叶尖区域,非定常得到的效率较小,而在25%~67%的叶中部分非定常的效率更大。研究了尾迹在这三个不同区域的传播规律。结果表明:在叶中,尾迹恢复效应在流动中起到了效率增益的作用;在叶尖,它会与机匣的附面层以及叶尖泄漏流作用,使得端壁附面层增厚和引起泄漏流熵增;在叶根,由端壁附面层产生的马蹄涡和分离涡会影响尾迹的传播,卷吸尾迹并导致尾迹破碎,引起掺混损失。     

多段翼型失速分离特性的计算

本文用高阶奇点分布板块法(Panel Method》及边界层理论,对多段翼型的失速性能进行计算。采用分离尾迹的模型,按尾迹是流线用迭代法求出分离尾迹边界(此为内迭代)。在尾迹上配置奇点,然后把尾迹考虑在内进行粘/位流迭代(此为外迭代)。收敛后可得多段翼型空气动力特性,且可计算到失速以后的情况。算例有单段翼型GA(W)-1,NACA 4412,和F111飞机的三段翼剖面,所得结果与实验值能很好地符合。     

H型叶栅通道流计算网格的改进

提出一种壁面正交、尾迹区网格节点分布逐渐均匀化叶栅通道网格生成方法。为了研究壁面处网格正交性对数值计算结果的影响，推导出壁面附近流动参数一阶导数离散误差与网格线交角的数学关系。采用N-S方程分别对平板紊流附面层流和叶栅通道流进行计算分析。结果表明：壁面处网格正交性对附面层内流动计算预测精度影响较大，而对势流区流动影响很小；尾迹区切向网格节点逐渐均匀化分布能更准确计算尾迹区流动。     

尾迹/势流干扰对跨声速风扇静叶附面层流动影响的数值研究

为揭示尾迹/势流干扰下压气机静叶附面层的非定常流动机理,采用数值方法对单级高负荷跨声速风扇中径处的非定常流场进行模拟,针对尾迹/势流干扰下的静叶附面层非定常流动特征进行研究。基于尾迹/势流与叶片附面层干扰模型,通过叶片壁面摩擦力、近壁面附面层湍动能和壁面静压脉动,详细分析了尾迹和势流干扰下尾迹对高负荷静叶附面层流动状态的影响。研究发现:在非定常条件下,尾迹干扰能够"刺穿"静叶附面层,使得静叶附面层发生跨越转捩现象。尾迹和受其诱导的势流碰撞静叶前缘产生的压力扰动波在压力面附面层内以声速传播,影响压力面整个弦长的静压和摩擦力分布,而在吸力面上,扰动波的传播仅局限在前缘区域部分弦长上,吸力面和压力面传播现象不同与逆压梯度和曲面凸、凹形状相关。     

基于流场计算的叶栅性能预测研究

通过对搜集到的损失和落后角模型进行优化筛选并修正,发展了一套基于S1流面计算结果的适用于亚音速叶型平面叶栅性能预测的组合模型。同时,考虑了轴向密流比对平面叶栅性能的影响,并验证了关于轴向密流比和落后角关联式的正确性。流场计算程序使用了S1流面无粘计算程序和附面层计算程序,更好地模拟了流场。通过计算与实验对比分析,验证了损失和落后角模型的有效性,为此类叶型的平面叶栅特性预测提供了数值计算技术。     

计算压气机特性的基元叶片法

本文介绍一种轴流压气机非设计性能的预测方法。该法特点是 ,在流线曲率法解径向平衡方程基础上 ,引入环壁附面层损失的计算 ,考虑叶型存在对气流熵增和总焓变化的影响 ,并计入二次流 /径向间隙流引起的损失 ,激波、尾迹损失等。结果与试验数据对比符合工程要求。     

透平跨音速叶栅的优化设计

将遗传算法用于透平跨音速叶栅的优化设计,提出了相应的数学模型和优化方法。优化算例表明,本文所提出的模型及算法能有效的优化叶型方案,避免了一般优化算法的局部最优或维数灾难问题,从而在新的领域内将优化计算在叶轮机械设计中的应用进行了有益的探索。      

跨音速粘流的计算

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

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

An inverse code based on a residual-correction method for wing design

A fast, efficient and reliable code for wing design has been developed at INTA coupling a residual-correction method by Bauer, McFadden and Garabedian to an inviscid solver for wing analysis in transonic flow. Smoothing procedures, including Bézier cubic splines, are used to avoid irregularities of the wing surface, as well as the twist distribution. A non-conservative full potential finite differences code is used as the flow solver. The redesign of two wings in transonic flow are presented with promising results. The computing time required to achieve convergence up to a good level of accuracy is less than 1 hour in a Workstation of 300 Mflops. The code is actually a powerful tool for preliminary wing design. Future work will be addressed to the addition of viscous effects and to studies on the optimization of target pressure distributions.     

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

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

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.     

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

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

跨声速涡轮叶栅三维流场的数值模拟

以具有试验测量数据的一个高压涡轮动叶中截面平面叶栅和一个高压涡轮导向器环形叶栅为对象,采用计算流体动力学(CFD)软件CFX-TASCflow,对两个涡轮叶栅的跨声速三维粘性流场进行了数值计算,验证了CFX-TASCflow软件在航空燃气涡轮设计计算中的可靠性,并对所研究的跨声速涡轮叶栅粘性流场进行了分析.      

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

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