Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “line search”

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 19 records
All sourcesNASA NTRSDOE OSTI
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 OSTI2025

A unified funnel restoration SQP algorithm

We consider nonlinearly constrained optimization problems and discuss a generic double-loop framework consisting of basic algorithmic ingredients that unifies a broad range of nonlinear optimization solvers. This framework has been implemented in the open-source solver Uno, a Swiss Army knife-like C++ optimization framework that unifies many nonlinearly constrained nonconvex optimization solvers. We illustrate the framework with a sequential quadratic programming (SQP) algorithm that maintains an acceptable upper bound on the constraint violation, called a funnel, that is monotonically decreased to control the feasibility of the iterates. Infeasible quadratic subproblems are handled by a feasibility restoration strategy. Globalization is controlled by a line search or a trust-region method. We prove global convergence of the trust-region funnel SQP method, building on known results from filter methods. We implement the algorithm in Uno, and we provide extensive test results for the trust-region line-search funnel SQP on small CUTEst instances.

Kiessling, David [Katholieke Univ. Leuven, Heverle↗
DOE OSTI2025

Optimizing the optimizer for physics-informed neural networks and Kolmogorov-Arnold networks

Physics-Informed Neural Networks (PINNs) have revolutionized the computation of PDE solutions by integrating partial differential equations (PDEs) into the neural network’s training process as soft constraints, becoming an important component of the scientific machine learning (SciML) ecosystem. More recently, physics-informed Kolmogorv-Arnold networks (PIKANs) have also shown to be effective and comparable in accuracy with PINNs. In their current implementation, both PINNs and PIKANs are mainly optimized using first-order methods like Adam, as well as quasi-Newton methods such as BFGS and its low-memory variant, L-BFGS. However, these optimizers often struggle with highly nonlinear and non-convex loss landscapes, leading to challenges such as slow convergence, local minima entrapment, and (non)degenerate saddle points. In this study, we investigate the performance of Self- Scaled BFGS (SSBFGS), Self-Scaled Broyden (SSBroyden) methods and other advanced quasi-Newton schemes, including BFGS and L-BFGS with different line search strategies. These methods dynamically rescale updates based on historical gradient information, thus enhancing training efficiency and accuracy. We systematically compare these optimizers – using both PINNs and PIKANs – on key challenging PDEs, including the Burgers, Allen-Cahn, Kuramoto-Sivashinsky, Ginzburg-Landau, and Stokes equations. Additionally, we evaluate the performance of SSBFGS and SSBroyden for Deep Operator Network (DeepONet) architectures, demonstrating their effectiveness for data-driven operator learning. Our findings provide state-of-the-art results with orders-of-magnitude accuracy improvements without the use of adaptive weights or any other enhancements typically employed in PINNs. More broadly, our work reveal insights into the effectiveness of quasi-Newton optimization strategies in significantly improving the convergence and accurate generalization of PINNs and PIKANs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗
DOE OSTI2024

O Corona, where art thou? eROSITA’s view of UV-optical-IR variability-selected massive black holes in low-mass galaxies

Finding massive black holes (MBHs,M BH ≈ 10 4 –10 7 M ⊙ ) in the nuclei of low-mass galaxies $\left( {{M_*}\mathop {\mathop < \limits_ }\limits_ {{10}^{10}}{M_ \odot }} \right)$ is crucial to constrain seeding and growth of black holes over cosmic time, but it is particularly challenging due to their low accretion luminosities. Variability selection via long-term photometric ultraviolet, optical, or infrared (UVOIR) light curves has proved effective and identifies lower-Eddington ratios compared to broad and narrow optical spectral lines searches. In the inefficient accretion regime, X-ray and radio searches are effective, but they have been limited to small samples. Therefore, differences between selection techniques have remained uncertain. Here, we present the first large systematic investigation of the X-ray properties of a sample of known MBH candidates in dwarf galaxies. We extracted X-ray photometry and spectra of a sample of ~200 UVOIR variability-selected MBHs and significantly detected 17 of them in the deepest available SRG/eROSITA image, of which four are newly discovered X-ray sources and two are new secure MBHs. This implies that tens to hundreds of LSST MBHs will have SRG/eROSITA counterparts, depending on the seeding model adopted. Surprisingly, the stacked X-ray images of the many non-detected MBHs are incompatible with standard disk-corona relations, typical of active galactic nuclei, inferred from both the optical and radio fluxes. They are instead compatible with the X-ray emission predicted for normal galaxies. After careful consideration of potential biases, we identified that this X-ray weakness needs a physical origin. A possibility is that a canonical X-ray corona might be lacking in the majority of this population of UVOIR-variability selected low-mass galaxies or that unusual accretion modes and spectral energy distributions are in place for MBHs in dwarf galaxies. This result reveals the potential for severe biases in occupation fractions derived from data from only one waveband combined with SEDs and scaling relations of more massive black holes and galaxies.

Astronomy & Astrophysics↗
DOE OSTI2026

Algebraic Multigrid with Filtering: An Efficient Preconditioner for Interior Point Methods in Large-Scale Contact Mechanics Optimization

Large-scale contact mechanics simulations are crucial in many engineering fields such as structural design and manufacturing. In the frictionless case, contact can be modeled by minimizing an energy functional; however, these problems are often nonlinear, nonconvex, and increasingly difficult to solve as mesh resolution increases. In this work, we employ a Newton-based interior-point (IP) filter line-search method, an effective approach for large-scale constrained optimization. While this method converges rapidly, each iteration requires solving a large saddle-point linear system that becomes ill-conditioned as the optimization process converges, largely due to IP treatment of the contact constraints. Such ill-conditioning can hinder solver scalability and increase iteration counts with mesh refinement. Here, to address this, we introduce a novel preconditioner, algebraic multigrid with filtering (AMGF), tailored to the Schur complement of the saddle-point system. Building on the classical AMG solver, commonly used for elasticity, we augment it with a specialized subspace correction that filters near null space components introduced by contact interface constraints. Through theoretical analysis and numerical experiments on a range of linear and nonlinear contact problems, we demonstrate that the proposed solver achieves mesh independent convergence and maintains robustness against the ill-conditioning that notoriously plagues IP methods. These results indicate that AMGF makes contact mechanics simulations more tractable and broadens the applicability of Newton-based IP methods in challenging engineering scenarios. More broadly, AMGF is well suited for problems, optimization or otherwise, where solver performance is limited by a low-dimensional subspace, such as those arising from localized constraints, interface conditions, or model heterogeneities. This makes the method widely applicable beyond contact mechanics and constrained optimization.

Mathematics and Computing↗
DOE OSTI2024

Spectroscopic search for optical emission lines from dark matter decay

We search for narrow-line optical emission from dark matter decay by stacking dark-sky spectra from the Dark Energy Spectroscopic Instrument (DESI) at the redshift of nearby galaxies from DESI’s Bright Galaxy and Luminous Red Galaxy samples. Our search uses regions separated by 5 to 20 arcsec from the centers of the galaxies, corresponding to an impact parameter of approximately 50 kpc. No unidentified spectral line shows up in the search, and we place a line flux limit of 10 −19 ergs/s/cm 2 /arcsec 2 on emissions in the wavelength range of 2000–9000⁢ Å. This places the tightest constraints yet on the two-photon decay of dark matter in the mass range of 5 to 12 eV, with a particle lifetime exceeding 3 × 10 25 s. Furthermore, this detection limit also implies that the line surface brightness contributed from all dark matter along the line of sight is at least 2 orders of magnitude lower than the measured extragalactic background light (EBL), ruling out the possibility that narrow optical-line emission from dark matter decay is a major source of the EBL.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Early results in the search for extreme coronal line emitters with the Dark Energy Spectroscopic Instrument

Here, we present the results of our search through the early data release (EDR) of the Dark Energy Spectroscopic Instrument (DESI) for extreme coronal line emitters (ECLEs) – a rare classification of galaxies displaying strong, high-ionization iron coronal emission lines within their spectra. With the requirement of a strong X-ray continuum to generate the coronal emission, ECLEs have been linked to both active galactic nuclei (AGNs) and tidal disruption events (TDEs). We focus our search on identifying TDE-linked ECLEs. We identify three such objects within the EDR sample, highlighting DESI’s effectiveness for discovering new nuclear transients, and determine a galaxy-normalized TDE-linked ECLE rate of RG=5 −3+5×10−6 galaxy−1 yr−1 at a median redshift of z=0.2 – broadly consistent with previous work. Additionally, we identify more than 200 AGNs displaying coronal emission lines, which serve as the primary astrophysical contaminants in searches for TDE-related events. We also include an outline of the custom python code developed for this search.

Clark, Peter [Southampton U.; Portsmouth U., ICG] ↗
DOE OSTI2024

Fixed lines in four fermion models in two dimensions

Motivated by conjectures about near-horizon dynamics in quantum gravity, we search for lines of perturbatively accessible fixed points emanating from models of N free fermions. Through two loops we find a new class of models, apart from the well-known Abelian Thirring models. Further study is needed to see whether these can lead to true conformal manifolds, or perhaps a new class of large- N fixed points. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

A simple introduction to the SiMPL method for density-based topology optimization

We introduce a novel method for solving density-based topology optimization problems: Sigmoidal Mirror descent with a Projected Latent variable (SiMPL). The SiMPL method (pronounced as “the simple method”) optimizes a design using only first-order derivative information of the objective function. The bound constraints on the density field are enforced with the help of the (negative) Fermi–Dirac entropy, which is also used to define a non-symmetric distance function called a Bregman divergence on the set of admissible designs. This Bregman divergence leads to a simple update rule that is further simplified with the help of a so-called latent variable. Because the SiMPL method involves discretizing the latent variable, it produces a sequence of pointwise-feasible iterates, even when high-order finite elements are used in the discretization. Numerical experiments demonstrate that the method outperforms other popular first-order optimization algorithms. In conclusion, to outline the general applicability of the technique, we include examples with (self-load) compliance minimization and compliant mechanism optimization problems.

Calculus of Variations and Optimization↗
DOE OSTI2024

Magnetothermopower of Nodal-Line Semimetals

The search for materials with large thermopower is of great practical interest. Dirac and Weyl semimetals have recently proven to exhibit superior thermoelectric properties, particularly when subjected to a quantizing magnetic field. Here, we consider whether a similar enhancement arises in nodal-line semimetals, for which the conduction and valence band meet at a line or ring in momentum space. We compute the Seebeck and Nernst coefficients for arbitrary temperature and magnetic field and we find a wealth of different scaling regimes. Most strikingly, when a sufficiently strong magnetic field is applied along the direction of a straight nodal line or in the plane of a nodal ring, the large degeneracy of states leads to a large linear-in- B thermopower that is temperature independent even at low temperatures. Our results suggest that nodal-line semimetals may offer significant opportunity for efficient low-temperature thermoelectrics. Published by the American Physical Society 2024

Chakraborty, Poulomi (ORCID:0000000213969473)↗
DOE OSTI2025

XRISM Constraints on Unidentified X-Ray Emission Lines, Including the 3.5 keV Line, in the Stacked Spectrum of 10 Galaxy Clusters

We stack 3.75 Ms of early XRISM Resolve observations of 10 galaxy clusters to search for unidentified spectral lines in the E = 2.5–15 keV band (rest frame), including the E = 3.5 keV line reported in earlier low spectral resolution studies of cluster samples. Such an emission line may originate from the decay of the sterile neutrino, a warm dark matter (DM) candidate. No unidentified lines are detected in our stacked cluster spectrum, with the 3σ upper limit on the m s ∼ 7.1 keV DM particle decay rate (which corresponds to an E = 3.55 keV emission line) of Γ ∼ 1.0 × 10 −27 s −1 . This upper limit is 3–4 times lower than the one derived by Hitomi Collaboration from the Perseus observation but still 5 times higher than the XMM-Newton detection reported by E. Bulbul et al. in the stacked cluster sample. XRISM Resolve, with its high spectral resolution but small field of view, may reach the sensitivity needed to test the XMM-Newton cluster sample detection by combining several years worth of future cluster observations.

Astronomy and AstroPhysics↗
DOE OSTI2025

Sensitivity of JWST to eV-Scale Decaying Axion Dark Matter

The recently launched James Webb Space Telescope can resolve eV-scale emission lines arising from dark matter decay. We forecast the end-of-mission sensitivity to the decay of axions, a leading dark matter candidate, in the Milky Way using the blank-sky observations expected during standard operations. Searching for unassociated emission lines will constrain axions in the mass range 0.18 to 2.6 eV with axion-photon couplings g a γ γ ≳ 5.5 × 10 − 12 GeV − 1 . In particular, these results will constrain nucleophobic QCD axions to masses ≲ 0.2 eV . Published by the American Physical Society 2025

Roy, Sandip (ORCID:0000000276387454)↗
DOE OSTI2026

A Search for Successful and Choked Jets in Nearby Broad-lined Type Ic Supernovae

The observational link between long gamma-ray bursts (GRBs) and broad-lined stripped-envelope core-collapse supernovae (SNe Ic-BL) is well established. Significant progress has been made in constraining what fraction of SNe Ic-BL may power high- or low-luminosity GRBs when viewed at small off-axis angles. However, the GRB–SN connection still lacks a complete understanding in the broader context of massive-star evolution and explosion physics. Models predict a continuum of outcomes for the fastest ejecta, from choked to ultrarelativistic jets, and observations from radio to X-rays are key to probing these scenarios across a range of viewing angles and velocities. Here, we present results from a coordinated radio-to-X-ray campaign targeting nearby (z ≲ 0.1) SNe Ic-BL designed to explore this diversity. With eight new radio-monitored events and updated data for one previously observed SN, we further tighten constraints on the fraction of SNe Ic-BL as relativistic as SN 1998bw (GRB 980425). We identify SN 2024rjw as a new radio-loud event likely powered by strong interaction with circumstellar material, and add evidence supporting a similar interpretation for SN 2020jqm. We also establish new limits on the properties of radio-emitting ejecta with velocities consistent with cocoons from choked jets, highlighting SN 2022xxf as a promising cocoon-dominated candidate. These results refine our understanding of the continuum linking ordinary SNe Ic-BL, engine-driven explosions, and GRBs, and contribute to building a sample that will inform future multimessenger searches for electromagnetic counterparts to high-energy neutrinos.

O’Dwyer, Tanner [Johns Hopkins Univ., Baltimore, M↗
DOE OSTI2024

Advancing Insights into Electrochemical Pre‐Treatments of Supported Nanoparticle Electrocatalysts by Combining a Design of Experiments Strategy with In Situ Characterization

Activation, break-in, and/or pre-treatment protocols are generally applied to energy conversion devices before regular operation to reach stable performance. There remains much to understand about the relationships among physical properties, performance, and electrochemical pre-treatments. Here, a design-of-experiments (DoE) strategy is employed to address this gap by demonstrating the influence of five pre-treatment parameters for carbon-supported Pt-nanoparticle catalysts on the electrocatalytic oxygen reduction reaction (ORR). A subset of pre-treatments, developed using a central composite design, are tested in a flow cell combined with an inductively-coupled plasma mass spectrometer (on-line ICP-MS). The DoE-based approach facilitates comprehensive insights from two orders of magnitude fewer experiments than a conventional grid search. The coupled on-line ICP-MS setup enables effective catalysis and real-time catalyst dissolution data. Leveraging insights from DoE for on-line ICP-MS and additional characterization, a model is built between the degradation of a multi-dimensional supported Pt surface, its performance, and applied electrochemical parameters. These investigations identify surface modifications, such as oxidation, and subsequent restructuring of Pt during pre-treatment as a primary cause of performance deterioration during ORR. By combining DoE with advanced characterization techniques, a powerful approach is demonstrated to gain a mechanistic understanding of pre-treatment protocols that can be broadly adapted to various reaction chemistries.

Platinum↗
DOE OSTI2025

Hunting Dark Matter Lines in the Infrared Background with the James Webb Space Telescope

Dark matter particles with a mass around 1 eV can decay into near-infrared photons. Utilising available public blank sky observations from the NIRSpec IFU on the James Webb Space Telescope (JWST), we search for a narrow emission line due to decaying dark matter and derive leading constraints in the mass range 0.8-3 eV on the decay rate to photons, and more specifcally, on the axion-photon coupling for the case of axion-like particles. We exclude τ < 3.5 x 10 26 s at m DM ≃ 0.8 eV and, in the case of axions, g aγγ > 1.3 x 10 -11 GeV -1 for m a ≃ 2.2 eV. Our results do not rely on dedicated observations, rather we use blank sky observations intended for sky subtraction, and thus our reach may be automatically strengthened as JWST continues to observe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

X-ray characterization of fully-depleted p-channel Skipper-CCDs for the DarkNESS mission

The Dark matter Nanosatellite Equipped with Skipper Sensors (DarkNESS) mission is a 6U CubeSat designed to search for X-ray lines from decaying dark matter using Skipper-CCDs. Thick, fully-depleted p-channel Skipper-CCDs provide low readout noise and high quantum efficiency for 1-10 keV X-rays, but their X-ray performance has not yet been demonstrated in the space environment. DarkNESS will operate in low-Earth orbit, where trapped protons induce displacement damage in the sensor that increases charge-transfer inefficiency and degrades the X-ray energy resolution. This work measures the X-ray line response of Skipper-CCDs before and after proton irradiation and quantifies the associated degradation. A sensor was exposed to 217 MeV protons at a fluence of 8.4 x 10^10 protons cm^-2, corresponding to a displacement-damage dose more than an order of magnitude above the three-year expectation for representative mid-inclination and Sun-synchronous low-Earth orbits. A 55Fe source was used to compare the energy resolution of the beam-exposed quadrant to adjacent unexposed quadrants and a non-irradiated reference sensor. These measurements provide a quantitative assessment of radiation-induced spectral degradation in Skipper-CCDs and enable an estimate of the end-of-life X-ray energy resolution expected for DarkNESS operation in low-Earth orbit.

Alpine, Phoenix [Illinois U., Urbana] (ORCID:00090↗
DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗