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.
Computer program has been developed as analytical tool to predict severity of effects of exhaust of rocket engines on adjacent spacecraft surfaces. Program computes forces, moments, pressures, and heating rates on surfaces immersed in or subjected to exhaust plume environments. Predictions will be useful in design of systems where such problems are anticipated.
An analysis of experimental nozzle, exhaust plume, and exhaust plume impingement data is presented. The data were obtained for subscale solid propellant motors with propellant Al loadings of 2, 10 and 15% exhausting to simulated altitudes of 50,000, 100,000 and 112,000 ft. Analytical predictions were made using a fully coupled two-phase method of characteristics numerical solution and a technique for defining thermal and pressure environments experienced by bodies immersed in two-phase exhaust plumes.
Human Factors Engineering support was provided for the 30% design review of the late Space Station Freedom Payload Control Area (PCA). The PCA was to be the payload operations control room, analogous to the Spacelab Payload Operations Control Center (POCC). This effort began with a systematic collection and refinement of the relevant requirements driving the spatial layout of the consoles and PCA. This information was used as input for specialized human factors analytical tools and techniques in the design and design analysis activities. Design concepts and configuration options were developed and reviewed using sketches, 2-D Computer-Aided Design (CAD) drawings, and immersive Virtual Reality (VR) mockups.
A series of computer codes has been developed to predict gas-particle flows and resulting impingement forces, moments and heating rates to surfaces immersed in the flow. The gas-particle flow solution is coupled via heat transfer and drag between the phases with chemical effects included in the gas phase. The flow solution and impingement calculations are discussed. Analytical results are compared with test data obtained to evaluate gas-particle effects on the Space Shuttle thermal protection system during the staging maneuver.
The fractionation of fine-grained aggregates into their major components is a problem in many scientific areas including earth and planetary science. Electrophoresis, the transport of electrically charged particles, immersed in a suspension medium, by a direct current field (Bier, 1959), was employed in this study as a means of separating simulated lunar soil into its constituent minerals. In these tests, conducted in a static analytical cylindrical microelectrophoresis apparatus, samples of simulated lunar soil and samples of pure mineral constituents were placed in the chamber; the electrophoretic mobilities of the lunar soil and the individual mineral constituents were measured. In most of the suspension buffers employed separability was indicated, on the basis of differences in mobility, for all the constituent mineral species except ilmenite and pyroxene, which were not efficiently separable in any of the buffers. Although only a few suspension media were employed, the success of this initial study suggests that electrophoresis may be an important mineral fractionation option in fine-grained aggregate processing.
An analytical model is fitted to experimental data to show that rod-to-coil (RTC) transformations in single strand polymers (SSPs) are close to true phase transitions. The RTC are known to occur when the temperature or pH of a solvent in which the SSP is immersed is varied abruptly. Long wavelength degrees of freedom in the SSPs induce long range interactions in the molecule which lead to a first order phase transition. The Landau formalism suggests that a long-range ordered state can appear at a low temperature. A microscopic model of the phonon degrees of freedom identifies the nonzero temperature threshold at which the phase transition will occur. The phenomenology described is similar to a previous model for transitions in double-strand polymers.
The capacitance of a grain immersed in a steady state plasma containing a size distribution of dust particles is studied. The grain charge is determined by assuming the equilibrium potential has been obtained by a simple balance of electron and ion collection currents. It is shown that the validity of the analytical treatment given here for the linearized Poisson equation is confined to a certain region of space. Within this region and starting at very small plasma Debye length lambda(D), the capacitance at first exhibits a monotonic increase with increasing lambda(D). The capacitance eventually reaches a maximum, followed by a monotonic decrease. The charge density of the dust in the plasma is found to be only a function of the lambda(D); there is no significant dependence on the interparticle spacing.
Infrared imaging systems can be used to measure temperatures of actively heated bodies immersed in an airstream. This monitoring of the convective heat transfer process, provides also information about the interaction between the body and the flow. The concept appeals to Nusselt/Reynolds numbers relations in order to produce data of interest from surface temperatures. Two test cases are presented and reference is made to analytical results: the mapping of a laminar jet and the temperature distribution along a constant power heated flat plate in laminar boundary layer regime. Although this research is currently focused on low speed aerodynamics, the extension to high speed aerodynamics, where the body undergoes frictional heating is of interest in this context, too.
The steady state interaction of two electrically biased parallel plates immersed in a flowing plasma characteristic of low earth orbit is studied numerically. Fluid equations are developed to describe the motion of the cold positively charged plasma ions, and are solved using finite-differences in two dimensions on a Cartesian grid. The behavior of the plasma electrons is assumed to be described by the Maxwell-Boltzmann distribution. Results are compared to an analytical and a particle simulation technique for a simplified flow geometry consisting of a single semi-infinite negatively biased plate. Comparison of the extent of the electrical disturbance into the flowing plasma and the magnitude of the current collected by the plate is very good. The interaction of two equally biased parallel plates is studied as a function of applied potential. The separation distance at which the current collected by either plate decreases by five and twenty percent is determined as a function of applied potential. The percent decreases were based on a non-interacting case. The decrease in overall current is caused by a decrease in ionic density in the region between the plates. As the separation between the plates decreases, the plates collect the ions at a faster rate than they are supplied to the middle region by the oncoming plasma flow. The docking of spacecraft in orbit is simulated by moving two plates of unequal potential toward one another in a quasi-static manner. One plate is held at a large negative potential while the other floats electrically in the resulting potential field. It is found that the floating plate does not charge continuously negative as it approaches the other more negatively biased plate. Instead, it charges more and then less negative as ionic current decreases and then increases respectively upon approach. When the two plates come into contact, it is expected that the electrically floating plate will charge rapidly negative to a potential near that of the other plate.
The results of four current analytical methods for predicting wing vortex strength and decay rate are compared with the results of a flight investigation of the wake characteristics of several large jet transport aircraft. An empirical expression defining the strength and decay rate of wake vortices is developed that best represents most of the flight-test data. However, the expression is not applicable to small aircraft that would be immersed in the vortex wake of large aircraft.
We examine the rate of accretion of mass and spin angular momentum by a spherical solid body ('planet') immersed in a differentially rotating disk of particles ('planetesimals'). If the planet travels on a circular orbit, the accretion process is described by two dimensionless parameters r and s, which measure the ratio of the planet's radius to its Hill radius and the ratio of the rms radial velocity of the planetesimals to the shear across the Hill radius, respectively. Using a combination of analytic arguments and numerical simulations, we derive the mass and angular momentum accretion rates and their scaling with r and s. By introducing an additional parameter, the effective mass of the planetesimals relative to the planet, we can derive the obliquities and spins arising from the stochastic nature of the accretion process. Our results are consistent with those of previous calculations by Lissauer and Kary (1991) wherever there is overlap in parameter space. In particular, we conclude that ordered accretion from a uniform or slowly varying disk of small bodies cannot result in rotation as rapid as that of earth or Mars for any value of the disk velocity dispersion s. The spin rates of these planets are most naturally explained as arising from one or a few 'giant' impacts by planetesimals.
Eleven samples from five partially degassed Mesozoic granitic xenoliths and two late Pleistocene host basalt samples from Sawmill Canyon in Sierra Nevada, California, have been analyzed using stepwise Ar-40/Ar-39 dating techniques. The petrography and analytic techniques are described and the results are presented. Isotopic compositions of Ar released from the xenoliths in several extraction steps at temperatures below about 900 C were colinear in Ar-36/Ar-40 versus Ar-39/Ar-40 diagrams and defined isochrons, giving a mean age of degassing of 119,000 + or - 7000 yr. Ar-40 extracted at higher temperatures included ancient radiogenic Ar-40 that never diffused from the xenoliths during immersion in the magma. This Ar-40 caused an increase in the apparent age for the high-temperature extractions. The high precision of the eruption age determined by this method is comparable to that obtained elsewhere by conventional K/Ar dating of sanidine.
Two- and three-lobed oscillations of a rotating liquid drop immersed in an immiscible fluid of comparable density and the same angular velocity were studied experimentally. Using acoustically suspended drops, it has been found that the relative change in the resonance frequencies of the axisymmetric drop-shape oscillations Delta omega-l/omega-l(0) is proportional to the square of the normalized angular speed (Omega/omega-l/0)-squared when omega-l greater than 2 Omega. This is in agreement with a recent analytical study of the same problem. Some preliminary results regarding the effect of rotation on the free-decay rate of the two-lobed oscillations are also presented.
As next-generation aircraft and vehicles continue to develop, so do their associated energy demands. Lithium metal batteries are a leading candidate to fulfill this energy requirement, but these batteries are prone to internal dendrite defects that can lead to catastrophic thermal runaway events. Current battery management systems are capable of mitigating such risks, but are unable to detect such defects until thermal runaway has already begun. Various nondestructive evaluation (NDE) techniques, particularly ultrasonic NDE, can directly monitor internal battery parameters giving them the potential to detect critical defects prior to catastrophic failure. However, most of the current battery NDE research has focused on improved battery state-of-charge (SOC) and stateof- health (SOH) monitoring with little emphasis on critical defect detection. Thus, a measurement technique sensitive to subtle battery defects is needed. In addition, the complex mechanics of ultrasound in porous, thin, multilayered batteries prompt the use of physics-based simulation to guide inspections. In this work, an ultrasonic NDE technique has been developed utilizing frequency domain analysis of local battery resonances to detect the presence of battery defects. This technique is a practical extension of local ultrasonic resonance spectroscopy (LURS) – which previously required non-contact laser ultrasonics – to measurements with piezoelectric contact and immersion scan transducers. To extend the technique to work with piezoelectric transducers, ultrasonic battery measurements were compared to a sans-battery calibration measurement. Then, a linear systems deconvolution was used to eliminate the transfer functions of extraneous factors such as the transducer and electronics, leaving only the frequency-dependent battery reflection coefficient. The LURS technique was first validated on stainless steel and aluminum plates, producing reflection coefficients in line with analytical and numerical finite element modeling (FEM) results. Functioning Li-metal pouch cells were then seeded with lithium chip defects prior to LURS measurements. The presence of these defects is shown to cause a measurable shift in the battery’s through-thickness local resonances. 2D, frequency-domain poroelastic models of ultrasonic propagation in a single-cell lithium metal pouch battery were created and corroborated these findings. Thus, this work has both extended and proven the feasibility of the LURS technique in the detection of local battery defects.
The hydrodynamic force on a cylinder oscillating transversely to its axis is a nonlinear function of the displacement amplitude x0. We report measurements and numerical calculations of the force at frequencies low enough that delta > R, where delta is the viscous penetration length and R is the cylinder radius. For small amplitudes, the numerically calculated Fourier transform of the force per unit length, F(sub small), agrees with Stokes' analytical calculation. For larger amplitudes, the force per unit length found by both calculation and measurement is F = F(sub small)C (x(sub 0)/delta,R/delta). The complex function C depends only weakly on R/delta, indicating that x0/delta is more appropriate as a scaling variable than the Keulegan-Carpenter number KC = pi*x(sub 0)/R. The measurements used a torsion oscillator driven at frequencies from 1 to 12 Hz while immersed in dense xenon. The oscillator comprised cylinders with an effective radius of R = 13.4 micron and oscillation amplitudes as large as x(sub 0)/delta = 4 (corresponding to KC as large as 71). The calculations used similar conditions except that the amplitudes were as large as x0/delta = 28.
The stability of a self-gravitating, nonrotating, isothermal gas layer threaded by a one-dimensional equipartition magnetic field, immersed in a rigid isothermal layer of stars, is considered with respect to waves with motions in and perpendicular to the X plane, where X is the equilibrium magnetic and gravitational field vectors. When motions are perpendicular to the X plane, the magnetic field hinders gravitational instability, increasing the minimum length necessary to produce instability by the factor (1 + alpha) to the 1/2 power, where alpha is the ratio of magnetic pressure to gas pressure. When motions are in the plane, no such simple analytical solution is found. However, the resulting system of equations reduces to a single fourth order differential equation in the perturbed gas potential that defines an eigenvalue problem for the marginally unstable state when the four appropriate boundary conditions are considered.
We analyze plasma, magnetic field, and electric field data for a flux transfer event (FTE) to highlight improvements in our understanding of these transient reconnection signatures resulting from high-resolution data. The approximate 20 s long, reverse FTE, which occurred south of the geomagnetic equator near dusk, was immersed in super-Alfvnic flow. The field line twist is illustrated by the behavior of flows parallel perpendicular to the magnetic field. Four-spacecraft timing and energetic particle pitch angle anisotropies indicate a flux rope (FR) connected to the Northern Hemisphere and moving southeast. The flow forces evidently overcame the magnetic tension. The high-speed flows inside the FR were different from those outside. The external flows were perpendicular to the field as expected for draping of the external field around the FR. Modeling the FR analytically, we adopt a non-force free approach since the current perpendicular to the field is nonzero. It reproduces many features of the observations.
The Fast Scattering Code (FSC) is a frequency domain noise prediction program developed at the NASA Langley Research Center (LaRC) to simulate the acoustic field produced by the interaction of known, time harmonic incident sound with bodies of arbitrary shape and surface impedance immersed in a potential flow. The code uses the equivalent source method (ESM) to solve an exterior 3-D Helmholtz boundary value problem (BVP) by expanding the scattered acoustic pressure field into a series of point sources distributed on a fictitious surface placed inside the actual scatterer. This work provides additional code validation studies and illustrates the range of code parameters that produce accurate results with minimal computational costs. Systematic noise prediction studies are presented in which monopole generated incident sound is scattered by simple geometric shapes - spheres (acoustically hard and soft surfaces), oblate spheroids, flat disk, and flat plates with various edge topologies. Comparisons between FSC simulations and analytical results and experimental data are presented.