极月轨道上太阳电池阵外热流规律分析

准确的月球表面温度分布模型对于开展月球探测具有重要意义. 目前有关月球表面温 度模型还缺乏对完整月球表面温度分布的计算方法研究. 本文建立了一套计算完整月球表面温度的方法, 其中月球阳面温度采用Racca模型直接计算得到; 对于月球阴面, 将其沿纬度方向划分为若干区域, 每个区域的地表土壤采用一维非稳态热传导模型, 根据嫦娥三号着陆器太阳电池阵在轨环月阶段的温度数据, 修正得到月球表面土壤导热系数、密度及比热容, 通过数值计算求解一维非稳态热传导方程, 得出任意时刻月球阴面表面温度随时间的变化. 嫦娥三号着陆器太阳电池阵环月阶段热分析结果与在轨温度符合较好, 初步说明本文建立的完整月球表面温度计算方法正确可行. 基于本文方法计算得到整个月球表面温度分布, 进一步研究了极月轨道太阳电池阵外热流变化规律.      

基于远程渲染的三维模型发布系统

当前三维模型在Internet上的应用越来越广,但先下载再显示的传统发布方式也一直受到安全性、网络带宽和客户端机器的图形性能3方面的制约.结合三维模型化简和远程渲染技术提出了一种新的三维模型发布系统.系统把高精度模型保存在渲染服务器上,客户端的三维浏览器只提供大幅精简的粗糙模型.用户可以在客户端对粗糙模型进行缩放、平移和旋转等操作,选择最佳的观察视点.当用户停止操作时,客户端向服务器发送一个渲染请求,该请求中包含用户当前的视点、方向和光照条件等信息.服务器根据客户端请求进行渲染并返回高精度模型的图像.因为敏感的高精度模型保存在服务器端,因此该系统具有先天的安全性优势;因为客户端和服务器之间只在用户停止操作时传输平均大小为30 kByte的图像,因此系统对网络带宽要求很低;此外客户端只需要渲染极为简单的精简版模型,所以不需要很高的图形性能.      

月球受光面上月尘静电浮扬特性分析

基于单粒子轨道理论及空间尘埃等离子体充电方程,建立了月球受光面上尘埃微粒的静态荷电模型.基于光电子能量Maxiwellian分布假设,确定了月面垂直空间电场强度和光电子鞘层内带电粒子密度的函数表达式.利用牛顿运动定律和静电场力表达式,构建了月球受光面上尘埃微粒的静电浮扬动力学模型,并进行月尘静态浮扬特性的数值计算.研究结果显示:太阳高度角与颗粒粒径是控制月尘静电浮扬发生及动力学特性的两个基本参量;月尘静电浮扬发生在月球的黎明和黄昏;随着粒径的减少,月尘颗粒的最大浮扬高度不断增加.     

最小数据丢失量的地月中继卫星任务调度研究

月球背面的探测器必须依靠地月中继卫星进行数据传输与通信。地月中继任务包括实时性任务和延迟容忍类任务,如数传任务。当探测器等待传输的数据量超出用户存储容量时,延迟容忍类任务会由于探测器本地存储资源不足和地月中继卫星天线资源受限而无法完成,导致任务数据丢失,所以需要设计一种合理的地月中继任务调度策略,提高地月中继卫星的资源利用率,减少数据的丢失。对地月中继卫星任务调度进行了研究,在分析地月中继卫星数传任务的特点及用户的存储限制的基础上,以最小化数据丢失量为优化目标,建立了地月中继卫星任务调度模型,并设计了一种基于离散烟花算法（DFWA）的地月中继卫星任务调度算法。仿真数据分析表明,基于离散烟花算法的地月中继卫星任务调度算法在求解结果上优于遗传算法,是一种合理、有效的调度方法。     

月球探测器垂直软着陆4D轨道全局优化

提出一种包含动态规划法与共轭梯度法的组合优化算法,求解月球探测器垂直软着陆问题. 其以动态规划法中求得的次优控制变量作为共轭梯度法的初始控制变量,求得更为精确的最优控制变量和最优轨道. 月球垂直软着陆的轨道可以分为两段,即制动段与着陆段. 以燃料消耗量最小为性能指标,采用该组合优化算法分别对两段轨道进行了4D全局优化. 数值仿真结果表明,该组合算法优化精度较高,收敛速度快,稳定性好,可用于机载计算机实时生成垂直软着陆4D轨道,同时还可推广到其他终端时间自由型两点边值问题.      

Effect of surface roughness on microwave brightness temperature from lunar surface: Numerical analysis with a hybrid method

A hybrid method, combining the radiative transfer theory and the method of moments (MoM), is proposed to study the potential effect of the lunar surface roughness on the microwave brightness temperature. The total upward emission reaching the lunar surface from below media is calculated by the radiative transfer theory, and then the brightness temperature is obtained by weighting the bidirectional transmission coefficients which is computed using the MoM. The method is validated by both flat and rough surface models with analytic solutions. With the hybrid method, brightness temperatures from simulated lunar model are calculated and compared to those from a flat layered model. The comparisons show that the effect of rough surface on brightness temperature cannot be ignored and also depends on many other factors, such as observation angle and polarizations. For vertical polarization, an optimal observation angle may exist to reduce the effect of surface roughness. These results indicate that the knowledge of lunar surface roughness is important in microwave remote sensing to the Moon and may probably provide a guide to lunar projects in future.     

Plasma electron temperature variability in lunar surface potential and in electric field under average solar wind conditions

The Moon is immersed in plasma environment. The most interesting challenge of the lunar plasma– field environment is that it is alternatively dominated by the extended but variable outer atmosphere of the Earth – the magnetosphere – and by the extended but highly variable solar atmosphere – the solar wind. Understanding the plasma environment and its interaction with the lunar surface will be beneficial to both manned and robotic surface exploration activities and to scientific investigations. Presented is a preliminary map of variations of lunar surface electric potential over the day side and night side using probe equations and a discussion on dust dynamics in this E-field structure using the data from Electron Reflectometer in Lunar Prospector spacecraft during 1998–1999. On the day side, potential is around 5 V and on the night side it reaches up to −82 V. On the night side region, only highly energetic electrons can overcome this large negative potential. The variation at electron temperature (T_e) strongly reflects in the surface potential. The potential reaches to a value of −82 V for T_e = 58 eV. Surface charging causes the electrostatic transport of charged dust grains. Dust grain size of 0.1 μm shows a levitation height of 4.92 m on lunar day side, 748 m on terminator region and 3.7 km on the night side. The radius of maximum sized grain to be lofted, R_max, peaks at the terminator region (R_max = 0.83 μm). At the transition region dust levitation is almost absent. This region is most suited for exploration activities as the region is free from hazards caused by lunar dust.     

月地转移轨道快速设计方法

月地转移轨道设计一般分为初步轨道设计和精确轨道设计.其中,初步轨道设计的准确性是确保后续精确轨道设计收敛的关键.提出了一种基于Lambert算法的月地转移轨道快速设计方法.以出月球影响球的时刻、位置和速度为中间变量,将轨道分为地心段和月心段分别进行计算.将探测器飞出月球影响球至指定再入点的地心段轨道简化为一个Lambert问题进行求解,提出了通过牛顿迭代法求解月地转移轨道Lambert问题的方法,避免了Lambert问题求解时大量的超几何函数和级数计算,提高了计算效率.在月心段轨道的快速计算中,提出了根据探测器出影响球速度矢量、月球停泊轨道倾角和近月点高度计算月心双曲线轨道根数的新方法.通过迭代计算,使得两段轨道在月球影响球处的位置和速度连续,从而获得一条完整的满足两端约束的双二体月地转移轨道.该方法计算速度快,精度相对较高.计算结果可以作为后续精确轨道设计的初值.      

Lunar gravity field modeling critical analysis and challenges

This paper summarizes and provides a critical analysis of the historical developments of lunar gravitational models from the earliest use of ground based tracking systems of the Lunar Orbiter to the Lunar Prospector mission. This encompasses a comprehensive and critical analysis of the various methods used in the estimation of the gravity coefficients and the processing of large batches of diverse measurements and data types. It has been shown that weakness exists in the current models of the lunar gravity field, which is primarily due to the lack of far side lunar tracking data information, which makes the lunar potential modeling difficult but expected to be overcome as data from SELENE satellite-to-satellite tracking becomes available. Comparisons of various lunar models reveal an agreement in the low order coefficients of the spherical harmonics. However, substantial differences in the models exist in the higher-order harmonics. A numerical comparison has been presented showing the performance of all the contemporary lunar gravitational models used within the astrodynamics community and available in public domain. Improvements to the current models are part of a continuing process and the recent model improvements and future possibilities in lunar gravity modeling are discussed. A brief review of the recent missions has been presented. It is hoped that this critical review will benefit the researchers by presenting the historical as well as state of the art in this field.     

Overview of the Latest Scientific Results of China's Lunar Exploration Program

China's Chang'E-4 probe successfully landed on 3 January 2019 in Von Kármán crater within the South Pole-Aitken (SPA) basin on the lunar far side. Based on the data acquired by the scientific payloads onboard the lander and the rover, the researchers obtained the related information such as the geologic and tectonic setting of the landing area, compositional characteristics of the landing surface materials, dielectric permittivity and density of the lunar soil. The experiments confirmed the existence of materials dominated by olivine and low-calcium pyroxene in the SPA basin on the lunar far side, which preliminary revealed the geological evolution history of the SPA basin and even that of the early time lunar crust, as well as the tectonic setting and formation mechanism of the materials in the lunar interior. The researchers also inves-tigated the particle radiation, Linear Energy Transaction (LET) spectrum, and so forth on the lunar surface. The low-frequency radio observations were carried out on the lunar far side for the first time as well. This article summarizes the latest scientific results in the past years, focusing on the Chang'E-4 mission. Key words CLEP, Chang'E-4, Scientific objectives, Scientific payloads, Scientific results      

月球天梯力学特性研究

月球天梯的构型对其力学特性和质量特性起决定性作用,直接关系到天梯建造材料的选取。目前对月球天梯的研究集中在组合横截面积构型上,未对其他构型天梯的特性展开深入研究。针对这一问题,文章进行了不同构型月球天梯的力学特性研究。建立了月球天梯动力学模型,分别将天梯缆绳线密度和应力作为常值,推导了等横截面积和变横截面积两种新构型月球天梯的动力学方程和质量特性方程,进行了不同构型天梯关键特性参数的计算和对比分析。研究结果表明,变横截面积天梯可将最大应力设计为任意恒定值,降低了对材料极限应力的要求,其总质量小于等横截面积和组合横截面积天梯,为最优构型方案。     

发动机羽流作用下的实时月尘运动仿真

月尘运动是月球探测器软着陆过程中不可或缺的重要环节,针对发动机羽流作用下月尘运动真实感不强和月尘颗粒运动模型过于简单的问题,提出了一种逼真的实时月尘运动仿真方法.通过计算流体动力学(CFD, Computational Fluid Dynamics)和二次谢别德插值(Quadratic Shepard)方法,分析和计算单个月尘颗粒的运动学模型,得到一定初始条件下粒子运动的二维轨迹曲线;通过分析粒子的数量、初始位置、初始速度、生命周期等参数对粒子运动学的影响和变化规律,抽象出月尘系统的粒子集;建立基于月尘粒子集的月尘运动模型.实验结果显示:该运动模型逼真的模拟了发动机羽流作用下月尘腾起、飞溅、弥漫、消散等运动过程,视觉真实感和实时性能良好,对研究真空环境中的月尘运动及月球软着陆等相关领域具有一定的参考意义.目前该方法已应用于北航虚拟现实国家重点实验室月球软着陆仿真系统.     

环月轨道交会的载人登月任务轨道与窗口规划

针对环月轨道(Low Lunar Orbit,LLO)共面交会支持的"人货分离"载人登月任务,提出了一种任务窗口与轨道一体化规划方法。分析了基于LLO共面交会的"人货分离"载人登月任务的基本流程和工程约束;针对任务各阶段窗口与轨道求解问题,提出了以动力下降时刻为迭代初值的窗口规划策略,并建立了高精度模型下的环月轨道、双二体模型下的人员和货物运输轨道规划模型。以载人月球探测中国科学家命名的环形山为假想背景,给出仿真实例,仿真结果验证了文章所提方法的正确性,为探月工程任务提供了一种有效的窗口与轨道设计工具。     

三层月壤模型的多通道微波辐射模拟与月壤厚度的反演

由月球表面数字高程试验性地构造了整个月球表面月壤厚度的分布.根据Clementine探月卫星的紫外-可见光光学数据,计算了整个月球表面月壤中FeO+TiO2含量分布,给出了整个月球表面月壤介电常数分布.由月球表层温度的观测结果以及月壤的导热特性,给出了月尘层与月壤层温度随纬度分布的经验公式.在这些条件的基础上,建立了月尘、月壤、月岩三层微波热辐射模型.由起伏逸散定理,模拟计算了该月球模型多通道辐射亮度温度.然后,以此辐射亮度温度模拟加随机噪声为理论观测值,按三层模型提出了月壤层厚度反演方法.由于高频通道穿透深度小,由高频通道的辐射亮度温度按照两层月尘-月壤微波热辐射模型反演月尘层与月壤层的物理温度,再由穿透深度较大的低频通道辐射亮度温度反演月壤层厚度.对于反演的相对误差也进行了讨论.      

Lunar geophysical parameters inversion based on gravity/topography admittance and particle swarm optimization

This paper presents a FORTRAN computer program. The program as code will be used for lunar parameter inversions based on gravity/topography admittance. This will be done by assuming that the lunar lithosphere is modeled as a thin elastic spherical shell. The parameters discussed here include; load ratio, crustal thickness, subsurface load depth, crustal density and elastic lithosphere thickness. The admittance of the best-fitting model can be found through automatically adjusting misfits between one theoretical admittance and an observed one. The results in this paper indicate that this research’s theoretical model is reasonable for exploring the best-fitting parameters. In addition, this code is not only able to automatically and simultaneously calculate the global optimum solution of the parameters studied, but also performs well in computational speed. The code can be easily modified to include more parameter inversions; such as the inversion for subsurface density anomaly and the case of considering infilling material in some lunar mare basins.     

Night side lunar surface potential in the Earth’s magnetosphere

In the Earth’s magnetotail, Japanese Moon orbiter Kaguya repeatedly encountered the plasmoid or plasma sheet. The encounters were characterized by the low energy ion signatures including lobe cold ions, cold ion acceleration in the plasma sheet-lobe boundaries, and hot plasma sheet ions or fast flowing ions associated with plasmoids. Different from the previous observations made in the magnetotail by the GEOTAIL spacecraft, the ions were affected by the existence of the Moon. On the dayside of the Moon, tailward flowing cold ions and their acceleration were observed. However, on the night side, tailward flowing cold ions could not be observed since the Moon blocked them. In stead, ion acceleration by the spacecraft potential and the electron beam accelerated by the potential difference between lunar surface and spacecraft were simultaneously observed. These electron and ion data enabled us to determine the night side lunar surface potential and spacecraft potential only from the observed data for the first time.     

Comparison analyses on the 150 × 150 lunar gravity field models by gravity/topography admittance,correlation and precision orbit determination

We analyzed the 150 × 150 lunar gravity field models, LP150Q, GLGM-3 and SGM150, using the power spectrum on the lunar nearside and farside, the lunar global and localized gravity/topography admittance and correlation, and Chang’E-2 precision orbit determination to investigate which model is a more effective tool to estimate geophysical parameters and determine the lunar satellite precision orbit. Results indicate that all gravity field models can be used to estimate the lunar geophysical parameters of the nearside of the Moon. However, SGM150 is better in such computation of the farside. Additionally, SGM150 is shown to be the most useful model for determining the lunar satellite orbit.     

Research on the evaluation method of lunar soil disturbance based on the fitting of Gaussian surface

It is crucial to maintain the layer sequence information on in situ lunar soil during the lunar soil sampling. In this paper, a discrete element model for direct push sampling of the lunar weathering layer is established, and the relationship between the perturbation rate of the lunar soil layer sequence information and the sampling depth is obtained. A quantitative evaluation method of the perturbation rate of the structural information on the lunar soil layer sequence based on the Gaussian curvature is proposed by reconstructing the center-of-mass point cloud model of the marked layer lunar soil particles in MATLAB. The proposed method can accurately describe the perturbation rate of lunar soil layer sequence information and can analyze the degree of lunar soil disturbance according to the displacement of the sampling tube and mechanical parameters in real-time inversion during the sampling process. When the sampling depth is 40 mm, the minimum and maximum perturbation rates of sampling along the penetration direction of the machine are 5.261 % and 27.649 %, respectively. The perturbation rate shows a trend of first decreasing and then increasing. The method proposed in this study could effectively solve the problem that the disturbance of lunar soil bedding information was difficult to evaluate in the sampling process of lunar regolith.     

Design and analysis of an extended mission of CE-2: From lunar orbit to Sun–Earth L2 region

Chang’E-2 (CE-2) has firstly successfully achieved the exploring mission from lunar orbit to Sun–Earth L2 region. In this paper, we discuss the design problem of transfer trajectory and at the same time analyze the visible segment of Tracking, Telemetry & Control (TT&C) system for this mission. Firstly, the four-body problem of Sun–Earth–Moon and Spacecraft can be decoupled in two different three-body problems (Sun–Earth + Moon Restricted Three-Body Problems (RTBPs) and Earth–Moon ephemeris model). Then, the transfer trajectory segments in different model are computed, respectively, and patched by Poincaré sections. The full-flight trajectory including transfer trajectory from lunar orbit to Sun–Earth L2 region and target Lissajous orbit is obtained by the differential correction method. Finally, the visibility of TT&C system at the key time is analyzed. Actual execution of CE-2 extended mission shows that the trajectory design of CE-2 mission is feasible.