Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 433 records · Page 24

Common workflows for computing material properties using different quantum engines

The prediction of material properties based on density-functional theory has become routinely common, thanks, in part, to the steady increase in the number and robustness of available simulation packages. This plurality of codes and methods is both a boon and a burden. While providing great opportunities for cross-verification, these packages adopt different methods, algorithms, and paradigms, making it challenging to choose, master, and efficiently use them. We demonstrate how developing common interfaces for workflows that automatically compute material properties greatly simplifies interoperability and cross-verification. We introduce design rules for reusable, code-agnostic, workflow interfaces to compute well-defined material properties, which we implement for eleven quantum engines and use to compute various material properties. Each implementation encodes carefully selected simulation parameters and workflow logic, making the implementer’s expertise of the quantum engine directly available to non-experts. All workflows are made available as open-source and full reproducibility of the workflows is guaranteed through the use of the AiiDA infrastructure.

36 MATERIALS SCIENCE↗

Simulation of the High Temperature Test Facility (HTTF) Core Using the 2D Ring Model with SAM

A two-dimensional ring model is developed with SAM to model the core of the High Temperature Test Facility (HTTF) at the system level. The ring model simplifies the complex structure of the HTTF core by converting the hexagonal rows of heaters and flow channels into layers of concentric annular rings. The ring model is first compared against a 3D-1D model where the solid structures are fully resolved in three dimensions while the fluid structures are modeled as one-dimensional flows. Comparison between the 3D-1D and the ring models shows that the latter can predict major parameters reasonably well under steady-state normal operating condition, but the heater temperatures are underpredicted. Adjustment is made to the effective thermal conductivity of the ceramic core of the ring model to improve the heater temperature predictions. Furthermore, the ring model is also used to simulate a transient Pressurized Conduction Cooldown (PCC) condition and is benchmarked with the experimental data from the HTTF Test PG-27. Good agreement is obtained between the experimental data and the predictions by the ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Improved modeling of in-ice particle showers for IceCube event reconstruction

The IceCube Neutrino Observatory relies on an array of photomultiplier tubes to detect Cherenkov light produced by charged particles in the South Pole ice. IceCube data analyses depend on an in-depth characterization of the glacial ice, and on novel approaches in event reconstruction that utilize fast approximations of photoelectron yields. Here, a more accurate model is derived for event reconstruction that better captures our current knowledge of ice optical properties. When evaluated on a Monte Carlo simulation set, the median angular resolution for in-ice particle showers improves by over a factor of three compared to a reconstruction based on a simplified model of the ice. The most substantial improvement is obtained when including effects of birefringence due to the polycrystalline structure of the ice. When evaluated on data classified as particle showers in the high-energy starting events sample, a significantly improved description of the events is observed.

47 OTHER INSTRUMENTATION↗

Shallow Cumulus Properties as Captured by Adiabatic Fraction in High-Resolution LES Simulations

Abstract Shallow convective clouds are important players in Earth’s energy budget and hydrological cycle, and are abundant in the tropical and subtropical belts. They greatly contribute to the uncertainty in climate predictions due to their unresolved, complex processes that include coupling between the dynamics and microphysics. Analysis of cloud structure can be simplified by considering cloud motions as a combination of moist adiabatic motions like adiabatic updrafts and turbulent motions leading to deviation from adiabaticity. In this work, we study the sizes and occurrence of adiabatic regions in shallow cumulus clouds during their growth and mature stages, and use the adiabatic fraction (AF) as a continuous metric to describe cloud processes and properties from the core to the edge. To do so, we simulate isolated trade wind cumulus clouds of different sizes using the System of Atmospheric Modeling (SAM) model in high resolution (10 m) with the Hebrew University spectral bin microphysics (SBM). The fine features in the clouds’ dynamics and microphysics, including small near-adiabatic volumes and a thin transition zone at the edge of the cloud (∼20–40 m in width), are captured. The AF is shown to be an efficient measure for analyzing cloud properties and key processes determining the droplet-size distribution formation and shape during the cloud evolution. Physical processes governing the properties of droplet size distributions at different cloud regions (e.g., core, edge) are analyzed in relation to AF. Significance Statement 1) This study investigates the evolution of cumulus clouds (Cu) using a 10-m-resolution LES model with spectral bin microphysics. 2) The study improves the understanding of the mutual effects of adiabatic updrafts and lateral entrainment and mixing. 3) The study demonstrates the existence of an adiabatic core in nonprecipitating Cu. 4) Shapes of the droplet size distributions are closely related to the adiabatic fraction values. 5) Utilization of high resolution reveals the existence of physically significant small features in the cloud structure, such as a narrow cloud interface zone and small adiabatic volumes.

54 ENVIRONMENTAL SCIENCES↗

Preliminary Implementation of Two-Dimensional Cartesian Solver in CTF-R

Sub-channel codes are one of the the modeling and simulation tools used for thermal-hydraulic analysis of nuclear reactors. A few examples of such sub-channel codes are the COolant Boiling in Rod Arrays (COBRA) family of codes. The approximations that are used to simplify the fluid conservation equations into sub-channel form, mainly that of axially-dominated flow, lead to noticeable limitations on sub-channels solvers for problems with significant flow in lateral directions. In this report, a two-dimensional Cartesian solver is developed and implemented within CTF-R, which is the residual solver in the North Carolina State University version of COBRA-TF (CTF). The new solver will enable CTF to simulate flow that is not axially-dominated. The appropriate Cartesian forms of the conservation equations are derived and implemented in the solver. Once the conservation equations are established, the process of constructing the matrix system was altered to solve a two-dimensional staggered grid system. A simple case was used to test that the two-dimensional Cartesian solver is accurate. The test problem does not include any source terms or flow in the lateral direction. The results show that the solver was able to run the simple case and converge to a steady-state solution. Future work will focus on testing existing capabilities by using test cases that include transients and equation cross-terms. Future work will also include adding additional capabilities such as enabling the solver to include cases with source terms and three dimensional cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fracture mechanics approach to TRISO fuel particle failure analysis

Weibull stress-based methods for failure probability assessment have been developed and analyzed to assess the integrity of tristructural isotropic (TRISO) fuel particles during fuel life cycles and accident operating conditions. While simple, these methods entail a number of drawbacks when stress concentrates near crack tips, including finite element mesh size dependency when the Weibull stress is averaged over the finite element domain. Fracture mechanics approaches involving the use of interaction integrals eliminate this lack of mesh convergence and produce consistent fracture predictions. In this work, we use an interaction integral approach to computing stress intensity factors for a crack in the inner pyrolytic carbon layer perpendicular to the silicon carbide layer, which is simplified representation of a failure mode in TRISO particles. Further, the interface between these two TRISO layers has been shown to become porous, which we simulate by considering a transition of mechanical properties over such porous length, i.e. the layers are modeled as a functionally graded material. Aspects such as porosity and thermal and irradiation eigenstrains are considered in computing the stress intensity factor from a fracture mechanics approach and compared with the known Weibull stress failure approach. The methodology introduced in this paper enables a more general fracture probability assessment in TRISO particles and eliminates the need to identify best suited parameters when using local or averaged stress-based failure criterion. Finally, the numerical sensitivity studies show how parameters such as the porous transition zone length, the material stiffness, and creep affect the probability of TRISO fuel particle failure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Structural Uncertainty in the Sensitivity of Urban Temperatures to Anthropogenic Heat Flux

One key source of uncertainty for weather and climate models is structural uncertainty arising from the fact that these models must simplify or approximate complex physical, chemical, and biological processes that occur in the real world. However, structural uncertainty is rarely examined in the context of simulated effects of anthropogenic heat flux in cities. Using the Weather Research and Forecasting (WRF) model coupled with a single-layer urban canopy model, it is found that the sensitivity of urban canopy air temperature to anthropogenic heat flux can differ by an order of magnitude depending on how anthropogenic heat flux is released to the urban environment. Moreover, varying model structures through changing the treatment of roof-air interaction and the parameterization of convective heat transfer between the canopy air and the atmosphere can affect the sensitivity of urban canopy air temperature by a factor of 4. Urban surface temperature and 2-m air temperature are less sensitive to the methods of anthropogenic heat flux release and the examined model structural variants than urban canopy air temperature, but their sensitivities to anthropogenic heat flux can still vary by as much as a factor of 4 for surface temperature and 2 for 2-m air temperature. Our study recommends using temperature sensitivity instead of temperature response to understand how various physical processes (and their representations in numerical models) modulate the simulated effects of anthropogenic heat flux.

54 ENVIRONMENTAL SCIENCES↗

Carbon nanotube porin diffusion in mixed composition supported lipid bilayers

Carbon nanotube porins (CNTPs), short pieces of carbon nanotubes capable of self-inserting into a lipid bilayer, represent a simplified model of biological membrane channels. We have used high-speed atomic force microscopy (HS-AFM) and all-atom molecular dynamics (MD) simulations to study the behavior of CNTPs in a mixed lipid membrane consisting of DOPC lipid with a variable percentage of DMPC lipid added to it. HS-AFM data reveal that the CNTPs undergo diffusive motion in the bilayer plane. Motion trajectories extracted from the HS-AFM movies indicate that CNTPs exhibit diffusion coefficient values broadly similar to values reported for membrane proteins in supported lipid bilayers. The data also indicate that increasing the percentage of DMPC leads to a marked slowing of CNTP diffusion. MD simulations reveal a CNTP-lipid assembly that diffuses in the membrane and show trends that are consistent with the experimental observations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

An extended Vlasov–Fokker–Planck approach for kinetic simulations of laser plasmas

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Furthermore, our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Photoelectron spectroscopy with entangled photons; enhanced spectrotemporal resolution

In this theoretical study, we show how photoelectron signals generated by time-energy entangled photon pairs can monitor ultrafast excited state dynamics of molecules with high joint spectral and temporal resolutions, not limited by the Fourier uncertainty of classical light. This technique scales linearly, rather than quadratically, with the pump intensity, allowing the study of fragile biological samples with low photon fluxes. Since the spectral resolution is achieved by electron detection and the temporal resolution by a variable phase delay, this technique does not require scanning the pump frequency and the entanglement times, which significantly simplifies the experimental setup, making it feasible with current instrumentation. Application is made to the photodissociation dynamics of pyrrole calculated by exact nonadiabatic wave packet simulations in a reduced two nuclear coordinate space. This study demonstrates the unique advantages of ultrafast quantum light spectroscopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Contribution of the Hydroxyl Group on Interfacial Heat Conduction of Monohydric Alcohols: A Molecular Dynamics Study

Monohydric alcohols have been used as promising phase change materials (PCMs) for low-temperature latent heat storage. However, the heat storage/retrieval rates are limited due to the low thermal conductivity of such alcohols. In this work, nonequilibrium molecular dynamics (NEMD) simulations were performed to study the microscopic heat conduction in example monohydric alcohols, i.e., 1-dodecanol (C12H26O), 1-tetradecanol (C14H30O), and 1-hexadecanol (C16H34O). A simplified ideal crystal model was proposed to exploit the potential for improving the thermal conductivity of monohydric alcohols. The effect of ideal crystalline structures, especially the contribution of the hydroxyl group, on the microscopic heat conduction process was analyzed. The thermal conductivity of the ideal crystals of the various monohydric alcohols was predicted to be more than twice as compared to that of their respective solids. The major thermal resistance in the ideal crystals was found around the molecular interfaces, as a result of the excellent heat conduction performance along the linear molecular chains. The calculated vibrational density of states (VDOS) and interfacial heat transfer were then investigated. When the interfaces are surrounded by hydroxyl groups as walls, strong hydrogen bond (HB) interactions were observed. The interfacial heat transfer coefficient of the ideal crystalline structures of 1-tetradecanol was found to reach up to ~735.6 MW/m2 W. It was elucidated that the high interfacial heat transfer rate is clearly related to the stronger intermolecular interactions.

Engineering↗

Conceptual and Mathematical Foundation for the FE/NETL CO 2 Prophet Model for Simulating CO 2 Enhanced Oil Recovery, Version 2

The United States (U.S.) Department of Energy (DOE) Office of Fossil Energy (FE) at the National Energy Technology Laboratory (NETL) has developed the FE/NETL CO 2 Prophet Model, which is Version 2 of CO 2 Prophet. This document describes the mathematical foundation for the FE/NETL CO 2 Prophet Model. CO 2 Prophet was originally developed in the 1990s by Texaco Exploration and Production and Technology Department for DOE. The FE/NETL CO 2 Prophet Model is an oil reservoir simulator that is suitable for simulating water floods and supercritical carbon dioxide (CO 2 ) enhanced oil recovery (EOR). The FE/NETL CO 2 Prophet Model uses a number of assumptions to simplify the equations describing the flow of oil, water (or brine), and CO 2 in the oil reservoir.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

User’s Manual for StrmtbFlow, the Stream Tube Multiphase Flow Part of the FE/NETL CO 2 Prophet Model, Version 2

The United States (U.S.) Department of Energy (DOE) Office of Fossil Energy (FE) at the National Energy Technology Laboratory (NETL) has developed the FE/NETL CO 2 Prophet Model, which is an updated version of the CO 2 Prophet. This document is a user’s manual for StrmtbFlow, which is one part of the FE/NETL CO 2 Prophet Model. CO 2 Prophet was originally developed in the 1990s by Texaco Exploration and Production Technology Department for DOE. The FE/NETL CO 2 Prophet Model is an oil reservoir simulator that is suitable for simulating water floods and supercritical carbon dioxide (CO 2 ) enhanced oil recovery (EOR). The FE/NETL CO 2 Prophet Model uses a number of assumptions to simplify the equations describing the flow of oil, water (or brine), and CO 2 in the oil reservoir.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

User’s Manual for StrmtbGen, the Stream Tube Generating Part of the FE/NETL CO 2 Prophet Model, Version 2

The United States (U.S.) Department of Energy (DOE) Office of Fossil Energy (FE) at the National Energy Technology Laboratory (NETL) has developed the FE/NETL CO 2 Prophet Model, which is Version 2 of CO 2 Prophet. This document is a User’s Manual for StrmtbGen, which is one part of the FE/NETL CO 2 Prophet Model. CO 2 Prophet was originally developed in the 1990s by Texaco Exploration and Production Technology Department for DOE. The FE/NETL CO 2 Prophet Model is an oil reservoir simulator that is suitable for simulating water floods and supercritical carbon dioxide (CO 2 ) enhanced oil recovery (EOR). The FE/NETL CO 2 Prophet Model uses a number of assumptions to simplify the equations describing the flow of oil, water (or brine), and CO 2 in the oil reservoir.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

PREFACE: A Search for Long-Lived Particles at the Large Hadron Collider

The Standard Model (SM) fails to explain many problems (neutrino masses, dark matter, and matter–antimatter asymmetry, among others) that may be resolved with new particles beyond the SM. No observation of such new particles may be explained either by their exceptionally high mass or by considerably small coupling to SM particles. The latter case implies relatively long lifetimes. Such long-lived particles (LLPs) then to have signatures different from those of SM particles. Searches in the “central region” are covered by the LHC general purpose experiments. The forward small angle region far from the interaction point (IP) is unexplored. Such particles are expected to have the energy as large as E = O(1 TeV) and Lorentz time dilation factor γ = E/m ≈ 10 2 –10 3 (with m the particle mass) hence long enough decay distances. A new class of specialized LHC detectors dedicated to LLP searches has been proposed for the forward regions. Among these experiments, FASER is already operational, and FACET is under consideration at a location 100 m from the LHC IP5 (the CMS detector intersection). However, some features of FACET require a specially enlarged beam pipe, which cannot be implemented for LHC Run 4. In this study, we explore a simplified version of the proposed detector PREFACE compatible with the standard LHC beam pipe in the HL-LHC Run 4. Realistic G EANT 4 simulations are performed and the background is evaluated. An initial analysis of the physics potential with the PREFACE geometry indicates that several significant channels could be accessible with sensitivities comparable to FACET and other LLP searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Setup and verification of a SCALE/KENO platform for generic FHR benchmark calculations

The work presented in this article is preliminary to downstream analysis of a generic fluoride salt-cooled high-temperature reactor (gFHR) core performed by the University of Tennessee in collaboration with Kairos Power (KP). A Monte Carlo transport model of the publicly available gFHR equilibrium core is developed in SCALE/KENO with multigroup energy treatment. Several output quantities of interest are used to verify the simulations against a benchmark model developed by KP using the continuous energy Monte Carlo code Serpent 2. Good agreement is seen in flux and fission rate profiles with a maximum relative difference of 1.4% and 2.8% respectively. Furthermore, an effective multiplication factor bias of 44 pcm was observed between the two simulations. The fuel temperature reactivity coefficient calculated with SCALE is within uncertainty to the reference model. This verification acts as a publicly reproducible benchmark for the gFHR in SCALE/KENO. A simplified depletion model is also presented where a single fuel pebble is depleted to discharge burnup through the equilibrium core while the equilibrium core is assumed to be invariant. This method produces results that intercept the equilibrium core concentrations in every case, however, an interesting artifact of this particular depletion model is uncovered. The phenomenon is shown to be a fundamental feature of the differential rate equations and inspires questions about how this system behaves when the time evolution of the equilibrium core is considered. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

2025 Advances in NekRS: Supporting improved performance for nuclear applications

This report presents several 2025 advancements in NekRS, a high-fidelity spectral element CFD code developed at Argonne National Laboratory to support the NEAMS thermal-hydraulics program. The forthcoming v25 release consolidates several of these advances, adding new features for portability across heterogeneous GPU architectures, real-time in situ visualization, improved turbulence modeling, and conjugate heat transfer coupling. Over the past year, NekRS has demonstrated strong scalability and performance on DOE’s leading exascale platforms, including Aurora and Frontier, confirming its readiness for some of the largest and most complex simulations attempted to date. These achievements provide a powerful new platform for high-fidelity data generation, which in turn supports the development and validation of advanced closure models critical for reactor safety and design. Significant algorithmic innovations have also been introduced. A new global runtime h-refinement capability simplifies workflows by reducing mesh preparation burdens and enabling coarse-to-fine restarts. Building on this, a novel multigrid strategy was implemented to accelerate pressure and transport solves at scale, addressing long-standing bottlenecks in exascale CFD. Together, these developments improve both the efficiency and accessibility of high-fidelity simulations for reactor-relevant problems. Collectively, these enhancements represent a major step forward in simulation technology, positioning NekRS as a cornerstone of NEAMS efforts to enable accurate, efficient, and scalable high-fidelity analysis of advanced nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗