Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “algorithms and data structure”

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 919 records · Page 51

Transfer-function-parameter estimation from frequency response data: A FORTRAN program

A FORTRAN computer program designed to fit a linear transfer function model to given frequency response magnitude and phase data is presented. A conjugate gradient search is used that minimizes the integral of the absolute value of the error squared between the model and the data. The search is constrained to insure model stability. A scaling of the model parameters by their own magnitude aids search convergence. Efficient computer algorithms result in a small and fast program suitable for a minicomputer. A sample problem with different model structures and parameter estimates is reported.

Seidel, R. C.↗

Food Quality and Phytoplankton Community Composition in San Francisco Bay using Imaging Spectroscopy Data from the California HyspIRI Airborne Campaign

The San Francisco Bay (SFB) is the largest estuary on the west coast of the United States. It is an important transition zone between marine, freshwater, and inland terrestrial watersheds. The SFB is an important region for the cycling of nutrients and pollutants and it supports nurseries of ecologically and commercially important fisheries, including some threatened species. Phytoplankton community structure influences food web dynamics, and the taxonomy of the phytoplankton may be more important in determining primary food quality than environmental factors. As such, estimating food quality from phytoplankton community composition can be a robust tool to understand trophic transfer of energy. Recent work explores phytoplankton food quality in SFB through the use of microscopy and phytoplankton chemotaxonomy to evaluate how changes in phytoplankton composition may have influenced the recent trophic collapse of pelagic fishes in the northern part of the SFB. The objective of this study is to determine if the approach can also be applied to imaging spectroscopy data in order to quantify phytoplankton food quality from space. Imaging spectroscopy data of SFB from the Airborne VisibleInfrared Imaging Spectrometer (AVIRIS) was collected during the Hyperspectral Infrared (HyspIRI) Airborne Campaign in California (2013 2015) and used in this study. Estimates of ocean chlorophyll and phytoplankton community structure were determined using standard ocean chlorophyll algorithms and the PHYtoplankton Detection with Optics (PHYDOTax) algorithms. These were validated using in situ observations of phytoplankton composition using microscopic cell counts and phytoplankton chemotaxonomy from the US Geological Surveys ship surveys of the SFB. The findings from this study may inform the use of future high spectral resolution satellite sensors with the spatial resolution appropriate for coastal systems (e.g., HyspIRI) to assess food quality from space.

HyspIRI↗

ICESat-2 Photon Classification: Finding Signal Photons in the ATL03 Geolocated Photon Data Product

ICESat-2 carries NASA's next-generation laser altimeter, ATLAS, (Advanced Topographic Laser Altimeter System), designed to measure changes in ice sheet height, sea ice freeboard, and vegetation canopy height. ATLAS contains a photon-counting lidar which transmits green (532-nm) pulses at 10kHz. Each pulse is split into 3 pairs of beams (one strong and one weak). Approximately 1014 photons per pulse travel from ATLAS through the atmosphere to reflect off the Earth's surface. Some return back into the ATLAS telescope where they are recorded. Photons from sunlight and instrument noise at the same wavelength are also recorded. The flight software time tags all photons within a 500m to 6 km range window and generates histograms. Using the histograms, it selects a telemetry window which varies from 20m over flat surfaces to hundreds of meters over rougher terrain. ATL03 contains the time, height (relative to the WGS-84 ellipsoid), latitude and longitude of every photon within the telemetry window. The basic challenge is to determine which of these photons were reflected off the surface. We have developed an algorithm that identifies these signal photons and assigns a confidence level (low, medium, or high) to each signal photon based on the signal to noise ratio. We present an overview of the signal identification algorithm and show the results on actual ICESat-2 data over ice sheet, sea ice, vegetated, and water surfaces. Higher level ATLAS products work with aggregations of the photons in order to determine the ellipsoidal height of the Earth, canopy height and structure, and other quantities of geophysical interest.

Remote Sensing↗

Understanding the Global Three-Dimensional Distribution of Precipitation Mean Particle Size with the Global Precipitation Measurement Mission

This study addresses the global distribution of precipitation mean particle size using data from the Global Precipitation Measurement (GPM) mission. The mass-weighted mean diameter Dm is a characteristic parameter of the precipitation particle size distribution (PSD), estimated from the GPM Combined Radar–Radiometer Algorithm (CORRA) using data from GPM’s dual-frequency precipitation radar and microwave imager. We examine Dm in individual precipitation systems in different climate regimes and investigate a 6-yr (2014–20) global climatology within 70°N–70°S. The vertical structure of Dm is demonstrated with cases of deep convection, frontal rain and snow, and stratocumulus light rain. The Dm values, detectable by GPM, range from ~0.7 mm in stratocumulus precipitation to >3.5 mm in the ice layers of intense convection. Within the constraint of the 12-dBZ detectability threshold, the smallest annual mean Dm (~0.8 mm) are found in the eastern oceans, and the largest values (~2 mm) occur above the melting levels in convection over land in summer. The standard deviation of the annual mean is generally <0.45 mm below 6 km. Climate regimes are characterized with Dm annual/seasonal variations, its convective/stratiform components, and vertical variabilities (2–10 km). The U.S. Central Plains and Argentina are associated with the largest Dm in a deep layer. Tropical Africa has larger Dm and standard deviation than Amazon. Large convective Dm occurs at high latitudes of Eurasia and North America in summer; the Southern Hemisphere high latitudes have shallower systems with smaller Dm. Oceanic storm tracks in both hemispheres have relatively large Dm, particularly for convective Dm in winter. Relatively small Dm occurs over tropical oceans, including ITCZ, requiring further investigation.

Mei Han↗

Architecture and data processing alternatives for the TSE computer. Volume 3: Execution of a parallel counting algorithm using array logic (Tse) devices

A parallel algorithm for counting the number of logic-l elements in a binary array or image developed during preliminary investigation of the Tse concept is described. The counting algorithm is implemented using a basic combinational structure. Modifications which improve the efficiency of the basic structure are also presented. A programmable Tse computer structure is proposed, along with a hardware control unit, Tse instruction set, and software program for execution of the counting algorithm. Finally, a comparison is made between the different structures in terms of their more important characteristics.

Metcalfe, A. G.↗

ED-cPSD: Fast Phase-Size Distribution via Sequential Erosion-Dilation

The Erosion-Dilation continuous Phase-Size Distribution, ED-cPSD, is an application for calculating continuous pore and particle-size distribution from digital reconstructions and/or image-based structural data. It is based on the erosion-dilation continuous phase-size distribution method. A continuous size distribution is a measure of the probability density of finding a particle or pore of a certain size. These distributions are of interest in any field of study involving porous media, including but not limited to electrochemistry, petroleum engineering, geology, and food science. The algorithm behind the software provides a computationally efficient way to calculate phase-size distributions for large domains. For a 3D battery electrode reconstruction with 1.3 x 10 8 voxels, the particle size distribution is derived in under 2 min on a desktop, while also retaining flexibility and computational efficiency for HPC-scale multi-threading. The software can handle structures with over 10 9 voxels. The algorithm is roughly 280 times faster than a previous version on the same task.

Characterization↗

Extraction of mechanical properties of materials through deep learning from instrumented indentation

Instrumented indentation has been developed and widely utilized as one of the most versatile and practical means of extracting mechanical properties of materials. This method is particularly desirable for those applications where it is difficult to experimentally determine the mechanical properties using stress–strain data obtained from coupon specimens. Such applications include material processing and manufacturing of small and large engineering components and structures involving the following: three-dimensional (3D) printing, thin-film and multilayered structures, and integrated manufacturing of materials for coupled mechanical and functional properties. Here, we utilize the latest developments in neural networks, including a multifidelity approach whereby deep-learning algorithms are trained to extract elastoplastic properties of metals and alloys from instrumented indentation results using multiple datasets for desired levels of improved accuracy. We have established algorithms for solving inverse problems by recourse to single, dual, and multiple indentation and demonstrate that these algorithms significantly outperform traditional brute force computations and function-fitting methods. Moreover, we present several multifidelity approaches specifically for solving the inverse indentation problem which 1) significantly reduce the number of high-fidelity datasets required to achieve a given level of accuracy, 2) utilize known physical and scaling laws to improve training efficiency and accuracy, and 3) integrate simulation and experimental data for training disparate datasets to learn and minimize systematic errors. The predictive capabilities and advantages of these multifidelity methods have been assessed by direct comparisons with experimental results for indentation for different commercial alloys, including two wrought aluminum alloys and several 3D printed titanium alloys.

36 MATERIALS SCIENCE↗

Multitiered computational methodology for extracting three-dimensional rotational diffusion coefficients from x-ray photon correlation spectroscopy data without structural information

X-ray photon correlation spectroscopy (XPCS) is a powerful technique for analyzing particle systems by investigating their dynamics in suspensions across a broad range of temporal and spatial scales. This is done by illuminating samples with coherent x-ray beams and calculating the correlation function of the obtained x-ray scattering images. XPCS is uniquely suited for studying Brownian dynamics, consisting of translational and rotational diffusion. While traditional XPCS image analysis techniques can extract translational diffusion components, they are unable to estimate rotational diffusion coefficients. Here, we introduce a methodology that combines the angular-temporal cross-correlation analysis and a algorithmic framework called Multi-Tiered Estimation for Correlation Spectroscopy in 3D for estimating three-dimensional rotational diffusion coefficients from XPCS images of three-dimensional particle systems. We demonstrate our methodology for extracting rotational diffusion coefficients from XPCS data by applying it to simulated noisy x-ray images of systems of crossing nanotubes and proteins that evolve under translational and rotational Brownian motion for different diffusion rates. Furthermore, our results show that our approach determines rotational diffusion coefficients within a few percent error.

97 MATHEMATICS AND COMPUTING↗

Geometric Interpretation of the Cluster Location Problem Part II: Application to the Pahala, Hawaii, Earthquake Sequence

In the companion “Theory” article, we presented a new framing of the seismic location problem in terms of differential geometry (Harris et al., 2025). From that viewpoint, we developed a “project and correct” approach for estimating the relative locations of earthquakes. Here, in this study, we use project and correct to estimate high-precision relative locations of events from an earthquake sequence beneath the town of Pahala, Hawaii, using high-precision correlation-derived picks. The sequence was active from 2020 through 2022 and produced many highly correlated signals at Hawaii Volcano Observatory (HVO) stations on the island of Hawaii. The data we inverted consisted of 2882 events with observations at 5 HVO stations. For comparison with the travel-time image, we also produced conventional hypocenter solutions using both the Bayesloc program (Myers et al., 2007, 2009) and a purpose-built double-difference code. There were obvious structural elements in the resulting image, the resolution of which we used to test the performance of the project and the correct algorithm. For the projection step, we first produced a 3D local basis using an singular value decomposition (SVD) of the 2882 groups of times. Projection of the travel-time vectors into this basis resulted in an image with structures similar to those produced by our conventional locators, but with distortion as predicted by theory. Removing the distortion requires an inverse operator generated from the metric tensor at the geometric centroid of the events. We compared two approaches to obtaining such an inverse operator. The first uses an estimate of the geographic centroid of the event cloud from the centroid of the travel-time data. The second approach uses the centroid of the conventionally produced locations. The first approach produces a corrected image very similar to the conventional results, but with a rotation. The corrected image produced using the conventionally derived centroid is a near-exact match to the conventional locations.

Dodge, Douglas A. [Lawrence Livermore National Lab↗

Discovering causal structure with reproducing-kernel Hilbert space ε -machines

We merge computational mechanics’ definition of causal states (predictively equivalent histories) with reproducing-kernel Hilbert space (RKHS) representation inference. The result is a widely applicable method that infers causal structure directly from observations of a system’s behaviors whether they are over discrete or continuous events or time. A structural representation—a finite- or infinite-state kernel ϵ-machine—is extracted by a reduced-dimension transform that gives an efficient representation of causal states and their topology. In this way, the system dynamics are represented by a stochastic (ordinary or partial) differential equation that acts on causal states. We introduce an algorithm to estimate the associated evolution operator. Paralleling the Fokker–Planck equation, it efficiently evolves causal-state distributions and makes predictions in the original data space via an RKHS functional mapping. We demonstrate these techniques, together with their predictive abilities, on discrete-time, discrete-value infinite Markov-order processes generated by finite-state hidden Markov models with (i) finite or (ii) uncountably infinite causal states and (iii) continuous-time, continuous-value processes generated by thermally driven chaotic flows. The method robustly estimates causal structure in the presence of varying external and measurement noise levels and for very high-dimensional data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Machine learning interatomic potential for silicon-nitride (Si 3 N 4 ) by active learning

Silicon nitride (Si 3 N 4 ) is an extensively used material in the automotive, aerospace, and semiconductor industries. However, its widespread use is in contrast to the scarce availability of reliable interatomic potentials that can be employed to study various aspects of this material on an atomistic scale, particularly its amorphous phase. In this work, we developed a machine learning interatomic potential, using an efficient active learning technique, combined with the Gaussian approximation potential (GAP) method. Our strategy is based on using an inexpensive empirical potential to generate an initial dataset of atomic configurations, for which energies and forces were recalculated with density functional theory (DFT); thereafter, a GAP was trained on these data and an iterative re-training algorithm was used to improve it by learning on-the-fly. When compared to DFT, our potential yielded a mean absolute error of 8 meV/atom in energy calculations for a variety of liquid and amorphous structures and a speed-up of molecular dynamics simulations by 3–4 orders of magnitude, while achieving a first-rate agreement with experimental results. Our potential is publicly available in an open-access repository.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis and compensation of an aircraft simulator control loading system with compliant linkage

A hydraulic control loading system for aircraft simulation was analyzed to find the causes of undesirable low frequency oscillations and loading effects in the output. The hypothesis of mechanical compliance in the control linkage was substantiated by comparing the behavior of a mathematical model of the system with previously obtained experimental data. A compensation scheme based on the minimum integral of the squared difference between desired and actual output was shown to be effective in reducing the undesirable output effects. The structure of the proposed compensation was computed by use of a dynamic programing algorithm and a linear state space model of the fixed elements in the system.

Johnson, P. R.↗

Preliminary data processing system for telemetry data from Interkosmos-8 type satellites

The preliminary data processing system (PDPS), a set of algorithms, programs, equipment, data storage and memory devices, and all devices controlling the processing, is described using the Intercosmos-8 satellite as an example. Each stage of the system is considered along with its content and sequence. The essential elements of the data processing structure are described and the concepts underlying the selection of these elements are discussed.

Vilchinskaya, A. S.↗

Geologic applications of thermal-inertia mapping from satellite

In the Powder River Basin, Wyo., narrow geologic units having thermal inertias which contrast with their surroundings can be discriminated in optimal images. A few subtle thermal inertia anomalies coincide with areas of helium leakage believed to be associated with deep oil and gas concentrations. The most important results involved delineation of tectonic framework elements some of which were not previously recognized. Thermal and thermal inertia images also permit mapping of geomorphic textural domains. A thermal lineament appears to reveal a basement discontinuity which involves the Homestake Mine in the Black Hill, a zone of Tertiary igneous activity and facies control in oil producing horizons. Applications of these data to the Cabeza Prieta, Ariz., area illustrate their potential for igneous rock type discrimination. Extension to Yellowstone National Park resulted in the detection of additional structural information but surface hydrothermal features could not be distinguished with any confidence. A thermal inertia mapping algorithm, a fast and accurate image registration technique, and an efficient topographic slope and elevation correction method were developed.

Offield, T. W.↗

Interaction of a cumulus cloud ensemble with the large-scale environment. IV - The discrete model

The Arakawa-Schubert (1974) parameterization is applied to a prognostic model of large-scale atmospheric circulations and used to analyze data in a general circulation model (GCM). The vertical structure of the large-scale model and the solution for the cloud subensemble thermodynamical properties are examined to choose cloud levels and representative regions. A mass flux distribution equation is adapted to formulate algorithms for calculating the large-scale forcing and the mass flux kernel, using either direct solution or linear programming. Finally, the feedback of the cumulus ensemble on the large-scale environment for a given subensemble mass flux is calculated. All cloud subensemble properties were determined from the conservation of mass, moist static energy, and total water.

Lord, S. J.↗

Fidelity of Satellite Soundings

Two pattern recognition procedures to provide improvements to first guess fields for satellite temperature retrievals were developed. The first technique where a radiometer measurement is used to select one or more historical radiosonde temperature profiles as analog estimates of ambient thermal structure. Radiosondes are the vertical scales of the analogs. The vertical resolving power of the satellite radiometer is relevant only to a decision process. The analog selection process is much more effective if implemented in an orthogonalized space of measurement information. The second procedure partitions a priori dependent data into shape coherent pattern libraries using structure information inherent in the data itself. It is shown that the pattern recognition techniques are able to reduce first guess profile errors by nearly 50%. The impact of pattern recognition on temperature retrieval error is assessed by regression and physical iterative retrieval algorithms. The influence of improved first guess fields is markedly different on these two types of algorithm. Pattern recognition was a strong positive impact on the physical iterative method.

Thompson, O.↗

Evolution and advanced technology

The NASREM architecture with its standard interfaces permits development and evolution of the Flight Telerobotic Servicer to greater autonomy. Technologies in control strategies for an arm with seven DOF, including a safety system containing skin sensors for obstacle avoidance, are being developed. Planning and robotic execution software includes symbolic task planning, world model data bases, and path planning algorithms. Research over the last five years has led to the development of laser scanning and ranging systems, which use coherent semiconductor laser diodes for short range sensing. The possibility of using a robot to autonomously assemble space structures is being investigated. A control framework compatible with NASREM is being developed that allows direct global control of the manipulator. Researchers are developing systems that permit an operator to quickly reconfigure the telerobot to do new tasks safely.

Ollendorf, Stanford↗

Model correlation and damage location for large space truss structures: Secant method development and evaluation

On-orbit testing of a large space structure will be required to complete the certification of any mathematical model for the structure dynamic response. The process of establishing a mathematical model that matches measured structure response is referred to as model correlation. Most model correlation approaches have an identification technique to determine structural characteristics from the measurements of the structure response. This problem is approached with one particular class of identification techniques - matrix adjustment methods - which use measured data to produce an optimal update of the structure property matrix, often the stiffness matrix. New methods were developed for identification to handle problems of the size and complexity expected for large space structures. Further development and refinement of these secant-method identification algorithms were undertaken. Also, evaluation of these techniques is an approach for model correlation and damage location was initiated.

Smith, Suzanne Weaver↗