应用压力间断解析解修正AUSM+-up格式

AUSM+-up格式是近年来发展的AUSM类格式,通过在对流通量中引入压力耗散,同时在压力通量中引入速度耗散,使之在各种马赫数下均有良好的收敛性和稳定性。然而,通过对格式的分析,我们发现当马赫数为零时,在压力间断位置会出现质量通量无穷大的非物理现象。为了克服这一缺陷,采用特征线法获得了压力间断问题的解析解;以此为基础,对对流通量中的压力耗散项进行修正,发展出新的AUSM+-up格式。以新的格式对不同速度驱动下的冲击波发展、不同强度压力间断问题和激波衍射问题进行了仿真,计算结果与理论解吻合得较好。     

迎风格式在接触间断的数值耗散及其诱导误差

在不同流场参数条件下对三种迎风格式在接触间断中的数值耗散问题进行了数值实验,并对数值耗散产生的机理进行了分析。数值计算结果和理论分析表明,矢通量分裂格式计算接触间断问题时,若流场静止或流场内存在亚声速区域,密度耗散的产生会诱导出以特征速度运动的数值扰动误差,该误差对数值耗散的大小无影响,但会影响流场的速度及压力参数分布,从而改变流场的结构。在二维问题中,诱导误差相互干扰会产生大量的复杂小尺度结构,给流场结构分析带来困难。同时,在密度参数线性分布的流场中,若空间离散格式重构的对象为对流通量,使用矢通量分裂格式计算流场会产生数值误差,使计算精度难以到达二阶。     

WENO格式权重分析与改进

为了优化WENO格式计算性能,在对Jiang和Shu的经典WENO格式(记为WENO-JS)加权方法分析的基础上,通过引入间接光滑指数,构造出一种新的WENO格式——WENO-E格式,取得减小间断区耗散的效果。理论分析表明,该格式与WENO-JS格式计算效率基本相同,可达到相同阶的计算精度;但在相同网格下,较之WENO-JS格式,该格式对光滑区域的求解有更小的截断误差,对间断的捕捉有更高的分辨率。与WENO-JS格式相比,采用WENO-E格式进行线性迁移方程、非线性Burgers方程、欧拉方程等相关问题的数值实验,均能取得更好的数值结果。     

半线性热传导方程Cauchy问题全局经典解的存在性与唯一性

本文证明了热传导方程Cauchy 问题 u_t=u_xx+u^a x∈R,t>0 u|_t=0=φ(x) 当实数α>3 时,只要初值φ(x)在某些Sobolev 空间中的范数充分小,就有唯一的全局经典解,且当t→+∞时,这个解具有一定的衰减性。本文所用的方法使得 Cauchy问题中的a的值同解与初值所在的空间紧密联系,a的值越大,解的性质越好。     

一种改进Van Leer格式的构造与分析

基于传统的Van Leer格式提出了一种改进的Van Leer(NVL)格式,解决了Van Leer格式在低马赫数下数值耗散过大,不能求解有粘流场等问题,同时保持Van Leer格式的计算效率.通过二维平板流动算例比较了NVL格式与传统Van Leer格式的数值耗散程度,并以NACA0012翼型的亚声速和跨声速流场以及M6机翼的跨声速流场求解为算例,与Roe格式和传统Van Leer格式进行比较,进一步验证了NVL格式的低耗散性和高效性.     

基于高精度WCNS格式的γ-Reθ转捩模型标定与应用*

在高精度数值风洞平台上,采用低速平板试验数据对基于高精度WCNS格式的γ-Re_θ转捩模型进行了标定,并在二维低速问题中进行了应用。计算结果与试验的对比表明,基于高精度WCNS格式的γ-Re_θ转捩模型可准确模拟自然转捩、旁路转捩及分离转捩的位置,并且具有较低的网格敏感性。同时也表明在中等雷诺数范围,层流区域长度和湍流区域有相同量级时,计算必须采用转捩模型才能准确模拟阻力系数。     

Zodiac算法新的不可能差分攻击

重新评估了Zodiac算法抵抗不可能差分攻击的能力。通过分析Zodiac算法的线性层,给出了Zodiac算法两条新的14轮不可能差分。利用新的不可能差分,结合Early-Abort技术对完整16轮的Zodiac算法进行了不可能差分攻击。攻击过程中一共恢复6个字节的密钥,其时间复杂度只有2_32.6次加密,数据复杂度约为2_85.6个明文,该攻击结果与已有最好的结果相比,时间复杂度降低了一个因子2₃₃。结果表明由于Zodiac算法线性层的扩散性差,使得该算法对不可能差分分析是不免疫的。     

Zodiac算法新的不可能差分攻击

#重新评估了Zodiac算法抵抗不可能差分攻击的能力。通过分析Zodiac算法的线性层,给出了Zodiac算法两条新的14轮不可能差分。利用新的不可能差分,结合Early-Abort技术对完整16轮的Zodiac算法进行了不可能差分攻击。攻击过程中一共恢复6个字节的密钥,其时间复杂度只有2_32.6次加密,数据复杂度约为2_85.6个明文,该攻击结果与已有最好的结果相比,时间复杂度降低了一个因子2₃₃。结果表明由于Zodiac算法线性层的扩散性差,使得该算法对不可能差分分析是不免疫的。     

流变抛光结合HF酸刻蚀提升熔石英元件抗激光损伤性能

为了提升紫外熔石英元件抗激光损伤性能,针对传统加工方法在加工过程中产生的破碎性缺陷和污染性缺陷,本文提出了使用磁流变抛光+HF刻蚀组合工艺提升紫外熔石英元件抗激光损伤性能的方法。磁流变抛光特有的剪切去除原理能够有效去除传统加工过程产生的破碎性缺陷,同时不产生新的破碎性缺陷。HF酸动态酸刻蚀能够有效减少加工过程中产生金属元素污染。实验结果表明：经过组合工艺处理的熔石英样品,在7J/cm₂.3ω激光通量辐照下损伤密度由0.2/mm₂降至0.008/mm₂,在8J/cm₂.3ω激光通量辐照下损伤密度由1/mm₂降至0.1/mm₂,对紫外熔石英元件损伤性能提升显著。     

高精度差分格式WNND的构造及数值实验

基于二阶NND格式,通过引入Jiang和Shu的加权思想以及具有TVD性质的三阶Runge Kutta方法,构造了一种时间、空间均达到三阶精度的WNND格式。分别以波动方程、一维Euler方程和三维全Navier Stokes方程为例,通过对WNND格式的数值结果分析表明,WNND格式引起的耗散和波动较小,并且能够高精度地分辨场间断。     

Investigation on energy output structure of explosives near-ground explosion

In order to give the energy output structure of typical explosives near-ground explosion in real ground conditions, the free-field shockwave, ground reflection shockwave and Mach wave overpressure time history of composition B explosive, RDX explosive and aluminized explosive were measured by air pressure sensors and ground pressure sensors. The shape of the free-field shock wave, ground reflection shock wave, and Mach wave and explosion flame were captured by high-speed camera. The experimental results show that, at the same horizontal distance from the initiation point, the peak overpressure of explosive shock wave of composition B explosive, both in the air and on the ground, is less than that of RDX and aluminized explosives. At a distance of 3.0 m from the initiation point, the peak overpressure of aluminized explosives is slightly less than that of RDX explosives. Owing to the exothermic effect of aluminum powder, the pressure drop of aluminized explosives is slower than that of RDX explosives. At 5.0 m from the initiation point, the peak overpressure of aluminized explosives is larger than that of RDX explosives. At the same position from the initiation point, among the three kinds of explosives, the impulse of aluminized explosives is the maximum and the impulse of composition B explosives is the minimum. With the increase of the horizontal distance from the initiation point, the height of Mach triple-points (Mach steam) of the three explosives increases gradually. At the same horizontal distance from the initiation point, there is poorly difference in the height of Mach triple-points between aluminized explosive and RDX explosive, and the height of Mach triple-points of composition B explosive is much smaller than that of other two explosives. The maximum diameter and duration of the fireball formed by aluminized explosives are the largest, followed by composition B explosive, and the maximum diameter and duration of the fireball formed by RDX explosive are the smallest.     

激波对密度间断层干扰及界面效应的数值模拟

运动激波在密度间断层传播时 ,引起密度分界面的变形 ,相应激波形状发生改变。以高密度层为例 ,采用TVD格式模拟这一过程 ,并与半解析解比较分析 ,初步讨论了激波对高密度层干扰作用引起高密度层的上扬和高密度层内激波与壁面的反射现象     

交叉激波湍流边界层干扰问题数值模拟

基于BANS方程,采用有限体积法离散计算区域,对流项使用Harten TVD格式进行流场求解,湍流部分应用Menter SST两方程湍流模型,建立了隐式LU-SGS方法的三维CFD计算程序,数值模拟了来流马赫数4.96条件下的18°×18°对称双楔产生的交叉激波/湍流边界层干扰问题,计算了不同剖面沿展向压强分布以及流场的详细结构,并将计算结果与相关文献的实验与计算结果进行了比较.     

改进的接触算法及其在光滑粒子流体动力学中的应用

采用以黎曼解描述粒子间相互作用的接触算法,模拟一维情况下包含间断的流场。在弱间断的条件下,采用弱波近似的黎曼解来描述粒子间的相互作用;在强间断的条件下,则引入非迭代黎曼解来描述粒子间的相互作用。并引入Taylor展开思想,提高了接触算法在自由边界的计算精度。该算法无需引入人为粘性,程序结构简洁。数值计算表明,改进后的接触算法能很好地描述强间断,并有效改善了原始接触算法在自由面的计算缺陷。     

Experimental study on the detonation process of a pulse detonation engine with ionized seeds

An experimental platform of a pulse detonation engine (PDE) was established to study the effect of different K₂CO₃ ionized seed mass contents on the detonation process. The pressure and ion concentration were detected in the detonation process of the PDE with different contents of ionized seeds. The initiation process of the PDE at different ignition frequencies was studied. The results show that the gas conductivity in the detonation process increased by adding ionized seeds to the PDE tube, and the conductivity increased with the increase in ionized seed mass content. With the increase in ionized seed mass content, the range of the conductivity decreased. The PDE was successfully ignited and formed a stable detonation wave at ignition frequencies of 5 Hz and 10 Hz, and the peak pressure of the stable detonation with the ignition frequency of 5 Hz was 17% higher than that with an ignition frequency of 10 Hz. The detonation wave intensity was weakened and degenerated to a shock wave that propagated in the tube without the fuel filled at the ignition frequency of 20 Hz.     

1D study of the detonation phenomenon and its influence on the interior ballistics of the combustion light gas gun

One-dimensional simulations with a detailed hydrogen-oxygen reaction mechanism have been performed to investigate detonation phenomenon in a combustion light gas gun (CLGG). Convection fluxes of the Navier-Stokes equations are solved using the WAF (weighted average flux) scheme HLLC Riemann solver. A high resolution fifth-order WENO scheme for the variable extrapolation at the volume interface and a fourth-order Runge-Kutta scheme for the time advancement are used. Validation tests of the stoichiometric hydrogen-oxygen deflagration to detonation transition process shows good agreement between the computed results and the analytical and documented solutions, demonstrating the reliability on the detonation simulation of the current scheme. Simulation results of the interior ballistic process of the CLGG show that the flame propagation experiences three distinct stages. The blast detonation wave causes a high-pressure shock and hazardous oscillations in the chamber and makes the projectile accelerate with fluctuations, but has only a small improvement to the muzzle velocity.     

含新型间断探测器的混合WCNS格式在间断无粘可压流的应用

为了让高精度数值格式在含间断和小尺度涡等复杂结构的超声速无粘可压缩流动情况下,仍能鲁棒地捕捉激波并快速得到流场高保真的模拟结果,研究了以子模板导数组合为基础的光滑度量算法,构造了精度与鲁棒性兼顾的新型间断探测器,使间断识别对小尺度涡也具有高分辨率;研究了混合加权紧致非线性格式(weighted compact nonlinear scheme, WCNS)方法,对流场中的光滑与间断区域分别使用线性与非线性加权格式求解,从而克服单一非线性格式在光滑区分辨率难以达到设计精度的问题。数值实验表明,使用新型间断探测器的混合WCNS格式对一维、二维Euler方程模拟结果良好,并且相比于在全流场使用局部特征分解的原始WCNS方法有计算效率的提高。     

Aerodynamics analysis of a hypersonic electromagnetic gun launched projectile

A hypersonic aerodynamics analysis of an electromagnetic gun (EM gun) launched projectile configuration is undertaken in order to ameliorate the basic aerodynamic characteristics in comparison with the regular projectile layout. Static margin and pendulum motion analysis models have been applied to evaluate the flight stability of a new airframe configuration. With a steady state computational fluid dynamics (CFD) simulation, the basic density, pressure and velocity contours of the EM gun projectile flow field at Mach number 5.0, 6.0 and 7.0 (angle of attack = 0°) have been analyzed. Furthermore, the static margin values are enhanced dramatically for the EM gun projectile with configuration optimization. Drag, lift and pitch property variations are all illustrated with the changes of Mach number and angle of attack. A particle ballistic calculation was completed for the pendulum analysis. The results show that the configuration optimized projectile, launched from the EM gun at Mach number 5.0 to 7.0, acts in a much more stable way than the projectiles with regular aerodynamic layout.     

基于水下脉冲放电的冲击波聚焦实验研究

冲击波的聚焦可以在聚焦区域形成局部较高的压力.基于能产生较强冲击波的水中脉冲放电声源和具有会聚作用的旋转椭球面反射罩建立了水下冲击波聚焦系统,开展了水下冲击波聚焦的实验研究,分析了水下冲击波的聚焦过程.对压力历史、轴线上的压力分布及峰值压力进行了分析,研究了旋转椭球面反射罩的聚焦特性.实验发现聚焦区域会产生负压,并导致局部空化.基于冲击动力学理论对这一现象进行了分析.实验结果表明,水下冲击波聚焦系统具有显著的聚焦效果.