用Ghost Fluid方法耦合流体和固体计算

在数值模拟中，欧拉方法适合计算大变形问题，但是物质界面的准确确定比较困难；拉氏方法计算采用物质网格，物质界面准确，但在大变形情况下，计算难以进行。在武器物理中爆炸产物与金属相互作用问题中，既存在多物质问题又存在大变形问题：爆炸产物和气体变形较大，金属相对变形较小。对于此类问题可以采用耦合拉氏欧拉方法(CEL)：流场大变形区域采用欧拉方法计算，在小变形区域采用拉氏方法计算，充分利用拉氏和欧拉方法各自的优点。

三维修正相容拉氏格式保球对称性质的证明和应用

计算格式的对称性保持性质是计算方法研究的一个重要内容。欧拉计算方法一般通过计算物理空间的坐标系选取来保持计算格式的对称性。如计算球对称性问题，一般采用球坐标系。在三维拉氏算法中，主要采用笛卡儿坐标系，球面放射状结构网格是保球对称性的必要条件。没有球对称性网格，严格的球对称性保持就失去了基础。但即使在球对称性网格的情况下，三维拉氏计算也不能保证一维球对称性。

三维多介质欧拉有限体积格式的混合算法

对于多介质欧拉方法，混合网格物理量的计算是其难点和关键点之一。这里提出的方法是运用Yonugs界面重构技术确定出混合网格内物质的界面，界面确定后，混合网格内每一部分可能是非规则的四面体、五面体、六面体或七面体，采用对非规则区域适应性很强的有限体积法对每一部分分别进行计算。这种方法虽然比较复杂，但是它兼有拉氏方法的优点，因此计算出的混合网格内每一部分物质的物理量比较精确。

拉氏流体动力学相容计算方法程序CHAP及其应用

对于所有的数值模拟方法而言，其基本宗旨是离散后的控制方程尽可能地忠实于原微分形式控制方程，准确地反映微分形式动力学控制方程所显含及隐含的基本运动规律。然而，在只考虑截断误差精度的离散过程中，往往失去了原微分形式方程所含有的重要物理性质，或者只是近似满足，这种近似性在某些情况下往往变得很差。例如，在能量方程的离散中，现有方法往往不考虑动能、内能和总能之间的严格物理关系，造成的结果是，对内能方程离散不能严格保证总能守恒。