Search NASA⌕ Search

SEARCH · Search NASA

Results for “quadratic”

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 109 records · Page 6

Leading nonlinear tidal effects and scattering amplitudes

We present the two-body Hamiltonian and associated eikonal phase, to leading post-Minkowskian order, for infinitely many tidal deformations described by operators with arbitrary powers of the curvature tensor. Scattering amplitudes in momentum and position space provide systematic complementary approaches. For the tidal operators quadratic in curvature, which describe the linear response to an external gravitational field, we work out the leading post-Minkowskian contributions using a basis of operators with arbitrary numbers of derivatives which are in one-to-one correspondence with the worldline multipole operators. Explicit examples are used to show that the same techniques apply to both bodies interacting tidally with a spinning particle, for which we find the leading contributions from quadratic in curvature tidal operators with an arbitrary number of derivatives, and to effective field theory extensions of general relativity. We also note that the leading post-Minkowskian order contributions from higher-dimension operators manifest double-copy relations. Finally, we comment on the structure of higher-order corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the planar limit of 3d $$ {\mathrm{T}}_{\rho}^{\sigma}\left[\mathrm{SU}\left(\mathrm{N}\right)\right] $$

A bstract We discuss a limit of 3d $$ {T}_{\rho}^{\sigma}\left[\mathrm{SU}(N)\right] $$ T ρ σ SU N quiver gauge theories in which the number of nodes is large and the ranks scale quadratically with the length of the quiver. The sphere free energies and topologically twisted indices are obtained using supersymmetric localization. Both scale quartically with the length of the quiver and quadratically with N , with trilogarithm functions depending on the quiver data as coefficients. The IR SCFTs have well-behaved supergravity duals in Type IIB, and the free energies match precisely with holographic results. Previously discussed theories with N 2 ln N scaling arise as limiting cases. Each balanced 3d quiver theory is linked to a 5d parent, whose matrix model is related and dominated by the same saddle point, leading to close relations between BPS observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A more natural composite Higgs model

Composite Higgs models provide an attractive solution to the hierarchy problem. However, many realistic models suffer from tuning problems in the Higgs potential. There are often large contributions from the UV dynamics of the composite resonances to the Higgs potential, and tuning between the quadratic term and the quartic term is required to separate the electroweak breaking scale and the compositeness scale. We consider a composite Higgs model based on the SU(6)/Sp(6) coset, where an enhanced symmetry on the fermion resonances can minimize the Higgs quadratic term. Moreover, a Higgs quartic term from the collective symmetry breaking of the little Higgs mechanism can be realized by the partial compositeness couplings between elementary Standard Model fermions and the composite operators, without introducing new elementary fields beyond the Standard Model and the composite sector. The model contains two Higgs doublets, as well as several additional pseudo-Nambu-Goldstone bosons. To avoid tuning, the extra Higgs bosons are expected to be relatively light and may be probed in the future LHC runs. The deviations of the Higgs couplings and the weak gauge boson couplings also provide important tests as they are expected to be close to the current limits in this model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

One-loop matrix elements of effective superstring interactions: α' -expanding loop integrands

In the low-energy effective action of string theories, non-abelian gauge interactions and supergravity are augmented by infinite towers of higher-mass-dimension operators. We propose a new method to construct one-loop matrix elements with insertions of operators D 2k F n and D 2k R n in the tree-level effective action of type-I and type-II superstrings. Inspired by ambitwistor string theories, our method is based on forward limits of moduli-space integrals using string tree-level amplitudes with two extra points, expanded in powers of the inverse string tension α'. Similar to one-loop ambitwistor computations, intermediate steps feature non-standard linearized Feynman propagators which eventually recombine to conventional quadratic propagators. With linearized propagators the loop integrand of the matrix elements obey one-loop versions of the monodromy and KLT relations. We express a variety of four- and five-point examples in terms of quadratic propagators and formulate a criterion on the underlying genus-one correlation functions that should make this recombination possible at all orders in α'. The ultraviolet divergences of the one-loop matrix elements are crosschecked against the non-separating degeneration of genus-one integrals in string amplitudes. Conversely, our results can be used as a constructive method to determine degenerations of elliptic multiple zeta values and modular graph forms at arbitrary weight.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

On the emerging potential of quantum annealing hardware for combinatorial optimization

Abstract Over the past decade, the usefulness of quantum annealing hardware for combinatorial optimization has been the subject of much debate. Thus far, experimental benchmarking studies have indicated that quantum annealing hardware does not provide an irrefutable performance gain over state-of-the-art optimization methods. However, as this hardware continues to evolve, each new iteration brings improved performance and warrants further benchmarking. To that end, this work conducts an optimization performance assessment of D-Wave Systems’ Advantage Performance Update computer, which can natively solve sparse unconstrained quadratic optimization problems with over 5,000 binary decision variables and 40,000 quadratic terms. We demonstrate that classes of contrived problems exist where this quantum annealer can provide run time benefits over a collection of established classical solution methods that represent the current state-of-the-art for benchmarking quantum annealing hardware. Although this work does not present strong evidence of an irrefutable performance benefit for this emerging optimization technology, it does exhibit encouraging progress, signaling the potential impacts on practical optimization tasks in the future.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Optical response of a dual membrane active–passive optomechanical cavity

Highlights: • The optical response is determined by competition between photon-tunneling, two-phonon process, gain effect. • The photon-tunneling adjusts supermode frequency. • The absorption coefficient can exceed 1. • Not only realize fast/slow light but conversion between fast and slow light. We investigate a dual membrane active–passive cavity where each mechanical membrane individually quadratically coupled to passive and active cavities via two-phonon process. Due to the fact that in the quadratically coupled optomechanical system mean-field approximation fails, hence to analyze the system completely, we switch to a more generalized out of equilibrium approach, namely Keldysh Green’s functional approach. We calculate transmission rate using predetermined full retarded Green’s function, and then numerically examine the effect of the various parameters on the transmission coefficient and discuss the features and physical mechanism behind them in detail. On the basis of the optical responsivity we further extend our study to fast and slow light phenomenon. The results show that our proposed system can not only realize ultra-fast/ultra-slow light under proper choice of cavity parameters but also a realization of the conversion between fast and slow light and vice versa has also been pointed out.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A fourth-order Lagrangian discontinuous Galerkin method using a hierarchical orthogonal basis on curvilinear grids

The existing high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic methods are restricted to using quadratic meshes with quadratic polynomials (P2), which in turn, yield up to third-order accuracy. Additionally, these existing DG hydrodynamic schemes, when extended to work with cubic meshes and cubic polynomials (P3), can be unstable on strong-shock problems. Therefore, this paper presents a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method to simulate compressible material dynamics ( e.g., gasses, fluids, and solids) with strong-shocks using cubic meshes and cubic polynomials, and delivers up to fourth-order accuracy on smooth flows. The stability on shock problems is achieved using new hierarchical orthogonal basis functions and a new subcell mesh stabilization (SMS) scheme for cubic meshes. The accuracy and robustness of the new high-order accurate Lagrangian DG hydrodynamic method is demonstrated by simulating a diverse suite of challenging test problems covering gas and solid dynamic problems on curvilinear meshes.

97 MATHEMATICS AND COMPUTING↗

Parallelized POD-based suboptimal economic model predictive control of a state-constrained Boussinesq approximation

Motivated by an energy efficient building application, we want to optimize a quadratic cost functional subject to the Boussinesq approximation of the Navier-Stokes equations and to bilateral state and control constraints. Since the computation of such an optimal solution is numerically costly, we design an efficient strategy to compute a sub-optimal (but applicationally acceptable) solution with significantly reduced computational effort. We employ an economic Model Predictive Control (MPC) strategy to obtain a feedback control. The MPC sub-problems are based on a linear-quadratic optimal control problem subjected to mixed control and state constraints and a convection-diffusion equation, reduced with proper orthogonal decomposition. Finally, to solve each sub-problem, we apply a primal-dual active set strategy. The method can be fully parallelized, which enables the solution of large problems with real-world parameters.

97 MATHEMATICS AND COMPUTING↗

Comparative Evaluation of Control-Oriented Heavy Duty Vehicle Air Drag Coefficient Models

Heavy-duty vehicles (HDVs) are a significant source of fuel consumption and greenhouse gas emissions, prompting solutions such as HDV platooning to mitigate these negative impacts through air drag reduction. The intervehicle distance in an HDV platoon needs to be carefully selected, such that the platoon-level energy efficiency and safety considerations can be well balanced. Underlying this problem lies in accurately modeling the relationship between HDV air drag coefficient and intervehicle distance. Through comprehensive evaluation and comparison, we analyze five control-oriented HDV air drag coefficient models, including the polynomial model, rational polynomial model, rational model, semi-quadratic model, and ridge model. Leveraging Scipy Curve-Fit toolbox and our previously compiled air drag coefficient datasets, we optimally identify the parameters inside each model. The calibrated models are then thoroughly evaluated via five complementary metrics. The comparison results reveal that the semi-quadratic model has the highest overall performance, while the widely adopted rational model only exhibits suboptimal performance.

Best, Micah↗

Reduced-order modeling on a near-term quantum computer

Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reducedorder modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. This reformulation was chosen as a means of facilitating implementation on a near-term quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and α = 15°. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗

OC7 phase I: Toward practical sea-state-dependent modeling of hydrodynamic viscous drag and damping

Here, this article presents a collaborative research campaign under the OC7 project on refining the engineering modeling approach for hydrodynamic viscous drag and damping on floating wind platforms, focusing on the adjustment of hydrodynamic drag and damping coefficients for different sea states. The participant simulation results show significant improvements over the previous OC6 project in predicting the low-frequency resonance motion under nonoperational conditions. The improvements are mainly due to enhanced modeling, including the adoption of wave stretching, and directly tuning the coefficients to measured platform motion in waves instead of free decay. For accurate predictions of mean- and slow-drift motion, the better performing models use a decreasing column splash zone drag coefficient and increasing surge damping/drag with increasing wave height. The model tuning for heave and pitch resonance shows less consistency. Generally, both quadratic drag and additional heave or pitch damping are needed for accurate predictions. Alternatively, a quadratic drag formulation with velocity filtering for the rectangular pontoons leads to improved predictions without additional damping. This model is also potentially more predictive, requiring minimal adjustment to its parameters for different conditions.

17 WIND ENERGY↗

Second-order wave excitation forces in WEC-Sim/MOST: Implementation, experimental validation, and code-to-code comparison

Accurate prediction of second-order hydrodynamic loads is essential for floating bodies, including floating offshore wind turbines, wave energy converters, and hybrid wind–wave platforms. These nonlinear effects, arising from both sum- and difference-frequency forcing, are critical for capturing key response characteristics but remain challenging to model efficiently. In this work, we extend the open-source Wave Energy Converter Simulator / MATLAB for Offshore Simulation Tool by implementing second-order wave excitation forces, supporting both the full Quadratic Transfer Function formulation and the Newman approximation. The full Quadratic Transfer Function method is used for all code-to-code comparisons and experimental validation, while the Newman approximation is provided as a computationally lighter alternative. To benchmark the new capability, we perform a code-to-code comparison with OpenFAST and OrcaFlex. We then validate the enhanced model using wave-tank measurements of a 1:96 scale DeepCwind semi-submersible, showing that second-order effects are required to reproduce platform motions. The implementation employs a computationally efficient pre-computation strategy for second-order wave excitation forces, reducing simulation cost while maintaining engineering accuracy. Overall, this work advances the tool as an open-source and versatile tool for modelling floating offshore renewable-energy systems requiring second-order hydrodynamic fidelity.

17 WIND ENERGY↗

On relaxations of the max k -cut problem formulations

Here, a tight continuous relaxation is a crucial factor in solving mixed integer formulations of many NP-hard combinatorial optimization problems. The (weighted) max k-cut problem is a fundamental combinatorial optimization problem with multiple notorious mixed integer optimization formulations. In this paper, we explore four existing mixed integer optimization formulations of the max k-cut problem. Specifically, we show that the continuous relaxation of a binary quadratic optimization formulation of the problem is: (i) stronger than the continuous relaxation of two mixed integer linear optimization formulations and (ii) at least as strong as the continuous relaxation of a mixed integer semidefinite optimization formulation. We also conduct a set of experiments on multiple sets of instances of the max k-cut problem using state-of-the-art solvers that empirically confirm the theoretical results in item (i). Furthermore, these numerical results illustrate the advances in the efficiency of global non-convex quadratic optimization solvers and more general mixed integer nonlinear optimization solvers. As a result, these solvers provide a promising option to solve combinatorial optimization problems. Our codes and data are available on GitHub.

97 MATHEMATICS AND COMPUTING↗

Lift & Learn: Physics-informed machine learning for large-scale nonlinear dynamical systems

In this work, we present Lift & Learn, a physics-informed method for learning low-dimensional models for large-scale dynamical systems. The method exploits knowledge of a system’s governing equations to identify a coordinate transformation in which the system dynamics have quadratic structure. This transformation is called a lifting map because it often adds auxiliary variables to the system state. The lifting map is applied to data obtained by evaluating a model for the original nonlinear system. This lifted data is projected onto its leading principal components, and low-dimensional linear and quadratic matrix operators are fit to the lifted reduced data using a least-squares operator inference procedure. Analysis of our method shows that the Lift & Learn models are able to capture the system physics in the lifted coordinates at least as accurately as traditional intrusive model reduction approaches. This preservation of system physics makes the Lift & Learn models robust to changes in inputs. Numerical experiments on the FitzHugh–Nagumo neuron activation model and the compressible Euler equations demonstrate the generalizability of our model.

97 MATHEMATICS AND COMPUTING↗

Studying nuclear medium modification using the Gerasimov-Drell-Hearn sum rule

The Gerasimov–Drell–Hearn sum rule is a generic relation that has been used to make significant contributions to research in hadronic physics. It connects the spin-dependent cross-section for photoproduction off a particle to the squared ratio of the particle’s anomalous magnetic moment and its mass, (κ/M)2. Thus, for a nucleon embedded in a nucleus, the sum rule relates the cross-section to κ/M averaged quadratically over the nucleons comprising the nucleus. This quadratic averaging can be used to constrain the mechanism responsible for the medium modification of the nucleon. We also point out that the global properties of the embedded nucleon like its axial charge, mass or magnetic moment are observables measurable through sum rules.

Deur, A. [Thomas Jefferson National Accelerator Fa↗

Reynolds-number scaling of wall-pressure–velocity correlations in wall-bounded turbulence

Wall-pressure fluctuations are a practically robust input for real-time control systems aimed at modifying wall-bounded turbulence. The scaling behaviour of the wall-pressure–velocity coupling requires investigation to properly design a controller with such input data so that it can actuate upon the desired turbulent structures. A comprehensive database from direct numerical simulations (DNS) of turbulent channel flow is used for this purpose, spanning a Reynolds-number range$Re_\tau \approx 550\unicode{x2013}5200$. Spectral analysis reveals that the streamwise velocity is most strongly coupled to the linear term of the wall pressure, at a Reynolds-number invariant distance-from-the-wall scaling of$\lambda _x/y \approx 14$(and$\lambda _x/y \approx 8$for the wall-normal velocity). When extending the analysis to both homogeneous directions in$x$and$y$, the peak coherence is centred at$\lambda _x/\lambda _z \approx 2$and$\lambda _x/\lambda _z \approx 1$for$p_w$and$u$, and$p_w$and$v$, respectively. A stronger coherence is retrieved when the quadratic term of the wall pressure is concerned, but there is only little evidence for a wall-attached-eddy type of scaling. An experimental dataset comprising simultaneous measurements of wall pressure and velocity complements the DNS-based findings at one value of$Re_\tau \approx 2$k, with ample evidence that the DNS-inferred correlations can be replicated with experimental pressure data subject to significant levels of (acoustic) facility noise. It is furthermore shown that velocity-state estimations can be achieved with good accuracy by including both the linear and quadratic terms of the wall pressure. An accuracy of up to 72 % in the binary state of the streamwise velocity fluctuations in the logarithmic region is achieved; this corresponds to a correlation coefficient of$\approx$0.6. This thus demonstrates that wall-pressure sensing for velocity-state estimation – e.g. for use in real-time control of wall-bounded turbulence – has merit in terms of its realization at a range of Reynolds numbers.

Mechanics↗

On the convergence of an IEQ-based first-order semi-discrete scheme for the Beris-Edwards system

We present a convergence analysis of an unconditionally energy-stable first-order semi-discrete numerical scheme designed for a hydrodynamic Q-tensor model, the so-called Beris-Edwards system, based on the Invariant Energy Quadratization Method (IEQ). The model consists of the Navier–Stokes equations for the fluid flow, coupled to the Q-tensor gradient flow describing the liquid crystal molecule alignment. By using the Invariant Energy Quadratization Method, we obtain a linearly implicit scheme, accelerating the computational speed. However, this introduces an auxiliary variable to replace the bulk potential energy and it is a priori unclear whether the reformulated system is equivalent to the Beris-Edward system. In this work, we prove stability properties of the scheme and show its convergence to a weak solution of the coupled liquid crystal system. We also demonstrate the equivalence of the reformulated and original systems in the weak sense.

Weber, Franziska↗