迎角直到60度旋转体涡流的多重线涡模型

本文用改进的多重线涡模型计算大迎角低速旋转体的对称和非对称涡流。通过给定物面上对称或非对称的分离线位置,现在提出的算法有效地解决了涡强度与自由涡线位置的迭代匹配问题,首次得到了迎角大到60度的涡流数值解。对一个切拱头体计算出的气动力特性与实验结果相符。     

不同迎角下脊形前体绕流数值模拟研究

脊形前体有较强的背风涡流场,不同的前体形状对前体涡流场和气动力有很大的影响。本文针对脊形前体飞行器大迎角湍流大分离流动计算的困难,采用IDDES混合湍流模型,以及与之匹配的非定常算法,研究了不同来流迎角下脊形前体的气动特性,以及背风涡非定常演化、破裂的细致流动结构。选取了不同脊形角,以及不同上、下高宽比的脊形前体进行计算。计算结果表明,在迎角较小时,随着迎角的增大,前体主涡会逐渐增强,在迎角较大时,前体主涡破裂;在相同迎角下,脊形角较小时,前体涡较强,涡升力也更大;对于相同脊形角的前体,当上半截面高宽比较小时,前体主涡强度较大,前体涡破裂临界迎角较小,即会提前破裂。     

大型飞机后体流动控制及减阻机理研究

通过数值模拟方法研究了大型飞机机身后体的流动特征和机理,并采用涡流发生器对机身后体流动分离进行控制,分析了其控制减阻的机理,讨论了不同参数对减阻效率的影响。计算结果表明:大型飞机机身后体流动分离是导致巡航阻力增加的一个重要因素;在后体底部附近安装涡流发生器,当来流迎角为负时减阻效果较明显,随着迎角的逐渐增大,减阻效果降低;涡流发生器的尺寸、周向距离、安装在机身后体上的前后位置、安装角度等参数对减阻效果具有显著影响。     

民用运输机短舱涡流片设计研究

采用已通过风洞试验验证的CFD流场数值模拟方法,从流动机理出发分析空间涡量的变化,研究了翼吊发动机短舱涡流片的安装位置对着陆构型升力特性的影响。数值模拟的结果表明:由于发动机短舱以及挂架的影响,缝翼下偏后会在主翼上留下的台阶以及缝翼的端面等,大迎角时诱导出许多空间涡系,容易引发主翼上表面的分离。短舱涡流片诱导的空间涡能够有效地抑制这些空间涡系和低速区,提高失速迎角和最大升力系数,对于文中着陆构型,失速迎角提高2°,最大升力系数提高0.15。短舱涡流片后移或者下移均会引起空间涡的下移,有利于抑制大迎角下主翼中段低速气流。但下移涡流片的同时会降低空间涡的强度,使其抑制作用减弱。因此,为了提高失速迎角和最大升力系数,在设计过程中需综合考虑短舱涡流片所诱导空间涡的强度和位置。针对某型民用运输机着陆构型中短舱涡流片因几何约束需对其位置进行重新设计的问题,根据上述研究结论,综合权衡空间涡的强度及位置,重新设计了短舱涡流片的位置。计算结果表明,重新设计涡流片的位置后,几何约束得到了满足,着陆构型的最大升力系数仅损失0.015,仍然能够满足设计指标。     

细长体大迎角非对称流动的高速PIV风洞试验研究

具有细长前体构型的飞行器在大迎角绕流中会出现明显的非对称涡系流动及其伴随而来的非对称力,该现象受多种因素影响,而其中对压缩性效应的研究相对较少。在0.6m亚跨超声速风洞中,采用PIV测量技术,对尖拱细长旋成体大迎角非对称流动开展了试验研究。试验M数范围为0.4~1.2,迎角为40°。试验结果表明:细长体模型在高速情况下仍然存在非对称多涡流动结构;Re数和压缩性均对非对称涡流动产生明显影响;模型头尖部人工微扰动与非对称涡之间存在确定的响应关系。     

二元翼型大迎角绕流的平衡态(解)的数值研究

基于k-ω的SST两方程湍流模型,在时间域求解雷诺平均Navier-Stokes方程(RANS)计算了翼型大迎角分离流动.通过给翼型施加不同形式的扰动,重点关注了不同扰动形式对翼型大迎角分离涡流场平衡态的影响.研究表明:在所给定的百万雷诺数条件下,当翼型厚度达到一定程度而且在一定迎角范围内时,翼型的分离涡流场存在多个平衡态,且在一定扰动条件下会发生平衡态转移现象.在不同振动形式上,沉浮和俯仰比前后振动更易激发平衡态的转移.正是这种分离涡流场平衡态的多值现象,可能引发风洞实验中大迎角气动力数据的分散.     

边条翼布局双垂尾抖振特性与机理风洞实验研究

对两种平面形状的边条翼布局模型分别作了双垂尾抖振实验和涡流场激光片光源显示实验研究。抖振实验测量了两种模型双垂尾的翼根弯矩响应和翼尖加速度响应,涡流场显示实验记录了两种模型上典型位置上的涡流场发展状态。通过边条涡流场随迎角的发展和破裂特性与模型垂尾抖振响应特性的对比分析发现:(1)垂尾翼根弯矩、翼尖加速度响应随迎角的变化均与边条涡的发展状态、是否破裂以及破裂程度密切相关;(2)主翼后掠角较大的情况下,机翼前缘涡与边条涡相互干扰,不但加快了涡的破裂使得双垂尾抖振起始迎角减小,而且使得垂尾的抖振响应较大。     

方形截面飞行器上仰机动对滚转特性影响的数值模拟

飞行器从中小迎角至大迎角范围内,由于背风区流动分离形态的演化,静态气动特性特别是横侧向气动特性也随迎角显著变化,可能诱发复杂的滚转运动。但飞行器一般是上仰机动时,才从平飞状态快速拉起至大迎角,此机动过程对横侧向气动特性和滚转运动可能产生较大影响。本文发展了刚体动力学方程和Navier-Stokes方程的松耦合求解技术,并通过数值模拟航天飞机脱落碎片的六自由度运动轨迹进行了验证。针对背风区涡流形态及横侧向气动特性复杂的方形截面飞行器,数值模拟研究了其不同迎角下的静态滚转气动特性、自由滚转运动特性,以及上仰机动时不同拉起速率对滚转运动特性的影响。结果表明,对于此飞行器,静态时存在临界迎角约为13°,当迎角小于临界迎角时,滚转方向是静不稳定的,诱发快速滚转运动;当迎角大于临界迎角时,滚转方向是静稳定的,其滚转运动是收敛的。但上仰机动时,滚转运动的形态还与拉起速率相关,即使拉起的终止迎角大于临界迎角,如果拉起速率较慢,也可能出现快速滚转运动。     

大迎角翼—身组合体涡流特性计算

论述了用非线性离散涡法来模拟大迎角下翼－身组合体涡流绕流的计算模型及计算方法；给出了一个典型翼－身组合体的涡流流态计算结果及非线性气动特性和截面压强分布。     

斜置平板的低速风洞实验研究

利用油流法和烟线法相结合的流动显示技术对不同斜掠角和迎角组合下的斜置平板进行了流场显示和和天平测力实验研究.油流显示结果揭示了斜置平板表面分离流动和再附流动的流动特征;烟线显示结果则揭示了斜置平板周围的涡系结构,从而解释了油流显示中所观察到的现象.通过天平测力,我们得到了斜置平板的滚转力矩特性和升力特性.实验结果表明:在小迎角下,随着模型斜掠角度的增加,一条展向发展的涡流动生成于模型的后掠段,该涡流动使滚转力矩保持为相对较小的正值;大迎角下,从斜置平板模型的前掠尖点生成了一股沿展向发展的涡流动,且随着模型斜掠角度的增加,涡流动的影响范围从模型的前掠段扩展到了后掠段,从而使滚转力矩先减小后增大.     

关于湍流问题的注记

指出N-S方程是研究湍流所遵循的基本方程, Reynolds方程反映了在湍流中随机性与确定性共存;Burgers-KdV方程与湍流机理无关, 不能作为湍流研究的“新方程”      

抛物化Navier-Stokes方程在求解高超声速流动稳定性分析 的基本流中的应用

通过在边界层内部使用边界层上缘的流向压力梯度代替该处的流向压力梯度,对传统抛物化Navier-Stokes(PNS)方程求解中的流向压力梯度的处理方法进行改进,从而使它可以求解用于稳定性分析的高超声速流动的基本流,并分别计算了平板、零攻角钝锥和小攻角钝锥三种典型算例.对于流向没有分离,且横向流动不强的流动,使用PNS方程计算高超声速流动稳定性分析用的基本流是可靠的.特别是改进后的方法求得的基本流的稳定性分析结果与直接求解Navier-Stokes(N-S)方程的结果非常吻合,但计算的时间消耗和空间占用都减少了一个数量级.     

Large-eddy simulation:Past,present and the future

Large-eddy simulation(LES) was originally proposed for simulating atmospheric flows in the 1960 s and has become one of the most promising and successful methodology for simulating turbulent flows with the improvement of computing power. It is now feasible to simulate complex engineering flows using LES. However, apart from the computing power, significant challenges still remain for LES to reach a level of maturity that brings this approach to the mainstream of engineering and industrial computations. This paper will describe briefly LES formalism first, present a quick glance at its history, review its current state focusing mainly on its applications in transitional flows and gas turbine combustor flows, discuss some major modelling and numerical challenges/issues that we are facing now and in the near future, and finish with the concluding remarks.     

高超声速等离子体流场混沌特性的实验分析

本文采用耗散系统的混沌动力学方法,分析处理了ＦＤ０４电弧风洞超高速等离子体流场实验中多道扫描探针的实验数据。分别给出了无锥体流场和有锥体尾迹等离子体流场的相关维数和Ｋｏｌｍｏｇｏｒｏｖ熵值的分布。该分布清晰地刻划了流场的湍动和分层结构。用Ｋｏｌｍｏｇｏｒｏｖ熵极小值方法,观测到了锥体尾迹流场的主激波结构。     

SST湍流模型在高超声速绕流中的改进

为模拟高超声速湍流问题,对剪切应力输运(SST)湍流模型系数进行了修正。数值格式采用改进的总变差递减(TVD)格式,并对湍流模型的负值强制项进行了隐式处理。在此基础上计算了绕平板以及具有分离、再附、激波/边界层干扰等复杂流动结构的压缩拐角的高超声速流动。计算结果与试验数据及半经验公式的对比表明:SST湍流模型引入的雷诺剪切应力与湍动能之比为常数(Bradshaw数)在高超声速绕流中并不成立。Bradshaw数修正后的SST湍流模型与原模型相比,所计算的壁面压力、摩擦阻力和壁面热流分布更接近试验结果。     

预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

Towards the use of large eddy simulation in engineering

This paper reviews some important aspects of Large Eddy Simulation (LES) as applied to engineering flows. We first summarize the present status of modeling in incompressible, compressible and reacting multi-phase flows, with a view towards the overall formalism instead of the intrinsic details of different subgrid models. On the basis of the assumed requirements on future LES, expected to handle full-scale flows and reacting flows with detailed chemistry, we discuss some potentially interesting LES methods for the future. These methods are exclusively based on multi-scale modeling, in which simplified equations are solved within each LES cell, instead of semi-empirical modeling based on the resolved flow scales only. After that we outline a few flows studied by LES at the Swedish Defense Research Agency—FOI, which form the basis for the subsequent discussion of validation and verification, and quality management, being of increasing importance for practical engineering flows. Next we summarize some practical aspects of LES of engineering applications, many of which being crucial to the successful use of LES, and being of increasing importance for engineering flows. Finally, a view of the future use of LES in engineering is presented, which is based primarily on the evolutionary use of LES during the last decade in the fields of hydrodynamics and combustion.     

一种基于动模态分解的翼型流动转捩预测新方法

考虑自由转捩的定常/非定常流动Navier-Stokes方程数值求解,对于翼型流动细节的精确模拟和气动力的精确预测均具有十分重要的意义。采用动模态分解(DMD)方法进行流动稳定性分析,再结合e N方法,提出了一套适用于翼型绕流的转捩预测新方法,称为DMD/e N方法。相比于传统的线性稳定性分析方法,DMD方法不需要求解附面层方程和线性稳定性方程,也没有引入平行流假设,具有更好的理论适用性和算法鲁棒性。开展了NLF0416、S809和SD7003等翼型的转捩预测数值验证研究,通过与实验结果以及与传统的基于线性稳定性分析的e N方法的比较,验证了本文所发展的转捩预测新方法在预测翼型的定常流动和非定常流动转捩方面的正确性,也表明了该方法具有解决含层流分离泡的翼型绕流转捩预测的能力。     

高超声速湍流直接数值模拟技术

概述了近年来国内外高超声速湍流直接数值模拟(DNS)技术方面最新的研究进展,主要集中在高精度、高鲁棒性数值方法方面,同时也介绍了近年来典型的高超声速湍流DNS算例。在数值方法方面,主要介绍了高精度激波捕捉格式以及保持计算稳定的数值技术,重点是WENO格式及高阶保单调格式的最新进展。在高超声速湍流DNS算例方面,介绍了压缩性影响、壁温影响、真实气体效应以及高超声速转捩等方面的DNS研究。此外,还简要介绍了作者开发的可压缩高精度计算流体力学软件OpenCFD。     

室内气流分布的大涡模拟

准确地模拟室内气流分布对优化室内通风设计及预测舒适和健康的室内环境有着重要的意义，传统的空内气流分布的数值模拟大都采用к-ε两方程湍流模型，该方法求解出的是流动变量的平均值，无法给出流场结构的详细信息，体现不出湍流流动的瞬时性的特点。本文采用大涡模拟（LES）法求解室内气流分布，使HSMAS法保证计算的稳定性和收敛性。通过对室内等温贴附射流和非等温受限射流两种类型的射流作用下的室内气流分布计算，以及导入к-ε两方程湍流模型对计算结果进行比较与分析表明，大涡模拟求解室内气流分布是一种确定可行和有效的方法。