SEARCH · Engineering Papers

Results for “Complex Constraints”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2024

CHMMPP: A c++ library for constrained Hidden Markov Models

SAND2024-13027O The CHMMPP: A c++ Library for Constrained Hidden Markov Models (HMM) software supports the analysis of multivariate time series data to detect patterns using HMM. Many applications involve the detection and characterization of hidden or latent states in a complex system using observable states and variables. This software supports inference of latent states integrating both an HMM and application-specific constraints that reflect known relationships in hidden states. The CHMMPP software supports application-specific and generic methods for constrained inference. This includes a framework for customized Viterbi methods, constrained inference of hidden states with A* and integer programming methods, and various constraint-informed methods for learning HMM model parameters. CHMMPP focuses on supporting generic methods that enable the agile expression of complex sets of constraints that naturally arise in many real-world applications.

Hart, William↗

DOE OSTI2025

A Smoothed Augmented Lagrangian Framework for Convex Optimization with Nonsmooth Constraints

Augmented Lagrangian (AL) methods have proven remarkably useful in solving optimization problems with complicated constraints. The last decade has seen the development of overall complexity guarantees for inexact AL variants. Yet, a crucial gap persists in addressing nonsmooth convex constraints. To this end, we present a smoothed augmented Lagrangian (AL) framework where nonsmooth terms are progressively smoothed with a smoothing parameter $\eta _k$ . The resulting AL subproblems are $\eta _k$ -smooth, allowing for leveraging accelerated schemes. By a careful selection of the inexactness level $\epsilon _k$ (for inexact subproblem resolution), the penalty parameter $\rho _k$ , and smoothing parameter $\eta _k$ at epoch k, we derive rate and complexity guarantees of $\tilde{\mathcal {O}}(1/{\varepsilon }^{3/2})$ and $\tilde{\mathcal {O}}(1/{\varepsilon })$ in convex and strongly convex regimes for computing an ${\varepsilon }$ -optimal solution, when $\rho _k$ increases at a geometric rate, a significant improvement over the best available guarantees for AL schemes for convex programs with nonsmooth constraints. Analogous guarantees are developed for settings with $\rho _k = \rho$ as well as $\eta _k = \eta$ . Preliminary numerics on a fused Lasso problem display promise.

augmented Lagrangian↗

DOE OSTI2024

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Structural constraint integration in a generative model for the discovery of quantum materials

Billions of organic molecules have been computationally generated, yet functional inorganic materials remain scarce due to limited data and structural complexity. Here, in this work, we introduce Structural Constraint Integration in a GENerative model (SCIGEN), a framework that enforces geometric constraints, such as honeycomb and kagome lattices, within diffusion-based generative models to discover stable quantum materials candidates. SCIGEN enables conditional sampling from the original distribution, preserving output validity while guiding structural motifs. This approach generates ten million inorganic compounds with Archimedean and Lieb lattices, over 10% of which pass multistage stability screening. High-throughput density functional theory calculations on 26,000 candidates shows over 95% convergence and 53% structural stability. A graph neural network classifier detects magnetic ordering in 41% of relaxed structures. Furthermore, we synthesize and characterize two predicted materials, TiPd 0.22 Bi 0.88 and Ti 0.5 Pd 1.5 Sb, which display paramagnetic and diamagnetic behaviour, respectively. Our results indicate that SCIGEN provides a scalable path for generating quantum materials guided by lattice geometry.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Directed Self-Assembly of the Organic Semiconductor C8-BTBT-C8 in Anodic Aluminum Oxide Nanopores

Controlling the self-assembly of organic semiconductors at the nanoscale is critical for advancing high-performance electronic and photonic devices, yet it remains challenging due to their intrinsic anisotropic crystallization and sensitivity to processing conditions. Here, in this study, we demonstrate that cylindrical nanoconfinement within anodic aluminum oxide membranes provides a versatile platform to precisely tune the molecular orientation and phase behavior of the prototypical organic semiconductor 2,7-dioctyl[1]benzothieno[3,2-b][1]benzothiophene (C8-BTBT-C8). Combining temperature-dependent high-resolution synchrotron X-ray scattering with optical birefringence measurements, we uncover that confinement geometries (pore diameters 25–180 nm) and surface chemistry govern the emergence of distinct smectic A textures, featuring molecular layers either parallel or perpendicular to the pore axis. The competition between axial and radial smectic layering is modulated by pore size, surface hydrophilicity, and thermal history, enabling reversible control over domain orientations and transitions between liquid crystalline and crystalline states. Notably, nanoconfinement stabilizes the smectic phase over an expanded temperature range compared to bulk, while inducing complex multidomain configurations owing to geometric constraints and anchoring conditions. Our results elucidate fundamental mechanisms by which anisotropic nanoscale confinement directs the self-organization of highly conjugated organic molecules, with implications for optimizing directional charge transport and anisotropic optical responses in organic–inorganic hybrid nanoarchitectures. This study establishes nanoconfinement as a powerful strategy to engineer morphology and functional properties in organic semiconducting materials with nanoscale precision.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Machine learning method for enforcing variable independence in background estimation with LHC data: ABCDisCoTEC

A novel solution is presented for the problem of estimating the backgrounds of a signal search using observed data while simultaneously maximizing the sensitivity of the search to the signal. The 'ABCD method' provides a reliable framework for background estimation by partitioning events into one signal-enhanced region (A) and three background-enhanced control regions (B, C, and D) via two smoothly varying, statistically independent variables. In practice, even slight correlations between the two variables can significantly undermine the method's performance. Thus, choosing appropriate variables by hand can present a formidable challenge, especially when background and signal differ only subtly. To address this issue, the ABCD with distance correlation (ABCDisCo) method was developed to construct two learned variables via a neural network trained to provide strong signal-background discrimination with small values of the distance correlation (DisCo) measure between the two learned variables. However, relying solely on minimizing the DisCo can result in learned variables that may not have distributions of background events that are smoothly varying and localized at extreme values, as necessary for the validity of the background estimation. The ABCDisCo training enhanced with closure (ABCDisCoTEC) method is introduced to solve this issue by directly minimizing the nonclosure, expressed as a dedicated differentiable loss term. This extended method is applied to a data set of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the CMS detector at the CERN Large Hadron Collider. Additionally, given the complexity of the minimization problem with constraints on multiple loss terms, the modified differential method of multipliers is applied and shown to greatly improve the stability and robustness of the ABCDisCoTEC method, compared to grid search hyperparameter optimization procedures.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

DOE OSTI2024

Energy Technology Proving Ground Program Plan

New methods of energy production and distribution are required to meet clean energy goals and demands across all U.S. energy sectors. Idaho National Laboratory’s (INL) Integrated Energy Systems (IES) initiative is enabling clean energy research, development, and demonstration (RD&D) activities. To date, IES demonstration programs have been limited by distributed infrastructure and a lack of large-scale facilities to accommodate industry-scale research of Technical Readiness Level (TRL) 6-8 technologies. The IES initiative plans to eliminate these constraints by establishing a new research complex at INL known as the Energy Technology Proving Ground (Proving Ground) to be led by the Energy and Environment Science and Technology Directorate. The Energy and Environment Science and Technology (EES&T) directorate, one of five Idaho National Laboratory (INL) RD&D organizations, focuses on clean energy technologies that anchor the industry-enabling research of the Proving Ground. The Proving Ground will combine diverse clean energy systems into lean integrated test bed of independent multiscale capabilities available to the government and commercial industries to perform research; and will enable INL’s goal of becoming a Net-Zero entity by 2031. This program encompasses existing and new research space at both the in-town Research and Education Campus (REC) and the Arco desert site (the Site). The Proving Ground will support the maturation of IES technologies from TRL 1 through 8 by providing the infrastructure and capabilities needed to sustain a continuum of RD&D from basic science to industry-scale. To establish The Proving Ground and meet INL’s net-zero goals by 2031, nine research program areas have been identified within the IES initiative that require expanded and new capital infrastructure. This program plan provides guidance for establishing the Proving Ground at the Site for plug-and-play pilot testing and proofing of integrated energy system functionality including fission and renewable energy sources for industry driven application platforms.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2025

Energy Technology Proving Ground FY-2026 Program Plan (Rev.1)

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

Search for light pseudoscalar boson pairs produced from Higgs boson decays using the 4τ and 2μ2τ final states in proton-proton collisions at $\sqrt{s}=13$ TeV

A search for a pair of light pseudoscalar bosons (a 1 ) produced in the decay of the 125 GeV Higgs boson is presented. The analysis examines decay modes where one a 1 decays into a pair of tau leptons and the other decays into either another pair of tau leptons or a pair of muons. The a1 boson mass probed in this study ranges from 4 to 15 GeV. The data sample was recorded by the CMS experiment in proton-proton collisions at a center-of-mass energy of 13 TeV and corresponds to an integrated luminosity of 138 fb −1 . No excess above standard model (SM) expectations is observed. The study combines the 4τ and 2μ2τ channels to set upper limits at 95% confidence level (CL) on the product of the Higgs boson production cross section and the branching fraction to the 4τ final state, relative to the Higgs boson production cross section predicted by the SM. In this interpretation, the a 1 boson is assumed to have Yukawa-like couplings to fermions, with coupling strengths proportional to the respective fermion masses. The observed (expected) upper limits range between 0.007 (0.011) and 0.079 (0.066) across the mass range considered. The results are also interpreted in the context of models with two Higgs doublets and an additional complex singlet field (2HD+S). The tightest constraints are obtained for the Type III 2HD+S model. In this case, assuming the Higgs boson production cross section equals the SM prediction, values of the branching ratio for the Higgs boson decay into a pair of a1 bosons exceeding 16% are excluded at 95% CL for a1 boson masses between 5 and 15 GeV and tan β > 2, with the exception of scenarios in which the a1 boson mixes with charm or bottom quark-antiquark bound states.

Higgs physics↗

DOE OSTI2025

The Atacama Cosmology Telescope: DR6 power spectrum foreground model and validation

We discuss the model of astrophysical emission at millimeter wavelengths used to characterize foregrounds in the multi-frequency power spectra of the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6), expanding on Louis et al. (2025) (2503.14452). We detail several tests to validate the capability of the DR6 parametric foreground model to describe current observations and complex simulations, and show that cosmological parameter constraints are robust against model extensions and variations. We demonstrate consistency of the model with pre-DR6 ACT data and observations from Planck and the South Pole Telescope. We evaluate the implications of using different foreground templates and extending the model with new components and/or free parameters. In all scenarios, the DR6 ΛCDM and ΛCDM+N eff cosmological parameters shift by less than 0.5σ relative to the baseline constraints. Some foreground parameters shift more; we estimate their systematic uncertainties associated with modeling choices. From our constraint on the kinematic Sunyaev-Zel'dovich power, we obtain a conservative limit on the duration of reionization of Δz rei < 4.4, assuming a reionization midpoint consistent with optical depth measurements and a minimal low-redshift contribution, with varying assumptions for this component leading to tighter limits. Finally, we analyze realistic non-Gaussian, correlated microwave sky simulations containing Galactic and extragalactic foreground fields, built independently of the DR6 parametric foreground model. Processing these simulations through the DR6 power spectrum and likelihood pipeline, we recover the input cosmological parameters of the underlying cosmic microwave background field, a new demonstration for small-scale CMB analysis. These tests validate the robustness of the ACT DR6 foreground model and cosmological parameter constraints.

CMBR experiments↗

DOE OSTI2025

Operator origin of anomalous dimensions in de Sitter space

The late-time limit of the power spectrum for heavy (principal series) fields in de Sitter (dS) space yields a series of polynomial terms with complex scaling dimensions. Such scaling behavior is expected to result from an associated operator with a complex dimension. In a free theory, these complex dimensions are known to match the constraints imposed by unitarity on the space of states. Yet, perturbative corrections to the scaling behavior of operators are naively inconsistent with unitary evolution of the quantum fields in dS space. This paper demonstrates how to compute one-loop corrections to the scaling dimensions that appear in the two-point function from the field theory description in terms of local operators. We first show how to evaluate these anomalous dimensions using Mellin space, which has the feature that it naturally accommodates a scaleless regulator. We then explore the consequences for the soft de Sitter effective theory (SdSET) description that emerges in the long wavelength limit. Carefully matching between the UV and SdSET descriptions requires the introduction of novel nondynamical “operators” in the effective theory. This is not only necessary to reproduce results extracted from the Källén-Lehmann representation (that use the space of unitary states directly), but it is also required by general arguments that invoke positivity. Published by the American Physical Society 2025

Cohen, Timothy↗

DOE OSTI2026

Design Status and Experimental Strategy for Initial Molten Salt Irradiation Experiments in HFIR

The development and deployment of molten salt reactor (MSR) technologies require experimental capabilities that can evaluate molten fuel salt behavior and structural material performance under representative irradiation conditions. Although modeling and separate effects testing provide important insight, there remains a critical lack of in-pile data that capture the coupled effects of neutron irradiation, temperature, salt chemistry, and time. Informed by lessons learned from historical MSR programs and recent international irradiation efforts, this report presents a structured approach to addressing existing gaps in molten salt irradiation testing using the High Flux Isotope Reactor (HFIR). A phased irradiation strategy is presented that incrementally increases experimental complexity while managing cost, risk, and facility constraints. The framework progresses from passive, nonfueled static capsule experiments to fueled, instrumented, and ultimately circulating molten salt systems, providing a pathway for materials screening, mechanistic understanding, and qualification-relevant testing. The report defines a near-term Phase 1 passive capsule concept, associated irradiation conditions, and a conceptual post-irradiation examination strategy. Ongoing modeling, out-of-pile testing, and safety analysis activities are described to support continued capability maturation and to enable future phases of molten salt irradiation testing in support of MSR technology deployment.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

How the Galaxy–Halo Connection Depends on Large-scale Environment

We investigate the connection between galaxies, dark matter halos, and their large-scale environments at z = 0 with Illustris TNG300 hydrodynamic simulation data. We predict stellar masses from subhalo properties to test two types of machine learning (ML) models: explainable boosting machines (EBMs) with simple galaxy environment features and E(3)-invariant graph neural networks (GNNs). The best-performing EBM models leverage spherically averaged overdensity features on 3 Mpc scales. Interpretations via SHapley Additive exPlanations also suggest that in the context of the TNG300 galaxy–halo connection, simple spherical overdensity on ∼3 Mpc scales is more important than cosmic web distance features measured using the DisPerSE algorithm. Meanwhile, a GNN with connectivity defined by a fixed linking length, L, outperforms the EBM models by a significant margin. As we increase the linking length scale, GNNs learn important environmental contributions up to the largest scales we probe (L = 10 Mpc). We conclude that 3 Mpc distance scales are most critical for describing the TNG galaxy–halo connection using the spherical overdensity parameterization, but that information on larger scales, which is not captured by simple environmental parameters or cosmic web features, can further augment these models. Our study highlights the benefits of using interpretable ML algorithms to explain models of astrophysical phenomena, and the power of using GNNs to flexibly learn complex relationships directly from data while imposing constraints from physical symmetries.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

X-ray Spectroscopy Characterization of Electronic Structure and Metal–Metal Bonding in Dicobalt Complexes

Developing multimetallic complexes with tunable metal–metal interactions has long been a target of synthetic inorganic chemistry efforts due to the unique properties that such compounds can exhibit. However, understanding relationships between metal–metal bonding and chemical properties is challenging due to system-dependent factors that influence metal–metal and metal–ligand interactions, including ligand identity, coordination geometry, and metal–metal distance. In this work, we apply X-ray absorption and emission spectroscopy and quantum chemical calculations to describe electronic structure and bonding in a series of dicobalt complexes. The compounds with silane ligands and pseudo-octahedral coordination geometry exhibit Co–Co σ and multicentered bonding character, which we characterize from both the occupied and vacant perspectives via their contributions to the Co X-ray emission and absorption spectra, respectively. In contrast, the dicobalt complexes with a pseudotetrahedral coordination environment do not exhibit Co–Co bonding due to symmetry constraints on orbital overlap. We extend these insights to the theoretical evaluation of related dicobalt complexes to explain how ligand coordination and symmetry dictate the presence or absence of a Co–Co bond. In conclusion, this work highlights how fundamental insights into electronic structure and bonding through X-ray spectroscopy uncover important factors governing metal–metal interactions and guide the rational design of multimetallic complexes with tunable metal–metal bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Robust wind farm layout optimization

Wake interactions in wind farms cause losses in annual energy production (AEP) on the order of 10%. Wind farm designers optimize the layout of the farm to mitigate wake losses, especially in the dominant site-specific wind directions. As wind turbines and wind farms grow in scale, optimization becomes more complex. Offshore wind farms regularly comprise more than 100 wind turbines and are characterized by complex boundaries due to shipping lanes, neighboring wind farms, and other constraints. Layout optimization methods are broadly split between gradient-based and gradient-free approaches. Gradient-based approaches can converge quickly and perform well for smaller, academic problems but are often sensitive to initial conditions and tuning parameters and require expert knowledge to use. On the other hand, gradient-free approaches can be more robust to problem complexities. We present a robust layout optimization approach based on a random search algorithm. The algorithm is intended for those who are not optimization experts and has few tuning parameters that need specification to achieve satisfactory results. Unlike off-the-shelf methods, which use generally available, non-domain-specific optimization routines that accept as inputs an optimization function and constraint definitions, this approach takes advantage of the relative computational costs of the different evaluations by evaluating cheaper computations first (boundary and minimum distance constraints) and running expensive AEP evaluations only if all other checks pass. Moreover, an outer genetic algorithm allows multiple solutions to evolve in parallel, enabling rapid solution development on high-performance computers. We discuss the relative ease of selecting necessary tuning parameters and demonstrate the efficacy of the genetic random search on a complex layout problem consisting of placing 70 turbines in a nonconvex and unconnected boundary region.

17 WIND ENERGY↗

DOE OSTI2025

Incorporating corrosion design constraints in desalination process optimization: A case study in mechanical vapor compression

Corrosion is an expensive and complex challenge for desalination, yet current design approaches do not explicitly account for corrosion mechanisms in process modeling and technoeconomic analysis. Here, to address this gap, we present a workflow for incorporating corrosion design constraints directly into desalination process optimization models. We develop surrogate models for general and localized corrosion metrics as functions of temperature, pH, salinity, dissolved oxygen, and material using data from OLI Systems’ Corrosion Analyzer. We then integrate these surrogates as corrosion design constraints in a cost-optimization MVC model that minimizes the levelized cost of water (LCOW). For a case study of mechanical vapor compression (MVC) treating seawater across a range of recoveries, we find dissolved oxygen (DO) is the dominant driver of localized corrosion, and thus of cost-optimal material choice and operating conditions. Reducing the DO from 8 mg/L to 0.5 mg/L reduces the LCOW by 15-35%, informing the breakeven costs for implementing DO removal or selecting highly corrosion-resistant alloys. This framework is broadly applicable across corrosion types, materials, and components and enables desalination process design that minimizes capital costs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Physicochemical properties of digital light processing 3D-Printed alumina and mullite ceramics

Alumina ceramics fabricated using conventional techniques such as uniaxial pressing or injection molding are popular due to their low density, excellent insulation, and mechanical properties. However, these methods often limit the fabrication of complex geometries with high dimensional accuracy due to tooling constraints and limited design freedom. To overcome these limitations, this study employed Digital Light Processing (DLP) additive manufacturing (AM), which enables the production of precise structures. In addition, the optimization of sintering parameters to enhance the densification and performance of alumina and mullite ceramics was investigated, with a specific focus on how varying sintering temperatures and hold times affect part shrinkage, geometric accuracy, and material integrity. Flexural strength of both alumina and mullite specimens was clearly influenced by the way the layers was stacked. When layers were arranged across the direction of the applied load (Z = 4; XZ), the strength was higher than when they were stacked along the same direction as the load (Z = 3; XY). Weibull analysis based on these results showed high modulus values across all samples, indicating good reliability in the flexural strength measurements. This reliability, combined with the clear influence of sintering parameters, highlights how processing conditions effect the material properties and mechanical performance of parts produced through DLP additive manufacturing. In conclusion, these findings open new avenues for ceramic AM across various applications, enhancing the potential for innovation in fields such as aerospace, biomedical engineering, and energy.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Physics-tailored machine learning reveals unexpected physics in dusty plasmas

Dusty plasma is a mixture of ions, electrons, and macroscopic charged particles that is commonly found in space and planetary environments. The particles interact through Coulomb forces mediated by the surrounding plasma, and as a result, the effective forces between particles can be nonconservative and nonreciprocal. Machine learning (ML) models are a promising route to learn these complex forces, yet their structure should match the underlying physical constraints to provide useful insight. Here, we demonstrate and experimentally validate an ML approach that incorporates physical intuition to infer force laws in a laboratory dusty plasma. Trained on 3D particle trajectories, the model accounts for inherent symmetries, nonidentical particles, and learns the effective nonreciprocal forces between particles with exquisite accuracy (R 2 > 0.99). We validate the model by inferring particle masses in two independent yet consistent ways. The model’s accuracy enables precise measurements of particle charge and screening length, identifying large deviations from common theoretical assumptions. Our ability to identify unknown physics from experimental data demonstrates how ML-powered approaches can guide new routes of scientific discovery in many-body systems. Furthermore, we anticipate our ML approach to be a starting point for inferring laws from dynamics in a wide range of many-body systems, from colloids to living organisms.

Science & Technology - Other Topics↗