Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite amplitude method”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Caravel: A C++ framework for the computation of multi-loop amplitudes with numerical unitarity

We present the first public version of Caravel, a C++17 framework for the computation of multi-loop scattering amplitudes in quantum field theory, based on the numerical unitarity method. Caravel is composed of modules for the D-dimensional decomposition of integrands of scattering amplitudes into master and surface terms, the computation of tree-level amplitudes in floating point or finite-field arithmetic, the numerical computation of one- and two-loop amplitudes in QCD and Einstein gravity, and functional reconstruction tools. Here, we provide programs that showcase Caravel's main functionalities and allow to compute selected one- and two-loop amplitudes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Two- and three-meson scattering amplitudes with physical quark masses from lattice QCD

We study systems of two and three mesons composed of pions and kaons at maximal isospin using four CLS ensembles with 𝑎 ≈ 0.063 fm, including one with approximately physical quark masses. Using the stochastic Laplacian-Heaviside method, we determine the energy spectrum of these systems including many levels in different momentum frames and irreducible representations. Using the relativistic two- and three-body finite-volume formalism, we constrain the two- and three-meson K matrices, including not only the leading 𝑠 wave, but also 𝑝 and 𝑑 waves. By solving the three-body integral equations, we determine, for the first time, the physical-point scattering amplitudes for 3⁢𝜋 + , 3⁢𝐾 + , 𝜋 + ⁢𝜋 + ⁢𝐾 + , and 𝐾 + ⁢𝐾 + ⁢𝜋 + systems. These are determined for total angular momentum 𝐽 𝑃 = 0 − , 1 + , and 2 − . We also obtain accurate results for 2⁢𝜋 + , 𝜋 + ⁢𝐾 + , and 2⁢𝐾 + phase shifts. We compare our results to chiral perturbation theory and to phenomenological fits.

FOS: Physical sciences↗
DOE OSTI2022

From force chains to nonclassical nonlinear dynamics in cemented granular materials

In this letter, we present evidence for a mechanism responsible for the nonclassical nonlinear dynamics observed in many cemented granular materials that are generally classified as mesoscopic nonlinear elastic materials. We demonstrate numerically that force chains are created within the complex grain-pore network of these materials when subjected to dynamic loading. The interface properties between grains along with the sharp and localized increase of the stress occurring at the grain-grain contacts leads to a reversible decrease of the elastic properties at macroscopic scale and peculiar effects on the propagation of elastic waves when grain boundary properties are appropriately considered. These effects are observed for relatively small amplitudes of the elastic waves, i.e., within tens of microstrain, and relatively large wavelengths, i.e., orders of magnitude larger than the material constituents. The mechanics are investigated numerically using the hybrid finite-discrete-element method and match those observed experimentally using nonlinear resonant ultrasound spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Time domain probabilistic seismic risk analysis using ground motion prediction equations of Fourier amplitude spectra

Modeling of Fourier amplitude spectra (FAS) of seismic motions has gained much attention in engineering seismology. In the past few years, several ground motion prediction equations (GMPEs) and inter-frequency correlation structure of FAS have been established. Due to many preferable characteristics of FAS, probabilistic seismic hazard/risk analysis is rapidly changing from ergodic, spectrum acceleration Sa(T 0 )-based approach to non-ergodic, site-specific, FAS-based approach. This paper presents time domain intrusive framework for probabilistic seismic risk analysis using GMPE of FAS. Herein, methodology for time domain stochastic ground motion modeling based on GMPEs of FAS is presented in some detail. The simulated uncertain motions are modeled as a random process and represented by polynomial chaos Karhunen-Loève expansion. The random process excitations are further propagated into the uncertain structural system using Galerkin stochastic finite element method (SFEM). Probabilistic evolution of structural response is solved, and such solution is used to develop seismic risk for any damage state. The presented framework is illustrated through seismic risk analysis of a four-story building subjected to possible earthquakes from two strike slip faults. The influences of the epistemic uncertainties in source stress drop Δσ and site attenuation κ0 on seismic risk are investigated. The need for non-ergodic seismic risk analysis with source-specific and site specific characterizations is emphasized.

58 GEOSCIENCES↗
DOE OSTI2020

Cosmic ray transport in mixed magnetic fields and their role on the observed anisotropies

ABSTRACT There is a growing set of observational data demonstrating that cosmic rays exhibit small-scale anisotropies (5°–30°) with amplitude deviations lying between 0.01–0.1 per cent that of the average cosmic ray flux. A broad range of models have been proposed to explain these anisotropies ranging from finite-scale magnetic field structures to dark matter annihilation. The standard diffusion transport methods used in cosmic ray propagation do not capture the transport physics in a medium with finite-scale or coherent magnetic field structures. Here, we present a Monte Carlo transport method, applying it to a series of finite-scale magnetic field structures to determine the requirements of such fields in explaining the observed cosmic ray, small-scale anisotropies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Here we developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Bosonic field digitization for quantum computers

Quantum simulation of quantum field theory is a flagship application of quantum computers that promises to deliver capabilities beyond classical computing. The realization of quantum advantage will require methods that can accurately predict error scaling as a function of the resolution and parameters of the model and that can be implemented efficiently on quantum hardware. In this paper, we address the representation of lattice bosonic fields in a discretized field amplitude basis, develop methods to predict error scaling, and present efficient qubit implementation strategies. A low-energy subspace of the bosonic Hilbert space, defined by a boson occupation number cutoff, can be represented with exponentially good accuracy by a low-energy subspace of a finite-size Hilbert space. The finite representation construction and the associated errors are directly related to the accuracy of the Nyquist-Shannon sampling and the finite Fourier transforms of the boson number states in the field and the conjugate-field bases. We analyze the relation between the boson mass, the discretization parameters used for wave function sampling, and the finite representation size. Numerical simulations of small size Φ 4 problems demonstrate that the boson mass optimizing the sampling of the ground state wave function is a good approximation to the optimal boson mass yielding the minimum low-energy subspace size. However, we find that accurate sampling of general wave functions does not necessarily result in accurate representation. Finally, we develop methods for validating and adjusting the discretization parameters to achieve more accurate simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Large Area Detection of Microstructural Defects with Multi-Mode Ultrasonic Signals

Cyclic loading or other stresses can lead to development of cracks and crack growth in mechanical structures, leading to eventual failure. While ultrasound imaging can be used for non-destructive testing of such structures, conventional ultrasound techniques are often limited by crack size, density, and areal coverage. An effective characterization of real-world, large-area structures is required at an early damage stage to prevent catastrophic failure and predict remaining life. In this study, a new nonlinear ultrasonic testing (NUT) method is proposed for large-area monitoring of practical structures with arbitrary complexity by using multiple-mode guided-wave ultrasonic signals. The proposed guided-wave NUT technique requires single-element transducers, simple electronics, and a mixed time-frequency domain signal processing. As a proof-of-concept demonstration, numerical simulations and experiments are performed on an A36 carbon steel beam assembly with previously formed microstructural defects that cause nonlinearities in ultrasonic response. The quadratic dependence of the nonlinear wave excitation on the input ultrasonic signal amplitude is shown by numerical simulations, and such a nonlinear ultrasonic response is experimentally observed in the zone with a high density of microstructural defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Adhesion of Self-Complementary, Sinusoidal Surfaces Fabricated Using Two-Photon Polymerization

Microscale, pick-and-place assembly is a non-lithographic assembly method poised to impact diverse fields including flexible electronics, microfluidics and robotics. However, a major technological challenge is the need to deterministically control adhesion between parts. Here, switchable adhesion involving 3D-printed, self-complementary surfaces is demonstrated. Mechanical properties of metasurfaces pressed against flat, rigid substrates are modeled using finite element methods. A series of flat slabs and metastructured slabs with 2D sinusoidal surfaces are printed using two-photon polymerization (2PP) of a shape-memory resin. The surface frequency of featured slabs was varied between $3.\bar3$ mm −1 and $26.\bar6$ mm −1 with similar amplitudes. Adhesion between printed metasurfaces and glass and between printed, self-complementary metasurfaces is studied above and below the cured resin’s glass transition temperature (∼45 °C). Simple heating of adhering surfaces to above 60 °C lowers adhesion, and compression of surfaces while above the glass transition temperature followed by cooling to room temperature elevates adhesion. The nominal adhesive strength between printed, self-complementary surfaces, as determined by the maximum observable pull-off stress, exceeds 3 MPa. Further tailoring complementary surfaces for adhesion control may facilitate microscale disassembly for recovery of components or precious metals.

amorphous materials↗
DOE OSTI2024

Integrating the full four-loop negative geometries and all-loop ladder-type negative geometries in ABJM theory

The decomposition of the four-point ABJM amplituhedron into negative geometries produces compact integrands of logarithmic of amplitudes such that the infrared divergence only comes from the last loop integration, from which we can compute the cusp anomalous dimension of the ABJM theory. In this note, we integrate L – 1 loop momenta of the L-loop negative geometries for all four-loop negative geometries and a special class of all-loop ladder-type negative geometries by a method based on Mellin transformation, and from these finite quantities we extract the corresponding contribution to the cusp anomalous dimension. We find that the infrared divergence of a box-type negative geometry at L = 4 is weaker than other negative geometries, then only tree-type negative geometries contribute to the cusp anomalous dimension at L = 4. For the all-loop ladder-type negative geometries, we prove and conjecture some recursive structures as integral equations in Mellin space and find that they cannot contribute zeta values like ζ 3 , ζ 5 to the cusp anomalous dimension.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A Dynamic Amplitude-Correcting Gradient Estimation Technique to Align X-ray Focusing Optics

High-brightness X-rays, as produced at synchrotrons and X-ray free electron laser (XFEL) facilities, are used to characterize materials in a variety of scientific experiments. In most cases, effective use of the high-energy light requires precisely-aligned focusing optics; one example being a compound refractive lens (CRL). To align a CRL, the position and rotation must be optimized along four axes. In practice, this is a labor-intensive, time-consuming manual process that can monopolize scarce experimental time at the necessary X-ray facilities. Models of the expected Xray transmission function suggest that this task can be automated; however, the temporally-varying intensity at X-ray free electron laser facilities preclude the direct use of standard implementations of optimization solvers such as steepest descent algorithms. In this paper, we propose a novel technique to estimate the gradient of noisy functions with temporally-varying amplitudes. We construct this dynamicamplitude correction by systematically sampling a fixed central location within the standard finite difference stencil, accounting for the observed changes in time, and normalizing the difference quotients against those fluctuations. In addition to a rigorous error analysis of the sampling technique, we demonstrate its efficacy in stochastic descent optimization methods. Further, we demonstrate how this approach may be implemented to optimize X-ray focusing optics at synchrotrons or XFEL facilities

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗
DOE OSTI2025

First Look at Quartic-in-Spin Binary Dynamics at Third Post-Minkowskian Order

We compute the conservative and radiation-reaction contributions to classical observables in the gravitational scattering between a spinning and a spinless black hole to the fourth order in spin and third order in the gravitational constant. The conservative results are obtained from two-loop amplitudes for the scattering process of a massive scalar with a massive spin-𝑠 field (𝑠=0, 1, 2) minimally coupled to gravity, employing the recently introduced spin interpolation method to resolve all spin-Casimir terms. The two-loop amplitude exhibits a spin-shift symmetry in both probe limits, which we conjecture to be a sign of yet unknown integrability of Kerr orbits through the quartic order in spin and to all orders in the gravitational constant. We obtain the radial action from the finite part of the amplitude and use it to compute classical observables, including the impulse and spin kick. This is done using the recently introduced covariant Dirac brackets, which allow for the computation of classical scattering observables for general (nonaligned) spin configurations. Finally, employing the radiation-reaction amplitude proposed by Alessio and Di Vecchia, together with the Dirac brackets, we obtain radiation-reaction contributions to observables at all orders in spin and beyond the aligned-spin limit. We find agreement with known results up to the quadratic order in spin for both conservative and radiation-reaction contributions. Our results advance the state of the art in the understanding of spinning binary dynamics in general relativity and demonstrate the power and simplicity of the Dirac bracket formalism for relating scattering amplitudes to classical observables.

general relativity↗
DOE OSTI2025

Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description

Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.

ab initio calculations↗
DOE OSTI2024

Theoretical Description of Pump-Probe Experiments in Charge-Density-Wave Materials out to Long Times

We describe coupled nonequilibrium electron-phonon systems semiclassically—Ehrenfest dynamics for the phonons and quantum mechanics for the electrons—using a classical Monte Carlo approach that determines the nonequilibrium response to a large pump field. The semiclassical approach is expected to be accurate, because the phonons are excited to average energies much higher than the phonon frequency, eliminating the need for a quantum description. The numerical efficiency of this method allows us to perform a self-consistent time evolution out to very long times (tens of picoseconds), enabling us to model pump-probe experiments of a charge-density-wave (CDW) material. Our system is a half-filled, one-dimensional (1D) Holstein chain that exhibits CDW ordering due to a Peierls transition. The chain is subjected to a time-dependent electromagnetic pump field that excites it out of equilibrium, and then a second probe pulse is applied after a time delay. By evolving the system to long times, we capture the complete process of lattice excitation and subsequent relaxation to a new equilibrium, due to an exchange of energy between the electrons and the lattice, leading to lattice relaxation at finite temperatures. We employ an indirect (impulsive) driving mechanism of the lattice by the pump pulse due to the direct driving of the electrons. We identify two driving regimes, where the pump can either cause small perturbations or completely invert the initial CDW order. Our work successfully describes the ringing of the amplitude mode in CDW systems that has long been seen in experiment but never successfully explained by microscopic theory. We also describe the fluence-dependent crossover that inverts the CDW order parameter and changes the phonon dynamics. Finally, we illustrate how this method can examine a number of different types of experiments including photoemission, x-ray diffraction, and two-dimensional (2D) spectroscopy. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

Stochastic Inversion of Gaussian Random Media Using Transverse Coherence Functions for Reflected Waves: Theory and Method

The transverse coherence functions (TCFs) of phase and amplitude fluctuations of a seismic wave are powerful to estimate the spatial distribution, length scales, and strength of random heterogeneities. However, TCFs have been formulated for transmitted waves only, not for reflected waves. In this paper, we derive reflection TCFs for Gaussian random media. Furthermore, we propose to invert for Gaussian random media using the reflection TCFs based on the grid search. We validate the new reflection TCF formulas using 2D finite-difference numerical experiments. The numerical example also illustrates the feasibility and efficiency of the inversion. The stochastic inversion using reflected waves can be used in both exploration and global seismology.

58 GEOSCIENCES↗
DOE OSTI2023

Primary observables for indirect searches at colliders

We consider the complete set of observables for collider searches for indirect effects of new heavy physics. They consist of SU(3) C × U(1) EM invariant interaction terms/operators that parameterize deviations from the Standard Model. We show that, un der very general assumptions, the leading deviations from the Standard Model are given by a finite number of ‘primary’ operators, with the remaining operators given by ‘Mandelstam descendants’ whose effects are suppressed by powers of Mandelstam variables divided by the mass scale M of the heavy physics. We explicitly determine all 3 and 4-point primary operators relevant for Higgs signals at colliders by using the correspondence between on-shell amplitudes and independent operators. We give a detailed discussion of the methods used to obtain this result, including a new analytical method for determining the independent operators. The results are checked using the Hilbert series that counts independent operators. We also give a rough sketch of the phenomenology, including unitarity bounds on the interaction strengths and rough estimates of their importance for Higgs decays at the HL-LHC. These results motivate further exploration of Higgs decays to $Z\overline{f}f$, $W\overline{f}f'$, $γ\overline{f}f$, and $Zγγ$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

AdS/CFT unitarity at higher loops: high-energy string scattering

What is the space of weakly-coupled, gravitational theories which contain massive, higher-spin particles? This class of theories is highly constrained and it is conjectured their ultraviolet completion must be string theory. We provide more evidence for this conjecture by studying the Regge limit in large N, 4d CFTs with single-trace operators of unbounded spin. We show that in the Regge limit, these theories have bulk scattering amplitudes which are consistent with the string theory prediction to all orders in 1/N for large, but finite, coupling. In the language of Regge theory, we show Pomeron exchange naturally exponentiates in the 1/N expansion. To do this, we solve the bootstrap equations at tree-level and then use the Lorentzian inversion formula to find the one-loop correlator in the Regge limit. This is a unitarity method for AdS/CFT which can be repeated iteratively to make all orders statements. We also explain under what conditions the tree-level result exponentiates in the 1/N expansion at arbitrary coupling. Finally, we comment on further inelastic effects and show they give subleading contributions at large coupling. As a consistency check, we recover results from bulk Einstein gravity in the limit where all higher-spin particles decouple.

1/N Expansion↗