重力插值方法研究

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

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

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

露天采石场边坡稳定性与加固

舟山群岛位于我国南北海运和长江内河航运干线的交汇点，近临上海、宁波等大都市，区域条件独特。区域内岛屿众多，深水港资源集中，渔业、旅游资源丰富，岛内以丘陵山地为主，地少人多。随着国民经济快速发展，船舶修造业建设、劈山造地、围填海涂及各类市政交通建设等日益增多。使舟山放炮开山的露天采石场仍然较多。     

市场经济下土地集约利用动因探析

土地是人类最宝贵的资源，土地资源的稀缺性、不可替代性和超常稳定性以及我国国情要求我们集约和合理利用每寸土地，摒弃过去那些粗放型利用方式，更新“花园式工厂”类的土地利用观念。当前土地集约利用情况并不乐观，产生的原因有内在的、外在的、体制上的、管理上的、观念上的等等。较为复杂，需要层层剖析，找出土地权利、制度之间的内在联系，究其根源加以解决。     

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

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

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

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

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

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

液核地球自转速率变化的带谐潮效应Ⅱ——数值计算与比较

地球自转速率的潮汐变化可由无量纲参数k／Cm)(k和Cm分别为壳幔的有效Love数和有效板转动惯量)来表征。对于一个具有弹性地幔、平衡海潮和核幔不耦合的地球k／Cm=0．944，且与潮波频率无关。海潮的非平衡扰动使k／Cm为复数，且与频率有关。大气对自转速率有效勒夫数的贡献约为△kat=0．0075。同时地幔滞弹性对勒夫数也产生扰动。利用本文得到的理论公式和最新的潮汐数据计算了地球自转速率的潮汐变化，及其有关地球物理机制的影响。