圆锥低超声速有攻角绕流非对称涡流的数值模拟与非线性分岔的研究

本文用数值方法研究了圆锥低超声速有攻角绕流的对称和非对称定常解，扰动响应以及在更大角时出现的准周期解问题。     

非对称体超音速无粘绕流流场数值计算

 本文给出了非对称再入飞行器的超音速三维无粘流场的数值计算方法。它归结为求解三个自变量—阶定常拟线性双曲型偏微分方程组的初边值问题。采用MacCormack二步显式差分格式。 为了说明计算方法的实际能力,对弯头钝锥和带控制翼的再入体作了计算。其精确度是令人满意的。     

弹道再入飞行器的无粘绕流及表面加热率的一种工程数值解法

RVSPHR是一种确定弹道再入飞行器无粘气动力和表面热流率的近似数值计算程序。在头部亚跨音速区内,由工程方法给出激波角及压力分布,并假定激波层内的物面法向压力和法向速度剖面,而相应的激波形状则按质量守恒条件迭代调整。下游超音速区仍由严格运动方程差分数值求解。在无粘流计算的基础上,用无需迭代的熵吞技术,计算表面热流率。所有计算结果与严格数值解及实验值作了比较,其一致性较好。     

几何微扰动对脱体涡非对称性影响的数值分析

利用拟压缩性方法，Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因式分解格式数值求解三维定常不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程采用Ｄｅｇａｎｉ－Ｓｃｈｉｆｆ修正的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，计算了绕尖头正切拱型旋成柱体的大迎角大雷诺数脱体涡流场，计算结果表明，在柱体头部施加微小的非对称几何扰动时，可以得到与实验相符的非对称脱体涡。     

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

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

球锥-柱-裙组合体超音速无粘流场的数值计算

本文叙述组合体超音速无粘流场的一种数值计算方法,跨越流场内部的间断面采用MacCormack二步显式差分格式(下称内间断“捕捉法”),在物面边界上混合使用两种计算方法。文中还论述了MacCormack二步格式的“异步”程序设计,对于低存贮量的计算机,可成倍提高计算速度。本文方法适用于钝锥、钝多锥,幂次体及钝锥-柱-裙组合体等轴对称物体,也可推广于身部非对称体的计算。大量计算表明,该方法简便、省机时、结果可靠、适用性强。     

对称驱动方腔定常流动的多解现象研究

对方腔内顶盖正弦速度驱动的定常流动进行了数值研究.采用非线性动力系统数值分岔分析方法,讨论了Re<2000时Navier-Stokes方程定常解的多解共存现象,计算了对称和非对称定常解的分岔曲线和相应的流动模式,给出了非对称解出现的临界参数Re=320,稳定的对称和非对称解同已有文献的结果一致,并获得了一组新的不稳定的非对称解.随着Re数的增加,侧壁出现流动分离并有小涡出现,随之发生流场拓扑结构分岔,但并不伴随N-S方程解的分岔.     

机翼大攻角气动力数值计算

本文描述了一个计算任意平面形状的机翼在无粘、亚音速、大攻角情况下的气动力数值方法,计算大攻角下气流分离时的定常非线性气动载荷。考虑翼面具有脱体涡分离,对旋涡强度和自由涡的位置进行松驰透代,取得收敛解。计算结果与实验数据很接近。     

细长锥体有攻角绕流的非对称分岔及其涡结构的数值模拟

本文用ＮＮＤ格式计算了具有复杂涡结构的细长锥体有攻角绕流问题，着重研究了非对称分离的分岔现象以及分岔前后涡结构的演变，数值计算结果表明，大约在ａ／（２β）＝１．０时，Ｎ－Ｓ方程出现分岔解，同时还给出了分岔前后涡结构及物理最分布的特征。     

航天飞机再入大气轨道分析

本文根据航天飞机再入大气轨道的运动方程,求出了升阻比为常值条件下的数值解。着重研究了改变升阻比和再入角对再入轨道的影响。并将精确数值解与近似解进行比较。此外,文中还估算了作用于航天飞机的最大热流。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

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

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

无人直升机组合外挂体气动特性研究

某型无人直升机机身左右两侧的组合外挂体由挂梁、挂架、导发架和导弹四部分组成。采用CFD计算方法对组合外挂体的气动特性进行了数值计算,计算状态包括不同的导弹发射安装角和机身侧滑角以及无挂载状态,计算过程中考虑了机身、主桨毂、尾梁和平垂尾等部件的气动干扰。结果表明导弹对组合外挂体气动特性的影响很大,组合外挂体的升力、俯仰和滚转力矩随导弹发射安装角的改变呈线性变化趋势;无人直升机转弯飞行或遇到侧风以及挂载导弹时,组合外挂体对无人直升机的气动特性影响更大。研究结果可为无人直升机选择合适的武器外挂形式与安装角以及估算外挂组合体的气动特性提供有用的参考。     

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

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

飞机亚音速非对称运动流场数值计算

本文以飞机亚音速对称运动（仅小攻角）流场数值计算为基础。采用涡格法理论开发飞机亚音速非对称运动（小攻角、小侧滑角、小角速度）流场特性数值计算。本文特点：从毕－萨公式出发推导出发推导的涡格法流场诱速公式比面涡法简便；采用多坐标系，详细推导了各坐标系转换公式；将复杂的干扰流场简化为一阶干扰和二阶干扰效应，并仔细推导了相应的流场诱速计算公式。已在ＩＢＭ－４３４１计算机上建立了流场计算实用性程序，既可单独     

飞机稳态尾旋平衡点的综合求解方法研究

为求解飞机稳态尾旋平衡点,推导了飞机稳态尾旋的动力学方程,并针对该方程给出了三种求解方法:作图法、数值求解法以及在前两种方法的基础上所提出的一种综合求解方法。综合求解方法首先使用作图法得出方程的所有近似解,并用稳态尾旋判据判断解是否为稳定解,然后将稳定近似解作为数值求解法的初值进行精确求解。运用综合求解方法对算例飞机进行了稳态尾旋平衡点的计算,算法收敛迅速,结果精确,表明该方法具有很好的工程应用价值。     

绕低频振动旋成体的跨音速非定常流核函数方法

本文对无粘跨音速流中作俯仰振动的细长旋成体,给出了一种计算表面非定常压力分布和动稳定性导数的近似方法,即用时间线化近似,在引入局部线化假设后,将跨音速小扰动方程化为一个积分方程,并用精确物面边界条件求数值解。本方法算例与其它理论和实验结果作了比较,表明用它可以得出较合理的结果。     

高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。     

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。