Gravity Field Recovery from GRACE: Unique Aspects of the High Precision Inter-Satellite Data and Analysis Methods

The very high accuracy of the Doppler and range measurements between the two low-flying and co-orbiting spacecraft of the GRACE mission, which will be at the μm/sec and ≈10 μm levels respectively, requires that special procedures be applied in the processing of these data. Parts of the existing orbit determination and gravity field parameters retrieval methods and software must be modified in order to fully benefit from the capabilities of this mission. This is being done in the following areas: (i) numerical integration of the equations of motion (summed form, accuracy of the predictor-corrector loop, Encke's formulation): (ii) special inter-satellite dynamical parameterization for very short arcs; (iii) accurate solution of large least-squares problems (normal equations vs. orthogonal decomposition of observation equations); (iv) handling the observation equations with high accuracy. Theoretical concepts and first tests of some of the newly implemented algorithms are presented. This revised version was published online in August 2006 with corrections to the Cover Date.     

求解航空发动机非线性方程组的变步长牛顿法

详细分析了航空涡轮发动机部件共同工作的非线性方程组的求解收敛性问题,并对计算中常见的迭代不收敛问题的机理进行了研究总结。针对常见的不收敛问题,提出了以传统牛顿迭代法为基础的变步长牛顿法加以解决。采用变步长牛顿法后,对某型定几何单轴涡喷发动机进行了一系列的验证计算,计算结果表明计算收敛性得到较大幅度的改善。变步长牛顿法对改善发动机总体性能计算收敛性具有很好的参考价值。     

线性拟合法建立航空发动机状态变量模型

针对非线性拟合法建立航空发动机状态变量模型过程中的初值问题,提出了一种能够避免该问题的线性拟合法,即在离散域情况下,利用已知的状态序列来构建系统矩阵与模型响应之间呈线性关系的目标函数,从而使非线性最小二乘问题转化为线性最小二乘问题.应用该方法建立某型涡扇发动机的小偏差状态变量模型,并与非线性拟合法进行了对比研究,结果表...     

基于Riemann问题的解在自适应无结构网格上求解Euler方程

给出了一种根据 Riemann问题的解 ,计算网格单元边界处的守恒量通量的方法。该方法不受网格单元形状的限制 ,具有较好的通用性。实际计算表明 ,采用这一方法在自适应无结构网格上求解 Euler方程 ,对于复杂的流动细节具有较高的分辨率 ,并且使用多重网格法能够加速收敛     

分拆大偶数为两个几乎相等的素数之和

我们称n =p+q 为偶数n的最平衡哥德巴赫分拆 ,其中 p与 q均为奇素数 ,p≤ q且 |p - q|取可能的最小值。根据哥德巴赫猜想 ,对于每个大于或等于 6的偶数 ,必存在其最平衡哥德巴赫分拆。本文通过对 10 7以内的偶数进行的数值计算 ,进一步提出如下猜想 :由所有 (≥ 6的 )偶数之最平衡哥德巴赫分拆产生的三组比值 p/q ,p/n 以及 (n/2 -p) /n的均值 ,分别具有极限 1,1/2及 0。     

定常幅值小推力登月飞行器轨道研究

 进行了基于平面三体模型的登月飞行器轨道控制方法的研究;研究了从近地低轨道到近月低轨道的飞行轨道;给出了在地球逃逸段、惯性漂移段和月球捕获段的运动轨迹和关键点的参数。提出使用“远地点可达”概念完成了地球逃逸段发动机推力终点的选择和使用飞行器相对月心能量完成了在月球捕获段止推发动机工作初始点的选择。     

基于圆型限制性三体问题模型的行星卫星逃逸能量研究

 基于由飞行器、行星及其卫星组成圆型限制性三体问题模型,通过庞加莱映射的方法,研究了飞行器从行星卫星附近逃逸的问题。在Jacobi常数确定的前提下,通过逆向积分,飞行器从L1或L2点附近返回近月点,得到近月点速度出发速度。研究结果表明绕飞L1点和L2点逃逸行星卫星需要的最低能量是不同的,从月球表面逃逸所需的速度脉冲分别比开普勒算法节省46.5m/s和42.3m/s,且均小于Villac等人在Hill模型下得到38.9m/s,从而改进了Villac等人的相关工作,同时也给出了从太阳系主要行星卫星表面逃逸所需的最小能量。     

新型高低速流统一算法及其叶轮机应用

 对于可压缩/不可压缩混合流动或高低速混合流动问题,目前没有较好的解决方法。基于Roe格式及预处理方法,提出了全速度Roe格式。全速度Roe格式是将Roe格式与低速Roe格式通过一个与本地马赫数相关的函数联合而得到的。与标准Roe格式相比,对亚音速流动,全速度Roe格式中声速对数值耗散的影响随着马赫数变化而变化。与预处理Roe格式相比,全速度Roe格式简单易于实现,而且没有计算稳定性的限制,不需要采用全局截断处理,从而在理论上对高低速流动具有一致精度。事实上,全速度Roe格式意味着可以通过直接改变Roe格式的非线性特征值,解决可压不可压流动统一计算精度的问题。这是对统一算法的新观点。数值模拟表明,全速度Roe格式具有与预处理Roe格式一样的激波捕获能力、极低速流动准确模拟能力。     

钉载作用下多裂纹板的裂纹闭合接触问题

 基于摩擦接触问题的数学规划解法,采用各向异性体平面弹性理论中的复势方法,建立了含裂纹群有限大各向异性板,在钉载作用下裂纹闭合或局部闭合问题的有效分析方法。忽略钉与孔间的摩擦,假设钉载沿孔边呈余弦分布,通过在可能闭合的裂纹边界引入互补变量函数并将其展成Fourier级数形式,以Faber级数为工具,应用保角映射技术和最小二乘边界配点法,导出裂纹面摩擦接触的线性互补模型,并通过算例验证了方法的有效性。数值结果表明,由于采用级数解描述板应力场和位移场,本方法具有较高的计算精度和效率,便于研究裂纹闭合对应力强度因子等断裂参数的影响。     

关于某个四阶非对称微分方程边值问题的唯一性

本文处理如下形式的四阶非对称微分方程边值问题的唯一性。 对于对称情形,传统的处理方法是利用Thichmarsh,E.C.的函数论方法与亏指标理论。而当方程为非对称时,上述方法会有些困难。但若通过构造某个特殊的与(?)λ有关的内积,并利用内插不等式及Opial不等式,则非对称情形能被容易地处理且得到许多唯一性的结果。