Oscillatory thermocapillary flows in simulated floating-zones with time-dependent temperature boundary conditions

This study deals with numerical simulations of the Maxus sounding rocket experiment on oscillatory Marangoni convection in liquid bridges. The problem is investigated through direct numerical solution of the non-linear, time-dependent, three-dimensional Navier-Stokes equations. In particular, a liquid bridge of silicon oil 2[cs] with a length L = 20 [mm] and a diameter D = 20 (mm) is considered. A temperature difference ΔT = 30 [K] is imposed between the supporting disks, by heating the top disk and cooling the bottom one with different rates of ramping. The results show that the oscillatory flow starts as an ‘axially running wave’, but after a transient time the instability is described by the dynamic model of a ‘standing wave’, with an azimuthal spatial distribution corresponding to m = 1 (where m is the critical wave number). After the transition, the disturbances become larger and the azimuthal velocity plays a more important role and the oscillatory field is characterized by a travelling wave. The characteristic times for the onset of the different flow regimes are computed for different rates of ramping.     

The effect of small vibrations on Marangoni convection and the free surface of a liquid bridge

The effects of small vibrations on Marangoni convection were investigated experimentally using a liquid bridge of 5 cSt silicone oil with a disk diameter of 7.0 mm, and an aspect ratio close to 0.5. Experiments were performed to determine the critical temperature difference data for no vibration case and with small vibrations applied. The experimental results have shown that the effect of small vibrations on the onset of oscillatory flow is small since the critical temperature difference data for different aspect ratios were not affected by the vibrations. To clarify the surface oscillation phenomena induced by external vibrations, a 3-D numerical simulation model was also developed using a level set algorithm to predict the surface oscillations of isothermal silicone oil bridges. By subjecting the liquid bridge to small vibrations, the surface oscillation characteristics were predicted numerically, and the numerical results compared well with the predictions of an analytical model proposed previously. Furthermore, the effect of small vibrations on the surface vibration amplitude of the liquid bridge is also discussed.     

高压液氧涡轮泵孔型/蜂窝阻尼密封的设计

针对某型火箭发动机高压液氧涡轮泵离心轮的前、后凸肩动密封,采用孔型/蜂窝阻尼密封代替原始迷宫密封方案,并对密封结构参数进行优化设计,旨在减小泵内动静间隙泄漏,提高泵容积效率.采用经水工质孔型密封泄漏量实验数据验证的,基于k-ε 湍流模型和三维定常Reynolds-Averaged Navier-Stokes(RANS)...     

variations of the atmospheric chemical composition in the earths’s polar regions after solar proton flare on July 14, 2000 (photochemical simulation)

Altitude—temporal cross-sections q(z, t) of atmospheric ionization rates by solar protons above the polar regions were calculated using the GOES-10 satellite data on solar proton fluxes for the period of solar proton flare (SPF) on July 14, 2000. The values of q(z, t) were used further in calculations of variations of the atmospheric chemical composition during the flare in the northern and southern polar regions (70°N and 70°S) by two different 1D photochemical models of the atmosphere (neutral and charged components). The calculation results have shown considerable variation of the ozone content after SPF: a decrease of [O3] was about 80% at altitudes of 65–75 km above northern and 25% in the layer of 55–65 km above the southern polar region. Such decrease of the ozone content is a result of reactions with [NO] and [OH] whose concentrations have grown substantially during SPF. According to calculations, the increase of electron concentration during SPF has reached 3–4 orders of magnitude at altitudes of 50–80 km. A comparison of the calculation results with the observational data on [NO], [NO2], and [O3] from the UARS and HALOE satellites for 70°N have shown a good qualitative correspondence, however, for variations of nitric oxides there are quantitative discrepancies.     

Thermocapillary simulation of single bubble dynamics in zero gravity

The lack of significant buoyancy effects in zero gravity conditions poses an issue with fluid transfer in a stagnant liquid. In this paper bubble movement in a stagnant liquid is analysed and presented numerically using a computational fluid dynamics (CFD) approach. The governing continuum conservation equations for two phase flow are solved using the commercial software package Ansys-Fluent v.13 and the Volume of Fluid (VOF) method is used to track the liquid/gas interface in 2D and 3D domains. The simulation results are in reasonable agreement with the earlier experimental observations, the VOF algorithm is found to be a valuable tool for studying the phenomena of gas–liquid interaction. The flow is driven via Marangoni influence induced by the temperature difference which in turn drives the bubble from the cold to the hot region. A range of thermal Reynolds (Re_T) and Marangoni numbers (Ma_T) are selected for the numerical simulations, specifically Re_T=13–658 and Ma_T=214–10,721 respectively. The results indicate that the inherent velocity of bubbles decreases with an increase of the Marangoni number, a result that is line with the results of previous space experiments (Kang et al., 2008) [1]. An expression for predicting the scaled velocity of bubble has been derived based on the data obtained in the present numerical study. Some three-dimensional simulations are also performed to compare and examine the results with two-dimensional simulations.     

Space-time structure of ion-cyclotron waves in the topside ionosphere as observed onboard the <Emphasis Type="Italic">ST</Emphasis>-<Emphasis Type="Italic">5</Emphasis> satellites

Results of observations of ion-cyclotron (IC) waves onboard the ST-5 satellites in the topside ionosphere (heights from a few hundred up to thousands of km) are presented. In this project, three identical micro-satellites were located during three months in 2006 in almost identical orbits with distances between them from first thousands to hundreds of km. All ion-cyclotron wave packets detected by two-three probes were observed at crossing one and the same latitude, which manifests their narrow localization in latitude with a characteristic scale from the first tens to 100 km. In no event IC waves were recorded with comparable amplitudes by all three satellites. At the same time, in the case of ST-5 flight near the ground-based induction magnetometer, a long emission in the same frequency range on the ground corresponded to a burst of IC waves in the topside ionosphere. This can indicate to the fact that an IC instability develops not continuously, but in the pulsing regime with a characteristic time of up to ∼10 min. A change in the rotation direction when a satellite crosses the wave structure is a characteristic feature of the polarization structure of registered transverse waves. The detected effects are discussed from the point of view of the existing models of generation and waveguide propagation of IC waves.     

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.     

Large-scale irregularities of the winter polar topside ionosphere according to data from <Emphasis Type="Italic">Swarm</Emphasis> satellites

An analysis of the electron density measurements (Ne) along the flyby trajectories over the high-latitude region of the Northern Hemisphere under winter conditions in 2014 and 2016 has shown that the main large-scale structure observed by Swarm satellites is the tongue of ionization (TOI). At the maximum of the solar cycle (F_10.7 = 160), the average value of Ne in the TOI region at an altitude of 500 km was 8 × 10⁴ cm^–3. Two years later, at F_10.7 = 100, Ne ~ 5 × 10⁴ cm^–3 and Ne ~2.5 × 10⁴ cm^–3 were observed at altitudes of 470 and 530 km, respectively. During the dominance of the azimuthal component of the interplanetary magnetic field, the TOI has been observed mainly on the dawn or dusk side depending on the sign of B _y . Simultaneous observations of the convective plasma drift velocity in the polar cap show the transpolar flow drift to the dawn (B_y < 0) or dusk side (B _y < 0). Observations and numerical simulation of the Ne distribution have confirmed the significant role of the electric field of the magnetospheric convection in the generation of large-scale irregularities in the polar ionosphere.     

Surface temperature distribution along a thin liquid layer due to thermocapillary convection

The surface temperature distributions due to thermocapillary convections in a thin liquid layer with heat fluxes imposed on the free surface are investigated. The nondimensional analysis predicts that, when convection is important, the characteristic length scale in the flow direction L, and the characteristic temperature difference ΔT₀, can be represented by and , respectively, where L_R and ΔT_T are the reference scales used in the conduction-dominant situations with A denoting the aspect ratio and Ma the Marangoni number. Having had L and ΔT₀ defined, the global surface-temperature gradient ( ), the global thermocapillary driving-force, and other interesting features can then be readily determined. Finally, numerical calculations involving a Gaussian heat flux distribution are presented to justify these two relations.     

Low Altitude Cusp–Cleft Region Signatures of Strong Geomagnetic Disturbances

The polar cusp being a region of the free access of the solar wind into the inner magnetosphere is also the site of turbulent plasma flow. The cusp area at low altitudes acts like a focus of a variety type of instability and disturbances from different regions of the Earth. Daily f ₀ F2 frequencies are discussed regarding the cusp position. The high time resolution wave measurements together with electron and ion energetic spectra measurements registered on the board the Freja satellite and Magion-3 and the electron density at the peak of f ₀ F2 layers collected from ground-based ionosonde measurements were used to study the response of ionospheric plasma within the cusp–cleft region to the strong geomagnetic storm. In this paper we present the response of the ionospheric plasma to the disturbed conditions seen in the topside wave measurements and in the ionospheric characteristics maps obtained from the ground-based VI network. The need of the cusp feature model for radio communication purposes is advocated.     

Numerical solution of nonlinear heat problem with moving boundary

Two phase gas–liquid flow in pipes is widely spread in space applications: bubble flows appear in cryogenic components transport through fuel/oxidant supply lines. Another important application is based on the fact that in liquid flows with small bubbles a close contact between the two phases occurs resulting in high rates of transfer between them. The compactness of a system makes it ideally suited to serve as a space-based two-phase bio-reactor which forms an important unit in environmental control and life support system deployed onboard. A numerical method was developed for solving a nonlinear problem of thermal interaction between a spherical gas bubble and surrounding liquid. The system of equations for describing this interaction was formulated. It includes ordinary and nonlinear partial differential equations. The problem was solved using finite-difference technique by dividing the system into spherical layers inside the bubble and employing the new variable which “freezes” the moving boundary of the bubble. A numerical solution is obtained for the problem of radial bubble motion induced by a sudden pressure change in the liquid—a situation which corresponds to the behavior of bubbles beyond a shock wave front when the latter enters a bubble curtain.     

Numerical investigation of the impact of asymmetric fuel injection on shock train characteristics

Numerical simulations are carried out to investigate the impact of asymmetric fuel injection on shock train characteristics using the commercial-code FLUENT. The asymmetry of fuel injection is examined by changing the fuel flow rates of the upper and lower wall fuel injectors. The numerical approach solves the two-dimensional Reynolds-averaged Navier–Stokes (RANS) equations, supplemented with a k-ω model of turbulence. As a result, different ways of fuel injections will always lead to shock train transitions, with the variations of shock train structure, strength and leading edge position. For symmetric fuel injection, the flowfield of the isolator is quite asymmetric with the boundary layer of the upper wall side developing much stronger than that of the lower wall, which is due to the heterogeneity of the incoming flow. Regarding to asymmetric fuel injection with more of lower wall side, though the pressures in the combustor are nearly the same, the first shock of the shock train converts between ‘Distinct symmetric X type shock’ and ‘Obscure and weaker asymmetric shock’ and the shock train leading edge moves upstream with the increase of the asymmetry level. With regard to asymmetric fuel injection with more of upper wall side, ‘incomplete asymmetric X type shock’ occurs and the shock train structures keep nearly the same with low level of fuel injection asymmetry. Unexpected results like unstart will happen when increasing the level of fuel injection asymmetry. And the isolator will come back to normal state by decreasing the differential of upper and lower wall sides fuel injections.     

The Fields of Magnetospheric Current Systems on the Earth's Surface in an Interval of the Magnetic Storm Studied under the International Program on Solar–Terrestrial Physics

Using the magnetic storm in January 1997 as an example, we examined the possibilities to employ the magnetospheric field T96 [1, 2] and the dynamic paraboloid model PM of the magnetosphere [3] for modeling the D _st variation. We have revealed the necessity to refine the results of normalizing the free parameters of the model T96 according to the solar wind parameters. The contributions to the D _st variation of magnetic fields of basic large-scale magnetospheric current systems (the field DCF on the magnetopause, the field DR of the ring current, and the field DT in the magnetotail) are estimated for different phases of the storm from model calculations. Possible causes of a discrepancy between the results of modeling D _st using the T96 and PM models are discussed. Special emphasis is made on the ratios of contributions into the D _st variation of the fields of the magnetotail and the ring current in the main phase of magnetic storms and on the contributions to D _st of the fields of various current systems at the recovery phase.     

不同高度通道内热薄燃料火焰传播特性的研究

为探究通道高度对逆风火焰传播特性的影响,采用实验、数值模拟、理论分析方法对微重力条件下不同高度通道内热薄燃料火焰传播特性进行了研究。通道高度分别为10mm, 14mm和60mm的短时落塔微重力实验结果表明,相同逆风气流速度时,随通道高度增加,火焰传播速度逐渐增大,且10mm, 14mm和16mm高度通道下的火焰长度比约为通道半高度比的二次方关系。数值模拟与理论分析表明,微重力条件下,随通道高度增加,火焰传播速度、火焰长度先增大再减小,呈非单调的变化趋势。燃烧热释放速率随高度变化呈非单调变化趋势,同火焰传播速度变化规律基本一致。     

Development Study on ATREX Engine

This is the status report of the development study on ATREX engine (Air Turbo Ramjet) that is now under way in the Institute of Space and Astronautical Science (ISAS) cooperation with the Ishikawajima Harima Heavy Industries (IHI), the Kawasaki Heavy Industries (KHI), the Mitsubishi Heavy Industries (MHI). ATREX engine will be applied for the propulsion system of fly-back booster of TSTO space plane. ATREX is the combined cycle (a fan-boosted ramjet) engine providing the effective thrust from sea level static to flight Mach number 6. ATREX is worked on the expander cycle with precooling the incoming air as shown in Fig. 1. ATREX employs the tip turbine configuration which allows the compactness and the light weight of turbo machinery and the variable geometry airintake and plugnozzle which allow the wide range operation conditions.From 1990 to 1992, “ ATREX-500“ has been tested at the sea level static conditions. ATREX-500 is the 1/4-scale model of which fan inlet diameter is 300 mm and overall length 2,200 mm. From 1992 have been performed the wind tunnel tests on the primary components of ATREX, the axisymmetric variable geometry airintakes, the precoolers and the variable geometry plug nozzles. In parallel to the windtunnel tests, the ram combusters have been tested simulating the hypersonic flight conditions and the application studies on advanced carbon-carbon composite for the tip-turbine and fan assembly has been proceeded.In 1994 initiated the flight test plan in which ATREX will be verified in the practical flight conditions by using an unmanned flying test bench.In 1995 will be tested ATREX-500 installing the precooler under the sea level static conditions to examine the engine performance and the icing on the precooler.The present paper addresses the high loading ram combuster experiment using the mixer with skewed lobes to generate swirl flow and the analytical studies and the designs on the precooler and the precooled ATREX engine and the flight test plan.     

Statistical Approach to the Evaluation of the Mean Correlation Length and the Propagation Velocity of Middle-Scale Plasma Structures in the Earth's Foreshock

Based on simultaneous measurements of ion fluxes made onboard the closely separated satellites Interball-1and Magion-4, the propagation velocity of middle-scale plasma structures in the Earth's foreshock relative to the solar wind flow is estimated. The derived value of this velocity allows these structures to be identified as a fast magnetosonic wave propagating upstream of the solar wind inflowing the Earth's bow shock. An evaluation is also made of the correlation length of these disturbances in the plane perpendicular to the Sun–Earth line. This length is approximately equal to 2R _E.     

旋转充液系统全飞行过程非线性动力学仿真

本文考虑一类刚认耦合的旋转充液系统的飞行动力学问题，在依据力学变分原理给出基本方程，导出全飞行过程运动状态方程的基础上，建立了相应的仿真模型。可能通过对Navier-Stokes方程求解和液体流动的惯性波振动特性，给出液体的反作用力矩。经过仿真实验，分析角运动的时间序列，给出了比较详尽的飞行不稳定机制分析结果。     

Numerical investigation on the shock wave transition in a three-dimensional scramjet isolator

The scramjet isolator, which is used to prevent the hypersonic inlet from disturbances that arise from the pressure rise in the scramjet combustor due to the intense turbulent combustion, is one of the most critical components in hypersonic airbreathing propulsion systems. Any engineering error that is possible in the design and manufacturing procedure of the experimental model, and the intense heat release in the scramjet combustor, may cause the performance of the isolator to decrease, leading to its lack of capability in supporting the back pressure. The coupled implicit Reynolds Averaged Navier–Stokes (RANS) equations and the two-equation standard k?ε turbulent model have been employed to numerically simulate the flow fields in a three-dimensional scramjet isolator. The effects of the divergent angle and the back pressure on the shock wave transition and the location of the leading edge of the shock wave train have been estimated and discussed. The obtained results show that the present numerical results are in very good agreement with the available experimental shadow-pictures, and the numerical method is more suitable for capturing the shock wave train and predicting the location of the leading edge of the shock wave train in the scramjet isolator than the present two-dimensional numerical methods. This is due to the small width-to-height ratio of the isolator and the intense three-dimensional flow structures. On increasing the divergent angle of the scramjet isolator, the static pressure along the central symmetrical line of the isolator decreases sharply. This is due to the strong expansion wave generated at the entrance of the isolator, and when the divergent angle of the isolator is sufficiently large, namely 1.5°, a zone of negative pressure is formed just ahead of the leading edge of the shock wave train. At the same time, the shock wave train varies from being oblique to being normal, and then back to oblique. With an increase in the prescribed back pressure at the exit of the scramjet isolator, the leading edge of the shock wave train moves forward towards the entrance of the isolator, and when the back pressure is sufficiently large, unstart conditions in the hypersonic inlet can take place if the shock train reaches the inlet.     

超音速复速级涡轮的气动设计改进

采用数值模拟方法对某液体火箭发动机超音速复速级涡轮进行了流场分析。根据分析结果对两列动叶的叶型进行了改进设计：第一列动叶栅的改进采用自由旋流法,通过等通道的叶栅流道设计,减弱了激波对附面层的干扰,有效抑制了流道内的流动分离;第二列动叶栅的改进采用参数化叶片造型法,型线用具有局部修改能力和保凸性较好的Bezier曲线表示,通过减小入口攻角降低了分离损失。数值模拟结果表明,改进后的复速级涡轮内部流动特性改善显著,分离损失明显减小,效率提高了5％以上。