A procedure for benchmarking laboratory exposures with 1 A GeV iron ions.

A new version of the HZETRN code capable of simulating HZE ions with either laboratory or space boundary conditions is under development. The computational model consists of combinations of physical perturbation expansions based on the scales of atomic interaction, multiple scattering, and nuclear reactive processes with use of asymptotic/Neumann expansions with non-perturbative corrections. The code contains energy loss with straggling, nuclear attenuation, nuclear fragmentation with energy dispersion and downshifts, and off-axis dispersion with multiple scattering under preparation. The present benchmark is for a broad directed beam for 1 A GeV iron ion beams with 2 A MeV width and four targets of polyethylene, polymethyl metachrylate, aluminum, and lead of varying thickness from 5 to 30 g/cm2. The benchmark quantities will be dose, track averaged LET, dose averaged LET, fraction of iron ion remaining, and fragment energy spectra after 23 g/cm2 of polymethyl metachrylate.     

An estimation of tropospheric corrections using GPS and synoptic data: Improving Urmia Lake water level time series from Jason-2 and SARAL/AltiKa satellite altimetry

Tropospheric correction is one of the most important corrections in satellite altimetry measurements. Tropospheric wet and dry path delays have strong dependence on temperature, pressure and humidity. Tropospheric layer has particularly high variability over coastal regions due to humidity, wind and temperature gradients. Depending on the extent of water body and wind conditions over an inland water, Wet Tropospheric Correction (WTC) is within the ranges from a few centimeters to tens of centimeters. Therefore, an extra care is needed to estimate tropospheric corrections on the altimetric measurements over inland waters. This study assesses the role of tropospheric correction on the altimetric measurements over the Urmia Lake in Iran. For this purpose, four types of tropospheric corrections have been used: (i) microwave radiometer (MWR) observations, (ii) tropospheric corrections computed from meteorological models, (iii) GPS observations and (iv) synoptic station data. They have been applied to Jason-2 track no. 133 and SARAL/AltiKa track no. 741 and 356 corresponding to 117–153 and the 23–34 cycles, respectively. In addition, the corresponding measurements of PISTACH and PEACHI, include new retracking method and an innovative wet tropospheric correction, have also been used. Our results show that GPS observation leads to the most accurate tropospheric correction. The results obtained from the PISTACH and PEACHI projects confirm those obtained with the standard SGDR, i.e., the role of GPS in improving the tropospheric corrections. It is inferred that the MWR data from Jason-2 mission is appropriate for the tropospheric corrections, however the SARAL/AltiKa one is not proper because Jason-2 possesses an enhanced WTC near the coast. Furthermore, virtual stations are defined for assessment of the results in terms of time series of Water Level Height (WLH). The results show that GPS tropospheric corrections lead to the most accurate WLH estimation for the selected virtual stations, which improves the accuracy of the obtained WLH time series by about 5%.     

对称面圆周角对轴对称降落伞流场特性的影响

根据降落伞的特点,通过3点假设(伞衣薄膜、伞的轴对称和流场定常),将三维复杂流动问题转化成二维轴对称问题,以节约计算时间.定义对称面圆周角,保持伞衣幅底部直径和顶孔直径不变,选取对称面圆周角在80°~140°范围内变化,建立一系列轴对称降落伞的计算模型.利用数值模拟手段,求解RNG (Renormalization Group)k-epsilon湍流模型下的N-S方程组,获得与有关单位试验相吻合的计算结果.分析发现对称面圆周角和伞衣幅高度对降落伞阻力影响很小.算例中阻力随对称面圆周角的变化在±0.28%以内.对称面圆周角的变化对轴对称降落伞尾流区流场的拓扑结构没有影响.对称轴上存在2个鞍点,随着对称面圆周角的增大,第1个鞍点的位置几乎不变,第2个鞍点的位置向尾流方向推移.      

卫星姿态测量信息的二阶非线性滤波技术研究

卫星应用任务要求卫星姿态在空间保持高精度定向 ,为此有必要研究精确的递推非线性姿态估计算法。文中针对星光 陀螺这种典型的三轴稳定卫星姿态确定系统模式 ,应用基于二阶泰勒级数近似得到的非线性滤波算法处理姿态测量信息 ,设计了二阶非线性姿态估计器的实现方案。在设计过程中 ,把QUEST法作为矢量观测数据压缩技术有效地结合进姿态估计器中 ,使得在多矢量观测情况下的滤波修正算法得到了简化。相同条件下的仿真测试结果证明 ,二阶非线性姿态估计器的滤波性能要优于采用扩展卡尔曼滤波技术的一阶线性姿态估计器     

Ionospheric Doppler correction model for satellite orbit determination

The use of satellites for navigation, geodesy and altimetry requires a precise knowledge of their orbits. Satellite orbits are commonly derived from Doppler counts of radio signals propagating to/from the satellites. However, corrections must be applied to these measurements to account for the presence of the ionosphere. The first order ionospheric perturbation term is usually eliminated by the use of dual frequency systems. However, the higher order terms have to be taken into account. In this paper, a correction model is described for the second and third order perturbation terms. The entries of the model are the estimates of the first order terms given by dual frequency systems. On the basis of computer simulations, it is shown that a significant reduction of the total ionospheric error can be achieved with this technique.     

An improved Green's function for ion beam transport.

Ion beam transport theory allows testing of material transmission properties in the laboratory environment generated by particle accelerators. This is a necessary step in materials development and evaluation for space use. The approximations used in solving the Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are the main emphasis of the present work. In consequence, an analytic solution of the linear Boltzmann equation is pursued in the form of a Green's function allowing flexibility in application to a broad range of boundary value problems. It has been established that simple solutions can be found for high charge and energy (HZE) ions by ignoring nuclear energy downshifts and dispersion. Such solutions were found to be supported by experimental evidence with HZE ion beams when multiple scattering was added. Lacking from the prior solutions were range and energy straggling and energy downshift with dispersion associated with nuclear events. Recently, we have found global solutions including these effects providing a broader class of HZE ion solutions.     

解Stokes问题的区域分解算法

对用非协调混合元离散的Ｓｔｏｋｅｓ问题，建立了延拓定理并分析了一个两子域不重叠区域分解算法，利用已有的对二阶椭圆型问题有效的预处理器，借助Ｂｒａｍｂｌｅ ＪＨ等的思想，将原鞍点问题变为一个在适当内积意义下的对称正定问题，并估计了它的条件数。     

链路统计分析中导出参数概率计算方法研究

为了适应深空探测任务需求,有效利用测控资源,提高链路计算精度,对天地链路统计分析基本方法进行了研究,对CCSDS标准建议的3种链路参数的概率密度函数特性进行了系统分析,提出了一种当链路参数不满足李雅普诺夫（Lyapunov）条件时,利用鞍点逼近估计尾概率函数确定导出参数取值概率的方法。通过理论分析和实例计算表明,采用该方法得到的导出参数的概率分布特性与真实值基本一致,进一步完善了天地链路统计计算方法,使天地链路统计分析方法更具有普遍性。     

流动对太阳黑子精细结构的影响

考虑磁扩散率为温度的函数以及绝热条件,线性扰动磁流力学方程组的色散方程是频率的三次方程。由于沿磁场流动的速度u₀小于声速u_s,即ε=u₀2/u_s2是一个小量,在零级近似下解色散方程,得到了磁扩散不稳定性。本文进一步讨论在一级近似下小流动对磁扩散不稳定性的影响。结果表明,小流动不破坏磁扩散不稳定性的存在,并且,它对波数空间中的不稳定区也影响不大。      

27-day variation in solar-terrestrial parameters: Global characteristics and an origin based approach of the signals

The Earth and the near interplanetary medium are affected by the Sun in different ways. Those processes generated in the Sun that induce perturbations into the Magnetosphere-Ionosphere system are called geoeffective processes and show a wide range of temporal variations, like the 11-year solar cycle (long term variations), the variation of ～27?days (recurrent variations), solar storms enduring for some days, particle acceleration events lasting for some hours, etc.In this article, the periodicity of ～27?days associated with the solar synodic rotation period is investigated. The work is mainly focused on studying the resulting 27-day periodic signal in the magnetic activity, by the analysis of the horizontal component of the magnetic field registered on a set of 103 magnetic observatories distributed around the world. For this a new method to isolate the periodicity of interest has been developed consisting of two main steps: the first one consists of removing the linear trend corresponding to every calendar year from the data series, and the second one of removing from the resulting series a smoothed version of it obtained by applying a 30-day moving average. The result at the end of this process is a data series in which all the signal with periods larger than 30?days are canceled.The most important characteristics observed in the resulting signals are two main amplitude modulations: the first and most prominent related to the 11-year solar cycle and the second one with a semiannual pattern. In addition, the amplitude of the signal shows a dependence on the geomagnetic latitude of the observatory with a significant discontinuity at approx. ±60°.The processing scheme was also applied to other parameters that are widely used to characterize the energy transfer from the Sun to the Earth: F10.7 and Mg II indices and the ionospheric vertical total electron content (vTEC) were considered for radiative interactions; and the solar wind velocity for the non-radiative interactions between the solar wind and the magnetosphere. The 27-day signal obtained in the magnetic activity was compared with the signals found in the other parameters resulting in a series of cross-correlations curves with maximum correlation between 3 and 5?days of delays for the radiative and between 0 and 1?days of delay for the non-radiative parameters. This result supports the idea that the physical process responsible for the 27-day signal in the magnetic activity is related to the solar wind and not to the solar electromagnetic radiation.     

角区三维分离流附着鞍点拓扑结构及其演化

实验研究角区层流边界层三维定常分离附着鞍点拓扑结构的存在性及其演化规律.二维粒子图像测速(PIV,Particle Image Velocimetry)激光片光源设置于给定平板和不同模型组成的角区对称面,采用微距镜头捕捉对称面奇点附近的局部流动拓扑结构.实验证实在一定流动参数条件下多种形状模型与平板形成的角区均存在有别于经典分离的附着鞍点拓扑结构;随着雷诺数增加或者随着模型钝度的变化,附着鞍点拓扑结构与经典分离鞍点拓扑结构之间存在着一定的演化规律;对附着鞍点拓扑结构与Lighthill经典三维分离模式的关系进行了分析,表明Lighthill经典三维分离模式也适用于分析附着鞍点拓扑结构.     

角区三维漩涡流动的新分离结构

分别用数值计算和PIV(Particle Image Velocimetry)实验方法研究角区层流边界层三维定常分离的流动结构,证实了角区确实存在有别于传统分离现象的附着鞍点结构,对称面上游的流线并非从壁面向上抬起从壁面分离,而是经由空间某个奇点向壁面附着.角区马蹄涡的传统分离鞍点结构和附着鞍点结构之间存在着一定的演化规律,影响参数包括模型头部钝度、边界层厚度和雷诺数.一定的边界层和雷诺数条件下随着模型头部钝度减小,角区马蹄涡将从传统的分离鞍点结构逐步过渡为附着鞍点结构;对一定的模型若雷诺数越大、边界层越厚则角区流动越趋向于传统的分离鞍点结构,反之则倾向于附着鞍点结构.     

卫星热网络模型修正技术进展及其改进方法研究

文章对目前卫星热网络模型修正技术的进展情况进行了总结分析 ,针对存在的问题提出了几种有关基于最小二乘法的热网络模型修正技术的改进思路和设想 ,包括参数化修正技术、分级修正技术和瞬态修正技术改进方法等 ,以解决目前工况数据不足、约束条件多、工程实现困难等问题 ,为今后卫星热网络模型修正技术研究提供技术参考     

平动点任务转移轨道中途修正研究

平动点在深空探测中具有重要作用. 由于发射中不可避免地会存在误差, 平动点卫星往往需要考虑转移轨道的中途修正. 本文从工程实用角度出发, 根据以往平动点任务的中途修正经验, 确定中途修正序列,综合以往文献对误差的处理, 给出了误差的分布参数; 将转移轨道中途修正问题转化为转移轨道设计问题, 使用微分修正算法求解转移轨道中途修正问题; 利用蒙特卡洛模拟给出统计性数据. 仿真结果验证了方法的有效性. 对于日地系 110天类型的转移轨道, 在95 %的置信水平下, 只需要78.0 m/s的机动速度, 即可以保证位置误差在66.1 km之内.      

线阵CCD平场校正及FPGA实现的研究

为了消除由于照明光源的非均匀性、镜头渐晕及CCD传感器噪声对图像质量造成的不良影响,首先从理论上分析了线阵CCD固定模式噪声和像素光电响应非均匀性的数学模型和产生的原因;对像素光电响应非均匀性这一主要影响因素,根据其光电响应模型提出了CCD像素光电响应非均匀性的两点校正方法,即求解出各个像素光电响应不均匀性的校正因子和暗偏置量系数,进而用以计算各个像素的校正量,实现逐点校正;针对校正算法大数据量、实时计算的要求,设计了基于现场可编程逻辑门阵列FPGA快速实现校正算法的硬件处理结构。最后,通过仿真和实验表明,两点校正方法可以有效地减小CCD像素光电响应不均匀性的图像噪声,提高图像质量,以便做后续精确处理。     

Energy requirements of a thermally processed ISRU radiation shield for a lunar habitat

The purpose of this study was to establish, on a first principles basis, the order of magnitude of energy requirements for a thermally processed, lunar regolith radiation shield constructed using an in-situ resource utilisation (ISRU) approach. This was done by developing a reference scenario habitat and using thermodynamic relationships and specific heat capacity expressions to determine the energy required to bring such a regolith volume up to sintering temperatures (c. 1,375 K). Once the energy requirements were developed some power system architectures were outlined conceptually and a nuclear power plant of c. 400 kW was suggested as a means to supply the necessary energy. This is well beyond current space nuclear applications. The study concludes that it is likely that the most efficient near-term solution is chemical processing of regolith, from an energy requirements perspective. The technology is also more mature and likely to be delivered on near term projects as it does not require such scaled-up power system architectures. Alternatively, bringing storm shelters up with the habitat to provide a means of weathering major solar events, and adding additional radiation protection to habitat quarters, possibly through a water blanket or similar mechanism, could provide a non-ISRU solution with current technology. However, in the longer term, the development of MW-scale power system architectures (fission, solar etc.), may permit a very large volume of material to be processed thermally for construction material, making a large, permanent human presence on the Moon more easily realisable.     

Combining frequency and time domain models to forecast space weather

Space weather series incorporate several distinct components, cycles at multiple frequencies, irregular trends, and nonlinear variability. The cycles are stochastic, i.e., the amplitude varies over time. Similarly, the trend is stochastic: the slope and direction of trending change repeatedly. This study sets out a combined model using both frequency and time domain methods, in two stages. In the first stage, a frequency domain algorithm is estimated and forecasted. In the second stage, the forecast is used as an input in a neural network. The combined model also includes a term enabling the model to react inversely to large deviations between the actual values and forecast. The models are evaluated using two data sets, the hemispheric power data obtained from the Polar Orbiting Environment satellites, and the Aa geomagnetic index. All the series are at a daily resolution. Forecasting experiments are run over horizons of 1–7 days. The models are estimated using a moving window or adaptive approach. The combined model consistently achieves the most accurate results. Among single equation methods, the frequency domain model is more accurate for the geomagnetic index because it is able to capture the underlying cycles more effectively. In the hemispheric power series, the cycles are less pronounced, so that time domain methods are more accurate, except at very short horizons. Nevertheless, in both data sets, the combined model works well because the frequency domain algorithm captures cyclical behavior, while the neural net is better able to capture short-term dependence and trending.