激波与SF_6梯形气柱相互作用的数值模拟

基于大涡模拟,结合五阶加权基本无振荡格式以及沉浸边界法对平面入射激波与两种SF_6梯形重气柱的相互作用过程进行了数值模拟,数值结果清晰地显示了激波诱导Richtmyer-Meshkov不稳定性所导致的两种梯形重气柱的变形过程,详细分析了入射激波在两种梯形重气柱界面发生反射、折射、绕射以及折射激波与透射激波在气柱内部来回反射的过程,并研究了该过程中所产生的复杂波系结构,对两种梯形气柱变形过程中与周围空气的混合过程进行了分析;通过记录气柱界面四个特征尺寸随时间的变化对两种梯形气柱界面的不同演化过程进行了定量分析。     

激波冲击R22重气柱所导致的射流与混合研究

基于大涡模拟, 结合五阶加权基本无振荡格式与沉浸边界法对激波自左向右与R22重气柱作用过程进行了数值模拟. 数值结果清晰地显示了激波诱导Richtmyer-Meshkov不稳定性所导致的重气柱变形过程, 并与Haas 和 Sturtevant 的实验结果符合. 另外, 结果还揭示了入射激波在气柱内右侧边界发生聚焦并诱导射流的过程, 以及在Kelvin-Helmhotz 次不稳定性作用下两个主涡滑移层形成次级涡的过程, 并分析了气柱变形过程中与周围空气的混合机理. 最后, 通过改变反射距离对反射激波与不同变形阶段的气柱的再次作用过程进行了研究. 结果表明: 当激波反射距离较长时, 反射激波与充分变形后的气柱作用, 使其在流向方向上进一步被压缩; 而当激波反射距离较短时, 反射激波会在气柱内发生马赫反射, 两个三波点附近产生两个高压区, 当其传播至气柱左侧边界时对气柱边界造成冲击加速, 诱导两道向左传播的反向射流. 关键词： Richtmyer-Meshkov不稳定性 R22重气柱 反射激波 射流     

激波诱导双层气柱演化的偏心效应研究

实验与数值研究了平面激波诱导下双层气柱的演化规律.利用肥皂膜技术生成了3种双层气柱,通过固定内、外层气柱半径,改变内层气柱在流向方向上的位置,研究了双层气柱偏心效应对流场演化的影响.结果表明,当内层气柱向上游偏移时,外层气柱的上游界面在后期会产生朝向上游的射流;当内层气柱向下游偏移时,下游两道界面会较早地耦合在一起...     

偏心对汇聚激波诱导的RM不稳定性影响的数值研究

数值研究汇聚激波与四种形状（圆形、小振幅单模、大振幅单模和正方形）的二维气柱界面相互作用,激波汇聚中心与界面同心和不同心（即偏心）时Richtmyer-Meshkov（RM）不稳定性的发展规律,重点考察界面中心的压力及混合区面积在两种情况下随时间的变化.数值方法使用VAS2D程序,该方法采用有限体积法结合网格自适应技术,能够达到时间和空间的二阶精度.结果表明,偏心情况下RM不稳定性是其在同心情况下的扰动和偏心小扰动叠加的结果.在本文采用的偏心程度下（20%）,偏心对于圆形无扰动界面发展的影响主要表现在后期界面出现微小扰动结构;而对于单模和正方形这种原本有扰动的界面,偏心使扰动结构呈现不对称及扭曲,同时也影响了界面中心压力和混合区面积,因而加剧了不稳定性的发展.     

激波诱导的液柱气液界面RM不稳定性的研究

为了进-步研究气/液界面上的Richtmyer-Meshkov(RM)不稳定性现象,在水平激波管中进行了激波与单、双排液柱相互作用的实验.通过使用驱动气体为氮气、被驱动气体为空气,获得了马赫数为1.10与1.25的激波;通过更换不同尺寸的试件获得了初始直径不同的液柱。使用高速摄影仪记录了激波作用之后液柱失稳变形的详细过程,测量分析了不稳定性发展过程中,尖钉高度、气泡深度及混合区域宽度随时间的发展规律,并与理论分析中的公式进行了比较。     

弱激波冲击无膜重气柱和气帘界面的实验研究

研制了一个改进的激波管设备,对马赫数为1.2的弱激波冲击作用下空气中SF6气柱和气帘界面的演变过程进行了初步的实验研究。通过设计激波管实验段、烟雾发生器、气体箱、进气吸气系统和激波管尾段,控制混合气体中SF6的峰值浓度和初始气流速度,建立了稳定、可重复的无膜气柱和气帘初始界面形成技术。利用高速摄影技术,在水平面内观测了气柱和气帘的初始界面图像,沿垂直方向观测了界面RM(Richtmyer-Meshkov)不稳定性的演变过程。气柱演变图像显示了典型的对涡结构,气帘演变图像显示了早期的多蘑菇形结构和后期的相邻波长干扰效应。图像后处理表明,气柱的高度和宽度、气帘的宽度均随时间单调增加,且宽度比高度增加快得多。从二维涡量动力学方程出发,对图像中涡的演变过程进行了初步解释。     

二维气/液界面不稳定性数值模拟

 以多介质的体积分数方法和三阶PPM（Piecewise Parabolic Method）方法为基础,给出了适用于多介质流体动力学数值模拟的计算方法和程序MFPPM。利用MFPPM程序对在高压气体冲击作用下的气体/液体交界面的Richtmyer-Meshkov（RM）不稳定性及其引起的流体混合现象进行了数值模拟研究。主要研究在不同的初始扰动情况下流体混合区的发展,并细致研究了流体混合区的宽度、气泡和尖钉高度随时间的增长情况及不同初始扰动对它们的影响;同时还研究了网格尺度不同时混合区、气泡以及尖钉的构型和高度的增长情况。通过对计算结果的分析得出,流体混合区、气泡以及尖钉的发展与初始扰动有密切的关系,特别是在后期影响更为显著;混合区宽度的变化过程和尖钉相似,而气泡高度的变化基本上呈线性增长趋势,且受初始扰动的影响比较小,但是其构型却有明显差别;网格的影响也主要体现在对混合区、气泡和尖钉的构型上。     

激波作用不同椭圆氦气柱过程中流动混合研究

在激波与气柱相互作用问题中,压力与密度间断不平行产生的斜压涡量会引起流动的不稳定性,从而促进物质间的混合.本文基于双通量模型,结合五阶加权基本无振荡(WENO)格式,求解多组分二维Navier-Stokes方程,分析激波作用面积相同结构不同的椭圆气柱所致的流动和混合.数值结果清晰地显示了激波诱导Richtmyer-Meshkov不稳定性引起的气柱界面变形和波系演化.同时定量地从界面运动、界面结构参数变化(长度和高度)、气柱体积压缩率、环量及混合率等角度分析激波诱导的流动混合机制,研究椭圆几何构型对氦气混合过程的影响.结果表明,界面及相关参数的演化与气柱初始形状密切相关.当激波沿椭圆长轴作用于气柱时,气柱前端出现空气射流结构,且射流不断增长并渗透到下游界面,致使气柱分离成两个独立涡团,离心率越大,射流发展越快;同时激波作用气柱后在界面处产生不规则反射现象.圆形气柱界面演化与这种作用情形类似.当激波沿椭圆短轴作用于气柱时,界面上游出现类平面结构,随后平面上下缘处产生涡旋,主导流动发展,激波在界面作用产生规则反射,离心率越大,这些现象越明显.界面高度、长度、体积压缩率也因此有所差异.对界面演化、环量和混合率的综合分析表明,激波沿长轴作用于气柱且离心率较大时,流动发展较快,不稳定性导致的流动越复杂,越有利于氦气与环境介质的混合.     

激波与SF6球形气泡相互作用的数值研究

本文基于大涡模拟方法, 采用高阶精度格式对平面入射激波以及不同反射距离条件下的反射激波与SF₆重气泡相互作用过程进行了三维数值模拟. 数值结果清晰地显示了SF₆重气泡在激波作用下诱导Richtmyer-Meshkov不稳定性过程, 揭示了入射激波以及反射激波在气泡界面聚焦诱导射流的过程, 详细分析了不同反射距离条件下反射激波与SF₆重气泡作用过程及流场结构.     

平面弱激波加载下球形气泡演化的实验研究

 利用高速摄影机直接观测了激波加载下球形气泡的演化过程。采用白色的烟雾颗粒对气泡内透明的测试气体进行染色,直接观察到了气泡界面。通过对平面弱激波加载下轻气体气泡和重气体气泡进行研究,成功验证了轻气体气泡和重气体气泡演化过程中出现的典型界面结构,获得了弱激波加载轻气体气泡中背风涡环的环状几何结构。运用相关方法分析了轻气体气泡实验初期的流场分布,得到的结果与理论分析结果吻合得很好,对于此类实验的数据处理具有一定的启发意义。同时,实验中使用的烟雾颗粒法为今后在球形气泡实验中引入更为精确的平面激光诱导荧光技术（PLIF）、激光粒子图像测速技术（PIV）等实验测试系统提供了布撒示踪粒子的范例。     

Discrete Alfvn eigenmodes in international thermonuclear experimental reactor operations with negative magnetic shear

Discrete Alfvén eigenmodes in steady-state operation scenarios with negative magnetic shear in the international thermonuclear experimental reactor are investigated in this paper. These magnetohydrodynamic eigenmodes are trapped by the α -induced potential wells along the magnetic field line. Here α =-q²Rdβ/d r with q being the safety factor, β the ratio between plasma and magnetic pressures, and R the major radius, and r the minor radius. Due to negligible continuum damping via wave energy tunneling, these Alfvén eigenmodes could be readily destabilized by energetic particles.     

Electromagnetic interactions in a cluster of cylinders

We present a rigorous multiple-scattering method to calculate the interaction energies due to electromagnetic field fluctuations in a cluster of parallel cylinders. Various kinds of cluster structure with both isotropic and anisotropic materials are considered. It is shown that Casimir interactions are sensitive to the positions of the cylinders in the cluster or the structure of the system. The sign of the interaction energy cannot only be changed by tuning the parameters of anisotropic materials, it can also be caused by changing the cylinders configuration of the system.     

The refined inviscid stability condition and cellular instability of viscous shock waves

Combining the work of Serre and Zumbrun, Benzoni-Gavage, Serre, and Zumbrun, and Texier and Zumbrun, we propose as a mechanism for the onset of cellular instability of viscous shock and detonation waves in a finite-cross-section duct, the violation of the refined planar stability condition of Zumbrun-Serre, a viscous correction of the inviscid planar stability condition of Majda. More precisely, we show for a model problem involving flow in a rectangular duct with artificial periodic boundary conditions that transition to multidimensional instability through violation of the refined stability condition of planar viscous shock waves on the whole space generically implies for a duct of sufficiently large cross-section, a cascade of Hopf bifurcations involving more and more complicated cellular instabilities. The refined condition is numerically calculable as described by Benzoni-Gavage-Serre-Zumbrun.     

Effect of external magnetic field on nonlinear propagation of cylindrical magnetoacoustic soliton waves in quantum plasma

The propagation of magnetoacoustic solitary waves is investigated in magnetized quantum plasma consisting of cold ions and hot electrons. By using the quantum magnetohydrodynamic (QMHD) model, the nonlinear characteristics of different features of solitary waves in an electron-ion quantum magnetoplasma are investigated. Magnetoacoustic solitary waves are stationary solutions of the equations composed of the nonlinear mass and momentum continuity, together with the Maxwell's equations. The important quantum-mechanical effects including the quantum statistical and diffraction are examined numerically on the profiles of the solitons. It is found that the non-dimensional characteristic of the quantum parameter plays a significant role in the formation of the solitons.     

磁场对螺旋波等离子体波和能量吸收影响的数值研究

本文考察在径向电子数密度呈抛物形分布的情况下,外加稳恒磁场,射频通过螺旋波天线在等离子体中激发电磁波的传播性质.采用线性扰动波假设,数值求解Maxwell方程组,得到80—800 G(1 G=10-4T)磁场条件下等离子体中径向电、磁场强度及能量沉积密度的分布情形.计算结果表明,磁场增大(80→800G)时,螺旋波受到的阻尼较小,可深入等离子体传播;Trivelpiece-Gould(TG)波受到的阻尼增大,在等离子体-真空边界处衰减增强;整体的能量吸收向边界集中.磁感应强度小于100 G时,TG波可深入主等离子体区传播,等离子体径向能量吸收相对均匀.     

同心旋转圆柱间隙流场中纤维取向的动态模拟

本文对两同心旋转圆柱间隙形成的流场以及处于流场中的纤维运动和取向进行了数值研究. 在贴体坐标网格下求解了流场控制方程, 得到了流场中的速度、压力等物理量. 研究了两同心圆柱同速反向旋转以及仅内层圆柱旋转这两种情况下的纤维运动和取向状态. 得到了处于这两种情况下的纤维在流场中从静止到开始运动和取向直至最终达到稳定状态的动态详细过程. 结果表明, 当两个圆柱同速反向旋转时, 纤维运动与取向也相应的呈现两层结构; 而仅内圆柱旋转时, 纤维运动与取向呈单层结构. 在两种情况下, 纤维均沿流线方向运动和取向. 讨论了纤维长径比对纤维取向的影响, 结果表明随着纤维长径比的增加, 纤维沿流线取向的取向度逐渐增强.     

Identification of the magnetic island structure in the HL-2A tokamak

The identification of the magnetic island structure in the HL-2A tokamak is presented. First, the perturbation current as a source for the perturbation flux can be determined by using Mirnov probe measurements. By superposing the perturbation flux and the equilibrium flux reconnected by equilibrium fitting, the structure and the width of the magnetic islands can be estimated. The method has been used in the HL-2A experiments.     

内孤立波沿缓坡地形传播特性的实验研究

以南中国海东北部海域底部缓坡地形为背景, 在大型重力式分层流水槽中模拟了下凹型内孤立波沿缓坡地形传播过程中的浅化、破碎、分裂等现象, 利用分层染色标识方法和多点组合探头阵列技术对内孤立波沿缓坡地形演化特征进行了定性分析和定量测量. 实验表明: 浅化效应使内孤立波传播速度减小, 对大振幅内孤立波具有抑制作用, 对小振幅波具有放大效应; 浅化效应可导致内孤立波的剪切失稳及破碎, 还可导致大振幅内孤立波的分裂. 利用Miles稳定性理论可定性描述内孤立波沿缓坡地形传播时发生不稳定状态的位置, 实验结果与理论分析相符合. 关键词： 分层流 缓坡地形 内孤立波 不稳定性