低纬子午环配备科学CCD的必要性和可行性

从基本天体测量学科的发展趋势和低纬子午环的课题目标出发，叙述了该仪器配备科学ＣＣＤ的必要性和对仪器性能的相应要求，并根据仪器的各种误差的测定方法和测定精度的分析，论述了该仪器配备科学ＣＣＤ后对天体位置作高精度绝对测定的可行性，文中还提出了对配备科学ＣＣＤ的总体要求。     

Perturbation of the radial and non-radial oscillations of a star by a magnetic field

A generalization of the perturbation method is applied to the problem of the radial and the non-radial oscillations of a gaseous star which is distorted by a magnetic field. An expression is derived for the perturbation of the oscillation frequencies due to the presence of a weak magnetic field when the equilibrium configuration is a spheroid. The particular application to the homogeneous model with a purely poloidal field inside, due to a current distribution proportional to the distance from the axis of symmetry, and a dipole type field outside is considered in detail. The main result is that the magnetic field has a large and almost stabilizing effect on unstableg-modes, particularly on higher order modes. With the considered magnetic field the surface layers appear to have a large weight.     

Monte Carlo Modeling of the Thermal Conductivity of Porous Cometary Ice

The basic structure of comet nuclei is an aggregation of grains, with a size distribution that extends over several orders of magnitude and a similar distribution of pores. Although attempts have been made to assess the effect of porosity on the thermal conductivity, the effect of pore size distribution has been ignored. Modeling a porous structure with a wide size distribution would require a very fine 3-D grid, so as to accommodate the smallest and largest voids. In order to circumvent this difficulty, we adopt a hierarchical procedure. Thus we assume a random and fractal porous structure and use a 3-D Monte Carlo model. The basic configuration is a cube made of unit cells of two types, (ice) filled and void, randomly distributed. Their relative number corresponds to a prescribed porosity. We solve the heat transport equation for this cube until a steady state is obtained, and from this solution the effective thermal conductivity is derived. The calculations are repeated for a range of porosities and temperatures, since the ice conductivity is temperature dependent. The basic cube serves as a unit filled cell in a larger cube, and in this way the hierarchical structure of the medium is built up.We find that the thermal conductivity is lowered by several orders of magnitude at high porosities. The correction factor, obtained as a fit to the results of our calculations, is expressed as a smooth function of the porosity, which tends to zero as the porosity approaches the percolation threshold of the solid. If only the porosity of the medium is known, this correction is not uniquely determined, but rather a range of values is possible. Only if the size distribution of the pores is known does the correction become uniquely determined.     

A new regularization of the restricted three-body problem and an application

A new regularizing transformation for the three-dimensional restricted three-body problem is constructed. It is explicitly derived and is equivalent to a simple rational map. Geometrically it is equivalent to a rotation of the 3-sphere. Unlike the KS map it is dimension preserving and is valid inn dimensions. This regularizing map is applied to the restricted problem in order to prove the existence of a family of periodic orbits which continue from a family of collision orbits.     

Propagation of auroral infrasonic waves studied by ray-tracings

On the basis of a ray tracing method the propagation and the attenuation of an auroral infrasonic wave are studied. Relations between the direct and reflected waves recorded at the Syowa Station, Antarctica, are clarified with regard to; (1) the delay time, (2) the intensity ratio, and (3) trace velocities. The time required for a wave to travel from the source to the ground is calculated as a function of a source altitude. The retardation time of the wave arrival behind the zenith crossing of the source current is deduced. A method is proposed for estimating the altitude of a source current from the retardation and a trace velocity of the wave. It is concluded that the existence of a supersonic equatorward motion of an electrojet which continues for a certain distance is necessary for the observation of auroral infrasonic waves. This distance must exceed at least 60 km equatorwards from the zenith to enable the direct wave to be observed and with total length of 930 km to enable the reflected wave to be observed. From these conditions it is also concluded that the infrasonic wave is not seen in mid latitudes and the reflected wave is a rare phenomenon.     

On the Influence of a Large-scale Magnetic Field on Turbulent Motions in an Electrically Conducting Medium

This paper deals with turbulent motions in a homogeneous incompressible electrically conducting medium in the presence of a magnetic field which is on average homogeneous and stationary. Using a model in which the turbulence is produced by a stochastic body force, and supposing a weak interaction between motion and magnetic field, a method is developed for calculating the pair correlation tensor of the velocity field from that occuring in a zero magnetic field. As an example, the pair-correlation tensor for a homogeneous stationary turbulence, which is isotropic and mirror-symmetric for zero magnetic field, is determined. With obvious assumptions on the correlation for zero field, two results are obtained. Firstly, the turbulent velocity is reduced by the magnetic field, the component parallel to the field, however, less than those perpendicular to it. Secondly, the correlation length parallel to the field turns out to be greater than the one perpendicular to it, indicating a tendency towards two-dimensional motion. Finally, the possibility of special situations is briefly discussed in which the turbulent velocity is enhanced by the magnetic field, and the anisotropies of the velocity components and the correlation lengths are opposite to those above.     

点在球面多边形内外的判定方法及应用

精确判断卫星进出南大西洋辐射异常区(SAA)对于保护卫星上重要仪器、延长卫星寿命和制定观测计划至关重要,这就是所谓点定位问题,简单说来就是判定球面上一点P是否位于指定的球面区域R内的问题。传统的研究方法是把球面简化为平面,研究点与平面多边形的关系。但SAA区域为球面上的区域,用平面简化得出的点定位误差较大。提出判断点在球面区域内外的算法,处理了可能出现的奇异情况,实现了点在球面区域内的准确定位。     

Self-consistent electromagnetic field in a nearly co-rotating magnetosphere

An example of the self-consistent solution which belongs to the non-trivial solution, obtained in a previous paper (Kaburaki, 1985), is found in a nearly co-rotating inner magnetosphere. Though the stellar wind is neglected there compared with the co-rotatinal velocity, drift motion around the magnetic axis, which is a manifestation of inertial effects, is determined self-consistently with the electromagnetic field. In this process, explicit expressions for the energy integral in the rotating frame and for the density distribution are also obtained. These expressions contain a fundamental length, which is to be evaluated according to physical conditions of a magnetosphere and determines the asymptotic-kinetic energy of a plasma particle at infinity. The electric current associated with the drift motion is too small to alter the original magnetic field, but the electric field is modified by the inertial effects even in the inner magnetosphere. The integrated Ohm's law is used to describe a force balance in the rotating frame, in the limits of weak and strong magnetic field.     

Resonant behaviour of MHD waves on magnetic flux tubes

Resonant absorption of MHD waves on a nonuniform flux tube is investigated as a driven problem for a 1D cylindrical equilibrium. The variation of the fractional absorption is studied as a function of the frequency and its relation to the eigenvalue problem of the MHD radiating eigenmodes of the nonuniform flux tube is established. The optimal frequencies producing maximal fractional absorption are determined and the condition for total absorption is obtained. This condition defines an impedance matching and is fulfilled for an equilibrium that is fine tuned with respect to the incoming wave. The variation of the spatial wave solutions with respect to the frequency is explained as due to the variation of the real and imaginary parts of the dispersion relation of the MHD radiating eigenmodes with respect to the real driving frequency.     

单位矢量法的数学模型(MMUVM)及其简化形式(PUVM1)的迭代法的收敛性

利用人造地球卫星观测资料测定有摄初轨的单位矢量法(PUVM1),已得到了非常广泛地实际应用.为了对单位矢量法作进一步地完善和改进,首先在考虑测量误差模型的基础上,建立单位矢量法所对应的数学模型MMUVM.它本质上就是一个非线性最优化问题.针对MMUVM,先分别使用多圈仿真数据和实测数据,形成了与之相对应的目标函数,再利用求解最优化问题的一种三对角二次插值模型的直接搜索方法,分别对其进行了数值处理.计算结果表明,所建立的优化模型MMUVM正确合理,所采用的直接搜索方法实用有效.其次,进一步指明了PUVM1和MMUVM之间关系,即:从本质上讲,PUVM1就是MMUVM的一种简化形式.从数学原理上,清楚地解释了利用PUVM1的准法化方程,只能使用单圈短弧段数据进行初始轨道确定,而不能使用长弧段多圈资料进行轨道确定或轨道改进的根本原因.最后,对PUVM1的迭代算法的收敛性问题进行了初步的理论分析,并给出了相应的数值验证实例,指出了PUVM1的迭代格式是条件收敛的,即:只有在满足一定条件后,才能收敛.这也就意味着:有的时候,尽管准法化方程是合情合理的,但是,此时该迭代法却是发散的,无法迭代求出所要的解.     

Mass motions in a heated flare filament

Energy transport in a hot flare plasma is examined with particular reference to the influence of fluid motion. On the basis of dimensional considerations the dynamical timescale of the flare plasma is shown to be comparable to the timescale for energy loss by conduction and radiation. It is argued that mass motion is likely to have a profound influence on the evolution of the flare.The detailed response of a flare filament to a localized injection of energy is then analyzed. Radiative, conductive and all dynamical terms are included in the energy equation. Apart from greatly enhancing the rate of propagation of the thermal disturbance through space, mass motion is found to be significant in transferring energy through the moving fluid.Finally the predicted thermal structure is discussed and it is concluded that the presence of mass motions in the flare may be inferred from the form of the soft X-ray differential emission measure.     

An analytic solution to the classical two-body problem with drag

An analytic solution to the two-body problem with a specific drag model is obtained. The model treats drag as a force proportional to the vector velocity and inversely proportional to the square of the distance to the center of attraction. The solution is expressed in terms of known functions and is of a simple and compact form. The time-of-flight is expressed as a quadrature in the ‘true anomaly’.     

The role of the Yoshizawa effect in the Archontis dynamo

The generation of mean magnetic fields is studied for a simple non-helical flow where a net cross-helicity of either sign can emerge. This flow, which is also known as the Archontis flow, is a generalization of the Arnold–Beltrami–Childress flow, but with the cosine terms omitted. The presence of cross-helicity leads to a mean-field dynamo effect that is known as the Yoshizawa effect. Direct numerical simulations of such flows demonstrate the presence of magnetic fields on scales larger than the scale of the flow. Contrary to earlier expectations, the Yoshizawa effect is found to be proportional to the mean magnetic field and can therefore lead to its exponential instead of just linear amplification for magnetic Reynolds numbers that exceed a certain critical value. Unlike α effect dynamos, it is found that the Yoshizawa effect is not notably constrained by the presence of a conservation law. It is argued that this is due to the presence of a forcing term in the momentum equation, which leads to a non-zero correlation with the magnetic field. Finally, the application to energy convergence in solar wind turbulence is discussed.     

Quantum and thermal phenomena in the matter-gravity system

The matter-gravity system is examined in a path integral approach for the case of conformal matter coupled to a Friedman-Robertson-Walker space time. In particular the case of gravitational potentials of interest in cosmology for which the universe tunnels from a small radius is examined. It is observed that in the presence of such gravitational horizons the universe evolves in a complex time and it is shown how a classical time and temperature emerge. Correspondingly, one will have compensating quantum and thermal fluctuations for the matter and gravity system and it is noted that the unstable mode of gravity corresponding to the universe tunneling into existence will be compensated by an analogous mode for matter corresponding to its creation. This last point is examined in a simple De Sitter model with conformal matter and a relation is found between the cosmological constant, the number of matter fields and the self coupling of matter responsable for its instability.     

Numerical simulation of slow shocks in the solar wind

The evolutionary state of slow forward shock waves is examined with the use of two MHD numerical codes. Our study is intended to be exploratory rather than a detailed parametric one. The first code is one-dimensional (with three components of velocity and magnetic field) which is used to follow a slow shock that propagates into a positive gradient of density versus distance. It is found that the slow shock evolves into an extraneous (intermediate) shock wave. The second code has a spherical, one-dimensional, planar geometry (with two velocity and magnetic field components) which is used to follow a spiral interplanetary magnetic field. It is found that a slow shock type perturbation can generate a forward slow shock; a fast forward shock is generated in the front of the slow shock; a contact discontinuity is formed behind the slow shock, and a compound nonlinear MHD wave is formed behind the contact discontinuity with a fast reverse shock formed further behind. Thus, we demonstrate that the evolution of a slow shock into (solely) a fast shock, as suggested by Whang (1987), is much more complicated.     

Substorm related changes in the geomagnetic tail: the growth phase

Changes in the configuration of the geomagnetic tail are known to play a fundamental role in magnetospheric substorms. Observations with the UCLA magnetometer on the eccentric orbiter OGO-5 indicate that the most pronounced changes in the midnight meridian occur in the cusp between 8 and 11 R_e. In order to organize the observations it is necessary to separate effects on the tail due to the solar wind magnetic field and effects due to substorms. Provided there are no changes in the solar wind there are two distinct phases of a substorm in the near tail: a growth phase and an expansion phase. During each phase the observations depend on the location of the satellite relative to the plasma sheet boundaries. Far behind the Earth is the pure tail region which consists of the lobe and the plasma sheet. In the lobe the field magnitude is characteristically enhanced relative to the dipole. Closer to the Earth is a region of transition. The field magnitude is close to that of the dipole but its orientation is distorted forming a cusp-like field line. Near the Earth is a region of depressed field. Here the field magnitude is much less than that of the dipole, but its orientation is similar. The growth phase of a substorm appears to be the direct consequence of the onset of a southward solar wind magnetic field. In the pure tail region the lobe field begins to increase in magnitude and the plasma sheet thins. The transition region moves earthward and the field lines become more tail-like as the field magnitude increases. In the inner region of depressed field, the field magnitude decreases rapidly. The onset of the expansion phase appears to be a process internal to the magnetosphere and independent of the solar wind. In the depressed field region there is a rapid, turbulent increase in field magnitude. In the transition region there is a sudden decrease in the field magnitude and a return to dipolar orientation. In the tail region the plasma sheet expands rapidly with the field becoming quite dipolar, decreasing slowly in the lobe of the tail.     

Electrostatic double layers and a plasma evacuation process

An evacuation process due to the growth of current driven instabilties in a plasma is discussed. The process, which leads to localized extreme density reductions, is related to the formation of electrostatic double layers. The initial linear phase is treated using the superposition of unstable plasma waves. In the long wavelength, non-dispersive limit, a density dip — which is initially present as a small disturbance — grows rapidly and remains localized in the plasma. The process works for a variety of plasma conditions provided a certain current density is exceeded. For a particular choice of plasma parameters the non-linear development is followed, by solving the coupled Vlasov-Poisson equations by finite difference methods. The evacuation process is found to work even more effectively in the non-linear phase and leads to an extreme density reduction within the dip. It is suggested that the growth of such structures produces weak points within the plasma that can lead to the formation of double layers.     

原子钟噪声的蒙特卡络模拟方法

讨论了原子钟噪声的蒙特卡络模拟方法．具体过程是：先产生模拟所需的正态分布随机数，从实际钟的阿仑方差估算出噪声电平，然后用１组递点函数产生高质量的模拟钟时差的随机序列．在对１台商品铯钟进行具体模拟后，给出了真实数据与模拟值的阿仑方差之问的相符程度，最后对实际问题进行了分析．     

Analytical results for solar sail optimal missions with modulated radial thrust

The aim of this paper is to analyze the optimal trajectories of a spacecraft subjected to a modulated radial thrust, whose magnitude is inversely proportional to the square of the distance from the primary body. This case is representative of a Sun-facing solar sail with a passive attitude control system. In this study the sailcraft is assumed to perform a finite number of reorientation maneuvers to set the propelling acceleration to zero and generate suitable coasting arcs along the trajectory. Accordingly, the resulting generalized orbit is a sequence of either propelled or ballistic conic arcs, whose main characteristics (in terms of semimajor axis, eccentricity, and perihelion radius) can be calculated in closed form. As a result, the sailcraft optimal performance can be studied using an analytical approach. In particular, some compact relationships are drawn and discussed that allow one to find the optimal sailcraft characteristics required to reach a prescribed final orbit.     

Radially accelerated optimal feedback orbits in central gravity field with linear drag

The solution of a feedback optimal control problem arising in orbital mechanics is addressed in this paper. The dynamics is that of a massless body moving in a central gravitational force field subject also to a drag and a radial modulated force. The drag is linearly proportional to the velocity and inversely proportional to the square of the distance from the center of attraction. The problem is tackled by exploiting the properties of a suitably devised linearizing map that transforms the nonlinear dynamics into an inhomogeneous linear system of differential equations supplemented by a quadratic objective function. The generating function method is then applied to this new system, and the solution is back transformed in the old variables. The proposed technique, in contrast to the classical optimal control problem, allows us to derive analytic closed-loop solutions without solving any two-point boundary value problem. Applications are discussed.