离子推力器羽流沉积对卫星热控影响研究

离子推力器的羽流是等离子体,等离子体的组成是带电粒子,这与传统的化学推进系统的羽流成分有很大不同,带电粒子有在卫星表面吸附的倾向,会形成羽流沉积污染。这种羽流沉积会改变卫星表面的吸收率和发射率,从而影响卫星的热控性能。为了预测离子推力器的羽流对卫星的热控性能的影响,建立了离子推力器羽流模型。所建模型采用了工程化离子推力器的在卫星上的布局位置和离子推力器的工作参数,模拟了离子推力器的正离子与中和电子束在工程化中分置的实际情况,使模型更为符合实际。通过数值模拟得到了离子、电子、中性粒子的空间分布,电场分布,得到了钼粒子在卫星表面的分布及沉积厚度,比较了模型计算的离子分布与实验获得的离子分布情况,说明了模型分析的正确性,给出了卫星表面热性能的变化及局部区域温升的最大包络可达二十多度的结果。     

Soliton-induced electric currents in plasma

This is a theoretical study of the nonequilibrium motion of charged particles in an electric field of solitons. We show that the self-consistent electric field of ion-acoustic and electron-acoustic solitons is characterized by one-way transfer of charged particles at a distance of several Debye radii. The dependence of relevant local currents on the amplitude of solitons is determined. We consider the practically important case of a moving cascade consisting of many solitons and show that the induced currents have a significant constant component. The kinetic energy acquired by charged particles in the soliton field is calculated. The temporal resolution required for the recording of soliton-induced currents is estimated. The calculations presented here can be used to interpret the results of experiments conducted to study solitons in the space plasma.     

卫星磁力矩器磁场效应仿真分析

磁力矩器是卫星上磁性较大的部件之一,局部的强磁场是影响卫星周围带电粒子分布的主要原因。为了定量分析磁力矩器与等离子体相互作用,文章基于磁场与运动粒子的耦合模型,利用有限元数值分析软件仿真了典型轨道航天器上磁力矩器工作时周围带电粒子入射轨迹变化及电荷(粒子数)分布。结果表明,受磁场作用,带电粒子在接近磁力矩器本体时会发生偏转,且偏转幅度随电子能量降低而增大;由磁场引起的这种偏转效应导致带电粒子不均匀入射到卫星内。     

磁力矩器对星外等离子体分布影响分析

文章通过计算某型号磁力矩器工作时产生的磁场,分析了等离子体在磁力矩器磁场中的运动状态,计算了电子在磁力矩器周围的密度分布。结果显示:在磁力矩器工作时,周围电子密度将增大,最大可为正常密度的10倍。密度增大的区域尺寸在轴向上略大于磁力矩器的长度,越靠近磁力矩器的两端电子密度越小,形成一个“凹陷”;径向上主要在半径小于4 m的范围内。在此区域应避免摆放对带电粒子敏感的部组件。     

Attitude dynamics of a pendulum-shaped charged satellite

In this paper, we investigate the possibility of the use of the Lorentz force, which acts on charged satellite when it is moving through the magnetic field, as a means of satellite attitude control. We first derive the equations of attitude motion of charged satellite and then investigate the stability of the motion. Finally we propose an attitude control method using the Lorentz force. Our method requires moderate charge level for future Lorentz-augmented satellite.     

The reflection model for propagation of solar cosmic rays

The principles of constructing the model are considered in this paper, together with the issues of its components and physical algorithm. The main emphasis is placed on the cause-and-effect relationships between the structure and dynamics of the solar and interplanetary media and the motion of solar cosmic ray particles. The specific features of the initial phase of motion of the particles after their acceleration in a solar flare are considered, as well as the subsequent phase of further particle propagation in the solar corona, and the phase of motion in the heliosphere. The processes of transport of charged particles in flare plasma ejections are considered in detail.     

Plasma irregularities in unstable outer ionosphere according to the Intercosmos-Bulgaria-1300 satellite data

The results of the satellite low-latitude and mid-latitude measurements of the disturbed plasma concentration, electron temperature, and quasi-stable electric field at heights of ～900 km after sunset are discussed. It is shown that the sharp fronts of changes in the electron temperature and plasma density observed in the experiment onboard the Intercosmos-Bulgaria-1300 satellite in the low-latitude (and equatorial) outer ionosphere can be related to damping of the oscillations of plasma electrons at local decreases of the plasma density (plasma “pits”) and formation of the vortex plasma structures at density and temperature gradients, which promotes conservation of ionosphere irregularities and makes the fronts of concentration variations steeper. Nonmonotonic variations in the plasma conductivity for the ionosphere currents in unstable plasma can be a cause of observed nonmonotonic disturbances of the vertical component of the “constant” electric field.     

Uncontrolled motion of the <Emphasis Type="Italic">Foton M-2</Emphasis> satellite and quasistatic microaccelerations on its board

The results of determining the uncontrolled rotational motion of the Foton M-2 satellite (in orbit from May 31 to June 16, 2005) are presented. The determination was made using the data of onboard measurements of the Earth’s magnetic field strength. Segments 270 min long (three orbits) were selected from these data covering the first two thirds of the flight. On each such segment the data were processed jointly by the least squares method using integration of the equations of attitude motion of the satellite. In processing, the initial conditions of motion and parameters of the used mathematical model were estimated. The thus obtained results gave a complete overview of the satellite motion. This motion, having started with a small angular velocity, gradually accelerated, and in two days became close to the regular Euler precession of an axisymmetric solid body. On June 09, 2005 (the last day of measurements) the angular velocity of the satellite relative to its lengthwise axis was about 1.1 deg/s, while the projection of the angular velocity onto a plane perpendicular to this axis had an absolute value of about 0.11 deg/s. Deviations of the lengthwise axis from a normal to the orbit plane did not exceed 60°. Based on the results of determination of the rotational motion of the satellite, calculations of quasi-static microaccelerations on its board are performed.     

A study of evolution of spin-up mode of a satellite in the plane of its orbit

We investigate the mode of spinning up a low-orbit satellite in the plane of its orbit. In this mode the satellite rotates around its principal central axis of the minimum moment of inertia which executes small oscillations with respect to the normal to the orbit plane; the angular velocity of the rotation around this axis several times exceeds the mean orbital motion. Gravitational and restoring aerodynamic moments are taken into account in the satellite’s equations of motion. A small parameter characterizing deviation of the satellite from a dynamically symmetric shape is introduced into the equations. A two-dimensional integral surface of the equations of motion, describing quasi-steady-state rotations of the satellite close to cylindrical precession of the corresponding symmetrical satellite in a gravitational field, has been studied by the method of small parameter and numerically. Such quasi-steady-state rotations are suggested to be considered as unperturbed motions of the satellite in the spin-up mode. Investigation of the integral surface is reduced to numerical solution of a periodic boundary value problem of a certain auxiliary system of differential equations and to calculation of quasi-steady-state rotations by the two-cycle method. A possibility is demonstrated to construct quasi-steady rotations by way of minimization of a special quadratic functional.     

Test-Particle Trajectories in a “Sheared” Stationary Field: Newton–Lorentz and First Order Drift Numerical Simulations1

We study a magnetic field distribution that is nonuniform and sheared like in tangential discontinuities. This distribution is an input parameter for the numerical integration of the equations of motion of the test-particle and of its guiding center. Two different electric field distributions are alternatively tested. In the first case, the electric field is uniform and constant like the electric field prescribed in the large-scale, steady-state reconnection models. The numerical solution shows that in this case the test-particle is trapped within the discontinuity into a region where (i) B goes to zero or (ii) the magnetic vector becomes exactly parallel to the electric field. In the second case, we consider an electric field, which is nonuniform. Its components are computed such that the zero order (or electric) drift is everywhere perpendicular to the discontinuity surface and its value is conserved throughout the simulation. In this case the numerically integrated trajectory of the test-particle penetrates the discontinuity for any angle of shear of B. Direct comparison between exact (Newton–Lorentz) and approximated (first order drift) numerical solutions shows that the mathematical singularities of the latter do not correspond to any physical singularity of the exact equation of motion of the particle.     

The use of refined model of aerodynamic moment in the problem of reconstruction of the <Emphasis Type="Italic">Foton</Emphasis> satellite rotational motion

A new mathematical model of the uncontrolled rotational motion of the Foton satellite is presented. The model is based on the Euler dynamic equations of rigid body motion and takes into account the action upon the satellite of four external mechanical moments: gravitational, restoring aerodynamic, moment with constant components in the satellite-fixed coordinate system, and moment arising due to interaction of the Earth’s magnetic field with the satellite’s proper magnetic moment. To calculate the aerodynamic moment a special geometrical model of the outer satellite shell is used. Detailed form of the formulas giving above-mentioned moments in the equations of satellite motion is agreed with the form of the considered motion. Model testing is performed by determining with its help the rotational motion of the Foton M-2 satellite (it was in orbit from May 31, 2005 to June 16, 2005) using the data of the onboard measurements of the Earth’s magnetic field strength. The use of the new model has led to a relatively small improvement in the accuracy of the motion determination, but allowed us to obtain physically real estimates of some parameters.     

Multispectral remote sensing of aerosols

Aerosols are common atmospheric constituents that occur both naturally (clouds, sea spray, dust, smoke, and volcanic emissions) and artificially (smog, smoke, certain hazes, detonation products, and industrial emissions). Some, like the great dust bowl storms in the U.S. in the 1930s, are a combination of natural and manmade agents. Most aerosols are difficult to model because they are composed of small, non-spherical particles whose optical constants and particle sizes are poorly known. Spectroscopic observations of aerosols in the thermal infrared atmospheric window between 8 and 13.5 μm offer the opportunity to detect aerosols both day and night down to very low column densities. Such observations can also identify the gross chemical composition of the particles and, in some cases, the actual sizes and shapes. In this paper, we discuss thermal i.r. observations of three types of aerosols: satellite measurements of volcanic dust, ground-based observations of airborne desert soil and both ground- and space-based measurements of cirrus clouds.     

20cm口径离子推力器性能建模与仿真

为了分析国内首台通过在轨飞行测试的20cm离子推力器栅极系统束流离子运行特性和推力器性能,针对该推力器栅极系统建立了束流引出二维数值仿真计算模型,利用PIC/MCC数值仿真计算方法,模拟束流引出过程中带电粒子在电场作用下的加速、聚焦与引出、带电粒子与中性原子之间的相互作用、电场和等离子体流场之间的相互耦合等过程。数值计算结果显示,屏栅截获的离子电流约为1.71×10 -4 A,加速栅截获的电流和CEX离子电流分别为0 A和9.11×10 -7 A,因此,加速栅电流的主要来源是冲击到其表面的CEX离子,证明了加速栅电流的主要来源是冲击到其表面的CEX离子,计算的加速栅截获电流与束流电流之比约为0.114%。试验测得推力器运行4000h期间,电子反流极限电压始终为75~90〖KG*9〗V,其变化幅度很小,这意味着中和器发射的电子在栅极系统中的反流不会导致其发生失效。理论计算结果与试验测试值相比,误差约为1.08%。〖JP〗     

基于粒子云方法研究地磁场对运动电子束的时空分布影响

空间电荷效应是空间带电粒子分布直接产生的相关效应,对于理解空间辐射环境有着重要意义。为研究地球磁场对空间电荷效应的影响,文章通过数值求解Poisson方程和相对论粒子运动方程,建立粒子云（particle-in-cell, PIC）自洽模型,并构建地球偶极子磁场模型作为内磁层背景磁场（空间分布L值为3～7）。在模型基础上,给出内磁层典型电子束结构（能量0.1～100 keV）的两种初始分布（Uniform和Gaussian分布）,模拟电子束在地磁场影响下的动力学过程。结果表明:电子束结构受地磁场影响出现纵向的整体偏转,并由于自场效应产生了纵向的拉伸和横向的发散;电子束整体的偏转角主要受地磁场强弱以及电子束运动与磁场夹角影响,而拉伸和发散效应则主要受电子束能量等参数所影响。     

Features of wave excitation of the electrostatic ion cyclotron type in the auroral ionosphere

Broadband electrostatic noise in the auroral ionosphere can be identified as a version of waves of an electrostatic ion cyclotron type, excited by plasma instability resulting from an inhomogeneous distribution of wave energy density. Broadband waves are generated due to both electric field inhomogeneities and plasma density inhomogeneities. The effect of the form of the distribution of electric field and plasma density inhomogeneities on the excitation of instabilities is studied. Also there is shown the role of the characteristic scale of inhomogeneities in the generation of electrostatic ion cyclotron waves due to the development of instability of this kind. The study of these issues, which are important for understanding the processes in the auroral region, is the subject of this paper. The work presents also a comparison of numerical results obtained using both satellite data and model approximations.     

Secular evolution of rotary motion of a charged satellite in a decaying orbit

An electrostatically charged Earth satellite whose orbit is decaying due to the Earths oblateness is considered. Secular perturbations of the orbit are taken into account: they are caused by the second zonal harmonic of the geopotential. These perturbations represent deviations of the longitude of the ascending node and perigee argument, the orbit form being invariable and the orbit inclination to the equatorial plane being constant. The attitude rotary motion of the satellite under the action of perturbing moments of the gravitational and Lorentz forces is studied. The magnetic field of the Earth is taken in a quadrupole approximation. The evolution of the satellites rotary motion is investigated on the basis of new differential equations in s-parameters specially constructed for this purpose. Using the method of averaging, basic regularities of the secular evolution of rotary motion of a screened satellite are revealed. It is found that the rotary motion of a charged satellite essentially depends on the quadrupole component of the geomagnetic potential.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 2, 2005, pp. 111–125.Original Russian Text Copyright © 2005 by Tikhonov.     

近赤道卫星绕飞编队的运动学分析及轨道设计研究

首先对一般运动学方法在进行近赤道卫星编队相对运动分析和轨道设计时存在的问题进行了分析 ;其次介绍了一种基于零倾角轨道变换的运动学新方法 ,用其对近赤道卫星编队中参考卫星轨道倾角对环绕卫星轨道根数的影响进行了研究 ;最后对零倾角卫星编队相对运动方程的线性化误差进行了理论分析 ,并利用数值仿真对两种运动学方法的相对运动方程线性化误差进行了比较分析。数值仿真的结果表明 :由基于零倾角轨道变换运动学方法得到的相对运动方程的线性化误差不随参考卫星轨道倾角改变 ,而由一般运动学得到的相对运动方程的线性化误差随着卫星编队接近赤道而呈非线性增大。     

Excitation and propagation of an electromagnetic pulse in magnetized plasma in the low-hybrid frequency range

We analyze specific features of slow quasi-static waves excited in the low-hybrid frequency band by a source whose dimensions are much smaller than electromagnetic wavelength. Main attention in this paper is given to an analysis of emission of harmonic signals by a dipole source in pulsed mode. First, the issue of influence of electromagnetic, thermal, and collision corrections in the dispersion equation on the field structure has been studied. Second, the structure of electric and magnetic fields near the resonant cone has been analyzed: in particular, effects of group delay and anomalous spreading of signals are considered. Thus constructed theory explains results of the OEDIPUS-C experiment, in which, for example, already at distances of order of ten wavelengths a signal delay of approximately 10^?4 s was observed. Finally, we have studied some aspects of the inverse problem of electrodynamics: the role played in field formation by smoothness of charge distribution over antenna is demonstrated, and a class of smooth distributions of charge on antenna is found that form a preset field structure.     

On the planar inertial motion of a system of three coupled bodies

Numerous papers are devoted to studying the motion of a system (coupling) of two bodies in the Earth’s satellite orbit ([1–4] and others). The problem on the planar inertial motion of three bodies, coupled by a non-extensible weightless string having the form of an unfastened chain, is considered in the paper. Such a configuration can be represented, for example, by a system of two coupled spacecraft rotating around their common center of mass (in order to simulate the gravity force) in long-term space missions, when the third body (the lift) is located on a connecting cable. The bodies are considered to be the material points (particles).     

One method of gravitational attitude control of a rotating satellite

The mode of spinning up a low-orbit satellite in the plane of its orbit is studied. In this mode, the satellite rotates around its longitudinal axis (principal central axis of the minimum moment of inertia), which executes small oscillations with respect to the normal to the orbit plane; the angular velocity of the rotation around the longitudinal axis is several tenths of a degree per second. Gravitational and restoring aerodynamic moments were taken into account in the equations of satellite’s motion, as well as a dissipative moment from eddy currents induced in the shell of the satellite by the Earth’s magnetic field. A small parameter characterizing deviation of the satellite from a dynamically symmetric shape and nongravitational external moments are introduced into the equations. A two-dimensional integral surface of the equations of motion, describing quasistationary rotations of the satellite close to cylindrical precession of the corresponding symmetrical satellite in a gravitational field, has been studied by the method of small parameter and numerically. We propose to consider such quasistationary rotations as unperturbed motions of the satellite in the spin-up mode.