Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel 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 415 records · Page 23

Application of the NASA Multiscale Analysis Tool: Multiscale Integration and Interoperability

The NASA Multiscale Analysis Tool (NASMAT) was developed recently to allow a wide variety of multiscale analysis problems to be effectively and efficiently solved. The architecture of NASMAT was established specifically to enable parallelized, “plug-and-play” functionality to reduce the complexity associated with adding new features to the code in the future and to allow end users to rapidly implement and evaluate user-defined capabilities. Additionally, the tool utilizes recursive data structures and subroutines to allow for an arbitrary number of length scales when performing multiscale analyses of heterogeneous materials. These features permit the rapid integration of user-defined capabilities (e.g., a material model, micromechanics approach, or failure theory) at all stages within a NASMAT calculation while leveraging built-in techniques where needed. Additionally, these features allow NASMAT to both be called from an external program as well as call an external program. This paper specifically focuses on the multiscale integration and interoperability of NASMAT with other analysis techniques through an illustrative, multiscale analysis of a 3D woven polymer matrix composite (PMC).

NASMAT↗

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↗

Parallel Tensor Decompositions for Massive Heterogeneous, Incomplete Data (Final Report: LDRD Project 199986).

We extended the fundamental capabilities of tensor decomposition to a broader range of problems - handling new data types and larger problems. This has implications for data analysis across a range of applications in sensor monitoring, cybersecurity, treaty verification, signal processing, etc. Identifies latent structure within data, enabling anomaly detection, process monitoring, scientific discovery.

97 MATHEMATICS AND COMPUTING↗

Soil and Water Chemistry and Trace Metal Extractability and Speciation in Wetland Soils from Illinois and South Carolina and Stream Sediments from Tennessee

Dataset revised on October 15, 2021. This revision adds sulfur and iron X-ray absorption near-edge structure spectra for the wetland soils and stream sediments from the field areas. It also renames the sample locations in a way that is more intuitive to readers of the companion paper that is under review. Finally, the data filenames and organization have been updated in their labeling to parallel the data sources in the associated paper. The abstract text and methods were also revised to reflect the data that was added to the dataset.Trace metals are essential for microbially-mediated biogeochemical processes occurring in anoxic wetland soils and stream bed sediments, such as denitrification, methanogenesis, and mercury methylation. Low availability of these elements may potentially inhibit key components of anaerobic carbon and nitrogen cycling and contaminant transformation. The solid-phase speciation of trace metals likely plays an important role in controlling their bioavailability. Metal speciation is well studied in contaminated soils and sediments as well as those naturally elevated in trace metals. However, less is known regarding the chemical forms of trace metals in systems having concentrations similar to geological background levels, the very settings where metal limitations may be most prevalent. We have investigated trace metal concentrations, extractability, and solid-phase speciation in three freshwater subsurface aquatic systems: marsh wetland soils, riparian wetland soils, and the sediments of a streambed.Data are provided for marsh wetland soils at Argonne National Laboratory, riparian wetland soils in the Tims Branch watershed at Savannah River National Laboratory, and stream bed sediments from East Fork Poplar Creek near Oak Ridge National Laboratory. Soil and sediment elemental abundances, mineralogy, and extractable nutrients as well as dissolved major elements, anions, trace metals, and nutrients in the overlying surface waters are provided. In addition, the results of sequential chemical extraction for the trace metals cobalt, nickel, copper, and zinc from the soils and sediment are reported as well as X-ray absorption near-edge structure (XANES) spectra in these materials are reported. To aid interpretation of these data, XANES spectra of sulfur in the soils and sediments as well as both XANES and extended X-ray absorption fine structure (EXAFS) spectra of iron in these materials are reported. The data package also includes the XANES spectra of reference standards and a potential interferent in the measurements. All data are provided in text-based CSV format with header sections indicating the data contained in each file and the corresponding units. Note that "u" is used in place of Greek lower case mu to indicate the micro prefix on units.

54 ENVIRONMENTAL SCIENCES↗

An experimental and theoretical study of adenine adsorption on Au(111)

A model study of adenine adsorption on the Au(111) surface is reported for molecular adlayers prepared by evaporation in vacuum and deposition from saturated aqueous solution. The electronic structure and adsorption geometry of the molecular films were studied experimentally by X-ray photoelectron spectroscopy and near edge X-ray absorption fine structure spectroscopy. Adsorption models are proposed for the adlayers arising from the different preparation methods. Density functional theory calculations were used to examine both parallel and upright adenine adsorption geometries, supply additional information on the bond strength, and identify which atom is involved in bonding to Au(111). In the case of deposition in vacuum, the adenine molecule is bound via van der Waals forces to Au(111) with the molecular plane parallel to the surface, consistent with the published scanning tunneling microscopy data on this system. The most stable parallel adenine configuration was found to have an adsorption energy of ca. -1.1 eV using the optB86b-vdW functional. For adenine deposition from aqueous solution, the adlayer is disordered, with molecules in an upright geometry, and with an adsorption energy of ca. -1.0 eV, coordinated via the imino N3 nitrogen atom. The present study contributes to the substantial literature of model studies of adenine on Au(111), complementing the existing knowledge with information on electronic structure, bonding geometry and adsorption energy of this system.

Acres, Robert G.↗

Proposed crystal structure of carbadox, C 11 H 10 N 4 O 4

A model for the crystal structure of carbadox has been generated and refined using synchrotron X-ray powder diffraction data and optimized using density functional theory techniques. Carbadox crystallizes in space groupP2 1 (#4) witha= 13.8155(3),b= 21.4662(1),c= 16.3297(3) Å,β= 110.0931(7)°,V= 4548.10(3) Å 3 , andZ= 16. The crystal structure is characterized by approximately parallel stacking of the eight independent carbadox molecules parallel to thebc-plane. There are two different molecular configurations of the eight carbadox molecules; five are in the lower-energy configuration and three are in a ~10% higher-energy configuration. This arrangement likely achieves the lowest-energy crystalline packing via hydrogen bonding. Hydrogen bonds link the molecules both within and between the planes. Each of the amino groups forms a N–H⋯O hydrogen bond to an oxygen atom of the 1,4-dioxidoquinoxaline ring system of another molecule. The result is four pairs of hydrogen-bonded molecules, which form rings with graph setR2,2(14). Variation in specimen preparation can affect the preferred orientation of particles considerably. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

Materials Science↗

Intermolecular correlations of liquid and glassy CS 2 studied by synchrotron radiation x-ray diffraction

How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. Here, we have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS 2 , which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16–25.7 Å -1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS 2 were also recorded in a wide temperature range of 160–300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS 2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗

Ultra-Structural Imaging Provides 3D Organization of 46 Chromosomes of a Human Lymphocyte Prophase Nucleus

Three dimensional (3D) ultra-structural imaging is an important tool for unraveling the organizational structure of individual chromosomes at various stages of the cell cycle. Performing hitherto uninvestigated ultra-structural analysis of the human genome at prophase, we used serial block-face scanning electron microscopy (SBFSEM) to understand chromosomal architectural organization within 3D nuclear space. Acquired images allowed us to segment, reconstruct, and extract quantitative 3D structural information about the prophase nucleus and the preserved, intact individual chromosomes within it. Our data demonstrate that each chromosome can be identified with its homolog and classified into respective cytogenetic groups. Thereby, we present the first 3D karyotype built from the compact axial structure seen on the core of all prophase chromosomes. The chromosomes display parallel-aligned sister chromatids with familiar chromosome morphologies with no crossovers. Furthermore, the spatial positions of all 46 chromosomes revealed a pattern showing a gene density-based correlation and a neighborhood map of individual chromosomes based on their relative spatial positioning. A comprehensive picture of 3D chromosomal organization at the nanometer level in a single human lymphocyte cell is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mechanisms and depths of atlantic transform earthquakes

Mechanisms and depths of 40 earthquakes on major transforms along the Mid-Atlantic Ridge are studied in order to identify events that deviate from the transform-parallel strike-slip motion. Long and short period waves and Rayleigh wave spectral amplitudes are used to analyze focal mechanisms, depths, source time functions, and seismic moments of earthquakes. The relationship between centroid depths and transform thermal structures is examined. The data reveal that transform earthquake centroid depths occur above the predicted 400 C isotherms and the oceanic intraplate depths extend to the 750 C isotherm. Slip rates inferred from seismic moment releases are compared to those predicted by plate motions and good correlation is detected. The difference in the centroid depths of transform and interplate seismicity indicates transforms are either weaker or higher temperatures than expected.

Engeln, J. F.↗

F-Nets and Software Cabling: Deriving a Formal Model and Language for Portable Parallel Programming

Parallel programming is still being based upon antiquated sequence-based definitions of the terms "algorithm" and "computation", resulting in programs which are architecture dependent and difficult to design and analyze. By focusing on obstacles inherent in existing practice, a more portable model is derived here, which is then formalized into a model called Soviets which utilizes a combination of imperative and functional styles. This formalization suggests more general notions of algorithm and computation, as well as insights into the meaning of structured programming in a parallel setting. To illustrate how these principles can be applied, a very-high-level graphical architecture-independent parallel language, called Software Cabling, is described, with many of the features normally expected from today's computer languages (e.g. data abstraction, data parallelism, and object-based programming constructs).

DiNucci, David C.↗

A parallel and performance portable implementation of a full-field crystal plasticity model

We have developed a parallel implementation of an Elasto-Viscoplastic Fast Fourier Transform-based (EVPFFT) micromechanical solver to enable computationally efficient crystal plasticity modeling for polycrystalline materials. Our primary focus lies in achieving performance portability, allowing a single EVPFFT implementation to run optimally on various homogeneous architectures, including multi-core Central Processing Units (CPUs), as well as on heterogeneous computer architectures comprising multi-core CPUs and Graphics Processing Units (GPUs) from different vendors. To accomplish this goal, we have leveraged MATAR, a C++ software library that simplifies the creation and utilization of multidimensional dense or sparse matrix and array data structures. These data structures are designed to be portable across diverse architectures through the use of Kokkos, a performance-portable library. Additionally, we have employed the Message Passing Interface (MPI) to efficiently distribute the computational workload among processors. The heFFTe (Highly Efficient FFT for Exascale) library is used to facilitate the performance portability of the fast Fourier transforms (FFTs) computation. The computational performance of EVPFFT is evaluated and presented in terms of parallel scalability and simulation runtime on different high-performance computing (HPC) architectures. As a result, the utility of the developed framework to efficiently simulate the micro-mechanical fields in polycrystalline microstructures in engineering applications is discussed.

36 MATERIALS SCIENCE↗

A calibration method for overconstrained spatial translational parallel manipulators

A calibration method for overconstrained spatial translational parallel mechanisms is proposed. The calibration of a Tri-pyramid overconstrained parallel robot is used as a vehicle to demonstrate the methodology and verify its feasibility. The proposed methodology consists of transforming the overconstrained mechanism of the robot into a non-overconstrained mechanism by replacing its redundant constraints with generalized forces (couples). The kinematic model of the non-overconstrained mechanism is derived using vector loops while its mechanical model is established by the energy method. An integrated kinematic model of the original overconstrained parallel robot is subsequently established by expressing the actual position of the end effector as the sum of the position derived by the kinematic model of the non-overconstrained mechanism and of the force-induced deformations at the end effector derived from the mechanical model. To assess the feasibility and the attainable improvements in accuracy through calibration, the positions of the end effector of the Tri-pyramid robot were measured by a method that uses a Double Ball Bar (DBB) and a 3-axes linear stage with micrometers. Based on the integrated kinematic model of the robot and the measured data, its structural parameters were identified by the Newton–Raphson iteration and least squares methods to complete the calibration. The calibration result shows that accuracy improvements on the order of 90% are achievable.

42 ENGINEERING↗

Convection beneath young oceanic lithosphere - Implications for thermal structure and gravity

Small-scale convection under the oceanic lithosphere which begins in the first 5 m.y. of cooling can produce a gravity signal with the amplitude and wavelength observed for large areas of the central Pacific and southern Indian oceans using Seasat altimeter data. The trend of the observed anomalies is parallel to the direction of plate motion as might be expected if they were produced by small-scale convection. Models predict that the wavelengths of gravity anomalies increase more rapidly with age than is observed. The persistence of short relatively uniform wavelength anomalies (less than 200 km) to crustal ages of 50 Ma may indicate that they were produced when the lithosphere was very young and thin and were 'frozen in' as cooling thickened the elastic lithosphere. Small-scale convection which begins under very young lithosphere does not violate other geophysical data such as the rate of seafloor subsidence and variations of geoid height with age. After convection has begun, the subsidence due to thermal contraction within the lithosphere varies linearly with age, in the absence of mantle heat sources, although the rate of change of these quantities is affected by convection. Much of the variation of the geoid height across fracture zones can be fit by a model which includes small-scale convection.

Buck, W. R.↗

Artic and subarctic environmental analyses utilizing ERTS-1 imagery. Discipline 8: Interpretation techniques development. Subdiscipline C: Classfication and pattern recognition

The author has identified the following significant results. Uncontrolled photo mosaics of ERTS-1 imagery using MSS band 5 and 7 at a scale of 1:1,000,000 were used to make a preliminary surficial geology map in northcentral Alaska. Seven distinct geologic units were recognized, defined, and mapped directly on a photo mosaic. Results are closely correlated with published surficial geology maps. Eight MSS images were examined to test utility of ERTS data in studies of coastal processes and stream hydrology, and in the identification and interpretation of geomorphic features throughout Alaska. The feasibility of using ERTS-1 data to map structural lineaments is well illustrated on a mosaic of 8, band 5 MSS images. Along the northern edge of the Brooks Range one lineament can be followed the entire width of the mosaic, a distance of 225 miles. Two nearly parallel lineaments can be seen running along the northern and southern edges of the Schwatka Mountains. About 135 miles south of these two lineaments another series located in the Chitanana River region can be followed for 45 miles. These lineaments appear to be faults, and it is interesting to note that the Yukon River parallels these and appears to be structurally controlled.

Anderson, D. M.↗

Effects of prolonged acceleration with or without clinostat rotation on seedlings of Arabidopsis thaliana (L.) Heynh

Three 21-day tests of the effects of chronic centrifugation were carried out on populations of Arabidopsis thaliana. In addition to 1 g the resultant g-forces tested were: 2,4,6,8,16, and 20 g. Observed end points included gross morphological characters such as size of plant organs and, at the other extreme, features of sub-cellular structure and ultrastructure. Plants were grown on banks of clinostats. The acceleration vector was directed either parallel with the plants' axes or transverse to the axes. Plant responses to chronic axial acceleration and to transverse acceleration with clinostated plants were determined. From the data obtained it was possible in some cases: (1) to determine the g-functions of specific plant developmental characters; (2) to extrapolate those functions to the hypothetical value at zero g in order to predict (tentatively) the morphology of a plant grown in space, (3) to describe morphological effects of clinostat rotation, (4) to determine which of those effects was influenced by the prevailing g-force, and (5) to put to direct test the assumption that clinostat rotation nullifies or compensates for the influence of gravity.

Brown, A. H.↗

Multipoint observations of planar interplanetary magnetic field structures

IMF data made on November 1, 1984, by three spatially well-separated spacecraft in the solar wind are presented. The IMF measured by each of the spacecraft is found to consist of a multiplicity of structures within which the magnetic field varies in parallel planes. The orientations of these planes at the three spacecraft locations are similar. The planes are inclined at a large angle to the ecliptic, and they lie almost perpendicular to the nominal Parker spiral direction in the ecliptic. Intercomparisons of the measurements at the various spacecraft show that the IMF features at one spacecraft are clearly reproduced at another, with time delays required for signal propagation. From these time delays and the mutual separations of the spacecraft, it is inferred that the structures are convecting with the ambient flow. Simultaneous observations made downstream of the bow shock in the magnetosheath reveal that the magnetosheath magnetic field, too, is planar.

Farrugia, C. J.↗

Observations of mirror waves and plasma depletion layer upstream of Saturn's magnetopause

The two inbound traversals of the Saturn's magnetosheath by Voyagers 1 and 2 have been studied using plasma and magnetic field data. In a great portion of the subsolar magnetosheath, large-amplitude compressional waves are observed at low frequency (approximately 0.1 f(sub p)) in a high-beta plasma regime. The fluctuations of the magnetic field magnitude and ion density are anticorrelated, as are those of the magnetic and thermal pressures. The normals to the structures are almost orthogonal to the background field, and the Doppler ratio is on the average small. Even though the data do not allow the determination of the ion thermal anisotropy, the observations are consistent with values of T(sub perpendicular)/T(sub parallel) greater than 1, producing the onset of the mirror instability. All the above features indicate that the waves should be most probably identified with mirror modes. One of the two magnetopause crossings is of the high-shear type and the above described waves are seen until the magnetopause. The other crossing is of the low-shear type and, similarly to what has been observed at Earth, a plasma depletion occurs close to the magnetopause. In this layer, waves with smaller amplitude, presumably of the mirror mode, are present together with higher-frequency waves showing a transverse component.

Violante, L.↗