弹塑性流体力学拉氏自适应网格和网格重分守恒重映数值模拟方法

拉氏自适应重分弹塑性流体力学有限元程序实现了网格完全自适应，具有良好、灵活的非结构自适应网格数据结构，实现了滑移界面两边(接触间断)网格动态调整，网格的细分和合并处理灵活，网格重分和网格自适应模块兼容、守恒重映，网格重分中采用多种方法控制新网格的质量，爆轰计算可采用Lee-Tarver的化学反应率模式。初步数值计算结果表明，弹塑性流体力学拉氏自适应重分数值模拟方法合理，计算结果正确，基本反映了流场的物理结构。     

强光光学元件波前检测算法

针对强光光学元件高精度检测中使用的两个新参数一——波前位相梯度和功率谱密度的数值计算方法进行了研究，并对目前使用的大口径干涉仪传递函数进行了实验标定。     

不同倾角的脉管制冷机内自然对流影响的研究

本文通过对脉管内自然对流的数值模拟，讨论了当脉管冷端处在不同方向时，脉管内的自然对流对脉管制冷机性能的影响．数值模拟结果表明，当冷端在下时获得的制冷温度最低；当热端在下并与重力方向成３０°时制冷温度最高．结论与实验测量结果定性一致．所做的理论分析对实现脉管制冷机的良好运行有指导意义。     

双组元姿控发动机反流特性的DSMC模拟

为准确获得双组元姿控发动机的羽流及其反流特性,采用直接模拟蒙特卡罗(DSMC)方法结合可变软球(VSS)模型数值模拟以N2O4/MMH为推进剂的真空羽流。喷口参数由喷管流动的DSMC法计算得到。与文献比较表明本文计算程序可信。模拟结果表明,分子质量越小的气体越容易进入反流区;本文的计算条件下,在偏离X轴正向116°~140°范围内,组分数密度和质量流率以及温度均出现了最大值。     

在大学物理教材中引入计算机数值解的尝试

本文从普通物理教学现代化的要求出发,通过新编《普通物理学》中若干例题的求解,论述了在大学物理教学中适当引入计算机数值模拟的研究方法将有助于物理概念的深化,以及学生综合能力的培养.     

注重培养研究能力的《微分方程数值解法》课程教学研究与实践

论述了注重培养研究能力的教学的重要性,数学专业基础课教学在研究能力培养中的作用,在《微分方程数值解法》课程教学中的体会,以及关于《微分方程数值解法》课程教学进一步改进的探讨.     

Electrical Potential Distribution around a Charged Colloidal Particle: Nonlinear Integral Equation

多电解质溶液中带电胶体粒子的电势分布由球形Ｐｏｉｓｓｏｎ－Ｂｏｌｔｚｍａｎｎ方程（ＰＢＥ）描述．ＰＢＥ是一个非线性的微分方程，且难以求得其解析解．本文采用非线性Ｐ－Ｂ积分方程，计算电势分布的数值解．首先，根据静电场和热力学系统中的物理定理，导出描述电势分布的Ｐ－Ｂ积分方程（ＰＢＩＥ）；其次，用迭代方法求ＰＢＩＥ的数值解．最后，计算了在３－１型电解液中无量纲半径κａ分别为０．１２和０．２２，无量纲表面电势ξ分别为１，２，４，６时球形胶体粒子外部的电势值．为了检验数值解的精度，计算了表面电荷密度，并与Ｌｏｅｂ（１９６１）和Ｏｓｈｉｍａ（１９９５）等人的结果比较，本文结果的相对误差小于１％，优于Ｏｓｈｉｍａ的结果．     

分子阻障内旋转的配分函数

本文用WIGNER-KIRKWOOD微扰展开法来推导分子阻障内旋转配分函数的近似表示式,准确到h[2]项,所得公式只包含两个参量(即分子自由旋转配分函数Qf和势垒高度UO与kT之比值)并可用BESSEL函数来表示,数值计算结果表明,在内旋转转动惯量较大,温度较高这一实用上最为重要的范围内,根据本文中所得公式,对任意的Qf值和U0/kT值,可直接算得阻障势能所产生的熵的修正值(Sir-Srir)/k,比用数值求解法较为简单而准确,如计及微扰展开的高次项,还可以在整个温度范围内取得良好的数值结果。不仅如此,本文所得出的解析公式容易推广到多维运动的情况中去,因此较PITZER的数值求解法具有更普遍的意义。