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.
Potential sweep method for organic and adsorption analyses, showing relation between galvanostatic and fast potential sweep transients, noting mathematical analysis of current-potential transients
A preliminary design tool for advanced propellers was developed combining a fast vortex lattice aerodynamic analysis, a fast subsonic point source noise analysis, and an optimization scheme using a conjugate directions method. Twist, chord and sweep distributions are optimized to simultaneously improve both the aerodynamic performance and the noise observed at a fixed relative position. The optimal noise/performance tradeoffs for straight and advanced concept blades are presented. The techniques used include increasing the blade number, blade sweep, reducing the rotational speed, shifting the spanwise loading and diameter changes.
A preliminary design tool for advanced propellers was developed combining a fast vortex lattice aerodynamic analysis, a fast subsonic point source noise analysis, and an optimization scheme using a conjugate directions method. Twist, chord and sweep distributions are optimized to simultaneously improve both the aerodynamic performance and the noise observed at a fixed relative position. The optimal noise/performance tradeoffs for straight and advanced concept blades are presented. The techniques used include increasing the blade number, blade sweep, reducing the rotational speed, shifting the spanwise loading and diameter changes.
Various methods are described for measuring heat generation in primary and secondary batteries as well as the specific heat of batteries and cell thermal conductance. Problems associated with determining heat generation in large batteries are examined. Special attention is given to monitoring temperature gradients in nickel cadmium cells, the use of auxiliary electrodes for conducting tests on battery charge control, evaluating the linear sweep of current from charge to discharge, and determining zero current voltage. The fast transient behavior of batteries in the microsecond range, and the electrical conductance of nickel sinters in the thickness direction are also considered. Mechanical problems experienced in the vibration of Ni-Cd batteries and tests to simulate cyclic fatigue of the steel table connecting the plates to the comb are considered. Methods of defining the distribution of forces when cells are compressed during battery packaging are also explored.
A general numerical strategy for solving partial differential equations and other functional problems by cycling between coarser and finer levels of discretization is described. Optimal discretization schemes are provided together with very fast general solvers. It is described in terms of finite element discretizations of general nonlinear minimization problems. The basic processes (relaxation sweeps, fine-grid-to-coarse-grid transfers of residuals, coarse-to-fine interpolations of corrections) are directly and naturally determined by the objective functional and the sequence of approximation spaces. The natural processes, however, are not always optimal. Concrete examples are given and some new techniques are reviewed. Including the local truncation extrapolation and a multilevel procedure for inexpensively solving chains of many boundary value problems, such as those arising in the solution of time-dependent problems.
Results are presented of a study conducted to analyze the potential aerodynamic improvements which can be obtained with aft-swept wing tips at subsonic speeds. A review and summary of previous swept-tip applications and studies are given. Additionally, the through evolution optimized tail-fin shapes of fast-cruising aquatic animals are examined as well as the wing planform shapes of sea birds. These superior long-distance travellers both display lifting-surface planform shapes with increased leading-edge sweep towards the tip. These observations and analytical results obtained with a nonlinear surface-panel method demonstrate that substantial induced-drag improvements can be obtained through the adoption of planform shapes with highly-tapered aft-swept tips.
Several techniques are discussed which can be used to explore laboratory photochemical processes and kinetics relevant to planetary atmospheres; these include time-resolved laser gain-versus-absorption spectroscopy and time-resolved Fourier transform infrared (FTIR) emission studies. The laser gain-versus-absorption method employed tunable diode and F-center lasers to determine the yields of excited photofragments and their kinetics. The time-resolved FTIR technique synchronizes the sweep of a commercial FTIR with a pulsed source of light to obtain emission spectra of novel transient species in the infrared. These methods are presently being employed to investigate molecular photodissociation, the yields of excited states of fragments, their subsequent reaction kinetics, Doppler velocity distributions, and velocity-changing collisions of translationally fast atoms. Such techniques may be employed in future investigations of planetary atmospheres, for example to study polycyclic aromatic hydrocarbons related to cometary emissions, to analyze acetylene decomposition products and reactions, and to determine spectral features in the near infrared and infrared wavelength regions for planetary molecules and clusters.
Numerical solutions of the Navier-Stokes equations using explicit schemes can be obtained at the expense of efficiency. Conventional implicit methods which often achieve fast convergence rates suffer high cost per iteration. A new implicit scheme based on lower-upper factorization and symmetric Gauss-Seidel relaxation offers very low cost per iteration as well as fast convergence. High efficiency is achieved by accomplishing the complete vectorizability of the algorithm on oblique planes of sweep in three dimensions.
Mass spectra are produced in most mass spectrometers by sweeping some parameter within the instrument as the sampled gases flow into the ion source. It is evident that any fluctuation in the gas during the sweep (mass scan) of the instrument causes the output spectrum to be skewed in its mass peak intensities. The time of flight mass spectrometer (TOFMS) with its fast, repetitive mode of operation produces spectra without skewing or varying instrument parameters and because all ion species are ejected from the ion source simultaneously, the spectra are inherently not skewed despite rapidly changing gas pressure or composition in the source. Methods of exploiting this feature by utilizing fast digital data acquisition systems, such as transient recorders and signal averagers which are commercially available are described. Applications of this technique are presented including TOFMS sampling of vapors produced by both pulsed and continuous laser heating of materials.
A direct simulation of turbulent flow in a channel is analyzed by the method of empirical eigenfunctions (Karhunen-Loeve procedure, proper orthogonal decomposition). This analysis reveals the presence of propagating plane waves in the turbulent flow. The velocity of propagation is determined by the flow velocity at the location of maximal Reynolds stress. The analysis further suggests that the interaction of these waves appears to be essential to the local production of turbulence via bursting or sweeping events in the turbulent boundary layer, with the additional suggestion that the fast acting plane waves act as triggers.
Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.
A variant of multigrid, based on zebra relaxation, and a new family of restriction/prolongation operators is described. Using zebra relaxation in combination with an operator-induced prolongation leads to fast convergence, since the coarse grid can correct all error components. The resulting algorithms are not only fast, but are also robust, in the sense that the convergence rate is insensitive to the mesh aspect ratio. This is true even though line relaxation is performed in only one direction. Multigrid becomes a direct method if an operator-induced prolongation is used, together with the induced coarse grid operators. Unfortunately, this approach leads to stencils which double in size on each coarser grid. The use of an implicit three point restriction can be used to factor these large stencils, in order to retain the usual five or nine point stencils, while still achieving fast convergence. This algorithm achieves a V-cycle convergence rate of 0.03 on Poisson's equation, using 1.5 zebra sweeps per level, while the convergence rate improves to 0.003 if optimal nine point stencils are used. Numerical results for two and three dimensional model problems are presented, together with a two level analysis explaining these results.
A variant of multigrid, based on zebra relaxation, and a new family of restriction/prolongation operators is described. Using zebra relaxation in combination with an operator-induced prolongation leads to fast convergence, since the coarse grid can correct all error components. The resulting algorithms are not only fast, but are also robust, in the sense that the convergence rate is insensitive to the mesh aspect ratio. This is true even though line relaxation is performed in only one direction. Multigrid becomes a direct method if an operator-induced prolongation is used, together with the induced coarse grid operators. Unfortunately, this approach leads to stencils which double in size on each coarser grid. The use of an implicit three point restriction can be used to factor these large stencils, in order to retain the usual five or nine point stencils, while still achieving fast convergence. This algorithm achieves a V-cycle convergence rate of 0.03 on Poisson's equation, using 1.5 zebra sweeps per level, while the convergence rate improves to 0.003 if optimal nine point stencils are used. Numerical results for two- and three-dimensional model problems are presented, together with a two level analysis explaining these results.
A new, fast, non-iterative version of the "Wall Pressure Signature Method" is described and used to determine blockage and angle-of-attack wind tunnel corrections for highly-powered jet-flap models. The correction method is complemented by the application of tangential blowing at the tunnel floor to suppress flow breakdown there, using feedback from measured floor pressures. This tangential blowing technique was substantiated by subsequent flow investigations using an LV. The basic tests on an unswept, knee-blown, jet flapped wing were supplemented to include the effects of slat-removal, sweep and the addition of unflapped tips. C sub mu values were varied from 0 to 10 free-air C sub l's in excess of 18 were measured in some cases. Application of the new methods yielded corrected data which agreed with corresponding large tunnel "free air" resuls to within the limits of experimental accuracy in almost all cases. A program listing is provided, with sample cases.
A major modification of and addition to existing Closed Brayton Cycle (CBC) space power system optimization codes was completed. These modifications relate to the global minimum mass search driver programs containing three nested iteration loops comprising iterations on cycle temperature ratio, and three separate pressure ratio iteration loops--one loop for maximizing thermodynamic efficiency, one for minimizing radiator area, and a final loop for minimizing overall power system mass. Using the method of steepest ascent, the code sweeps through the pressure ratio space repeatedly, each time with smaller iteration step sizes, so that the three optimum pressure ratios can be obtained to any desired accuracy for each of the objective functions referred to above (i.e., maximum thermodynamic efficiency, minimum radiator area, and minimum system mass). Two separate options for the power system heat source are available: 1. A nuclear fission reactor can be used. It is provided with a radiation shield 1. (composed of a lithium hydride (LiH) neutron shield and tungsten (W) gamma shield). Suboptions can be used to select the type of reactor (i.e., fast spectrum liquid metal cooled or epithermal high-temperature gas reactor (HTGR)). 2. A solar heat source can be used. This option includes a parabolic concentrator and heat receiver for raising the temperature of the recirculating working fluid. A useful feature of the code modifications is that key cycle parameters are displayed, including the overall system specific mass in kilograms per kilowatt and the system specific power in watts per kilogram, as the results for each temperature ratio are computed. As the minimum mass temperature ratio is encountered, a message is printed out. Several levels of detailed information on cycle state points, subsystem mass results, and radiator temperature profiles are stored for this temperature ratio condition and can be displayed or printed by users.
A production rate of 0.065 + or - 0.003 Ar-37 atom/kg min of K-39 at 2-mwe depth below sea level was measured by sweeping argon from potassium solutions. This rate is unaffected by surrounding the solution by paraffin and is attributed to negative muon captures and the electromagnetic interaction of fast muons, and not to nucleonic cosmic ray component. The Ar-37 yield from K-39 by the stopping of negative muons in a muon beam of a synchrocyclotron was measured to be 8.5 + or - 1.7%. The stopping rate of negative cosmic ray muons at 2-mwe depth below sea level from these measurements and an estimated 17% electromagnetic production is 0.63 + or - 0.13 muon(-)/kg min. Previous measurements on the muon stopping rate vary by a factor of 5. Our value is slightly higher but is consistent with two previous high values. The sensitivity of the Ar-37 radiochemical method for the detection of muons is considerably higher than that of the previous radiochemical methods and could be used to measure the negative muon capture rates at greater depths.
The scientific target of NASA's Magnetospheric Multiscale (MMS) mission is to study the fundamentally important phenomenon of magnetic reconnection. Theoretical models of this process predict a small size, on the order of hundred kilometers, for the ion diffusion region where ions are demagnetized at the dayside magnetopause. This region may typically sweep over the spacecraft at relatively high speeds of 50 km/s, requiring the fast plasma investigation (FPI) instrument suite to have an extremely high time resolution for measurements of the 3D particle distribution functions. As part of the FPI on MMS, the 16 dual ion spectrometers (DIS) will provide fast (150 ms) 3D ion velocity distributions, from 10 to 30,000 eV/q, by combining the measurements from four dual spectrometers on each of four MMS spacecraft. For any multispacecraft mission, the response uniformity among the spectrometer set assumes an enhanced importance. Due to these demanding instrument requirements and the effort of calibrating more than 32 sensors (16 × 2) within a tight schedule, a highly systematic and precise calibration was required for measurement repeatability. To illustrate how this challenge was met, a brief overview of the FPI DIS was presented with a detailed discussion of the calibration method of approach and implementation. Finally, a discussion of DIS performance results, their unit-to-unit variation, and the lessons learned from this calibration effort are presented.
The Fast Plasma Investigation (FPI) on the Magnetospheric Multiscale mission (MMS) combines data from eight spectrometers, each with four deflection states, into a single map of the sky. Any systematic discontinuity, artifact, noise source, etc. present in this map may be incorrectly interpreted as legitimate data and incorrect conclusions reached. For this reason it is desirable to have all spectrometers return the same output for a given input, and for this output to be low in noise sources or other errors. While many missions use statistical analyses of data to calibrate instruments in flight, this process is difficult with FPI for two reasons: 1. Only a small fraction of high resolution data is downloaded to the ground due to bandwidth limitations and 2: The data that is downloaded is, by definition, scientifically interesting and therefore not ideal for calibration. FPI uses a suite of new tools to calibrate in flight. A new method for detection system ground calibration has been developed involving sweeping the detection threshold to fully define the pulse height distribution. This method has now been extended for use in flight as a means to calibrate MCP voltage and threshold (together forming the operating point) of the Dual Electron Spectrometers (DES) and Dual Ion Spectrometers (DIS). A method of comparing higher energy data (which has low fractional voltage error) to lower energy data (which has a higher fractional voltage error) will be used to calibrate the high voltage outputs. Finally, a comparison of pitch angle distributions will be used to find remaining discrepancies among sensors.