间断Galerkin有限元隐式算法GPU并行化研究

为了提高间断伽辽金（discontinuous Galerkin, DG）有限元方法的计算效率,围绕求解Euler方程,构建了基于图形处理器（graphics processing unit, GPU）并行加速的隐式DG算法。算法结合Roe格式进行空间离散,采用人工黏性法处理激波等间断问题,时间推进选用下上对称高斯-赛德尔（lower-upper symmetric Gauss-Seidel, LU-SGS）隐式格式。为了克服传统隐式格式固有的数据关联依赖问题,借助于本文提出的面向任意网格的单元着色分组技术,先给出了LUSGS隐式格式的并行化改造,使得隐式时间推进能按颜色组别依次并行,由于同一颜色组内算法已不存在数据关联,可以据此实现并行化。在此基础上,再结合DG算法局部紧致等特点,基于统一计算设备架构（compute unified device architecture, CUDA）编程模型,设计了依据单元的核函数,并构建了对应的线程与数据结构,给出了DG有限元隐式GPU并行算法。最后,发展的算法通过了多个二维和三维典型流动算例考核与性能测试,展示出隐式算法GPU加速的效果,且获得...     

轴对称Boltzmann模型方程统一算法与空天飞行环境喷管流动模拟

拦截机动飞行器周围大范围区域存在主喷/侧喷流/羽流影响,而传统的纳维-斯托克斯（Navier-Stokes, N-S）方程不能很好模拟发动机喷管扩张段出口附近流动情况,需要一种新的方法来处理这种全流域流动问题。为解决该问题,针对特定轴对称喷管内流动,本文通过数学推导确立描述不同克努森数稀薄环境条件下的轴对称喷管内流动Boltzmann模型方程,初步建立适于该模型方程的数值格式与气体动理论统一算法。通过开展同轴圆筒间的定常/非定常旋转流动以及轴对称喷管内流动数值计算研究,发现统一算法计算流场与其他途径得到的结果吻合较好,验证了统一算法在全局克努森数喷管流动模拟的适应性和可靠性。通过与低密度风洞实验对比,喷管出口核心区羽流结构一致,羽流轴线压力分布一致,表明统一算法可以有效解决喷管入口压缩段到扩张段多流域混合,尤其是出口附近稀薄气体真空低压环境流动问题。     

吉赫横电磁室（GTEM）的波分析

吉赫横电磁室的主体是一种角锥状双导体同轴结构。为得到内部场分布，应对Ｈｅｌｍｈｏｌｔｚ方程求解。实现了求解过程，并应用了横向谐振绕射（ＴＲＤ）法进行处理。用自行编制的软件程序实现了场分布的三维图形显示，给出了高阶模截止频率的实算数据。     

报知通信系统抗干扰性研究

在野战地域通信网中,通信设备的通信能力和抗干扰能力至关重要,对战争的胜负具有决定性的作用。主要介绍了影响干扰能力的因素、干扰功率与干扰距离之间的关系,并结合EC 130H"罗盘呼叫"通信干扰飞机分析了通信设备的抗干扰能力。     

架空线自阻尼的测试方法

自阻尼是导线的基本特性之一。在高压架空导线防振设计中,自阻尼是一个基础数据。本文是介绍对导线自阻尼的三种测试方法。其中的“剪力法”是我们自己研创的新方法。文中不仅介绍了测量技术,还附有在我国第一次测出的 LGJQ-400型导线的自阻尼数据。     

吸气式内外流一体化飞行器动导数数值模拟

动导数是飞行器轨道及姿态控制系统设计时的重要参数,对飞行器开环系统受到扰动时振荡的敛散特性起重要作用。在标模验证的基础上开展吸气式内外流一体化飞行器WR-A动态特性分析,采用N-S方程模拟强迫简谐振动的非定常流场,获得了飞行器的直接阻尼导数、加速度导数和旋转导数,并对WR-A外形大振幅振动时的进气道性能参数进行了分析。研究表明,在某些情况下反映流动时滞效应的加速度导数占直接阻尼导数的比例最高可以达到40%以上,而在某些情况与直接阻尼导数符号相反。此外,在动导数研究基础上,采用参数化的动导数概念初步开展了非定常气动参数建模研究。针对WR-A大振幅强迫振动的研究显示,采用所建立的模型预测的气动参数与完全非定常计算吻合较好。     

基于时间约束的无人机避障研究

对无人机避开障碍物这一热点问题展开了研究。在极坐标系下,基于无人机与障碍物之间的几何关系,建立了无人机与障碍物之间的运动学方程。通过设计滑模变结构有限时间收敛制导律,使连接无人机与避障点的视线角速率快速收敛到零,相对速度方向收敛到期望的避障方向,保证了无人机能够顺利避开运动障碍物。通过有限时间收敛分析,得到了相对速度收敛到期望的避障方向时间与制导律参数的表达式。通过选择合适的参数,可使收敛时间小于到达避障点的时间,保证了避障的完成,也确定了制导律参数的取值范围。最后对设计的避障算法进行了仿真,仿真结果验证了算法的有效性。     

TA7钛合金高温流变行为研究

通过热压缩试验研究了TA7钛合金在变形温度850~1000℃、应变速率0.001~0.1s-1条件下的流变应力变化规律,计算并建立了描述TA7钛合金高温变形特性的本构方程。结果表明:变形温度和应变速率对TA7钛合金流变应力影响显著,随变形温度升高和变形速率的降低,相同变形程度下合金的流变应力显著降低,并且在较低的应变下合金即可达到稳态流变状态。     

基于代理模型的飞翼多目标气动优化设计

针对飞翼布局飞行器,采用雷诺平均N-S方程（ RANS）计算流场,使用基于代理模型的多目标优化方法进行了同时考虑起飞性能和巡航性能的多点多目标气动优化设计。在设计过程中,将飞翼的平面形状、剖面形状及扭转角同时作为设计变量（共58个设计变量）,将提高起飞时的升力系数和提高巡航升阻比为设计目标,以起飞状态和巡航状态的力矩系数作为气动约束,并以飞翼平面面积不减和剖面厚度不减作为几何约束。通过采用基于Kriging模型的多目标优化方法,以较小的计算花费得到了较好的Perato前沿。取Pareto前沿中一个最优解与基准外形的性能进行了对比,结果显示,优化外形的性能较基准外形的气动性能得到全面大幅提高且所有约束得到严格满足。     

波浪前缘静子叶片对高速轴流风扇单音噪声的影响

航空发动机压气机噪声与大型风洞压缩机噪声问题日益凸显,相关研究机构迫切需求新的降噪手段以指导大型叶轮机械降噪设计。为了探索波浪前缘静子叶片在大型叶轮机械降噪中的应用前景,采用非定常雷诺平均Navier-Stoke（URANS）方程与FW-H方程混合方法对基准静子叶片和3种波浪前缘静子叶片的降噪效果进行了数值模拟,研究对象静子来流平均马赫数约为0.49,基于静子叶片弦长的雷诺数约为1 040 000。数值预测结果显示：波浪前缘静子叶片可以显著降低高速轴流风扇单音噪声,但会对风扇的气动性能产生少许不利影响;相较于基准静子叶片,3种波浪前缘静子叶片可以在1BPF时降低风扇入口声功率级0.97~1.5 dB,2BPF时降低风扇入口声功率级2.89~4.9 dB,3BPF时降低风扇入口声功率级3.32~4.72 dB;同时,总压比降低0.1%~0.8%,等熵效率降低0.1%~0.3%。进一步研究表明：不同频率下声源振幅和相位关系是风扇单音噪声强度的主要影响因素,总的来说,幅值的增加会降低声源强度,然而通过改变声源相位关系的降噪方式则需要兼顾径向模态与波长两个方面。