加权梯度限制器在节点DG格式中的应用

在使用节点间断Galerkin (DG)方法求解存在激波的流动问题时,为抑制Gibbs效应导致的非物理振荡,对一种加权梯度限制器进行研究,在二维非结构网格上对定常欧拉方程求解.限制器以原始变量作为限制对象,对其梯度进行限制,以加权形式进行重构.对翼型的跨音速流场、超燃冲压发动机内流场和前台阶模型流场进行模拟,结果表明限制器能很好抑制间断解附近的非物理振荡,有效避免计算中断,体现了DG格式精度高、对间断捕捉敏感的特点.     

大型航天器再入解体气动力热特性模拟的直接模拟蒙特卡洛方法研究

为模拟大型航天器离轨再入近连续过渡流区高超声速气动力/热绕流特征,构建了基于直接模拟蒙特卡洛法碰撞限制器技术的混合方法,发展了基于密度梯度的动态自适应混合网格处理技术与变时间步长计算方案。利用当地流动梯度的克努森数作为判断连续流失效的参数,将流场划分为不同区域,在连续流区采用碰撞限制器以及大网格尺度和大时间步长,在流场的大梯度区域——包括激波和壁面边界层区域——采用基于当地密度梯度的动态自适应碰撞网格和取样网格处理技术。为保证整个流场范围每个碰撞网格内的模拟粒子数分布更加均匀,采用变时间步长计算方案,并固定当地时间步长与粒子权重的比值,避免了因分子穿越网格界面产生的复制或消失。通过计算类天宫飞行器低密度风洞试验状态的气动力系数,并与试验数据对比,验证了上述算法的高精度模拟能力与可靠性。同时模拟分析了带太阳电池帆板的类天宫飞行器再入85 km高超声速复杂气动力热,及头部对接台与板舱非规则物形绕流所致激波/边界层干扰、流动分离与强气动力热致太阳电池帆板毁坏发生首次解体机制。     

基于二维结构化网格的间断Galerkin方法研究

在二维结构化网格上,对求解Euler方程的间断Galerkin方法进行了研究。应用二阶、三阶精度的问断Galerkin方法对绕NACA0012翼型的跨音速无粘流动进行了数值模拟。在整体逼近精度基本保持不变的前提下,采用moment限制器以消除数值解在激波附近的伪振荡。数值结果表明,间断Galerkin方法能够很好地模拟流场,准确捕捉激波;moment限制器成功抑制了激波附近的伪振荡;高阶格式具有比低阶格式更强的激波捕捉能力。     

MUSCL格式在跨声速转子流场计算中应用

以工程应用为背景,采用Spalart-Allmaras湍流模型和LU-SGS隐式解法,开发了基于MUSCL格式的CFD程序,并将程序应用到压气机跨声速转子内部流场的数值模拟中,研究了5种差分格式限制器的计算效果.与Rotor 67和Rotor 37转子试验数据的对比表明:当程序采用VanLeer限制器,且入口湍流粘性取值约为3.6倍的分子粘性时,收敛效果最好,能较准确地模拟出流场的总性能参数和激波结构.     

非结构网格、混合网格下计算方法研究

为了提高计算精度，本文研究了计算格式和网格技术。通过比较非结构网格下的NND格式、MUSCL格式、ENO格式及最小二乘格式的计算精度和波动，选择利用NND方法中的守恒变量法来建立非结构网格下的计算格式，通过改进NND方法中的限制器来提高计算精度和分辨率，计算结果比较显示出该格式具有很高的分辨率和计算精度。本文也研究了非定常网格自适应技术，计算给出了双马赫反射问题的新旧限制器及网格自适应计算结果比较，新限制器网格自适应后计算得到的滑移线出现了剪切涡现象，这表明，数值计算中取足够密的网格和构造高精度格式是非常必要的。     

基于非结构动网格的非定常激波装配法

从ALE方程出发,通过修正体积的新思路构建离散几何守恒律,基于非结构动网格的信息传递方法得到网格重构以后新网格的流场参数,采用一阶格式捕捉法得到的流场确定装配法的初始激波位置,成功地把激波装配法应用到非定常流动的模拟。模拟超声速飞行环境下整流罩分离的非定常流固耦合过程,在马赫数6条件下装配法和捕捉法的计算结果完全一致,随着马赫数增加,二阶精度格式的捕捉法需要进行调整限制器参数、时间步长等人工干预,但是装配法直到马赫数20计算过程没有出现异常。该基于非结构动网格技术的激波装配方法,与传统的激波装配方法相比,能够用于复杂外形产生的非规则形状激波和非定常流动产生的运动激波,解决二阶精度有限体积法用于高马赫数飞行条件下多体分离过程的流动模拟时遇到的稳定性问题。     

涡扇发动机加力燃烧室扩压器流场的数值计算

 本文用数值计算的方法得到了某型涡扇发动机加力燃烧室扩压器内的流场,与台架状态测得的实验数据基本吻合。研究了不同进口气流参数分布对扩压器流场的影响,计算了台架不同转速状态,不同飞行工况下扩压器内的流场。     

结冰区流场特征一致的混合缩比翼型设计

提出了一种基于翼型前缘范围流场相似的多控制点混合缩比翼型优化设计方法:以前缘局部范围内的空间流场特征相同为目标,采用多目标遗传优化算法结合翼型参数化、多控制点、网格自适应技术及CFD,进行了不同工况下的混合缩比翼型设计,分别通过低湍流度风洞及结冰风洞对翼型的气动性能及积冰过程进行了实验。设计及实验结果表明:该混合缩比翼型与对应的全尺寸翼型在前缘15%弦长范围内的压力系数偏差均在5%以内,翼型绕流相似;在相同时间历程内两者的冰型增长及气动特性均具有相同的变化规律。该设计方法可以应用于不同实验状态下的混合缩比翼型的设计,提高多状态下结冰实验效率。     

基于深度神经网络的超声速隔离段湍流涡黏性系数辨识

雷诺平均Navier-Stokes(RANS)方程由于计算成本较低，当前仍然在工程设计领域广泛应用。为了进一步提升计算精度和减少时间，应用深度神经网络（deep neural networks, DNN）方法自适应辨识稳态湍流涡黏性系数。以隔离段的激波串前缘位置检测流场生成为例，使用Wilcox-2006 k-ω湍流模型进行模拟。在不同的背压情况下，产生稳定状态的湍流涡黏性流场作为模型学习的训练数据集。最后在不同背压条件下开展了测试。结果表明：提出的DNN方法能够快速预测湍流涡黏性系数的值，预测值与计算流体力学数值模拟计算的参数值相比，方均根误差较小，可决系数大于99%，预测的流场结果与真实流场基本一致，进一步验证了深度学习技术在湍流模型参数辨识的可行性。     

非结构网格通量重构算法下三种紧致WENO限制器对比研究

在通量重构算法框架下对比研究了3种最新提出的用于间断伽辽金算法的紧致型WENO(weighted essentially non-oscillatory)限制器，即简单WENO限制器、p阶加权WENO限制器和多精度WENO限制器。这3种WENO限制器能够较高精度地模拟流场并捕捉激波，其紧致性在于问题网格单元中的重构只涉及问题网格单元及其相邻网格单元。同时，使用保精度的保正限制器用以避免可能出现的非物理的密度和压力负值。最后，采用非结构四边形网格在双马赫反射、激波与涡相互作用、激波反射、黏性弓形激波和激波与层流边界层相互作用等多个二维算例中，对这3种WENO限制器进行分析比较。结果显示：多精度WENO限制器与p阶加权WENO限制器能够高精度模拟流场并捕捉激波，同时一定程度抑制通量重构方法本身在激波附近的数值伪振荡；p阶加权WENO限制器与多精度WENO限制器相比，其稳态收敛性相对更好；简单WENO限制器则性能较差。根据以上研究，提出亟须发展能够使高阶WENO限制器在稳态问题中收敛的间断探测器。     

迎角过载边界限制器的设计方法

介绍了迎角过载边界限制器的意义、设计要求和一些有代表性的国内外飞机的边界限制器设计方法。边界限制经历了硬限制、软限制、反馈限制和比较逻辑反馈限制几个阶段。分析了这些设计方法的优缺点,并对比较逻辑反馈限制器的设计方法进了理论推导,提出了这类边界限制器重要参数的确定方法。     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

Analysis and Design Considerations of Hard Limiters for LF and VLF Navaid Receivers

Hard limiters, followed by a D-type flip-flop as a digital-signal-polarity detector, are very effective receiver/phase detectors for low frequency (LF) and very low frequency (VLF) navigation receivers. However the performance not only depends on the signal quality, but also on the specifications of the hard limiter and the flip-flop. Analysis of the tracking accuracy is given as a function of the dc offsets of the limiter and the flip-flop, the linear gain of the limiter, the signal-to-noise ratios of one or more input signals, and the power consumption of the limiter. The results are presented in formulas and graphs for application by circuit designers. A design example of a low-power, high-gain limiter is given.     

Carrier and intermodulation performance of limiters excited bymulticarriers

A general formula, using a Fourier-series representation, is derived for the input-output characteristic of ideal hard limiters, biased hard limiters and biased soft limiters. Using this formula one can obtain an expression, involving Bessel functions for calculating the amplitudes of the output products of a limiter excited by a multicarrier signal plus Gaussian noise. In contrast with previously published results, the implementation of the present expression requires the computation of functions that are available in most mathematical-library subroutines and can even be implemented by using programmable pocket calculators. Results obtained by using the present expression are in excellent agreement with the previously published results. Additionally, the expression is used to generate a soft- and a hard-limiter signal-suppression database for the case in which a single large carrier and a number of equal-power small carriers exist in a noisy environment     

带通限幅器对锁相环噪声性能的影响

锁相环路中,在鉴相器前加上AGC（自动增益控制）,可以使环路性能在输入信噪比发生变化时保持稳定,带通限幅器也可以起到类似作用。尽管有一些文章对限幅器的性能作了研究,但是在PLL（锁相环）闭环条件下,鲜有文章分析其对环路噪声性能的影响。因此,文章以基本PLL环路的数学模型为基础,推导出PLL中分别加入相干AGC和带通限幅器后新的环路数学模型,以一阶PLL为例,计算和比较不同控制方式下环路相位误差均方值的变化,最后说明应当根据飞行任务的不同选择对环路更有利的控制方式。     

任意单元间断Galerkin有限元计算方法研究

基于龙格库塔间断Galerkin(RKDG)有限元法的构造思想,通过局部坐标变换,发展了非正交单元DG有限元计算方法;借鉴非结构网格有限体积隐式计算方法,发展了适应于DG有限元方法的隐式计算方法;借鉴一维和二维(三角形单元)DG有限元限制器构造方法,提出了非正交三棱柱单元限制器方法;利用上述方法数值模拟了球头和双椭球的高超声速粘性绕流,得到了清晰的流场结构,并得到了较好的压力和热流分布,表明该方法在复杂高超声速流动的数值模拟方面具有广阔的应用前景。     

具有迎角／过载限制器的俯仰增稳系统控制律的研究

本文以JJ-6飞机为例,研究了采用具有迎角/过载限制器的俯仰增稳系统对该机在大迎角下稳定性和操纵性的影响.并根据该机在大迎角下的偏离/失速运动机理,探讨系统中限制器边界值确定的原则。研究结果表明,该系统具有抗偏离/失速的能力,进而可以防止尾旋发生。但在大迎角下过大的横侧操纵,飞机仍能进入尾旋运动。     

基于自适应网格的钝头体激波诱导燃烧现象数值模拟研究

钝头体激波诱导燃烧是爆震研究的一个基本问题。针对化学恰当量比的H2/Air预混气在Ma=4.79和Ma=6.46时的激波诱导燃烧现象开展数值模拟研究，采用基于有限体积法的块结构自适应网格加密程序AMROC对带化学反应源项的轴对称Euler方程解耦求解，考察了数值模拟中不同形式的MUSCL重构格式、限制器类型以及化学反应机理等重要因素对模拟结果的影响。结果表明，程序能够根据设定的加密判据较好地实现网格自适应加密，减小总网格量，实现高效数值模拟。通过与实验数据的对比，表明非定常激波诱导燃烧算例的准确程度不仅取决于化学反应机理，也取决于限制器类型，而采用两种不同形式的MUSCL重构格式获得的振荡频率则几乎一致，与试验结果的误差分别为1.17%和0.97%。模拟对比经典的Jachimowski机理和近年来新发展的几种包含压力相关反应步的氢/氧反应机理，模拟结果表明：对于Ma=4.79时的非定常激波诱导燃烧模拟，经典的Jachimowski机理仍然是能够给出与实验结果最接近的反应机理；而对于Ma=6.46时的定常激波诱导燃烧模拟，几种反应机理均能给出与实验吻合较好的结果。     

星敏感器低频误差与陀螺漂移离线校正方法

获取高精度事后姿态数据是提高遥感平台成像质量的必要条件之一,离线处理可有效降低敏感器测量误差,从而获得更高的姿态确定精度。基于滤波的校正方法中,星敏感器低频误差(LFE)与陀螺漂移将产生耦合影响导致校正精度低,本文针对该问题推导了耦合误差的数学模型,并设计了一种两步双向平滑事后处理算法,将陀螺漂移与低频误差分两步校正,通过反复滤波剥离陀螺漂移与低频误差。同时,针对低频误差参数收敛速度慢、噪声参数调节困难的问题,利用一种基于极大似然估计(MLE)的固定窗口自适应双向滤波算法进行处理以获得更好的噪声估计,提高了收敛速度和收敛精度。文中仿真工况下,离线姿态确定精度可达到0.8″(3σ),低频误差参数完全收敛时间不超过4个轨道周期。     

不完全数据回归分析方法

基于不完全数据顺序统计量,提出一种不完全数据回归分析方法。建立不完全数据回归方程,给出回归系数的最佳无偏整体估计及其协方差矩阵。文中详细讨论了工程中常见的极值分布、Weibull分布和正态分布的等尺度和非等尺度(或异方差)线性回归分析。该方法可以将各个不同状态下的不完全数据作为一个整体进行统计分析,具有信息量大、精度高的特点。