SEARCH · Engineering Papers

Results for “Shock Physics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗

NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗

NASA NTRS2023

The Faraday Effect Tracker of Coronal and Heliospheric Structure (FETCH) Instrument

There continue to be open questions regarding the solar wind and coronal mass ejections (CMEs). For example: how do magnetic fields within CMEs and corotating/stream interaction regions (CIRs/SIRs) evolve in the inner heliosphere? What is the radially distributed magnetic profile of shock-driving CMEs? What is the internal magnetic structure of CMEs that cause magnetic storms? It is clear that these questions involve the magnetic configurations of solar wind and transient interplanetary plasma structures, for which we have limited knowledge. In order to better understand the origin of the magnetic field variability in steady-state structures and transient events, it is necessary to probe the magnetic field in Earth-directed structures/disturbances. This is the goal of the Multiview Observatory for Solar Terrestrial Science (MOST) mission (Gopalswamy et al., 2022). For MOST to answer the aforementioned questions, we propose the instrument concept of the Faraday Effect Tracker of Coronal and Heliospheric structures (FETCH), a simultaneous quad-line-of-sight polarization radio remote-sensing instrument. With FETCH, spacecraft radio beams passing through the Sun–Earth line offer the possibility of obtaining information of plasma conditions via analysis of radio propagation effects such as Faraday rotation and wave dispersion, which provide information of the magnetic field and total electron content (TEC). This is the goal of the FETCH instrument, one of ten instruments proposed to be hosted on the MOST mission. The MOST mission will provide an unprecedented opportunity to achieve NASA’s heliophysics science goal to “explore and characterize the physical processes in the space environment from the Sun” (Gopalswamy et al., 2022).

Solar corona↗

NASA NTRS1994

Fine structure in plasma waves and radiation near the plasma frequency in Earth's foreshock

Novel observations are presented of intrunsic fine structure in the frequency spectrum of electomagnetic (EM) radiation and plasma waves near the electron plasma frequency f(sub p) during a period of unusually high interplanetary magnetic field strength. Measured using the wideband receiver on the International Sun-Earth Explorer (ISEE) 1 spacecraft, fine-structured emissions are observed both in the solar wind and the foreshock, The fine structure is shown to correspond to emissions spaced above f(sub p) near half harmonies of the electon cyclotron frequency f(sub ce), i.e., near f(sub p) + nf(sub ce)/2. These appear to be the first space physics observations of emissions spaced by f(sub ce)/2. Indirect but strong arguments are used to discriminate between EM and electrostatic (ES) signals, to identify whether ISEE 1 is in the solar wind or the foreshock, and to determine the relative frequencies of the emissions and the local f(sub p). The data are consistent with generation of the ES and EM emissions in the foreshock, with subsequent propagation of the EM emissions into the solar wind. It remains possible that some emissions currently identified as ES have significant EM character. The ES and EM emisions often merge into one another with minimal changes in frequency, arguing that their source regions and generation mechanisms are related and imposing significant constraints on theories. The f(sub ce)/2 ES and EM fine structures observed may be intrinsic to the emission mechanisms or to superposition of two series of signals with f(sub ce) spacing that differ in starting frequency by f(sub ce)/2. Present theories for nonlinear wave coupling processes, cyclotron maser emission, and other linear instability processes are all unable to explain multiple EM and/or ES components spaced by approximately f(sub ce)/2 above f(sub p) for f(sub p)/f(sub ce) much greater than 1 and typical for shock beams parameters. Suitable avenues for further theoretical research are identified. Empirically, the observed fine structures appear very similar to those in split bnad and multiple-lane type II solar radio bursts; interpretation of both these type II fine structures in terms of f(sub ce)/2 splitting is suggested, thereby supporting and generalizing a suggestion by Wild (1950). A possible application to continuum radiation is mentioned. The ubiquity of these fine structures in the Earth's f(sub p) radiation and foreshock waves remains unknown. Only the ISEE 1 wideband receiver has sufficient frequency resolution (approximately less than or equal to 100 Hz) to perform a dedicated search. Further study of the ubiquity of these fine structures, of how reliably the splitting corresponds to f(sub ce)/2, and of the other interpretations above is necessary.

Cairns, Iver H.↗

NASA NTRS1994

Wind Prelaunch Mission Operations Report (MOR)

The National Aeronautics and Space Administration (NASA) Wind mission is the first mission of the Global Geospace Science (GGS) initiative. The Wind laboratory will study the properties of particles and waves in the region between the Earth and the Sun. Using the Moon s gravity to save fuel, dual lunar swing-by orbits enable the spacecraft to sample regions close to and far from the Earth. During the three year mission, Wind will pass through the bow shock of Earth's magnetosphere to begin a thorough investigation of the solar wind. Mission objectives require spacecraft measurements in two orbits: lunar swing- by ellipses out to distances of 250 Earth radii (RE) and a small orbit around the Lagrangian point L-l that remains between the Earth and the Sun. Wind will be placed into an initial orbit for approximately 2 years. It will then be maneuvered into a transition orbit and ultimately into a halo orbit at the Earth-Sun L-l point where it will operate for the remainder of its lifetime. The Wind satellite development was managed by NASA's Goddard Space Flight Center with the Martin Marietta Corporation, Astro-Space Division serving as the prime contractor. Overall programmatic direction was provided by NASA Headquarters, Office of Space Science. The spacecraft will be launched under a launch service contract with the McDonnell Douglas Corporation on a Delta II Expendable Launch Vehicle (ELV) within a November l-l4, 1994 launch window. The Wind spacecraft carries six U.S. instruments, one French instrument, and the first Russian instrument ever to fly on an American satellite. The Wind and Polar missions are the two components of the GGS Program. Wind is also the second mission of the International Solar Terrestrial Physics (ISTP) Program. The first ISTP mission, Geotail, is a joint project of the Institute of Space and Astronautical Science of Japan and NASA which launched in 1992. The Wind mission is planned to overlap Geotail by six months and Polar by one year. The Wind and Polar missions, together with the Geotail mission (launched on July 24, 1992) and supporting equatorial measurements, will provide simultaneous data to enable the study of solar wind input to the magnetosphere and key elements of the magnetospheric response: ring current energy storage, geomagnetic tail energy storage, and ionospheric energy input.

Source record↗

NASA NTRS2000

A Collaborative Analysis Tool for Integrated Hypersonic Aerodynamics, Thermal Protection Systems, and RBCC Engine Performance for Single Stage to Orbit Vehicles

Presented is a computer-based tool that connects several disciplines that are needed in the complex and integrated design of high performance reusable single stage to orbit (SSTO) vehicles. Every system is linked to every other system, as is the case of SSTO vehicles with air breathing propulsion, which is currently being studied by NASA. An RBCC propulsion system integrates airbreathing and rocket propulsion into a single engine assembly enclosed within a cowl or duct. A typical RBCC propulsion system operates as a ducted rocket up to approximately Mach 3. Then there is a transition to a ramjet mode for supersonic-to-hypersonic acceleration. Around Mach 8 the engine transitions to a scramjet mode. During the ramjet and scramjet modes, the integral rockets operate as fuel injectors. Around Mach 10-12 (the actual value depends on vehicle and mission requirements), the inlet is physically closed and the engine transitions to an integral rocket mode for orbit insertion. A common feature of RBCC propelled vehicles is the high degree of integration between the propulsion system and airframe. At high speeds the vehicle forebody is fundamentally part of the engine inlet, providing a compression surface for air flowing into the engine. The compressed air is mixed with fuel and burned. The combusted mixture must be expanded to an area larger than the incoming stream to provide thrust. Since a conventional nozzle would be too large, the entire lower after body of the vehicle is used as an expansion surface. Because of the high external temperatures seen during atmospheric flight, the design of an airbreathing SSTO vehicle requires delicate tradeoffs between engine design, vehicle shape, and thermal protection system (TPS) sizing in order to produce an optimum system in terms of weight (and cost) and maximum performance. To adequately determine the performance of the engine/vehicle, the Hypersonic Flight Inlet Model (HYFIM) module was designed to interface with the RBCC engine model. HYFIM performs the aerodynamic analysis of forebodies and inlet characteristics of RBCC powered SSTO launch vehicles. HYFIM is applicable to the analysis of the ramjet/scramjet engine operations modes (Mach 3-12), and provides estimates of parameters such as air capture area, shock-on-lip Mach number, design Mach number, compression ratio, etc., based on a basic geometry routine for modeling axisymmetric cones, 2-D wedge geometries. HYFIM also estimates the variation of shock layer properties normal to the forebody surface. The thermal protection system (TPS) is directly linked to determination of the vehicle moldline and the shaping of the trajectory. Thermal protection systems to maintain the structural integrity of the vehicle must be able to mitigate the heat transfer to the structure and be lightweight. Herein lies the interdependency, in that as the vehicle's speed increases, the TPS requirements are increased. And as TPS masses increase the effect on the propulsion system and all other systems is compounded. The need to analyze vehicle forebody and engine inlet is critical to be able to design the RBCC vehicle. To adequately determine insulation masses for an RBCC vehicle, the hypersonic aerodynamic environment and aeroheating loads must be calculated and the TPS thicknesses must be calculated for the entire vehicle. To accomplish this an ascent or reentry trajectory is obtained using the computer code Program to Optimize Simulated Trajectories (POST). The trajectory is then used to calculate the convective heat rates on several locations on the vehicles using the Miniature Version of the JA70 Aerodynamic Heating Computer Program (MINIVER). Once the heat rates are defined for each body point on the vehicle, then insulation thicknesses that are required to maintain the vehicle within structural limits are calculated using Systems Improved Numerical Differencing Analyzer (SINDA) models. If the TPS masses are too heavy for the performance of the vehicle the process may be repeated altering the trajectory or some other input to reduce the TPS mass. E-PSURBCC is an "engine performance" model and requires the specification of inlet air static temperature and pressure as well as Mach number (which it pulls from the HYFIM and POST trajectory files), and calculates the corresponding stagnation properties. The engine air flow path geometry includes inlet, a constant area section where the rocket is positioned, a subsonic diffuser, a constant area afterburner, and either a converging nozzle or a converging-diverging nozzle. The current capabilities of E-PSURBCC ejector and ramjet mode treatment indicated that various complex flow phenomena including multiple choking and internal shocks can occur for combinations of geometry/flow conditions. For a given input deck defining geometry/flow conditions, the program first goes through a series of checks to establish whether the input parameters are sound in terms of a solution path. If the vehicle/engine performance fails mission goals, the engineer is able to collaboratively alter the vehicle moldline to change aerodynamics, or trajectory, or some other input to achieve orbit. The problem described is an example of the need for collaborative design and analysis. RECIPE is a cross-platform application capable of hosting a number of engineers and designers across the Internet for distributed and collaborative engineering environments. Such integrated system design environments allow for collaborative team design analysis for performing individual or reduced team studies. To facilitate the larger number of potential runs that may need to be made, RECIPE connects the computer codes that calculate the trajectory data, aerodynamic data based on vehicle geometry, heat rate data, TPS masses, and vehicle and engine performance, so that the output from each tool is easily transferred to the model input files that need it.

Stanley, Thomas Troy↗

NASA NTRS1996

CFD Study of Turbo-Ramjet Interactions in Hypersonic Airbreathing Propulsion System

Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually involve turbo-ramjet configurations. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. In the first phase of our study, an over/under nozzle configuration was analyzed. The two plumes from the turbojet and ramjet interact at the end of a common 2-D cowl, where they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. For the problem analyzed, the turbojet engine operates at a higher nozzle pressure ratio than the ramjet, causes the turbojet plume overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data were used to compare with a computational fluid dynamics (CFD) study using the PARC2D code. The CFD results were in good agreement with both static pressure distributions on the cowl separator and on nozzle walls. The thrust coefficients were also in reasonable agreement. In addition, inviscid relationships were developed around the confluence point, where the two exhaust jets meet, and these results compared favorably with the CFD results. In the second phase of our study, a 3-D CFD solution was generated to compare with the 2-D solution. The major difference between the 2-D and 3-D solutions was the interaction of the shock waves, generated by the plume interactions, on the sidewall. When a shock wave interacts with a sidewall and sidewall boundary layer, it is called a glancing shock sidewall interaction. These interactions entrain boundary layer flow down the shockline into a vortical flow pattern. The 3-D plots show the streamlines being entrained down the shockline. The pressure of the flow also decreases slightly as the sidewall is approached. Other difference between the 2-D and 3-D solutions were a lowering of the nozzle thrust coefficient value from 0.9850 (2-D) to 0.9807 (3-D), where the experimental value was 0.9790. In the third phase of our study, a different turbo-ramjet configuration was analyzed. The confluence of a supersonic turbojet and a subsonic ramjet in the turbine based combined-cycle (TBCC) propulsion system was studied by a 2-D CFD code. In the analysis, Mach 1.4 primary turbojet was mixed with the subsonic ramjet secondary flow in an ejector mode operation. Reasonable agreements were obtained with the supplied I-D TBCC solutions. For low downstream backpressure, the Fabri choke condition (Break-Point condition) was observed in the secondary flow within mixing zone. For sufficient high downstream backpressure, the Fabri choke no longer exist, the ramjet flow was reduced and the ejector flow became backpressure dependent. Highly non-uniform flow at ejector exit were observed, indicated that for smooth downstream combustion, the mixing of the two streams probably required some physical devices.

Chang, Ing↗

NASA NTRS1994

What we have learned about Mars from SNC meteorites

The SNC meteorites are thought to be igneous martian rocks, based on their young crystallization ages and a close match between the composition of gases implanted in them during shock and the atmosphere of Mars. A related meteorite, ALH84001, may be older and thus may represent ancient martian crust. These petrologically diverse basalts and ultramafic rocks are mostly cumulates, but their parent magmas share geochemical and radiogenic isotopic characteristics that suggest they may have formed by remelting the same mantle source region at different times. Information and inferences about martian geology drawn from these samples include the following: Planetary differentiation occured early at approximately 4.5 GA, probably concurrently with accretion. The martian mantle contains different abundances of moderately volatile and siderophile elements and is more Fe-rich than that of the Earth, which has implications for its mineralogy, density, and origin. The estimated core composition has a S abundance near the threshold value for inner core solidification. The former presence of a core dynamo may be suggested by remanent magnetization in Shergottite-Nakhlite-Chassignite (SNC) meteorites, although these rocks may have been magnetized during shock. The mineralogy of martian surface units, inferred from reflectance spectra, matches that of basaltic shergottites, but SNC lithologies thought to have crystallized in the subsurface are not presently recognized. The rheological properties of martian magmas are more accurately derived form these metorites than from observations of martian flow morphology, although the sampled range of magma compositions islimited. Estimates of planetary water abundance and the amount of outgassed water based on these meteorites are contridictory but overlap estimates based on geological observations and atmospheric measurements. Stable isotope measurements indicate that the martian hydrosphere experienced only limited exchange with the lithosphere, but it is in isotopic equilibrium with the atmosphere and has been since 1.3 Ga. The isotopically heavy atmosphere/hydrosphere composition deduced from these rocks reflects a loss process more severe than current atmospheric evolution models, and the occurence of carbonates in SNC meteorites suggest that they, rather than scapolite or hydrous carbonates, are the major crustal sink for CO2. Weathering products in SNC meteorites support the idea of limited alteration of the lithosphere by small volumes of saline, CO2-bearing water. Atmospheric composition and evolution are further constrained by noble gases in these meteorites, although Xe and Kr isotopes suggest different origins for the atmosphere. Planetary ejection of these rocks has promoted an advance in the understanding of impact physics, which has been accomplished by a model involving spallation during large cratering events. Ejection of all the SNC meteorites (except ALH84001) in one or two events may provide a plausible solution to most constraints imposed by chronology, geochemistry, and cosmic ray exposure, although problems remain with this scenario; ALH84001 may represent older martian crust sampled during a separate impact.

Mcsween, Harry Y., Jr.↗

NASA NTRS2012

Multi-Sensors Observations of Pre-Earthquake Signals. What We Learned from the Great Tohoku Earthquake?

The lessons learned from the Great Tohoku EQ (Japan, 2011) will affect our future observations and an analysis is the main focus of this presentation. Multi-sensors observations and multidisciplinary research is presented in our study of the phenomena preceding major earthquakes Our approach is based on a systematic analysis of several physical and environmental parameters, which been reported by others in connections with earthquake processes: thermal infrared radiation; temperature; concentration of electrons in the ionosphere; radon/ion activities; and atmospheric temperature/humidity [Ouzounov et al, 2011]. We used the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model, one of several possible paradigms [Pulinets and Ouzounov, 2011] to interpret our observations. We retrospectively analyzed the temporal and spatial variations of three different physical parameters characterizing the state of the atmosphere, ionosphere the ground surface several days before the March 11, 2011 M9 Tohoku earthquake Namely: (i) Outgoing Long wave Radiation (OLR) measured at the top of the atmosphere; (ii) Anomalous variations of ionospheric parameters revealed by multi-sensors observations; and (iii) The change in the foreshock sequence (rate, space and time); Our results show that on March 8th, 2011 a rapid increase of emitted infrared radiation was observed and an anomaly developed near the epicenter with largest value occurring on March 11 at 07.30 LT. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. Starting from this day in the lower ionosphere there was also observed an abnormal TEC variation over the epicenter. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. We use the Japanese GPS network stations and method of Radio Tomography to study the spatiotemporal structure of ionospheric perturbations, and to distinguish ionospheric responses to processes of EQ preparation against the effects of other factors. The 2-D snapshots of the electron density over Japan showed abnormal increase over the maximum stress during the night, a few hours before the main shock. Our results from recording atmospheric and ionospheric conditions during the earthquake indicate the presence of anomalies in the atmosphere and ionospheres occurring consistently over regions of maximum stress near the epicenter. Due to their long duration (hours and days) and spatial appearance (only over the Sendai region) these results do not appear to be caused by meteorological or magnetic activity. They reveal the existence of atmospheric and ionospheric phenomena occurring prior to the earthquake, which indicates new evidence of a distinct coupling between the lithosphere and atmosphere/ionosphere. Similar results have been reported before the catastrophic events in Chile (M8.8, 2010), Italy (M6.3, 2009) and Sumatra (M9.3, 2004).

Ouzonounov, D.↗

NASA NTRS2022

The X-Ray Polarization View of Mrk 421 in an Average Flux State as Observed By the Imaging X-Ray Polarimetry Explorer

Particle acceleration mechanisms in supermassive black hole jets, such as shock acceleration, magnetic reconnection, and turbulence, are expected to have observable signatures in the multiwavelength polarization properties of blazars. The recent launch of the Imaging X-Ray Polarimetry Explorer (IXPE) enables us, for the first time, to use polarization in the X-ray band (2–8 keV) to probe the properties of the jet synchrotron emission in high-synchrotron-peaked BL Lac objects (HSPs). We report the discovery of X-ray linear polarization (degree Π x = 15% ± 2% and electric vector position angle ψ x = 35° ± 4°) from the jet of the HSP Mrk 421 in an average X-ray flux state. At the same time, the degree of polarization at optical, infrared, and millimeter wavelengths was found to be lower by at least a factor of 3. During the IXPE pointing, the X-ray flux of the source increased by a factor of 2.2, while the polarization behavior was consistent with no variability. The higher level of Π x compared to longer wavelengths, and the absence of significant polarization variability, suggest a shock is the most likely X-ray emission site in the jet of Mrk 421 during the observation. The multiwavelength polarization properties are consistent with an energy-stratified electron population, where the particles emitting at longer wavelengths are located farther from the acceleration site, where they experience a more disordered magnetic field.

acceleration of particles↗

NASA NTRS2010

Effects of Temperature on the Performance and Stability of Recent COTS Silicon Oscillators

Silicon oscillators have lately emerged to serve as potential replacement for crystal and ceramic resonators to provide timing and clock signals in electronic systems. These semiconductor-based devices, including those that are based on MEMS technology, are reported to be resistant to vibration and shock (an important criteria for systems to be deployed in space), immune to EMI, consume very low current, require few or no external components, and cover a wide range of frequency for analog and digital circuits. In this work, the performance of five recently-developed COTS silicon oscillator chips from different manufacturers was determined within a temperature range that extended beyond the individual specified range of operation. In addition, restart capability at extreme temperatures, i.e. power switched on while the device was soaking at extreme (hot or cold) temperature, and the effects of thermal cycling under a wide temperature range on the operation of these silicon oscillators were also investigated. Performance characterization of each oscillator was obtained in terms of its output frequency, duty cycle, rise and fall times, and supply current at specific test temperatures. The five different oscillators tested operated beyond their specified temperature region, with some displaying excellent stability throughout the whole test temperature range. Others experienced some instability at certain temperature test points as evidenced by fluctuation in the output frequency. Recovery from temperature-induced changes took place when excessive temperatures were removed. It should also be pointed out that all oscillators were able to restart at the extreme test temperatures and to withstand the limited thermal cycling without undergoing any significant changes in their characteristics. In addition, no physical damage was observed in the packaging material of any of these silicon oscillators due to extreme temperature exposure and thermal cycling. It is recommended that additional and more comprehensive testing under long term cycling be carried out to fully establish the reliability of these devices and to determine their suitability for use in space exploration missions under extreme temperature conditions.

Patterson, Richard L.↗

NASA NTRS1994

Measurement of a surface heat flux and temperature

The Heat Flux Microsensor is a new sensor which was recently patented by Virginia Tech and is just starting to be marketed by Vatell Corp. The sensor is made using the thin-film microfabrication techniques directly on the material that is to be measured. It consists of several thin-film layers forming a differential thermopile across a thermal resistance layer. The measured heat flux q is proportional to the temperature difference across the resistance layer q= k(sub g)/delta(sub g) x (t(sub 1) - T(sub 2)), where k(sub g) is the thermal conductivity and delta (sub g) is the thickness of the thermal resistance layer. Because the gages are sputter coated directly onto the surface, their total thickness is less than 2 micrometers, which is two orders of magnitude thinner than previous gages. The resulting temperature difference across the thermal resistance layer (delta is less than 1 micrometer) is very small even at high heat fluxes. To generate a measurable signal many thermocouple pairs are put in series to form a differential thermopile. The combination of series thermocouple junctions and thin-film design creates a gage with very attractive characteristics. It is not only physically non-intrusive to the flow, but also causes minimal disruption of the surface temperature. Because it is so thin, the response time is less than 20 microsec. Consequently, the frequency response is flat from 0 to over 50 kHz. Moreover, the signal of the Heat Flux Microsensor is directly proportional to the heat flux. Therefore, it can easily be used in both steady and transient flows, and it measures both the steady and unsteady components of the surface heat flux. A version of the Heat Flux Microsensor has been developed to meet the harsh demands of combustion environments. These gages use platinum and platinum-10 percent rhodium as the thermoelectric materials. The thermal resistance layer is silicon monoxide and a protective coating of Al2O3 is deposited on top of the sensor. The superimposed thin-film pattern of all six layers is presented. The large pads are for connection with pins used to bring the signal out the back of the ceramic. In addition to the heat flux measurement, the surface temperature is measured with a platinum resistance layer (RTS). The resistance of this layer increases with increasing temperature. Therefore, these gages simultaneously measure the surface temperature and heat flux. The demonstrated applications include rocket nozzles, SCRAM jet engines, gas turbine engines, boiling heat transfer, flame experiments, basic fluid heat transfer, hypersonic flight, and shock tube testing. The laboratory involves using one of these sensors in a small combustion flame. The sensor is made on a 2.5 cm diameter piece of aluminum nitride ceramic.

Davis, R. M.↗

NASA NTRS2018

Solar Energetic Particle Events Observed by the PAMELA Mission

Despite the significant progress achieved in recent years, the physical mechanisms underlying the origin of solar energetic particles (SEPs) are still a matter of debate. The complex nature of both particle acceleration and transport poses challenges to developing a universal picture of SEP events that encompasses both the low-energy (from tens of keV to a few hundreds of MeV)observations made by space-based instruments and the GeV particles detected by the worldwide network of neutron monitors in ground-level enhancements (GLEs). The high-precision data collected by the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) satellite experiment offer a unique opportunity to study the SEP fluxes between ∼80 MeV and a few GeV, significantly improving the characterization of the most energetic events. In particular, PAMELA can measure for the first time with good accuracy the spectral features at moderate and high energies, providing important constraints for current SEP models. In addition, the PAMELA observations allow the relationship between low and high-energy particles to be investigated, enabling a clearer view of the SEP origin. No qualitative distinction between the spectral shapes of GLE, sub-GLE and non-GLE events is observed, suggesting that GLEs are not a separate class, but are the subset of a continuous distribution of SEP events that are more intense at high energies. While the spectral forms found are to be consistent with diffusive shock acceleration theory, which predicts spectral rollovers at high energies that are attributed to particles escaping the shock region during acceleration, further work is required to explore the relative influences of acceleration and transport processes on SEP spectra.

Bruno, A.↗

NASA NTRS2014

NEQAIRv14.0 Release Notes: Nonequilibrium and Equilibrium Radiative Transport Spectra Program

NEQAIR v14.0 is the first parallelized version of NEQAIR. Starting from the last version of the code that went through the internal software release process at NASA Ames (NEQAIR 2008), there have been significant updates to the physics in the code and the computational efficiency. NEQAIR v14.0 supersedes NEQAIR v13.2, v13.1 and the suite of NEQAIR2009 versions. These updates have predominantly been performed by Brett Cruden and Aaron Brandis from ERC Inc at NASA Ames Research Center in 2013 and 2014. A new naming convention is being adopted with this current release. The current and future versions of the code will be named NEQAIR vY.X. The Y will refer to a major release increment. Minor revisions and update releases will involve incrementing X. This is to keep NEQAIR more in line with common software release practices. NEQAIR v14.0 is a standalone software tool for line-by-line spectral computation of radiative intensities and/or radiative heat flux, with one-dimensional transport of radiation. In order to accomplish this, NEQAIR v14.0, as in previous versions, requires the specification of distances (in cm), temperatures (in K) and number densities (in parts/cc) of constituent species along lines of sight. Therefore, it is assumed that flow quantities have been extracted from flow fields computed using other tools, such as CFD codes like DPLR or LAURA, and that lines of sight have been constructed and written out in the format required by NEQAIR v14.0. There are two principal modes for running NEQAIR v14.0. In the first mode NEQAIR v14.0 is used as a tool for creating synthetic spectra of any desired resolution (including convolution with a specified instrument/slit function). The first mode is typically exercised in simulating/interpreting spectroscopic measurements of different sources (e.g. shock tube data, plasma torches, etc.). In the second mode, NEQAIR v14.0 is used as a radiative heat flux prediction tool for flight projects. Correspondingly, NEQAIR has also been used to simulate the radiance measured on previous flight missions. This report summarizes the database updates, corrections that have been made to the code, changes to input files, parallelization, the current usage recommendations, including test cases, and an indication of the performance enhancements achieved.

Radiation Solver↗

NASA NTRSDate not supplied

Application of an Affine Nonlinear Galerkin Reduced-order Model to Compressible Fluid Flows

Galerkin reduced-order models (ROMs) often struggle to accurately capture multiscale fluid physics in challenging flow regimes such as flows experiencing compressibility effects. This in part stems from the global nature of both the basis construction problem and the spectral formulation itself. In this work, a multi-basis ROM is developed in an affine space based on proper orthogonal decomposition (POD) by projecting the full Navier-Stokes equations expressed in terms of the specific volume, velocity, and pressure primitive variables. The model is applied to high-fidelity numerical simulation datasets obtained for a canonical compressible flow configuration: the flow over a backward facing step at different subsonic Mach numbers. It is observed that application of an eigenvalue reassignment (ER) stabilization method is required to avoid early divergence of the ROM predictions for this configuration in the three Mach numbers tested. The sensitivity of the POD-ROM results to the choice of parameters in the stabilization algorithm is discussed.

reduced-order model↗

NASA NTRSDate not supplied

Application of an Affine Nonlinear Galerkin Reduced-order Model to Compressible Fluid Flows

reduced-order model↗

NASA NTRS1995

A molecular line study of NGC 1333/IRAS 4

Molecular line surveys and fully sampled spectral line maps at 1.3 and 0.87 mm are used to examine the physical and chemical characteristics of the extreme Class I sources IRAS 4A and 4B in the L1450/NGC 1333 molecular cloud complex. A very well collimated, jetlike molecular flow emanates from IRAS 4A, with a dynamical age of a few thousand years. Symmetric, clumpy structure along the outflow lobes suggests that there is considerable variability in the mass-loss rate or wind velocity even at this young age. Molecular emission lines toward IRAS 4A and 4B are observed to be weak in the velocity range corresponding to quiescent material surrounding the young stellar objects (YSOs). Depletion factors of 10-20 are observed for all molecules, including CO, even for very conservative mass estimates from the measured millimeter and submillimeter dust continuum. However, abundance scaled with respect to CO are similar to other dark molecular cloud cores. Such depletions could be mimicked by high dust and optical depths or increased grain emissivities at the observing frequencies of 230 and 435 GHz, but the millimeter and submillimeter spectral energy distributions suggest that this is unlikely over the single-dish size scales of 5000-10,000 AU. Dense, outflowing gas is found to be kinematically, but not spatially, distinct from the quiescent material on these size scales. If CO is used as a chemical standard for the high-velocity gas, we find substantial enhancements in abundances of several molecules in outflowing material, most notably CS, SIO, and CH3OH. The SiO emission is kinematically well displaced from the bulk cloud velocity and likely arises from directly shocked material. As is the case for CO, however, the outflow features from more volatile species are centered near the cloud velocity and are often characterized by quite low rotational temperatures. We suggest that grain-grain collisions induced by velocity shear zones surrounding the outflow axes transiently desorb the grain mantles, resulting in large abundance enhancements of selected species. Similar results have recently been obtained in several other low-mass YSOs, where outflowing gas is often both kinematically and spatailly distinct, and are illustrative of the ability of accretion and outflow processes to simultaneously modify the composition of the gas and dust surrounding young stars.

Blake, Geoffrey A.↗

NASA NTRS1999

Cosmic Ray Telescope Experiment (CRT), Pioneer 10/11 Program

In May 1996 the electrical power on Pioneer 10 was no longer adequate to support the Cosmic Ray Telescope Experiment (CRT). In March 1997 ground based operations for the mission were terminated. The 25 years of Pioneer 10 were a remarkable voyage of discovery as it ventured into a vast new unexplored region of space. These observations by the Pioneer experiments led in the development of a new discipline of space science - heliospheric physics. For cosmic ray studies it was an incredible era, leading to the identification of new energetic particle populations and processes and initiating the study of the dynamics and large-scale structure and began the study of the dynamics and large-scale structure of the outer heliosphere. A summary of some of the principal scientific findings of the CRT experiment over this period is given in the next section. Even with the cessation of data from the Pioneer 10/11 mission there remained a great deal of scientific analysis and data archiving that required on the order of an additional three years of effort on the part of the CRT team. We have tried to select those tasks where the P10/11 CRT data have an important role to play. Special emphasis is placed on long-term synoptic studies that make use of the extended temporal and spatial coverage of the missions and the full capabilities of the CRT experiment. The major scientific objectives of this study are: (1) A phenomenological study of the modulation process in the heliosphere, thereby laying the foundation for developing a unified model of cosmic ray modulation over a complete 22 year heliomagnetic cycle; (2) These modulation studies and the temporal changes in the cosmic ray intensities also provide information on the large-scale structure and dynamics of the outer heliosphere; (3) use the galactic and anomalous cosmic ray data from Pioneer 10 to obtain a more accurate estimate of the distance to the modulation boundary and to the termination shock; (4) Use the results from (2) and from temporal variations to determine whether significant modulation occurs in the region of the heliosphere; (5)study the acceleration and transport of low energy solar interplanetary energetic particles and their relation to solar activity and interplanetary phenomena.

McDonald, Frank B.↗