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

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

非对称射流形成的几何理论(二)

将平板非对称正碰撞形成射流满足的定常方程组封闭条件推广到非正碰情况，提出了平面二维条件下非对称射流形成的一般理论，给出了射流形成的临界条件和出流射流的理论公式。对斜碰撞情形下非对称射流的出流方向、厚度随碰撞角和来流初始位形的变化规律进行了理论预测和数值模拟，理论解与二维欧拉程序的仿真结果符合较好。     

强动载作用下孔洞汇合对延性金属层裂损伤演化过程的影响

针对强动载作用下延性金属的层裂问题，在分析孔洞之间几何关联的基础上，定义了一个新的耦合损伤及孔洞几何信息的孔洞汇合判定方法，同时，基于能量守恒原理，解析了孔洞汇合对损伤快速增长影响的物理机理.通过分析数值计算结果和对比相关文献的实验可知:孔洞汇合后不仅引起损伤增长，而且导致了损伤材料内部微孔洞数目的减少、孔洞平均尺寸的增加。     

聚能装药侵彻混凝土靶板的数值模拟

将MEPH2Y程序和LTZ-2D程序衔接在一起,应用于分析射流、射弹对混凝土靶板的侵彻问题.模拟射流对混凝土靶板的侵彻以及靶板的损伤破坏情况,计算结果与实验结果符合得较好.     

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

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

多组份计算中封闭性模型研究

在使用单速度多组份方法计算多介质混合网格中的物理量时,需要给出一种封闭模型来使控制方程封闭。本文分析了压力增量相等封闭性模型存在基础遭到破坏的原因,之后在保证流场守恒性的前提下构造可以保证假设基础的可解方程组对压力增量相等封闭性模型进行修正,给出了一种改进方法。使用改进的模型进行了爆轰驱动飞片问题的计算,计算结果表明,改进的模型能很好的追踪物质界面和冲击波位置,同时较好的抑制界面处压力、内能震荡。     

三角波强加载下延性金属多次层裂破坏问题

基于应力瞬时断裂判据和Tuler-Butcher损伤累积判据，分析了没有升压的简单三角波强加载下延性金属的多次层裂破坏问题。分析结果显示:层裂片的厚度随着冲击波宽度与强度的比值的增大而增大，材料的破坏深度小于冲击波宽度的一半；在强冲击加载下，材料的破坏深度约为冲击波宽度的一半。     

缺陷形状对铝材料微喷射的影响

通过二维弹塑性流体力学欧拉程序MEPH,对包含不同形状不同角度缺陷铝材料的微射流现象进行数值模拟,并对其最大喷射速度以及喷射系数进行数值统计,得到不同形状、不同角度缺陷最大喷射速度及喷射系数的变化规律.并将计算结果与实验结果进行比较,基本吻合,能够反映试验规律.     

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

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