Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Elaborated thermochemical treatment of HF, CO, N 2 , and H 2 O: Insight into HEAT and its extensions

Here, empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N 2 , and H 2 O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm -1 (0.12 kJ mol -1 ), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol—where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure—this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

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

Solving the electronic structure problem for over 100000 atoms in real space

Using a real-space high-order finite-difference approach, we investigate the electronic structure of large spherical silicon nanoclusters. Within Kohn-Sham density functional theory and using pseudopotentials, we report the self-consistent field convergence of a system with over 100000 atoms: a Si 107,641 ⁢H 9,084 nanocluster with a diameter of 16 nm. Our approach uses Chebyshev-filtered subspace iteration to speed up the convergence of the eigenspace, and blockwise Hilbert space-filling curves to speed up sparse matrix-vector multiplications, all of which are implemented in the parsec code. For the largest system, we utilized 2048 nodes (114 688 cores) on the Frontera machine in the Texas Advanced Computing Center. Our quantitative analysis of the electronic structure shows how it gradually approaches its bulk counterpart as a function of nanocluster size. The band gap is enlarged due to quantum confinement in nanoclusters, but decreases as the system size increases, as expected. In conclusion, our work serves as a proof of concept for the capacity of the real-space approach in efficiently parallelizing very large calculations using high-performance computer platforms, which can straightforwardly be replicated in other systems with more than 10 5 atoms.

0-dimensional systems↗
DOE OSTI2020

Towards Robustness Guarantees for Feedback-Based Optimization

Feedback-based online optimization algorithms have gained traction in recent years because of their simple implementation, their ability to reject disturbances in real time, and their increased robustness to model mismatch. While the robustness properties have been observed both in simulation and experimental results, the theoretical analysis in the literature is mostly limited to nominal conditions. In this work, we propose a framework to systematically assess the robust stability of feedback-based online optimization algorithms. We leverage tools from monotone operator theory, variational inequalities and classical robust control to obtain tractable numerical tests that guarantee robust convergence properties of online algorithms in feedback with a physical system, even in the presence of disturbances and model uncertainty. The results are illustrated via an academic example and a case study of a power distribution system.

approximation algorithms↗
DOE OSTI2021

IEA Wind Energy Task 37 - System Engineering - Aerodynamic Optimization Case Study

This paper presents the results from the first IEA Wind Task 37 aerodynamic optimization case study. 8 participants applied their optimization tools to a purely aerodynamic problem and the results were compared. Overall, the different tools produced widely different designs, while there was better agreement in the improvement achieved. This highlights the fact that further investigation is needed to try to understand these differences and develop best practices. There were too many differences between the different analysis and optimization codes to determine the sources of these differences. However, several potential sources of discrepancy were identified for further investigation. One hypothesis is that one potential source of discrepancy is that the design problem itself is relatively flat in design directions of constant loading. This flatness would impact the convergence in the optimization, while at the same time mean that discrepancies in the design itself would be less severe.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2023

LensWatch. I. Resolved HST Observations and Constraints on the Strongly Lensed Type Ia Supernova 2022qmx (“SN Zwicky”)

Supernovae (SNe) that have been multiply imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply imaged Type Ia SN 2022qmx (aka, "SN Zwicky") at z = 0.3544. SN Zwicky was discovered by the Zwicky Transient Facility in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multicycle "LensWatch (www.lenswatch.org)" program. We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, and F814W) but unresolved with WFC3/IR F160W, and present an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between lens-model-predicted time delays (≲1 day), and delays estimated with the single epoch of Hubble Space Telescope colors (≲3.5 days), including the uncertainty from chromatic microlensing (~1–1.5 days). Our lens models converge to an Einstein radius of ${\theta }_{{\rm{E}}}=({0.168}_{-0.005}^{+0.009})^{\prime\prime} $, the smallest yet seen in a lensed SN system. The "standard candle" nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter than predicted by $\sim {1.7}_{-0.6}^{+0.8}$ mag and $\sim {0.9}_{-0.6}^{+0.8}$ mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

IEA Wind Energy Task 37 System Engineering Aerodynamic Optimization Case Study: Preprint

This paper presents the results from the first IEA Wind Task 37 aerodynamic optimization case study. Eight participants applied their optimization tools to a purely aerodynamic problem and the results were compared. Overall, the different tools produced widely different designs, while there was better agreement in the improvement achieved. This highlights the fact that further investigation is needed to try to understand these differences and develop best practices. There were too many differences between the different analysis and optimization codes to determine the sources of these differences. However, several potential sources of discrepancy were identified for further investigation. One hypothesis is that one potential source of discrepancy is that the design problem itself is relatively flat in design directions of constant loading. This flatness would impact the convergence in the optimization, and at the same time mean that discrepancies in the design itself would be less severe.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2021

An embedded variable step IMEX scheme for the incompressible Navier–Stokes equations

This report presents a series of implicit–explicit (IMEX) variable stepsize algorithms for the incompressible Navier–Stokes equations (NSE). Here, IMEX means the nonlinear term is treated fully explicitly while the remaining terms are treated implicitly. With the advent of new computer architectures there has been growing demand for low memory solvers of this type. The addition of time adaptivity improves the accuracy and greatly enhances the efficiency of the algorithm. We prove energy stability of an embedded first–second order IMEX pair. For the first order member of the pair, we prove stability for variable stepsizes, and analyze convergence. We believe this to be the first proof of this type for a variable stepsize IMEX scheme for the incompressible NSE. We then define and test a variable stepsize, variable order IMEX scheme using these methods. Our work contributes several firsts for IMEX NSE schemes, including an energy argument and error analysis of a two-step, variable stepsize method, and embedded error estimation for an IMEX multi-step method.

42 ENGINEERING↗
DOE OSTI2025

Stochastic Adaptive Droop Control in Frequency Regulation of Power Systems With Intermittent Generators

Modern power systems (MPSs), including microgrids (MGs), are increasingly incorporating multiple renewable energy sources (RESs) such as wind and solar power, as well as battery storage and controllable loads. While environmentally beneficial, these sources pose challenges for control and management due to their intermittent and stochastic nature, especially in maintaining frequency stability with multiple interconnected generators of varying capacities. Traditional droop control methods are effective in systems with generators that are dispatchable and have fixed generation capacities, but they fall short when applied to systems with RESs, where generation capacities are dynamic and affected by unpredictable environmental conditions. To address these challenges, this paper introduces a novel stochastic adaptive droop control (SADC) method for load frequency control (LFC). The proposed method adapts droop coefficients in real time, based on the measured stochastic data of power generation capacities, enabling more effective frequency regulation in systems with variable and intermittent power generation. Unlike traditional adaptive control methods, which assume constant or slowly-varying system parameters, this approach accounts for stochastic processes by modeling them as Markov chains, enabling robust performance under highly dynamic and unpredictable conditions. The key contributions of this work include the development of real-time droop coefficient adaptation algorithms, derivation of their stability and convergence properties, and the demonstration of the advantages of the method through simulations. Case studies highlight the improved performance of frequency regulation, particularly in addressing the impact of stochastic weather conditions and the benefits of reducing dependence on battery reserves in dealing with intermittency of RESs. Finally, this paper provides a comprehensive analysis of the theoretical foundations of the method, as well as practical implementation insights for future power systems with high penetration of RESs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Effects of threshold resummation for large-x PDF in large momentum effective theory

Parton distribution functions (PDFs) at large x are challenging to extract from experimental data, yet they are essential for understanding hadron structure and searching for new physics beyond the Standard Model. Within the framework of the large momentum P z expansion of lattice quasi-PDFs, we investigate large x PDFs, where the matching coefficient is factorized into the hard kernel, related to the active quark momentum xP z , and the threshold soft function, associated with the spectator momentum (1 − x)P z . The renormalization group equation of the soft function enables the resummation of the threshold double logarithms α k ln 2k (1 − x), which is crucial for a reliable and controllable calculation of large x PDFs. Our analysis with pion valence PDFs indicates that perturbative matching breaks down when the spectator momentum (1 − x)P z approaches Λ QCD , but remains valid when both xP z and (1 − x)P z are much larger than Λ QCD . Additionally, we incorporate leading renormalon resummation within the threshold framework, demonstrating good perturbative convergence in the region where both spectator and active quark momenta are perturbative scales.

factorization↗
DOE OSTI2026

Coexisting kagome and heavy fermion flat bands in YbCr 6 Ge 6

Flat bands, electronic states with nearly dispersionless energy-momentum structure, provide fertile ground for unconventional quantum phases. Recent observations of flat bands at the Fermi level in kagome metals open the possibility of unifying topology and correlation-driven heavy-fermion physics. Here we show that topology and heavy-fermion correlations coexist in the layered kagome metal YbCr 6 Ge 6 . At high temperatures, an intrinsic kagome flat band—arising from frustrated hopping on the kagome lattice—dominates the Fermi level. Upon cooling, localized Yb 4f-states hybridize with the topological kagome flat bands, transforming this state into momentum-independent Kondo resonance states across the entire Brillouin zone. Topological analysis of the hybridization gaps reveals filling-tunable weak and strong topological Kondo-insulating regimes, and identifies a topological Dirac–Kondo semimetal. Taken together, these results identify YbCr 6 Ge 6 as a prototype of a topological heavy-fermion system and a platform where geometric frustration, strong correlations, and topology converge, with broad implications for correlated quantum matter.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

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

Dark Energy Survey Year 3 results: Imprints of cosmic voids and superclusters in the Planck CMB lensing map

ABSTRACT The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set (∼5000 deg2), expanding on previous measurements that used Y1 catalogues (∼1300 deg2). Given the increased statistical power compared to Y1 data, we report a 6.6σ detection of negative CMB convergence (κ) imprints using approximately 3600 voids detected from a redMaGiC luminous red galaxy sample. However, the measured signal is lower than expected from the MICE N-body simulation that is based on the ΛCDM model (parameters Ωm = 0.25, σ8 = 0.8), and the discrepancy is associated mostly with the void centre region. Considering the full void lensing profile, we fit an amplitude $A_{\kappa }=\kappa _{{\rm DES}}/\kappa _{{\rm MICE}}$ to a simulation-based template with fixed shape and found a moderate 2σ deviation in the signal with Aκ ≈ 0.79 ± 0.12. We also examined the WebSky simulation that is based on a Planck 2018 ΛCDM cosmology, but the results were even less consistent given the slightly higher matter density fluctuations than in MICE. We then identified superclusters in the DES and the MICE catalogues, and detected their imprints at the 8.4σ level; again with a lower-than-expected Aκ = 0.84 ± 0.10 amplitude. The combination of voids and superclusters yields a 10.3σ detection with an Aκ = 0.82 ± 0.08 constraint on the CMB lensing amplitude, thus the overall signal is 2.3σ weaker than expected from MICE.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Phenomenological Model for Cavitation

A phenomenological model of cavitation is presented, based on observations that both large relative negative pressures and large negative time derivatives of pressure are required for cavitation onset. We simulated two cavitation experiments to generate cavitation scaling parameters for relative pressure drop and rate of pressure drop. Our results show the model, while simple, is effective at reproducing results from laboratory experiments of cavitation. The parameters were then used in conjunction with a human surrogate computational model to predict, at any position within the head, the probability of intracranial cavitation caused by exposure to a blast event. The results suggest that the magnitude of blast overpressure observed in field data is sufficient to cause intracranial cavitation. Our analysis indicates that the helmeted head, when compared to the unhelmeted head configuration, results in a decrease but not elimination of cavitation exposure. When density functions of cavitation probability versus cumulative brain volume are combined with an injury severity model, the results show helmet efficacy at low and moderate risk levels. However, the convergence of unhelmeted and helmeted probability density functions at high-to-excessive risk thresholds indicates the helmet offers diminishing protection at elevated exposure levels, relative to the unhelmeted baseline. Future investigation and collaboration with neuroscience subject matter experts are needed to contextualize the current computational results. While the present work contributes specific and quantified predictions of intracranial cavitation location and severity, more research is required to apply our results to clinical settings with population-based brain injury subjects and controls. The relationship between our intracranial cavitation predictions with their anticipated clinical sequelae remains a topic in need of exploration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Blockchain for Fault-Tolerant Grid Operations

Radial topology and vast geographic coverage make distribution systems prone to widespread power outages upon the failure of a single (or multiple) upstream component. Fault-handling algorithms depend heavily on correct estimations of the system’s state to effectively isolate the affected area and reduce the number of affected customers while maintaining operational safety. The work described here leverages the core features of distributed, consensus-based decision-making processes and the immutability of blockchain, and demonstrates their value in improving fault-tolerant grid operations. In this work, blockchain was used to create a trusted data-sharing platform that enables independent actors to reconstruct the system state; this enables distributed resources to make intelligent decisions with limited knowledge. Although the process requires data sharing, its algorithms have been designed to limit the amount of private information that is exchanged, which helps preserve business-sensitive data and maintain customer privacy. In addition, by reducing the information that must be shared, the communication requirements are also reduced; (however, an in-depth analysis of the communication requirements is beyond the scope of this project). The proposed use cases are intended to represent a foundational basis for third parties to develop functional solutions that can eventually be deployed in the field. To further provide guidance, the envisioned use cases have incorporated design requirements that consider the blockchain characteristics and a need to limit information from surrounding resources, which preserve the assumption and the possibility that such resources could belong to different entities. This report presents a detailed design of the three use cases with the tools needed to enable the analysis being tested. The implemented gross error detection method can detect mismatches when the error exceeds 3.8 times the sensor’s rated accuracy. Detection of the circuit breaker state successfully identified the correct states across all simulation tests. A distribution-system power-flow solution in the simulator OpenDSS generally possesses a convergency tolerance of 0.01% on the voltage magnitude. The evaluation of possible reconnection using voltage magnitude—preserving the data ownership—has a voltage magnitude difference smaller than 0.001% from the OpenDSS result. The results preserving data ownership have a difference within the expected power flow tolerance with full knowledge of the system, which surpasses expectations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Optimizing processing conditions for additively reinforced thermoforming (ART) in convergent manufacturing

This study utilized additively reinforced thermoforming (ART) to enhance the thermomechanical properties of polyethylene terephthalate glycol (PETG) sheet. ART materials were produced by overprinting PETG/carbon fiber filament (PETG/CF) on neat PETG sheets at varying conditions. The mechanical properties of the PETG sheet, PETG/CF, and ART materials were assessed, showing that ART exhibited superior tensile strength and modulus of elasticity. The tensile strength and modulus in the x-direction for ART at 265°C were 57.32 ± 2.9 MPa and 3.41 ± 0.4 GPa, respectively, compared to 49.1 ± 0.5 MPa and 1.92 ± 0.09 GPa for neat PETG. Microstructural analysis revealed strong interfacial adhesion between layers, while thermogravimetric analysis (TGA), differential scanning calorimetry (DSC), and heat deflection temperature analysis provided insights into the ART material's thermoforming behavior, aiding design optimization for enhanced stiffness, reduced necking, and improved customization. In conclusion, this information can be used to design for the thermoforming operation.

Additive reinforcement↗
DOE OSTI2025

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2026

Nanosecond Structure of Radical Pair Intermediates from High-Frequency Quantum Oscillations: Insight into the Q A •– to Q B Electron Transfer Step in Purple Bacterial Photosynthesis

We demonstrate the validity of our approach to deduce, from the anisotropy of quantum oscillations, the geometry of short-lived radical pair intermediates in photosynthesis. A global fit of a two-dimensional W-band (94 GHz) electron paramagnetic resonance (EPR) experiment provides the same global minimum values for the geometry of the A-side radical pair P 700 •+ A 1A •− in photosystem I (PSI) as observed in a previous Q-band (34 GHz) EPR study, yet with a significantly increased convergence rate of 62%. This demonstrates that the global fit yields the correct radical pair geometry even at Q-band frequencies. With this information, we revisit our previous Q-band study of the cofactor arrangement of P 865 •+ Q A •− , the stabilized charge-separated state in purple bacterial reaction centers (RCs). Analysis of calculated two-dimensional data sets of P 865 •+ Q A •− reveals that the quantum oscillation technique is unaffected by a mirror ambiguity in disordered solids and thus can provide unambiguous solutions for all five Euler angles of the radical pair geometry. This enables us to elucidate the Q A •− to Q B electron transfer step in purple bacterial photosynthesis, the subject of controversial discussions for more than 25 years. Our results show that this electron transfer step involves a gating mechanism requiring a 60° rotation of the headgroup of Q A •− in its binding pocket.

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

Large-Angle Rocking Beam Electron Diffraction of Large Unit Cell Crystals Using Direct Electron Detector

We report a large-angle rocking beam electron diffraction (LARBED) technique for electron diffraction analysis. Diffraction patterns are recorded in a scanning transmission electron microscope (STEM) using a direct electron detector with large dynamical range and fast readout. We use a nanobeam for diffraction and perform the beam double rocking by synchronizing the detector with the STEM scan coils for the recording. Using this approach, large-angle convergent beam electron diffraction (LACBED) patterns of different reflections are obtained simultaneously. By using a nanobeam, instead of a focused beam, the LARBED technique can be applied to beam-sensitive crystals as well as crystals with large unit cells. Here, this paper describes the implementation of LARBED and evaluates the performance using silicon and gadolinium gallium garnet crystals as test samples. We demonstrate that our method provides an effective and robust way for recording LARBED patterns and paves the way for quantitative electron diffraction of large unit cell and beam-sensitive crystals.

4DSTEM↗