Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “F 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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Constraints on cosmology and baryonic feedback with joint analysis of Dark Energy Survey Year 3 lensing data and ACT DR6 thermal Sunyaev-Zel'dovich effect observations

We present a joint analysis of weak gravitational lensing (shear) data obtained from the first three years of observations by the Dark Energy Survey and thermal Sunyaev-Zel'dovich (tSZ) effect measurements from a combination of Atacama Cosmology Telescope (ACT) and Planck data. A combined analysis of shear (which traces the projected mass) with the tSZ effect (which traces the projected gas pressure) can jointly probe both the distribution of matter and the thermodynamic state of the gas, accounting for the correlated effects of baryonic feedback on both observables. We detect the shear$~\times~$tSZ cross-correlation at a 21$\sigma$ significance, the highest to date, after minimizing the bias from cosmic infrared background leakage in the tSZ map. By jointly modeling the small-scale shear auto-correlation and the shear$~\times~$tSZ cross-correlation, we obtain $S_8 = 0.811^{+0.015}_{-0.012}$ and $\Omega_{\rm m} = 0.263^{+0.023}_{-0.030}$, results consistent with primary CMB analyses from Planck and P-ACT. We find evidence for reduced thermal gas pressure in dark matter halos with masses $M < 10^{14} \, M_{\odot}/h$, supporting predictions of enhanced feedback from active galactic nuclei on gas thermodynamics. A comparison of the inferred matter power suppression reveals a $2-4\sigma$ tension with hydrodynamical simulations that implement mild baryonic feedback, as our constraints prefer a stronger suppression. Finally, we investigate biases from cosmic infrared background leakage in the tSZ-shear cross-correlation measurements, employing mitigation techniques to ensure a robust inference. Our code is publicly available on GitHub.

Pandey, S. [Johns Hopkins U.; Columbia U.] (ORCID:↗
DOE OSTI2024

3-center and 4-center 2-particle Gaussian AO integrals on modern accelerated processors

We report an implementation of the McMurchie–Davidson (MD) algorithm for 3-center and 4-center 2-particle integrals over Gaussian atomic orbitals (AOs) with low and high angular momenta l and varying degrees of contraction for graphical processing units (GPUs). This work builds upon our recent implementation of a matrix form of the MD algorithm that is efficient for GPU evaluation of 4-center 2-particle integrals over Gaussian AOs of high angular momenta (l ≥ 4) [A. Asadchev and E. F. Valeev, J. Phys. Chem. A 127, 10889–10895 (2023)]. The use of unconventional data layouts and three variants of the MD algorithm allow for the evaluation of integrals with double precision and sustained performance between 25% and 70% of the theoretical hardware peak. Performance assessment includes integrals over AOs with l ≤ 6 (a higher l is supported). Preliminary implementation of the Hartree–Fock exchange operator is presented and assessed for computations with up to a quadruple-zeta basis and more than 20 000 AOs. The corresponding C++ code is part of the experimental open-source LibintX library available at https://github.com/ValeevGroup/libintx.

Chemistry↗
DOE OSTI2022

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

MSTAR2019 Code Description and User's Manual

This report describes an ideal cascade model of uranium enrichment that includes side feed and side product streams and flexible input options. It is based on the original MSTAR model developed by Ed Von Halle and implemented in a Visual Basic code. The current version allows the user to specify integer numbers of stages or the stage numbers of external flows instead of assays in those flows. The computational engine is written in FORTRAN-90 and is invoked by a user-friendly GUI written in C++. This version of MSTAR has been demonstrated to operate on Linux, Mac, and Windows platforms, and has undergone significant testing and quality analysis. A number of examples are presented to illustrate the operation of the code and guide the user. The code is extremely fast, and results are returned immediately. The input and output have been specially configured for analysis by the environmental sampling team of the International Atomic Energy Agency (IAEA). A number of additional output features are included to assist the user in visualizing computational results and downloading data to files for use in other analysis software.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

NCBI’s Virus Discovery Codeathon: Building “FIVE” —The Federated Index of Viral Experiments API Index

Viruses represent important test cases for data federation due to their genome size and the rapid increase in sequence data in publicly available databases. However, some consequences of previously decentralized (unfederated) data are lack of consensus or comparisons between feature annotations. Unifying or displaying alternative annotations should be a priority both for communities with robust entry representation and for nascent communities with burgeoning data sources. To this end, during this three-day continuation of the Virus Hunting Toolkit codeathon series (VHT-2), a new integrated and federated viral index was elaborated. This Federated Index of Viral Experiments (FIVE) integrates pre-existing and novel functional and taxonomy annotations and virus–host pairings. Variability in the context of viral genomic diversity is often overlooked in virus databases. As a proof-of-concept, FIVE was the first attempt to include viral genome variation for HIV, the most well-studied human pathogen, through viral genome diversity graphs. As per the publication of this manuscript, FIVE is the first implementation of a virus-specific federated index of such scope. FIVE is coded in BigQuery for optimal access of large quantities of data and is publicly accessible. Many projects of database or index federation fail to provide easier alternatives to access or query information. To this end, a Python API query system was developed to enhance the accessibility of FIVE.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

SOLPS-ITER validation with TCV L-mode discharges

This work presents a quantitative test of SOLPS-ITER simulations against tokamak a configuration variable (TCV) L-mode experiments. These simulations account for drifts, currents, kinetic neutrals, and carbon impurities providing the most complete edge transport simulations for TCV to date. The comparison is performed on nominally identical discharges carried out to assess the effectiveness of TCV’s divertor baffles in the framework of the European Plasma Exhaust program and employs numerous edge diagnostics providing a detailed code-experiment benchmark for TCV. The simulations show a qualitative consistency, but the quantitative differences remain, which are assessed herein. It is found that, for a given separatrix density, the simulations most notably yield a colder, and denser, divertor state with a higher divertor neutral pressure than measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Details of the neutral energy distribution and ionization source using spectrally resolved Balmer-alpha measurements on DIII-D

Spectrally resolved passive Balmer- α (D- α, H- α) measurements from the DIII-D 16 channel edge main-ion charge exchange recombination system confirm the presence of higher energy neutrals (“thermal” neutrals) in addition to the cold neutrals that recycle off the walls in the edge region of DIII-D plasmas. Charge exchange between thermal ions and edge neutrals transfers energy and momentum between the populations giving rise to thermal neutrals with energies approximating the ions in the pedestal region. Multiple charge exchange events in succession allow an electron to effectively take a random walk, transferring from ion to ion, providing a pathway of increasing energy and velocity, permitting a neutral to get deeper into the plasma before a final ionization event that contributes to the ion and electron particle fueling. Spectrally resolved measurements provide information about the density and velocity distribution of these neutrals, which has been historically valuable for validating Monte Carlo neutral models, which include the multi stage charge exchange dynamics. Here, in this study, a multi-channel set of such measurements is used to specifically isolate the details of the thermal neutrals that are responsible for fueling inside the pedestal top. Being able to separate the thermal from the cold emission overcomes several challenges associated with optical filter-based neutral density measurements. The neutral dynamics, deeper fueling by the thermal neutrals, and spectral measurement are modeled with the FIDASIM Monte Carlo collisional radiative code, which also produces synthetic spectra with a shape that is in close agreement with the measurements. By scaling the number of neutrals in the simulation to match the intensity of the thermal emission, we show it is possible to obtain local neutral densities and ionization source rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of particle transport, neutrals and radiation in magnetically-confined plasmas with Aurora

In this work, we present Aurora, an open-source package for particle transport, neutrals and radiation modeling in magnetic confinement fusion plasmas. Aurora's modern multi-language interface enables simulations of 1.5D impurity transport within high-performance computing frameworks, particularly for the inference of particle transport coefficients. A user-friendly Python library allows simple interaction with atomic rates from the Atomic Data and Atomic Structure database as well as other sources. This enables a range of radiation predictions, both for power balance and spectroscopic analysis. We discuss here the superstaging approximation for complex ions, as a way to group charge states and reduce computational cost, demonstrating its wide applicability within the Aurora forward model and beyond. Aurora also facilitates neutral particle analysis, both from experimental spectroscopic data and other simulation codes. Leveraging Aurora's capabilities to interface SOLPS-ITER results, we demonstrate that charge exchange is unlikely to affect the total radiated power from the ITER core during high performance operation. Finally, we describe the ImpRad module in the one modeling framework for integrated task framework, developed to enable experimental analysis and transport inferences on multiple devices using Aurora.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Understanding detachment of the W7-X island divertor

The fundamental behavior of the W7-X island divertor under detached conditions, which has been theoretically predicted with the EMC3-Eirene code, is re-examined here under the experimental conditions achieved so far and compared with the first experimental results. Both simulations and experiments cover a range of divertor configurations and plasma parameters, and show the following common trends: (1) with rising impurity radiation, the target heat load decreases 'uniformly' over the entire target surface in the sense that both the peak and average heat loads can drop by an order of magnitude. Impurity radiation (mainly from intrinsic carbon) occurs primarily at the plasma edge and the resulting negative impact on the stored energy is less than 10%. (2) When the total radiation exceeds a critical level, the target particle flux (the recycling flux Γrecy) begins to fall and can drop by a factor of 3–5 at high radiation levels without an obvious indication of significant volume recombination. (3) While Γ recy decreases, the divertor neutral pressure continues to build up and reaches a maximum, at which point Γ recy has declined significantly. (4) During detachment, the electron temperature at the last closed flux surface falls in a way that is not quantitatively understandable from parallel classical heat conduction processes. This paper presents a physical explanation of the numerical/experimental results described above. Furthermore, using the EMC3-Eirene code as a diagnostic tool, we are able, apparently for the first time, to provide a full quantitative analysis of each transport channel in the island divertor, aiming to clarify how the island divertor plasma self-regulates to maintain particle, energy, and momentum balance under detached conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Atmospheric neutrinos in next-generation xenon and argon dark matter experiments

In this work, we study the sensitivity of future xenon- and argon-based dark matter and neutrino detection experiments to low-energy atmospheric neutrinos. Not accounting for experimental backgrounds, the primary obstacle for identifying nuclear recoils induced by atmospheric neutrinos in xenon is the tail of the electron recoil distribution due to pp solar neutrinos. We use the NEST code to model the solar and atmospheric neutrino signals in a xenon detector and find that an exposure of 700 tonne-years will produce a 5σ detection of atmospheric neutrinos. We explore the effect of different detector properties and find that a sufficiently long electron lifetime is essential to the success of such a measurement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Performance of an Astrophysical Radiation Hydrodynamics Code under Scalable Vector Extension Optimization

We present results of a performance study of an astrophysical radiation hydrodynamics code, V2D, on the Arm-based A64FX processor developed by Fujitsu. The code solves sparse linear systems, a task for which the A64FX architecture should be well suited. Here, we performed the performance analysis study on Ookami, an Apollo 80 platform utilizing the A64FX processor. We explored several compilers and performance anal-ysis packages and found the code did not perform as expected under scalable vector extension optimization, suggesting that a “deeper dive” into analyzing the code is worthwhile. However, a simple driver program that exercised basic sparse linear algebra routines used by V2D did show significant speedup with the use of the scalable vector extension optimization. We present the initial results from the study which used V2D on a relatively simple test problem that emphasized the repeated solution of sparse linear systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Implementation of a real-time MSE system

Motional Stark effect polarimetry is a key diagnostic for plasma fusion research since its usage on PBX-M. The MSE diagnostic measures the radial magnetic pitch angle profile in a plasma from a neutral beam by observation of Stark split D-alpha emission from atoms excited by collision with ions and electrons in the plasma. The pitch angle measurement is used with equilibrium reconstruction codes to determine the q-profile for studies of plasma stability, confinement, and transport. Historically, the algorithm was used in a post-processing fashion. The goal of our work was to apply this method in real time and pass the results to the plasma control system computer for real-time equilibrium reconstruction and control.

Instruments & Instrumentation↗
DOE OSTI2024

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Bayesian inference applied to x-ray spectroscopy data analysis enables uncertainty quantification necessary to rigorously test theoretical models. However, when comparing to data, detailed atomic physics and radiation transfer calculations of x-ray emission from non-uniform plasma conditions are typically too slow to be performed in line with statistical sampling methods, such as Markov Chain Monte Carlo sampling. Furthermore, differences in transition energies and x-ray opacities often make direct comparisons between simulated and measured spectra unreliable. Here, we present a spectral decomposition method that allows for corrections to line positions and bound–bound opacities to best fit experimental data, with the goal of providing quantitative feedback to improve the underlying theoretical models and guide future experiments. In this work, we use a neural network (NN) surrogate model to replace spectral calculations of isobaric hot-spots created in Kr-doped implosions at the National Ignition Facility. The NN was trained on calculations of x-ray spectra using an isobaric hot-spot model post-processed with Cretin, a multi-species atomic kinetics and radiation code. The speedup provided by the NN model to generate x-ray emission spectra enables statistical analysis of parameterized models with sufficient detail to accurately represent the physical system and extract the plasma parameters of interest.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Code-to-code comparison between FLASH and HYDRA in gas-puff Z-pinch modeling

The numerical modeling of gas-puff Z pinches involves the nonlinear coupling of a broad range of complex, multi-physics phenomena that makes such simulations challenging. The challenge is further compounded by nonlinear instabilities that can impact the dynamics of imploding gas-puff Z pinches, such as the magneto Rayleigh–Taylor instability (MRTI). If the growth rate and amplitude of the latter is comparable to the relevant timescales and properties of the imploding plasma, the MRTI can dramatically alter implosion dynamics, dictate pinch stability, and govern the plasma properties achievable in pulsed-power-driven laboratory experiments. National Laboratories and academic teams have developed numerical tools that can accurately model Z-pinch configurations and provide reliable design capabilities that can guide experimental choices and assist in interpreting experimental results. Most such tools, however, are not broadly available. Here, we present newly developed Z-pinch simulation capabilities of the publicly available FLASH code, applied in the study of MRTI growth and dynamical effects in gas-puff implosions. To verify the new implementations, we perform a comparison of FLASH gas-puff implosion simulations with previously published calculations with the HYDRA code from Lawrence Livermore National Laboratory, which have been validated with experimental data from the CESZAR pulsed-power driver at the University of California, San Diego. The experiments involved double- and triple-nozzle configurations, in an experimental attempt to stabilize the pinch to the MRTI. The code-to-code comparison shows similar results between the FLASH and HYDRA simulations, supporting the use of FLASH in the modeling of future gas-puff Z-pinch experiments at CESZAR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit

Abstract Summary Bayesian inference in biological modeling commonly relies on Markov chain Monte Carlo (MCMC) sampling of a multidimensional and non-Gaussian posterior distribution that is not analytically tractable. Here, we present the implementation of a practical MCMC method in the open-source software package PyBioNetFit (PyBNF), which is designed to support parameterization of mathematical models for biological systems. The new MCMC method, am, incorporates an adaptive move proposal distribution. For warm starts, sampling can be initiated at a specified location in parameter space and with a multivariate Gaussian proposal distribution defined initially by a specified covariance matrix. Multiple chains can be generated in parallel using a computer cluster. We demonstrate that am can be used to successfully solve real-world Bayesian inference problems, including forecasting of new Coronavirus Disease 2019 case detection with Bayesian quantification of forecast uncertainty. Availability and implementation PyBNF version 1.1.9, the first stable release with am, is available at PyPI and can be installed using the pip package-management system on platforms that have a working installation of Python 3. PyBNF relies on libRoadRunner and BioNetGen for simulations (e.g. numerical integration of ordinary differential equations defined in SBML or BNGL files) and Dask.Distributed for task scheduling on Linux computer clusters. The Python source code can be freely downloaded/cloned from GitHub and used and modified under terms of the BSD-3 license (https://github.com/lanl/pybnf). Online documentation covering installation/usage is available (https://pybnf.readthedocs.io/en/latest/). A tutorial video is available on YouTube (https://www.youtube.com/watch?v=2aRqpqFOiS4&t=63s). Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Observation of synergy between lower hybrid waves at two frequencies in EAST

Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (n_e,lin ≈ 2.0 × 10^19m^−3), simultaneous injection of a 65/35 mix of 2.45 GHz/4.6 GHz power achieved an LHCD efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker-Planck solver CQL3D. For each discharge, profiles of current and hard x-ray from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.

Lower hybrid synergy↗
DOE OSTI2021

PETSc Users Manual (Revision 3.15)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Observation of synergy between lower hybrid waves at two frequencies in EAST

Synergistic effects between two frequencies of lower hybrid (LH) waves—operating at 2.45 and 4.6 GHz—were observed in experiment on EAST for the first time. At low density (ne,lin≈2.0×1019 m−3), simultaneous injection of a 65/35 mix of 2.45/4.6 GHz power achieved an lower hybrid current drive efficiency that was 25% higher than what should be expected from the linear combination of the two sources. The experiment was interpreted with time-dependent simulations, using the equilibrium and transport solver TRANSP, coupled with the ray-tracing code GENRAY and the Fokker–Planck solver CQL3D. For each discharge, profiles of current and hard x-rays from simulation and measurement agree within uncertainties. An examination of the electron distribution function indicates that the LH synergy is supported by the increased width of the LH resonance plateau in the simultaneous injection case compared to independent injection.

Choi, W. (ORCID:0000000268912300)↗