基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

一种新型的直角网格生成方法及应用研究

提出了一种基于背景网格的直角网格新型生成方法。背景网格是用来提供网格加密的信息．通过网格尺度与所得计算参数的比较，定量判断确定网格是否需要加密；而在求解Euler方程计算中，采用的是Jameson的有限体积法，它可以适合任意形状的网格单元，以及四步Runge—Kutta时间推进。本文对NACA0012翼型和复杂多段翼型等问题进行了数值实验。结果表明，这种网格生成方法易于推广到三维情况中去，具有网格生成时间短．收敛速度快，易于处理复杂边界的优点。     

适用于气动学问题的DRP类有限差分方法

研究了一类低耗散低色散的高阶精度有限差分方法,目的是直接计算非定常欧拉方程用于气动声学问题.采用的数值方法是空间四阶、时间三阶精度的色散关系保持(DRP:dispersion-relation-preserving)类有限差分格式,通过波动方程算例验证了格式模拟波动问题的能力;采用加变幅值的高波数人工耗散项来抑制高频数值振荡,使得该格式可推广用于含激波、初始间断或非线性波动的问题.通过含激波、初始间断条件的算例验证了所用变幅值人工耗散的有效性,最后作为初步应用计算了二维亚、超音速均匀流中简单点声源辐射问题,得到了很好的结果,表明此类方法在气动声学问题计算上深有潜力.     

地下电缆传热分析及其内部温度场数值计算

高压电力输送在城市一般都采用地下电缆的形式.为确保电缆能够安全正常的工作,采用精确方法预计地下电缆的载流量和热特性是很有意义的.采用有限差分法,以单芯电缆为例,进行传热分析,最终可以确定电缆允许的载流量.计算与实验结果表明,完全可以满足工程需要.     

一种降低数字式FM/FM解调系统分路滤波器阶数的方法

FM/ FM遥测体制有着广泛的应用 ,研究数字式 FM/ FM解调系统具有重要意义。分路滤波器 (通常是带通滤波器 )是解调系统的重要组成部分 ,有限冲激响应 (FIR)数字滤波器具有平坦的通带幅度响应和线性相位特性 ,是数字式解调系统的最佳选择。但是 ,符合 IRIG标准的分路滤波器通常是窄带滤波器 ,需要很高的阶数 ,这给系统的实现带来了很大的困难。文中根据变采样率多级实现 FIR滤波器的原理和带通信号采样定理 ,并结合 FM/ FM系统 IRIG标准中各路信号的频谱特性提出了先按频谱分割进行预滤波 ,再降低采样率滤波的方法 ,有效地降低了 FM/ FM系统分路滤波器的阶数。     

利用耦合实现相控阵平顶单元方向图的方法

提出了一种用于低轨卫星有限扫描相控阵的新型波导辐射单元。该辐射单元利用缝隙耦合实现平顶的单元方向图以解决有限扫描相控阵由于采用大间距设计时引起的扫描边缘增益跌落和栅瓣的问题。分析采用多模散射矩阵的方法,计算速度快、精度高。数值计算和测试结果表明这种辐射单元可以用于扫描角度宽达 ±25° 的有限扫描相控阵的设计。计算结果同时表明该辐射单元至少可以满足10％的频带宽度。     

纳秒级脉冲照射下隐身飞行器的隐身效果分析

为研究冲激雷达的反隐身潜力,首先用时域有限差分方法(FDTD)和解析方法计算了电磁波垂直入射时涂敷窄带雷达吸波材料(RAM)和涂敷等离子体的金属平板的反射率,然后用FDTD仿真了在纳秒级双高斯脉冲照射下某金属缩比隐身飞机模型,及其涂敷窄带RAM及涂敷等离子体的散射特性。结果表明,在低频区外形隐身的隐身效果不明显,RAM和等离子体的吸波效果较差。这证实了发射纳秒级脉冲信号的冲激雷达具有反隐身的潜力。     

有限变形塑性的指数算法

基于不同的应力测量方式，将得到不同形式的有限变形塑性的算法列式。本文基于柯西应力来表示屈服准则，给出了以柯西应力表示的有限变形塑性的指数算法列式。利用对数应变，在主轴下建立了返回映射算法，得到了一种非对称的算法切线模量。它保持了在无穷小理论中返回映射的算法结构。算例表明，对于静水压力相关的塑性材料，应力测量形式的选择对数值结果有较大的影响，柯西应力列式得到了合理的结果。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。