Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “M 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The accuracy of multi-group models for nonlocal electron transport in magnetized plasmas

In the extreme conditions of inertial confinement fusion experiments, heat flow plays a vital role, but local diffusive models frequently break down and overestimate the heat flow. The situation becomes more complicated again in the significant magnetic fields generated during laser–plasma interactions or in magnetized fusion schemes. Accurate non-local and magnetized heat flow computations can be carried out using Vlasov–Fokker–Planck (VFP) simulations, but these are computationally expensive. There is, therefore, significant interest in using faster multi-group models to accurately calculate the non-local heat flow in magnetized plasmas. We benchmark two such multi-group models for calculating the heat flow, M1 and hybrid-AWBS-BGK, against diffusive models and full VFP simulations, before applying the models to realistic example test cases, both magnetized and unmagnetized. We find that the multi-group models generally perform very well for moderate non-localities up to kλmfp∼0.01, but the computational cost increases dramatically. hybrid-AWBS-BGK performs more effectively than M1 at high non-localities, up to kλmfp∼1, due to its adaptive solver and robust P1 closure, but tends to fail in very strong magnetic fields. Both codes are much faster than VFP simulations but are still slow in steep temperature gradients.

Arran, C. (ORCID:0000000286448118)↗
DOE OSTI2023

Inference of main ion particle transport coefficients with experimentally constrained neutral ionization during edge localized mode recovery on DIII-D

Abstract The plasma and neutral density dynamics after an edge localized mode are investigated and utilized to infer the plasma transport coefficients for the density pedestal. The Lyman-Alpha Measurement Apparatus (LLAMA) diagnostic provides sub-millisecond profile measurements of the ionization and neutral density and shows significant poloidal asymmetries in both. Exploiting the absolute calibration of the LLAMA diagnostic allows quantitative comparison to the electron and main ion density profiles determined by charge-exchange recombination, Thomson scattering and interferometry. Separation of diffusion and convection contributions to the density pedestal transport are investigated through flux gradient methods and time-dependent forward modeling with Bayesian inference by adaptation of the Aurora transport code and IMPRAD framework to main ion particle transport. Both methods suggest time-dependent transport coefficients and are consistent with an inward particle pinch on the order of 1 m s −1 and diffusion coefficient of 0.05 m 2 s −1 in the steep density gradient region of the pedestal. While it is possible to recreate the experimentally observed phenomena with no pinch in the pedestal, low diffusion in the core and high outward convection in the near scrape-off layer are required without an inward pedestal pinch.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Turbulent drifts of impurity ions as an explanation for anomalous radial transport in the far-SOL of DIII-D

Abstract Successful fusion reactor operation relies on minimal core contamination by impurities, otherwise too much power may be radiated and harm performance. This requires reliable predictions of impurity transport from the scrape-off layer (SOL) into the core, beyond the traditional ‘anomalous’ diffusion approach. We report a set of far-SOL tungsten transport simulations that demonstrate the role of turbulent drifts on radial impurity transport. A turbulent plasma background is simulated using the gyrokinetic SOL code Gkeyll. Tungsten ions are followed within the plasma background using only their drifts. We find that tungsten tends to travel radially outwards with velocities between v r = 300–1200 m s −1 primarily due to polarization drift. We also extract an anomalous radial diffusion coefficient that varies from D r anom = 5–20 m 2 s −1 . These results are compared to and agree with previous interpretive modeling results. We also show how the turbulent polarization drift can transport some tungsten ions from the wall inwards with effective pinch velocities up to 10 000 m s −1 . We conclude that turbulent drifts are a likely explanation for historically anomalous radial impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

SetBERT: the deep learning platform for contextualized embeddings and explainable predictions from high-throughput sequencing

MOTIVATION: High-throughput sequencing (HTS) is a modern sequencing technology used to profile microbiomes by sequencing thousands of short genomic fragments from the microorganisms within a given sample. This technology presents a unique opportunity for artificial intelligence to comprehend the underlying functional relationships of microbial communities. However, due to the unstructured nature of HTS data, nearly all computational models are limited to processing DNA sequences individually. This limitation causes them to miss out on key interactions between microorganisms, significantly hindering our understanding of how these interactions influence the microbial communities as a whole. Furthermore, most computational methods rely on post-processing of samples which could inadvertently introduce unintentional protocol-specific bias. RESULTS: Addressing these concerns, we present SetBERT, a robust pre-training methodology for creating generalized deep learning models for processing HTS data to produce contextualized embeddings and be fine-tuned for downstream tasks with explainable predictions. By leveraging sequence interactions, we show that SetBERT significantly outperforms other models in taxonomic classification with genus-level classification accuracy of 95%. Furthermore, we demonstrate that SetBERT is able to accurately explain its predictions autonomously by confirming the biological-relevance of taxa identified by the model. AVAILABILITY AND IMPLEMENTATION: All source code is available at https://github.com/DLii-Research/setbert. SetBERT may be used through the q2-deepdna QIIME 2 plugin whose source code is available at https://github.com/DLii-Research/q2-deepdna.

Ludwig, David W↗
DOE OSTI2025

Accelerating Control Systems with GitOps: A Path to Automation and Reliability

GitOps is a foundational approach for modernizing infrastructure by leveraging Git as the single source of truth for declarative configurations. The poster explores how GitOps transforms traditional control system infrastructure, services and applications by enabling fully automated, auditable, and version-controlled infrastructure management. Cloud-native and containerized environments are shifting the ecosystem not only in the IT industry but also within the computational science field, as is the case of CERN and Diamond Light Source among other Accelerator/Science facilities which are slowly shifting towards modern software and infrastructure paradigms. The ACORN project, which aims to modernize Fermilab’s control system infrastructure and software is implementing proven best-practices and cutting-edge technology standards including GitOps, containerization, infrastructure as code and modern data pipelines for control system data acquisition and the inclusion of AI/ML in our accelerator complex.

Gonzalez, M. [Fermilab]↗
DOE OSTI2021

Mini minutes for a mini CSEWG Virtual mini-CSEWG Meeting 16-19 August, 2021

This meeting is a collaboration meeting of the Cross Section Evaluation Working Group (CSEWG), a collaboration including members from various DOE NNSA and non-NNSA laboratories as well as several university groups that support these laboratory missions. The focus of this meeting will be preparations for the next major release of the ENDF/B nuclear data library. The ENDF/B library is the main neutronics library used in various US neutronics codes such as MCNP and SCALE for advance reactor work, safeguards, stockpile stewardship, etc.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Universal Utility Data Exchange (UUDEX) – Security and Administration: Cybersecurity of Energy Delivery Systems (CEDS) Research and Development

A critical component of the Universal Utility Data Exchange (UUDEX) approach is the integrated security contained within its processing. This document describes how that security is designed and expected to be implemented by UUDEX Implementations (U-Implementations), including the UUDEX Server (U-Server) and UUDEX Clients (U-Clients). The UUDEX security hierarchy consists of three levels: 1. The UUDEX Instance (U-Instance) itself, which sits at the top of the hierarchy and contains the U-Server, the UUDEX Identity Authority (U-Identity Authority), and the UUDEX Administrator (U-Administrator) functions; 2. A group of one or more UUDEX Participants (U-Participants) that present “organizations” that participate in the U-Instance and contains the UUDEX Administrator Participant (U-U-Administrator Participant) function; 3. A group of one or more UUDEX Endpoints (U-Endpoints) that represent the individual UUDEX Publish Clients (U-Publish Client) responsible for supplying data to the U-Instance that is consumed by UUDEX Subscriber Clients (U-Subscriber Clients). U-Endpoints can be either autonomous devices that publish and subscribe data such as data exchange servers found in supervisory control and data acquisition and energy management systems, or they can be tied to users of applications that, for example, submit DOE OE-417 disturbance reports. U-Participants and U-Endpoints can be organized into UUDEX Groups (U-Groups). Any number of U-Participants or U-Endpoints can be members of a U-Group. A given U-Participant or U-Endpoint can be a member of multiple U-Groups, but a U-Group cannot contain other U-Groups. For example, a U-Group could be created to contain all U-Participant Transmission Operators within the purview of a Reliability Coordinator, and another U-Group could be created to contain all U-Participant Generator Operators within the purview of a Reliability Coordinator. U-Participants that are both Transmission Operators and Generator Operators would be members of both U-Groups. U-Participants, U-Endpoints, and U-Groups are used in the access control structures to provide access to individual UUDEX Subjects (U-Subjects). U-Groups are created by the U-Administrator and are managed by the U-Administrator or the designated U-Group Managers. U-Endpoints can be assigned UUDEX Roles (U-Roles) that can be used to further restrict access. U-Roles are assigned to individual U-Endpoints. For example, a U-Role of “Security Analyst” could be used to restrict which U-Endpoints can publish or subscribe security incident reports and vulnerability notifications, while a U-Role of “Transmission Planner” can be used to restrict which U-Endpoints can publish power system model updates. U-Role definitions are created by the U-Administrator, but the U-Roles are assigned to U-Endpoints by their respective UUDEX Participant Administrators (U-Participant Administrator). Because all information required to make security decisions is either included within the U-Endpoint’s X.509 digital certificate or stored in a datastore on the U-Server, all security decisions are performed and enforced within the U-Server. This reduces the complexity of the U-Client code and minimizes the chance for compromise of the integrity of the UUDEX security features.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Modeling of the Jupiter Lase Facility Using Virtual Beamline ++ (Vbl++) Code

The Jupiter Laser Facility (JLF) at Lawrence Livermore National Laboratory (LLNL) is a mid-scale laser facility, with three different laser platforms: Janus (Target Area 1), Titan, and COMET. To help the facility and its users better understand the laser performance, a model of Janus and Titan has been developed using the Virtual Beamline ++ code (VBL++), which was initially developed by LLNL to model the National Ignition Facility beamlines and predict their performance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

“Development of Interactive Chemistry in E3SM v2” Final Report (3/1/2023-2/29/2024)

We are pleased to report that the proposed milestones and tasks under this subcontract, or those modified subsequently by the Technical Representative based on project needs, have been met. In general, the accomplishments under this subcontract exceeded the originally proposed milestones and effort, and went from code development to scientific results and high-profile publications, to mission critical software development for the E3SM project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Dielectric screening and vacancy formation for large neutral and charged Si n H m ( n > 1500 ) nanocrystals using real-space pseudopotentials

A commonly used procedure for computing the properties of defects in crystalline materials is to consider a large supercell that includes the defect of interest. This is a straightforward technique as standard energy band codes can be used for such computations. For neutral defects, the only impediment of such an approach is to avoid defect-defect interactions between adjoining cells. However, this procedure can be complex if the defect of interest is charged as the system at large contains Coulombic divergences. Moreover, some have recently argued that the conventional definition of formation energies for charged defects cannot be reconciled with statistical mechanics. Here, we focus on an alternative approach. We consider large nanocrystals wherein a charged defect can be placed. Since the system is confined, a charged defect within the nanocrystal does not result in a Coulombic divergence. The chief impediment is computational, i.e., while no defect-defect or Coulombic divergences are present, the nanocrystal must be sufficiently large to allow the system to properly replicate a bulklike configuration. With the development of new algorithms and hardware advances, computations for systems of sufficient size to address this issue are feasible. In particular, we solve the Kohn-Sham equation in real space using pseudopotential-density-functional theory for large silicon nanocrystals, which contain thousands of atoms. Further, we focus on (i) the screening of a point charge and (ii) the formation of a charged vacancy in hydrogen-terminated silicon nanocrystals. This approach allows us to examine the role of quantum confinement in addition to exploring the bulk limit. Comparisons to other methods confirm the viability of this approach.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Validation and Independent Uncertainty Analysis of the MIX-SOL-THERM-003 ICSBEP Benchmark

The International Criticality Safety Benchmark Evaluation Project (ICSBEP) was started in 1992 by the United States Department of Energy and later in 1995 became an international project with contributions from 22 countries. The project is now organized by the OECD (Organisation for Economic Co-operation and Development) Nuclear Energy Agency (NEA). In its most recent iteration, the ICSBEP handbook contains over five thousand evaluations of critical, near-critical, and subcritical experiments conducted in facilities all around the world. These benchmarks serve as valuable information for criticality safety engineers who can use them to validate calculation techniques and establish minimum subcritical margins for operations with fissionable materials. The benchmarks in the handbook are categorized by their fissile material composition, material form (oxide, solution, or metal), and fission energy spectra. This is especially useful for those looking for benchmarks similar to a system they are working on to compare methods and identify trends. The ICSBEP Handbook Uncertainty Guide is document outlining recommended practices and methods for determining uncertainties in these benchmarks. Quantifying these uncertainties thoroughly is crucial as it allows a higher degree of confidence that data used from them is valid and relevant. The guide stresses the importance of a thorough and well documented uncertainty analysis when evaluating an experiment. All measured values of a system, whether they be dimensions or material compositions, have a certain amount of uncertainty associated with them and can be analyzed one by one to determine their effects on the system. Many evaluated benchmarks in the handbook present this in detail, however some do not, mostly earlier evaluations performed in the 1990’s and early 2000’s. Recently at Los Alamos National Laboratory (LANL), the Nuclear Criticality Safety Division (NCSD) of LANL has been validating MCNP6.2 ® input files of criticality benchmarks for use by Whisper, a criticality safety code developed at LANL. This effort is also part of the OECD NEA Working Party on International Nuclear Data Evaluation Co-operation (WPEC) Subgroup 45, also known as Validation of Nuclear Data Libraries (VaNDaL). The goal of VaNDaL is to compile a set of validated simulation inputs for use in validating nuclear data and simulation codes. As part of these efforts, one of the benchmarks reviewed was the MIX SOL-THERM-003 ICSBEP benchmark. This paper provides an independent uncertainty analysis of this benchmark experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

ANTARES: a new IRSN computational chain for safety analysis

Since some years, the IRSN has developed a new project to simulate accurately the interaction between neutronic and thermal hydraulic phenomena. This project called ANTARES (Advanced Neutronics and Thermal-hydraulic for the Analysis of Reactor Safety) couples the thermal hydraulic code CATHARE-3 with the 3D neutronic nodal code PARCS. The purpose of this project is to take the advantage of the complementary of these 2 codes, by the knowledge of the 3D core power distribution in CATHARE-3 (computed by PARCS) and improving the thermal hydraulic boundaries conditions in PARCS. In this way, the IRSN diversifies and enlarges its available tools to perform safety analysis with improved accuracy. The current technique usually adopted in France for the safety demonstrations for Reactivity Insertion Accident (RIA) is the so-called 'conservative' approach, which consists of reducing all the feedback (Doppler and moderator effects) and in modifying some physical quantities in such a way to increase a power peak. For this reason, these facilities ('penalties') have been implemented in ANTARES. In this paper we will give two examples of accidental transients that can be simulated with ANTARES: a Main Steam Line Break (MSLB) and an inadvertent boron dilution event. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Ancestral archaea expanded the genetic code with pyrrolysine

The pyrrolysyl-tRNA synthetase (PylRS) facilitates the cotranslational installation of the 22nd amino acid pyrrolysine. Owing to its tolerance for diverse amino acid substrates, and its orthogonality in multiple organisms, PylRS has emerged as a major route to install noncanonical amino acids into proteins in living cells. Recently, a novel class of PylRS enzymes was identified in a subset of methanogenic archaea. Enzymes within this class (ΔPylSn) lack the N-terminal tRNA-binding domain that is widely conserved amongst PylRS enzymes, yet remain active and orthogonal in bacteria and eukaryotes. In this study, we use biochemical and in vivo UAG-readthrough assays to characterize the aminoacylation efficiency and substrate spectrum of a ΔPylSn class PylRS from the archaeon Candidatus Methanomethylophilus alvus. We show that, compared with the full-length enzyme from Methanosarcina mazei, the Ca. M. alvus PylRS displays reduced aminoacylation efficiency but an expanded amino acid substrate spectrum. To gain insight into the evolution of ΔPylSn enzymes, we performed molecular phylogeny using 156 PylRS and 105 pyrrolysine tRNA (tRNA Pyl ) sequences from diverse archaea and bacteria. This analysis suggests that the PylRS•tRNA Pyl pair diverged before the evolution of the three domains of life, placing an early limit on the evolution of the Pyl-decoding trait. Furthermore, our results document the coevolutionary history of PylRS and tRNA Pyl and reveal the emergence of tRNA Pyl sequences with unique A73 and U73 discriminator bases. The orthogonality of these tRNA Pyl species with the more common G73-containing tRNA Pyl will enable future efforts to engineer PylRS systems for further genetic code expansion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

IC reports (2021)

We requested HPC support to continue research on the seismic waves generated by impacts. We used the following codes: (1) the Hybrid Optimization Software Suite (HOSS), developed at LANL. HOSS is based on a combined Finite and Discrete Element Method (FDEM). New material models are developed for the sedimentary rocks. HOSS has been recently benchmarked to iSale and FLAG codes. SPECFEM3D is an open-source code developed since the last 90s. It won the Gordon Bell award for best performance in 2003, was finalist again in 2008 for a run at 0.16 petaflops on 149,784 cores on the ‘Jaguar’ Cray system at Oak Ridge National Laboratory. It also won the BULL Joseph Fourier supercomputing award in 2010.; SW4 is a 4th-order finite difference code developed at LLNL which is currently actively developed to handle complex 3D models and to be ported on future exascale platforms. We assessed our need to a total of 3.9M CPU-hrs for year 1 and and 3.1 M for year 2.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Implementation of a Drift-Flux Model in SAM for Modeling of Passively Transported GAS in Molten Salt Reactors

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is further developing the Multiphysics Object Oriented Simulation Environment (MOOSE) system thermal hydraulics (T/H) code, System Analysis Module (SAM), to include modeling of molten salt reactor (MSR) designs. MSR designs in which the fuel is dissolved in the salt coolant poses unique modeling and simulation (M&S) challenges because fission products—some of which will be noncondensable gasses—travel with the coolant, thus allowing for deposition and interaction throughout the entire system flow loop. This behavior must be modeled to successfully understand MSR system behavior during normal and off-normal operation events. In past work, SAM has been extended to include a passive species transport capability for modeling of MSR fission products that travel with the salt coolant. This work extends the model to include a drift-flux model to capture the tendency of the gas phase to move at a different velocity than the bulk liquid. Additional constitutive models have also been added for calculation of gas volume fraction, bubble radius, and interfacial area, which will be important for calculation of species phase migration in future coupling activities. Testing has been performed to ensure that the new model agrees with expected values for several simple models, but additional future work will be required to test the model for more complex geometries that will be needed for full-scale MSR simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Computation and Comparison of the Stable Northeastern US Marine Boundary Layer

In this work we investigate the behavior of stable marine boundary layers located near the coast of the Northeastern United States. Using the ExaWind large eddy simulation (LES) codes, three stable atmospheric conditions were chosen to match the Cape Wind measurements of Archer et al. with wind speeds of 5 m/s, 10 m/s, and 15 m/s at the 20 m measurement height. The behavior of the stable boundary layers, including mean flow quantities and turbulent statistics, are examined and compared to previous computations of the neutral and unstable offshore boundary layer at the same location. This study also examines the domain and mesh requirements necessary to capture the turbulent scales for the stable offshore boundary layers. Finally, we compare solutions computed using both Nalu-Wind and AMR-Wind solvers, and compare their predicted solutions and performance in this study.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2022

Farpoint: A High-resolution Cosmology Simulation at the Gigaparsec Scale

Abstract In this paper we introduce the Farpoint simulation, the latest member of the Hardware/Hybrid Accelerated Cosmology Code (HACC) gravity-only simulation family. The domain covers a volume of (1000 h −1 Mpc) 3 and evolves close to two trillion particles, corresponding to a mass resolution of m p ∼ 4.6 × 10 7 h −1 M ⊙ . These specifications enable comprehensive investigations of the galaxy–halo connection, capturing halos down to small masses. Further, the large volume resolves scales typical of modern surveys with good statistical coverage of high-mass halos. The simulation was carried out on the GPU-accelerated system Summit, one of the fastest supercomputers currently available. We provide specifics about the Farpoint run and present an initial set of results. The high mass resolution facilitates precise measurements of important global statistics, such as the halo concentration–mass relation and the correlation function down to small scales. Selected subsets of the simulation data products are publicly available via the HACC Simulation Data Portal.

79 ASTRONOMY AND ASTROPHYSICS↗