Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrodynamic radius”

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 NTRS1994

The heating of nova ejecta by radioactive decays of the beta-unstable nuclei

Recent nucleosynthesis and hydrodynamic calculations of the consequences of accretion onto massive ONeMg white dwarf stars show that under certain circumstances significant amounts of the beta-unstable nuclei can be produced and ejected by the resulting explosion. We use these calculations as a guide in order to obtain the conditions under which the heating of the ejected material by the nonthermal electrons and positrons produced by the decays of the beta-unstable nuclei is sufficient to overcome the cooling from adiabatic expansion and lead to the production of X-ray-emitting coronal gas. These conditions are as follows: (1) a mass fraction for Na-22 of the order of 10(exp -3) or greater, (2) an expansion velocity in the range approximately 10(exp 2) - 10(exp 3) km/s, (3) a photospheric radius of approximately 10(exp 14) cm, (4) if the density distribution in the atmosphere satisfies a power law, then the exponent must be less than 3 for heating to overcome adiabatic cooling. Both the simulations of the outburst and the model atmosphere fits to the observed energy distributions, however, imply that the exponent is greater than or = 3 during the early phases of the outburst. Nevertheless, for a value of the exponent of 2, we predict the time when hot coronal gas can form during the expansion phases of the envelope.

Pistinner, Shlomi↗
DOE OSTI2025

Spatially and temporally resolved plasma parameter estimations of laser heated MagLIF relevant gas pipes at NIF

The ability to control laser pre-heat is an integral part of the inertial confinement fusion concept known as Magnetized Liner Inertial Fusion. This process is studied at the National Ignition Facility (NIF) where 4 of the 192 laser beams are propagated through a 1-cm long gas cell where they deposit >20 kJ of energy into the gaseous fuel via inverse bremsstrahlung absorption. This process ionizes the gas, producing a plasma that follows behind the laser front and expands over the radius of the cell. Emission from this plasma, as viewed by a gated x-ray detector, can be used to build spatially and temporally resolved estimations of the pre-heat plasma's density and temperature profiles. This can then be used to estimate the plasma pressure, internal energy, and radiation losses. Estimations show the evolution of the plasma in magnetized and unmagnetized gas cells filled with ambient temperature neopentane (C5H12) +1% Ar, as well as unmagnetized cryogenically cooled (32 K) deuterium +1% Ne filled targets. This analysis shows the effects of initial gas-fill density, composition, and axial magnetization on the time-dependent plasma parameters. Previously, these parameters at the NIF had not been experimentally characterized, and these estimations provided a potential new means of testing radiation magneto-hydrodynamic predictive capability models. Results in unmagnetized targets have strong agreement with simulations. However, in targets with a 19 T applied axial magnetic field, this method yields electron temperatures up to 100% hotter than those predicted by HYDRA codes.

Bremsstrahlung↗
NASA NTRS2021

VLTI-MATISSE Chromatic Aperture-Synthesis Imaging of η Carinae’s Stellar Wind Across the Brα Line: Periastron Passage Observations in February 2020

Context: Eta Carinae is a highly eccentric, massive binary system (semimajor axis∼15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, ηCar A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims: Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods: Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results: We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar wind in several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R∼960). Our observations were performed close to periastron passage in February 2020 (orbital phase∼14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (∼14 au). The radius of the faintest outer wind regions is∼26 mas (∼60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84±0.06 mas (6.54±0.14 au). We compare the derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

G Weigelt↗
DOE OSTI2025

Testing the thermal Sunyaev-Zel’dovich power spectrum of a halo model using hydrodynamical simulations

Statistical properties of large-scale cosmological structures serve as powerful tools for constraining the cosmological properties of our Universe. Tracing the gas pressure, the thermal Sunyaev-Zel’dovich (tSZ) effect is a biased probe of mass distribution and, hence, can be used to test the physics of feedback or cosmological models. Therefore, it is crucial to develop robust modelling of hot gas pressure for applications to tSZ surveys. Since gas collapses into bound structures, it is expected that most of the tSZ signal is within halos produced by cosmic accretion shocks. Hence, simple empirical halo models can be used to predict the tSZ power spectra. In this study, we employed the HMx halo model to compare the tSZ power spectra with those of several hydrodynamical simulations: the Horizon suite and the Magneticum simulation. We examine various contributions to the tSZ power spectrum across different redshifts, including the one- and two-halo term decomposition, the amount of bound gas, the importance of different masses, and the electron pressure profiles. Our comparison of the tSZ power spectrum reveals discrepancies between the halo model and cosmological simulations that increase with redshift. We find a 20% to 50% difference between the measured and predicted tSZ angular power spectrum over the multipole range ℓ = 10 3 − 10 4 . Our analysis reveals that these differences are driven by the excess of power in the predicted two-halo term at low k and in the one-halo term at high k . At higher redshifts ( z ∼ 3), simulations indicate that more power comes from outside the virial radius than from inside, suggesting a limitation in the applicability of the halo model. We also observe differences in the pressure profiles, despite the fair level of agreement on the tSZ power spectrum at low redshift with the default calibration of the halo model. In conclusion, our study suggests that the properties of the halo model need to be carefully controlled against real or mock data to be proven useful for cosmological purposes.

Ayçoberry, Emma (ORCID:0000000292351195)↗
NASA NTRS1995

Gaseous Vortices in Barred Spiral Galaxies

During the course of examining many two-dimensional, as well as a smaller sample of three-dimensional, models of gas flows in barred spiral galaxies, we have been impressed by the ubiquitous presence fo vortex pairs, oriented roughly perpendicular to their bars, with one vortex on each side. The vortices are obvious only when viewed in the bar frame, and the centers of their velocity fields usually are near Lagrangian points L(sub 4,5). In all models that we have studied, the vortices form on essentially the same time scale as that for the development of gaseous spiral arms, typically two bar rotations. Usually the corotation radius, r(sub c), lies slightly beyond the end of the bar. Depending upon the mass distributions of the various components, gas spirals either into, or out of, the vortices: In the former case, the vortices become regions of high density, whereas the opposite is true if the gas spirals out of a vortex. The models described in this paper have low-density vortices, as do most of the models we have studied. Moreover, usually the vortex centers lie approximately within +/- 15 deg of L(sub 4,5). In the stellar dynamic limit, when pressure and viscous forces are absent, short-period orbits exist, centered on L(sub 4,5). These orbits need not cross and therefore their morphology is that of gas streamlines, that is, vortices. We believe that the gas vortices in our models are hydrodynamic analogues of closed, short-period, libration orbits centered on L(sub 4,5).

England, Martin N.↗
DOE OSTI2023

Theory of an Automatic Seepage Meter and Ramifications for Applications

A new approach for measuring fluxes across surface water—groundwater interfaces was recently proposed. The Automatic Seepage Meter (ASM) is equipped with a precise water level sensor and digital memory that analyzes water level time series in a vertical tube inserted into a streambed. The ability to infer flux values with high temporal resolution relies on an accurate interpretation of water level dynamics inside the tube. Here, we reduce the three-dimensional hydrodynamic problem that describes the ASM water level in a variety of field conditions to a single ordinary differential equation. This novel general analytical solution for estimating ASM responses is more comprehensive and flexible than previous approaches and is applicable to the entire range of field conditions, including steady or transient stream stages, evaporation, rainfall, and noise. For example, our analysis determines the timing of the nonmonotonic ASM response to a monotonic linear stream stage variation and explains previously used empirical parabolic approximation for estimating fluxes. We present algorithms for simultaneous inference of vertical interface flux and hydraulic conductivity values together with an example code. We quantify how the accuracy of parameter estimation depends on test duration and noise amplitude and propose how our analysis can be used to optimize field test protocols. On this basis, changing the ASM geometry by increasing the radius and decreasing tube insertion depth may enable ASM field test protocols that estimate interface flux and hydraulic conductivity faster while maintaining desired accuracy. Potential applications of joint parameter estimation are suggested.

47 OTHER INSTRUMENTATION↗
NASA NTRS2003

Cellular Instabilities and Self-Acceleration of Expanding Spherical Flames

In the present investigation we aim to provide experimental information on and thereby understanding of the generation and propagation of spark-ignited, outwardly propagating cellular flames, with three major focuses. The first is to unambiguously demonstrate the influence of the four most important parameters in inducing hydrodynamic and diffusional-thermal cellularities, namely thermal expansion, flame thickness, non-unity Lewis number, and global activation energy. The second is to investigate the critical state for the onset of cellularity for the stretch-affected, expanding flame. The third is to identify and consequently quantify the phenomena of self-acceleration and possibly auto-turbulization of cellular flames. Due to space limitation the effects of activation energy and the critical state for the onset of cellularity will not be discussed herein. Experiments were conducted using C3H8-air and H2-O2-N2 mixtures for their opposite influences of non-equidiffusivity. The additional system parameters varied were the chamber pressure (p) and the mixture composition including the equivalence ratio (phi). From a sequence of the flame images we can assess the propensity of cell formation, and determine the instantaneous flame radius (R), the flame propagation rate, the global stretch rate experienced by the flame, the critical flame radius at which cells start to grow, and the average cell size.

Law, C. K.↗
DOE OSTI2021

Hohlraum x-ray preheat asymmetry measurement at the ICF capsule via Mo ball fluorescence imaging

In inertial confinement fusion, penetrating asymmetric hohlraum preheat radiation (>1.8 keV, which includes high temperature coronal M-band emission from laser spots) can lead to asymmetric ablation front and ablator–fuel interface hydrodynamic instability growth in the imploding capsule. First experiments to infer the preheat asymmetries at the capsule were performed on the National Ignition Facility for high density carbon (HDC) capsules in low density fill (0.3 mg/cc 4He) Au hohlraums by time resolved imaging of 2.3 keV fluorescence emission of a smaller Mo sphere placed inside the capsule. Measured Mo emission is pole hot (P2 > 0) since M-band is generated mainly by the outer laser beams as their irradiance at the hohlraum wall is 5× higher than for the inner beams. P2 has a large swing vs time, giving insight into the laser heated hohlraum dynamics. P4 asymmetry is small at the sphere due to efficient geometric smoothing of hohlraum P4 asymmetries at large hohlraum-to-capsule radii ratios. The asymmetry at the HDC capsule is inferred from the Mo emission asymmetry accounting for the Mo/HDC radius difference and HDC capsule opacity.

Dewald, E. L. (ORCID:0000000340546967)↗
NASA NTRS2013

The 0.5-2.22 micrometer Scattered Light Spectrum of the Disk around TW Hya: Detection of a Partially Filled Disk Gap at 80 AU*

We present a 0.5-2.2 micrometer scattered light spectrum of the circumstellar disk around TW Hya from a combination of spatially resolved Hubble Space Telescope STIS spectroscopy and NICMOS coronagraphic images of the disk. We investigate the morphology of the disk at distances greater than 40 AU over this wide range of wavelengths, and identify the presence of a depression in surface brightness at approximately 80 AU that could be caused by a gap in the disk. Additionally, we quantify the surface brightness, azimuthal symmetry, and spectral character of the disk as a function of radius. Our analysis shows that the scattering efficiency of the dust is largely neutral to blue over the observed wavelengths. We model the disk as a steady α-disk with an ad hoc gap structure. The thermal properties of the disk are selfconsistently calculated using a three-dimensional radiative transfer code that uses ray tracing to model the heating of the disk interior and scattered light images. We find a good fit to the data over a wide range of distances from the star if we use a model disk with a partially filled gap of 30% depth at 80 AU and with a self-similar truncation knee at 100 AU. The origin of the gap is unclear, but it could arise from a transition in the nature of the disk's dust composition or the presence of a planetary companion. Based on scalings to previous hydrodynamic simulations of gap-opening criteria for embedded proto-planets, we estimate that a planetary companion forming the gap could have a mass between 6 and 28 solar mass.

radiative transfer↗
NASA NTRS2002

Hydrodynamic Models of Line-Driven Accretion Disk Winds III: Local Ionization Equilibrium

We present time-dependent numerical hydrodynamic models of line-driven accretion disk winds in cataclysmic variable systems and calculate wind mass-loss rates and terminal velocities. The models are 2.5-dimensional, include an energy balance condition with radiative heating and cooling processes, and includes local ionization equilibrium introducing time dependence and spatial dependence on the line radiation force parameters. The radiation field is assumed to originate in an optically thick accretion disk. Wind ion populations are calculated under the assumption that local ionization equilibrium is determined by photoionization and radiative recombination, similar to a photoionized nebula. We find a steady wind flowing from the accretion disk. Radiative heating tends to maintain the temperature in the higher density wind regions near the disk surface, rather than cooling adiabatically. For a disk luminosity L (sub disk) = solar luminosity, white dwarf mass M(sub wd) = 0.6 solar mass, and white dwarf radii R(sub wd) = 0.01 solar radius, we obtain a wind mass-loss rate of M(sub wind) = 4 x 10(exp -12) solar mass yr(exp -1) and a terminal velocity of approximately 3000 km per second. These results confirm the general velocity and density structures found in our earlier constant ionization equilibrium adiabatic CV wind models. Further we establish here 2.5D numerical models that can be extended to QSO/AGN winds where the local ionization equilibrium will play a crucial role in the overall dynamics.

Pereyra, Nicolas Antonio↗
DOE OSTI2022

Shocks in the stacked Sunyaev-Zel’dovich profiles of clusters II: Measurements from SPT-SZ + Planck Compton- y map

ABSTRACT We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg2 SPT-SZ survey on the locations of clusters identified in that same data set. The sample contains 516 clusters with mean mass $\langle M_{\rm 200m}\rangle = 10^{14.9} \, {\rm M}_\odot$ and redshift 〈z〉 = 0.55. We analyse in parallel a set of zoom-in hydrodynamical simulations from the three hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = 1.08 ± 0.09, measured at 3.1σ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = 4.58 ± 1.24 at 2.0σ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyse clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius ($R_{\rm sh,\, acc}$) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, $R_{\rm sh,\, acc}/R_{\rm sp}\gt 2.16 \pm 0.59$.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

Fractionation of the Early Terrestrial Atmospheres: Dynamical Escape

Hydrodynamic escape may have played a significant role in the early fractionation of the atmospheres of the terrestrial planets. This possibility has been demonstrated in the last two decades by numerous models that show radial, transonic flow of hydrogen can occur in the presence of sufficient solar EUV Hydrodynamic escape may have played a significant role in the early fractionation of the atmospheres of the terrestrial planets. This possibility has been demonstrated in the last two decades by numerous models that show radial, transonic flow of hydrogen can occur in the presence of sufficient solar EUV flux, thought to exist in the first 500 My. The models show that the larger the solar flux the greater the mass of the fractionating species, which are accelerated to escape speeds by the hydrogen wind through drag processes. As the atmospheres evolve and the solar EUV flux wanes, the maximum mass of flowing gas constituents decreases until all gases become static. We show that fractionation can continue beyond this point when non-radial flow and dynamically enhanced Jeans escape are considered. For example, the early terrestrial atmospheres are thought to have had large hydrogen contents, resulting in exobase altitudes of a planetary radius or more. In this case, rotational speeds at the exobases of Earth and Mars would be large enough so that light constituents would "spin" off and fractionate, especially at equatorial latitudes. Also, in the presence of transonic flow of hydrogen only, non-radial expansion throws heavier gases to high altitudes in the exosphere, accompanied by strong bulk speeds at the exobase, which results in enhanced thermal escape fluxes and fractionation. flux, thought to exist in the first 500 My. The models show that the larger the solar flux the greater the mass of the fractionating species, which are accelerated to escape speeds by the hydrogen wind through drag processes. As the atmospheres evolve and the solar EUV flux wanes, the maximum mass of flowing gas constituents decreases until all gases become static. We show that fractionation can continue beyond this point when non-radial flow and dynamically enhanced Jeans escape are considered. For example, the early terrestrial atmospheres are thought to have had large hydrogen contents, resulting in exobase altitudes of a planetary radius or more. In this case, rotational speeds at the exobases of Earth and Mars would be large enough so that light constituents would "spin" off and fractionate, especially at equatorial latitudes. Also, in the presence of transonic flow of hydrogen only, non-radial expansion throws heavier gases to high altitudes in the exosphere, accompanied by strong bulk speeds at the exobase, which results in enhanced thermal escape fluxes and fractionation.

Hartle, Richard E.↗
NASA NTRS1960

NASA Research on the Hydrodynamics of the Gaseous Vortex Reactor

The experimental and analytical results to date of a study of a two-component gaseous vortex system are presented in this paper. Analytical expressions for tangential velocity and static-pressure profiles in a turbulent vortex show good agreement with experimental data. Airflow rates from 0.075 to 0.14 pound per second and corresponding tangential velocities from 160 to 440 feet per second are correlated by turbulent Reynolds numbers from 1.95 to 2.4. An analysis of an air-bromine gas mixture in a turbulent vortex indicates that a boundary value of bromine-to-air radial velocity ratio (u(2)/u(1)) of 0.999 gives essentially no bromine buildup, while a value of 0.833 results in considerable separation. For a constant value of (u(2)/u(1))(0) the bromine buildup increases as (1) the tangential velocity increases, (2) the air-to-bromine weight-flow ratio decreases, (3) the airflow rate decreases, (4) the temperature decreases, and (5) the turbulence decreases. Analytical temperature, pressure, and tangential-velocity profiles are also presented. Preliminary experimental results indicate that the flow of an air-bromine mixture through a vortex field results in a bromine density increase to a maximum value; followed by a decrease; the air density exhibits a uniform decrease from the outer vortex radius to the exhaust-nozzle radius.

Ragsdale, Robert G.↗
DOE OSTI2020

Limitations to the ‘basic’ HOD model and beyond

ABSTRACT We use the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based halo occupation distribution (HOD) approach to modelling the galaxy–halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard (basic) HOD model, we find that, on average, the ‘basic’ HOD model underpredicts the real-space correlation function in the TNG300 box by ∼15 per cent on scales of $1 \,\,\lt\,\, r \,\,\lt\,\, 20 \ {\rm Mpc}\, h^{-1}$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of ‘secondary’ parameters with which to augment the ‘basic’ HOD, we find that the local environment of the halo, the velocity dispersion anisotropy, β, and the product of the half-mass radius and the velocity dispersion, σ2Rhalfmass, are the three most effective measures of assembly bias that help reconcile the ‘basic’ HOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch, and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy–halo connection can be constructed if we combine both mass and local environment information about the halo.

42 ENGINEERING↗
NASA NTRS1994

X-ray clusters from a high-resolution hydrodynamic PPM simulation of the cold dark matter universe

A new three-dimensional hydrodynamic code based on the piecewise parabolic method (PPM) is utilized to compute the distribution of hot gas in the standard Cosmic Background Explorer (COBE)-normalized cold dark matter (CDM) universe. Utilizing periodic boundary conditions, a box with size 85 h(exp-1) Mpc, having cell size 0.31 h(exp-1) Mpc, is followed in a simulation with 270(exp 3)=10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, Sigma(sub 8)=1.05, Omega(sub b)=0.06, we find the X-ray-emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. The results, which are compared with those obtained in the preceding paper (Kang et al. 1994a), may be used in conjuction with ROSAT and other observational data sets. Overall, the results of the two computations are qualitatively very similar with regard to the trends of cluster properties, i.e., how the number density, radius, and temeprature depend on luminosity and redshift. The total luminosity from clusters is approximately a factor of 2 higher using the PPM code (as compared to the 'total variation diminishing' (TVD) code used in the previous paper) with the number of bright clusters higher by a similar factor. The primary conclusions of the prior paper, with regard to the power spectrum of the primeval density perturbations, are strengthened: the standard CDM model, normalized to the COBE microwave detection, predicts too many bright X-ray emitting clusters, by a factor probably in excess of 5. The comparison between observations and theoretical predictions for the evolution of cluster properties, luminosity functions, and size and temperature distributions should provide an important discriminator among competing scenarios for the development of structure in the universe.

Bryan, Greg L.↗
DOE OSTI2023

Proto-neutron star convection and the neutrino-driven wind: implications for the r-process

The neutrino-driven wind from proto-neutron stars is a proposed site for r-process nucleosynthesis, although most previous work has found that a wind heated only by neutrinos cannot produce the third r-process peak. However, several groups have noted that introducing a secondary heating source within the wind can change the hydrodynamic conditions sufficiently for a strong r-process to proceed. One possible secondary heating source is gravito-acoustic waves, generated by convection inside the proto-neutron star. As these waves propagate into the wind, they can both accelerate the wind and shock and deposit energy into the wind. Additionally, the acceleration of the wind by these waves can reduce the total number of neutrino captures and thereby reduce the final electron fraction of the wind. In neutron rich conditions, all of these effects can make conditions more favourable for r-process nucleosynthesis. Here, we present a systematic investigation of the impact of these convection-generated gravito-acoustic waves within the wind on potential nucleosynthesis. We find that wave effects in the wind can generate conditions favourable for a strong r-process, even when the energy flux in the waves is a factor of 10 –4 smaller than the total neutrino energy flux and the wind is marginally neutron rich. Nevertheless, this depends strongly on the radius at which the waves become non-linear and form shocks. We also find that both entropy production after shock formation and the acceleration of the wind due to stresses produced by the waves prior to shock formation impact the structure and nucleosynthesis of these winds.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Eigenmodes of Ducted Flows With Radially-Dependent Axial and Swirl Velocity Components

This report characterizes the sets of small disturbances possible in cylindrical and annular ducts with mean flow whose axial and tangential components vary arbitrarily with radius. The linearized equations of motion are presented and discussed, and then exponential forms for the axial, circumferential, and time dependencies of any unsteady disturbances are assumed. The resultant equations form a generalized eigenvalue problem, the solution of which yields the axial wavenumbers and radial mode shapes of the unsteady disturbances. Two numerical discretizations are applied to the system of equations: (1) a spectral collocation technique based on Chebyshev polynomial expansions on the Gauss-Lobatto points, and (2) second and fourth order finite differences on uniform grids. The discretized equations are solved using a standard eigensystem package employing the QR algorithm. The eigenvalues fall into two primary categories: a discrete set (analogous to the acoustic modes found in uniform mean flows) and a continuous band (analogous to convected disturbances in uniform mean flows) where the phase velocities of the disturbances correspond to the local mean flow velocities. Sample mode shapes and eigensystem distributions are presented for both sheared axial and swirling flows. The physics of swirling flows is examined with reference to hydrodynamic stability and completeness of the eigensystem expansions. The effect of assuming exponential dependence in the axial direction is discussed.

Kousen, Kenneth A.↗
NASA NTRS1994

Comet Shoemaker-Levy 9: Impact on Jupiter and plume evolution

The impact of fragments of Comet Shoemaker-Levy 9 on Jupiter and the resulting vapor plume expansion are investigated by conducting three-dimensional numerical simulations using the smoothed particle hydrodynamics (SPH) method. An icy body, representing the cometary fragments, with a velocity of 60 km/sec and a diameter of 2 km can penetrate to 350 km below the 1-bar pressure level in the atmosphere. Most of the initial kinetic energy of the fragment is transferred to the atmosphere between 50 km and 300 km below the 1-bar pressure level. The shock-heated atmospheric gas in the wake is totally dissociated and partially ionized. Scaling our SPH results to other sizes indicates that fragments larger than approximately 100 m in diameter can penetrate to below the visible cloud decks. The energy deposited in the atmosphere is explosively released in the upward expansion of the resulting plume. The plume preferentially expands upward rather than horizontally due to the density gradient of the ambient atmosphere. It rises greater than or equal to 10(exp 2) km in approximately 10(exp 2) sec. Eventually the total atmospheric mass ejected to above 1 bar is greater than or equal to 40 times the initial mass of the impactor. The plume temperature at a radius approximately 10(exp 3) km is greater than 10(exp 3) K for 10(exp 3) sec for a 2-km fragment. We predict that impact-induced plumes will be observable with the remote sensing instruments of the Galileo spacecraft. As the impact site rotates into the view of Earth some 20 min after the impact, the plume expansion will be observable using the Hubble Space Telescope (HST) and from visible and infrared instruments on groundbased telescopes. The rising plume reaches approximately 3000 km altitude in approximately 10 min and will be visible from Earth.

Takata, Toshiko↗