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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Thermal Analysis of Unusual Local-scale Features on the Surface of Vesta

At 525 km in mean diameter, Vesta is the second-most massive object in the main asteroid belt of our Solar System. At all scales, pyroxene absorptions are the most prominent spectral features on Vesta and overall, Vesta mineralogy indicates a complex magmatic evolution that led to a differentiated crust and mantle [1]. The thermal behavior of areas of unusual albedo seen on the surface at the local scale can be related to physical properties that can provide information about the origin of those materials. Dawn's Visible and Infrared Mapping Spectrometer (VIR) [2] hyperspectral images are routinely used, by means of temperature-retrieval algorithms, to compute surface temperatures along with spectral emissivities. Here we present temperature maps of several local-scale features of Vesta that were observed by Dawn under different illumination conditions and different local solar times.

Tosi, F.↗
NASA NTRS1977

Impact on multilayered composite plates

Stress wave propagation in a multilayer composite plate due to impact was examined by means of the anisotropic elasticity theory. The plate was modelled as a number of identical anisotropic layers and the approximate plate theory of Mindlin was then applied to each layer to obtain a set of difference-differential equations of motion. Dispersion relations for harmonic waves and correction factors were found. The governing equations were reduced to difference equations via integral transforms. With given impact boundary conditions these equations were solved for an arbitrary number of layers in the plate and the transient propagation of waves was calculated by means of a Fast Fourier Transform algorithm. The multilayered plate problem was extended to examine the effect of damping layers present between two elastic layers. A reduction of the interlaminar normal stress was significant when the thickness of damping layer was increased but the effect was mostly due to the softness of the damping layer. Finally, the problem of a composite plate with a crack on the interlaminar boundary was formulated.

Kim, B. S.↗
NASA NTRSDate not supplied

LIPA: Lunar Ice Perception Algorithm

Introduction: The highest concentration of Lu-nar water-ice stores exists within the Permanently Shadowed Regions (PSRs) of the Lunar South Pole [1-3]. As such, the ability to locate in situ water-ice stores in an accurate, systematic, and safe manner will prove vital for future Lunar activities which rely on hydrogen-based resources. Here we show how the strong absorptive properties of ice can be exploited so that surface ice located in PSRs can be easily differentiated from the surrounding frozen regolith. Testbeds which simulate an icy lunar landscape were created and then systematically imaged using a mid-wave infrared (MWIR) camera system. Testbeds were imaged under two filter modes (1) high-absorption (high-abs) mode: whereby imagery captured were confined to a single central wave-length (CWL) of 3.15 ± 0.03 μm and (2) low-absorption (low-abs) mode: whereby imagery cap-tured were confined to a single CWL of 3.80 ± 0.04 μm (Figure 1). High- and low-absorption modes are related to the absorptive properties of ice at each selected wavelength, respectively. Corresponding images from each filter mode were differenced (i.e., pixels were subtracted) to enhance contrast between ice-bearing and non-ice-bearing pixels, and then fed into a semantic segmen-tation model. The model was trained to detect and differentiate between water, ice, shadows, and lunar regolith. Results: Modeling results accurately discrimi-nated ice from other materials (such as frozen lunar regolith) and were used to visually resolve the spa-tial extent of surface ice. Further, outputs produced through semantic segmentation were used to estimate water-ice contents in collected imagery [(Pixels with Class = “Water Ice”)/(Sum of Pixels)*100]. Summary: These works prove promising for future in situ resource utilization (ISRU) missions which employ robotics in combination with infrared camera systems to advance science objectives (e.g., locate water-ice in frozen regolith) on the lunar sur-face. References: [1] Cannon K. M., Deutsch A. N., Head J. W., and Britt D. T. (2020) Geophysical Re-search Letters, 46, e2020GL088920. [2] Honniball C. I. et al. (2021) Nature Astronomy 5, no. 2, 121-127. [3] Li S. et al. (2018) Proceedings of the National Academy of Sciences, 115(36), 8907-8912.

A. Slabic↗
NASA NTRS1993

An algorithm for solving the perturbed gas dynamic equations

The present application of a compact, higher-order central-difference approximation to the linearized Euler equations illustrates the multimodal character of these equations by means of computations for acoustic, vortical, and entropy waves. Such dissipationless central-difference methods are shown to propagate waves exhibiting excellent phase and amplitude resolution on the basis of relatively large time-steps; they can be applied to wave problems governed by systems of first-order partial differential equations.

Davis, Sanford↗
NASA NTRS1991

A Large-Grain Mapping Approach for Multiprocessor Systems Through Data Flow Model Ph.D. Thesis

A large-grain level mapping method is presented of numerical oriented applications onto multiprocessor systems. The method is based on the large-grain data flow representation of the input application and it assumes a general interconnection topology of the multiprocessor system. The large-grain data flow model was used because such representation best exhibits inherited parallelism in many important applications, e.g., CFD models based on partial differential equations can be presented in large-grain data flow format, very effectively. A generalized interconnection topology of the multiprocessor architecture is considered, including such architectural issues as interprocessor communication cost, with the aim to identify the 'best matching' between the application and the multiprocessor structure. The objective is to minimize the total execution time of the input algorithm running on the target system. The mapping strategy consists of the following: (1) large-grain data flow graph generation from the input application using compilation techniques; (2) data flow graph partitioning into basic computation blocks; and (3) physical mapping onto the target multiprocessor using a priority allocation scheme for the computation blocks.

Kim, Hwa-Soo↗
NASA NTRS1993

Automatic partitioning of unstructured meshes for the parallel solution of problems in computational mechanics

Most of the recently proposed computational methods for solving partial differential equations on multiprocessor architectures stem from the 'divide and conquer' paradigm and involve some form of domain decomposition. For those methods which also require grids of points or patches of elements, it is often necessary to explicitly partition the underlying mesh, especially when working with local memory parallel processors. In this paper, a family of cost-effective algorithms for the automatic partitioning of arbitrary two- and three-dimensional finite element and finite difference meshes is presented and discussed in view of a domain decomposed solution procedure and parallel processing. The influence of the algorithmic aspects of a solution method (implicit/explicit computations), and the architectural specifics of a multiprocessor (SIMD/MIMD, startup/transmission time), on the design of a mesh partitioning algorithm are discussed. The impact of the partitioning strategy on load balancing, operation count, operator conditioning, rate of convergence and processor mapping is also addressed. Finally, the proposed mesh decomposition algorithms are demonstrated with realistic examples of finite element, finite volume, and finite difference meshes associated with the parallel solution of solid and fluid mechanics problems on the iPSC/2 and iPSC/860 multiprocessors.

Farhat, Charbel↗
NASA NTRS1975

Simulation of multicorrelated random processes using the FFT algorithm

A technique for the digital simulation of multicorrelated Gaussian random processes is described. This technique is based upon generating discrete frequency functions which correspond to the Fourier transform of the desired random processes, and then using the fast Fourier transform (FFT) algorithm to obtain the actual random processes. The main advantage of this method of simulation over other methods is computation time; it appears to be more than an order of magnitude faster than present methods of simulation. One of the main uses of multicorrelated simulated random processes is in solving nonlinear random vibration problems by numerical integration of the governing differential equations. The response of a nonlinear string to a distributed noise input is presented as an example.

Wittig, L. E.↗
NASA NTRS1983

Advanced methods for preparation and characterization of infrared detector materials

Crystals were prepared by the Bridgman-Stockbarger method with a wide range of crystal growth rates and temperature gradients adequate to prevent constitutional supercooling under diffusion-limited, steady-state, growth conditions. The longitudinal compositional gradients for different growth conditions and alloy compositions were calculated and compared with experimental data to develop a quantitative model of solute redistribution during the crystal growth of the alloys. Measurements were performed to ascertain the effect of growth conditions on radial compositional gradients. The pseudobinary HgTe-CdTe constitutional phase diagram was determined by precision differential-thermal-analysis measurements and used to calculate the segregation coefficient of Cd as a function of x and interface temperature. Experiments were conducted to determine the ternary phase equilibria in selected regions of the Hg-Cd-Te constitutional phase diagram. Electron and hole mobilities as functions of temperature were analyzed to establish charge-carrier scattering probabilities. Computer algorithms specific to Hg(1-x)CdxTe were developed for calculations of the charge-carrier concentration, charge-carrier mobilities, Hall coefficient, and Dermi Fermi energy as functions of x, temperature, ionized donor and acceptor concentrations, and neutral defect concentrations.

Broerman, J. G.↗
NASA NTRS1986

Constitutive modeling of cyclic plasticity and creep, using an internal time concept

Using the concept of an internal time as related to plastic strains, a differential stress-strain relation for elastoplasticity is rederived, such that (1) the concept of a yield-surface is retained; (2) the definitions of elastic and plastic processes are analogous to those in classical plasticity theory; and (3) its computational implementation, via a 'tangent-stiffness' finite element method and a 'generalized-midpoint-radial-return' stress-integration algorithm, is simple and efficient. Also, using the concept of an internal time, as related to both the inelastic strains as well as the Newtonian time, a constitutive model for creep-plasticity interaction, is discussed. The problem of modeling experimental data for plasticity and creep, by the present analytical relations, as accurately as desired, is discussed. Numerical examples which illustrate the validity of the present relations are presented for the cases of cyclic plasticity and creep.

Watanabe, O.↗
NASA NTRS1987

Numerical studies of the thermal design sensitivity calculation for a reaction-diffusion system with discontinuous derivatives

The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.

Hou, Jean W.↗
NASA NTRS1992

Using adaptive grid in modeling rocket nozzle flow

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which cannot be solved analytically. However, this system of equations called the Navier-Stokes equations can be solved numerically. The accuracy and the convergence of the solution of the system of equations will depend largely on how precisely the sharp gradients in the domain of interest can be resolved. With the advances in computer technology, more sophisticated algorithms are available to improve the accuracy and convergence of the solutions. An adaptive grid generation is one of the schemes which can be incorporated into the algorithm to enhance the capability of numerical modeling. It is equivalent to putting intelligence into the algorithm to optimize the use of computer memory. With this scheme, the finite difference domain of the flow field called the grid does neither have to be very fine nor strategically placed at the location of sharp gradients. The grid is self adapting as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzles by taking the refinement part of grid generation out of the hands of computational fluid dynamics (CFD) specialists and place it into the computer algorithm itself.

Chow, Alan S.↗
NASA NTRSDate not supplied

LIPA: Lunar Ice Perception Algorithm

The highest concentration of Lunar water-ice stores exists within the Permanently Shadowed Regions (PSRs) of the Lunar South Pole. As such, the ability to locate in situ water-ice stores in an accurate, systematic, and safe manner will prove vital for future Lunar activities which rely on hydrogen-based resources. Here we show how the strong absorptive properties of ice can be exploited (by coupling robotics, infrared imaging techniques, and machine learning) so that surface frost located in PSRs can be easily differentiated from the surrounding frozen regolith. Testbeds which simulate an icy lunar landscape were created and then imaged using a mid-wave infrared (MWIR) camera system mounted to a robotic arm (UR5e). Testbeds were imaged under two filter modes (1) wide band mode: whereby imagery captured filled a spectral range of 3.0 - 5.0 μm and (2) narrow band mode: whereby imagery captured were confined to a single central wavelength (CWL) of 3.15 ± 0.03 μm. A CWL of 3.15 μm was chosen due to the highly absorptive nature of ice at that specific wavelength. Images produced under both camera modes were processed in MATLAB. Narrow band images (NB) were subtracted from their wide band (WB) counterparts to produce differenced images (DI) which clearly demonstrated the spatial extent of ice (e.g., WB – NB = DI). Differenced images were used to train a Microsoft Azure model to discriminate between frozen regolith which did and did not contain ice. These works prove promising for future in situ resource utilization (ISRU) missions which employ robotics in combination with camera systems to advance science objectives (e.g., locate water-ice in frozen regolith) on the lunar surface.

Ane Slabic↗
NASA NTRS1986

Tectonics and evolution of the Juan Fernandez microplate at the Pacific-Nazca-Antarctic triple junction

Magnetic and bathymetric profiles derived from the R/V Endeavor survey and focal mechanism studies for earthquakes on two of the Juan Fernandez microplate boundaries are analyzed. It is observed that the Nazca-Juan Fernandez pole is in the northern end of the microplate since the magnetic lineation along the East Ridge of the microplate fans to the south. The calculation of the relative motion of the Juan Fernandez-Pacific-Nazca-Antarctic four-plate system using the algorithm of Minster et al. (1974) is described. The development of tectonic and evolutionary models of the region is examined. The tectonic model reveals that the northern boundary of the Juan Fernandez microplate is a zone of compression and that the West Ridge and southwestern boundary are spreading obliquely; the evolutionary model relates the formation of the Juan Fernandez microplate to differential spreading rates at the triple junction.

Anderson-Fontana, S.↗
NASA NTRS1996

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗
NASA NTRS2014

A New Hybrid Spatio-temporal Model for Estimating Daily Multi-year PM2.5 Concentrations Across Northeastern USA Using High Resolution Aerosol Optical Depth Data

The use of satellite-based aerosol optical depth (AOD) to estimate fine particulate matter PM(sub 2.5) for epidemiology studies has increased substantially over the past few years. These recent studies often report moderate predictive power, which can generate downward bias in effect estimates. In addition, AOD measurements have only moderate spatial resolution, and have substantial missing data. We make use of recent advances in MODIS satellite data processing algorithms (Multi-Angle Implementation of Atmospheric Correction (MAIAC), which allow us to use 1 km (versus currently available 10 km) resolution AOD data.We developed and cross validated models to predict daily PM(sub 2.5) at a 1X 1 km resolution across the northeastern USA (New England, New York and New Jersey) for the years 2003-2011, allowing us to better differentiate daily and long term exposure between urban, suburban, and rural areas. Additionally, we developed an approach that allows us to generate daily high-resolution 200 m localized predictions representing deviations from the area 1 X 1 km grid predictions. We used mixed models regressing PM(sub 2.5) measurements against day-specific random intercepts, and fixed and random AOD and temperature slopes. We then use generalized additive mixed models with spatial smoothing to generate grid cell predictions when AOD was missing. Finally, to get 200 m localized predictions, we regressed the residuals from the final model for each monitor against the local spatial and temporal variables at each monitoring site. Our model performance was excellent (mean out-of-sample R(sup 2) = 0.88). The spatial and temporal components of the out-of-sample results also presented very good fits to the withheld data (R(sup 2) = 0.87, R(sup)2 = 0.87). In addition, our results revealed very little bias in the predicted concentrations (Slope of predictions versus withheld observations = 0.99). Our daily model results show high predictive accuracy at high spatial resolutions and will be useful in reconstructing exposure histories for epidemiological studies across this region.

Air pollution↗
NASA NTRS2012

Optimal Force Control of Vibro-Impact Systems for Autonomous Drilling Applications

The need to maintain optimal energy efficiency is critical during the drilling operations performed on future and current planetary rover missions (see figure). Specifically, this innovation seeks to solve the following problem. Given a spring-loaded percussive drill driven by a voice-coil motor, one needs to determine the optimal input voltage waveform (periodic function) and the optimal hammering period that minimizes the dissipated energy, while ensuring that the hammer-to-rock impacts are made with sufficient (user-defined) impact velocity (or impact energy). To solve this problem, it was first observed that when voice-coil-actuated percussive drills are driven at high power, it is of paramount importance to ensure that the electrical current of the device remains in phase with the velocity of the hammer. Otherwise, negative work is performed and the drill experiences a loss of performance (i.e., reduced impact energy) and an increase in Joule heating (i.e., reduction in energy efficiency). This observation has motivated many drilling products to incorporate the standard bang-bang control approach for driving their percussive drills. However, the bang-bang control approach is significantly less efficient than the optimal energy-efficient control approach solved herein. To obtain this solution, the standard tools of classical optimal control theory were applied. It is worth noting that these tools inherently require the solution of a two-point boundary value problem (TPBVP), i.e., a system of differential equations where half the equations have unknown boundary conditions. Typically, the TPBVP is impossible to solve analytically for high-dimensional dynamic systems. However, for the case of the spring-loaded vibro-impactor, this approach yields the exact optimal control solution as the sum of four analytic functions whose coefficients are determined using a simple, easy-to-implement algorithm. Once the optimal control waveform is determined, it can be used optimally in the context of both open-loop and closed-loop control modes (using standard realtime control hardware).

Aldrich, Jack B.↗
NASA NTRS1984

Numerical modeling of turbulent flow

Three dimensional combustor calculations are currently stretching the computer hardware capabilities and the computing budgets of gas turbine manufacturers. One of the main reasons for this relates to the large number of complex physical processes occurring in the combustor. Airflow, fuel spray, reaction kinetics, flame radiation, and not the least of which, turbulence must be modeled and the related differential equations solved. Discussions in this conference will address methods to improve the accuracy of combustor flow field calculations and methods to speed the convergence of the modeled equations. This report will focus on aspects of merging these two new technologies. The improved accuracy discretization schemes have a negative impact on the speed of convergence of the modeled equations that the improved solution algorithms may not overcome. A description of the causes of this problem and potential solutions will be examined.

Claus, R. W.↗
NASA NTRS1987

Analysis of SO2 signals in SBUV/TOMS data

Absorption bands between 300 nm and 315 nm were observed in spectral scans of the atmospheric albedo made by the Solar Backscatter Ultraviolet (SBUV) spectrometer following the eruption of El Chichon. It is shown that those bands coincide with peaks in the absorption coefficient spectrum of SO sub 2. The magnitude of the absorption is used to the column content of present SO sub 2. These observations confirm that the differential absorption between 312.5 nm and 317.5 nm in Total Ozone Mapping Spectrometer (TOMS) data can be used to measure SO sub 2 with high spatial resolution. A maximum concentration of SO sub 2 of 15 matm-cm was observed by SBUV on April 15, 1982. The measurement based on direct measurement of band intensity can be used to calibrate the TOMS algorithm for deriving SO sub 2 amounts.

Mcpeters, Richard↗