SEARCH · Search NASA
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.
The resistance of water, even at low speeds, may be such as to render it impossible to take off. For this type of aircraft, it is therefore especially important to determine the take-off characteristics. A new and simple graphical method is given and appears to be the only one for the complete solution of the problem.
A hydrodynamic model was developed of a spherically symmetric dusty gas flow in a cometary atmosphere assuming a single fluid, inviscid, perfect gas. The hydrodynamics for gas and dust, which involves the gas drag force (momentum transfer), heat exchange between gas and dust, photodissociation energy for H2O gas, and radiative heating and cooling terms for dust particles, are solved using the Gear method for stiff, coupled differential equations. Calculations were done with a dust size distribution for radii alpha = 0.01 micron to 10 cm with densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g/cu cm and a total dust-to-gas mass ratio chi = 1 was adopted. There are indications from in situ observations that dust particle fragment into smaller ones. Fragmentation of dust particles was incorporated into the model. This is done by adding source and sink terms in the continuity equations for the dust. Lifetimes for the decay of dust particles were assumed as a function of particle size. It is also assumed that dust particles always fragment only into the next smaller size.
Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3–4 μm), M (4.5–5 μm), and N (8–13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes OPTIM3D and RADMC3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N-band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of VX Sgr’s stellar surface and circumstellar environment across a very large spectral domain for the first time.
Weakly nonlinear Magneto Hydrodynamic (MHD) stability of the Halley cometosheath determined by the balance between the outward ion-neutral drag force and the inward Lorentz force is investigated including the transverse plasma motion as observed in the flanks with the help of the method of multiple scales. The eigenvalues and the eigenfunctions are obtained for the linear problem and the time evolution of the amplitude is obtained using the solvability condition for the solution of the second order problem. The diamagnetic cavity boundary and the adjacent layer of about 100 km thickness is found unstable for the travelling waves of certain wave numbers. Halley ionopause has been observed to have strong ripples with a wavelength of several hundred kilometers. It is found that nonlinear effects have stabilizing effect.
The purpose of this paper is to study the dynamic behavior of a three-wave journal bearing using a transient approach. The transient analysis permits the determination of the rotor behavior after the fractional frequency whirl appears. The journal trajectory is determined by solving a set of nonlinear equations of motion using the Runge-Katta method. The fluid film forces are computed by integrating the transient Reynolds equation at each time step location of the shaft with respect to the bearing. Because of the large values of the rotational speeds, turbulent effects were included in the computations. The influence of the temperature on the viscosity was also considered. Numerical results were compared to experimenta1 results obtained at the NASA Glenn Research Center. Comparisons of the theoretical results with experimental data were found to be in good agreement. The numerical and experimental results showed that the fluid film of a three-wave journal bearing having a diameter of 30 mm, a length of 27 mm, and a wave amplitude ratio greater than 0.15 is stable even at rotational speeds of 60,000 RPM. For lower wave amplitude ratios, the threshold speed at which the fluid film becomes unstable depends on the wave amplitude and on the supply pocket pressure. Even if the fluid film is unstable, the wave bearing maintains the whirl orbit inside the bearing clearance.
The problem of acoustic radiation generated by instability waves of a compressible plane turbulent shear layer is solved. The solution provided is valid up to the acoustic far-field region. It represents a significant improvement over the solution obtained by classical hydrodynamic-stability theory which is essentially a local solution with the acoustic radiation suppressed. The basic instability-wave solution which is valid in the shear layer and the near-field region is constructed in terms of an asymptotic expansion using the method of multiple scales. This solution accounts for the effects of the slightly divergent mean flow. It is shown that the multiple-scales asymptotic expansion is not uniformly valid far from the shear layer. Continuation of this solution into the entire upper half-plane is described. The extended solution enables the near- and far-field pressure fluctuations associated with the instability wave to be determined. Numerical results show that the directivity pattern of acoustic radiation into the stationary medium peaks at 20 degrees to the axis of the shear layer in the downstream direction for supersonic flows. This agrees qualitatively with the observed noise-directivity patterns of supersonic jets.
The author suggests that since clouds are born co-moving in a turbulent intra-cluster medium (ICM), the allowed parameter space can now be opened up to a more acceptable range. Large-scale motions can be driven in the central parts of cooling flows by a number of mechanisms including the motion of the central and other galaxies, and the dissipation of advected, focussed rotational and magnetic energy. In addition to the velocity width paradox, two other paradoxes (Heckman et al. 1989) can be solved if the ICM is turbulent. Firstly, the heating source for the emission line regions has always been puzzling - line luminosities are extremely high for a given (optical or radio) galaxy luminosity compared to those in non-cooling flow galaxies, therefore a mechanism peculiar to cooling flows must be at work. However most, if not all, previously suggested heating mechanisms either fail to provide enough ionization or give the wrong line ratios, or both. The kinetic energy in the turbulence provides a natural energy source if it can be efficiently converted to cloud heat. Researchers suggest that this can be done via magneto-hydrodynamic waves through plasma slip. Secondly, while the x ray observations indicate extended mass deposition, the optical line emission is more centrally concentrated. Since many of the turbulence-inducing mechanisms are strongest in the central regions of the ICM, so is the method of heating. In other words material is dropping out everywhere but only being lit up in the center.
Due to complexity of the program which used a hydrodynamic, axially and radially deformable disk-ring model and the resulting long computing time only the output gap was investigated. Results from independent studies were used to initiate the starting conditions for the electrons and the RF voltage using our program. Although this method of computation is less exact than processing the entire klystron interaction in 3-Dimensions it is shown that, for a confined flow focused throughout the penultimate cavity, radial velocities remain very small and the beam is highly laminar. It is concluded that possible errors resulting from treating only the output cavity in 3-D would remain small.
For aerodynamic and hydrodynamic vehicles, it is highly desirable to reduce drag and noise levels. A reduction in drag leads to fuel savings. In particular for submersible vehicles, a decrease in noise levels inhibits detection. A suggested means to obtain these reduction goals is by delaying the transition from laminar to turbulent flow in external boundary layers. For hydrodynamic applications, a passive device which shows promise for transition delays is the compliant coating. In previous studies with a simple mechanical model representing the compliant wall, coatings were found that provided transition delays as predicted from the semi-empirical e(sup n) method. Those studies were concerned with the linear stage of transition where the instability of concern is referred to as the primary instability. For the flat-plate boundary layer, the Tollmien-Schlichting (TS) wave is the primary instability. In one of those studies, it was shown that three-dimensional (3-D) primary instabilities, or oblique waves, could dominate transition over the coatings considered. From the primary instability, the stretching and tilting of vorticity in the shear flow leads to a secondary instability mechanism. This has been theoretical described by Herbert based on Floquet theory. In the present study, Herbert's theory is used to predict the development of secondary instabilities over isotropic and non-isotropic compliant walls. Since oblique waves may be dominant over compliant walls, a secondary theory extention is made to allow for these 3-D primary instabilities. The effect of variations in primary amplitude, spanwise wavenumber, and Reynolds number on the secondary instabilities are examined. As in the rigid wall case, over compliant walls the subharmonic mode of secondary instability dominates for low-amplitude primary disturbances. Both isotropic and non-isotropic compliant walls lead to reduced secondary growth rates compared to the rigid wall results. For high frequencies, the non-isotropic wall suppresses the amplification of the secondary instabilities, while instabilities over the isotropic wall may grow with an explosive rate similar to the rigid wall results. For the more important lower frequencies, both isotropic and non-isotropic compliant walls suppress the amplification of secondary instabilities compared to the rigid wall results. The twofold major discovery and demonstration of the present investigation are: (1) the use of passive devices, such as compliant walls, can lead to significant reductions in the secondary instability growth rates and amplification; (2) suppressing the primary growth rates and subsequent amplification enable delays in the growth of the explosive secondary instability mechanism.
Preliminary analysis showed that the doubly curved triangular shell elements were too stiff for these shell structures. The doubly curved quadrilateral shell elements were found to give much improved results. A total of six load cases were analyzed in this study. The load cases were either those resulting from a static test using reaction straps to simulate the drop conditions or under assumed hydrodynamic conditions resulting from a drop test. The latter hydrodynamic conditions were obtained through an emperical fit of available data. Results obtained from a linear analysis were found to be consistent with results obtained elsewhere with NASTRAN and BOSOR. The nonlinear analysis showed that the originally assumed loads would result in failure of the shell structures. The nonlinear analysis also showed that it was useful to apply internal pressure as a stabilizing influence on collapse. A final analysis with an updated estimate of load conditions resulted in linear behavior up to full load.
Observations of the interplanetary medium are considered together with the structure of the solar atmosphere and the kinetic properties of the interplanetary plasma. Methods of gas kinetics and continuum flows are examined, taking into account microscopic and macroscopic representations, the kinetic flow equations, the formal solution of the Vlasov equation, and the continuum flow equations. The collective particle behavior of that interplanetary gas is discussed along with the hydrodynamic coronal expansion, a wave-pump problem, the free expansion phenomenon, and a generalized free expansion problem.
This project investigates the origin of giant planets, both in the Solar System and around other stars. It is assumed that the planets form by the core accretion process: small solid particles in a disk surrounding a young star gradually coagulate into objects of a few kilometers in size, known as planetesimals, which then accumulate into solid protoplanetary cores. Once the cores have become large enough, they are able to attract gas from the surrounding disk to form the deep gaseous envelope of the giant planet. Our code simulates giant planet growth in a spherical approximation, and it has been quite successful in addressing a number of basic planetary properties. Further improvements to the code have been made to achieve a more realistic understanding of planetary formation. The computations of the models were based on an earlier version of our code and were stopped at the onset of runaway gas accretion. Now, improved boundary conditions have been incorporated into the code to allow for hydrodynamic inflow of gas and to handle the late stages of evolution when the planet evolves at constant mass. These changes were made to the version of the code that uses a constant accretion rate and to the version that uses a self-consistent method for calculating both the solid and gas accretion rates. The equation of state has been updated to incorporate the detailed tables of Saumon, Chabrier, and Van Horn. The opacities were updated to include the results of Alexander and Ferguson. The outer boundary conditions were modified. During the accretion phase when the planet's radius is between the accretion radius and the tidal radius, we set the outer boundary at a 'modified' accretion radius, which is the point where thermal energy is enough to bring gas to the edge of the Hill sphere.
Scattering from the ocean surface at microwave frequencies is often described by the Bragg scattering model or the two-scale model. Here, the ability of these models to predict cross section modulation is investigated by comparing Bragg and two-scale cross-sections with cross sections calculated by a numerical electromagnetic approach. The numerical method is used as the standard of comparison. For slightly rough surfaces, the Bragg model is found to be very accurate, but for much rougher surfaces, the two-scale model is superior and should be used for calculating both the absolute cross section and the change in cross section due to spatial or temporal modulation of the ocean wave height spectrum. These results are relevant to the interpretation of SAR images of ocean surface phenomena which perturb the wave height spectrum at ripple wavelengths.
A brief overview of some theoretical and computational studies of the stability of streamwise vortices is given. The local induction model and classical hydrodynamic vortex stability theories are discussed in some detail. The importance of the three-dimensionality of the mean velocity profile to the results of stability calculations is discussed briefly. The mean velocity profile is provided by employing the similarity solution of Donaldson and Sullivan. The global method of Bridges and Morris was chosen for the spatial stability calculations for the nonlinear eigenvalue problem. In order to test the numerical method, a second order accurate central difference scheme was used to obtain the coefficient matrices. It was shown that a second order finite difference method lacks the required accuracy for global eigenvalue calculations. Finally the problem was formulated using spectral methods and a truncated Chebyshev series.
A three-dimensional hydrodynamic computer program is used to study the growth of nonaxisymmetric structures in rapidly rotating, self-gravitating polytropes. Models with polytropic index n = 0.8, 1.0, 1.3, 1.5, and 1.8 are studied. The initially axisymmetric equilibria are constructed by the Ostriker-Mark self-consistent-field method. The nonaxisymmetric pattern that develops out of low-amplitude random noise is a two-armed spiral with a well-defined pattern speed and growth rate which closely match properties of the toroidal mode predicted from the linear, second-order tensor-virial equation. A Fourier analysis of each polytrope's azimuthal density distribution shows that, even in the linear amplitude regime, higher-order angular patterns also develop exponentially in time. The higher-order patterns ultimately move in synchronization with the broad two-armed spiral, creating a narrow two-armed spiral. As the polytropic index is decreased, a more open and centrally more barlike pattern develops.
We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids, and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest, such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.
Context. We investigate iPTF13bvn, a core-collapse (CC) supernova (SN) in the nearby spiral galaxy NGC 5806. This object was discovered by the intermediate Palomar Transient Factory (iPTF) very close to the estimated explosion date and was classified as a stripped-envelope CC SN, likely of Type Ib. Furthermore, a possible progenitor detection in pre-explosion Hubble Space Telescope (HST) images was reported, making this the only SN Ib with such an identification. Based on the luminosity and color of the progenitor candidate, as well as on early-time spectra and photometry of the SN, it was argued that the progenitor candidate is consistent with a single, massive Wolf-Rayet (WR) star. Aims. We aim to confirm the progenitor detection, to robustly classify the SN using additional spectroscopy, and to investigate if our follow-up photometric and spectroscopic data on iPTF13bvn are consistent with a single-star WR progenitor scenario. Methods. We present a large set of observational data, consisting of multi-band light curves (UBVRI, g′r′i′z′) and optical spectra. We perform standard spectral line analysis to track the evolution of the SN ejecta. We also construct a bolometric light curve and perform hydrodynamical calculations to model this light curve to constrain the synthesized radioactive nickel mass and the total ejecta mass of the SN. Late-time photometry is analyzed to constrain the amount of oxygen. Furthermore, image registration of pre- and post-explosion HST images is performed. Results. Our HST astrometry confirms the location of the progenitor candidate of iPTF13bvn, and follow-up spectra securely classify this as a SN Ib. We use our hydrodynamical model to fit the observed bolometric light curve, estimating the total ejecta mass to be 1.9 solar mass and the radioactive nickel mass to be 0.05 solar mass. The model fit requires the nickel synthesized in the explosion to be highly mixed out in the ejecta. We also find that the late-time nebular r′-band luminosity is not consistent with predictions based on the expected oxygen nucleosynthesis in very massive stars. Conclusions. We find that our bolometric light curve of iPTF13bvn is not consistent with the previously proposed single massive WR-star progenitor scenario. The total ejecta mass and, in particular, the late-time oxygen emission are both significantly lower than what would be expected from a single WR progenitor with a main-sequence mass of at least 30 solar mass.
In this paper we construct higher order Godunov schemes for isothermal flow. Isothermal hydrodynamics serves as a good representation for several systems of astrophysical interest. The schemes designed here have second-order accuracy in space and time and some are third-order accurate for advection. Moreover, several ingredients of these schemes are essential components of even higher order. The methods designed here have excellent ability to represent smooth flow yet capture shocks with high resolution. Several test problems are presented. The algorithms presented here are compared with other algorithms having a comparable formal order of accuracy.