Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock initiation”

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.

At least 1,297 records · Page 72

The initial stage of development of type 4 radio bursts and the relation to expanding magnetic bottles

Using observed data for wide-band type IV solar radio bursts, the onset time differences between the microwave and metric frequencies and the peak flux intensities of the metric component were analyzed as a function of the longitudinal position of associated flares on the solar disk. It is shown that this time difference is dependent on the position of associated flares and that the peak flux intensity reaches maximum when a flare occurs in the region 10 to 40 deg west of the central meridian of the solar disk. These results are explained by taking into account the eastward expansion of magnetic bottles which trap mildly relativistic electrons responsible for type IV bursts. Discussion is also made of the relation between these magnetic bottles and shock waves which excite type II radio bursts.

Sakurai, K.↗

High-energy emission from a magnetar giant flare in the Sculptor galaxy

Magnetars are the most highly magnetized neutron stars in the cosmos (with magnetic field 10 13 –10 15 G). Giant flares from magnetars are rare, short-duration (about 0.1 s) bursts of hard X-rays and soft γ rays. Owing to the limited sensitivity and energy coverage of previous telescopes, no magnetar giant flare has been detected at gigaelectronvolt (GeV) energies. Here, we report the discovery of GeV emission from a magnetar giant flare on 15 April 2020 (refs. 3,4 and A. J. Castro-Tirado et al., manuscript in preparation). The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope detected GeV γ rays from 19 s until 284 s after the initial detection of a signal in the megaelectronvolt (MeV) band. Our analysis shows that these γ rays are spatially associated with the nearby (3.5 megaparsecs) Sculptor galaxy and are unlikely to originate from a cosmological γ-ray burst. Thus, we infer that the γ rays originated with the magnetar giant flare in Sculptor. Furthermore, we suggest that the GeV signal is generated by an ultra-relativistic outflow that first radiates the prompt MeV-band photons, and then deposits its energy far from the stellar magnetosphere. After a propagation delay, the outflow interacts with environmental gas and produces shock waves that accelerate electrons to very high energies; these electrons then emit GeV γ rays as optically thin synchrotron radiation. This observation implies that a relativistic outflow is associated with the magnetar giant flare, and suggests the possibility that magnetars can power some short γ-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

Bridging the Bondi and Event Horizon Scales: 3D GRMHD Simulations Reveal X-shaped Radio Galaxy Morphology

X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up ≲10% of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole (SMBH), SMBH merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the SMBH sphere of influence of radius RB to the SMBH of gravitational radius Rg at the largest scale separation, RB/Rg = 103, to date. Initially, our axisymmetric system of constant-density hot gas contains a weak vertical magnetic field and rotates in the equatorial plane of a rapidly spinning SMBH. We seed the gas with small-scale 2% level pressure perturbations. Infalling gas forms an accretion disk, and the SMBH launches relativistically magnetized collimated jets reaching well outside RB. Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the SMBH, becomes dynamically important, and leads to a magnetically arrested disk state. The SMBH accretes at 2% of the Bondi rate ($\dot{M}\simeq 2.4\times {10}^{-3}{M}_{\odot }\,{\mathrm{yr}}^{-1}$ for M87*) and launches twin jets at η = 150% efficiency. These jets are powerful enough (P jets ≃ 2 × 10 44 erg s -1 ) to escape along the SMBH spin axis and end the short-lived intermittent jet state, whose transient nature can account for the rarity of XRGs.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical solution of the unsteady Navier-Stokes equation

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws are discussed. These schemes share many desirable properties with total variation diminishing schemes, but TVD schemes have at most first-order accuracy, in the sense of truncation error, at extrema of the solution. In this paper a uniformly second-order approximation is constructed, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Osher, Stanley J.↗

Ultra-high Temperature Ceramic Composites for Leading Edges

Ultra-high temperature ceramics (UHTC) have performed unreliably due to material flaws and attachment design. These deficiencies are brought to the fore by the low fracture toughness and thermal shock resistance of UHTC. If these deficiencies are overcome, we are still faced with poor oxidation resistance as a limitation on UHT applicability to reusable launch vehicles. We have been addressing the deficiencies of UHTC for the past two years via a small task at GRC that is in the Airframe part of the Next Generation Launch Technology Program. Our focus is on composite constructions and functional grading to address the mechanical issues and on composition modification to address the oxidation issue. The progress on approaches to improving oxidation resistance by alloying and functional grading will be reported. In particular, initial tests of tantalum additions have shown potential for major improvement. Less promising results of additional tests at higher temperatures will be presented. Two fabrication approaches are being pursued to produce carbon fiber reinforced UHTC composites: prepregging and rigid perform infiltration. Fabrication procedures and microstructures for composites will be reported.

Levine, Stanley R.↗

Ultra-High Temperature Ceramic Composites for Leading Edges

Ultra-high temperature ceramics (UHTC) have performed unreliably due to material flaws and attachment design. These deficiencies are brought to the fore by the low fracture toughness and thermal shock resistance of UHTC. If these deficiencies are overcome, we are still faced with poor oxidation resistance as a limitation on UHTC applicability to reusable launch vehicles. We have been addressing the deficiencies of UHTC for the past two years via a small task at GRC that is in the Airframe part of the Next Generation Launch Technology Program. Our focus is on composite constructions and functional grading to address the mechanical issues and on composition modification to address the oxidation issue. The progress on approaches to improving oxidation resistance by alloying and functional grading will be reported. In particular, initial tests of tantalum additions have shown potential for major improvement. Results for additional tests at higher temperatures will be presented. These oxidation improvements are being incorporated in the composites approaches. Two fabrication approaches are being persued to produce carbon fiber reinforced UHTC composites: prepregging and rigid perform infiltration. Fabrication procedures, microstructures, and initial mechanical property and oxidation results for composites will be reported.

Levine, Stanley R.↗

Flood basalts and mass extinctions

There appears to be a correlation between the times of flood basalts and mass-extinction events. There is a correlation of flood basalts and hotspot tracks--flood basalts appear to mark the beginning of a new hotspot. Perhaps there is an initial instability in the mantle that bursts forth as a flood basalt but then becomes a steady trickle that persists for many tens of millions of years. Suppose that flood basalts and not impacts cause the environmental changes that lead to mass-extinctions. This is a very testable hypothesis: it predicts that the ages of the flows should agree exactly with the times of extinctions. The Deccan and K-T ages agree with this hypothesis; An iridium anomaly at extinction boundaries apparently can be explained by a scaled-up eruption of the Hawaiian type; the occurrence of shocked-quartz is more of a problem. However if the flood basalts are all well dated and their ages indeed agree with extinction times, then surely some mechanism to appropriately produce shocked-quartz will be found.

Morgan, W. Jason↗

Time-Dependent Photodissociation Regions

We present theoretical models of the time-dependent thermal and chemical structure of molecular gas suddenly exposed to far-ultraviolet (FUV) (6 eV less than hv less than 13.6 eV) radiation fields and the consequent time- dependent infrared emission of the gas. We focus on the response of molecular hydrogen for cloud densities ranging from n = 10(exp 3) to 10(exp 6)/cu cm and FUV fluxes G(sub 0) = 10(exp 3)-10(exp 6) times the local FUV interstellar flux. For G(sub 0)/n greater than 10(exp -2) cu cm, the emergent H(sub 2) vibrational line intensities are initially larger than the final equilibrium values. The H(sub 2) lines are excited by FUV fluorescence and by collisional excitation in warm gas. Most of the H(sub 2) intensity is generated at a characteristic hydrogen column density of N approximately 10(exp 21)/sq cm, which corresponds to an FUV optical depth of unity caused by dust opacity. The time dependence of the H(sub 2) intensities arises because the initial abundances of H(sub 2) at these depths is much higher than the equilibrium values, so that H(sub 2) initially competes more effectively with dust in absorbing FUV photons. Considerable column densities of warm (T approximately 1000) K H(sub 2) gas can be produced by the FUV pumping of H(sub 2) vibrational levels followed by collisional de-excitation, which transfers the energy to heat. In dense (n greater than or approximately 10(exp 5)/cu cm) gas exposed to high (G(sub 0) greater than or approximately 10(exp 4)) fluxes, this warm gas produces a 2-1 S(1)/1-0 S(l) H(sub 2) line ratio of approximately 0.1, which mimics the ratio found in shocked gas. In lower density regions, the FUV pumping produces a pure-fluorescent ratio of approximately 0.5. We also present calculations of the time dependence of the atomic hydrogen column densities and of the intensities of 0 I 6300 A, S II 6730 A, Fe II 1.64 microns, and rotational OH and H20 emission. Potential applications include star-forming regions, clouds near active galactic nuclei, and planetary nebulae. We apply our models to five planetary nebulae and conclude that only BD +30deg3639 shows evidence of enhanced H(sub 2) emission due to (high) nonequilibrium H(sub 2) abundances.

Hollenbach, David↗

The Late Peaking Afterglow of GR8 100418A

GRB 100418A is a long Gamma-Ray Burst at redshift z=0.6235 discovered with the Swift Gamma-Ray Burst Explorer with unusual optical and X-ray light curves ' After an initial short-lived, rapid decline in X-rays, the optical and X-ray light curves observed with Swift are approximately flat or rising slightly out to at least approx.7 ks after the trigger, peak at approx.50 ks, and then follow an approximately power-law decay. Such a long optical plateau and late peaking is rarely seen in 6R8 afterglows. Observations with REM during a gap in the Swift coverage indicate a bright optical flare at approx.25 ks, The long plateau phase of the afterglow is interpreted using either a model with continuous injection of energy into the forward shock of the burst or a model in which the 'et of the burst is viewed off-axis. In both models the isotropic kinetic energy in the late afterglow after the plateau phase is >100 times the 10(exp 51) erg of the prompt isotropic gamma-ray energy release. The energy injection model is favored because the off-axis 'et model would require the intrinsic $T f801$ for the GR8 'et viewed on-axis to be very short, approx.10 ms, and the intrinsic isotropic gamma-ray energy release and the true jet energy to be much higher than the typical values of known short GRBs^ The non-detection of a 'et break up to approx.2 Ms indicates a jet half-opening angle of at least 14 degrees, and a relatively high collimation-corrected 'et energy of at least 10(exp 52) erg.

Marshall, Frank↗

Plasma Clouds and Snowplows: Bulk Plasma Escape from Mars Observed by MAVEN

We present initial Mars Atmosphere and Volatile EvolutioN (MAVEN) observations and preliminary interpretation of bulk plasma loss from Mars. MAVEN particle and field measurements show that planetary heavy ions derived from the Martian atmosphere can escape in the form of discrete coherent structures or "clouds." The ions in these clouds are unmagnetized or weakly magnetized, have velocities well above the escape speed, and lie directly downstream from magnetic field amplifications, suggesting a "snowplow" effect. This postulated escape process, similar to that successfully used to explain the dynamics of active gas releases in the solar wind and terrestrial magnetosheath, relies on momentum transfer from the shocked solar wind protons to the planetary heavy ions, with the electrons and magnetic field acting as intermediaries. Fluxes of planetary ions on the order of 10(exp 7)/sq cm/s can escape by this process, and if it operates regularly, it could contribute 10-20% of the current ion escape from Mars.

Halekas, J. S.↗

Data Mining for Vortices on the Earth's Magnetosphere - Algorithm Application for Detection and Analysis

Unsteady processes in the solar wind– magnetosphere interaction, such as vortices developed at the magnetopause boundary by the Kelvin–Helmholtz instability, may contribute to the process of mass, momentum and energy transfer into the Earth’s magnetosphere. The research described in this paper validates an algorithm to automatically detect and characterize vortices based on velocity data from simulations. The vortex identification algorithm (VIA) systematically searches the 3-D velocity fields to identify critical points where the magnitude of the velocity vector vanishes. The velocity gradient tensor is computed and its invariants are used to assess vortex structure in the flow field. We use the Community Coordinated Modeling Center (CCMC) Runs on Request capability to create a series of model runs initialized from the conditions observed by the Cluster mission in the Hwang et al. (2011) analysis of Kelvin–Helmholtz vortices observed during southward interplanetary magnetic field (IMF) conditions. We analyze further the properties of the vortices found in the runs, including the velocity changes within their motion across the magnetosheath. We also demonstrate the potential of our tool to identify and characterize other transient features (e.g., flux transfer events, FTEs) with vortical internal structures. We find that the vortices are associated with flows on the magnetosheath side of the magnetopause that reach speeds greater than the solar wind speed at the bow shock.

Collado-Vega, Yaireska M.↗

Initial Results on the Extraterrestrial Component of New Sediment Cores Containing Deposits of the Eltanin Impact Event

Background: The impact of the Eltanin asteroid into the Bellingshausen Sea (2.15 Ma) is the only known impact in a deep-ocean (approx. 5 km) basin. In 1995, Polarstern expedition ANT XII/4 made the first geological survey of the suspected impact region. Three sediment cores sampled around the San Martin seamounts (approx. 57.5 S, 91 W) contained well-preserved impact deposits. Sediments of Eocene age and younger were ripped up and redeposited by the impact. The depositional sequence produced by the impact has three units: a chaotic assemblage of sediment fragments up to 50 cm, followed by laminated sands deposited as a turbulent flow, and finally silts and clays that accumulated from dispersed sediments in the water column. The meteoritic impact ejecta, which is composed of shock-melted asteroidal materials and unmelted meteorites, settled through the water column and concentrated near the top of the laminated sands.

Kyte, Frank T.↗

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

Chondrites and the Protoplanetary Disk, Part 1

The papers discussed the following: The Formation Process of Adhering and Consorting Compound Chondrules Inferred Their Petrology and Major-Element Composition. The Prospect of High-Precision Pb Isotopic Dating of Meteorites. Evolution of UV-Irradiated Protoplanetary Disks. A Model for the Formation of E Chondrites. Oxygen Isotopic Diffusion and Exchange Experiments on Olivine and Chondrule Melts: Preliminary Results. Shock Heating: Origin of Shock Waves in the Protoplanetary Disk. Thermal Structures of Protoplanetary Disks. Meteoritical Astrophysics: A New Subdiscipline. Origin and Thermal History of FeNi-Metal in Primitive Chondrites. The Collisions of Chondrules Behind Shock Waves. Primary Signatures of the Nebular Dust Preserved in Accretionary Rims and Matrices of CV Chondrites. History of Thermally Processed Solids in the Protoplanetary Disk: Reconciling Theoretical Models and Meteoritical. Evidence Evaporation and Condensation During CAI and Chondrule Formation. Shock Heating: Effects on Chondritic Material. Rhounite-bearing Inclusions E201 and E202 from Efremovka: Constraints from Trace. Element Measurements Element Mapping in Anhydrous IDPs: Identification of the Host Phases of Major/Minor Elements as a Test of Nebula Condensation Models. Theoretical Studies of Disk Evolution Around Solar Mass Stars. Chemical Effects of High-Temperature Processing of Silicates. I-Xe and the Chronology of the Early Solar System. The Effects of X-Rays on the Gas and Dust in Young Stellar Objects. Origin of Short-lived Radionuclides in the Early Solar System. On Early Solar System Chronology: Implications of an Initially Heterogeneous Distribution of Short-lived Radionuclides. The Origin of Short-lived Radionuclides and Early Solar System Irradiation. Disequilibrium Melting and Oxygen Isotope Exchange of CAIs and Chondrules in the Solar Nebula. Mineralogy and Chemistry of Fine-grained Matrices, Rims, and Dark Inclusions in the CR Carbonaceous Chondrites Acfer/El Djouf 001 and the Ungrouped Carbonaceous Chondrites Acfer 094 and Adelaide. Oxygen Isotopes of Aluminum-rich Chondrules from Unequilibrated Enstatite Chondrites.

Source record↗

Merging neutron stars. 1. Initial results for coalescence of noncorotating systems

We present three-dimensional Newtonian simulations of the coalescence of two neutron stars, using a smoothed particle hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes (gamma = 2.4) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates result in greater mass loss and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of 2.4 solar mass, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of approximately 10 MeV, while the disk is heated by shocks to a temperature of 2 to 4 MeV.

Davies, M. B.↗

Scramjet Performance Assessment Using Water Absorption Diagnostics (U)

Simultaneous multiple path measurements of temperature and H2O concentration will be presented for the AIMHYE test entries in the NASA Ames 16-Inch Shock Tunnel. Monitoring the progress of high temperature chemical reactions that define scramjet combustor efficiencies is a task uniquely suited to nonintrusive optical diagnostics. One application strategy to overcome the many challenges and limitations of nonintrusive measurements is to use laser absorption spectroscopy coupled with optical fibers. Absorption spectroscopic techniques with rapidly tunable lasers are capable of making simultaneous measurements of mole fraction, temperature, pressure, and velocity. The scramjet water absorption diagnostic was used to measure combustor efficiency and was compared to thrust measurements using a nozzle force balance and integrated nozzle pressures to develop a direct technique for evaluating integrated scramjet performance. Tests were initially performed with a diode laser tuning over a water absorption feature at 1391.7 nm. A second diode laser later became available at a wavelength near 1343.3 nm covering an additional water absorption feature and was incorporated in the system for a two-wavelength technique. Both temperature and mole fraction can be inferred from the lineshape analysis using this approach. Additional high temperature spectroscopy research was conducted to reduce uncertainties in the scramjet application. The lasers are optical fiber coupled to ports at the combustor exit and in the nozzle region. The output from the two diode lasers were combined in a single fiber, and the resultant two-wavelength beam was subsequently split into four legs. Each leg was directed through 60 meters of optical fiber to four combustor exit locations for measurement of beam intensity after absorption by the water within the flow. Absorption results will be compared to 1D combustor analysis using RJPA and nozzle CFD computations as well as to data from a nozzle metric balance measuring thrust and integrated pressure measurements along the length of the nozzle. Assessment of its value as a combustor performance evaluation tool will be conducted.

Cavolowsky, John A.↗

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Here, using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Dynamics of Computing a Chemically Relaxed Nonequilibrium Flow

In order to gain insights into the strong dependence of numerical solutions on initial data for finite time steps, a set of nonlinear test problems rich enough to capture the behavior of difference schemes were recently identified and the numerical basins of attraction for these problems were computed using commonly used time discretizations in CFD. Our study revealed a wealth of surprisingly nonlinear behavior of numerical schemes that were not observed before, in particular for the implicit time discretizations that are commonly used in CFD. The goal of this work is to apply these tools to study a practical model from non-equilibrium flowfield relaxation. This type of problem arises in chemically nonequilibrium hypersonic flows such as in a shock tube experiment or an expanding nozzle. Here we consider a reacting mixture of (N2, N) for an inviscid one-dimensional steady model. Preliminary numerical results indicate that, aside from the possibility of spurious numerical solutions being introduced by the time discretizations, limitations on the model for physical or accurate solutions may also play a part in the dynamics observed.

Sweby, P. K.↗