保质量的自适应坐标变换方法

使用Lagrange方法求解流体力学问题时，常常会遇到网格扭曲变形的现象。这种扭曲常由流场（也是网格运动场）的旋度而引起。当旋度的计算不能保持精确时，就会引起非物理的网格扭曲，并导致计算失准甚至中止。     

近似保角的自适应坐标变换方法

自适应坐标变换方法是为解决多介质和大变形问题而提出的一类网格生成方法，该方法中的一种为近似保持网格夹角不变，保持物质界面为拉氏描述，并要求网格速度在最小二乘意义下尽量靠近流体运动速度。这里所讨论的坐标变换的自适应性，指的是新坐标系自动适应流体流场的一些重要特性（接近流体速度）以及保持网格的几何特性（保角）。为了处理多介质情况，网格方程应在子区域的所有边界上给出边界条件。     

RESEARCH ANNOUNCEMENTS High order Zumerical Solution of Integral Transport Equation in slab Geometry

There are some common numerical methods for solving neutron transport equation, which including the well-known discrete ordinates method, PN approximation and integral transport methods[1]. There exists certain singularities in the solution of transport equation near the boundary and interface[2]. It gives rise to the difficulty in the construction of high order accurate numerical methods. The numerical solution obtained by now can not attain the second order convergent accuracy[3,4].     

二维浅水波方程的数值激波不稳定性

研究二维浅水波方程的数值激波不稳定性问题.线性稳定性分析和数值实验表明,格式的临界稳定性与数值激波的不稳定现象有重要的联系.基于扰动量的增长矩阵分析,本文将高分辨率的数值格式和HLL格式进行特定的加权,设计一类新的混合型数值格式.其中可以调节非线性波速的HLLC与HLL的混合格式,数值试验展示了消除浅水波方程激波不稳定现象的有效性和鲁棒性.     

自适应坐标变换方法

提出求解流体力学Euler方程的自适应坐标变换方法,它们保持网格夹角近似不变、保持物质界面为Lagrange描述,且使得网格速度与流体速度的差在最小二乘意义下达到极小.所提出的新坐标系能自适应于流体流场的若干重要特性.     

格子Boltzmann方法求解Burgers方程

众所周知,格子方法(包括格子气和格子Boltzmann方法)在计算物理领域取得巨大进展。与之形成鲜明对比,格子方法的数学理论始终处于停滞不前的状况。为求解Burgem方程。一类带有BGK模型格子方法被构造出来,经过变量替换,发现他们属于三层非线性差分方法。使用极值原理,给出此类格式稳定性的严格证明。最后,从数值实验中可以看出,使用LBM得到的结果,与经典二阶守恒差分方法的结果符合得非常好。     

格子Boltzmann方法求解Burgers方程

众所周知,格子方法（包括格子气和格子Ｂｏｌｔｚｍａｎｎ方法）在计算物理领域取得巨大进展。与之形成鲜明对比,格子方法的数学理论始终处于停滞前的状况。为求解Ｂｕｒｇｅｒｓ方程,一类带有ＢＧＫ模型格子方法被构造出来,经过变量替换,发现他们属于三层非性差分方法。使用极值原理,给出此类格式稳定性的严格证明,最后,从数值实验中可以看出,使用ＬＢＭ得到的结果,与经典二阶守恒差分方法的结果符合得非常好。     

SPH方法中的Riemann解与人工粘性

本文描述了光滑粒子动力学方法的人工粘性和Riemann解方法,分析了Godunov方法与传统人工粘性方法的耗散项．给定一种人工粘性,总可以找到一种相应的Riemann解法器,使得它们的耗散项在形式上几乎相同．本文利用各种近似Riemann解构造了相应的新的人工粘性．这些新的粘性,无需人工调节粘性系数．本文完成了多个数值试验,比较了使用传统人工粘性方法与Riemann解方法的不同．采用新的人工粘性,辅助热通量粘性,可以获得令人满意的计算结果．     

Burgers方程的声速点故障问题

研究了在Burgers方程跨声速稀疏波计算中遇到的sonic point glitch问题,对它产生的原因及其与数值格式熵条件的关系进行了分析.对若干著名格式,按照是否满足熵条件进行了分类.为了消除sonic point glitch现象,提出了一种新的两步分裂方法,并用这种方法改进了一系列典型格式.数值试验表明这是一种很好的消除sonic point glitch的方法.     

电子与高荷电离子碰撞激发的扭曲波截面

采用我们的“快速”扭曲波方法和程序MCDW(9),计算了类氖Ce~(22+)27个精细能级之间351个过程的激发截面,并且与其它理论值进行了比较.结果表明:碰撞强度大的过程一般符合在10%左右.但是由于组态混合导致跃迁矩阵元严重抵消而变成为相对小的过程以及激发能很小的光学禁戒(电偶极振子强度为零)过程,其截面(或碰撞强度)对模型细节敏感,不同作者计算差别较大. 关键词：