高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

开放式气动力数值模拟系统研究

气动力数值模拟系统是ＣＦＤ流场解算技术、网格生成技术、数据可视化技术和网络技术相结合的产物。描述了一种运行于计算机网格环境下多任务、多用户、交互式的数值模拟软件环境－ 气动力数值模拟系统,该环境能计算模拟亚、跨、超音速绕流的气动力,也能方便地集成模拟不同复杂气动布局飞机流场的解算程序。介绍了其网格生成、数据可视化、系统接口等核心部分的结构和功能。     

前体边界层状态对高超声速进气道流动结构及性能的影响

采用CFD方法研究了前体边界层状态对轴对称高超声速进气道的波系结构、激波/边界层干扰特性以及流量捕获能力、总压恢复系数等的影响,获得了不同马赫数下进气道性能参数随一级锥边界层强制转捩位置的变化规律,并对采用二维粗糙体作为强制转捩措施给流场所带来的干扰进行了分析,研究中还对比了具体进气道设计方案间的差别.      

NNW-HyFLOW高超声速流动模拟软件框架设计

NNW-HyFLOW在国家数值风洞（NNW）工程支持下,由风雷开源软件提供基础框架支撑,将建设成为基于结构/非结构混合网格,面向高超声速应用领域的国产自主工业CFD软件,具备高温气体热化学非平衡效应模拟及其相关气动力、气动热和气动物理特性计算分析等主要特色功能。本文从框架设计、数据结构、耦合方法、并行计算方法以及接口设计等方面对软件设计思想和框架特点进行了介绍,给出了求解器采用的理论模型、核心数值算法及其实现方法,结合HEG风洞试验、RAM-C飞行试验、Electre飞行试验以及航天飞机OV102飞行试验等典型算例开展了数值模拟。研究表明：NNW-HyFLOW具有底层代码复用、功能兼容性好、扩展能力强和接口灵活等优点,其当前测试版本已经具备了较好的高超声速非平衡流动数值模拟能力,在热化学非平衡效应及其影响的气动力特性、气动热环境和等离子体分布特性预测与评估方面,具有较高的数值计算精准度,初步满足了高超声速复杂飞行器高温非平衡流动数值模拟的需求。     

多斜孔壁冷却方式小孔内对流换热的数值模拟

采用数值计算的方法研究了多斜孔壁冷却方式小孔内对流换热的特性.数值分析采用了只考虑流体区域对流换热和考虑流体与固体间流固耦合传热共两种方法.计算结果表明:流固耦合传热法更能反映小孔内实际的换热特性;小孔进口处的Nu数最大,且随着孔内雷诺数的变化,孔进口区对孔内换热的影响也不同.      

Pattern control for large-scale spacecraft swarms in elliptic orbits via density fields

Space swarms, enabled by the miniaturization of spacecraft, have the potential capability to lower costs, increase efficiencies, and broaden the horizons of space missions. The formation control problem of large-scale spacecraft swarms flying around an elliptic orbit is considered. The objective is to drive the entire formation to produce a specified spatial pattern. The relative motion between agents becomes complicated as the number of agents increases. Hence, a density-based method is adopted...     

基于广义Richardson外插方法的阻力预测精度分析

利用计算流体力学(CFD)方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA 0012翼型无黏绕流和平板湍流边界层两个算例中验证了NSAWET程序和广义Richardson外插方法,分别得到数值算法预测压差阻力和摩擦阻力能达到的名义精度。进而模拟三维通用研究模型(CRM)翼身组合体绕流,得到的阻力名义精确值在DPW 5的统计误差带范围之内;综合DPW 5的计算结果来看,不同CFD解算器的结果之间存在一定差别,阻力预测精度总体上不符合二阶。可见,标准Richardson方法采用的二阶精度假设难以普遍适用,有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度,采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性,可以进一步降低网格对阻力预测的影响。     

高超声速飞行器表面横缝旋涡结构及气动热环境数值模拟

针对高超声速飞行器表面横缝内部流动,通过求解可压缩Navier-Stokes方程,自主研发了一套能够较好模拟缝隙流动特性的CFD软件,利用该软件对横缝进行了数值模拟。研究表明,对于缝隙内黏性起主导作用的低速流动,通常需要密度足够充分的网格才能获得较好的计算结果。收敛后的壁面热流与文献实验热流吻合较好;缝隙内部旋涡个数近似与缝隙深宽比成比例,与文献中的有关结论完全一致。缝隙壁面热流分布受旋涡结构影响,会出现波浪式变化。因此,合理捕捉旋涡结构对于缝隙壁面热流的数值模拟有着重要意义,本文CFD软件对高超声速飞行器表面横缝的计算结果具有较高可信度。     

再入飞行器高超声速绕流的一步和二步有限元计算

本文针对再入飞行器的空气动力问题，采用一步和二步有限元格式求解可缩守恒型Ｅｕｌｅｒ方程，并结合不同的人工粘性模型来捕捉激波，给出了高超声速二维钝头绕流的结果。     

From geometry to CFD grids—An automated approach for conceptual design

The CEASIOM software developed in the EU-funded collaborative research project SimSAC generates stability and control data for preliminary aircraft design using different methods of varying fidelity. In order to obtain the aerodynamic derivatives by CFD, the aircraft geometry must be defined, computational meshes generated, and numerical parameters set for the flow solvers. An approach to automation of the process is discussed, involving geometry generation and mesh generation for inviscid as well as RANS flow models.     

气动热数值模拟中的网格相关性及收敛

针对气动热数值模拟中的网格相关性和收敛问题,采用计算流体力学(CFD)计算方法,以典型的钝头为算例,进行了数值模拟研究。研究结果表明:网格是气动热数值模拟中的关键因素,壁面附近法向网格间距最为敏感,热流结果随网格间距不同会出现数倍乃至数量级的差异,在得出正确的壁面压力结果的网格上不能保证得到正确的热流结果。分析认为壁面附近客观存在的巨大温度梯度是网格要求敏感的原因。计算表明:气动热计算的收敛过程比压力和流场收敛慢得多;在压力和流场完全收敛时,热流与最终的迭代收敛值偏差达20%以上;以方程残值的下降量级和压力收敛作为热流收敛的判别标准是不合适的。在气动热计算中,应直接观察热流数据的收敛,确保得到真正的收敛解。     

Some vortical-flow flight experiments on slender aircraft that impacted the advancement of aeronautics

This paper highlights the three aerodynamic pillars of aeronautics; namely, theory/CFD, wind-tunnel experiments and flight tests, and notes that at any given time these three are not necessarily at the same level of maturity. After an initial history of these three pillars, the focus narrows to a brief history of some vortical-flow flight experiments on slender aircraft that have impacted the advancement of aeronautics in recent decades. They include the F-106, Concorde, SR-71, light-weight fighters (F-16, F/A-18), and F-16XL. These aircraft share in common the utilization of vortical flow and have flown at transonic speeds during a part of the flight envelope. Due to the vast amount of information from flight and CFD that has recently become available for the F-16XL, this aircraft is highlighted and its results detailed. Lastly, it is interesting to note that, though complicated, vortical flows over the F-16XL aircraft at subsonic speeds can be reliably and generally well-predicted with the current CFD flow solvers. However, these solvers still have some problems in matching flight pressure data at transonic speeds. That this problem has been highlighted is both an advancement in aeronautics and a tempting prize to those who would seek its solution.     

High-enthalpy hypersonic project at ONERA

This article deals with the research performed in the framework of the internal ONERA Hyperenthalpic Hypersonic Project. This project involved fifteen researchers specialised in different domains, which are the following: fundamental studies on shock/wave and shock/boundary layers interactions; study of the laminar/turbulent transition; real gas effects and associated physical modelling; the development of non-intrusive optical diagnostic methods usable in cold and hot wind tunnels; flow rebuilding in a hot shot wind tunnel such as F4; conception, realisation and validation of computational fluid dynamics (CFD) solvers.     

超声速粘性剪切流空间稳定性对称紧致差分数值分析

用一类高精度对称紧致差分格式数值离散一阶改型三维可压粘性扰动方程,对导出的非线性离散特征值问题采用二阶修正Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ边值迭代局部解法,实现了超声速剪切流的线性空间稳定性分析。基于对流Ｍ数、Ｒｅ数、速度比、密度比等参数研究,讨论了压缩性效应、粘性效应、超声速扰动快／慢模态等,结果显示超声速粘性剪切流的弱不稳定性同多种影响因素密切相关。     

高超声速飞行器热环境与结构传热的多场耦合数值研究

为了准确预测高超声速飞行器面临的严峻气动热/力环境以及结构的热力响应,发展了高超声速流动与结构传热耦合框架。采用分区求解方法,通过耦合界面的实时数据传递,实现了基于Navier-Stokes方程的高超声速化学非平衡计算流体力学（CFD）求解器与结构的热力全耦合有限元法（FEM）求解器的多场耦合计算,建立了高超声速飞行器的多场耦合数值分析方法。首先对经典高超声速圆柱绕流实验进行了耦合计算,结果与实验值吻合良好。然后针对典型的超高温陶瓷（UHTC）材料的耦合传热问题进行了数值研究,考虑热传导效应对气动热环境和结构热响应预测的影响,结果表明对于复杂外形且热导率相对较高的UHTC材料,结构内部热传导对热环境和表面温度分布的影响不可忽略。最后针对UHTC材料热物性（比热和热导率）非线性对高超声速流动传热过程的影响进行了研究,结果表明当比热和热导率处于合理的误差范围内时,材料表面温度响应对其变化并不敏感。     

雷诺数对涡轮叶尖流场影响的数值研究

采用基于压力修正的三维计算流体力学程序，结合雷诺应力湍流模型和剪切应力传输湍流模型加壁面函数的方法，对某一轴流涡轮转子叶尖间隙流场进行了数值计算研究，详细计算了不同涡轮叶尖间隙高度和来流湍流度条件下雷诺数对涡轮转子间隙流场的影响，最后计算了转子效率。结果表明：当泄漏流流经叶尖时因为叶尖剪切力做功有块总压增大区；雷诺数带来的影响比湍流度和叶尖间隙高度带来的影响要大，湍流度的变化对流场影响不大；雷诺数对泄漏涡尺寸的影响不大，但低雷诺数会引起流动分离，带来损失，当雷诺数在文中的范围内变化时，效率下降近八个百分点。     

高超声速升力体气动力气动热数值计算

升力体气动布局，由于其良好的高超声速气动特性、有效的内部空间和有益的防热特性，越来越受到气动工作者的重视，本文应用张涵信院士的NND格式思想，用数值模拟的方法对类似于美国X－33的升力体外形进行了气动力、气动热研究。特别是在粘性项的离散和表面热流率的计算中应用了散度和梯度的积分定义，避免了数值计算的奇异性，保证了通量守恒，提高了计算效率。结果表明气动热计算结果与美国同类飞行器结果接近，驻点热流率与经典理论结果一致。     

轨迹优化的直接数值解法综述

对飞行器轨迹优化的数值解法进行了系统的归纳与分析。从不同角度对数值解法的分支——直接法进行了概述;重点阐述了直接配点法的新进展——伪谱离散化方法;最后,介绍了相应的非线性规划求解器和实际应用中的多种轨迹优化软件包。     

高性能计算支持的静气弹优化设计方法

为扩大国产超级计算机在国内航空航天工业领域的应用,面向大展弦比机翼气动结构综合优化设计的实际问题,基于具有大量工业应用和自主知识产权的CFD软件TRIP和CSD软件SiPESC,应用改进的三维类别/形状变换（CST）参数化方法、差分敏度并行计算策略和移动渐近线（MMA）等大规模优化算法,通过TRIP+SiPESC耦合气弹分析软件预测气动载荷,发展了一套机翼型架外形的优化设计方法。以CRM标模内插机翼为例,针对无弯扭的初始型架构型,以气动外形展向几何扭转分布和结构有限元模型一万个壳单元厚度为设计变量,实现1 g载荷气动优化和2.5 g载荷结构优化的综合;在不降升阻比、满足最大应力/位移的约束下,俯仰力矩系数改善近40%,针对本文假设材料和初定的尺寸参数,机翼减重近9%,验证了方法和软件的有效性。该方法直接采用产品级高保真的自主静气弹分析软件预测气动力/载荷,将计算资源耗费大的静气弹优化和设计变量/约束规模大的结构优化分治,有较好工程应用前景。     

三维离散矢量场中的涡核识别

基于临界点理论，使用特征矢量方法在三维离散矢量场中进行了涡核寻找的研究。该方法能够快速准确地找到三维矢量场中的涡核，为CFD工作提供了一个有效的数值分析工具，搜寻结果与其它方法（如流线法等）进行了比较。