Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

General overview of multiphysics modeling and simulation at CEA/DES/IRESNE supporting safety, operating and development of nuclear reactors and facilities

This paper presents a general overview of the innovative multiphysics research/development activities carried on by the CEA/DES/IRESNE institute. Eight main different multiphysics fields addressed by IRESNE have been listed (among others: fuel fabrication, behavior of fuels, of reactors, of corium, radionuclide transfer to environment,...) in relation with the involved coupled physics: neutron transport, fluid mechanics, electro-magnetism, heat transfer, solid mechanics and (physico-)chemistry. The major mission of IRESNE institute is to study current and future nuclear systems integrated into a low carbon energy system. To fulfill this mission, IRESNE's activities rely on the development and the implementation of multiphysics modeling, advanced coupling scheme, data transfer techniques and dedicated multiphysics experimentations. Indeed, such multiphysics simulation approach of nuclear systems is today required to deal with the complex intrinsic features of studied phenomena. This approach aims at taking into account physical phenomena interactions, at each interest modeling scale, in order to reach an accurate and predictive behavior representation of a system component or of the whole system. The final objective is to contribute to the development of digital reactor, core and fuel at CEA, by progressing on applied mathematics, high-performance computing and physical mechanism understanding.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

PDG: A Composite Method Based on the Resolution of the Identity

The Gaussian-3 (G3) composite approach for thermochemical properties is revisited in light of the enhanced computational efficiency and reduced memory costs by applying the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs) to the computationally demanding component methods in the G3 model: the energy and gradient computations via the second-order Møller–Plesset perturbation theory (MP2) and the energy computations using the coupled-cluster singles–doubles method augmented with noniterative triples corrections [CCSD(T)]. Efficient implementation of the RI-based methods is achieved by employing a hybrid distributed/shared memory model based on MPI and OpenMP. The new variant of the G3 composite approach based on the RI approximation is termed the RI-G3 scheme, or alternatively the PDG method. The accuracy of the new RI-G3/PDG scheme is compared to the “standard” G3 composite approach that employs the memory-expensive four-center ERIs in the MP2 and CCSD(T) calculations. Taking the computation of the heats of formation of the closed-shell molecules in the G3/99 test set as a test case, it is demonstrated that the RI approximation introduces negligible changes to the mean absolute errors relative to the standard G3 model (less than 0.1 kcal/mol), while the standard deviations remain unaltered. Furthermore, the efficiency and memory requirements for the RI-MP2 and RI-CCSD(T) methods are compared to the standard MP2 and CCSD(T) approaches, respectively. The hybrid MPI/OpenMP-based RI-MP2 energy plus gradient computation is found to attain a 7.5× speedup over the standard MP2 calculations. For the most demanding CCSD(T) calculations, the application of the RI approximation is found to nearly halve the memory demand, confer about a 4–5× speedup for the CCSD iterations, and reduce the computational time for the compute-intensive triples correction step by several hours.

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

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (N k ) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, N ISDF ) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of N ISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. In conclusion, we perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

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

Magnetic-Field-Driven Electron Dynamics in Graphene

Graphene exhibits unique optoelectronic properties originating from the band structure at the Dirac points. It is an ideal model structure to study the electronic and optical properties under the influence of the applied magnetic field. In graphene, electric field, laser pulse, and voltage can create electron dynamics which is influenced by momentum dispersion. However, computational modeling of momentum-influenced electron dynamics under the applied magnetic field remains challenging. In this work, we perform computational modeling of the photoexcited electron dynamics achieved in graphene under an applied magnetic field. Our results show that magnetic field leads to local deviation from momentum conservation for charge carriers. With the increasing magnetic field, the delocalization of electron probability distribution increases and forms a cyclotron-like trajectory. Our work facilitates understanding of momentum resolved magnetic field effect on non-equilibrium properties of graphene, which is critical for optoelectronic and photovoltaic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Seven open problems in applied combinatorics

We present and discuss seven different open problems in applied combinatorics. Additionally, the application areas relevant to this compilation include quantum computing, algorithmic differentiation, topological data analysis, iterative methods, hypergraph cut algorithms, and power systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Gaussian process autoregressive model capturing microstructure evolution paths in a Ni–Mo–Nb alloy

Additive manufacturing is increasingly being employed to produce components of complex geometries in structural alloys because of the expected energy savings associated with the near-net-shape capability and the ability to build in novel internal features that are not possible with many conventional manufacturing approaches. However, because of the extreme thermal conditions encountered, the non-equilibrium microstructures produced during powder bed-based additive manufacturing processes must be subjected to custom post-heat treatment processes to recover the target mechanical properties. Phase-field models and simulation techniques have matured to a state where the microstructure evolution paths, and the morphologies of the resulting precipitate phases can be predicted reasonably accurately, considering alloy-specific thermodynamic and kinetic aspects of the nucleation and growth processes. However, phase-field simulations are computationally intensive, which precludes the ability to apply the simulations directly to the length scale of the entire component. Therefore, it is highly desirable to develop low-computational-cost surrogate models that effectively capture the physics at the microstructural length scale, while facilitating the design of optimized processing conditions resulting in location-specific targeted microstructures at the component scale. The work presented here demonstrates the application of the materials knowledge system framework to develop a surrogate model that effectively captures the microstructural path during annealing of a Ni–Mo–Nb alloy containing different Mo and Nb compositions known to segregate during solidification under additive manufacturing conditions. Specifically, the surrogate model built in this work is based on a Gaussian process autoregressive model informed by statistical representation of simulated microstructures using two-point correlations and dimensionality reduction through principal component analysis. In conclusion, this surrogate model is shown to capture the bifurcation of the microstructural path during precipitation, which yields a microstructure dominated by the $\gamma^{\prime\prime}$ phase at high Nb concentrations and the $\delta$ phase at low Nb concentrations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Construction of 3D MHD pressure drop correlation and flow characterization in the contraction region of a fusion blanket manifold

Inlet and outlet manifolds are typical components of liquid metal (LM) blanket designs of a fusion power reactor to be used to distribute the LM flow into breeding channels and collect it at the exit of the blanket. High pressure loss in the magnetohydrodynamic (MHD) flows featuring abrupt geometrical changes is one of the main feasibility issues of such designs. Recently, optimization studies were conducted to construct 3D MHD pressure drop correlations for a LM flow in an electrically insulating manifold with gradual expansion. Here, the 3D computational approach developed in that study is applied to the outlet manifold featuring gradual contraction. A systematic analysis was performed with a total number of 135 flow cases computed with COMSOL Multiphysics for Hartmann numbers 1000 < Ha < 10,000, Reynolds numbers 100 < Re < 12,000, and contraction angles 45° < θ < 75° for a fixed contraction ratio of 4. The effects of Ha, Re and θ on the flow recirculation, development length and the total pressure drop were carefully examined. A linear regression analysis was used to determine the power rule of pressure drop coefficient k related to Ha and Re, demonstrating a good match with the Ludford layer theory. Eventually, a correlation for the 3D MHD pressure drop coefficient was constructed as a function of Ha, Re and θ. Further, the results were compared against the inlet manifold. It was found that the flow in the inlet manifold exhibits larger recirculation zones. In the investigated range of Ha, Re and θ, the pressure drop coefficient k of the LM MHD flow in the gradual contraction is only slightly lower (< 8 %) than that in the gradual expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Mechanistic insights into metal, nitrogen doped carbon catalysts for oxygen reduction: progress in computational modeling

We report metal and nitrogen doped carbon materials (denoted as M–N–C) synthesized through high-temperature pyrolysis have been found to exhibit activity for oxygen reduction reaction (ORR) approaching that of Pt and electrochemical stability higher than previous MN 4 -containing macrocyclic molecular catalysts. Tremendous efforts have thus been devoted to the advancement of M–N–C catalysts as an economical alternative to Pt-based catalysts for proton exchange membrane fuel cell cathodes with a focus on simultaneously improving activity and stability. To this end, novel computational modeling techniques have been developed and applied to acquire knowledge crucial for accelerating the pace of M–N–C catalyst development. In this review, recent progress in computational method development, as well as the predictions of chemical structure of active sites, reaction pathways, ORR kinetics, and catalyst stability in electrochemical environments, are critically surveyed. Moreover, the crucial role of computational modeling to elucidate the functional mechanism of M–N–C catalysts for ORR in acid media and enable rational design of M–N–C catalysts is discussed with a visionary outlook for the field.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Reinforcement learning in discrete action space applied to inverse defect design

Abstract Reinforcement learning (RL) algorithms that include Monte Carlo Tree Search (MCTS) have found tremendous success in computer games such as Go, Shiga and Chess. Such learning algorithms have demonstrated super-human capabilities in navigating through an exhaustive discrete action search space. Motivated by their success in computer games, we demonstrate that RL can be applied to inverse materials design problems. We deploy RL for a representative case of the optimal atomic scale inverse design of extended defects via rearrangement of chalcogen (e.g. S) vacancies in 2D transition metal dichalcogenides (e.g. MoS 2 ). These defect rearrangements and their dynamics are important from the perspective of tunable phase transition in 2D materials i.e. 2H (semi-conducting) to 1T (metallic) in MoS 2 . We demonstrate the ability of MCTS interfaced with a reactive molecular dynamics simulator to efficiently sample the defect phase space and perform inverse design—starting from randomly distributed S vacancies, the optimal defect rearrangement of defects corresponds a line defect of S vacancies. We compare MCTS performance with evolutionary optimization i.e. genetic algorithms and show that MCTS converges to a better optimal solution (lower objective) and in fewer evaluations compared to GA. We also comprehensively evaluate and discuss the effect of MCTS hyperparameters on the convergence to solution. Overall, our study demonstrates the effectives of using RL approaches that operate in discrete action space for inverse defect design problems.

42 ENGINEERING↗
DOE OSTI2026

Smokescreen: A Python package for data vector blinding and encryption in cosmological analyses

Smokescreen is an open-source Python library for data-vector concealment (blinding) in cosmological analyses. Data-vector blinding works by applying cosmology-dependent shifts to the observed data vector, moving it away from the true cosmological signal without affecting its statistical properties, so that analysts cannot infer the true result until the analysis is frozen and the blinding is lifted. The package computes these shifts using Firecrown likelihoods applied to data vectors stored in the SACC format, ensuring that the theoretical model used for blinding is identical to that used for inference whilst remaining agnostic to the specific observable being blinded. To prevent accidental unblinding, the original SACC file, containing the true cosmology, is encrypted. Although developed for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), Smokescreen is applicable to any experiment using Firecrown likelihoods and the SACC data format.

Loureiro, Arthur [Stockholm U., OKC; Imperial Coll↗
DOE OSTI2021

Optimizing multigrid reduction-in-time and Parareal coarse-grid operators for linear advection

Parallel-in-time methods, such as multigrid reduction-in-time (MGRIT) and Parareal, provide an attractive option for increasing concurrency when simulating time-dependent partial differential equations (PDEs) in modern high-performance computing environments. While these techniques have been very successful for parabolic equations, it has often been observed that their performance suffers dramatically when applied to advection-dominated problems or purely hyperbolic PDEs using standard rediscretization approaches on coarse grids. In this paper, we apply MGRIT or Parareal to the constant-coefficient linear advection equation, appealing to existing convergence theory to provide insight into the typically nonscalable or even divergent behavior of these solvers for this problem. To overcome these failings, we replace rediscretization on coarse grids with improved coarse-grid operators that are computed by applying optimization techniques to approximately minimize error estimates from the convergence theory. Therefore, one of our main findings is that, in order to obtain fast convergence as for parabolic problems, coarse-grid operators should take into account the behavior of the hyperbolic problem by tracking the characteristic curves. Our approach is tested for schemes of various orders using explicit or implicit Runge–Kutta methods combined with upwind-finite-difference spatial discretizations. In all cases, we obtain scalable convergence in just a handful of iterations, with parallel tests also showing significant speed-ups over sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Smoky Mountain Data Challenge 2021: An Open Call to Solve Scientific Data Challenges Using Advanced Data Analytics and Edge Computing

The 2021 Smoky Mountains Computational Sciences and Engineering Conference enlists scientists from across Oak Ridge National Laboratory (ORNL) and industry to be data sponsors and help create data analytics and edge computing challenges for eminent datasets in a variety of scientific domains. This work describes the significance of each of the eight datasets and their associated challenge questions. The challenge questions for each dataset were required to cover multiple difficulty levels. An international call for participation was sent to students, asking them to form teams of up to six people and apply novel data analytics and edge computing methods to solve these challenges.

Devineni, Pravallika↗
DOE OSTI2024

Fluid-structure interaction solver for transient dynamics of fracturing media

Computer implemented methods include providing a numerical simulation solid mesh comprising a plurality of solid nodes and a numerical simulation fluid mesh comprising a plurality of fluid nodes, and performing a numerical simulation with the solid mesh and fluid mesh including: determining a fluid-solid interaction between the fluid mesh and the solid mesh at a time step by (i) immersing boundaries of the solid and fluid meshes, (ii) computing current corrected fluid and solid velocities using an interaction parameter applied to an intermediate fluid-solid relative velocity, and (iii) computing a fluid-solid interaction force based on the current corrected fluid and solid velocities.

Rougier, Esteban↗
DOE OSTI2020

Utilizing Essential Symmetry Breaking in Auxiliary-Field Quantum Monte Carlo: Application to the Spin Gaps of the C 36 Fullerene and an Iron Porphyrin Model Complex

We present three distinct examples where phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) can be reliably performed with a single-determinant trial wave function with essential symmetry breaking. Essential symmetry breaking was first introduced by Lee and Head-Gordon [ Phys. Chem. Chem. Phys. 2019, 21, 4763-4778, 10.1039/C8CP07613H]. We utilized essential complex and time-reversal symmetry breaking with ph-AFQMC to compute the triplet–singlet energy gap in the TS12 set. We found statistically better performance of ph-AFQMC with complex-restricted orbitals than with spin-unrestricted orbitals. We then showed the utilization of essential spin symmetry breaking when computing the singlet–triplet gap of a known biradicaloid, C 36 . ph-AFQMC with spin-unrestricted Hartree–Fock (ph-AFQMC+UHF) fails catastrophically even with spin-projection and predicts no biradicaloid character. With approximate Brueckner orbitals obtained from regularized orbital-optimized second-order Møller–Plesset perturbation theory (κ-OOMP2), ph-AFQMC quantitatively captures strong biradicaloid character of C 36 . Lastly, we applied ph-AFQMC to the computation of the quintet–triplet gap in a model iron porphyrin complex where brute-force methods with a small active space fail to capture the triplet ground state. We show unambiguously that neither triplet nor quintet is strongly correlated using UHF, κ-OOMP2, and coupled-cluster with singles and doubles (CCSD) performed on UHF and κ-OOMP2 orbitals. There is no essential symmetry breaking in this problem. By virtue of this, we were able to perform UHF+ph-AFQMC reliably with a cc-pVTZ basis set and predicted a triplet ground state for this model geometry. The largest ph-AFQMC in this work correlated 186 electrons in 956 orbitals. Our work highlights the utility, scalability, and accuracy of ph-AFQMC with a single-determinant trial wave function with essential symmetry breaking for systems mainly dominated by dynamical correlation with little static correlation.

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

Knowledge gaps for neuromorphic ionic computing

BACKGROUND Neuromorphic computing, inspired by the human brain’s ability to process information efficiently, represents a transformative approach to computation. In this Review, we explore the emerging field of neuromorphic ionic computing, which leverages ionic conduction and coupling to mimic neural processes, and identify critical knowledge gaps that must be addressed to realize its full potential. A central theme of the discussion is energy efficiency, a challenge that is both a limitation and an opportunity for this technology. Although complementary metal-oxide semiconductor (CMOS)–based neuromorphic technologies have made strides in scaling to billions of neurons and are increasingly applied in artificial intelligence and numerical computing, they remain orders of magnitude behind the human brain in terms of connectivity and energy efficiency. Neuromorphic ionic computing promises to overcome these limitations by leveraging the distinct architectural and operational principles of the brain. Our brains achieve this energy efficiency by combining several key features: using the same network elements to store and process information; using an incredibly complex and massively interconnected three-dimensional (3D) network of locally active elements that enables sparsity, robustness in the presence of noise, adaptation, and life-long learning; computing at comparatively low voltage and frequency; and last, taking advantage of a plethora of ions and small molecules as information carriers. Here, we propose that ionic computing systems can take advantage of similar features to achieve substantial gains in energy efficiency. ADVANCES Since the first reports of neuromorphic ionic behavior in nanofluidic channels, we have witnessed an explosion of reports that used ionic devices to produce synaptomimetic behaviors. However, achieving the goals of ionic computing requires not only implementation of much more sophisticated device functionality but also overcoming fundamental barriers in materials science, device architecture, and system integration. Current ionic devices, even those incorporating state-of-the-art materials, still suffer from limited functionality and stability, which restrict their performance and increase energy demands. Developing new materials with enhanced ionic properties is essential to overcome these limitations. Similarly, the design of neuromorphic devices must evolve to leverage the particular advantages of ionic processes. Existing architectures often follow a single-information-carrier logic of conventional electronics or are constructed of mesoscale fluidics, failing to capitalize on the energy-efficient mechanisms inherent to ionic systems or implement the multiple-information-carrier paradigm. Current neuromorphic chips focus on large-scale networks of analog memory elements based on mechanisms such as charge trap (flash), filamentary, phase change, or spin, which are built on top of a network of artificial CMOS neurons. Although such prototype networks have achieved impressive performance, it is difficult to envision how they can implement the key features such as massive connectivity, sophisticated plasticity, adaptability, sparsity, and “multichromatic” computing. Although small-scale devices have demonstrated promising results, integrating them, maintaining energy efficiency, and implementing temperature control as systems grow in complexity and size to computationally relevant scale remain major hurdles. Furthermore, interfacing neuromorphic ionic devices with existing computing technologies presents technical and conceptual challenges that will require innovative approaches that combine insights from neuroscience, materials science, and engineering. OUTLOOK Despite these challenges, the potential impact of neuromorphic ionic computing is profound with potential applications ranging from artificial intelligence to robotics and beyond. We also argue that neuromorphic ionic computing systems should not, at least in the beginning, compete with CMOS technologies but rather should focus on applications that require extreme energy efficiency with chemical and/or biological compatibility, such as biomedical applications (for example, brain-computer interfaces), environmental monitoring, and agricultural and food applications. Ultimately, this Review highlights the crucial role of interdisciplinary collaboration in advancing the field. Neuromorphic ionic computing is not merely a technological innovation; it represents a substantial step toward sustainable computation, aligning with the growing demand for energy-conscious solutions in a world that is increasingly reliant on data and computation.

Neuromorphic↗
DOE OSTI2023

I/O Access Patterns in HPC Applications: A 360-Degree Survey

The high-performance computing I/O stack has been complex due to multiple software layers, the inter-dependencies among these layers, and the different performance tuning options for each layer. In this complex stack, the definition of an “I/O access pattern” has been reappropriated to describe what an application is doing to write or read data from the perspective of different layers of the stack, often comprising a different set of features. It has become common to have to redefine what is meant when discussing a pattern in every new study, as no assumption can be made. This survey aims to propose a baseline taxonomy, harnessing the I/O community’s knowledge over the past 20 years. This definition can serve as a common ground for high-performance computing I/O researchers and developers to apply known I/O tuning strategies and design new strategies for improving I/O performance. We seek to summarize and bring a consensus to the multiple ways to describe a pattern based on common features already used by the community over the years.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Machine learning for arbitrary single-qubit rotations on an embedded device

Here, in this study, we present a technique for using machine learning (ML) for single-qubit gate synthesis on field-programmable logic for a superconducting transmon-based quantum computer based on simulated studies. Our approach is multi-stage. We first “bootstrap” a model based on simulation with access to the full state vector for measuring gate fidelity. We next present an algorithm, named adapted randomized benchmarking (ARB), for fine-tuning the gate on hardware based on measurements of the devices. We also present techniques for deploying the model on programmable devices with care to reduce the required resources. While the techniques here are applied to a transmon-based computer, many of them are portable to other architectures.

97 MATHEMATICS AND COMPUTING↗