超临界压力下航空煤油RP-3在水平细圆管内对流换热特性实验研究

为探究超临界压力下碳氢燃料在水平管内的对流换热规律,文章针对超临界条件下航空煤油RP-3在水平细圆管内的对流换热,分析了热流密度、进口雷诺数及浮升力对对流换热的影响。研究表明:沿流动方向,管内表面传热系数随热流密度的增大先减小后增大;在低进口温度及低进口雷诺数情况下,管内换热均出现先恶化后强化的现象,而随着进口温度和雷诺数的增加,此现象消失;浮升力对换热的影响随热流密度的增加而增加;浮升力对下表面换热的加强使得入口效应的影响在下表面先于上表面结束;受浮升力影响,上下壁最大温差可达50 K;质量流速的增加会抑制浮升力对换热的影响;准则数Grq/Grth可以很好地反映浮升力的变化趋势。以上研究结果可为采用碳氢燃料作冷却介质的各类飞行器主动热防护技术方案提供技术支撑。     

Sources and Losses of Energetic Electrons at Low Altitudes from Satellite Data

Fluxes and differential spectra of 0.04–2.0 MeV electrons at low altitudes in the Earth's magnetosphere are considered in comparison with the model spectra AE-8 MAX and AE-8 MIN. Possible causes of the discrepancy between the observational and model spectra are discussed. The coefficients of radial diffusion at various L-shells are estimated for the maximum of solar activity (using the Interkosmos-19 data) and for the minimum of solar activity (using the Kosmos-1686 data) and are derived from the model AE-8. A quantitative evaluation of the electron yield from radial diffusion at low L shells is derived. Ionization losses, Coulomb angle scattering, and resonant wave–particle interaction are considered as the loss mechanisms. A calculation of these losses at the low L-shells is given. The electron distribution at low L-shells is best fitted by a combination of dissipative terms from different models: Coulomb scattering dominates at the lower L-shells (L = 1.2–1.4) and the resonant wave–particle interaction controls the radiation belt maximum and the gap (L = 1.4–2.0).     

Determination of the Turbulent Diffusion Coefficient in the Plasma Sheet Using the Project INTERBALL Data

The results of an analysis of velocity fluctuations in the plasma sheet of the Earth's magnetotail measured onboard INTERBALL Tail Probe satellite are presented. The hodographs of the velocity in directions (Y, Z) and correlation functions are presented for a number of passages when the satellite was in the plasma sheet for a long time. The turbulent diffusion coefficients are calculated. A comparison of the obtained diffusion coefficients with those predicted theoretically in [1] is carried out. It is shown that the results of observations confirm theoretical predictions.     

Mass transfer interaction of moving particles in a reactive dispersive system

Mass transfer from a system of particles moving linearly in a viscous fluid is considered. A chemical reaction on the particle surfaces is taken into account. The analysis is based on the investigation of concentration and velocity fields in a particle wake and of the distortion of these fields due to particles interaction. The values of both local and average Sherwood numbers are calculated.The results obtained are applied to the system of drops and solid spherical particles moving in line along the axis of uniform Stokes flow. It is shown, that if the distance 1 between the particles satisfies the condition (a is the sphere radius and P is the Peclet number) then the total diffusion flux to the surface of kth sphere is as follows

Full-size image (1K)

. Here k is the particle number along the stream (for the leading particle k = 1), I₁ is the diffusion flux to the leading particle: n = 1 for drops, n = 2 for solid particles.     

Verification of the Accuracy of the Magnetopause Empirical Models Using the INTERBALL-1 Satellite Data

A statistical analysis of the shape and location of the magnetopause according to the INTERBALL-1 satellite data for the period 1995–1997 is carried out. The instants of crossing the magnetosphere boundaries obtained by the plasma and magnetic data are compared with computations based on three empirical models, namely, Petrinec and Russel, 1996; Shue et al., 1997; and Shue et al., 1998. The state of the interplanetary medium (dynamic pressure of the solar wind plasma P _d and the B _z component of the interplanetary magnetic field) was determined by the measurements onboard the WIND spacecraft. We estimate the accuracy of the considered models for different groups of boundary crossings: single, multiple with small duration (less than 40 min), and multiple with large duration (more than 40 min). It is demonstrated that the small-scale motions of the boundary (<1R _E) are observed more often in the dayside magnetosphere, especially near the cusp region. Large-scale boundary oscillations (>1R _E) are more common in the tail region of the magnetosphere, namely, its flanks. Various models give similar results: about 50% of all events have deviations by more than 1R _E from the model locations. In some cases, the deviation of the measured location of the magnetosphere boundary from the model prediction may be as large as 5–6R _E for all three models considered, the actual boundary being more often located nearer to the Earth than the result of model computations. The best model is that of Shue et al., 1998, but it does not differ significantly from the other models.     

Turbulence in the Plasma Sheet during Substorms: A Case Study for Three Events Observed by the INTERBALL Tail Probe

Variations of turbulence properties of the plasma sheet during geomagnetic substorms are investigated using observations of the INTERBALL Tail Probe satellite. The periods are chosen when the satellite was inside the plasma sheet. Fluctuations of the plasma bulk velocity across the plasma sheet are studied for the growth, expansion, and recovery phases of geomagnetic substorms on October 14, 1997; October 30, 1997; and December 16, 1998. It was demonstrated that the level of turbulence increases considerably after the onset of the substorm expansion phase and slowly decreases nearly to the presubstorm level later. The correlation times of plasma fluctuations in the Z-direction are estimated, and diffusion coefficients in the Z-direction are calculated.     

A study of detonation in methane/air clouds

In the present paper, we study the problem of detonation in unconfined, gaseous mixtures of methane/oxygen/nitrogen. A numerical simulation approach is employed in which we use a one-dimensional (spherical symmetry), time-dependent computer model to simulate the coupled compressible fluid dynamics-chemical kinetics processes which occur upon direct initiation of detonation. We establish the magnitude of explosive yield of tetryl required to initiate detonation in mixtures of CH₄ + 2O₂ with varying degrees of nitrogen dilution, up to and including stoichiometric . The numerical simulations illustrate the features of direct initiation observed in many experimental investigations, e.g. shock-wave breakaway followed by detonation reestablishment via a quasi-steady, oscillatory flow regime which occurs before the establishment of a steadily propagating spherical detonation. Our results compare well with recent experimental data obtained by Bull et al. (1976) over the range of tetryl masses studied by them. We find that tetryl explosive masses in excess of 10⁷ grams would be required to initiate detonation in an unconfined, stoichiometric mixture.     

一种半解析方法火星再入探测器一体化设计模型

提出了一种火星再入探测器的半解析法一体化模型,为未来火星探测器的初步设计提供包括弹道、气动参数、防热系统在内的设计依据。在任务初步设计阶段,由于涉及多个学科,火星大气与地球大气模型差异很大,需要考虑的因素很多,因此需要一个半解析的,能够实时进行弹道、气动参数和防热层温度分布计算的设计工具。该工具应该包括气动外形、弹道、防热为一体的设计模型,且能满足计算负担要求。基于此,提出了一种普遍适应的火星探测器外形定义方程,并对其中的一种特殊情况进行了解析公式的推导。在合理假设的前提下,给出了计算防热层内部温度分布的半解析计算公式,并以MER\|B (Mars Exploration Rover)探测器为例进行了弹道和气动受热仿真。计算结果表明,该一体化计算模型的仿真计算落点与实际落点很接近,达到了火星探测器初步设计阶段一体化的计算要求。
     

Evaporation and ignition of droplets in combustion chambers modeling and simulation

Computer simulation of liquid fuel jet injection into heated atmosphere of combustion chamber, mixture formation, ignition and combustion need adequate modeling of evaporation, which is extremely important for the curved surfaces in the presence of strong heat and mass diffusion fluxes. Combustion of most widely spread hydrocarbon fuels takes place in a gas-phase regime. Thus, evaporation of fuel from the surface of droplets turns to be one of the limiting factors of the process as well. The problems of fuel droplets atomization, evaporation being the key factors for heterogeneous reacting mixtures, the non-equilibrium effects in droplets atomization and phase transitions will be taken into account in describing thermal and mechanical interaction of droplets with streaming flows. In the present paper processes of non-equilibrium evaporation of small droplets will be discussed. As it was shown before, accounting for non-equilibrium effects in evaporation for many types of widely used liquids is crucial for droplet diameters less than 100 μm, while the surface tension effects essentially manifest only for droplets below 0.1 μm. Investigating the behavior of individual droplets in a heated air flow allowed to distinguish two scenarios for droplet heating and evaporation. Small droplets undergo successively heating, then cooling due to heat losses for evaporation, and then rapid heating till the end of their lifetime. Larger droplets could directly be heated up to a critical temperature and then evaporate rapidly. Droplet atomization interferes the heating, evaporation and combustion scenario. The scenario of fuel spray injection and self-ignition in a heated air inside combustion chamber has three characteristic stages. At first stage of jet injection droplets evaporate very rapidly thus cooling the gas at injection point, the liquid jet is very short and changes for a vapor jet. At second stage liquid jet is becoming longer, because evaporation rate decreases due to decrease of temperature. But combustion of fuel vapor begins which brings to increase of heat flux to droplets and accelerates evaporation. The length of the liquid jet decreases again and remains constant slightly oscillating.     

A Diffusion Model of Radial Distributions of Energetic Protons in Saturn's Magnetosphere

A simple diffusion model of Saturn's radiation belts with a constant internal source and particle losses at the boundaries is considered. The radial diffusion coefficient is supposed to be D _LL L ⁴. The analytical solution to the stationary diffusion equation describes well the observed radial profiles of the distribution function of energetic protons in Saturn's magnetosphere.     

The turbulent combustion of a propane-air mixture

This paper models the combustion of a turbulent homogeneous mixture of propane and air within a duct having a stationary one-dimensional mean flow. The Bray-Moss model is applied to the closure of the chemical production terms, using a probability density function (pdf) of the temperature which is chosen as the characteristic variable. Under the conditions chosen for the study, chemical kinetic factors are important and the conventional assumption, that heat release is controlled by turbulent mixing, is not valid. The chemical model of Edelman and Fortune for the combustion of hydrocarbons is used and simplifying assumptions are made which reduce the systems of unknowns to that of the temperature alone. This leads to the introduction of two chemical production terms which are defined respectively in a “delay zone”, where the heat release is modest, and a “combustion zone”. The required equations for the Favre-averaged temperature, turbulence kinetic energy and the mean square fluctuation of the temperature are solved numerically. In the delay zone, a comparison is made between a second order Borghi type closure and the pdf closure. Good agreement is found in the case of relatively small turbulence intensity. It is shown that the pdf formulation does not require the two zones to be spatially distinct. Differing chemical source terms can be discriminated instantaneously by the reaction progress variable and contributions to the average production terms appropriately apportioned by its pdf. Predictions are made of the profiles of mean temperature and mean square fluctuation under different initial turbulence levels.     

吸波外热流模拟装置等效热分析模型

吸波材料可用于微波载荷无线功能测试。在微波载荷的真空热试验条件下,容器壁安装由大量四棱锥尖劈阵列而成的吸波结构,其兼具外热流模拟功能。为简化热分析以指导工程优化设计,文章基于体积平均思想,为吸波结构建立等效热分析模型,并与真实结构模型的分析结果进行对比,研究吸波装置的瞬态温度变化和外热流模拟精度以及对于复杂结构的温度场控制精度,验证了等效模型的精度与工程实用性。     

Prediction of turbulent diffusion flames with a four-equation turbulence model

For turbulent flame prediction a four-equation turbulence model consisting of the turbulent kinetic energy, the dissipation rate and the density-velocity correlations using unweighted statistics is suggested. The fluctuating concentrations are taken into account by using Spalding's (1971) equation for the variance of the mixture fraction f and by defining the probability density function of f as beta-function. The thermodynamic model for the combustion of H₂ with air assumes infinitely fast reaction in a global single step. Preliminary results are shown for the calculation of Kent and Bilger (1972) turbulent H₂-air diffusion flame. For this flame with fluctuating density the inclusion of the velocity-density correlations appears to be essential, because they have a strong effect especially on the development of the turbulent viscosity and the mean velocity and mass fraction profiles. Consequently the mean concentrations are in good agreement with the experimental results.     

氧气/煤油发动机燃气热物理参数及输运系数计算

为了给氧气/煤油发动机设计和热防护设计提供必要的设计参数,针对氧气/煤油燃气进行热力学计算。运用吉布斯最小自由焓计算模型得到燃气平衡组成,通过拟合公式的方法得到燃气的热物理参数及输运系数。通过计算,得到氧气/煤油燃气的组分及比焓、密度、比熵、粘性系数等热物理参数和输运系数随温度和压力的变化特性。分析结果表明:水离解对氧气/煤油燃气组分变化存在显著影响,压力增大会导致水离解起始温度升高;氧气/煤油燃气比焓、比熵、定压比热、粘性系数、热传导系数变化在温度较低时受压力影响较小,当水开始离解后,压力的影响显著增强;组分在燃气中的扩散系数同时受到了温度和组分摩尔分数的影响;燃气普朗特数变化受热传导系数变化的影响较大,水离解后,热传导系数的迅速增大使燃气的普朗特数迅速减小。     

真空热试验中冷能利用的思考

在大型航天器的真空热试验中,仅从过冷器排除的液氮和低温气氮的混合物就达约5×104kg,其中蕴涵了108kJ量级的冷能.文章首先介绍了自然冷能和LNG冷能利用的研究现状,并结合航天器真空热试验的实际情况,研究了试验过程中排气造成的冷能损失,对它的再利用可行性进行了讨论,提出可以采用联合法冷能回收流程,并提出了冷能利用中的几个关键问题.认为:需要解决液氮和低温气氮供应的连续性及找到恰当的二次冷媒,是真空热试验中冷能利用的关键.     

卫星稳态热模型参数修正方法研究

为了研究卫星稳态热模型参数的修正方法,应用Thermal Desktop和Sinda／Fluint软件对一 卫星热分析模型的表面涂层热光学性质和接触换热系数进行了参数修正。研究结果显示,在 对上述两种参数进行同时修正时,尽管采用多种方法,但 所得出的修正结果精度都不高。针对上述问题,提出了对不确定参数进行分层修正的方 法,即先修正对设备温度影响较大、与设备温度相关性较高的全局关键参数,在此基础上再 进一步修正局部关键参数。分析结果表明,采用分层修正时,参数修正值与真值之间最大相 对误差为3.33％,设备温度计算值与真值之间的最大误差为2.60％,和常规修正的方法 相比,分层修正的修正精度大大提高。
     

A Series of Proton Flares in July 1985 from Observations by the Prognoz-10 Satellite and Vega-1, 2 Spacecraft

The enhancements of solar protons with an energy of more than 30 MeV, originating from flares in one active region and observed simultaneously aboard three spacecraft Vega-1, 2 and Prognoz-10 in July 1985, are analyzed and approximated in this work by using the reflection model [2]. The numerical values of several physical parameters (the parameters of distribution of solar cosmic ray (SCR) particles inside heliosphere traps, the coefficients of trap transparency) that are not observed by direct measurements are estimated.     

Autogeneration of an Instability of a Plasma Composition Equilibrium in the Plasmaspheric Transition Region during a Postsunset High-Altitude Launch of the Vertikal'-10 Rocket

A cause of the instability of equilibrium of plasma ion composition is discussed and exemplified by the data on a burst of amplitudes of small-scale plasma irregularities in the plasmaspheric transition region during an evening launch of the Vertikal'-10 rocket. This burst is accompanied by a simultaneous decrease in the average plasma density at altitudes of 700–1100 km. Specific features of the observed events are compared to postsunset incoherent scatter radar observations of the ion diffusion flux density. It is demonstrated that the instability is caused by peculiarities of the protonosphere–ionosphere interaction associated with a sharp difference between thermal conditions of the ionospheric and protonospheric air shortly after sunset. The induced nonuniformity of postsunset cooling of the protonospheric–ionospheric plasma causes density irregularities in ion diffusion fluxes and generates local bunches of heavy ions, which are usually only a minor impurity to lighter ions. As a result, conditions are created that are favorable for the nondissipative accumulation of potential energy for the mutual opposition of two or more groups of ions with different masses and for the subsequent release of this energy by a threshold excitation of impurity-driven plasma instabilities.