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 37 records · Page 2

A numerical hydrodynamic study of coalescence in head-on collisions of identical stars.

The specific results reported refer to head-on collisions between identical polytropes of index 3 having solar mass and radius. If the polytropes were initially at rest at infinity, then about 5% of the combined mass is lost by ejection following collision. The volatilized mass fraction rises to about 18% for an initial relative collision velocity of 1000 km/sec at infinite separation, and to about 60% for the 2000 km/sec case. Since the initial kinetic and gravitational energies balance for a relative velocity of 1512 km/sec at infinity, it may be seen that net coalescence persists to velocities somewhat in excess of this figure. Mass ejection takes place in two ways simultaneously: (1) by a rapid sideward expulsion of fluid in a massive lateral sheet normal to the collision axis, and (2) as a result of two recoil shocks which lead momentum flows backward along this axis. The lateral effect has similarities to the expansion of gas into a vacuum i.e., shocks are not involved. However, the ejection of material from the rear colliding hemisphere due to the recoil shocks predominates at low collision velocities. As the velocity increases, both effects strengthen, but the lateral expulsion intensifies more rapidly than the recoil shocks.

Seidl, F. G. P.↗

Hydrodynamic study of Supernova 1987A: The phase of a wave of cooling and recombination

The dependence of the bolometric light curve and of the effective temperature on the density distribution in a progenitor, its chemical composition, mass, radius, and explosion energy is studied. It is shown that, just before the supernova 1987 A outburst, the outside layers of the blue supergiant Sk-69.202 deg had a density distribution similar to that of the polytropic model with an index of n equals 3, a chemical composition with a mass fraction of hydrogen of the order of 0.1, and a relative helium abundance of about 0.9.

Utrobin, V. P.↗

Elemental Fractionation During Rapid Accretion of the Moon Triggered by a Giant Impact

Recently, Ida et al. made an N-body simulation of lunar accretion from a protolunar disk formed by a giant impact. One of their important conclusions is that the accretion time of the Moon is as short as one month. Such rapid accretion is a necessary consequence of the high surface density of a lunar mass disk accreting just beyond the Roche limit (about 3Re); the Safronov accretion time (a few days) is even shorter. The energy of accretion always exceeds the gravitational binding energy of newly arriving matter. Hence, without an energy sink, the accreting body is thermally unstable. For the Earth and other planets, radiation acts as the sink. However, in such a short accretion time, the Moon cannot radiate the accretional energy. Even radiating at a silicate cloudtop temperature of roughly 2000 K, it would take more than 100 yr to radiatively cool the Moon. The plausible alternative heat sinks are heat capacity, latent heat of vaporization, and thermal escape of the gas to space (i.e., hydrodynamic blowoff). The latter becomes plausible for the Moon because the scale height at 2000 K (about 300 km) is a significant fraction of the lunar radius. The early stages of lunar (or "lunatesimal") growth release relatively little energy and can occur simply by heating the material, especially if the accreting material is originally cold. However, the material is unlikely to be cold, because the disk itself is hot and cooling time is long, while the lunar accretion time iss very short. Therefore, the moon is likely to accrete condensed material just after it condenses. Accordingly, the newly accreted material will be on the verge of vaporization and will have very little heat capacity to spare. The immediate heat sink is the latent heat of vaporization. Most of the vapor will escape from the moon, because the thermal energy in the gas can be used to drive escape. However, vaporization is generally incomplete. the latent heat of vaporization exceeds the energy of accretion. Viewed globally, the accretional energy is about half the energy required to vaporize the entire Moon. Thus to first approximation, half of the Moon-forming material can be vaporized and lost during accretion. During this process, we would expect preferential loss of relatively volatile elements. Escape will retard the rate of accretion. To test these ideas, we computed detailed models of the thermal state of the Moon during accretion. We pay special attention to the structure of the silicate atmosphere and its loss rate by calculating the chemical species at equilibrium. We used the PHEQ program which includes 12 elements (H,O,C,Mg,Si,Fe,Ca, Al, Na,Ti, and N.) and 272 compounds (including ionic compounds). Because of the large heats of vaporization and ionization, the adiabatic atmosphere is nearly isothermal and massive escape is expected. The pressure of the atmosphere is determined by the balance between vaporization of a accreting material and escape. If the accretion time is one month, a 0.3 bar atmosphere is expected. Elemental fractionation depends strongly on the temperature of the accreting material. The initial temperature of the material can be estimated from the condition of gravitational instability in the protolunar disk. As shown by Ida et al, accretion starts when gravitational instability occurs when more than 99% of the material condenses. At this point, all of Ca, Al, Si, Mg, and Fe, and 95% of Na (probably K also), are in condensed phases. If the moon is formed from the accretion of such material, volatile elements such as Na, and K are retained by the moon only early in accretion. At later times, K and Na are lost and a fraction of the MG, Si and Fe is lost. However, refractory elements such as Ca and Al are retained and so achieve a mild degree (factor 2) of superabundance.

Abe, Y.↗

Two-dimensional viscous accretion disk models. I - On meridional circulations in radiative regions

It has been suggested that thin Keplerian disks, like rotating stars, cannot attain hydrostatic and thermal equilibrium simultaneously if the angular velocity is a function of radius only. To address the problem of vertical structure and stability of thin accretion disks, we performed 2D axisymmetric calculations by solving the set of fully nonlinear hydrodynamic equations, including radiation transport. The results indicate that the variation of the rotational velocity with height is too small to sustain hydrostatic equilibrium. Instead, the variation of the radial (inflow) velocity with height suffices to establish a stationary state. In particular, we find that for low values of the viscosity parameter alpha, there is mass outflow in the central parts of the disk close to the equatorial plane, and inflow only near the surface. Only for modestly high values of alpha, the flow throughout the disk is directed inward. Thus, the flow within the disk is the result of viscous radial inflow plus a circulatory flow directed outward in the midplane and inward near the surface.

Kley, W.↗

Curvilinear grids for sinuous river channels

In order to effectively analyze the flow in sinuous river channels, a curvilinear grid system was developed for use in the appropriate hydrodynamic code. The CENTERLINE program was designed to generate a two dimensional grid for this purpose. The Cartesian coordinates of a series of points along the boundaries of the sinuous channel represent the primary input to CENTERLINE. The program calculates the location of the river centerline, the distance downstream along the centerline, and both radius of curvature and channel width as a function of such distance downstream. These parameters form the basis for the generation of the curvilinear grid. Based on input values for longitudinal and lateral grid spacing, the corresponding grid system is generated and a file is created containing the appropriate parameters for use in the associated explicit finite difference hydrodynamic programs. Because of the option for a nonuniform grid, grid spacing can be concentrated in areas containing the largest flow gradients.

Tatom, F. B.↗

Galactic scale gas flows in colliding galaxies: 3-dimensional, N-body/hydrodynamics experiments

We present some results from three dimensional computer simulations of collisions between models of equal mass galaxies, one of which is a rotating, disk galaxy containing both gas and stars and the other is an elliptical containing stars only. We use fully self consistent models in which the halo mass is 2.5 times that of the disk. In the experiments we have varied the impact parameter between zero (head on) and 0.9R (where R is the radius of the disk), for impacts perpendicular to the disk plane. The calculations were performed on a Cray 2 computer using a combined N-body/smooth particle hydrodynamics (SPH) program. The results show the development of complicated flows and shock structures in the direction perpendicular to the plane of the disk and the propagation outwards of a density wave in both the stars and the gas. The collisional nature of the gas results in a sharper ring than obtained for the star particles, and the development of high volume densities and shocks.

Lamb, Susan A.↗

Tidal disruption of Periodic Comet Shoemaker-Levy 9 and a constraint on its mean density

The apparent tidal disruption of Periodic Comet Shoemaker-Levy 9 (1993e) during a close encounter within approximately 1.62 planetary radii of Jupiter can be used along with theoretical models of tidal disruption to place an upper bound on the density of the predisruption body. Depending on the theoretical model used, these upper bounds range from rho(sub c) less than 0.702 +/- 0.080 g/cu cm for a simple analytical model calibrated by numerical smoothed particle hydrodynamics (SPH) simulations to rho(sub c) less than 1.50 +/- 0.17 g/cu cm for a detailed semianalytical model. The quoted uncertainties stem from an assumed uncertainty in the perijove radius. However, the uncertainty introduced by the different theoretical models is the major source of error; this uncertainty could be eliminated by future SPH simulations specialized to cometary disruptions, including the effects of initially prolate, spinning comets. If the SPH-based upper bound turns out to be most appropriate, it would be consistent with the predisruption body being a comet with a relatively low density and porous structure, as has been asserted previously based on observations of cometary outgassing. Regardless of which upper bound is preferable, the models all agree that the predisruption body could not have been a relatively high-density body, such as an asteroid with rho approximately = 2 g/cu cm.

Boss, Alan P.↗

On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds

The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.

Klein, Richard I.↗

Cooling flows and the stability of radio jets

This paper considers what happens to radio jets which, propagating out from the nucleus of an elliptical galaxy, encounters a point at which the ambient gas pressure increases very rapidly. The pressure variation within the sonic radius of a cooling flow is calculated, and a semianalytical calculation of the stability of a jet which encounters a rapidly increasing pressure gradient is given. The results of numerical hydrodynamic simulations which test the results of the analytic calculations are shown. The case where the jet encounters the sonic region at an oblique angle is considered, and it is suggested that this may provide an explanation for jets which are bent on small scales. These results are applied to radio sources associated with central dominant cluster galaxies and with more normal elliptical galaxies.

Soker, Noam↗

Near-critical spherical accretion by neutron stars - General relativistic treatment

Time-independent spherical accretion by a neutron star is studied using general relativistic radiation hydrodynamics. Numerical integrations of the flow equations are presented. These show that when the luminosity is sufficiently close to (but below) the Eddington limit, the flow velocity increases with decreasing radius far from the neutron star, reaches a maximum at an intermediate radius, and decreases at small radii. A large fraction of the binding energy of the flow is transferred to the radiation through scattering before the flow strikes the surface of the neutron star. Following Miller's treatment of accretion at luminosites near the Eddington limit (which neglected general relativistic effects), analytic approximations for the decelerating phase of the flow's velocity profile are derived. The dependence of the solutions on the variable Eddington factor prescription chosen to close the radiation moment equations is also examined.

Park, Myeong-Gu↗

Latitudinal Dependence of the Radial IMF Component: Coronal Imprint

Measurements by Ulysses have confirmed that there is no significant gradient with respect to heliomagnetic latitude in the radial component, B(sub r,) of the interplanetary magnetic field. In the corona, the plasma, beta is much less than 1, except directly above streamers, so longitudinal and latitudinal gradients in field strength will relax due to the transverse magnetic pressure gradient force as the solar wind carries magnetic flux away from the Sun. This happens quickly enough so that the field is essentially uniform by 5 - 10 solar radius, apparently remaining so as it is carried to beyond 1 AU. Here, we illustrate the coronal relaxation with a qualitative physical argument and by reference to a detailed Magneto HydroDynamics (MHD) simulation.

Suess, S. T.↗

Squeeze effects in radial face seals

Squeeze effects in a liquid lubricated radial face seal are analyzed. The analysis considers face misalignment with both axial and angular vibrations of the primary seal ring. Translational, rotational, and cross-coupled damping coefficients of the fluid film are derived analytically from a solution of the Reynolds equation utilizing the narrow seal approximation. Results are given for a wide range of practical radius ratios. At each radius ratio, the complete range of angular misalignment - from parallel faces to touch down - is covered. It is shown that squeeze effects in face seals are usually larger than the more familiar hydrodynamic effects. These effects play an important role in the seal's mechanism of operation and therefore have to be considered in any realistic seal model.

Etsion, I.↗

The formation phase of the solar nebula

Hydrodynamical calculations of the collapse of an axisymmetric, rotating protostellar cloud, with radiation transport and without magnetic fields, are presented. The collapse is assumed to start from a centrally condensed sphere of radius 5 x 10 to the 15th cm, a mean density of 4 x 10 to the -15th g/cu cm, a total mass of 1 solar mass, and a total angular momentum of 10 to the 53rd g sq cm per sec. The numerical grid is chosen to resolve the region of disk formation between 1 and 60 AU from the center. Frequency-dependent radiative transfer calculations show how the emergent spectrum of the structure depends upon viewing angle with respect to the rotation axis and how the observed isophotal contours should depend on wavelength and viewing angle. The central part of the protostar, interior to 1 AU, is not resolved numerically but is modeled approximately. At the end of the calculation, this region is found to have a mass of 0.6 solar mass and a ratio of rotational to gravitational energy of about 0.4, sufficiently large to be unstable to nonaxisymmetric perturbations. Although the disk is gravitationally stable according to the local Toomre criterion, the nonaxisymmetric structure in the center is likely to lead to angular momentum transport.

Bodenheimer, Peter↗

Propagation of solar disturbances in interplanetary space

Time-dependent solutions of a one-fluid model of the interplanetary medium are investigated. This set of unsteady hydrodynamic equations has been written in conservation form in order to apply the Lax-Wendroff method for the solution of this problem. The initial condition is specified by a pulse at 1 solar radius. The equilibrium condition is chosen to be the steady solution of a quiet solar wind. The specified solar disturbances in this calculation are allowed to be both sub- and supersonic by the present theoretical formulation. The results are presented in terms of density, velocity, and temperature profiles of the interplanetary gas flow at heliocentric distances up to about 10 AU at any particular time. The trajectories of disturbances for various initial pulses are shown. Some 1972 solar-flare observational data are compared with these theoretical calculations. From these calculations, the effects on the interplanetary environment, due to the propagation of solar disturbances, can be determined.

Wu, S. T.↗

Negative Ion Chemistry in the Coma of Comet 1P/Halley

Negative ions (anions) were identified in the coma of comet 1P/Halley from in-situ measurements performed by the Giotto spacecraft in 1986. These anions were detected with masses in the range 7-110 amu, but with insufficient mass resolution to permit unambiguous identification. We present details of a new chemical-hydrodynamic model for the coma of comet Halley that includes - for the first time - atomic and molecular anions, in addition to a comprehensive hydrocarbon chemistry. Anion number densities arc calculated as a function of radius in the coma, and compared with the Giotto results. Important anion production mechanisms arc found to include radiative electron attachment, polar photodissociation, dissociative electron attachment, and proton transfer. The polyyne anions C4H(-) and C6H(-) arc found to be likely candidates to explain the Giotto anion mass spectrum in the range 49-73 amu. Thc CN(-) anion probably makes a significant contribution to the mass spectrum at 26 amu. Larger carbon-chain anions such as C8H(1) can explain the peak near 100 amu provided there is a source of large carbon-chain-bearing molecules from the cometary nucleus.

Cordiner, M. A.↗

Cosmological implications of ROSAT observations of groups and clusters of galaxies

We have combined ROSAT Position Sensitive Proportional Counter (PSPC) and optical observations of a sample of groups and clusters of galaxies to determine the fundamental parameters of these systems (e.g., the dark matter distribution, gas mass fraction, baryon mass fraction, mass-to-light ratio, and the ratio of total-to-luminous mass). Imaging X-ray spectroscopy of groups and clusters show that the gas is essentially isothermal beyond the central region, indicating that the total mass density (mostly dark matter) scales as rho(sub dark) varies as 1/r squared. The density profile of the hot X-ray emitting gas is fairly flat in groups with rho(sub gas) varies as 1/r and becomes progressively steeper in hotter richer systems, with rho(sub gas) varies as 1/r squared in the richest clusters. These results show, that in general, the hot X-ray-emitting gas is the most extended mass component in groups and clusters, the galaxies are the most centrally concentrated component, and the dark matter is intermediate between the two. The flatter density rofile of the hot gas compared to the dark matter produces a gas mass fraction that increases with radius within each object. There is also a clear trend of increasing gas mass fraction (from 2% to 30%) between elliptical galaxies and rich clusters due to the greater detectable extent of the X-ray emission in richer systems. For the few systems in which the X-ray emission can be traced to the virial radius (where the overdensity delta is approximately equal 200), the gas mass fraction (essentially the baryon mass fraction) approaches a roughly constant value of 30%, suggesting that this is the true primordial value. Based on standard big bang nucleosynthesis, the large baryon mass fraction implies that Omega = 0.1 - 0.2. The antibiased gas distribution suggests that feedback from galaxy formation and hydrodynamics play important roles in the formation of structure on the scale of galaxies to rich clusters. All the groups and clusters in our sample have mass-to-light ratios of M/L(sub V) approximately 100 - 150 solar mass/solar luminosity, which strongly contrasts with the traditional view that the mass-to-light ratio of rich clusters is significantly greater than individual galaxies or groups with M/L(sub V) approximately 250 - 300 solar mass/solar luminosity. We also show that M/L(sub V is essentially constant within the virial radius of clusters (where delta is greater than or approximately 200), which is consistent with the peaks formalism of biased galaxy formation. While the mass-to-light ratios of groups and clusters are comparable (indicating a constant mass fraction of optically luminous material), the ratio of the total mass-to-luminous mass (gas plus stars) monotonically decreases between galaxies and clusters. The decrease in M(sub total)/M(sub lum) arises from two factors: (1) the composition of baryonic matter varies from a predominance of optically luminous material (stars) on the scale of galaxies (approximately 10 kpc) to a predominance of X-ray luminous material (hot gas) on the scale of rich clusters (approximately 1 Mpc), and (2) the hot gas has a more extended spatial distribution than the gravitating matter. The observed decrease M(sub total)/M(sub lum) between galaxies and clusters indicates that the universe actually becomes `brighter' on mass scales between 10(exp 12) and 10(exp 15) solar mass, in the sense that a greater fraction of the gravitating mass is observable.

David, Laurence P.↗

The clumpy circumstellar medium around young supernova remnants

Each of the youngest supernova remnants known in the Milky Way, Cas A, Kepler's SNR, and Tycho's SNR, shows a different morphological structure caused by different conditions in the progenitor stars and their surroundings. In all three cases, however, the observed shells have a thickness of about 1/4 the radius, a sharp outer edge, and significant brightness irregularities. These features require that the circumstellar medium be highly clumped. To investigate the phenomenon, models of the expansion have been constructed using a one-dimensional spherical hydrodynamic code. As a supernova shock moves down the external density gradient of the star, material behind the shock begins to go into free expansion. Then as surrounding material is encountered a reverse shock moving back into the ejectum will be formed. Until the expansion has swept up about eight times the ejected mass when the situation can be considered as a point explosion in its surroundings, the dynamics are controlled by conditions between the shocks. The region is also where the synchrotron radio emission from relativistic electrons trapped in magnetic fields arises. Initial particles and fields are accelerated and amplified by eddy motion at the interface between the ejected and swept-up material and at the boundaries of clumps. Polarimetry shows that these SNR have a net radial orientation of their magnetic fields apparently from stretching by Rayleigh-Taylor instabilities at the contact surfaces. Without clumps the observed shell is much too narrow and steep on the inside.

Dickel, J. R.↗

Numerical simulation of nonaxisymmetric adiabatic accretion flow

The hydrodynamics of gas flow past a finite-sized gravitating central object is studied in the adiabatic approximation in two spatial dimensions. The flow morphology is shown to depend on the ratio (epsilon) of the axisymmetric accretion radius to the length scale of the density variation. For small to intermediate epsilon values the flows are highly time dependent, and for large epsilon values a quasi-steady state is obtained.

Fryxell, B. A.↗