Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate computing”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

REYNO: A Reactive Hydrodynamics Modeling Suite

REYNO is a Python library that is primarily designed as a platform for the agile development and testing of novel reactive flow models. Equations of state and multi-step rate laws can be implemented using the provided class hierarchy with different thermodynamic closure conditions. The library includes a one-dimensional Lagrangian hydrodynamic solver to simulate multi-material reactive or inert problems in planar, cylindrical, or spherical coordinates. Custom time-dependent boundary conditions are also available for problems such as ramp compression. This document describes the different numerical methods and algorithms that are included in REYNO.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Novel Structural Motif To Promote Mg-Ion Mobility: Investigating ABO 4 Zircons as Magnesium Intercalation Cathodes

There is an increasing need for sustainable energy storage solutions as fossil fuels are replaced by renewable energy sources. Multivalent batteries, specifically Mg batteries, are one energy storage technology that researchers continue to develop with hopes to surpass the performance of Li-ion batteries. However, the limited energy density and transport properties of Mg cathodes remain critical challenges preventing the realization of high-performance multivalent batteries. In this work, ABO 4 zircon materials (A = Y, Eu and B = V, Cr) are computationally and experimentally evaluated as Mg intercalation cathodes. Remarkably good Mg-ion transport properties were predicted and Mg-ion intercalation was experimentally verified in sol-gel synthesized zircon YVO 4 , EuVO 4 , and EuCrO 4 . Among them, EuVO 4 exhibited the best electrochemical performance and demonstrated repeated reversible cycling. While we believe that the one-dimensional diffusion channels and redox-active species tetragonal coordination limit the value of many zircons as high-performance cathodes, their unique structural motif of overlapping polyhedra along the diffusion pathway appears instrumental for promoting good Mg-ion mobility. The motif results in a favorable "6-5-4" change in coordination that avoids unfavorable sites with lower coordination along the diffusion pathway and a structural design metric for future Mg cathode development.

25 ENERGY STORAGE↗
DOE OSTI2020

DESC: A stellarator equilibrium solver

In this paper, the new code DESC is presented to solve for fixed-boundary ideal magnetohydrodynamic equilibria in stellarators. The approach directly solves the equilibrium force balance as a system of nonlinear equations in the form f(x) = 0. The independent variables x represent nested magnetic flux surfaces expressed in the inverse representation with toroidal flux coordinates, and the equations f(x) quantify equilibrium force balance errors at discrete points in real space. Discretizing with global Fourier–Zernike basis functions properly treats the magnetic axis and minimizes the number of coefficients needed to describe the flux surfaces. The pseudospectral method provides great flexibility in where the errors are evaluated, and the system of equations is efficiently solved with a Newton–Raphson iteration. Equilibria are computed and compared against VMEC for both axisymmetric and non-axisymmetric examples. The results show fast convergence rates and solutions with low errors throughout the plasma volume.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Real-Time Regional PV Spinning Reserve Estimator with AGC Look-Ahead Windows

Curtailed PV generation is a zero-marginal cost spinning reserve that can be used for a number of active power control services. However, unlike the traditional spinning reserve providers, i.e., fossil-fueled generators, who have well-defined operating characteristics, e.g., available headroom or potential high limit (PHL), PV plants have by nature variable and uncertain operating characteristics. To ensure the effective coordination between PV plants and the system operator during an active power control event, accurate forecasts of the PV PHL are essential. A novel reference-control grouping based scaling method has been proposed by NREL to estimate the PV PHL in real-time. This work further enhances the methodology by: 1) improving the model accuracy through machine learning; 2) considering look-ahead windows introduced by the computation and communication latencies; 3) applying the method to regional spinning reserve estimation. A significant performance improvement, over 99% of estimation error reduction, has been observed based on real-world data collected by CAISO and PV plant operators.

potential high limit↗
DOE OSTI2023

Real-Time Regional PV Spinning Reserve Estimator with AGC Look-Ahead Windows: Preprint

Curtailed PV generation is a zero-marginal cost spinning reserve that can be used for a number of active power control services. However, unlike the traditional spinning reserve providers, i.e., fossil-fueled generators, who have well-defined operating characteristics, e.g., available headroom or potential high limit (PHL), PV plants have by nature variable and uncertain operating characteristics. To ensure the effective coordination between PV plants and the system operator during an active power control event, accurate forecasts of the PV PHL are essential. A novel reference-control grouping based scaling method has been proposed by NREL to estimate the PV PHL in real-time. This work further enhances the methodology by: 1) improving the model accuracy through machine learning; 2) considering look-ahead windows introduced by the computation and communication latencies; 3) applying the method to regional spinning reserve estimation. A significant performance improvement, over 1.6% of estimation error reduction, has been observed based on real-world data collected by CAISO and PV plant operators.

potential high limit↗
DOE OSTI2020

Impacts of hydrogen bonding interactions with Np(v/vi)O 2 Cl 4 complexes: vibrational spectroscopy, redox behavior, and computational analysis

The neptunyl (Np(V)O 2 + /Np(VI)O 2 2+ ) cation is the dominant form of 237 Np in acidic aqueous solutions and the stability of the Np(V) and Np(VI) species is driven by the specific chemical constituents present in the system. Hydrogen bonding with the oxo group may impact the stability of these species, but there is limited understanding of how these intermolecular interactions influence the behavior of both solution and solid-state species. In the current study, we systematically evaluate the interactions between the neptunyl tetrachloride species and hydrogen donors in coordination complexes and in the related aqueous solutions. Both Np(V) compounds (N 2 C 4 H 12 ) 2 [Np(V)O 2 Cl 4 ]Cl (Np(V)pipz) and (NOC 4 H 10 ) 3 [Np(V)O 2 Cl 4 ] (Np(V)morph) exhibit directional hydrogen bonding to the neptunyl oxo group while Np(VI) compounds (NC 5 H 6 ) 2 [Np(VI)O 2 Cl 4 ] (Np(VI)pyr) and (NOC 4 H 10 ) 4 [Np(VI)O 2 Cl 4 ]·2Cl (Np(VI)morph) assemble via halogen interactions. Here, the Raman spectra of the solid-state phases indicate the activation of vibrational bands when there is asymmetry of the neptunyl bond, while these spectral features are not observed within the related solution phase spectra. Density functional theory calculations of the Np(V)pipz system suggest that activation of the ν 3 asymmetric stretch and other combination modes lead to additional complexity within the solid-state spectra. Electrochemical analyses of complexes in the solution phases are consistent with the results of the crystallization experiments as the voltammetric potentials of Np(V)/Np(VI) complexes in the presence of protonated heterocycles differ from the potentials of pure Np(V) and may correlate with the hydrogen bonding interactions.

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

Integrate FARM with PID controllers: IES Simulation Ecosystem Control System Development

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM was designed to support the HERON software module in the solution of the optimal dispatch problem for IES units. As the result of HERON-FARM dispatch simulation, the set-point trajectories are optimized to meet constraints on both the production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and the process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.) at a coarse time resolution (every 10 or 100 seconds) over long time horizons (several days or weeks). In case the operational constraints need to be met at finer time resolution, the computational burden of HERON-FARM would linearly increase with the sampling rate, and sub-optimal solutions might be obtained. System responses characterized by overshoots and damped oscillations temporarily violating the imposed constraints might occur during abrupt power transients. In this report, a hierarchical control system architecture for the operation of the Dynamic Energy Technology and Integration Laboratory (DETAIL) facility constructed at INL was proposed. First, the preliminary studies on the proposed control strategy for operating the facility and the designed PI controllers were reviewed. In particular, the current approach for generating the set-point trajectories was studied, and its limits were identified. To this aim, the inclusion of a Supervisory Control layer embedding a modified version of the FARM algorithm for preserving the system safe operation over both long and real-time horizons was proposed. In this way, FARM would be applied twice, i.e., the original version (“FARM Validator”) aiding the solution of the power dispatch problem, and the modified version (“FARM Supervisory” coordinating the PI controllers to address the real-time control tasks. Despite the kernel of the two modules is the same algorithm, their roles, tasks, and capabilities are quite different. A detailed description of the role of FARM at addressing low-level control tasks is provided, along with tentative operational procedures for training the models embedded into the algorithm by using the collected experimental data.

42 ENGINEERING↗
DOE OSTI2025

Exploring the influence of transition metals on f-element bonding

The United States aims to triple its nuclear energy production by 2050, which will result in increased uranium usage and spent nuclear fuel generation. This highlights the need for a comprehensive understanding of actinide coordination chemistry, which is crucial for the extraction, processing, purification, and fabrication of uranium-based fuels. Moreover, it is vital for the reprocessing or safe disposal of nuclear waste and effective remediation efforts. Achieving this successfully requires an in-depth understanding of f-electron behavior, as “the role of 5f electrons in bond formation remains a fundamental topic in actinide chemistry”. Introducing a second metal into the system can increase structural dimensionality and diversify structural architecture. Heterometallic systems can also alter material properties, such as magnetic and spectroscopic characteristics, luminescence, and actinide mobility. Additionally, secondary transition metals, even when present only in the second coordination sphere and not directly coordinated, can influence the electron density at the actinyl metal center. In this study, uranium heterometallic single crystals were synthesized by incorporating transition metals such as iron(III), iron(II), nickel(II), manganese(II), copper(I), and cobalt(II). The crystals were formed using 2,6-pyridine dicarboxylic acid and other structurally similar ligands with varying functional groups. The synthesized crystals were analyzed using an extensive array of analytical and computational characterization techniques, including single crystal X-ray diffraction, Raman and infrared spectroscopy, and density functional theory calculations.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2021

Effects of Al:Si and (Al+Na):Si Ratios on the Properties of the International Simple Glass, Part II: Structure

High-alumina containing high-level waste (HLW) will be vitrified at the Waste Treatment Plant at the Hanford Site. The resulting glasses, high in alumina, will have distinct composition-structure-property (C-S-P) relationships compared to previously studied HLW glasses. These C-S-P relationships determine the processability and product durability of glasses and therefore must be understood. The main purpose of this study is to understand the detailed structural changes caused by Al:Si and (Al+Na):Si substitutions in a simplified nuclear waste model glass (ISG, International Simple Glass) by combining experimental structural characterizations and molecular dynamics (MD) simulations. The structures of these two series of glasses were characterized by neutron total scattering and 27Al, 23Na, 29Si, and 11B solid-state nuclear magnetic resonance (NMR) spectroscopy. Additionally, MD simulations were used to generate atomistic structural models of the borosilicate glasses and simulation results were validated by the experimental structural data. Short-range (e.g., bond distance, coordination number, etc.) and medium-range (e.g., oxygen speciation, network connectivity, polyhedral linkages) structural features of the borosilicate glasses were systematically investigated as a function of the degree of substitution. The results show that bond distance and coordination number of the cation-oxygen pairs are relatively insensitive to Al:Si and (Al+Na):Si substitutions with the exception of the B-O pair. Additionally, the Al:Si substitution results in an increase of tri-bridging oxygen species, while (Al+Na):Si substitution creates non-bridging oxygen species. Charge compensator preferences were found for Si-[NBO] (Na+), [3]B-[NBO] (Na+), [4]B (mostly Ca2+), [4]Al (nearly equally split Na+ and Ca2+), and [6]Zr (mostly Ca2+). The network former-BO-network former linkages preferences were also tabulated; Si-O-Al and Al-O-Al were preferred at the expense of lower Si-O-[3]B and [3]B-O-[3]B linkages. These results provide insights on the structural origins of property changes such as glass transition temperature caused by the substitutions, providing a basis for future improvements of theoretical and computer simulation models.

Lu, Xiaonan↗
DOE OSTI2023

Spatio-Temporal Surrogates for Interaction of a Jet with High Explosives: Part II - Clustering Extremely High-Dimensional Grid-Based Data

Building an accurate surrogate model for the spatio-temporal outputs of a computer simulation is a challenging task. A simple approach to improve the accuracy of the surrogate is to cluster the outputs based on similarity and build a separate surrogate model for each cluster. This clustering is relatively straightforward when the output at each time step is of moderate size. However, when the spatial domain is represented by a large number of grid points, numbering in the millions, the clustering of the data becomes more challenging. In this report, we consider output data from simulations of a jet interacting with high explosives. These data are available on spatial domains of different sizes, at grid points that vary in their spatial coordinates, and in a format that distributes the output across multiple files at each time step of the simulation. We first describe how we bring these data into a consistent format prior to clustering. Borrowing the idea of random projections from data mining, we reduce the dimension of our data by a factor of thousand, making it possible to use the iterative k-means method for clustering. We show how we can use the randomness of both the random projections, and the choice of initial centroids in k-means clustering, to determine the number of clusters in our data set. Our approach makes clustering of extremely high dimensional data tractable, generating meaningful cluster assignments for our problem, despite the approximation introduced in the random projections.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Significant acceleration of solid-state NMR simulations via three-angle powder averaging

The anisotropic frequency shifts imparted onto the NMR resonance frequency depend on the spherical angular coordinates that describe the orientations of the NMR interaction tensors with respect to the applied magnetic field direction. Experiments performed using magic-angle spinning, however, gain a dependence on a third angle: the rotor phase γ. Traditionally, a carousel average is performed to integrate over γ, which leads to a slow convergence of intensities without contributing to the underlying powder patterns. Herein, we show an order of magnitude acceleration in computation time may be obtained by including the γ-averaging into the main powder average to eliminate redundant calculation of resonance frequencies.

magic-angle spinning↗
DOE OSTI2025

Derivative-free stochastic optimization via adaptive sampling strategies

In this paper, we present a novel derivative-free framework for solving unconstrained stochastic optimization problems. Many problems in fields ranging from simulation optimization to reinforcement learning to quantum computing involve settings where only stochastic function values are obtained via a zeroth-order oracle, which has no available gradient information and necessitates the usage of derivative-free optimization methodologies. Our approach includes estimating gradients using stochastic function evaluations and integrating adaptive sampling techniques to control the accuracy in these stochastic approximations. Our framework encapsulates several gradient estimation techniques, including standard finite-difference, Gaussian smoothing, sphere smoothing, randomized coordinate finite-difference, and randomized subspace finite-difference methods. We provide theoretical convergence guarantees for our framework and analyze the worst-case iteration and sample complexities associated with each gradient estimation method. Finally, we demonstrate the empirical performance of the methods on logistic regression and nonlinear least squares problems.

Adaptive sampling↗
DOE OSTI2022

BETO 2021 Peer Review - 3.4.2.302 - Process Scale Up for Production Environments

This project supports the Department of Energy's Bioenergy Technology Office mission of transitioning bioenergy technologies to market by de-risking integration, developing scaling relations and modeling that de-risks scale up, and providing a feedback loop with lower TRL projects to understand process and technology fundamentals. Two objectives were pursued during this merit cycle. First, the Thermal and Catalytic Process Development Unit (TCPDU), a 0.5 ton/day pilot plant, was used to validate kinetic models developed by the Consortium for Computational Physics and Chemistry (CCPC). Validating the models accurately predict product yields and composition is a first step toward de-risking a common industry failure - taking too large of steps between scales. In addition, these models help with design and troubleshooting of new unit operations for the TCPDU. The second objective was to conduct the FY22 verification campaign around ex situ CFP. The primary challenge of this effort was the coordination of multiple national labs and projects. This was overcome by a dedicated leadership role and effective communication strategy between the projects. Currently the project is undergoing a pivot based on a stagegate decision to not conduct the verification campaign in the TCPDU and it is focused on closing out the CFP technology at a smaller scale by the end of the current fiscal year.

BIOMASS FUELS↗
DOE OSTI2024

Understanding Europium and Terbium Speciation and Ion Pairing in Carbonate Complexes Using Advanced Spectroscopy Techniques

Lanthanide (Ln) elements are critical materials that are typically extracted/mined together. Their separation by solvent extraction from acidic media is well known; however, there are few studies in basic media with carbonate anions. We investigated the complexation of Eu(III) and Tb(III) carbonates as solids and solutions in alkaline K2CO3, wherein we sought to access a Tb(IV) carbonate complex through ozonolysis. L3-edge XANES of Eu and Tb carbonate solids, colorless solutions, and a red-hued Tb solution (obtained by ozonolysis) all showed Ln(III) cations. The absence of evidence for a Tb(IV) complex was confirmed through XAS and EPR analyses, despite the solution exhibiting a deep red color. For solids and solutions, EXAFS results indicate molecular Ln(III)-carbonato anions. In terms of the Eu(III) carbonate coordination number, the coordination does not change upon dissolution of the solid sample. Furthermore, EXAFS for the solutions revealed evidence for the association of potassium cations with the Ln(III)-carbonato anions. Furthermore, this direct observation of contact ion pairing by EXAFS at room temperature is rare. The insights into Ln(III) carbonate complexation and solution speciation afforded by XANES-EXAFS, FT-IR, and EPR provides perspectives that serve as benchmarks for future computational and experimental efforts focused on caustic-side solvent extraction of Ln(III) ions.

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

Computing Potential of the Mean Force Profiles for Ion Permeation Through Channelrhodopsin Chimera, C1C2

Umbrella sampling, coupled with a weighted histogram analysis method (US-WHAM), can be used to construct potentials of mean force (PMFs) for studying the complex ion permeation pathways of membrane transport proteins. Despite the widespread use of US-WHAM, obtaining a physically meaningful PMF can be challenging. Here, we provide a protocol to resolve that issue. Then, we apply that protocol to compute a meaningful PMF for sodium ion permeation through channelrhodopsin chimera, C1C2, for illustration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Optimal Orbital Selection for Full Configuration Interaction (OptOrbFCI): Pursuing the Basis Set Limit under a Budget

Full configuration interaction (FCI) solvers are limited to small basis sets due to their expensive computational costs. An optimal orbital selection for FCI (OptOrbFCI) is proposed to boost the power of existing FCI solvers to pursue the basis set limit under a computational budget. The optimization problem coincides with that of the complete active space SCF method (CASSCF), while OptOrbFCI is algorithmically quite different. OptOrbFCI effectively finds an optimal rotation matrix via solving a constrained optimization problem directly to compress the orbitals of large basis sets to one with a manageable size, conducts FCI calculations only on rotated orbital sets, and produces a variational ground-state energy and its wave function. Coupled with coordinate descent full configuration interaction (CDFCI), we demonstrate the efficiency and accuracy of the method on the carbon dimer and nitrogen dimer under basis sets up to cc-pV5Z. We also benchmark the binding curve of the nitrogen dimer under the cc-pVQZ basis set with 28 selected orbitals, which provide consistently lower ground-state energies than the FCI results under the cc-pVDZ basis set. Furthermore, the dissociation energy in this case is found to be of higher accuracy.

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

Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product Report

Objective analysis (OA) is a method of mapping unstructured data to a structured grid. In the context of scanning radar data, OA is used to interpolate data in antenna coordinates (range, azimuth, and elevation) onto a regularly spaced Cartesian grid (Trapp and Doswell 2000). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product (VAP) uses the Python ARM Radar Toolkit (Py-ART), a data model-driven interactive architecture for working with weather radar data, to map the data to a Cartesian grid (Helmus and Collis 2016). MMCG, with Py-ART built in, has the ability to take radar data in antenna coordinates and map the gates to a Cartesian grid using inverse distance weight functions such as Cressman (square) and Barnes (exponential), but also can filter the data during the interpolation. MMCG also allows arbitrary formulations for the radius of influence, which are matched to particular radar scanning strategies. This creates a complex parameter space for optimizing the retention of storm structure detail while minimizing artifacts. MMCG takes data processed with ARM’s Corrected Precipitation Radar Moments in Antenna Coordinates (CMAC) VAP and maps it to a Cartesian grid as the output product. A variety of fields that have been mapped to the Cartesian grid are then saved to plots to complement each grid file.

54 ENVIRONMENTAL SCIENCES↗