A Class of TVD Type Combined Numerical Scheme for MHD Equations With a Survey About Numerical Methods in Solar Wind Simulations

It has been believed that three-dimensional, numerical, magnetohydrodynamic (MHD) modelling must play a crucial role in a seamless forecasting system. This system refers to space weather originating on the sun; propagation of disturbances through the solar wind and interplanetary magnetic field (IMF), and thence, transmission into the magnetosphere, ionosphere, and thermosphere. This role comes as no surprise to numerical modelers that participate in the numerical modelling of atmospheric environments as well as the meteorological conditions at Earth. Space scientists have paid great attention to operational numerical space weather prediction models. To this purpose practical progress has been made in the past years. Here first is reviewed the progress of the numerical methods in solar wind modelling. Then, based on our discussion, a new numerical scheme of total variation diminishing (TVD) type for magnetohydrodynamic equations in spherical coordinates is proposed by taking into account convergence, stability and 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 (PPM). Finally, some conclusions are made. This revised version was published online in August 2006 with corrections to the Cover Date.     

Discrete energy release in microwave emission in the preflare stage

Examples of solar microwave precursors and estimates of their energetic characteristics are discussed based on spectral investigations of the preflare activity.     

Design-to-Cost Applied Within an Integrated System Model

In an environment of declining financial budgets for space projects, new approaches - such as Design-To-Cost - are being implemented to improve today's satellite design processes. Using an example of a current mission (the power subsystem of the Solar Probe spacecraft) under study at NASA's Jet Propulsion Laboratory, the main part of the paper discusses an Integrated System Model (structured into a performance model, a cost model, and an effectiveness model) that is part of a model-based design process used to perform cost-effectiveness trades. A simulation tool is used during the first step to size the components of the power subsystem, and then simulate its performance during operation. The determined dimensions are transferred into an EXCEL^TM-spreadsheet and linked to the components' costs. With a cost accounting tool that combines cost estimating relationships with the Work Breakdown Structure of the power subsystem, the life-cycle cost of each alternative design concept is computed. To determine the cost-risk of the different design alternatives for each component, cost probability distributions are introduced. By performing Monte-Carlo simulations, cost sensitivities are revealed. In the next step of the trade study process, the effectiveness of the alternatives is analyzed. Having determined cost and effectiveness, estimates can be made for where the different alternatives lie in the design space. The last part of the paper identifies the main drivers for the spacecraft's performance and cost. Finally it is shown how the mission design benefited from the Integrated System Model and from the application of Design-To-Cost.     

Properties of Interplanetary Coronal Mass Ejections

Interplanetary coronal mass ejections (ICMEs) originating from closed field regions on the Sun are the most energetic phenomenon in the heliosphere. They cause intense geomagnetic storms and drive fast mode shocks that accelerate charged particles. ICMEs are the interplanetary manifestations of CMEs typically remote-sensed by coronagraphs. This paper summarizes the observational properties of ICMEs with reference to the ordinary solar wind and the progenitor CMEs.     

面向深空探测的电磁帆推进技术研究进展

电磁帆是一类利用太阳风驱动的无工质或少工质新型推进技术,在需要推力器长时间工作的深空探测任务中具有非常广阔的应用前景。首先介绍了电磁帆的基本概念和主要分类,其次分别介绍了电帆、纯磁帆和磁等离子体帆三类推进方式的系统组成和工作原理,重点介绍了各自的研究现状,并梳理出了相关关键技术,同时介绍了等离子体磁罩技术。最后对电磁帆推进技术研究进展进行了总结,为国内开展此方面研究提供了研究方向和发展思路。     

Modeling and analysis of mover gaps in tubular moving-magnet linear oscillating motors

A tubular moving-magnet linear oscillating motor (TMMLOM) has merits of high efficiency and excellent dynamic capability. To enhance the thrust performance, quasi-Halbach permanent magnet (PM) arrays are arranged on its mover in the application of a linear electro-hydrostatic actuator in more electric aircraft. The arrays are assembled by several individual segments, which lead to gaps between them inevitably. To investigate the effects of the gaps on the radial magnetic flux density and the machine thrust in this paper, an analytical model is built considering both axial and radial gaps. The model is validated by finite element simulations and experimental results. Distributions of the magnetic flux are described in condition of different sizes of radial and axial gaps. Besides, the output force is also discussed in normal and end windings. Finally, the model has demonstrated that both kinds of gaps have a negative effect on the thrust, and the linear motor is more sensitive to radial ones.     

椭圆振动切削加工的同步控制研究

利用椭圆振动切削加工的回转体表面微织构是在机床主轴旋转、平台进给和装置椭圆振动的联合运动下得到的,机床主轴转动与装置位移输出如果不能保持完全的实时同步,会造成微织构加工的偏差。为了解决这一问题,提出了一种新的基于转角的同步控制方法,该方法在主轴的对位基准处实现装置正弦信号的复位,以修正转速漂移造成的加工位置误差,并通过判断主轴角度信号实时输出对应幅值电压来实时同步信号输出和机床旋转运动。基于LabVIEW FPGA搭建同步控制模块,经试验验证,利用该方法加工的矩形阵列微织构轴向排列与中轴线平行,与周向排列的夹角和仿真结果的偏差不超过0.1°,织构深度随转角的不同发生预设的变化,因此,同步控制方法具有实时同步控制椭圆振动切削装置的能力。     

碳纳米管阵列推力器放电特性的数值模拟研究

基于直接蒙特卡洛(DSMC)算法、浸入式有限元算法(IFE)、粒子云及蒙特卡洛碰撞(PIC-MCC)算法以及Ammosov-Delone-Krainov(ADK)隧穿电离模型建立了碳纳米管阵列推力器工质气体和CNT场电离过程的三维仿真模型。利用文献中的实验数据对本文仿真模型的有效性进行了验证,结果表明本文仿真模型与实验结果基本一致。针对碳纳米管阵列推力器的放电特性进行了研究,并重点分析了离子与原子碰撞过程对碳纳米管阵列推力器放电特性的影响。结果表明CNT可以极大地增强局部电场,稳态时场增强因子约为1308;与原子的碰撞使得离子在CNT附近更加集中,降低了局部电场,进而降低了其场增强能力。     

Direct position determination using single moving rotating linear array: Noncoherent and coherent processing

To improve the resolution and accuracy of Direct Position Determination (DPD), this paper investigates the problem of positioning multiple emitters directly with a single moving Rotating Linear Array (RLA). Firstly, the geometry of the RLA is formulated and analysed. According to its geometry, the intercepted noncoherent signals in multiple interception intervals are modeled. Correspondingly, the MUltiple SIgnal Classification (MUSIC) based noncoherent DPD approach is proposed. Secondly, the synchronous coherent pulse signals are individually considered and formulated. And the coherent DPD approach which aims for localizing this special type of signal is presented by stacking all array responses at different interception intervals. Besides, we also derive the constrained Cramér-Rao Lower Bound (CRLB) expression for both noncoherent and coherent DPD with RLA under the constraint that the altitudes of the emitters are known. At last, computer simulations are included to examine the performance of the proposed approach. The results demonstrate that the localization accuracy and resolution of DPD with single moving linear array can be significantly improved by the array rotation. In addition, coherent DPD with RLA further improves the resolution and increases the maximum emitter number that can be localized compared with the noncoherent DPD with RLA.     

Direction finding for two-dimensional incoherently distributed sources with Hadamard shift invariance in non-uniform orthogonal arrays

This paper proposes a novel algorithm for Two-Dimensional (2D) central Direction-of-Arrival (DOA) estimation of incoherently distributed sources. In particular, an orthogonal array structure consisting of two Non-uniform Linear Arrays (NLAs) is considered. Based on first-order Taylor series approximation, the Generalized Array Manifold (GAM) model can first be established to separate the central DOAs from the original array manifold. Then, the Hadamard rotational invariance relationships inside the GAMs of two NLAs are identified. With the aid of such relationships, the central elevation and azimuth DOAs can be estimated through a search-free polynomial rooting method. Additionally, a simple parameter pairing of the estimated 2D angular parameters is also accomplished via the Hadamard rotational invariance relationship inside the GAM of the whole array. A secondary but important result is a derivation of closed-form expressions of the Cramer-Rao lower bound. The simulation results show that the proposed algorithm can achieve a remarkably higher precision at less complexity increment compared with the existing low-complexity methods, which benefits from the larger array aperture of the NLAs. Moreover, it requires no priori information about the angular distributed function.