Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bounded objective 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Diagnosing Open Cluster Stock 2: Member Candidates and Mass Distribution with Gaia DR2 and LAMOST

We identify 1325 member candidates of the open cluster (OC) Stock 2 using data from Gaia DR2. We use the algorithms Clusterix 2.0 and HDBSCAN to select cluster candidates and further refine the final cluster membership by defining neighbors in 5D phase space (X {sub cp}, Y {sub cp}, Z{sub cp},κ⋅μ{sub α}{sup ∗}/ϖ, κ · μ {sub δ}/ϖ). Among these candidates, less than half have G, G {sub BP}, and G {sub RP} extinctions from Gaia. When Gaia extinctions are unavailable, we compute extiction using empirical formulas and E(B − V) = 0.350. We analyze the spatial distribution and mass profile of Stock 2. Our results reveal Stock 2 is still a bound OC and we find evidence of mass segregation. By comparing initial mass functions, the present-day mass function indicates that Stock 2 is a massive stellar cluster with a mass of 4000 M {sub ⊙}. The core radius and tidal radius, calculated via the radial density profile and total mass, are 3.97 pc and 22.65 pc, respectively. Common stars between our selected member candidates and the Large Sky Area Multi-Object Fiber Spectroscopic Telescope DR7 medium-resolution catalog give a metalliclity of [Fe/H] = −0.040 ± 0.147.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗
DOE OSTI2023

Convex Q-Learning in Continuous Time with Application to Dispatch of Distributed Energy Resources

Convex Q-learning is a recent approach to reinforcement learning, motivated by the possibility of a firmer theory for convergence, and the possibility of making use of greater a priori knowledge regarding policy or value function structure. This paper explores algorithm design in the continuous time domain, with a finite-horizon optimal control objective. The main contributions are (i) The new Q-ODE: a model-free characterization of the Hamilton-Jacobi-Bellman equation. (ii) A formulation of Convex Q-learning that avoids approximations appearing in prior work. The Bellman error used in the algorithm is defined by filtered measurements, which is necessary in the presence of measurement noise. (iii) Convex Q-learning with linear function approximation is a convex program. It is shown that the constraint region is bounded, subject to an exploration condition on the training input. (iv) The theory is illustrated in application to resource allocation for distributed energy resources, for which the theory is ideally suited.

Lu, Fan↗
DOE OSTI2025

Hybrid epoxy–acrylate resins for wavelength-selective multimaterial 3D printing

Structures in nature have evolved to combine hard and soft materials in precise 3D arrangements, which imbues bulk properties and functionality that remain elusive to mimic synthetically. However, the potential for biomimetic analogs to seamlessly interface hard materials with soft surfaces for applications ranging from robotics and sealants to medical devices (e.g., prosthetics and wearable health monitors) has driven the demand for innovative chemistries and manufacturing approaches. Herein, we unveil a liquid resin for rapid, high resolution digital light processing (DLP) 3D printing of multimaterial objects with an unprecedented combination of strength, elasticity, and resistance to aging. Two enabling discoveries are the use of a covalently bound (hybrid) epoxy-acrylate monomer that precludes plasticization of soft domains and a wavelength-selective photosensitizer that accelerates cationic curing for hard domains. Using dual projection for multicolor DLP 3D printing (UV and violet light), several bioinspired metamaterial structures are fabricated, including one with hard springs embedded in a soft cylinder to adjust compressive behavior and a detailed knee joint featuring “bones” and “ligaments” for smooth motion. Lastly, the application of this system to facilitate selective stretching for electronic devices is demonstrated with a proof-of-concept device.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Analysis of the SiMPL Method for Density-Based Topology Optimization

We present a rigorous convergence analysis of a new method for density-based topology optimization that provides pointwise bound-preserving design updates and faster convergence than other popular first-order topology optimization methods. Due to its strong bound preservation, the method is exceptionally robust, as demonstrated in numerous examples here and in the companion article [D. Kim et al., Struct. Multidiscip. Optim., 68 (2025), 74]. Furthermore, it is easy to implement with clear structure and analytical expressions for the updates. Our analysis covers two versions of the method, characterized by the employed line search strategies. We consider a modified Armijo backtracking line search and a Bregman backtracking line search. For both line search algorithms, our algorithm delivers a strict monotone decrease in the objective function and further intuitive convergence properties, e.g., strong and pointwise convergence of the density variables on the active sets, norm convergence to zero of the increments, convergence of the Lagrange multipliers, and more. In addition, the numerical experiments demonstrate apparent mesh-independent convergence of the algorithm. Here, we refer to the new algorithm as the SiMPL method (pronounced “simple”), which stands for Sigmoidal Mirror descent with a Projected Latent variable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Possibility of 14 C cluster as a building block of medium-mass nuclei

The possibility of the 14 C cluster being a basic building block of medium-mass nuclei is discussed. Although α cluster structures have been widely discussed in the light N ≈ Z mass region, the neutron-to-proton ratio deviates from unity in the nuclei near β-stability line and in neutron-rich nuclei. Thus, more neutron-rich objects with N > Z could become the building blocks of cluster structures in such nuclei. The 14 C nucleus is strongly bound and can be regarded as such a candidate. In addition, the path to the lowest shell-model configuration at short relative distances is closed for the 14 C+ 14 C structure, contrary to the case of the 12 C+ 12 C structure; this allows an appreciable separation distance between the 14 C clusters. The recent development of the antisymmetrized quasicluster model (AQCM) allows us to utilize a j j-coupling shell-model wave function for each cluster in a simplified way. The AQCM results for the 14 C+ 14 C structure in 28 Mg are compared with the ones of cranked relativistic mean field (CRMF) calculations. Although theoretical frameworks of these two models are quite different, they give similar results for the nucleonic densities and rotational properties of the structure under investigation. The existence of a linear chain three- 14 C cluster structure in 42 Ar has also been predicted in AQCM. Furthermore, these results confirm the role of the 14 C cluster as a possible building block of cluster structures in medium-mass nuclei.

20 ≤ A ≤ 38↗
DOE OSTI2023

Search for Z' bosons decaying to pairs of heavy Majorana neutrinos in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

A search for the production of pairs of heavy Majorana neutrinos (N ℓ ) from the decays of Z' bosons is performed using the CMS detector at the LHC. The data were collected in proton-proton collisions at a center-of-mass energy of $\sqrt s$ = 13 TeV, with an integrated luminosity of 138 fb –1 . The signature for the search is an excess in the invariant mass distribution of the final-state objects, two same-flavor leptons (e or μ) and at least two jets. No significant excess of events beyond the expected background is observed. Upper limits at 95% confidence level are set on the product of the Z' production cross section and its branching fraction to a pair of N ℓ , as functions of N ℓ and Z' boson masses (m N ℓ and m Z' , respectively) for m Z' from 0.4 to 4.6 TeV and m N ℓ from 0.1 TeV to m Z' /2. In the theoretical framework of a left-right symmetric model, exclusion bounds in the m N ℓ -m Z' plane are presented in both the electron and muon channels. The observed upper limit on m Z' reaches up to 4.42 TeV. These are the most restrictive limits to date on the mass of N ℓ as a function of the Z' boson mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

IDAES-PSE 2.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.0.0 Release Highlights Removal of deprecated features from IDAES v1 Update to Pyomo v6.5 – this required a number of updates to support the new NL solver writer and to address some changes in Pyomo Creation of new testing suite for backward compatibility, model robustness and verification More general implementation of the Helmholtz EoS. This brings some new features like standard property diagrams, choice of mass or mole basis, and new state variable options Standardizing names in Heat Exchanger models (breaking change from v2.0.0a2): Control Volumes named hot_side and cold_side Ports names hot_side_inlet, hot_side_outlet, cold_side_inlet and cold_side_outlet Config Blocks names hot_side_config and cold_side_config Config arguments for user provided names for each side: hot_side_name and cold_side_name. Updating Keras surrogate tool to use v1.1 of OMLT New prototype API for model initialization (idaes.core.initialization) The new API uses "Model Initializer" objects instead of class methods, allowing for the definition of multiple initialization routines for a single model A number of common, model agnostic initialization routines have also been defined, including initialization from data, block-decomposition and a general hierarchical approach equivalent to the existing method for common unit models New metadata for thermophysical properties – valid_range This can be used to record the range of values over which a property value can be trusted, such as the range of experimental data used to regress parameters A number of new utility functions have been added to check for properties with values outside the valid range and to set bounds based on this metadata Updated construction of balance expressions in Control Volumes to remove unneeded terms In the past, unneeded terms were added as a constant 0 term, however they will now be dropped entirely from the expression This was necessary due to more strict unit checking in the new Pyomo solver writer which no longer ignores 0 terms Updates to metadata for thermophysical properties to better define known properties and units of measurement This results in more strict enforcement of standard naming for thermophysical and reaction properties Users can still define custom properties, but these must be done explicitly using the define_custom_properties() method instead of being implicitly created by add_property() Updated convergence tester utility tool to support definition of benchmark files (JSON format) and comparison of performance to benchmarks Set default iteration limit for IPOPT in IDAES config to 200 iterations Update scaling of example models to work with new Pyomo NL solver writer Improve testing of extensions and examples infrastructure to avoid need for downloading files Updated distillation column to centralize common functionality and remove a number of Pyomo warnings

IDAES↗
DOE OSTI2025

QCD Constraints on Isospin-Dense Matter and the Nuclear Equation of State

Understanding the behavior of dense hadronic matter is a central goal in nuclear physics as it governs the nature and dynamics of astrophysical objects such as supernovae and neutron stars. Because of the nonperturbative nature of quantum chromodynamics (QCD), little is known rigorously about hadronic matter in these extreme conditions. Here, lattice QCD calculations are used to compute thermodynamic quantities and the equation of state of QCD over a wide range of isospin chemical potentials with controlled systematic uncertainties. Agreement is seen with chiral perturbation theory when the chemical potential is small. Comparison to perturbative QCD at large chemical potential allows for an estimate of the gap in the superconducting phase, and this quantity is seen to agree with perturbative determinations. Since the partition function for an isospin chemical potential μ I bounds the partition function for a baryon chemical potential μ B = 3 μ I / 2 , these calculations also provide rigorous nonperturbative QCD bounds on the symmetric nuclear matter equation of state over a wide range of baryon densities for the first time. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Augmenting subspace optimization methods with linear bandits

In this work, we consider the framework of methods for unconstrained minimization that are, in each iteration, restricted to a model that is only a valid approximation to the objective function on some affine subspace containing an incumbent point. These methods are of practical interest in computational settings where derivative information is either expensive or impossible to obtain. Recent attention has been paid in the literature to employing randomized matrix sketching for generating the affine subspaces within this framework. We consider a relatively straightforward, deterministic augmentation of such a generic subspace optimization method. In particular, we consider a sequential optimization framework where actions consist of one-dimensional linear subspaces and rewards consist of (approximations to) the magnitudes of directional derivatives computed in the direction of the action subspace. Reward maximization in this context is consistent with maximizing lower bounds on descent guaranteed by first-order Taylor models. This sequential optimization problem can be analysed through the lens of dynamic regret. We modify an existing linear upper confidence bound (UCB) bandit method and prove sublinear dynamic regret in the subspace optimization setting. We demonstrate the efficacy of employing this linear UCB method in a setting where forward-mode algorithmic differentiation can provide directional derivatives in arbitrary directions and in a derivative-free setting. For the derivative-free setting, we propose SS-POUNDers, an extension of the derivative-free optimization method POUNDers that employs the linear UCB mechanism to identify promising subspaces. Our numerical experiments suggest a preference, in either computational setting, for employing a linear UCB mechanism within a subspace optimization method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evaluating Space Object Conjunction Probabilities Using Characteristic Function Inversion

This report discusses an approach to computing the probability of a conjunction between two space objects in the short-term encounter scenario. A conjunction is defined here as an event where the miss distance between the objects is less than some specified value. The scenario assumptions are that the motion of the objects is linear, their positions are Gaussian distributed, and their velocities are known and constant. Under these assumptions, the squared-miss distance is shown to have the generalized chi-square distribution. An established statistical technique called characteristic function inversion is employed to evaluate the distribution and obtain conjunction probabilities. The method is closely related to a recent approach based on moment generating function inversion, and a qualitative comparison of the approaches is provided. Last, the method is tested on several benchmark test cases where it agrees with numerical integration on the cases with conjunction probabilities above 10 –12 . However, the exact probability in these cases is usually not needed and this probability can be bounded above using an independent Gaussian approximation. Overall, the report shows how to compute conjunction probabilities using a standard statistical method, though numerical integration seems to perform equally well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Scalable, In-situ Data Clustering Data Analysis for Extreme Scale Scientific Computing (Final Report)

The objective of this project is to address challenges in the design and development of scalable in-situ data clustering and analytics algorithms and software. Our goal is to develop parallel software consisting of a set of spatio-temporal data clustering and anomaly detection functions, both of which are very important for large-scale analysis and have wide applicability for in-situ runs as well as post-processing analysis. Our design principles for in-situ analysis consider the following: (1) identify parts of the computation can be done close to the data within the nodes, while it is still in memory; (2) extract analysis components can (and should) be performed in remote staging and analysis nodes; (3) develop error-bound approximation methods for applications tolerable for small errors; (4) identify the type of derived distributions and statistics, for spatio-temporal data, that can be kept locally in order to both accelerate computations and meet energy constraints in subsequent iterations and phases; (5) use a self-describing data format so that data can be consistent and understood among local storage (memory and SSDs) and at staging and analysis nodes, thereby providing portability and flexibility; (6) develop service-oriented functions that can schedule in-situ and post-hoc analysis tasks based on the dynamic requirements of applications. Our development focus is to produce the parallel data analysis software/library that will be scalable, reusable, extensible, and generic for applications in different disciplines. The software will be able to run in-situ with the simulations as well as post-hoc analysis. This approach will satisfy many synergistic requirements for data intensive applications executed on data coming from instruments and experiments. In particular, the proposed multilevel approach is directly applicable to perform design tradeoffs for running part of the algorithms near the instruments and the rest on remote (analysis) systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

gC1qR/C1qBP/HABP-1: Structural Analysis of the Trimeric Core Region, Interactions With a Novel Panel of Monoclonal Antibodies, and Their Influence on Binding to FXII

The protein gC1qR/C1qBP/HABP-1 plays an essential role in mitochondrial biogenesis, but becomes localized at the cellular surface in numerous pathophysiological states. When this occurs on endothelial cells, surface-exposed gC1qR activates the classical pathway of complement. It also promotes assembly of a multi-protein complex comprised of coagulation factor XII (FXII), pre-kallikrein (PK), and high-molecular weight kininogen (HMWK) that activates the contact system and the kinin-generating system. Since surface-exposed gC1qR triggers intravascular inflammatory pathways, there is interest in identifying molecules that block gC1qR function. Here we further that objective by reporting the outcome of a structure/function investigation of gC1qR, its interactions with FXII, and the impact of a panel of monoclonal anti-gC1qR antibodies on FXII binding to gC1qR. Although deletion mutants have been used extensively to assess gC1qR function, none of these proteins have been characterized structurally. To that end, we determined a 2.2 Å resolution crystal structure of a gC1qR mutant lacking both of its acidic loops, but which retained nanomolar-affinity binding to FXII and FXIIa. This structure revealed that the trimeric gC1qR assembly was maintained despite loss of roughly thirty residues. Characterization of a novel panel of anti-gC1qR monoclonal antibodies identified several with biochemical properties distinct from previously described antibodies, as well as one which bound to the first acidic loop of gC1qR. Intriguingly, we found that each of these antibodies could partly inhibit binding of FXII and FXIIa to gC1qR. Based on these results and previously published studies, we offer new perspectives for developing gC1qR inhibitors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Multistage economic MPC for systems with a cyclic steady state: A gas network case study

Multistage model predictive control (MPC) provides a robust control strategy for dynamic systems with uncertainties and a setpoint tracking objective. Moreover, extending MPC to minimize an economic cost instead of tracking a pre-calculated optimal setpoint improves controller performance. This paper presents a novel multistage economic nonlinear model predictive control (E-NMPC) framework for dynamic systems operating under uncertainty, with specific application to natural gas transmission networks. A key innovation lies in the integration of cyclic steady-state (CSS) constraints within the multistage MPC formulation, enabling the controller to manage periodic operating conditions commonly observed in energy systems. A Lyapunov-based descent condition is enforced to ensure robust stability of the controller. The multistage economic MPC framework is validated on two gas pipeline case studies, where it successfully minimizes net energy consumption, respects operational constraints under uncertain demand profiles, and guides the network to optimal cyclic operation. The Lyapunov function remains bounded in both case studies, validating the robust stability of multistage E-NMPC.

03 NATURAL GAS↗
DOE OSTI2025

Five-point functions and the permutation group 𝑆 5

Five-point functions and five-body wave functions play an important role in many areas of nuclear and particle physics, e.g., in 2 →3 scattering processes, in the five-gluon vertex, or in the study of pentaquarks. In this work we consider the permutation group 𝑆 5 to facilitate the description of such objects. We work out the multiplets transforming under irreducible representations of 𝑆 5 and provide compact formulas allowing one to cast the permutations of an object 𝑓 12345 into combinations with definite permutation symmetry. We also give the explicit expressions for the irreducible multiplet products. We consider several practical applications as examples: We arrange the four-momenta and Lorentz invariants of a five-point function into the multiplet structure, we work out the color tensors of the five-gluon vertex in the multiplet notation, and we discuss applications for five-body wave functions like those of pentaquarks.

Bethe-Salpeter equation↗
DOE OSTI2022

Unraveling small-scale defects in irradiated ThO 2 using kinetic Monte Carlo simulations

Point defects and their clusters generated through irradiation can have significant impact on the physical and mechanical properties of materials. However, direct experimental visualization of these small-scale defects using high-resolution scanning transmission electron microscopy remains a challenging task. In this study, using thorium dioxide (ThO 2 ) with the fluorite structure as a model system, we demonstrate the use of ab initio basin-hopping simulations in synergy with object kinetic Monte Carlo simulations as a powerful tool for identifying small defect complexes in irradiated materials. In addition to providing quantitative insights into defect evolution in ThO 2 under irradiation, our study reveals an unexpected role of bound anti-Schottky defect clusters in mediating defect transport. Remarkably, despite their poor thermal stability against dissociation at high temperatures, the transient formation of bound anti-Schottky defects under irradiation and their subsequent migration provide the dominant mechanism for the growth of large interstitial loops that have been experimentally observed in ThO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reliable and Efficient Machine Learning (Final Technical Report)

Modern scientific experiments generate massive amounts of data at a pace much faster than humans can manually analyze. While machine learning has revolutionized commercial data analysis (such as recommending movies or recognizing faces), applying these tools to complex scientific discovery is challenging because scientific answers must be precise, interpretable, and adhere to physical laws. The research under this project aims to develop new mathematical tools and computer algorithms specifically designed for scientific applications. Major progress has been made in automatically cleaning and deconstructing messy experimental data, analyzing the visual information of physical phenomena, determining the underlying physical variables, and providing rig orous mathematical analysis of interesting algorithms and concepts widely used in machine learning. This project addressed the critical gap between our ability to generate massive scientific data and our ability to extract interpretable information from it. We established mathematical foundations for Scientific Machine Learning (SciML) aimed at effective data analytics and automated discovery. Our work focused on three core objectives: (1) developing reliable feature extraction methods for dynamic high-dimensional data, (2) establishing mathematical foundations for discovering dynamics via neural networks, and (3) creating rigorous optimization techniques for these models. Key outcomes come from two fronts. On the practical side, they include the development of algorithms that significantly enhance the extraction of signals from field data, as well as the capability to handle situations that exhibit smooth variations or physical stretching due to temperature changes. They also include the creation of an automated framework for discovering fundamental state variables from raw experimental data, demonstrating the ability to identify intrinsic physical dimensions without prior knowledge of the governing laws. On the theoretical front, the research results in theoretical advances in Optimal Transport, a widely used notion in SciML, specifically regarding functions with fixed-size nodal sets, provide sharp bounds relevant to uncertainty quantification. Meanwhile, the outcomes also include the establishment of convergence theories for nonlocal gradient descent methods, enabling robust optimization with noisy data in high-dimensional settings commonly encountered in scientific modeling. The project also helps creating opportunities to train the next generation of researchers, equipping them with the necessary technical skills for today’s workplace and preparing them for future advances.

97 MATHEMATICS AND COMPUTING↗