Search NASA⌕ Search

SEARCH · Search NASA

Results for “B CODES”

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 343 records · Page 19

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗

tapir: A tool for topologies, amplitudes, partial fraction decomposition and input for reductions

The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluation of Polyethylene and Blank Pulsed Sphere Experiments Using Deuteron Transport Feature in COG

From the late 1960s to about 1985, LLNL undertook a series of ‘pulsed sphere’ experiments for 31 materials involving 148 different experiments using 75 different spheres. The purpose was to measure the neutron leakage spectra for different materials from 14-MeV neutrons generated by 3 H(d,n) 4 He reactions induced by an incident D+ beam from the Insulated Core Transfer (ICT) accelerator at LLNL. This deuteron beam was focused to impinge upon tritium loaded onto a titanium substrate within a low mass target assembly placed within the center of a spherical shell constructed from the material of interest. COG is a three-dimensional, continuous energy, Monte Carlo code, developed by Lawrence Livermore National Laboratory (LLNL). Recent additions to COG include alpha and deuteron particle transport in order to directly calculate (α,n) and 3 H(d,n) 4 He reactions, and resulting neutron emission spectra. Results of the direct deuteron transport calculation for the 0.8 MFP (mean free path) polyethylene sphere with PILOT B detector was previously published. The purpose of this paper is to present expanded COG results showing the neutron peak from deuterium buildup ( 2 H(d,n) 3 He) due to tritium decay in the target. Cross section data libraries used are ENDF/B-VIII.0 for neutrons and ENDF/B-VII.1 for deuterons, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

QSpace - An open-source tensor library for Abelian and non-Abelian symmetries

This is the documentation for the tensor library QSpace (v4.0), a toolbox to exploit ‘quan tum symmetry spaces’ in tensor network states in the quantum many-body context. QSpace permits arbitrary combinations of symmetries including the abelian symmetries $\mathbb{Z}_n$ and U(1), as well as all non-abelian symmetries based on the semisimple classical Lie algebras: A n , B n , C n , and D n , or respectively, the special unitary group SU(n), the odd orthogonal group SO(2n+1), the symplectic group Sp(2n), and the even orthogonal group SO(2n). The code (C++ embedded via the MEX interface into Matlab) is available open source as of QSpace v4.0 on bitbucket under the Apache 2.0 license. QSpace is designed as a bottom-up approach for non-abelian symmetries. It starts from the defining representation and the respective Lie algebra. By explicitly comput ing and tabulating generalized Clebsch-Gordan coefficient tensors, QSpace is versatile in the type of operations that it can perform across all symmetries. At the level of an ap plication, much of the symmetry-related details are hidden within the QSpace C++ core libraries. Hence when developing tensor network algorithms with QSpace, these can be coded (nearly) as if there are no symmetries at all, despite being able to fully exploit general non-abelian symmetries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Synthesis and Crystal Structure of the Zintl Phases Na2CaCdSb2, Na2SrCdSb2 and Na2EuCdSb2

This work details the synthesis and the crystal structures of the quaternary Zintl phases Na2CaCdSb2, Na2SrCdSb2 and Na2EuCdSb2. They are isostructural and their noncentrosymmetric structure is with the space group Pmc21 (Pearson code oP12). All structural work is carried out via single-crystal X-ray diffraction methods. The structure features [CdSb2]4– layers of corner-shared CdSb4 tetrahedra, which are stacked along the b-crystallographic axis and are separated by cations. The results from the structure refinements suggest that in addition to full cation ordering, which is typical for this structure, there also exists a possibility for an accommodation of a small degree of cation disorder.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser transport and backscatter in low-density SiO 2 and Ta2O 5 foams

Experiments using a single 527 nm wavelength beam interacting with sub- and supercritical density SiO 2 and Ta 2 O 5 foams examined laser propagation and backscatter from laser–plasma instabilities such as Stimulated Brillouin Scattering (SBS). Two densities of each material were examined, and multiple diagnostics were used to characterize the propagation and backscatter. For 5 mg/cc SiO 2 (n e /n c = 0.375), the laser propagation distance was well approximated by treating the foam as a gas. However, for the 2 mg/cc SiO 2 foam (n e /n c = 0.15), the same model over-predicts the propagation distance by ~40%. Furthermore, existing analytical theories on propagation through subcritical foams were able to account for this difference. The laser heat wave propagated ~1/2 as far in Ta 2 O 5 than SiO 2 foams with similar electron density. We showed that this difference is due to the increased radiation losses in the higher Z foam. The fraction of backscattered light scales linearly with incident laser intensity for the range of intensities examined. Ta 2 O 5 foams had significantly lower levels of backscatter (1–3%) than the SiO 2 (4–8%), which is consistent with estimates of large Landau damping due to the presence of the oxygen atoms. The measured fraction of SBS backscattered laser energy for a 2 mg/cc SiO 2 foam shot was ~4 times lower than predicted by simulations assuming a gas-like foam. We found that we needed to assume increased ion heating such that T i /T e ~1.2–1.5 in the plasma to agree with the measured SBS reflectivity. Analytical models of laser-heated foams predict preferential heating of the ions as has been observed in previous experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale thermal hydraulic coupling methods for boiling water reactor simulation

During the last two years, the Virtual Environment for Reactor Applications (VERA) has been extended to simulate boiling water reactors (BWRs). The thermal hydraulic effects present in BWRs are far more complex than those in pressurized water reactors (PWRs). Therefore, the runtime is significantly increase and convergence behavior worse compared to PWRs. Most of the additional runtime is spent during the first few coupled iterations, when the power shape is still rapidly evolving, dramatically affecting the thermal hydraulics (TH). To alleviate the increased computational expense, a multiscale TH coupling approach was developed in VERA in which a highly efficient, simplified TH model is solved for several coupled iterations until the power shape is partially converged. The simplified solution is used to inform the assembly-wise flow distribution in the high-fidelity TH solver CTF, reducing the amount of work required to properly balance the pressure drop in each channel and reducing runtime for the couple calculation. Because the final TH calculations are performed with CTF, there is ultimately no impact on the accuracy of the converged solution. This paper presents the details of this multiscale TH coupling approach, along with results for a 4x4 array of GE-14 fuel bundles and whole-core coupled simulations of the Hatch core. These two cases show that the multiscale approach dramatically reduces the runtime of coupled BWR simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

S/U Comparison Study with a Focus on USLs

Under a DOE Nuclear Criticality Safety Program (NCSP) task involving Analytical Methods, three Laboratories collaborated in a comparison of results obtained from Sensitivity/Uncertainty (S/U) packages relevant to validation of transport codes. The task involves Institut de Radioprotection et de Sûreté Nucléaire (IRSN), Los Alamos National Laboratory (LANL), and Oak Ridge National Laboratory (ORNL) comparing results of MORET 5/MACSENS V3.0, MCNP6.2/Whisper-1.1, and SCALE 6.2.3/TSUNAMI/USLSTATS respectively. All Monte Carlo transport code results utilize ENDF/B-VII.1. Four cases from the International Handbook of Evaluated Criticality Safety Benchmark Experiments (ICSBEP Handbook) were selected as application models: HEU-MET-FAST-013-001, HEU-SOL-THERM 001-008, PU-MET-FAST-022-001, and PU-SOL-THERM 001-001. Ultimately, comparison is made between Upper Subcritical Limits (USLs) obtained using each code package for each application case. Since differences exist in whether packages take into account margin of subcriticality (MOS), the USL may be computed using bias and bias uncertainty, also known as the calculational margin (CM) in ANSI/ANS 8.24. Application of portions of MOS to the USL for nuclear data uncertainty of and potential code margin is referred to as USL herein. In either case, additional MOS is considered for actual application cases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sensitivity/Uncertainty Comparison Study Involving IRSN, LANL, and ORNL Tools to Support Validation [Abstract]

Under a DOE Nuclear Criticality Safety Program (NCSP) task involving Analytical Methods, three Laboratories collaborated in a comparison of results obtained from Sensitivity/Uncertainty (S/U) packages relevant to validation of transport codes. The task involves Institut de Radioprotection et de Sûreté Nucléaire (IRSN), Los Alamos National Laboratory (LANL), and Oak Ridge National Laboratory (ORNL) comparing results of MORET 5/MACSENS V3.0, MCNP6.2/Whisper-1.1, and SCALE 6.2.3/TSUNAMI/USLSTATS respectively. All Monte Carlo transport code results utilize nuclear data from ENDF/B-VII.1 evaluation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Generalized Bayesian Framework for Evaluation of Integral Benchmark Experiments

A recently published generalized Bayesian optimization framework has provided a way to retract any or all of the three common assumptions underlying the conventional Generalized Linear Least Squares (GLLS) optimization method based on the concepts introduced in reference two. These assumptions are: 1. Perfection: The model used for data evaluation and the prior probability distribution function (PDF) of generalized* data are perfect; 2. Normality: The prior and posterior PDF are normal; and 3. Linearity: The model is linear. In this work we outline how the framework in 1 could be directly adopted for improved evaluation of nuclear criticality integral benchmark experiments (IBEs) by: 1. Removing the first assumption alone by utilizing the concept of imperfections introduced in 1 to enable evaluation in the presence of discrepancies between the data and model or of missing covariance information by a GLLS method that will be seen as a generalization of the conventional GLLS method employed by the TSURFER code, and by 2. Removing the remaining two assumptions by implement- ing a Markov Chain Monte Carlo method for computation of the posterior PDF in the SAMPLER code, where TSURFER and SAMPLER are the uncertainty quantification (UQ) codes for IBEs in the SCALE code system based on the GLLS, and the stochastic method, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neutrino-dark sector equilibration and primordial element abundances

After neutrinos decouple from the photon bath, they can populate a thermal dark sector. If this neutrino cooling occurs at a temperature above ∼ 100 keV , it can impact light element abundances. We calculate the changes to the Y P and D / H abundances in this scenario with a suitably customized version of the public big bang nucleosynthesis code PRIMAT. If electron neutrinos significantly participate in the cooling, we find limits down to temperatures as low as 100 keV. And if the dark sector undergoes a “step” in N eff , there can be additional, ω b -dependent constraints. Published by the American Physical Society 2025

Giovanetti, Cara (ORCID:0000000316113379)↗

Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)

Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.

97 MATHEMATICS AND COMPUTING↗

Design and physics basis for the upcoming DIII-D SAS-VW campaign to quantify tungsten leakage and transport in a new slot divertor geometry

A set of experiments are planned to exploit the high SOL collisionality enabled by a tightly baffled slot divertor geometry to suppress tungsten leakage in DIII-D. A toroidal row of graphite tiles from the Small Angle Slot (SAS) divertor is being coated with 10–15 μ m of tungsten. New spectroscopic viewing chords with in-vacuo optics will measure the W gross erosion source from the divertor surface with high spatial and temporal resolution. In parallel, the bottom of the SAS divertor is changed from a flat to a ‘V’ shape. New SOLPS-ITER/DIVIMP simulations conducted with drifts using the planned ‘V’ shape predict a substantial reduction in W sourcing and SOL accumulation in either B × ∇B direction relative to either the old SAS divertor shape or the open, lower divertor. Dedicated studies are planned to carefully characterize the level of W sourcing, leakage, and scrape-off-layer (SOL) accumulation in DIII-D over a wide range of plasma scenarios. Various actuators will be assessed for their efficacy in further reducing high-Z impurity sources and leakage from the slot divertor geometry. Finally, this coupled code-experiment validation effort will be used to stress-test physics models and build confidence in extrapolations to advanced, high-Z divertor geometries for next-step devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The exceptional form and function of the giant bacterium Ca. Epulopiscium viviparus revolves around its sodium motive force

Epulopiscium spp. are the largest known heterotrophic bacteria; a large cigar-shaped individual is a million times the volume of Escherichia coli. To better understand the metabolic potential and relationship of Epulopiscium sp. type B with its host Naso tonganus, we generated a high-quality draft genome from a population of cells taken from a single fish. We propose the name Candidatus Epulopiscium viviparus to describe populations of this best-characterized Epulopiscium species. Metabolic reconstruction reveals more than 5% of the genome codes for carbohydrate active enzymes, which likely degrade recalcitrant host-diet algal polysaccharides into substrates that may be fermented to acetate, the most abundant short-chain fatty acid in the intestinal tract. Moreover, transcriptome analyses and the concentration of sodium ions in the host intestinal tract suggest that the use of a sodium motive force (SMF) to drive ATP synthesis and flagellar rotation is integral to symbiont metabolism and cellular biology. In natural populations, genes encoding both F-type and V-type ATPases and SMF generation via oxaloacetate decarboxylation are among the most highly expressed, suggesting that ATPases synthesize ATP and balance ion concentrations across the cell membrane. High expression of these and other integral membrane proteins may allow for the growth of its extensive intracellular membrane system. Further, complementary metabolism between microbe and host is implied with the potential provision of nitrogen and B vitamins to reinforce this nutritional symbiosis. The few features shared by all bacterial behemoths include extreme polyploidy, polyphosphate synthesis, and thus far, they have all resisted cultivation in the lab.

59 BASIC BIOLOGICAL SCIENCES↗

A weighted state redistribution algorithm for embedded boundary grids

State redistribution is an algorithm that stabilizes cut cells for embedded boundary grid methods. This work extends the earlier algorithm in several important ways. First, state redistribution is extended to three spatial dimensions. Second, we discuss several algorithmic changes and improvements motivated by the more complicated cut cell geometries that can occur in higher dimensions. In particular, we introduce a weighted version with less dissipation in an easily generalizable framework. Third, we demonstrate that state redistribution can also stabilize a solution update that includes both advective and diffusive contributions. Notably, the stabilization algorithm is shown to be effective for incompressible as well as compressible reacting flows. Finally, we discuss the implementation of the algorithm for several exascale-ready simulation codes based on AMReX, demonstrating ease of use in combination with domain decomposition, hybrid parallelism and complex physics.

97 MATHEMATICS AND COMPUTING↗

Multinucleon transfer mechanism in Cf 250 + Th 232 collisions using the quantal transport description based on the stochastic mean-field approach

Production cross sections of heavy neutron-rich isotopes are calculated by employing quantal transport description in 250 Cf + 232 Th collisions. This quantal transport description is based on the stochastic mean-field approach, and it provides a microscopic approach beyond time-dependent Hartree-Fock theory to include mean-field fluctuations. Deexcitation of primary fragments is determined by employing the statistical GEMINI++ code. Here, calculations provide predictions for production cross sections of neutron rich transfermium isotopes without any adjustable parameters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗