Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave equation”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Baryon spectrum from superconformal quantum mechanics and its light-front holographic embedding

We describe the observed light-baryon spectrum by extending superconformal quantum mechanics to the light front and its embedding in AdS space. This procedure uniquely determines the confinement potential for arbitrary half-integer spin. To this end, we show that fermionic wave equations in AdS space are dual to light-front supersymmetric quantum-mechanical bound-state equations in physical space-time. The specific breaking of conformal invariance explains hadronic properties common to light mesons and baryons, such as the observed mass pattern in the radial and orbital excitations, from the spectrum generating algebra. Lastly, the holographic embedding in AdS also explains distinctive and systematic features, such as the spin-J degeneracy for states with the same orbital angular momentum, observed in the light-baryon spectrum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Decay of the de Sitter vacuum

The decay rate of the Bunch-Davies state of a massive scalar field in the expanding flat spatial sections of de Sitter space is determined by an analysis of the particle pair creation process in real time. The Feynman definition of particle and antiparticle Fourier mode solutions of the scalar wave equation and their adiabatic phase analytically continued to the complexified time domain show conclusively that the Bunch-Davies state is not the vacuum state at late times. The closely analogous creation of charged particle pairs in a uniform electric field is reviewed and Schwinger’s result for the vacuum decay rate is recovered by this same real time analysis. The vacuum decay rate in each case is also calculated by switching the background field on adiabatically, allowing it to act for a very long time, and then adiabatically switching it off again. In both the uniform electric field and de Sitter cases, the particles created while the field is switched on are verified to be real, in the sense that they persist in the final asymptotic flat zero-field region. In the de Sitter case, there is an interesting residual dependence of the rate on how the de Sitter phase is ended, indicating a greater sensitivity to spatial boundary conditions. The electric current of the created particles in the E-field case and their energy density and pressure in the de Sitter case are also computed, and the magnitude of their backreaction effects on the background field estimated. Finally, possible consequences of the Hubble scale instability of the de Sitter vacuum for cosmology, vacuum dark energy, and the cosmological “constant” problem are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Positive Energy Functional for Massless Scalars in Rotating Black Hole Backgrounds of Maximal Ungauged Supergravity

Here, we outline a proof of the stability of a massless neutral scalar field ψ in the background of a wide class of four dimensional asymptotically flat rotating and “electrically charged” solutions of supergravity, and the low energy limit of string theory, known as STU metrics. Despite their complexity, we find it possible to circumvent the difficulties presented by the existence of ergo regions and the related phenomenon of superradiance in the original metrics by following a strategy due to Whiting, and passing to an auxiliary metric admitting an everywhere lightlike Killing field and constructing a scalar field ψ (related to a possible unstable mode ψ by a nonlocal transformation) which satisfies the massless wave equation with respect to the auxiliary metric. By contrast with the case for ψ, the associated energy density of ψ is not only conserved but is also non-negative.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A New Model of Frequency Delay in Power Systems

Observations of slow frequency signals in power systems have been explained by continuum wave equations that approximate the system's electromechanical dynamics. These continuum models predict that the signal speed is a function of the rotational inertia of the generators and impedances of the transmission lines. Here, we show that frequency signals with speeds modulated by these factors are predicted by the swing and Kirchhoff equations, without resorting to a continuum model. Our analysis also shows that the precision of the sensors used to measure frequency has a significant impact on the apparent speed of the frequency signal.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

Feasibility of RACT for 3D dose measurement and range verification in a water phantom

Purpose: The objective of this study is to establish the feasibility of using radiation-induced acoustics to measure the range and Bragg peak dose from a pulsed proton beam. Simulation studies implementing a prototype scanner design based on computed tomographic methods were performed to investigate the sensitivity to proton range and integral dose. Methods: Derived from thermodynamic wave equation, the pressure signals generated from the dose deposited from a pulsed proton beam with a 1 cm lateral beam width and a range of 16, 20, and 27 cm in water using Monte Carlo methods were simulated. The resulting dosimetric images were reconstructed implementing a 3D filtered backprojection algorithm and the pressure signals acquired from a 71-transducer array with a cylindrical geometry (30 × 40 cm) rotated over 2π about its central axis. Dependencies on the detector bandwidth and proton beam pulse width were performed, after which, different noise levels were added to the detector signals (using 1 μs pulse width and a 0.5 MHz cutoff frequency/hydrophone) to investigate the statistical and systematic errors in the proton range (at 20 cm) and Bragg peak dose (of 1 cGy). Results: The reconstructed radioacoustic computed tomographic image intensity was shown to be linearly correlated to the dose within the Bragg peak. And, based on noise dependent studies, a detector sensitivity of 38 mPa was necessary to determine the proton range to within 1.0 mm (full-width at half-maximum) (systematic error < 150 μm) for a 1 cGy Bragg peak dose, where the integral dose within the Bragg peak was measured to within 2%. For existing hydrophone detector sensitivities, a Bragg peak dose of 1.6 cGy is possible. Conclusions: This study demonstrates that computed tomographic scanner based on ionizing radiation-induced acoustics can be used to verify dose distribution and proton range with centi-Gray sensitivity. Realizing this technology into the clinic has the potential to significantly impact beam commissioning, treatment verification during particle beam therapy and image guided techniques.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

SU-D-210-03: Limited-View Multi-Source Quantitative Photoacoustic Tomography

Purpose: This work is to investigate a novel limited-view multi-source acquisition scheme for the direct and simultaneous reconstruction of optical coefficients in quantitative photoacoustic tomography (QPAT), which has potentially improved signal-to-noise ratio and reduced data acquisition time. Methods: Conventional QPAT is often considered in two steps: first to reconstruct the initial acoustic pressure from the full-view ultrasonic data after each optical illumination, and then to quantitatively reconstruct optical coefficients (e.g., absorption and scattering coefficients) from the initial acoustic pressure, using multi-source or multi-wavelength scheme.Based on a novel limited-view multi-source scheme here, We have to consider the direct reconstruction of optical coefficients from the ultrasonic data, since the initial acoustic pressure can no longer be reconstructed as an intermediate variable due to the incomplete acoustic data in the proposed limited-view scheme. In this work, based on a coupled photo-acoustic forward model combining diffusion approximation and wave equation, we develop a limited-memory Quasi-Newton method (LBFGS) for image reconstruction that utilizes the adjoint forward problem for fast computation of gradients. Furthermore, the tensor framelet sparsity is utilized to improve the image reconstruction which is solved by Alternative Direction Method of Multipliers (ADMM). Results: The simulation was performed on a modified Shepp-Logan phantom to validate the feasibility of the proposed limited-view scheme and its corresponding image reconstruction algorithms. Conclusion: A limited-view multi-source QPAT scheme is proposed, i.e., the partial-view acoustic data acquisition accompanying each optical illumination, and then the simultaneous rotations of both optical sources and ultrasonic detectors for next optical illumination. Moreover, LBFGS and ADMM algorithms are developed for the direct reconstruction of optical coefficients from the acoustic data. Jing Feng and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000) and the Shanghai Pujiang Talent Program (#14PJ1404500)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Multilevel Convergence Analysis of Multigrid-Reduction-in-Time

This study presents a multilevel convergence framework for multigrid-reduction-in-time (MGRIT) as a generalization of previous two-grid estimates. The framework provides a priori upper bounds on the convergence of MGRIT V- and F-cycles, with different relaxation schemes, by deriving the respective residual and error propagation operators. The residual and error operators are functions of the time-stepping operator, analyzed directly and bounded in the norm, both numerically and analytically. We present various upper bounds of different computational cost and varying sharpness. These upper bounds are complemented by proposing analytic formulae for the approximate convergence factor of V-cycle algorithms that take the number of fine grid time points, the temporal coarsening factors, and the eigenvalues of the time-stepping operator as parameters. The paper concludes with supporting numerical investigations of parabolic (anisotropic diffusion) and hyperbolic (wave equation) model problems. We assess the sharpness of the bounds and the quality of the approximate convergence factors. Observations from these numerical investigations demonstrate the value of the proposed multilevel convergence framework for estimating MGRIT convergence a priori and for the design of a convergent algorithm. We further highlight that observations in the literature are captured by the theory, including that two-level Parareal and multilevel MGRIT with F-relaxation do not yield scalable algorithms and the benefit of a stronger relaxation scheme. An important observation is that with increasing numbers of levels MGRIT convergence deteriorates for the hyperbolic model problem, while constant convergence factors can be achieved for the diffusion equation. The theory also indicates that L-stable Runge--Kutta schemes are more amendable to multilevel parallel-in-time integration with MGRIT than A-stable Runge--Kutta schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗
DOE OSTI2026

From the Great Wave of Translation to the Force between Quarks

Here, the chance observation of a novel traveling wave in a canal led over time to the formulation of a nonlinear wave equation—the Korteweg–de Vries equation—that describes strikingly robust disturbances now called solitons. The figure of an isolated soliton corresponds to a reflectionless potential that supports a single bound state in the one-dimensional Schrödinger equation. An appropriate combination of individual solitons yields a symmetric reflectionless potential that supports multiple bound states. Thus, the KdV equation opens the path to solving the inverse scattering problem for a collection of bound states. Applied to the quarkonium spectra, this formalism allows the construction of reflectionless approximations to the confining potentials that account for the force between quarks, and to tests of the flavor-independence of the interquark interaction.

Inverse Scattering↗
DOE OSTI2023

Sparse Approximate Multifrontal Factorization with Composite Compression Methods

This article presents a fast and approximate multifrontal solver for large sparse linear systems. In a recent work by Liu et al., we showed the efficiency of a multifrontal solver leveraging the butterfly algorithm and its hierarchical matrix extension, HODBF (hierarchical off-diagonal butterfly) compression to compress large frontal matrices. The resulting multifrontal solver can attain quasi-linear computation and memory complexity when applied to sparse linear systems arising from spatial discretization of high-frequency wave equations. To further reduce the overall number of operations and especially the factorization memory usage to scale to larger problem sizes, in this article we develop a composite multifrontal solver that employs the HODBF format for large-sized fronts, a reduced-memory version of the nonhierarchical block low-rank format for medium-sized fronts, and a lossy compression format for small-sized fronts. This allows us to solve sparse linear systems of dimension up to 2.7 × larger than before and leads to a memory consumption that is reduced by 70% while ensuring the same execution time. The code is made publicly available in GitHub.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Seismic Waves, 4th order accurate

SW4 is a program for simulating seismic wave propagation on parallel computers. SW4 colves the seismic wave equations in Cartesian corrdinates. It is therefore appropriate for regional simulations, where the curvature of the earth can be neglected. SW4 implements a free surface boundary condition on a realistic topography, absorbing super-grid conditions on the far-field boundaries, and a kinematic source model consisting of point force and/or point moment tensor source terms. SW4 supports a fully 3-D heterogeneous material model that can be specified in several formats. SW4 can output synthetic seismograms in an ASCII test format, or in the SAC finary format. It can also present simulation information as GMT scripts, whixh can be used to create annotated maps. Furthermore, SW4 can output the solution as well as the material model along 2-D grid planes.

Source record↗
DOE OSTI2025

Atmospheric Profile Builder

APBuilder is a software tool designed to generate atmospheric profiles for use with AC2Dr by downloading weather model data and transforming it into 1D or 2D binary profiles. These binary profiles serve as required inputs for AC2Dr, an LLNL-developed, open source, two-dimensional numerical solver for the acoustic wave equation based on the finite difference method.

Kim, Keehoon [Lawrence Livermore National Laborato↗
DOE OSTI2026

NLIN-BST

Nonlinear scalar wave equation and Einstein equations solver with higher dimensional black hole and black string spacetime

Lim, Hyun [Los Alamos National Laboratory]↗
DOE OSTI2024

A graphics processing unit accelerated sparse direct solver and preconditioner with block low rank compression

We present the GPU implementation efforts and challenges of the sparse solver package STRUMPACK. The code is made publicly available on github with a permissive BSD license. STRUMPACK implements an approximate multifrontal solver, a sparse LU factorization which makes use of compression methods to accelerate time to solution and reduce memory usage. Multiple compression schemes based on rank-structured and hierarchical matrix approximations are supported, including hierarchically semi-separable, hierarchically off-diagonal butterfly, and block low rank. Here, in this paper, we present the GPU implementation of the block low rank (BLR) compression method within a multifrontal solver. Our GPU implementation relies on highly optimized vendor libraries such as cuBLAS and cuSOLVER for NVIDIA GPUs, rocBLAS and rocSOLVER for AMD GPUs and the Intel oneAPI Math Kernel Library (oneMKL) for Intel GPUs. Additionally, we rely on external open source libraries such as SLATE (Software for Linear Algebra Targeting Exascale), MAGMA (Matrix Algebra on GPU and Multi-core Architectures), and KBLAS (KAUST BLAS). SLATE is used as a GPU-capable ScaLAPACK replacement. From MAGMA we use variable sized batched dense linear algebra operations such as GEMM, TRSM and LU with partial pivoting. KBLAS provides efficient (batched) low rank matrix compression for NVIDIA GPUs using an adaptive randomized sampling scheme. The resulting sparse solver and preconditioner runs on NVIDIA, AMD and Intel GPUs. Interfaces are available from PETSc, Trilinos and MFEM, or the solver can be used directly in user code. We report results for a range of benchmark applications, using the Perlmutter system from NERSC, Frontier from ORNL, and Aurora from ALCF. For a high frequency wave equation on a regular mesh, using 32 Perlmutter compute nodes, the factorization phase of the exact GPU solver is about 6.5× faster compared to the CPU-only solver. The BLR-enabled GPU solver is about 13.8× faster than the CPU exact solver. For a collection of SuiteSparse matrices, the STRUMPACK exact factorization on a single GPU is on average 1.9× faster than NVIDIA’s cuDSS solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Compensating attenuation effects in full-waveform inversion with dissipation-dispersion decoupling

Seismic attenuation poses challenges to velocity model building from acoustic/elastic full-waveform inversion (FWI). In particular, when constructing the FWI gradient for velocity inversion, recent studies have indicated that the high-attenuation structure could distort the gradient by damping the amplitudes and shifting the kinematic phases, resulting in an imbalanced update and thus unreliable velocity model. These Q effects are particularly significant in the reflection acquisition geometry due to the “double-damping” issue. Here, we develop a Q-compensated FWI algorithm for constructing a Q-free FWI gradient. By using a recently developed viscoacoustic wave propagator, this compensation can be done conveniently by keeping the dispersion term and flipping the dissipation term in the wave equation when we simulate the forward and time-reversed adjoint wavefields. The resultant gradient obtained by interacting these two wavefields has correct kinematics and Q-free amplitudes. This Q-compensated FWI can balance the update between Q- and no- Q-areas, which we determine using synthetic examples. In addition, we illustrate how to take advantage of the dissipation-dispersion decoupling to determine the anomalous Q value in the Q model building workflow via Q-compensated reverse time migration.

Geochemistry & Geophysics↗
DOE OSTI2017

Polarized optical scattering by inhomogeneities and surface roughness in an anisotropic thin film

We extend the theory for scattering by oblique columnar structure thin films to include the induced form birefringence and the propagation of radiation in those films. We generalize the 4 × 4 matrix theory to include arbitrary sources in the layer, which are necessary to determine the Green function for the inhomogeneous wave equation. We further extend first-order vector perturbation theory for scattering by roughness in the smooth surface limit, when the layer is anisotropic. Scattering by an inhomogeneous medium is approximated by a distorted Born approximation, where effective medium theory is used to determine the effective properties of the medium and strong fluctuation theory is used to determine the inhomogeneous sources. In this manner, we develop a model for scattering by inhomogeneous films, with anisotropic correlation functions. Here, the results are compared to Mueller matrix bidirectional scattering distribution function measurements for a glancing-angle deposition (GLAD) film. While the results are applied to the GLAD film example, the development of the theory is general enough that it can guide simulations for scattering in other anisotropic thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗
DOE OSTI2019

Thermal decay without information loss in horizonless microstate geometries

We develop a new hybrid WKB technique to compute boundary-to-boundary scalar Green functions in asymptotically-AdS backgrounds in which the scalar wave equation is separable and is explicitly solvable in the asymptotic region. We apply this technique to a family of six-dimensional \frac{1}{8} 1 8 -BPS asymptotically AdS _3\,\times\, 3 × S ^3 3 horizonless geometries that have the same charges and angular momenta as a D1-D5-P black hole with a large horizon area. At large and intermediate distances, these geometries very closely approximate the extremal-BTZ \,\times\, × S ^3 3 geometry of the black hole, but instead of having an event horizon, these geometries have a smooth highly-redshifted global-AdS _3\,\times\, 3 × S ^3 3 cap in the IR. We show that the response function of a scalar probe, in momentum space, is essentially given by the pole structure of the highly-redshifted global-AdS _3 3 modulated by the BTZ response function. In position space, this translates into a sharp exponential black-hole-like decay for times shorter than N_1 N_5 N 1 N 5 , followed by the emergence of evenly spaced “echoes from the cap,” with period \sim N_1 N_5 ∼ N 1 N 5 . Our result shows that horizonless microstate geometries can have the same thermal decay as black holes without the associated information loss.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗