Engineering Papers⌕ Search

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.

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

At least 325 records · Page 18

On Greenstadt's binary index criterion.

The Greenstadt-type criterion for distinguishing between laminar and 'pulsating' (i.e., irregular) bow shock crossings can be expressed as a critical angle between the shock normal and the magnetic field. When the field is more tangent to the shock plane than this critical angle, the shock is laminar. More importantly, the requirement of Galilean invariance of the physics underlying such a criterion reveals that the Greenstadt p index must vary inversely as the solar wind velocity.

Northrop, T. G.↗

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗

Physics of magnetospheric boundary layers

This final report was concerned with the ideas that: (1) magnetospheric boundary layers link disparate regions of the magnetosphere-solar wind system together; and (2) global behavior of the magnetosphere can be understood only by understanding its internal linking mechanisms and those with the solar wind. The research project involved simultaneous research on the global-, meso-, and micro-scale physics of the magnetosphere and its boundary layers, which included the bow shock, the magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical, and simulation projects were performed on these subjects, as well as comparisons of theoretical results with observational data. Other related activity included in the research included: (1) prediction of geomagnetic activity; (2) global MHD (magnetohydrodynamic) simulations; (3) Alfven resonance heating; and (4) Critical Ionization Velocity (CIV) effect. In the appendixes are list of personnel involved, list of papers published; and reprints or photocopies of papers produced for this report.

Cairns, Iver H.↗

Experimental and Computational Investigation of the Tip Clearance Flow in a Transonic Axial Compressor Rotor

Experimental and computational techniques are used to investigate tip clearance flows in a transonic axial compressor rotor at design and part speed conditions. Laser anemometer data acquired in the endwall region are presented for operating conditions near peak efficiency and near stall at 100% design speed and at near peak efficiency at 60% design speed. The role of the passage shock/leakage vortex interaction in generating endwall blockage is discussed. As a result of the shock/vortex interaction at design speed, the radial influence of the tip clearance flow extends to 20 times the physical tip clearance height. At part speed, in the absence of the shock, the radial extent is only 5 times the tip clearance height. Both measurements and analysis indicate that under part-speed operating conditions a second vortex, which does not originate from the tip leakage flow, forms in the endwall region within the blade passage and exits the passage near midpitch. Mixing of the leakage vortex with primary flow downstream of the rotor at both design and part speed conditions is also discussed.

Suder, Kenneth L.↗

Pair cascades in extragalactic jets. 1: Gamma rays

A model of the approximately 0.1-10 GeV gamma-ray jets observed by the EGRET instrument on the Compton Gamma Ray Observatory (CGRO) is developed. It is shown that the soft X-ray background in an active galactic nuclei (AGN) contributes an opacity to pair production and that a gamma-ray photosphere or 'gamma-sphere' can be defined whose radius increases with gamma-ray energy E(sub gamma). It is proposed that the observed gamma-ray emission is due to inverse Compton scattering of the ambient soft X-rays by relativistic pairs accelerated in situ by shock fronts in a relativistic jet. For a wide range of assumed physical conditions, the emission at a given E(sub gamma) originates from near the associated gamma-spheres; emission from below the gamma-sphere initiates a cascade down to the energy where the gamma-rays can escape freely. In this model, the slope of the emergent gamma-ray spectrum is determined by the scattered, soft X-ray spectrum and the variation of the particle acceleration rate with jet radius. In general it is expected that the variation in the gamma-ray flux will be either slower or later at higher energy. It is also shown that the efficiency of conversion of energy from injected high-energy pairs to 0.1-10 GeV gamma-rays is typically high so that the models are radiatively efficient. It is argued that the observed gamma-ray jets are likely to be particle-dominated, though magnetically confined. The gamma-ray spectrum should continue down to an energy approximately 5 MeV emitted from an annihilation radius within which the pair content of the jet is limited by annihilation. This is probably the site of the beamed hard X-ray emission. It is speculated that the relativistic jets associated with radio-loud AGNs are powered electromagnetically by a spinning black hole and that they are collimated by an encircling MHD wind leaving the accretion disk at a slower speed. Powerful FR2 radio sources are formed when the hole spins rapidly and the relativistic core accelerates the MHD sheath; low-power FR1 sources ensue when the opposite occurs. Finally, it is suggested that the key factor which determines whether or not a given active nucleus can form a jet and a radio to gamma-ray nonthermal continuum is the central density of mass-losing stars which, when large, precludes the formation of a super-Alfvenic, collimating wind.

Blandford, R. D.↗

Automated Fluid Feature Extraction from Transient Simulations

In the past, feature extraction and identification were interesting concepts, but not required to understand the underlying physics of a steady flow field. This is because the results of the more traditional tools like iso-surfaces, cuts and streamlines were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of much interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one "snap-shot" of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for co-processing environments like pV3). And methods must be developed to abstract the feature and display it in a manner that physically makes sense. The following is a list of the important physical phenomena found in transient (and steady-state) fluid flow: (1) Shocks, (2) Vortex cores, (3) Regions of recirculation, (4) Boundary layers, (5) Wakes. Three papers and an initial specification for the (The Fluid eXtraction tool kit) FX Programmer's guide were included. The papers, submitted to the AIAA Computational Fluid Dynamics Conference, are entitled : (1) Using Residence Time for the Extraction of Recirculation Regions, (2) Shock Detection from Computational Fluid Dynamics results and (3) On the Velocity Gradient Tensor and Fluid Feature Extraction.

Haimes, Robert↗

Flash Spectroscopy: Emission Lines From the Ionized Circumstellar Material Around 10-Day-Old Type II Supernovae

Supernovae (SNe) embedded in dense circumstellar material (CSM) may show prominent emission lines in their early-time spectra (< or =10 days after the explosion), owing to recombination of the CSM ionized by the shock-breakout flash. From such spectra ("flash spectroscopy"), we can measure various physical properties of the CSM, as well as the mass-loss rate of the progenitor during the year prior to its explosion. Searching through the Palomar Transient Factory (PTF and iPTF) SN spectroscopy databases from 2009 through 2014, we found 12 SNe II showing flash-ionized (FI) signatures in their first spectra. All are younger than 10 days. These events constitute 14% of all 84 SNe in our sample having a spectrum within 10 days from explosion, and 18% of SNe II observed at ages <5 days, thereby setting lower limits on the fraction of FI events. We classified as "blue/featureless" (BF) those events having a first spectrum that is similar to that of a blackbody, without any emission or absorption signatures. It is possible that some BF events had FI signatures at an earlier phase than observed, or that they lack dense CSM around the progenitor. Within 2 days after explosion, 8 out of 11 SNe in our sample are either BF events or show FI signatures. Interestingly, we found that 19 out of 21 SNe brighter than an absolute magnitude M(sub R) = −18.2 belong to the FI or BF groups, and that all FI events peaked above M(sub R) = −17.6 mag, significantly brighter than average SNe II.

supernovae: general↗

Electron Energy Partition Across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

Solar wind↗

An adaptive pseudospectral method for discontinuous problems

The accuracy of adaptively chosen, mapped polynomial approximations is studied for functions with steep gradients or discontinuities. It is shown that, for steep gradient functions, one can obtain spectral accuracy in the original coordinate system by using polynomial approximations in a transformed coordinate system with substantially fewer collocation points than are necessary using polynomial expansion directly in the original, physical, coordinate system. It is also shown that one can avoid the usual Gibbs oscillation associated with steep gradient solutions of hyperbolic pde's by approximation in suitably chosen coordinate systems. Continuous, high gradient solutions are computed with spectral accuracy (as measured in the physical coordinate system). Discontinuous solutions associated with nonlinear hyperbolic equations can be accurately computed by using an artificial viscosity chosen to smooth out the solution in the mapped, computational domain. Thus, shocks can be effectively resolved on a scale that is subgrid to the resolution available with collocation only in the physical domain. Examples with Fourier and Chebyshev collocation are given.

Augenbaum, Jeffrey M.↗

An adaptive pseudospectral method for discontinuous problems

The accuracy of adaptively chosen, mapped polynomial approximations is studied for functions with steep gradients or discontinuities. It is shown that, for steep gradient functions, one can obtain spectral accuracy in the original coordinate system by using polynomial approximations in a transformed coordinate system with substantially fewer collocation points than are necessary using polynomial expansion directly in the original, physical, coordinate system. It is also shown that one can avoid the usual Gibbs oscillation associated with steep gradient solutions of hyperbolic pde's by approximation in suitably chosen coordinate systems. Continuous, high gradient solutions are computed with spectral accuracy (as measured in the physical coordinate system). Discontinuous solutions associated with nonlinear hyperbolic equations can be accurately computed by using an artificial viscosity chosen to smooth out the solution in the mapped, computational domain. Thus, shocks can be effectively resolved on a scale that is subgrid to the resolution available with collocation only in the physical domain. Examples with Fourier and Chebyshev collocation are given.

Augenbaum, J. M.↗

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗

Physics of Boundaries and their Interactions in Space Plasmas

This final report describes a brief summary of our accomplishments during the complete contract period. Traditionally, due to computational limitations, it has been impossible to obtain a global view of the magnetosphere on ion time and spatial scales. As a result, kinetic-simulations have concentrated on the local structure of different magnetospheric discontinuities and boundaries. However, due to the emergence of low cost supercomputers, as well as by taking full advantage of latest advances in data mining and visualization technology, we were able to bypass our planned (proposed) regional simulations and proceed to large-scale 3-D and 2-D global hybrid simulations of the magnetosphere. As a result, although we are only finishing the second year of the proposed activity, much of the original scientific objectives have been surpassed and new avenues of investigation have been opened. Such simulations have led us to possible explanations of some long-standing issues in magnetospheric physics. They have also enables us to make a number of important discoveries predictions, which need to be looked for in satellite data. Examples include the finding that the bow shock can become unstable to the Kelvin-Helmholtz (KH), (2) the discovery of a mechanism for intermittent reconnection due to ion physics which may be relevant to the explanation of the recurrence rate of flux transfer events (FTEs), and (3) this finding that the current sheet in the near-Earth magnetotail region can become unstable to KH with detectable, unique ionospheric signatures. Further, we demonstrated a viable mechanism for the onset of reconnection at the magnetopause, examined the detailed structure of the boundary layer incorporating curvature effects, and provided an explanation for the large core fields observed within FTEs as well as flux ropes in the magnetotail.

Omidi, Nojan↗

Physics of Boundaries and their Interactions in Space Plasmas

This final report describes a brief summary of our accomplishments during the complete contract period. Traditionally, due to computational limitations, it has been impossible to obtain a global view of the magnetosphere on ion time and spatial scales. As a result, kinetic simulations have concentrated on the local structure of different magnetospheric discontinuities and boundaries. However, due to the emergence of low cost desktop superconductors, as well as by taking full advantage of latest advances in data mining and visualization technology, we were able to bypass our planned (proposed) regional simulations and proceed to large-scale 3-D and 2-D global hybrid simulations of the magnetosphere. As a result, although we are only finishing the second year of the proposed activity, much of the original scientific objectives have been surpassed and new avenues of investigation have been opened. Such simulations have led us to possible explanations of some long-standing issues in magnetospheric physics. They have also enabled us to make a number of important discoveries/predictions, which need to be looked for in satellite data. Examples include: (1) the finding that the bow shock can become unstable to the Kelvin-Helmholtz (KH;) (2) the discovery of a mechanism for intermittent reconnection due to ion physics which may be relevant to the explanation of the recurrence rate of flux transfer events (FTEs;) and (3) the finding that the current sheet in the near-Earth magnetotail region can become unstable to KH with detectable, unique ionospheric signatures. Further, we demonstrated a viable mechanism for the onset of reconnection at the magnetopause, examined the detailed structure of the boundary layer incorporating curvature effects, and provided an explanation for the large core fields observed within FTEs as well as flux ropes in the magnetotail.

Omidi, Nojan↗

Multiple transonic solutions and a new class of shock transitions in solar and stellar winds

The steady isothermal solar wind equations are shown to admit, under certain circumstances, mutliple transonic solutions when, for example, momentum deposition gives rise to multiplee critical points in the flow. These multiple solutions consist of a continuous solution and solutions which involve shock transitions between critical solutions. The ambiguity arising from the multiplicity of the solutions can be resolved by following the time evolution of a wind profile with one critical point. Results of the numerical integration of the time-dependent equations with momentum addition show that each of these multiple solutions is physically accessible and depends on the rate of change of momentum deposition. These results suggest that standing shocks are likely to be present in the inner solar wind flow.

Habbal, S. R.↗

CFD Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature

An axisymmetric full Navier-Stokes computational fluid dynamics study is conducted to examine nozzle exhaust jet plume effects on the sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock research airplane, is considered. The computational fluid dynamics code is validated using available wind-tunnel sonic boom experimental data. The effects of grid size, spatial order of accuracy, grid type, and flow viscosity on the accuracy of the predicted sonic boom pressure signature are quantified. Grid lines parallel to the Mach wave direction are found to give the best results. Second-order accurate upwind methods are required as a minimum for accurate sonic boom simulations. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature is observed for the highly underexpanded nozzle flow. Axisymmetric computational fluid dynamics simulations show the flow physics inside the F-15 nozzle to be nonisentropic and complex. Although the one-dimensional isentropic nozzle plume results look reasonable, they fail to capture the sonic boom shock train in the highly underexpanded nozzle flow.

Bui, Trong T.↗

A unique coannular plug nozzle jet noise prediction procedure

The work presents a method for predicting jet mixing and shock noise spectra from coannular plug nozzles operated in the inverted velocity profile mode. It is shown that the jet mixing noise prediction method is based on the physics of the flow and its noise radiating characteristics, modeling the source spectrum, eddy convection, and fluid shrouding in a simplistic fashion. Further, the shock noise prediction procedure is described, noting that it is an extension of the Fisher-Harper-Bourne method for conical nozzles. Finally, it is concluded that the method has the potential of being extended to suppressors of complex geometries.

Bhutiani, P. K.↗

Aerothermodynamic heating due to shock wave/laminar boundary-layer interactions in high-enthalpy hypersonic flow

The interaction between a swept shock wave and a laminar boundary layer was investigated experimentally in high-enthalpy hypersonic flow. The effect of high-temperature, real gas physics on the interaction was examined by conducting tests in air and helium. Heat transfer measurements were made on the surface of a flat plate and a shock-generating fin using thin-film resistance sensors for fin incidence angles of 0, 5, and 10 deg at Mach numbers of 6.9 in air and 7.2 in helium. The experiments were conducted in the NASA HYPULSE expansion tube, an impulse-type facility capable of generating high-enthalpy, high-velocity flow with freestream levels of dissociated species that are particularly low. The measurements indicate that the swept shock wave creates high local heat transfer levels in the interaction region, with the highest heating found in the strongest interaction. The maximum measured heating rates in the interaction are order of magnitude greater than laminar flat plate boundary layer heating levels at the same location.

Hackett, Charles M.↗

Physical aspects of computing the flow of a viscous fluid

One of the main themes in fluid dynamics at present and in the future is going to be computational fluid dynamics with the primary focus on the determination of drag, flow separation, vortex flows, and unsteady flows. A computation of the flow of a viscous fluid requires an understanding and consideration of the physical aspects of the flow. This is done by identifying the flow regimes and the scales of fluid motion, and the sources of vorticity. Discussions of flow regimes deal with conditions of incompressibility, transitional and turbulent flows, Navier-Stokes and non-Navier-Stokes regimes, shock waves, and strain fields. Discussions of the scales of fluid motion consider transitional and turbulent flows, thin- and slender-shear layers, triple- and four-deck regions, viscous-inviscid interactions, shock waves, strain rates, and temporal scales. In addition, the significance and generation of vorticity are discussed. These physical aspects mainly guide computations of the flow of a viscous fluid.

Mehta, U. B.↗