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.
Nonlinear combustion instability in liquid propellant rocket motors, and research on damping devices, reaction kinetics, and cold flow tests
Complexity of the phenomena of panel flutter instability has resulted in the necessity of developing separate design criteria for a variety of flow conditions and panel configurations. Vehicle panel configurations with low aspect ratios are of interest in low supersonic flow, where boundary layer effects are important.
This paper describes a generalized dynamic model which has been developed for use in compression component aerodynamic stability studies. The model is a one-dimensional, pitch-line, blade row, lumped volume system. Arbitrary placement of blade free volumes upstream, within, and downstream of the compression component as well as the removal of bleed flow from the exit of any rotor or stator are model options. The model has been applied to a two-stage fan and an eight-stage compressor. The clean inlet pressure ratio/flow maps and the surge line have been reproduced, thereby validating the capability of the dynamic model to reproduce the steady-flow characteristics of the compression component. A method for determining the onset of an aerodynamic instability which is associated with surge is described. Sinusoidally time-varying inlet and exit boundary conditions have been applied to the eight stage compressor as examples of the manner in which this model may be used for stability studies.
To investigate the cause of flares and quiescent filament eruptions the quasi-static evolution of a magnetohydrostatic (MHS) model was studied. The results lead to a proposal that: the sudden disruption of an active-region filament field configuration and the accompanying flare result from the lack of a neighboring equilibrium state as magnetic shear is increased above the critical value; and a quiescent filament eruption is due to an ideal MHD kink instability of a highly twisted detached flux tube formed by the increase of plasma current flowing along the length of the filament. A numerical solution was developed for the 2-D MHS equation for the self-consistent equilibrium of a filament and overlying coronal magnetic field. Increase of the poloidal current causes increase of magnetic shear. As shear increases past a critical point, there is a discontinuous topological change in the equilibrium configuration. It was proposed that the lack of a neighboring equilibrium triggers a flare. Increase of the axial current results in a detached tube with enough helical twist to be unstable to ideal MHD kink modes. It was proposed that this is the condition for the eruption of a quiescent filament.
The three-dimensional time-dependent, compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation algorithm to study the stability of a Mach 4.5 flow over a flat plate. Several nonlinear direct simulations suggest the existence of a secondary instability which might provide a possible route to transition. Pertinent differences in the energy content of the various Fourier modes between this instability and the more common incompressible K-type instabilities are pointed out.
The spectral evolution of accretion disks in X-ray binaries containing black holes is studied, based on the disk instability model. The thermal transition of the outer portions of the disk controls the mass flow rate into the inner portions of the disk, thus modulating the soft X-ray flux which is thought to arise from the inner disk. Calculated soft X-ray spectra are consistent with the observations of the X-ray transient A0620 - 00 and especially ASM 2000 + 25, the soft X-ray spectra of which are well fitted by blackbody radiation with a fixed inner edge of the disk, Rin, and with monotonically decreasing temperature at Rin with time. Since the gas pressure is always dominant over the radiation pressure during the decay in these models, a two-temperature region is difficult to create. Instead, it is suggested that hard X-rays are generated in a hot (kT greater than 10 keV) accretion disk corona above the cool (kT less than 1 keV) disk.
The interaction of current sheets embedded in the upstream flow with a shock, relevant to the study of the earth's bow shock, is examined. It is shown that a hot flow anomaly (HFA) can be generated by a direct method that does not involve an instability. The HFA is shown to be due to the interaction of reflected ions with the current sheets. An important part of the interaction, pointed out by Burgess (1989), is that for a class of current sheets reflected ions are always focused toward the current sheet by the motional electric field, while for another class of current sheets the electric field defocuses reflected ions away from the current sheet. In addition, not just the behavior of reflected ions upstream of the shock but also the behavior of those behind the shock front is related to HFA formation.
The present numerical simulation explores a thermal-convective mechanism for oscillatory thermocapillary convection in a shallow Cartesian cavity for a Prandtl number 6.78 fluid. The computer program developed for this simulation integrates the two-dimensional, time-dependent Navier-Stokes equations and the energy equation by a time-accurate method on a stretched, staggered mesh. Flat free surfaces are assumed. The instability is shown to depend upon temporal coupling between large scale thermal structures within the flow field and the temperature sensitive free surface. A primary result of this study is the development of a stability diagram presenting the critical Marangoni number separating steady from the time-dependent flow states as a function of aspect ratio for the range of values between 2.3 and 3.8. Within this range, a minimum critical aspect ratio near 2.3 and a minimum critical Marangoni number near 20,000 are predicted below which steady convection is found.
Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz instability for a magnetopauselike configuration. The unidirectional magnetic field is essentially transverse to the plasma flow velocity, which is the most unstable case according to linear theory and models the flow dynamics in the subsolar region of the magnetopause for northward interplanetary magnetic field. We recover effects analogous to those found in MHD simulations, including a mode cascade to longer wavelengths. The boundary layer consists of coherent structure and is not well described by a diffusive process. Isolated structures on the order of the ion gyroradius are formed which can cross the boundary in either direction. We describe how the time evolution of these structures represents transport across boundary layers, and we consider the possible connection of these entities to flux transfer events and other structure seen in the low-latitude boundary layer at the Earth's magnetopause as well as to flux ropes commonly observed near the ionopause of Venus. We also discuss the relation of the hybrid calculations to previous MHD simulations and to observations.
Convective fluid dynamics of immiscible double and triple liquid layers are considered. First results on multilayer convective flow, in preparation for spaceflight experiment aboard IML-2 (International Microgravity Laboratory), are discussed. Convective flow in liquid layers with one or two horizontal interfaces with heat flow applied parallel to them is one of the systems investigated. The second system comprises two horizontally layered immiscible liquids heated from below and cooled from above, that is, heat flow orthogonal to the interface. In this system convection results due to the classical Benard instability.
The instability of an incompressible three-dimensional boundary layer (that is, one with cross-flow) is considered theoretically and computationally in the context of vortex/wave interactions. Specifically the work centers on two low amplitude, lower-branch Tollmien-Schlichting waves which mutually interact to induce a weak longitudinal vortex flow; the vortex motion, in turn, gives rise to significant wave-modulation via wall-shear forcing. The characteristic Reynolds number is taken as a large parameter and, as a consequence, the waves' and the vortex motion are governed primarily by triple-deck theory. The nonlinear interaction is captured by a viscous partial-differential system for the vortex coupled with a pair of amplitude equations for each wave pressure. Three distinct possibilities were found to emerge for the nonlinear behavior of the flow solution downstream - an algebraic finite-distance singularity, far downstream saturation or far-downstream wave-decay (leaving pure vortex flow) - depending on the input conditions, the wave angles, and the size of the cross-flow.
The influence of time-dependent heat addition on the linear instablity of shear layers is of considerable interest in understanding the dynamic behavior of reacting flows and combustion-turbulence interactions. The approach is based upon the Bernoulli enthalpy aeroacoustics theory, which utilizes the specific enthalpy and specific entropy as the primary thermodynamic variables. In addition, velocity oscillations are split into Helmoholtz decomposition theorem.
An experimental study was conducted to determine the performance of a two-dimensional, mixed-compression bifurcated duct inlet system designed for a free-stream Mach number of 2.7. Thirty percent of the supersonic area contraction occurred internally. A movable ramp was used to vary the contraction ratio for off-design operation. Boundary layer bleed regions were located on the cowl, centerbody, and sidewall surfaces. There were also provisions for vortex generators on the cowl and centerbody of the subsonic diffuser. Data were obtained over the Mach number range of 2.0 to 2.8 and at angles of yaw from 0 deg. to the maximum value prior to inlet un-start. The test at Mach 2.8 was to obtain data for an over- speed condition. The Reynolds number varied from 2.5 to 2.3 million/ft for Mach numbers above 2.5. At Mach numbers of 2.5 and lower, the Reynolds number was set at 2.5 million/ft. Bleed patterns, vortex generator patterns, and ramp position were varied, and three inlet configurations were selected for more extensive study. Two of these configurations had self-starting capability. The self-starting configuration that was developed produced 89 percent total pressure recovery at the compressor face station with 6.8 percent total bleed. The compressor face distortion was about 16 percent. Vortex generators were extremely effective in re-distributing flow but were not as effective in reducing distortion. Excellent flow symmetry was achieved between the separated halves of the inlet, and twin-duct instability was not observed. The ramp tip shock was steeper than expected. This caused the cowl lip shock to be reflected from the ramp instead of being cancelled at the shoulder. However, peak recovery at the throat was still obtained with the ramp near the design position.
Io has very high surface heat flow and an abundance of volcanic activity, which are thought to be driven by nonuniform tidal heating in its interior. This nonuniform heat is transported to the base of the lithosphere by very vigorous convection in Io's silicate mantle, the form of which is presumably responsible for the distribution of surface features such as volcanoes and mountains. We here present three-dimensional spherical calculations of mantle convection in Io, in order to ascertain the likely form of this convection and the resulting distribution of heat flow at the surface and core-mantle boundary. Different models of tidal dissipation are considered: the endmember scenarios (identified by M. N. Ross and G. Schubert) of dissipation in the entire mantle, or dissipation in a thin (approximately 100-km-thick) asthenosphere, as well as the 'preferred' distribution of M. N. Ross et al. comprising 1/3 mantle and 2/3 asthenosphere heating. The thermal structure of Io's mantle and asthenosphere is found to be strongly dependent on tidal heating mode, as well as whether the mantle-asthenosphere boundary is permeable or impermeable. Results indicate a large-scale flow pattern dominated by the distribution of tidal heating, with superimposed small-scale asthenospheric instabilities that become more pronounced with increasing Rayleigh number. These small-scale instabilities spread out the surface heat flux, resulting in smaller heat flux variations with increasing Rayleigh number. Scaled to Io's Rayleigh number of O(10(exp 12)) variations of order a few percent are expected. This small but significant variation in surface heat flux may be compatible with the observed distributions of volcanic centers and mountains, which appear fairly uniform at first sight but display a discernible distribution when suitably processed. The observed distribution of volcanic centers is similar to the asthenosphere heating distribution, implying that most of the tidal heating in Io occurs in an asthenosphere.
The basic kinetic and thermodynamic concepts of crystal growth will be revisited in view of recent AFM and interferometric findings. These concepts are as follows: 1) The Kossel crystal model that allows only one kink type on the crystal surface. The modern theory is developed overwhelmingly for the Kessel model; 2) Presumption that intensive step fluctuations maintain kink density sufficiently high to allow applicability of Gibbs-Thomson law; 3) Common experience that unlimited step bunching (morphological instability) during layer growth from solutions and supercooled melts always takes place if the step flow direction coincides with that of the fluid.
Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.
The effects of forward- and backward-facing steps on the receptivity and stability of three-dimensional supersonic boundary layers over a swept wing with a blunt leading edge are numerically investigated for a freestream Mach number of 3 and a sweep angle of 30 degrees. The flow fields are obtained by solving the full Navier-Stokes equations. The evolution of instability waves generated by surface roughness is simulated with and without the forward- and backward-facing steps. The separation bubble lengths are about 5-10 step heights for the forward-facing step and are about 10 for the backward-facing step. The linear stability calculations show very strong instability in the separated region with a large frequency domain. The simulation results show that the presence of backward-facing steps decreases the amplitude of the stationary crossflow vortices with longer spanwise wavelengths by about fifty percent and the presence of forward-facing steps does not modify the amplitudes noticeably across the steps. The waves with the shorter wavelengths grow substantially downstream of the step in agreement with the linear stability prediction.
We develop an iterative automated method to predict transition locations in boundary-layer flows by using the FUN3D solver to perform flow simulations and the LASTRAC code for linear stability computations. The coupling of FUN3D and LASTRAC allows for a robust physics-based approach to model boundary-layer transition by analyzing the growth of different instability waves and then using that information to iteratively update the resulting transition location. There is no user involvement during the iterative computations. We apply this automated method to subsonic flow over a flat plate with a sharp leading edge. The final solution has regions of laminar and turbulent flow with a transition onset location that agrees with experiments and stability-based correlations. This iterative automated method is also applied to an NLF(1)-0416 airfoil for conditions with and without a separation bubble. Along with predicting transition locations, we compare the streamwise distributions of surface-pressure and skin-friction coefficients to a transport-equation-based model. We consider a 6:1 prolate spheroid at three angles of attack, namely, five, ten, and fifteen degrees, where mixed-mode transition occurs due to both Tollmien-Schlichting and crossflow instabilities. The skin-friction contours and transition fronts at every angle of attack from our iterative automated method show good agreement with past experimental and computational results in the literature for the 6:1 prolate spheroid.