N-S方程数值模拟弹体不可压脱体涡绕流

应用在连续方程中增加压力对时间导数项的拟压缩性方法 ,数值求解三维定常不可压N- S方程。采用隐式 Beam- Warming近似因式分解格式求解。以修正的 Baldwin- Lomax代数模型模拟湍流流动。对绕正切拱形旋成柱体的大迎角脱体涡流动进行了数值模拟 ,并分析了层流粘性和湍流粘性对脱体涡及二次分离涡的影响。在头部施加微小的非对称射流扰动可得到与实验相符的非对称脱体涡     

前体涡非对称分离机理及前缘吹气控制研究

通过设计对称性算法,求解层流Navier-Stokes方程,数值模拟了细长体在低超声速情况下前体背风涡随攻角演化的规律,在此基础上,进一步研究了前缘吹气对前体涡演化和侧向力特性的控制机理.根据数值模拟结果及分析,倾向于支持在层流框架内,前体涡的非对称失稳是一种对流不稳定机制,要想根据需要产生对称或不对称的前体涡,就必须外加持续的扰动.在约16°～48°攻角区间内,前缘吹气可产生规律性较好的侧向力,有可能直接利用前体涡进行横侧向控制.为工程实用化,需提高前缘吹气的激励收益效费比.     

空间发展混合流的大涡模拟

采用大涡模拟方法数值模拟了空间发展混合流，通过在平均流中引入小扰动，研究了混合流大涡结构的产生和演化过程，捕捉了展向涡的卷起、配对、合并，以及二次流向涡的出现等大尺度的三维拟序结构。与实验及直接数值模拟的结果相比，吻合程度较好，表明本文所采用的LES方法对湍流的数值模拟是切实可行的。     

细长旋成体背风面非对称流场的DES模拟

利用脱体涡模拟DES方法对细长旋成体非对称绕流进行数值模拟,比较了基于S-A湍流模型的DES与RANS方法对背风面非对称分离涡的数值模拟能力,发现在大迎角非对称涡未破裂情况下,DES与RANS均能模拟出与实验相符的非对称绕流;当迎角增大,背风面非对称涡发生破裂时,RANS无法准确捕捉到背风面流场的非定常性,而DES能准确预测非对称涡的飘起与破裂,并与实验值接近。计算结果表明:与RANS方法相比,DES方法具有更好的模拟大攻角,大分离流动的能力,尤其在非对称流场的分离涡破裂模拟方面具有明显优势,能够更真实模拟出细长旋成体背风面分离涡破裂之后的非定常流动特征。     

头部钝度对大迎角非对称多涡流动特性的影响

在西北工业大学NF 3风洞通过表面测压和物面油流实验以及在北京航空航天大学1.2m水洞利用染色线显示和激光片光技术,对尖头和钝头拱形细长体大迎角绕流结构和气动特性进行了研究。结果表明,尖头拱形细长体在大迎角下(45°≤α≤60°)表现出非对称多涡的复杂涡系流动现象,相应的截面侧向力CZ沿轴向呈现出类似正弦曲线形式的减幅振荡,并且侧向力幅值较大(1相似文献   

一种强螺旋流现场的数值试验研究

螺旋流一般通过切向进流、安装导流片或旋转管道三种方式产生，但Horii等人通过实验发现有三种装置也产生了非常稳定的螺旋流。本文利用数值模拟对其中的一种螺旋流形成过程与机制进行了研究，计算结果证实了双涡结构向单涡结构演化现象并揭示了收缩管道的优点以及不同倾向进流的区别。在数值模拟时本文分别选取了线性、二及三阶涡粘性涡流模式，但数值实验表明线性与二阶模式的计算结果差别不大，说晨在数值模拟湍流螺旋流时二阶涡粘性模式作用不大。而三阶模式则表现出具有预测对称性进流向非对称结构转化的能力。     

细长体空间涡系的观测

本文用流态显示技术和空间总压测量,对于钝头和尖拱形头部两个细长旋成体模型在中等和大迎角状态的复杂背涡系进行了实验研究。侧重于“二次分离区”的流动观察和测量。实验表明:在中等迎角时细长体背风面上存在有一对二次涡,它们的旋转方向相反;靠外侧的一个涡尺度较大,其旋转方向与同侧主涡相同,当同侧主涡破裂时,它也发生破裂。随迎角进一步增加此涡呈现出从物面上间断地形成和脱落。     

叶栅二次流旋涡结构与损失分析

采用三维粘性程序对某型动力涡轮的第一级进行了数值模拟, 模拟结果捕捉到了该涡轮级叶栅的内部流的流动细节, 展示了涡轮叶栅端壁和型面流动及叶栅通道内的三维流动结构.通过对叶栅中的二次流现象和流动损失机理的分析, 揭示了该涡轮级叶栅通道内二次流旋涡结构(马蹄涡、通道涡、壁角涡、尾迹涡、泄漏涡等)的演变过程, 以及旋涡结构对损失分布的影响.      

旋转飞行器固体火箭发动机内涡核流动的数值分析

为分析旋转飞行器固体火箭发动机内涡核流动,采用数值计算的方法对发动机流场进行模拟。计算结果显示发动机的旋转增强了轴向涡核流动,随着涡核的不断延伸,在喷管入口处,由于潜入式喷管结构以及喷喉直径较小,振荡运动随之加强。涡核在流经喷管时的不断增强,在喷管区域形成了非对称的压力形式,并且产生了实际中的扰动侧向力,最终导致了飞行器的章动不稳定性。     

高压涡轮导向器非轴对称端壁设计

为进一步降低高性能涡轮端壁二次损失,根据其内端壁二次流形成和发展机理,提出一种基于多段三角函数的非轴对称端壁正交造型方法。根据此造型方法,通过优化的方式对高压涡轮导向器进行了三维端壁设计,并对此具有工程约束的寻优设计结果进行了全三维N-S方程数值模拟,分析和讨论了采用的非轴对称端壁造型方法对二次流损失的影响。结果表明:采用多段三角函数的非轴对称端壁正交造型方法,能有效控制马蹄涡压力侧分支的流动和通道涡的形成,并能有效降低高压涡轮导向器的端壁二次流损失。     

消除小展弦比尾舵大舵偏非对称流动现象的方法研究

在前期的风洞试验结果中已经证实,跨声速条件下,俯仰舵偏角增大到一定程度时,呈开口布局的一对尾舵形成的翼尖涡会相互干扰,使得两侧的翼尖涡强度不同,形成非对称流动现象,进而使全弹产生较大的横侧向气动力,对尾舵的控制能力提出了挑战,需要开展消除非对称流动的方法研究。本文基于数值模拟方法对非对称流场随舵偏角的变化进行了研究,并对该非对称流动现象产生的机理进行了分析,基于此,在保持舵面效率不降低的情况下开展了消除非对称流动现象的方法研究。研究结果表明,随着舵偏角的增加,呈开口布局一对尾舵的翼尖涡距离逐渐减小,翼尖涡由弱的不对称性,逐渐发展直到其中的一侧翼尖涡消失,研究结果还表明通过尾舵前缘局部切角和降低尾舵根弦长的方法在保持尾舵效率不降低的情况下,能够有效消除非对称流动现象,对类似布局设计具有指导意义。     

涵道螺旋桨桨叶涡强度分布计算

本文依据开式螺旋桨理论 ,涵道后缘的库塔条件以及微元涡丝的诱导方程 ,建立了涵道螺旋桨叶片涡强度分布的计算模型并进行了数值计算。结果显示 ,涵道螺旋桨桨叶涡强度在桨尖处达到最大     

细长旋成体大迎角绕流的背涡结构与压力脉动特性

在亚临界流动范围内,通过尖拱细长旋成体在无侧滑状态下的脉动压力测量,对大迎角下细长旋成体非对称背涡结构沿轴向的压力脉动特征进行了研究分析,揭示了压力的脉动幅度与频率沿旋成体轴向的变化规律及其与背涡结构之间的对应关系,并对压力脉动特征沿轴变化的原因进行了初步探讨.     

侧滑边条翼旋涡和气动特性实验研究

本文综合一组小展弦比边条翼实验研究结果,给出不同展弦比机翼前缘脱体涡系变化,并分析涡绞合及破裂对气动特性的影响。研究表明,侧滑使机翼迎风侧涡绞合推迟,涡提前破裂,背风侧正好相反。侧滑时不对称涡破裂对气动特性影响显著。     

Excitation of instabilities in the wake of an airfoil by means of active winglets

Strong wake vortices that develop behind every aircraft as a byproduct of lift production pose a threat where aircraft fly in close staggering such as in the vicinity of airports. One approach to alleviate these vortex wakes is the use of high lift systems or control surfaces of the wing to create an unstable vortex system. The inherent instability of this vortex system shall then lead to an accelerated decay of the vortex wake, triggered for example by a periodic motion of the control surfaces. In the work presented here a simple wing model with winglets able to produce a vortex system of up to six distinct vortices is investigated in towing tank experiments. Theoretical studies show that these vortex systems potentially have a high degree of instability. By means of active oscillation of rudders integrated into the winglets, these vortex systems are to be excited to initiate an accelerated decay of the vortices. It is shown that configurations exist which exhibit strong instabilities, that lead to a significantly lower hazard level behind the vortex generating wing, even when not actively excited. However, an additional oscillation does not seem to accelerate decay of these vortex systems in relation to the statical reference case.     

Effects of wing locations on wing rock induced by forebody vortices

Previous studies have shown that asymmetric vortex wakes over slender bodies exhibit a multi-vortex structure with an alternate arrangement along a body axis at high angle of attack. In this investigation, the effects of wing locations along a body axis on wing rock induced by forebody vortices was studied experimentally at a subcritical Reynolds number based on a body diameter. An artificial perturbation was added onto the nose tip to fix the orientations of forebody vortices. Par-ticle image velocimetry was used to identify flow patterns of forebody vortices in static situations, and time histories of wing rock were obtained using a free-to-roll rig. The results show that the wing locations can affect significantly the motion patterns of wing rock owing to the variation of multi-vortex patterns of forebody vortices. As the wing locations make the forebody vortices a two-vortex pattern, the wing body exhibits regularly divergence and fixed-point motion with azimuthal varia-tions of the tip perturbation. If a three-vortex pattern exists over the wing, however, the wing-rock patterns depend on the impact of the highest vortex and newborn vortex. As the three vortices together influence the wing flow, wing-rock patterns exhibit regularly fixed-points and limit-cycled oscillations. With the wing moving backwards, the newborn vortex becomes stronger, and wing-rock patterns become fixed-points, chaotic oscillations, and limit-cycled oscillations. With fur-ther backward movement of wings, the vortices are far away from the upper surface of wings, and the motions exhibit divergence, limit-cycled oscillations and fixed-points. For the rearmost location of the wing, the wing body exhibits stochastic oscillations and fixed-points.     

波瓣混合器涡系结构及射流掺混机理的数值研究

借助流体力学软件ANSYS CFX,对波瓣混合器射流掺混流场进行了全三维定常数值模拟,研究了流场中各涡系结构的形成机理及发展过程,并详细探讨了其加速射流掺混过程的作用机制.结果表明:基于SST(shear stress transport)模型的封闭N-S方程能较好地模拟波瓣混合器射流掺混过程,波瓣特殊几何外形诱导产生的流向涡主要通过扭曲内外涵交界面的间接方式加速射流掺混过程,波瓣下游剪切层中K-H(Kelvin-Helmholtz)不稳定性发展而成的正交涡是直接加速射流掺混的关键因素,波谷附近二次流之间的相互作用所产生的通道涡对该区域内的射流掺混有明显的加速作用,受波瓣前缘切割的边界层在径向压力梯度作用下沿波瓣表面卷起而形成的马蹄涡对射流掺混的影响不是特别明显.      

Equilibrium state of trailing vortices: two-dimensional model of a far-field vortex wake

The equilibrium statistical mechanics of a system composed of a large number of two-dimensional point vortices is employed to describe the vortex system shed from aircraft wings. According to this theory, these higher energy states of the vortex system can only be achieved by segregating the point vortices of like kind into two clusters that descend with a constant velocity. The solution is given in terms of the integral constraints for each cluster: total circulation, center of inertia, and kinetic energy. The negative non-dimensional inverse temperature of the system and the length scale related to angular momentum of a single trailing vortex are obtained versus initial interaction energy of the vortex system. Comparison of the theoretical results with available experimental data shows good agreement between the calculated tangential velocity distribution in the trailing vortex and the data. The flow characteristics for three different wing loads are also compared to emphasize the effect of initial circulation distribution along a lifting wing on the vorticity distribution in the equilibrium trailing vortices.     

压气机叶栅内二次流的数值研究

本文应用Beam-Warming近似隐式因子分解格式以及MML代数湍流模型,采用拟压缩性方法求解雷诺平均拟压缩N-S方程组,对压气机叶栅内三维粘性流场进行了数值研究。结果发现,在给定的流动条件下,叶栅内有多涡系组成的二次流旋涡结构。这些涡系的产生、发展、破裂以及相互作用显示了叶栅通道内三维流动的主要特征并对流动损失具有重要影响。计算结果与实验结果比较,两者吻合较好。      

翼身组合体大迎角前体非对称涡流动性态研究

在亚临界流动范围内,对于带有鸭翼、机翼的翼身组合体,在其头尖部带有确定扰动的条件下,研究模型大迎角下的非对称背涡结构及其气动力特性随扰动周向角的演化规律。通过对模型表面的压力分布和侧向力分布分析,结合流场显示结果,表明翼身组合体绕流中鸭翼前各截面均处于非对称二涡区,头部截面侧向力分布随头尖部滚转而呈现出双稳态特性,鸭翼和机翼上方的流动在大迎角下处于完全分离流动状态,从而使得模型上鸭翼之后的截面侧向力接近为零。