单喷管火箭自由喷流噪声数值模拟

利用非定常燃气流场数值模拟方法，研究超声速喷流核心区及其外围亚声速区燃气流动态分布和流动特性，在此基础上利用欧拉方程条件的伽辽金有限元方法以及FWH方法，实现并完成喷流噪声传播特性、辐射特性数值模拟。喷流噪声数值模拟结果显示：燃气流推进初期，强喷流噪声区域紧随燃气流前锋；燃气流场相对稳定后，强喷流噪声区域主要位于燃气流等能区末稍，一些小尺度试验的燃气流激波系附近也存在较强喷流噪声。这些强喷流噪声主要由燃气流前锋带动的大涡动态卷吸、燃气流强湍流脉动以及激波扰动引起。受数值模拟网格分辨率影响，当前仅能保持中低频段声压级数值模拟结果与实测结果总体接近。     

Plasma Jet Motion Across the Geomagnetic Field in the “North Star” Active Geophysical Experiment

The active geophysical rocket experiment North Star was carried out in the auroral ionosphere on January 22, 1999, at the Poker Flat Research Range (Alaska, USA) using the American research rocket Black Brant XII with explosive plasma generators on board. Separable modules with scientific equipment were located at distances of from 170 to 1595 m from the plasma source. The experiment continued the series of the Russian–American joint experiments started by the Fluxus experiment in 1997. Two injections of aluminum plasma across the magnetic field were conducted in the North Star experiment. They were different, since in the first injection a neutral gas cloud was formed in order to increase the plasma ionization due to the interaction of neutrals of the jet and cloud. The first and second injections were conducted at heights of 360 and 280 km, respectively. The measurements have shown that the charged particle density was two orders of magnitude higher in the experiment with the gas release. The magnetic field in the first injection was completely expelled by the dense plasma of the jet. The displacement of the magnetic field in the second injection was negligible. The plasma jet velocity in both injections decreased gradually due to its interaction with the geomagnetic field. One of the most interesting results of the experiment was the conservation of high plasma density during the propagation of the divergent jet to considerable distances. This fact can be explained by the action of the critical ionization velocity mechanism.     

数值模拟侧向超声速单喷流干扰流场特性

采用数值方法研究了平板上超/高超声速来流与超声速横向喷流相撞引起的复杂干扰流场特性。所建立的单介质冷喷流数值模拟方法,经过了表面多方位压力分布测量结果、纹影显示的激波结构以及表面油流图谱表现的表面分离范围的实验验证。根据数值模拟与实验对比的结果,合理地描述了喷流干扰流场压力分布以及表面、空间结构特性,并分析了压力比对流场结构和特性的影响。     

Integrated Studies of Electric Propulsion Engines during Flights in the Earth’s Ionosphere

Fifty years ago, on October 1, 1966, the first Yantar satellite laboratory with a gas plasma–ion electric propulsion was launched into orbit as part of the Yantar Soviet space program. In 1966–1971, the program launched a total of four laboratories with thrusters operating on argon, nitrogen, and air with jet velocities of 40, 120, and 140 km/s, respectively. These space experiments were the first to demonstrate the long-term stable operation of these thrusters, which exceed chemical rocket engines in specific impulse by an order of magnitude and provide effective jet charge compensation, under the conditions of a real flight at altitudes of 100–400 km. In this article, we have analyzed the potential modern applications of the scientific results obtained by the Yantar space program for the development of air-breathing electric propulsion that ensure the longterm operation of spacecraft in very low orbits.     

Hypervelocity fragment formation technology for ground-based laboratory tests

Cumulative jet formation was regarded aimed at the formation of hypervelocity fragments up to 14 km/s for the investigation of space debris effects on shielding screens. The basic physical problems of jet formation process are analyzed. It is shown that in process of realization of hyper-cumulation conditions for jet formation without complete stagnation point involving formation of the inner zone of constant pressure (dead zone), the flow mass is always greater than slug mass. That opens wide the possibilities for increasing fragment mass in laboratory tests.     

导弹轨控喷流气动干扰全空域插值方法研究

导弹侧向喷流气动干扰受海拔高度影响显著,主要原因在于雷诺数和压力比会随海拔高度的变化而变化。受大气参数的影响,雷诺数和压力比的变化规律并不一致。为能在现有气动模型基础上既保证插值精度又减少节点数,通过数值模拟方法研究了雷诺数和压力比对轨控喷流气动干扰的影响。结果表明:海拔高度对喷流气动干扰的影响主要由压力比变化造成,雷诺数变化产生的影响相对较小,且喷流干扰与压力比间存在一定的线性关系。因此,提出了一种基于压力比的喷流气动干扰插值方法。相比于传统的基于高度的插值方法,该方法可显著提高插值精度,并减少节点数,适用于全空域轨控喷流气动建模。     

三维横向流喷与超声速主流干扰流场的数值模拟

本文对三维隅角机翼声速喷流与超声速主流的干扰流扬进行了数值模拟。三维欧拉方程的求解采用非结构网络有限体积伽辽金法（Finite Volume Galerkin Method）。引入了总体结点积分域的概念,简化了从单元矩阵到总体矩阵的汇总过程。通量的分裂采用Osher格式,通过外差使其由一阶精度上升为二阶精度。发展了一种基于线化流量的逆风非结构网格隐式有限元格式以提高求解精度及效率。最后给出了三维隅角机翼流场的算例。     

侧向喷流干扰流场建立与消退过程数值模拟

针对高超声速轨控侧向喷流的非定常干扰效应问题,应用非定常数值模拟方法和多重网格加速收敛技术,研究了锥-柱-裙外形轨控侧向喷流干扰流场的建立和消退过程,获得了详细的喷流瞬时干扰流场结构特性,分析了法向力放大系数、干扰力矩系数、法向力系数及俯仰力矩系数随时间的变化特性。研究结果表明：在侧向喷流干扰流场建立和消退过程中气动力变化较大,存在峰谷值;法向力放大系数及干扰力矩系数的定常值和非定常时均值之间存在明显差异。     

类火星环境中火箭起降过程喷流影响研究

采用尺度自适应模拟（Scale Adaptive Simulation, SAS）方法对火箭在起降过程中喷流的地面效应进行了数值研究，对比分析了在海平面和类火星环境下的喷流冲击倾斜地面和水平地面时的流场结构以及喷管的侧向载荷。在海平面环境下，喷管内部出现流动分离，喷管内壁压力先减小，在分离点处压力达到最低点，然后压力增加至环境压力。喷流冲击倾斜地面时具有更强烈的流动不稳定性，喷流出现明显的时间演化现象。在类火星环境下，喷管内部不发生流动分离，内壁的压力沿着流向逐渐减小。喷流冲击倾斜地面时，在喷管出口的外侧形成了上下移动的马赫盘；喷流冲击水平地面时，喷流和地面反弹气流相互作用，阻碍了马赫盘的形成。研究结果还表明在类火星环境和海平面环境下，喷流冲击倾斜地面时喷管内壁温度均上升。研究结果将为火箭回收和火星着陆的方案选取提供参考。     

Formation of a Polarization Jet during the Expansion Phase of a Substorm: Results of Ground-Based Measurements

The ground-based polarization jet measurements at the Yakutsk ionosphere station (L= 3.0) for the years 1989–1991 (110 events) are compared with variations of the AE-index and with parameters of the local magnetic activity. It is shown that polarization jet development in the near midnight sector can occur during a period of no longer than 10 min on the expansion phase of a substorm. The formation of the polarization jet is accompanied by a specific magnetic field variation corresponding in shape to a fast passage of the Harang Discontinuity above the station. Statistical data are given on ground level observations of the polarization jet, which are close to those measured from satellites. The mean delay (averaged over the full data bank) between the onset of a substorm with AE 500 nT and the moment of the polarization jet appearance at L= 3.0 is equal to 0.5 h near midnight and to 1.0 – 1.5 h in the evening sector. Estimations show that the duration of the polarization jet formation when energetic ions are injected into the Harang Discontinuity region above the ground station can last for about 10 min, and during this time the Harang Discontinuity can be shifted to the west. This is in qualitative agreement with the described observations.     

高层大气模型对空间站轨道漂移和寿命的影响分析

本文以轨道摄动分析方法一阶理论为基础，其中大气阻力摄动采用数值积分方法，给出一种可利用各种大气模型进行轨道摄动分析的计算方法，并利用三种高层大气模型（ＣＩＲＡ—７２，ＣＩＲＡ—８６和ＤＴＭ）和三个太阳活动水平（Ｆ１０．７＝１００，１５０和２００）分析比较了大气阻力振动对高度为４００ｋｍ的空间站轨道漂移和寿命的影响，以及估算修正轨道漂移所需的能量。给出的定量分析结果将为空间站或航天飞行器的轨道设计和能量估算提供依据。     

两种物态方程在超高速碰撞数值模拟中的比较

材料在超高速碰撞下将发生熔化甚至气化相变。为反映超高速碰撞中材料相变带来的影响,文章采用GRAY三相物态方程与Tillotson物态方程,对超高速碰撞进行数值模拟对比研究。研究表明:当碰撞速度在3 km/s以下时,除了靶板穿孔直径外,两种物态方程所给出的碎片云的结果基本一致;但当碰撞速度在3 km/s以上时,两种物态方程给出的数值模拟结果有较大差异,这说明在超高速碰撞中相变的产生对碎片云形状参数有较大影响。     

超声速横向气流中喷雾的数值模拟

对超声速横向气流中的喷雾过程进行了数值模拟,采用二维N-S方程计算气相,应用一次雾化模型和二次雾化模型模拟了喷雾雾化过程,并与实验测量结果进行了对比。研究了湍流度和附面层厚度对液雾穿透深度的影响,发现湍流度和附面层厚度并不是主要的影响因素,认为雾化模型是影响液雾穿透深度的关键因素。     

Shock–Acoustic Waves Generated during Rocket Launches and Earthquakes

Shock–acoustic waves generated during rocket launches and earthquakes are investigated by a method developed earlier for processing data from a global network of receivers of the GPS navigation system. Disturbances of the total electron content in the ionosphere accompanying the launches of the Proton, Soyuz, and Space Shuttle space vehicles from the Baikonur cosmodrome and Kennedy Space Center launch site in 1998–2000, as well as the earthquakes in Turkey on August 17 and November 12, 1999, were analyzed. It was shown that, regardless of the source type, the impulsive disturbance has the character of an N-wave with a period of 200–360 s and an amplitude exceeding background fluctuations under moderate geomagnetic conditions by a factor of 2–5 as a minimum. The elevation angle of the disturbance wave vector varies from 25° to 65°, and the phase velocity (900–1200 m/s) approaches the speed of sound at heights of the ionospheric F-region maximum. The source location corresponds to a segment of the booster trajectories at a distance of no less than 500–1000 km from the start position and to a flight altitude of no less than 100 km. In the case of earthquakes the source location approximately coincides with the epicenter.     

有限空间轴对称射流流场的数值模拟

本文在к－ε紊流模型下，采用罚有限元法对有限空间轴对称射流流场进行数值模拟。引用罚函数优化迎风格式和退化积分，建立有限元解法的数值矩阵，在不同截面上求得工，压力和紊流参数分布，计算结果与实验数据十分接近。本文研究结果可应用于火箭发动机的喷射设计和射流量论的研究。     

Changes in the Earth’s ozonosphere due to ionization of high-latitude atmosphere by solar protons in October, 2003

Using the data of the Russian KORONAS-F satellite and American GOES spacecraft on solar cosmic ray fluxes associated with powerful events which occurred on the Sun at the end of October - the beginning of November, 2003, calculations of ionization of high-latitude (70° N) atmosphere were carried out. The calculations have shown that the maximum values of ionization for the chosen latitude lie in the range of 50–70 km. The largest ionization was caused by the flare on November 28, 2003. Based on a numerical photochemical simulation it is shown that, as a result of intensification of catalytic cycles with participation of ozone-destroying NO and OH, the concentration of ozone decreased by 30% at ionization maximum altitudes.Translated from Kosmicheskie Issledovaniya, Vol. 42, No. 6, 2004, pp. 653–662.Original Russian Text Copyright © 2004 by Krivolutsky, Kuminov, Vyushkova, Kuznetsov, Myagkova.     

SMC模式下RBCC发动机4Ma工况性能仿真

为研究SMC模式下火箭混合比对RBCC发动机性能的影响规律,完成了氢/氧火箭推力室中心布局、二元定几何结构模型发动机飞行马赫数Ma0=4、高度H=17 km弹道点流场仿真,获得了不同火箭混合比(MR=2、3、4、5、6、8)及燃烧室长度的推力、比冲性能.研究表明:在火箭燃气富燃条件下(MR<8),产生了正的火箭推力增益...     

快速绕飞卫星空间圆编队设计方法

连续小推力条件下,针对圆参考轨道卫星,推导了满足快速绕飞条件的空间圆编队动力学模型。分析了一个绕飞周期内的燃耗情况。在绕飞周期确定的条件下,给出了最小燃耗的计算方法,并利用数值方法验证了所推导公式的正确性。计算结果表明,在高度为500km的圆轨参考道卫星实现周期为10分钟、半径100m的空间圆绕飞,一个绕飞周期所需最小速度增量约为4.77m/s。文章提出的基本设计原理并不局限于空间圆编队,也可用于其它快速绕飞编队的设计工作。
     

High speed gas flows in explosions of cavitated explosives

A possibility of attaining steady flow of detonation products with specific energy much larger than the specific chemical energy of explosive is demonstrated in the case when a cylindrical charge of explosive is fitted with an evacuated cavity. Simple estimates and results of numerical analysis of the process are presented. Steady process may be considered to occur under the following assumptions: (1) effects arising due to jet interaction with cavity walls are negligible; (2) the detonation process is steady. In the case of limited explosive lengths these assumptions have been shown to be correct.When the cavity is filled with gas or liquid, a variety of steady and non-steady flow regimes is possible, depending on the density of the filling medium. One well-known case is that of flow with irregular reflection of shock waves at the cavity axis accompanied with the formation of Mach intersections. Another interesting flow regime is observed to occur in the case of low density filling medium (liquid hydrogen, for example). In this case the filling medium is driven by a “detonation piston” at constant velocity, equal to the velocity of detonation, forming a uniform growing column of hot shock-compressed matter, specific energy of which exceeds by one order of magnitude the specific energy of the explosive. Obviously, the walls of the vessel containing hydrogen must be able to withstand radial loads for a sufficiently long time (20 μ sec).The relative merits of these methods in comparison to others in high speed gas-dynamics is discussed.