欧拉流体力学数值方法中的界面计算

本文所提算法适用于二维和三维多介质流体力学两步欧拉数值方法中输运计算的混合网格(包括自由面网格)界面处理。在一个混合网格中，界面被近似地看作直线(二维)或平面(三维)。整个方法分为三步：(1)第一步，用混合网格周围的八个网格的介质面积份额(二维)或二十六个网格的介质体积份额(三维)确定界面的法线方向；第二步，用混合网格的本身的介质面积份额(二维)或体积份额(三维)确定界面的方程(位置)；第三步，用此直线方程求出通过网格边界的流以及下一时刻网格的面积份额(二维)或体积份额(三维)。最后给出了用此方法所做的一些数值计算及与SLIC算法的比较。     

二维多介质欧拉自适应网格方法研究

要对聚能和流场的剧烈变化进行模拟，欧拉方法具有明显的优势。而在这些方面的研究中，所涉及的流场十分复杂，为达到所需的计算精度，必须采用很密的网格才能以较好的分辨率去模拟流场的剧烈变化部分和介质的界面，特别是大空间尺度流场局部细节的数值模拟，有些问题如果用统一网格计算，即使最快的计算机也不能提供足够的分辨率。所以目前计算机的内存和速度限制了整个计算区域的网格细分程度，对计算区域作局部的网格自适应细分、用大的动态空间分辨率划分流场是必要的和迫切的。     

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

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

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

在采用Youngs界面重构技术的基础上，对三维欧拉方法混合网格的计算格式进行了研究。运用Youngs技术确定界面后，混合网格内每一部分物质一般不再是正规的六面体结构，可能是非规则的四面体、五面体、六面体或七面体。本文采用对非规则网格适应性很强的有限体积法对每一部分分别进行计算，给出了混合网格内每种物质的压力、速度、能量等的计算公式，比较有效地解决了混合网格的计算问题。     

矩形网格上二阶精度流体体积界面追踪方法

欧拉方法模拟多流体问题的难点在于混合网格中界面的处理及物理量的确定。因此，界面处理技术是欧拉方法的主要研究内容之一。     

二维高速碰撞问题欧拉数值模拟的混合网格计算

提出了适用于二维平面或轴对称多介质流体力学两步欧拉数值方法中输运计算的混合网格界面处理.在一个混合网格中,将界面近似看作直线.整个方法分为3步:①用混合网格周围的8个网格的介质面积份额确定界面的法线方向;②用混合网格的介质面积份额或体积份额确定界面的直线方程;③用此直线方程求出通过网格边界的流.给出了用此方法所做的测试、数值计算及与其它算法的比较.     

多介质流体力学两步欧拉方法的模型封闭性方法

描述多介质流体力学两步欧拉方法拉氏步中方程的不封闭问题.分析已存在的模型封闭性方法的优缺点,给出根据声速特征的自适应的模型封闭性方法.数值例子表明方法合理可靠,并说明压力弛豫过程在数值计算中的作用.     

金属圆筒内爆轰波对碰效应的欧拉数值模拟

运用基于Youngs技术的多介质欧拉方法及应用该方法研制的二维多介质弹塑性流体动力学欧拉程序meph2d对金属圆筒内爆轰波相互作用过程进行了数值模拟并获得成功，计算图像与实验照片吻合，计算结果与实验测量结果符合较好。说明采用多介质欧拉方法可以对此类爆轰波相互作用及爆轰驱动问题进行有效的数值模拟。     

高速碰撞的三维欧拉数值模拟方法

介绍一种可用于高速碰撞问题数值模拟的三维弹塑性流体动力学数值计算方法,以及应用该方法研制的应用软件ＭＥＰＨ３Ｄ,并给出了一些高速碰撞问题的数值计算结果,计算结果表明该程序具有计算三维高速碰撞问题的能力。     

混凝土动态冲击问题的一种欧拉数值方法

针对混凝土的动态冲击问题给出了混凝土的Holmquist Johnson Cook本构关系 ,本构关系中的等效强度作为压力、应变率和损伤的函数 ;压力表示为体积应变的函数 ,并考虑了永久破碎的效应 ;积累损伤作为等效塑性应变、塑性体积应变和压力的函数。结合三维弹塑性流体力学欧拉数值方法 ,提出了混凝土的Holmquist Johnson Cook本构关系与欧拉程序相结合的计算方法。介绍了带混凝土Holmquist Johnson Cook本构关系的三维弹塑性流体力学欧拉程序。