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

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

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

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

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

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

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

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

炸药爆轰驱动超高速击波管的数值摸拟

 高能炸药（HE）柱壳沿轴向传播的稳定态爆轰，在置于炸药柱壳内部的泡沫塑料（工作流体）中产生很强的轴向击波，达到极高的压力和能量密度。适当设计，可使工作流体的流动接近一稳定态，与炸药爆轰以相同的相速度前进。击波阵面后的工作流体可用Nozzle流动方程很好地描述^[1]，Nozzle收缩段中的流动不受散开段流动的影响。利用这一性质，可以设计出达到极高压力和能量密度的超高速击波管，并且不受管壁物质强度的限制。数值模拟计算给出了上述稳定态物理图象，并显示出开始阶段轴向击波的形成过程及其后对稳定态的逼近，计算结果及物理图象与理论分析符合得很好。     

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

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

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

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

局部采用标志质点的流体网格法

本文提出了一个解具有大形变的二维非定常的多物质流体力学问题的数值方法,在采用流体网格法思想的同时采用了标志质点,在不同的物质的界面附近,借助于标志质点计算流量;对连续流的区域,用流体网格法计算流量。计算检验表明,与质点网格法相比,减少了数值波动。     

水下爆炸特性的二维数值模拟研究

 对带有3 mm左右厚铝壳或钢壳的柱形TNT、RS-211、T/Γ（35/65）以及ROT901装药在水下采用一端点起爆后,水中冲击波超压分布进行了二维数值模拟研究。对该四种炸药的每一种,分别考虑了药量为1 kg、3 kg、5 kg和8 kg四种TNT当量的柱形装药情况。柱形装药的长径比为L/D=1.5左右。计算使用的程序为二维欧拉多流体网格法流体力学计算程序MFIC。给出了水中冲击波超压的空间和时间分布;对铝壳和钢壳两种情况的结果进行了分析、比较;同时给出了缩比模型与1∶1模型之间相应物理量之间的缩比关系。