Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

一种基于约束框架的棱柱网格生成方法

提出了一种基于约束框架的棱柱网格生成方法,即利用物面网格节点与背景结构框架网格之间的单一映射关系线性插值生成棱柱网格。背景结构框架网格基于TFI方法生成。物面网格节点与背景框架网格之间的映射关系通过一种线性投影方法建立。然后通过V型槽、圆球、M6机翼和F6翼身组合体四种不同外形对该棱柱网格生成方法进行了测试,测试结果表明该方法能够避免棱柱网格生成中的网格交叉现象,保证棱柱网格的成功生成。该方法在棱柱网格质量控制、棱柱网格局部修改等方面具有较好的效率和优势,但背景框架网格的生成依然是该方法的一个瓶颈,需要进一步的研究。     

基于超网格的重叠网格守恒插值方法

保证重叠网格边界数据插值的守恒性是计算流体力学面临的一大难题。基于超网格方法与格心格式有限体积法,发展了一种新型混合重叠网格守恒插值方法。详细研究了用网格切割技术在重叠网格边界构造局部超网格,用网格求交算法合理地扩大贡献单元模板,建立了一种适用于任意网格类型的隐式并行重叠网格守恒插值方法,以超网格为媒介可实现重叠网格二阶精度的守恒插值。数值结果表明,本文方法对二阶分布的流场变量具有严格的守恒性,相比三线性插值方法和逆向距离权插值方法,本文方法减小了数值误差,提高了重叠网格边界的插值精度,加快了计算收敛速度,改善了重叠区域网格尺度相差较大时流场的光滑性和连续性。     

基于误差估计的伴随网格自适应方法

网格自适应是一种有效减少解与网格相关性,改善计算结果的重要途径。给出了一种基于误差估计的伴随网格自适应计算技术,定义了基于误差估计的自适应准则,给出了非结构四面体网格单元细化剖分规则,以及一般的网格自适应流程,介绍了一种基于Rechardson外推的网格收敛技术,运用本文提供的方法对NLR7301翼型进行了数值计算研究,数值试验表明给出的网格自适应具有较高的可有效改善计算结果。     

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

一种用于二维燃烧室流场计算的非结构网格的生成方法

本文依照Ｄｅｌａｕｎａｙ 准则构造二维非结构网格,并提出了一种与之相应的基于阵面推进法的自动生成新的内部节点的方法。文中给出了利用这一方法对边界形状复杂的燃烧室流场进行网格划分时,所得的非结构网格,及网格相应的质量统计特性。从这些算例可以看出,这一方法能够较好地适应复杂的几何外形,得到分布合理,并且大多数单元形状接近正三角形的非结构网格     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

非结构聚合多重网格法流场数值模拟研究

 在非结构网格上应用聚合多重网格技术来加速 Euler方程的收敛过程。设计了一种高效率的网格聚合方法。采用四重网格,在每一层网格上采用有限体积法进行计算。通过对多段翼型绕流的数值求解验证了多重网格加速收敛的高效性。     

面向映射单调性的TC4初生α相晶粒尺寸超声评价方法

随着被检对象微观组织结构复杂信息的增加,以误差最小为目标来建立有效的评价曲线越来越难,针对这一问题提出了一种以单调性为目标的TC4钛合金初生α相晶粒尺寸超声评价方法。依据相关性度量从多个超声参数中选取有效参数,经映射函数降成单维参数并归一化处理后与晶粒尺寸进行一次项拟合,构建以单调性（特征参数样本点依次序列差值同时为正或为负）最大为目标的优化问题并结合自适应差分进化（SADE）算法进行求解,寻找最佳的映射函数和拟合函数系数,最终建立映射单调性的多参数超声评价模型。实验结果表明,因考虑了单调性的重要性且与误差为目标的评价模型相比,所建模型的评价效应更显优势。与单一参数的声速模型、衰减系数模型和非线性系数模型对比分析,所建模型误差小、单调性好、性能稳定且有着评价结果精度高等特点。     

非结构网格限制器研究

研究了6种限制器形式及其参数计算方法,其中新的自适应VenkatakrishnanM限制器能够根据当地流场特征自动调节参数从而实现不同区域不同特性.数值试验表明:非结构网格限制器对邻居集合选择方法较敏感,其中共面法比共点法精度高约15%;高马赫数时压缩型限制器比耗散型限制器计算精度高约10%~40%;参变量算法中,方向梯度法得到的流场比增量法平滑;自适应VenkatakrishnanM限制器能够模拟复杂超声速流动.      

高超声速压缩拐角湍流计算的收敛性分析

应用基于shear stress transport(SST)湍流模型的隐式紧耦合方法对高超声速压缩拐角流动进行了数值计算.考察了差分格式、限制器以及湍流时间步长等7个因素对湍流计算过程稳定性和收敛性的影响.结果表明:差分格式和限制器对收敛性有很大影响,总体而言耗散性越大的差分格式和限制器收敛性越好;增加湍流时间步长对收敛性影响很小;适当降低湍流变量最低限制值有利于湍流方程的收敛;增大Courant-Friedrichs-Lewy(CFL)数和隐式内迭代次数可显著加速收敛,但均存在最优值,超过最优值则加速效果不明显.      

An Analysis of the Performance of the Oscillating Limiter Driven by FM Signals Corrupted by Noise

The behavior of the oscillating limiter (OL) driven by FM signals is surveyed, and its performance with signal corrupted by noise is investigated. For high values of the carrier-to-noise ratio (CNR), if the frequency deviation of the signal is small in comparison with the locking range of the OL, it is calculated, and experimentally verified, that a system OL discriminator is equivalent to a system bandpass limiter discriminator followed by a linear network whose frequency response has been specified. When the frequency deviation is not so small, the baseband noise power increases with it; a formula is given that allows the calculation of this power when the signal is such that the circuit operates in quasistationary fashion. For low values of the CNR, a mathematical analysis presents unsurmountable difficulties. However, heuristic argumentation leads to an interpretation of the operation of the OL in the threshold region, which is substantiated by an experimental investigation. The results of this paper enable a comparative evaluation of a system OL discriminator and a system bandpass limiter discriminator, to which the former reduces when the feedback path in the OL is open.     

Analysis of a Class of Nonlinear Sampled-Data Systems with a Gaussian Input Signal

The problem treated here is the analysis of a class of nonlinear sampled-data systems with a Gaussian input signal. The nonlinear element is a symmetrical limiter (memoryless device) that feeds the linear element. Preceding the limiter is a sampler and zero-order hold unit for data reconstruction. The analysis is predicated on finding an equivalent gain for the nonlinear element such that the mean-square error is a minimum. The nonlinearity is thus replaced by a linear component and the system is then analyzed by conventional techniques. It is also shown that for appropriate sampling rates the sampled-data system is equivalent to a continuous system due to the action of the sampler-hold combination. All calculations pertaining to the limiter input signal require the use of a transfer function that is interior to the feedback loop. However, in developing this transfer function, the overall system is considered so that the system output may be evaluated after the value of equivalent gain has been found.     

Limiting the Impulsive Noise in a PLL

This paper analyzes the operation of a phase-locked loop (PLL) preceded by a bandpass limiter (BPL) in the presence of impulsive noise. It is shown that the effect of the limiter consists essentially in a change of the statistics of the pulse strengths of the noise, so that the behavior of a BPL + PLL can be deduced from that of a simple PLL by suitable adjustment of the noise model. It is found that the limiter greatly enhances the PLL performance by reducing both the phase-error variance in the loop and the probability of cycle slippage. Finally, the design of the filter of the BPL is discussed, resulting in the conclusion that the best results are obtained by using a singletuned RLC circuit.     

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

基于密切原理的Bump进气道外压缩鼓包逆向设计

基于密切乘波理论提出一种Bump进气道外压缩鼓包的设计方法,可根据制定的激波形状及其曲率中心分布来逆向求解外压缩鼓包型面。通过引入曲率中心分布这一变量,可以控制横截面激波形状并调节外压缩鼓包的三维外型及其表面横向压力分布,进而提高外压缩鼓包的附面层排移能力。同时,发展了一种Bump进气道的流量系数快速估算法,能够在设计初期以不超过2%的误差快速给出进气道的流量系数。结果表明:基于密切原理的外压缩鼓包设计有利于改进Bump进气道的流量捕获和附面层排移能力。算例中,较锥导鼓包模型,密切鼓包方法设计的Bump进气道流量系数提升4.03%,附面层排移能力提升2.12%。      

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

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

大迎角限制器设计与飞机敏捷性

在介绍高机动性飞机大迎角限制器和飞机敏捷性要求的基础上，以典型战半机为例，计算了飞机的敏捷性尺度，根据计算结果进一步分析大迎角限制器对飞机敏捷性的影响，并以空战周期时间为例，说明迎角限制器参数的选择应兼顾飞机的敏捷性与飞行安全两方面。