Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer program”

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 55 records · Page 3

7th World Congress on Integrated Computational Materials Engineering (ICME 2023) (Final Technical Report)

Integrated Computational Materials Engineering (ICME) has received international attention due to its potential to shorten product development time, while lowering cost and improving design and manufacturing outcomes. ICME is an approach to designing materials solutions for specific applications that use computer modeling programs to predict the behavior of materials and integrate this information into the overall materials, processing, and manufacturing design cycle. The 7th World Congress on Integrated Computational Materials Engineering (ICME 2023) was held in Orlando, Florida from May 21–25, 2023 with the goal to convene stakeholders from across all areas of modeling and simulation, experimental specialization, and design, as well as from across academia, government, and industry, to address ICME tools and techniques and their integration, as well as to examine their application in engineering. This atmosphere facilitated rich interactions between the experimentalists, modelers, and computational and design, from academia, government, and industry, to discuss ICME tools and techniques and their application in engineering.

36 MATERIALS SCIENCE↗

COLUMBUS─An Efficient and General Program Package for Ground and Excited State Computations Including Spin–Orbit Couplings and Dynamics

The COLUMBUS program system provides the tools for performing high-level multireference (MR) computations, including the multireference configuration interaction (MRCI) method and its multireference averaged quadratic coupled cluster (MR-AQCC) extension, allowing computations on a wide range of fascinating atomic and molecular systems, including the treatment of open-shells and complicated excited state phenomena. The inclusion of spin−orbit coupling (SOC) directly within the MRCI step enables the description of systems containing heavy elements, such as lanthanides and actinides, whose properties are strongly influenced by SOC. Analytic energy gradients and nonadiabatic couplings at the correlated MRCI level provide the foundation for a variety of dynamics studies, giving insight into ultrafast photochemistry. New and ongoing method developments in COLUMBUS include the computation of spin densities, improved descriptions of ionic states, enhancements to the AQCC method, and the porting of COLUMBUS to graphical processing units (GPUs). New external interfaces enable an enhanced description of electronic resonances and molecules in strong laser fields. This work highlights these new developments while providing a detailed account of the diverse applications of COLUMBUS in recent years.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficiently Computable Limits on EPR Pair Generation in Quantum Broadcast Channels

We investigate the generation of EPR pairs between three observers in a general causally structured setting, where communication occurs via a noisy quantum broadcast channel. The most general quantum codes for this setup take the form of tripartite quantum channels. Since the receivers are constrained by causal ordering, additional temporal relationships naturally emerge between the parties. These causal constraints enforce intrinsic no-signalling conditions on any tripartite operation, ensuring that it constitutes a physically realizable quantum code for a quantum broadcast channel. We analyze these constraints and, more broadly, characterize the most general quantum codes for communication over such channels. We examine the capabilities of codes that are fully no-signalling among the three parties, positive partial transpose (PPT)-preserving, or both, and derive simple semidefinite programs to compute the achievable entanglement fidelity. We then establish a hierarchy of semidefinite programming converse bounds -- both weak and strong -- for the capacity of quantum broadcast channels for EPR pair generation, in both one-shot and asymptotic regimes. Notably, in the special case of a point-to-point channel, our strong converse bound recovers and strengthens existing results. Finally, we demonstrate how the PPT-preserving codes we develop can be leveraged to construct PPT-preserving entanglement combing schemes, and vice versa.

FOS: Physical sciences↗

From zonal to nodal capacity expansion planning: Spatial aggregation impacts on a realistic test-case

Solving power system capacity expansion planning (CEP) problems at realistic spatial resolutions is computationally challenging. Thus, a common practice is to solve CEP over zonal models with low spatial resolution rather than over full-scale nodal power networks. Due to improvements in solving large-scale stochastic mixed integer programs, these computational limitations are becoming less relevant, and the assumption that zonal models are realistic and useful approximations of nodal CEP is worth revisiting. Here, this work is the first to conduct a systematic computational study on the assumption that spatial aggregation can reasonably be used for ISO-scale CEP. By considering a realistic, large-scale test network based on the state of California with over 8000 buses, we find that well-designed small spatial aggregations can yield good approximations but that coarser zonal models may result in large distortions of investment decisions, e.g., capacity under-investment of up to 41% for the lowest resolution model considered.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Future Generation High Performance Computing Center (FG-HPCC): RFI Technical Considerations

Lawrence Livermore National Security, LLC (LLNS) is interested in receiving information about technologies that could be available in the 2029-2030 timeframe that may serve to enable the vision for a Future Generation High Performance Computing (HPC) Center (FG-HPCC) described in this document. The future HPC Center vision has been conceived to meet the future mission needs of the Advanced Simulation and Computing (ASC) Program within the National Nuclear Security Administration (NNSA). LLNS envisions a center composed not of many independent clusters, but of heterogeneous elements accessible to users as a single system. The capabilities will be integrated to create a scalable, flexible, yet tightly coupled computing center capable of integrated HPC, AI, and cloud-like workloads.

97 MATHEMATICS AND COMPUTING↗

HPC Center of the Future: R&D Acquisition Intent

This document contains intended technical requirements for an anticipated future procurement for Lawrence Livermore National Laboratory (LLNL), hereafter referred to as “LLNL” or “the Laboratory”, which seeks to fund a few new and innovative proposals for developing concepts and performing studies that would lead to products available in the 2029-2030 timeframe. Our interest is to incentivize explorations and evaluations of technical capabilities in support of our vision for a High Performance Computing (HPC) center of the future at Livermore Computing (LC), the HPC center at LLNL. To realize this vision a new approach to our approach to acquisition is being considered by LC. The approach is designed to decrease risks for the vendors, while encouraging a broader set of vendors able and we hope willing to partner with LLNL and LC. The future HPC Center vision and acquisition approach was conceived to meet the future mission needs of the Advanced Simulation and Computing (ASC) Program within the National Nuclear Security Administration (NNSA).

97 MATHEMATICS AND COMPUTING↗

Tailoring Carbide Dispersed Steels: A Path to Increased Strength and Hydrogen Tolerance

The use of transition metal carbides is reported for use as a hydrogen trapping mechanism for ferritic and austenitic steel materials. The program combined computational modeling and simulations to guide experiments towards candidate metal carbide traps, both for interfacial and interior trapping. It was found that interfacial trapping is less effective than interior trapping, with the group IVB transition metal carbides being the most effect internal traps with a loss of carbon. The sub-stoichiometric rocksalt structure accommodate the hydrogen atoms in its octahedral interstices. Using percolation theory, carbon loss of approximately 25% or more was sufficient to ensure an interconnected network of vacancies for such trapping from the surface to the internal sites within the carbide. Using this as a guide, the program developed a means to provide a uniform dispersion of ZrC nanoparticles with either Fe or 304L micron-scale powders which was then consolidated by direct current sintering. Electrolytic hydrogen diffusivity studies confirmed the reduction of hydrogen diffusion in the matrix with increasing ZrC content, which was a linear response over the sample range studied (0.01 to 1.0 wt.%). The consolidated material was micro-tensile tested in either a non-hydrogen or hydrogen charge condition and compared to a control with no carbides. Additions up to 0.05 wt.% ZrC increased the yield strength with no loss in ductility in either the non-hydrogen or hydrogen tested condition. ZrC concentrations above this amount further increased the yield strength at the expense of ductility. While these samples had a lower absolute ductility value prior to failure, the relative change in ductility between the non-hydrogen and hydrogen charge states was less for the carbides than that of the control. Metal-rich ZrC nanoparticles were fabricated through a conformal coating process yielding ZrC0.66 particles that were then incorporated into a metal matrix. Notch fatigue testing in a hydrogen environment was conducted where the number of cycles to failure was found to be less in the control than that of the carbide addition. However, the spread in experimental data and the number of samples tested limits a conclusive outcome based on defects noticed in the gauge section of all the powder processed samples. The collective outcomes of this report provide further insight into the mechanisms by which carbides act as hydrogen traps; a means to process such carbides through powder metallurgy; and their associated mechanical performance in either a non-hydrogen or hydrogen-charged condition.

08 HYDROGEN↗

Reflected entropy in random tensor networks. Part III. Triway cuts

For general random tensor network states at large bond dimension, we prove that the integer Rényi reflected entropies (away from phase transitions) are determined by minimal triway cuts through the network. This generalizes the minimal cut description of bipartite entanglement for these states. A natural extrapolation away from integer Rényi parameters, suggested by the triway cut problem, implies the holographic conjecture S R = 2EW, where S R is the reflected entropy and EW is the entanglement wedge cross-section. Minimal triway cuts can be formulated as integer programs which cannot be relaxed to find a dual maximal flow/bit-thread description. This sheds light on the gap between the existence of tripartite entanglement in holographic states and the bipartite entanglement structure motivated by bit-threads. In particular, we prove that the Markov gap that measures tripartite entanglement is lower bounded by the integrality gap of the integer program that computes the triway cut.

AdS-CFT correspondence↗

Measurement of early time outer laser beam reflection inside a cylindrical hohlraum at the National Ignition Facility

In a National Ignition Facility (NIF) cylindrical hohlraum, any specular reflection (“glint”) off the wall from the outer cone (incident angle of 50° and 44° relative to the surface of the wall) laser beam will irradiate the capsule poles. If the glint power is sufficiently large during the picket (early time) of the laser pulse, it may seed high-mode perturbations on the capsule surface that can grow during the implosion. To quantify the glint power on the capsule during the picket by the outer beams, we performed dedicated experiments on NIF using a flat witness foil as a surrogate for the capsule in a half-hohlraum target. We found that the measured glint power is approximately comparable to simulations using a nominal electron conduction flux limiter of f = 0.15, and over an order of magnitude lower than that predicted using f = 0.03 in the wall. Based on our current understanding, we conclude that the glint from the outer beams plays an insignificant role in capsule drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and analysis of dudded fuel experiments at the National Ignition Facility

Recent experiments conducted at the National Ignition Facility (NIF) within the past 2 years have achieved the burning plasma state and exceeded the Lawson criterion for the first time in the laboratory. Here, we report on a set of experiments where the deuterium and tritium (DT) ice layers were replaced with dudded tritium, hydrogen, and deuterium (THD) fuel mixtures to remove the influence of alpha-heating on hot spot dynamics. The hot spot compression and yield in the absence of alpha particle self-heating were measured to assess the proximity of NIF implosions toward the ignition cliff. We find that the “burn-off” Lawson parameters χnoα inferred from the THD experiments are in good agreement with the inferences from postshot simulations of the DT-layered implosions. The THD for burning plasma shot N210307 yielded χnoα≈0.88±0.03 while the THD for ignition shot N210808 yielded χnoα≈1.04±0.04. These results also provide important context for the observed variability in the repeat attempts of ignition shot N210808 since implosions on the ignition cliff are expected to exhibit very large variations in the fusion yield from small changes in the initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MHD-induced beta limits in the Large Helical Device

Using the extended-magnetohydrodynamics code, M3D-C1, we perform a systematic numerical study of the effect of externally applied heating on the achievable plasma beta in a ten field-period heliotron. Heat sources of varying intensity are applied to a vacuum magnetic field that is representative of the standard configuration of the Large Helical Device, with R 0 = 3.66 m, where R 0 is the radial position of the magnetic axis in vacuum. As the system is driven to a state that is unstable to low-n magnetohydrodynamic (MHD) modes, nonlinear mode interactions lead to the formation of chaotic magnetic fields. With sufficiently strong heating, a collapse of the electron temperature profile is observed. This demonstrates the necessity of simulating the self-consistent evolution of plasma profiles, without imposing assumptions on the structure of the magnetic field, to accurately determine transport properties in stellarator plasmas. It also highlights the value of these advanced simulation capabilities for accelerating the development of high-performance stellarator operating scenarios.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatial description of dislocation nucleation in the shock response of single-crystal aluminum

Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.

36 MATERIALS SCIENCE↗

Evaluating nonlocal heat transport in directly driven chromium spheres using x-ray spectroscopy

We report on experiments investigating heat transport in laser-generated plasmas using directly driven chromium spheres. The spheres are fielded at the OMEGA laser facility and are driven with laser intensities of 5×10 14 Wcm −2 . Plasma conditions in the corona and scattered light are measured experimentally and compared against predictions from two-dimensional (2D) radiation-hydrodynamic simulations using different heat transport models. Spectroscopic analysis of x-ray self-emission is used as an additional diagnostic. X-ray emission is integrated over a large region of the plasma, probing regions that are not observed by localized optical Thomson scattering. In particular, x-ray emission peaks near the plasma critical density, so emission from optically thin lines provides information on plasma conditions where nonlocal transport is most likely to be significant. Three common heat transport models are considered: local transport with flux limiters f = 0.15 and f = 0.03, and the nonlocal Schurtz–Nicolai–Busquet (SNB) model. Consistent with previous work, both the high-flux (f = 0.15) and SNB models show good agreement with experimentally measured plasma conditions in the corona despite overpredicting laser absorption, whereas the low-flux (f = 0.03) model fails to match any experimental data. Conditions inferred from x-ray self-emission line ratios support this conclusion during the period of laser peak power, although synthetic spectra for all models fail to match the experiment during the transient portions of the pulse. For these reasons, the low-flux model is again rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Denoising Seismograms in the Time Domain Using a Deep Learning Model

Deep learning has emerged as a transformative tool for enhancing the extraction of reliable information from seismograms, addressing the increasing demand for precise and efficient seismic data analysis. We introduce an innovative encoder–decoder deep learning model, named WaveDenoiser, designed for noise reduction in the time domain, thereby eliminating the need for spectrogram computations that have been used for existing deep learning tools and significantly improving processing speed. Utilizing the benchmark dataset that is Stanford Earthquake Dataset, we developed three models of varying sizes: base, medium, and large. Notably, the large (referred to as WaveDenoiser) model demonstrated superior performance, achieving a median signal‐to‐noise ratio improvement of 8.8 dB on in‐distribution unseen data (in the same geographic region) and 7.7 dB on out‐distribution unseen data (in a new geographic region), outpacing both the base and medium models. Further evaluation of the WaveDenoiser model revealed a reduction in median arrival‐time errors by 0.02 s for P waves and 0.01 s for S waves when processing waveforms prior to phase picking using PhaseNet on in‐distribution unseen data. When tested on out‐distribution unseen data, the model also effectively reduced the P‐wave median arrival‐time error by 0.02 and 0.01 s in median arrival‐time error for S waves. Importantly, the application of WaveDenoiser resulted in a significant reduction of phase picking outliers by 1.1% to 3.6% for both P and S waves. In addition, we achieved over five times acceleration in processing speed compared with the seisBench implementation of DeepDenoiser. Our findings underscore the potential of WaveDenoiser as a powerful tool for improving seismic data analysis and processing efficiency.

P-waves↗

Cardinal: Seismic and Geoacoustic Array Processing

Data collected via seismic and infrasound array deployments are leveraged in the geosciences to detect and characterize a myriad of natural and anthropogenic sources. These deployments consist of numerous sensors placed in a predetermined configuration to amplify signal strength and improve the efficacy of array processing techniques used to measure signal directionality and waveform coherence. High‐fidelity feature extraction is often predicated on interstation distance as well as the frequency content and wavelength of an incident signal. Numerous array processing softwares analyze data in sequential frequency bands to obtain a more detailed characterization of a signal. However, current algorithms are limited in their ability to determine optimal array configuration for each band. We introduce an open‐source Python code, called Cardinal, to process seismic and infrasound array data in discretized time–frequency space with the option of applying an adaptive array design to determine optimal subarray configuration for each frequency band. To reduce computational time, the array processing step can be run in parallel using multithreading. Furthermore, the software has the capability to aggregate array processing results from different time–frequency pixels to produce separate sets of detections, or families, with added utility via the application of an adaptive semblance threshold, which aids in isolating signals‐of‐interest from coherent background noise. Upon appropriate configuration, Cardinal exhibits the potential to combine distinct seismic and infrasound phases into separate families.

Adaptive Array↗

SAGIPS: A scalable Framework for scidac quantom

As part of the Scientific Discovery through Advanced Computing (SciDAC) program, the Quantum Chromodynamics Nuclear Tomography (QuantOM) project aims to analyze data from Deep Inelastic Scattering (DIS) experiments conducted at Thomas Jefferson National Accelerator Facility and the upcoming Electron Ion Collider. The DIS data analysis is performed on an event level by taking into leveraging nuclear theory models and accounting for experimental conditions. In order to efficiently run multiple analyses under varying conditions, a composable workflow was designed where each section (theory, experiment, objective minimization, etc.) has its own dedicated module. This presentation gives an overview over of the current status of this workflow, highlights present and future challenges, and highlights possible extensions to other projects with similar requirements.

Lersch, Daniel [Thomas Jefferson National Accelera↗

Accelerating discoveries at DIII-D with the Integrated Research Infrastructure

DIII-D research is being accelerated by leveraging high performance computing (HPC) and data resources available through the National Energy Research Scientific Computing Center (NERSC) Superfacility initiative. As part of this initiative, a high-resolution, fully automated, whole discharge kinetic equilibrium reconstruction workflow was developed that runs at the NERSC for most DIII-D shots in under 20 min. This has eliminated a long-standing research barrier and opened the door to more sophisticated analyses, including plasma transport and stability. These capabilities would benefit from being automated and executed within the larger Department of Energy Advanced Scientific Computing Research program’s Integrated Research Infrastructure (IRI) framework. The goal of IRI is to empower researchers to meld DOE’s world-class research tools, infrastructure, and user facilities seamlessly and securely in novel ways to radically accelerate discovery and innovation. For transport, we are looking at producing flux matched profiles and also using particle tracing to predict fast ion heat deposition from neutral beam injection before a shot takes place. Our starting point for evaluating plasma stability focuses on the pedestal limits that must be navigated to achieve better confinement. This information is meant to help operators run more effective experiments, so it needs to be available rapidly inside the DIII-D control room. So far this has been achieved by ensuring the data is available with existing tools, but as more novel results are produced new visualization tools must be developed. In addition, all of the high-quality data we have generated has been collected into databases that can unlock even deeper insights. This has already been leveraged for model and code validation studies as well as for developing AI/ML surrogates. The workflows developed for this project are intended to serve as prototypes that can be replicated on other experiments and can be run to provide timely and essential information for ITER, as well as next stage fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Radiological Safety Analysis Computer (RSAC 7.2) Presentation

This presentation will be for the graduate seminar class for the University of New Mexico. We will discuss the Radiological Safety Analysis Computer (RSAC) version 7.2, a critical tool in radiological consequence analysis within the nuclear safety domain. We will begin with an introduction to radiological consequence analysis and the importance of dose modeling in nuclear safety, referencing key regulatory frameworks such as DOE-STD-3009 and 10 CFR 20. The presentation will provide an overview of RSAC, including its history, primary functions, and the inputs and outputs involved in its calculations. A detailed examination of RSAC's atmospheric dispersion model, based on the Gaussian plume methodology, will be presented, highlighting how meteorological factors influence the results. We will discuss the regulatory alignment of RSAC with DOE and NRC guidelines, emphasizing the assumptions required for compliance and the balance between conservative and realistic inputs. The presentation will also address the limitations and assumptions inherent in RSAC, comparing it with other consequence analysis codes like MACCS, HotSpot, and RASCAL. An example RSAC run will be demonstrated, showcasing a hypothetical accident scenario and interpreting the resulting dose calculations. Finally, we will introduce new ideas to enhance the current RSAC framework. This future work aims to provide multi-industry adaptability and integration with new dispersion models, addressing the evolving needs of DOE and NRC. The presentation will conclude with a discussion on RSAC's ongoing role in nuclear safety and the opportunities for future innovation.

computer program↗