Dual attitude and parameter estimation of passively magnetically stabilized nano satellites

The attitude determination capability of a nano satellite is limited by a lack of traditional high performance attitude sensors, a result of having small budgets for mass and power. Attitude determination can still be performed on a nano satellite with low fidelity sensors, but an accurate model of the spacecraft attitude dynamics is required. The passive magnetic stabilization systems commonly employed in nano satellites are known to introduce uncertainties in the parameters of the attitude dynamics model that cannot easily be resolved prior to launch. In this paper, a batch estimation problem is formulated that simultaneously solves for the attitude of the spacecraft and performs parameter estimation on the magnetic properties of the magnetic materials using only a measurement of the solar vector. The estimation technique is applied to data from NASA Ames Research Center's O/OREOS nano satellite and the University of Michigan's RAX-1 nano satellite, where clear differences are detected between the magnetic properties as measured before launch and those that fit the observed data. To date this is the first known on-orbit verification of the attitude dynamics model of a passively magnetically stabilized spacecraft.     

On Lambert’s problem and the elliptic time of flight equation: A simple semi-analytical inversion method

In the two-body model, time of flight between two positions can be expressed as a single-variable function and a variety of formulations exist. Lambert’s problem can be solved by inverting such a function. In this article, a method which inverts Lagrange’s flight time equation and supports the problematic 180°$180°$ transfer is proposed. This method relies on a Householder algorithm of variable order. However, unlike other iterative methods, it is semi-analytical in the sense that flight time functions are derived analytically to second order vs. first order finite differences. The author investigated the profile of Lagrange’s elliptic flight time equation and its derivatives with a special focus on their significance to the behaviour of the proposed method and the stated goal of guaranteed convergence. Possible numerical deficiencies were identified and dealt with. As a test, 28 scenarios of variable difficulty were designed to cover a wide variety of geometries. The context of this research being the orbit determination of artificial satellites and debris, the scenarios are representative of typical such objects in Low-Earth, Geostationary and Geostationary Transfer Orbits. An analysis of the computational impact of the quality of the initial guess vs. that of the order of the method was also done, providing clues for further research and optimisations (e.g. asteroids, long period comets, multi-revolution cases). The results indicate fast to very fast convergence in all test cases, they validate the numerical safeguards and also give a quantitative assessment of the importance of the initial guess.     

A robust semi-analytic method to solve the minimum ΔV two-impulse rendezvous problem

A novel semi-analytic approach is developed to determine the minimum $\mathrm{\Delta }V$ for a two-impulse rendezvous and validated both empirically and analytically. A previously published closed-form $\mathrm{\Delta }V$ estimate and the Lambert minimum energy transfer is used to establish upper and lower bounds of the minimum $\mathrm{\Delta }V$ transfer between two orbits. These bounds, in conjunction with the bisection method, operate on a nonlinear radical cost function to guarantee linear convergence. This approach has several real world applications including a low earth orbit (LEO) to highly elliptical orbit (HEO), and a HEO to retrograde geosynchronous orbit transfer. The minimum $\mathrm{\Delta }V$ estimates are better than those reported in the existing literature, while run times improved as much as two orders of magnitude over a fixed time Lambert solver. All singularity cases were addressed such that any orbital geometry, including Hohmann and radial elliptic transfers, converged to the global minimum $\mathrm{\Delta }V$. This approach will work for both coplanar and non-coplanar 3D geometries for any orbit type.     

红外热像仪在固体火箭发动机羽焰测温中的应用

利用红外热像仪测量固体火箭发动机尾焰温度场。考察了红外发射率影响测量精度的主要因素在空间的分布,通过在尾焰区域建立二维直角坐标的方法,成功实现了热图分析位置和空间实际位置之间的双向互换,在某种程度上弥补了红外热图分析软件的不足。     

发动机喷管喉衬烧蚀及热结构工程计算

给出了一种理论模型和基本计算方法。采用有限元法对某固体火箭发动机喉衬结构的烧蚀量、瞬态温度场和应力场进行了理论计算,计算结果与实测值比较一致,喉部直径实测值与计算值相对误差约为2％。     

固体火箭发动机磁悬浮试验架系统动态性能分析

在磁悬浮系统数学模型的基础上,将磁悬浮试验架对发动机的悬浮支承等效为质量一弹簧．阻尼系统。计算了磁悬浮试验架的刚度、阻尼及固有频率,计算结果得到了试验的验证。表明除系统机械设计参数外,控制方式及控制器参数是影响系统动态特性的关键因素,并通过仿真说明了控制器参数对系统动态性能的影响。为设计适合发动机试验的磁悬浮系统提供参考。     

磁通门磁强计在深空探测中的应用

磁场测量是深空探测的重要任务之一,通过磁场可以遥感行星内部、研究行星演化历史、认知太阳系天体空间环境。基于法拉第电磁感应原理的磁通门磁强计,因空间适应性强、技术成熟度高、可靠性高等特点,是深空磁场测量最为常用的载荷。简要描述了磁通门磁强计的基本测量原理,探讨了地面和在轨标定的原理和实施方法,并介绍了磁强计在空间任务中的应用方式。目前,我国已经具备了星载高精度磁通门磁强计的研制能力。在不久的将来,磁通门磁强计有望在深空探测任务中发挥重要作用。     

空间太阳能电站微波能量传输验证方案设计

微波能量传输技术作为空间太阳能电站（SpaceSolarPowerStation,SSPS）的关键技术之一,目前的研究和验证工作均集中在各单项技术的突破和验证,缺乏针对SSPS系统特点的全面优化设计。文章根据SSPS的工作模式给出了全面验证空间太阳能电站微波能量传输的验证系统方案设计,对收发天线进行了一体化设计,利用了幅度近似高斯分布的发射阵列场分布设计和低反射的接收整流阵列设计,以高精度来波方向测量和高精度移相控制为波束指向控制的技术途径。对验证系统的波束收集效率进行了分析,收集效率可达94.2%,比传统均匀分布系统高出17.6%。验证系统可从系统规模缩比、波束扫描范围、发射天线口径场分布、整流天线处功率密度、反向波束控制方法等方面模拟SSPS微波能量传输工作模式,推动SSPS系统技术的发展。     

阴阳极进气量对附加场磁等离子体推力器性能的影响

针对附加场磁等离子体推力器阴阳极供气量对推力器性能的影响进行了研究。采用试验方法测量了稳态AF MPDT的束流参数以及放电电压、放电电流、附加磁感应强度等工况参数,测量计算了推力、效率等性能参数。通过改变阴阳极供气比例,分析研究了初始气体分布对推力器性能及物理机制的影响。研究结果显示,阴极供气比重增大时,推力器推力及效率相应增大,但增大到一定比例,推力器出现工作不稳定导致性能下降,说明阴阳极存在最佳的供气比例,表明供气比例对推力器性能具有重要影响。     

槽道内涡波流场展向涡的分布特征

通过2DPIV测量了槽道内涡波流场的瞬态速度矢量场,利用雷诺分解和旋转强度的λci准则进行了展向涡的识别;并从识别的逆时针和顺时针展向涡的尺度属性、力学属性和运动属性方面分析了展向涡的分布特征。结果表明,流场中的顺时针展向涡的数量大于逆时针展向涡的数量;在法向位置0.05~0.45范围内,逆时针展向涡的旋转强度呈抛物线形状变化,顺时针展向涡整体呈下降趋势;逆时针展向涡的平均直径和椭圆平均长轴长度随着法向位置的增加而缓慢增加,而顺时针展向涡趋于下降;逆时针展向涡的椭圆长轴倾角最大值为67.28°,顺时针展向涡的等效椭圆的离心率呈下降趋势,似圆性较好;逆时针展向涡对全局平均的〈vωz〉分量在法向位置0.15处出现最大值61%;展向涡对流场雷诺应力的影响较小;逆时针展向涡的数量密度在0~0.08法向位置范围内呈急剧上升趋势,并达到最大值0.36,顺时针展向涡数量密度整体呈倒U形分布;不同旋向的展向涡占总展向涡数量密度的比值都大于0.3;展向涡的流向对流速度整体上小于流场的平均速度,且在槽道中心区域差值较大;展向涡的法向对流速度随着法向距离增加,从负值上升到正值,在中心主流两侧分别发生喷射和扫掠2种现象。