A theory of axial-mode shock-wave oscillations in a solid-rocket combustor
Axial mode shock wave oscillations in end-burning solid rocket motor, obtaining periodic solutions for flow properties, solid temperature and burning rate
SEARCH · Engineering Papers
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.
Axial mode shock wave oscillations in end-burning solid rocket motor, obtaining periodic solutions for flow properties, solid temperature and burning rate
Experimental observations of instability of cavitating inducers were made for two different inducers operating at different flow coefficients. In general, instability occurred just before head breakdown. Auto-oscillation and rotating cavitation were observed. Analysis of small-amplitude behavior of the inducer and hydraulic system is carried out, and analytical predictions of stability limits were compared with experiment.
A wind tunnel technique for obtaining gust frequency response functions for use in predicting the response of flexible aircraft to atmospheric turbulence is evaluated. The tunnel test results for a dynamically scaled cable supported aeroelastic model are compared with analytical and flight data. The wind tunnel technique, which employs oscillating vanes in the tunnel throat section to generate a sinusoidally varying flow field around the model, was evaluated by use of a 1/30 scale model of the B-52E airplane. Correlation between the wind tunnel results, flight test results, and analytical predictions for response in the short period and wing first elastic modes of motion are presented.
The concept of the precessing accretion disk in HZ Her/Her X-1 in its varied forms, to account for the 35 day periodicity in the X-ray flux, has met many objections from a number of workers on various grounds, but it is still being invoked in current publications. These objections are reviewed and additional arguments are presented against the precessing accretion disk model. The implausibility of the disk models is demonstrated. An alternate clock mechanism, based on nonlinear oscillations in the normal star, which provides the modulation of the mass flow is discussed.
The effects of alternate-blade structural detuning and adjacent-blade alternate-circumferential-spacing aerodynamic detuning on the supersonic unstalled torsional flutter stability of a turbomachine rotor are investigated analytically. An unsteady aerodynamic model employing influence coefficients is constructed for the case of a flat-plate-airfoil cascade in torsion-mode harmonic oscillation in a supersonic inviscid isentropic adiabatic irrotational perfect-gas inlet flow with a subsonic leading-edge locus. The influence coefficients and equations of motion are derived; the model is verified by applying it to the 12-blade cascade-B flow geometry studied by Verdon and McCune (1975); and the results are presented graphically. It is found that the rotor can be stabilized over the entire reduced frequency range by applying a combination of structural and aerodynamic detuning as the passive flutter-control mechanism.
A brief review is given of non-flare investigations using data obtained by the Ultraviolet Spectrometer and Polarimeter on the Solar Maximum Mission. The major topics described are sunspot research including magnetic field measurements, oscillations, and models; mass motions in quiet and active regions including steady flows and acoustic waves; and prominence research including physical conditions, dynamics, and mass motions around prominences. Also discussed are studies of UV bursts, the formation of the Cl I line at 1351 A, ozone in the terrestrial atmosphere, and active regions using correlated observations from other instruments on the spacecraft or on the ground.
The model wing consists of a set of fiberglass panels mounted on a steel spar that spans the 8 ft test section of the UTRC Large Subsonic Wind Tunnel. The first use of this system was to measure surface pressures and flow conditions for a series of constant pitch rate ramps and sinusoidal oscillations a Mach number range, a Reynolds number range, and a pitch angle range. It is concluded that an increased pitch rate causes stall events to be delayed, strengthening of the stall vortex, increase in vortex propagation, and increase in unsteady airloads. The Mach number range causes a supersonic zone near the leading edge, stall vortex to be weaker, and a reduction of unsteady airloads.
Topics discussed include formation of a filament around a magnetic region, evolution of fine structures in a filament, the spatial distribution of prominence threads, high resolution analysis of quiescent prominences at NSO/Sacramento Peak Observatory, small-scale Doppler velocities in a quiescent prominence, Doppler velocity oscillations in quiescent prominences, oscillatory relaxation of an eruptive prominence, and matter flow velocities in an active region emission loop observed in H-alpha. Attention is also given to an automated procedure for measurement of prominence transverse velocities, the nonlinear evolution of magnetized filaments, thermal equilibrium of coronal loops and prominence formation, thermal instability in planar coronal strucutres, radiative transfer in cylindrical prominence threads, numerical simulation of a catastrophe model for prominence eruptions, and the law of evolution and destruction of solar prominences.
The goal of this research was the understanding of the various oscillatory, transient, and quasi-steady motions in sunspots and the basic structure of a sunspot. The research involved both theoretical modeling (based on thermohydrodynamic theory) and observations of dynamical phenomena in sunspots. The principal topics of the research were sunspot seismology (the interaction of solar p-modes with a sunspot as a probe of the subsurface structure of a sunspot); three minute umbral oscillations and their relation to the structure of the umbral atmosphere; siphon flows in isolated magnetic flux tubes and their relation to the photospheric Evershed flow and to intense magnetic elements outside of sunspots; and more general theoretical work on magneto-atmospheric waves. Here, a summary of results is given.
The flow of an incompressible, stratified fluid over an obstacle will produce an oscillation in which buoyancy is the restoring force, called a gravity wave. For disturbances of this scale, the atmosphere may be treated as incompressible; and even the linear approximation will explain many of the phenomena observed in the lee of mountains. However, nonlinearities arise in two ways: (1) through the large (scaled) size of the mountain, and (2) from dynamically singular levels in the fluid field. These produce a complicated array of phenomena that present hazards to aircraft and to lee surface areas. If there is no dynamic barrier, these waves can penetrate vertically into the middle atmosphere (30-100 km attitude), where recent observations show them to be of a length scale that must involve the Coriolis force in any modeling. At these altitudes, the amplitude of the waves is very large, and the waves are studied with a view to their potential impact on the projected National Aerospace Plane. This paper presents the results of analyses and state-of-the-art numerical simulations, validated where possible by observational data.
The potential of future data sets from global oscillations network group (GONG) and solar oscillations investigation (SOI) for resolving long-lived azimuthal jets and shearing flows, is investigated. Various artificial data sets are constructed, containing noise resembling that of a one-year observation run. These are inverted using a two dimensional regularized least squares inversion. The ability of this method to form well localized averages of the rotation rate, as measured by the averaging kernels, is investigated using an extensive mode set and subsets. It is shown that it is possible to keep the noise in the solution down to a few nHz in much of the solar interior, while obtaining a reasonable resolution for a GONG-like data set. At low latitudes in the middle of the convection zone, an angular resolution of less than 10 deg and a radial resolution of about 0.04 solar radii, are obtained. The averaging kernels depend on the mode set, and a reduction in the number of modes tends to introduce small-scale near surface structures into the averaging kernels which would adversely affect the inferred rotation rate.
A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.
The unsteady, three-dimensional, full Navier-Stokes (NS) equations and the Euler equations of rigid-body dynamics are sequentially solved to simulate the natural rolling response of slender delta wings of zero thickness at moderate to high angles of attack, to transonic and subsonic flows. The governing equations of fluid flow and dynamics of the present multi-disciplinary problem are solved using the time-accurate solution of the NS equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The main focus is to analyze the effect of Mach number and angle of attack on the leading edge vortices and their breakdown, the resultant rolling motion, and overall aerodynamic response of the wing. Three cases demonstrate the natural response of a 65 deg swept, cropped delta wing in a transonic flow with breakdown of the leading edge vortices and an 80 deg swept delta wing in a subsonic flow undergoing either damped or self-excited limit-cycle rolling oscillations as a function of angle of attack. Comparisons with an experimental investigation completes this study, validating the analysis and illustrating the complex details afforded by computational investigations.
NASA's Advanced Subsonic Technology (AST) program is developing new technologies to increase the fuel efficiency of commercial aircraft engines, improve the safety of engine operation, and reduce engine emissions and noise. With the development of new designs for ducted fans, compressors, and turbines to achieve these goals, a basic aeroelastic requirement is that there should be no flutter or high resonant blade stresses in the operating regime. To verify the aeroelastic soundness of these designs, we need an accurate prediction and analysis code. Such a two-dimensional viscous propulsion aeroelastic code, named TURBO-AE, is being developed at the NASA Lewis Research Center. The TURBO-AE aeroelastic code is based on a three-dimensional unsteady aerodynamic Euler/Navier-Stokes turbomachinery code TURBO, developed under a grant from NASA Lewis. TURBO-AE can model viscous flow effects that play an important role in certain aeroelastic problems, such as flutter with flow separation (or stall flutter) and flutter in the presence of shock and boundary-layer interaction. The structural dynamics representation of the blade in the TURBO-AE code is based on a normal mode representation. A finite element analysis code, such as NASTRAN, is used to calculate in-vacuum vibration modes and the associated natural frequency. A work-per-cycle approach is used to determine aeroelastic (flutter) stability. With this approach, the motion of the blade is prescribed to be a harmonic vibration in a specified in vacuum normal mode. The aerodynamic forces acting on the vibrating blade and the work done by these forces on the vibrating blade during a cycle of vibration are calculated. If positive work is being done on the blade by the aerodynamic forces, the blade is dynamically unstable, since it will extract energy from the flow, leading to an increase in the amplitude of the blade's oscillation. Initial calculations have been done for a configuration representative of the Energy Efficient Engine fan rotor. The accompanying figure shows the work-per-cycle after each cycle of vibration. It can be seen that the work-per-cycle does not vary much after the fourth cycle. The negative sign of the converged work-per-cycle shows that the fan blade is dynamically stable and will not flutter. TURBO-AE will provide a useful aeroelastic prediction/analysis capability for engine manufacturers. It will reduce design cycle times by allowing new blade designs to be verified for aeroelastic soundness before blades are built and tested. With this prediction capability, it will be possible to build thinner, lighter, and faster rotating blades without encountering aeroelastic problems like stall flutter and high-cycle fatigue due to forced vibrations.
This paper discusses Spectral Proper Orthogonal Decomposition (SPOD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements in recent NASA Ames wind tunnel test. The uPSP measurements were collected using Innovative Scientific Solutions, Inc. (ISSI) porous, fast-response pressure-sensitive paint, 40 ISSI four-inch air-cooled Light-Emitting Diodes, and 8 Phantom v2512 high-speed cameras at 10,000 frames per second in the uPSP Launch Vehicle Demonstration Test (LVDT) of the Space Launch System (SLS) vehicle in the 11-by 11-foot transonic test section of the Unitary Plan Wind Tunnel at NASA Ames Research Center in April 2024. SPOD is derived from a space-time proper orthogonal decomposition problem for statistically stationary flows. SPOD modes are determined in the frequency domain. Each SPOD mode oscillates at a single frequency. SPOD can be viewed as an extension of the Discrete Fourier Transform composition and the Dynamic Mode Decomposition. In this paper, the outputs of SPOD of the uPSP measurements in the uPSP LVDT are presented and the effectiveness of SPOD in the identification, diagnosis and analysis of the aerodynamic and aeroacoustic phenomena is demonstrated. The unsteady and dynamic property of the pressure field on the surface of the SLS Block 1B crew vehicle is presented with the visualization of the SPOD modes of the uPSP measurements in the tests of a Mach sweep run of the uPSP LVDT. The SPOD outputs were generated with the execution in parallel of a code in Python, with the library of Message Passing Interface for parallel processing, on the NASA Pleiades supercomputer. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities Portfolio Office.
Despite the fact that the low-level jet of the southern Great Plains (the GPLLJ) of the U.S. is primarily a nocturnal phenomenon that virtually vanishes during the daylight hours, it is one of the most persistent and stable climatological features of the low-level continental flow during the warm-season months, May through August. We have used significant-level data to validate the skill of the GEOS-1 Data Assimilation System (DAS) in realistically detecting this jet and inferring its structure and evolution. We have then carried out a 15-year reanalysis with the GEOS-1 DAS to determine its climatology and mean diurnal cycle and to study its interannual variability. Interannual anomalies of the meridional flow associated with the GPLLJ are much smaller than the mean diurnal fluctuations, than random intraseasonal anomalies, and than the mean wind itself. There are three maxima of low-level meridional flow variance over the Great Plains and the Gulf of Mexico: a 1.2 m2 s-2 peak over the southeast Texas, to the east and south of the mean velocity peak, a 1.0 m2 s-2 peak over the western Gulf of Mexico, and a .8 m2 s-2 peak over the upper Great Plains (UGP), near the Nebraska/South Dakota border. Each of the three variance maxima corresponds to a spatially coherent, jet-like pattern of low-level flow interannual variability. There are also three dominant modes of interannual variability corresponding to the three variance maxima, but not in a simple one-to-one relationship. Cross-sectional profiles of mean southerly wind over Texas remain relatively stable and recognizable from year to year with only its eastward flank showing significant variability. This variability, however, exhibits a distinct, biennial oscillation during the first six to seven years of the reanalysis period and only then. This intermittent biennial oscillation (IBO, one of the three modes discussed in the previous paragraph) in the lowlevel flow is restricted to the region surrounding eastern Texas and is also evident in the NCEP/NCAR reanalysis data set from about 1978 to 1985 or 1986 and again from 1995 to 2000. It is evident as well in surface pressure in both the GEOS-1 and NCEP/NCAR sets. The interannual anomalies do not necessarily persist uniformly throughout an entire season, but can fluctuate from one part of the season to the next. To estimate the characteristic sub-seasonal time scales for coherence of these fluctuations, we have taken the weekly anomaly of low-level wind at each point of the domain from the climatological average for that given point and that given week of the season and computed the covariance of its fluctuations over all weeks and over all years with the weekly climatological anomaly of the meridional wind at each of the three reference points discussed above. The typical duration of a coherent interannual anomaly within a given warm season increases with decreasing latitude from 2 to 3 weeks over the UGP, to 6 to 7 weeks over eastern Texas. Coherence over the western Gulf of Mexico is intermediate between the two with a typical duration of 4 to 5 weeks. There appears to be evidence that the interannual anomalies over Texas the Gulf propagate to the UGP after a week and those over the Gulf propagate there after 2 to 3 weeks. There also appears to be some reverse propagation of interannual anomalies over the UGP to Texas and to the Gulf after a period of about one week. The interannual anomalies in southerly flow over eastern Texas seem to correlate well with interannual anomalies of surface temperature and (negative) ground wetness and over western Texas.
This study was a theoretical and experimental investigation of thin liquid sheet flows in vacuum. A sheet flow created by a narrow slit of width, W, coalesces to a point at a distance, L, as a result of surface tension forces acting at the sheet edges. As the flow coalesces, the fluid accumulates in the sheet edges. The observed triangular shape of the sheet agrees with the calculated triangular result. Experimental results for L/W as a function of Weber number, We, agree with the calculated result, L/W = the sq. root of 8We. The edge cross sectional shape is found to oscillate from elliptic to 'cigar' like to 'peanut' like and then back to elliptic in the flow direction. A theoretical one-dimensional model was developed that yielded only elliptic solutions for the edge cross section. At the points where the elliptic shapes occur, there is agreement between theory and experiment.
This paper considers transonic flutter mechanisms of the Benchmark Supercritical Wing, a model under study in the Aeroelastic Prediction Workshop series. Flutter boundaries are mapped out across an angle of attack sweep at Mach 0.8, utilizing both time-domain and linearized frequency-domain solvers, manual meshes and adapted meshes, and various governing equations. With increased angle of attack, linearized and finite amplitude flutter predictions exhibit differences above 3◦ as the flow begins to separate; the latter predictions are found to be driven by subcritical limit cycle oscillations whose strength increases with angle of attack. Moderate perturbation values provide a stability boundary at 5◦ which matches the experimental data, but it is not clear how the experimental perturbation, from one test condition to the next, can be reasonably characterized.