Attitude stability criteria of axisymmetric solar sail

Passive attitude stability criteria of a solar sail whose membrane surface is axisymmetric are studied in this paper under a general SRP model. This paper proves that arbitrary attitude equilibrium position can be designed through adjusting the deviation between the pressure center and the mass center of the sail. The linearized method is applied to inspect analytically the stability of the equilibrium point from two different points of views. The results show that the attitude stability depends on the membrane surface shape and area. The results of simulation with full dynamic equations confirm that the two stability criteria are effective in judging the attitude stability for axisymmetric solar sail. Several possible applications of the study are also mentioned.     

Influence of curved thin-film device on deformation of a solar sail

A spinning solar sail IKAROS’s membrane is estimated to unexpectedly deform into an inverted pyramid shape due to thin-film devices with curvature, such as thin-film solar cells and steering devices on the membrane. It is important to investigate the deformation caused by the curved thin-film devices and predict the sail shape because the out-of-plane deformation greatly affects solar radiation pressure (SRP) and SRP torque. The purpose of this paper is to clarify the relationship between the global shape and orientation and position of curved thin-film devices and to evaluate SRP torque on the global shape using finite element analysis. The global shape is evaluated based on the out-of-plane displacement and the SRP torque. When the curved thin-film devices make the membrane shrink in the circumferential, diagonal, and radial direction, the sail deforms into a pyramid shape, an inverted pyramid one, and a saddle one, respectively. The saddle shape is more desirable for solar sails than the inverted pyramid shape and the pyramid one from the viewpoint of shape stability to SRP and control of SRP torque in the normal direction of the sail (windmill torque). The position of the thin-film device tends to increase the absolute value of windmill torque when it is biased circumferentially from the petal central axis. The suggested design principles for the arrangement of thin-film devices is that the curved thin-film devices should be directed so that the sail shrinks in the radial direction in order to deform the sail into a saddle shape with high shape stability, and the position of the thin-film devices should be biased in the circumferential direction paying attention to the absolute value of windmill torque to determine the direction of windmill torque.     

基于目标特性的航天器太阳光压面积计算方法

太阳光压是影响深空探测航天器轨道确定与预报精度最主要的摄动力.针对实际任务需求,采用了一种基于目标特性的光压面积建模与计算方法,根据航天器形状、尺寸、表面材料以及材料光学特性等信息,实现了分析型光压模型的建立与求解,提高了计算效率和精度,可快速计算目标在光照方向上的光压面积、投影面积以及光压比例因子等参数.通过长方体光压面积理论值与仿真值的对比,验证了该方法的准确性和有效性.针对复杂结构探测器开展了光压面积计算,可为深空探测航天器精密定轨中的光压模型解算、定轨及预报提供参考.      

Attitude control actuator scaling laws for orbiting solar reflectors

The rapid evolution of in-orbit manufacturing will enable the fabrication of low-cost, large-scale space structures. In particular, the use of 3D printing technologies will remove traditional payload constraints associated with launch vehicles, due to fairing size and launch loads, thus allowing the construction of larger and lighter structures, such as orbiting solar reflectors. These structures will require efficient attitude control systems, able to provide the necessary torque for maneuvers and to counteract perturbations, such as gravity gradient and solar radiation pressure. In this paper, a top-level overview of actuator performances for orbiting solar reflectors is provided, and scaling laws associated with the required actuator mass and input power are developed. For each class of actuator, upper bounds on the maximum size of the structure which can be effectively controlled are presented. The results can also be extended to other classes of large planar Earth-pointing structures such as solar power satellites, solar sails, or large antennae.     

通气位置对潜射航行体流体动力特性影响分析

基于均质平衡流理论,通过求解混合介质的RANS方程、SST湍流输运方程和各相之间的质量输运方程,开展了通气位置对潜射航行体流体动力特性影响的三维数值模拟研究,对比分析了不同通气位置条件下空泡形态特性、表面压力分布以及阻力变化特性.结果表明:当两个通气口间距增加到一定值后,在两个通气口之间区域空泡发生断裂,在断裂区域迎流面和背流面表面压力值升高,且在空泡断裂闭合位置出现压力峰值;在相同通气量条件下,随着通气口间距的增大,压差阻力系数和粘性阻力系数均呈减小趋势变化.     

日冕扰动传播的数值模拟——(Ⅱ)热传导和辐射的影响

利用包含完全能量方程的磁流体力学方程组,数值研究了磁静力平衡基态下由于日面冷物质径向喷射所引起的日冕动力学响应。这些结果可以解释日珥的形成及伴随的日冕扰动传播。结果还表明,热传导和辐射损失对扰动传播的基本性质影响很小;对于日珥的形成,也未显出重要作用。      

单元构架式抛物面天线张紧绳索多层设计方法

星载抛物面天线在轨运行时由于姿态调整或环境因素变化,可能产生振动,从而降低其工作性能。为了减弱振动的影响,需采取加装张紧绳索的方法来提升天线刚度。提出了一种新型张紧绳索多层设计方法,通过将抛物面天线划分为多层,调整各层绳索的张紧力,从而提高天线结构刚度,同时尽量降低张紧力引起的形面误差。为了验证所提方法的有效性,建立了天线的有限元模型。在此基础上进行了有限元分析,研究了不同张紧力参数配置下的结构刚度和形面误差,并以提升结构刚度的同时降低张紧力引起的形面误差为目标,对张紧力参数进行了优化设计。为了提高计算效率,采用响应面法建立代理模型参与优化迭代计算。采用非劣排序遗传算法（NSGA-Ⅱ）完成优化迭代计算,优化后的张紧力参数使天线的性能得到了进一步的提升,为构架式可展开抛物面天线的设计提供理论指导。      

The feasibility of using an L1 positioned dust cloud as a method of space-based geoengineering

In this paper a method of geoengineering is proposed involving clouds of dust placed in the vicinity of the L₁ point as an alternative to the use of thin film reflectors. The aim of this scheme is to reduce the manufacturing requirement for space-based geoengineering. It has been concluded that the mass requirement for a cloud placed at the classical L₁ point, to create an average solar insolation reduction of 1.7%, is 7.60 × 10¹⁰ kg yr⁻¹ whilst a cloud placed at a displaced equilibrium point created by the inclusion of the effect of solar radiation pressure is 1.87 × 10¹⁰ kg yr⁻¹. These mass ejection rates are considerably less than the mass required in other unprocessed dust cloud methods proposed and are comparable to thin film reflector geoengineering requirements. Importantly, unprocessed dust sourced in-situ is seen as an attractive scheme compared to highly engineered thin film reflectors. It is envisaged that the required mass of dust can be extracted from captured near Earth asteroids, whilst stabilised in the required position using the impulse provided by solar collectors or mass drivers used to eject material from the asteroid surface.     

Plasma polarization electric field derived from radio sounding of solar wind acceleration region with spacecraft signals

Presented is the analytical approximation of averaged solar wind velocity radial dependence in the solar wind acceleration region at heliolatitudes below 60° under low and moderate solar activity. This empirical approximation is based on the data of radio sounding of the solar corona with radio signals from various spacecraft. Deduced is an equation connecting the solar wind velocity radial dependence and the radial dependence of solar wind plasma polarization electric field intensity. This allows constructing a semi-empirical radial dependence of plasma polarization electric field corresponding to the empirical radial dependence of solar wind velocity. Main properties of the semi-empirical dependence, which is based on radio sounding data, are described.     

基于流动控制的水下航行体流体动力特性分析

利用CFX对水下航行体流体动力特性进行了数值模拟研究.求解了混合介质的剪切压力输运湍流模型控制方程、雷诺-平均奈维尔-斯托克斯方程以及各相间的质量输运方程,采用三维数值模拟方法对比分析了不同流动控制方案的水下航行体空泡形态特征、表面压力分布和阻力系数变化特性,讨论了不同参数对减阻效果的影响.结果表明,对于有横流作用的水下垂直发射航行体,多相流动控制不仅可以降低通气空泡的不对称性和航行体阻力,同时可以均匀迎流面和背流面压力,从而实现航行体多相流空泡形态及表面压力特性等的控制.此外,通气口的位置对减阻效果具有显著影响.      

太阳帆飞行器轨道动力学分析

通过分析对轨道要素影响最大的加速度分量,使太阳帆总是位于光压力沿这个加速度方向的分量最大的方位。通过分析轨道要素调整时的相互影响关系,提出了同时修正多种轨道偏差的控制方案。对处于地球静止轨道上太阳帆飞行器的轨道调整进行了数值仿真。结果证明利用太阳光压力进行轨道调整是可行的,而且有利的太阳方位是进行快速有效的轨道调整的必要条件。     

空泡摆动对超空泡航行体尾拍影响分析

针对超空泡航行体存在尾拍运动与空泡形态变化互相耦合作用的特点,建立了超空泡航行体尾拍耦合运动方程,并对超空泡摆动对超空泡航行体尾拍相互作用过程进行了数值研究,对比分析了不同速度和角速度的超空泡航行体运动过程及空泡摆动对尾拍的影响规律.研究结果表明:空泡摆动使航行体转角、整体运动周期、尾拍升力和阻力增大;在不同速度下,空泡摆动对航行体转角、整体运动周期及尾拍周期的影响随速度增大而逐渐减小;不同初始角速度下,空泡摆动对航行体转角、整体运动周期、航行体尾拍升力和阻力的影响随角速度增大而逐渐增大.     

Sailing with solar and planetary radiation pressure

Literature on solar sailing has thus far mostly considered solar radiation pressure (SRP) as the only contribution to sail force. However, considering a sail in a planetary mission scenario, a new contribution can be added. Since the planet itself emits radiation, this generates a radial planetary radiation pressure (PRP) that is also exerted on the sail. Hence, this work studies the combined effects of both SRP and PRP on a sail for two case studies, i.e. Earth and Venus. In proximity of the Earth, the effect of PRP can be significant under specific conditions. Around Venus, instead, PRP is by far the dominating contribution. These combined effects have been studied for single- and double-sided reflective coating and including eclipse. Results show potential increase in the net acceleration and a change in the optimal attitude to maximise the acceleration in a given direction. Moreover, an increasing semi-major axis manoeuvre is shown with and without PRP, to quantify the difference on a real-case scenario.     

Active shape control for flexible space structures using an optimal gyricity distribution

An active shape control approach for a circular flexible space structure is studied. The shape of these flexible structures may need active shape control due to some particular mission objectives. For example, antenna reflectors may have requirements for their shape accuracy to guarantee communication performance; and some solar sails may change their shape to realize different commissions during interplanetary flight. This paper investigates an active shape control process executed by the gyricity (stored angular momentum) distribution existing in the circular flexible space structure. The shape of the circular flexible space structure is expected to be reformed from a plane plate to a paraboloidal plate. An optimal gyricity distribution is derived based on the optimal control theory for systems described by partial differential equations. Distributed forces can be generated by the optimal gyricity distribution in a rotating field to implement the active shape control. The dynamic model of the circular flexible space structure is established for the active shape control through the finite element method. Numerical simulation validates the functionality of the shape control method, and illustrates the shape-adjusting process of the circular flexible space structure.     

航天器在轨道运行中太阳帆板的温度问题

本文应用能量平衡的基本原理, 获得了航天器在几种典型热工况中太阳帆板的能量平衡的基本方程以及宇宙空间的热辐射角系数基本方程.利用数值积分法, 得到了角系数的数值以及太阳帆板在不同的轨道所承受的最大温度波动, 在考虑这一温度效应后, 用有限元法分析计算了太阳能电池帆板的固有频率特性.计算结果表明, 航天器在空间轨道运行时, 通过阴区和阳区所产生的温度变化对太阳帆板的动力学特性影响是不容忽略的问题.      

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

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

嫦娥三号着陆区月表辐射特性测算

月球表面辐射特性研究是月球探索的一个重要领域,对于勘测月球表面形态及进行资源勘探等具有重要意义,但是目前在月球表面直接进行近距离测量的数据非常少.本文建立了基于嫦娥三号着陆器地形地貌相机和缓释电机在轨遥测数据研究月球表面辐射特性的方法,初步测算得到嫦娥三号着陆区月球表面的太阳辐射光谱反射率为0.105,发射率为0.866.所得结果能够为进一步空间探测任务提供原始测算数据,以降低设计难度及提高探测能力.      

Computer modelling of a penetrator thermal sensor

The Philae lander is part of the Rosetta mission to investigate comet 67P/Churyumov-Gerasimenko. It will use a harpoon like device to anchor itself onto the surface. The anchor will perhaps reach depths of 1–2 m. In the anchor is a temperature sensor that will measure the boundary temperature as part of the MUPUS experiment. As the anchor attains thermal equilibrium with the comet ice it may be possible to extract the thermal properties of the surrounding ice, such as the thermal diffusivity, by using the temperature sensor data. The anchor is not an optimal shape for a thermal probe and application of analytical solutions to the heat equation is inappropriate. We prepare a numerical model to fit temperature sensor data and extract the thermal diffusivity. Penetrator probes mechanically compact the material immediately surrounding them as they enter the target. If the thermal properties, composition and dimensions of the penetrator are known, then the thermal properties of this pristine material may be recovered although this will be a challenging measurement. We report on investigations, using a numerical thermal model, to simulate a variety of scenarios that the anchor may encounter and how they will affect the measurement.     

Solar atmosphere wave dynamics generated by solar global oscillating eigenmodes

The solar atmosphere exhibits a diverse range of wave phenomena, where one of the earliest discovered was the five-minute global acoustic oscillation, also referred to as the p-mode. The analysis of wave propagation in the solar atmosphere may be used as a diagnostic tool to estimate accurately the physical characteristics of the Sun’s atmospheric layers.In this paper, we investigate the dynamics and upward propagation of waves which are generated by the solar global eigenmodes. We report on a series of hydrodynamic simulations of a realistically stratified model of the solar atmosphere representing its lower region from the photosphere to low corona. With the objective of modelling atmospheric perturbations, propagating from the photosphere into the chromosphere, transition region and low corona, generated by the photospheric global oscillations the simulations use photospheric drivers mimicking the solar p-modes. The drivers are spatially structured harmonics across the computational box parallel to the solar surface. The drivers perturb the atmosphere at 0.5?Mm above the bottom boundary of the model and are placed coincident with the location of the temperature minimum. A combination of the VALIIIC and McWhirter solar atmospheres are used as the background equilibrium model.We report how synthetic photospheric oscillations may manifest in a magnetic field free model of the quiet Sun. To carry out the simulations, we employed the magnetohydrodynamics code, SMAUG (Sheffield MHD Accelerated Using GPUs).Our results show that the amount of energy propagating into the solar atmosphere is consistent with a model of solar global oscillations described by Taroyan and Erdélyi (2008) using the Klein-Gordon equation. The computed results indicate a power law which is compared to observations reported by Ireland et al. (2015) using data from the Solar Dynamics Observatory/Atmospheric Imaging Assembly.     

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

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