随机选取法计算化学反应流爆轰过程引起的参数涨落

利用ＲＣＭ（随机选取法）解反应流爆轰问题时，掌握影响流动参数随机涨落的主要因素是选择适当计算参数、保证得到正确结果的关键。各种流动参数的随机涨落相互关联，只需要搞清楚炸药未反应质量分数的随机涨落并做出分析，其它流动参数涨落的情况也就清楚了。给出了影响未反应质量分数涨落的因数，同时还给出了具体算例，检验所得结论的正确性。     

三维动力学分析软件MSC／DYTRAN功能及其应用

简要介绍了三维动力学分析软件ＭＳＣ／ＤＹＴＲＡＮ的功能，基本特征，以及在数值计算领域如：浅滩水下爆破，弹丸高速碰撞，成型装药等问题中的应用，并给出了部分算例。     

气体动理学方法对斜界面不稳定性实验的数值模拟

介绍了基于BGK模型(Bhatnagar-Gross-Krook)的气体动理学数值方法(Gas Kinetic Scheme,GKS),通过计算多组分激波管问题对以GKS方法实现的MBCD程序进行了验证,表明了在接触间断处产生的震荡非常小,网格收敛性好。在多组分问题求解时各组分满足同温和同速的假设条件下,利用该程序开展了Air/SF6斜界面不稳定性实验的数值模拟,得出了不同时刻SF6的密度分布图,并得到交界面宽度、斜界面左右端相对位置在激波穿过斜界面后的变化情况。通过与实验结果的比较可知,GKS方法和程序在网格数为1024×200和512×100时得到的交界面宽度与实验值的最大误差分别为6.1%和7.3%,可用于对界面不稳定性问题的计算。     

激波冲击下air/SF6界面的Richtmyer-Meshkov不稳定性

实验研究了低马赫数(1.27)激波作用air/SF6界面的RM不稳定性问题.air/SF6初始正弦界面由厚度为1～2 μm的薄膜相隔得到,用阴影法测试界面演化过程.实验结果表明:由于不稳定性重流体(SF6)向轻流体(air)演化成“尖钉”结构,而轻流体演化为“气泡”结构;由于界面切向速度差的Kelvin-Helm-ho...     

钝感炸药爆轰波传播和驱动的数值模拟研究

为实现对非理想爆轰波传播和爆轰驱动问题的精确计算,研究了与贴体坐标下Level Set方法相结合的爆轰冲击波动力学(DSD)计算方法;根据Hamilton-Jacobi方程的Godunov差分格式,提出了非正交的贴体坐标系中Level Set函数方程的差分格式及其相应的数值方法;并将该DSD方法与二维流体动力学程序TD...     

激波管实验和果冻实验界面不稳定性数值计算

运用算子分裂技术,增加二阶空间中心差方法和两步Rung-Kutta时间推进方法计算动力学粘性以及热流部分对流场的影响,将可压缩多介质流体动力学高精度欧拉计算方法MFPPM发展到适用于NS(Nayier Stokes)方程的可压缩多介质粘性流体计算方法MVPPM.通过与界面不稳定性实验结果的比较,来检验计算方法的正确性,并验证计算程序的有效性.主要包括一个激波管实验和两个果冻实验,即英国AWE(Atomic Weapons Establishment)激波管实验和LLNL实验室的果冻环实验以及中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室进行的爆轰驱动下单层果冻界面不稳定性实验.研究结果表明:数值模拟结果与实验测量结果以及对应时刻的实验图象均吻合较好.     

二维流体力学有限差分和有限元程序的耦合计算

根据压力和法向速度连续准则，将Euler方法为基础的MFPPM（Piecewise-Parabolic Method）程序和Lagrange方法为基础的DEFEL（2-D Finite ElementsCode，二维流体弹塑性动力有限元）程序进行耦合，发展了基于levelset的GEL（Ghost—fluid Euler-Lagrange）方法。该方法在处理大变形流场与小变形结构以及复杂流动与多物体相互作用等问题具有优越性。通过二维算例的计算结果与文献比较，检验了GEL方法和耦合程序的正确性，并对水下爆炸形成的流场对多物体作用过程进行了数值模拟。     

Euler-Lagrange耦合计算的GEL方法

研究了一种Euler-Lagrange耦合数值方法——GEL（Ghost-fluid Euler-Lagrange）方法，编写了GEL二维计算程序。     

重新初始化的LS方法跟踪二维可压缩多介质流界面运动

 用非耦合求解方法计算Level Set函数方程与流体力学方程组，应用重新初始化的Level Set函数确保距离函数性质，流体力学方程组采用二阶精度多介质流波传播差分格式计算，重新初始化方程采用五阶WENO格式计算。并给出了二维可压缩多介质流界面运动的计算结果。     

准等熵压缩流场反演技术研究

 与一般流体力学计算中将流体力学方程组进行空间离散、再进行时间推进计算不同,采用Hayes提出的流场反演方法,将流体力学方程组进行时间离散,然后将VISAR测量得到的样品后自由面速度历史作为输入数据,再进行空间的反演计算。与一般流体力学计算进行了比较,并与柏劲松设计的Pillow飞片撞击铜样品实验的计算结果作了比较,该方法反演出的样品加载面处的压力历史与其计算结果基本一致,通过参数优化方法多次进行反积分计算调试出的等熵压缩线具有较高的可信度。