基于多普勒测速技术的JB-9014炸药反应区结构研究

为了解TATB基JB-9014炸药的爆轰过程,利用火炮驱动飞片加载,采用光子多普勒测速技术,对JB-9014炸药的爆轰反应区结构进行了实验研究。实验中利用火炮发射高速蓝宝石飞片冲击起爆被测炸药,在炸药后表面安装镀膜氟化锂（LiF）窗口测量炸药爆轰时的界面粒子速度,测试过程的时间分辨率小于2 ns。将粒子速度剖面对时间进行一阶求导,通过一阶导数的拐点来确定炸药反应区宽度、反应时间。研究结果表明,钝感炸药JB-9014的反应时间为（0.26±0.02）μs,对应的化学反应区宽度为（1.5±0.2）mm,反应结束点处的压力为27.3 GPa,von Neumann峰处压力为40.3 GPa。     

起爆高密度TATB炸药的飞片速度阈值

采用实验与数值模拟相合的方法研究了飞片起爆钝感TATB炸药的性能。设计了飞片可靠起爆钝感TATB炸药的起爆序列,利用全光纤位移干涉测速系统分别测出飞片可靠起爆和未起爆TATB炸药的速度,初步确定飞片起爆钝感TATB炸药的起爆速度阈值。采用DYNA2D程序对飞片起爆TATB炸药过程进行数值模拟,模拟结果与实验结果相符合。     

波后慢放能过程对钝感炸药爆轰过程的影响

采用唯象反应率模型与慢放能函数相结合方法,研究钝感炸药波后缓慢放能过程对爆轰过程的影响,计算结果与实验对比分析表明,慢放能反应率模型可更好地描述钝感炸药的爆轰物理过程.     

炸药TATB/粘结剂的短脉冲冲击起爆阈值测量

利用金属箔电爆炸驱动聚酯薄膜飞片产生短脉冲冲击波的加载技术(电炮),依据DRM(Delayed Robbins-Monro)实验程序,研究了炸药TATB/粘结剂在各种激励条件下的短脉冲冲击起爆特性,获得了其50%起爆概率条件下的短脉冲冲击起爆阈值,及实验条件下其短脉冲冲击起爆判据为:lnp=0.205-0.87lnτ。并与相同条件下TATB/HMX为基的高聚物粘结炸药的短脉冲冲击起爆阈值进行了比较。结果表明飞片面积和压力脉宽均对炸药的短脉冲冲击起爆阈值有重要影响,与TATB/HMX为基的高聚物粘结炸药相比,TATB/粘结剂更钝感。     

用三项式点火增长反应速率模型计算JB-9014的反应区

在一维流体动力学编码SSS程序中,利用三项式点火增长反应速率模型对钝感炸药JB-9014的反应区进行了计算。计算中,未反应炸药采用固态HOM状态方程,产物采用气态JWL状态方程,计算得到了钝感炸药JB-9014化学反应区的峰值压力、CJ压力、反应区宽度和反应时间等参数,计算结果与实验结果符合较好。分析表明,三项式点火增长反应速率模型可用于研究钝感炸药的反应区结构。     

多级炸药爆轰高速驱动技术的Gurney模型优化分析

利用多级炸药 飞片装置可得到速度在10km/s以上、直径数厘米的末级平面飞片,在高压物态方程和高速碰撞等实验研究中有重要应用前景。改进了爆炸驱动的Gurney模型,推广用于典型的二级爆轰驱动(前级飞片 炸药层 后级飞片)装置的分析,得到前、后级飞片终速度的计算公式,与实验结果符合较好。根据对炸药层厚度优化计算结果的讨论可知,对于综合利用前级飞片动能和炸药层能量最有利的后级飞片速度处于2～3倍Gurney速度的范围。前级飞片速度很高时,后级飞片速度趋于两者直接弹性碰撞的结果。     

钝感炸药点火增长模型的欧拉数值模拟

在自主研发的二维多介质欧拉弹塑性流体力学程序中,通过引入点火增长的反应率模型以及炸药 减敏模型,借助网格自适应技术,研究钝感炸药的冲击点火、直径效应以及死区形成等爆轰现象。数值模拟结 果表明,该程序能够正确模拟平面爆轰波的爆速、CJ状态、vonNeumann尖点等爆轰参数;并能够较好模拟炸 药的直径效应。另外,通过引入考虑减敏效应的反应率模型,能较好地模拟钝感炸药的死区形成过程。     

钝感炸药的超压爆轰与冲击起爆过程数值模拟

采用Hybrid反应率结合修正的JWL方程,研究了LX-17、超细TATB等钝感炸药的冲击起爆（SDT）过程,并计算了爆轰波的对碰现象。结果表明,该方法计算钝感炸药的冲击起爆过程与实验数据符合较好;计算爆轰波对碰区的峰值压力提高了10％。     

一种以TATB/HMX为基炸药的到爆轰距离

利用光纤探针/光电转换器/示波器技术和楔形炸药块方法,采取电炮驱动薄飞片冲击加载产生高压短脉冲激励的装置,测量了一种以TATB/HMX为基炸药的到爆轰距离,研究了到爆轰距离与初始冲击波压力幅值和脉宽的关系,给出了该炸药的POP曲线和表达式。实验结果有助于了解压力脉宽对钝感炸药的冲击起爆和爆轰成长的影响。对于短脉冲冲击加载,入射压力脉宽对炸药的到爆轰距离影响明显,相同实验条件下,压力脉宽越长,炸药的到爆轰距离越短;相同压力脉宽下,加载压力越高,炸药的到爆轰距离越短。     

含铝炸药爆轰驱动的非线性特征线模型

针对含铝炸药爆轰的非理想特性,提出了含铝炸药爆轰产物膨胀的局部等熵假设,建立含铝炸药爆轰驱动的非线性特征线模型,为研究含铝炸药爆轰产物的非等熵流动和膨胀做功提供了一种新的理论分析方法。设计了5、50 μm含铝炸药和含LiF炸药驱动0.5、1 mm厚金属板实验,通过激光位移干涉仪测试金属板运动的速度历程,再通过实验结果计算得到铝粉在爆轰产物中的反应度变化规律,结合含铝炸药爆轰产物的非线性特征线模型,理论计算了含铝炸药驱动金属板的速度历程。对比理论与实验结果,理论方法能够很好地描述铝粉二次反应对炸药做功能力的贡献,同时验证了含铝炸药爆轰驱动的非线性特征线模型的正确性。     

RDX基含铝炸药爆轰波结构实验研究

为了获得含铝炸药爆轰反应区附近铝粉的反应情况,对两种RDX/Al炸药和一种RDX/LiF炸药的爆轰波结构进行了测量。实验过程中,利用火炮加载产生一维平面波,通过光子多普勒测速仪测量炸药/LiF窗口的界面粒子速度。结果表明:含铝炸药爆轰波的结构与理想炸药的差异较大,其界面粒子速度曲线没有明显的拐点;反应初期,由于气相产物与添加物之间温度的非平衡性,RDX/Al界面的粒子速度低于RDX/LiF炸药的;随后,由于铝粉反应放能,RDX/Al界面的粒子速度高于RDX/LiF炸药的;微米尺度铝粉在CJ面前几乎不发生反应;2、10 μm等两种粒度铝粉的反应延滞时间小于0.8 μs;在本文中,两种粒度铝粉的反应度为16%~31%。     

基于深度神经网络的横流转捩预测

基于深度神经网络DNN构建了从层流流场无量纲速度梯度、流向涡强度等物理量到横流转捩模态下间歇因子间的映射关系,获得一种新的数据驱动转捩模型.通过将数据驱动转捩模型与SST k-ω湍流模型耦合,有效简化了转捩模型输运方程求解,实现高效、准确的亚音速三维边界层横流转捩流场计算. DNN训练数据来自变雷诺数的NLF(2)-0415无限展长后掠翼计算结果,并以两种工况进行测试,数据驱动转捩模型预测精度与γ-Reθ转捩模型近似.将数据驱动转捩模型用于其他典型横流转捩算例的计算,以验证其泛化能力.对于变后掠角的NLF(2)-0415后掠翼,数据驱动转捩模型与γ-Reθt-CF模型预测的转捩位置几乎一致,并且能够预测出后掠角从45°增长到65°的过程中,转捩位置先向前再向后移动的现象;对于标准椭球体,使用低分辨率网格进行计算,数据驱动转捩模型依然能够实现转捩位置预测,对椭球体表面Cf的计算结果与多个平台的横流转捩模型、实验结果基本一致.研究表明,以横流转捩相关物理量作为输入对DNN进行训练,并将获得的数据驱动转捩模型与SST k-ω湍流模型耦合,可以实现对横流转捩的有效预测,且具有较强的泛化能力.数...     

格子玻尔兹曼法多块网格的交界结构优化研究

基于格子玻尔兹曼方法LBM(Lattice Boltzmann Method)对多块网格方法(Multi-Block)的粗细网格交界结构进行了研究,提出了一种新的优化处理方案。解决了原有网格交界结构存在的三个问题,即两套插值运算造成的程序结构复杂的问题,存储前几个时间步的节点流场数据以备插值运算造成内存浪费的问题和基于时间插值结果进行空间插值计算造成插值误差积累的问题。用一次多点二维空间插值的方式,将原方法的空间和时间双插值,简并成一次空间插值。通过对经典的非定常的圆柱绕流算例和定常的标准顶盖方腔驱动流算例的仿真模拟,验证了交界面处质量、动量及应力的连续性以及网格交界面数据过渡的流畅度,最终验证了改进方法的正确性。数值模拟结果表明,改进后多块算法可实现局部网格细化,进一步推动LBM方法在实际工程问题中的应用。     

Dual Mesh Method for Upscaling in Waterflood Simulation

Detailed geological models typically contain many more cells than can be accommodated by reservoir simulation due to computer time and memory constraints. However, recovery predictions performed on a coarser upscaled mesh are inevitably less accurate than those performed on the initial fine mesh. Recent studies have shown how to use both coarse and fine mesh information during waterflooding simulations. In this paper, we present an extension of the dual mesh method (Verdière and Guérillot, 1996) which simulates water flooding injection using both the coarse and the original fine mesh information. The pressure field is first calculated on the coarse mesh. This information is used to estimate the pressure field within each coarse cell and then phase saturations are updated on the fine mesh. This method avoids the most time consuming step of reservoir simulation, namely solving for the pressure field on the fine grid. A conventional finite difference IMPES scheme is used considering a two phase fluid with gravity and vertical wells. Two upscaling methodologies are used and compared for averaging the coarse grid properties: geometric average and the pressure solve method. A series of test cases show that the method provides predictions similar to those of full fine grid simulations but using less computer time.     

点引爆圆板装药驱动飞板的实验及二维数值模拟

本文报导点引爆圆板装药驱动飞板的实验及二维数值模拟研究结果。实验测定了有机玻璃、铝、钢和铜四种飞板在爆轰产物驱动下的碰靶时间及碰靶波形,并将二维数值模拟结果与测量结果进行了比较,二者符合较好。文中还就数值模拟中存在的一些问题进行了简要分析、讨论。     

Flexible Unstructured Gridding for 2D Reservoir Coarse Scale Modeling using Fine Scale Permeability Filtering

Accurate up-scaling is an essential part of creating a valid reservoir coarse scale dynamic model. In this article, unstructured discretization of spatial domain is accompanied by numerical permeability up-scaling in order to construct an accurate coarse scale model. A new technique for generating a course scale triangular mesh is presented in which the density of elements in key flow regions is kept high to capture accuracy. The fine scale permeability map is investigated using image processing techniques, especially steerable filters, and the results are converted into a high-resolution element size map. This element size map will be refined by the integration of other important factors such as well-position effects and used to construct a coarse triangular mesh. The combination of flux-continuous pressure approximation and mass conservative, total variation diminishing finite volume schemes have been considered to solve two phase flow equations on the control volume finite element mesh. Fine scale simulations results are compared with the coarse scale ones for a series of water flooding examples to investigate the efficiency and accuracy of the presented gridding methodology. This method is developed for 2D cases, but can be easily extended to 3D problems.     

A simplified adaptive Cartesian grid system for solving the 2D shallow water equations

This paper presents a new simplified grid system that provides local refinement and dynamic adaptation for solving the 2D shallow water equations (SWEs). Local refinement is realized by simply specifying different subdivision levels to the cells on a background uniform coarse grid that covers the computational domain. On such a non‐uniform grid, the structured property of a regular Cartesian mesh is maintained and neighbor information is determined by simple algebraic relationships, i.e. data structure becomes unnecessary. Dynamic grid adaptation is achieved by changing the subdivision level of a background cell. Therefore, grid generation and adaptation is greatly simplified and straightforward to implement. The new adaptive grid‐based SWE solver is tested by applying it to simulate three idealized test cases and promising results are obtained. The new grid system offers a simplified alternative to the existing approaches for providing adaptive mesh refinement in computational fluid dynamics. Copyright © 2011 John Wiley & Sons, Ltd.