Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient methods”

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 1,297 records · Page 72

Thermal Design and Analysis for the Cryogenic MIDAS Experiment

The Materials In Devices As Superconductors (MIDAS) spaceflight experiment is a NASA payload which launched in September 1996 on the Shuttle, and was transferred to the Mir Space Station for several months of operation. MIDAS was developed and built at NASA Langley Research Center (LaRC). The primary objective of the experiment was to determine the effects of microgravity and spaceflight on the electrical properties of high-temperature superconductive (HTS) materials. The thermal challenge on MIDAS was to maintain the superconductive specimens at or below 80 K for the entire operation of the experiment, including all ground testing and 90 days of spaceflight operation. Cooling was provided by a small tactical cryocooler. The superconductive specimens and the coldfinger of the cryocooler were mounted in a vacuum chamber, with vacuum levels maintained by an ion pump. The entire experiment was mounted for operation in a stowage locker inside Mir, with the only heat dissipation capability provided by a cooling fan exhausting to the habitable compartment. The thermal environment on Mir can potentially vary over the range 5 to 40 C; this was the range used in testing, and this wide range adds to the difficulty in managing the power dissipated from the experiment's active components. Many issues in the thermal design are discussed, including: thermal isolation methods for the cryogenic samples; design for cooling to cryogenic temperatures; cryogenic epoxy bonds; management of ambient temperature components self-heating; and fan cooling of the enclosed locker. Results of the design are also considered, including the thermal gradients across the HTS samples and cryogenic thermal strap, electronics and thermal sensor cryogenic performance, and differences between ground and flight performance. Modeling was performed in both SINDA-85 and MSC/PATRAN (with direct geometry import from the CAD design tool Pro/Engineer). Advantages of both types of models are discussed. Correlation of several models to ground testing and flight data (where available) is presented. Both SINDA and PATRAN models predicted the actual thermal performance of the experiment well, even without post-flight correlation adjustments of the models.

Amundsen, Ruth McElroy↗

Application of the generalized vertical coordinate ocean model for better representing satellite data

It is found that two adaptive parametric functions can be introduced into the basic ocean equations for utilizing the optimal or hybrid features of commonly used z-level, terrain- following, isopycnal, and pressure coordinates in numerical ocean models. The two parametric functions are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara (1 974), the Jacobian pressure gradient formulation of Song (1 998), and a newly developed metric factor that permits both compressible (non-Boussinesq) and incompressible (Boussinesq) approximations. Based on the new formulation, an adaptive modeling strategy is proposed and a staggered finite volume method is designed to ensure conservation of important physical properties and numerical accuracy. Implementation of the combined techniques to SCRUM (Song and Haidvogel1994) shows that the adaptive modeling strategy can be applied to any existing ocean model without incurring computational expense or altering the original numerical schemes. Such a generalized coordinate model is expected to benefit diverse ocean modelers for easily choosing optimal vertical structures and sharing modeling resources based on a common model platform. Several representing oceanographic problems with different scales and characteristics, such as coastal canyons, basin-scale circulation, and global ocean circulation, are used to demonstrate the model's capability for multiple applications. New results show that the model is capable of simultaneously resolving both Boussinesq and non-Boussinesq, and both small- and large-scale processes well. This talk will focus on its applications of multiple satellite sensing data in eddy-resolving simulations of Asian Marginal Sea and Kurosio. Attention will be given to how Topex/Poseidon SSH, TRMM SST; and GRACE ocean bottom pressure can be correctly represented in a non- Boussinesq model.

numerical ocean model generalized coordinate syste↗

Bulk Crystal Growth of Piezoelectric PMN-PT Crystals Using Gradient Freeze Technique for Improved SHM Sensors

There has been a growing interest in recent years in lead based perovskite ferroelectric and relaxor ferroelectric solid solutions because of their excellent dielectric, piezoelectric and electrostrictive properties that make them very attractive for various sensing, actuating and structural health monitoring (SHM) applications. We are interested in the development of highly sensitive and efficient PMN-PT sensors based on large single crystals for the structural health monitoring of composite materials that may be used in future spacecrafts. Highly sensitive sensors are needed for detection of defects in these materials because they often tend to fail by distributed and interacting damage modes and much of the damage occurs beneath the top surface of the laminate and not detectable by visual inspection. Research is being carried out for various combinations of solid solutions for PMN-PT piezoelectric materials and bigger size crystals are being sought for improved sensor applications. Single crystals of this material are of interest for sensor applications because of their high piezoelectric coefficient (d33 greater than 1700 pC/N) and electromechanical coefficients (k33 greater than 0.90). For comparison, the commonly used piezoelectric ceramic lead zirconate titanate (PZT) has a d33 of about 600 pC/N and electromechanical coefficients k33 of about 0.75. At the present time, these piezoelectric relaxor crystals are grown by high temperature flux growth method and the size of these crystals are rather small (~3x4x5 mm(exp 3). In the present paper, we have attempted to grow bulk single crystals of PMN-PT in a 2 inch diameter platinum crucible and successfully grown a large size crystal of 67%PMN-33%PT using the vertical gradient freeze technique with no flux. Piezoelectric properties of the grown crystals are investigated. PMN-PT plates show excellent piezoelectric properties. Samples were poled under an applied electric field of 5 kV/cm. Dielectric properties at a frequency of 1 kHz are examined. The grown PMN-PT crystals show typical relaxor dielectric properties. Additionally, the thermal properties of the sample are tested. The results are in good agreement with those found in the literature and some are reported for the first time.

Aggarwal, Mohan D.↗

Stability of the Dual-Frequency Radar Equations and a New Method Applied to the GPM’s Dual-Frequency Precipitation Radar (DPR) Data

A new algorithm is proposed that estimates two parameters of the particle size distribution (PSD) at each range bin from the Global Precipitation Measurement’s (GPM’s) Dual-frequency Precipitation Radar data. The equation that expresses the relationship of the PSD parameters between adjacent range bins is derived. By including the attenuation effect within the bin in the discretized equation, the new algorithm alleviates the double-solution problem when attenuation within the bin is sufficiently large. The stability of the solutions to the equation depends on the value of the mean diameter D(sub m) and its gradient with respect to range in the case of liquid precipitation. If the critical diameter above which the dual-frequency ratio (DFR) of radar reflectivity factors becomes a monotonically increasing function of D(sub m), the backward processing of the equation provides a stable or moderately diverging solution, unlike the forward processing that often gives unstable solutions. To provide a set of initial conditions without using the surface reference technique (SRT) in the backward processing, an initialization method using the Hitschfeld-Bordan attenuation correction method is proposed and tested. The proposed algorithm may provide a tool for investigating the assumptions used in various algorithms.

GPM↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗

Application of a self-adaptive grid method to complex flows

A directional-split modular user-friendly grid-point distribution code is applied to several test problems. Grid-point spacing is determined by user-specified constants denoting maximum and minimum grid spacings and relating the relative influence of smoothness and orthogonality. Estimates of truncation error, in terms of flowfield gradients and/or geometric features, are used to determine the point distribution.

Deiwert, G. S.↗

Vibrational Power Flow Analysis of Rods and Beams

A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.

Wohlever, James Christopher↗

Poisson equations of rotational motion for a rigid triaxial body with application to a tumbling artificial satellite

A mathematical model is developed for studying the effects of gravity gradient torque on the attitude stability of a tumbling triaxial rigid satellite. Poisson equations are used to investigate the rotation of the satellite (which is in elliptical orbit about an attracting point mass) about its center of mass. An averaging method is employed to obtain an intermediate set of differential equations for the nonresonant, secular behavior of the osculating elements which describe the rotational motions of the satellite, and the averaged equations are then integrated to obtain long-term secular solutions for the osculating elements.

Liu, J. J. F.↗

A technique for determining cloud free vs cloud contaminated pixels in satellite imagery

Since the first earth orbiting satellite sent pictures of the earth back to them, atmospheric scientists have been focused on the possibilities of using that information as both a forecasting tool and as a meteorological research tool. With the latest generation of Geostationary Operational Environmental Satellites (GOES) now entering service, that view of the earth yields views at a frequency and resolution never before available. These satellites have imagers with a five band multi-spectral capability with high spatial resolution. In addition, the sounder has eighteen thermal infrared (IR) channels plus one low-resolution visible band. With a resolution as small as one kilometer, GOES provides scientists with a powerful eye on the atmosphere. Menzel and Purdom (1994) detail both the imager and sounder capability as well as other systems on the GOES satellites. Immediately apparent in the visible channel are the patterns of clouds swirling over both oceans and continents. These clouds range in size from huge planetary systems covering thousands of kilometers to puffy fair weather cumulus clouds on the order of half a kilometer in size. With the IR sensors temperature patterns are observed. High clouds appear very cold, while low stratus field show temperatures near that of the surface. The surface, in turn, generally appears warmer than the clouds. It would seem then a simple manner to determine cloud and surface temperature from the imagery, but such is not the case. While most of the atmospheric constituents are well mixed and homogeneous, water vapor is not. The water molecule, because of its unique structure and vibration modes, affects the transmittance of the atmosphere most notably in the infrared regions. There are regions of the IR spectrum where water vapor acts as a strong absorber, and at others it is nearly transparent. The transparent wavelengths are called windows, and one such window occurs at 11.2 microns. Adjacent to this window at 12.7 microns which is strongly absorbed by water vapor. These two wavelengths form what is known as a split window, the utility of which was used. Using the linearized form of the radiative transfer equation, they were able to use the split window to determine the amount of water vapor present in the atmosphere. Jedlovec developed the physical split-window (PSW) technique which determines the integrated water content (IWC). The PSW method using Visible Infrared Spin Scan Radiometer (VISSR) Atmospheric Sounder (VAS) found on the older versions of the GOES satellites was used. Recently, Jedlovec and colleagues have been attempting to apply the PSW method using full disk IR imagery obtained by the new generation of GOES satellites. IWC is essential for improved analysis and prediction of convective storms which have been observed to develop in regions of both strong and rapidly evolving moisture gradients. It has also been used in the prediction of clouds and precipitation.

Wohlman, Richard A.↗

Advanced Global Atmospheric Gases Experiment (AGAGE): MIT Contribution

We describe in detail the instrumentation and calibrations used in the ALE, GAGE and AGAGE experiments and present a history of the majority of the anthropogenic ozone- depleting and climate-forcing gases in air based on these experiments. Beginning in 1978, these three successive automated high frequency in-situ experiments have documented the long-term behavior of the measured concentrations of these gases over the past twenty years, and show both the evolution of latitudinal gradients and the high frequency variability due to sources and circulation. We provide estimates of the long-term trends in total chlorine contained in long- lived halocarbons involved in ozone depletion. We summarize interpretations of these measurements using inverse methods to determine trace gas lifetimes and emissions. Finally, we provide a combined observational and modeled reconstruction of the evolution of chlorocarbons by latitude in the atmosphere over the past sixty years which can be used as boundary conditions for interpreting trapped air in glaciers and oceanic measurements of chlorocarbon tracers of the deep oceanic circulation. Some specific conclusions are: (a) International compliance with the Montreal Protocol is so far resulting in chlorofluorocarbon and chlorocarbon mole fractions comparable to target levels, (b) Mole fractions of total chlorine contained in long-lived halocarbons (CCl2F2, CCl3F, CH3CCl3, CCl4, CHClF2, CCl2FCClF2, CH3Cl, CH2Cl2, CHCl3, CCl2=CCl2) in the lower troposphere reached maximum values of about 3.6 ppb in 1993 and are beginning to slowly decrease in the global lower atmosphere, (c) The chlorofluorocarbons have atmospheric lifetimes consistent with destruction in the stratosphere being their principal removal mechanism, (d) Multi-annual variations in chlorofluorocarbon and chlorocarbon emissions deduced from ALUGAGWAGAGE data are consistent approximately with variations estimated independently from industrial production and sales data where available (CCl2F2 (CFC-12) and CCl2FCClF2 (CFC-113) show the greatest discrepancies), (e) The mole fractions of the hydrochlorofluorocarbons and hydrofluorocarbons, which are replacing the regulated halocarbons, are rising very rapidly in the atmosphere but, with the exception of the much longer manufactured CHClF2 (HCFC-22), they are not yet at levels sufficient to contribute significantly to atmospheric chlorine loading. These replacement species could in the future provide independent estimates of the global weighted-average OH concentration provided their industrial emissions are accurately documented, (f) In the future, analysis of pollution events measured using high frequency in-situ measurements of chlorofluorocarbons and their replacements may enable emission estimates at the regional level which, together with industrial end-use data, are of sufficient accuracy to be capable of identifying regional non-compliance with the Montreal Protocol.

Kurylo, Michael↗

Exploring the Science Trade Space with the JPL Innovation Foundry A-Team

The Jet Propulsion Laboratory Innovation Foundry has established a new approach for exploring, developing, and evaluating early concepts with a group called the Architecture Team. The Architecture Team combines innovative collaborative methods and facilitated sessions with subject matter experts and analysis tools to help mature mission concepts. Science, implementation, and programmatic elements are all considered during an A-Team study. In these studies, Concept Maturity Levels are used to group methods. These levels include idea generation and capture (Concept Maturity Level 1), initial feasibility assessment (Concept Maturity Level 2), and trade space exploration (Concept Maturity Level 3). Methods used for exploring the science objectives, feasibility, and scope will be described including the use of a new technique for understanding the most compelling science, called a Science Return Diagram. In the process of developing the Science Return Diagram, gradients in the science trade space are uncovered along with their implications for implementation and mission architecture. Special attention is paid toward developing complete investigations, establishing a series of logical claims that lead to the natural selection of a measurement approach. Over 20 science-focused A-Team studies have used these techniques to help science teams refine their mission objectives, make implementation decisions, and reveal the mission concept's most compelling science. This article will describe the A-Team process for exploring the mission concept's science trade space and the Science Return Diagram technique.

science↗

Validation of Cryogenic Propellant Tank Self-Pressurization by Leveraging Reduced Order Modeling within Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank self-pressurization was performed using a hybrid Computational Fluid Dynamics (CFD) and reduced order modeling methodology. Data from a liquid hydrogen ground test conducted at the K-site facility at the National Aeronautics and Space Administration (NASA) Glenn Research Center was used for the validation effort. Liquid phase dynamics were explicitly resolved with a CFD tool. Vapor phase dynamics were modeled as a point mass that communicated heat from the tank wall to the liquid phase via a boundary condition used at the gas-liquid interface. The method proved to be more accurate, robust, and efficient than explicit resolution of the dynamics using a standard Volume of Fluid (VOF) methodology. The subject pressurization process was found to be heavily dependent upon both the relatively high liquid temperature gradient near the gas-liquid interface and the natural convection flow path. Modeling the gas-liquid interface as an immovable surface eliminated temperature gradient destroying gas-liquid interface velocities observed in VOF simulations, and correspondingly enabled more rapid simulation since interface advection was not allowed. The single phase computational domain also facilitated the ability to demonstrate spatial resolution convergence of natural convection cells within the liquid which significantly impacted the tank pressurization rate. This work was used to demonstrate the critical physics for tank self-pressurization and numerical methodologies that may be used to best resolve those physics. The findings informed development and operation of production level CFD tools used in the Fluid Dynamics Branch at NASA Marshall Space Flight Center.

J. M. Brodnick↗

Experimental examination of similarity measures and preprocessing methods used for image registration

The criterion used to measure the similarity between images and thus find the position where the images are registered is examined. The three similarity measures considered are the correlation coefficient, the sum of the absolute differences, and the correlation function. Three basic types of preprocessing are then discussed: taking the magnitude of the gradient of the images, thresholding the images at their medians, and thresholding the magnitude of the gradient of the images at an arbitrary level to be determined experimentally. These multitemporal registration techniques are applied to remote imagery of agricultural areas.

Svedlow, M.↗

Preventing Damaging Pressure Gradients at the Walls of an Inflatable Space System

An inflatable structural system to deploy a space system such as a solar shield, an antenna or another similar instrument requires a stiffening element after it is extended by the inflated gas pressure. The stiffening element has to be packaged in folded configuration before the deployment. It must be relatively small, lightweight, non-damaging to the inflated system and be able to become stiff in a short time. One stiffening method is to use a flexible material inserted in the deployable system, which, upon a temperature curing, can become stiff and is capable of supporting the entire structure. There are two conditions during the space operations when the inflated volume could be damaged: during the transonic region of the launch phase and when the curing of the rigidizing element occurs. In both cases, an excess of pressure within the volume containing the rigid element could burst the walls of the low-pressure gas inflated portion of the system. This paper investigates those two conditions and indicates the vents, which will prevent those damaging overpressures. Vent openings at the non-inflated volumes have been calculated for the conditions existing during the launch. Those vents allow the initially folded volume to exhaust the trapped atmospheric gas at approximately the same rate as the ambient pressure drops. That will prevent pressure gradients across the container walls which otherwise could be as high as 14.7 psi. The other condition occurring during the curing of the stiffening element has been investigated. This has required the testing of the element to obtain the gas generation during the curing and the transformation from a pliable material to a rigid on The tested material is a composite graphite/epoxy weave. The outgassing of the uncured sample at 121 deg Celcius was carried with the Cahn Microbalance and with other outgassing facilities including the micro-CVCM ASTM E-595 facility. The test provided the mass of gas evolved during the test. That data, including the chemical nature of the evolved gas, provided the data for the calculation of the pressure produced within the volume. The evaluation of the areas of the vents that would prevent excessive pressures and provide a rapid release of the gas away from contamination sensitive surfaces has be carried out. The pressure decay with time has been indicated.

Scialdone, John J.↗

Electrophoretic separation of proteins in space

Commercially available and synthetic wide range and short range ampholytes used in the isoelectric focusing of proteins was analyzed by ion exchange chromatography. A pH gradient over the pH range 3.8 to 11.0 was used to elute the ampholytes from a column of a sulfonated polystyrene resin. The wide range ampholytes were resolved into some 60 to 70 ninhydrin positive components. The recovery obtained with the method was quantitative. Acid short range ampholytes have approximately 35 components which elute readily from the ion exchange resin. Basic short range ampholytes gave about 50 components, most of which eluted at alkaline pH.

Brown, R. K.↗

Application of Gauss-Seidel multilevel control to a single axis torsional model

An approach to the application of multilevel control techniques for large space structures is presented. Gauss-Seidel second level controls and an extension of standard linear quadratic regulator techniques are applied to a model consisting of a flexible space vehicle comprising three rigid bodies rotating about a common axis. The method incorporates second order derivatives with respect to time, entailing the use of a two level hierarchy of subsystems. Gauss-Seidel second level control formulations were chosen to avoid the necessity of having as many controls as constraints or using gradient techniques in the choice of the Hamiltonian. The modular nature of the resulting control system allows applications for spacecraft with larger numbers of modules.

Chichester, F. D.↗

Engineering calculations for communications satellite systems planning

Observed solution times were analyzed for the extended gradient and cyclic coordinate search procedures. The times used in the analysis come from computer runs made during a previously-reported experiment conducted to assess the quality of the solutions to a BSS synthesis problem found by the two search methods. The results of a second experiment with a Fixed Satellite Service (FSS) test problem are also presented. Computational results are summarized for mixed integer programming approaches for solving FSS synthesis problems. A promising heuristic algorithm is described. A synthesis model is discussed for orbital arc allotment optimization. Research plans for the near future are also presented.

Reilly, C. H.↗

Optimal cooperative time-fixed impulsive rendezvous

A method has been developed for determining optimal, i.e., minimum fuel, trajectories for the fixed-time cooperative rendezvous of two spacecraft. The method presently assumes that the vehicles perform a total of three impulsive maneuvers with each vehicle being active, that is, making at least one maneuver. The cost of a feasible 'reference' trajectory is improved by an optimizer which uses an analytical gradient developed using primer vector theory and a new solution for the optimal terminal (rendezvous) maneuver. Results are presented for a large number of cases in which the initial orbits of both vehicles are circular but in which the initial positions of the vehicles and the allotted time for rendezvous are varied. In general, the cost of the cooperative rendezvous is less than that of rendezvous with one vehicle passive. Further improvement in cost may be obtained in the future when additional, i.e., midcourse, impulses are allowed and inserted as indicated for some cases by the primer vector histories which are generated by the program.

Mirfakhraie, Koorosh↗