Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock hydrodynamics”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Instability of a supersonic shock-free elliptic jet

A comparison of the measured and the computed spatial stability properties of an aspect ratio 2, supersonic, shock free, elliptic jet is presented. The shock free nature of the elliptic jet furnishes a perfect test of the validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with the linear hydrodynamic stability theory. Aerodynamic and acoustic data are measured to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. The measured mean velocity profiles are employed to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles utilized in the stability analysis. It is shown that the measured frequency associated with peak amplitude noise radiation agrees with that frequency predicted to dominate the large scale structure near the end of the potential core.

Baty, Roy S.↗
NASA NTRS2001

Jet Stability and the Generation of Superluminal and Stationary Components

We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. The increase in velocity of the trailing components is an indirect consequence of the acceleration of the expanding fluid, which is assumed to be relativistically hot; if observed, such accelerations would therefore favor an electron-positron (as opposed to proton rest mass) dominated jet.

Agudo, Ivan↗
NASA NTRS2015

A Coordinated X-Ray and Optical Campaign of the Nearest Massive Eclipsing Binary, Delta Orionis Aa. II. X-Ray Variability

We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral data set of the delta Ori Aa binary system. The four observations, obtained with Chandra ACIS (Advanced CCD Imaging Spectrometer) HETGS (High Energy Transmission Grating), have a total exposure time approximately equal to 479 kiloseconds and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range of 5-25 angstroms is confirmed, with a maximum amplitude of about plus or minus15 percent within a single approximately equal to125 kiloseconds observation. Periods of 4.76 and 2.04 days are found in the total X-ray flux, as well as an apparent overall increase in the flux level throughout the nine-day observational campaign. Using 40 kiloseconds contiguous spectra derived from the original observations, we investigate the variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S (sub XV), Si (sub XIII), and Ne (sub IX). For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phi = 0.0 when the secondary delta Ori Aa2 is at the inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.

X-ray flux↗
NASA NTRS1986

A model of inner cometary ionospheres

A hydrodynamical model to describe the movement of the thermalized charged components in the inner ionosphere of comet Halley is presented. Photoelectrons are included by applying a two-stream-type approach. The numerical scheme describes shock transitions in a natural way. Solutions are obtained for a number of different assumptions concerning electron heating rates but all show that the electron temperature increases sharply where the collisional electron neutral coupling becomes unimportant. This temperature increase is accompanied by an increase in the plasma pressure and in its associated polarization electric field, and causes the plasma flow to go subsonic. In certain cases this transition occurs as an inner shock which may explain the observed ion pile-up in Halley's comet.

Korosmezey, A.↗
NASA NTRS1995

A perturbation calculation of the dynamical effects of cosmic rays and interstellar gas on the solar wind

The modulation of galactic cosmic rays and other energetic particles by the solar wind produces a gradient in their pressure, which in turn influences the wind dynamics. The basic equations describing this interaction are presented in the 'hydrodynamic' approximation. A perturbation solution of the equations is presented for the case in which both the galactic cosmic ray pressure, and the pressure of the anomalous cosmic ray component accelerated at the wind termination shock, are small compared with the wind ram pressure. Analytical expressions for the deceleration of the wind and the modification of the shock and its location are derived in this case. These effects are estimated to have a relative magnitude of several percent in the solar wind. The interstellar neutral gas which penetrates the heliosphere is ionized, predominantly by photoionization and charge exchange with the wind, and may also have a significant dynamical effect on the wind. The basic equations describing this interaction are also presented, A perturbation solution is presented under the assumption that the mass loading and momentum-loading of the wind by the interstellar pickup ions is small compared with the wind ram pressure. Analytical expressions for the deceleration of the wind, the contribution of the pickup ions to the wind pressure, and the modification of the termination shock location are derived. Again, with the exception of pickup ion pressure which is large compared with solar wind thermal pressure, the effects are estimated to be several percent in relative magnitude in the solar wind.

Lee, M. A.↗
NASA NTRS2012

Atomic Physics of Shocked Plasma in Winds of Massive Stars

High resolution diffraction grating spectra of X-ray emission from massive stars obtained with Chandra and XMM-Newton have revolutionized our understanding of their powerful, radiation-driven winds. Emission line shapes and line ratios provide diagnostics on a number of key wind parameters. Modeling of resolved emission line velocity profiles allows us to derive independent constraints on stellar mass-loss rates, leading to downward revisions of a factor of a few from previous measurements. Line ratios in He-like ions strongly constrain the spatial distribution of Xray emitting plasma, confirming the expectations of radiation hydrodynamic simulations that X-ray emission begins moderately close to the stellar surface and extends throughout the wind. Some outstanding questions remain, including the possibility of large optical depths in resonance lines, which is hinted at by differences in line shapes of resonance and intercombination lines from the same ion. Resonance scattering leads to nontrivial radiative transfer effects, and modeling it allows us to place constraints on shock size, density, and velocity structure

Physics↗
NASA NTRS1994

The CO/SiO radiative instability in cool star atmospheres revisited

We revisit the formation of radiative instabilities in cool star atmospheres and compare our results with those given by Muchmore, Nuth, & Stencel. We have considered the combined influence of CO and SiO molecules and have computed models for a grid of effective temperatures and geometrical dilution factors for the stellar radiation. Our results are based on the analysis of the energy balance of gas elements with prescribed thermodynamic properties. Our results show that radiative instabilities are most likely primarily caused by CO, whereas SiO is expected to play only a minor role, except when the CO density is reduced compared to Local Thermodynamic Equilibrium (LTE) values or the CO band can be assumed to be optically thick. The onset of radiative instabilities is expected to be strongly modified when dynamic phenomena such as stochastic shocks are present. Our results provide strong evidence that dust formation can most likely occur via a radiative instability alone. Therefore, we present a revised version of the Muchmore et al. dust formation paradigm, which also considers hydrodynamic cooling. The new paradigm is particularly relevant in cases where dust is formed relatively close to the stellar photosphere.

Cuntz, M.↗
NASA NTRS2000

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗
NASA NTRS1992

Flux-limited diffusion in hydrodynamics

The radiation transport equation, correct to order v/c, is used to derive a flux limiter valid in high-velocity hydrodynamic flows. In the limit v = 0 the widely used result of Levermore and Pomraning (LP) (1982) is recovered. In the case of homologous flow (v varies as R) the LP result remains valid, but in the presence of shock waves and other nonhomologous flows differences may be large. We calculate the Eddington factors, which are significantly altered, and discuss the importance of having them possess the correct form in the presence of a shock wave, e.g., on the neutrino transport during the formation of a supernova.

Cooperstein, J.↗
NASA NTRS2021

Boundary Layer Instabilities Over a Cone-Cylinder-Flare Model at Mach 6

Computations are performed to investigate the boundary-layer instabilities over a sharp cone-cylinder-flare model at zero degrees angle of attack. The model geometry and the flow conditions are selected to match the experiments conducted in the Boeing/AFOSR Mach 6 Quiet Tunnel (BAM6QT) at Purdue University. The geometry consists of a nominally sharp 5-degree half-angle cone, followed by a cylindrical segment and then a 10-degree flare. An axisymmetric separation bubble is generated as a result of the laminar shock/boundary-layer interaction in the cylinder-flare region. The comparison of the laminar flow solution and the schlieren images shows a remarkable agreement between the respective locations of both the boundarylayer edge and the reattachment shock. The predicted heat flux distribution is also in agreement with the measured values downstream of the reattachment location. The analysis of convective and global instabilities is performed for flare half angles equal to 8, 10, and 12 degrees and nosetip radii equal to 0.1, 1, and 5 mm. The linear amplification of first and second Mack mode instabilities that begin to amplify in the cone region are computed with a combination of the parabolized stability equations (PSE) and the harmonic linearized Navier-Stokes equations (HLNSE). The predicted frequency spectra of the surface pressure fluctuations associated with both planar and oblique instability waves are compared with the measured spectra at the various locations of the PCB and Kulite sensors. The comparison shows that the computational analysis captures the distinct lobes within the disturbance amplification spectra measured in the experiments, but some differences in amplification characteristics are noted at low frequencies. Overall, the oblique disturbances are found to be more amplified than the planar disturbances. To our knowledge, this represents the first successful comparison between convective instability analysis and measured surface pressure fluctuations for a hypersonic configuration with a separation bubble. Finally, the global instability analysis shows that the laminar flow becomes supercritical for flare half angles larger than 8 degrees. The unstable global mode for the experimental configuration of a 10 degrees flare and a sharp nosetip cone corresponds to a stationary three-dimensional disturbance that is concentrated in the recirculation region and achieves its maximum growth rate for an azimuthal wavenumber of 5.

Boundary layer transition↗
NASA NTRS1984

Magnetic field amplification and generation in hypervelocity meteoroid impacts with application to lunar paleomagnetism

A one-dimensional numerical model for the expansion of impact-produced vapor clouds is used to investigate magnetic field generation mechanisms in events such as meteor collisions with the moon. The resulting cloud properties, such as ionization fraction, electrical conductivity, radial expansion velocity, mass density, and energy density are estimated. The model is initiated with the peak shock states and pressure thresholds for incipient and complete vaporization of anorthosite lunar surface materials by iron and GA composition meteorites. The expansion of the spherical gas cloud into a vacuum was traced with a one-dimensional explicit lagrangian hydrodynamic code. The hypervelocity impact plasmas produced are found to be significant in the amplitudes and orientations of the magnetic fields generated. An ambient magnetic field could have been provided by the core dynamo, which would have interacted with the expanding plasmas and formed induced paleomagnetic fields. Several other field-contribution mechanisms are discussed and discarded as potential remanent magnetism contributors.

Hood, L. L.↗
NASA NTRS1991

Energy Models for One-Carrier Transport in Semiconductor Devices

Moment models of carrier transport, derived from the Boltzmann equation, made possible the simulation of certain key effects through such realistic assumptions as energy dependent mobility functions. This type of global dependence permits the observation of velocity overshoot in the vicinity of device junctions, not discerned via classical drift-diffusion models, which are primarily local in nature. It was found that a critical role is played in the hydrodynamic model by the heat conduction term. When ignored, the overshoot is inappropriately damped. When the standard choice of the Wiedemann-Franz law is made for the conductivity, spurious overshoot is observed. Agreement with Monte-Carlo simulation in this regime required empirical modification of this law, or nonstandard choices. Simulations of the hydrodynamic model in one and two dimensions, as well as simulations of a newly developed energy model, the RT model, are presented. The RT model, intermediate between the hydrodynamic and drift-diffusion model, was developed to eliminate the parabolic energy band and Maxwellian distribution assumptions, and to reduce the spurious overshoot with physically consistent assumptions. The algorithms employed for both models are the essentially non-oscillatory shock capturing algorithms. Some mathematical results are presented and contrasted with the highly developed state of the drift-diffusion model.

Jerome, Joseph W.↗
NASA NTRS1982

Supernovae and the formation of galaxies

Attention is given to the cumulative effect of many SNs in a galaxy, especially a young galaxy, as well as the effect of the blast wave resulting from the combination of all individual SNs on the intergalactic medium (IGM). It is suggested that the IGM can be shock-heated by this process, so that in the late stages of an individual expanding galaxy's shock a shell may form whose mass is many times greater than that of the seed galaxy. If this shell fragments, due to gravitational instability, it will give rise to a new generation of galaxies. SN explosions may thereby provide a hydrodynamic amplifier of some significance in the process of galaxy formation.

Ostriker, J. P.↗
NASA NTRS2020

The Dynamics, Destruction, and Survival of Supernova-Formed Dust Grains

Observations have demonstrated that supernovae efficiently produce dust. This is consistent with the hypothesis that supernovae and asymptotic giant branch stars are the primary producers of dust in the universe. However, there has been a longstanding question of how much of the dust detected in the interiors of young supernova remnants can escape into the interstellar medium. We present new hydrodynamical calculations of the evolution of dust grains that were formed in dense ejecta clumps within a Cas A–like remnant. We follow the dynamics of the grains as they decouple from the gas after their clump is hit by the reverse shock. They are subsequently subject to destruction by thermal and kinetic sputtering as they traverse the remnant. Grains that are large enough (∼0.25 μm for silicates and ∼0.1 μm for carbonaceous grains) escape into the interstellar medium while smaller grains get trapped and destroyed. However, grains that reach the interstellar medium still have high velocities, and are subject to further destruction as they are slowed down. We find that for initial grain size distributions that include large (∼0.25–0.5 μm) grains, 10%–20% of silicate grains can survive, while 30–50% of carbonaceous grains survive even when the initial size distribution cuts off at smaller (0.25 μm) sizes. For a 19 Mꙩ star similar to the progenitor of Cas A, up to 0.1 Mꙩ of dust can survive if the dust grains formed are large. Thus we show that supernovae under the right conditions can be significant sources of interstellar dust.

Jonathan D. Slavin↗
NASA NTRS1993

When clusters collide - A numerical Hydro/N-body simulation of merging galaxy clusters

A 3D numerical simulation of two merging clusters of galaxies, using a hybrid Hydro/N-body code, is presented. The hydrodynamics of the code is solved by an Eulerian finite difference method. Initial results disclose that the X-ray emission of the dominant cluster becomes elongated and broadened; heating occurs at the core of the dominant cluster as a result of multiple shocks, and high velocity gas motions within the intracluster medium. It is predicted that clusters which have undergone recent mergers and do not have cooling flows will have high peculiar gas velocities and that the shocks and turbulence generated during the merger may power cluster-wide radio halos. Prolonged high-velocity gas motions through the dominant cluster core possibly play a major role in the formation and shaping of wide-angle tailed radio sources associated with central dominant galaxies. The N-body component of the simulation reveals the subcluster to be dispersed as it passes through the dominant cluster.

Roettiger, Kurt↗
NASA NTRS1993

Numerical simulation of plasma processes driven by transverse ion heating

The plasma processes driven by transverse ion heating in a diverging flux tube are investigated with numerical simulation. The heating is found to drive a host of plasma processes, in addition to the well-known phenomenon of ion conics. The downward electric field near the reverse shock generates a doublestreaming situation consisting of two upflowing ion populations with different average flow velocities. The electric field in the reverse shock region is modulated by the ion-ion instability driven by the multistreaming ions. The oscillating fields in this region have the possibility of heating electrons. These results from the simulations are compared with results from a previous study based on a hydrodynamical model. Effects of spatial resolutions provided by simulations on the evolution of the plasma are discussed.

Singh, Nagendra↗
NASA NTRS1993

Protostellar hydrodynamics: Constructing and testing a spacially and temporally second-order accurate method. 2: Cartesian coordinates

In Boss & Myhill (1992) we described the derivation and testing of a spherical coordinate-based scheme for solving the hydrodynamic equations governing the gravitational collapse of nonisothermal, nonmagnetic, inviscid, radiative, three-dimensional protostellar clouds. Here we discuss a Cartesian coordinate-based scheme based on the same set of hydrodynamic equations. As with the spherical coorrdinate-based code, the Cartesian coordinate-based scheme employs explicit Eulerian methods which are both spatially and temporally second-order accurate. We begin by describing the hydrodynamic equations in Cartesian coordinates and the numerical methods used in this particular code. Following Finn & Hawley (1989), we pay special attention to the proper implementations of high-order accuracy, finite difference methods. We evaluate the ability of the Cartesian scheme to handle shock propagation problems, and through convergence testing, we show that the code is indeed second-order accurate. To compare the Cartesian scheme discussed here with the spherical coordinate-based scheme discussed in Boss & Myhill (1992), the two codes are used to calculate the standard isothermal collapse test case described by Bodenheimer & Boss (1981). We find that with the improved codes, the intermediate bar-configuration found previously disappears, and the cloud fragments directly into a binary protostellar system. Finally, we present the results from both codes of a new test for nonisothermal protostellar collapse.

Myhill, Elizabeth A.↗
NASA NTRS1985

MHD modelling of the earth's magnetosphere

Attention is given to the shock capturing technique, nonuniform grid system, and multiple time scale problem aspects of a global MHD model of the earth's magnetosphere that uses numerical methods. Because of the fast Alfven waves near the earth, the model is more difficult to solve than problems of hydrodynamic flow past bodies. The Rusanov (1962) scheme, which is a variant of the Lax (1954) method that maintains viscosity at a minimum value except where a large one is required, is used. The magnetosphere shape obtained by these means is noted.

Wu, C. C.↗