Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave propagation”

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 OSTI2023

Towards the usage of quarry blasts as an active seismic source for subsurface characterization and imaging

SUMMARY Active or man-made seismic sources are often used to delineate subsurface geological structures via seismic imaging techniques. Although conventional controlled seismic sources have provided high-resolution subsurface images, the high cost of data acquisition necessitates improved use of cheaper alternative seismic sources for subsurface imaging, such as traffic noise, quarry blasts and earthquakes. Here, we explore the potential use of quarry blasts as an active seismic source for subsurface imaging. The goal of our study is to demonstrate the utility of quarry blasts as seismic sources for investigating subsurface structure with a field experiment in Mill Creek, southern Oklahoma. We deployed a 7455-m-long receiver line of 72 receivers (105 m spacing) extending away from a quarry, with a recording time of 28 d. We analyse similarities between blasts, wave types, origin time and excitation location of quarry blasts, which is important information for using such signal as an active source. Given a relative homogeneous geological condition in the study area, we could not observe clear reflected waves. However, the use of P- and S-wave velocities led to the delineation of a major basement damage zone in the study area, which further demonstrates the viability of quarry blasts as a useful seismic source for subsurface geological imaging.

58 GEOSCIENCES↗
DOE OSTI2024

GLAD-M35: a joint P and S global tomographic model with uncertainty quantification

We present our third and final generation joint P and S global adjoint tomography (GLAD) model, GLAD-M35, and quantify its uncertainty based on a low-rank approximation of the inverse Hessian. Starting from our second-generation model, GLAD-M25, we added 680 new earthquakes to the database for a total of 2160 events. New P-wave categories are included to compensate for the imbalance between P- and S-wave measurements, and we enhanced the window selection algorithm to include more major-arc phases, providing better constraints on the structure of the deep mantle and more than doubling the number of measurement windows to 40 million. Two stages of a Broyden–Fletcher–Goldfarb–Shanno (BFGS) quasi-Newton inversion were performed, each comprising five iterations. With this BFGS update history, we determine the model’s standard deviation and resolution length through randomized singular value decomposition.

58 GEOSCIENCES↗
DOE OSTI2025

Estimating an airborne dipole source using 3-D wavefield simulations and seismic receivers on the ground

Acoustic signals can couple to the ground, providing an opportunity to use seismic stations to investigate airborne sources. The study of Bishop et al. used wavefield simulations in a fluid-solid medium to quantify the role of topography on the seismic (ground) recordings of a monopole source in the air. We build upon this study by linking wavefield forward modelling with the source estimation code MTUQ, which can accommodate point forces or moment tensors in a solid medium, as well as sources in the air (or water) if they are enabled by the forward-modelling solver. We perform a series of synthetic numerical experiments to demonstrate that a dipole airborne source can be estimated using ground-based receivers, even within the presence of realistic topography. We investigate the influence of receiver coverage, topography, and assumed source location on the estimated results. The established capabilities raise the prospects for future efforts to estimate dipole sources in 3-D models that include heterogeneity in the air and the earth in addition to topography.

58 GEOSCIENCES↗
DOE OSTI2025

Data-informed grid refinement to improve traveltime accuracy in the regional seismic traveltime (RSTT) model

The regional seismic traveltime (RSTT) model predicts traveltimes of regional seismic phases accounting for 3-D structure of the crust and the upper mantle on a global scale. Previous versions of the RSTT model have been implemented using nodes separated by ∼1° spacing across the globe. A regional-scale study using regional Pn and Pg traveltimes across Israel and the Middle East demonstrated that data driven, systematic grid refinement reduces traveltime residuals and enhances resolution of smaller tectonic features in regions having dense ray coverage. High density Pn ray coverage in the western US, Europe, Middle East and East Asia can likewise provide the resolution that allows systematic global grid refinement of the RSTT model. In this study, we use a large number of Pn ray paths originating from events located with an epicentral location uncertainty of 25 km (GT25) or better. We conduct targeted grid refinements at 1.0°, 0.5°, 0.25° and 0.125° on a global scale, producing a refined RSTT model that yields a 21.6 per cent reduction in median event location error in Europe and the Middle East, when compared with the original global RSTT model presented in Begnaud et al. The new model also resolves finer tectonic structures in regions with high Pn ray density.

58 GEOSCIENCES↗
DOE OSTI2024

Unified wavefront singularity characterization of three-dimensional elastodynamic time-domain half-space Green's function under impulsive boundary and internal loads

Founded on a novel analytical formulation that led to a rigorous yet compact path-integral representation of the time-domain elastodynamic half-space Green's function, a unified analysis of the possible occurrence of different singular wavefront behaviour in the response under arbitrary impulsive internal or surface point loads at arbitrary source-receiver locations is presented. With the decomposition of the general solution into distinct initiating and reflected wave group integrals that share a common factored format and simple contour definitions, the mathematical framework is shown to allow a straightforward identification of the specific conditions and the particular wave groups that are responsible for the singular wavefront phenomena without resorting to advanced analytic function theories or asymptotic methods. Analytic characterizations of the nature, strength and direction of all intrinsic singular wavefront behaviours of the three-dimensional Green's function in three canonical cases of source-receiver configurations are given in a dual integral-closed form format to facilitate their theoretical understanding as well as computational applications. Graphical illustrations of their variation with the source-receiver configuration and the medium's Poisson's ratio together with relevant comparison and clarifications of some classical treatments are included.

Green’s functions↗
DOE OSTI2020

Compression-induced solidification of shock-melted cerium

Compression induced solidification has been observed in cerium on nanosecond timescales. For this study, a series of experiments were conducted on Sandia National Laboratories' Z Facility in which cerium was shock melted and subsequently shocklessly, or ramp, loaded across the melt line inducing solidification. The signature of solidification manifested in the recovery of material strength and the propagation of waves at the local elastic sound velocity. Density functional theory simulations of cerium along the experimental phase-space path exhibit spontaneous freezing to a tetragonal phase at the same pressure and closely predict the observed physical properties of solid and liquid cerium near melt.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Perturbative quantum evolution of the gravitational state and dressing in general backgrounds

This paper sets up a perturbative treatment of the evolving quantum state of a gravitational system, in a Schrödinger-like picture, working about a general background. This connects gauge symmetry, the constraints, gravitational dressing, and evolution. Starting with a general time slicing, we give a simple derivation of the relation between the constraints, the Hamiltonian, and its well-known boundary term. Among different approaches to quantization with constraints, we focus on a “gauge-invariant canonical quantization,” which is developed perturbatively in the gravitational coupling. The leading-order solution of the constraints (including the Wheeler-DeWitt equation) for perturbations about the background is given in terms of an explicit construction of gravitational dressings built using certain generalized Green’s functions; different such dressings corresponding to adding propagating gravitational waves to a particular solution of the constraints. Dressed operators commute with the constraints, expressing their gauge invariance, and have an algebraic structure differing significantly from the undressed operators of the underlying field theory. These operators can act on the vacuum to create dressed states, and evolution of general such states is then generated by the boundary Hamiltonian, and alternately may be characterized using other relational observables. This provides a concrete approach to studying perturbative time evolution, including the leading gravitational backreaction, of quantum states of black holes with flat or anti–de-Sitter asymptotics, for example on horizon-crossing slices. This description of evolution in turn provides a starting point for investigating possibly important corrections to quantum evolution, that go beyond quantized general relativity. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Phonon spectra of a two-dimensional solid dusty plasma modified by two-dimensional periodic substrates

Here, phonon spectra of a two-dimensional (2D) solid dusty plasma modulated by 2D square and triangular periodic substrates are investigated using Langevin dynamical simulations. The commensurability ratio, i.e., the ratio of the number of particles to the number of potential well minima, is set to 1 or 2. The resulting phonon spectra show that propagation of waves is always suppressed due to the confinement of particles by the applied 2D periodic substrates. For a commensurability ratio of 1, the spectra indicate that all particles mainly oscillate at one specific frequency, corresponding to the harmonic oscillation frequency of one single particle inside one potential well. At a commensurability ratio of 2, the substrate allows two particles to sit inside the bottom of each potential well, and the resulting longitudinal and transverse spectra exhibit four branches in total. We find that the two moderate branches come from the harmonic oscillations of one single particle and two combined particles in the potential well. The other two branches correspond to the relative motion of the two-body structure in each potential well in the radial and azimuthal directions. The difference in the spectra between the square and triangular substrates is attributed to the anisotropy of the substrates and the resulting alignment directions of the two-body structure in each potential well.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Absorption of charged particles in perfectly matched layers by optimal damping of the deposited current

Perfectly matched layers (PMLs) are widely used in particle-in-cell simulations, in order to absorb electromagnetic waves that propagate out of the simulation domain. However, when charged particles cross the interface between the simulation domain and the PMLs, a number of numerical artifacts can arise. Here, in order to mitigate these artifacts, we introduce a PML algorithm whereby the current deposited by the macroparticles in the PML is damped by an analytically derived optimal coefficient. The benefits of this algorithm are illustrated in practical simulations. In particular, it is shown that this algorithm is well suited for particles exiting the box in near-normal incidence, in the sense that the fields behave as if the exiting particle is propagating in an infinite vacuum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Microscopic Mechanisms of Glasslike Lattice Thermal Transport in Cubic Cu 12 Sb 4 S 13 Tetrahedrites

Materials based on cubic tetrahedrites (Cu 12 Sb 4 S 13 ) are useful thermoelectrics with unusual thermal and electrical transport properties, such as very low and nearly temperature-independent lattice thermal conductivity (κ L ). We explain the microscopic origin of the glasslike κ L in Cu 12 Sb 4 S 13 by explicitly treating anharmonicity up to quartic terms for both phonon energies and phonon scattering rates. We show that the strongly unstable phonon modes associated with trigonally coordinated Cu atoms are anharmonically stabilized above approximately 100 K and continue hardening with increasing temperature in accord with experimental data. This temperature-induced hardening effect reduces scattering of heat carrying acoustic modes by reducing the available phase space for three-phonon processes, thereby balancing the conventional ∝ T increase in scattering due to phonon population and yielding nearly temperature independent κ L . Furthermore, we find that very strong phonon broadening leads to a qualitative breakdown of the conventional phonon-gas model and modify the dominant heat transport mechanism from the particlelike phonon wave packet propagation to incoherent contributions described by the off-diagonal terms in the heat-flux operator, which are typically prevailing in glasses and disordered crystals. Furthermore, our work paves the way to a deeper understanding of glasslike thermal conductivity in complex crystals with strong anharmonicity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Time Refraction of Spin Waves

Here, we present an experimental study of time refraction of spin waves (SWs) propagating in microscopic waveguides under the influence of time-varying magnetic fields. Using space- and time-resolved Brillouin light scattering microscopy, we demonstrate that the broken translational symmetry along the time coordinate results in a loss of energy conservation for SWs and thus allows for a broadband and controllable shift of the SW frequency. With an integrated design of SW waveguide and microscopic current line for the generation of strong, nanosecond-long, magnetic field pulses, a conversion efficiency up to 39% of the carrier SW frequency is achieved, significantly larger compared to photonic systems. Given the strength of the magnetic field pulses and its strong impact on the SW dispersion relation, the effect of time refraction can be quantified on a length scale comparable to the SW wavelength. Furthermore, we utilize time refraction to excite SW bursts with pulse durations in the nanosecond range and a frequency shift depending on the pulse polarity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Ytterbium ion trap quantum computing: The current state-of-the-art

In this work, we present an overview of contemporary quantum computing with ytterbium ion traps, placing the emphasis on industry implementations. We provide brief, concrete descriptions of various key features, such as trap loading, electronic structure, qubit function, gates, error analysis, and benchmarking. We focus on the underlying science and current technologies to provide readers with a holistic picture of available techniques for using ytterbium in contemporary ion trap designs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Generalized analytical and numerical modeling of optical second harmonic generation in anisotropic crystals and complex heterostructures using #SHAARP package

Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. The closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, here we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP) for single interface (si) and in multilayers (ml) for homogeneous crystals. The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials. SHAARP.si and SHAARP.ml are available through GitHub https://github.com/Rui-Zu/SHAARP and https://github.com/bzw133/SHAARP.ml, respectively.

complex systems↗
DOE OSTI2022

Dynamics of active liquid interfaces

Controlling interfaces of phase-separating fluid mixtures is key to the creation of diverse functional soft materials. Traditionally, this is accomplished with surface-modifying chemical agents. Using experiment and theory, we studied how mechanical activity shapes soft interfaces that separate an active and a passive fluid. Chaotic flows in the active fluid give rise to giant interfacial fluctuations and noninertial propagating active waves. At high activities, stresses disrupt interface continuity and drive droplet generation, producing an emulsion-like active state composed of finite-sized droplets. When in contact with a solid boundary, active interfaces exhibit nonequilibrium wetting transitions, in which the fluid climbs the wall against gravity. These results demonstrate the promise of mechanically driven interfaces for creating a new class of soft active matter.

Science & Technology - Other Topics↗
DOE OSTI2021

Rock Fracturing Using High-Pressure Ethylene/Nitrous Oxide Detonations

The present work investigates high initial pressure detonations of a stoichiometric mixture of ethylene and nitrous oxide (C2H4 + 6N2O) as a method of fracturing rock beneath the ground surface. These tests were conducted at a test site operated by the Energetic Materials Research and Testing Center (EMRTC), Socorro, New Mexico. The volume under the surface used for testing (called the Down Hole Assembly) consists of a 0.438 in. ID x 50 ft. long stainless-steel tube running down from the test site to a well bore which is 3 in. ID x 10 ft. long and the rock in the well bore is exposed to the propagating combustion wave. The testing carried out at Zucrow Laboratories in the smaller, alloy steel combustion vessel provided a scaling of pressures expected in the well bore. The combustion is initiated by energizing an EBW (Exploding Bridge Wire) above the ground surface. The experimental setup accommodates one high pressure (100,000 psia) transducer to measure the pressure peak and is placed approximately 5 ft. above the ground surface and 5 ft. downstream of the EBW. The focus of this series of experiments is to investigate the dependence of fracture to the rock beneath the surface on initial pressures of the mixture of ethylene and nitrous oxide. Experiments were carried out at initial pressures varying between 125 psia and 300 psi. The transducer recorded elevated pressures, which were 2.3 to 2.6 times in excess of the CJ values. The experimental results are discussed and explained in this report.

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

SURFplus model calibration for PBX 9012

PBX 9012 is a plastic bonded explosive composed of 90.2 wt % HMX (cyclo-tetramethylene-tetranitramine) and 9.8 % Viton A (vinylidene-hexafluoropropylene copolymer). It has nearly the same composition as LX-07 (90% HMX and 10% Viton A) but with a slightly lower density. Hence it has a higher porosity and is slightly more shock sensitive. Here we describe the SURFplus model calibration for PBX 9012 using data from five shock-to-detonation transition (SDT) gas gun experiments by Burns and Chiquete, detonation front shape data from rate stick experiments and cylinder test experiment reported in Anderson et al., Jackson et al. The SDT experiments provide Pop plot data points (run distance-to-detonation) which characterize shock initiation, Lagrangian velocity time histories from embedded magnetic velocity gauges which give information on the reaction behind the lead shock, and Hugoniot data points for the reactants EOS. The rate stick experiments provide data on the curvature effect which characterize propagating detonation waves. The cylinder test experiments provide data on the detonation release isentrope used to calibrate the products EOS. The HE model uses a Davis reactants EOS with parameters that Burns and Chiquete [2020] fit to Hugoniot data, and Davis products EOS with parameters that Anderson et al. [2020] fit to cylinder test data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2021

Oregon Roadway CSZ Liquefaction and Landslide Impact Screening Analysis

The roadway transportation routes and the surrounding environments in Oregon will be subject to ground acceleration and potential associated deformation due to the propagation of waves of energy during a Cascadia Subduction Zone (CSZ) earthquake. Linear features such as highways may be subject to partial or complete disruption of their function resulting from ground deformations at and along the roadway. A potential source of these disruptions associated with a CSZ event are landslides and liquefaction in the soils along and beneath the highways and local roadways. The construction of the Oregon roadway system and urban development have occurred largely along the banks of the rivers—those along the Pacific coast, and along the valleys throughout the Oregon Coast Range and Cascade Mountains. These locations also represent areas where landslides and liquefaction have occurred historically, or may potentially occur in the future, presenting a significant vulnerability for the roadway system and dependent post-CSZ earthquake emergency response and recovery effort. Many of the Oregon areas with landslide- and liquefaction-prone soils have the potential to impact roadway transportation routes after a major earthquake causing delays in establishing routes to transport emergency response personnel and resources to impacted areas and communities. This report presents methods used to project the impacts that CSZ earthquake-induced landslides and liquefaction may have on roadways across Oregon. First, this report presents an approach for estimating the risk of landslides to Oregon roadways, and determining approximate reopening times. Second, this report presents a similar approach for estimating the risk of liquefaction-induced ground displacements to roadway pavements, also resulting in approximate reopening times for roadway segments. The outcomes of these two ground failure analyses serve as inputs to a broader analysis of CSZ earthquake impacts to Oregon’s surface transportation system and its ability to support post-disaster response and recovery activities, as part of the Oregon Transportation Systems Regional Resiliency Assessment Program project (CISA 2021).

58 GEOSCIENCES↗