Modeling Research of the 27-day Forecast of 10.7 cm Solar Radio Flux (I)

Adopting the autoregressive method for time-series modeling, we have made a study on the medium-term forecast of solar 10.7 cm radio flux (F10.7). The result of forecast experiments and the error analysis indicate that when the solar activity is at a rather low level and the 27-day periodicity of F10.7 is apparent, the autoregressive forecast method has a high accuracy and relatively ideal effectiveness, but when a large active region appears or disappears on the solar dusk, the forecast effectiveness is not ideal. This means that the autoregressive method for the time-series modeling can reflect well the 27-day periodicity of F10.7, and that it has certain applicability for building a mediumterm forecast model of F10.7. By comparing the forecast results in the period from 21th September 2005 to 7th June 2007, it is demonstrated that the accuracy of the autoregressive forecast method is equivalent to that of the forecast made by the American Air Force.     

科学计算程序库WHIGG F90 LIB及其地学应用

WHIGG F90LIB(WFL)是一款用FORTRAN90语言编写的用于地球学科的计算程序库。应用WFL可以减少重复性编程,方便地解决科学计算中数据的预处理问题,结合多种分析方法,对科学问题进行定量分析。利用WFL在地球重力场、气候涛动和数字信号处理3个方面的应用实例,以验证WFL在科学计算中应用的有效性。     

Symplectic adaptive algorithm for solving nonlinear two-point boundary value problems in Astrodynamics

In this paper, from a Hamiltonian point of view, the nonlinear optimal control problems are transformed into nonlinear two-point boundary value problems, and a symplectic adaptive algorithm based on the dual variational principle is proposed for solving the nonlinear two-point boundary value problem. The state and the costate variables within a time interval are approximated by using the Lagrange polynomial and the costate variables at two ends of the time interval are taken as independent variables. Then, based on the dual variational principle, the nonlinear two-point boundary value problems are replaced by a system of nonlinear equations which can preserve the symplectic structure of the nonlinear optimal control problem. Furthermore, the computational efficiency of the proposed symplectic algorithm is improved by using the adaptive multi-level iteration idea. The performance of the proposed algorithm is tested by the problems of Astrodynamics, such as the optimal orbital rendezvous problem and the optimal orbit transfer between halo orbits.     

Mutual eclipses of J2 Europa by J1 Io observed at Yunnan Observatory in 2009

Mutual events between natural satellites include mutual occultation and mutual eclipse. Mutual eclipse is another kind of mutual occultation as viewed from the center of the Sun instead of the Earth. Two mutual eclipses of J2 Europa by J1 Io (2009 Aug. 28 and Sept. 12) were observed at Yunnan Observatory during the PHEMU09 international campaign. We will calculate the astrometric data of these Galilean satellites by analyzing and fitting the light curves we obtained. The limb-darkening was considered during...     

从永夏陈四楼煤矿水文地质特征看岩溶陷落柱并不存在

根据区域和矿区地质资料 ,结合岩溶陷落柱的形成地质条件进行认真的分析研究 ,得出了否定的认识 ,对正确认识永夏矿区的水文地质条件 ,有效地开展矿井防治水工作 ,正确指导生产有重要意义。     

基于OpenGL的GNSS仿真软件中的飞行模拟可视化研究

探讨利用OpenGL实现GNSS仿真软件中的飞行模拟可视化方法,给出根据仿真输出的俯仰角、偏航角和速度信息实现的三维可视化结果。     

Fiber Fine Structures Superposed on the Solar Continuum Emission Near 3 GHz

On April 21, 2002, a broadband solar radio burst was observed at about 01:00 – 03:00 UT with the digital spectrometers of National Astronomical Observatories of China (NAOC). Also many fiber bursts superposed on the continuum bursts were detected in the frequency range of 2.6 – 3.8 GHz during the time interval. After data processing, some parameters of the fibers such as frequency drift rate, duration, bandwidth, and relative bandwidth were determined. The mean value of the frequency drift was in the range of 42.3 – 87.4 MHz s⁻¹ (negative). A theoretical interpretation for the fibers was presented based upon a model of the velocity of Alfvén solitons. In this model, the source of the fiber emission was considered as the ducting of the solitons within the magnetic-mirror loop. Then the magnetic field strength of the fiber source was estimated to be about 130 ≤ B₀ ≤ 270 G. Also a comparison of the magnetic field estimation was made with another model of whistler group velocity.     

Solar Stokes Spectrum Observation at Yunnan Astronomical Observatory

1　ＳｉｇｎｉｆｉｃａｎｃｅｏｆＳｏｌａｒＳｔｏｋｅｓＳｐｅｃｔｒｕｍＯｂｓｅｒｖａｔｉｏｎＤｕｅｔｏｔｈｅｄｅｖｅｌｏｐｍｅｎｔｉｎｏｂｓｅｒｖｉｎｇｉｎｓｔｒｕｍｅｎｔ,ｔｈｅｃｏｎｃｅｐｔｏｆｓｏｌａｒｓｐｅｃｔｒｕｍｓｈｏｕｌｄｂｅｇｅｎｅｒａｌｉｚｅｄ .Ｉｎｔｈｅｐａｓｔ,ｉｔｉｓｏｎｌｙａｓｓｉｇｎｅｄｔｏｔｈｅｏｒｄｉｎａｒｙｏｒｕｎｐｏｌａｒｉｚｅｄｓｐｅｃｔｒｕｍｒｅｃｏｒｄｉｎｇｔｈｅｗａｖｅｌｅｎｇｔｈ ｄｉｓｔｒｉｂｕｔｅｄｉｎｔｅｎｓｉｔｙｏｆｔｈｅｃｏｎｔｉｎｕｕｍａｎｄ…     

Excitation of Annual Polar Wobble by Global Oceans

A reasonable and quantitative result on the variation of polar wobble excited by the oceans is not available at present. Numerous researches have shown that atmospheric motion is the greatest excitation source for the seasonal variations in the polar wobble and that oceanic motion is one of the main remaining excitation sources. The excitation of variation in the annual polar wobble caused by oceans from 1992 to 2004 both globally and in latitude dependence, have been studied in depth by means of the new generation of SODA oceanic data assimilation (SODA-1.4.2 and SODA-1.4.3) and the ECCO oceanic data assimilation. The result shows that the variation in the seasonal polar wobble excited by the SODA oceans is very close to that of the residual after the action of the atmosphere and land water is deducted from the geodesic excitation function for a large part of the investigated time interval, and that there is overall agreement between the two as regards the annual amplitude and phase. In addition, in comparison with the result of early SODA-Bata 7, the new generation of SODA oceanic excitation has achieved obvious improvements. The latitude distributions of the excitations of the annual polar wobble by the SODA and ECCO oceans are consistent in the Greenwich direction, while having obvious differences in the direction of 90° E.     

AMR Simulations of Magnetohydrodynamic Problems by the CESE Method in Curvilinear Coordinates

The objective of this paper is to present new extensions of the space – time conservation element and solution element (CESE) method for simulations of magnetohydrodynamic (MHD) problems in general curvilinear coordinates by using an adaptive mesh refinement (AMR) grid system. By transforming the governing MHD equations from the physical space (x,y,z) to the computational space (ξ,η,ζ) while retaining the form of conservation, the CESE method is established for MHD in the curvilinear coordinates. Utilizing the parallel AMR package PARAMESH, we present the first implementation of applying the AMR CESE method for MHD (AMR-CESE-MHD) in both Cartesian and curvilinear coordinates. To show the validity and capabilities of the AMR-CESE-MHD code, a suite of numerical tests in two and three dimensions including ideal MHD and resistive MHD are carried out, with two of them in both Cartesian and curvilinear coordinates. Numerical tests show that our results are highly consistent with those obtained previously by other authors, and the results under both coordinate systems confirm each other very well.