2002 Vol. 22, No. 4

Display Method:
COSMIC PLASMA DYNAMO
YANG Zhiliang
2002, 22(4): 289-299. doi: 10.11728/cjss2002.04.20020401
Abstract(2246) PDF 867KB(1039)
Abstract:
A new dynamo model based on the polarization of plasma is presented in this paper. From the Maxwell equations in a moving medium, a magnetization vector can be caused with Rongon current. The steady solar magnetic field is solved from the equations. On the assumption that the meridianal flow is ignored, the distribution of magnetic field is put out. In the model, there is no additional parameter considered. The intensity of magnetic field inside the sun ranges from 1-6 T. The surface magnetic field around the pole is in the order of 1 ×10-3 T, at low latitude the calculated surface magnetic field has the order of 1×10-2 T. The maximum magnetic field is around 30° in latitude.
A CLASS OF TVD TYPE COMBINED NUMERICAL SCHEME FOR MHD EQUATIONS AND ITS APPLICATION TO MHD NUMERICAL SIMULATION
FENG Xueshang, S. T. WU, FAN Quanlin, WEI Fengsi
2002, 22(4): 300-308. doi: 10.11728/cjss2002.04.20020402
Abstract(2461) PDF 804KB(759)
Abstract:
In the past three decades, solar-terrestrial physicists have introduced many kinds of numerical schemes in computational fluid mechanics to magnetohydrody-namic system in order to simulate various phenomena of solar-terrestrial physics. Along with the recent advance of space weather, great attention has been paid to the development of quick convergence and high resolution numerical schemes of MHD system for the purpose of solar wind storm simulation or establishing numerical prediction methods of space weather. The numerical study of solar wind has undergone a transit from its early simulation problems in supersonic and superAlfvenic domain to the recent works from solar surface, interplanetary space to the interaction between solar wind and the earth's magnetosphere. In this paper, according to the characteristics of numerically modeling solar wind, a new numerical scheme of TVD type for magnetohydrodynamic equations in spherical coordinates is proposed by taking into account of the quality such as convergence, stability, resolution. This new MHD model is established by solving the fluid equations of MHD system with a modified Lax-Friedrichs scheme and the magnetic induction equations with MacCormack II scheme for the purpose of developing a combined scheme of quick convergence as well as of TVD property. To verify the validation of the scheme, the propagation of one-dimensional MHD fast and slow shock problem is discussed with the numerical results conforming to the existing results obtained by the piece-wise parabolic method. Under typical physical parameters on the solar surface, using a dipole as initiation, the steady state numerical results for the solar wind flow is reached by time-relaxation approach. This shows that this numerical model has potential application in modeling solar wind of complex magnetic field and realistic solar-interplanetary storms.
THE VARIATIONAL PROPERTIES OF GALATIC COSMIC RAYS’ INTENSITIES OBSERVED AT GUANGZHOU STATION DURING JAN.7—10,1997
LE Guiming, YE Zonghai, YU Shaohua, LI Runsen, XIAO Shaoyu, GONG Juhong
2002, 22(4): 309-313. doi: 10.11728/cjss2002.04.20020403
Abstract(2551) PDF 458KB(861)
Abstract:
The paper presents the galactic cosmic rays intensities in vertical direction, south direction and north direction respectively observed by Guangzhou muon telescope (113°18'E, 23°6'N, 21m, 16 GV) before the January 1997 CME/magnetic cloud reached the earth. The enhancement in daily variation amplitude of galactic cosmic rays' intensity in vertical direction was obvious, while the north-south anisotropy was small. Based on standard transport theory and data observed in interplanetary space, the paper concludes that the enhancement of daily variation amplitude was due to the small angle between particle velocity direction and the magnetic force line, B ×▽n drift was another reason for enhanced daily variation amplitude. The small north-south anisotropy was mainly due to the fact that the B ×▽n drift with ▽n around 16 GV was small.
THE COLLAPSES OF ELECTRIC FIELD AND FORMATION OF DENSITY CAVITONS IN THE WAKE OF MOVING BODY IN SPACE
MA Shanjun, CHENG Xiaojin
2002, 22(4): 314-320. doi: 10.11728/cjss2002.04.20020404
Abstract(2259) PDF 723KB(1004)
Abstract:
In this paper, using the method of numerical simulation, the nonsteady and nonlinear interaction between the moving body and plasma in wake region was studied in detail. Starting from the two fluid equations and the collisionless Boltz-mann equation for the ion distribution, under the condition of a static limit, for the low-frequency disturbance, the nonlinear coupling equations have been obtained. In order to get the value of Δn, the fast Fourier translation has been used and the initial condition of wavepacket form that the envelope is slowly varying function over time has been chosen in Cartesian coordinates system. For the purpose of easy to calculation the initial condition has the form of two dimensions with three field components. The results of whole wake region can be obtained by rotating the results in a two dimensional space. Of course, the choice of initial condition is not unique. In the simulation, the natural boundary conditions in x direction i.e. the field quantities are zero when x →∞, and a periodic boundary conditions in z direction have been used. The method of numerical simulation is FTCS (time forward-difference and space central-difference method). Finally the distributions of the electromagnetic soliton and density disturbance have been obtained. The basic parameters in calculation are v0 = 106 cm/s, R0 = 100cm, ne = 105/cm3, Te = 3000 K, |Ê|max2(τ= 0) = 2.42. The results show that in the wake region of the moving body, there are the formation of the electromagnetic soliton and density disturbance which can be detected due to the collapse of electric field. As the field and density increasing, the condition W < 1 is not valid, the numerical simulation has to stop. We can trace out the moving body by means of observing the structure and intensity of the density caviton and electromagnetic soliton although the body may be a stealthy vehicle.
SOME TYPICAL RECORDS OF IONOSPHERIC DOPPLER SHIFT AND THEIR SIGNIFICANCE IN THE STUDY OF IONOSPHERIC MORPHOLOGY
XIAO Zuo, LIU Kaijun, ZHANG Donghe
2002, 22(4): 321-329. doi: 10.11728/cjss2002.04.20020405
Abstract(2449) PDF 751KB(1010)
Abstract:
The measurements of HF Doppler shift in the ionosphere have long been a powerful tool in monitoring the ionospheric disturbances due to its features of continuity in time and simplicity. In this paper, based on the analysis of data from the HF Doppler station at Peking University, some typical records with special morphological characteristics of the ionosphere were shown. These records are recognized as the manifestations of sporadic-E, irregularities, acoustic-gravity waves in different scales and response to the sudden change of electron density due to solar flares. Using IRI model of the ionosphere, responses of Doppler shift to solar flare and ionospheric irregularities are simulated. The results are in agreements with observations so that the mechanics of such physical processes are explained. A special kind of records——"S" shaped trace is introduced and programs designed to simulate such records. It is shown that the "S" shaped traces are produced by horizontally propagated acoustic-gravity waves under proper conditions. Relations between shapes and AGW's parameters are also given. These results are helpful in deducing some properties of AGWs from Doppler shift records.
A NUMERICAL STUDY ON THE NONLINEAR PROPAGATION OF GRAVITY-WAVE PACKETS EXCITED BY TEMPERATURE DISTERBANCE
HUANG Chunming, ZHANG Shaodong, YI Fan
2002, 22(4): 330-338. doi: 10.11728/cjss2002.04.20020406
Abstract:
By using the spectral allocation method, a fully-nonlinear dynamical numerical model for the propagation of disturbance in the mesosphere is found. Applying the newly founded numerical model, the evolution and propagation of an initially given Gassian wave-like temperature disturbance in a compressible atmosphere were simulated. The simulation results show that the temperature disturbance excites two gravity-wave packets soon, one propagates upward and another propagates downward. The energy propagation paths of these two packets are slightly different from the ray paths predicted by the linear theory. A further numerical analysis indicates that after 3 hours' evolution of the initially given temperature disturbance, the disturbance converts entirely to wave motion, and during the conversion, only 79% disturbance energy is converted to be wave-associated energy.
SIMULATION AND RESEARCH OF SUDDEN SODIUM LAYER GUIDED BY GRAVITY WAVES
TAN Hui, WAN Weixing
2002, 22(4): 339-345. doi: 10.11728/cjss2002.04.20020407
Abstract(2287) PDF 638KB(1027)
Abstract:
The causes of formation of sudden sodium layer are not confirmed. The main viewpoints include effects of meteor ablation, dust release and atmosphere dynamics. To research into the influence of atmosphere dynamics to formation of sudden sodium layer, we simulate sudden sodium layer by applying the continuity equations for both atomic sodium and the background neutral atmosphere. As far as vertical winds are concerned, we adopt the model of a sine travelling wave which approaches the real one. The results reflect the process of sudden sodium layer's formation and development. The instant of sudden sodium layer's formation is from 5-minute to 15-minute and sudden sodium layer can last to 30-minute. The altitude of sudden sodium layer is from 90-kilometer to 100-kilometer. The width of sudden sodium layer is from 0.5 kilometers to 2 kilometers. All these are in agreement with the data obtained. And this paper also researches the tendency of sudden sodium layer's development when the parameters of gravity waves alter, including vertical wave length, wave period, velocity of the vertical winds and half-width of normal sodium layer. In all these factors, velocity of the vertical winds and half-width of normal sodium layer influence the formation of sudden sodium layer greatly. To validate the reliability of the simulation, this paper provides two cases of observation. Comparing the result of the simulation with the observation, we find that they are close to each other.
A PHYSICAL EXPLANATION OF MICROGRAVITY EFFECTS AND ITS APPLICATIONS
XI Risheng
2002, 22(4): 346-356. doi: 10.11728/cjss2002.04.20020408
Abstract(2635) PDF 992KB(1067)
Abstract:
The effects and its explanation in the point of view of molecular physics are discussed in this paper Accumulative and macroscopic features of universal gravitation are determined by its' monopole characteristics and it has almost no effects on free molecules. However, inside liquid, Especially in the liquor, obvious effects from gravity on molecular groups exist in the form of macroscopically static pressure, as gravity gets smaller and smaller, the static pressure inside liquid will disappear gradually and the secondary phenomena, such as sedimentation and floating in different phases and floating convection in same phases will also disappear. In this paper, it is pointed out that Bossinesq approximation is substantially a hypothesis of linear thermodynamics or quasi-equilibrium state, so it is possible to express the nondimensional numbers in the critical point with some physical property parameters in the initial equilibrium state. An idea is put forward that by eliminating gravity to induce a transform from non-equilibrium state to equilibrium state. The expression of Rayleigh number Ra derived in this paper shows that as gravity acceleration disappearing, the value of Ra will change from a turbulent state far from equilibrium to a ray Leigh-Benard convection cell in a stable flow, even to an equilibrium state without macroscopic flow. Finally, this theoretical explanation model is applied to estimate the microgravity level necessary for space materials processing and the interchangeability of low vacuum and microgravity in simulation of the heat conduction state in natural convections of gases in space devices.
ON-LINE MINIMUM VARIANCE TRACKING FILTER FOR THRUST ACCELERATION ESTIMATION
LI Hengnian, LI Jisheng, HUANG Yongxuan
2002, 22(4): 357-362. doi: 10.11728/cjss2002.04.20020409
Abstract(2345) PDF 476KB(1272)
Abstract:
The uncertainty of thrust acceleration will cause enormous error to model and track the accelerated motion of maneuvered target during orbit transfer process, an on-line minimum variance tracking filter for thrust acceleration estimation is developed in this paper. The propellant mass-flow rate and exhaust velocity are selected as state variants, which are estimated by employing an integral state model and Expended Kalman Filter. The variation equations for measurement vector to state variants have been established to linearize the discrete-time measurement equations. The algorithm has applied successfully to maneuver process in commanding satellite into geo-stationary orbit. The results show the algorithm developed here can estimate thrust acceleration precisely and quickly during orbit transfer process.
RESEARCH ON SECOND-ORDER NONLINEAR FILTER ALGORITHM FOR GYROLESS SATELLITE ATTITUDE DETERMINATION USING STAR SENSORS
LIN Yurong, DENG Zhenglong
2002, 22(4): 363-371. doi: 10.11728/cjss2002.04.20020410
Abstract(2138) PDF 736KB(945)
Abstract:
A second-order nonlinear filter which can estimate the attitude of a three-axis stabilized gyroless satellite from star sensor vector observations is developed using the nonlinear filter technigue based on second-order Taylor series approximation. Considering the characteristics of vector observation, the methods to dispose vector observation information in the case of one vector and multi-vector are discussed separately. In addition, the QUEST algorithm is efficiently combined into the attitude estimator as the compressing technique of multi-vector observation information in order to simplify the update implementation. Performance comparison made by computer simulation demonstrates that the second-order nonlinear attitude estimator consistently performs better than the traditional extended Kalman attitude estimator.
FREQUENCY CHARACTERISTICS AND EXPERIMENT INVESTIGATION ON A HUB-BEAM SYSTEM
YANG Hui, HONG Jiazhen, YU Zhengyue
2002, 22(4): 372-379. doi: 10.11728/cjss2002.04.20020411
Abstract(2156) PDF 853KB(1014)
Abstract:
The frequency characteristics of a hub-beam system are investigated in this paper. This rigid-flexible coupling system consists of a rigid hub, a flexible beam and a torsional spring. A finite-element model for this system is derived and all the couplings between the system rigid and flexible degrees of freedom are included. The frequency characteristics of this system are related to two ratios: the ratio J' of the beam and rigid hub inertia, the ratio K' of the torsional spring and beam bending stiffness. The effects of two ratios on the system frequencies are analyzed by numerical simulation and experimental test, respectively. The experiment is conducted in a single-axis air-bearing testbed. Both the theoretical results and experimental results show that the two ratios, i.e. K' and J', have significant influences on the system frequency characteristics. The dominant frequency of this system in steady responses can transit between the rigid and flexible frequencies in case of different values of K' and J'.