Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “algorithmic differentiation”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Band-Moment Compression of AVIRIS Hyperspectral Data and its Use in the Detection of Vegetation Stress

A remote sensing campaign was conducted over a U.S. Department of Agriculture test farm at Shelton, Nebraska. An experimental field was set off in plots that were differentially treated with anhydrous ammonia. Four replicates of 0-kg/ha to 200-kg/ha plots, in 50-kg/ha increments, were set out in a random block design. Low-altitude (GSD of 3 m) Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) hyperspectral data were collected over the site in 224 bands. Simultaneously, ground data were collected to support the airborne imagery. In an effort to reduce data load while maintaining or enhancing algorithm performance for vegetation stress detection, band-moment compression and analysis was applied to the AVIRIS image cube. The results indicated that band-moment techniques compress the AVIRIS dataset significantly while retaining the capability of detecting environmentally induced vegetation stress.

Estep, L.↗
NASA NTRS1994

The formation of small grains in shocks in the ISM

Carbonaceous and silicate grains swept up, and betatron accelerated, by supernova-generated shock waves in the interstellar medium are exposed to grain destructive processing. The degree of grain destruction is determined by the differential gas-grain and grain-grain velocities, which lead to sputtering of the grain surface and grain core disruption (deformation, vaporization and shattering), respectively. The threshold pressure for grain shattering in grain-grain collisions (100 k bar) is considerably lower than that for vaporization (approximately 5 M bar). Therefore, collisions between grains shatter large grains into smaller fragments (i.e., small grains and PAH's). Using a new algorithms for the destructive processes, it was possible to model the formation fo small grain fragments in grain-grain collisions in the warm phase of the interstellar medium. It was found that in one cycle through the warm medium (approximately 3 x 10(sup 6) years) of order 1-2% of the total grain mass is shattered into particles with radii of less than 50 A.

Jones, Anthony P.↗
NASA NTRS2012

Airborne Lidar Measurements of Atmospheric Pressure Made Using the Oxygen A-Band

Accurate measurements of greenhouse gas mixing ratios on a global scale are currently needed to gain a better understanding of climate change and its possible impact on our planet. In order to remotely measure greenhouse gas concentrations in the atmosphere with regard to dry air, the air number density in the atmosphere is also needed in deriving the greenhouse gas concentrations. Since oxygen is stable and uniformly mixed in the atmosphere at 20.95%, the measurement of an oxygen absorption in the atmosphere can be used to infer the dry air density and used to calculate the dry air mixing ratio of a greenhouse gas, such as carbon dioxide or methane. OUT technique of measuring Oxygen uses integrated path differential absorption (IPDA) with an Erbium Doped Fiber Amplifier (EDF A) laser system and single photon counting module (SPCM). It measures the absorbance of several on- and off-line wavelengths tuned to an O2 absorption line in the A-band at 764.7 nm. The choice of wavelengths allows us to maximize the pressure sensitivity using the trough between two absorptions in the Oxygen A-band. Our retrieval algorithm uses ancillary meteorological and aircraft altitude information to fit the experimentally obtained lidar O2 line shapes to a model atmosphere and derives the pressure from the profiles of the two lines. We have demonstrated O2 measurements from the ground and from an airborne platform. In this paper we will report on our airborne measurements during our 2011 campaign for the ASCENDS program.

Riris, Haris↗
NASA NTRSDate not supplied

Obtaining Accurate Derivatives Across Adaptive-Step Integration

Obtaining derivatives through the integration of a set of ordinary differential equations is of interest to a number of fields of research. Approaches for doing so that are available in the literature, however, do not account for details associated with variable step integrator that cause inaccuracies in the total derivatives. In this paper, we explore the causes of the accumulation of derivative error in adaptive-step Runge-Kutta integration schemes and propose a solution to eliminate them. The new algorithm proposed is a transformation of each step of a Runge-Kutta integration into a set of residual equations. This work demonstrates the improvement in accuracy of total derivatives when the new algorithm is utilized.

Kaushik Ponnapalli↗
NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗
NASA NTRS2014

A Preliminary Analysis of Precipitation Properties and Processes during NASA GPM IFloodS

The Iowa Flood Studies (IFloodS) is a NASA Global Precipitation Measurement (GPM) ground measurement campaign, which took place in eastern Iowa from May 1 to June 15, 2013. The goals of the field campaign were to collect detailed measurements of surface precipitation using ground instruments and advanced weather radars while simultaneously collecting data from satellites passing overhead. Data collected by the radars and other ground instruments, such as disdrometers and rain gauges, will be used to characterize precipitation properties throughout the vertical column, including the precipitation type (e.g., rain, graupel, hail, aggregates, ice crystals), precipitation amounts (e.g., rain rate), and the size and shape of raindrops. The impact of physical processes, such as aggregation, melting, breakup and coalescence on the measured liquid and ice precipitation properties will be investigated. These ground observations will ultimately be used to improve rainfall estimates from satellites and in particular the algorithms that interpret raw data for the upcoming GPM mission's Core Observatory satellite, which launches in 2014. The various precipitation data collected will eventually be used as input to flood forecasting models in an effort to improve capabilities and test the utility and limitations of satellite precipitation data for flood forecasting. In this preliminary study, the focus will be on analysis of NASA NPOL (S‐band, polarimetric) radar (e.g., radar reflectivity, differential reflectivity, differential phase, correlation coefficient) and NASA 2D Video Disdrometers (2DVDs) measurements. Quality control and processing of the radar and disdrometer data sets will be outlined. In analyzing preliminary cases, particular emphasis will be placed on 1) documenting the evolution of the rain drop size distribution (DSD) as a function of column melting processes and 2) assessing the impact of range on ground‐based polarimetric radar estimates of DSD properties.

Carey, Lawrence↗
NASA NTRS1989

A vectorized algorithm for 3D dynamics of a tethered satellite

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Wilson, Howard B.↗
NASA NTRS2013

SMACK: A New Algorithm for Modeling Collisions and Dynamics of Planetesimals in Debris Disks

We present the Superparticle Model/Algorithm for Collisions in Kuiper belts and debris disks (SMACK), a new method for simultaneously modeling, in 3-D, the collisional and dynamical evolution of planetesimals in a debris disk with planets. SMACK can simulate azimuthal asymmetries and how these asymmetries evolve over time. We show that SMACK is stable to numerical viscosity and numerical heating over 10(exp 7) yr, and that it can reproduce analytic models of disk evolution. We use SMACK to model the evolution of a debris ring containing a planet on an eccentric orbit. Differential precession creates a spiral structure as the ring evolves, but collisions subsequently break up the spiral, leaving a narrower eccentric ring.

circumstellar matter↗
NASA NTRS1995

Two-dimensional mesh embedding for Galerkin B-spline methods

A number of advantages result from using B-splines as basis functions in a Galerkin method for solving partial differential equations. Among them are arbitrary order of accuracy and high resolution similar to that of compact schemes but without the aliasing error. This work develops another property, namely, the ability to treat semi-structured embedded or zonal meshes for two-dimensional geometries. This can drastically reduce the number of grid points in many applications. Both integer and non-integer refinement ratios are allowed. The report begins by developing an algorithm for choosing basis functions that yield the desired mesh resolution. These functions are suitable products of one-dimensional B-splines. Finally, test cases for linear scalar equations such as the Poisson and advection equation are presented. The scheme is conservative and has uniformly high order of accuracy throughout the domain.

Shariff, Karim↗
NASA NTRS2014

A Walsh Function Module Users' Manual

The solution of partial differential equations (PDEs) with Walsh functions offers new opportunities to simulate many challenging problems in mathematical physics. The approach was developed to better simulate hypersonic flows with shocks on unstructured grids. It is unique in that integrals and derivatives are computed using simple matrix multiplication of series representations of functions without the need for divided differences. The product of any two Walsh functions is another Walsh function - a feature that radically changes an algorithm for solving PDEs. A FORTRAN module for supporting Walsh function simulations is documented. A FORTRAN code is also documented with options for solving time-dependent problems: an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the usage of the Walsh function module including such features as operator overloading, Fast Walsh Transforms in multi-dimensions, and a Fast Walsh reciprocal.

Gnoffo, Peter A.↗
NASA NTRS1986

Computer-aided design and distributed system technology development for large space structures

Proposed large space structures have many characteristics that make them difficult to analyze and control. They are highly flexible, with components mathematically modeled by partial differential equations or very large systems of ordinary differential equations. They have many resonant frequencies, typically low and closely spaced. Natural damping may be low and/or improperly modeled. Coupled with stringent operational requirements of orientation, shape control, and vibration suppression, and the inability to perform adequate ground testing, these characteristics present an unconventional identification and control design problem to the systems theorist. Some of the research underway within Langley's Spacecraft Control Branch, Guidance and Control Division aimed at developing theory and algorithms to treat large space structures systems identification and control problems is described. The research areas to be considered are computer-aided design algorithms, and systems identification and control of distributed systems.

Armstrong, Ernest S.↗
NASA NTRS1984

On initial orbit determination

The classical methods of initial orbit determination are brought together within a larger viewpoint. This new synthesis stresses that all such techniques follow one of three approaches. Either they seek to compute the orbital element set, or its equivalent, by attacking the differential equations of motion (Laplace), the first integrals of the equations of motion (Taff), or the solution itself (Gauss). The particular technique pursued within a given type of approach should depend upon the nature of the observational data, the amount of a priori information one is willing to presume, and the object of the exercise. This might be a binary star system, a moon, a minor planet, or an artificial satellite. The efficacy of some algorithms for each approach is discussed briefly. Unfortunately, none of them work very well. Extensions of these techniques to radars or laser radars are trivial and have provided no new insights into the overall problem.

Taff, L. G.↗
NASA NTRS2015

Technology Advancement for Active Remote Sensing of Carbon Dioxide from Space Using the ASCENDS CarbonHawk Experiment Simulator: First Results

The ASCENDS CarbonHawk Experiment Simulator (ACES) is a newly developed lidar developed at NASA Langley Research Center and funded by NASA's Earth Science Technology Office (ESTO) Instrument Incubator Program (IIP) that seeks to advance technologies critical to measuring atmospheric column carbon dioxide (CO2) mixing ratios in support of the NASA Active Sensing of CO2 Emissions over Nights, Days, and Seasons (ASCENDS) mission. The technology advancements targeted include: (1) increasing the power-aperture product to approach ASCENDS mission requirements by implementing multi-aperture telescopes and multiple co-aligned laser transmitters; (2) incorporating high-efficiency, high-power Erbium-Doped Fiber Amplifiers (EDFAs); (3) developing and incorporating a high-bandwidth, low-noise HgCdTe detector and transimpedence amplifier (TIA) subsystem capable of long-duration autonomous operation on Global Hawk aircraft, and (4) advancing algorithms for cloud and aerosol discrimination. The ACES instrument architecture is being developed for operation on high-altitude aircraft and will be directly scalable to meet the ASCENDS mission requirements. These technologies are critical towards developing not only spaceborne instruments but also their airborne simulators, with lower platform requirements for size, mass, and power, and with improved instrument performance for the ASCENDS mission. ACES transmits five laser beams: three from commercial EDFAs operating near 1.57 microns, and two from the Exelis oxygen (O2) Raman fiber laser amplifier system operating near 1.26 microns. The three EDFAs are capable of transmitting up to 10 watts average optical output power each and are seeded by compact, low noise, stable, narrow-linewidth laser sources stabilized with respect to a CO2 absorption line using a multi-pass gas absorption cell. The Integrated-Path Differential Absorption (IPDA) lidar approach is used at both wavelengths to independently measure the CO2 and O2 column number densities and retrieve the average column CO2 mixing ratio. The ACES receiver uses three fiber-coupled 17.8-cm diameter athermal telescopes. The transmitter assembly consists of five fiber-coupled laser collimators and an associated Risley prism pair for each laser to co-align the outgoing laser beams and to align them with the telescope field of view. The backscattered return signals collected by the three telescopes are combined in a fiber bundle and sent to a single low noise detector. The detector/TIA development has improved the existing detector subsystem by increasing its bandwidth to 4.7 MHz from 500 kHz and increasing the duration of autonomous, service-free operation periods from 4 hours to >24 hours. The new detector subsystem enables the utilization of higher laser modulation rates, which provides greater flexibility for implementing advanced thin-cloud discrimination algorithms as well as improving range-determination resolution and error reduction. The cloud/aerosol discrimination algorithm development by Langley and Exelis features a new suite of algorithms for the minimization/elimination of bias errors in the return signal induced by the presence of intervening thin clouds. Multiple laser modulation schemes are being tested in an effort to significantly mitigate the effects of thin clouds on the retrieved CO2 column amounts. Full instrument development concluded in the spring of 2014. After ground range tests of the instrument, ACES successfully completed six test flights on the Langley Hu-25 aircraft in July, 2014, and recorded data at multiple altitudes over land and ocean surfaces with and without intervening clouds. Preliminary results from these test flights will be presented in this paper.

Obland, Michael D.↗
NASA NTRS1972

A numerical investigation of film cooling as an alternate thermal protection system for spacecraft during atmospheric reentry

Research was undertaken to determine the mass required to protect a surface during atmospheric reentry by injecting a cold gas tangentially along the surface. A consideration of this problem revealed that the boundary layer equations could be employed to approximate the most important flow configurations which occur in film cooling. However, even with this simplification, the equations to be solved are complex, nonlinear, coupled partial differential equations. Due to the relative weakness of classical methods of boundary layer analysis, finite difference methods were chosen as the best approach to the formidable problem under consideration. Instead of blindly adopting a proposed technique and trying to force it to the problem, a more systematic study was deemed to be not only desirable but necessary. A relatively basic study of the application of finite difference methods was undertaken with the objective of developing an algorithm capable of obtaining accurate solutions to film cooling problems.

Clausing, A. M.↗
NASA NTRS1987

Characterizing forest stands with multi-incidence angle and multi-polarized SAR data

The potential for using HH-polarized L-band SAR data obtained at different incidence angles from satellite altitudes to identify and map different forest cover types and stand density classes is studied. Reasonably accurate results are obtained if the speckle characteristics of the data are suppressed by low-pass spatial filters and a contextual classification algorithm. Multipolarized L-band SAR data obtained from aircraft altitudes over the same test site are also analyzed to assess the relationships between polarization and forest stand characteristics. It is found that incidence angle controls, to a very large extent, the characteristics of the data and the type of information that can be obtained from L-band, HH-polarized satellite SAR data. Cross-polarization of L-band SAR data enhances and differentiates various forest stand characteristics which cannot be defined using only the like-polarized data, and vice-versa.

Hoffer, R. M.↗
NASA NTRS1992

Solder creep-fatigue interactions with flexible leaded parts

With flexible leaded parts, the solder-joint failure process involves a complex interplay of creep and fatigue mechanisms. To better understand the role of creep in typical multi-hour cyclic loading conditions, a specialized non-linear finite-element creep simulation computer program has been formulated. The numerical algorithm includes the complete part-lead-solder-PWB system, accounting for strain-rate dependence of creep on applied stress and temperature, and the role of the part-lead dimensions and flexibility that determine the total creep deflection (solder strain range) during stress relaxation. The computer program has been used to explore the effects of various solder creep-fatigue parameters such as lead height and stiffness, thermal-cycle test profile, and part/board differential thermal expansion properties. One of the most interesting findings is the strong presence of unidirectional creep-ratcheting that occurs during thermal cycling due to temperature dominated strain-rate effects. To corroborate the solder fatigue model predictions, a number of carefully controlled thermal-cycle tests have been conducted using special bimetallic test boards.

Ross, R. G., Jr.↗
NASA NTRSDate not supplied

Q-Law for Rapid Assessment of Low Thrust Cislunar Trajectories Via Automatic Differentiation

Q-Law is a Lyapunov-based control law used to determine optimal controls for a low thrust trajectory. One major issue with its use is the difficult derivatives re-quired for calculating optimal controls at a given time. In this paper, an implementation of Q-Law with automatic differentiation via a Python package called JAX is applied. With automatic differentiation, the difficult derivatives for Q-Law’s optimal controls are calculated with ease, and derivatives of final states with respect to Q-Law’s weights are found enabling gradient-based optimization of Q-Law for the first time. Different search and optimization methods for finding optimal weights are then compared using the LEO to GEO problem, and it was found that gradient-free methods like design of experiments and genetic algorithm produced the best results, but they took the longest time to get a solution, while the gradient-based method found a locally optimal result in a much faster time. Overall, the run time for a single propagation is manageable and well-suited for a mission designer to use as an initial guess generator for trajectory optimization, or for simple orbit transfer analysis.

Nathan Steffen↗
NASA NTRS2003

First Year Wilkinson Microwave Anisotropy Probe(WMAP) Observations: Data Processing Methods and Systematic Errors Limits

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight using the dipole modulation of the CMB due to the observatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter to the time-ordered data to remove a low level of l/f noise. We infer and correct for a small (approx. 1 %) transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel covariance matrix.

Hinshaw, G.↗