高超声速再入体表面热流计算

运用非结构网格计算外部无粘流场,并结合边界层内粘性主导区域的工程算法,计算高超声速飞行器的气动加热。通过求解三维Euler方程确定复杂外形飞行器的边界层外缘参数,在理论与经验公式的基础上,利用局部相似性解的方法计算了钝锥和钝双锥外形有攻角再入的表面热流,并与国内外文献的NS方程数值计算结果和风洞试验结果进行了对比,三者的结果吻合较好。     

高超声速飞行器关键部位气动热计算

运用快速算法对高超声速飞行器外表面的一些关键部位经受的气动热环境进行计算分析。在理论和经验公式的基础上,利用轴对称比拟法考虑攻角影响,采用局部相似性解及参考焓等方法确定飞行器有攻角再入的表面气动加热,发展了一套高超声速飞行器关键部位气动热的计算方法。以钝锥为算例对计算方法进行了验证,结果表明,本文所述方法具有较高的效率和精度。     

高速双锥绕流中热化学与输运模型影响研究

激波与边界层之间相互作用是高超声速飞行中的常见现象,对飞行器气动性能与飞行安全至关重要.对于高焓来流,流场中通常存在复杂的物理化学现象,此时准确模拟流场中激波边界层相互作用的难度大,相关物理化学建模仍有待进一步考察和研究.本文针对最近文献中纯净空气高超声速双锥绕流实验开展数值研究,分别研究了不同热化学模型与输运模型对壁面压力与热流的影响.热力学模型包括完全气体、热力学平衡和非平衡模型,化学模型包括冻结和非平衡化学模型,输运模型包括经典的Wilke/Blottner/Eucken模型与更加复杂的Gupta/SCEBD模型,以及考虑壁面催化/非催化影响的模型.计算了6个不同算例,涵盖了低焓至高焓来流等不同工况.壁面压力与热流的数值计算结果与实验结果符合较好;对于低焓来流,计算结果主要受到分子内能分布的影响,输运模型对计算结果的影响不大;对于高焓来流,一方面计算结果受到化学反应与壁面催化的影响较大,另一方面不同输运模型对计算结果的影响也更加明显.     

高速双锥绕流中热化学与输运模型影响研究

激波与边界层之间相互作用是高超声速飞行中的常见现象,对飞行器气动性能与飞行安全至关重要.对于高焓来流,流场中通常存在复杂的物理化学现象,此时准确模拟流场中激波边界层相互作用的难度大,相关物理化学建模仍有待进一步考察和研究.本文针对最近文献中纯净空气高超声速双锥绕流实验开展数值研究,分别研究了不同热化学模型与输运模型对壁面压力与热流的影响.热力学模型包括完全气体、热力学平衡和非平衡模型,化学模型包括冻结和非平衡化学模型,输运模型包括经典的Wilke/Blottner/Eucken模型与更加复杂的Gupta/SCEBD模型,以及考虑壁面催化/非催化影响的模型.计算了6个不同算例,涵盖了低焓至高焓来流等不同工况.壁面压力与热流的数值计算结果与实验结果符合较好;对于低焓来流,计算结果主要受到分子内能分布的影响,输运模型对计算结果的影响不大;对于高焓来流,一方面计算结果受到化学反应与壁面催化的影响较大,另一方面不同输运模型对计算结果的影响也更加明显.      

局部凸起在V形钝前缘模型中的降热特性研究

三维内转式进气道的唇口结构通常存在复杂的激波干扰及严酷的气动热载荷,严重威胁高超声速飞行器的性能与安全.在6.0马赫的高超声速流动中,以V形钝前缘模型为研究对象,设计了局部凸起的被动流动控制降热方案.采用数值模拟手段,首先研究了局部凸起方案的降热能力以及降热原理,然后初步优化了局部凸起的位置、高度以及宽度等关键设计参数,最后分析了优化后的局部凸起方案的攻角、侧滑角及马赫数的适用性.研究结果表明:上游凸起边缘形成的斜激波与主马赫反射结构形成的透射激波发生干扰,能够减弱其冲击壁面的强度,实现降热的目的;驻点凸起通过改变超声速射流的对撞角度,能够降低其对撞的强度,实现降热的目的.原始方案的降热能力约为37.75%,在对局部凸起的关键设计参数进行初步优化后,优化方案的降热能力将提升至44.60%.设计工况下的优化方案具有良好的攻角适用性,而高度可变的优化方案可以较好地适用于有侧滑角及高马赫数的流动.在研究范围内,高度可变的优化局部凸起方案的降热能力均高于20%.     

高超声速飞行器气动防热新概念研究

传统乘波构型的高超声速飞行器尖锐的前缘存在严重的气动加热问题,而简单的前缘钝化气动防热方法由于造成很大的升阻比损失,难以发挥实质性作用. 引入``人工钝前缘(ABLE)'概念,拟以一种新的思路解决这一矛盾. 通过定义ABLE构型的外形参数,并采用CFD数值计算方法研究了各参数对气动力和气动热特性的影响规律,在流场分析的基础上进行了外形优化,最终得到令人满意的新型高超声速飞行器头部外形,总结了运用ABLE概念进行气动防热的相关设计原则和规律.      

多维迎风方法在非结构/混合网格热流计算中的应用研究

非结构/混合网格具有极强的几何灵活性,在复杂外形飞行器的气动力特性数值模拟中已得到广泛应用,但目前还难以准确地预测气动热环境。本文从非结构/混合网格热流计算的三个需求出发,选取了多维迎风方法,并与其他方法进行了对比研究。以二维圆柱高超声速绕流这一Benchmark典型问题为例,对比研究了多维迎风方法和几种广泛使用的无粘通量格式(Roe格式、Van Leer格式和AUSMDV格式)对混合网格热流计算精度的影响。结果表明,多维迎风方法在热流计算精度、鲁棒性以及收敛性方面表现良好。最后,将多维迎风方法应用于常规混合网格上的圆柱和钝双锥绕流问题,均得到了较好的热流计算结果,为非结构/混合网格热流计算在复杂高超飞行器中的应用奠定了基础。     

乘波体前体/进气道一体化设计研究

为了研究乘波体几何外形参数和飞行参数对前体/进气道一体化设计的影响,采用理论分析和数值模拟相结合的方法,以马赫数Ma=6和攻角α=0为设计状态、进气道总压恢复系数和前体阻力系数为目标函数,对乘波体前体/进气道进行了优化设计,并在此基础上研究了攻角、马赫数、前缘半径、前体宽度对气动参数的影响。结果表明:该乘波体前体/进气道构型具有良好的攻角特性,总压恢复系数比基准构型提高17.79%,阻力系数比基准构型降低78.5%,符合高超声速飞行器高升力、低阻力的要求,且非常适合小攻角高超声速巡航飞行;为了得到较高升阻比的前体,在前缘半径R≤2mm的范围内进行流场反设计时,可以将设计马赫数的取值比预期低一些。     

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

高超声速湍流边界层研究现状

一、高超声速湍流边界层研究的重要性随着再入导弹武器从惯性弹道导弹发展到可作机动飞行的多弹头分导弹道导弹,以及航天飞机的出现,高超声速再入飞行器的气动外形变得更复杂了。由于出现了多个激波的相互干扰,激波与边界层的相互影响,以及边界层的分离(入射激波和后台阶产生      

Experimental demonstration of a new concept of drag reduction and thermal protection for hypersonic vehicles

A new idea of drag reduction and thermal protection for hypersonic vehicles is proposed based on the combination of a physical spike and lateral jets for shockreconstruction. The spike recasts the bow shock in front of a blunt body into a conical shock, and the lateral jets work to protect the spike tip from overheating and to push the conical shock away from the blunt body when a pitching angle exists during flight. Experiments are conducted in a hypersonic wind tunnel at a nominal Mach number of 6. It is demonstrated that the shock/shock interaction on the blunt body is avoided due to injection and the peak pressure at the reattachment point is reduced by 70% under a 4° attack angle.     

Shock tunnel study of spiked aerodynamic bodies flying at hypersonic Mach numbers

A three-component accelerometer balance system is used to study the drag reduction effect of an aerodisc on large angle blunt cones flying at hypersonic Mach numbers. Measurements in a hypersonic shock tunnel at a freestream Mach number of 5.75 indicate more than 50% reduction in the drag coefficient for a 120° apex angle blunt cone with a forward facing aerospike having a flat faced aerodisc at moderate angles of attack. Enhancement of drag has been observed for higher angles of attack due to the impingement of the flow separation shock on the windward side of the cone. The flowfields around the large angle blunt cone with aerospike assembly flying at hypersonic Mach numbers are also simulated numerically using a commercial CFD code. The pressure and density levels on the model surface, which is under the aerodynamic shadow of the flat disc tipped spike, are found very low and a drag reduction of 64.34% has been deduced numerically.     

A new method for characterizing patterns of neural spike trains and its application

A method for characterizing and identifying firing patterns of neural spike trains is presented. Based on the time evolution of a neural spike train, the counting process is constructed as a time-dependent stair-like function. Three characteristic variables defined at sequential moments, including two formal derivatives and the integration of the counting process, are introduced to reflect the temporal patterns of a spike train. The reconstruction of a spike train with these variables verify the validity of this method. And a model of cold receptor is used as an example to investigate the temporal patterns under different temperature conditions. The most important contribution of our method is that it not only can reflect the features of spike train patterns clearly by using their geometrical properties, but also it can reflect the trait of time, especially the change of bursting of spike train. So it is a useful complementarity to conventional method of averaging.     

智能物面对非定常分离流的最优自适应控制

以低雷诺数二维大攻角翼型绕流为研究对象, 将非定常动边界计算流体力学方法与 最优控制方法有机结合, 研究二维不可压非定常流智能物面最优自适应流 动控制的理论与算法, 并将其用于固定攻角和俯仰振荡翼型绕流. 结果表明: 在给定合适的最优控制目标函数下, 智能物面可最优地实时改变形状, 得到能显著提高翼型性能的最优翼型. 最优翼型在非设计工况下的气动性能也比对照翼型有 所提高.     

无网格法在三维复杂超音速流场模拟中的应用研究

对最小二乘无网格方法在含复杂外形三维超音速流场中的应用进行了研究.选用分解法求解采用最小二乘法得到的对称方程组,针对最小二乘无网格方法的计算特点生成近似正交均匀分布的离散点,对B1AC2R常规导弹超音速流场采用最小二乘无网格方法进行了无粘数值模拟,计算了B1AC2R常规导弹在不同攻角下的轴向力、法向力及俯仰力矩系数,并将数值结果与实验结果进行了比较.结果表明,最小二乘无网格方法在求解含复杂外形超音速流场时具有较高的准确度,将其应用于三维含复杂外形超音速流场的模拟是完全可行的.     

Pattern discrimination and firing behavior of cold receptor model aided with noise

A method for characterizing and identifying firing patterns of neural spike trains is presented. Three characteristic variables defined at sequential moments, including two formal derivatives and the integration of the counting process, are introduced to reflect the temporal patterns of a spike train. This paper also examines how noise interacts with encoding mechanisms of neuronal stimulus in a cold receptor. From ISI series and bifurcation diagrams it is shown that there are considerable differences in interval distributions and impulse patterns caused by purely deterministic simulations and noisy simulations. The ISI-distance can be used as an effective and powerful way to measure the noise effects on spike trains quantitatively. It is found that spike trains observed in cold receptors can be more strongly affected by noise for low temperatures than for high temperatures in some aspects; meanwhile, spike train has greater variability with the noise intensity increasing.     

Computation of 3D vortex flows past a flat plate at incidence through a variational approach of the full steady euler equations

A variational method for solving directly the full steady Euler equations is presented. This method is based on both Newton's linearization and a least squares formulation. The validity of the Euler model and boundary conditions to capture the vortex sheet is discussed. A finite element approximation of the groups of conservative variables is described and results are given for 3D subsonic flows as well as supersonic flows past a flat plate at high angle of attack.     

基于当地流活塞理论的气动弹性计算方法研究

发展了一种高效、高精度的超音速、高超音速非定常气动力计算 方法------基于定常CFD技术的当地流活塞理论. 运用当地流活塞理论计算非定常 气动力，耦合结构运动方程，实现超音速、高超音速气动弹性的时域模拟. 运用这 种方法计算了一系列非定常气动力算例和颤振算例，并和原始活塞理论、非定 常Euler方程结果作了比较. 由于局部地使用活塞理论假设，这种方法大大地克服 了原始活塞理论对飞行马赫数、翼型厚度和飞行迎角的 限制. 与非定常Euler方程方法相比，当地流活塞理论的效率很高.     

Wavelets for uncertainty estimates of propagating shock and detonation waves

The propagation speed of a shock or detonation wave in a shock or detonation tube is usually determined by a time-of-flight method by dividing the distance between two transducers with the propagation time of the disturbance signal. Some arbitrariness is inherent in determining the propagation time by this method. A new method based on Haar and Morlet wavelet transforms is reported. The method was applied to shock and detonation waves representing a step and a decaying spike discontinuity. The wavelet methods can be applied to the step discontinuity provided that the SNR ratio is good. The wavelet methods worked well for a decaying spike in the presence of noise.     

CFD analysis of an oscillating wing at various reduced frequencies

The effect of various reduced frequencies has been examined for an oscillating aspect ratio 10 NACA 0015 wing. An unsteady, compressible three‐dimensional (3D) Navier–Stokes code based on Beam and Warming algorithm with the Baldwin–Lomax turbulence model has been used. The code is validated for the study against published experimental data. The 3D unsteady flow field is simulated for reduced frequency values of 0.1, 0.2 and 0.3 for a fixed mean angle of attack position and fixed amplitude. The type of motion is sinusoidal harmonic. The force coefficients, pressure distributions and flow visualization show that at the given conditions the flow remains attached to the wing surface even at high angles of attack with no clear separation or typical light‐to‐deep category of dynamic stall. Increased magnitude of hysteresis and higher gradients are seen at higher reduced frequencies. The 3D effects are even found at midspan locations. In addition, the rate of decrease in lift near the wing tips compared with the wing root is not much like in the static cases. Copyright © 2008 John Wiley & Sons, Ltd.