SEARCH · Engineering Papers
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.
The agenda for the presentation is: introduction and overview; method description; Hydrodynamics problems; LED multi-physics problems; HED multi-physics problems; Kokkos port for ATS-1 and ATS-2; Summary and future work; and, Wrap-up.
ABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperature-dependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. ABINIT has strong links with many other software projects that are briefly mentioned.
This report highlights the development of advanced analytical capabilities at Oak Ridge National Laboratory to characterize PuO 2 powder to support the Multi-Lab Plutonium Process Signatures Campaign in FY 2022. Two samples from a batch of 76 statistical precipitation runs were packaged and shipped from Pacific Northwest National Laboratory to Oak Ridge National Laboratory at the beginning of FY 2022. Techniques were used to characterize the sample, including powder X-ray diffraction, Raman spectroscopy, and diffuse reflectance spectroscopy. A brief description of each technique, the results and discussion, and future work are described in this report. Opportunities to expand current capabilities and include additional analytical capabilities for the next campaign will also be discussed.
This full-scale pilot of a machine learning based nutrient controller/digital twin is being done as part of The Water Research Foundation (WRF) project 5121: Development of Innovative Predictive Control Strategies for Nutrient Removal. The overall goal of this project is to develop and full-scale test a hybrid (machine learning + mechanistic model) nutrient management controller at four different water resource recovery facilities (WRRFs). The controller demonstrates both short-term optimization functions and predictive capabilities. This project includes the design of the controller, its performance at the WRRFs, and a description of the future improvements.
The Materials Protection, Accounting, and Control Technologies (MPACT) campaign, within the U.S. Department of Energy Office of Nuclear Energy, has developed a Virtual Facility Distributed Test Bed for safeguards and security design for future nuclear fuel cycle facilities. The purpose of the Virtual Test Bed is to bring together experimental and modeling capabilities across the U.S. national laboratory and university complex to provide a one-stop-shop for advanced Safeguards and Security by Design (SSBD). Experimental testing alone of safeguards and security technologies would be cost prohibitive, but testbeds and laboratory processing facilities with safeguards measurement opportunities, coupled with modeling and simulation, provide the ability to generate modern, efficient safeguards and security systems for new facilities. This Virtual Test Bed concept has been demonstrated using a generic electrochemical reprocessing facility as an example, but the concept can be extended to other facilities. While much of the recent work in the MPACT program has focused on electrochemical safeguards and security technologies, the laboratory capabilities have been applied to other facilities in the past (including aqueous reprocessing, fuel fabrication, and molten salt reactors as examples). This paper provides an overview of the Virtual Test Bed concept, a description of the design process, and a baseline safeguards and security design for the example facility. Parallel papers in this issue go into more detail on the various technologies, experimental testing, modeling capabilities, and performance testing.
AS4A/SASSYS-1 (SAS) is a simulation tool used to perform deterministic analyses of anticipated events as well as design basis and beyond design basis accidents for advanced liquid-metal-cooled nuclear reactors. With its origin as SAS1A in the late 1960s, the SAS series of codes has been under continuous use and development for over sixty years and represents a critical investment in safety analysis capabilities for the U.S. Department of Energy. To support the dedication effort, this report has been generated to provide a detailed description of the available verification testing. The verification testing presented in this report captures functionality testing, focusing mainly on the testing of specific functions and algorithms for accuracy and precision of output, and interface testing, focusing mainly on the testing of critical input parameters and their valid ranges. Although SAS was developed to support the analysis of any liquid-metal-cooled nuclear reactor, the testing described in this document primarily focuses on the verification of SAS capabilities as they relate to a generic pool-type Sodium Fast Reactor (SFR).
Many important physical processes in inertial confinement fusion (ICF) and dense Z-pinch (DZP) experiments require a kinetic (velocity-space-dependent) description. Conventional particle-in-cell (PIC) methods are poorly suited for high-energy-density (HED) plasmas, due to restrictive time-step constraints and the inability to conserve energy. In a previous LDRD (21-FS-048), we demonstrated that a fully implicit PIC formulation overcomes these limitations: it conserves energy even when coupled with Coulomb collision models and can be solved efficiently with large grid cells and large time steps. Thus, it is feasible to use this method to study kinetic effects in ICF and DZP plasmas on hydro-like time and spatial scales. In this follow-on LDRD, we advanced this methodology into a high-fidelity tool for production-scale simulations and used it to answer key questions relevant to ICF and DZP experiment.
The adequate odorization of natural gas is critical to identify gas leaks and to reduce accidents. To ensure odorization, natural gas utility companies collect samples to be processed at core facilities or a trained human technician smells a diluted natural gas sample. In this work, we report a detection platform that addresses the lack of mobile solutions capable of providing quantitative analysis of mercaptans, a class of compounds used to odorize natural gas. Detailed description of the platform hardware and software components is provided. Designed to be portable, the platform hardware facilitates extraction of mercaptans from natural gas, separation of individual mercaptan species, and quantification of odorant concentration, with results reported at point-of-sampling. The software was developed to accommodate skilled users as well as minimally trained operators. Detection and quantification of six commonly used mercaptan compounds (ethyl mercaptan, dimethyl sulfide, n-propylmercaptan, isopropyl mercaptan, tert–butyl mercaptan, and tetrahydrothiophene) at typical odorizing concentrations of 0.1–5 ppm was performed using the device. We demonstrate the potential of this technology to ensure natural gas odorizing concentrations throughout distribution systems.
This mini-review describes developments over the last ~30 years in characterizing the nucleation & growth of epitaxially-supported metal nanoclusters (NCs) or islands during vapor deposition, as well as their post-deposition coarsening. A beyond-mean-field treatment for homogeneous nucleation & growth corrects the deficiencies of traditional treatments in describing, e.g., the island size distribution, but also necessitates consideration of the spatial distribution of islands and their capture zones. We discuss advances in modeling capabilities, including those based upon on an ab-initio level treatment of periphery diffusion kinetics, for description of the non-equilibrium growth shapes of these NCs, focusing on 2D NCs. For post-deposition coarsening of arrays of NCs, there is generally a competition between Ostwald Ripening (OR) and Smoluchowski Ripening (SR). SR is also known as Particle Migration & Coalescence. For 2D NCs in homoepitaxial systems, conventional OR is observed on pristine fcc(111) surfaces, dramatically enhanced OR in the presence of even trace amounts of chalcogens for Cu(111) and Ag(111), and anomalous OR on anisotropic fcc(110) surfaces. The unexpected discovery of SR for fcc(100) homoepitaxial systems prompted extensive analysis of the underlying diffusivities of 2D NCs as a function of size, as well as of NC coalescence dynamics. A comprehensive understanding of these processes is now available. Self-assembly of 3D NCs during deposition, issues related to heterogeneous nucleation, directed assembly, NC growth structure selection, and coarsening are addressed. Finally, for SR of 3D epitaxial NCs, recent insights into the size-dependence of diffusivity are described.
We analyze the potential of rare charm decays as probes of new physics at a high-luminosity flavor facility operating at the Z pole, such as the FCC-ee or CEPC. In particular, we identify clean null-test observables in $D^0 → π^+ π^- v\bar{v}$ and in polarized $Λ^+_c → p ℓ^+ ℓ^-$ decays with $ℓ = e, μ$. Complementarity with the LHC and HL-LHC flavor programs arises from the characteristic features of a Tera-Z environment: the capability to study missing-energy modes and charm production with significant polarization. We improve the theoretical description of $D^0 → π^+ π^- v\bar{v}$ decays and work out the phenomenology of polarization-induced null-test observables in $Λ^+_c → p ℓ^+ ℓ^-$ decays. In regions of dilepton mass near the φ resonance, polarization asymmetries can reach $O$(5%) for muons and $O$(14%) for electrons times the $Λ^+_c$ polarization. We also point out synergies between the dineutrino and the dilepton modes using the SMEFT framework of heavy new physics. Using the IDEA detector concept at FCC-ee, we find in simulation studies that dineutrino branching fractions as low as ~ 2 x 10 -7 can be probed, which reaches well into the parameter space of new physics, and also allows for discrimination of lepton flavor structures. Furthermore, the measurement of asymmetries in $Λ^+_c → p ℓ^+ ℓ^-$ at $O$(1%) will be possible. Similar sensitivities are expected for dielectron final states, although robust predictions will require further dedicated studies.
This paper introduces an innovative method for characterizing, implementing, and validating both three-phase and single-phase induction motor models, accompanied by a variable speed drive model. The primary goal is to investigate interactions between the electrical power grid and other dynamic domains (e.g., thermofluidic) that impact motor/load drive behavior. Our approach involves establishing a mechanical interface based on a physically meaningful equation linking motor torque/speed to the electrical model in the phasor domain. This allows seamless integration of diverse domain subsystems into a unified multi-domain model using Modelica v4.0.0 and the OpenIPSL library v3.0.1, overcoming co-simulation limitations. The proposed model, which requires only one Modelica-compliant tool for simulation, introduces additional dynamics through the mechanical interface, enabling explicit simulation of load disturbances based on constitutive physics. This deepens our understanding of dynamic interactions between the electrical power domain and other subsystems connected through the motor. We detail the modeled components using mathematical equations and textual descriptions, emphasizing the Modelica modeling approach. Simulation examples validate the implementation, demonstrating the multi-domain modeling capabilities of the newly developed components.
FECM/NETL Unconventional Shale Well Economic Model (UShWEM) is an Excel-based model that evaluates the economics of an unconventional shale well on a per-well and per-pad basis. This document serves as the user’s manual for the model with descriptions of the procedures the user must follow to run the model. This document also describes the capabilities of the model and provides the equations that are used by the model to calculate technical quantities and key model outputs including net cash flow, internal rate of return (IRR), net present value (NPV), earnings before interest, taxes, depreciation, and amortization (EBITDA), payout month and year, and breakeven price (for either oil- or gas-wells).
Understanding the emergence of hadron mass (EHM) represents one of the most challenging and still open problems in the Standard Model. Results of EHM studies from the nucleon resonance electroexcitation amplitudes (gvpN* electrocouplings) determined from exclusive meson electroproduction data measured with the CLAS detector will be presented in this talk. These CLAS data have provided the first and only available results on the evolution of the ?vpN* electrocouplings with photon virtuality Q2 up to 5 GeV2 for most nucleon resonances in the mass range up to 1.8 GeV. A successful description of the CLAS results on the Q2-evolution of the ?(1232)3/2+, N(1440)1/2+, and ?(1600)3/2+ electrocouplings achieved within the CSM approach by employing the same dressed quark mass function inferred from the QCD Lagrangian, which also provided a successful description of the pion elastic electromagnetic form factor and PDF, and nucleon elastic form factor, conclusively demonstrated the capability of gaining insight into EHM. The expected results from the CLAS12 detector on the Q2-evolution of the gvpN* electrocouplings for all prominent resonances in the mass range up to 2 GeV at the still almost unexplored range of Q2>5 GeV2 will allow us to place tight constraints on the dressed quark mass function in the range of quark momenta where the substantial part of hadron mass ~50% is expected to be generated. Consistent results on the momentum dependence of the dressed quark mass obtained from analyses of the ?vpN* electrocouplings of nucleon resonances of different structure will validate insight into EHM in a nearly model-independent way, allowing us to address key open problems in hadron physics on the emergence of hadron mass and nucleon resonance structure from QCD.
Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.
Thermodynamic modeling of the molten salt systems is critical in the design of advanced nuclear power plants, materials recycling, and many other important clean energy applications. As one of the most capable computational thermodynamic approaches, the CALPHAD (Calculation of Phase Diagrams) method has been widely used for multicomponent thermodynamic database development and further support more physics-based modeling. Although considerable CALPHAD modeling efforts for molten salt systems have been made, we found that no evaluation has been performed regarding the unary thermodynamic models for the molten salt systems despite the fact that more than one pure substance database is available for the CALPHAD modeling. Therefore, the thermodynamic descriptions of the pure salts should be critically evaluated to ensure a high-fidelity molten salt thermodynamic database. In this work, we comprehensively analyzed some selected molten salts using different available thermodynamic descriptions from two databases as an example. One is the SSUB database released by the SGTE (Scientific Group Thermodata Europe), and the other is the FactPS database from the FactSage. Finally, the significant discrepancies observed in the thermodynamic description between the existing CALPHAD pure substance databases and experiments call for an urgent effort to improve.
Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.
This report documents the efforts to model the conceptual design of the STS natural circulation flow loop using TRACE. The results, summarized in Section 2.1, support the conclusion that the conceptual STS natural circulation flow loop design provides sufficient cooling of the tungsten target and exhibits predictable flow behavior, provided that certain design features and/or operating conditions are used. In particular: A sufficient initial liquid volume fraction is required to avoid dryout of the target flow channels; Inclusion of a pressurizer or bellows volume at the top of the loop is recommended and is capable of achieving steady, non-oscillatory flow and temperature conditions at sufficiently high loop pressure. Section 1 provides a detailed description of the geometry and assumptions used in the model, and Section 2 presents the results for the current flow loop design as well as some investigation of potential design modifications with the goal of enhancing the thermal hydraulic stability of the loop.
Abstract Coarse-grained descriptions of dislocation motion in crystalline metals inherently represent a loss of information regarding dislocation-dislocation interactions. In the present work, we consider a coarse-graining framework capable of re-capturing these interactions by means of the dislocation-dislocation correlation functions. The framework depends on a convolution length to define slip-system-specific dislocation densities. Following a statistical definition of this coarse-graining process, we define a spatial correlation function which will allow the arrangement of the discrete line system at two points—and thus the strength of their interactions at short range—to be recaptured into a mean field description of dislocation dynamics. Through a statistical homogeneity argument, we present a method of evaluating this correlation function from discrete dislocation dynamics simulations. Finally, results of this evaluation are shown in the form of the correlation of dislocation densities on the same slip-system. These correlation functions are seen to depend weakly on plastic strain, and in turn, the dislocation density, but are seen to depend strongly on the convolution length. Implications of these correlation functions in regard to continuum dislocation dynamics as well as future directions of investigation are also discussed.