Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 505 records · Page 28

Validation of Thermal Scattering Laws for Light Water at Elevated Temperatures with Diffusion Experiments

Light water is the most important neutron moderator in many reactor applications. Thermal neutron scattering kernels for H bound in H 2 O are represented by ENDF-format thermal scattering laws (TSLs) evaluated at discrete temperatures T . Nuclear data evaluations are conventionally validated utilizing a combination of critical benchmarks and experimental cross section data. Existing public critical benchmarks may be inadequate for testing water TSLs over the full range of T of interest in reactor applications, and experimental thermal scattering cross section data for water is sparse at elevated T . In this work, MC21 is used to simulate the decay of the equilibrium thermal neutron flux in light-water spheres of several radii. The fundamental-mode time eigenvalue $\mathcal{a}$ is calculated at 22 °C and 227 °C (at saturation pressure) for each sphere using three different H-H 2 O TSLs. Fitting $\mathcal{a}$ to a polynomial function of geometric buckling allows calculation of the thermal neutron diffusion length L . Extrapolation lengths are treated as a function of the transport mean free path and geometry. Experimental L data from 25 publications at 49 temperatures (from 10 °C to 295 °C), computed by several different time-dependent and space-dependent decay methods, is used to develop an empirical fit of L vs. T . The MC21-calculated L for the TSLs tested are within ±1% of the predicted values at 22 °C and 227 °C. This validation approach, which may be repeated at arbitrary T , constitutes an integral benchmark specific to and characterizing the detailed physics of the thermal scattering kernel applied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A fast platform for simulating semi-flexible fiber suspensions applied to cell mechanics

We present a novel platform for the large-scale simulation of three-dimensional fibrous structures immersed in a Stokesian fluid and evolving under confinement or in free-space in three dimensions. One of the main motivations for this work is to study the dynamics of fiber assemblies within biological cells. For this, we also incorporate the key biophysical elements that determine the dynamics of these assemblies, which include the polymerization and depolymerization kinetics of fibers, their interactions with molecular motors and other objects, their flexibility, and hydrodynamic coupling. This work, to our knowledge, is the first technique to include many-body hydrodynamic interactions (HIs), and the resulting fluid flows, in cellular assemblies of flexible fibers. We use non-local slender body theory to compute the fluid–structure interactions of the fibers and a second-kind boundary integral formulation for other rigid bodies and the confining boundary. A kernel-independent implementation of the fast multipole method is utilized for efficient evaluation of HIs. The deformation of the fibers is described by nonlinear Euler–Bernoulli beam theory and their polymerization is modeled by the reparametrization of the dynamic equations in the appropriate non-Lagrangian frame. We use a pseudo-spectral representation of fiber positions and implicit time-stepping to resolve large fiber deformations, and to allow time-steps not excessively constrained by temporal stiffness or fiber–fiber interactions. The entire computational scheme is parallelized, which enables simulating assemblies of thousands of fibers. We use our method to investigate two important questions in the mechanics of cell division: (i) the effect of confinement on the hydrodynamic mobility of microtubule asters; and (ii) the dynamics of the positioning of mitotic spindle in complex cell geometries. Finally to demonstrate the general applicability of the method, we simulate the sedimentation of a cloud of semi-flexible fibers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Taking Lessons Learned from a Proxy Application to a Full Application for SNAP and PARTISN

SNAP is a proxy application which simulates the computational motion of a neutral particle transport code, PARTISN. Here in this work, we have adapted parts of SNAP separately; we have re-implemented the iterative shell of SNAP in the task-model runtime Legion, showing an improvement to the original schedule, and we have created multiple Kokkos implementations of the computational kernel of SNAP, displaying similar performance to the native Fortran. We then translate our Kokkos experiments in SNAP to PARTISN, necessitating engineering development, regression testing, and further thought.

97 MATHEMATICS AND COMPUTING↗

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Balar: A SST GPU Component for Performance Modeling and Profiling

Programmable accelerators have become commonplace in modern computing systems. Advances in programming models and the availability of massive amounts of data have created a space for massively parallel accelerators capable of maintaining context for thousands of concurrent threads resident on-chip. These threads are grouped and interleaved on a cycle-by-cycle basis among several massively parallel computing cores. One path for the design of future supercomputers relies on an ability to model the performance of these massively parallel cores at scale. The SST framework has been proven to scale up to run simulations containing tens of thousands of nodes. A previous report described the initial integration of the open-source, execution-driven GPU simulator, GPGPU-Sim, into the SST framework. This report discusses the results of the integration and how to use the new GPU component in SST. It also provides examples of what it can be used to analyze and a correlation study showing how closely the execution matches that of a Nvidia V100 GPU when running kernels and mini-apps.

97 MATHEMATICS AND COMPUTING↗

SST-GPU: A Scalable SST GPU Component for Performance Modeling and Profiling

Programmable accelerators have become commonplace in modern computing systems. Advances in programming models and the availability of unprecedented amounts of data have created a space for massively parallel accelerators capable of maintaining context for thousands of concurrent threads resident on-chip. These threads are grouped and interleaved on a cycle-by-cycle basis among several massively parallel computing cores. One path for the design of future supercomputers relies on an ability to model the performance of these massively parallel cores at scale. The SST framework has been proven to scale up to run simulations containing tens of thousands of nodes. A previous report described the initial integration of the open-source, execution-driven GPU simulator, GPGPU-Sim, into the SST framework. This report discusses the results of the integration and how to use the new GPU component in SST. It also provides examples of what it can be used to analyze and a correlation study showing how closely the execution matches that of a Nvidia V100 GPU when running kernels and mini-apps.

97 MATHEMATICS AND COMPUTING↗

Computational Benefit of GPU Optimization for the Atmospheric Chemistry Modeling

Abstract Global chemistry‐climate models are computationally burdened as the chemical mechanisms become more complex and realistic. Optimization for graphics processing units (GPU) may make longer global simulation with regional detail possible, but limited study has been done to explore the potential benefit for the atmospheric chemistry modeling. Hence, in this study, the second‐order Rosenbrock solver of the chemistry module of CAM4‐Chem is ported to the GPU to gauge potential speed‐up. We find that on the CPU, the fastest performance is achieved using the Intel compiler with a block interleaved memory layout. Different combinations of compiler and memory layout lead to ~11.02× difference in the computational time. In contrast, the GPU version performs the best when using a combination of fully interleaved memory layout with block size equal to the warp size, CUDA streams for independent kernels, and constant memory. Moreover, the most efficient data transfer between CPU and GPU is gained by allocating the memory contiguously during the data initialization on the GPU. Compared to one CPU core, the speed‐up of using one GPU alone reaches a factor of ~11.7× for the computation alone and ~3.82× when the data transfer between CPU and GPU is considered. Using one GPU alone is also generally faster than the multithreaded implementation for 16 CPU cores in a compute node and the single‐source solution (OpenACC). The best performance is achieved by the implementation of the hybrid CPU/GPU version, but rescheduling the workload among the CPU cores is required before the practical CAM4‐Chem simulation.

54 ENVIRONMENTAL SCIENCES↗

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Prediction of thermal conductivity in dielectrics using fast, spectrally-resolved phonon transport simulations

We present a new method for predicting effective thermal conductivity in materials, informed by ab initio material property simulations. Using the Boltzmann transport equation in a self-adjoint angular flux formulation, we performed simulations in silicon at room temperatures over length scales varying from 10 nm to 10 micron and report temperature distributions, spectral heat flux and thermal conductivity. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. In this method, a closure term is introduced to the transport equation which acts as a redistribution kernel for the total energy bath of the system. This term is an effective indicator of the degree of disorder between the spectral phonon radiance and the angular phonon intensity of the transport system. In conclusion, we employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations.

36 MATERIALS SCIENCE↗

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗

A Cold Moderator For Sub-Thermal Neutron Flux Enhancement at the RPI-LINAC

The importance of high-quality thermal neutron scattering kernels coupled with a lack of accurate sub-thermal total cross section data prompted researchers at Rensselaer Polytechnic Institute (RPI) to design a cold moderation system to enhance sub-thermal neutron flux (defined as flux below 5 meV), for nuclear data measurements at the RPI Gaerttner LINAC facility. Using a series of neutronic and heat transfer simulations obtained with MCNP and COMSOL, an optimized final design was produced. The final design exceeded current target sub-thermal neutron producing capabilities by a factor of 8 in the sub-thermal region and is easily configurable with existing neutron producing targets.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Genetic Algorithm for Hyperparameter Optimization in Gaussian Process Modeling

A genetic algorithm is developed and applied to optimize hyperparameters of convolutional recursively determined dual neural network-Gaussian process (NNGP) kernels. As a specific application of the combined GPNN-GA algorithm, it is applied to image classification in publicly available data of Hyper Suprime-Cam Subaru Strategic Program. Matthews correlation coefficient is calculated based on results of binary star-galaxy classification and used as a fitting function of the GA module of the algorithm. The simulation results confirm significant improvement of the classification accuracy with optimized hyperparameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗

Computational diagnostics and characterization of combustion recession in diesel sprays

While low-temperature combustion (LTC) strategies have been found to mitigate nitrogen oxides and particulate matter emissions in diesel engines, studies have also associated LTC with an increase in unburnt hydrocarbons. With more recent studies on diesel after end-of-injection (AEOI), combustion recession is identified as a phenomenon where at near nozzle region, high-temperature ignition (HTI) combustion can propagate back to the nozzle tip consuming the unburnt hydrocarbons AEOI. Current literature has suggested that combustion recession is controlled by auto-ignition. However, high-fidelity simulations and detailed analysis of such a mechanism are missing. In this study, comprehensive Large Eddy Simulations of a reacting spray at “Spray A” conditions are performed, where detailed analysis of combustion recession concerning flame morphology and propagation modes are included. In particular, this study demonstrated for the first time that while combustion recession is mainly auto-ignition dominated (consistent with the literature), a cool flame was found to deflagrate towards the richer regions of the mixture, promoting mixing and increasing the mixture temperature. This leads to HTI kernels, which then grow and develop as deflagrative waves, therefore sustaining the combustion recession process. The study also detailed the extinction mechanism of combustion: the entrainment wave will overlean the near-nozzle mixtures, rendering it unable to support HTI, which leads to the extinction of the upstream flame AEOI in lower reactivity mixtures. Combustion recession is also observed to be contingent on the chemical and diffusion processes, even at low scalar dissipation rates. Finally, a new criterion for combustion recession based on chemical explosive mode is proposed and validated with previous combustion recession index to quantify the extent of HTI in near-nozzle mixtures AEOI. The newly developed metric combined with a previous experimentally-based metric can provide simple but valuable measurements of the degree and propensity of the upstream flame AEOI.

CEMA↗

Assessing Radiative Feedbacks and Their Contribution to the Arctic Amplification Measured by Various Metrics

Arctic amplification (AA), characterized by a more rapid surface air temperature (SAT) warming in the Arctic than the global average, is a major feature of global climate warming. Various metrics have been used to quantify AA based on SAT anomalies, trends, or variability, and they can yield quite different conclusions regarding the magnitude and temporal patterns of AA. This study examines and compares various AA metrics for their temporal consistency in the region north of 70°N from the early twentieth to the early 21st century using observational data and reanalysis products. We also quantify contributions of different radiative feedback mechanisms to AA based on short-term climate variability in reanalysis and model data using the Kernel-Gregory approach. Albedo and lapse rate feedbacks are positive and comparable, with albedo feedback being the leading contributor for all AA metrics. The net cloud feedback, which has large uncertainties, depends strongly on the data sets and AA metrics used. By quantifying the influence of internal variability on AA and related feedbacks based on global climate model ensemble simulations, we find that water vapor and cloud feedbacks are most heavily affected by internal variability.

54 ENVIRONMENTAL SCIENCES↗

Predicting Cycle-to-cycle Variations in a Spark-ignition Engine using Multi-cycle Large Eddy Simulation

Multi-cycle large eddy simulation (LES) of a spark-ignition engine has been performed to predict cycle-to-cycle variations (CCV). To describe spark-ignited premixed flame propagation, the front propagation formulation (FPF) method is combined with a laminar-to-turbulent flame transition model that describes the non-equilibrium flame development processes during the initial stage of the kernel growth. In the laminar-to-turbulent flame transition model, the effects of turbulence on the variations in the laminar-to-turbulent flame transition time are modeled based on the Kolmogorov energy cascade. 15 consecutive cycle LES of in-cylinder processes in a spark-ignition four-valve single cylinder engine is conducted. The predicted in-cylinder pressure is compared with the experimental measurements. The range of CCV observed in the experiment is reproduced to an acceptable degree. The predicted CCV is found to be closely related to the difference in the turbu- lence field near the spark plug around the time of spark discharging, which is consistent with the previous LES studies. The current results also suggest the importance of capturing the variations of the laminar-to-turbulent flame transition in predicting CCV.

Su, Yunde↗

Predicting Cycle-to-cycle Variations in a Spark-ignition Engine using Multi-cycle Large Eddy Simulation

Multi-cycle large eddy simulation (LES) of a spark-ignition engine has been performed to predict cycle-to-cycle variations (CCV). To describe spark-ignited premixed flame propagation, the front propagation formulation (FPF) method is combined with a laminar-to-turbulent flame transition model that describes the non-equilibrium flame development processes during the initial stage of the kernel growth. In the laminar-to-turbulent flame transition model, the effects of turbulence on the variations in the laminar-to-turbulent flame transition time are modeled based on the Kolmogorov energy cascade. 15 consecutive cycle LES of in-cylinder processes in a spark-ignition four-valve single cylinder engine is conducted. The predicted in-cylinder pressure is compared with the experimental measurements. The range of CCV observed in the experiment is reproduced to an acceptable degree. The predicted CCV is found to be closely related to the difference in the turbulence field near the spark plug around the time of spark discharging, which is consistent with the previous LES studies. The current results also suggest the importance of capturing the variations of the laminar-to-turbulent flame transition in predicting CCV.

Yu, Yunde↗

Demonstrate FARM supervisory capabilities for a thermal energy storage problem for the DETAIL facility: IES Simulation Ecosystem Control System Development

The goal of the power dispatch problem for an Integrated Energy System (IES) is to adjust the power output and the heat flow of each component to maximize the profitability of the whole unit. Facilities that can integrate real-time digital signals, mock nuclear power, thermal energy storage and industrial heat use via high-temperature electrolysis were constructed at INL to support the research activities. The Dynamic Energy Technology and Integration Laboratory (DETAIL) houses the Microreactor Agile Non-nuclear Experimental Test Bed (MAGNET) and the Thermal Energy Distribution System (TEDS). In this report, the hierarchical control system architecture proposed in June 2023 milestone for the flexible operation of DETAIL facility is finalized and demonstrated. A brief description of the components and the corresponding Dymola models from the HYRBID repository is first provided. Then, the current control strategy is presented. In particular, the approach for generating the set-point trajectories to be fed to the PI controllers is analyzed, and its limits were identified. To preserve safe operation over both long-time and real-time horizons, the integration of a Supervisory Control layer embedding a modified version of FARM (Feasible Actuator Range Modifier) module is proposed. FARM is a component of the RAVEN-based FORCE framework designed to support HERON module at optimizing the operation of IES units. The proposed control system for DETAIL foresees FARM to be applied twice, i.e., the original version (“FARM-Validator”) aiding the solution of the power dispatch problem, and a modified version (“FARM-Supervisory”) coordinating the PID controllers. Despite the kernel of the two modules is the same, their tasks are quite different. The former intervenes at the beginning of each hour to prevent constraint violations over long time periods, the latter addresses real-time control tasks and monitors the response of constrained variables at a much finer time resolution. A tentative procedure for training the embedded Digital Twins with the experimental data is also proposed. Finally, the capabilities of the designed architecture and the impact of the added Supervisory Control layer are demonstrated by simulating a representative power dispatch scenario.

25 ENERGY STORAGE↗