Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 595 records · Page 33

Adaptive computational methods for aerothermal heating analysis

The development of adaptive gridding techniques for finite-element analysis of fluid dynamics equations is described. The developmental work was done with the Euler equations with concentration on shock and inviscid flow field capturing. Ultimately this methodology is to be applied to a viscous analysis for the purpose of predicting accurate aerothermal loads on complex shapes subjected to high speed flow environments. The development of local error estimate strategies as a basis for refinement strategies is discussed, as well as the refinement strategies themselves. The application of the strategies to triangular elements and a finite-element flux-corrected-transport numerical scheme are presented. The implementation of these strategies in the GIM/PAGE code for 2-D and 3-D applications is documented and demonstrated.

Price, John M.↗

CFD application to supersonic/hypersonic inlet airframe integration

Supersonic external compression inlets are introduced, and the computational fluid dynamics (CFD) codes and tests needed to study flow associated with these inlets are outlined. Normal shock wave turbulent boundary layer interaction is discussed. Boundary layer control is considered. Glancing sidewall shock interaction is treated. The CFD validation of hypersonic inlet configurations is explained. Scramjet inlet modules are shown.

Benson, Thomas J.↗

Aerodynamic tailoring of the Learjet Model 60 wing

The wing of the Learjet Model 60 was tailored for improved aerodynamic characteristics using the TRANAIR transonic full-potential computational fluid dynamics (CFD) code. A root leading edge glove and wing tip fairing were shaped to reduce shock strength, improve cruise drag and extend the buffet limit. The aerodynamic design was validated by wind tunnel test and flight test data.

Chandrasekharan, Reuben M.↗

Overview of Results from the Cassini Magnetospheric Imaging Instrument (MIMI) During the First Year of Operations

The MIMI investigation comprises three sensors covering the indicated energy ranges: the Ion and Neutral Camera (INCA) -- 7 keV/nuc 3 Mev) and protons (1.6 < E < 160 Mev) from the back end of the dual field-of-view telescope. The Saturn observation sequences began in January, 2004 and culminated in Saturn Orbit Insertion on July 1, 2004. The MIMI sensors observed substantial activity in interplanetary space for several months prior to SOI, including several interplanetary shocks associated with corotating interaction regions, numerous increases most likely originating from particle streams in the vicinity of the Saturnian bow shock and energetic neutral atoms (ENA) emanating from Saturn s magnetosphere. Results following SOI revealed: a dynamical magnetosphere with a day-night asymmetry and an 11-hour periodicity; several water-product ions (O+, OH+, H2O+), but little N+; inferred quantities of neutral gas sufficient to cause major losses in the trapped ions and electrons in the middle and inner magnetosphere; a Titan exosphere that is a copious source of ENA; INCA imaging through ENA has also revealed a previously unknown radiation belt residing inward of the D-ring that is most likely the result of double charge-exchange between the main radiation belt and the upper layers of Saturn s exosphere. Finally, there is ample evidence for the presence of substorm-like injections of plasma that subsequently corotates for a number of days before dissipating on the night-side magnetotail. The observations will be presented and discussed in the context of current theoretical models.

Krimigis, S. M.↗

Application of Micro-ramp Flow Control Devices to an Oblique Shock Interaction

Tests are planned in the 15cm x 15cm supersonic wind tunnel at NASA Glenn to demonstrate the applicability of micro-ramp flow control to the management of shock wave boundary layer interactions. These tests will be used as a database for computational fluid dynamics (CFD) validation and Design of Experiments (DoE) design information. Micro-ramps show potential for mechanically simple and fail-safe boundary layer control.

Hirt, Stefanie↗

When a Standard Candle Flickers

The Crab is the only bright steady source in the X-ray sky. The Crab consists of a pulsar wind nebula, a synchrotron nebula, and a cloud of expanding ejecta. On small scales, the Crab is extremely complex and turbulent. X-ray astronomers have often used the Crab as a standard candle to calibrate instruments, assuming its spectrum and overall flux remains constant over time. Four instruments (Fermi/GBM, RXTE/PCA, Swift/BAT, INTEGRAL/ISGRI) show a approx.5% (50 m Crab) decline in the Crab from 2008-2010. This decline appears to be larger with increasing energy and is not present in the pulsed flux, implying changes in the shock acceleration, electron population or magnetic field in the nebula. The Crab is known to be dynamic on small scales, so it is not too surprising that its total flux varies as well. Caution should be taken when using the Crab for in-orbit calibrations.

Wilson-Hodge, Colleen A.↗

Type II Radio Bursts and Their Association With Coronal Mass Ejections in Solar Cycles 23 and 24

Context. Meter-wavelength type II solar radio bursts are thought to be the signatures of shock-accelerated electrons in the corona. Studying these bursts can give information about the initial kinematics, dynamics, and energetics of coronal mass ejections (CMEs) in the absence of white-light observations. Aims. We investigate the occurrence of type II bursts in solar cycles 23 and 24 and their association with CMEs. We also explore whether type II bursts might occur in the absence of a CME. Methods. We performed a statistical analysis of type II bursts that occurred between 200 and 25 MHz in solar cycles 23 and 24 and determined the temporal association of these radio bursts with CMEs. We categorized the CMEs based on their linear speed and angular width and studied the distribution of type II bursts with fast (≥500 km s−1 ), slow (<500 km s−1), wide (≥60◦), and narrow (<60◦) CMEs. We explored the dependence of type II bursts occurrence on the phases of the solar cycle. Results. Our analysis shows that during solar cycles 23 and 24, 768 and 435 type II bursts occurred, respectively. Of these, 79% were associated with CMEs in solar cycle 23, and 95% were associated with CMEs in solar cycle 24. However, only 4% and 3% of the total number of CMEs were accompanied by type II bursts in solar cycle 23 and 24, respectively. Most of the type II bursts in both cycles were related to fast and wide CMEs (48%). We also determined the typical drift rate and duration for type II bursts, which is 0.06 MHz s−1 and 9 min. Our results suggest that type II bursts dominate at heights ≈1.7−2.3±0.3 R . A clear majority have an onset height around 1.7 ± 0.3 R assuming the four-fold Newkirk model. Conclusions. The results indicate that most of the type II bursts had a white-light CME counterpart, but a few type II bursts lacked a clear CME association. There were more CMEs in cycle 24 than in cycle 23. However, cycle 24 contained fewer type II radio bursts than cycle 23. The onset heights of type II bursts and their association with wide CMEs reported in this study indicate that the early lateral expansion of CMEs may play a key role in the generation of these radio bursts.

Sun: activity↗

Shock capturing schemes for multidimensional flow

Progress made in the development of a genuinely multidimensional finite volume algorithm for problems in inviscid gas dynamics is presented. The approach entails: (1) the reconstruction of flowfield data using a planar wave pattern in which the strengths and orientations of the component waves are derived independently of the mesh geometry; and (2) the development of a flux formula which provides a numerical approximation to the flux at a finite volume cell face during the passage of waves which are in general oblique to the face. Several algorithms are outlined, and the most recent developments are included. The results of several numerical cases are also included.

Parpia, Ijaz H.↗

The Three-Dimensional Bow Shock of Mars as Observed by MAVEN

The Martian magnetosphere is a product of the interaction of Mars with the interplanetary magnetic field and the supersonic solar wind. The location of the bow shock has been previously modeled as conic sections using data from spacecraft such as Phobos 2, Mars Global Surveyor, and Mars Express. The Mars Atmosphere and Volatile EvolutioN (MAVEN) mission spacecraft arrived in orbit about Mars in November 2014 resulting in thousands of crossings to date. We identify over 1,000 bow shock crossings. We model the bow shock as a three-dimensional surface accommodating asymmetry caused by crustal magnetic fields. By separating MAVEN's bow shock encounters based on solar condition, we also investigate the variability of the surface. We find that the shock surface varies in shape and location in response to changes in the solar radiation, the solar wind Mach number, dynamic pressure of the solar wind, and the relative local time location of the strong crustal magnetic fields (i.e., whether they are on the dayside or on the nightside).

Supersonic↗

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

This paper presents the results of a computational study on porous-plug nozzles for their eventual application in reducing takeoff noise levels for supersonic civil transport. The study explores a wide range of porous-plug geometries, encompassing variations in hole size, plug length and porosity. Computational fluid dynamics (CFD) analyses are employed in concert with an advanced automated mesh refinement (AMR) scheme, Sketch-to-Solution (S2S), to investigate porous plug nozzle design strategies that may diminish shock strength, broad-band shock noise (BBSN), flow separation, and other undesirable characteristics associated with high takeoff noise levels. A selected subset of the porous-plug nozzles that were analyzed computationally were manufactured at subscale using stereolithography 3D printing. These nozzles were tested experimentally, obtaining far-field noise spectral data and schlieren flow visualization pictures. Detailed comparisons between CFD predictions and experimental results were conducted. These comparisons focus on gaining a better understanding of the underlying physical mechanisms responsible for noise reduction with the porous plugs relative to corresponding solid plugs.

computational fluid dynamics↗

Second coefficient of viscosity in air

Acoustic attenuation measurements in air were analyzed in order to estimate the second coefficient of viscosity. Data over a temperature range of 11 C to 50 C and at relative humidities between 6 percent and 91 percent were used. This analysis showed that the second coefficient of viscosity varied between 1900 and 20,000 times larger than the dynamic or first coefficient of viscosity over the temperature and humidity range of the data. In addition, the data showed that the molecular relaxation effects, which are responsible for the magnitude of the second coefficient of viscosity, place severe limits on the use of time-independent, thermodynamic equations of state. Compressible flows containing large streamwise velocity gradients, like shock waves, which cause significant changes in particle properties to occur during time intervals shorter than hundredths of seconds, must be modeled using dynamic equations of state. The dynamic model approach is described briefly.

Ash, Robert L.↗

Temporal and Spatial Dynamics of Nitric Oxide Production at High Latitudes During an ICME-driven Storm on Dec. 14, 2006

During geomagnetic storms driven by interplanetary coronal mass ejections (ICMEs), Earth’s upper atmosphere neutral mass density is directly driven by the geomagnetic storm cycle from intense driving and subsequent recovery. The neutral density heating and cooling times are shorter for stronger storms and longer for weaker storms. The influx of energy (Joule Heating and precipitation) into Earth’s upper atmosphere during the storm main phase allows for the enhanced production of nitric oxide (NO) at high latitudes, which in turn cools the thermosphere by radiating away excess energy. As a result, greater NO production results in quicker thermospheric cooling. Yet while the production of NO on a global scale has been examined in detail, the space and time evolution of NO is not well-known. In this study, we investigate the effects of a particular ICME-driven storm, which occurred on Dec. 14th, 2006, on the production of NO at high latitudes. We compare NO measurements from the Thermosphere, Ionosphere, Mesosphere Dynamics (TIMED) spacecraft to neutral mass density measurements from the Challenging Minisatellite Payload (CHAMP) spacecraft and find that the impact of the shock prior to the storm, not just the onset of the storm itself, is responsible for an increase in NO production. We also observe a dawn-dusk asymmetry in high-latitude NO production and discuss how the geometry of the solar wind driver may play a role in creating this asymmetry. Finally, we discuss how the improved hemispheric and local time coverage provided by the GDC and DYNAMIC missions will enable more robust studies of the spatiotemporal evolution of NO during storms. Since NO is an important factor in regulating the neutral density of the thermosphere, such studies will be relevant for more accurate orbital drag calculations.

Kevin Delano↗

Electron Dynamics in a Subproton-Gyroscale Magnetic Hole

Magnetic holes are ubiquitous in space plasmas, occurring in the solar wind, downstream of planetary bow shocks, and inside the magnetosphere. Recently, kinetic-scale magnetic holes have been observed near Earth's central plasma sheet. The Fast Plasma Investigation on NASA's Magnetospheric Multiscale (MMS) mission enables measurement of both ions and electrons with 2 orders of magnitude increased temporal resolution over previous magnetospheric instruments. Here we present data from MMS taken in Earth's nightside plasma sheet and use high-resolution particle and magnetometer data to characterize the structure of a subproton-scale magnetic hole. Electrons with gyroradii above the thermal gyroradius but below the current layer thickness carry a current sufficient to account for a 10-20 depression in magnetic field magnitude. These observations suggest that the size and magnetic depth of kinetic-scale magnetic holes is strongly dependent on the background plasma conditions.

Gershman, Daniel J.↗

The bow shock structure of IRS 7 - Wind-wind collision near the Galactic center

New structural details of ionized gas associated with IRS 7, a mass-losing supergiant near the Galactic center, are reported, and radio continuum and theoretical maps are compared. High-resolution multiconfiguration VLA observations of IRS 7 reveal ionized gas with a bow shock structure at a wavelength of 2 cm. The apex of the bow shock is facing more toward the cluster of blue stellar objects, known as IRS 16, than Sgr A*, the compact nonthermal source near the Galaxy's dynamical center. It is inferred that the shape of the mass-losing envelope of IRS 7 is influenced by Galactic winds from IRS 16. It is concluded from a comparison of the morphology of the ionized envelope of IRS 7 with simple theoretical modeling of the bow shock structure that the shape of the shock results from the collision between the Galactic center and stellar winds rather than from the motion of IRS 7 through the interstellar medium.

Yusef-Zadeh, Farhad↗

A survey of grid generation techniques in computational fluid dynamics

The use of numerically-generated body-conforming coordinate systems has become well-established in computational fluid dynamics, allowing general codes to be written which can treat boundaries of arbitrary shape. There is now a rapidly increasing interest in the dynamic coupling of the grid with the physical solution so that the grid lines continually move to achieve concentration in regions of strong variation of the physical solution and alignment with shocks, flame fronts, etc. Such adaptive coordinate systems have been developed using a variety of techniques for sensing and reacting to the developing regions of concentration and/or alignment. This paper reviews the techniques and applications of these dynamically adaptive systems and other systems.

Thompson, J. F.↗

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗

Micrometeorite dynamic pyrometamorphism: Nonstoichiometric clinoenstatite (CLEN)

Polymorphs of enstatite are common phases in many meteorites. They contain clues on their formation and the thermal evolution of their host rock which includes shock metamorphism. Rare, micron-sized, CLEN whiskers and thin platelets in chondritic porous micrometeorites were interpreted as solar nebula condensates that remained unaffected during atmospheric entry flash-heating. This CLEN formed by protoenstatite (PEN) inversion whereby the surface energy of the micron-sized PEN crystals aided the OREN-CLEN transformation or by metastable growth. Ca-poor, Mg,Fe-pyroxene with unequilibrated, intraparticle, Mg/(Mg+Fe) distributions occur in most chondritic micrometeorites. These distributions are a parent body signature that survived dynamic pyrometamorphism because the duration of the thermal spike during atmospheric entry is too short but this conclusion does not consider the ultrafine grain size of micrometeorites. The maximum temperature and duration of the heating event will depend on the kinetic energy and entry angle of the incoming micrometeorite. But lacking detailed petrological data for an individual particle, its thermal profile during atmospheric entry can not be deduced from its mass alone as a function of entry angle. In order to constrain dynamic pyrometamorphism in unmelted micrometeorites, I determined the petrological composition and silicate mineralogy in non-chondritic micrometeorites L2005T13, L2005E40, and L2006A28.

Rietmeijer, Frans J. M.↗