Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrodynamic stability”

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.

281 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗
NASA NTRS2009

Liquid-Propellant Rocket Engine Throttling: A Comprehensive Review

Liquid-Propellant Rocket Engines (LREs) are capable of on-command variable thrust or thrust modulation, an operability advantage that has been studied intermittently since the late 1930s. Throttleable LREs can be used for planetary entry and descent, space rendezvous, orbital maneuvering including orientation and stabilization in space, and hovering and hazard avoidance during planetary landing. Other applications have included control of aircraft rocket engines, limiting of vehicle acceleration or velocity using retrograde rockets, and ballistic missile defense trajectory control. Throttleable LREs can also continuously follow the most economical thrust curve in a given situation, compared to discrete throttling changes over a few select operating points. The effects of variable thrust on the mechanics and dynamics of an LRE as well as difficulties and issues surrounding the throttling process are important aspects of throttling behavior. This review provides a detailed survey of LRE throttling centered around engines from the United States. Several LRE throttling methods are discussed, including high-pressure-drop systems, dual-injector manifolds, gas injection, multiple chambers, pulse modulation, throat throttling, movable injector components, and hydrodynamically dissipative injectors. Concerns and issues surrounding each method are examined, and the advantages and shortcomings compared.

Casiano, Matthew↗
NASA NTRS1994

A waved journal bearing concept with improved steady-state and dynamic performance

Analysis of the waved journal bearing concept featuring a waved inner bearing diameter for use with a compressible lubricant (gas) is presented. A three wave, waved journal bearing geometry is used to show the geometry of this concept. The performance of generic waved bearings having either three, four, six, or eight waves is predicted for air lubricated bearings. Steady-state performance is discussed in terms of bearing load capacity, while the dynamic performance is discussed in terms of dynamic coefficients and fluid film stability. It was found that the bearing wave amplitude has an important influence on both steady-state and dynamic performance of the waved journal bearing. For a fixed eccentricity ratio, the bearing steady-state load capacity and direct dynamic stiffness coefficient increase as the wave amplitude increases. Also, the waved bearing becomes more stable as the wave amplitude increases. In addition, increasing the number of waves reduces the waved bearing's sensitivity to the direction of the applied load relative to the wave. However, the range in which the bearing performance can be varied decreases as the number of waves increases. Therefore, both the number and the amplitude of the waves must be properly selected to optimize the waved bearing design for a specific application. It is concluded that the stiffness of an air bearing, due to the hydrodynamic effect, could be doubled and made to run stably by using a six or eight wave geometry with a wave amplitude approximately half of the bearing radial clearance.

Dimofte, Florin↗
NASA NTRS1994

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗
NASA NTRS2003

Mixing and Transport in the Solar Nebula

Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.

Boss, Alan P.↗
NASA NTRS1993

The formation of protostellar disks. I - 1 M(solar)

Hydrodynamical calculations of the collapse of an axisymmetric, rotating one solar mass protostellar cloud, including the effects of radiative transfer and radiative acceleration but without magnetic fields, are presented. The results include calculations of infrared protostellar spectra as a function of time and viewing angle. A numerical algorithm involving explicit nested grids is used to resolve the region of initial disk formation and at the same time to include the outer regions in the calculation. The central part of the protostar is modeled approximately. Initial conditions are systematically varied to investigate their influence on the evolution and final configuration of central star plus circumstellar disk. The initial state for the standard case is a centrally condensed molecular cloud core of one solar mass with a mean density of 8 x 10 exp -18 g/cu cm and a specific angular momentum at the outer edge of 7 x 10 exp 20 sq cm/s. The collapse is followed for 8 x 10 exp 4 yr, at which point 0.45 solar mass is contained in a rapidly rotating central object and most of the remainder in a surrounding equilibrium disk. The stability of this final structure is qualitatively analyzed.

Yorke, Harold W.↗
NASA NTRS2022

Hypersonic Boundary-Layer Instabilities over Ogive-Cylinder Models

Computational investigations of an ogive-cylinder geometry with varying nosetips at zero degrees angle of attack are presented. The model geometry and conditions are selected to match experiments conducted in the Air Force Research Laboratory (AFRL) Mach 6 Ludwieg Tube. Five nosetips of interest were selected for the computational studies herein: two sharp ogives, two blunt ogives, and one hemispherical nosetip. Computations are performed at a freestream Reynolds number of 7.01 × 10 6 m -1 . Each sharp and blunt tip ogives had a 14 and 28 degree version for the tip angles. A cylindrical section follows the nosetip, resulting in a meter long model such as in the experiments. The laminar flow solutions are analyzed. The boundary-layer-edge properties and the velocity and temperature profiles are compared across streamwise locations aft of the ogive-cylinder junction. The blunt nosetips induce an entropy layer that envelopes the boundary-layer profiles. Modal stability analysis identifies most amplified frequencies corresponding to Mack’s second modes that agree with experimental results for the sharp tips. Similar to the experimental measurements based on wall-mounted pressure sensors, no unstable modes are found for the blunt models. Nonmodal analysis revealed a broadband set of disturbances present for the blunter tips, in agreement with experimental observations. Flow perturbation contours of most amplified planar and oblique disturbances are shown to qualitatively match wind tunnel Schlieren images, with a switch from rope-like to elongated structures, i.e., from high frequency Mack’s second modes to low frequency Mack’s first modes, as the nosetip angle is increased for the sharp tip, and from boundary-layer to entropy-layer disturbances as the bluntness is increased.

Boundary layer transition↗
NASA NTRS2022

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

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

Hydrodynamic instability↗
NASA NTRS1994

A hydrodynamic approach to cosmology: The mixed dark matter cosmological scenario

We compute the evolution of spatially flat, mixed cold and hot dark matter models containing both baryonic matter and two kinds of dark matter. Hydrodynamics is treated with a highly developed Eulerian hydrodynamic code (see Cen 1992). A standard particle-mesh (PM) code is also used in parallel to calculate the motion of the dark matter components. We adopt the following parameters: h equivalent to (sub 0)/100 km/s Mpc(exp -1) = 0.5, OMEGA(sub C) = 0.3, and OMEGA(sub B) = 0.06, with amplitude of the perturbation spectrum fixed by the Cosmic Background Explorer Satellite (COBE) Dark Matter Radiation (DMR) measurements (Smoot et al. 1992) being sigma (sub 8) = 0.67. Four different boxes are simulated with box sizes of L = (64, 16, 4, 1) h(exp -1) Mpc, respectively, the two small boxes providing good resolution but little valid information due to absence of large-scale power. We use 128(exp 3) approximate 10(exp 6.3) baryonic cells, 128(exp .3) cold dark matter particles, and 2 x 128(exp 3) hot dark matter particles. In addition to the dark matter we follow separately six baryonic species (H, H(+), He, He(+), He(++), e(-)) with allowance for both (nonequilibrium) collisional and radiative ionization in every cell. The background radiation field is also followed in detail with allowance made for both continuum and line processes, to allow nonequilibrium heating and cooling processes to be followed in detail. The mean final Zeldovich-Sunyaev y parameter is estimated to be y Bar = (5.4 + or - 2.7) x 10(exp -7) below currently attainable observations, with a rms fluctuation of approximately delta bar y = (0.6 + or - 3.0) x 10(exp -7) on arcminute scales. The rate of galaxy formation peaks at an even later epoch (z approximate 0.3) than in the standard (OMEGA = 1, sigma sub 8 = 0.67) cold dark matter (CDM) model (z approximate 0.5) and, at a redshift of z = 4, is nearly a factor of 100 lower than for the CDM model with the same value of sigma sub 8. With regard to mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum has a broad peak in the vicinity of M(sub B) = 10(exp 9.5) solar mass with a reasonable fit to the Schechter luminosity function if the ratio of baryon mass to blue light is approximately 4. In addition, one very large PM simulation was made in a box with size (320 h(exp - 1) Mpc) containing 3 x 200(exp 3) = 10(exp 7.4) particles. Utilizing this simulation we find that the model yields a cluster mass function which is about a factor of 4 higher than observed, but a cluster-cluster correlation length marginally lower than observed, but that both are closer to observations than in the (COBE) normalized CDM model. The one-dimensional pairwise velocity dispersion is 605 + or - 8 km/s at 1/h separation, lower than that of the DCM model normalized to COBE, but still significant higher than observations (Davis & Peebles 1983). A plausible velocity bias b(sub v) = 0.8 + or - 0.1 on this scale will reduce but not remove the discrepancy. The velocity auto-correlat ion function has a coherence length of 40/h Mpc, which is somewhat lower than the observed counterpart. In all these respects the model would be improved by decreasing the cold fraction of the dark OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDB). But formation of galaxies and clusters of galaxies is much later in this model than in COBE-normalized CDM, perhaps too late. To improve on these constraints a larger ratio of OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDM)) is required than the value of 0.67 adopted here. It does not seem possible to find a value for this ratio which would satisfy all tests. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of sigma(sub B), but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter, significantly improves the ability of the COBE-normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same sigma(sub 8) and is still probably unsatisfactory with regard to several critical tests.

Cen, Renyue↗
NASA NTRS2002

Containerless Ripple Turbulence

One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Putterman, Seth↗
NASA NTRSDate not supplied

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗