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

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

非线性抛物型方程计算方法 值此周毓麟先生九十寿辰之际, 谨以本文表达对他的崇高敬意和真诚祝福

本文简要回顾非线性抛物型方程差分方法若干研究工作, 包括周毓麟先生在该研究方向取得的部分研究成果, 并对近年来相关的部分研究进展进行综述, 展望拟开展的研究工作.     

格子Boltzmann方法求解Burgers方程

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

格子Boltzmann方法求解Burgers方程

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

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

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

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描述,且使得网格速度与流体速度的差在最小二乘意义下达到极小.所提出的新坐标系能自适应于流体流场的若干重要特性.     

辐射扩散方程的自适应坐标变换方法

借鉴自适应坐标变换流体力学方法和其他一些有益工作，研究辐射扩散方程的自适应坐标变换方法，并考虑辐射扩散和流体力学耦合的辐射流体力学问题。对辐射扩散方程和一般二维对流扩散方程，从积分和微分形式出发，导出自适应坐标变换下的等价形式，对数值离散格式进行初步讨论。对二维三温辐射扩散方程，采用时空有限体积方法，进行数值离散；利用局部函数近似离散扩散算子。利用合理的进程分裂，简化耦合问题的计算。