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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Classical Non-Markovian Noise in Symmetry-Preserving Quantum Dynamics

In quantum dynamics, symmetries are vital for identifying and assessing conserved quantities that govern the evolution of a quantum system. When promoted to the open quantum system setting, dynamical symmetries can be negatively altered by system-environment interactions, thus, complicating their analysis. Previous work on noisy symmetric quantum dynamics has focused on the Markovian setting, despite the ubiquity of non-Markovian noise in a number of widely used quantum technologies. Here, in this Letter, we develop a framework for quantifying the impact of non-Markovian noise on symmetric quantum evolution via root space decompositions and the filter function formalism. We demonstrate analytically that symmetry-preserving noise maintains the symmetric subspace, while nonsymmetric noise leads to highly specific leakage errors that are block diagonal in the symmetry representation. We support our findings with numerical studies of a transverse-field Ising model and a quantum error detecting code subject to spatiotemporally correlated multiaxis noise. Our results are broadly applicable, providing new analytic insights into the control and characterization of open quantum system dynamics.

decoherence↗
DOE OSTI2020

Spitfire

Spitfire is a Python and C++ code for solving canonical flame structure problems to build tabulated chemistry models for the simulation of large-scale fires. Spitfire is a modular software with performance critical code written in C++, abstract numerical methods provided in Python, and chemistry tabulation techniques made convenient for high-level use by non-expert users. It facilitates both high-level usage and novel research in tabulated chemistry approaches. SAND2019-11791 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Armstrong, Elizabeth↗
DOE OSTI2021

A method for achieving complete microbial genomes and improving bins from metagenomics data

Metagenomics facilitates the study of the genetic information from uncultured microbes and complex microbial communities. Assembling complete genomes from metagenomics data is difficult because most samples have high organismal complexity and strain diversity. Some studies have attempted to extract complete bacterial, archaeal, and viral genomes and often focus on species with circular genomes so they can help confirm completeness with circularity. However, less than 100 circularized bacterial and archaeal genomes have been assembled and published from metagenomics data despite the thousands of datasets that are available. Circularized genomes are important for (1) building a reference collection as scaffolds for future assemblies, (2) providing complete gene content of a genome, (3) confirming little or no contamination of a genome, (4) studying the genomic context and synteny of genes, and (5) linking protein coding genes to ribosomal RNA genes to aid metabolic inference in 16S rRNA gene sequencing studies. We developed a semi-automated method called Jorg to help circularize small bacterial, archaeal, and viral genomes using iterative assembly, binning, and read mapping. In addition, this method exposes potential misassemblies from k-mer based assemblies. We chose species of the Candidate Phyla Radiation (CPR) to focus our initial efforts because they have small genomes and are only known to have one ribosomal RNA operon. In addition to 34 circular CPR genomes, we present one circular Margulisbacteria genome, one circular Chloroflexi genome, and two circular megaphage genomes from 19 public and published datasets. We demonstrate findings that would likely be difficult without circularizing genomes, including that ribosomal genes are likely not operonic in the majority of CPR, and that some CPR harbor diverged forms of RNase P RNA. Code and a tutorial for this method is available at https://github.com/lmlui/Jorg and is available on the DOE Systems Biology KnowledgeBase as a beta app.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Understanding Power and Energy Utilization in Large Scale Production Physics Simulation Codes

Power is an often-cited reason for moving to advanced architectures on the path to Exascale computing. This is due to the practical concern of delivering enough power to successfully site and operate these machines, as well as concerns over energy usage while running large simulations. Since accurate power measurements can be difficult to obtain, processor thermal design power (TDP) is a possible surrogate due to its simplicity and availability. However, TDP is not indicative of typical power usage while running simulations. Using commodity and advance technology systems at Lawrence Livermore National Laboratory (LLNL) and Sandia National Laboratory, we performed a series of experiments to measure power and energy usage in running simulation codes. These experiments indicate that large scale LLNL simulation codes are significantly more efficient than a simple processor TDP model might suggest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

An Extension of the Athena++ Code Framework for Radiation-magnetohydrodynamics in General Relativity Using a Finite-solid-angle Discretization

We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical procedure in detail, verifying its correctness with a suite of tests. Motivated in particular by black hole accretion in the high-accretion-rate, thin-disk regime, we demonstrate the application of the method to this problem. With excellent scaling on flagship computing clusters, the port of the algorithm to the GPU-enabled AthenaK code now allows the simulation of many previously intractable radiation-GRMHD systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Methods for Incorporating Model Uncertainty into Exoplanet Atmospheric Analysis

A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Land-based wind turbines with flexible rail-transportable blades – Part 2: 3D finite element design optimization of the rotor blades

Abstract. Increasing growth in land-based wind turbine blades to enable higher machine capacities and capacity factors is creating challenges in design, manufacturing, logistics, and operation. Enabling further blade growth will require technology innovation. An emerging solution to overcome logistics constraints is to segment the blades spanwise and chordwise, which is effective, but the additional field-assembled joints result in added mass and loads, as well as increased reliability concerns in operation. An alternative to this methodology is to design slender flexible blades that can be shipped on rail lines by flexing during transport. However, the increased flexibility is challenging to accommodate with a typical glass-fiber, upwind design. In a two-part paper series, several design options are evaluated to enable slender flexible blades: downwind machines, optimized carbon fiber, and active aerodynamic controls. Part 1 presents the system-level optimization of the rotor variants as compared to conventional and segmented baselines, with a low-fidelity representation of the blades. The present work, Part 2, supplements the system-level optimization in Part 1 with high-fidelity blade structural optimization to ensure that the designs are at feasible optima with respect to material strength and fatigue limits, as well as global stability and structural dynamics constraints. To accommodate the requirements of the design process, a new version of the Numerical Manufacturing And Design (NuMAD) code has been developed and released. The code now supports laminate-level blade optimization and an interface to the International Energy Agency Wind Task 37 blade ontology. Transporting long, flexible blades via controlled flapwise bending is found to be a viable approach for blades of up to 100 m. The results confirm that blade mass can be substantially reduced by going either to a downwind design or to a highly coned and tilted upwind design. A discussion of active and inactive constraints consisting of material rupture, fatigue damage, buckling, deflection, and resonant frequencies is presented. An analysis of driving load cases revealed that the downwind designs are dominated by loads from sudden, abrupt events like gusts rather than fatigue. Finally, an analysis of carbon fiber spar caps for downwind machines finds that, compared to typical carbon fibers, the use of a new heavy-tow carbon fiber in the spar caps is found to yield between 9 % and 13 % cost savings.

17 WIND ENERGY↗
DOE OSTI2023

ClimSim: A large multi-scale dataset for hybrid physics-ML climate emulation

Modern climate projections lack adequate spatial and temporal resolution due to computational constraints. A consequence is inaccurate and imprecise predictions of critical processes such as storms. Hybrid methods that combine physics with machine learning (ML) have introduced a new generation of higher fidelity climate simulators that can sidestep Moore’s Law by outsourcing compute-hungry, short, high-resolution simulations to ML emulators. However, this hybrid ML-physics simulation approach requires domain-specific treatment and has been inaccessible to MLexperts because of lack of training data and relevant, easy-to-use workflows. Wepresent ClimSim, the largest-ever dataset designed for hybrid ML-physics research. It comprises multi-scale climate simulations, developed by a consortium of climate scientists and ML researchers. It consists of 5.7 billion pairs of multivariate input and output vectors that isolate the influence of locally-nested, high-resolution, high-fidelity physics on a host climate simulator’s macro-scale physical state. The dataset is global in coverage, spans multiple years at high sampling frequency, and is designed such that resulting emulators are compatible with downstream coupling into operational climate simulators. We implement a range of deterministic and stochastic regression baselines to highlight the ML challenges and their scoring. The data (https://huggingface.co/datasets/LEAP/ClimSim_high-res2) and code(https://leap-stc.github.io/ClimSim)arereleasedopenlytosupport the development of hybrid ML-physics and high-fidelity climate simulations for the benefit of science and society.

artificial intelligence, machine learning↗
DOE OSTI2023

Isotope dependence of transport in ST40 hot ion mode plasmas

The ST40 compact, high-field spherical tokamak, operating at 2.1 T and with 1.8 MW of neutral beam heating, achieved central carbon impurity ion temperatures in excess of 10 keV, surpassing their business milestone of 100 M °C (8.6 keV). The high temperature discharges were in the hot ion mode, with $T_i \gg T_e$, and they were achieved in both hydrogen and deuterium plasmas with deuterium neutral beam injection. In order to achieve the high temperature scenarios, careful wall conditioning and scenario optimization was carried out in $H^0 \rightarrow H^+$, $D^0 \rightarrow H^+$, and finally in $D^0 \rightarrow D^+$ plasmas. The TRANSP transport code was employed to study the dependence of confinement and transport on isotopic mass, and the conditions that led to the high measured central ion temperatures in the $D^0 \rightarrow H^+$ and $D^0 \rightarrow D^+$ plasmas. The kinetic profiles input into TRANSP were inferred from line-of-sight and limited radial measurements as well as consistency with a number of other experimental constraints. The TRANSP results first showed that the main species central ion temperature was only 0.5–1 keV lower than the measured carbon impurity temperature, and thus in the high performance plasmas also surpassed the 100 M °C level. TRANSP also showed that while the electron thermal conduction loss was dominant, reductions in central ion transport, and an effective decoupling of the ions from the electrons, led to an increase in central ion temperatures with increasing plasma mass. In fact, the associated confinement times exhibited a strong dependence on the isotopic mass of the thermal plasma. These preliminary results are foundations for dedicated experiments with full kinetic profile measurements on ST40 in the next run campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

High-performance and high-fidelity Monte Carlo solutions to the BEAVRS benchmark

The BEAVRS (Benchmark for Evaluation and Validation of Reactor Simulation) benchmark is solved by PRAGMA, the GPU-based continuous energy Monte Carlo code. The resulting solutions are comprised of the detailed simulation results of two cycles, each of which consists of the zero power physics test (ZPPT) and the core depletion calculations. The ZPPT consists of characteristic parameters, such as critical boron concentration (CBC), control rod bank worth, isothermal temperature coefficients, and assembly-wise detector signal, which are compared with measured data provided by the benchmark administration. The core depletion calculations were performed for both the hot full power and the load follow modes, and the comparison was made with the measured or deduced CBCs and assembly-wise detector signals. In the load follow calculations, the operating power history was approximately applied to simulate the real operation as closely as possible. PRAGMA performed the various calculations with a tremendous number of histories ranging up to hundreds of millions per cycle, exploiting GPUs' massively parallel performance. The load follow run time was shorter than 16 hours on a single rack of computing nodes mounded with 24 gaming GPUs with a remarkable agreement with the measurements for most comparisons. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Synthesis of magnesiowüstite nanocrystallites embedded in an amorphous silicate matrix via low energy multiple ion implantations

The synthesis process is presented for experimentally simulating modifications in cosmic dust grains using sequential ion implantations or irradiations followed by thermal annealing. Cosmic silicate dust analogues were prepared via implantation of 20–80 keV Fe - , Mg - , and O - ions into commercially available p-type silicon (100)wafers. The as-implanted analogues are amorphous with a Mg/(Fe+Mg) ratio of 0.5 tailored to match theoretical abundances in circumstellar dusts. Before the ion implantations were performed, Monte-Carlo-based ion-solid interaction codes were used to model the dynamic redistribution of the implanted atoms in the silicon substrate. 600 keV helium ion irradiation was performed on one of the samples before thermal annealing. Two samples were thermally annealed at a temperature appropriate for an M-class stellar wind, 1000 K, for 8.3 h in a vacuum chamber with a pressure of 1 x 10 -7 torr. The elemental depth profiles were extracted utilizing Rutherford Backscattering Spectrometry (RBS) in the samples before and after thermal annealing. X-ray diffraction (XRD)analysis was employed for the identification of various phases in crystalline minerals in the annealed analogues. Transmission electron microscopy (TEM) analysis was utilized to identify specific crystal structures. RBS analysis shows redistribution of the implanted Fe, Mg, and O after thermal annealing due to incorporation into the crystal structures for each sample type. XRD patterns along with TEM analysis showed nanocrystalline Mg and Fe oxides with possible incorporation of additional silicate minerals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Simulations of divertor plasmas with inverse sheaths

The effect of strong electron emission from material surfaces has been proposed to form an “inverse sheath”: a region with a positive potential relative to the near-wall plasma which prevents the flow of ions to the wall [M. D. Campanell, “Negative plasma potential relative to electronemitting surfaces,” Phys. Rev. E. 88, 033103 (2013); M. D. Campanell and M. V. Umansky, “Strongly emitting surfaces unable to float below plasma potential,” Phys. Rev. Lett. 116, 1–5 (2016); M. D. Campanell and G. R. Johnson, “Thermionic cooling of the target plasma to a sub-ev temperature,” Phys. Rev. Lett. 122, 1–5 (2019)]. We assess the viability of this regime in a tokamak device using the 2D edge plasma transport code UEDGE [T. Rognlien et al ., “A fully implicit, time dependent 2-D fluid code for modeling tokamak edge plasmas,” J. Nucl. Mater. 196–198, 347–351 (1992)]. Since the UEDGE code does not consider the sheath region directly, we apply boundary conditions at the divertor targets which emulate the physics of both “standard” and “inverse” sheath regimes [R. Masline et al ., “Influence of the inverse sheath on divertor plasma performance in tokamak edge plasma simulations,” Contrib. Plasma Phys. 60, e201900097 (2020)]. Using these boundary conditions, we perform scoping studies to assess plasma parameters near the target by varying the density at the core-edge interface. We observe a smooth transition in the resultant profiles of plasma parameters for the standard sheath, and a bifurcation across the simulation set for plasmas with an inverse sheath. The cause of this bifurcation is assessed by performing the parameter scan both with and without impurity radiation; we observe that the bifurcation persists in both cases, indicating that this bifurcation is caused by plasma recombination.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effect of magnetic field ripple on parallel electron transport during microwave plasma heating in the Proto-MPEX linear plasma device

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a linear plasma device located at Oak Ridge National Laboratory to develop the plasma source concept for the Materials Plasma Exposure eXperiment (MPEX). Recent experiments have demonstrated the heating of electrons in a high-density deuterium plasma (n e ~ 5 × 10 19 m -3 ) via electron Bernstein waves and source plasma production via helicon waves in this linear configuration. Moreover, experimental observations suggest that the magnetic ripple adversely affects the parallel transport of the heated electrons toward the target where material samples are to be exposed to the plasma. In order to understand the transport process during microwave application, a test-particle Monte-Carlo (MC) code has been developed that incorporates the effects of Coulomb collisions, magnetic mirror adiabatic trapping and electron cyclotron interaction via a quasilinear RF heating operator. MC calculations indicate that the absence of 2nd turning points along the trajectory of cyclotron heated electrons significantly reduces adiabatic trapping and maximizes power transport to the Target surface. Thermalization of cyclotron heating electrons is analyzed with the MC code and it is found that for conditions relevant to ProtoMPEX, a significant fraction of the absorbed heating power is coupled to the Target surface via ballistic fast electrons; however, under certain conditions, up to 60% of the absorbed RF power can be dissipated via thermalization of the fast electrons on the background plasma. Results are compared with experimental measurements. Finally, the implications of the results are discussed in the context of the upcoming MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

CHARM-SYCL & IRIS: A Tool Chain for Performance Portability on Extremely Heterogeneous Systems

Performance portability is becoming crucial as high-performance computing systems become increasingly heterogeneous. We have many options for CPUs and accelerators (e.g., GPUs) but also for non-Von Neumann architectures such as field-programmable gate arrays. This paper presents the CHARM-SYCL unified programming environment for multiple accelerator types as a performance-portable programming environment. It uses the IRIS library developed at Oak Ridge National Laboratory as the back end accelerator runtime. IRIS has a high-performance scheduler to distribute tasks across accelerators. This design allows us to run an application from the same source on multiple systems with multiple configurations. We provide three types of portability with CHARM-SYCL: Portable Workflow, Compiler and Runtime Portability, and Application and Performance Portability. We implement a Monte Carlo simulation benchmark code on the CHARM-SYCL execution environment and demonstrate that our programming environment can accommodate extremely heterogeneous systems.

Fujita, Norihisa↗
DOE OSTI2024

SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations

Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world’s largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach is used ubiquitously across scientific domains, and in many cases leads to a set of ordinary differential equations (ODEs) that need to be solved as part of a larger “outer-loop” time-stepping approach. The SUNDIALS library provides a plethora of robust time integration algorithms for solving ODEs, and the U.S. Department of Energy Exascale Computing Project (ECP) has supported its extension to applications on exascale-capable computing hardware. In this paper, we highlight some SUNDIALS capabilities and its deployment in combustion and cosmology application codes (Pele and Nyx, respectively) where operator splitting gives rise to numerous, small ODE systems that must be solved concurrently.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

TNSL improvements in FUDGE [Slides]

FUDGE TNSL processing capabilities are improving. New development focuses on reading and writing GNDS, but FUDGE also supports writing to ACE and ENDL. Work provides a chance to compare in-depth with other codes. TNSL processing is mostly complete, but MCGIDI sampling has room for improvement. There are plans to revisit interpolation and to improve strategies for sampling coherent/incoherent elastic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Spacecraft surface charging as a function of material properties

Spacecraft material behavior plays a very important role in space missions. Spacecraft immersed in plasma get charged by absorbing plasma particles and by emitting electrons from spacecraft surfaces via photoelectron and secondary electron emission. Spacecraft charging depends heavily on material properties such as work function, secondary electron yield, dielectric constant, and electric conductivity among other. Material properties are typically assumed to be static in charging models. However, it is well known that this is not the case in space. This makes spacecraft charging predictions very challenging. Material properties are well characterized before the spacecraft is put in orbit through characterization in the lab under clean conditions. However, once in space, material properties change due to the harsh and very dynamic space environment. We present a new capability to predict material behavior in space from first-principles modeling. The ongoing effort seeks to couple material models, density functional theory (DFT) and molecular dynamic (MD) codes, with environment models, plasma kinetic codes. This preliminary study will show results of surface charging as a function of material work function, dielectric constant, and conductivity.

36 MATERIALS SCIENCE↗