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 397 records · Page 22

Observations of ultraviolet variability in RV Tauri stars

An IUE program to monitor the ultraviolet variability in RV Tauri stars was initiated. The Mg II region was investigated as a potential probe of atmospheric shocks, which are believed to be associated with the pulsational variability of this class of objects. Observations, a description of the spectra, and findings for V Vul and AC Her are presented. The Mg II emission does vary significantly during the cycle; major changes in the emission line strength occur on a time scale much less than 0.2 in phase; and as the UV (and optical) continuum flux increases, the Mg II lines decrease and increased emission may be seen at 2823, 2844, and 2900 A.

Brugel, Edward W.↗

Star-forming galaxy models: Blending star formation into TREESPH

We have incorporated star-formation algorithms into a hybrid N-body/smoothed particle hydrodynamics code (TREESPH) in order to describe the star forming properties of disk galaxies over timescales of a few billion years. The models employ a Schmidt law of index n approximately 1.5 to calculate star-formation rates, and explicitly include the energy and metallicity feedback into the Interstellar Medium (ISM). Modeling the newly formed stellar population is achieved through the use of hybrid SPH/young star particles which gradually convert from gaseous to collisionless particles, avoiding the computational difficulties involved in creating new particles. The models are shown to reproduce well the star-forming properties of disk galaxies, such as the morphology, rate of star formation, and evolution of the global star-formation rate and disk gas content. As an example of the technique, we model an encounter between a disk galaxy and a small companion which gives rise to a ring galaxy reminiscent of the Cartwheel (AM 0035-35). The primary galaxy in this encounter experiences two phases of star forming activity: an initial period during the expansion of the ring, and a delayed phase as shocked material in the ring falls back into the central regions.

Mihos, J. Christopher↗

The location of planetary bow shocks: A critical overview of theory and observations

A bow shock (BS has been observed in the collisionless solar wind upstream of every planet except Pluto, which has yet to be visited by a spacecraft. They are all of similar character, but their size relative to the planet varies widely, e.g., the planeto-centric distance to the BS nose ranges from about 1.4 R(sub V) for Venus to 88 R(sub J) or more for Jupiter. Comparisons are reviewed that show its location may be represented satisfactorily by a gasdynamic (GD) model, provided the properties of the solar wind and planetary magnetic field and ionosphere are known and used as input in the application. Factors that determine the location are discussed, and examples are presented to illustrate effects of their variation, including which part of a BS is influenced by a local variation of the magneto/ionopause (MIP) shape. The interplanetary magnetic field (IMF) has no influence on the BS location in the GD model, but is shown to have a small effect in corresponding solutions of the basic MHD model from which the GD model is derived as the limit for weak IMF. Nearly all GD and MHD solutions are for steady flow, but a solution for unsteady flow associated with the passage of an interplanetary shock is also presented. It shows that the BS moves rapidly from its initial to final location, e.g., in about minute for the earth. Since many changes in the solar wind occur over longer intervals, these results help explain the success of quasi-stationary solutions in modeling the BS in time-varying solar wind flows.

Spreiter, J. R.↗

Calculation of H2-He Flow with Nonequilibrium Ionization and Radiation: an Interim Report

The nonequilibrium ionization process in hydrogen-helium mixture behind a strong shock wave is studied numerically using the detailed ionization rate model developed recently by Park which accounts for emission and absorption of Lyman lines. The study finds that, once the avalanche ionization is started, the Lyman line is self-absorbed. The intensity variation of the radiation at 5145 Angstroms found by Leibowitz in a shock tube experiment can be numerically reproduced by assuming that ionization behind the shock wave prior to the onset of avalanche ionization is 1.3%. Because 1.3% initial ionization is highly unlikely, Leibowitz s experimental data is deemed questionable. By varying the initial electron density value in the calculation, the calculated ionization equilibration time is shown to increase approximately as inverse square-root of the initial electron density value. The true ionization equilibration time is most likely much longer than the value found by Leibowitz.

Furudate, Michiko↗

The influence of the computational mesh on accuracy for initial value problems with discontinuous or nonunique solutions

Discontinuous, or weak, solutions of the wave equation, the inviscid form of Burgers equation, and the time-dependent, two-dimensional Euler equations are studied. A numerical method of second-order accuracy in two forms, differential and integral, is used to calculate the weak solutions of these equations for several initial value problems, including supersonic flow past a wedge, a double symmetric wedge, and a sphere. The effect of the computational mesh on the accuracy of computed weak solutions including shock waves and expansion phenomena is studied. Modifications to the finite-difference method are presented which aid in obtaining desired solutions for initial value problems in which the solutions are nonunique.

Maccormack, R. W.↗

The X-ray spectra of the Vela and Puppis supernova remnants and the shock-wave model of supernova remnants

The structure and spectrum of the Vela and Puppis supernova remnants (SNRs) were observed in soft (0.1 to 1.5 keV) X rays from a rocket. The spectral data of both objects below 1.5 keV are consistent with thermal radiation from a hot plasma and are not consistent with a simple power-law function. Fitting the data to a Tucker-Koren model results in values of T = (4.3 plus or minus 0.3) x 1,000,000 K for Vela and T = (7 plus or minus 1) x 1,000,000 K for Puppis. The present data, earlier data from the Cygnus Loop, and results on Cas A and Tycho reported by others are considered within the context of a simple model in which X rays are thermally produced as a result of a shock wave expanding into the interstellar medium. X-ray data on the temperature, intrinsic diameter, and intrinsic luminosity of the five SNRs are used to compute the age, energy release, and initial interstellar density for each one. From measurements of X-ray absorption along the line of sight, and also from shock-wave model interpretations of temperature-diameter-luminosity relations, it is concluded that the average density of interstellar material is about 0.1 per cu cm between the sun and the nearest X-ray emitting SNR, Vela.

Gorenstein, P.↗

What's wrong with the NOL card gap test?

The card gap test, originally developed for explosives, was adopted for determining the shock sensitivity of propellants. Testing solid propellants resulted in the discovery of several deficiencies in the test. The effects of critical diameter, donor pulse duration, and initiator variability are discussed. Reproducibility of results is shown to be worse for solid propellants than is reported for explosives. The calibration to express test results as shock pressure is of questionable accuracy. Use of the test in its present form leads to erroneous conclusions.

Brown, B.↗

Shock propagation and the generation of magnetohydrodynamic wave fields in inhomogeneous molecular clouds

We develop a simple one-dimensional model for the interaction of a steady, thin, planar shock wave with a nonrigid cloud which may be in motion relative to the surrounding medium, and we apply the model to shocks impinging on, and propagating through, molecular clouds. Both 'adiabatic' (gamma = 5/3) and radiative (gamma = 1) shocks are considered and we allow for the presence of a uniform magnetic field directed either parallel or perpendicular to the shock normal. The former field orientation is equivalent to the hydrodynamic case, and the latter involves only fast MHD shocks. We focus on the manner in which such shocks can generate internal kinetic motions in the cloud on a range of size and density scales through the direct acceleration of cores and clumps by shocks transmitted into them and through the generation of an MHD wavefield via the reflection of the incident shock at clump boundaries. We find that stronger incident Mach numbers and smaller density contrasts lead to more efficient cloud acceleration, as do isothermal intercloud shocks and small intercloud magnetic field strengths. The acceleration efficiency is insensitive to the adiabatic index and the magnetic field strength in the cloud itself. For typical parameter choices, the direct acceleration of clouds and clumps by strong shocks is found to be substantial and could at least in part account for their observed velocity dispersions. If the shocks are moderately weak, the final velocity of the cloud is linearly related to its initial velocity, with higher acceleration giving shallower slopes (i.e., final velocity distributions which are less sensitive to the initial distribution). Compared to the kinetic energy of the postshock cloud, the energy given to the wavefield at each encounter is small, and the heating of the interclump medium by the dissipation of this wavefield is found to be insufficient to balance the cooling rate in the cloud as a whole (although it may be important in particular regions), even if this medium is warm, unless it is also extremely tenuous (n approximately less than 0.1/cu cm). Nevertheless, the correction for the velocity imparted to the cloud leads to a substantial increase in the critical incident Mach numner for wave emission over that reported by Spitzer for the rigid case. The implications of our model for shock-induced star formation are discussed briefly.

Miesch, Mark S.↗

Shock wave compression of iron-silicate garnet.

Shock wave compression data to over 650 kb are presented for single-crystal almandine garnet. The data indicate the initiation of a phase transformation near 200 kb. Total transition to the high-pressure polymorph occurs at approximately 300 kb. The elastic properties of the high-pressure phase are calculated from the metastable Hugoniot data by using the linear shock velocity-particle velocity relationships. The overall results obtained strongly suggest that upper mantle minerals are likely to occur in the ilmenite structure over a substantial part of the lower mantle.

Graham, E. K.↗

Temporal evolution of the solar wind and the formation of a standing shock

The temporal evolution of the solar wind from one steady state to another is explored when momentum deposition produces multiple critical points in the flow. It is shown that the wind always evolves in time to a new steady state compatible with the solution of the steady state equation of motion. However, for the same initial state and identical asymptotic momentum deposition rate, the temporal evolution pattern of the wind depends on the detailed time history of momentum addition and is therefore not unique. This feature plays an important role in the particular case when multiple (three in this study) steady states exist for identical boundary conditions; each one of these solutions is thus shown to be physically accessible. The details of the temporal evolution pattern of the wind reveal the formation of a shock discontinuity whenever the flow becomes supersonic at a critical point upstream from the initial critical point. If the flow remains supersonic at that inner critical point, the shock can become a standing one, depending on the strength and the temporal history of momentum addition. The results of this study indicate that the time scale required for the solar wind to evolve between steady states is of the order of 30-60 hours. Furthermore, the results also reveal the interesting and novel phenomenon that a standing shock is likely to develop in the inner solar wind flow within this time frame, in particular, in coronal hole regions with rapidly diverging geometries.

Habbal, S. R.↗

Coordinate perturbation and multiple scale in gasdynamics.

A simple way of applying the coordinate perturbation technique is developed which bypasses the usual transformations of the equations to perturbed coordinates, and which provides the appropriate coordinate stretching by simple algebraic manipulations on the nonuniformly valid expansion obtained by a straightforward expansion from the original equations. Interesting results are obtained by applying the procedure to two gasdynamical problems. In the first, the flowfield around a supersonic two-dimensional wing is determined up to third order, including a uniformly valid representation of the front shock shape. The second concerns the oscillations in a closed tube following an arbitrary initial disturbance, both when the two ends are closed, and when one of the two ends contains an oscillating piston.

Crocco, L.↗

Combined energy and pitch angle diffusion of pickup ions at Comet Halley

It is well known that cometary pickup ions, e.g., H2O(+), OH(+), O(+), CO(+), H(+), initially form a ring-beam distribution in the solar wind reference frame, which is highly unstable to the growth of MHD waves (such as ion-cyclotron waves). The low-frequency magnetic fluctuations (or waves), which were observed upstream of Comet Halley, cannot only pitch-angle-scatter the pickup ions so that the distribution becomes at least partially isotropized, but also stochastically accelerate the ions, resulting in the energetic ion populations observed in the vicinity of Comet Halley. Here, numerical solutions of the quasi-linear diffusion equation were used to investigate the cometary ion pickup process at Comet Halley. Both pitch angle and energy diffusion are taken into account. Many quasi-linear models of cometary pickup ions exist which involve one type of diffusion or the other but not both types at once. It is found that the pitch angle scattering occurs faster than the energy diffusion, as expected. Moreover, the results demonstrate that the distribution of accelerated energetic ions is more isotropic than that of ions which have just been picked up. In fact, the ion distribution function on the initial pickup shell is quite anisotropic, even close to the Comet-Halley bow shock.

Ye, Gang↗

The Galactic center wind

The combined effect of winds from a cluster of stars in the central 0.8 pc of the Galaxy is modeled as uniform power and mass input over the central region. The flow becomes supersonic outside the central region, and the expected decrease in pressure is in approximate accord with observations. The pressure variations on a larger scale suggest that the Galactic center wind passes through a shock front at a radius of a few pc, leading to a shocked wind bubble on a scale of 100 pc. The tangential magnetic field can come to dominate the pressure in the shocked wind flow even if the energy density of the magnetic field in the initial wind is only 0.1 percent of the wind kinetic energy density. The magnetic region produced in this way may be related to some of the apparently magnetized structures observed in the central region of the Galaxy.

Chevalier, Roger A.↗

Collisional and dynamical history of Gaspra

Interpretation of the impact record on Gaspra requires understanding of the effects of collisions on a target body of Gaspra's size and shape, recognition of impact features that may have different morphologies from craters on larger planets, and models of the geological processes that erase and modify impact features. Crater counts on the 140 sq km of Gaspra imaged at highest resolution by the Galileo spacecraft show a steep size-frequency distribution (cumulative power-law index near -3.5) from the smallest resolvable size (150 m diameter) up through the large feature (1.5 km diameter crater) of familiar crater-like morphology. In addition, there appear to be as many as eight roughly circular concavities with diameters greater than 3 km visible on the asteroid. If we restrict our crater counts to features with traditionally recognized crater morphologies, these concavities would not be included. However, if we define craters to include any concave structures that may represent local or regional damage at an impact size, then the larger features on Gaspra are candidates for consideration. Acceptance of the multi-km features as craters has been cautious for several reasons. First, scaling laws (the physically plausible algorithms for extrapolating from experimental data) indicate that Gaspra could not have sustained such large-crater-forming impacts without being disrupted; second, aside from concavity, the larger structures have no other features (e.g. rims) that can be identified with known impact craters; and third, extrapolation of the power-law size distribution for smaller craters predicts no craters larger than 3 km over the entire surface. On the other hand, recent hydrocode modeling of impacts shows that for given impact (albeit into a sphere), the crater size is much larger than given by scaling laws. Gaspra-size bodies can sustain formation of up to 8-km craters without disruption. Besides allowing larger impact craters, this result doubles the lifetime since the last catastrophic fragmentation event up to one billion years. Events that create multi-km craters also globally damage the material structure, such that regolith is produced, whether or not Gaspra 'initially' had a regolith, contrary to other models in which initial regolith is required in order to allow current regolith. Because the globally destructive shock wave precedes basin formation, crater size is closer to the large size extrapolated from gravity-scaling rather than the strength-scaling that had earlier been assumed for such small bodies. This mechanism may also help explain the existence of Stickney on Phobos. Moreover, rejection of the large concavities as craters based on unfamiliar morphology would be premature, because (aside from Stickney) we have no other data on such large impact structures on such a small, irregular body. The eight candidate concavities cover an area greater than that counted for smaller craters, because they are most apparent where small craters cannot be seen: on low resolution images and at the limb on high resolution images. We estimate that there are at least two with diameter greater than 4 km per 140 sq km, which would have to be accounted for in any model that claims these are impact craters.

Greenberg, R.↗

A Comparison of Computer Codes for the Propagation of Sonic Booms Through Realistic Atmospheres Utilizing Actual Acoustic Signatures

A numerical exercise to compare computer codes for the propagation of sonic booms through the atmosphere is reported. For the initial portion of the comparison, artificial, yet realistic, waveforms were numerically propagated through identical atmospheres. In addition to this comparison, one of these codes has been used to make preliminary predictions of the boom generated from a recent SR-71 flight. For the initial comparison, ground waveforms are calculated using four different codes or algorithms: (1) weak shock theory, an analytical prediction, (2) SHOCKN, a mixed time and frequency domain code developed at the University of Mississippi, (3) ZEPHYRUS, another mixed time and frequency code developed at the University of Texas, and (4) THOR, a pure time domain code recently developed at the University of Texas. The codes are described and their differences noted.

Chambers, James P.↗

Analytical Modeling of Pressure Wall Hole Size and Maximum Tip-to-Tip Crack Length for Perforating Normal and Oblique Orbital Debris Impacts

This report presents the results of a study whose objective was to develop first-principles-based models of hole size and maximum tip-to-tip crack length for a spacecraft module pressure wall that has been perforated in an orbital debris particle impact. The hole size and crack length models are developed by sequentially characterizing the phenomena comprising the orbital debris impact event, including the initial impact, the creation and motion of a debris cloud within the dual-wall system, the impact of the debris cloud on the pressure wall, the deformation of the pressure wall due to debris cloud impact loading prior to crack formation, pressure wall crack initiation, propagation, and arrest, and finally pressure wall deformation following crack initiation and growth. The model development has been accomplished through the application of elementary shock physics and thermodynamic theory, as well as the principles of mass, momentum, and energy conservation. The predictions of the model developed herein are compared against the predictions of empirically-based equations for hole diameters and maximum tip-to-tip crack length for three International Space Station wall configurations. The ISS wall systems considered are the baseline U.S. Lab Cylinder, the enhanced U.S. Lab Cylinder, and the U.S. Lab Endcone. The empirical predictor equations were derived from experimentally obtained hole diameters and crack length data. The original model predictions did not compare favorably with the experimental data, especially for cases in which pressure wall petalling did not occur. Several modifications were made to the original model to bring its predictions closer in line with the experimental results. Following the adjustment of several empirical constants, the predictions of the modified analytical model were in much closer agreement with the experimental results.

Schonberg, William P.↗

Summary of Transition Results From the F-16XL-2 Supersonic Laminar Flow Control Experiment

A variable-porosity suction glove has been flown on the F-16XL-2 aircraft to demonstrate the feasibility of this technology for the proposed High-Speed Civil Transport. Boundary-layer transition data on the titanium glove primarily have been obtained at speeds of Mach 2.0 and altitudes of 15,240-16,764 m (50,000-55,000 ft). The objectives of this flight experiment have been to achieve 0.50-0.60 chord laminar flow on a highly swept wing at supersonic speeds and to provide data to validate codes and suction design. The most successful laminar flow results have not been obtained at the glove design point, a speed of Mach 1.9 at an altitude of 15,240 m (50,000 ft); but rather at a speed of Mach 2.0 and an altitude of 16,154 m (53,000 ft). Laminar flow has been obtained to more than 0.46 wing chord at a Reynolds number of 22.7 x 10(exp 6). A turbulence diverter has been used to initially obtain a laminar boundary layer at the attachment line. A lower-surface shock fence was required to block an inlet shock from the wing leading edge. This paper discusses research variables that directly impact the ability to obtain laminar flow and techniques to correct for these variables.

Marshall, Laurie A.↗

Evaluation of Thermo-Mechanical Stability of COTS Dual-Axis MEMS Accelerometers for Space Applications

Microelectromechanical systems in MEMS is one of the fastest growing technologies in microelectronics, and is of great interest for military and aerospace applications. Accelerometers are the earliest and most developed representatives of MEMS. First demonstrated in 1979, micromachined accelerometers were used in automobile industry for air bag crash- sensing applications since 1990. In 1999, N4EMS accelerometers were used in NASA-JPL Mars Microprobe. The most developed accelerometers for airbag crash- sensing are rated for a full range of +/- 50 G. The range of sensitivity for accelerometers required for military or aerospace applications is much larger, varying from 20,000 G (to measure acceleration during gun and ballistic munition launches), and to 10(exp -6) G, when used as guidance sensors (to measure attitude and position of a spacecraft). The presence of moving parts on the surface of chip is specific to MEMS, and particularly, to accelerometers. This characteristic brings new reliability issues to micromachined accelerometers, including cyclic fatigue cracking of polysilicon cantilevers and springs, mechanical stresses that are caused by packaging and contamination in the internal cavity of the package. Studies of fatigue cracks initiation and growth in polysilicon showed that the fatigue damage may influence MEMS device performance, and the presence of water vapor significantly enhances crack initiation and growth. Environmentally induced failures, particularly, failures due to thermal cycling and mechanical shock are considered as one of major reliability concerns in MEMS. These environmental conditions are also critical for space applications of the parts. For example, the Mars pathfinder mission had experienced 80 mechanical shock events during the pyrotechnic separation processes.

Sharma, Ashok K.↗