The shock expansion tube and its application as a sonic boom simulator
Shock expansion tube characteristics and use as sonic boom simulator
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Shock expansion tube characteristics and use as sonic boom simulator
Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.
The Richtmyer–Meshkov instability (RMI) of a twice-shocked gas interface is researched using high-speed planar-laser induced fluorescence (PLIF) in the Wisconsin Shock Tube Laboratory’s vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium-acetone mixture and argon. This IC is accelerated by a shock of nominal strength M = 1.8, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An estimate of the light gas mole fraction is extracted from high-speed imaging using an iterative process that accounts for the nonlinear temperature dependence of the acetone’s fluorescence quantum yield (FQY) and absorption cross-section. A vorticity deposition model for the initial growth rate after reshock is compared with the Mikaelian model for re-shock. Previously used in literature, the number of generations is shown to naturally arise from a normalisation of the scalar transport equation. A self-similar analysis is then further performed using the mole fraction data to explore the evolution of the RMI after reshock and the higher order moments of the light gas mole fraction are compared with a proposed model.
The Richtmyer–Meshkov instability of a twice-shocked gas interface is studied using high-speed planar laser-induced fluorescence in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. Furthermore, this initial condition is accelerated by a shock of nominal strength M = 1.9, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. Three individual experiments are analyzed, the energy spectrum and the structure functions of the light gas mole fraction field are calculated and compared.
The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.
Radiative enhancement factors for the CN violet and N2(+) first negative band systems caused by nonequilibrium thermochemistry in the shock layer of a blunt-nosed vehicle during earth entry are reported. The results are based on radiometric measurements obtained with the aid of a combustion-driven shock tube. The technique of converting the shock-tube measurements into predictions of the enhancement factors for the blunt-body case is described, showing it to be useful for similar applications of other shock-tube measurements.
Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.
Subsonic Ludwieg /simple expansion/ tube as shock tube for aerodynamic testing, examining flow characteristics
Present test facilities are seriously limited regarding investigations involving high Reynolds numbers due to financial considerations. Quasi-steady testing facilities offer a practical immediate solution to the problem of high-Re testing. A familiar example is the blowdown wind tunnel, but even more flexibility and economy may be provided by using shock-tube devices. The Ludwieg tube is the shock-tube device most often proposed as a means of generating high-Re flows. Two-stage nozzles may be used with a Ludwieg tube. Quasi-steady facilities will be useful only if the available test time exceeds that required to establish steady flow.
A brief review of laser elements is given. Flowing gas lasers are represented to have the best potential for high average power. The background of shock-tube researchers and the shock tube itself are alleged to be ideally suited for the development of such lasers. Three types - the electric discharge, the gasdynamic, and the chemical laser - are discussed briefly. A legion number of possible gas lasers is enumerated. With the development of their potential for higher power and efficiency, many additional and important uses of lasers are predicted, even beaming power through space for long distances, up to 1 AU. A few details of some current high-power gasdynamic laser devices are given.
The experimental program for validating real gas hypersonic flow codes at NASA Ames is described. Ground based test facilities used include ballistic ranges, shock tubes and shock tunnels, arcjet facilities and heated air hypersonic wind tunnels. Also included are large scale computer systems for kinetic theory simulations and benchmark code solutions. Flight tests consist of the Aeroassist Flight Experiment, the Space Shuttle, Project Fire 2, and planetary probes such as Galileo, Pioneer Venus and PAET.
The experimental program for validating real gas hypersonic flow codes at NASA Ames Rsearch Center is described. Ground-based test facilities used include ballistic ranges, shock tubes and shock tunnels, arc jet facilities and heated-air hypersonic wind tunnels. Also included are large-scale computer systems for kinetic theory simulations and benchmark code solutions. Flight tests consist of the Aeroassist Flight Experiment, the Space Shuttle, Project Fire 2, and planetary probes such as Galileo, Pioneer Venus, and PAET.
We synthesize profiles of the infrared line Fe I 15648.5 A (g = 3) for a recently developed theoretical model of siphon flows along photospheric magnetic loops. The synthesized line profiles are compared with the observations from which Rueedi et al. (1992) deduced the presence of such flows across the neutral line of an active region plage. This comparison supports the interpretation of Rueedi et al. (1992). It also suggests that the average footpoint separation of the observed loops carrying the siphon flow is 8-15 sec and that the siphon flow experiences a standing tube shock in the downstream leg near the top of the arch.
The Entry Systems Modeling project (ESM) is supported by both the NASA Space Technology and the Science Mission Directorates and focuses on developing simulation tools and validated models for characterizing the performance of entry systems tailored to planetary destinations across the Solar System. ESM is organized into six technical capability areas that together address all relevant factors related to spacecraft entry, as well as some aspects of descent: Thermal Protection System (TPS) Materials; Aerothermodynamics; Entry & Descent Vehicle Dynamics; Guidance, Navigation, and Control; Vehicle Systems Analysis; and Advanced Tools and Numerical Methods. Development within the capability areas is undertaken explicitly with a focus on transition and infusion to science missions, human exploration missions, and commercial space activities. The present talk details developments that specifically impact science missions, including simulation tool capabilities that aid in mission design and model development to understand entry system performance at a given destination. Examples of the successful infusion and transition of such project outcomes to science missions also are provided. Several simulation tool development efforts within ESM have resulted in new design capabilities for missions. One such outcome is improved toolsets for mission trajectory and concept of operations design. Specifically, an initiative to couple a leading tool for entry, ascent/descent, and orbital trajectory optimization (Program to Optimize Simulated Trajectories II or POST2) to those used within the Agency for interplanetary trajectory optimization (Copernicus and Monte) has made substantial progress, with the outcomes to date promising to allow efficient trajectory optimization across mission phases. Additionally, toolchains for the evaluation of vehicle performance during entry and descent have been developed that allow assessment of multi-dimensional aeroheating on detailed vehicle geometries, characterization of deployment and inflation of parachutes, and assessment of vehicle dynamic stability during descent. These capabilities are achieved by coupling diverse sets of physics together – material response, computational fluid dynamics, radiation, and vehicle dynamics – to suitably describe complex entry and descent phenomena. Several model development and validation efforts for specific destinations and entry regimes also are underway within the ESM project. For instance, new experimental capabilities to validate radiation models at low densities/high altitudes recently have been established with project support, specifically the Low-Density Shock Tube (LDST) at the NASA Ames Research Center Electric Arc Shock Tube (EAST) facility. The LDST is being leveraged to develop improved models of shock layer kinetics and radiation in Titan atmospheres, while future studies will be conducted in the LDST and the existing high velocity shock tube to provide validation data for radiation models of Venus, Ice Giants, and Mars atmospheres. Models describing the aerothermal and thermo-structural performance of Thermal Protection System (TPS) materials has been another focus, with multiscale modeling activities on-going for the two leading TPS materials applicable to a range of entry conditions and science missions: the Phenolic-Impregnated Carbon Ablator (PICA) and woven materials like 3D Mid-Density Carbon Phenolic (3MDCP). A continual effort is made to infuse and transition outcomes from ESM simulation tool and model development activities into relevant science missions. Significant progress has been made on this front, with missions such as Dragonfly, DAVINCI, and Mars Missions benefitting from project outcomes. The groundwork also is being laid to provide insights into forward looking missions to Gas/Ice Giants as well as for potential sample returns.
The ignition, oxidation, and pyrolysis chemistry of methyl propyl ether (MPE) was probed experimentally at several different conditions, and a comprehensive chemical kinetic model was constructed to help understand the observations, with many of the key parameters computed using quantum chemistry and transition state theory. Experiments were carried out in a shock tube measuring time variation of CO concentrations, in a flow tube measuring product concentrations, and in a rapid compression machine (RCM) measuring ignition delay times. The detailed reaction mechanism was constructed using the Reaction Mechanism Generator software. Sensitivity and flux analyses were used to identify key rate and thermochemical parameters, which were then computed using quantum chemistry to improve the mechanism. Validation of the final model against the 1-20 bar 600-1500 K experimental data is presented with a discussion of the kinetics. The model is in excellent agreement with most of the shock tube and RCM data. Strong non-monotonic variation in conversion and product distribution is observed in the flow-tube experiments as the temperature is increased, and unusually strong pressure dependence and significant heat release during the compression stroke is observed in the RCM experiments. These observations are largely explained by a close competition between radical decomposition and addition to O2 at different sites in MPE; this causes small shifts in conditions to lead to big shifts in the dominant reaction pathways. The validated mechanism was used to study the chemistry occurring during ignition in a diesel engine, simulated using Ignition Quality Test (IQT) conditions. At the IQT conditions, where the MPE concentration is higher, bimolecular reactions of peroxy radicals are much more important than in the RCM.
Detailed spectrally and spatially resolved radiance has been measured in the Electric Arc Shock Tube for conditions relevant to Titan entry, varying atmospheric composition, free-stream density (equivalent to altitude) and shock velocity. Permutations in atmospheric composition include 1.1, 2, 5 and 8.6 CH4 by mole with a balance of N2 and 1.5 CH4 0.5 Ar 98 N2 by mole, which is consistent with the current understanding of Titan's atmosphere. The effect of gas impurities identified in previous shock tube studies were also examined by testing in pure N2 and deliberate addition of air to the CH4N2 mixtures. The test campaign measured radiation at velocities from 4.7 kms to 8 kms and free-stream pressures from 0.1 to 0.47 Torr. These conditions cover a range of potential trajectories for flight missions, including a direct ballistic trajectory, a fly by or an extremely high speed entry. Radiances measured in this work are substantially larger compared to that reported both in past EAST test campaigns and other shock tube facilities. Depending on the metric used for comparison, the discrepancy can be as high as an order of magnitude. Potential causes for the discrepancy, such as the effect of oxygen due to Air leakage, gas composition and purity are discussed. The present work provides a new benchmark set of data to replace those published in previous studies.
Results of research on uranium and tungsten spectra are summarized. Measurements of visible line spectra and opacities were carried out on shock tube plasmas which, prior to shock compression, were mixtures of rare gases and UF6 or WF6. Opacities were compared to theoretical predictions. Feasibility of light source methods other than the shock tube was explored for future applications in the spectroscopy of heavy metals and ions.
Flame modeling studies have highlighted the role of branching in the kinetics for the prompt NO switch reaction, NCN + H. In this reaction, there is strong competition between the CH + N2 (R1) and HCN + N (R2) product channels. Increased branching towards reaction R2 promotes the subsequent formation of Fenimore or Prompt NO. Recent direct shock tube studies on this complex reaction conclude that branching to reaction R2 is much more favorable than the predicted branching from theoretical studies. The experimentally predicted prompt NO switch temperature (TS), i.e., the temperature at which k2/(k1+k2) = 0.5, is at 1670 K, in contrast to theory that predicts Ts > 3200 K. The present theory/modeling work was initiated to highlight potential causes for this significant discrepancy between experiment and theory. Our analysis of the simulations reported in the experimental studies indicate that including a fall-off representation for C2H5 radical dissociation (with C2H5 produced from C2H5I dissociation and acting as an H-atom source) has a noticeable influence on the simulations for NCN decay and HCN formation in the shock tube studies. With the inclusion of a theory-based pressure fall-off representation for C2H5 dissociation, we show that this radical persists for much longer timescales even at the high-temperatures in the shock tube studies. This persistence then facilitates the rapid conversion of reactive H-atoms to CH3 radicals at the shock tube conditions via the H + C2H5 → CH3 + CH3 reaction. An updated theoretical analysis for this classic addition-elimination reaction is also provided in this work. The rapid formation of CH3 radicals necessitates the inclusion of additional reactions to adequately simulate the NCN and HCN data. Our simulations conclude that, with the inclusion of these additional reactions, the NCN and HCN profiles can be reasonably well simulated by the most recent theoretical predictions.