SEARCH · Engineering Papers

Results for “Point Optimization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Controlling Ligand Excimer Formation with Dipole Changes in Emissive Rare-Earth/Phosphonic Acid Complexes

The interactions between substituted arylvinyl phosphonic acid (AVPA) ligands within a Eu-AVPA complex are shown to influence the outcomes of excited state evolution after photoexcitation. Compared with unfunctionalized AVPAs, pairs of ligands functionalized with CF 3 in the para position preassociate in the ground state of complexes with Eu 3+ according to calculated geometry optimizations. The CF 3 -substituted AVPA complexes show evidence of red-shifted optical absorption and undergo more efficient excimer formation, as revealed by transient absorption spectroscopy. We rationalize this behavior through simulations of excited-state geometry optimizations that reveal evolution toward interligand phenyl-phenyl planarity for specific excited states. Emission from complexed Eu 3+ after energy transfer from the ligand is found to be weaker with CF 3 substitution, which we hypothesize is due to intracomplex, interligand aggregates with excimer-promoting geometries. These observations point to the need to consider ground-state geometries as well as dynamic excited-state processes to understand the flow of energy in rare earth coordination complexes.

09 BIOMASS FUELS↗

DOE OSTI2023

The role of electron correlations in the electronic structure of putative Chern magnet TbMn6Sn6

Abstract A member of the RMn 6 Sn 6 rare-earth family materials, TbMn 6 Sn 6 , recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard U (DFT + U ) calculations to examine the electronic structure of TbMn 6 Sn 6 . To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration–interaction (CI) methods. We find that DFT + U and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT + U with Mn magnetic moments adjusted to the experiment predict the Dirac crossing in bulk to be close to the Fermi level, within ~120 meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Machine learning changes the rules for flux limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations has seen recent interest for reducing the computational load for fluid simulations. Here, we focus on determining a flux-limiter for shock capturing methods. Focusing on flux limiters provides a specific plug-and-play component for existing numerical methods. Since their introduction, an array of flux limiters has been designed. Using the coarse-grained Burgers' equation, we show that flux-limiters may be rank-ordered in terms of their log-error relative to high-resolution data. We then develop a theory to find an optimal flux-limiter and present flux-limiters that outperform others tested for integrating Burgers' equation on lattices with [Formula: see text], and 2x, 3x, 4x, and 8x coarse-grainings. We train a continuous piecewise linear limiter by minimizing the mean-squared misfit to six-grid point segments of high-resolution data, averaged over all segments. While flux limiters are generally designed to have an output of φ(r) = 1 at a flux ratio of r = 1, our limiters are not bound by this rule and yet produce a smaller error than standard limiters. Here we find that our machine learned limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we use our theory to learn flux-limiters that outperform standard limiters across a range of values (as opposed to at a specific fixed value) of coarse-graining, number of discretized bins, and diffusion parameter. This demonstrates the ability to produce flux limiters that should be more broadly useful than standard limiters for general applications.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Chemical control of polarization in thin strained films of a multiaxial ferroelectric: Phase diagrams and polarization rotation

The emergent behaviors in thin films of a multiaxial ferroelectric (FE) due to electrochemical coupling between the rotating polarization and surface ions are explored within the framework of the 2–4 Landau-Ginzburg-Devonshire (LGD) thermodynamic potential combined with the Stephenson-Highland (SH) approach. The combined LGD-SH approach allows us to describe the electrochemical switching and rotation of a polarization vector in a multiaxial ferroelectric film covered by surface ions with a charge density defined by the oxygen pressure. In this work, we calculate phase diagrams, analyze the dependence of polarization components on the applied voltage, and discuss the peculiarities of quasistatic ferroelectric, dielectric, and piezoelectric hysteresis loops in thin strained multiaxial ferroelectric films. The nonlinear surface screening by oxygen ions makes the diagrams very different from the known diagrams of, e.g., strained BaTiO 3 films. Quite unexpectedly, we predict the appearance of ferroelectric reentrant phases. The obtained results point to the possibility to control the appearance and features of ferroelectric, dielectric, and piezoelectric hysteresis in multiaxial FE films covered with surface ions by varying their concentration via the partial oxygen pressure. The LGD-SH description of a multiaxial FE film can be further implemented within the Bayesian optimization framework, paving the way toward predictive materials optimization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

A Multi-Function AAA Algorithm Applied to Frequency Dependent Line Modeling

Modeling of power distribution system components that are valid for a wide range of frequencies are crucial for highly accurate modeling of electromagnetic transient (EMT) events. This has recently become of interest due to the improvements needed for the resilient operation of distribution systems. Vector fitting (VF) is a very popular and commonly used algorithm for wide band representations of power system components in EMT simulations. In this research, we present a new multi-input rational approximation algorithm (MIAAA) and illustrate its advantages with respect to VF using examples of approximations of admittance matrices discussed in the literature. We show that MIAAA not only outperforms VF in terms of achieving better accuracy using lesser number of poles, but also has no numerical issues achieving convergence. In contrast to VF, MIAAA is not sensitive to the location of input sample points and it does not require good estimates for the location of the desired approximation poles. The novelty of this research work is the use of recent mathematical results to solve existing challenges in distribution system modeling and to develop rational approximations for power system models that intend to be optimal in terms of accuracy and performance.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2023

Savannah River National Laboratory – General Fusion 2023 INFUSE Report (Rev.1)

This report describes the results from an INFUSE research project, where Savannah River National Laboratory (SRNL) in collaboration with General Fusion (GF) used process modeling to understand and optimize commercial power plant (CPP) fuel cycle designs based on parameters provided by GF. The study primarily focused on two candidate fuel cycles with different blanket materials, one with a lead lithium eutectic (LLE) blanket and the other with a pure lithium (Li) blanket. LLE benefits from a low melting point, favorable neutronics, and lower reactivity, but liquid lithium has the potential for higher tritium breeding ratios (TBR) and does not poison the plasma as a high Z contaminant. It was found that the main differences between LLE and Li designs are the extraction technologies required to remove tritium from the blanket and the amount of tritium and its distribution within the facility. More than 80% of the in-process tritium inventory for the LLE design is contained in the isotope separation system, while for the Li design, over 60% of the in-process tritium inventory is contained within the blanket material. This is due to significant tritium retention by Li. For the Li blanket, the burden of tritium processing rests on the blanket extraction technology rather than the traditional exhaust processing route. Thus, the blanket extraction technology is a main driver of tritium inventory in the Li system and determines the subsequent interface with the tritium processing plant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A MIMO AAA Agorithm for Frequency Dependent Line Modeling

AAA algorithm↗

DOE OSTI2024

Stochastic Trust-Region Algorithm in Random Subspaces with Convergence and Expected Complexity Analyses

Here, this work proposes a framework for large-scale stochastic derivative-free optimization (DFO) by introducing STARS, a trust-region method based on iterative minimization in random subspaces. This framework is both an algorithmic and theoretical extension of a random subspace derivative-free optimization (RSDFO) framework, and an algorithm for stochastic optimization with random models (STORM). Moreover, like RSDFO, STARS achieves scalability by minimizing interpolation models that approximate the objective in low-dimensional affine subspaces, thus significantly reducing per-iteration costs in terms of function evaluations and yielding strong performance on largescale stochastic DFO problems. The user-determined dimension of these subspaces, when the latter are defined, for example, by the columns of so-called Johnson-Lindenstrauss transforms, turns out to be independent of the dimension of the problem. For convergence purposes, inspired by the analyses of RSDFO and STORM, both a particular quality of the subspace and the accuracies of random function estimates and models are required to hold with sufficiently high, but fixed, probabilities. Using martingale theory under the latter assumptions, an almost sure global convergence of STARS to a first-order stationary point is shown, and the expected number of iterations required to reach a desired first-order accuracy is proved to be similar to that of STORM and other stochastic DFO algorithms, up to constants.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

An asymptotically compatible approach for Neumann-type boundary condition on nonlocal problems

In this paper we consider 2D nonlocal diffusion models with a finite nonlocal horizon parameter δ characterizing the range of nonlocal interactions, and consider the treatment of Neumann-like boundary conditions that have proven challenging for discretizations of nonlocal models. We propose a new generalization of classical local Neumann conditions by converting the local flux to a correction term in the nonlocal model, which provides an estimate for the nonlocal interactions of each point with points outside the domain. While existing 2D nonlocal flux boundary conditions have been shown to exhibit at most first order convergence to the local counter part as δ → 0, the proposed Neumann-type boundary formulation recovers the local case as O(δ 2 ) in the L∞(Ω) norm, which is optimal considering the O(δ 2 ) convergence of the nonlocal equation to its local limit away from the boundary. We analyze the application of this new boundary treatment to the nonlocal diffusion problem, and present conditions under which the solution of the nonlocal boundary value problem converges to the solution of the corresponding local Neumann problem as the horizon is reduced. To demonstrate the applicability of this nonlocal flux boundary condition to more complicated scenarios, we extend the approach to less regular domains, numerically verifying that we preserve second-order convergence for non-convex domains with corners. Finally, based on the new formulation for nonlocal boundary condition, we develop an asymptotically compatible meshfree discretization, obtaining a solution to the nonlocal diffusion equation with mixed boundary conditions that converges with O(δ 2 ) convergence.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗

DOE OSTI2024

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Phase-based velocity extraction method for photonic Doppler velocimetry with potential higher time resolution

We present an extension of the [Takeda et al., J. Opt. Soc. Am. 72, 156 (1982)] phase extraction method to heterodyne photonic Doppler velocimetry applications. The method yields results equivalent to those obtained by the short-time Fourier transform (STFT), while offering potential improvements in time resolution. Unlike STFT, which relies on window functions, such as the Hamming window, that emphasize central data points and diminish the influence of edges, the extended Takeda method utilizes all data uniformly. This uniform treatment allows for the derivation of empirical equations that directly relate velocity error to the actual time resolution rather than to the local analysis duration. The established equation provides a useful metric for both optimizing hardware configuration and guiding data analysis. Simulation and experimental results confirm that, for a given dataset, specifying a target time resolution yields consistent velocity errors for both methods. These findings underscore the Takeda method’s advantages, particularly its potential higher time resolution and reduced computational burden, making it a valuable tool for high-throughput applications such as laser dynamic compression experiments.

Computer simulation↗

DOE OSTI2020

A Physical Model for Understanding the Activation of MoS 2 Basal-plane Sulfur Atoms for the Hydrogen Evolution Reaction

All the DFT calculations are done with the Vienna Ab Initio Simulation Package (VASP) using the projector augmented wave method. The Bayesian error estimation exchange-correlation functionals (BEEF) with van der Waals interactions are employed. A plane-wave cutoff energy of 400 eV is used together with PAW-PBE potentials where semi core p states are treated as valence. All the calculations allow for spin-polarization. The structures are relaxed until the force is converged to < 0.01 eV/Å. The lattice parameter of MoS 2 unit cell, optimized with this functional, is 3.19 Å. A (4×4) supercell is used to model all the transition metal doped MoS 2 systems studied here, including those with S vacancies. For calculations in the initial dopant structure exploration, the Brillouin zone is sampled with a 3x3x1 Monkhorst-Pack k-point mesh. A 6×6×1 Monkhorst-Pack k-point mesh is used for the H binding energy and density of states calculations. In all calculations, the vacuum layer is set as 15 Å to eliminate periodic interaction perpendicular to the basal plane.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Review of TEAMER Awards for WEC-Sim Support: Preprint

Testing Expertise and Access for Marine Energy Research (TEAMER) is a U.S. Department of Energy Water Power Technologies Office sponsored program, overseen by the Pacific Ocean Energy Trust, which aims to advance the state of marine energy technologies. The program connects technology developers with experts at U.S. facilities, including numerical modeling and analysis facilities, like WEC-Sim. The WEC-Sim facility is supported by the WEC-Sim development team at Sandia National Laboratories and the National Renewable Energy Laboratory. WEC-Sim (Wave Energy Converter SIMulator) is an open-source software for simulating wave energy converters. WEC-Sim can model the multi-body dynamics of devices comprised of bodies, joints, power take-off systems, and mooring systems. Since TEAMER's first round of support in 2020, there have been eighteen TEAMER awards focused on numerical model development in WEC-Sim. TEAMER awards for WEC-Sim support have modeled a wide range of wave energy converter archetypes, including point absorbers, attenuators, oscillating water columns, and many other novel architectures. A wide variety of studies have been conducted, leading to important insights for TEAMER partners and software improvements for WECSim. This article highlights several successful WEC-Sim TEAMER awards. The awards described herein include TEAMER recipients Ocean Motion Technologies, AquaHarmonics, iProTech, East Carolina University, Virginia Tech, Maiden Wave Energy, and the University of Massachusetts Dartmouth. The awards of these seven partners contain a wide range of investigations and cover the creation of baseline hydrodynamic models, PTO modeling, geometry optimization in both boundary element methods and WECSim, and model tuning and validation.

industry support↗

DOE OSTI2026

Using Filter Methods to Guide Convergence for ADMM, with Applications to Nonnegative Matrix Factorization Problems

Nonconvex, nonlinear optimization problems arise naturally in parameter fitting and machine learning. While augmented Lagrangian methods have demonstrated robust convergence for classes of these problems, their convergence for block updates has been relatively unexplored outside of the context of the alternating direction method of multipliers (ADMM). ADMM has seen extensive use in these applications, but may exhibit uncertain convergence behavior in many practical nonconvex settings, and struggles with general nonlinear constraints. In contrast, filter methods have proved effective in enforcing convergence for sequential quadratic programming methods and interior point methods with feasibility criteria. We develop an ADMM-filter method for highly nonlinear and nonconvex problems. Here, we show convergence under mild assumptions for several types of coordinate descent schemes, and demonstrate our algorithm on nonnegative matrix factorization and completion problems in imaging and chemical spectrum analysis.

Nonconvex optimization↗

DOE OSTI2021

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗

DOE OSTI2021

Anisotropic colloidal interactions & assembly in AC electric fields

In this work, we match experimental and simulated configurations of anisotropic epoxy colloidal particles in high frequency AC electric fields by identifying analytical potentials for dipole-field and dipole-dipole interactions. We report an inverse Monte Carlo simulation algorithm to determine optimal fits of analytical potentials by matching simulated and experimental distribution functions for non-uniform liquid, liquid crystal, and crystal microstructures in varying amplitude electric fields. Two potentials that include accurate particle volume and dimensions along with a concentration dependent prefactor quantitatively capture experimental observations. At low concentrations, an effective ellipsoidal point dipole potential works well, whereas a novel stretched point dipole potential is found to be suitable at all concentrations, field amplitudes, and degrees of ordering. The simplicity, accuracy, and adjustability of the stretched point dipole potential suggest it can be applied to model field mediated microstructures and assembly of systematically varying anisotropic particle shapes.

42 ENGINEERING↗

DOE OSTI2023

Dissolved gas recovery from water using a sidestream hollow-fiber membrane module: First principles model synthesis and steady-state validation

This paper presents a first-principles model for the recovery of dissolved gases from liquids using a sidestream hollow-fiber membrane module. The model avoids the use of new empirical coefficients, thus providing a parametric understanding of the process behavior for future design and optimization of membrane modules. This type of first-principles model could be particularly useful when gas recovery is beneficial to biological or chemical reactions of interest, such as the acetogenesis reactions in two-stage anaerobic digesters. The steady-state behavior of the model was validated against both new experimental data for the recovery of H 2 , CH 4 and H 2 –CH 4 mixtures from pure water, as well as existing published data. The modeled gas recovery predictions agreed with experimental data to an absolute average error of 13%, and an average R value of 0.98. Parametric analysis of mixed-gas recovery suggests possible key transition points in the composition of the recovered gases. For example, at 40 °C, increasing trans-membrane pressure while keeping hydraulic residence time (HRT) under 0.5 s will result in an increase in the ratio of H 2 to CH 4 recovered. Otherwise, increasing trans-membrane pressure will instead decrease the ratio of H 2 to CH 4 recovered. The model has potential to be extended to transient analysis, but has yet to be validated with transient experimental data. Furthermore, this model was successfully implemented in both Python and MATLAB, and provides valuable insights for future net-energy optimization for anaerobic digestion systems with in-situ gas recovery.

Anaerobic Digestion↗