Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energy minimization”

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

A FETI approach to domain decomposition for meshfree discretizations of nonlocal problems

We propose a domain decomposition method for the efficient simulation of nonlocal problems. Our approach is based on a multi-domain formulation of a nonlocal diffusion problem where the subdomains share “nonlocal” interfaces of the size of the nonlocal horizon. This system of nonlocal equations is first rewritten in terms of minimization of a nonlocal energy, then discretized with a meshfree approximation and finally solved via a Lagrange multiplier approach in a way that resembles the finite element tearing and interconnect method. Specifically, we propose a distributed projected gradient algorithm for the solution of the Lagrange multiplier system, whose unknowns determine the nonlocal interface conditions between subdomains. Several two-dimensional numerical tests on problems as large as 191 million unknowns illustrate the strong and the weak scalability of our algorithm, which outperforms the standard approach to the distributed numerical solution of the problem. Finally, this work is the first rigorous numerical study in a two-dimensional multi-domain setting for nonlocal operators with finite horizon and, as such, it is a fundamental step towards increasing the use of nonlocal models in large scale simulations.

42 ENGINEERING↗
DOE OSTI2022

Q-POP-Thermo: A general-purpose thermodynamics solver for ferroelectric materials

We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

On the convexity of phase-field fracture formulations: Analytical study and comparison of various degradation functions

Efficient and accurate fracture modeling is of great importance in applications where catastrophic outcomes under extreme scenarios are possible. The phase-field (PF) approach to fracture received significant attention over the past decade, due to its capability to capture complicated fracture patterns (e.g., crack merging and branching). Specifically, crack initiation and propagation are modeled via minimization of the total energy functional, which is regularized with the aid of a phase field. Despite the promising results and modeling capabilities of the PF method in many applications, the solution of fracture problems remains computationally challenging mainly due to the non-convexity of the total energy functional with respect to the combined unknown (phase field and displacement) fields. Understanding the effects of their coupling on convexity is crucial in order to address frequently encountered hurdles in fracture modeling (e.g., inefficient solvers and non-physical crack nucleation). In this paper, we develop convexity criteria for a wide class of PF fracture formulations. For this class of formulations, the second variation of the total energy functional is expressed in terms of Hessian matrices (evaluated at individual material points). Depending on the choice of geometric crack functions and degradation functions, we classify the formulations into three categories and analytically study each one separately. To study the sign of the second variation, we derive inequalities which are satisfied at material points when the Hessian matrix is locally positive semi-definite. These inequalities provide objective criteria for comparing degradation functions. Finally, the applicability of the proposed convexity criteria is demonstrated in the context of a one-dimensional problem, solved using a conventional monolithic solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Effect of static local distortions vs. dynamic motions on the stability and band gaps of cubic oxide and halide perovskites

Ternary ABX 3 perovskites made of corner-sharing BX 6 octahedra have long featured prominently in solid-state chemistry and condensed matter physics. Still, the joint understanding of their two main subgroups—halides and oxides—has not been fully developed. Indeed, unlike the case in simpler compounds having a single, robust repeated motif (“monomorphous”), certain cubic perovskites can manifest a non-thermal (= intrinsic) distribution of local motifs (“polymorphous networks”). Such static deformations can include positional degrees of freedom (e.g., atomic displacements and octahedral tilting) or magnetic moment degrees of freedom in paramagnets. Unlike thermal motion, such static distortions do not time-average to zero, being an expression of the intrinsic symmetry breaking preference of the chemical bonding. Here, the present study compares electronic structure features of oxide and halide perovskites starting from the static polymorphous distribution of motifs described by Density Functional Theory (DFT) minimization of the internal energy, continuing to finite temperature thermal disorder modeled via finite temperature DFT molecular dynamics. We find that (i) different oxide vs. halide ABX 3 compounds adopt different energy-lowering symmetry-breaking modes. The calculated pair distribution function (PDF) of SrTiO 3 from the first-principles agrees with recently measured PDF. (ii) In both oxides and halides, such static distortions lead to band gap blueshifts with respect to undistorted cubic Pm-3m structure. (iii) For oxide perovskites, high-temperature molecular dynamics simulations initiated from the statically distorted polymorphous structures reveal that the thermally-induced distortions can lead to a band gap redshift. (iv) In contrast, for cubic halide perovskite CsPbI 3 , both the intrinsic distortions and the thermal distortions contribute in tandem to band gap blueshift, the former, intrinsic effect being dominant. (v) In the oxide SrTiO 3 and CaTiO 3 (but not in halide) perovskites, octahedral tilting leads to the emergence of a distinct Γ–Γ direct band gap component as a secondary valley minimum to the well-known indirect R–Γ gap. Understanding such intrinsic vs. thermal effects on oxide vs. halide perovskites holds the potential for designing target electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Shape Control of Colloidal Semiconductor Nanocrystals through Thermodynamically Driven Aggregative Growth

The optoelectronic properties of colloidal semiconductor nanocrystals (NCs) can be manipulated by changing their geometric shapes. The precise synthetic control over particle morphologies, however, has remained elusive. Conventional growth techniques rely on the kinetic assembly of atomic units, where supersaturation and precipitation processes can lead to a broad distribution of particle shapes. In this paper, we demonstrate that replacing atomic precursors with small-size nanocrystals as building blocks for larger colloids offers an easier, more predictive control over nanoparticle shape evolution. The reported growth strategy is illustrated via shape-selective syntheses of CdSe and CdS NC cubes, spheres, rods, as well as unprecedented “donut” and ring-like structures. Growth of desired particle morphologies was achieved by choosing an appropriate reaction solvent that minimizes the surface free energy through a unique shape formation. This so-called "aggregative growth mechanism" is explained using a thermodynamic model for interacting viscous colloids. We expect that this novel growth strategy will offer a viable platform for shape–selective synthesis of many inorganic colloids, providing pathways to optimizable optoelectronic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Improving the Accuracy of Variational Quantum Eigensolvers with Fewer Qubits Using Orbital Optimization

Near-term quantum computers will be limited in the number of qubits on which they can process information as well as the depth of the circuits that they can coherently carry out. To date, experimental demonstrations of algorithms such as the Variational Quantum Eigensolver (VQE) have been limited to small molecules using minimal basis sets for this reason. In this work we propose incorporating an orbital optimization scheme into quantum eigensolvers wherein a parametrized partial unitary transformation is applied to the basis functions set in order to reduce the number of qubits required for a given problem. The optimal transformation is found by minimizing the ground state energy with respect to this partial unitary matrix. Through numerical simulations of small molecules up to 16 spin orbitals, we demonstrate that this method has the ability to greatly extend the capabilities of near-term quantum computers with regard to the electronic structure problem. Finally, we find that VQE paired with orbital optimization consistently achieves lower ground state energies than traditional VQE when using the same number of qubits and even frequently achieves lower ground state energies than VQE methods using more qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Thermal Weight Determination and Interstate Coupling in State-Averaged ADAPT-VQE

Characterizing electronic thermal states at low temperatures is an important but challenging task in quantum chemistry and condensed matter physics, making it a prime candidate for a useful application in quantum computing. One of the most successful methods for state preparation on quantum computers is the Adaptive, Problem-Tailored (ADAPT) Variational Quantum Eigensolver (VQE), which has recently been generalized to treat excited states within a state-averaged framework as well as Gibbs states. In this work, we introduce Helmholtz-Optimized Thermal (HOT) ADAPT-VQE, an ancilla-free strategy for preparing Gibbs states that directly minimizes the Helmholtz free energy by targeting the dominant eigenstates of the thermal ensemble. We demonstrate the usefulness of HOT-ADAPT-VQE by predicting the free energy of two model systems with strongly correlated ground states: (1) the Fe 2+ cation in a magnetic field and (2) a [Cu 2 O 7 ] 10– fragment of the Mott insulator La 2 CuO 4 . Our results demonstrate that HOT-ADAPT-VQE significantly improves upon Gibbs-state estimates from multistate variants of ADAPT-VQE, often with substantially shallower quantum circuits, making it a promising candidate for thermal-state calculations.

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

How Well Does Kohn–Sham Regularizer Work for Weakly Correlated Systems?

Kohn–Sham regularizer (KSR) is a differentiable machine learning approach to finding the exchange-correlation functional in Kohn–Sham density functional theory that works for strongly correlated systems. Here we test KSR for a weak correlation. We propose spin-adapted KSR (sKSR) with trainable local, semilocal, and nonlocal approximations found by minimizing density and total energy loss. We assess the atoms-to-molecules generalizability by training on one-dimensional (1D) H, He, Li, Be, and Be 2+ and testing on 1D hydrogen chains, LiH, BeH 2 , and helium hydride complexes. The generalization error from our semilocal approximation is comparable to other differentiable approaches, but our nonlocal functional outperforms any existing machine learning functionals, predicting ground-state energies of test systems with a mean absolute error of 2.7 mH.

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

Evaluating the Structural Response of Amphiphilic Monolayers to Environmental Stimuli

Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.

functional groups↗
DOE OSTI2023

Plasma Properties in the Earth's Magnetosheath Near the Subsolar Magnetopause: Implications for Geocoronal Density Estimates

Combined in situ ion measurements and remote sensing of energetic neutral atoms are used to determine the geocoronal Hydrogen density at large (∼10 R E ) distances from the Earth. This method for determining the geocoronal density requires global magnetospheric modeling. Observations in the Earth's subsolar magnetosheath from the Magnetospheric Multiscale mission are used to determine the accuracy of using global models to predict the geocoronal density. On average, gas dynamic and magnetohydrodynamic (MHD) models and observations are in reasonable agreement, with differences <25%. In addition, the MHD model subsolar magnetopause is about 0.5 R E sunward of the observed location. However, variations around averages are large (up to a factor of 2), indicating that global models introduce relatively large uncertainties in geocoronal density estimates. Finally, the critical ion flux in the Interstellar Boundary Explorer IBEX‐Hi energy range is often minimally affected by fluctuations of a factor of 2 in the density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Carbon monoxide oxidation expands the known metabolic capacity in anaerobic methanotrophic consortia

Consortia of anaerobic methane-oxidizing archaea (ANME-2) and sulphate-reducing bacteria (SRB) represent globally relevant syntrophic associations capable of growing with minimal amounts of free energy and can persist when methane becomes limiting. Carbon monoxide (CO) has been reported in seep environments and represents a thermodynamically favourable alternative electron donor due to its low reduction potential. Here, we show that environmental ANME-SRB consortia can oxidize CO in the absence of methane, in anoxic microcosm experiments using a combination of stable isotope geochemical tracers, metatranscriptomics, and single cell activity measurements (FISH–nanoSIMS). The oxidation of CO was coupled with sulphate-reduction by syntrophic consortia, and, in the absence of sulphate, through CO 2 reduction to methane by ANME-2. Under these conditions, the production of methane was one ninth the rate of methanotrophy coupled to sulphate-reduction. Paired single cell FISH-nanoSIMS analysis of anabolic activity indicates that CO respiration appears to support cell maintenance rather than active growth, consistent with the observed down-regulation of energy generating complexes in ANME (e.g., mtr, rnf, etc.). The versatile capability of CO oxidation by anaerobic methanotrophic consortia broadens our understanding of carbon cycling in methane seeps and highlights potential mechanisms of resilience by methanotrophic archaea under changing geochemical regimes.

03 NATURAL GAS↗
DOE OSTI2023

A j eff = 1/2 Kitaev material on the triangular lattice: the case of NaRuO 2

Motivated by recent reports of a quantum-disordered ground state in the triangular lattice compound NaRuO 2 , we derive a j eff = 1/2 magnetic model for this system by means of first-principles calculations. The pseudospin Hamiltonian is dominated by bond-dependent off-diagonal Γ interactions, complemented by a ferromagnetic Heisenberg exchange and a notably antiferromagnetic Kitaev term. In addition to bilinear interactions, we find a sizable four-spin ring exchange contribution with a strongly anisotropic character, which has been so far overlooked when modeling Kitaev materials. The analysis of the magnetic model, based on the minimization of the classical energy and exact diagonalization of the quantum Hamiltonian, points toward the existence of a rather robust easy-plane ferromagnetic order, which cannot be easily destabilized by physically relevant perturbations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Downward quantum learning from element 118: Automated generation of Fermi–Löwdin orbitals for all atoms

A new algorithm based on a rigorous theorem and quantum data computationally mined from element 118 guarantees automated construction of initial Fermi–Löwdin-Orbital (FLO) starting points for all elements in the Periodic Table. It defines a means for constructing a small library of scalable FLOs for universal use in molecular and solid-state calculations. The method can be systematically improved for greater efficiency and for applications to excited states such as x-ray excitations and optically silent excitations. FLOs were introduced to recast the Perdew–Zunger self-interaction correction (PZSIC) into an explicit unitarily invariant form. The FLOs are generated from a set of N quasi-classical electron positions, referred to as Fermi-Orbital descriptors (FODs), and a set of N-orthonormal single-electron orbitals. FOD positions, when optimized, minimize the PZSIC total energy. However, creating sets of starting FODs that lead to a positive definite Fermi orbital overlap matrix has proven to be challenging for systems composed of open-shell atoms and ions. The proof herein guarantees the existence of a FLOSIC solution and further guarantees that if a solution for N electrons is found, it can be used to generate a minimum of N – 1 and a maximum of 2 N – 2 initial starting points for systems composed of a smaller number of electrons. As a result, applications to heavy and super-heavy atoms are presented. All starting solutions reported here were obtained from a solution for element 118, Oganesson.

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

Origin of medium-range atomic correlation in simple liquids: Density wave theory

The atomic pair-distribution function of simple liquid and glass shows exponentially decaying oscillations beyond the first peak, representing the medium-range order (MRO). The structural coherence length that characterizes the exponential decay increases with decreasing temperature and freezes at the glass transition. Conventionally, the structure of liquid and glass is elucidated by focusing on a center atom and its neighboring atom shell characterized by the short-range order (SRO) and describing the global structure in terms of overlapping local clusters of atoms as building units. However, this local bottom-up approach fails to explain the strong drive to form the MRO, which is different in nature from the SRO. We propose to add an alternative top-down approach based upon the density wave theory. In this approach, one starts with a high-density gas state and seeks to minimize the global potential energy in reciprocal space through density waves using the pseudopotential. The local bottom-up and global top-down driving forces are not mutually compatible, and the competition and compromise between them result in a final structure with the MRO. This even-handed approach provides a more intuitive explanation of the structure of simple liquid and glass.

42 ENGINEERING↗
DOE OSTI2022

On-policy learning-based deep reinforcement learning assessment for building control efficiency and stability

Artificial intelligence technologies have emerged as a game changer not only in specific applications such as image recognition and machine translation but also in many scientific domains. In particular, as deep reinforcement learning (DRL) has shown great success in complex control problems, DRL-based control has been considered as a potential solution to efficiently control and manage building systems. However, broad assessment of DRL-based building control is still required to characterize their pros and cons in comparison with conventional building control methods (e.g., rule-based feedback controls). In this paper, we assessed DRL-based controls with on-policy learning-based algorithms and continuous control actions for cooling control of large office buildings in the summer season to minimize whole-building energy use and occupant discomfort. We compared DRL-based control methods with two baseline control methods: (1) a pre-determined schedule with supply temperature and static pressure setpoints, and (2) advanced reset method that adjusts setpoints based on heuristic rules, i.e., ASHRAE Guideline 36. We also tested the DRL algorithms to evaluate their performances in multiple climate locations. We found that DRL-based control methods outperformed the baseline control methods in terms of energy savings while maintaining a thermal comfort. DRL reduced energy use between ~4%–22% on average compared to the baseline methods, depending on climate location. We also evaluated DRL-based control in terms of control stability and showed that DRL-based methods should address the span of hardware lifetimes in practical operations.

control stability↗
DOE OSTI2020

Equilibrium shapes of facetted 3D metal nanoclusters intercalated near the surface of layered materials

Experimental studies indicate that 3D crystalline metal nanoclusters (NCs) intercalated under the surface of graphite have flat-topped equilibrated shapes. In this work, we characterize the shapes of these facetted NCs sandwiched between a blanketing graphene layer and the underlying graphite substrate. Specifically, we focus on the cases of fcc Cu and hcp Fe NCs. Furthermore, the analysis involves numerical minimization of the system energy for a specified NC volume and NC height, the latter corresponding to the separation between parallel top and bottom facets. Our numerical analysis quantifies how the distance of the side facet planes from center of the nanocluster varies linearly with a natural characteristic linear dimension of the nanocluster. Calculated shapes of fcc Cu and hcp Fe NCs are consistent with the hexagonal footprints observed in scanning tunneling microscopy studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Flux trapping in NbTiN strips

We use scanning superconducting quantum interference device (SQUID) microscopy to image vortices in superconducting structures fabricated from NbTiN thin films. We repeatedly cool superconducting strips with different width in an applied magnetic field and image the individual vortices. From these images we determine the threshold field at which the first vortex enters a strip and the number and configuration of vortices beyond this threshold field. We model the behavior of the vortices with and without considering the effect of pinning using numerical methods to minimize the Gibbs free energy of vortices in the strips. Our measurements provide a first benchmark to investigate the flux trapping properties of NbTiN thin films directly relevant to NbTiN-based superconducting digital circuits.

Bai, Ruiheng [Cornell University] (ORCID:000000025↗
DOE OSTI2020

Automated Design of Neuromorphic Networks for Scientific Applications at the Edge

Designing spiking neural networks for neuromorphic deployment is a non-trivial task. It is further complicated when there are resource constraints for the neuromorphic implementation, such as size or power constraints, that may be present in edge applications. In this work, we utilize a previously presented approach, EONS, to design spiking neural networks for a memristive neuromorphic implementation for scientific data applications. We specifically use a multi-objective approach in EONS to maximize network accuracy on the scientific data application task, but also to minimize network size and energy. We illustrate that EONS determines both the network structure and the parameters, removing the burden from the user on determining the appropriate spiking neural network structure, and we show that the resulting networks are very different from the layered structure of typical neural networks. Finally, we show that the multi-objective approach produces smaller, more energy efficient networks than the original EONS approach and produces comparable accuracy to a back-propagation style training approach.

Schuman, Catherine↗