Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cost function”

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 271 records · Page 15

Mechanical properties of zeolite-templated carbons from approximate density functional theory calculations

Zeolite-templated carbon (ZTC) is a unique porous carbonaceous material whose structure is ordered at the nanometre scale, enabling a representative periodic description at the atomistic level. Utilizing an existing, well-defined reference model for ZTCs, a structural library of varying compositions was developed by refinement using density-functional tight-binding (DFTB) potentials parameterized for materials science applications. We first determined the quantum chemical-refined structures of models with CH, CHO, CHON, CHOB, and CHOBN compositions with various degrees of heteroatom substitution. These structural models comprise the characteristic morphological features of highly porous carbon materials, such as open-blade surfaces, edges, saddles, and closed-strut formations, spanning a range of curvatures and characteristic sizes. Second, we carried out alternating compression and expansion of the CHO model unit cell to determine the lowest energy structure as well as to obtain its bulk modulus in order to demonstrate a close connection between macroscopic observations and atomic-scale structures. Further, the agreement between experimental measurements and the computational model is remarkable and demonstrates the power of approximate density functional theory as a cost-effective computational tool with chemical accuracy for the investigation of structure/property relationships in real-world carbon-based solids.

03 NATURAL GAS↗

Design and dispatch optimization of packaged ice storage systems within a connected community

The traditional implementation of cool thermal energy storage (CTES) must be reimagined within the context of a dynamic grid and smart buildings operating as connected communities. As most buildings do not operate central chillers or connect to district cooling loops, this necessitates a broader use of packaged CTES. Our objective is to begin answering the question of how such packaged CTES should be implemented within a connected community. We do so by presenting a simulation–optimization workflow employing building energy modeling software and a mixed-integer linear program to design and dispatch a packaged CTES technology to achieve minimum total annual cost. We demonstrate this methodology on a seven-building case study using current utility rates and find that total annual cooling energy costs can be reduced by 17.8% compared to baseline, after accounting for the cost of storage. We perform three parametric sensitivity studies to evaluate modeling assumptions and obtain the prioritization of storage procurement as a function of annualized life-cycle cost of storage. We find that a community optimization approach provides significantly different results than individual building optimizations and provides greater savings compared to baseline.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗

A deformation-dependent coupled Lagrangian/semi-Lagrangian meshfree hydromechanical formulation for landslide modeling

Abstract The numerical modelling of natural disasters such as landslides presents several challenges for conventional mesh-based methods such as the finite element method (FEM) due to the presence of numerically challenging phenomena such as severe material deformation and fragmentation. In contrast, meshfree methods such as the reproducing kernel particle method (RKPM) possess unique features conducive to modelling extreme events such as the absence of a structured mesh and the ease of adaptive refinement, among others. While the semi-Lagrangian reproducing kernel (SL-RK) shape functions of RKPM defined in the current configuration have proven to be effective in extreme event modelling, the computational cost for the re-evaluation of the shape functions at every time step is costly. In this work, a deformation-dependent coupling of the Lagrangian reproducing kernel (L-RK) and SL-RK approximations is proposed for the solution of a hydro-mechanical formulation for effective simulations of landslides. The ramp function is constructed based on an equivalent plastic strain as a deformation-dependent transition from L-RK shape functions to SL-RK ones as the deformation progresses. The particular focus of the paper will be on modelling seepage-induced landslides with a mixed $$u$$ u – $$p$$ p formulation to couple the solid and fluid phases. Examples are presented to examine the effectiveness of this coupled Lagrangian/semi-Lagrangian reproducing kernel (L–SL RK) formulation and to highlight its performance in landslide modelling.

Baek, Jonghyuk↗

Efficient Quantum Gate Discovery with Optimal Control

Optimal control theory provides a framework for numerical discovery of device controls that implement quantum logic gates, but common objective functions used for optimization often assign arbitrarily high costs to otherwise useful controls. We propose a framework for designing objective functions that permit novel gate designs such as echo pulses or locally-equivalent gates. We use numerical simulations to demonstrate the efficacy of the new objective functions by designing microwave-only pulses that act as entangling gates for superconducting transmon architectures. We observe that the proposed objective functions lead to higher fidelity controls in fewer optimization iterations than obtainable by traditional objective functions.

Kairys, Paul↗

Cyclic olefin copolymer-based reinforced anion exchange membranes for water electrolyzers

Anion exchange membranes (AEMs) have emerged as a promising technology for water electrolysis in hydrogen production since they offer significant cost reduction in choices of electrocatalysts and bipolar plates. However, AEMs satisfying multiple requirements of high ionic conductivity, good chemical stability, robust mechanical properties, scalable synthesis, and low manufacturing costs are rare. Herein, we introduce quaternary ammonium functionalized cyclic olefin copolymers (COCs) as a new class of chemically stable and low-cost AEM materials. To further enhance the mechanical robustness, we prepared reinforced composite AEMs by impregnating the ionically functionalized COC into a mechanically robust matrix. The resulting reinforced composite membrane exhibits a high hydroxide conductivity of 127 mS cm −1 and excellent mechanical strength. In water electrolyzers, the MEA demonstrated outstanding performance, achieving a current density of 2.24 A cm −2 at 1.8 V, attributable to high conductivity, enhanced mechanical properties, and good alkaline stability of the composite membrane. These results indicate that the COC-based AEMs demonstrate good potential for application in AEM electrolyzers.

08 HYDROGEN↗

A deep potential model with long-range electrostatic interactions

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Structure prediction of epitaxial inorganic interfaces by lattice and surface matching with Ogre

We present a new version of the Ogre open source Python package with the capability to perform structure prediction of epitaxial inorganic interfaces by lattice and surface matching. In the lattice matching step, a scan over combinations of substrate and film Miller indices is performed to identify the domain-matched interfaces with the lowest mismatch. Subsequently, surface matching is conducted by Bayesian optimization to find the optimal interfacial distance and in-plane registry between the substrate and the film. For the objective function, a geometric score function is proposed based on the overlap and empty space between atomic spheres at the interface. The score function reproduces the results of density functional theory (DFT) at a fraction of the computational cost. The optimized interfaces are pre-ranked using a score function based on the similarity of the atomic environment at the interface to the bulk environment. Final ranking of the top candidate structures is performed with DFT. Ogre streamlines DFT calculations of interface energies and electronic properties by automating the construction of interface models. The application of Ogre is demonstrated for two interfaces of interest for quantum computing and spintronics, Al/InAs and Fe/InSb.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of Tailored Hydrocarbon-Based Pentablock Copolymer Membranes for Sodium-Polysulfide Flow Batteries

Long-duration energy storage (LDES) technologies are pivotal for the adoption of renewables like wind and solar. Non-aqueous redox flow batteries (NARFBs) with a sodium-polysulfide hybrid system feature high energy density independent of power density, yet face challenges with polysulfide shuttling. This study investigates a hydrocarbon-based penta-block copolymer membrane, Nexar, to mitigate crossover effects by balancing TFSI conversion and their crosslink density. The membranes are annealed to induce crosslinking for reducing electrolyte uptake and enhancing mechanical stability while demonstrating excellent ionic conductivity. The hydrocarbon-based membranes address environmental concerns associated with perfluoroalkyl substances and improve the performance and durability of NARFBs. In conclusion, our findings suggest that annealed Nexar membranes with tailored TFSI functionality offer a scalable, cost-effective solution for enhancing the efficiency of high-capacity energy storage systems, pivotal for grid integration of renewable sources.

25 ENERGY STORAGE↗

Apollo Next Generation Sample Analysis (ANGSA): an Apollo Participating Scientist Program to Prepare the Lunar Sample Community for Artemis

As a first step in preparing for the return of samples from the Moon by the Artemis Program, NASA initiated the Apollo Next Generation Sample Analysis Program (ANGSA). ANGSA was designed to function as a low-cost sample return mission and involved the curation and analysis of samples previously returned by the Apollo 17 mission that remained unopened or stored under unique conditions for 50 years. These samples include the lower portion of a double drive tube previously sealed on the lunar surface, the upper portion of that drive tube that had remained unopened, and a variety of Apollo 17 samples that had remained stored at -27 °C for approximately 50 years. ANGSA constitutes the first preliminary examination phase of a lunar “sample return mission” in over 50 years. It also mimics that same phase of an Artemis surface exploration mission, its design included placing samples within the context of local and regional geology through new orbital observations collected since Apollo and additional new “boots-on-the-ground” observations, data synthesis, and interpretations provided by Apollo 17 astronaut Harrison Schmitt. ANGSA used new curation techniques to prepare, document, and allocate these new lunar samples, developed new tools to open and extract gases from their containers, and applied new analytical instrumentation previously unavailable during the Apollo Program to reveal new information about these samples. Most of the 90 scientists, engineers, and curators involved in this mission were not alive during the Apollo Program, and it had been 30 years since the last Apollo core sample was processed in the Apollo curation facility at NASA JSC. There are many firsts associated with ANGSA that have direct relevance to Artemis. ANGSA is the first to open a core sample previously sealed on the surface of the Moon, the first to extract and analyze lunar gases collected in situ, the first to examine a core that penetrated a lunar landslide deposit, and the first to process pristine Apollo samples in a glovebox at -20 °C. All the ANGSA activities have helped to prepare the Artemis generation for what is to come. The timing of this program, the composition of the team, and the preservation of unopened Apollo samples facilitated this generational handoff from Apollo to Artemis that sets up Artemis and the lunar sample science community for additional successes.

79 ASTRONOMY AND ASTROPHYSICS↗

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Agent-Supervisor Coordination for Decentralized Event-Triggered Optimization

This letter proposes decentralized resource-aware coordination schemes for solving network optimization problems defined by objective functions that combine locally evaluable costs with network-wide coupling components. These methods are well suited for a group of supervised agents trying to solve an optimization problem under mild coordination requirements. Each agent has information on its local cost and coordinates with the network supervisor for information about the coupling term of the cost. The proposed approach is feedback-based and asynchronous by design, guarantees anytime feasibility, and ensures the asymptotic convergence of the network state to the desired optimizer. Numerical simulations on a power system example illustrate our results.

decentralized algorithms↗

Science Applications of Phased Array Radars

Phased array radars (PARs) are a promising observing technology, at the cusp of being available to the broader meteorological community. PARs offer near instantaneous sampling of the atmosphere with flexible beam forming, multi-functionality, low operational and maintenance costs, and without mechanical inertia limitations. These PAR features are transformative compared to those offered by our current reflector-based meteorological radars. The integration of PARs into meteorological research has the potential to revolutionize the way we observe the atmosphere. The rate of adoption of PARs in research will depend on many factors including i) the need to continue educating the scientific community on the full technical capabilities and trade-offs of PARs through an engaging dialogue with the science and engineering communities and ii) the need to communicate the breadth of scientific bottlenecks that PARs can overcome in atmospheric measurements and the new research avenues that are now possible using PARs in concert with other measurement systems. The former is the subject of a companion article that focuses on PAR technology while the latter is the objective here.

54 ENVIRONMENTAL SCIENCES↗

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Cooperative Research and Development Agreement between National Energy Technology Laboratory and VariGrid Explorations Inc. [Final Report]

The scope of this project is to advance the National Energy Technology Laboratory’s (NETL) Variable Grid Method (VGM) technology to commercial status thought enhancements upon prior work to produce a VGM application that can rapidly quantify multiple types of data and uncertainty commonly encountered by the subsurface energy industry. Along with development, the project will also include technology demonstrations, through applications using common industry data types and uncertainty metrics, to stress the value the VGM can offer the subsurface energy industry by providing an improved understanding of the strengths and limitations of their data, models, and analytical results. The resulting commercially ready VGM will offer interoperability with common industry software applications to ensure broad ‘plug and play’ functionality that will help keep costs relatively low and foster more rapid industry adoption.

97 MATHEMATICS AND COMPUTING↗

Recovery of Rare Earths, Precious Metals and other Critical Materials from Geothermal Waters with Advanced Sorbent Structures - CRADA 355 (Abstract)

The ability to recover valuable trace level minerals from geothermal brines using high-performance solid-phase sorbents will be explored and developed. A compressive range of sorbent materials will be screened for application to metal extraction from geothermal brines. Preferred sorbents from extraction of trace levels of rare earths (REs), precious metals (PMs), and other critical/strategically valuable materials (CMs) such as Zn, Mn, Te, Sc, Se and U from geothermal brines will be identified. For the preferred sorbents PNNL will determine the volumes they are capable of providing efficient extraction from. The thermal and chemical limits (including, acid, sulfur, salt) for performance of the preferred sorbent materials will be determined; with a target of at least 125°C and perhaps as high as 400° C. Sorbent form factors (including, packed bed, membrane, mats) that can function efficiently and be installed cost effectively in geothermal energy plants will be assessed for chemical and economic viability. Material regeneration and cyclic utilization will be demonstrated, targeting hundreds to thousands of cycles. Options for recovery and purification (including, selective separation of heavy REs) of collected materials will be explored. A techno-economic analysis (TEA) will be performed to assess the best approach to provide a value-added extraction process for geothermal energy systems. The sorbent materials and engineering analysis will be applicable to other industrial processes in which secondary recovery of valuable materials could provide economic benefit.

15 GEOTHERMAL ENERGY↗

Integrated Off-gas System: A Preconceptual Design of an Integrated Off-Gas Treatment System

The U.S. has a vested interest in the advancement of nuclear energy to achieve aggressive net-zero goals, with reprocessing and recycling of used nuclear fuels (UNF) playing a vital role. It will not be possible to meet U.S. regulatory requirements without robust off-gas treatment, so it is crucial to advance treatment technologies to facilitate the design of future reprocessing facilities. For many years, teams of researchers across the U.S. Department of Energy (DOE) National Laboratory complex have been investigating off-gas treatment technologies for the capture and removal of volatile radionuclides (i.e., 85 Kr, Xe, 14 C, and 129 I) and oxides of nitrogen (NO X ) that are produced from reprocessing. These investigations have been focused on developing individual technologies for the capture of Kr, Xe, iodine, and CO 2 . Capture technologies for each constituent were tested independently from one another by utilizing nonradioactive surrogates to simulate simplified off-gas streams. The tests have been relatively small, laboratory-scale experiments of up to approximately 1 L/minute total gas flow rate. To increase the readiness of these technologies for deployment, an integrated test system with a larger-scale capacity is needed to bridge the gap between promising bench scale and fully scalable UNF reprocessing off-gas treatment. This document contains the goals, design basis, functional requirements, preconceptual design, and cost estimates for an integrated off-gas demonstration system for the capture and removal of NO x , Kr, Xe, CO 2 , and iodine at 10× higher throughput than earlier laboratory studies. The order-of-magnitude cost estimate for this system is approximately $\$$886,000. Next phases include conceptual design, detailed design, fabrication, and commissioning.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Discontinuous Galerkin discretization for quantum simulation of chemistry

All-electron electronic structure methods based on the linear combination of atomic orbitals method with Gaussian basis set discretization offer a well established, compact representation that forms much of the foundation of modern correlated quantum chemistry calculations—on both classical and quantum computers. Despite their ability to describe essential physics with relatively few basis functions, these representations can suffer from a quartic growth of the number of integrals. Recent results have shown that, for some quantum and classical algorithms, moving to representations with diagonal two-body operators can result in dramatically lower asymptotic costs, even if the number of functions required increases significantly. We introduce a way to interpolate between the two regimes in a systematic and controllable manner, such that the number of functions is minimized while maintaining a block-diagonal structure of the two-body operator and desirable properties of an original, primitive basis. Techniques are analyzed for leveraging the structure of this new representation on quantum computers. Empirical results for hydrogen chains suggest a scaling improvement from O(N 4.5 ) in molecular orbital representations to O(N 2.6 ) in our representation for quantum evolution in a fault-tolerant setting, and exhibit a constant factor crossover at 15 to 20 atoms. Moreover, we test these methods using modern density matrix renormalization group methods classically, and achieve excellent accuracy with respect to the complete basis set limit with a speedup of 1–2 orders of magnitude with respect to using the primitive or Gaussian basis sets alone. These results suggest our representation provides significant cost reductions while maintaining accuracy relative to molecular orbital or strictly diagonal approaches for modest-sized systems in both classical and quantum computation for correlated systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗