大密度变化非连续介质气体流场的数值仿真

针对大密度变化非连续介质气体流场的仿真难题,实现了分区可变时间步长方法,为了保证宏观物理量通量在界面上连续,针对穿过两个时间步长不相等的区之间的界面的分子执行剩余时间缩放处理。提出了一种高效的分区界面和壁面处理算法,将每个界面设定为只属于一个区,将跨过多个区的壁面分解成分别只属于一个区的多个壁面,并规定仿真分子只能与本区的界面和壁面相互作用。对喷管流场的模拟表明,此方法可以大大提高计算效率。     

基于粘性笛卡尔网格的N-S方程计算

本文采用基于叉树数据结构的粘性笛卡尔网格求解N-S方程。为了使网格的分布更为合理,采用基于物体几何外形的方法对网格进行自适应加密。在物面处采用光顺、投影、分割加密的方法保证了物面处网格的正交性与贴体性,从而生成高质量粘性笛卡尔网格。将有限体积法应用于粘性笛卡尔网格求解N-S方程,采用了SA紊流模型。算例计算结果与实验数据符合良好。     

网格结构对Euler方程解的影响

 采用2种不同对接分区网格技术讨论了不同的网格结构对绕偏副翼的三维机翼流动的Euler方程解的影响.一种是部分不连续的网格(PDG);另一种是完全连续的网格(CCG),分别在PDG和CCG中采用Van Leer的矢通量分裂格式和Jameson的中心差分格式求解了Euler方程.算例表明,2种网格各有其优缺点,且都能得到有用的结果,但为了高效地取得更好的计算结果,必须根据流动和几何外形选取恰当的计算网格.     

桥梁断面的气动导数和颤振临界风速的数值计算

选取有代表性的大海带东桥、南京二桥流线型闭口箱梁断面和荆沙桥钝头开口板梁断面,用动网格法考虑流体和结构的耦合作用,用有限单元法求解原始变量二维不可压粘性流体的Ｎ－Ｓ方程,计算其气动导数,并用半逆解法确定颤振临界风速。计算结果和实验值相当吻合。     

高速旋转冲压增程弹用进气道复杂流场数值模拟

采用块结构网格与二阶精度流场分区求解技术,对高速旋转、含侧向支柱冲压增程弹丸进气道内外复杂流场进行了数值模拟,得到了高速旋转工况下对应于不同来流马赫数和攻角,以及临界工况时超音速进气道内外流场复杂的波系结构。在旋转工况下,进气道流场结构与不旋转时相似,但旋转速度以及攻角的存在对冲压发动机进气道的总体性能产生了负面影响,进气道总压恢复系数和动能系数均有所降低,而流场畸变指数则显著增大;冲压发动机进气道在较低马赫数工作时的综合性能优于在较高马赫数工作时的综合性能。     

基于k-ωSST模型的同心筒发射装置流场数值模拟

围绕同心筒发射装置流场数值模拟及结构优化展开研究,流场数值模拟基于轴对称Navier-Stokes(N-S)方程、k-ω shear stress transfer(SST)湍流模型,并运用域动分层动网格方法,以同心筒初始方案为参照,研究了同心筒发射过程的动态流场,在此基础上,分别对内筒收敛段的收敛角度以及筒口扩张段的扩张角度进行了优化研究,研究结果表明:内筒收敛段使得燃气排导更为通畅,筒内的热环境有所改善,并且内筒收敛段起到了一定的防喷作用,收敛角越大,筒底热环境越好;筒口扩张段对发射过程动态流场也有明显影响,收敛角为15°,扩张角为20°时,导弹从开始发射到完全出筒过程中热环境最为良好。      

基于Euler方程解摄动的超声速旋转导数计算方法

基于笛卡尔网格和采用Euler方程基本流场解,通过Riemann不变量摄动,对旋转扰动引起的物面法向速度构建摄动关系,发展了一种超声速旋转导数快速计算方法。文中给出了具体的推导过程,可以看出由此得到的超声速旋转导数,仅与定常欧拉解基本流相关。因此,一旦CFD计算得到定常欧拉解,就可一次性地快速计算出感兴趣的超声速旋转导数。为了验证该计算方法,本文选用国际动导数标准模型Basic Finner算例进行了计算,计算结果与试验数据、文献数据具有很好的一致性,展示出本文方法计算耗时少、精度高、工程实用性强等优点,且该算法基于笛卡尔网格,特别适合复杂气动外形飞行器超声速旋转导数的快速计算。     

采用改进的CE/SE方法模拟方管中氢氧爆轰波的稳定传播结构

采用改进的时空守恒元/解元（CE/SE）格式和新型二步化学反应模型，对小扰动情况下方管通道中的氢氧气相爆轰传播过程进行了数值模拟，对不同截面尺寸和两种扰动模式下的爆轰波三波线结构进行了详细计算和分析。数值模拟花费了相对较少的计算代价，清晰地得到了方形截面通道中3种相对稳定的爆轰波传播结构，即直角模式、对角模式和螺旋模式。随着截面尺寸的减小，爆轰波胞格图样与波阵面结构随之产生适应性的变形，乃至出现胞格数目的减少和部分三波线的消失；最终，直角模式和对角模式都会转化为顺时针或逆时针的螺旋模式；临界尺寸下，得到了转换过程中新的过渡结构和压力脉冲振幅的变化，并以此初步讨论了直角和对角传播模式的稳定性。     

格序群的极子群在凸l－子群格中生成的闭子格

设Ｇ是一个格序群，Ｐ（Ｇ）是Ｇ的极子群集合，Ｃ（Ｇ）是Ｇ的凸ｌ－子群集合，Ω（Ｇ）是Ｐ（Ｇ）在Ｃ（Ｇ）中生成的闭子格，本文通过对Ω（Ｇ）的研究得到了一系列重要结果，对有限值格序群给予了一系列新的刻画，证明了ｇ是基元当且仅当ｇ－是ｇ的唯一的Ω－值等结论，实际上对Ｐ．Ｃｏｎｒａｄ提出的问题＂Ｄｅｓｃｒｉｂｅｔｈｅｅｌｅｍｅｎｔｓｉｎ＂进行了富有意义的探索。     

基于SOVA的并行级联空时格码的译码算法

并行级联空时格码调制(Turbo—STCM)将并行迭代译码结构(Turbo)应用于空时编码，可以极大地提高空时格码的编码增益。本基于一种改进的(Soft—output Viterbi a1gorithm，SOVA)译码算法，提出了针对Turbo—STCM的并行迭代译码方案。仿真结果表明，在保持分集增益的前提下，Turbo—STCM和原始的空时格码调制相比，编码增益有了极大的提高．同时采用该SOVA算法译码复杂度低，速度快，更利于实际系统的开发实现．