气相爆轰波马赫反射影响因素的实验研究北大核心CSCD

在矩形截面的爆轰管道中,对C2H2+2.5O2+8.17Ar和C2H2+5N2O在CJ爆轰状态下经过不同楔面所发生的马赫反射的影响因素进行了实验研究。实验中,由烟膜记录爆轰波马赫反射的胞格结构转变过程;采用纹影技术捕捉爆轰波马赫反射波阵面的不稳定性及波后流场分布。实验结果表明:两种实验气体在爆轰波马赫反射过程中均存在由CJ区域向过驱区域转变的胞格结构;初始压力对楔面与马赫反射三波点轨迹线之间的夹角(χ)影响明显,楔角θw对χ的影响随θw的增大而增大;根据实验测得的θw+χ与θw之间的关系,可知爆轰波马赫反射三波点轨迹线的斜率随着θw的增大而增大,与CJ区域内胞格轨迹线的相交距离也更短,使马赫杆后的过驱度升高。另外,不稳定气体C2H2+5N2O的不稳定性高于稳定气体C2H2+2.5O2+8.17Ar,导致二者的爆轰波马赫反射行为存在较大的差异。     

柱面Love波频散分析与SH波场的数值计算

对贴井壁环型剪切源在柱状双层弹性介质中激发的SH波场进行了理论求解,导出了柱面Love波频散方程,讨论了柱面Love波存在的条件及其区域.通过数值计算考察了柱面Love波的频散特性和激发强度,发现最低阶柱面Love波具有截止频率,这与平面半空间双层弹性介质模型下的Love波无截止频率的特征不同.渐近分析与数值考察都表明,井径r₁→∞时,柱面Love波频散方程趋向平面双层半空间的Love波方程,柱面Love波的截止频率趋于零.全波计算还显示用激发SH波来探测侵入带外原状地层的横波信息是一个十分简洁的途径..     

铝板中Lamb波检测的实验研究

Lamb波在传播过程中具有频散及多模态特性,若相关参数选择不当,会导致在实际应用中信号相互叠加而无法识别。本文基于Lamb波的频散曲线是其频散方程实数解分布的特点,采用二分法绘制了铝板中Lamb波的频散曲线、波结构曲线和入射角曲线。根据曲线选择S0模态的Lamb波对1mm厚铝板中不同类型的缺陷进行检测。实验结果表明,S0模态的Lamb波对裂纹型缺陷和贯穿型缺陷十分敏感,但对于裂纹型缺陷,其幅值变化并不与缺陷大小成线性关系,并且S0模态Lamb波的声场指向性十分集中,在偏离声束轴线时无法检测到缺陷。     

联合简正波水波和底波频散特性的贝叶斯地声参数反演

大多数基于浅海简正波模态频散数据的地声参数反演方法无法对深层底质声学参数进行可靠估计,究其原因是仅利用了简正波水波频散特征,忽略了与深层底质声学参数密切相关的底波频散特征,因此,本文在分析了包含水波和底波的浅海宽带数据的基础上,结合底波频散特征对深层底质声学参数变化更加敏感的物理特性,实现了基于完整简正波频散特性的贝叶斯地声参数反演,并针对简正波宽带声场模型计算复杂度较高的现实问题,利用变分贝叶斯蒙特卡罗方法的推断优势,完成了未知参数的可靠估计和快速后验分析.仿真和海上实验结果表明:联合简正波水波和底波频散特征数据的贝叶斯地声参数反演,不仅可以有效估计深层底质声学参数,而且降低了其他相关环境参数的估计不确定性.     

曲率相关的爆轰波在非结构四边形网格上的数值方法

假设爆轰波阵面的法向速度是曲率的线性函数,在非结构四边形网格上采用水平集方法模拟爆轰波阵面的运动过程.水平集方程的曲率无关项采用正格式离散,曲率项采用伽辽金等参有限元方法空间离散,时间离散采用半隐格式.在笛卡儿网格和随机网格上,含曲率的水平集方程的离散格式为强一阶精度,重新初始化方程的离散格式精度为近似一阶精度.曲率收缩的不光滑界面和多个爆轰波阵面相互作用的算例说明格式可有效地模拟爆轰波与曲率相关的运动.     

一种捕捉多向刚性爆轰波的激波定位随机投影方法

提出一种激波定位随机投影方法,利用当地随机投影,引入激波指示器并提出恰当的判定当地投影区域的规则,使爆轰波多向推进的问题得以解决。数值实验表明,激波定位随机投影方法不仅适用于爆轰波的多向推进,还能高效求解刚性反应方程,易于推广到三维情况。     

一维爆轰波在扰动来流中传播的数值研究

爆轰波在静止气体或定常来流中的传播得到了广泛研究, 然而在扰动来流中的传播研究较少。这方面的研究不仅是爆轰传播机制的重要组成部分, 还可为爆轰发动机的应用提供参考。文章基于两步诱导-放热总包反应模型, 开展了一维爆轰波在正弦密度扰动来流中的传播数值模拟。通过对数值结果分析, 获得了放热反应控制参数与爆轰波内在不稳定性的关系, 并在此基础上研究了扰动波长和幅值对一维爆轰波动力学过程的影响。研究发现, 在波前施加连续扰动会诱导爆轰波表现出更复杂的动力学行为, 且影响过程与爆轰波的内在不稳定性相关。对于稳定爆轰波, 扰动只在特定波长范围内引起前导激波后的压力振荡。对于不稳定爆轰波, 扰动会进一步强化其内在不稳定性。扰动幅值越大, 对爆轰波动力学过程的影响越显著。      

基于静动态重构的高阶DG/FV混合格式在二维非结构网格中的推广

通过比较间断Galerkin有限元方法(DGM)和有限体积方法(FVM),提出"静态重构"和"动态重构"的概念,进一步建立基于静动态"混合重构"算法的三阶DG/FV混合格式.在DG/FV混合格式中,单元平均值和一阶导数由DGM方法"动态重构",二阶导数利用FVM方法"静态重构";在此基础上,构造高阶多项式插值函数,得到...     

一种基于简正波模态频散的远距离宽带海底参数反演方法

在浅海环境中, 海底环境参数对声传播有着重要的影响. 由于利用单个宽带声源进行海底参数反演时, 随着距离的增大, 误差变大, 本文提出利用warping变换对在浅海波导中传播的, 不同距离上的两个宽带爆炸声源进行简正波的有效分离, 实现了宽带爆炸声源的远距离海底参数反演. 采用全局寻优遗传算法对提取出的模态频散到达时间差与理论计算的模态频散到达时间差进行匹配处理, 并结合随距离连续变化的声传播损失, 实现了利用单水听器进行海底参数的反演. 实验结果表明: 运用反演出的海底参数提取模态频散时间差和实测数据提取出的模态频散时间差吻合得较好; 而通过传播损失反演得到的海底衰减系数与频率呈指数关系. 最后, 对反演结果进行了后验概率分析, 并将本组爆炸声源的反演结果用于另一组不同距离上爆炸声源时仍然有效, 来评价反演结果的有效性.     

一种基于简正波模态消频散变换的声源距离深度估计方法

针对浅海环境中传播的低频宽带水声脉冲信号,基于简正波水平波数差和波导不变量之间的关系,本文提出了一种利用距离-频散参数二维平面聚焦测距与匹配模态能量定深的目标声源定位方法.首先,通过将由频散参数和波导不变量表示的前几阶模态相速度与由环境模型计算的相速度进行对比分析,从而估计出前几阶模态的频散参数和环境的波导不变量.其次,利用估计出的频散参数值和波导不变量对接收信号进行消频散变换处理,只有当接收信号的距离参数等于目标声源距离时,各号简正波的幅度均达到最大值,在距离-频散参数二维平面上,出现声压聚焦的现象,利用此现象可以估计目标声源的距离.不仅如此,消频散变换后的接收信号,前几阶模态在时域上明显地分离开来,可以准确地估计出前几阶模态的能量,采用多模态能量匹配的方式,可以估计出目标声源的深度.最后,通过对仿真和冬季获得的气枪信号数据处理结果验证了本文方法的有效性.     

CTVD格式数值计算非均质炸药爆轰问题

将高分辨率激波捕捉格式CTVD格式拓展应用到非均质炸药爆轰的数值模拟问题.增加了化学反应率控制方程,引入Lee-Tarver点火成长模型,未反应的固体炸药和化学反应气体产物都使用JWL形式状态方程.数值模拟了非均质固体炸药PBX-9404和TATB的冲击起爆问题.获得了较高的爆轰波分辨率和光滑解区的数值精度,对具有复杂物态方程形式的固体炸药爆轰问题,CTVD格式具有简单实用、高效和高分辨的特点.     

一种基于特征理论模拟凝聚炸药爆轰的单元中心型Lagrange方法

提出一种数值模拟凝聚炸药爆轰问题的单元中心型Lagrange方法.利用有限体积离散爆轰反应流动方程组,基于双曲型偏微分方程组的特征理论获得离散网格节点的速度与压力,获得的网格节点速度与压力用于更新网格节点位置以及计算网格单元边的数值通量.以这种方式获得的网格节点解是一种"真正多维"的理论解,是一维Godunov格式在二维Riemann问题的推广.有限体积离散得到的爆轰反应流动的半离散系统使用一种显-隐Runge-Kutta格式来离散求解：显式格式处理对流项,隐式格式处理化学反应刚性源项.算例表明,提出的单元中心型Lagrange方法能够较好地模拟凝聚炸药的爆轰反应流动.     

耗散紧致格式的频谱特性研究与应用

采用Fourier分析方法,给出线性耗散紧致格式和基于m级Runge-Kutta时间积分方法的全离散格式的频谱特性,并应用五阶耗散紧致格式模拟典型高频波传播和超声速平面Couette流动的特征值问题及其稳定性边值问题,展示耗散紧致差分格式良好的频谱特性.     

A One-Dimensional Discrete Boltzmann Model for Detonation and an Abnormal Detonation Phenomenon

A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem,and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model,the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition.The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.     

黏性激波管的有限体积WENO格式数值模拟

本文采用OpenFOAM软件下实现的一种可实现任意阶数,可应用于非结构网格的有限体积WENO格式对黏性激波管问题进行模拟.模拟中对流项离散采用3阶精度、4阶精度该类WENO格式,网格形式采用结构网格和三角形非结构网格.结果表明,采用该类格式,三角形非结构网格的算精度、效率优于结构网格,3阶精度格式计算效率优于4阶精度....     

非结构四边形网格下的一类保对称有限体元格式

针对定常扩散问题,在非结构四边形网格下,通过选取特殊的控制体和有限体元空间,建立两种保对称有限体元格式,在拟一致网格剖分下,当扩散系数光滑时,证明有限体元解函数在L²和H¹范数下均具有饱和误差阶.数值实验验证理论结果的正确性,同时表明新格式对扭曲大变形四边形网格、间断系数问题具有较强的适应性.在正交网格下,第二种格式对流(flux)函数在单元中心点的值还具有超逼近性.     

基于连续旋转爆震的推进技术研究进展

基于爆震燃烧的推进技术是未来空间技术的重要发展趋势,特别是可实现结构简单化设计和高热力学效率.针对火箭式连续旋转爆震发动机、吸气式爆震发动机的实验测试和数值仿真,文章综述了其国内外研究进展,分别总结了不同燃料、燃烧室结构、喷注方式以及工作方式等对连续旋转爆震波的传播规律和发动机的特性影响规律.虽然上述探索性研究得到了诸多有益的结论,但是由于连续旋转爆震燃烧技术涉及的流动、物理化学过程十分复杂,对旋转爆震燃烧的机理研究仍然有待进一步深入开展.      

铝粉-空气混合物爆燃转爆轰过程及爆轰波结构

 在长为32.4 m、内径为0.199 m的大型长直水平管道中,对铝粉-空气两相流的燃烧转爆轰（DDT）过程及爆轰波结构进行了实验研究。对铝粉-空气混合物弱点火条件下DDT过程不同阶段的特征进行了分析,实验结果显示混合物经历了缓慢反应压缩阶段、压缩波加速冲击波形成阶段、冲击反应过渡阶段、冲击反应向过压爆轰过渡阶段和爆轰阶段,得到了混合物各阶段的DDT参数,由此进一步分析了DDT浓度的上、下限。在1.4 m爆轰测试段的4个截面的环向上各均匀安装8个传感器,对爆轰波结果进行测试,并对铝粉-空气混合物爆轰波的单头结构进行了分析。     

An Improved CE/SE Scheme and Its Application to Detonation Propagation

An improved two-dimensional space-time conservation element and solution element （CE/SE） method with second-order accuracy is proposed, examined and extended to simulate the detonation propagations using detailed chemical reaction models. The numerical results of planar and cellular detonation are compared with corresponding results by the Chapman-Jouguet theory and experiments, and prove that the method is a new reliable way for numerical simulations of detonation propagation.     

爆轰波在弯管内传播过程数值分析

应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波在弯管内的传播过程进行了数值模拟.计算中氢氧混合物化学反应采用了8种组分20个反应方程式.在处理化学反应引起的刚性问题时采用了时间算子分裂的方法.计算结果表明,在弯管小曲率半径壁面附近,由于膨胀稀疏作用,爆轰波强度减弱,在局部出现前导激波与放热反应区的解藕以及二次起爆现象;在弯管大曲率半径壁面上爆轰波在马赫反射和正规反射之间相互转变,使爆轰波加强.弯管内的爆轰现象与弯管曲率半径有关.