激波管实验和果冻实验界面不稳定性数值计算

运用算子分裂技术,增加二阶空间中心差方法和两步Rung-Kutta时间推进方法计算动力学粘性以及热流部分对流场的影响,将可压缩多介质流体动力学高精度欧拉计算方法MFPPM发展到适用于NS(Nayier Stokes)方程的可压缩多介质粘性流体计算方法MVPPM.通过与界面不稳定性实验结果的比较,来检验计算方法的正确性,并验证计算程序的有效性.主要包括一个激波管实验和两个果冻实验,即英国AWE(Atomic Weapons Establishment)激波管实验和LLNL实验室的果冻环实验以及中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室进行的爆轰驱动下单层果冻界面不稳定性实验.研究结果表明:数值模拟结果与实验测量结果以及对应时刻的实验图象均吻合较好.     

用激波管研究超音速气固两相流

在激波管里进行了可压缩性气固两相流的实验研究。测量了激波通过颗粒群时的压力的衰减过程。用纹影仪拍摄了激波与颗粒群相互干涉的照片。试验了颗粒群的不同构造对压力衰减的影响。指出了激波反射、聚焦等非线性气动因素是可压缩性气固两相流的关键问题。     

Impulse Loading of Plates using a Diverging Shock Tube

Experimental Mechanics - Impulsive blast loading studies typically require facilities with sub-millisecond decay time simulation capability. Using a conventional shock tube for such studies,...     

Focusing of Strong Shocks in an Annular Shock Tube

Focusing of strong shock waves in a gas-filled thin convergence chamber with various forms of the reflector boundary is investigated experimentally and numerically. The convergence chamber is mounted at the end of the horizontal co-axial shock tube. The construction of the convergence chamber allows the assembly of the outer chamber boundaries of various shapes. Boundaries with three different shapes have been used in the present investigation—a circle, an octagon and a smooth pentagon. The shock tube in the current study was able to produce annular shocks with the initial Mach number in the range M _s = 2.3 − 3.6. The influence of the shape of the boundary on the shape and properties of the converging and reflected shock waves in the chamber has then been investigated both experimentally and numerically. It was found that the form of the converging shock is initially governed by the shape of the reflector and the nonlinear interaction between the shape of the shock and velocity of shock propagation. Very close to the center of convergence the shock obtains a square-like form in case of a circular and octagonal reflector boundary. This is believed to stem from the instability of the converging shock front triggered by the disturbances in the flow field. The outgoing, reflected shocks were also observed to be influenced by the shape of the boundary through the flow ahead as created by the converging shocks.     

Analytical Solution for Isothermal Flow in a Shock Tube Containing Rigid Granular Material

Analytical solution of shock wave propagation in pure gas in a shock tube is usually addressed in gas dynamics. However, such a solution for granular media is complex due to the inclusion of parameters relating to particles configuration within the medium, which affect the balance equations. In this article, an analytical solution for isothermal shock wave propagation in an isotropic homogenous rigid granular material is presented, and a closed-form solution is obtained for the case of weak shock waves. Fluid mass and momentum equations are first written in wave and (mathematical) non-conservation forms. Afterwards by redefining the sound speed of the gas flowing inside the pores, an analytical solution is obtained using the classical method of characteristics, followed by Taylor’s series expansion based on the assumption of weak flow which finally led to explicit functions for velocity, density and pressure. The solution enables plotting gas velocity, density and pressure variations in the porous medium, which is of high interest in the design of granular shock isolators.     

High-Precision Digital Image Correlation for Investigation of Fluid-Structure Interactions in a Shock Tube

Experimental Mechanics - Background: Structural response measurements are challenging in aerodynamic testing environments due to high-speed requirements, facility vibrations, and the desire for...     

Gas and Particle Dynamics of a Contoured Shock Tube for Pre-clinical Microparticle Drug Delivery

We investigate the gas-particle dynamics of a device designed for biological pre-clinical experiments. The device uses transonic/supersonic gas flow to accelerate microparticles such that they penetrate the outer skin layers. By using a shock tube coupled to a correctly expanded nozzle, a quasi-one-dimensional, quasi-steady flow (QSF) is produced to uniformly accelerate the microparticles. The system utilises a microparticle “cassette” (a diaphragm sealed container) that incorporates a jet mixing mechanism to stir the particles prior to diaphragm rupture. Pressure measurements reveal that a QSF exit period – suitable for uniformly accelerating microparticles – exists between 155 and 220 mus after diaphragm rupture. Immediately preceding the QSF period, a starting process secondary shock was shown to form with its (x,t) trajectory comparing well to theoretical estimates. To characterise the microparticle, flow particle image velocimetry experiments were conducted at the nozzle exit, using particle payloads with varying diameter (2.7–48 μm), density (600–16,800 kg/m³) and mass (0.25–10 mg). The resultant microparticle velocities were temporally uniform. The experiments also show that the starting process does not significantly influence the microparticle nozzle exit velocities. The velocity distribution across the nozzle exit was also uniform for the majority of microparticle types tested. For payload masses typically used in pre-clinical drug and vaccine applications (≤ 1 mg), it was demonstrated that payload scaling does not affect the microparticle exit velocities. These characteristics show that the microparticle exit conditions are well controlled and are in agreement with ideal theory. These features combined with an attention to the practical requirements of a pre-clinical system make the device suitable for investigating microparticle penetration into the skin for drug delivery.     

激波管实验研究非均匀流场R-M不稳定性

利用激波管装置及马赫数为1.27的弱入射激波实验研究了SF6非均匀流场的R-M不稳定性。Air/SF6初始正弦界面由厚度为0.5μm的薄膜相隔得到,由阴影方法记录界面演化过程。实验结果表明:由于不稳定性,重流体(SF6)向轻流体(Air)演化成"尖钉"结构,而轻流体演化为"气泡"结构;由于界面切向速度差的Kelvin-Helmholtz不稳定性,"尖钉"头部翻转成蘑菇头形状;由于流场密度分布不均,低密度区流场扰动增长较快,扰动振幅发展的实验结果与PPM数值计算的结果较吻合。     

Shock Wave Pressure and Velocity Measuring Using a Novel Optic Sensor in a Newly Designed Diaphragm-less Shock Tube

Experimental Techniques - This paper suggests a new optical system for measuring shock wave characteristic that is installed in a shock tube. This shock tube has an innovative pneumatic rapid...     

A High Vacuum Fracture Facility for Molecular Interactions

This paper describes the development of a high vacuum fracture rig for delaminating functionalized silicon surfaces. The main focus here is on examining molecular interactions by functionalizing the silicon surfaces with carboxyl and diamine terminated self-assembled monolayers (SAMs). However, many other interactions can be considered. The crack front location and normal crack opening displacements (NCOD) are measured by infrared crack opening interferometry (IR-COI). This allows the normal traction-separation relation of the interactions to be determined. Some mixed-mode fracture experiments were conducted on silicon/carboxyl/diamine/silicon sandwich specimens in ambient and high vacuum. Interesting differences in behavior were noted.     

一种激波驱动的新型固粒冲蚀试验系统

针对目前大部分固体颗粒冲蚀试验系统存在反弹颗粒干扰、颗粒质量计算误差大等问题,采用激波驱动气固两相流的方式,并结合动态压力传感器、高速摄影仪对激波及颗粒测速的方法,设计了1种新型冲蚀磨损试验系统,具有结构简单、操作方便、颗粒质量计算准确等优点.该装置最大激波马赫数可达2.3,能将颗粒加速到200 m/s以上,适用不同颗粒、材料及表面处理试件进行15～90°冲击角下的磨损特性试验.通过SiO2颗粒对低碳钢的冲蚀试验,证实了该装置测得数据与经典试验数据规律一致.因此,本试验系统有望作为1种新的冲蚀磨损试验方法,应用于耐固粒冲蚀的优化设计和寿命预测等领域.     

激波风洞内超燃冲压发动机三面压缩进气道流场实验观测

主要进行了超燃冲压发动机三面压缩进气道的实验观测。利用来流马赫数4.5的直通式激波风洞,考察了三组具有不同压缩角度的进气道模型内部的流场情况。实验观测手段为油流法、丝线法和高速纹影,同时,辅以数值模拟以有助于流场细节分析。纹影照片展示了进气道内部以激波边界层相互作用为主要影响因素的流场复杂结构,数值模拟也显示了相近的结果。油流技术与丝线法显示了近壁面处的流动图像,照片中可见激波、分离线、再附线等分界线位置。根据实验结果,可以推测唇口激波与进气道内边界层的相互作用及其引起的壁面分离是影响进气道内流动的主要因素。同时,尝试了利用抽吸方法减弱激波与边界层相互作用诱发的壁面流动分离,并取得一定结果。     

Measuring Dynamic Pressure in Shock Tube and Musical Instruments with WIKA A-10 Sensor*

International Applied Mechanics - Devices used in scientific research related to musical performance for recording dynamic mechanical processes are described. Recorder LMS SCADAS Mobile, dynamic...     

A Measurement Technique for Shock Wave-Loaded Structures and Its Applications

In engineering problems it is necessary to predict the deformations of structural elements subjected to shock waves. In the literature a wide range of structural theories, constitutive equations and simulation techniques is available in order to simulate the occurring deformations. However, an objective statement about the accuracy of calculated structural deformations is only possible by comparing these results to experiments. In the present work a measurement technique with shock tubes is introduced which was especially developed to measure fast deections of plates during the impulse duration. This technique provides a possibility to validate and to improve constitutive and structural theories. Furthermore, very precise measurements can be performed in order to observe limit states and buckling of repeatedly loaded plates. These applications are shown in this study.     

Shock Hugoniot for Biological Materials

This paper presents shock Hugoniot compression data for several bio-related materials by using flat plate impact experiments. Shock pressure covered in this study ranges at least up to 1 GPa. It is emphasized here that spatial and temporal uniformity of pressure distribution behind a shock wave front is very important and it can be realized in these materials by the impact method, and their precise shock propagation characteristics have been obtained by the application of the procedure developed previously by our group. Hugoniot measurements for different and systematic data for various samples are compared with shock Hugoniot curve for water. Samples used in the experiment include water gel of gelatin, NaCl aqueous solution, and finally chicken breast meat. Several samples with different initial density were used for gelatin and NaCl solution. Shock Hugoniot function for all the specimens tested can be summarized as u _s =A+B u _p , B ~ 2Value of the intercept of the relationship, A, which has the meaning of the sound velocity, is apparently dependent on the material and ambient temperature. Physical meanings of the obtained results have been discussed.     

Shock conditions for hypoelastic materials

The equations governing the motion of hypoelastic materials (and related models of non-Newtonian fluids) are not in conservation form. Hence there is no obvious formulation of Rankine-Hugoniot jump conditions across a shock. In this paper we demonstrate that a viscosity criterion can be used to obtain meaningful shock conditions. In particular, we discuss shocks of small amplitude. The shock conditions obtained will in general depend on the form of the viscosity term.Michael Renardy's research was supported by the National Science Foundation under Grant DMS-9008497.     

Weighted Compact Scheme for Shock Capturing

In this paper a new class of finite difference schemes - the Weighted Compact Schemes are proposed. According to the idea of the WENO schemes, the Weighted Compact Scheme is constructed by a combination of the approximations of derivatives on candidate stencils with properly assigned weights so that the non-oscillatory property is achieve when discontinuities appear. The primitive function reconstruction method of ENO schemes is applied to obtain the conservative form of the Weighted Compact Scheme. This new scheme not only preserves the characteristic of standard compact schemes and achieves high order accuracy and high resolution using a compact stencil, but also can accurately capture shock waves and discontinuities without oscillation. Numerical examples show the new scheme is very promising and successful.     

Nonlinear Strong Shock Interactions: A Shock-Fitted Approach

This paper addresses nonlinear effects which result from the interaction of shock waves with vortices. A series of experiments are carried out, which involve the interaction of a strong shock wave with a single plane vorticity wave and a randomly distributed wave system. These experiments are first conducted in the linear regime to obtain a mutual verification of theory and computation. They are subsequently extended into the nonlinear regime. A systematic study of the interaction of a plane shock wave and a single vortex is then conducted. Specifically, we investigate the conditions under which nonlinear effects become important, both as a function of shock Mach number, M ₁, and incident vortex strength (characterized by its circulation Γ). The shock Mach number is varied from 2 to 8, while the circulation of the vortex is varied from infinitesimally small values (linear theory) to unity. Budgets of vorticity, dilatation, and pressure are obtained. They indicate that nonlinear effects become more significant as both the shock Mach number and the circulation increase. For Mach numbers equal to 5 and above, the dilatation in the vortex core grows quadratically with circulation. An acoustic wave propagates radially outward from the vortex center. As circulation increases, its upstream-facing front steepens at low Mach numbers, and its downstream-facing front steepens at high Mach numbers. A high Mach number asymptotic expansion of the Rankine--Hugoniot conditions reveals that nonlinear effects dominate both the shock motion and the downstream flow for ΓM ₁ > 1. Received 28 June 1997 and accepted 25 November 1997     

Nonlinear Stability for Multidimensional Fourth-Order Shock Fronts

We consider the question of stability for planar wave solutions that arise in multidimensional conservation laws with only fourth-order regularization. Such equations arise, for example, in the study of thin films, for which planar waves correspond to fluid coating a pre-wetted surface. An interesting feature of these equations is that both compressive, and undercompressive, planar waves arise as solutions (compressive or undercompressive with respect to asymptotic behavior relative to the un-regularized hyperbolic system), and numerical investigation by Bertozzi, Münch, and Shearer indicates that undercompressive waves can be nonlinearly stable. Proceeding with pointwise estimates on the Green's function for the linear fourth-order convection–regularization equation that arises upon linearization of the conservation law about the planar wave solution, we establish that under general spectral conditions, such as appear to hold for shock fronts arising in our motivating thin films equations, compressive waves are stable for all dimensions d≧2 and undercompressive waves are stable for dimensions d≧3. (In the special case d=1, compressive waves are stable under a very general spectral condition.) We also consider an alternative spectral criterion (valid, for example, in the case of constant-coefficient regularization), for which we can establish nonlinear stability for compressive waves in dimensions d≧3 and undercompressive waves in dimensions d≧5. The case of stability for undercompressive waves in the thin films equations for the critical dimensions d=1 and d=2 remains an interesting open problem.