Engineering PapersSearch

SEARCH · Engineering Papers

Results for “formulation generation”

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.

At least 37 records · Page 2

Modelling the Shuttle Remote Manipulator System: Another flexible model

High fidelity elastic system modeling algorithms are discussed. The particular system studied is the Space Shuttle Remote Manipulator System (RMS) undergoing full articulated motion. The model incorporates flexibility via a methodology the author has been developing. The technique is based in variational principles, so rigorous boundary condition generation and weak formulations for the associated partial differential equations are realized, yet the analyst need not integrate by parts. The methodology is formulated using vector-dyad notation with minimal use of tensor notation, therefore the technique is believed to be affable to practicing engineers. The objectives of this work are as follows: (1) determine the efficacy of the modeling method; and (2) determine if the method affords an analyst advantages in the overall modeling and simulation task. Generated out of necessity were Mathematica algorithms that quasi-automate the modeling procedure and simulation development. The project was divided into sections as follows: (1) model development of a simplified manipulator; (2) model development of the full-freedom RMS including a flexible movable base on a six degree of freedom orbiter (a rigid-body is attached to the manipulator end-effector); (3) simulation development for item 2; and (4) comparison to the currently used model of the flexible RMS in the Structures and Mechanics Division of NASA JSC. At the time of the writing of this report, items 3 and 4 above were not complete.

Barhorst, Alan A.

A Preliminary Study on the Vapor/Mist Phase Lubrication of a Spur Gearbox

Organophosphates have been the primary compounds used in vapor/mist phase lubrication studies involving ferrous bearing material. Experimental results have indicated that the initial formation of an iron phosphate film on a rubbing ferrous surface, followed by the growth (by cationic diffusion) of a lubricious pyrophosphate-type coating over the iron phosphate, is the reason organophosphates work well as vapor/mist phase lubricants. Recent work, however, has shown that this mechanism leads to the depletion of surface iron atoms and to eventual lubrication failure. A new organophosphate formulation was developed which circumvents surface iron depletion. This formulation was tested by generating an iron phosphate coating on an aluminum surface. The new formulation was then used to vapor/mist phase lubricate a spur gearbox in a preliminary study.

Morales, Wilfredo

Overview of Proposed ISRU Technology Development

ISRU involves any hardware or operation that harnesses and utilizes in-situ resources (natural and discarded) to create products and services for robotic and human exploration: Assessment of physical, mineral chemical, and volatile water resources, terrain, geology, and environment (orbital and local). Production of replacement parts, complex products, machines, and integrated systems from feedstock derived from one or more processed resources. Civil engineering, infrastructure emplacement, and structure construction using materials produced from in situ resources. Radiation shields, landing pads, roads, berms, habitats, etc. Generation and storage of electrical, thermal, and chemical energy with in situ derived materials. Solar arrays, thermal wadis, chemical batteries, etc. ISRU is a disruptive capability: Enables more affordable exploration than todays paradigm. Allows more sustainable architectures to be developed. Understand the ripple effect in the other Exploration Elements: MAV: propellant selection, higher rendezvous altitude (higher DV capable with ISRU propellants). EDL: significantly reduces required landed mass. Life Support: reduce amount of ECLSS closure, reduce trash mass carried through propulsive maneuvers. Power: ISRU drives electrical requirements, reactant and regeneration for fuel cells for landers, rovers, and habitat backup. Every Exploration Element except ISRU has some flight heritage (power, propulsion, habitats, landers, life support, etc.) ISRU will require a flight demonstration mission on Mars before it will be included in the critical path. Mission needs to be concluded at least 10 years before first human landed mission to ensure lessons learned can be incorporated into final design. ISRU Formulation team has generated a (still incomplete) list of over 75 technical questions on more than 40 components and subsystems that need to be answered before the right ISRU system will be ready for this flight demo.

Mars exploration

Stochastic sea state for SRB studies

Ocean surface characteristics at two locations were studied for a space shuttle solid rocket booster ocean impact and recovery analysis. Probability distributions of wave heights, wave slopes, wave velocities, ocean currents, and 1-kilometer altitude winds were formulated. Procedures for generating ocean wave simulations are described.

Perlmutter, M.

Application of second-order turbulent modeling to the prediction of radiated aerodynamic sound

The Ribner formulation of the generation of aerodynamic sound is coupled with predictions of second-order velocity correlations and integral scale to estimate the sound radiated from several complicated jet flows. In particular, it is shown that the sound radiated from a cold swirling jet is greater than from its nonswirling equal thrust counterpart. The noise radiated from the flow field of a multitube suppressor was estimated and compared with an equal thrust diameter Gaussian jet. It is shown that the multitube concept is indeed quieter.

Bilanin, A. J.

The CCSDS Lossless Data Compression Algorithm for Space Applications

In the late 80's, when the author started working at the Goddard Space Flight Center (GSFC) for the National Aeronautics and Space Administration (NASA), several scientists there were in the process of formulating the next generation of Earth viewing science instruments, the Moderate Resolution Imaging Spectroradiometer (MODIS). The instrument would have over thirty spectral bands and would transmit enormous data through the communications channel. This was when the author was assigned the task of investigating lossless compression algorithms for space implementation to compress science data in order to reduce the requirement on bandwidth and storage.

Yeh, Pen-Shu

Ambipolar Electric Field, Photoelectrons, and Their Role in Atmospheric Escape From Hot Jupiters

Atmospheric mass loss from Hot Jupiters can be large due to the close proximity of these planets to their host star and the strong radiation the planetary atmosphere receives. On Earth, a major contribution to the acceleration of atmospheric ions comes from the vertical separation of ions and electrons, and the generation of the ambipolar electric field. This process, known as the "polar wind," is responsible for the transport of ionospheric constituents to Earth's magnetosphere, where they are well observed. The polar wind can also be enhanced by a relatively small fraction of super-thermal electrons (photoelectrons) generated by photoionization.We formulate a simplified calculation of the effect of the ambipolar electric field and the photoelectrons on the ion scale height in a generalized manner. We find that the ion scale height can be increased by a factor of 2-15 due to the polar wind effects. We also estimate a lower limit of an order of magnitude increase of the ion density and the atmospheric mass-loss rate when polar wind effects are included.

atmospheres

Second Generation Prototype Design and Testing for a High Altitude Venus Balloon

This paper describes the development of a second generation prototype balloon intended for flight in the upper atmosphere of Venus. The design of this new prototype incorporates lessons learned from the construction and testing of the first generation prototype, including finite element analyses of the balloon stresses and deformations, measured leak performance after handling and packaging, permeability and optical property measurements on material samples, and sulfuric acid testing. An improved design for the second generation prototype was formulated based on these results, although the spherical shape and 5.5 m diameter size were retained. The resulting balloon has a volume of 87 cubic meters and is capable of carrying a 45 kg payload at a 55 km altitude at Venus. The design and fabrication of the new prototype is described, along with test data for inflation and leakage performance.

mobility

The electron Boltzmann equation in a plasma generated by fission fragments

A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material show that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux but increases sharply in the presence of a sustainer electric field.

Hassan, H. A.

Electron distribution function in a plasma generated by fission fragments

A Boltzmann equation formulation is presented for the determination of the electron distribution function in a plasma generated by fission fragments. The formulation takes into consideration ambipolar diffusion, elastic and inelastic collisions, recombination and ionization, and allows for the fact that the primary electrons are not monoenergetic. Calculations for He in a tube coated with fissionable material shows that, over a wide pressure and neutron flux range, the distribution function is non-Maxwellian, but the electrons are essentially thermal. Moreover, about a third of the energy of the primary electrons is transferred into the inelastic levels of He. This fraction of energy transfer is almost independent of pressure and neutron flux.

Hassan, H. A.

Element Library for Three-Dimensional Stress Analysis by the Integrated Force Method

The Integrated Force Method, a recently developed method for analyzing structures, is extended in this paper to three-dimensional structural analysis. First, a general formulation is developed to generate the stress interpolation matrix in terms of complete polynomials of the required order. The formulation is based on definitions of the stress tensor components in term of stress functions. The stress functions are written as complete polynomials and substituted into expressions for stress components. Then elimination of the dependent coefficients leaves the stress components expressed as complete polynomials whose coefficients are defined as generalized independent forces. Such derived components of the stress tensor identically satisfy homogenous Navier equations of equilibrium. The resulting element matrices are invariant with respect to coordinate transformation and are free of spurious zero-energy modes. The formulation provides a rational way to calculate the exact number of independent forces necessary to arrive at an approximation of the required order for complete polynomials. The influence of reducing the number of independent forces on the accuracy of the response is also analyzed. The stress fields derived are used to develop a comprehensive finite element library for three-dimensional structural analysis by the Integrated Force Method. Both tetrahedral- and hexahedral-shaped elements capable of modeling arbitrary geometric configurations are developed. A number of examples with known analytical solutions are solved by using the developments presented herein. The results are in good agreement with the analytical solutions. The responses obtained with the Integrated Force Method are also compared with those generated by the standard displacement method. In most cases, the performance of the Integrated Force Method is better overall.

Kaljevic, Igor

Generating Global Leaf Area Index from Landsat: Algorithm Formulation and Demonstration

This paper summarizes the implementation of a physically based algorithm for the retrieval of vegetation green Leaf Area Index (LAI) from Landsat surface reflectance data. The algorithm is based on the canopy spectral invariants theory and provides a computationally efficient way of parameterizing the Bidirectional Reflectance Factor (BRF) as a function of spatial resolution and wavelength. LAI retrievals from the application of this algorithm to aggregated Landsat surface reflectances are consistent with those of MODIS for homogeneous sites represented by different herbaceous and forest cover types. Example results illustrating the physics and performance of the algorithm suggest three key factors that influence the LAI retrieval process: 1) the atmospheric correction procedures to estimate surface reflectances; 2) the proximity of Landsatobserved surface reflectance and corresponding reflectances as characterized by the model simulation; and 3) the quality of the input land cover type in accurately delineating pure vegetated components as opposed to mixed pixels. Accounting for these factors, a pilot implementation of the LAI retrieval algorithm was demonstrated for the state of California utilizing the Global Land Survey (GLS) 2005 Landsat data archive. In a separate exercise, the performance of the LAI algorithm over California was evaluated by using the short-wave infrared band in addition to the red and near-infrared bands. Results show that the algorithm, while ingesting the short-wave infrared band, has the ability to delineate open canopies with understory effects and may provide useful information compared to a more traditional two-band retrieval. Future research will involve implementation of this algorithm at continental scales and a validation exercise will be performed in evaluating the accuracy of the 30-m LAI products at several field sites. ©

Landsat

A Corrected Formulation of the Multilayer Model (MLM) for Inferring Gaseous Dry Deposition to Vegetated Surfaces

The Multilayer Model (MLM) has been used for many years to infer dry deposition fluxes from measured trace species concentrations and standard meteorological measurements for national networks in the U.S., including the U.S. Environmental Protection Agency's Clean Air Status and Trends Network (CASTNet). MLM utilizes a resistance analogy to calculate deposition velocities appropriate for whole vegetative canopies, while employing a multilayer integration to account for vertically varying meteorology, canopy morphology and radiative transfer within the canopy. However, the MLM formulation, as it was originally presented and as it has been subsequently employed, contains a non-physical representation related to the leaf-level quasi-laminar boundary layer resistance that affects the calculation of the total canopy resistance. In this note, the non-physical representation of the canopy resistance as originally formulated in MLM is discussed and a revised, physically consistent, formulation is suggested as a replacement. The revised canopy resistance formulation reduces estimates of HNO3 deposition velocities by as much as 38% during mid-day as compared to values generated by the original formulation. Inferred deposition velocities for SO2 and O3 are not significantly altered by the change in formulation (less than 3%). Inferred deposition loadings of oxidized and total nitrogen from CASTNet data may be reduced by 10-20% and 5-10%, respectively, for the Eastern U. S. when employing the revised formulation of MLM as compared to the original formulation.

vegetated surfaces

Aerodynamic shape optimization of wing and wing-body configurations using control theory

This paper describes the implementation of optimization techniques based on control theory for wing and wing-body design. In previous studies it was shown that control theory could be used to devise an effective optimization procedure for airfoils and wings in which the shape and the surrounding body-fitted mesh are both generated analytically, and the control is the mapping function. Recently, the method has been implemented for both potential flows and flows governed by the Euler equations using an alternative formulation which employs numerically generated grids, so that it can more easily be extended to treat general configurations. Here results are presented both for the optimization of a swept wing using an analytic mapping, and for the optimization of wing and wing-body configurations using a general mesh.

Reuther, James

A New Recycling Method to Generate Turbulent Inflow Profiles

The accuracy of the scale-resolving simulations for practical geometries strongly depends on the inflow boundary conditions. Imposing experimentally observed turbulent inflow profiles for the numerical simulations is a major challenge. Existing methods available in the literature assume self-similar behavior, which is not true for most of the experiments. In the present work, we formulate the turbulent inflow profile generation technique as an optimization problem. An adjoint technique is exploited to evaluate the sensitivities of multiple input parameters for the present problem. The present formulation is then tested to generate a laminar boundary layer profile, turbulent boundary layer profile, and turbulent jet profile.

Garai, Anirban

Numerical Study of Boundary Layer Interaction with Shocks: Method Improvement and Test Computation

The objective is the development of a high-order and high-resolution method for the direct numerical simulation of shock turbulent-boundary-layer interaction. Details concerning the spatial discretization of the convective terms can be found in Adams and Shariff (1995). The computer code based on this method as introduced in Adams (1994) was formulated in Cartesian coordinates and thus has been limited to simple rectangular domains. For more general two-dimensional geometries, as a compression corner, an extension to generalized coordinates is necessary. To keep the requirements or limitations for grid generation low, the extended formulation should allow for non-orthogonal grids. Still, for simplicity and cost efficiency, periodicity can be assumed in one cross-flow direction. For easy vectorization, the compact-ENO coupling algorithm as used in Adams (1994) treated whole planes normal to the derivative direction with the ENO scheme whenever at least one point of this plane satisfied the detection criterion. This is apparently too restrictive for more general geometries and more complex shock patterns. Here we introduce a localized compact-ENO coupling algorithm, which is efficient as long as the overall number of grid points treated by the ENO scheme is small compared to the total number of grid points. Validation and test computations with the final code are performed to assess the efficiency and suitability of the computer code for the problems of interest. We define a set of parameters where a direct numerical simulation of a turbulent boundary layer along a compression corner with reasonably fine resolution is affordable.

Adams, N. A.

On the Ice Nucleation Spectrum

This work presents a novel formulation of the ice nucleation spectrum, i.e. the function relating the ice crystal concentration to cloud formation conditions and aerosol properties. The new formulation is physically-based and explicitly accounts for the dependency of the ice crystal concentration on temperature, supersaturation, cooling rate, and particle size, surface area and composition. This is achieved by introducing the concepts of ice nucleation coefficient (the number of ice germs present in a particle) and nucleation probability dispersion function (the distribution of ice nucleation coefficients within the aerosol population). The new formulation is used to generate ice nucleation parameterizations for the homogeneous freezing of cloud droplets and the heterogeneous deposition ice nucleation on dust and soot ice nuclei. For homogeneous freezing, it was found that by increasing the dispersion in the droplet volume distribution the fraction of supercooled droplets in the population increases. For heterogeneous ice nucleation the new formulation consistently describes singular and stochastic behavior within a single framework. Using a fundamentally stochastic approach, both cooling rate independence and constancy of the ice nucleation fraction over time, features typically associated with singular behavior, were reproduced. Analysis of the temporal dependency of the ice nucleation spectrum suggested that experimental methods that measure the ice nucleation fraction over few seconds would tend to underestimate the ice nuclei concentration. It is shown that inferring the aerosol heterogeneous ice nucleation properties from measurements of the onset supersaturation and temperature may carry significant error as the variability in ice nucleation properties within the aerosol population is not accounted for. This work provides a simple and rigorous ice nucleation framework where theoretical predictions, laboratory measurements and field campaign data can be reconciled, and that is suitable for application in atmospheric modeling studies.

Barahona, D.

Hyperbolic Methods for Surface and Field Grid Generation

This chapter describes the use of hyperbolic partial differential equation methods for structured surface grid generation and field grid generation. While the surface grid generation equations are inherently three dimensional, the field grid generation equations can be formulated in two or three dimensions. The governing equations are derived from orthogonality relations and cell area/volume constraints; and are solved numerically by marching from an initial curve or surface. The marching step size and marching distance can be prescribedly the user. Exact specifications of the side and outer boundaries are not possible with a one sweep marching scheme but limited control is achievable. Excellent orthogonality and grid clustering characteristics are provided by hyperbolic methods with one to two orders of magnitude savings in time over typical elliptic methods. Since hyperbolic grid generation methods do not require the exact specifications of the side and outer boundaries of a grid, these methods are particularly well suited for the overlapping grid approach for solving problems on complex configurations. Grid generation software based on hyperbolic methods and their applications on several complex configurations will be described.

Chan, William M.