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.
Hydrodynamic differential and difference equations and computer algorithm for nonlinear combustion instability in liquid propellant rocket engines, and Fortran IV listing
In the early stages of the J-2X upper stage engine program, various inducer configurations proposed for use in the fuel turbopump (FTP) were tested in water. The primary objectives of this test effort were twofold. First, to obtain a more comprehensive data set than that which existed in the Pratt & Whitney Rocketdyne (PWR) historical archives from the original J-2S program, and second, to supplement that data set with information regarding the cavitation induced vibrations for both the historical J-2S configuration as well as those tested for the J-2X program. The J-2X FTP inducer, which actually consists of an inducer stage mechanically attached to a kicker stage, underwent 4 primary iterations utilizing sub-scaled test articles manufactured and tested in PWR's Engineering Development Laboratory (EDL). The kicker remained unchanged throughout the test series. The four inducer configurations tested retained many of the basic design features of the J-2S inducer, but also included variations on leading edge blade thickness and blade angle distribution, primarily aimed at improving suction performance at higher flow coefficients. From these data sets, the effects of the tested design variables on hydrodynamic performance and cavitation instabilities were discerned. A limited comparison of impact to the inducer efficiency was determined as well.
Instead of computing a large number of coarsely zoned hydrodynamic models covering the entire atmospheric instability strip, the author computed a single model as well as computer limitations allow. The implicit hydrodynamic code of Kutter and Sparks was modified to include radiative transfer effects in optically thin zones.
Calculations with a two-dimensional hydrodynamic simulation show that a generic Raleigh-Taylor-like instability occurs in the mantles of nascent neutron stars, that it is possibly violent, and that the standard spherically symmetric models of neutron star birth and supernova explosion may be inadequate. Whether this 'convective' instability is pivotal to the supernova mechanism, pulsar nagnetic fields, or a host of other important issues that attend stellar collapse remains to be seen, but its existence promises to modify all questions concerning this most energetic of astronomical phenomena.
The evolution of young supernova remnants has been modeled using a one-dimensional hydrodynamics code. Rayleigh-Taylor instabilities and the growth of Rayleigh-Taylor fingers have been included in the code. Turbulent dynamo amplification of magnetic fields and both turbulent and shock acceleration of relativistic electrons have been included macroscopically. From this, the distribution of synchrotron luminosity in the remnant has been calculated. It is found that the radio morphology of model remnants expanding into a homogeneous medium does not agree with observations. Expansion into a circumstellar medium with many small cloudlets does produce radio shells which agree with observations. It is suggested that supernova remnants reflect the interaction of ejected matter with a cloudy circumstellar medium.
Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.
The Ingard-Myers impedance boundary condition is widely recognized to be prone to hydrodynamic Kelvin-Helmholtz-type instability, primarily due to its use of a vortex sheet to model the flow at the boundary. Recently, a stabilized approximation of the Ingard-Myers condition, called the Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC), has been introduced. This approximation offers a good representation of the Ingard-Myers condition for grazing flows characterized by low to mid subsonic Mach numbers. This paper explores an extension of the TIMIBC, referred to as TIMIBC-ext, designed for liners in grazing flows with high subsonic Mach numbers. The TIMIBC-ext introduces a tunable parameter and demonstrates that by judiciously selecting the parameter value, the accuracy of the TIMIBC can be improved for either upstream or downstream propagating waves in grazing flows with high subsonic Mach numbers. Consequently, the TIMIBC-ext represents an enhancement to the TIMIBC when the location of the source relative to the liner (upstream or downstream) is known. Time domain implementations of the TIMIBC-ext are also provided, employing a multipole expansion model for the impedance function. A numerical example is provided to assess the performance of the TIMIBC-ext, comparing computational results with measurements from a recent set of measurements in the NASA Langley Grazing Flow Impedance Tube.
The paper develops the theory of the hydrodynamic stability of an electrolytic cell as a function of the imposed electric current. A new electrohydrodynamic instability is encountered when the current is forced to exceed the Nernst limit. The convection is driven by the volume force exerted by the electric field on space charges in the electrolyte. This intrinsic instability is found to be easily masked by extrinsic convection sources such as gravity or stirring. A linear stability analysis is performed and a dimensionless number Le is derived whose value determines the convection pattern.
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.
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.
The instability of a neutron star at the minimum mass is investigated using a hydrodynamical model to follow the evolution of the unstable star. A detailed analysis of the linear stability of the equilibrium configurations near the minimum mass is performed by solving the radial eigenvalue problem for the fundamental mode. The minimum mass configuration M(mmc) of 0.196 solar mass is found to be stable to small perturbations. Stability to radial perturbations is lost only at a lower critical mass of 0.16 solar mass, corresponding to about 0.8 M(mmc). It is shown that a simple dynamical model constructed using a 3-polytrope equation of state for hot dense matter reproduces the key dynamical features of the instability in the explosion phase.
The linear inviscid hydrodynamic stability of slightly curved free mixing layers is studied in this paper. The disturbance equation is solved numerically using a shooting technique. Two mean velocity profiles that represent stably and unstably curved free mixing layers are considered. Results are shown for cases of five curvature Richardson numbers. The stability characteristics of the shear layer are found to vary significantly with the introduction of the curvature effects. The results also indicate that, in a manner similar to the Goertler vortices observed in a boundary layer along a concave wall, instability modes of spatially developing streamwise vortex pairs may appear in unstable curved mixing layers.
The linear inviscid hydrodynamic stability of slightly curved free mixing layers is studied in this paper. The disturbance equation is solved numerically using a shooting technique. Two mean velocity profiles that represent stably and unstably curved free mixing layers are considered. Results are shown for cases of five curvature Richardson numbers. The stability characteristics of the shear layer are found to vary significantly with the introduction of the curvature effects. The results also indicate that, in a manner similar to the Goertler vortices observed in a boundary layer along a concave wall, instability modes of spatially developing streamwise vortex pairs may appear in centrifugally unstable curved mixing layers.
Hydrodynamic models for effects of thermocapillary instability, film thickness, and heat flux on breakdown of falling liquid films
Combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system. The problem of combustion-acoustic instability is a concern in many devices for various reasons, as each device may have a unique mechanism causing unsteady heat release rates and many have unique boundary conditions. To accurately predict and quantify combustion-acoustic stabilities, the unsteady heat release rate and boundary conditions need to be accurately determined. The present review brings together work performed on a variety of practical combustion devices. Many theoretical and experimental investigations of the unsteady heat release rate have been performed, some based on perturbations in the fuel delivery system particularly for rocket instabilities, while others are based on hydrodynamic processes as in ramjet dump combustors. The boundary conditions for rocket engines have been analyzed and measured extensively. However, less work has been done to measure acoustic boundary conditions in many other combustion systems.
We present the numerical techniques and tests used to construct and validate a computer code designed to study the multidimensional nonlinear hydrodynamics of large-scale sheet structures in the universe, especially the fragmentation of such structures under various instabilities. This code is composed of two codes, the hydrodynamical code ZEUS-2D and a particle-mesh code. The ZEUS-2D code solves the hydrodynamical equations in two dimensions using explicit Eulerian finite-difference techniques, with modifications made to incorporate the expansion of the universe and the gas cooling due to Compton scattering, bremsstrahlung, and hydrogen and helium cooling. The particle-mesh code solves the equation of motion for the collisionless dark matter. The code uses two-dimensional Cartesian coordinates with a nonuniform grid in one direction to provide high resolution for the sheet structures. A series of one-dimensional and two-dimensional linear perturbation tests are presented which are designed to test the hydro solver and the Poisson solver with and without the expansion of the universe. We also present a radiative shock wave test which is designed to ensure the code's capability to handle radiative cooling properly. And finally a series of one-dimensional Zel'dovich pancake tests used to test the dark matter code and the hydro solver in the nonlinear regime are discussed and compared with the results of Bond et al. (1984) and Shapiro & Struck-Marcell (1985). Overall, the code is shown to produce accurate and stable results, which provide us a powerful tool to further our studies.
Circumplanetary discs can be linearly unstable to the growth of disc tilt in the tidal potential of the star–planet system. We use3D hydrodynamical simulations to characterize the disc conditions needed for instability, together with its long-term evolution.Tilt growth occurs for disc aspect ratios, evaluated near the disc outer edge, ofH/r0.05, with a weak dependence on viscosityin the wave-like regime of warp propagation. Lower mass giant planets are more likely to have circumplanetary discs thatsatisfy the conditions for instability. We show that the tilt instability can excite the inclination to above the threshold where thecircumplanetary disc becomes unstable to Kozai–Lidov (KL) oscillations. Dissipation in the KL unstable regime caps furthertilt growth, but the disc experiences large oscillations in both inclination and eccentricity. Planetary accretion occurs in episodicaccretion events. We discuss implications of the joint tilt–KL instability for the detectability of circumplanetary discs, for theobliquity evolution of forming giant planets, and for the formation of satellite systems.
A multifluid hydrodynamic approximation allowing for the relative motion along the magnetic field of the newly created ions and the original fluid is used to treat the ion-pickup process. Due to the processes characterized by these means, the ion tail of a comet may not be antisolar; the derivation from radial is anticipated to be largest for oxygen due to its ionization at the greatest distances. Other ions, created nearer the comet where flow speed is lower, should have smaller transverse velocities.