SEARCH · Engineering Papers
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.
A two-dimensional hydrodynamic numerical model has been used to study the motion of the termination shock in response to an 11 year variation in the solar wind ram pressure. For a total variation in the ram pressure by a factor of 2, a termination shock at 89 AU moves inward and outward about +/- 8 percent of its distance with a typical velocity of 12 km/s.
The trajectories of ions reflected or leaked upstream from the earth's bow shock and subject solely to the Lorentz force in a steady interplanetary magnetic field B and the V x B electric field are studied theoretically. Expressions are obtained for the guiding center motion and gyromotion in a frame (the Hoffman-Teller frame) moving parallel to the shock surface with sufficient speed to transform the incident solar wind velocity into motion entirely along the interplanetary magnetic field. Equations are derived which transform these motions back to the observer's frame. The predicted upstream motions for four different source models for upstream ions are compared using these expressions: magnetic moment-conserving reflection of solar wind ions, specular reflection of solar wind ions, magnetic moment-conserving leakage of magnetosheath ions, and leakage of magnetosheath ions parallel to the shock normal.
Here, this paper overviews experimental and numerical results on acceleration of narrow energy spread ion beams by an electrostatic collisionless shockwave driven by 1 μm (Omega EP) and 10 μm (UCLA Neptune Laboratory) lasers in near critical density CH and He plasmas, respectively. Shock waves in CH targets produced high-energy ~50 MeV protons (ΔE/E of ≤30%) and 314 MeV C 6+ ions (ΔE/E of ≤10%). Observation of acceleration of both protons and carbon ions to similar velocities is consistent with reflection of particles off the moving potential of a shock front. For shocks driven by CO 2 laser in a gas jet, ~30 MeV peak in He ion spectrum was detected. Particle-in-cell simulations indicate that regardless of the target further control over its density profile is needed for optimization of accelerated ion beams in part of energy spread, yield and maximum kinetic energy.
A large set of bow shock crossings (i.e., 1392) observed by 17 spacecraft has been used to explore the three-dimensional shape and location of the Earth's bow shock and its dependence on solar wind and interplanetary magnetic field (IMF) conditions. This study investigates deviations from gas dynamic flow models associated with the magnetic terms in the magnetohydrodynamic (MHD) equations. Empirical models predicting the statistical position and shape of the bow shock for arbitrary values of the solar wind pressure, IMF, and Alfvenic Mach number (M(sub A)) have been derived. The resulting data set has been used to fit three-dimensional bow shock surfaces and to explore the variations in these surfaces with sonic (M(sub S)), Alfvenic (M(sub A)) and magnetosonic (M(sub MS)) Mach numbers. Analysis reveals that among the three Mach numbers, M(sub A) provides the best ordering of the least square bow shock curves. The subsolar shock is observed to move Earthward while the flanks flare outward in response to decreasing M(sub A); the net change represents a 6-10% effect. Variations due to changes in the IMF orientation were investigated by rotating the crossings into geocentric interplanetary medium coordinates. Past studies have suggested that the north-south extent of the bow shock surface exceeds the east-west dimension due to asymmetries in the fast mode Mach cone. This study confirms such a north-south versus east-west asymmetry and quantifies its variation with M(sub S), M(sub A), M(sub MS), and IMF orientation. A 2-7% effect is measured, with the asymmetry being more pronounced at low Mach numbers. Combining the bow shock models with the magnetopause model of Roelof and Sibeck (1993), variations in the magnetosheath thickness at different local times are explored. The ratio of the bow shock size to the magnetopause size at the subpolar point is found to be 1.46; at dawn and dusk, the ratios are found to be 1.89 and 1.93, respectively. The subsolar magnetosheath thickness is used to derive the polytropic index gamma according to the empirical relation of Spreiter et al. (1966). The resulting gamma = 2.3 suggests the empirical formula is inadequate to describe the MHD interaction between the solar wind and the magnetosphere.
The effect of area ratio variation on the performance of a supersonic Hall MHD duct is investigated. Results indicate that for a given combustion pressure there exists an area ratio below which the power generating region of the duct is shock free and the power output increases linearly with the square of the magnetic field. For area ratios greater than this, a shock forms in the power generating region which moves upstream with increasing magnetic field strength resulting in a less rapid raise in the power output. The shock can be moved downstream by either increasing the combustion pressure of decreasing the exhaust pressure. The influence of these effects upon duct performance is presented.
The effect of area ratio variation on the performance of a supersonic Hall MHD duct showed that for a given combustion pressure there exists an area ratio below which the power generating region of the duct is shock free and the power output increases linearly with the square of the magnetic field. For area ratios greater than this, a shock forms in the power generating region which moves upstream with increasing magnetic field strength resulting in a less rapid raise in the power output. The shock can be moved downstream by either increasing the combustion pressure or decreasing the exhaust pressure. The influence of these effects upon duct performance is presented in this paper.
In this paper, we present an experiment that explores the plasma dynamics of a 7μm diameter carbon wire after being irradiated with a near-relativistic-intensity short pulse laser. Using an x-ray free electron laser pulse to measure the small angle x-ray scattering signal, we observe that the scattering surface is bent and prone to instability over tens of picoseconds. The dynamics of this process are consistent with the presence of a sharp, propagating shock front inside the wire, moving at a speed close to the hole boring velocity or that expected from a thermal shock at a few tens of Mbar.
We present the first images of the pulsar wind nebula (PWN) MSH 15−52 in the hard X-ray band (8 keV), as measured with the Nuclear Spectroscopic Telescope Array (NuSTAR). Overall, the morphology of the PWN as measured by NuSTAR in the 3-7 keV band is similar to that seen in Chandra high-resolution imaging. However, the spatial extent decreases with energy, which we attribute to synchrotron energy losses as the particles move away from the shock. The hard-band maps show a relative deficit of counts in the northern region toward the RCW 89 thermal remnant, with significant asymmetry. We find that the integrated PWN spectra measured with NuSTAR and Chandra suggest that there is a spectral break at 6 keV, which may be explained by a break in the synchrotron emitting electron distribution at approximately 200 TeV and/or imperfect cross calibration. We also measure spatially resolved spectra, showing that the spectrum of the PWN softens away from the central pulsar B1509−58, and that there exists a roughly sinusoidal variation of spectral hardness in the azimuthal direction. We discuss the results using particle flow models. We find non-monotonic structure in the variation with distance of spectral hardness within 50 of the pulsar moving in the jet direction, which may imply particle and magnetic-field compression by magnetic hoop stress as previously suggested for this source. We also present two-dimensional maps of spectral parameters and find an interesting shell-like structure in the N(sub H) map. We discuss possible origins of the shell-like structure and their implications.
Experimental studies are conducted to examine the utilization of transpiration cooling to reduce the peak-heating loads in areas of shock/shock interaction. Smooth and transpiration-cooled nosetip models, 12 inches in diameter, were employed in these studies, which focused on defining the pressure distributions and heat transfer in type III and IV interaction areas. Transpiration cooling was determined to significantly increase the size of the shock layer and to move the peak-heating point around the body. A transpiration-cooling rate of more than 30 percent of the freestream maximum flux did not lower the peak-heating level more than 10 percent, but the integrated heating loads were reduced.
We report on the radio-emission characteristics of 222 interplanetary (IP) shocks detected by spacecraft at Sun-Earth L1 during solar cycle 23 (1996 to 2006, inclusive). A surprisingly large fraction of the IP shocks (approximately 34%) was radio quiet (RQ; i.e., the shocks lacked type II radio bursts). We examined the properties of coronal mass ejections (CMEs) and soft X-ray flares associated with such RQ shocks and compared them with those of the radio-loud (RL) shocks. The CMEs associated with the RQ shocks were generally slow (average speed approximately 535 km/s) and only approximately 40% of the CMEs were halos. The corresponding numbers for CMEs associated with RL shocks were 1237 km/s and 72%, respectively. Thus, the CME kinetic energy seems to be the deciding factor in the radio-emission properties of shocks. The lower kinetic energy of CMEs associated with RQ shocks is also suggested by the lower peak soft X-ray flux of the associated flares (C3.4 versus M4.7 for RL shocks). CMEs associated with RQ CMEs were generally accelerating within the coronagraph field of view (average acceleration approximately +6.8 m/s (exp 2)), while those associated with RL shocks were decelerating (average acceleration approximately 3.5 m/s (exp 2)). This suggests that many of the RQ shocks formed at large distances from the Sun, typically beyond 10 Rs, consistent with the absence of metric and decameter-hectometric (DH) type II radio bursts. A small fraction of RL shocks had type II radio emission solely in the kilometric (km) wavelength domain. Interestingly, the kinematics of the CMEs associated with the km type II bursts is similar to those of RQ shocks, except that the former are slightly more energetic. Comparison of the shock Mach numbers at 1 AU shows that the RQ shocks are mostly subcritical, suggesting that they were not efficient in accelerating electrons. The Mach number values also indicate that most of these are quasi-perpendicular shocks. The radio-quietness is predominant in the rise phase and decreases through the maximum and declining phases of solar cycle 23. About 18% of the IP shocks do not have discernible ejecta behind them. These shocks are due to CMEs moving at large angles from the Sun-Earth line and hence are not blast waves. The solar sources of the shock-driving CMEs follow the sunspot butterfly diagram, consistent with the higher-energy requirement for driving shocks.
Errors found in the original edition are corrected. Refinement was made in procedures for solving the conservation relations for an incident (moving), standing, and reflected normal shock, as well as in computational methods for determining thermochemical-equilibrium hydrogen properties. A six-species hydrogen model replaces the original four-species model, and the heat of formation and spectroscopic constants used in this six-species model are listed in appendix A. In appendix B, comparisons are made between a number of methods for determining equilibrium thermodynamic properties for hydrogen for several values of pressure and temperatures to 50000 K. A comparison is also performed between the present method and a second method for determining thermodynamic properties and flow velocity behind an incident shock into pure hydrogen and behind a reflected shock.
Equations are presented for the growth or decay rates of moving and deforming vortex sheets (gust fronts) and shock waves. The equations exhibit a cross-coupling between these two types of surface discontinuity.
The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of an axisymmetric external compression supersonic inlet is being developed. The model utilizes compressible flow computational fluid dynamics to model the internal inlet segment as well as the external inlet portion between the cowl lip and normal shock, and compressible flow relations with flow propagation delay to model the oblique shocks upstream of the normal shock. The external compression portion between the cowl-lip and the normal shock is also modeled with leaking fluxes crossing the sonic boundary, with a moving CFD domain at the normal shock boundary. This model has been verified in steady state against tunnel inlet test data and it s a first attempt towards developing a more comprehensive model for inlet dynamics.
Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.
It is proposed that in loop coronal transients an erupting loop moves away from the solar surface, with a velocity exceeding the local Alfven speed, pushing against the overlying magnetic fields and driving a shock in the front of the moving part of the loop. Lower hybrid waves are excited at the shock front and propagate radially toward the center of the loop with phase velocity along the magnetic field that exceeds the thermal velocity. The lower hybrid waves stochastically accelerate the tail of the electron distribution inside the loop. The manner in which the accelerated electrons are trapped in the moving loop are discussed, and their radiation signature is estimated. It is suggested that plasma radiation can explain the power observed in stationary and moving type IV bursts.
In connection with investigations conducted to determine the influence of water vapor on experiments in wind tunnels, the question arose as to what changes due to vapor condensation might be expected in airfoil measurements. Density measurements on circular-arc airfoils aided by an interferometer in choked tunnels with parallel walls show that increasing humidity produces increasing changes in the flow field. The flow becomes nonstationary at high humidity. At the airfoil, however, the influence of the condensation is only felt, inasmuch as the shock bounding the local supersonic region moves upstream with increasing humidity while its intensity decreases. The density distribution upstream of the shock remains unchanged. Even if the flow becomes nonstationary in the vicinity of the airfoil, no changes occur at the airfoil.
An examination is presented of the association of solar radio bursts of spectral type II and coronal shocks with solar flare ejecta observed in H-alpha, the green coronal line, and white-light coronagraphs. It is suggested that fast-moving optical coronal transients should for the most part be identified with piston-type phenomena well behind the outward-traveling shock waves that generate type II radio bursts. A general model is presented which relates type II radio bursts and coronal shocks to optically observed ejecta and consists of three main velocity regimes: (1) a quasi-hemispherical shock wave moving outward from the flare at speeds of 1000-2000 km/sec and Alfven Mach number of about 1.5; (2) the velocity of the piston driving the shock, on the order of 0.8 that of the shock; and (3) the regime of the slower-moving H-alpha ejecta, with velocities of 300-500 km/sec.
A shock that is continuously driven by a moving exciter will propagate at a speed that depends on the exciter speed. We obtained this dependence experimentally, in a strongly coupled dusty plasma that was prepared as a single two-dimensional layer of charged microparticles. Attaining this result required an experimental advance, developing a method of driving a shock continuously, which we did using an exciter moving at a constant supersonic speed, analogous to a piston in a cylinder. The resulting compressional pulse was a shock that propagated steadily without weakening, ahead of the moving exciter. In this work, we compare our experimental results to an empirical form $M_{\text{shock}} = 1 + sM_{\text{exciter}}$, and to the prediction of a recent simulation.