Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed data structures”

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 937 records · Page 52

Bridging the Gap between Cosmological Simulations with Graph Neural Networks and Domain Adaptation

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross dataset tasks (up to 28% better relative error and up to almost an order of magnitude better χ 2 ). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

97 MATHEMATICS AND COMPUTING↗

Potentially adaptive SARS-CoV-2 mutations discovered with novel spatiotemporal and explainable AI models

Abstract Background A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic. Results Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp 614 Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro 323 Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus. Conclusions These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

59 BASIC BIOLOGICAL SCIENCES↗

Interface properties and built-in potential profile of a LaCr O 3 / SrTi O 3 superlattice determined by standing-wave excited photoemission spectroscopy

LaCrO 3 (LCO)/SrTiO 3 (STO) heterojunctions are intriguing due to a polar discontinuity along [001], exhibiting two distinct and controllable charged interface structures [(LaO) + /(TiO 2 ) 0 and (SrO) 0 /(CrO 2 ) - ] with induced polarization, and a resulting depth-dependent potential. In this study, we have used soft- and hard-x-ray standing-wave excited photoemission spectroscopy (SW-XPS) to quantitatively determine the elemental depth profile, interface properties, and depth distribution of the polarization-induced built-in potentials. We observe an alternating charged interface configuration: a positively charged (LaO) + /(TiO 2 ) 0 intermediate layer at the LCO top /STO bottom interface and a negatively charged (SrO) 0 /(CrO 2 ) - intermediate layer at the STO top /LCO bottom interface. Using core-level SW data, we have determined the depth distribution of species, including through the interfaces, and these results are in excellent agreement with scanning transmission electron microscopy and electron energy-loss spectroscopy mapping of local structure and composition. SW-XPS also enabled deconvolution of the LCO and STO contributions to the valence-band (VB) spectra. Using a two-step analytical approach involving first SW-induced core-level binding-energy shifts and then VB modeling, the variation in potential across the complete superlattice is determined in detail. As a result, this potential is in excellent agreement with density functional theory models, confirming this method as a generally useful tool for interface studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Large Impact Basins on Mercury: Global Distribution, Characteristics, and Modification History from MESSENGER Orbital Data

The formation of large impact basins (diameter D greater than or equal to 300 km) was an important process in the early evolution of Mercury and influenced the planet's topography, stratigraphy, and crustal structure. We catalog and characterize this basin population on Mercury from global observations by the MESSENGER spacecraft, and we use the new data to evaluate basins suggested on the basis of the Mariner 10 flybys. Forty-two certain or probable impact basins are recognized a few additional basins that may have been degraded to the point of ambiguity are plausible on the basis of new data but are classified as uncertain. The spatial density of large basins (D greater than or equal to 500 km) on Mercury is lower than that on the Moon. Morphological characteristics of basins on Mercury suggest that on average they are more degraded than lunar basins. These observations are consistent with more efficient modification, degradation, and obliteration of the largest basins on Mercury than on the Moon. This distinction may be a result of differences in the basin formation process (producing fewer rings), greater relaxation of topography after basin formation (subduing relief), and/or higher rates of volcanism during the period of heavy bombardment on Mercury compared to the Moon (burying basin rings and interiors).

Fassett, Caleb I.↗

Domain Adaptive Graph Neural Networks for Constraining Cosmological Parameters Across Multiple Data Sets

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross-dataset tasks (up to $28\%$ better relative error and up to almost an order of magnitude better $\chi^2$). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

Vertical Distribution of Aersols and Water Vapor Using CRISM Limb Observations

Near-infrared spectra taken in a limb-viewing geometry by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) on-board the Mars Reconnaissance Orbiter (MRO) provide a useful tool for probing atmospheric structure. Specifically, the observed radiance as a function of wavelength and height above the limb allows the vertical distribution of both dust and ice aerosols to be retrieved. These data serve as an important supplement to the aerosol profiling provided by the MRO/MCS instrument allowing independent validation and giving additional information on particle physical and scattering properties through multi-wavelength studies. A total of at least ten CRISM limb observations have been taken so far covering a full Martian year. Each set of limb observations nominally contains about four dozen scans across the limb giving pole-to-pole coverage for two orbits at roughly 100 and 290 W longitude over the Tharsis and Syrtis/Hellas regions, respectively. At each longitude, limb scans are spaced roughly 10 degrees apart in latitude, with a vertical spatial resolution on the limb of roughly 800 m. Radiative transfer modeling is used to model the observations. We compute synthetic CRISM limb spectra using a discrete-ordinates radiative transfer code that accounts for multiple scattering from aerosols and accounts for spherical geometry of the limb observations by integrating the source functions along curved paths in that coordinate system. Retrieved are 14-point vertical profiles for dust and water ice aerosols with resolution of 0.4 scale heights between one and six scale heights above the surface. After the aerosol retrieval is completed, the abundances of C02 (or surface pressure) and H20 gas are retrieved by matching the depth of absorption bands at 2000 nm for carbon dioxide and at 2600 run for water vapor. In addition to the column abundance of water vapor, limited information on its vertical structure can also be retrieved depending on the signal available from aerosol scattering.

Smith, Michael D.↗

Field-Programmable Gate Array Computer in Structural Analysis: An Initial Exploration

This paper reports on an initial assessment of using a Field-Programmable Gate Array (FPGA) computational device as a new tool for solving structural mechanics problems. A FPGA is an assemblage of binary gates arranged in logical blocks that are interconnected via software in a manner dependent on the algorithm being implemented and can be reprogrammed thousands of times per second. In effect, this creates a computer specialized for the problem that automatically exploits all the potential for parallel computing intrinsic in an algorithm. This inherent parallelism is the most important feature of the FPGA computational environment. It is therefore important that if a problem offers a choice of different solution algorithms, an algorithm of a higher degree of inherent parallelism should be selected. It is found that in structural analysis, an 'analog computer' style of programming, which solves problems by direct simulation of the terms in the governing differential equations, yields a more favorable solution algorithm than current solution methods. This style of programming is facilitated by a 'drag-and-drop' graphic programming language that is supplied with the particular type of FPGA computer reported in this paper. Simple examples in structural dynamics and statics illustrate the solution approach used. The FPGA system also allows linear scalability in computing capability. As the problem grows, the number of FPGA chips can be increased with no loss of computing efficiency due to data flow or algorithmic latency that occurs when a single problem is distributed among many conventional processors that operate in parallel. This initial assessment finds the FPGA hardware and software to be in their infancy in regard to the user conveniences; however, they have enormous potential for shrinking the elapsed time of structural analysis solutions if programmed with algorithms that exhibit inherent parallelism and linear scalability. This potential warrants further development of FPGA-tailored algorithms for structural analysis.

Singleterry, Robert C., Jr.↗

Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations

This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it-an improved and generalized version of Bayesian Blocks [Scargle 1998]-that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piece- wise linear and piecewise exponential representations, multivariate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by [Arias-Castro, Donoho and Huo 2003]. In the spirit of Reproducible Research [Donoho et al. (2008)] all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.

signal detection↗

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗

Wettability Effects on Primary Drainage Mechanisms and NAPL Distribution: A Pore-Scale Study

The pore-scale processes governing water drainage behavior in porous media have implications for geoscience multiphase scenarios including carbon capture and storage, contaminant site remediation, oil recovery, and vadose zone processes. However, few studies report directly observed pore-scale water drainage phenomena in 3-D soils. This knowledge gap limits our ability to verify assumptions underlying existing models and develop optimal solutions. As such, this paper utilizes synchrotron X-ray microtomography to present an experimental pore-scale examination of nonaqueous phase liquid (NAPL)/water distribution along a primary drainage front as dense NAPL was injected upward into water wetting (WW) and intermediate wetting (IW) sand-packed columns. Pore-network structures were extracted from imaged data sets and mapped onto segmented NAPL/water data sets which allowed quantitative examinations of wettability impacts on (a) the extent to which NAPL fills individual pore bodies and (b) relationships between pore size and the phase occupying the pore, with both considered as a function of distance (and capillary pressure) relative to the NAPL front. These results revealed that several hypotheses treating IW sand similarly to WW sands are simplistic. IW systems exhibited a sequence of pore filling that deviated from traditional capillary pressure-based model predictions: NAPL invades smaller pores, while larger, adjacent pores are bypassed leaving multipore residual water ganglia. NAPL pore saturations were close to 1 and did not change with capillary pressure in IW systems. Overall, the results illustrate how a relatively small change in operative contact angle alters NAPL distribution during water drainage, with important implications for geoscience multiphase flow scenarios.

58 GEOSCIENCES↗

Solving Stochastic Inverse Problems for Property–Structure Linkages Using Data-Consistent Inversion and Machine Learning

Determining process–structure–property linkages is one of the key objectives in material science, and uncertainty quantification plays a critical role in understanding both process–structure and structure–property linkages. In this work, we seek to learn a distribution of microstructure parameters that are consistent in the sense that the forward propagation of this distribution through a crystal plasticity finite element model matches a target distribution on materials properties. This stochastic inversion formulation infers a distribution of acceptable/consistent microstructures, as opposed to a deterministic solution, which expands the range of feasible designs in a probabilistic manner. Furthermore, to solve this stochastic inverse problem, we employ a recently developed uncertainty quantification framework based on push-forward probability measures, which combines techniques from measure theory and Bayes’ rule to define a unique and numerically stable solution. This approach requires making an initial prediction using an initial guess for the distribution on model inputs and solving a stochastic forward problem. To reduce the computational burden in solving both stochastic forward and stochastic inverse problems, we combine this approach with a machine learning Bayesian regression model based on Gaussian processes and demonstrate the proposed methodology on two representative case studies in structure–property linkages.

36 MATERIALS SCIENCE↗

Quasi-brittle porous material: Simulated effect of stochastic air void structure on compressive strength

Highlights: • A simulation procedure to study the effect of porosity on the compressive strength of a quasi-brittle material is introduced • Based on the simulated results, the relationship between compressive strength and air void fraction (porosity) is presented. • A linear relationship was achieved between the cubic root of porosity and the simulated compressive strength. The effect of porosity comprised of spherical air voids on the compressive strength of quasi-brittle material was studied via simulations. The simulated porous structures were based on pore size distributions of two mortar samples measured by X-ray microtomography. While the simulation method set practical limits on the size of sample, the base of the statistics was established by simulating 128 small structures generated by sampling from pore structures of two mortars. By studying the application of the classical strength-porosity formulas to the simulated data, a new simple model was formed. A linear relationship was achieved between the cubic root of air void fraction (porosity) and the simulated compressive strength. The reasons for scattering of simulated strength around fitted trend remained unresolved in this study; no clear dependence on pore number or other distribution properties was observed. With the presented simulation approach, the dependence of compressive strength on porosity is achieved independently of disturbances that occur in experimental studies creating understanding of compressional behavior of low porosity materials.

36 MATERIALS SCIENCE↗

Composition and structure of the atmosphere of Venus

Some of the main questions regarding the composition, structure, and origin of the atmosphere of Venus are posed. These questions are (1) the distribution of the constituents of the lower atmosphere, (2) cloud composition, (3) the planet's surface and interior as revealed by atmospheric data, (4) the state property profiles and their variation over the planet, (5) the reason for the high temperatures of the lower atmosphere, (6) composition and temperature profiles of the upper atmosphere and location of the homopause, (7) spatial and temporal variations in the upper atmosphere, (8) the cause of the stability of CO2 - global circulation or local turbulence, (9) influence of neutral composition on the thermal structure, (10) response of upper atmosphere to change in solar EUV and solar wind, (11) the source and destination of the atmosphere, and (12) the location of Venus's water. The main parameters to be measured which will aid in resolving these problems are enumerated.

Hoffman, J. H.↗

Dual polarization VLBI observations of stellar binary systems at 5 GHz

The binary systems Sigma CrB, SZ Psc, HR 1099, UX Ari, HR 5110, Algol, and Cyg X-1 are characterized in terms of marcsec-scale radio structure on the basis of dual-polarization observations at 4.969-4.997 GHz, obtained with an international VLBI array comprising the 100-m Effelsberg, 36-m Haystack, 42-m Green Bank, 26-m Fort Davis, 130-m-equivalent VLA, and 40-m Owens Valley telescopes (system temperatures 40-70 K) at 22-14 h UT on July 26-27, 1983. The data are presented in graphs and tables and discussed in detail. Single Gaussian brightness distributions with theta(FWHM) = 0.5-2.0 marcsec and brightness temperatures (2-200) x 10 to the 8th K are determined for four sources; core-halo structures are found for UX Ari and Algol; and an upper limit of 50 microarcsec is established for the angular separation of the regions of left and right circular polarization. The gross features seen in RS CVn systems are explained in terms of a simple expanding-coronal-loop model of flare events.

Mutel, R. L.↗

X-Ray Scattering Data Processing and Analysis

Semicrystalline polymers have a degree of crystallinity that depends on chain structure and crystallization conditions. Thermomechanical properties are influenced by the degree of crystallinity and crystalline morphology. A common method used to understand crystalline content, crystal and amorphous layer thickness, and layer thickness distributions, is X-ray scattering (XS). This report uses XS to first investigate crystalline orientation in high density polyethylene (HDPE) samples. Then, a single orientation is analyzed via XS in a temperature study to investigate crystalline behavior under changing temperatures. Lastly, block copolymer filaments extruded at varying speeds are analyzed. The X-ray data are used to gain a baseline understanding of crystalline behavior in these unique filaments.

36 MATERIALS SCIENCE↗

Investigation of the impact of supplemental reflective films to improve micro-light climate within tomato plant canopy in solar greenhouses

The non-uniform growth and development of crops within Chinese Solar Greenhouses (CSG) is directly related to the micro-light climate within canopy. In practice, reflective films are used to improve micro-light climate within plant canopy by homogenizing light distribution and so increasing total plant light interception. However, as to our knowledge, the contributions to light distribution within canopy have not been investigated for passive reflector like reflective films. Field experiments dealing with light conditions and growth behavior over time, are complicated to carry out, time-consuming and hard to control, while however, accurate measurements of how reflective films influence the micro-light climate of canopy are an essential step to improve the growth conditions for any crop. Here, we propose a supplementary light strategy using reflective films to improve light distribution within plant canopy. Based on the example of CSG, a 3D greenhouse model including a detailed 3D tomato canopy structure was constructed to simulate the influence of supplementary reflective films to improve micro-light climate. Comparison of measured solar radiation intensity with predicted model data demonstrated that the model could precisely predict light radiation intensity over time with different time points and positions in the greenhouse. A series of reflective film configurations were investigated based on features analysis of light distribution in the tomato canopy on sunny days using the proposed model. The reflective film configuration scheme with the highest impact significantly improved the evenness of horizontal and vertical light distribution in tomato canopy. The strategy provided here can be used to configure reflective films that will enhance light conditions in CSG, which can be applied and extended in different scenarios.

Liu, Anhua↗

First-Principle Calculations of Structural, Elastic, and Electronic Properties of Intermetallic Rare Earth R{sub 2}Ni{sub 2}Pb (R = Ho, Lu, and Sm) Compounds

The structural, elastic, and electronic properties of rare earth intermetallic R{sub 2}Ni{sub 2}Pb (where R = Ho, Lu, and Sm) compounds were investigated with the density functional theory (DFT) calculations. The calculations are performed using the full potential-linearized augmented plane wave (FP-LAPW) method within the framework of local density approximation (LDA). The calculated values of the equilibrium lattice constants were in agreement with the available experimental values. The elastic constants (C{sub ij}) were also calculated to understand the mechanical properties and structural stability of the compounds. Furthermore, the density of states and the charge density distributions of the compounds were calculated to understand the nature of the bonding in the material. The calculated results are in accordance with the available data in the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗