Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bounds”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Stability Analysis of Inline ZFP Compression for Floating-Point Data in Iterative Methods

Currently, the dominating constraint in many high performance computing applications is data capacity and bandwidth, in both internode communications and even moreso in intranode data motion. A new approach to address this limitation is to make use of data compression in the form of a compressed data array. Storing data in a compressed data array and converting to standard IEEE-754 types as needed during a computation can reduce the pressure on bandwidth and storage. However, repeated conversions (lossy compression and decompression) introduce additional approximation errors, which need to be shown to not significantly affect the simulation results. Here, we extend recent work that analyzed the error of a single use of compression and decompression of the ZFP compressed data array representation to the case of time-stepping and iterative schemes, where an advancement operator is repeatedly applied in addition to the conversions. We show that the accumulated error for iterative methods involving fixed-point and time evolving iterations is bounded under standard constraints. An upper bound is established on the number of additional iterations required for the convergence of stationary fixed-point iterations. An additional analysis of traditional forward and backward error of stationary iterative methods using ZFP compressed arrays is also presented. The results of several 1D, 2D, and 3D test problems are provided to demonstrate the correctness of the theoretical bounds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Polyhedral Relaxations for Optimal Pump Scheduling of Potable Water Distribution Networks

The classic pump scheduling or optimal water flow (OWF) problem for water distribution networks (WDNs) minimizes the cost of power consumption for a given WDN over a fixed time horizon. In its exact form, the OWF is a computationally challenging mixed-integer nonlinear program (MINLP). It is complicated by nonlinear equality constraints that model network physics, discrete variables that model operational controls, and intertemporal constraints that model changes to storage devices. To address the computational challenges of the OWF, this paper develops tight polyhedral relaxations of the original MINLP, derives novel valid inequalities (or cuts) using duality theory, and implements novel optimization-based bound tightening and cut generation procedures. The efficacy of each new method is rigorously evaluated by measuring empirical improvements in OWF primal and dual bounds over 45 literature instances. The evaluation suggests that our relaxation improvements, model strengthening techniques, and a thoughtfully selected polyhedral relaxation partitioning scheme can substantially improve OWF primal and dual bounds, especially when compared with similar relaxation-based techniques that do not leverage these new methods.

bound tightening↗
DOE OSTI2025

Metasurfaces with Multipolar Resonances and Enhanced Light–Matter Interaction

Metasurfaces, composed of engineered nanoantennas, enable unprecedented control over electromagnetic waves by leveraging multipolar resonances to tailor light–matter interactions. This review explores key physical mechanisms that govern their optical properties, including the role of multipolar resonances in shaping metasurface responses, the emergence of bound states in the continuum (BICs) that support high-quality factor modes, and the Purcell effect, which enhances spontaneous emission rates at the nanoscale. These effects collectively underpin the design of advanced photonic devices with tailored spectral, angular, and polarization-dependent properties. This review discusses recent advances in metasurfaces and applications based on them, highlighting research that employs full-wave numerical simulations, analytical and semi-analytic techniques, multipolar decomposition, nanofabrication, and experimental characterization to explore the interplay of multipolar resonances, bound and quasi-bound states, and enhanced light–matter interactions. A particular focus is given to metasurface-enhanced photodetectors, where structured nanoantennas improve light absorption, spectral selectivity, and quantum efficiency. By integrating metasurfaces with conventional photodetector architectures, it is possible to enhance responsivity, engineer photocarrier generation rates, and even enable functionalities such as polarization-sensitive detection. The interplay between multipolar resonances, BICs, and emission control mechanisms provides a unified framework for designing next-generation optoelectronic devices. This review consolidates recent progress in these areas, emphasizing the potential of metasurface-based approaches for high-performance sensing, imaging, and energy-harvesting applications.

Kerker effect↗
DOE OSTI2020

The C‐Terminal Domain of α‐Synuclein Confers Steric Stabilization on Synaptic Vesicle‐Like Surfaces

Abstract While α‐synuclein, an intrinsically disordered protein linked to Parkinson's disease, has been shown to associate with membrane organelles, its overall cellular function remains nebulous. α‐Synuclein binds to membranes through its amino‐terminal domain (first ≈100 residues), but there is no consensus on the biophysical function of the carboxyl‐terminal domain (last ≈40 residues) due, in part, to its lack of strong interaction partners and persisting intrinsic disorder even when membrane bound. Here, by directly applying force on α‐synuclein bound to spherical nanoparticle‐supported lipid bilayers (SSLBs) and tracking higher‐order structural changes through small‐angle X‐ray scattering, strong evidence is presented that α‐synuclein sterically stabilizes membrane surfaces through its carboxyl‐terminal domain. Full‐length α‐synuclein dramatically increases the critical osmotic pressure at which SSLBs cluster ( P C ≈ 1.3 × 10 5 Pa) compared to α‐synuclein without the carboxyl‐terminal domain ( P C ≈ 1.9 × 10 4 Pa) at physiological salt and temperature conditions. This clustering of α‐synuclein‐bound SSLBs is shown to be reversible and sensitive to monovalent/divalent salt, both features of grafted polyelectrolyte‐mediated steric stabilization. In elucidating the biophysical function of α‐synuclein in the framework of polymer science, it is demonstrated that the carboxyl‐terminal domain can potentially utilize its persisting intrinsic disorder to functionalize membrane surfaces.

Chung, Peter J.↗
DOE OSTI2021

Machine‐learning‐based construction of barrier functions and models for safe model predictive control

Abstract In this paper, we propose a control Lyapunov‐barrier function‐based model predictive control method utilizing a feed‐forward neural network specified control barrier function (CBF) and a recurrent neural network (RNN) predictive model to stabilize nonlinear processes with input constraints, and to guarantee that safety requirements are met for all times. The nonlinear system is first modeled using RNN techniques, and a CBF is characterized by constructing a feed‐forward neural network (FNN) model with unique structures and properties. The FNN model for the CBF is trained based on data samples collected from safe and unsafe operating regions, and the resulting FNN model is verified to demonstrate that the safety properties of the CBF are satisfied. Given sufficiently small bounded modeling errors for both the FNN and the RNN models, the proposed control system is able to guarantee closed‐loop stability while preventing the closed‐loop states from entering unsafe regions in state‐space under sample‐and‐hold control action implementation. We provide the theoretical analysis for bounded unsafe sets in state‐space, and demonstrate the effectiveness of the proposed control strategy using a nonlinear chemical process example with a bounded unsafe region.

Chen, Scarlett↗
DOE OSTI2024

Anderson acceleration with approximate calculations: Applications to scientific computing

Here we provide rigorous theoretical bounds for Anderson acceleration (AA) that allow for approximate calculations when applied to solve linear problems. We show that, when the approximate calculations satisfy the provided error bounds, the convergence of AA is maintained while the computational time could be reduced. We also provide computable heuristic quantities, guided by the theoretical error bounds, which can be used to automate the tuning of accuracy while performing approximate calculations. For linear problems, the use of heuristics to monitor the error introduced by approximate calculations, combined with the check on monotonicity of the residual, ensures the convergence of the numerical scheme within a prescribed residual tolerance. Motivated by the theoretical studies, we propose a reduced variant of AA, which consists in projecting the least-squares used to compute the Anderson mixing onto a subspace of reduced dimension. The dimensionality of this subspace adapts dynamically at each iteration as prescribed by the computable heuristic quantities. We numerically show and assess the performance of AA with approximate calculations on: (i) linear deterministic fixed-point iterations arising from the Richardson's scheme to solve linear systems with open-source benchmark matrices with various preconditioners and (ii) non-linear deterministic fixed-point iterations arising from non-linear time-dependent Boltzmann equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Preconditioned least‐squares Petrov–Galerkin reduced order models

Abstract In this article, we introduce a methodology for improving the accuracy and efficiency of reduced order models (ROMs) constructed using the least‐squares Petrov–Galerkin (LSPG) projection method through the introduction of preconditioning. Unlike prior related work, which focuses on preconditioning the linear systems arising within the ROM numerical solution procedure to improve linear solver performance, our approach leverages a preconditioning matrix directly within the minimization problem underlying the LSPG formulation. Applying preconditioning in this way has the potential to improve ROM accuracy for several reasons. First, preconditioning the LSPG formulation changes the norm defining the residual minimization, which can improve the residual‐based stability constant bounding the ROM solution's error. The incorporation of a preconditioner into the LSPG formulation can have the additional effect of scaling the components of the residual being minimized to make them roughly of the same magnitude, which can be beneficial when applying the LSPG method to problems with disparate scales (e.g., dimensional equations, multi‐physics problems). Importantly, we demonstrate that an “ideal preconditioned” LSPG ROM (a ROM in which the preconditioner is the inverse of the Jacobian of its corresponding full order model) emulates projection of the full order model solution increment onto the reduced basis. This quantity defines a lower bound on the error of a ROM solution for a given reduced basis. By designing preconditioners that approximate the Jacobian inverse—as is common in designing preconditioners for solving linear systems—it is possible to obtain a ROM whose error approaches this lower bound. The proposed approach is evaluated on several mechanical and thermo‐mechanical problems implemented within the Albany HPC code and run in the predictive regime, with prediction across material parameter space. We demonstrate numerically that the introduction of simple Jacobi, Gauss‐Seidel, and ILU preconditioners into the proper orthogonal decomposition/LSPG formulation reduces significantly the ROM solution error, the reduced Jacobian condition number, the number of nonlinear iterations required to reach convergence, and the wall time (thereby improving efficiency). Moreover, our numerical results reveal that the introduction of preconditioning can deliver a robust and accurate solution for test cases in which the unpreconditioned LSPG method fails to converge.

Lindsay, Payton↗
DOE OSTI2024

Ternary complexes of isopentenyl phosphate kinase from Thermococcus paralvinellae reveal molecular determinants of non‐natural substrate specificity

Abstract Isopentenyl phosphate kinases (IPKs) have recently garnered attention for their central role in biocatalytic “isoprenol pathways,” which seek to reduce the synthesis of the isoprenoid precursors to two enzymatic steps. Furthermore, the natural promiscuity of IPKs toward non‐natural alkyl‐monophosphates (alkyl‐Ps) as substrates has hinted at the isoprenol pathways' potential to access novel isoprenoids with potentially useful activities. However, only a handful of IPK crystal structures have been solved to date, and even fewer of these contain non‐natural substrates bound in the active site. The current study sought to elucidate additional ternary complexes bound to non‐natural substrates using the IPK homolog fromThermococcus paralvinellae(TcpIPK). Four such structures were solved, each bound to a different non‐natural alkyl‐P and the phosphoryl donor substrate/product adenosine triphosphate (ATP)/adenosine diphosphate (ADP). As expected, the quaternary, tertiary, and secondary structures ofTcpIPK closely resembled those of IPKs published previously, and kinetic analysis of a novel alkyl‐P substrate highlighted the potentially dramatic effects of altering the core scaffold of the natural substrate. Even more interesting, though, was the discovery of a trend correlating the position of two α helices in the active site with the magnitude of an IPK homolog's reaction rate for the natural reaction. Overall, the current structures ofTcpIPK highlight the importance of continued structural analysis of the IPKs to better understand and optimize their activity with both natural and non‐natural substrates.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Filtering adaptive output feedback control for multivariable nonlinear systems with mismatched uncertainties and unmodeled dynamics

This article synthesizes a filtering adaptive output feedback controller for multivariable nonlinear systems with mismatched uncertainties and unmodeled dynamics. The multivariable nonlinear systems under consideration have both matched and mismatched uncertainties, which satisfy the semiglobal Lipschitz condition. The unmodeled dynamics are bounded-input bounded-output stable. By adopting an estimation/cancellation strategy, a piecewise constant adaptive law drives the estimation error to zero at every time instant, which yields the adaptive parameters; a disturbance rejection control law is designed to compensate the nonlinear uncertainties within the bandwidth of low-pass filters. The matched uncertainties are cancelled directly by adopting their opposite in the control signal, while a dynamic inversion of the system is required to eliminate the effect of the mismatched uncertainties on the output. A feedforward control law is designed to track a given command. Since the virtual reference system defines the best performance that can be achieved by the closed-loop system, the uniform performance bounds are derived for the states and control signals via comparison. Both numerical and practical examples are provided to illustrate the effectiveness of the proposed filtering adaptive output feedback control architecture, comparisons with the model reference adaptive control demonstrates the superiority of the proposed control method.

filtering adaptive control, mismatched uncertainit↗
DOE OSTI2023

Veneziano variations: how unique are string amplitudes?

String theory offers an elegant and concrete realization of how to consistently couple states of arbitrarily high spin. But how unique is this construction? In this paper we derive a novel, multi-parameter family of four-point scattering amplitudes exhibiting i) polynomially bounded high-energy behavior and ii) exchange of an infinite tower of high-spin modes, albeit with a finite number of states at each resonance. These amplitudes take an infinite-product form and, depending on parameters, exhibit mass spectra that are either unbounded or bounded, thus corresponding to generalizations of the Veneziano and Coon amplitudes, respectively. For the bounded case, masses converge to an accumulation point, a peculiar feature seen in the Coon amplitude but more recently understood to arise naturally in string theory [1]. Importantly, our amplitudes contain free parameters allowing for the customization of the slope and offset of the spin-dependence in the Regge trajectory. We compute all partial waves for this multi-parameter class of amplitudes and identify unitary regions of parameter space. For the unbounded case, we apply similar methods to derive new deformations of the Veneziano and Virasoro-Shapiro amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Bootstrapping 2d ϕ 4 theory with Hamiltonian truncation data

We combine the methods of Hamiltonian Truncation and the recently proposed generalisation of the S-matrix bootstrap that includes local operators to determine the two-particle scattering amplitude and the two-particle form factor of the stress tensor at s > 0 in the 2d Φ 4 theory. We use the form factor of the stress tensor at s ≤ 0 and its spectral density computed using Lightcone Conformal Truncation (LCT), and inject them into the generalized S-matrix bootstrap set-up. The obtained results for the scattering amplitude and the form factor are fully reliable only in the elastic regime. We independently construct the “pure” S-matrix bootstrap bounds (bootstrap without including matrix elements of local operators), and find that the sinh-Gordon model and its analytic continuation the “staircase model” saturate these bounds. Surprisingly, the Φ 4 two-particle scattering amplitude also very nearly saturates these bounds, and moreover is extremely close to that of the sinh-Gordon/staircase model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

BPS fivebrane stars. Part I. Expectation values of observables

We study ensembles of 1/2-BPS bound states of fundamental strings and NS-fivebranes (NS5-F1 states) in the AdS decoupling limit. We revisit a solution corresponding to an ensemble average of these bound states, and find that the appropriate duality frame for describing the near-source structure is the T-dual NS5-P frame, where the bound state is a collection of momentum waves on the fivebranes. We find that the fivebranes are generically well-separated; this property results in the applicability of perturbative string theory. The geometry sourced by the typical microstate is not close to that of an extremal non-rotating black hole; instead the fivebranes occupy a ball whose radius is parametrically much larger than the “stretched horizon” scale of the corresponding black hole. These microstates are thus better characterized as BPS fivebrane stars than as small black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The geometry of the modular bootstrap

Abstract We explore the geometry behind the modular bootstrap and its image in the space of Taylor coefficients of the torus partition function. In the first part, we identify the geometry as an intersection of planes with the convex hull of moment curves onR + ⊗ℤ, with boundaries characterized by the total positivity of generalized Hankel matrices. We phrase the Hankel constraints as a semi-definite program, which has several advantages, such as the validity of bounds irrespective of spin truncation. We derive bounds on the gap, twist-gap, and the space of Taylor coefficients themselves. We find that if the gap is above$$ {\Delta }_{\textrm{gap}}^{\ast } $$ ∆ gap ∗ , where$$ \frac{c-1}{12}<{\Delta}_{\textrm{gap}}^{\ast }<\frac{c}{12} $$ c − 1 12 < Δ gap ∗ < c 12 , all coefficients become bounded on both sides and kinks develop in the space. In the second part, we propose an analytic method of imposing the integrality condition for the degeneracy number in the spinless bootstrap, which leads to a non-convex geometry. We find that even at very low derivative order this condition rules out regions otherwise allowed by bootstraps at high derivative order.

Physics↗
DOE OSTI2026

Testing real WIMPs with CTAO

We forecast the reach of the upcoming Cherenkov Telescope Array Observatory (CTAO) to the full set of real representations within the paradigm of minimal dark matter. We employ effective field theory techniques to compute the annihilation cross section and photon spectrum that results when fermionic dark matter is the neutral component of an arbitrary odd and real representation of SU(2), including the Sommerfeld enhancement, next-to-leading log resummation of the relevant electroweak effects, and the contribution from bound states. We also compute the corresponding signals for scalar dark matter, with the exception of the bound state contribution. Results are presented for all real representations from the ~ 3 TeV triplet (or wino), a 3 of SU(2), to the ~ 300 TeV tredecuplet, a 13 of SU(2) that is at the threshold of the unitarity bound. Using these results, we forecast that with 500 hrs of Galactic Center observations and assuming background systematics are controlled at the level of $\mathcal{(O)}$(1%), then should no signal emerge, CTAO could exclude all representations up to the 11 of SU(2) in even the most conservative models for the dark-matter density in the inner galaxy, in both the fermionic and scalar dark matter cases. Assuming the default CTAO configuration, the tredecuplet will marginally escape exclusion, although we outline steps that CTAO could take to test even this scenario. In summary, CTAO appears poised to make a definitive statement on whether real WIMPs constitute the dark matter of our universe.

Models for Dark Matter↗
DOE OSTI2020

The Higgs trilinear coupling and the scale of new physics

We consider modifications of the Higgs potential due to new physics at high energy scales. These upset delicate cancellations predicted by the Standard Model for processes involving Higgs bosons and longitudinal gauge bosons, and lead to a breakdown of the theory at high energies. We focus on modifications of the Higgs trilinear coupling and use the violation of tree-level unitarity as an estimate of the scale where the theory breaks down. We obtain a completely model-independent bound of $ \underset{\sim }{<} $ 13 TeV for an order-1 modification of the trilinear. We argue that this bound can be saturated only in fine-tuned models, and the scale of new physics is likely to be much lower. The most stringent bounds are obtained from amplitudes involving multiparticle states that are not conventional scattering states. Our results show that a future determination of the Higgs cubic coupling can point to a well-defined scale of new physics that can be targeted and explored at future colliders.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

New tools for dissecting the general $\mathrm{2HDM}$

Two Higgs doublet models (2HDM) provide the low energy effective theory (EFT) description in many well motivated extensions of the Standard Model. It is therefore relevant to study their properties, as well as the theoretical constraints on these models. In this article we concentrate on three relevant requirements for the validity of the 2HDM framework, namely the perturbative unitarity bounds, the bounded from below constraints, and the vacuum stability constraints. In this study, we concentrate on the most general renormalizable version of the 2HDM — without imposing any parity symmetry, which may be violated in many UV extensions. We derive novel analytical expressions that generalize those previously obtained in more restrictive scenarios to the most general case. We also discuss the phenomenological implications of these bounds, focusing on CP violation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Kink-antikink scattering in a quantum vacuum

We study kink-antikink scattering in the sine-Gordon model in the presence of interactions with an additional scalar field, ψ, that is in its quantum vacuum. In contrast to the classical scattering, now there is quantum radiation of ψ quanta and the kink-antikink may form bound states that resemble breathers of the sine-Gordon model. We quantify the rate of radiation and map the parameters for which bound states are formed. Even these bound states radiate and decay, and eventually there is a transition into long-lived oscillons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Positivity from Cosmological Correlators

Abstract Effective field theories in flat space and in anti-de Sitter space are constrained by causality and unitarity, often in the form of positivity bounds. Similar bounds have been harder to demonstrate in cosmological backgrounds, where the roles of unitarity and causality are more obscure. Fortunately, the expansion of the universe ensures that late-time cosmological correlators are effectively classical and the role of unitarity is played by classical statistical inequalities. For multi-field inflation, the resulting positivity constraints have long been known in terms of the Suyama-Yamaguchi inequality. In this paper, we demonstrate that similar statistical bounds imply nontrivial constraints for massive fields in the early universe. We show that any real anomalous dimensions for principal series fields in de Sitter space must be positive. We also derive a limit on the amplitude of particular oscillatory signals from inflation, including those arising in cosmological collider physics. Finally, we demonstrate that these constraints manifest themselves directly in the two-point statistics of matter and galaxies that will be measured in upcoming surveys.

Physics↗