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 397 records · Page 22

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Efficiently evaluating loop integrals in the EFTofLSS using QFT integrals with massive propagators

We develop a new way to analytically calculate loop integrals in the Effective Field Theory of Large Scale-Structure. Previous available methods show severe limitations beyond the one-loop power spectrum due to analytical challenges and computational and memory costs. Our new method is based on fitting the linear power spectrum with cosmology-independent functions that resemble integer powers of quantum field theory massive propagators with complex masses. A remarkably small number of them is sufficient to reach enough accuracy. Similarly to former approaches, the cosmology dependence is encoded in the coordinate vector of the expansion of the linear power spectrum in our basis. We first produce cosmology-independent tensors where each entry is the loop integral evaluated on a given combination of basis vectors. For each cosmology, the evaluation of a loop integral amounts to contracting this tensor with the coordinate vector of the linear power spectrum. The 3-dimensional loop integrals for our basis functions can be evaluated using techniques familiar to particle physics, such as recursion relations and Feynman parametrization. We apply our formalism to evaluate the one-loop bispectrum of galaxies in redshift space. The final analytical expressions are quite simple and can be evaluated with little computational and memory cost. We show that the same expressions resolve the integration of all one-loop N-point function in the EFTofLSS. This method, which is originally presented here, has already been applied in the first one-loop bispectrum analysis of the BOSS data to constraint ΛCDM parameters and primordial non-Gaussianities [1, 2].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Double Diffusion Maps and their Latent Harmonics for scientific computations in latent space

In this work, we introduce a data-driven approach to building reduced dynamical models through manifold learning; the reduced latent space is discovered using Diffusion Maps (a manifold learning technique) on time series data. A second round of Diffusion Maps on those latent coordinates allows the approximation of the reduced dynamical models. This second round enables mapping the latent space coordinates back to the full ambient space (what is called lifting); it also enables the approximation of full state functions of interest in terms of the reduced coordinates. In our work, we develop and test three different reduced numerical simulation methodologies, either through pre-tabulation in the latent space and integration on the fly or by going back and forth between the ambient space and the latent space. The data-driven latent space simulation results, based on the three different approaches, are validated through (a) the latent space observation of the full simulation through the Nyström Extension formula, or through (b) lifting the reduced trajectory back to the full ambient space, via Latent Harmonics. Latent space modeling often involves additional regularization to favor certain properties of the space over others, and the mapping back to the ambient space is then constructed mostly independently from these properties; here, we use the same data-driven approach to construct the latent space and then map back to the ambient space.

97 MATHEMATICS AND COMPUTING↗

Theoretical Assessment of the Transition Between Electron Emission Mechanisms for Nonplanar Diodes

Theoretically and computationally describing the operation of nanodiodes requires characterizing the transitions between multiple electron emission mechanisms for nanodiodes with complicated geometries. This motivates our development of techniques to determine when simplified theories for individual mechanisms suffice compared to more complete, but more computationally expensive, models. Leveraging recent theories that define a canonical gap distance to translate planar theory to nonplanar diodes, we derive the conditions for the transitions among thermal emission, field emission, and space-charge-limited current density (SCLCD) in vacuum and with collisions for non-Cartesian coordinate systems, including spherical, cylindrical, and prolate spheroidal coordinate systems. Particle-in-cell (PIC) simulations of the current density as a function of applied voltage for a tip-to-plate geometry in vacuum agreed qualitatively with the asymptotes for thermal emission at low voltage and SCLCD at higher voltage using the canonical gap distance. As a result, this demonstrates the utility of this approach for guiding system design and suggests future extensions to save simulation time for more realistic geometries that are more computationally expensive.

Conformal mapping↗

Mobile Sensing for Wind Field Estimation in Wind Farms: Preprint

This paper introduces a novel approach for estimating the wind field over an entire wind farm using a mobile sensor to collect limited amounts of data. The proposed method estimates the boundary conditions of a simplified turbine wake model by computing the model sensitivity matrix and using a recursive least-squares algorithm to recover the model parameters from the wind field measurements. To address the fact that it is not practical to take measurements across the entire wind farm, the proposed method classifies each area on the map based on its sensitivity to parameter variations. This classification is then used to generate a suitable path for a mobile sensor, which is charged with collecting data for the recursive least-squares algorithm. The proposed framework can successfully estimate the model boundary conditions using just the measurements collected along the path of the mobile sensor. This preliminary result paves the way for using real-time wind field estimates for the coordinated control of all the turbines within a wind farm.

estimation↗

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Enabling Modular Autonomous Feedback‐Loops in Materials Science through Hierarchical Experimental Laboratory Automation and Orchestration

Abstract Materials acceleration platforms (MAPs) operate on the paradigm of integrating combinatorial synthesis, high‐throughput characterization, automatic analysis, and machine learning. Within a MAP, one or multiple autonomous feedback loops may aim to optimize materials for certain functional properties or to generate new insights. The scope of a given experiment campaign is defined by the range of experiment and analysis actions that are integrated into the experiment framework. Herein, the authors present a method for integrating many actions within a hierarchical experimental laboratory automation and orchestration (HELAO) framework. They demonstrate the capability of orchestrating distributed research instruments that can incorporate data from experiments, simulations, and databases. HELAO interfaces laboratory hardware and software distributed across several computers and operating systems for executing experiments, data analysis, provenance tracking, and autonomous planning. Parallelization is an effective approach for accelerating knowledge generation provided that multiple instruments can be effectively coordinated, which the authors demonstrate with parallel electrochemistry experiments orchestrated by HELAO. Efficient implementation of autonomous research strategies requires device sharing, asynchronous multithreading, and full integration of data management in experimental orchestration, which to the best of the authors’ knowledge, is demonstrated for the first time herein.

36 MATERIALS SCIENCE↗

Realizing Wide-Gamut Human-Centric Display Lighting with K 3 AlP 3 O 9 N:Eu 2+

Abstract Computers, televisions, and smartphones are revolutionized by the invention of InGaN blue light-emitting diode (LED) backlighting. Yet, continual exposure to the intense blue LED emission from these modern displays can cause insomnia and mood disorders. Developing “human-centric” backlighting that uses a violet-emitting LED chip and a trichromatic phosphor mixture to generate color images is one approach that addresses this problem. The challenge is finding a blue-emitting phosphor that possesses a sufficiently small Stokes’ shift to efficiently down-convert violet LED light and produce a narrow blue emission. This work reports a new oxynitride phosphor that meets this demand. K 3 AlP 3 O 9 N:Eu 2+ exhibits an unexpectedly narrow (45 nm, 2206 cm -1 ), thermally robust, and efficient blue photoluminescence upon violet excitation. Computational modeling and temperature-dependent optical property measurements reveal that the narrow emission arises from a rare combination of preferential excitation and site-selective quenching. The resulting chromaticity coordinates of K 3 AlP 3 O 9 N:Eu 2+ lie closer to the vertex of the Rec. 2020 than a blue LED chip and provides access to ≈10% more colors than a commercial tablet when combined with commercial red- and green-emitting phosphors. Alongside the wide gamut, tuning the emission from the violet LED and phosphor blend can reduce blue light emissions to produce next-generation, human-centric displays.

36 MATERIALS SCIENCE↗

Grain Boundary‐Derived Cu + /Cu 0 Interfaces in CuO Nanosheets for Low Overpotential Carbon Dioxide Electroreduction to Ethylene

Abstract Electrochemical CO 2 reduction reaction can be used to produce value‐added hydrocarbon fuels and chemicals by coupling with clean electrical energy. However, highly active, selective, and energy‐efficient CO 2 conversion to multicarbon hydrocarbons, such as C 2 H 4 , remains challenging because of the lack of efficient catalysts. Herein, an ultrasonication‐assisted electrodeposition strategy to synthesize CuO nanosheets for low‐overpotential CO 2 electroreduction to C 2 H 4 is reported. A high C 2 H 4 Faradaic efficiency of 62.5% is achieved over the CuO nanosheets at a small potential of −0.52 V versus a reversible hydrogen electrode, corresponding to a record high half‐cell cathodic energy efficiency of 41%. The selectivity toward C 2 H 4 is maintained for over 60 h of continuous operation. The CuO nanosheets are prone to in situ restructuring during CO 2 reduction, forming abundant grain boundaries (GBs). Stable Cu + /Cu 0 interfaces are derived from the low‐coordinated Cu atoms in the reconstructed GB regions and act as highly active sites for CO 2 reduction at low overpotentials. In situ Raman spectroscopic analysis and density functional theory computation reveal that the Cu + /Cu 0 interfaces offer high *CO surface coverage and lower the activation energy barrier for *CO dimerization, which, in synergy, facilitates CO 2 reduction to C 2 H 4 at low overpotentials.

36 MATERIALS SCIENCE↗

Surface-Controlled TiO 2 Nanocrystals with Catalytically Active Single-Site Co Incorporation for the Oxygen Evolution Reaction

The design of advanced electrocatalysts is often hindered by uncertainties in identifying and controlling the active surfaces and catalytic centers within heterogeneous materials. Here we present the synthesis of single-site Co catalysts, substitutionally doped into surface-controlled TiO 2 anatase nanocrystals, aimed at enhancing the oxygen evolution reaction (OER). Grand canonical quantum mechanics calculations reveal that the kinetics of the OER, following an adsorbate evolution mechanism, is markedly influenced by the coordination environment of Co. The simulations suggest significantly higher turnover frequencies when Co is doped into the (001) surface of TiO 2 compared to the (101) surface. Consistent with the computational findings, experimental results show that Co-doped TiO 2 (Co-TiO 2 ) nanoplates with selectively exposed {001} surfaces exhibit enhanced current densities and turnover frequencies compared to Co-TiO 2 nanobipyramids with {101} surfaces. This study highlights the synergy between theoretical calculations and precision synthesis in the development of more effective catalysts.

25 ENERGY STORAGE↗

Symmetry Determining Equations of the Euler Equations

This purpose of this document is to provide a clear derivation of the symmetry determining equations (SDEs) of the Euler equations to supplement the lecture document. The Euler equations are a coupled system of three first-order quasilinear hyperbolic partial differential equations (PDEs). The Euler equations are a system of partial differential equations (PDEs) for fluid velocity v(t, r), fluid density ρ(t, r), and fluid pressure P(t, r), where t is time and r is the spatial coordinate.

97 MATHEMATICS AND COMPUTING↗

Exascale workflow applications and middleware: An ExaWorks retrospective

Exascale computers offer transformative capabilities to combine data-driven and learning-based approaches with traditional simulation applications to accelerate scientific discovery and insight. However, these software combinations and integrations are difficult to achieve due to the challenges of coordinating and deploying heterogeneous software components on diverse and massive platforms. Here, we present the ExaWorks project, which addresses many of these challenges. We developed a workflow Software Development Toolkit (SDK), a curated collection of workflow technologies that can be composed and interoperated through a common interface, engineered following current best practices, and specifically designed to work on HPC platforms. ExaWorks also developed PSI/J, a job management abstraction API, to simplify the construction of portable software components and applications that can be used over various HPC schedulers. The PSI/J API is a minimal interface for submitting and monitoring jobs and their execution state across multiple and commonly used HPC schedulers. We also describe several leading and innovative workflow examples of ExaWorks tools used on DOE leadership platforms. Furthermore, we discuss how our project is working with the workflow community, large computing facilities, and HPC platform vendors to address the requirements of workflows sustainably at the exascale.

97 MATHEMATICS AND COMPUTING↗

Machine Learning-Based PV Reserve Determination Strategy for Frequency Control on the WECC System

Frequency control from photovoltaic (PV) power plants has great potential to address the frequency response challenge of the power system with high penetrations of renewable generation. Using model-based approaches to determine the optimal PV headroom reserve, however, requires significant online computation and is intractable for an interconnection level system. This paper proposes a machine learning based strategy, that is suitable for real-time operation, to determine the optimal PV reserve for frequency control. The proposed machine learning algorithm is trained and tested on 1,987 offline simulations of a 60% renewable penetration Western Electricity Coordinating Council (WECC) system. Furthermore, the proposed reserve determination strategy is applied on a realistic 1-day operation profile of the WECC system and demonstrates a savings of more than 40% PV headroom compared to a conservative approach. It is evident that the proposed strategy can efficiently and effectively determine the optimal PV frequency control reserve for realistic interconnection systems.

frequency control↗

Validation of a Tapered Impact Damper for Traffic Signal Structure Vibrations Using the Method of Harmonic Balance

Significant vibration amplitudes and cycles can be produced when traffic signal structures with low inherent damping are excited near one of their natural frequencies. For the mitigation of wind-induced vibrations, dynamic vibration absorbers coupled to the structure are often used. Here, this research investigates the performance of a tapered impact damper, consisting of a hanging spring-mass oscillator inside a housing capable of reducing vibration amplitude over a broader frequency range than the conventional tuned mass damper. A nonlinear, two degree-of-freedom model is developed with coordinates representing the traffic structure and the tapered impact damper. This research focuses on the application of the harmonic balance method to approximate the periodic solutions of the nonlinear equations to compute the nonlinear dynamics of the damped traffic signal structure. After designing and manufacturing a tapered impact damper, the traffic signal structure is tested with and without the damper using free vibration snapback tests. The experimental frequency and damping backbone curves are used to validate the analytical model, and the effectiveness of the damper is discussed.

42 ENGINEERING↗

Spatially separating redox centers on 2D carbon nitride with cobalt single atom for photocatalytic H 2 O 2 production

Redox cocatalysts play crucial roles in photosynthetic reactions, yet simultaneous loading of oxidative and reductive cocatalysts often leads to enhanced charge recombination that is detrimental to photosynthesis. This study introduces an approach to simultaneously load two redox cocatalysts, atomically dispersed cobalt for improving oxidation activity and anthraquinone for improving reduction selectivity, onto graphitic carbon nitride (C 3 N 4 ) nanosheets for photocatalytic H 2 O 2 production. Spatial separation of oxidative and reductive cocatalysts was achieved on a two-dimensional (2D) photocatalyst, by coordinating cobalt single atom above the void center of C 3 N 4 and anchoring anthraquinone at the edges of C 3 N 4 nanosheets. Such spatial separation, experimentally confirmed and computationally simulated, was found to be critical for enhancing surface charge separation and achieving efficient H 2 O 2 production. This center/edge strategy for spatial separation of cocatalysts may be applied on other 2D photocatalysts that are increasingly studied in photosynthetic reactions.

2D photocatalyst↗

Simulating water adsorption in metal–organic frameworks with open metal sites using the 12-6-4 Lennard–Jones potential

Metal-organic frameworks (MOFs) with coordinatively unsaturated open metal sites (OMSs) are promising sorbent materials in atmospheric water harvesting (AWH) systems at low relative humidity (RH) due to their strong interactions with water molecules. However, accurate computational modelling of water adsorption in those materials are challenging, as standard force fields (FFs) based on the 12-6 Lennard-Jones (L-J) potential cannot properly describe the water-OMS interactions. In biomolecular simulations, the 12-6-4 L-J potential has been successfully used to model metal ion-water interactions in which the extra 1/r 4 -term is for charge-induced dipole electrostatics. In this work, we adopted this strategy and used the 12-6-4 L-J potential to model water-OMS interactions in MOFs for low RH AWH applications. Notably, without any modifications, the parameters used for aqueous metal ions were able to greatly improve the accuracy in predicting water adsorption isotherms in MOF-74, which highlights the simplicity and transferability of this method.

74 ATOMIC AND MOLECULAR PHYSICS↗

An Efficient Distributed Reinforcement Learning for Enhanced Multi-Microgrid Management

Economic dispatch in multi-microgrid (MMG) systems requires coordinating distributed energy resources (DERs) of different microgrids, which leads to a significant increase in the number of states for energy management. In these cases, traditional reinforcement learning (RL) approaches become computationally expensive or output a solution that causes extra-operating costs for the system. This paper proposes an RL approach that employs local learning agents to interact with microgrid environments in a distributed manner and aggregates the outcomes to train the global agent to learn the policy for the MMG system. This distributed exploration and aggregation process provides an effective solution and guides the global agent to learn the dispatch policy efficiently. Case studies are performed on a system with three microgrids with different types of DERs. Results obtained using the proposed RL and comparisons with conventional methods substantiate the effectiveness of the proposed approach in terms of operation costs, computation time, and peak-to-average ratio.

Das, Avijit↗

Positive-Sequence Modeling of Droop-Controlled Grid-Forming Inverters for Transient Stability Simulation of Transmission Systems

Here, this paper describes a positive-sequence model to represent two widely reported droop-controlled grid-forming inverters for bulk power systems simulation study. Methods of how to develop the equivalent voltage source behind impedance to represent inverters with and without inner control loops, modeling of P-f and Q-V droop controls, active and reactive power limiting modeling, and modeling of fault current limiting controls are described in detail. The model is implemented in commercially available positive-sequence simulation tools and compared with detailed electromagnetic transient models in both a single-grid-forming-inverter infinite-bus system and a modified IEEE 39-bus system. Finally, the model is tested on the full US Western Electricity Coordinating Council (WECC) system. Study results show that the model has a good level of accuracy compared to the detailed EMT models, and at the same time also achieves a high computational efficiency, which is suitable for the bulk power system simulation study.

42 ENGINEERING↗