模拟高度复杂流动的自适应网格算法

提出一种基于对等弧分布律进行扭转修正的自适应网格算法。该方法将三维网格系统分裂为３个单方向自适应网格子系统,在计算中与主控方程的求解格式相耦合,交替地使用。应用所构造的算法结合ＭａｃＣｏｒｍａｃｋ显式格式对喷流／绕流干扰流场进行了雷诺平均Ｎ－Ｓ方程模拟,给出了清晰的剪切层与激波波系结构。计算结果表明,该自适应网格算法可较大程度地提高ＭａｃＣｏｒｍａｃｋ格式的分辨率。     

二元矢量喷管内流特性的数值模拟

应用多重网格ＭａｃＣｏｒｍａｃｋ显示格式对不同矢量角下二元喷管的流动进行了模拟分析。计算时,在细网络上采用了雷诺平均Ｎ－Ｓ方程和Ｋ－ε湍流模型方程;在粗网格上则采用Ｅｕｌｅｒ方程。数值计算给出了较清晰的流谱。壁面压力分布的计算值也与实验结果吻合。通过计算结果分析了喷管几何参数对流场的影响。所得结果显示出多重网格算法高效可靠,可用于开展参数研究     

一种求解S形进气道三维粘性流场的多重网格算法

提出了一种基于多重网格技术,可快速求解Ｓ形进气道粘性内流场的实用工程算法：从Ｂｒａｎｄｔ的多重网格理论出发,采用隐式Ｂｅａｍ－Ｗａｒｍｉｎｇ时间推进格式,结合“锯齿式”迭代技术交替在粗、细两种网格上对雷诺平均三维非定常Ｎ－Ｓ方程进行求解。用此算法数值模拟了某型号Ｓ形进气道的内流场,与试验结果对比显示：计算得到的流场符合一般Ｓ形管道流动特性,所捕获的二次流图谱清晰合理,进气道出口截面处总压恢复系数与试验值较为接近。与普通单层网格计算方法相比,计算结果精度高,收敛速度快。     

可压缩流动模拟的多重网格-扰动域推进方法

针对可压缩流动数值模拟的加速需求，在作者前期提出的扰动域推进计算框架下，融入多重网格加速技术，提出了多重网格-扰动域推进方法。建立了粗、细网格上动态计算域的更新逻辑与高效算法，并针对多重网格方法的特点，对粗网格初始化、增大对流动态域等步骤进行了改进。设计了新型粗网格生成策略，可消除结构网格粗网格生成对单元数目的限制。典型算例验证表明，多重网格-扰动域推进方法可同时降低粗、细网格上的计算量并减少总迭代步数；通过在粗、细网格上同时采用随数值扰动传播而增大、随求解收敛而缩小的动态计算域，多重网格-扰动域推进方法可在传统多重网格方法加速的基础上再将计算效率提升57.9%。     

高阶DG格式多重网格方法构造中的不同隐式策略影响研究

DG方法是一种非常具有潜力的高精度方法,但其在对复杂外形的数值模拟方面仍存在内存需求量大、计算量巨大等不足.为了进一步提高DG方法求解Euler方程的效率,在传统p型多重网格的基础上,结合LU-SGS和GMRES两种隐式迭代方法,研究其整体加速性能.p型多重网格方法通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛的目的.高阶近似(p＞0)使用显式龙格库塔格式,最低阶近似(p=0)使用隐式格式.对NACA0012翼型和ONERA M6机翼跨音速无粘流动进行数值模拟,结果表明:与显式TVD-RKDG时间格式相比,DG(p0)层上采用LU-SGS和GMRES的p型多重网格方法收敛速度均得到明显提高,且GMRES迭代法性能最佳,LU-SGS迭代法次之.     

双时间步长加权ENO-强紧致高分辨率格式及在叶轮机械非定常流动中的应用

在非结构网格下,针对非定常三维N-S方程组发展了一种双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间上将r=3的加权ENO格式与强紧致格式相结合去处理N-S方程中的对流项以及离散方程的右端项,并用四阶精度的紧致格式去计算N-S方程中的粘性项。典型的3个算例从不同侧面对格式进行了考核。计算表明:该算法具有高效率与高分辨率的特征,所得的计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下模拟非定常流动的物理过程。      

p型多重网格间断Galekin有限元方法研究

在二维非结构网格上,使用p型多重网格间断Galerkin方法求解定常可压缩欧拉方程。p型多重网格方法主要特征是通过对不同阶次多项式的近似解进行递归迭代求解。文中高阶近似(p0)上使用显式格式,在最低阶近似(p=0)上选用隐式格式,而非显式格式,从而在保证精度和占用较小内存的情况下加速收敛到定常解。运用该方法对NACA0012跨音速无粘流动进行数值模拟,数值结果表明:p型多重网格方法同单重显式Runge-Kutta方法相比,收敛速度能够提高6倍左右,并且精度保持不变。     

不同矢量偏角下二元喷管内流特性的数值研究

应用多重网格MacCormack显式格式对不同矢量角下二元喷管的内部流动进行了模拟分析。计算时, 在细网格上采用了雷诺平均N-S方程与k-ε湍流模型方程;在粗网格上则采用Euler方程。数值计算给出了较清晰的流谱, 壁面压力分布的计算也与实验结果相吻合。通过计算结果分析了矢量喷管几何参数对流场的影响, 并对矢量喷管的推力特性进行了数值研究。      

局部网格加密方法计算不可压N—S方程

本文发展了二维非定常不可压N-S方程的局部网格加密方法(LMR):需加密区域预先指定;仅对空间方向加密,分粗细二层网格,细网格覆盖在粗网格上;粗和细网格上分别用显式和隐式差分格式。结合压力修正法类的余量型差分格式,恰当地处理了粗细网格之间的信息传递,使得粗细网格交接面上质量守恒,即满足连续性条件。粗网格通过插值给细网格边界值影响细网格;细网格通过粗细网格压力Poisson方程耦合求解影响粗网格,并且由于压力Poisson方程从动量方程形成,在细网格覆盖下的粗网格上,压力Poisson方程是细网格方程的组合,从而耦合求解时粗细网格压力Poisson方程不需迭代。 本文中计算了二维方腔内的自然对流模型问题,Rayleigh数是10~6。粗网格用显式ULWC格式,细网格用隐式余量型近似因式分解格式。Poisson方程用快速直接算法PO1STG(在FISHPACK中)。粗网格距△=1/16,细网格距△=1/64,加密1/4区域,计算到定常局部网格加密方法所需CPU时间比同等均匀网格(△=1/64)节约一倍以上。     

动能BGK格式模拟高雷诺数定常流的收敛性研究

为克服动能BGK格式计算耗时多、收敛慢的缺点,在格式中引入当地时间步长、隐式LU-SGS方法和多重网格技术,基于RAE2822翼型粘性跨声速绕流的数值模拟,对其收敛性的改善进行了研究。结果表明:上述3种加速手段对求解定常问题时的收敛性有明显的改善,达到定常解所需的迭代步数随着CFL数的提高而减少;因而求解某类问题时,隐式格式比显式格式更具优势;多重网格既适用于显式格式,也适用于隐式格式,加速性更加突出。     

有翼高超声速再入飞行器气动设计难点问题

有翼高超声速再入飞行器是近年来的研究热点,气动设计是飞行器设计的关键。为了更清楚地认识有翼高超声速再入飞行器气动设计的难点问题,对有翼高超声速再入飞行器的发展、优势及总体任务剖面进行了介绍,从5个方面详细介绍了该类飞行器气动设计的难点问题,包括多约束复杂面对称气动布局设计、高温真实气体效应对气动特性影响、天地差异与天地换算方法、反作用控制系统(RCS)喷流干扰对气动特性的影响以及气动数据不确定度等,简要阐明了这些难点问题对总体设计的重要性以及初步的解决思路,为有翼高超声速再入飞行器气动设计提供了一些参考。     

Simulation and Analysis of Crashworthiness of Fuel Tank for Helicopters

Crashworthiness requirement of fuel tanks is one of the important requirements in helicopter designs. The relations among the protection frame, textile layer and rubber layer of the fuel tank are introduced. Two appropriate FE models are established, one is for an uncovered helicopter fuel tank without protection frame, and the other is for fuel tank with protection frame. The dynamic responses of the two types of fuel tanks impinging on the ground with velocities of 17.3 m/s are numerically simulated for the purpose of analyzing energy-absorbing capabilities of the textile layer and protection frame. The feasibility of the current crashworthiness design of the fuel tank is examined though comparing the dynamic response behaviors of the two fuel tanks.     

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

突破技术趋同，波音再现竞争优势

 波音公司面临着来自空客公司的巨大挑战,企业战略性创新才是公司成功的关键。为此波音公司的全部战略性研究集中在扩大产品的差异性上,体现在3个方面：电子化（e-Enabled）运营环境、整体复合材料机身部件的制造技术和支持波音787客机的全球协同环境（GCE）。     

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

基于弯曲激波压缩系统的高超声速进气道反设计研究进展

总结了近十年来弯曲激波压缩研究的主要成果。提出了弯曲激波压缩系统的新概念,即利用特殊设计的楔形弯曲压缩面或空间弯曲压缩面,产生一系列与前缘弱激波相互交汇或叠加的压缩波系,从而使前缘激波弯曲,形成特殊的弯曲激波,它与波后的等熵压缩波来共同完成对气流的压缩。在此基础上,实现了由给定出口气动参数的超声速内流道反设计,实现了由给定压缩面压力分布和给定压缩面马赫数分布要求的型面反设计,实现了由给定激波波面的压缩型面反设计。研究证明,弯曲压缩面-弯曲激波压缩系统具有良好的综合气动性能,为高性能高超声速进气系统的气动设计提供了一种全新的设计方法。     

临近空间飞行器测控与信息传输系统频段选择

 临近空间飞行器是高性能信息化武器平台,测控（TT&;C）与信息传输系统是其信息保障的核心,而选择合理、可行的频段是展开系统设计的前提和基础。频段选择影响到整个技术方案的制定,是一个需综合考虑、影响深远并具有战略意义的关键问题,从国际电联（ITU）国际标准、高速数传、接收信噪比（SNR）、“三抗”、超视距中继、黑障、雨衰以及设备研制成熟度8个方面全面、细致论证了近空间平台测控系统的频段选择问题,最终得出在视距链路中以Ka频段为宜,在超视距链路中以Ku/Ka双频段为宜的结论。     

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

Experimental Study on Effects of Fuel Injection on Scramjet Combustor Performance

In order to investigate the effects of fuel injection distribution on the scramjet combustor performance, there are conducted three sets of test on a hydrocarbon fueled direct-connect scramjet test facility. The results of Test A, whose fuel injection is carried out with injectors located on the top-wall and the bottom-wall, show that the fuel injection with an appropriate close-front and centralized distribution would be of much help to optimize combustor performances. The results of Test B, whose fuel injection is performed at the optimal injection locations found in Test A, with a given equivalence ratio and different injection proportions for each injector, show that this injection mode is of little benefit to improve combustor performances. The results of Test C with a circumferential fuel injection distribution displaies the possibility of ameliorating combustor performance. By analyzing the effects of injection location parameters on combustor performances on the base of the data of Test C, it is clear that the injector location has strong coupled influences on combus- tor performances. In addition, an inner-force synthesis specific impulse is used to reduce the errors caused by the disturbance of fuel supply and working state of air heater while assessing combustor performances.     

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。