分子动力学模拟的主要技术

综述了分子动力学模拟技术的发展，介绍了分子动力学的基本 原理、有关的有限差分技术、势函数的发展、初始条件和边界条件的 选取、平衡态系综及其调控、感兴趣量的提取、分子动力学的特别用 途以及与其他计算方法的结合. 最后还指出了分子动力学模拟方法本 身进一步研究的方向.  

摄动有限差分方法研究进展

振动有限差分（ＰＦＤ）方法，既离散徽商项也离散非微商项（包括微商系数），在微商用直接差分近似的前提下提高差分格式的精度和分辨率．ＰＦＤ方法包括局部线化微分方程的摄动精确数值解（ＰＥＮＳ）方法和摄动数值解（ＰＮＳ）方法以及考虑非线性近似的摄动高精度差分（ＰＨＤ）方法。论述了这些方法的基本思想、具体技巧、若干方程（对流扩散方程、对流扩散反应方程、双曲方程、抛物方程和ＫｄＶ方程）的ＰＥＮＳ、ＰＮＳ和ＰＨＤ格式，它们的性质及数值实验．并与有关的数值方法作了必要的比较．最后提出值得进一步研究的一些课题．  

爆炸荷载作用下影响RC梁破坏形态的主要因素分析

在冲击、爆炸等强荷载作用下，钢筋混凝土结构通常是发生弯曲破坏。室内外试验结果表明，当动荷载峰值较大、作用时间较短（脉冲荷载）时，钢筋混凝土结构易发生脆性的剪切破坏，其主要原因是脉冲荷载的高频成份丰富和加载速率高等特点容易激发结构构件中的剪切变形和剪 应力。本文基于Timoshenko梁理论，提出了一种爆炸荷载作用下钢筋混凝土梁破坏形态的有限差分预报方法，利用该方法分析了荷载加载速率、截面高度、混凝土强度、配筋率等对钢筋混凝土破坏形态的影响规律，并指出，随着荷载升压段加载速率、截面高度、主筋配筋率的增加或混凝土强度的降低，梁的破坏形态从弯曲破坏转变为剪切破坏。  

迷宫密封中流场的有限差分模拟

本文用k-ε紊流模型和可压缩流体压力修正方程以及SIMPLEC算法对时均形式的二维N-S方程进行了求解,得到了单腔室迷宫密封中的流场特性,其结果可以用来改进迷宫密封转子动力特性的分析模型.  

多孔介质中的非达西自然对流的分岔研究

利用分岔理论研究了多孔介质底部加热所引起的非达西自然对流。用有限差分方法计算了对流的分岔；确定了Beta数与临界瑞利数的关系。结果表明：随着Be从0增大到1，出现分岔的单胞对流的临界瑞利数Rac从39．35单调地增大到41．15。双胞对流亦有类似的趋势。这说明惯性－湍流效应有使对流稳定性增强的趋势。  

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。  

变温场中具损伤粘弹性矩形板的非线性动力响应分析

基于热粘弹性理论、Von Karman板理论和连续损伤力学,导出了二维状态下各向同性材料的变温粘弹性本构方程,建立了含损伤效应的各向同性粘弹性矩形板在变温场中的非线性运动控制方程,且应用有限差分法对问题进行求解.算例中,讨论了损伤演化及温度场等因素对粘弹性矩形板非线性动力学行为的影响,得出一些有意义的结论.  

非线性动力学常微分方程组高精度数值积分方法

建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．  

聚能射流三维数值模拟

研究三维多介质界面处理及数值模拟问题。采用拉格朗日法、欧拉法相结合的方法在矩形网格上离散差分基本方程组；在欧拉步中引入模糊方法处理界面，计算各输运量；编写了数值模拟程序，并对线型装药金属罩聚能射流模型进行模拟。证明模糊界面描述和模糊输运计算有效、可行。