Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation kernel”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Surface enrichment dictates block copolymer orientation

Orientation of block copolymer (BCP) morphology in thin films is critical to applications as nanostructured coatings. Despite being well-studied, the ability to control BCP orientation across all possible block constituents remains challenging. Here, in this study, we deploy coarse-grained molecular dynamics simulations to study diblock copolymer ordering in thin films, focusing on chain makeup, substrate surface energy, and surface tension disparity between the two constituent blocks. We explore the multi-dimensional parameter space of ordering using a machine-learning approach, where an autonomous loop using a Gaussian process (GP) control algorithm iteratively selects high-value simulations to compute. The GP kernel was engineered to capture known symmetries. The trained GP model serves as both a complete map of system response, and a robust means of extracting material knowledge. We demonstrate that the vertical orientation of BCP phases depends on several counter-balancing energetic contributions, including entropic and enthalpic material enrichment at interfaces, distortion of morphological objects through the film depth, and of course interfacial energies. BCP lamellae are found more resistant to these effects, and thus more robustly form vertical orientations across a broad range of conditions; while BCP cylinders are found to be highly sensitive to surface tension disparity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Post-irradiation 1600°C Heating Test of AGR-1 UCO Fuel Kernels

Five bare kernels were selected from the U.S. Advanced Gas Reactor (AGR) 1 irradiation experiment Compact 5-3-1 to perform post-irradiation safety tests. The safety test involved heating the kernels in the Fuel Accident Condition Simulator (FACS) furnace in the inert atmosphere to a peak temperature of 1600°C and isothermally holding at this temperature for about 47 hours while collecting fission products released. This test was to assess the retention of fission products in bare kernels without the effects of the other TRISO layers (buffer, IPyC, SiC, and OPyC) or the graphitic matrix material. The bare kernels released nearly 100% of cesium and antimony, while they were able to maintain about 30% of europium, 50% of strontium, and the majority of cerium and ruthenium. In addition, about 48% of the calculated inventory of Kr-85 released during the test, indicating that kernels were capable of retaining a considerable fraction of fission gas K-85 during irradiation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Computational Fluid Dynamics Modeling of Low Temperature Ignition Processes From a Nanosecond Pulsed Discharge at Quiescent Conditions

Recent interest in nonequilibrium plasma discharges as sources of ignition for the automotive industry has not yet been accompanied by the availability of dedicated models to perform this task in computational fluid dynamics (CFD) engine simulations. The need for a low-temperature plasma (LTP) ignition model has motivated much work in simulating these discharges from first principles. Most ignition models assume that an equilibrium plasma comprises the bulk of discharge kernels. LTP discharges, however, exhibit highly nonequilibrium behavior. In this work, a method to determine a consistent initialization of LTP discharge kernels for use in engine CFD codes like CONVERGE is proposed. The method utilizes first principles discharge simulations. Such an LTP kernel is introduced in a flammable mixture of air and fuel, and the subsequent plasma expansion and ignition simulation is carried out using a reacting flow solver with detailed chemistry. Finally, the proposed numerical approach is shown to produce results that agree with experimental observations regarding the ignitability of methane-air and ethylene-air mixtures by LTP discharges.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

STEPS: A Portable Numerical Simulation Toolkit for Electrical Power System Dynamic Studies

Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.

42 ENGINEERING↗
DOE OSTI2022

Using kernel-based statistical distance to study the dynamics of charged particle beams in particle-based simulation codes

Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. Here, the focus is on kernel-based methods such as maximum mean discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Using Numba for GPU acceleration of Neutron Beamline Digital Twins

Digital twins of neutron instruments using Monte Carlo ray tracing have proven to be useful in neutron data analysis and verifying instrument and sample designs. However, these simulations can become quite complex and computationally demanding with tens of billions of neutrons. In this paper, we present a GPU accelerated version of MCViNE using Python and Numba to balance user extensibility with performance. Numba is an open-source just-in-time (JIT) compiler for Python using LLVM to generate efficient machine code for CPUs and GPUs with NVIDIA CUDA. The JIT nature of Numba allowed complex instrument kernels to be generated easily. Initial simulations have shown a speedup between 200-1000x over the original CPU implementation. The performance gain with Numba enables more sophisticated data analysis and impacts neutron scattering science and instrument design.

Kendrick, Coleman↗
DOE OSTI2025

Physics and Components syntax to enable a systems-based approach to multiphysics

Simulations in MOOSE have traditionally used kernel and boundary condition classes to describe the equations. Downstream applications leveraged a system called Actions to define a pre-packaged discretization of the equations they solve. Unfortunately, the Action base class was very limited, and most applications implemented the same concepts in their Actions. This led to an increased maintenance burden and a reduction in coupling opportunities, save for the use of MultiApps which renders each input mostly independent. With the introduction of multi-system capabilities in MOOSE, there is growing interest in defining entire simulations of complex multiphysics systems in a single input file. By introducing a new Physics system, with its dedicated syntax and a new base class providing wide-ranging capabilities, we are now able to define multiple equations in a single input file in a compact and user-friendly way. With new interactions between Physics and the Component system, these equations can be defined on each component of a complex system. In this talk, we will present the capabilities of these new systems, their interactions, and how to define complex systems multiphysics simulations with Physics and Components.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Evaluating Climate Model Simulations of the Radiative Forcing and Radiative Response at Earth’s Surface

In this work, we analyze the radiative forcing and radiative response at Earth’s surface, where perturbations in the radiation budget regulate the atmospheric hydrological cycle. By applying a radiative kernel-regression technique to CMIP5 climate model simulations where CO 2 is instantaneously quadrupled, we evaluate the intermodel spread in surface instantaneous radiative forcing, radiative adjustments to this forcing, and radiative responses to surface warming. The cloud radiative adjustment to CO 2 forcing and the temperature-mediated cloud radiative response exhibit significant intermodel spread. In contrast to its counterpart at the top of the atmosphere, the temperature-mediated cloud radiative response at the surface is found to be positive in some models and negative in others. Also, the compensation between the temperature-mediated lapse rate and water vapor radiative responses found in top-of-atmosphere calculations is not present for surface radiative flux changes. Instantaneous radiative forcing at the surface is rarely reported for model simulations; as a result, intermodel differences have not previously been evaluated in global climate models. We demonstrate that the instantaneous radiative forcing is the largest contributor to intermodel spread in effective radiative forcing at the surface. We also find evidence of differences in radiative parameterizations in current models and argue that this is a significant, but largely overlooked, source of bias in climate change simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

ORNL Analysis of Leach-Burn-Leach Round-Robin Test Samples

An international round-robin test to examine the consistency in leach-burn-leach (LBL) analysis of tristructural-isotropic- (TRISO-) coated particle fuel was conducted by three research organizations from the Generation IV International Forum member countries of the People’s Republic of China, the Republic of Korea, and the United States of America. Two sets of round-robin test samples were exchanged for analysis. One set of samples consisted of a series of nonuranium-bearing, TRISO-coated zirconium dioxide particles seeded with up to four depleted uranium-bearing, TRISO-coated uranium dioxide (UO 2 ) particles, which had intentionally damaged coating layers to simulate either particles with either exposed-kernel defects (i.e., particles with a cracked TRISO coating that should be detected during preburn leaching) or particles with silicon carbide (SiC) defects (i.e., particles with an intact pyrocarbon coating and a hole in the SiC layer that should be detected during postburn leaching). These simulated samples also contained added powder with known quantities of impurities from a coal standard. The other sample set consisted of representative sublots of UO 2 -TRISO particles fabricated in a production-scale coater, except they all contained depleted uranium instead of enriched uranium. In this report, the methodology used at Oak Ridge National Laboratory to conduct LBL analysis of the round-robin samples is presented, and the general results are summarized.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

QASMBench: A Low-Level Quantum Benchmark Suite for NISQ Evaluation and Simulation

The rapid development of quantum computing (QC) in the NISQ era urgently demands a low-level benchmark suite and insightful evaluation metrics for characterizing the properties of prototype NISQ devices, the efficiency of QC programming compilers, schedulers and assemblers, and the capability of quantum system simulators in a classical computer. In this work, we fill this gap by proposing a low-level, easy-to-use benchmark suite called QASMBench based on the OpenQASM assembly representation. It consolidates commonly used quantum routines and kernels from a variety of domains including chemistry, simulation, linear algebra, searching, optimization, arithmetic, machine learning, fault tolerance, cryptography, and so on, trading-off between generality and usability. To analyze these kernels in terms of NISQ device execution, in addition to circuit width and depth, we propose four circuit metrics including gate density, retention lifespan, measurement density, and entanglement variance, to extract more insights about the execution efficiency, the susceptibility to NISQ error, and the potential gain from machine-specific optimizations. Applications in QASMBench can be launched and verified on several NISQ platforms, including IBM-Q, Rigetti, IonQ and Quantinuum. For evaluation, we measure the execution fidelity of a subset of QASMBench applications on 12 IBM-Q machines through density matrix state tomography, comprising 25K circuit evaluations. In addition we also compare the fidelity of executions among the IBM-Q machines, the IonQ QPU and the Rigetti Aspen M-1 system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A modified approach for simulating electronically nonadiabatic dynamics via the generalized quantum master equation

We present a modified approach for simulating electronically nonadiabatic dynamics based on the Nakajima-Zwanzig generalized quantum master equation (GQME). The modified approach utilizes the fact that the Nakajima-Zwanzig formalism does not require casting the overall Hamiltonian in system-bath form, which is arguably neither natural nor convenient in the case of the Hamiltonian that governs nonadiabatic dynamics. Within the modified approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density operator is fully captured by a memory kernel super-operator. A methodology for calculating the memory kernel from projection-free inputs is developed. Simulating the electronic dynamics via the modified approach, with a memory kernel obtained using exact or approximate methods, can be more cost effective and/or lead to more accurate results than direct application of those methods. The modified approach is compared to previously proposed GQME-based approaches, and its robustness and accuracy are demonstrated on a benchmark spin-boson model with a memory kernel which is calculated within the Ehrenfest method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Hardware-Accelerated Ray Tracing of CAD-Based Geometry for Monte Carlo Radiation Transport

Monte Carlo radiation transport (MCRT) methods have been used to simulate radiation environments for many decades by tracking individual particles through a model to accumulate statistical information. MCRT geometry is historically formed using the constructive solid geometry (CSG). Recently, significant work has been performed to support simulations using computer-aided design (CAD)-based tessellated surfaces to support highly complex geometries. Ray tracing acceleration data structures from the rendering and visualization community are applied to accelerate particle tracking in CAD-based models. Despite these efforts, CSG representations provide the superior performance in surface intersection operations during particle flight. Concurrently, pseudo Monte Carlo methods have become prevalent in rendering applications to support more realistic models for scattering media, motivating innovations that are advantageous for MCRT simulations. Finally, the authors’ work extends these innovations by employing Intel’s Embree ray tracing kernel within a geometry toolkit for Monte Carlo to improve the simulation performance using CAD-based models by factors of 1.5 to 2.

42 ENGINEERING↗
DOE OSTI2023

Simulation of Mechanical Fractionation of Chopped Whole-Plant Corn (WPC) Using Discrete Element Method (DEM)

Fractionating whole-plant corn (WPC) in a single-pass harvesting system requires studies on the WPC-to-equipment interaction for improved property control, as well as mechanical and air-driven separation processes compared to the traditional multi-pass grain and stover harvesting system. The discrete element method (DEM) technique has the potential to simulate WPC mechanical fractionation and support simulation-based design of WPC separation processes. In this study, methods to develop DEM particle models of WPC (kernel, cob, stalk, and husk) and their material properties for simulating mass fractionation using the ASABE standard mechanical shaker were proposed. Measurement was done on the axial dimensions (major, intermediate, and minor) and mass of each WPC type (mean sample size is 56), sampled from single-pass harvesting. Applying gaussian multivariate regression and bootstrapping re-sampling techniques, a DEM particle approximate to each WPC was developed. Sensitivity analysis of the DEM Young‘s modulus, Poisson‘s ratio, and interaction parameters of coefficient of restitution, coefficient of rolling friction, and coefficient of static friction on mass fraction was performed after 156 ASABE sieve-shaking DEM simulation runs, generated using Latin Hypercube Design (LHD) design of experiment (DOE) from 19 DEM material parameters. DEM simulation using Hertz-Mindlin with flexible bond contact laws and DOE optimized material properties successfully reproduced the mass fractions retained in ASABE sieves at 9.8% mean relative error and a coefficient of determination of R2 = 0.87. Here, the DEM methodology developed for mechanical WPC mass fractionation could be deployed to perform virtual design of feedstock handling equipment and performance analysis of mechanical fraction systems.

09 BIOMASS FUELS↗
DOE OSTI2019

A comprehensive modeling procedure for estimating statistical properties of forced ignition

In this work, a comprehensive modeling procedure for estimating the probability of ignition with application to high altitude relights of aircraft combustors is developed. In these configurations, an ignitor is used to introduce high-enthalpy discharge into a fuel-laden but stratified flow. Because of the inherent variabilities inflow conditions, kernel discharge process, and the chaotic turbulent flow, ignition cannot be described deterministically, but only as a probabilistic measure. The proposed modeling framework consists of three components. The turbulent flow is represented using the large eddy simulation (LES) framework. The ignition process is modeled using a manifold approach, where the initial kernel evolution is represented using a homogeneous reactor while the latter part of the evolution is represented as a competition between diffusion and chemical reactions using a flamelet-type mapping. A combined lookup table that can track the evolution of the kernel through these two distinct reaction stages is developed. The table lookup variables are solved within the LES framework. The resulting simulation tool is then embedded within an uncertainty quantification approach, where variations in the turbulent flow, as well as operating and kernel properties, are simulated using a Monte-Carlo-based sampling approach. Techniques to reduce computational cost are used to obtain a robust, numerically accurate, and physically representative model for engine relight. The method is validated using experimental data for ignition of methane/air mixtures. Due to the comprehensive nature of the modeling procedure, it is found that the simulation tool reproduces experimentally observed ignition probabilities over a wide range of operating conditions.

42 ENGINEERING↗
DOE OSTI2024

DNS of ignition and flame stabilization in a simplified gas turbine premixer

With the increasing need for fuel flexibility, mitigation of auto-ignition (AI) inside gas turbine (GT) premixers becomes crucial. They must be designed to yield a sufficiently homogeneous fuel-air mixture to achieve low emissions while at the same time avoiding the occurrence of AI and subsequent flame stabilization. This challenge requires a detailed understanding of turbulent mixing and chemistry interactions. In the present work, a direct numerical simulation (DNS) of an array of jets in crossflow (JICF), representative of an industrial GT premixer, is reported to shed light on these complex phenomena. It is found that AI kernels form in the aft part of the premixer and coalesce into a flame front that then propagates upstream, mainly through the boundary layer, and successively engulfs the jets. This, therefore, suggests a significant role of the jet array pattern on the flame stabilization. It is noted that AI kernels continue to form independently during the whole time of the simulation. To clarify the contribution of AI and diffusion in the ignition kernels and the main flame, chemical explosive mode analysis (CEMA) is employed jointly with a kernel tracking algorithm. It is found that during the initial formation of the flame, many ignition kernels form in mixtures with low scalar dissipation rate and large contribution from AI mode. As they quickly grow, they merge into a single flame front that becomes increasingly more diffusion-assisted over time, balancing the AI mode. Turbulence is shown to have a significant enhancing effect in lean premixed flames, but further analysis is required to fully characterize it. These findings are relevant for the industrial premixer studied, and also for novel micromix concepts that may be used in the next generation of GT combustion systems.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

MONTE CARLO CROSS SECTION LOOKUP KERNEL FOR THE CEREBRAS WSE-2 IN CSL

This is a small kernel that was used to collect data for an upcoming paper. We would like to have the code be open source so that the reviewers (and then readers) of the paper can see the whole code, and can reproduce/verify our results. This is not a fully featured application, it cannot produce any useful simulation results, it just executes a small abstracted kernel using synthetic data. The purpose of the kernel is to understand the basic performance characteristics of an HPC kernel on novel AI accelerator architectures. The main kernel is written in the CSL coding language for use with the Cerebras WSE-2 AI accelerator. The kernel represented is the Monte Carlo cross section lookup kernel, which is a small kernel used by the Monte Carlo neutral particle transport algorithm. There is also a baseline kernel written in CUDA that we will include in the repository to form a basis for comparing the WSE-2 to GPU.

Tramm, John↗
DOE OSTI2021

Time-dependent orbital-free density functional theory: Background and Pauli kernel approximations

Time-dependent orbital-free density functional theory (DFT) is an efficient method for calculating the dynamic properties of large-scale quantum systems due to the low computational cost compared to standard time-dependent DFT. In this work, we formalize this method by mapping the real system of interacting fermions onto a fictitious system of noninteracting bosons. The dynamic Pauli potential and associated kernel emerge as key ingredients of time-dependent orbital-free DFT. Using the uniform electron gas as a model system, we derive an approximate frequency-dependent Pauli kernel. Pilot calculations suggest that space nonlocality is a key feature for this kernel. Nonlocal terms arise already in the second-order expansion with respect to unitless frequency and reciprocal space variable ($\frac{ω}{qk_F}$ and $\frac{q}{2k_F}$, respectively). Given the encouraging performance of the proposed kernel, we expect it will lead to more accurate orbital-free DFT simulations of nanoscale systems out of equilibrium. Additionally, the proposed path to formulate nonadiabatic Pauli kernels presents several avenues for further improvements which can be exploited in future work to improve the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗