基于非结构网格有限差分法的扎染算法

应用有限差分法时遵循贴体坐标系下离散等价方程的离散准则,格式中仅用到当地网格点的坐标变换系数,据此可以构建新的非结构网格有限差分法,数值算例表明在空间二维情况下用连接离散点的3条网格线构造出来的一阶迎风格式可以稳定收敛。把这种非结构网格有限差分法推广到三维,在4条网格线基础上进行离散计算,然后提出一种局部区域多重网格、重复计算的新算法。首先在包含物体的整个区域采用直角坐标系下的均匀网格进行计算,然后消除物体内部点对物体外部流场区域的影响,最后在物面和直角网格之间很小的局部区域填补非结构网格进行计算,计算过程类似于中国传统的扎染工艺。处理流场内存在曲面物体边界问题时,相较于常规结构网格差分算法,扎染算法尽管在计算中包含需要特殊处理的无用网格点,也存在重复计算的局部区域,但是以这些计算能力的浪费为代价,换取了编程简便、内存量小、网格生成快捷、易于扩展网格规模等优势。定性分析表明,这种扎染算法非常适合研制大规模并行计算,建立的计算流体力学应用软件可以充分发挥数十万至百万处理器核心的超级计算机的效能。     

多局部节点异步抗差航迹关联算法

为解决航迹异步与系统误差并存情况下的多局部节点航迹关联问题,提出一种基于区间序列离散度的多局部节点异步抗差航迹关联算法。定义区间型数据集的离散信息度量,给出系统误差下航迹序列区间化方法,通过累次积分计算离散度,结合多维分配进行关联判定。针对多局部节点上报目标不完全一致现象,设置零号航迹管理关联质量。与传统算法相比,无需时域配准,可在系统误差下对异步航迹直接关联。仿真结果表明,算法能在局部节点上报目标不完全一致场景下实现有效关联,且正确关联率随局部节点数目的增加或目标密集程度的增大而提高。     

基于离散伴随的流场反演在湍流模拟中的应用

精确模拟湍流流动是学术界和工业界均普遍关注的问题。采用数据驱动湍流建模的思路,建立了基于离散伴随方法的流场反演框架。通过为SA模型涡黏性输运方程的生成项乘以非均匀分布的系数,并利用有限的观测数据对该系数进行推断,实现对SA模型的修正。为了提高带有物理约束的离散伴随优化的效率,发展了约束增广的伴随方法,其高效性在本文得到了验证。选取了结冰翼型和周期山2个分离算例进行分析,所得结果在2个算例中均能以很高的精度拟合观测数据,并能借助湍流模型的修正将有限的观测信息泛化到整个流场。分析表明,流场反演所推断出的修正区域具有较为明确的物理意义,能够指导湍流模型的进一步改进。     

基于气体动理学方法的多尺度稳态热辐射计算

随着航空发动机燃烧室内燃气温度的提升,高温区域的热辐射问题将会变得更为重要,同时由于高温气体辐射导致辐射特性的剧烈变化,将会呈现出多尺度的现象,增加了辐射计算的难度。文章通过稳态离散统一气体动力学格式（SDUGKS,Steady Discrete Unified Gas Kinetics Scheme）的思路实现稳态辐射问题的求解,SDUGKS格式通过特征线差分离散的方法实现单元界面的重构,通过某种迭代格式来实现单元数据的更新,这一过程实现了对网格内部的辐射特性的有效模拟,该过程可适用于任何具体的辐射尺度。因此就可以在任意网格条件下实现多尺度问题的计算。文章引入了修正的延迟修正法（Deferred-correction,DC）实现单元数据更新,对单一尺度和多尺度问题进行了计算,在单一尺度问题中,验证了SDUGKS格式求解稳态辐射问题的正确性,进一步构造多尺度问题进行求解,论证了SDUGKS在多尺度问题计算中的渐近保持性质和准确性。     

基于DPM的试验台导流槽喷水冷却数值研究

采用计算流体动力学(CFD)技术,对某型运载火箭动力系统试验时导流槽的热环境进行了分析。计算采用氢氧(H-O)单步燃烧反应模型,标准k-ε湍流模型获得燃烧流场,并采用离散相模型(DPM)模拟了试验台的喷水降噪系统和导流槽底部喷水冷却系统。随后,用经过校验的模型对某型号发动机试验台导流槽的设计参数进行了对比研究。结果表明:试验台的喷水降噪系统及导流槽底部的喷水冷却系统均对导流槽底部有冷却效果,两者同时工作时对导流槽底部的冷却有叠加增强效果;导流槽深度加大时或者喷水量加大时会对底部热环境的改善同样有增强作用。      

基于离散数据的可变流量固冲发动机性能分析软件

应用Matlab的GUI平台开发了基于离散数据的固冲发动机性能分析软件PERFORMANCE SHOW,它以发动机性能计算获得的离散数据为输入,采用多维插值、多项式拟合、多项式求根、等值面计算、逻辑运算等方法,能够获得描述发动机工作特性的曲线、曲面和实体,从而可对发动机性能进行多维度深入分析。PERFORMANCE SHOW软件界面友好,使用简便,适合于各类冲压发动机的性能分析。     

离散式粗糙元诱导翼型边界层转捩的数值和实验研究

在Ma=0.6来流条件下,针对三维离散式圆柱粗糙元诱导NACA0012翼型边界层转捩问题开展了直接数值模拟(Direct numerical simulation,DNS)和油膜干涉风洞实验研究,分析了粗糙元诱导转捩的机理及粗糙元高度和相邻粗糙元间距对转捩位置的影响。结果表明,翼型表面粗糙元能够通过诱导出三维Λ涡和马蹄涡促进边界层转捩,达到转捩控制的效果。粗糙元高度和相邻粗糙元间距对边界层转捩有影响,且增加粗糙元高度和减小相邻粗糙元间距能够促进转捩。粗糙元高度对转捩的影响大于粗糙元间距,且对粗糙元后方区域影响大,对相邻粗糙元中间区域影响小。     

一种抵抗几何变换攻击的数字图像水印方法

在水印嵌入的过程中,首先在原图像经过对数极坐标变换后的图像中寻找水印嵌入的特定区域,并在该特定区域所映射的原图像经过离散傅立叶变换后的图像区域中随机选择像素点嵌入水印,保证了水印嵌入位置的不重复和不可侦测性。在水印提取的过程中,提出了获得图像旋转角度与缩放比例的方法,简化了水印的提取过程。实验证明,将对数极坐标变换和离散傅离叶变换结合后,在图像指定区域选择随机嵌入点嵌入的水印对于几何攻击具有良好的鲁棒性。     

低速长航时无人机垂直突风过载分析

低速中空长航时无人机的升力非线性较强,利用标准离散突风模型和连续紊流模型计算得到的突风过载差别较大.针对该问题,本文采用改进的离散突风模型求解突风载荷,即直接求解突风条件下原始的运动学方程,并在运动历程中考虑非线性升力的影响.结果表明:改进的离散突风模型计算得到的突风载荷明显低于标准离散突风模型所得的载荷,更贴近无人机的实际情况,且改进的离散突风模型与连续紊流模型的计算结果也相近.