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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A Higher Order, Stable Partitioned Scheme for Fluid-Structure Interaction Problems

Although still a very active area of research with many open questions, there exist highly-accurate and efficient algorithms for numerically estimating solutions to the equations of fluid motion. Similarly, great strides have been made in numerically modeling the motion of solids. However, there exist many practical applications where one must solve both fluid and structure in tandem. There is a gap in the literature for high order unconditionally stable partitioned algorithms for problems of this type. We investigate a new partitioned scheme which shows promise in filling this gap. We use a finite-element-based model with an added mass approach on problems coupling the Euler beam equation with the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Ternary Bismuthide SrPtBi 2 : Computation and Experiment in Synergism to Explore Solid-State Materials

A combination of theoretical calculation and the experimental synthesis to explore the new ternary compound is demonstrated in the Sr–Pt–Bi system. Because Pt–Bi is considered as a new critical charge-transfer pair for superconductivity, it inspired us to investigate the Sr–Pt–Bi system. With a thorough calculation of all the known stable/metastable compounds in the Sr–Pt–Bi system and crystal structure predictions, the thermodynamic stability of hypothetical stoichiometry, SrPtBi2, is determined. Following the high-temperature synthesis and crystallographic analysis, the first ternary bismuthide in Sr–Pt–Bi, SrPtBi2 was prepared, and the stoichiometry was confirmed experimentally. SrPtBi 2 crystallizes in the space group Pnma (S.G. 62, Pearson Symbol oP48), which matches well with theoretical prediction using an adaptive genetic algorithm. Using first-principles calculations, we demonstrate that the orthorhombic structure has lower formation energies than other 112 structure types, such as tetragonal BaMnBi 2 (CuSmP 2 ) and LaAuBi 2 (CuHfSi 2 ) structure types. The bonding analysis indicates that the Pt–Bi interactions play a critical role in structural stability. The physical property measurements show the metallic properties at the low temperature, which agrees with the electronic structure assessment.

36 MATERIALS SCIENCE↗

DOE OSTI2015

A revised partiality model and post-refinement algorithm for X-ray free-electron laser data

An updated partiality model and post-refinement algorithm for XFEL snapshot diffraction data is presented and confirmed by observing anomalous density for S atoms at an X-ray wavelength of 1.3 Å. Research towards using X-ray free-electron laser (XFEL) data to solve structures using experimental phasing methods such as sulfur single-wavelength anomalous dispersion (SAD) has been hampered by shortcomings in the diffraction models for X-ray diffraction from FELs. Owing to errors in the orientation matrix and overly simple partiality models, researchers have required large numbers of images to converge to reliable estimates for the structure-factor amplitudes, which may not be feasible for all biological systems. Here, data for cytoplasmic polyhedrosis virus type 17 (CPV17) collected at 1.3 Å wavelength at the Linac Coherent Light Source (LCLS) are revisited. A previously published definition of a partiality model for reflections illuminated by self-amplified spontaneous emission (SASE) pulses is built upon, which defines a fraction between 0 and 1 based on the intersection of a reflection with a spread of Ewald spheres modelled by a super-Gaussian wavelength distribution in the X-ray beam. A method of post-refinement to refine the parameters of this model is suggested. This has generated a merged data set with an overall discrepancy (by calculating the R{sub split} value) of 3.15% to 1.46 Å resolution from a 7225-image data set. The atomic numbers of C, N and O atoms in the structure are distinguishable in the electron-density map. There are 13 S atoms within the 237 residues of CPV17, excluding the initial disordered methionine. These only possess 0.42 anomalous scattering electrons each at 1.3 Å wavelength, but the 12 that have single predominant positions are easily detectable in the anomalous difference Fourier map. It is hoped that these improvements will lead towards XFEL experimental phase determination and structure determination by sulfur SAD and will generally increase the utility of the method for difficult cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Feasibility Study on Implementing a Staggered-Grid Finite Volume Method for System Analysis Code Development Under the MOOSE Framework

Here, this work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key to the test bed is the implementation of high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. The test bed utilized a more flexible code structure to enable the finite volume method implementation and direct interacting with the solver package, instead of using the natively supported finite element method by the framework. Using a suite of selected test problems with different problem sizes and levels of complexity, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. For a complex reactor model, transient simulation was performed using the newly developed finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development.

MOOSE↗

DOE OSTI2013

First-principles simulation of solvation structure and deprotonation reactions of ore forming metal ions in very nonideal solutions

Advances in algorithms and computational performance achieved in this grant period have allowed the atomic level dynamical simulation of complex nanoscale materials using interparticle forces calculated directly from an accurate density functional solution to the electronic Schrödinger equation (ab-initio molecular dynamics, AIMD). Focus of this program was on the prediction and analysis of the properties of environmentally important ions in aqueous solutions. AIMD methods have provided chemical interpretations of these very complex systems with an unprecedented level of accuracy and detail. The structure of the solvation region neighboring a highly charged metal ion (e.g., 3+) in an aqueous solution is very different from that of bulk water. The many-body behaviors (polarization, charge transfer, etc.) of the ion-water and water-water interactions in this region are difficult to capture with conventional empirical potentials. However, a large numbers of waters (up to 128 waters) are required to fully describe chemical events in the extended hydrations shells and long simulation times are needed to reliably sample the system. Taken together this makes simulation at the 1 st principles level a very large computational problem. Our AIMD simulation results using these methods agree with the measured octahedral structure of the 1 st solvation shell of Al 3+ at the 1 st shell boundary and a calculated radius of 1.937 (exp. 1.9). Our calculated average 2 nd shell radius agrees remarkably well with the measured radius, 4.093 Å calculated vs. the measured value of 4.0-4.15 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Unsupervised discovery of extreme weather events using universal representations of emergent organization

Spontaneous self-organization is ubiquitous in systems far from thermodynamic equilibrium. While organized structures that emerge dominate transport properties, universal representations that identify and describe these key objects remain elusive. Here, we introduce a theoretically grounded framework for describing emergent organization that, via data-driven algorithms, is constructive in practice. Its building blocks are spacetime lightcones that embody how information propagates across a system through local interactions. We show that predictive equivalence classes of lightcones—local causal states—capture organized behaviors in complex spatiotemporal systems. Employing an unsupervised physics-informed machine learning algorithm and a high-performance computing implementation, we demonstrate automatically discovering organized structures in two real-world domain science problems. We show that local causal states identify vortices and track their power-law decay behavior in two-dimensional fluid turbulence. We then show how to detect and track familiar extreme weather events—hurricanes and atmospheric rivers—and discover other novel structures associated with precipitation extremes in high-resolution climate data at the grid-cell level.

Rupe, Adam [Pacific Northwest National Laboratory ↗

DOE OSTI2025

Dissipative ground state preparation in ab initio electronic structure theory

Dissipative engineering is a powerful tool for quantum state preparation, and has drawn significant attention in quantum algorithms and quantum many-body physics in recent years. In this work, we introduce a novel approach using the Lindblad dynamics to efficiently prepare the ground state for general ab initio electronic structure problems on quantum computers, without variational parameters. These problems often involve Hamiltonians that lack geometric locality or sparsity structures, which we address by proposing two generic types of jump operators for the Lindblad dynamics. Type-I jump operators break the particle number symmetry and should be simulated in the Fock space. Type-II jump operators preserves the particle number symmetry and can be simulated more efficiently in the full configuration interaction space. For both types of jump operators, we prove that in a simplified Hartree-Fock framework, the spectral gap of our Lindbladian is lower bounded by a universal constant. For physical observables such as energy and reduced density matrices, the convergence rate of our Lindblad dynamics with Type-I jump operators remains universal, while the convergence rate with Type-II jump operators only depends on coarse grained information such as the number of orbitals and the number of electrons. To validate our approach, we employ a Monte Carlo trajectory-based algorithm for simulating the Lindblad dynamics for full ab initio Hamiltonians, demonstrating its effectiveness on molecular systems amenable to exact wavefunction treatment.

Quantum chemistry↗

DOE OSTI2019

Structure retrieval from four-dimensional scanning transmission electron microscopy: Statistical analysis of potential pitfalls in high-dimensional data

Four-dimensional (4D) scanning transmission electron microscopy is one of the most rapidly growing modes of electron microscopy imaging. The advent of fast pixelated cameras and the associated data infrastructure have greatly accelerated this process. Yet conversion of the 4D datasets into physically meaningful structure images in real space remains an open issue. In this work, we demonstrate that it is possible to systematically create filters that will affect the apparent resolution or even qualitative features of the real-space structure image, reconstructing artificially generated patterns. As initial efforts, we explore statistical model selection algorithms, aiming for robustness and reliability of estimated filters. This statistical model selection analysis demonstrates the need for regularization and cross validation of inversion methods to robustly recover structure from high-dimensional diffraction datasets.

36 MATERIALS SCIENCE↗

DOE OSTI2022

SDynPy: A Structural Dynamics Python Library

SAND2023-11957O SDynPy software can be used in digital signal processing, modal analysis, and geometry algorithms that are available in open literature, specifically the Synthesize Modes and Correlate and polynomial-based multiple reference modal fitters. This software can be used to perform structural dynamic testing and analysis. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Rohe, Daniel↗

DOE OSTI2025

A physical basis for cosmological correlators from cuts

Significant progress has been made in our understanding of the analytic structure of FRW wavefunction coefficients, facilitated by the development of efficient algorithms to derive the differential equations they satisfy. Moreover, recent findings indicate that the twisted cohomology of the associated hyperplane arrangement defining FRW integrals overestimates the number of integrals required to define differential equations for the wave-function coefficient. We demonstrate that the associated dual cohomology is automatically organized in a way that is ideal for understanding and exploiting the cut/residue structure of FRW integrals. Utilizing this understanding, we develop a systematic approach to organize compatible sequential residues, which dictates the physical subspace of FRW integrals for any n -site, ℓ-loop graph. In particular, the physical subspace of tree-level FRW wavefunction coefficients is populated by differential forms associated to cuts/residues that factorize the integrand of the wavefunction coefficient into only flat space amplitudes. After demonstrating the validity of our construction using intersection theory, we develop simple graphical rules for cut tubings that enumerate the space of physical cuts and, consequently, differential forms without any calculation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Model-Based Reconstruction for Collimated Beam Ultrasound Systems

Collimated beam ultrasound systems are a novel technology for imaging inside multi-layered structures such as geothermal wells. Such systems include a transmitter and multiple receivers to capture reflected signals. Common algorithms for ultrasound reconstruction use delay-and-sum (DAS) approaches; these have low computational complexity but produce inaccurate images in the presence of complex structures and specialized geometries such as collimated beams.In this paper, we propose a multi-layer, ultrasonic, model-based iterative reconstruction algorithm designed for collimated beam systems. We introduce a physics-based forward model to accurately ac-count for the propagation of a collimated ultrasonic beam in multi-layer media and describe an efficient implementation using binary search. We model direct arrival signals, detector noise, and a spatially varying image prior, then cast the reconstruction as a maximum a posteriori estimation problem. Using simulated and experimental data we obtain significantly fewer artifacts relative to DAS while running in near real time using commodity compute resources.

Alanazi, Abdulrahman M.↗

DOE OSTI2019

Computationally Driven Discovery of a Family of Layered LiNiB Polymorphs

Abstract Two novel lithium nickel boride polymorphs, RT ‐LiNiB and HT ‐LiNiB, with layered crystal structures are reported. This family of compounds was theoretically predicted by using the adaptive genetic algorithm (AGA) and subsequently synthesized by a hydride route with LiH as the lithium source. Unique among the known ternary transition‐metal borides, the LiNiB structures feature Li layers alternating with nearly planar [NiB] layers composed of Ni hexagonal rings with a B–B pair at the center. A comprehensive study using a combination of single crystal/synchrotron powder X‐ray diffraction, solid‐state 7 Li and 11 B NMR spectroscopy, scanning transmission electron microscopy, quantum‐chemical calculations, and magnetism has shed light on the intrinsic features of these polymorphic compounds. The unique layered structures of LiNiB compounds make them ultimate precursors for exfoliation studies, thus paving a way toward two‐dimensional transition‐metal borides, MBenes.

Gvozdetskyi, Volodymyr↗

DOE OSTI2019

Computationally Driven Discovery of a Family of Layered LiNiB Polymorphs

Gvozdetskyi, Volodymyr↗

DOE OSTI2021

Light-Front Field Theory on Current Quantum Computers

We present a quantum algorithm for simulation of quantum field theory in the light-front formulation and demonstrate how existing quantum devices can be used to study the structure of bound states in relativistic nuclear physics. Specifically, we apply the Variational Quantum Eigensolver algorithm to find the ground state of the light-front Hamiltonian obtained within the Basis Light-Front Quantization (BLFQ) framework. The BLFQ formulation of quantum field theory allows one to readily import techniques developed for digital quantum simulation of quantum chemistry. This provides a method that can be scaled up to simulation of full, relativistic quantum field theories in the quantum advantage regime. As an illustration, we calculate the mass, mass radius, decay constant, electromagnetic form factor, and charge radius of the pion on the IBM Vigo chip. This is the first time that the light-front approach to quantum field theory has been used to enable simulation of a real physical system on a quantum computer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Experimental demonstration of coupled multi-peak Bragg coherent diffraction imaging with genetic algorithms

Bragg coherent diffraction imaging has the potential to provide significant insight into the structure-properties relationship for crystalline materials by imaging, with nanoscale resolution, three-dimensional strain fields within individual grains and nanoparticles. The capability of present-day synchrotrons to locate and measure a multiplicity of Bragg reflections from a single grain makes it possible to recover the full strain tensor with nanometer resolution. Recent methods for coupling reconstructions from several peaks to determine the strain tensor have been developed and applied to synthetic data, but have not been applied to experimental data. Here, using a coupled genetic reconstruction algorithm, we reconstruct an experimental data set and demonstrate improvements in the ability to resolve vector-valued displacement fields internal to the particle as compared to what is achieved with a noncoupled approach. The coupled approach developed in this work was also validated on simulated data sets. In both simulated and experimental data, reconstructions from our coupled Bragg peak algorithm show improvements over the noncoupled independent reconstruction method of 5% in terms of accuracy and 53% in terms of consistency.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2016

A Retrieval of Tropical Latent Heating Using the 3D Structure of Precipitation Features

Radar-based latent heating retrievals typically apply a lookup table (LUT) derived from model output to surface rain amounts and rain type to determine the vertical structure of heating. In this study, a method has been developed that uses the size characteristics of precipitating systems (i.e., area and mean echo-top height) instead of rain amount to estimate latent heating profiles from radar observations. This technique [named the convective–stratiform area (CSA) algorithm] leverages the relationship between the organization of convective systems and the structure of latent heating profiles and avoids pitfalls associated with retrieving accurate rainfall information from radars and models. The CSA LUTs are based on a high-resolution regional model simulation over the equatorial Indian Ocean. The CSA LUTs show that convective latent heating increases in magnitude and height as area and echo-top heights grow, with a congestus signature of midlevel cooling for less vertically extensive convective systems. Stratiform latent heating varies weakly in vertical structure, but its magnitude is strongly linked to area and mean echo-top heights. The CSA LUT was applied to radar observations collected during the DYNAMO/Cooperative Indian Ocean Experiment on Intraseasonal Variability in the Year 2011 (CINDY2011)/ARM MJO Investigation Experiment (AMIE) field campaign, and the CSA heating retrieval was generally consistent with other measures of heating profiles. Finally, the impact of resolution and spatial mismatch between the model and radar grids is addressed, and unrealistic latent heating profiles in the stratiform LUT, namely, a low-level heating peak, an elevated melting layer, and net column cooling, were identified. These issues highlight the need for accurate convective–stratiform separations and improvement in PBL and microphysical parameterizations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

3D Experimental Measurements of Evolution of Force Chains in Natural Silica Sand

The mechanisms of force transmission in granular materials is a classic physics problem that has been addressed since the 19th century, when Heinrich Rudolf Hertz investigated the interaction between two similar objects that were in contact under compression. However, the study of force transmission mechanisms in assemblies of more particles has proven to be a formidable problem due to the complex nature of granular materials. In recent years, synchrotron microcomputed tomography (SMT) and three-dimensional X-ray diffraction microscopy (3DXRD) have been employed to study the mechanics of granular materials experimentally. Combining SMT and 3DXRD offers unique three-dimensional (3D) experimental measurements of the internal structure, kinematics (such as rotation and translation), and lattice strains of individual sand particles. In this paper, in situ SMT and 3DXRD scans were acquired at multiple load steps for a specimen composed of 2,705 natural Ottawa sand particles that were subjected to one-dimensional (1D) confined compression. An algorithm was developed to combine SMT images and 3DXRD lattice strain measurements and used to characterize the constitutive behavior of sand particles. The results were used to identify the crystal structure and the evolution of the stresses and lattice strains of individual sand particles. Another algorithm was developed to characterize the force structures within the specimen. Force structures were identified, and their properties (such as length) and evolution through the experiment were examined. The contact number of particles is a particle-scale property that affects the mechanics of granular materials. The effect of the contact number of the sand particles on the onset and evolution of the force structures was also investigated and discussed.

3D x-ray diffraction microscopy↗

DOE OSTI2023

ReaLigands: A Ligand Library Cultivated from Experiment and Intended for Molecular Computational Catalyst Design

Computational catalyst design requires identification of a metal and ligand that together result in the desired reaction reactivity and/or selectivity. A major impediment to translating computational designs to experiments is evaluating ligands that are likely to be synthesized. Here we provide a solution to this impediment with our ReaLigands library that contains >30,000 monodentate, bidentate (didentate), tridentate, and larger ligands cultivated by dismantling experimentally reported crystal structures. Individual ligands from mononuclear crystal structures were identified using a modified depth-first search algorithm and charge was assigned using a machine learning model based on quantum-chemical calculated features. In the library ligands are sorted based on direct ligand-to-metal atomic connections and on denticity. Representative principal component analysis (PCA) and uniform manifold approximation and projection (UMAP) analyses were used to analyze several tridentate ligand categories, which revealed both the diversity of ligands and connections between ligand categories. Furthermore, we also demonstrated the utility of this library by implementing it with our building and optimization tools, which resulted in the very rapid generation of barriers for 750 bidentate ligands for Rh-hydride ethylene migratory insertion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗