团树断层场地地震反应的数值模拟

团树断层场地在多次地震中出现了烈度异常。为了弄清烈度异常的原因，应用显式有限元方法，数值模拟了在脉冲地震波斜入射情形下断层场地的地震反应。结果表明：在基岩山脚附近土层场地地震动放大比基岩山区平均高4倍、比远离基岩山脚土层场地平均高2倍，与历史震害分布有较好的对应。     

青藏高原东侧一次β中尺度对流系统的数值模拟

利用PSU／NCAR的高分辨率中尺度非静力数值模式MM5，模拟了2001年9月18日发生在青藏高原东侧的绵阳大暴雨过程。结果表明，高分辨率数值模式对触发本次降水过程的β中尺度对流系统具有较好的模拟能力。高分辨率模拟输出显示“9．18”绵阳大暴雨与边界层内一个中尺度辐合扰动的发展和移动相伴，且中尺度辐合扰动还诱发了一个时间尺度约为6h的β中尺度涡旋和强烈发展的降水雨团。模拟还显示受高原地形影响，该β中尺度对流系统具有独特的动力和热力结构。其动力特征是强上升运动和超强散度柱与对流层低层强涡度互耦发展，热力特征是对流层低层具有对流不稳定能量，中层具有斜压不稳定能量。中尺度对流系统具有两支上升人流和两支下沉出流，低层人流(东南气流)触发低层对流不稳定能量释放，中层人流(高原近地层偏西暖湿气流)触发中层大气斜压不稳定能量释放，两种不稳定能量共同作用促使对流雨团强烈发展，形成绵阳大暴雨。     

磁场重联中离子轨道的混合模拟研究

在使用2.5维混合模拟方法研究了Petschek模型稳态驱动磁场 重联演化的基础上， 本文考察了计算域内各典型区域中粒子分布函数的变化，描绘了重联区不同位置几种类型的 非Maxwell分布函数. 结果表明，磁场重联会将重联区少部分粒子加速到很高的能量，不同 加速程度的粒子将形成球壳状的速度分布. 粒子的轨道特征表明，在重联区中出流的粒子， 有一部分被磁镜捕获，其回旋半径大于重联区宽度，并构成整个流体速度的低速部分. 另外 ，在X中性点附近进入重联区的粒子沿磁力线向出流区以三种形式漂移，分别为：沿磁力线 逃逸、捕获在磁镜中随流体运动、横越磁力线漂移，其比例分别约为70%，20%和10%。     

黄海、东海海域出海气旋发展过程中尺度数值模拟

利用MM5中尺度模式对1999年6月两个出海气旋发展过程进行 数值模拟. 数值模拟的气旋出海后移动路径与实际情况基本一致. 在数值模拟基础上重点讨 论了出海气旋发展过程潜热通量和感热通量的分布及其演变情况. 气旋出海后在气旋中心区 南方和东方存在负潜热通量和感热通量区. 出海气旋的东移和发展，其前方强大正热通量区 的存在可能是重要原因之一.     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

热层风场的理论模拟与分析

通过求解中性大气Navier Stokes动量方程建立了一个时变的三维风场理论模式，利用目前新版的中性大气模式NRLMSISE 00及国际电离层参考模式IRI2000作为输入参数给出热层风场. 基于该模式，计算得到中等太阳活动年磁静日风场的变化形态及其受电场和离子曳力的影响. 同时，将Navier Stokes动量方程作不同形式的简化，并利用简化模式与本文的模式计算结果的对比，分析中性大气Navier Stokes动量方程中黏性项以及非线性项(U·Δ)U的作用. 结果表明，本文所建立伪三维风场模式给出的结果更为合理，而简化模式在某些地区尤其在低纬和赤道区不适用，黏性项及非线性项的作用不可忽略. 本文所建立的风场模式将对研究电离层动力学过程、电离层与热层的耦合过程以及空间天气学研究都有着重要意义.     

高斯平面坐标换算到测区平均高程面上的方法

依据高斯投影原理，提出了一种在高斯平面与任决高程平面间进行坐标转换的方法，并设计了相应的计算软件。     

空间直角坐标的协因数阵转换到高斯平面上的计算公式

导出了由空间直角坐标的协因数阵向大地坐标的协因数阵；由大地坐标的协因数阵向高斯平面直角坐标的协因数阵转换的计算公式。在此基础上，导出了直接由空间角坐标的协因数阵计算高斯平面直角坐标的因数阵的简便公式。     

高斯投影正反算的递推算法

本文通过引入等量坐标概念，讨论高斯坐标和等量坐标以及大地坐标之间的关系，推导高斯投影坐标正反算的递推算法公式。