Shock Waves最新文献

Normal detonation reflection generates a shock wave that exhibits complicated dynamics as it propagates through the incident detonation and post-detonation flow. Ideal models have historically neglected the influence of a finite detonation thickness on the reflected shock due to its small size relative to laboratory scales. However, one-dimensional numerical simulations show that the reflected shock accelerates to a large shock speed not predicted by ideal theory as it propagates through the incident detonation. Analysis with a derived shock-change equation identifies the principal role of the highly nonuniform upstream flow on producing the large shock acceleration. Simulations of detonation reflection show how a finite detonation thickness affects the entire trajectory of the reflected shock.

Computational fluid dynamics shows that a shock wave can detach from the sharp leading edge of a curved wedge at a wedge angle smaller than the classical maximum flow deflection as well as the sonic wedge angle. This is attributed to the inability of the sonic flow, at the wedge trailing edge, to pass as much mass flow as is being admitted through the shock wave attached at the leading edge. At this condition, the flow is unsteady, causing both the sonic surface and the shock to make adjustments in their shapes and positions to achieve a steady state with mass-flow balance. As a result, the shock wave becomes detached. Time-accurate CFD calculations show the gasdynamic details of the adjustment where the flow and the detached shock assume a steady state as the mass-flow imbalance gradually decreases to zero. This mechanism of shock detachment, occurring near the leading edge, is called local choking to distinguish it from shock detachment due to global choking that occurs because of flow choking at the exit of a convergent duct and to distinguish it as well from detachment due to an excessive leading-edge deflection. The local choking mechanism has been postulated to be a cause of shock detachment from doubly curved wedges. An analysis, based on curved shock theory and confirmed by CFD, shows that local choking and shock detachment from a doubly curved leading edge are dependent on Mach number, wedge angle, wedge curvature (both streamwise and cross-stream), and wedge length.

The Hugoniot shock wave velocity ((U_{textrm{S}}))–particle velocity ((u_{textrm{p}})) curve of polymethyl methacrylate (PMMA) was measured in an experiment using only 2.5 g of a high explosive. The thickness of the plate was varied to accurately determine (U_{textrm{S}}) at an arbitrary position in the PMMA. Image analysis was conducted to obtain the x–t diagram of shock wave propagating in PMMA along the axis of the explosive, and its derivative was used to obtain the on-axis (U_{textrm{S}}) at an arbitrary location. Using the pressure measurement results and (U_{textrm{S}}) values, the Hugoniot (U_{textrm{S}})–(u_{textrm{p}}) curve of PMMA was obtained by calculating (u_{textrm{p}}) from the momentum conservation law. The results are in very good agreement with the reported values for flat-plate impact experiments conducted using an impact gun. It was found that the Hugoniot (U_{textrm{S}})–(u_{textrm{p}}) curve of PMMA on the low-pressure side ((u_{textrm{p}} < 0.5, {textrm{km}}/{textrm{s}})) can be evaluated with high accuracy using a simple measurement method that does not use plane waves.

Turbulence modeling has the potential to revolutionize high-speed vehicle design by serving as a co-equal partner to costly and challenging ground and flight testing. However, the fundamental assumptions that make turbulence modeling such an appealing alternative to its scale-resolved counterparts also degrade its accuracy for practical high-speed configurations, especially when fully 3D flows are considered. The current investigation develops a methodology to improve the performance of turbulence modeling for a complex Mach 8.3, 3D shock boundary layer interaction (SBLI) in a double fin geometry.A representative two-equation model, with low-Reynolds-number terms, is used as a test bed. Deficiencies in the baseline model are first elucidated using benchmark test cases involving a Mach 11.1 zero pressure gradient boundary layer and a Mach 6.17 flow over an axisymmetric compression corner. From among different possibilities, two coefficients are introduced to inhibit the non-physical over-amplification of (i) turbulence production and (ii) turbulence length-scale downstream of a shock wave. The coefficients rely on terms already present in the original model, which simplifies implementation and maintains computational costs. The values of the coefficients are predicated on the distribution of turbulence quantities upstream of the shock; this ensures that the modifications do not degrade the model predictions in simpler situations such as attached boundary layers, where they are unnecessary. The effects of the modifications are shown to result in significant improvements in surface pressure and wall heat flux for the 3D SBLI test case, which contains numerous features not observed in 2D situations, such as 3D separation, skewed boundary layers, and centerline vortices. Considerations on the inflow values of turbulence variables and mesh resolution are provided.

Full-field direct simulation of sonic boom has only been applied to the analysis of axisymmetric geometries. In this work, a more realistic analysis of complex geometries over buildings is achieved by employing a combination of the following four numerical approaches: (i) a hierarchical structured adaptive mesh refinement method, (ii) a ghost fluid method for incorporating the immersed boundary conditions on the solid–fluid interfaces, (iii) a well-balanced finite volume method to allow stable stratification of the atmosphere, and (iv) a segmentation method of the computational domain to increase the efficiency of the computations. The three-dimensional Euler equations with a gravitational source term are solved over a stratified atmosphere. The simulation is split into two stages. Firstly, the entire flow field that involves a delta wing body is solved without buildings. Thereafter, the flow behaviors near the ground are recomputed considering rectangular and L-type buildings. Computational results show that the near- and far-fields waveforms are comparable to those from the wind tunnel experiment and the waveform parameter method, respectively. The waveform shape behind the shock waves is spiked due to the diffracted waves around buildings, with the spiking effect in L-type buildings being stronger than that in rectangular buildings. The pressure rises for rectangular and L-type buildings are significantly amplified due to double and triple reflections, respectively, each with an amplification factor comparable to the theoretical value. These results indicate that full-field simulation is promising for analyzing three-dimensional characteristics of sonic boom emanating from complex geometries passing over buildings.

Time-resolved schlieren photography was used to visualise mixing zones induced by the Richtmyer–Meshkov instability. These were initiated with four different initial conditions: three of them with monotonic, single-mode shapes and one with a non-monotonic, multi-mode shape. These initial conditions were generated by an innovative experimental concept, the Micro Rotating Shutter System. The results of this experimental campaign reveal that the shape of the initial air–helium interface influences the subsequent development of the resulting mixing zone. Over the measurement time range, the width of the mixing zone induced by this instability is correctly fitted by a power law. Its growth exponent depends on the monotonicity of the initial air–helium interface: while mixing widths originating from single-mode initial conditions are almost superimposed, a lesser growth exponent is found for the multi-mode initial condition. The Reynolds number based on the width of the mixing zone suggests that both flows initiated with single- and multi-mode initial conditions reach a fully turbulent state after the interaction with the reflected shock wave (reshock). The schlieren photography visualisations presented here also allow to illustrate the structure of the induced mixing and highlight the effect of the initial conditions on the large-scale structures of the Richtmyer–Meshkov instability-induced mixing.

An ongoing rotating detonation rocket engine program is investigating the influence of combustor annulus radii on RDRE operating characteristics with flat-faced impinging injectors. To facilitate the isolation of all but the radius of curvature effects in the experiments, the annular gap was kept constant at 5 mm in combustors having either 25-mm or 51-mm outer diameter. The mixing processes were kept similar by utilizing injectors with the same net injector-to-annular gap area ratio (AR = 0.11), same radial separation distance of the orifices, and same center-of-gap impingement distance from the front-end wall. The wave dynamics, plenum pressure, and axial pressure profiles in these RDREs were compared over the mass flux and equivalence ratio ranges of (80{-}400,text {kg/s/m}^{2}) and 0.26(-)2.6, respectively, with gaseous methane–oxygen propellant. Experiments showed that stable one-wave operation would occur in the 25-mm RDRE at most mass fluxes where stable two-wave operation was established in the 51-mm RDRE. Stable one-wave operation with a single counter-rotating wave was maintained in the 51-mm RDRE at mass fluxes of (240,text {kg/s/m}^{2}) and below. Under these fueling conditions in the 25-mm RDRE, a counter-rotating wave also appeared while it operated with a single dominant wave. The wave spin speeds were typically 20–40% less than the Chapman–Jouguet detonation speed of the propellant and depended only on mass flux and wave number rather than the annulus diameter.

We studied the real gas effect on the ignition characteristics in chemical reactors with one-step irreversible reaction. The real gas effects were characterized by the inter-molecular attraction term ((alpha )) and the finite molecular volume term ((beta )). The Noble-Abel and van der Waals equations of state were employed to derive non-dimensional reactor models. In addition to ideal reactors, i.e., constant volume and constant pressure, non-ideal reactors that account for the non-ideal pressure variation in shock tube and rapid compression machine were also considered. For all reactors, low value of (alpha /beta ) and high value of (beta ) (approximately (alpha /beta <{{1.0}}) and (beta >{{0.1}})) induce a decrease of the ignition delay-time, while high value of both (alpha /beta ) and (beta ) (approximately (alpha /beta >{{2.0}}) and (beta >{{0.1}})) induces an increase of the ignition delay-time. The variations of the ignition delay-time induced by real gas effects are mainly related to the change of the fugacity coefficient with (alpha ) and (beta ). Additional contributions are due to the real gas heat capacity at constant pressure when considering a constant pressure reactor and to non-ideal volume variation when considering non-ideal reactors. The impact of various parameters was also investigated, including the heat capacity ratio of perfect gas, the reduced activation energy of the one-step reaction, and the heat content of the mixtures. Comparison with simulation performed with detailed reaction mechanisms and considering real gas models demonstrates that the present approach constitutes a rapid and simple, yet qualitatively or even quantitatively accurate method to assess the need of accounting for real gas effects to model chemical kinetics under high-pressure conditions.

To study the shock wave initiation characteristics of 2,4,6-trinitrotoluene (TNT) under different charging types, the shock wave pressure and shock wave attenuation of standard Pentolite explosives under different diaphragm thicknesses were quantitatively studied using the ion probe method. The gap tests of three explosives were carried out, including pressed TNT without restraint, pressed TNT with steel pipe restraint, and cast TNT with steel pipe restraint. The shock wave initiation pressures of TNT under the three different conditions were compared. Moreover, combined with the numerical simulation technology, the critical initiation pressure and the pressure cloud diagram of the gap test of TNT were obtained, and the dynamic change process of the shock wave in the diaphragm was acquired, which was difficult to measure in the experiments. The results showed that the critical initiation pressure of pressed TNT was significantly lower than that of cast TNT and that restraint can reduce the measured critical initiation pressure of TNT under certain conditions. Therefore, the research results may provide a basis for the damage range of TNTs with different charging types and the determination of the safety protection distance of shock wave initiation.

The blast loading from a detonation of a high explosive charge at high altitude is quite different from that at sea level. Due to diminished ambient pressure, the damage caused by the blast load may be more minor at high altitude. However, the shock wave parameters at diminished ambient pressure have not yet been thoroughly studied. In this research, experiments were carried out to study the relation between ambient air pressure and shock wave parameters. The explosion experiments were carried out in a sealed explosion chamber with an initial pressure of 95 kPa, 74 kPa, and 57 kPa. For these three atmospheric conditions, the history profiles of incident shock wave pressure generated by TNT charges of 106 g and 292 g were recorded. The influence of ambient pressure and temperature on the shock wave parameters was analyzed through numerical simulations. By analyzing the experimental and numerical data, it was found that ambient pressure is the main factor affecting the shock wave parameters, while the effect of temperature is not so obvious. Furthermore, based on the analysis of experimental data, formulas for evaluating shock wave overpressure, specific impulse, and arrival time using the Sachs variables are given, and the shock wave parameters at an altitude of 5000 m are calculated using these formulas. The observed maximum reduction in the shock wave overpressure was 23%, in specific impulse 27%, and in arrival time 12%, compared to the results calculated at sea level. The results can be applied to blast-resistant analyses of buildings in low-pressure environment.