Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Tool”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quantum simulation of hyperbolic space with circuit quantum electrodynamics: From graphs to geometry

We show how quantum many-body systems on hyperbolic lattices with nearest-neighbor hopping and local interactions can be mapped onto quantum field theories in continuous negatively curved space. Here, the underlying lattices have recently been realized experimentally with superconducting resonators and therefore allow for a table-top quantum simulation of quantum physics in curved background. Our mapping provides a computational tool to determine observables of the discrete system even for large lattices, where exact diagonalization fails. As an application and proof of principle we quantitatively reproduce the ground state energy, spectral gap, and correlation functions of the noninteracting lattice system by means of analytic formulas on the Poincaré disk, and show how conformal symmetry emerges for large lattices. This sets the stage for studying interactions and disorder on hyperbolic graphs in the future. Importantly, our analysis reveals that even relatively small discrete hyperbolic lattices emulate the continuous geometry of negatively curved space, and thus can be used to experimentally resolve fundamental open problems at the interface of interacting many-body systems, quantum field theory in curved space, and quantum gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Optimization of algorithmic errors in analog quantum simulations

Due to rapidly improving quantum computing hardware, Hamiltonian simulations of relativistic lattice field theories have seen a resurgence of attention. Furthermore, this computational tool requires turning the formally infinite-dimensional Hilbert space of the full theory into a finite-dimensional one. For gauge theories, a widely used basis for the Hilbert space relies on the representations induced by the underlying gauge group, with a truncation that keeps only a set of the lowest dimensional representations. This works well at large bare gauge coupling, but becomes less efficient at small coupling, which is required for the continuum limit of the lattice theory. In this work, we develop a new basis suitable for the simulation of an SU(2) lattice gauge theory in the maximal tree gauge. In particular, we show how to perform a Hamiltonian truncation so that the eigenvalues of both the magnetic and electric gauge-fixed Hamiltonian are mostly preserved, which allows for this basis to be used at all values of the coupling. Little prior knowledge is assumed, so this may also be used as an introduction to the subject of Hamiltonian formulations of lattice gauge theories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multimessenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two monoenergetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results in this work intimate the promise of this “variational annealing” methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Handing off the outcome of binary neutron star mergers for accurate and long-term postmerger simulations

In this work, we perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (~1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the HandOff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to harm3d, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the HandOff approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

New basis for Hamiltonian SU(2) simulations

Due to rapidly improving quantum computing hardware, Hamiltonian simulations of relativistic lattice field theories have seen a resurgence of attention. This computational tool requires turning the formally infinite-dimensional Hilbert space of the full theory into a finite-dimensional one. For gauge theories, a widely used basis for the Hilbert space relies on the representations induced by the underlying gauge group, with a truncation that keeps only a set of the lowest dimensional representations. This works well at large bare gauge coupling, but becomes less efficient at small coupling, which is required for the continuum limit of the lattice theory. In this work, we develop a new basis suitable for the simulation of an SU(2) lattice gauge theory in the maximal tree gauge. In particular, we show how to perform a Hamiltonian truncation so that the eigenvalues of both the magnetic and electric gauge-fixed Hamiltonian are mostly preserved, which allows for this basis to be used at all values of the coupling. Little prior knowledge is assumed, so this may also be used as an introduction to the subject of Hamiltonian formulations of lattice gauge theories. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Green's functions applied to the theory of spectroscopy

Green’s functions are a powerful analytical and computational tool for ab initio calculations of X-ray spectra. For example, Green’s functions provide an efficient means for calculations over broad energy ranges since many-body effects can be incorporated naturally in terms of the electron self-energy. Here, their role in the theory of X-ray absorption and related spectroscopies is discussed, with particular focus on many-body effects such as quasi-particle energy shifts and lifetimes, core-hole interactions and electron–phonon interactions. Additionally, the cumulant expansion for the Green’s function is reviewed and compared with the usual GW Dyson equation approach for including multielectron excitations that produce satellite features in X-ray spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Radiation-damage investigation of a DNA 16-mer

In macromolecular crystallography, a great deal of effort has been invested in understanding radiation-damage progression. While the sensitivity of protein crystals has been well characterized, crystals of DNA and of DNA–protein complexes have not thus far been studied as thoroughly. Here, a systematic investigation of radiation damage to a crystal of a DNA 16-mer diffracting to 1.8 Å resolution and held at 100 K, up to an absorbed dose of 45 MGy, is reported. The RIDL ( Radiation-Induced Density Loss ) automated computational tool was used for electron-density analysis. Both the global and specific damage to the DNA crystal as a function of dose were monitored, following careful calibration of the X-ray flux and beam profile. The DNA crystal was found to be fairly radiation insensitive to both global and specific damage, with half of the initial diffraction intensity being lost at an absorbed average diffraction-weighted dose, D 1/2 , of 19 MGy, compared with 9 MGy for chicken egg-white lysozyme crystals under the same beam conditions but at the higher resolution of 1.4 Å. The coefficient of sensitivity of the DNA crystal was 0.014 Å 2 MGy -1 , which is similar to that observed for proteins. These results imply that the significantly greater radiation hardness of DNA and RNA compared with protein observed in a DNA–protein complex and an RNA–protein complex could be due to scavenging action by the protein, thereby protecting the DNA and RNA in these studies. In terms of specific damage, the regions of DNA that were found to be sensitive were those associated with some of the bound calcium ions sequestered from the crystallization buffer. In contrast, moieties farther from these sites showed only small changes even at higher doses.

Instruments & Instrumentation↗
DOE OSTI2020

Experimental Study on the Effect of Localized Blockages on the Friction Factor of a 61-Pin Wire-Wrapped Bundle

Abstract The thermal-hydraulic behavior of the flow in rod bundles has motivated numerous experimental and computational investigations. Previous studies have identified potential for accumulation of debris within the small subchannels of typical wire-wrapped assemblies with subsequent total or partial blockage of subchannel coolant flow. A test campaign was conducted to study the effects of localized blockages on the bundle averaged friction factor of a tightly packed wire-wrapped rod bundle. Blockages were installed within the bundle, and fluid pressure drop was measured across one wire pitch for a Reynolds number range of 500–17,200. The Darcy–Weisbach friction factor of the perturbed rod bundle geometry was compared with that of the unblocked bundle, as well as with the predictions of a well-established friction factor correlation. Differing effects based on blockage size and location for various flow regimes were studied. A number of conclusions can be made about the effects of the blockages on the friction factor, such as an increasing effect of the blockage on friction factor with an increase in Reynolds number, a change in flow behavior in the turbulent transition flow regime near Reynolds number 3000, differences in effect on friction factor for different types of subchannel blockage, and a nonlinear trend in friction factor variation with flow area impeded for edge subchannels. To this end, all data and quantified uncertainty produced in this study are made available for comparison and validation of advanced computational tools.

Engineering↗
DOE OSTI2021

Numerical Simulation of Container Breach and Airborne Release of Solids Due to Mechanical Insults

Throughout U.S. Department of Energy (DOE) complexes, safety engineers employ the five-factor formula to calculate the source term (ST) that includes parameters of airborne release fraction (ARF), respirable fraction (RF) and damage ratio (DR). Limited experimental data on fragmentation of solids, such as ceramic pellets (i.e., PuO2), and container breach due to mechanical insults (i.e., drop and forklift impact), can be supplemented by modeling and simulation using high fidelity computational tools to estimate these parameters. This paper presents the use of Sandia National Laboratories' SIERRA solid mechanics (SM) finite element code to investigate the behavior of the widely utilized waste container (such as 7A Drum) subject to a range of free fall impact and puncture scenarios. The resulting behavior of the container is assessed, and the estimates are presented for bounding DRs from calculated breach areas for the various accident conditions considered. This paper also describes a novel multiscale constitutive model recently implemented in SIERRA/SM that simulates the fracture of brittle materials such as PuO 2 and determines ARF during the fracture process. Furthermore, comparisons are made between model predictions and simple bench-top experiments.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Network Uncertainty Quantification for Analysis of Multi-Component Systems

To impact physical mechanical system design decisions and realize the full promise of high-fidelity computational tools, simulation results must be integrated at the earliest stages of the design process. This is particularly challenging when dealing with uncertainty and optimizing for system-level performance metrics, as full-system models (often notoriously expensive and time-consuming to develop) are generally required to propagate uncertainties to system-level quantities of interest. Methods for propagating parameter and boundary condition uncertainty in networks of interconnected components hold promise for enabling design under uncertainty in real-world applications. These methods avoid the need for time consuming mesh generation of full-system geometries when changes are made to components or subassemblies. Additionally, they explicitly tie full-system model predictions to component/subassembly validation data which is valuable for qualification. These methods work by leveraging the fact that many engineered systems are inherently modular, being comprised of a hierarchy of components and subassemblies that are individually modified or replaced to define new system designs. By doing so, these methods enable rapid model development and the incorporation of uncertainty quantification earlier in the design process. The resulting formulation of the uncertainty propagation problem is iterative. We express the system model as a network of interconnected component models, which exchange solution information at component boundaries. We present a pair of approaches for propagating uncertainty in this type of decomposed system and provide implementations in the form of an open-source software library. We demonstrate these tools on a variety of applications and demonstrate the impact of problem-specific details on the performance and accuracy of the resulting UQ analysis. This work represents the most comprehensive investigation of these network uncertainty propagation methods to date.

42 ENGINEERING↗
DOE OSTI2023

Attosecond probing and control of charge migration in carbon-chain molecule

In quantum mechanics, when an electron is quickly ripped off from a molecule, a superposition of new eigenstates of the cation creates an electron wave packet that governs the charge flow inside, which has been called charge migration (CM). Experimentally, extracting such dynamics at its natural (attosecond) timescale is quite difficult. We report the first such experiment in a linear carbon-chain molecule, butadiyne (C 4 H 2 ), via high-harmonic spectroscopy (HHS). By employing advanced theoretical and computational tools, we showed that the wave packet and the CM of a single molecule are reconstructed from the harmonic spectra for each fixed-in-space angle of the molecule. For this one-dimensional molecule, we calculate the center of charge $\langle$ x $\rangle$ ( t ) to obtain v cm , to quantify the migration speed and how it depends on the orientation angle. The findings also uncover how the electron dynamics at the first few tens to hundreds of attoseconds depends on molecular structure. The method can be extended to other molecules where the HHS technique can be employed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Performance simulation of the soft gamma-ray concentrator

The soft gamma-ray concentrator is a telescope mission concept utilizing a suitable arrangement of bent multilayer structures of alternating low- and high-density materials. This lens is able to channel gamma-ray photons via total external reflection and concentrate the incident radiation to a point. The channeling technique offers the potential for concentrating gamma rays with focal lengths <10 m and energies >100 keV, beyond the reach of current grazing-incidence hard x-ray mirrors. For the performance estimation of such an instrument, we have developed a flexible set of computer modeling tools to compute the optical properties of multilayer structures, predict the channeling efficiency for a given multilayer configuration, and aid in the optimization of potential gamma-ray concentrator-based telescope designs. This modeling includes the multilayer optical properties calculated by the IMD software, the ray tracing using an IDL code, and the focal plane detector simulation by MEGAlib. We illustrate the potential of this approach by presenting simulated astronomical observations from a balloon-borne platform. The final result, including simulated effective area, instrument sensitivity, and polarization performance, shows that the gamma-ray concentrator will provide greatly increased sensitivity for next-generation soft gamma-ray missions with modest cost and complexity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Formation of a constructed microbial community in a nutrient-rich environment indicates bacterial interspecific competition

ABSTRACT Understanding the organizational principles of microbial communities is essential for interpreting ecosystem stability. Previous studies have investigated the formation of bacterial communities under nutrient-poor conditions or obligate relationships to observe cooperative interactions among different species. How microorganisms form stabilized communities in nutrient-rich environments, without obligate metabolic interdependency for growth, is still not fully disclosed. In this study, three bacterial strains isolated from the Populus deltoides rhizosphere were co-cultured in complex medium, and their growth behavior was tracked. These strains co-exist in mixed culture over serial transfer for multiple growth-dilution cycles. Competition is proposed as an emergent interaction relationship among the three bacteria based on their significantly decreased growth levels. The effects of different initial inoculum ratios, up to three orders of magnitude, on community structure were investigated, and the final compositions of the mixed communities with various starting composition indicate that community structure is not dependent on the initial inoculum ratio. Furthermore, the competitive relationships within the community were not altered by different initial inoculum ratios. The community structure was simulated by generalized Lotka-Volterra and dynamic flux balance analysis to provide mechanistic predictions into emergence of community structure under a nutrient-rich environment. Metaproteomic analyses provide support for the metabolite exchanges predicted by computational modeling and for highly altered physiologies when microbes are grown in co-culture. These findings broaden our understanding of bacterial community dynamics and metabolic diversity in higher-order interactions and could be significant in the management of rhizospheric bacterial communities. IMPORTANCE Bacteria naturally co-exist in multispecies consortia, and the ability to engineer such systems can be useful in biotechnology. Despite this, few studies have been performed to understand how bacteria form a stable community and interact with each other under nutrient-rich conditions. In this study, we investigated the effects of initial inoculum ratios on bacterial community structure using a complex medium and found that the initial inoculum ratio has no significant impact on resultant community structure or on interaction patterns between community members. The microbial population profiles were simulated using computational tools in order to understand intermicrobial relationships and to identify potential metabolic exchanges that occur during stabilization of the bacterial community. Studying microbial community assembly processes is essential for understanding fundamental ecological principles in microbial ecosystems and can be critical in predicting microbial community structure and function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism allows the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In current research using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For 13 C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Relativistic Magnetic Reconnection in Astrophysical Plasmas: A Powerful Mechanism of Nonthermal Emission

Magnetic reconnection—a fundamental plasma physics process, where magnetic field lines of opposite polarity annihilate—is invoked in astrophysical plasmas as a powerful mechanism of nonthermal particle acceleration, able to explain fast-evolving, bright high-energy flares. Near black holes and neutron stars, reconnection occurs in the relativistic regime, in which the mean magnetic energy per particle exceeds the rest mass energy. This review reports recent advances in our understanding of the kinetic physics of relativistic reconnection (RR): ▪ Kinetic simulations have elucidated the physics of plasma heating and nonthermal particle acceleration in RR. ▪ The physics of radiative RR, with its self-consistent interplay between photons and reconnection-accelerated particles—a peculiarity of luminous, high-energy astrophysical sources—is the new frontier of research. ▪ RR plays a key role in global models of high-energy sources, in terms of both global-scale layers and reconnection sites generated as a by-product of local magnetohydrodynamic instabilities. We summarize themes of active investigation and future directions, emphasizing the role of upcoming observational capabilities, laboratory experiments, and new computational tools.

acceleration of particles↗
DOE OSTI2025

Modeling Electrodeposition in 3D Porous Architectures for Solid-State Li-Metal Batteries

Li-metal storage in three-dimensional (3D) electrodes is considered a potential dendrite-mitigation strategy. The large surface area and high porosity of these electrodes result in reduced local Li-plating current densities. The porous topology provides a scaffold for Li-deposition and stripping, maintaining both mechanical integrity and Li accessibility. The goal of this study is to understand how characteristics, such as geometry and material properties, affect the current distribution and deposition pattern. To this end, we developed a computational method to track material growth driven by electrodeposition within a complex geometry. This method ensures that the finite-element discretization remains conforming to the moving boundary while preserving an adequate mesh quality, and thus maintains solution accuracy. Using this new computational tool, we analyze the conditions under which porous anode architectures effectively expand the surface area of the charge-transfer interface, and self-regulate current density and dendrite growth.

3D electrode architectures↗
DOE OSTI2021

Omics Mock Generator Library (OMG) v0.1.1

The Omics Mock Generator (OMG) library is used to provide the synthetic multiomics data needed to test the computational tools for bioengineering metabolic models. Since experimental multiomics data is expensive and non-trivial to produce, OMG provides a quick and easy way to produce large amounts of multiomics data that are based on plausible metabolic assumptions.

Garcia Martin, Hector↗
DOE OSTI2019

PathParser [SWR-19-54]

PathParser is a python-based computational tool for the thermodynamics and robustness analysis of native and designed metabolic pathways. The following functionalities are provided: 1 Max-min driving force (MDF) optimization. The Gibbs free energy change of the least favorable reaction will be maximized to evaluate thermodynamic feasibility of the entire pathway. 2 Protein cost estimation. The smallest enzyme investment will be assessed by calculating the minimal total enzyme cost supporting a given pathway flux. 3 Robustness analysis. Specifically, an ensemble of models is generated to simulate the system response to enzyme perturbation based on bifurcation theory and a continuation method. Probability of system failure and flux fold change against enzyme perturbation are estimated as well as a flux control index.

Xiong, Wei↗