Proposal of application of same beam VLBI measurements in precision orbit determination of lunar orbiter and return capsule and lunar gravity field simulation in Chang’E-3 mission

High accuracy differenced phase delay can be obtained by observing multiple point frequencies of two spacecraft using the same beam Very Long Baseline Interferometry (VLBI) technology. Its contribution in lunar spacecraft precision orbit determination has been performed during the Japanese lunar exploration mission SELENE. In consideration that there will be an orbiter and a return capsule flying around the moon during the Chinese lunar exploration future mission Chang’E-3, the contributions of the same beam VLBI in spacecraft precision orbit determination and lunar gravity field solution have been investigated. Our results show that the accuracy of precision orbit determination can be improved more than one order of magnitude after including the same beam VLBI measurements. There are significant improvements in accuracy of low and medium degree coefficients of lunar gravity field model obtained from combination of two way range and Doppler and the same beam VLBI measurements than the one that only uses two way range and Doppler data, and the accuracy of precision orbit determination can reach meter level.     

Using cross correlation to estimate Doppler frequency

Doppler frequency estimation with high accuracy is very important for precise orbit determination and scientific studies in deep space exploration. In this paper, we propose a new method which implements cross correlation of only received signals from one single station to estimate Doppler frequency. The algorithm is relatively simple and it can achieve high accuracy even when signal to noise ratio (SNR) is very low. Moreover, the accuracy can be further improved by implementing frequency compensation, especially when the frequency dynamic range is high. Simulations were performed and results proved the validity of this method. X-band observations on MEX (Mars Express) and Juno were performed by a 13 m telescope of Wuhan University, which was equipped by two back ends, i.e. CDAS (Chinese Data Acquisition System) and RSR (Radio Science Receiver). Raw data recorded by the CDAS were processed by using the proposed method, and the estimated frequencies were in good agreement with the real time frequency estimation values by the RSR. What’s more, the accuracy by using cross correlation was higher than the one of RSR.     

Identification of window frequency at different places of India

Millimeterwave window frequencies between the two successive absorption maxima of 60?GHz and 120?GHz respectively, play a significant role in the context of radar and communication applications. Atmospheric parameters like temperature, water vapour and oxygen play major roles for determination of window frequencies which are latitude dependent. Radiosonde data were analyzed to identify a frequency at which minimum signal attenuation occurs in the millimeterwave band, between 60?GHz and 120?GHz, at various places in India. The data are taken from British Atmospheric Data Centre (BADC, U.K). Analysis shows that water vapour plays a dominant role for determining window frequency during the monsoon season. On the other hand, temperature dominates water vapour in shifting window frequency during the winter and summer seasons. The results obtained also show that minimum attenuation occurs at 73?GHz and maximum at 83?GHz over the chosen places in India during January to December depending on the latitudinal position. Another salient result obtained from our analyses is that water vapour is mainly responsible for lowering the window frequency from its conventionally accepted value, over certain places in India. Hence, these climatological parameters play a major role in determining window frequency over certain places of choice in India throughout the year.     

Variations of the solar acoustic high-degree mode frequencies over solar cycle 23

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present variations of the solar acoustic mode frequencies caused by the solar activity cycle. High-degree (100 < ? < 900) solar acoustic modes were analyzed using global helioseismology analysis techniques over most of solar cycle 23. We followed the methodology described in details in [Korzennik, S.G., Rabello-Soares, M.C., Schou, J. On the determination of Michelson Doppler Imager high-degree mode frequencies. ApJ 602, 481–515, 2004] to infer unbiased estimates of high-degree mode parameters ([see also Rabello-Soares, M.C., Korzennik, S.G., Schou, J. High-degree mode frequencies: changes with solar cycle. ESA SP-624, 2006]). We have removed most of the known instrumental and observational effects that affect specifically high-degree modes. We show that the high-degree changes are in good agreement with the medium-degree results, except for years when the instrument was highly defocused. We analyzed and discuss the effect of defocusing on high-degree estimation. Our results for high-degree modes confirm that the frequency shift scaled by the relative mode inertia is a function of frequency and it is independent of degree.     

失重时球腔内液体晃动特性的研究

本文研究了在失重状态下球形腔体内液体的晃动问题.首先分析了静液面的形状, 并用多项式构造了一类特征函数, 进而求解了球腔内液体晃动的力学特性, 并通过数字计算得到了一些晃动参数:如晃动频率、动力系数、阻尼比等.      

星载可展开天线研究

海洋二号卫星微波校正辐射计天线(HY-2 CRA)是高增益、低副瓣、高波束效率、低损耗的天线子系统. 其高主波束效率、低损耗等要求对天线展开精度、形面精度、测试等提出了极大挑战. 由于卫星发射整流罩结构限制, 天线采用可展开机构形式, 并用正交模耦合器实现频率复用. 天线频率分别达到18.7 GHz, 23.8 GHz 和37 GHz.为了减小天线形面变形, 保证天线波束指向, 天线反射面(包括反射面支架、馈源支架)采用碳纤维复合材料, 确保天线在宇宙空间变形尽可能小. 本文就HY-2 CRA的设计、加工、测试等进行了论述.      

Single-track effects and new directions in GCR risk assessment.

Light flashes in the eye as recorded by astronauts on missions outside the geomagnetosphere are presumably caused by single particle traversals of galactic cosmic rays traversing the retina. Although these flashes are not considered to have deleterious short- or long-term effects on vision, they are testimony that the body can detect single particle traversals. The frequencies of the flashes implicate ions in the charge range of 6 to 8 (i.e., carbon and/or oxygen ions). Other particles with higher charge and causing more ionization are present at lower frequencies. The possibility of the importance of such single-track effects in radiation carcinogenesis and other late effects suggest that a risk assessment system based on particle fluence rather than absorbed dose might be useful for assessing risk on long-term space missions. Such a system based on the concept of a risk cross section is described. Human cancer risk cross sections obtained from recently compiled A-bomb survival data are presented, and problems involving the determination of the LET-dependence of such cross sections are discussed.     

一枚Chaff火箭观测到的湍流谱

在一个异常强水平风速和极端大风切变层内, 垂直速度扰动谱表明几个相当重要特征:(1)垂直速度扰动谱有二个区域, 一个在低频区, 有谱斜率一1.65, 这个谱斜率接近于惯性子区内-5/3谱斜率, 另一个在高频区, 有谱斜率-7.11, 它接近于粘性子区内-7谱斜率.(2)垂直速度扰动谱有一个崩溃点, 约位于70m, 这与计算的内尺度有好的一致.(3)由垂直速度扰动谱计算的湍流能量耗散率是0.087m2·s-3, 这比MAP/WINE试验中得到的值约大一个数量级.      

Wave measurements in active experiments on plasma beam injection

The paper presents the results and discussion of VLF wave measurements carried out in the course of two rocket experiments on injection of a dense cesium ion beam into the ionosphere at the heights of 165–240 km. The injection was accompanied by enhancement of the broad-band noise by several orders of magnitude. The wave measurements in the frequency range of 1–11 kHz revealed two pronounced frequencies somewhat exceeding that of the lower-hybrid resonance oscillations in the background plasma and the cesium beam. The oscillations were modulated by frequencies close to the ion-cyclotron frequencies of the background plasma.     

Analysis of the code and phase between-receiver inter-system biases of the overlapping frequencies for GPS/Galileo/BDS

The overlapping-frequency signals from different GNSS constellations are interoperable and can be integrated as one constellation in multi-GNSS positioning when inter-system bias (ISB) is properly disposed. The look-up table method for ISB calibration can enhance the model strength, maximize the number of integer-estimable ambiguities, and thus is preferred. However, the characteristics and magnitudes of the receiver code ISB and phase fractional ISB (F-ISB) are not well known and the wrong values of the biases can seriously degrade the positioning results. In this contribution, we first estimate the between-receiver code ISB and phase F-ISB of hundreds of the baselines up to around 25km in the European Permanent GNSS Network (EPN) and the Multi-GNSS Experiment (MGEX) for the overlapping frequencies L1-E1 (L1), L5-E5a (L5) and E5b-B2b (L7). The data collected from 1st January 2016 to 1st January 2019. Second, the receiver-type and firmware-version combinations for the receivers of Trimble, Leica, Javad, Septentrio and NovAtel are carefully classified. Results show that the Septentrio receivers have consistent code and phase ISB values for the three overlapping frequencies i.e. only one value for each frequency and no receivers are different. The Leica, Trimble and Javad receivers have two or more ISB values for at least one of the three frequencies. A few receivers with biases to the groups are also found and listed. Third, the code ISB and phase F-ISB of the groups are adjusted by the least-squares method. The root mean square errors (RMSE) of the least square adjustment are 0.240 m, 0.250 m and 0.200 m for code of L1, L5 and L7 frequencies, respectively, and are 0.0009 m, 0.0015 m and 0.0031 m for phase of L1, L5 and L7 frequencies, respectively. Finally, the effects of code ISB errors on code positing are investigated with the zero-baseline MAT1_MATZ. The distance root mean square error (DRMS) of L1-E1 code positioning can be reduced by 48.2% with 5 GPS and Galileo satellites and the DRMS degrades quickly when the code ISB error is larger.     

Seismological studies of the sun and other stars

Observations of a large number of different oscillation frequencies in the Sun provide an opportunity for detailed testing of the theory of stellar structure and evolution. At present highly significant discrepancies remain between observed and computed frequencies, and so our models of the solar interior have to be modified. With further improvements in the observations it might become possible to make a direct empirical determination of the density structure throughout the Sun.Similar oscillations have so far not been detected in other stars, but attempts to do so are under way. Theoretical estimates indicate that amplitudes somewhat greater than for the Sun might be expected for early F stars on the main sequence, and that the amplitude increases rapidly with decreasing gravity. Observation of such oscillations would enable investigations of the structure of these stars, and would in addition provide valuable information about the excitation mechanism of the oscillations.     

基于自适应滤波的编队卫星实时相对定轨

提出基于自适应滤波的编队卫星实时相对定轨算法,利用2005-12-09—10两颗GRACE（Gravity Recovery and Climate Experiment）卫星的GPS（Global Positioning System）实测数据进行实时相对定轨试验计算,采用JPL（Jet Propulsion Laboratory）轨道对试验结果外部检核,结果表明:①自适应滤波相对定轨通过自适应因子,可以较好地平衡编队卫星的观测信息和相对动力学信息,其相对定轨结果精度优于Kalman滤波相对定轨结果;②自适应滤波相对定轨结果随着星间基线缩短而精度提高;③两颗GRACE卫星采用单频伪距和广播星历进行自适应滤波相对定轨,可以得到精度优于6cm的星间基线。     

“龙江2号”月球轨道微卫星定轨分析

微纳卫星深空探测任务中,通常所分配的测控资源有限,因此有必要对有限测控资源条件下微纳卫星的定轨精度进行分析。以微纳卫星深空探测为背景,采用"龙江2号"微卫星的轨道测量数据对其定轨精度进行了分析。"龙江2号"微卫星只有USB轨道测量数据,且环月段测控资源相对紧张,每天有两站跟踪,共约3～4 h的轨道测量数据。首先介绍了"龙江2号"微卫星飞行任务及其飞行过程中影响测定轨的因素;其次给出了定轨的动力学模型,对微卫星地月转移段的定轨精度进行了分析;最后通过分析摄动力、动量轮卸载以及数据弧段长度的影响,给出了微卫星环月阶段所使用的定轨策略,并通过重叠弧段比较的模式,给出了微卫星环月段的定轨精度。研究结论可以为后续微纳卫星深空探测任务提供有益参考。     

超空泡水下航行体的结构动力响应特性

通过综合有关文献的理论,得到作用于超空泡水下航行体尾部的冲击载荷与航行速度的关系,然后利用有限元法研究了在不同速度条件下航行体的结构动力响应和变形特性.得到了加速度响应的主频率及其取值区间、航行体内的最大应变值及其分布位置;并给出了加速度响应主频率以及航行体内最大应变值与运动速度的关系.计算结果表明:存在2个加速度响应主频率,当航行体航速相对较低时,应采用第1个主频率作为结构设计时的参考频率,速度较高时应采用第2个主频率;计算还表明:航行体内最大应变值随着航速的增加而显著增加.结果对超空泡水下航行体的结构强度分析和结构优化设计有着指导意义.     

Recent results on solar activity at submillimeter wavelengths

Since the installation of the Solar Submillimeter Telescope (SST) in 1999 in the Complejo Astronómico El Leoncito (CASLEO, Argentina), the almost unexplored solar emissions at frequencies >100 GHz started to reveal new insights about thermal and non-thermal processes in active regions. SST operates at the frequencies of 212 and 405 GHz providing the unique opportunity to distinguish and investigate emission mechanisms. We present a review of the most relevant findings obtained. An statistical study made with observations of a selected sample of active regions shows that their flux density spectra increase with frequency. Rapid brightenings (pulses) are always observed both at 212 and 405 GHz in association to solar flares lasting for some tens to hundreds of milliseconds. They are well correlated between the two frequencies and have flux spectra either flat or increasing with frequency. The flux of submillimeter wave pulses remain within the same order of magnitude for different bursts, ranging typically 100–300 s.f.u. at 212 GHz and 500–1000 s.f.u. at 405 GHz. The time evolution of the pulse occurrence rate usually reproduces the time profile of the X-rays/γ-rays emission, and the bulk emission at submillimeter waves, when the latter is observable. There are examples of good correlation between individual pulses at submillimeter waves and hard X-rays/γ-rays. Submillimeter pulses are not restricted to flare events, but appear to be a general phenomenon that occurs over active regions as well. The starting time of the rapid submillimeter wave pulses is coincident or precedes the projected launch time of the coronal mass ejections. SST observations of the November 4, 2003 large flare revealed a new and yet unknown spectral component with intensities increasing towards even higher frequencies, appearing along with, but separated from the well-known microwave component.     

RTK model and positioning performance analysis using Galileo four-frequency observations

This paper proposes a real-time kinematic (RTK) model that uses one common reference satellite for the Galileo system with four frequency observations. In the proposed model, the double-differenced (DD) pseudorange and carrier phase biases among the different frequencies are estimated as unknown parameters to recover the integer features of the DD ambiguities among the different frequencies for ambiguity resolution and precise positioning. Analysis results show that the E5a, E5b, and E5 frequencies have virtually the same performance in terms of the positioning accuracy, observation residuals, and ratio values of ambiguity resolution. However, the E1 frequency performs worse than the E5a, E5b, and E5 frequencies. The RTK results for the combination of multiple frequencies are much better than those for a single-frequency observation, the coordinates’ standard deviation is improved about 20–30%, and the ambiguity fix time is improved about 10%.     

Schumann resonance frequency and conductivity in the nighttime ionosphere of Mars: A source for lightning

We have solved the Maxwellian equations of electromagnetic waves which oscillate within the cavity formed in the lower ionosphere of Mars between 0 and 70?km. The electrical conductivity and Schumann Resonance (SR) frequencies are calculated in the lower ionosphere of Mars, in the presence of a major dust storm that occurred in Martian Year (MY) 25 at low latitude region (25°–35°S). It is found that the atmospheric conductivity reduced by one to two orders of magnitude in the presence of a dust storm. It represents a small dust layer at about 25–30?km altitudes where lightning can occur. We also found that the SR frequencies peak at?～18?km with values 19.9, 34.5 and 48.8?Hz for the modes l?=?1, 2 and 3, respectively, in the non-homogeneous medium. Our results indicate that practical or measurable values of SR are dependent on the altitudes.     

基于预设载波的随机开关频率调制方法

现有的随机开关频率调制方法未考虑载波切换频率与载波频率间的关系,因此会使逆变器输出交流电产生一定的谐波畸变.为了克服这一不足,通过分析载波切换频率对随机载波信号的影响,得出预置载波随机切换时,前后两条载波的相位不匹配是造成逆变器输出电压产生畸变的根本原因,并提出采用随机切换预设载波的方法实现开关频率的随机化.仿真及实验结果表明,随机脉冲宽度调制(PWM,Pulse Width Modulation)方法不但可以有效抑制输出电压的高次谐波而且在波形正弦度方面优于常规随机PWM方法.