超声速钝体逆向喷流减阻的数值模拟研究

为研究逆向喷流技术对超声速钝体减阻的影响,采用标准k-ε湍流模型,通过求解二维Navier-Stokes方程对超声速球头体逆向冷喷流流场进行了数值模拟,并着重分析了喷口总压、喷口尺寸对流场模态和减阻效果的影响。计算结果显示:随着喷流总压的变化,流场可出现两种流动模态,即长射流穿透模态和短射流穿透模态;喷流能使球头体受到的阻力明显减小;存在最大减阻临界喷流总压值(在所研究参数范围内最大减阻可达51.1%);在其它喷流物理参数不变时,随着喷口尺寸的增大,同一流动模态下的减阻效果下降。本文的研究对超声速钝体减阻技术在工程上的应用具有一定的参考价值。     

圆球诱发斜爆轰波的数值研究

斜爆轰发动机是飞行器在高马赫数飞行条件下的一种新型发动机,具有结构简单、成本低和比冲高等优点.但是斜爆轰发动机的来流马赫数范围广,来流条件复杂,为实现斜爆轰波的迅速、可靠引发,采用钝头体来诱发.利用Euler方程和氢氧基元反应模型,对超声速氢气/空气混合气体中圆球诱导的斜爆轰流场进行了数值研究.不同于楔面诱发的斜爆轰波,球体首先会在驻点附近诱发正激波/爆轰波,然后在稀疏波作用下发展为斜激波/爆轰波.模拟结果显示,经过钝头体压缩的预混气体达到自燃温度后,会出现两种流场:当马赫数较低时,由于稀疏波的影响,燃烧熄灭,钝头体下游不会出现燃烧情况;而当马赫数较高时,燃烧阵面能传到下游.分析表明,当钝头体的尺度较小时,驻点附近的能量不足以诱发爆轰波,只会形成明显的燃烧带与激波非耦合结构;当钝头体的尺度较大时,流场中不会出现燃烧带与激波的非耦合现象,且这一特征与马赫数无关.通过调整球体直径,获得了激波和燃烧带部分耦合的燃烧流场结构,这一流场结构在楔面诱发的斜爆轰波中并不存在,说明稀疏波与爆轰波面的相互作用是决定圆球诱发斜爆轰波的关键.     

介质与静压对激波管校准压阻式绝压传感器动态特性的影响

激波管通常用于动态压力传感器的校准,压阻式绝压传感器在激波管校准过程当中,会出现谐振频率等动态性能指标随着激波管静态压力环境、气体介质变化而改变的情况,影响传感器动态特性的校准。基于压阻式传感器的工作原理,对传感器的敏感膜片结构进行了机理分析,建立了膜片结构与校准环境中介质和静压关系的动态模型;通过ANSYS与SIMULINK软件开展了数值模拟验证工作,模拟结果与理论推导一致。通过激波管校准实验验证了气体介质与静压的影响关系,结果表明:传感器的谐振频率与静压间存在非线性关系,并且随着敏感膜片径厚比的增大而显著增大;系统阻尼比大小与气体介质有关,随着气体密度的降低而升高;传感器的灵敏度与气体介质和静压无太大直接关系。在使用激波管校准压阻式绝压传感器时,应当考虑介质与静压参数对校准结果的影响。     

功能梯度锥-柱连接壳的环向自由振动分析

论文旨在分析功能梯度锥-柱连接壳的环向自由振动,以提高其结构的振动性能和稳定性.采用Voigt模型和四参数幂函数体积分数描述功能梯度材料属性,基于Donnell薄壳理论推导出锥壳和柱壳的位移与应变关系,分别得出锥壳和柱壳的能量表达式.引入人工弹簧模拟边界和壳体间的连接条件,依据Chebyshev多项式构造位移函数,基于Rayleigh-Ritz法求解FGMs锥-柱连接壳模态频率,分析梯度指数、边界条件和几何参数对模态频率的影响.结果表明：增加陶瓷体积分数能有效提高结构的模态频率,而增大梯度指数则会降低结构的模态频率;边界约束条件越强,FGMs锥-柱连接壳的模态频率越高;随着环向波数的增大,边界条件对结构模态频率的影响越来越弱,边界约束效果作用于圆柱壳明显强于圆锥壳;当环向波数大于3时,随着壳体厚度增大,结构的模态频率呈线性提高,而增大锥柱壳长度比会降低结构模态频率;在锥柱壳长度比一定时,随着锥角的增大会使结构的模态频率先增加到峰值后减小.     

自旋挠性航天器非约束模态动力学特性研究

在研究挠性航天器动力学问题时,关注的问题是挠性航天器系统的刚柔耦合作用问题,即航天器挠性附件的振动可能会造成航天器运动失稳。针对中心刚体-双侧大挠性结构的自旋航天器,提出了航天器帆板结构的梁式简化模型,建立了一种非约束模态动力学模型。本研究考虑受到万有引力作用,探讨自旋挠性航天器非约束模态的动力学建模及动态特性。首先利用欧拉方程和哈密顿原理建立了自旋挠性航天器动力学方程,方程解释了刚性模态和弹性模态之间的耦合;然后进行了模态离散化,分别在约束模态和非约束模态下对特征值问题开展研究,对频率和相关振型进行了定量比较;最后进行了数值仿真,求解了自旋挠性航天器非约束模态特征值问题,比较约束模态与非约束模态之间的差异,并用有限元进行验证,得到了随着梁长度的增加,即刚柔惯量比、质量比的减小,非约束模态比约束模态更加准确的结论。     

弹性支撑功能梯度微圆柱壳模态频率的研究

在建立弹性支撑功能梯度薄壁微圆柱壳模型的基础上,基于修正的偶应力理论和一阶剪切变形理论,推导了微圆柱壳的模态频率方程,讨论了弹性支撑、尺寸效应、温度梯度、材料组分指数、孔隙以及几何尺寸等参数对微圆柱壳模态频率的影响。结果表明：微尺度下,弹性刚度系数在0～10⁵ N/m³范围内对微圆柱壳的模态频率基本无影响,剪切刚度系数在0～5×10⁴ N/m范围内对模态频率的影响较大,且增大剪切刚度系数有益于提高微圆柱壳的模态频率;由修正的偶应力理论得到的模态频率大于由经典连续体理论得到的模态频率;在弹性支撑和尺寸效应有无考虑的4种组合下,模态频率随温度梯度和微圆柱壳长度的增大而减小,随陶瓷体积分数指数的增大而增大,随孔隙体积分数和微圆柱壳厚度的变化规律不同;温度梯度对考虑尺寸效应或弹性基础的微圆柱壳模态频率影响较大,而孔隙调节具弹性支撑微圆柱壳的模态频率尤其显著。     

高超声速自适应激波针数值研究

针对传统的与钝体轴线共线安装的固定式激波针方法在有攻角状态所存在的问题, 在前人工作基础上得到一种新型高超声速飞行器减阻/降热方法------自适应激波针方法. 将该方法应用于三维高超声速轴对称钝锥外形以及扁平楔外形, 并采用数值模拟的方法对其进行了概念验证. 在0○~120○攻角范围内, 对不同L/D参数的激波针外形流场以及前缘壁面的压力、热流分布等进行了对比分析. 结果表明, 这种新型自适应激波针方法无论在无攻角还是有攻角状态, 均可有效降低高超声速飞行器头部壁面的压力和热流, 可以有效解决传统激波针方法在较大攻角情况状态下失效的问题.      

支杆-钝头体带攻角流场和“军刺”挡板作用研究

通过三维N-S方程的数值求解, 研究了支杆-钝头体结构在10o攻角M∞=6.0飞行条件下的流场结构和特点, 指出其气动力特性恶化的原因, 提出采用``军刺'挡板改善流场和气动力特性, 并通过对比两种不同挡板作用下的流场和气动力特性变化分析其作用机理, 发现``军刺'挡板结构分割流场抑制三维效应形成的周向流动, 迎风面形成稳定的回流区和剪切层结构, 将迎风面锥激波推离轴线, 降低钝头体肩部流动结构相互作用强度, 并在一定程度上缓解背风面流动干扰, 明显改善支杆-钝头体带攻角飞行时的气动力特性.      

自然超空泡与尾喷流相互作用的数值模拟

认识带尾喷流和自然超空泡的水下高速航行体流体动力特性并发展其预报与控制方法仍是水动力学领域极具挑战性的课题．本研究采用CFD方法对尾喷流和自然超空泡之间的相互作用进行了数值研究．针对发动机欠膨胀超音速喷流,采用现有实验结果验证了基于两方程湍流模型的二维轴对称流动数值模型的可靠性．尾喷流在空气和蒸汽环境中流动的数值计算结果表明,由于蒸汽环境中背压较低,欠膨胀尾喷流无法及时形成压缩波,使得蒸汽环境中尾喷流的过膨胀区和气相扩散区的体积比空气中大;尾喷流很难形成规则的激波格栅,波系结构相对简单．针对携尾喷流的平头航行体超空泡流状态的数值模拟结果表明,尾喷流注入超空泡后迅速充满航行体周围的空腔区域;尾喷流与超空泡尾迹区域形成的回射流相互作用最终导致超空泡断裂,断裂过程中伴随着燃气泡的下泄现象;受空泡壁面约束,尾喷流难以在狭窄的超空泡空腔内完全膨胀,尾喷流的激波波系结构有显著的变化：在喷嘴附近受到压缩,在远离喷嘴区域受到超空泡水汽掺混的破坏;空泡内压强基本维持在饱和蒸汽压附近,没有显著增加．     

矢量喷流对细长体大迎角非对称流动影响研究

采用测压方法研究了矢量喷流对细长旋成体大迎角非对称流动的影响特性.实验结果表明:矢量喷流对细长旋成体大迎角非对称侧向力有明显的抑制作用,该抑制作用是通过喷流诱导作用,改变其空间绕流涡系结构的分布来实现的,但是矢量喷流的存在并不能改变大迎角机身空间绕流涡系的本质结构;随着迎角的增大,矢量喷流对细长旋成体大迎角非对称流动的影响区域不断前移,甚至影响到头部;随着喷流落压比的增加,矢量喷流对细长旋成体大迎角非对称侧向力的抑制作用加强,但当喷流落压比达到临界落压比后(即喷管出口处达到设计马赫数时),喷流影响作用将不会随喷流落压比的增加而改变.     

Numerical heat transfer study around a spiked blunt-nose body at Mach 6

An aerospike attached to a blunt body significantly alters its flowfield and influences aerodynamic drag at high speeds. The dynamic pressure in the recirculation area is highly reduced and this leads to the decrease in the aerodynamic drag. Consequently, the geometry of the aerospike has to be simulated in order to obtain a large conical recirculation region in front of the blunt body to get beneficial drag reduction. Axisymmetric compressible Navier–Stokes equations are solved using a finite volume discretization in conjunction with a multistage Runge–Kutta time stepping scheme. The effect of the various types of aerospike configurations on the reduction of aerodynamic drag is evaluated numerically at a length to diameter ratio of 0.5, at Mach 6 and at a zero angle of incidence. The computed density contours are showing satisfactory agreement with the schlieren pictures. The calculated pressure distribution on the blunt body compares well with the measured pressure data on the blunt body. Flowfield features such as formation of shock waves, separation region and reattachment point are examined for the flat-disc spike and on the hemispherical disc spike attached to the blunt body. One of the critical heating areas is at the stagnation point of a blunt body, where the incoming hypersonic flow is brought to rest by a normal shock and adiabatic compression. Therefore, the problem of computing the heat transfer rate near the stagnation point needs a solution of the entire flowfield from the shock to the spike body. The shock distance ahead of the hemisphere and the flat-disc is compared with the analytical solution and a good agreement is found between them. The influence of the shock wave generated from the spike is used to analyze the pressure distribution, the coefficient of skin friction and the wall heat flux facing the spike surface to the flow direction.     

Shock unsteadiness creation and propagation: experimental analysis

The possibility of creating unsteady distortions of the tip shock by waves emitted from an aircraft is assessed experimentally. The model chosen is a cylindrical fore body equipped with a spike. This configuration is known for generating an important level of unsteadiness around the spike in supersonic regime. The wind tunnel Mach number is equal to 2. The experiments show that waves emitted from this source propagate along the tip shock and interact with it. It is then assessed that this interaction produces a periodic distortion of the shock that propagates to the external flow. Unsteady pressure sensors, high speed schlieren films, hot wire probing and laser Doppler velocimetry are used as complementary experimental means. The final result is a coherent representation of the complex mechanism of wave propagation that has been evidenced. The principle of creating unsteady shock deformation by onboard equipments could be examined as a possibly promising method of sonic boom control.     

Characteristics of unsteady type IV shock/shock interaction

Characteristics of the unsteady type IV shock/shock interaction of hypersonic blunt body flows are investigated by solving the Navier–Stokes equations with high-order numerical methods. The intrinsic relations of flow structures to shear, compression, and heating processes are studied and the physical mechanisms of the unsteady flow evolution are revealed. It is found that the instantaneous surface-heating peak is caused by the fluid in the “hot spot” generated by an oscillating and deforming jet bow shock (JBS) just ahead of the body surface. The features of local shock/boundary layer interaction and vortex/boundary layer interaction are clarified. Based on the analysis of flow evolution, it is identified that the upstream-propagating compression waves are associated with the interaction of the JBS and the shear layers formed by a supersonic impinging jet, and then the interaction of the freestream bow shocks and the compression waves results in entropy and vortical waves propagating to the body surface. Further, the feedback mechanism of the inherent unsteadiness of the flow field is revealed to be related to the impinging jet. A feedback model is proposed to reliably predict the dominant frequency of flow evolution. The results obtained in this study provide physical insight into the understanding of the mechanisms relevant to this complex flow.     

逆向喷流流场模态分析及减阻特性研究

逆向喷流减阻的基本原理是利用逆向高速喷流与飞行器绕流的相互作用，使飞行器周围的流场结构发生变化，致使飞行器的气动特性发生改变，从而改善飞行器的气动性能。利用数值模拟方法对轴对称球头、截锥的逆向喷流流场开展了研究，考虑了高温非平衡化学反应对流场的影响。模拟了球头和截锥在不同总压比时流场不同的模态：长穿透流模态（LPM）和短穿透流模态（SPM），得到了不同模态下钝体表面压力、气动力系数和不同模态之间转换的瞬态效应．简单分析了喷流在减阻方面的应用，给出了几个喷口参数与减阻效率之间的关系，提出了喷流减阻工程应用时应考虑的主要因素。     

Numerical analysis of pressure oscillations over axisymmetric spiked blunt bodies at Mach 6.80

The pressure oscillations over a forward facing spike attached to an axisymmetric blunt body are simulated by solving time-dependent compressible Navier–Stokes equations. The governing fluid flow equations are discretized in spatial coordinates employing a finite volume approach which reduces the equations to semidiscretized ordinary differential equations. Temporal integration is performed using the two-stage Runge–Kutta time stepping scheme. A global time step is used to obtain a time-accurate numerical solution. The numerical computation is carried out for a freestream Mach number of 6.80 and for spike length to hemispherical diameter ratios of 0.5, 1.0 and 2.0. The flow features around the spiked blunt body are characterized by a conical shock wave emanating from the spike tip, a region of separated flow in front of the hemispherical cap, and the resulting reattachment shock wave. Comparisons of the numerical results are made with the available experimental results, such as schlieren pictures and the surface pressure distribution along the spiked blunt body. They are found to be in good agreement. Spectral analysis of the computed pressure oscillations are performed employing fast Fourier transforms. The surface pressure oscillations over the spike and phase plots exhibit a behaviour analogous to that of the Van der Pol equation for a self-sustained oscillatory flow. Received 28 February 2001 / Accepted 17 January 2002     

Disturbance evolution in the leading-edge region of a blunt flat plate in supersonic flow

The spatio-temporal dynamics of small disturbances in viscous supersonic flow over a blunt flat plate at freestream Mach number M_∞=2.5 is numerically simulated using a spectral approximation to the Navier–Stokes equations. The unsteady solutions are computed by imposing weak acoustic waves onto the steady base flow. In addition, the unsteady response of the flow to velocity perturbations introduced by local suction and blowing through a slot in the body surface is investigated. The results indicate distinct disturbance/shock-wave interactions in the subsonic region around the leading edge for both types of forcing. While the disturbance amplitudes on the wall retain a constant level for the acoustic perturbation, those generated by local suction and blowing experience a strong decay downstream of the slot. Furthermore, the results prove the importance of the shock in the distribution of perturbations, which have their origin in the leading-edge region. These disturbance waves may enter the boundary layer further downstream to excite instability modes.     

Supersonic motion of bodies in a gas with shock waves

The authors consider the problem of supersonic unsteady flow of an inviscid stream containing shock waves round blunt shaped bodies. Various approaches are possible for solving this problem. The parameters in the shock layer on the axis of symmetry have been determined in [1, 2] by using one-dimensional theory. The authors of [3, 4] studied shock wave diffraction on a moving end plane and wedge, respectively, by the through calculation method. This method for studying flow around a wedge with attached shock was also used in [5]. But that study, unlike [4], used self-similar variables, and so was able to obtain a clearer picture of the interaction. The present study gives results of research into the diffraction of a plane shock wave on a body in supersonic motion with the separation of a bow shock. The solution to the problem was based on the grid characteristic method [6], which has been used successfully to solve steady and unsteady problems [7–10]. However a modification of the method was developed in order to improve the calculation of flows with internal discontinuities; this consisted of adopting the velocity of sound and entropy in place of enthalpy and pressure as the unknown thermodynamic parameters. Numerical calculations have shown how effective this procedure is in solving the present problem. The results are given for flow round bodies with spherical and flat (end plane) ends for various different values of the velocities of the bodies and the shock waves intersected by them. The collision and overtaking interactions are considered, and there is a comparison with the experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 141–147, September–October, 1984.     

Numerical investigation of supersonic flow past blunt bodies with protruding spikes

A procedure for the calculation of a supersonic flow of ideal gas near axisymmetric blunt bodies with protruding spikes is developed. The flow past a frustum of a cone with a protruding spherically blunt cylindrical spike as a dependence on the ratio K of the spike length1 to the diameter D of the flat end of the body and the Mach number M of the oncoming flow is studied. Several steady flow regimes are obtained, including the formation of circulation zones and internal shock waves in the shock layer. It is shown that mounting a spike in front of the frustum of a cone can lead to a 40–50% reduction in its drag. A full investigation of the variation of the drag coefficient as a dependence on K is carried out for M = 3.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 119–127, May–June, 1986.The authors express their gratitude to V. A. Levin for the formulation of the problem and his constant attention to the work.     

Interaction between an external compression shock and a blunt body in a hypersonic stream

A complex shock configuration with two triple points can occur during the interaction between an external oblique compression shock and the detached shock ahead of a blunt body (for instance, ahead of a wing or stabilizer edge). This results in the formation of a high-pressure, low-entropy supersonic gas jet [1–6]. Here two flow modes are possible [1], which differ substantially in the intensity of the thermal and dynamic effects of the stream on the blunt body: mode I corresponds to the impact of a supersonic jet [2–6], while the supersonic jet in mode II does not reach the body surface in the domain of shock interaction because of curvature under the effect of a pressure drop. Conditions for the realization of the above-mentioned flow modes are investigated experimentally and theoretically, and an approximate method is proposed to determine the magnitude of the compression shock standoff in the interaction domain. Blunt bodies with plane and cylindrical leading edges are examined. The results of a computation agree satisfactorily with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 97–103, January–February, 1976.The author is grateful to V. V. Lunev for discussing the research and for useful remarks.