带制退器的膛口流场数值模拟

带制退器的膛口流场是一种具有高度瞬变性、强烈方向性的复杂波系结构流场.采用分块网格划分的整体运动处理方法和基于动态层变方法的结构动网格技术处理弹丸高速运动造成的网格变化,模拟某型带炮口制退器的火炮发射过程中的非定常燃气射流流场,根据数值结果,详细讨论了初始流场、火药燃气流场的形成与发展以及与弹丸的耦合和相互作用过程.数值计算结果可以为炮口制退器的设计与实验提供一定理论指导,为深入研究奠定良好基础.     

非结构动网格分区并行计算方法

针对包含运动边界的非定常流动问题,基于弹簧近似法设计了一套非结构动网格分区并行计算流程,引入了METIS软件进行非结构网格二次剖分,研究了对接边界通信类型及其处理方法,基于消息传递并行编程模式实现了非结构动网格求解器的并行化。测试结果表明:该求解器具有较高的并行效率,网格变形模块对非结构动网格并行计算性能具有重要影响。     

基于新型非结构有限体积解耦算法的激波诱导燃烧数值模拟

把基于结构网格有限差分方法建立的化学非平衡流动模拟的新型解耦算法推广到基于非结构网格的有限体积法,从而可用于模拟复杂几何构型下的化学反应流动.对H2/ Air预混气体中激波诱导振荡燃烧的Lehr实验进行了数值模拟,计算得到的振荡频率与实验结果符合很好,表明计算方法具有时间和空间二阶精度.通过对不同几何外形发射体的计算发...     

运载火箭芯级与助推级分离的数值模拟

采用动态重叠网格方法和双欧法,耦合求解非定常流动控制NS方程和6DOF运动方程,数值模拟运载火箭助推级和芯级的自由分离过程和附加外力作用下的强迫分离过程,并对芯级、助推级分离的主要影响因素进行分析。数值结果与试验结果吻合,表明动态重叠网格方法能正确预测火箭芯级与助推级间复杂的流动现象。根据数值模拟结果,给出助推级在不同分离工况下的分离参数及分离规律。     

掠飞末敏弹耦合运动动导数计算

以新型掠飞末敏弹为研究对象,提出了基于欧拉转动定理和滑移网格技术的复杂角运动模拟方法,利用著名的罗德里格斯转换矩阵插值求得弹箭在每个时间步的角速度修正值,并指定给球形滑移网格区。通过对非定常气动参数进行求解辨识,分析了不同马赫数下掠飞末敏弹滚转运动对其俯仰组合动导数和升力系数动导数的影响规律。结果表明:所提角运动模拟方法可有效消除姿态角计算的累积误差,实现对弹箭任意给定角运动的准确模拟;弹箭滚转运动对俯仰组合导数和升力系数动导数的辨识结果均存在显著影响,在进行弹箭动导数计算和稳定性分析时需充分考虑俯仰耦合效应的影响。     

激波诱导异质气体界面失稳的数值模拟

基于二维N-S方程,利用有限差分数值离散方法,对激波诱导异质气体界面失稳的现象进行了数值模拟,与文献中实验结果和计算结果进行了定性比较,并进一步分析了整个流动的非定常动态变化特性和非线性特征.研究表明,本文数值模拟的非定常流场图谱与文献中的实验结果和数值结果吻合较好;数值结果捕捉到了六氟化硫界面的演变过程及流场中复杂的波系结构.     

平面运动激波绕障碍物流动的非结构Euler解

用三角形非结构网格及有限体积法求解Euler方程 ,对运动激波跨越不同截面形状的二维障碍物绕流进行了数值模拟。计算结果给出了不同时刻各流场的复杂结构 ,该结构与实验结果一致。计算表明用非结构算法模拟复杂外形绕流具有简单、准确的特点 ,并能避免尖角处非物理结果的产生     

高超声速飞行器前缘流-热-固一体化计算

针对高超声速流动气动加热与结构传热的复杂耦合问题,探索和研究基于有限体积法的高超声速流-热-固一体化求解方法,将流场与结构温度场进行统一建模与数值模拟。该方法避开了传统气动加热和结构传热耦合求解方法在时间域内进行流场与结构温度场耦合交替迭代计算所带来的大量数据交换与计算,将流场与结构温度场作为一个物理场,采用统一的控制方程进行求解。采用典型高超声速绕流二维圆管稳态或非稳态流-热-固耦合算例对该一体化方法进行验证,稳态时圆管驻点温度最高达到648 K,非稳态时的热流密度和结构温度与参考文献和实验值吻合较好,由此证明了该方法的可靠性和正确性。与耦合计算方法的对比分析结果表明:该一体化求解方法所得计算结果更接近实验值,并且计算量和网格依赖性都相对较小,具有更好的稳定性和计算精度,能为高超声速飞行器一体化热防护设计提供有效的理论和技术支撑。     

潜艇应急上浮六自由度运动及黏性流场数值模拟

为揭示潜艇应急上浮过程中的强非线性运动及水动力作用规律,基于RANS方程及流体体积模型,针对Suboff标模建立了潜艇应急上浮数值预报方法。通过开展不同纵倾攻角及斜航漂角下全附体模型的水动力计算,验证了该计算方法的有效性,确定了其适用范围,进而结合整体动网格技术模拟了潜艇的应急上浮过程,获得了艇体六自由度运动的时历参数及流场细节信息。数值模拟结果表明该方法能够较合理地描述潜艇上浮的运动规律,也证明了该方法在潜艇上浮多自由度运动及水动力性能研究中的潜力和适用性。     

机体/推进一体化飞行器动导数频域计算方法

运用Etkin非定常气动力模型,将谐波平衡法应用于复杂外形动导数辨识。计算结果表明:谐波平衡法和时域方法的动导数辨识结果一致,与实验值吻合,验证了程序的正确性和谐波平衡法的可靠性。采用谐波平衡法对类X-51高超声速机体/推进一体化飞行器WR-A进行动导数辨识,并与时域方法的计算结果比较,考察了该方法在复杂外形下数值模拟非定常流动和动导数辨识的能力。实验结果表明:对于WR-A这种复杂外形飞行器,谐波平衡法也能在保证精度的同时具有较高的计算效率。     

钝体高超声速气动加热与结构热传递耦合的数值计算

气动加热与结构热传递耦合问题在航天和工程应用领域非常重要。分别采用松耦合与紧耦合方法,数值模拟了高超声速二维圆管绕流的流场与结构传热耦合的非定常过程。在紧耦合方法中,流场部分采用基于Navier Stokes方程的有限体积法,将AUSM+格式与时间方向的显式多步Runge Kutta法结合;结构传热部分采用基于二维热传导方程的Galerkin有限元法。流场与结构区通过交界面的热流和温度边界条件实现耦合。计算结果分别与实验、文献做了对比,结构内部温度变化关系以及壁面的热流分布均较好地吻合。两种耦合方式的计算结果对比表明,对于流场特征时间远小于结构传热特征时间的问题,松耦合方法计算效率高,精度与紧耦合方法接近。     

高效计算网格壁面距离的KD树方法

计算流体动力学的一些领域中经常要用到流场点到壁面的最小距离,为了提高壁面距离计算的精度和效率,将KD树应用于非结构网格的壁面距离计算,计算了若干三维外形的壁面距离,结果表明:基于KD树的网格壁面距离计算方法鲁棒性好,计算效率和结果精度高,适应复杂外形的能力强;算法的通用性好,可以拓展应用到多种网格类型。     

管道气液固三相流的耦合瞬变模型研究

在管输水力瞬变过程中,管道与流体之间存在耦合互动作用,需要用流固耦合模型来进行描述.针对工程中的多相流动问题,以最常见的气液固三相流为对象,建立了考虑流固耦合效应的水力瞬变模型.通过改变流体中固相或气相的体积分数,模型可适用于气液、固液两相流以及纯液体单相流的耦合水力瞬变,是单相流耦合瞬变模型向多相流耦合瞬变模型的拓展...     

非定常Euler方程数值计算中高精度格式应用

采用五阶精度加权紧致非线性格式(WCNS)和非定常“双时间步”方法求解非定常Euler方程,模拟NACA0012翼型强迫俯仰振动流场,研究了高精度格式应用到非定常计算时“双时间步”方法物理时间步长、子迭代收敛判据、子迭代步数以及物理时间导数离散方法对计算精度和计算效率的影响.     

基于OpenFOAM的某型炮弹初始段气动仿真研究

为研究某型炮弹初始段的气动特性并获得相关气动参数,运用结构化网格划分方法建立该型炮弹初始段仿真模型;依托OpenFOAM软件平台,应用可压缩流动求解器对初始段气动特性进行仿真,并提取可用于外弹道设计中的相关气动参数;将仿真结果与风洞试验值进行对比,发现可压缩流动求解器的仿真结果具有很高的精度,完全满足该型炮弹高精度射表的编制要求．     

基于单元变形的网格间流场信息传递方法

针对单元网格发生变形重构中新旧网格之间的数据插值问题,提出一种基于格点格式有限体积法的流场数据传递方法。利用非结构动网格技术,将旧网格单元移动到新网格单元,同时时间推进求解流场控制方程,获得移动后旧网格单元物理量,并将其赋给新网格单元,以此来实现两套网格间的信息传递。计算结果表明,该方法在信息传递过程中没有引入插值误差,理论上流场求解方法的精度就是信息传递精度,验证结果表明其显著优于二阶插值精度。     

液管道系统动力响应的研究进展

综述了自70年代以来关于输液管道系统的动态响应问题的研究进展情况.根据管路结构的特点,将问题分为两端支撑直管、悬臂管、不稳定流直管、圆柱壳管和曲管5个部分分别进行了阐述.并在此基础上,概括了应用于该领域的一些计算方法.最后,在总结研究成果的同时分析了这一领域的研究趋向.     

凝结水泵三维流场数值模拟与分析

为深入理解和分析电站凝结水泵性能,采用计算流体力学方法(CFD)对其内部复杂的三维流场进行了数值模拟与分析.首先,建立了某三级凝结水泵各部件几何模型,并对其进行网格划分.通过求解RANS方程对凝结水泵内部三维流场进行数值求解,湍流流动采用SST模型进行模拟.计算得到该泵若干工况下的扬程和功率等性能参数和泵内部流动特征,计算结果与试验数据吻合很好.模拟结果及分析表明:采用CFD方法模拟结构复杂的凝结水泵内部流动和进行性能预报是可行的.通过内流场分析还发现导叶体等部件的设计有可改进之处.     

水下应用栅格翼动态展开参数预示方法

基于水下应用栅格翼动态展开过程动力学模型,根据展开特征角度定常水动力方法,拟合获得展开全程受到的流体力矩,并引入考虑相对速度影响的修正因子,形成水下应用栅格翼动态展开过程参数工程预示方法,以典型展开时序点航行体运动参数为设计输入,对栅格翼展开过程运动参数进行预示。通过与水下应用栅格翼非定常流场仿真计算数据以及水下航行体弹射试验数据对比,验证了上述预示方法的正确性及工程适用性,为水下应用栅格翼方案设计的优化及展开不同步性分析提供设计参考。