多种观测数据驱动的三维行星际太阳风MHD模拟

三维磁流体力学（MHD）数值模拟是行星际太阳风研究的重要手段.本文发展了一种由多种观测数据驱动的三维行星际太阳风MHD数值模型.模型的计算区域为0.1AU到1AU附近，使用Lax-Friedrich差分格式在六片网格系统中进行数值求解.边界条件中磁场使用GONG台站观测的光球磁图外推获得，密度通过LASCO观测的白光偏振亮度反演得到，速度根据以上两种观测数据并利用一种基于人工神经网络技术（ANN）的方法得到，温度通过自洽方法根据磁场和密度导出.利用该模型模拟了第2062卡灵顿周（CR2062）时期的行星际太阳风，模拟结果显示出丰富的观测特征，并与OMNI以及Ulysses的实际观测值符合得较好.该模型可用于提供接近真实的行星际太阳风，有助于提高空间天气预报的精度.      

高频Alfven波供能的日冕及太阳风模型

空间探测的最新结果表明，太阳大气中广泛存在小尺度喷发事件，这些事件被认为是由新浮现磁通量导致的小尺度磁重联的表现．小尺度磁重联同时会产生高频Alfven波．本文给出了建立在上述思想基础上的描述日冕加热及太阳风加速的时变二元流体模型的一组数值解，显示解的时间演化、收敛特性及守恒特性．在对由过渡区发出的高频Alfven波频谱的适当假设下，数值解能够解释日冕及太阳风等多方面的观测结果，包括日冕加热，极冕洞密度径向变化，太阳风近日冕底部加速，及在03AU观测到的高速流特征．     

WSA太阳风经验模型及其应用

Wang-Sheely-Arge （WSA）模型是对准稳态太阳风的经验和物理相结合的描述,其利用观测的日面磁图作为输入,可以提前3到4天预测L1点处的太阳风速度和行星际磁场极性.WSA模型是在WS模型的基础上经过若干改进形成的实时预报模式,之后又借鉴Distance from the Coronal Hole Boundary （DCHB）模型的参数,进一步改进了太阳风速度关系式,形成了目前常见的形式.WSA经验模型由日冕磁场模型、太阳风速度关系式和一维运动学模型三部分组成.在实际应用过程中,基本步骤包括观测磁图预处理、日冕三维磁场反演、计算日冕磁场参数、计算太阳风的速度分布和将太阳风映射到1AU等环节.在发展过程中,WSA模型经历了一些细节上的调整变化,例如观测磁图数据的来源、日冕磁场模型的类型、经验速度关系中自由系数的取值等.许多研究对如何改善模型的预报效果进行了探索.      

两冕流间CME事件数值模拟的改进

通过偶极子场和六极子场适当叠加，改进猜解磁场，使猜解磁场在太阳南北极符号相反，然后采用理想磁流体力学方程组(MHD)，由猜解磁场与太阳风流动相互作用计算出稳态自洽解，得到定性上与观测比较接近的具有两个冕流的背景结构．在两个冕流间采用具有同心圆磁场位形的触发模型触发CME事件，研究CME的日冕传播特征．模拟结果表明，CME被约束在两冕流间传播，CME闭磁场位形和磁云横截面磁场位形相似，可以解释1AU处观测磁云的部分特征；在CME附近，存在压力和Lorentz力起主要作用的区域，这可以为分析1AU处CME事件的观测数据提供帮助．     

1998年5月2日日冕亮度观测图的数值研究

采用数值手段模拟了1998年5月2日日冕亮度观测图，计算模式改进为球坐标下系下特殊的二维理想磁流体（MHD）模型，即把（r,φ）坐标建立在SOHO观测日冕亮度的子午面上，消除了子午面极区的几何奇异性，根据SOHD日冕观测布置磁极子得到初始磁场位形，内边界条件采用自治的投影特征边界条件，计算迭代出稳态的多 磁场结构，得得到了与观测基本一致的亮度图，计算结果表明在太阳表面附近磁场位形对太阳凤等离子参数分布起控制作用。     

无振荡、无自由参数格式在三维太阳风流动数值模拟中的应用

从快速、有效地进行空间天气数值预报的需要出发,针对1998年5月份行星际太阳风暴事件,采用全三维流体力学(HD)模型进行数值试验,考察了无振荡、无自由参数(NND)格式在三维太阳风流动数值模拟中的应用.计算中内边界的密度分布由观测K日冕数据来确定,由此得到的源表面密度分布具有和源表面电流片相似的结构.数值试验表明: 虽然该格式不需要加入人工粘性,但是具有很好的计算稳定性.在扰动计算时激波间断耗散小,在三维任何方向上间断所占的网格数比较少,没有数值色散现象.      

基于CPU/GPU的日冕偏振亮度并行计算模型

三维磁流体力学(MHD)数值模拟是用来研究日冕和太阳风最常用的方法之一, 其中将计算得到的日冕电子数密度转化为日冕偏振亮度(Polarization Brightness, PB)是与观测对比的重要方法. 由于待转换电子数据网格密度、PB数据网格密度和计算模型的复杂度, 使得日冕偏振亮度的计算比较耗时, 利用单CPU计算无法达到近实时转换日冕偏振亮度的要求, 从而影响了数值模拟的验证效率. 本文在CPU/GPU环境下, 利用CUDA编程技术, 提出了一个日冕偏振亮度并行计算模型. 实验结果表明, 该模型比CPU上的串行模型计算速度提高了31.86倍, 达到了近实时模拟与观测数据比对的计算要求.      

一个基于观测的太阳风背景模型

建立由太阳光球磁场和日冕偏振亮度等观测约束的单流体太阳风模型,包括日冕和太阳风的等离子体密度、速度和磁场,温度还有待于以后处理.这里采用高山观测台(HAO)MKⅢ的日冕偏振亮度(pB)在1.36Rs上的观测概图,根据Guhathakurta在1996年发展的日冕电子密度反演模型确定日冕的电子密度分布.同时采用Wilcox太阳观测台(WSO)的光球磁场视向分量的观测概图作为底部边界,根据Zhao等在1994年发展的水平电流-电流片(HCCS)模型得到全球磁场.Phillips在1995年及McComas在2003年分别用Ulysses第一次和第二次跨极飞行的观测发现,归一化到1 AU的太阳风动量流密度除了在10°～30°的纬度范围内略低以外几乎不变.根据这一结论,结合已经得到的密度数据,就可以得到日冕和太阳风的速度.将上面的模型应用于1918卡林顿自转周稳态太阳风的研究,结果与太阳活动极小期的观测基本相符,但是与观测相比较低速高密度区偏大,因此密度模型还有待改进.      

二维磁流体动力学平衡日冕(Ⅱ)

基于解析和数值相结合的方法,进一步讨论了非均匀引力场中日冕的二维磁流体动力学平衡。对临界点进行了比较仔细的处理。得到了包含闭场区、中性片和开场区的大尺度日冕磁场位形,闭场区和中性片构成冕流结构。在高纬和低纬地区几个太阳半径之外,等离子体径向流动速度超过了局地声速和局地Alfvén速度。在1AU处,太阳风速度可达到400kms_-1以上      

一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点，考虑到算法的质量（收敛速度、稳定性、精度等），结合磁流体数值计算的特性，对三维球坐标磁流体动力学（MHD）方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacComack格式，发展了一类快捷的具有TVD特性的组合数值新方法，作为格式的检验，在一维情况下，将其与PPM格式进行了比较，对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果，然后将其诮到定态太阳风的数值模拟上，在不同等离子体β情形下，可得到理想的太阳风定态结构，为今后将此数值模式应用到具有复杂磁场位型或三维直实太阳风暴的数值模拟研究奠定了基础。     

Learning in an immersive digital theater

The Houston Museum of Natural Science, in collaboration with Rice University has an outreach program taking portable digital theaters to schools and community sites for over five years and has conducted research on student learning in this immersive environment. By using an external independent evaluator, the effectiveness of NASA-funded Education and Public Outreach (EPO) projects can be assessed. This paper documents interactive techniques and learning strategies in full-dome digital theaters. The presentation is divided into Evaluation Strategies and Results and Interactivity Strategies and Results. All learners from grades 3–12 showed statistically significant short-term increase in knowledge of basic Earth science concepts after a single 22-min show. Improvements were more significant on items that were taught using more than one modality of instruction: hearing, seeing, discussion, and immersion. Thus immersive theater can be an effective as well as engaging teaching method for Earth and Space science concepts, particularly those that are intrinsically three-dimensional and thus most effectively taught in an immersive environment. The portable system allows taking the educational experience to rural and tribal sites where the underserved students could not afford the time or expense to travel to museums.     

Stochastic particle acceleration by the forced interaction of relativistic current sheets

We address the problem of interacting relativistic current sheets in self-consistent kinetic plasma simulations within the framework of the Particle-In-Cell model. The interaction is enforced in head-on collisions of up to 10 current sheets at relativistic bulk speeds. The simulations are motivated by the general problem of Poynting flux dissipation in ‘striped wind’ configurations presumably governing the relativistic outflows pervasive in pulsar winds and gamma-ray bursts. We identify the generation of non-thermal particles and formation of a stable power-law shape in the particle energy distributions f(γ) dγ ∝ γ^−s dγ. In 1D, a spectral index s ∼ 2 is observed and attributed to a stochastic Fermi-type acceleration mechanism. In 2D, the generic index of s ∼ 3–4 is retained as in previous simulations of individual current sheets. Whereas in 2D the high energy cut-off is constrained by the limited dissipation of magnetic energy, in 1D the process converts the bulk motion of current sheets towards directed particle momentum of an exclusive class of non-thermal particles.     

Scaling analysis and model estimation of solar corona index

A monthly average solar green coronal index time series for the period from January 1939 to December 2008 collected from NOAA (The National Oceanic and Atmospheric Administration) has been analysed in this paper in perspective of scaling analysis and modelling. Smoothed and de-noising have been done using suitable mother wavelet as a pre-requisite. The Finite Variance Scaling Method (FVSM), Higuchi method, rescaled range (R/S) and a generalized method have been applied to calculate the scaling exponents and fractal dimensions of the time series. Autocorrelation function (ACF) is used to find autoregressive (AR) process and Partial autocorrelation function (PACF) has been used to get the order of AR model. Finally a best fit model has been proposed using Yule-Walker Method with supporting results of goodness of fit and wavelet spectrum. The results reveal an anti-persistent, Short Range Dependent (SRD), self-similar property with signatures of non-causality, non-stationarity and nonlinearity in the data series. The model shows the best fit to the data under observation.     

Development of a gas chromatography compatible Sample Processing System (SPS) for the in-situ analysis of refractory organic matter in martian soil: preliminary results

In the frame of the 2009 Mars Science Laboratory (MSL) mission a new sample preparation system (SPS) compatible with gas chromatography–mass spectrometry (GC–MS) has been developed for the in situ analysis of complex organic molecules in the Martian soil. The goal is to detect, if they exist, some of the key compounds that play an important role in life on Earth including carboxylic acids, amino acids and nucleobases.     

基于LuGre模型的电液加载系统摩擦补偿

为了提高电液加载系统控制精度,针对摩擦问题提出了基于LuGre摩擦模型的前馈补偿方法.建立了用于摩擦仿真分析和补偿器设计的电液加载系统数学模型;通过实验获取并分析了相关的摩擦数据;基于实验数据进行LuGre模型参数辨识,把LuGre模型和辨识结果引入电液加载系统数学模型,并进行仿真结果与实际摩擦数据的对照,证明了LuGre摩擦模型的准确性.设计前馈补偿器,进行了实验对比,实验结果表明前馈补偿器可将摩擦产生的控制误差有效地降至未补偿时的30%左右.      

CoRoT and Kepler results: Solar-like oscillators

The space-borne observatories CoRoT (Convection Rotation and planetary Transits) and Kepler have provided photometric time series data of unprecedented precision for large numbers of stars. These data have revolutionized the fields of transiting exoplanets and asteroseismology. In this review some important asteroseismic results obtained using data from the CoRoT and Kepler space missions concerning stars that show solar-like oscillations are discussed. These results comprise, among others, measurements of the location of the base of the convection zone and helium second-ionization zone in main-sequence stars, the presence (or not) of core-helium burning in red-giant stars, as well as differential rotation in these stars.     

Radioactivity of the Moon and planets

In this report the main results of the study of radioactivity of the solar sistem bodies are considered. The radioactivity of the Moon and planets was measured by means of vehicles in situ. The radioactivity of the lunar samples, brought to the Earth was studied with laboratory equipment.     

Return to Europa: Overview of the Jupiter Europa orbiter mission

Missions to explore Europa have been imagined ever since the Voyager mission first suggested that Europa was geologically very young. Subsequently, the Galileo spacecraft supplied fascinating new insights into this satellite of Jupiter. Now, an international team is proposing a return to the Jupiter system and Europa with the Europa Jupiter System Mission (EJSM). Currently, NASA and ESA are designing two orbiters that would explore the Jovian system and then each would settle into orbit around one of Jupiter’s icy satellites, Europa and Ganymede. In addition, the Japanese Aerospace eXploration Agency (JAXA) is considering a Jupiter magnetospheric orbiter and the Russian Space Agency is investigating a Europa lander.     

基于最小值原理的冗余度机器人动力学控制

为了解冗余度机器人全局法优化中数值求解的困难,本文讨论了动力学方程的建立,无约束和有约束最优控制问题之间的内在联系,重点分析了求解最优控制问题的数值方法,文中提出了双向异步积分迭代的求解正则方程组的直接迭代法,较好解决了状态方程和协态方程稳定相逆给求解两点这值问题带来的困难。     

基于GPS载波相位的微卫星姿态算法研究

GPS姿态系统是利用GPS载波相位测量来确定载体的航向和姿态角。本文对微卫星的GPS姿态系统进行了研究。重点解决短基线的GPS状态算法，讨论了姿态价格函数的了小化方法。测试结果表明，所提的算法对小于1m的基线是有效的。