凝聚介质中斜激波的反射

结合数值模拟和理论分析,研究凝聚介质中斜激波反射.采用龙格-库塔控制体积间断有限元方法,数值求解"刚性气体"状态方程形式下的欧拉方程组;理论分析凝聚介质中斜激波反射模式;运用激波极曲线理论,给出典型激波强度下正规反射向非正规反射过渡的临界角及波后状态.比较数值模拟结果和激波极曲线理论分析结果,得到典型弱、强斜激波的反射图像.     

爆轰波在炸药-金属界面上的折射分析

针对爆轰波在炸药-金属界面上折射时由实验获得的金属折射冲击波压力与经典爆轰波极曲线理论预测的压力存在显著差异这一问题, 本文展开了进一步的理论和数值模拟分析研究. 首先通过分析指出经典爆轰波极曲线理论的缺陷, 并对爆轰波极曲线理论进行了改进, 改进爆轰波极曲线理论给出了炸药爆轰波折射类型以及折射冲击作用点处的压力值. 然后发展了一个基于次特征理论来数值求解爆轰反应流动控制方程的二阶中心型Lagrange方法, 并数值模拟了一个典型的炸药爆轰波折射实验. 改进爆轰波极曲线理论和数值模拟分析结果表明, 爆轰波折射类型有三种:反射冲击波的正规折射、带Mach反射的非正规折射、无反射波的正规折射, 并且金属折射冲击波压力值随入射角增大而单调减小.     

凝聚炸药爆轰波在高声速材料界面上的折射现象分析

凝聚炸药爆轰在边界高声速材料约束下传播时,爆轰波会在约束材料界面上产生复杂的折射现象.本文针对凝聚炸药爆轰波在高声速材料界面上的折射现象展开理论和数值模拟分析.首先通过建立在爆轰ZND模型上的改进爆轰波极曲线理论给出爆轰波折射类型,然后发展一种求解爆轰反应流动方程的基于特征理论的二阶单元中心型Lagrange计算方法来数值模拟典型的爆轰波折射过程.从改进爆轰波极曲线理论和二阶Lagrange方法数值模拟给出的结果看出,凝聚炸药爆轰波在高声速材料界面上的折射类型有四种:反射冲击波的正规折射、带束缚前驱波的非正规折射、带双Mach反射的非正规折射、带λ波结构的非正规折射.     

二维平面激波折射的理论计算方法

围绕二维平面激波与物质界面作用,建立了流场波系结构的理论计算方法.首先分析了界面两侧激波独自沿界面传播过程,识别了两类正规折射和三类非正规折射;然后,根据两侧扰动沿界面传播的相对快慢,得到了三种不同的作用情形.与已有的Catherasoo方法相比,此方法:1)考虑了激波波后紧邻波阵面流场中扰动对界面附近折射的影响,区分了扰动的强弱及其是否可以追赶上激波阵面; 2)除了少数情形外,大都基于Euler方程精确解描述流场中波系相互作用.利用该方法对马赫数1.17的激波与空气/SF6界面作用进行计算,得到的波系结构与数值结果、已有实验数据基本一致,透射波阵面与水平方向夹角吻合程度优于Catherasoo结果;给出了反射波参数、界面偏折角度等,而Catherasoo方法不能计算这些参数.对于马赫数2.00的激波,通过比较不同介质密度比和界面倾角下的理论与数值结果,发现此方法对波系类型的识别比Catherasoo方法更为准确,且识别了波后强扰动追赶上激波形成三波结构的折射类型,后者无法识别此类型.上述结果表明,新方法具有良好的适用性,对波系类型的识别较已有方法的准确度更高,可获取更为丰富的波系结构...     

激波聚焦反射的实验和数值研究

 采用高速摄影技术和数值模拟方法,对入射激波在两种不同形状的抛物形反射器表面聚焦和反射的过程进行了研究,得到了激波聚焦反射过程的波系结构的实验阴影照片和数值计算结果,两者符合得很好。对激波聚焦形成的气体动力学焦点的特性进行了分析,结果表明,入射激波在两种反射器反射后聚焦所形成的气体动力学焦点均是由三波点在轴心处的会聚所导致的,气体动力学焦点位于相应的反射器壁面的几何焦点附近。不同的反射器中,激波聚焦前后的波系结构也不同,对较浅的反射器,入射激波反射前在反射器壁面形成了弓形激波,反射之后需要相对较长的时间完成聚焦,形成聚焦反射激波之后弓形激波仍未相交;对较深的反射器,入射激波反射后在更短的时间内聚焦,聚焦时弓形激波已经相交,聚焦反射激波之后的流场波系结构更加复杂。     

超音速气膜冷却及其受斜激波的影响

选取了一组接近超燃发动机内部实际工况的边界条件,对二维平行缝槽形式的超音速气膜冷却进行了数值计算,比较了吹风比、缝槽高度和隔板厚度等因素对壁面冷却效率的影响,并考虑了不同强度的斜激波入射气膜冷却边界层对壁面冷却效果和压力分布的影响.通过回归分析,不受激波影响的超音速气膜冷却的壁面冷却效率可以整理成统一的关联式.计算结果表明,斜激波入射将使入射点及其下游的壁面冷却效率降低,但是这个影响没有向上游传播.     

激波对超音速平板发汗冷却的影响

本文分析了激波影响下,不同激波强度,冷却气体和注入率对发汗冷却效果的影响。马赫数3的超音速主流遇到流道内的楔形激波发生器产生斜激波入射到多孔平板表面。楔块楔角φ分别为0°、4°、8°、12°模拟不同的入射激波强度,冷却介质分别为空气、甲烷与氢气。计算结果表明,激波使多孔区域出口表面静压上升,冷却流体流出受阻而破坏冷却效果,且冷却效率随激波强度增强而下降,随冷却气体分子量增大,冷却效率下降幅度减小;但在激波强度较强时,激波在多孔表面形成逆压梯度,迫使冷却流体流向入射点上游区域而使该区域冷却效果得到恢复。提高冷却剂注入率可以减弱激波对冷却效率的破坏作用,但使壁面温度不均匀性增加。     

理想磁流体中激波与矩形密度界面相互作用的数值研究

发展一套采用三阶WENO格式和混合GLM方法的理想磁流体数值方法,并对激波与矩形密度界面相互作用进行数值研究.通过圆极化阿尔芬波和旋转激波管问题对数值方法的稳定性和可靠性进行验证.在入射激波马赫数为10,界面内外气体密度比为10的情况,对比不同磁场中矩形密度界面的演变过程.结果表明,磁场能够减少界面上涡量的生成从而抑制界面不稳定性,并且磁场对界面的加速过程以及界面内外气体混合率有影响;而界面的存在将会使波后部分区域磁场增强;由于尖角的存在,矩形界面的发展与圆形界面不同.     

收缩楔腔内激波汇聚诱导点火实验

文章提出了一种采用圆柱形汇聚激波实现可燃气体点火特性研究的新方法.通过采用激波动力学理论合理地设计壁面型线, 将激波管中产生的平面运动激波近乎连续地转变为扇形区内圆柱形汇聚激波.以氢氧预混气体为考察对象, 开展了相关激波管实验, 实现了可控圆柱面激波汇聚诱导点火.实验发现两种点火现象:强点火和弱点火.在强点火过程中, 点火由入射激波直接诱导产生; 而在弱点火过程中, 点火则是在波后气流经历热压缩过程后发生.      

入射和反射激波诱导重气泡变形和失稳的三维数值研究

 采用大涡模拟方法,对入射激波及其反射激波诱导球形重气泡的变形失稳过程进行了三维数值模拟,利用已有实验验证了计算模型的可靠性,重点考察了反射激波与已经失稳的气泡界面的再次作用,讨论了涡环的形成及其三维失稳的过程。研究结果显示：入射和反射激波与球形重气泡作用产生斜压效应,会在流场中产生旋转方向截然相反的多个涡环;反射激波诱导的涡环具有较小的强度,故更加容易失稳,甚至能完全形成具有流向涡量的复杂小尺度涡结构。     

粘性物质中正激波稳定性分析

用线性稳定性理论,分析粘性物质中的正激波稳定性问题.粘性物质中任意强度的一维激波,其稳定性问题可归结为处理复数范围内的特征值问题,该特征值问题由两个一阶常微分方程及一个二阶常微分方程构成.这些常微分方程的系数依赖于流动的基本流场的物理量及其梯度.所获得的特征值问题由一个四阶精度的有限差分离散求解.分析考虑物质粘性的金属铝中的正激波稳定性,可以看出,正激波运动是稳定的,并且激波速度对波前和波后的小扰动量的衰减有相反的作用,而物质粘性有致稳的作用.     

冲击波在空气和熔石英元件界面传播的超快诊断

基于球面波在弹性介质界面上的反射和折射特性的基本理论,采用超快时间分辨的光学诊断技术,研究了532nm纳秒激光辐照熔石英元件表面产生的冲击波在空气和样品界面的传播特性,获得了冲击波在材料内部传输以及在空气与样品界面反射的时间分辨图像。结果表明:激光脉冲与材料作用在前后表面产生了向体内传输的冲击波,且产生的冲击波在玻璃与空气界面处反射为两个波,即反射波和反射剪切波;反射波和反射剪切波的强度与入射冲击波的入射角有关。     

并联铁电体冲击去极化输出电流实验和计算

 为分析不同数量铁电体的输出电流特征，进行了多片并联PZT95/5铁电体冲击去极化实验，测量了电阻负载下铁电体的输出电流；建立了并联铁电体输出电流的计算模型，计算了并联铁电体的输出电流。通过比较计算和实验结果，分析了铁电体数量和冲击波状态对铁电体输出电流的影响。结果表明：若进入铁电体的冲击波为平面波，铁电体输出电流的波形接近方波，电流峰值随铁电体并联数量的增加成比例增加，冲击波波形和侧向稀疏波对铁电体输出电流波形和去极化率有很大的影响。     

Anisotropic sound and shock waves in dipolar Bose-Einstein condensate

We study the propagation of anisotropic sound and shock waves in dipolar Bose-Einstein condensate in three dimensions (3D) as well as in quasi-two (2D, disk shape) and quasi-one (1D, cigar shape) dimensions using the mean-field approach. In 3D, the propagation of sound and shock waves are distinct in directions parallel and perpendicular to dipole axis with the appearance of instability above a critical value corresponding to attraction. Similar instability appears in 1D and not in 2D. The numerical anisotropic Mach angle agrees with theoretical prediction. The numerical sound velocity in all cases agrees with that calculated from Bogoliubov theory. A movie of the anisotropic wave propagation in a dipolar condensate is made available as supplementary material.     

Delta shock waves in conservation laws with impulsive moving source: some results obtained by multiplying distributions

The present paper concerns the study of a Riemann problem for the conservation law u_t + [? (u)]_x = kδ (x ? vt) where x, t, k, v and u = u(x,t) are real numbers. We consider ? an entire function taking real values on the real axis and δ stands for the Dirac measure. Within a convenient space of distributions we will explicitly see the possible emergence of waves with the shape of shock waves, delta waves and delta shock waves. For this purpose, we define a rigorous concept of a solution which extends both the classical solution concept and a weak solution concept. All this framework is developed in the setting of a distributional product that is not constructed by approximation. We include the main ideas of this product for the reader’s convenience. Recall that delta shock waves are relevant physical phenomena which may be interpreted as processes of concentration of mass or even as processes of formation of galaxies in the universe.     

Measurement of gas flow parameters by the laser produced breakdown method

A method of measuring gas flow parameters, based on testing the gas flow by shock waves, is described. The shock waves in the examined flow were produced by a pulsed laser beam focused on to the surface of an object. To illustrate the method, experimental data are presented dealing with gas parameters on the surface of a free-flying object.     

Colliding Wave Solutions in a Symmetric Non-Metric Theory

A method is given to generate the non-linear interaction (collision) of linearly polarized gravity coupled torsion waves in a non-metric theory. Explicit examples are given in which strong mutual focussing of gravitational waves containing impulsive and shock components coupled with torsion waves does not result in a curvature singularity. However, the collision of purely torsion waves displays a curvature singularity in the region of interaction.