二维激波与剪切层相互作用的直接数值模拟研究

采用五阶weighed esseritially non-oscillatory (WENO) 格式和三阶total variation diminishing (TVD) Runge-Kutta 格式, 通过求解二维非定常Navier-Stokes 方程, 直接数值模拟了激波与剪切层相互作用, 目的在于揭示激波与剪切层相互作用过程中噪声产生的机理. 研究发现:(1) 当入射激波穿过剪切层时, 剪切层中心位置向下层区域偏移;(2) 入射激波穿过剪切层产生小激波, 在小激波与剪切层接触点处产生声波并向外辐射;(3) 反射激波穿过剪切层后形成了分段弧状激波;(4) 当反射激波穿过剪切层时, 激波在鞍点处泄漏并向外辐射声波, 这是一种激波泄漏机制.     

激波诱导下的气体界面不稳定性实验研究

界面不稳定性, 特别是Richtmyer–Meshkov (RM) 不稳定性, 是流体
力学中一项重要的研究内容, 无论在学术研究领域还是工程应用领域都有着
重要的研究价值和应用背景. RM 不稳定性问题自提出以来, 得到了学术界
广泛的关注, 其研究无论是在实验方法、数值模拟还是在理论分析方面都取
得了很大的进展. 在激波管中开展激波与界面相互作用的实验研究, 即研究
界面初始扰动在激波诱导下的演化规律, 是目前研究RM 不稳定性的重要手
段. RM 不稳定性实验研究包括3 个部分, 分别是激波的产生、界面的形成
以及流场的观测. 综述了RM 不稳定性的实验研究进展, 并针对目前研究的
局限性提出了RM 不稳定性今后实验研究的重点和方向: 汇聚激波作用下界
面不稳定性的发展规律; 激波冲击下多种形状及大振幅界面的演化机理; 三
维界面的RM 不稳定性发展规律; 可压缩湍流的形成与混合机理.      

Numerical modeling of the initial stage of the generation of unsteady vortices from sharp corner in plane compressible flow

ＮＵＭＥＲＩＣＡＬＭＯＤＥＬＩＮＧＯＦＴＨＥＩＮＩＴＩＡＬＳＴＡＧＥＯＦＴＨＥＧＥＮＥＲＡＴＩＯＮＯＦＵＮＳＴＥＡＤＹＶＯＲＴＩＣＥＳＦＲＯＭＳＨＡＲＰＣＯＲＮＥＲＩＮＰＬＡＮＥＣＯＭＰＲＥＳＳＩＢＬＥＦＬＯＷＨｕａｎｇＤｕｎ（黄敦）（Ｄｅｐａｒｔ．...     

Shock wave deformation in shock-vortex interactions

An initially planar shock wave can undergo significant distortion to its shape along with changes in its strength during the period of its interaction with a compressible vortex. This phenomenon is studied by numerically simulating the shock wave-vortex interaction with a high resolution shock-capturing scheme. Incident shock waves of various Mach numbers are made to interact with a compressible vortex and the dependence of the shock wave distortion on the strength of the incident shock wave is studied in detail. It is known that the type of complex shock structure formed in the later stages of a compressible vortex-shock wave interaction is dependent on the Mach number of the incident shock wave. A simple physical model based on the principle of shock wave reflection is proposed to explain this complex shock structure formation and its dependence on the relative strengths of the interacting vortex and shock wave. Received July 28, 1997 / Accepted November 17, 1997     

A lattice boltzmann model for compressible perfect gas

A new lattice Boltzmann model for compressible perfect gas is proposed. The numerical example shows that it can be used to simulate shock wave and contact discontinuity. The results are comparable with those obtained by traditional methods. The ratio of specific heats γ may be chosen according to the requirement of problems. The project supported by the National Natural Science Foundation of China     

可压缩流向涡与反向运动激波相互作用的实验

对可压缩流向涡与反向运动激波相互作用的现象进行了实验研究．实验在９４ｍｍ×９４ｍｍ的方截面激波管中进行．在实验段上游安装了一个有限翼展平直机翼．当入射激波通过机翼后，波后２区气流在模型翼尖诱导出一条流向涡．入射激波在激波管端壁反射后，形成的反射激波在观察窗处和流向涡发生作用．实验中拍摄了激波与流向涡作用全过程的纹影照片，观察到了一些和定常激波与旋涡相互作用不同的现象，并与数值计算结果进行了初步比较     

Shock wave propagation and interactions

激波的传播特性既取决于激波的产生条件,也与所处的传播环境密切相关.驱动条件、几何边界、介质的物理化学属性等发生变化时,都会引起激波传播特性的改变,而激波的变化反过来又会对其波及的流场产生影响.尽管激波传播及其干扰现象广泛存在于自然界和人类科技活动之中,其复杂机理的认识、规律的描述乃至应用潜力的挖掘仍有漫长的路要走.本文根据气体中激波传播和干扰现象以及与之相关的理论描述特征,在对激波传播以及反射、折射等基本现象进行简要阐述的基础上,重点围绕目前的热点问题,包括激波/激波干扰、激波/边界层干扰、激波与湍流作用、激波的聚焦与点火以及激波作用下气体界面不稳定性等研究进行了介绍和讨论,旨在对近年来该领域的进展及获得的成果做一个概述和归纳,期望对将来的深入研究有一个鉴借意义.     

Stability of Contact Discontinuities for the 1-D Compressible Navier-Stokes Equations

In this paper, we study the large-time asymptotic behavior of solutions of the one-dimensional compressible Navier-Stokes system toward a contact discontinuity, which is one of the basic wave patterns for the compressible Euler equations. It is proved that such a weak contact discontinuity is a metastable wave pattern, in the sense introduced in [24], for the 1-D compressible Navier-Stokes system for polytropic fluid by showing that a viscous contact wave, which approximates the contact discontinuity on any finite-time interval for small heat conduction and then runs away from it for large time, is nonlinearly stable with a uniform convergence rate provided that the initial excess mass is zero. This result is proved by an elaborate combination of elementary energy estimates with a weighted characteristic energy estimate, which makes full use of the underlying structure of the viscous contact wave.     

Vorticity produced by shock wave diffraction

Numerical simulations with a monotonicity preserving flow solver have been performed to study shock diffraction phenomena and shock wave generated vorticity. The computations were performed using the conservative Finite Element Method-Flux Corrected Transport (FEM-FCT) scheme, which has been shown to have an excellent predictive capability for various compressible flows with both strong and weak shocks. An adaptive unstructured methodology based on adapting to high density and entropy gradients was used in conjunction with a conservative shock-capturing scheme to adequately resolve strong and weak flowfield gradients. The chief interest was the formation of vorticity arising from shock wave propagation over a sharp corner and the high accuracy and resolution of the interacting compressible wave features. Numerical simulations were compared with previous experimental results and exhibited remarkably good agreement in terms of compressible wave propagation, as well as vorticity development and transport. The computations also allowed insight into the fundamental fluid dynamics, specifically shock diffraction, vortex convection and shock-vortex interactions.     

ANTI-BLAST ANALYSIS OF CELLULAR SACRIFICIAL CLADDING

运用一维冲击波模型和三维细观有限元模型分析了多胞牺牲层的抗爆炸行为. 基于刚性-塑性硬化(R-PH)的多胞材料模型,建立了一维冲击波模型,得到了多胞牺牲层中冲击波传播的控制方程. 揭示了冲击波在多胞牺牲层中的传播特性,并阐述了附加质量和爆炸载荷强度两个参数对牺牲层设计的重要影响. 比较了基于刚性-理想塑性-锁定(R-PP-L) 模型和基于刚性-塑性硬化(R-PH) 模型的多胞牺牲层的结构设计,指出了两种模型的适用范围. 通过基于三维Voronoi 技术的细观有限元方法验证了基于R-PH 模型的多胞牺牲层结构的设计准则.     

激波驱动的气固两相流的动态压力测量

利用水平圆柱形激波管对激波驱动的可压缩性气固两相流进行了试验研究.利用压电式压力传感器、电荷放大器、示波器及计算机组成的压力信号测试系统, 对激波与颗粒作用前后的气相参数进行测量及分析. 试验中测得了激波在管中的传播速度, 波后气流的压力, 反射激波、透射激波的压力和速度等. 分别考察颗粒、装载比、驱动气源以及入射激波马赫数等因素的差异对气相参数的影响.试验结果表明: 激波与颗粒群相互作用时, 会产生反射激波和透射激波,其强度与驱动气源、颗粒大小、颗粒装载比等参数有关;激波衰减率随着装载比、马赫数的增大而减小. 研究指出,在颗粒群被激波加速的初始阶段, 颗粒间的弹性碰撞起着重要的作用.      

A STUDY ON JET PHENOMENON OF R22 GAS CYLINDER UNDER THE IMPACT OF SHOCK WAVE

重物质气泡在平面激波作用下出现射流是一个非常有趣的物理现象, 它与气体属性、激波强度密切相关, 气泡内外波系结构复杂, 值得进行深入探讨. 论文采用水平集方法结合真实虚拟流体方法捕捉物质界面, 对R₂₂圆柱气泡在平面激波作用下的射流现象进行了数值模拟, 详细阐述了波系结构在气泡内外的演变过程, 总结了导致气泡射流现象形成的典型波系结构特征, 从波系结构角度解释了射流形成的原因.论文也对不同强度的激波冲击R₂₂圆柱气泡的情况进行了归纳, 发现激波强度对射流现象的出现以及射流结构的尺寸特征有重要影响, 得出了 R₂₂圆柱气泡能否出现射流结构的临界激波强度, 以期对实验结果起到指导与验证作用.     

Dynamics of stress wave propagation in a chain of photoelastic discs impacted by a planar shock wave; Part I, experimental investigation

The propagation of stress waves through a chain of discs has been studied experimentally. Optically transparent 20-mm diameter discs, made of epoxy, were loaded dynamically by head-on collision with an incident planar shock wave. The loading was done in a vertical shock tube. The head-on collision between the punch-plate, placed on top of the chain of discs, and the incident shock wave resulted in a head-on reflected shock wave inducing behind it a fairly uniform step-wise pressure pulse having duration of about 6 ms. The recorded fringe patterns of the stress field, in the discs-chain, show that the input pressure pulse was broken into several oscillating cycles. The back and forth bouncing of stress waves gave rise to two different modes of the contact stress oscillations, which continued until the overall stress reaches equilibrium with the input conditions. The registered propagation velocity of the stress wave was significantly lower than the appropriate speed of sound in the material from which the discs were made.      

Transient phenomena in one-dimensional compressible gas–particle flows

The transient behavior of compressible gas– particle flows produced in shock tubes with particle-laden driver section is studied. Particular attention is focused on the time scales with which the solution approaches the equilibrium state. Theoretical estimates indicate that the gas and particle contact surfaces equilibrate first, followed by the shock wave, and finally by the expansion fan. The estimates are in good agreement with numerical simulations. The simulations also show that the approach to equilibrium condition of the shock speed is non-monotonic (monotonic) if the mass fraction of particles initially located in the driver section is below (above) a particle-diameter dependent critical value. For the speed of the particle contact surface, the reverse trends are observed.      

Shock wave diffraction by a square cavity filled with dusty gas

A simple two-dimensional square cavity model is used to study shock attenuating effects of dust suspension in air. The GRP scheme for compressible flows was extended to simulate the fluid dynamics of dilute dust suspensions, employing the conventional two-phase approximation. A planar shock of constant intensity propagated in pure air over flat ground and diffracted into a square cavity filled with a dusty quiescent suspension. Shock intensities were and , dust loading ratios were and , and particle diameters were and {\rm \mu}$m. It was found that the diffraction patterns in the cavity were decisively attenuated by the dust suspension, particularly for the higher loading ratio. The particle size has a pronounced effect on the flow and wave pattern developed inside the cavity. Wall pressure histories were recorded for each of the three cavity walls, showing a clear attenuating effect of the dust suspension. Received 15 November 1999 / Accepted 25 October 2000     

A Free-Lagrange method for unsteady compressible flow: simulation of a confined cylindrical blast wave

A Free-Lagrange numerical procedure for the simulation of two-dimensional inviscid compressible flow is described in detail. The unsteady Euler equations are solved on an unstructured Lagrangian grid based on a density-weighted Voronoi mesh. The flow solver is of the Godunov type, utilising either the HLLE (2 wave) approximate Riemann solver or the more recent HLLC (3 wave) variant, each adapted to the Lagrangian frame. Within each mesh cell, conserved properties are treated as piece-wise linear, and a slope limiter of the MUSCL type is used to give non-oscillatory behaviour with nominal second order accuracy in space. The solver is first order accurate in time. Modifications to the slope limiter to minimise grid and coordinate dependent effects are described. The performances of the HLLE and HLLC solvers are compared for two test problems; a one-dimensional shock tube and a two-dimensional blast wave confined within a rigid cylinder. The blast wave is initiated by impulsive heating of a gas column whose centreline is parallel to, and one half of the cylinder radius from, the axis of the cylinder. For the shock tube problem, both solvers predict shock and expansion waves in good agreement with theory. For the HLLE solver, contact resolution is poor, especially in the blast wave problem. The HLLC solver achieves near-exact contact capture in both problems. Received May 25, 1995 / Accepted September 11, 1995     

分数阶理论下实心球体受热冲击作用的渐进分析

基于分数阶广义热弹性理论,针对实心球体在外表面受均匀热冲击作用下的一维广义热弹性问题进行研究分析. 利用热冲击的瞬时特征,借助于Laplace 正、反变换技术及柱函数的渐进性质,推导了热冲击作用周期内位移场、温度场和应力场的渐进表达式. 通过计算,得到了不同传热能力下受热冲击作用时热波、热弹性的传播规律以及位移场、温度场及应力场的分布规律. 结果表明:分数阶参数取值的不同,热波、热弹性波的传播以及各物理场的分布均有所不同,分数阶参数可视为延迟时间的影响因子,通过改变延迟效应对热弹性行为的影响来改变热冲击的作用效果.      

激波作用下三维涡的形成及其界面捕捉计算

针对三维多介质可压缩流体,给出了可压缩多介质流体三维高精度数值计算方法,以及界面捕捉方程和带重新初始化的三维LevelSet方法,对初始压力间断和密度间断条件形成激波、接触间断以及稀疏波的三维复杂流场相互作用情况进行数值计算,给出流场中涡的形成过程和界面位置。并对计算方法进行理论验证。     

A high‐resolution scheme for compressible multicomponent flows with shock waves

A simple methodology for a high‐resolution scheme to be applied to compressible multicomponent flows with shock waves is investigated. The method is intended for use with direct numerical simulation or large eddy simulation of compressible multicomponent flows. The method dynamically adds non‐linear artificial diffusivity locally in space to capture different types of discontinuities such as a shock wave, contact surface or material interface while a high‐order compact differencing scheme resolves a broad range of scales in flows. The method is successfully applied to several one‐dimensional and two‐dimensional compressible multicomponent flow problems with shock waves. The results are in good agreement with experiments and earlier computations qualitatively and quantitatively. The method captures unsteady shock and material discontinuities without significant spurious oscillations if initial start‐up errors are properly avoided. Comparisons between the present numerical scheme and high‐order weighted essentially non‐oscillatory (WENO) schemes illustrate the advantage of the present method for resolving a broad range of scales of turbulence while capturing shock waves and material interfaces. Also the present method is expected to require less computational cost than popular high‐order upwind‐biased schemes such as WENO schemes. The mass conservation for each species is satisfied due to the strong conservation form of governing equations employed in the method. Copyright © 2010 John Wiley & Sons, Ltd.     

Uniaxial wave propagation in a nonlinear thermoviscoelastic medium with temperature dependent properties

The problem of finite wave propagation in a nonlinearly thermoviscoelastic thin rod whose viscoelastic properties are temperature dependent is considered. The rod is subjected to mechanical or thermal time-dependent loading. The coupled equations of motion and heat conduction are based on a constitutive theory of nonisothermal nonlinear viscoelasticity which is described by single-integral terms only. This theory is reformulated here for the uniaxial motion of a compressible rubbery material. The solution of the field equations is obtained by a numerical procedure which is developed for the present case and is able to handle successfully shock waves whenever they built up in the nonlinear material.