任务轨道设计

本文首先说明太空任务与轨道设计的关系 ,接着介绍轨道的基本性质。从地球重力势的观点看各种常用的绕地轨道 ,包括地球和太阳同步轨道及Molniya轨道。从扰动的观点看常用的星际轨道 ,包括LISA、ASTROD、SOHO轨道。最后对星际轨道设计 ,说明二点边界值问题的数值解法、飞掠星体的应用、最佳化的考虑 ,并用以设计 2 0 1 5年发射的ASTROD初步任务轨道。     

富氏谱分析法模拟地壳垂直运动的地学依据及其策应力模型

应用地壳波浪与镶嵌构造学说对富氏谱分析法提取地壳垂直形变信息的科学性做了地质学意义上的阐释 ,并提出了根据多期形变资料提取特定波段上构造策应力的数学模型     

绝热非径向太阳振动的数值计算方法

本采用改进后的平行Ｓｈｏｏｔｉｎｇ方法，即连续正交归一化方法，通过增加系统的秩和把微分方程从线性变成非线性的方法，使得病态方程变成非病态方程，克服了平行Ｓｈｏｏｔｉｎｇ方法（不连续正交归一化方法）的主要困难：一是选择径向网格点数目的问题，如果被检验振动模的径向节点数大或小，所需要的网格点数则多或少；二是在远离阈值的时候，尚未有合适的方法提高计算精度的问题。数值计算出了标准太阳模型的绝热非径向太阳     

褶皱岩层初始挠曲的分叉特性

基于拖带坐标法建立了岩层大变形的一维非线性微分方程，然后用解析与数值方法研究了变刚度、可伸缩岩层在侧向压力及自重作用下的分叉特性，初步揭示了褶皱岩层初始挠曲的形成机理。     

重力插值方法研究

为了根据离散观测数据构制连续空间重力变化图像，分析和讨论了３种数值插值方法，计算结果表明多面函数方法插值精度最高。由于逐步回归分析筛选核函数中心点的计算繁琐，文中提出根据分形理论和Ｓｈａｎｎｏｎ取样定理来确定核函数中心点。对滇西试验场进行模拟试算，插值精度可达到４～５（１０－８ｍｓ－２）。     

关于非球形粒子光散射的T-矩阵数值计算方法的研究

首先介绍了无规取向轴对称非球形粒子光散射的T－矩阵数值方法。然后，进一步研究了T－矩阵计算的收敛问题，提出了一种新的关于无规取向无吸收非线形散射物体的物理收敛方法，我们研究了椭球粒子的收敛问题，结果表明我们的数值方法和T－矩阵计算程序是有效的， 收敛精度与粒子的尺度和形状有很强的依赖性，在一定的条件下，我们的物理收敛速度优于NASA Mishchenko的数学收敛速序。     

三类密近双星自转同步性参量的理论计算和统计分析

利用双星自转同步性理论给出了６９个三类密近双星系统中９３个子星的临界同步自转参量和临界自转周期．并把利用临界自转同步参量所计算的临界自转周期与由气体星自转不稳定理论所计算的临界自转周期做了比较，其结果是两者均属同一量级．     

在重力分层的太阳大气中磁流浮现的数值解

为了解释太阳观测中磁流浮现形成活动区的主要特征——AFS弧顶的上升运动和弧脚的下沉运动,谱斑,Ellerman炸弹等,我们采用了MHD-Lagrangian数值模拟法计算了重力分层(压力密度指数递减),太阳大气中同基底磁场反向和同向的磁偶极子从光球下浮出时产生的动力学演化过程.计算证明:反向的偶极子会导致流体产生会聚和下降运动,为了能让反向偶极子浮出必须适当增压,演化的结果形成电流片,因而有可能引起如Ellerman炸弹等小型爆发现象.还证明了:在浮现处一定产生加热增温,即导致谱斑发亮现象.在弧顶气体是上升运动,速度小;弧脚处气体作下降运动,速度较大,这主要是重力分层气体压力密度非均匀引起的,和Brunt-Vaisala振动有同一的起因.我们还模拟了通常活动区形成时的初期现象:1500高斯强磁偶极子浮进100高斯的弱光球色球基底场时产生的动力学现象.计算表明:为了不致产生激波必须减压浮出,即类似于黑子降温式的浮现,浮现时间应是小时的量级,浮现后确实在色球里导致强的电流片,当这电流片升至低日冕后,会成为Heyvaerts,Priest等人期望的耀斑发生源,因而浮现过程有可能是耀斑储能和触发机理之一.另外我们利用计算结果试解释了北京天文台怀柔站的磁流浮现资料,即文[1]和文[2]的一些分析结果。     

介质横向不均匀对震后地表位移场的影响

利用三维有限元模拟两层地球模型(包括弹性层和粘弹性层)介质的横向不均匀和断层错动的复杂性对震后地表位移场的影响，其数值结果包括弹性解和粘弹性解。结果表明，震后位移场的分布特征不仅与粘弹性层的流变特性有关，还与弹性层介质的剪切模量以及断层的同震错动有关。