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.
We discuss considerations that can be used to formulate recommendations for initiating a lepton collider project that would provide precision studies of the Higgs boson and related electroweak phenomena.
This project adapts and applies the gyrokinetic Eulerian code COGENT to an emerging area of highly-magnetized Z-pinch plasma systems. Simulation efforts for these plasmas currently use either magneto-hydrodynamic (MHD) fluid codes, which suffer from fidelity being restricted to narrow plasma conditions, or fully kinetic particle-in-cell (PIC) methods, which suffer from the run-time needed to resolve the very small spatial and temporal scales required for the plasma conditions. The progress with achieving detailed theoretical understanding of Z-pinch plasma systems can be greatly expedited using some combination of gyrokinetic and drift-ordered MHD models, developed in the tokamak community and implemented in the COGENT code. These formulations are applicable to a much larger parameter space than standard MHD models, while providing a substantial speed-up in run-time over fully kinetic codes. In particular, the gyrokinetic formalism is successfully applied to a flow-stabilized Z pinch (FSZP) system – a promising fusion concept where a sheared axial flow is used to provide stability against large-scale ideal-MHD modes. Results of gyrokinetic simulations with COGENT agree well with those from fully kinetic simulations, while requiring ~500x less run-time.
This project is a collaborative effort between InnoSense LLC (ISL) and Oak Ridge National Laboratory (ORNL) aimed at mitigating the mercury (Hg) contamination problem at the Y-12 National Security Complex (NSC). The report summarizes the research and technology development work of the Phase IIB Project, led by Dr. Uma Sampathkumaran at ISL for the period of July 31, 2017–October 31, 2021. The Y-12 site cleanup includes technologies for safe and cost-effective facility demolition, soil remediation, and on-site treatment and disposal of mercury bearing wastes. One of the technology demonstration goals was a reactive strippable coating for targeted pre-demolition and de-contamination to limit Hg mobilization and the volume of debris requiring treatment for disposal. Toward this, ISL has developed a water-based reactive strippable coating (Trap & See-Hg™) to sequester elemental Hg(0) and mercuric Hg(II) species in the coating. The coating provides multiple benefits including (1) Hg vapor suppression, (2) sealing in the hazardous contaminants, (3) visually locating the contaminant through color changes, and (4) remediating the contaminated surfaces upon coating removal. Through collaboration between ORNL and ISL in Phase II work, the team successfully demonstrated the Trap & See-Hg technology to (1) visually locate the presence of Hg on concrete, soil and drywall surfaces, (2) effectively sequester all mercury species (elemental, soluble and insoluble Hg) from lab-generated and ex-situ soil and brick/rubble samples collected from the Y 12 site with 80 to 95% efficiency to capture/reduce Hg vapor concentration, (3) be applied by manual and robotic spray applicators and cure at room ambient, (4) remain a cohesive and peelable coating upon drying, and (5) sequester Hg0 from walls in the alpha-2 building basement at Y-12 during a small scale in-situ field trial. In Phase IIB, ISL’s effort was focused on refining the coating formulations with additives containing passive and reactive sorbents, colorimetric indicator, and thixotropic agents to achieve the desired performance. The paints were then evaluated at ORNL for applicability to metal and concrete surfaces of different sizes and shapes, and soil rubble, Hg-sorption capacity and achieving sag-free coating on surfaces when applied by a manual spray applicator. Ambient cure of the coating resulted in a dry thickness of 16–20 mil. The coatings were removed from smooth or rough surfaces with low to moderate force. Both lab-generated Hg-contaminated samples and soil and brick/rubble samples containing a range of Hg concentrations (collected from the Y-12 site by UCOR) were used in this project. Optimized formulations were applied to these samples and evaluated for ability to suppress Hg vapors, sequester Hg 0 and Hg(II) species and color mapping of mercury species on the test surface. The ORNL team initially identified three to four potential sites at Y-12 and ORNL for potential demonstration. However, the pandemic closed any possible avenues for on-site demonstration of the robotic spray applicator. The metal and concrete samples with low to high Hg content provided by UCOR from an ongoing demolition were subsequently used for ex-situ sample characterization and analysis at ORNL over a period of 25 days. The coatings suppressed Hg vapors with an efficiency of 95–99% from heavily contaminated metal and 60% from concrete substrate. Analysis of total Hg in the peeled paints indicates sorption capacities ranging from 12 µg Hg/g of paint from concrete to 1.2 mg Hg/g of paint from the metal substrates while peeled paints from soil samples could capture a total Hg from 2.3 mg Hg/g to 66 mg/g paint. Toxicity characteristic leaching procedure (TCLP) tests recovered total mercury leached at 0.15±0.0008 µg/L for concrete to 126.81±0.67 µg/L and 377.25±5.67 µg/L for the two metal substrates, to allow for informed decisions on downstream disposal considerations. Together, this Phase IIB demonstrated the potential value of the Trap & See-Hg for Hg uptake, vapor suppression, and color mapping, and may thus be used to enhance worker safety during site cleanup, facility demolition, and downstream waste disposal at Y-12 NSC.
The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).
Demonstrate how to use project PARETO framework for screening different produced water networks and identify critical minerals recovery opportunities. More specifically, demonstrate how different mathematical formulations can have an impact on the solution time.
The residential sector consumes a significant portion of the electricity sold in the United States. Above 60% of the energy used in the sector is used to operate heating, ventilation, and air conditioning (HVAC) systems and water heating (WH) systems. With the increase of intelligence in the grid and the new decision and control options enabled by the Internet of Things; control of these devices can be used to support the grid. Therefore, this article presents a scalable multiagent system for optimizing HVAC and WH systems while maintaining comfort. It allows a utility to orchestrate the shifting of energy from critical periods without direct control, but instead by using a price signal. The architecture, optimization formulation, implementation strategy and results from an implementation project are discussed.
Recently, we introduced the “Newman-Penrose map”, a novel correspondence between a certain class of solutions of Einstein’s equations and self-dual solutions of the vacuum Maxwell equations, which we showed was closely related to the classical double copy. Here, we give an alternative definition of this correspondence in terms of quantities that are defined naturally on twistor space, and a shear-free null geodesic congruence on Minkowski space whose twistorial character is articulated by the Kerr theorem. The advantage of this reformulation is that it is purely geometrical in nature, being manifestly invariant under both spacetime diffeomorphisms and projective transformations on twistor space. While the original formulation of the map may be more convenient for most explicit calculations, the twistorial formulation we present here may be of greater theoretical utility.
Elastomeric rubbers serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose-liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. Cyclic (de)pressurization is known to degrade these materials through the process of cavitation. This readily visible failure mode occurs as a fracture or rupture of the material and is due to the oversaturated gas localizing to form gas bubbles. Computational modeling in the Hydrogen Materials Compatibility Program (H-Mat), co-led by Sandia National Laboratories and Pacific Northwest National Laboratory, employs multi-scale simulation efforts to build a predictive understanding of hydrogen-induced damage in materials. Modeling efforts within the project aim to provide insight into how to formulate materials that are less sensitive to high-pressure hydrogen-induced failure. In this document, we summarize results from atomistic molecular dynamics simulations, which make predictive assessments of the effects of compositional variations in the commonly used elastomer, ethylene propylene diene monomer (EPDM).
In this project, we sought to develop new cloud and rain microphysics frameworks within the Energy Exascale Earth System Model (E3SM). This work encompassed two primary avenues of research: 1) Further development of a Bayesian-based scheme called BOSS (Bayesian Observationally-constrained Statistical-physical Scheme) to represent cloud and rain microphysics, testing it in realistic high-resolution cloud models, and implementing it in E3SM; 2) Development of a methodology utilizing machine learning to enable computationally tractable use of tractable use of Markov chain Monte Carlo sampling for Bayesian parameter estimation in Earth system and cloud models. In this project, we adapted the BOSS microphysics scheme, originally formulated for rain-only, to include all liquid-phase microphysical processes for cloud and rain, in particular the processes that mediate between these two categories, for example the conversion from cloud to rain through collision and coalescence of drops. We constrained the scheme via comparison and testing against a detailed model that explicitly represents the evolution of cloud and rain particles, called a bin microphysics scheme.
A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.
We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.
This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.
The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.
Chris Saunders and three technologists are in high demand from Sandia’s deep learning teams, and they’re kept busy by building new clusters of computer nodes for researchers who need the power of supercomputing on a smaller scale. Sandia researchers working on Laboratory Directed Research & Development (LDRD) projects, or innovative ideas for solutions on short timeframes, formulate new ideas on old themes and frequently rely on smaller cluster machines to help solve problems before introducing their code to larger HPC resources. These research teams need an agile hardware and software environment where nascent ideas can be tested and cultivated on a smaller scale.
Solid surfaces act as catalysts for many important chemical processes, including commodity chemicals production, energy conversion and pollution mitigation. The outcome of catalytic processes is mainly governed by interactions of adsorbates with active sites on the catalyst surface. These interactions can be dramatically affected by the changes in the local characteristics of the active site, which can be accomplished by an introduction of various promoters and poisons, or by alloying with other metal atoms. In fact, many commercial heterogeneous catalysts contain metal nanoparticles promoted by various additives. Our goal is to develop fundamental understanding about the impact of perturbations in the local geometry of a catalytic site on the local surface reactivity. These insights are crucial for a more systematic and rational catalyst design. The project aims to advance the field of catalysis by formulating reliable structure/performance relationships and utilizing these relationships to identify novel catalytic materials. Benefits for DOE and our society are in the development of fundamental theories that guide the process of discovery of energy-efficient and environmentally friendly heterogeneous catalysis.
This ARPA-E project developed a machine learning tool to use in formulation design of cementitious binders for concrete having 50% less embodied CO 2 and possessing twice the durability compared to concrete based on ordinary portland cement (OPC) binders. The technical focus was on limestone/calcined clay cement (LC3), the leading replacement for OPC. Here, hierarchical machine learning (HML) was used to model the flowability, set time, strength, and durability of LC3 concrete. This methodology identifies latent variables derived from domain knowledge and empirical models that develop an accurate model for a response surface from small datasets. For the flowability metric, particle packing was a dominant factor, while strength and durability were both strongly determined by the fraction of metakaolin and the water:solids ratio. Under constraints of water:binder ratio, material performance metrics, embodied CO 2 , and cost per tonne of OPC, multi-objective optimization was used to design binders parameterized by the mineral composition replacing OPC, particle size distributions, and water:solids ratio. The trained algorithm was able to predict multiple mixes met these performance criteria, and experimental testing validated the predictions. The model demonstrated here is relevant for North America, where pure kaolin deposits are found broadly. The approach is being taken forward into commercial application by Ansatz AI, a materials informatics company founded by PI Washburn and co-PI Poczos. Through collaborations with the cement and concrete industry, and funding from SBIR programs, a commercial software will be developed in future research.
We report the spin decomposition of a spin-contaminated different-orbitals-for-different-spins (DODS) wave function is formulated in terms of the Krylov space of the $\hat{S}^2$ ||operator. The number of determinants that contribute to the various projected spin eigenfunctions vertical | $\psi ; S, M \rangle$ depends on the sparsity of the orbital overlap matrix X between the α and β spatial orbitals. The cosine-sine decomposition (CSD) procedure may be applied to this overlap matrix, and the resulting redundant orbital transformations may be applied to the α and β spatial orbitals. This produces a sparse X' matrix in the transformed basis in which each row or column has either one or two nonzero elements. This sparse X' matrix simplifies the spin-decomposition procedure in three ways:1) it reduces the number of contributing determinants within the Krylov basis functions and within the projected spin functions, 2) it reduces the effective orbital dimension through elimination of the frozen core and frozen virtual orbitals, and 3) it simplifies the spin-decomposition procedure in both the original and transformed bases by limiting the Krylov space dimension. These simplifications all reduce the computational effort for the spin-decomposition process. This procedure is implemented within a string-based DODS determinant formulation and applied to spin projection of unrestricted Hartree-Fock determinants.
The advantage of using Technology Performance Level (TPL) in conjunction with Technology Readiness Level (TRL) assessments in guiding technology development trajectories to successful outcomes in less time, at less overall cost, and with less encountered risk has been well articulated in the literature. In partnership with industry and international collaborators, a TPL assessment methodology for grid-connected applications has been developed through the application of the systems engineering approach. Metrics under seven different categories have been developed, weighted based on their relative relevance, and combined to yield a composite score. The methodology has been implemented in a spreadsheet tool plus a web application specifically aimed at assessing early stage (TRL 1-3) concepts. The target use cases are (a) technology developers improving their design, to find fatal flaws early, to get feedback on current design, to identify areas of improvement that will yield the highest return on investment, (b) reviewers assessing technologies in competitions or for making funding decisions, (c) investor or project developer doing due diligence, (d) policy makers landscaping the technology domain for formulating R&D strategy. The methodology and the tools are undergoing continuous improvement based on the experience and lessons learnt from applying it to internal and external marine energy technology development projects. The methodology is also being adapted for assessing WECs servicing markets outside the continental grid - broadly categorized as Powering the Blue Economy (PBE) applications. Such applications have vastly different functional requirements entailing a modification of the methodology to account for their higher risk tolerance, reduced price sensitivities, lower power needs, different permitting protocols, etc. This paper presents the latest status of the TPL assessment methodology and tools, describes its adaptation to select PBE markets, and explores its extension to other domains where it could provide a comprehensive and holistic measure of a nascent or disruptive technology's technoeconomic performance potential.