Search NASA⌕ Search

SEARCH · Search NASA

Results for “algorithmic 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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Monte Carlo simulation of a biological object with optical coherent tomography structural images using a voxel-based geometry of a medium

We describe a Monte Carlo algorithm for simulating an interference signal and constructing a structural image of a biological object using optical coherent tomography (OCT). The algorithm is specific in that the object geometry is described by using a three-dimensional array of voxels combined into parallelograms, and the border cross checking is carried out by using an improved Smith algorithm. The geometry of the simulated object is reconstructed based on the structure of real biological tissues obtained by OCT. The algorithm efficiency is estimated by comparing experimental and corresponding simulated images of a subcutaneous blood vessel and upper layer of human skin in vivo . (biophotonics)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Studying Transient Phenomena in Thin Films with Reinforcement Learning

Neutron reflectometry has long been a powerful tool to study the interfacial properties of energy materials. Recently, time-resolved neutron reflectometry has been used to better understand transient phenomena in electrochemical systems. Those measurements often comprise a large number of reflectivity curves acquired over a narrow q range, with each individual curve having lower information content compared to a typical steady-state measurement. In this work, we present an approach that leverages existing reinforcement learning tools to model time-resolved data to extract the time evolution of structure parameters. Further, by mapping the reflectivity curves taken at different times as individual states, we use the Soft Actor-Critic algorithm to optimize the time series of structure parameters that best represent the evolution of an electrochemical system. We show that this approach constitutes an elegant solution to the modeling of time-resolved neutron reflectometry data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Scalable implementation of polynomial filtering for density functional theory calculation in PARSEC

In this work, we present an efficient implementation of polynomial filtering methods in PARSEC, a real-space pseudopotential based Kohn–Sham density functional theory solver. The implementation described here improves upon a Chebyshev-filtered subspace iteration algorithm used in the previous version of PARSEC. We present a hybrid polynomial filtering scheme that combines Chebyshev-filtered subspace iteration and a spectrum slicing method that partitions the spectrum into several spectral slices and uses bandpass-filtered subspace iteration to compute approximate eigenpairs within each interior slice simultaneously. We describe a procedure to partition a spectrum and construct polynomial filters. We also discuss a number of practical issues such as the use of appropriate data layouts for carrying out the computation on a two-dimensional process grid and how to achieve good load balance by allocating an appropriate number of process groups to each spectral slice. Numerical examples are presented to demonstrate the effectiveness of the hybrid polynomial filtering method as well as the superior parallel scalability of spectrum slicing in comparison to that of Chebyshev-filtered subspace iteration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

The accurate particle tracer code

The Accurate Particle Tracer (APT) code is designed for systematic large-scale applications of geometric algorithms for particle dynamical simulations. Based on a large variety of advanced geometric algorithms, APT possesses long-term numerical accuracy and stability, which are critical for solving multi-scale and nonlinear problems. To provide a flexible and convenient I/O interface, the libraries of Lua and Hdf5 are used. Following a three-step procedure, users can efficiently extend the libraries of electromagnetic configurations, external non-electromagnetic forces, particle pushers, and initialization approaches by use of the extendible module. APT has been used in simulations of key physical problems, such as runaway electrons in tokamaks and energetic particles in Van Allen belt. As an important realization, the APT-SW version has been successfully distributed on the world’s fastest computer, the Sunway TaihuLight supercomputer, by supporting master–slave architecture of Sunway many-core processors. Here, based on large-scale simulations of a runaway beam under parameters of the ITER tokamak, it is revealed that the magnetic ripple field can disperse the pitch-angle distribution significantly and improve the confinement of energetic runaway beam on the same time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Understanding nanoscale structural distortions in Pb(Zr 0.2 Ti 0.8 )O 3 by utilizing X-ray nanodiffraction and clustering algorithm analysis

Hard X-ray nanodiffraction provides a unique nondestructive technique to quantify local strain and structural inhomogeneities at nanometer length scales. However, sample mosaicity and phase separation can result in a complex diffraction pattern that can make it challenging to quantify nanoscale structural distortions. In this work, a k-means clustering algorithm was utilized to identify local maxima of intensity by partitioning diffraction data in a three-dimensional feature space of detector coordinates and intensity. This technique has been applied to X-ray nanodiffraction measurements of a patterned ferroelectric PbZr 0.2 Ti 0.8 O 3 sample. The analysis reveals the presence of two phases in the sample with different lattice parameters. A highly heterogeneous distribution of lattice parameters with a variation of 0.02 Å was also observed within one ferroelectric domain. This approach provides a nanoscale survey of subtle structural distortions as well as phase separation in ferroelectric domains in a patterned sample.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Decomposition of Feynman integrals by multivariate intersection numbers

We present a detailed description of the recent idea for a direct decomposition of Feynman integrals onto a basis of master integrals by projections, as well as a direct derivation of the differential equations satisfied by the master integrals, employing multivariate intersection numbers. We discuss a recursive algorithm for the computation of multivariate intersection numbers, and provide three different approaches for a direct decomposition of Feynman integrals, which we dub the straight decomposition, the bottom-up decomposition, and the top-down decomposition. These algorithms exploit the unitarity structure of Feynman integrals by computing intersection numbers supported on cuts, in various orders, thus showing the synthesis of the intersection-theory concepts with unitarity-based methods and integrand decomposition. We perform explicit computations to exemplify all of these approaches applied to Feynman integrals, paving a way towards potential applications to generic multi-loop integrals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Derivative-free optimization of a rapid-cycling synchrotron

Here, we develop and solve a constrained optimization model to identify an integrable optics rapid-cycling synchrotron lattice design that performs well in several capacities. Our model encodes the design criteria into 78 linear and nonlinear constraints, as well as a single nonsmooth objective, where the objective and some constraints are defined from the output of Synergia, an accelerator simulator. We detail the difficulties of the 23-dimensional simulation-constrained decision space and establish that the space is nonempty. We use a derivative-free manifold sampling algorithm to account for structured nondifferentiability in the objective function. Our numerical results quantify the dependence of solutions on constraint parameters and the effect of the form of objective function.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

A Graph Theory-Based Automated Twin Recognition Technique for Electron Backscatter Diffraction Analysis

Abstract The present article introduces a new software, Microstructure Evaluation Tool for Interface Statistics (METIS), that performs high-throughput microstructure statistical analysis from electron backscatter diffraction maps. Emphasis is placed on the detection of twin domains in hexagonal close-packed metals. The numerical framework on which METIS is built leverages graph theory, group structures, and associated numerical algorithms to automatically detect twins and unravel both their intrinsic characteristics features and those pertaining to their interactions. The proposed graphical interface allows for the detection and correction of unlikely twin/parent associations rendering the approach applicable to highly deformed microstructures. Twin statistics and microstructural data are classified and saved in a relational database that can be interrogated via either GUI or SQL requests to reveal a wide spectrum of features of the microstructure. Illustration of the approach is performed in the case of zirconium.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Acceleration of the Particle Swarm Optimization for Peierls–Nabarro modeling of dislocations in conventional and high-entropy alloys

Dislocations are among the most important defects in determining the mechanical properties of both conventional alloys and high-entropy alloys. The Peierls-Nabarro model supplies an efficient pathway to their geometries and mobility. The difficulty in solving the integro-differential Peierls-Nabarro equation is how to effectively avoid the local minima in the energy landscape of a dislocation core. Among the other methods to optimize the dislocation core structures, we choose the algorithm of Particle Swarm Optimization, an algorithm that simulates the social behaviors of organisms. By employing more particles (bigger swarm) and more iterative steps (allowing them to explore for longer time), the local minima can be effectively avoided. But this would require more computational cost. The advantage of this algorithm is that it is readily parallelized in modern high computing architecture. We demonstrate the performance of our parallelized algorithm scales linearly with the number of employed cores.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Implementing a neural network interatomic model with performance portability for emerging exascale architectures

The two main thrusts of computational science are increasingly accurate predictions and faster calculations; to this end, the zeitgeist in molecular dynamics (MD) simulations is pursuing machine learned and data driven interatomic models, e.g. neural network potentials, and novel hardware architectures, e.g. GPUs. Current implementations of neural network potentials are orders of magnitude slower than traditional interatomic models and while looming exascale computing offers the ability to run large, accurate simulations with these models, achieving portable performance for MD with new and varied exascale hardware requires rethinking traditional algorithms, using novel data structures, and library solutions. We re-implement a neural network interatomic model in CabanaMD, an MD proxy application, built on libraries developed for performance portability. Our implementation shows significantly improved thread scaling in this complex kernel as compared to a current LAMMPS implementation, across both strong and weak scaling. Our single-source solution enables simulations up to 20 million atoms on a single CPU node and 4 million atoms with improved performance on a single GPU. Furthermore, we also explore parallelism and data layout choices (using flexible data structures called AoSoAs) and their effect on performance, seeing up to ~50% and ~5% improvements in performance on a GPU by choosing the right level of parallelism and data layout respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mesoscale Modeling of the Effects of Accelerated Burnup on UO2 Microstructural Evolution

Accelerating the nuclear fuel qualification process will rely on some combination of advanced modeling and simulation techniques with accelerated irradiation testing and separate effects experiments to enable the development of new fuel concepts in a shorter time frame. One of the key challenges to successfully leveraging accelerated irradiation tests will be understanding the artifacts that may be introduced with accelerated accumulation of dose and/or burnup. This work presents phase field (MARMOT) simulations of the evolution of representative 2D UO2 microstructures up to 40 MWd/kgU. Simulations were performed under both commercial light water reactor fuel conditions as well as those that would be expected for highly accelerated (~10x) burnup conditions similar to those used in the MiniFuel irradiations in Oak Ridge National Laboratory’s High Flux Isotope Reactor. The phase field model was coupled with a discrete nucleation algorithm to model re- structuring at high burnup. The effect of the different fission rates in both microstructures was investigated at two temperatures: 650?C and 800?C. The lower temperature simulations both showed an onset of restructuring at nearly 60 MWd/kgU. More extensive restructuring was obtained in the MiniFuel microstructure compared with that of the PWR fuel. At 800?C, no restructuring was obtained as a result of the thermally activated diffusion of Xe atoms and U vacancies to fission gas bubbles, which reduces the nucleation driving force. These results highlight the importance of using modeling and simulation tools to inform the environmental conditions during targeted accelerated irradiation tests to extract the most useful fuel performance data.

accelerated fuel qualification, Phase Field, Restr↗
DOE OSTI2015

A fast method for optical simulation of flood maps of light-sharing detector modules

Optical simulation of the detector module level is highly desired for Position Emission Tomography (PET) system design. Commonly used simulation toolkits such as GATE are not efficient in the optical simulation of detector modules with complicated light-sharing configurations, where a vast amount of photons need to be tracked. Here, we present a fast approach based on a simplified specular reflectance model and a structured light-tracking algorithm to speed up the photon tracking in detector modules constructed with polished finish and specular reflector materials. We also simulated conventional block detector designs with different slotted light guide patterns using the new approach and compared the outcomes with those from GATE simulations. And while the two approaches generated comparable flood maps, the new approach was more than 200–600 times faster. The new approach has also been validated by constructing a prototype detector and comparing the simulated flood map with the experimental flood map. The experimental flood map has nearly uniformly distributed spots similar to those in the simulated flood map. In conclusion, the new approach provides a fast and reliable simulation tool for assisting in the development of light-sharing-based detector modules with a polished surface finish and using specular reflector materials.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Artificial Intelligence for Conjugated Polymers

Conjugated polymers have garnered significant attention due to their diverse applications in electronics, photonics, and energy storage. However, realizing their full potential poses a formidable challenge, as their design has historically relied on iterative adjustments and continuous inspiration from researchers. Traditional methods often struggle to efficiently navigate their vast chemical landscape. In this work, the application of artificial intelligence (AI), specifically machine learning (ML), needs to be discussed in the realm of conjugated polymers. Our paper emphasizes the importance of understanding the structure–property relationships of these polymers and how ML can facilitate property prediction and inverse-design. We delve into various chemical fingerprints, structural descriptors, and ML algorithms, showcasing their utility across a spectrum of applications, including simulations, glass transition temperature determination, photovoltaics, reorganization energy for charge transport, photocatalysts, and sensors. Finally, we give some outlooks in this filed and propose unexplored areas within the field that hold the potential to benefit from ML techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Discovering type I cis-AT polyketides through computational mass spectrometry and genome mining with Seq2PKS

Type 1 polyketides are a major class of natural products used as antiviral, antibiotic, antifungal, antiparasitic, immunosuppressive, and antitumor drugs. Analysis of public microbial genomes leads to the discovery of over sixty thousand type 1 polyketide gene clusters. However, the molecular products of only about a hundred of these clusters are characterized, leaving most metabolites unknown. Characterizing polyketides relies on bioactivity-guided purification, which is expensive and time-consuming. To address this, we present Seq2PKS, a machine learning algorithm that predicts chemical structures derived from Type 1 polyketide synthases. Seq2PKS predicts numerous putative structures for each gene cluster to enhance accuracy. The correct structure is identified using a variable mass spectral database search. Benchmarks show that Seq2PKS outperforms existing methods. Applying Seq2PKS to Actinobacteria datasets, we discover biosynthetic gene clusters for monazomycin, oasomycin A, and 2-aminobenzamide-actiphenol.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Magnetism of new metastable cobalt-nitride compounds

The search for new magnetic materials with high magnetization and magnetocrystalline anisotropy is important for a wide range of applications including information and energy processing. There is only a limited number of naturally occurring magnetic compounds that are suitable. This situation stimulates an exploration of new phases that occur far from thermal-equilibrium conditions, but their stabilization is generally inhibited due to high positive formation energies. Here a nanocluster-deposition method has enabled the discovery of a set of new non-equilibrium Co–N intermetallic compounds. The experimental search was assisted by computational methods including adaptive-genetic-algorithm and electronic-structure calculations. Conventional wisdom is that the interstitial or substitutional solubility of N in Co is much lower than that in Fe and that N in Co in equilibrium alloys does not produce materials with significant magnetization and anisotropy. By contrast, our experiments identify new Co–N compounds with favorable magnetic properties including hexagonal Co 3 N nanoparticles with a high saturation magnetic polarization ($J_s$ = 1.28 T or 12.8 kG) and an appreciable uniaxial magnetocrystalline anisotropy ($K_1$ = 1.01 MJ m -3 or 10.1 Mergs per cm 3 ). This research provides a pathway for uncovering new magnetic compounds with computational efficiency beyond the existing materials database, which is significant for future technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Automated laser-based barely visible impact damage detection in honeycomb sandwich composite structures

Nondestructive evaluation (NDE) for detection and quantification of damage in composite materials is fundamental in the assessment of the overall structural integrity of modern aerospace systems. Conventional NDE systems have been extensively used to detect the location and size of damages by propagating ultrasonic waves normal to the surface. However they usually require physical contact with the structure and are time consuming and labor intensive. An automated, contactless laser ultrasonic imaging system for barely visible impact damage (BVID) detection in advanced composite structures has been developed to overcome these limitations. Lamb waves are generated by a Q-switched Nd:YAG laser, raster scanned by a set of galvano-mirrors over the damaged area. The out-of-plane vibrations are measured through a laser Doppler Vibrometer (LDV) that is stationary at a point on the corner of the grid. The ultrasonic wave field of the scanned area is reconstructed in polar coordinates and analyzed for high resolution characterization of impact damage in the composite honeycomb panel. Two methodologies are used for ultrasonic wave-field analysis: scattered wave field analysis (SWA) and standing wave energy analysis (SWEA) in the frequency domain. The SWA is employed for processing the wave field and estimate spatially dependent wavenumber values, related to discontinuities in the structural domain. The SWEA algorithm extracts standing waves trapped within damaged areas and, by studying the spectrum of the standing wave field, returns high fidelity damage imaging. While the SWA can be used to locate the impact damage in the honeycomb panel, the SWEA produces damage images in good agreement with X-ray computed tomographic (X-ray CT) scans. The results obtained prove that the laser-based nondestructive system is an effective alternative to overcome limitations of conventional NDI technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An ultra-incompressible Mn3N compound predicted by first-principles genetic algorithm

Using genetic algorithms for an unbiased structure search and first-principles total-energy calculations, a stable manganese nitride, Mn3N, is discovered. Mn3N is a nonmagnetic metal and isostructural to superhard Re3N. Mn3N exhibits a large bulk modulus and incompressibility comparable to that of the ultra-incompressible OsB. We show that the large bulk modulus can be attributed to the strong covalent bonding in this system. Phonon calculations and analysis confirm the dynamical stability of the Mn3N compound. We also show that weak electron–phonon coupling leads to a small superconducting transition temperature for Mn3N.

Zhang, Chao (ORCID:0000000259572287)↗
DOE OSTI2022

Coupled Cluster Downfolding Methods: the effect of double commutator terms on the accuracy of ground-state energies

Downfolding coupled cluster (CC) techniques have recently been introduced into quantum chemistry as a tool for the dimensionality reduction of the many-body quantum problem. As opposed to earlier formulations in physics and chemistry based on the concept of effective Hamiltonians, the appearance of the downfolded Hamiltonians is a natural consequence of the single-reference exponential parametrization of the wave function. In this paper, we discuss the impact of higher-order terms originating in double commutators. In analogy to previous studies, we consider the case when only one- and two-body interactions are included in the downfolded Hamiltonians. We demonstrate the efficiency of the many-body expansions involving single and double commutators for the unitary extension of the downfolded Hamiltonians on the example of the beryllium atom, and bond-breaking processes in the Li2 and H2O molecules. For the H2O system, we also analyze energies obtained with downfolding procedures as functions of the active space size.

quantum computing, quantum algorithms, electronic ↗