Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained optimization”

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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

On the constrained Chebyshev approximation problem on ellipses

Constrained Chebyshev approximation problems of the type with minimum (p is an element of Pi(sub n):p(c)=1) and maximum (z is an element of E) with /p(z)/ are considered. Here Pi(sub n) denotes the set of all complex polynomials of degree at most n, E is any ellipse in the complex plane, and c is an element of C/E. Such approximation problems arise in the context of optimizing semi-iterative methods for the solution of large, sparse systems of linear equations Ax=b with complex non-Hermitian coefficient matrices A. The problem of obtaining optimal polynomial preconditioners for conjugate gradient type methods for Ax=b also leads to problems of this type. A new family of polynomials -- q(sub n)(z;c), n is an element of N, and c is an element of C/E -- are introduced as the polynomials which are optimal for a modified version of the Chebyshev approximation problem with Pi(sub n) replaced by a certain subfamily. Some simple properties of q(sub n) are also listed. A necessary and sufficient condition for q(sub n) to be the extremal polynomial for the approximation problem is then derived. Finally, it is shown that q(sub n) is indeed optimal for the problem for all fixed n whenever the distance between c and E is sufficiently large. Results of some numerical tests are presented.

Freund, Roland↗
NASA NTRS2021

A sampling-based optimization approach to handling environmental uncertainty for a planetary lander

Planning for unknown environments presents a number of technical challenges. The planner must ensure robustness to unknown phenomena and manage unpredictable variation in execution, all while operating in a capacity that maximizes its objective. Productivity in the face of these challenges re-quires an integrated approach to planning and execution that is capable of accomplishing goals, reacting to variation, and maximizing overall utility. We examine this problem in the context of a Europa Lander concept mission. We model the problem as a hierarchical task network, framing it as a utility maximization problem constrained on a depletable energy resource. We propose an uncertainty–sensitive deterministic planning framework that utilizes periodic replanning to better handle model uncertainty and variable execution. We demonstrate the efficacy of our framework through simulations of a Europa Lander concept mission in which our algorithm out-performs several baseline approaches in both utility maximization and robustness

Zilberstein, Shlomo↗
NASA NTRS2021

A Sampling-Based Optimization Approach to Handling Environmental Uncertainty for a Planetary Lander

Planning for unknown environments presents a number of technical challenges. The planner must ensure robustness to unknown phenomena and manage unpredictable variation in execution, all while operating in a capacity that maximizes its objective. Productivity in the face of these challenges re-quires an integrated approach to planning and execution that is capable of accomplishing goals, reacting to variation, and maximizing overall utility. We examine this problem in the context of a Europa Lander concept mission. We model the problem as a hierarchical task network, framing it as a utility maximization problem constrained on a depletable energy resource. We propose an uncertainty–sensitive deterministic planning framework that utilizes periodic replanning to better handle model uncertainty and variable execution. We demonstrate the efficacy of our framework through simulations of a Europa Lander concept mission in which our algorithm out-performs several baseline approaches in both utility maximization and robustness

Chien, Steve↗
DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Characterization and Optimization of Skipper CCDs for the SOAR Integral Field Spectrograph

We present results from the characterization and optimization of Skipper charge-coupled devices (CCDs) for use in a focal plane prototype for the Southern Astrophysical Research Integral Field Spectrograph (SIFS). We tested eight Skipper CCDs and selected six for SIFS based on performance results. The Skipper CCDs are 6k × 1k, 15 μm pixels, thick, fully depleted, p-channel devices that have been thinned to ~250 μm, backside processed, and treated with an anti-reflective coating. We demonstrate a single-sample readout noise of <4.3 e - rms pixel -1 in all amplifiers. We optimize the readout sequence timing to achieve a readout noise of 0.5 e - rms pixel -1 after 74 non-destructive measurements, which can be accomplished in a region covering 5% of the detector area in a readout time of <4 minutes. We demonstrate single-photon-counting in all 24 amplifiers (four amplifiers per detector) with a readnoise of σ N ~ 0.18 e - rms pixel -1 after N samp = 400 samples, and we constrain the degree of nonlinearity to be ≲1% at low signal levels (0 e - to 50 e - ). Clock-induced charge (CIC) remains an important issue when the Skipper CCD is configured to provide a large full-well capacity. We achieve a CIC rate of <1.45 × 10 -3 e - pixel -1 frame -1 for a full-well capacity of ~900 e - , which increases to a CIC rate of ~3 e - pixel -1 frame -1 for full-well capacities ~40,000-65,000 e - . We also perform conventional CCD characterization measurements such as charge transfer inefficiency (3.44 × 10 -7 on average), dark current (~2 × 10 -4 e - pixel -1 s -1 ), photon transfer curves, cosmetic defects (<0.45% “bad” pixels), and charge diffusion (point-spread function < 7.5 μm) to verify that these properties are consistent with expectations from conventional p-channel CCDs used for astronomy. Furthermore, we provide the first measurements of the brighter-fatter effect and absolute quantum efficiency (≳80% between 450 and 980 nm; ≳90% between 600 and 900 nm) using Skipper CCDs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Search for the nonresonant production of Higgs boson pairs via gluon fusion and vector-boson fusion in the $b\bar{b}$τ + τ – final state in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for the nonresonant production of Higgs boson pairs in the HH → $b\overline{b}$$\tau^{+}\tau^{-}$ channel is performed using 140 fb -1 of proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. The analysis strategy is optimized to probe anomalous values of the Higgs boson self-coupling modifier $\kappa_{\tau}$ and of the quartic HHVV (v = W,Z) coupling modifier $\kappa$ 2⁢V . No significant excess above the expected background from Standard Model processes is observed. An observed (expected) upper limit μ HH < 5.9⁢(3.3) is set at 95% confidence-level on the Higgs boson pair production cross section normalized to its Standard Model prediction. The coupling modifiers are constrained to an observed (expected) 95% confidence interval of -3.1 <$\kappa_{\tau}$ < 9.0 (-2.5 < $\kappa_{\tau}$ < 9.3) and -0.5 <$\kappa$ 2⁢V < 2.7 (-0.2 <$\kappa$ 2⁢V < 2.4), assuming all other Higgs boson couplings are fixed to the Standard Model prediction. The results are also interpreted in the context of effective field theories via constraints on anomalous Higgs boson couplings and Higgs boson pair production cross sections assuming different kinematic benchmark scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Fast Adaptive Neural Control of Resonant Extraction at Fermilab

We present the development of a machine learning (ML) based regulation system for third-order resonant beam extraction in the Mu2e experiment at Fermilab. Classical and ML-based controllers have been optimized using semi-analytic simulations and evaluated in terms of regulation performance and training efficiency. We compare several controller architectures and discuss the integration of neural control into an adaptive framework. We also present progress on surrogate models that predict the controller response given a spill intensity and controller action history. To enable real-time deployment, we report progress on implementing low-latency, edge-based inference suitable for hardware-constrained environments. Our results demonstrate the feasibility and advantages of ML-based control in managing complex, time-varying physical systems, with broader implications for accelerator operations and other domains requiring fast, adaptive regulation.

Berlioz, Jose Rene [Fermilab]↗
DOE OSTI2024

Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument

The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^{2}$ over 5 yr to constrain the cosmic expansion history through precise measurements of baryon acoustic oscillations (BAO). The scientific program for DESI was evaluated during a 5 month survey validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar Milky Way Survey (MWS), Bright Galaxy Survey (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the 5 yr program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a One-Percent Survey conducted at the conclusion of SV covering 140 deg$^{2}$ using the final target selection algorithms with exposures of a depth typical of the main survey. The SV indicates that DESI will be able to complete the full 14,000 deg$^{2}$ program with spectroscopically confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval z < 1.1, 0.39% over the redshift interval 1.1 < z < 1.9, and 0.46% over the redshift interval 1.9 < z < 3.5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Using a surrogate-assisted Bayesian framework to calibrate the runoff-generation scheme in the Energy Exascale Earth System Model (E3SM) v1

Abstract. Runoff is a critical component of the terrestrial water cycle, and Earth system models (ESMs) are essential tools to study its spatiotemporal variability. Runoff schemes in ESMs typically include many parameters so that model calibration is necessary to improve the accuracy of simulated runoff. However, runoff calibration at a global scale is challenging because of the high computational cost and the lack of reliable observational datasets. In this study, we calibrated 11 runoff relevant parameters in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) using a surrogate-assisted Bayesian framework. First, the polynomial chaos expansion machinery with Bayesian compressed sensing is used to construct computationally inexpensive surrogate models for ELM-simulated runoff at 0.5∘ × 0.5∘ for 1991–2010. The error metric between the ELM simulations and the benchmark data is selected to construct the surrogates, which facilitates efficient calibration and avoids the more conventional, but challenging, construction of high-dimensional surrogates for the ELM simulated runoff. Second, the Sobol' index sensitivity analysis is performed using the surrogate models to identify the most sensitive parameters, and our results show that, in most regions, ELM-simulated runoff is strongly sensitive to 3 of the 11 uncertain parameters. Third, a Bayesian method is used to infer the optimal values of the most sensitive parameters using an observation-based global runoff dataset as the benchmark. Our results show that model performance is significantly improved with the inferred parameter values. Although the parametric uncertainty of simulated runoff is reduced after the parameter inference, it remains comparable to the multimodel ensemble uncertainty represented by the global hydrological models in ISMIP2a. Additionally, the annual global runoff trend during the simulation period is not well constrained by the inferred parameter values, suggesting the importance of including parametric uncertainty in future runoff projections.

58 GEOSCIENCES↗
DOE OSTI2026

Applications of Multibody for Everybody (M4E) in Marine Energy

Multibody for Everybody (M4E) is an open-source symbolic dynamics modeling framework designed to automate the derivation of equations of motion and simulation of constrained multibody systems using the joint coordinate formulation.

16 TIDAL AND WAVE POWER↗
NASA NTRS1991

Red light-regulated growth. I. Changes in the abundance of indoleacetic acid and a 22-kilodalton auxin-binding protein in the maize mesocotyl

We examined the changes in the levels of indoleacetic acid (IAA), IAA esters, and a 22-kilodalton subunit auxin-binding protein (ABP1) in apical mesocotyl tissue of maize (Zea mays L.) during continuous red light (R) irradiation. These changes were compared with the kinetics of R-induced growth inhibition in the same tissue. Upon the onset of continuous irradiation, growth decreased in a continuous manner following a brief lag period. The decrease in growth continued for 5 hours, then remained constant at 25% of the dark rate. The abundance of ABP1 and the level of free IAA both decreased in the mesocotyl. Only the kinetics of the decrease in IAA within the apical mesocotyl correlated with the initial change in growth, although growth continued to decrease even after IAA content reached its final level, 50% of the dark control. This decrease in IAA within the mesocotyl probably occurs primarily by a change in its transport within the shoot since auxin applied as a pulse move basipetally in R-irradiated tissue at the same rate but with half the area as dark control tissue. In situ localization of auxin in etiolated maize shoots revealed that R-irradiated shoots contained less auxin in the epidermis than the dark controls. Irradiated mesocotyl grew 50% less than the dark controls even when incubated in an optimal level of auxin. However, irradiated and dark tissue contained essentially the same amount of radioactivity after incubation in [14C]IAA indicating that the light treatment does not affect the uptake into the tissue through the cut end, although it is possible that a small subset of cells within the mesocotyl is affected. These observations support the hypothesis that R causes a decrease in the level of auxin in epidermal cells of the mesocotyl, consequently constraining the growth of the entire mesocotyl.

Non-NASA Center↗
NASA NTRS2005

Single-Scale Retinex Using Digital Signal Processors

The Retinex is an image enhancement algorithm that improves the brightness, contrast and sharpness of an image. It performs a non-linear spatial/spectral transform that provides simultaneous dynamic range compression and color constancy. It has been used for a wide variety of applications ranging from aviation safety to general purpose photography. Many potential applications require the use of Retinex processing at video frame rates. This is difficult to achieve with general purpose processors because the algorithm contains a large number of complex computations and data transfers. In addition, many of these applications also constrain the potential architectures to embedded processors to save power, weight and cost. Thus we have focused on digital signal processors (DSPs) and field programmable gate arrays (FPGAs) as potential solutions for real-time Retinex processing. In previous efforts we attained a 21 (full) frame per second (fps) processing rate for the single-scale monochromatic Retinex with a TMS320C6711 DSP operating at 150 MHz. This was achieved after several significant code improvements and optimizations. Since then we have migrated our design to the slightly more powerful TMS320C6713 DSP and the fixed point TMS320DM642 DSP. In this paper we briefly discuss the Retinex algorithm, the performance of the algorithm executing on the TMS320C6713 and the TMS320DM642, and compare the results with the TMS320C6711.

Hines, Glenn↗
NASA NTRS2012

Deriving Plasma Densities and Elemental Abundances from SERTS Differential Emission Measure Analysis

We use high-resolution spectral emission line data obtained by the SERTS instrument during three rocket flights to demonstrate a new approach for constraining electron densities of solar active region plasma.We apply differential emission measure (DEM) forward-fitting techniques to characterize the multithermal solar plasma producing the observed EUV spectra, with constraints on the high-temperature plasma from the Yohkoh Soft X-ray Telescope. In this iterative process, we compare line intensities predicted by an input source distribution to observed line intensities for multiple iron ion species, and search a broad range of densities to optimize chi-square simultaneously for the many available density-sensitive lines. This produces a density weighted by the DEM, which appears to be useful for characterizing the bulk of the emitting plasma over a significant range of temperature. This "DEM-weighted density" technique is complementary to the use of density-sensitive line ratios and less affected by uncertainties in atomic data and ionization fraction for any specific line. Once the DEM shape and the DEM-weighted density have been established from the iron lines, the relative elemental abundances can be determined for other lines in the spectrum. We have also identified spectral lines in the SERTS wavelength range that may be problematic

corona↗
DOE OSTI2026

Quantum Routing and Entanglement Dynamics Through Bottlenecks

To implement arbitrary quantum circuits in architectures with restricted interactions, one may effectively simulate all-to-all connectivity by routing quantum information. We consider the entanglement dynamics and routing between two regions only connected through an intermediate “bottleneck” region with few qubits. In such systems, where the entanglement rate is restricted by a vertex boundary rather than an edge boundary of the underlying interaction graph, existing results such as the small incremental entangling theorem give only a trivial constant lower bound on the routing time (the minimum time to perform an arbitrary permutation). We significantly improve the lower bound on the routing time in systems with a vertex bottleneck. Specifically, for any system with two regions 𝐿,𝑅 with 𝑁 𝐿 ,𝑁 𝑅 qubits, respectively, coupled only through an intermediate region 𝐶 with 𝑁 𝐶 qubits, for any 𝛿 > 0 we show a lower bound of Ω⁢(𝑁$^{1−𝛿}_{𝑅}$/√𝑁 𝐿⁢ 𝑁 𝐶 ) on the Hamiltonian quantum routing time when using piecewise time-independent Hamiltonians, or time-dependent Hamiltonians subject to a smoothness condition. We also prove an upper bound on the average amount of bipartite entanglement between 𝐿 and 𝐶,𝑅 that can be generated in time 𝑡 by such architecture-respecting Hamiltonians in systems constrained by vertex bottlenecks, improving the scaling in the system size from 𝑂⁡(𝑁 𝐿⁢ 𝑡) to 𝑂⁡(√𝑁 𝐿⁢ 𝑡). As a special case, when applied to the star graph (i.e., one vertex connected to 𝑁 leaves), we obtain an Ω⁡(√𝑁 1−𝛿 ) lower bound on the routing time and on the time to prepare 𝑁/2 Bell pairs between the vertices. We also show that, in systems of free particles, we can route optimally on the star graph in time Θ⁡(√𝑁) using Hamiltonian quantum routing, obtaining a speedup over gate-based routing, which takes time Θ⁡(𝑁).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Data-Driven Compression of Electron-Phonon Interactions

First-principles calculations of electron interactions in materials have seen rapid progress in recent years, with electron-phonon ( e − ph ) interactions being a prime example. However, these techniques use large matrices encoding the interactions on dense momentum grids, which reduces computational efficiency and obscures interpretability. For e − ph interactions, existing interpolation techniques leverage locality in real space, but the high dimensionality of the data remains a bottleneck to balance cost and accuracy. Here we show an efficient way to compress e − ph interactions based on singular value decomposition (SVD), a widely used matrix and image compression technique. Leveraging (un)constrained SVD methods, we accurately predict material properties related to e − ph interactions—including charge mobility, spin relaxation times, band renormalization, and superconducting critical temperature—while using only a small fraction (1%–2%) of the interaction data. These findings unveil the hidden low-dimensional nature of e − ph interactions. Furthermore, they accelerate state-of-the-art first-principles e − ph calculations by about 2 orders of magnitude without sacrificing accuracy. Our Pareto-optimal parametrization of e − ph interactions can be readily generalized to electron-electron and electron-defect interactions, as well as to other couplings, advancing quantitative studies of condensed matter. Published by the American Physical Society 2024

Physics↗
NASA NTRS2016

Constraining the Sulfur Dioxide Degassing Flux from Turrialba Volcano, Costa Rica Using Unmanned Aerial System Measurements

Observed sulfur dioxide (SO2)mixing ratios onboard unmanned aerial systems (UAS) duringMarch 11-13, 2013 are used to constrain the three-day averaged SO2 degassing flux fromTurrialba volcanowithin a Bayesian inverse modeling framework. A mesoscale model coupled with Lagrangian stochastic particle backward trajectories is used to quantify the source-receptor relationships at very high spatial resolutions (i.e., b1 km). The model shows better performance in reproducing the near-surface meteorological properties and observed SO2 variations when using a first-order closure non-local planetary boundary layer (PBL) scheme. The optimized SO2 degassing fluxes vary from 0.59 +/- 0.37 to 0.83 +/- 0.33 kt d−1 depending on the PBL scheme used. These fluxes are in good agreement with ground-based gas flux measurements, and correspond to corrective scale factors of 8-12 to the posteruptive SO2 degassing rate in the AeroCom emission inventory. The maximum a posteriori solution for the SO2 flux is highly sensitive to the specification of prior and observational errors, and relatively insensitive to the SO2 loss term and temporal averaging of observations. Our results indicate relatively low degassing activity but sustained sulfur emissions from Turrialba volcano to the troposphere during March 2013. This study demonstrates the utility of low-cost small UAS platforms for volcanic gas composition and flux analysis.

Xi, Xin↗
DOE OSTI2025

Probing the Higgs boson CP properties in vector-boson fusion production in the H → τ + τ − channel with the ATLAS detector

The CP properties of the Higgs boson are studied in the vector-boson fusion production mode. The analysis exploits the decay mode of the Higgs boson into two τ-leptons using 140 fb −1 of proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. Results are obtained using the Optimal Observable method. CP-violating interactions between the Higgs boson and electroweak gauge bosons are considered in the effective field theory framework, with the interaction strength described in the HISZ basis by $\tilde{d}$, and in the Warsaw basis by c H$\tilde{W}$ , c H$\tilde{B}$ , and c H$\tilde{W}$B . No deviations relative to the Standard Model are observed, and limits are obtained on the strength parameters. The $\tilde{d}$ parameter is constrained to the interval [−0.012, 0.044] at the 95% confidence level while c H$\tilde{W}$ is constrained to [−0.24, 0.83], when considering both linear and quadratic effects of physics beyond the Standard Model.

Hadron-Hadron Scattering↗
DOE OSTI2025

Optimizing lossy state preparation for quantum sensing using Hamiltonian engineering

One of the most prominent platforms for demonstrating quantum sensing below the standard quantum limit is the spinor Bose–Einstein condensate. While a quantum advantage using several tens of thousands of atoms has been demonstrated in this platform, it faces an important challenge: atom loss. Atom loss is a Markovian error process modeled by Lindblad jump operators, and a no-go theorem, which we also show here, states that the loss of atoms in all spin components reduces the quantum advantage to a constant factor. Here, we show that this no-go theorem can be circumvented if we constrain atom losses to a single spin component. Moreover, we show that in this case, the maximum quantum Fisher information with N atoms scales as N 3/2 , establishing that a scalable quantum advantage can be achieved despite atom loss. Although Lindblad jump operators are generally non-Hermitian and non-invertible, we use their Moore–Penrose inverse to develop a framework for constructing several states with this scaling of Fisher information in the presence of losses. We use Hamiltonian engineering with realistic Hamiltonians to develop experimental protocols for preparing these states. Finally, we discuss possible experimental techniques to constrain the losses to a single spin mode.

74 ATOMIC AND MOLECULAR PHYSICS↗