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

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

Richtmyer-Meshkov不稳定性湍流混合数值模拟

采用PPM方法数值求解Euler方程；采用Shyue提出的考虑压力平衡的混合网格状态方程的处理方法，完成R-M不稳定性问题后期混合的数值模拟。界面不稳定性后期混合具有明显的三维特征，二维计算不能分辨后期混合流体团之间三维扭曲拉伸现象，因此要求三维数值模拟。另一方面，界面不稳定性后期，通过非线性作用，小尺度运动被充分激发，必须模拟从大尺度到小尺度的级串现象，因此数值模拟要求很高的空间分辨率，要求大规模数值计算。由此我们采用MPI、应用区域分解方法完成程序并行化，并行程序具有较好的可扩展性。     

炸药爆轰产物Jones-Wilkins-Lee状态方程不确定参数

Jones-Wilkins-Lee (JWL)状态方程是一种不显含化学反应、由实验方法确定参数的半经验状态方程, 能比较精确地描述爆轰产物的膨胀驱动做功过程. 在JWL状态方程中有多个未知(不确定)参数需要确定. 传统的确定JWL状态方程参数的方法是“调参数”, 人为因素影响较大, 无法给出参数的不确定性信息. 本文利用贝叶斯分析方法研究了炸药的不确定参数, 该方法能够基于以往的认识、实验和模拟数据标定(calibration)不确定参数. 在本文结果中, 不确定参数的后验分布均值与文献结果相符合, 基于参数标定结果的数值模拟90%置信区间完全包含实验数据. 数值标定结果说明贝叶斯参数标定适用于确定样品炸药的JWL状态方程参数. 特别是, 在本文JWL状态方程参数标定过程中极大地减少了人为因素的影响.     

含尘气体Rayleigh-Taylor不稳定性诱导混合的多相浮阻力模型研究

文章首先采用单相浮阻力模型对不同加速度下Rayleigh-Taylor不稳定性诱发的物质渗透边界的演化过程进行了计算, 揭示了该混合在常加速度和变加速度情况下不同的发展规律, 并通过与实验结果的比较分析, 验证了该模型的适用性.在此基础上, 发展了多相浮阻力模型, 采用该模型对常加速度情况下含尘气体中的Rayleigh-Taylor不稳定性诱导混合进行了研究, 发现混合区宽度随着颗粒数密度和颗粒尺寸的增大而减小, 揭示了气体中所含杂质抑制混合发展的规律.      

适应强扭曲网格的扩散方程时空二阶精度计算

以全局支撑算子方法为基础,通过引入面通量,构造了具有局部模板点的时空二阶精度格式。对于大变形扭曲网格,格式采用法向修正技术和合理的单元角体积计算方法,可以保持通量的精确性。算例表明该格式在非凸网格上能够精确获得线性解; 在非光滑网格上可以达到时空二阶精度; 能够较好地保持对称性; 并适合于三维非结构网格上的求解。     

双层Kidder自相似解及其Rayleigh-Taylor不稳定性研究

将单层Kidder自相似解推广到双层,使得两层壳体的交界面两侧存在密度跳跃,使得轻流体向重流体加速产生Rayleigh-Taylor不稳定性;通过采用Lagrange坐标下的Godunov方法进行一维直接数值模拟,将模拟解与双层Kidder自相似基本解进行比较,验证了双层Kidder自相似解的可靠性;最后,通过编制球形内爆的三维扰动的线性稳定性分析程序,对双层Kidder自相似解的Rayleigh-Taylor不稳定性进行了分析计算.计算结果表明:初始扰动越集中于交界面,会造成后期扰动增长得越快,越不稳定;扰动波数越大,扰动增长得越快,越不稳定;从扰动在空间上的发展来看,可压缩性研究表明内外壳体的可压缩性对扰动增长起着相反的作用,外层壳体的可压缩性对Rayleigh-Taylor不稳定起失稳作用,而内层壳体的可压缩性对Rayleigh-Taylor不稳定起致稳作用.     

Linear Rayleigh-Taylor instability analysis of double-shell Kidder''s self-similar implosion solution

This paper generalizes the single-shell Kidder's self-similar solution to the double-shell one with a discontinuity in density across the interface. An isentropic implosion model is constructed to study the Rayleigh-Taylor instability for the implosion compression, A Godunov-type method in the Lagrangian coordinates is used to compute the one-dimensional Euler equation with the initial and boundary conditions for the double-shell Kidder's self-similar solution in spherical geometry. Numerical results are ob...     

冲击波加载速率对微射流影响的数值模拟

材料表面在发生熔化前，微射流可能是微喷射的主要物理机制之一。曾鉴荣等在纯铅的实验中发现，当靶板中出现三波结构（即弹性先驱波、相变波和塑性波）时，测得峰值压力为22GPa时纯铅样品的微喷射量比峰值压力为20GPa的单次冲击加载喷射量几乎减少了1／2。Asay在铝平面样品的微喷射实验中，也发现随着冲击波加载速率的减小（上升沿宽度增加），喷射量大致按指数规律减小。对于自由面上缺陷平均尺度为5lain的样品，在冲击加载变到35ns波阵面宽度的加载条件时，喷射量约降低了2个数量级。     

材料特性对爆轰驱动飞层对碰凸起影响的数值模拟

应用多介质弹塑性L-R两步欧拉计算方法,对炸药两端起爆情况下柱面爆轰驱动飞层对碰凸起现象进行了数值模拟。在飞层为铅材料情况下,给出了铅飞层熔化区域及断裂图像,将对碰凸起计算结果与实验X光图像进行比较的基础上,进一步探索凸起现象的形成机制。从材料的熔点、弹塑性、密度、声速等方面对金属飞层对碰区凸起的影响进行了分析。综合分析几个计算模型的数值模拟结果,初步认为:材料强度大、密度大、熔点高能使飞层对碰凸起钝化,而材料声速较小会造成对碰区波系关系发生变化且正压作用时间增加,从而形成较严重的对碰凸起,材料声速大小对对碰凸起的影响较强度、密度等因素更明显。     

FeMnNi合金冲击相变与卸载逆相变历程研究

采用VISAR测试技术,将正向加载和逆向加载实验相结合,研究了不同加载压力下FeMnNi合金的冲击相变和卸载逆相变历程.结果表明:冲击加载压力与冲击相变阈值之差小于逆相变阈值时,即0＜σp-σt＜σn,FeMnNi合金发生α→ε相变,加载段呈现弹性波、塑性波和相变波三波结构,塑性波在自由面反射稀疏波使FeMnNi合金完...