SEARCH · Search NASA
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.
A three dimensional direct field solver for the electric and magnetic fields using the Finite Difference Time Domain (FDTD) method is described. Example simulations using periodic and absorbing boundary conditions in a vacuum medium for a plane wave propagating in three dimensions are explored and are compared to theory.
These are the raw ultrasonic waveform files and nanoindentation measurements for three additively manufactured 316L stainless steel components with different fabrication parameters. Each part's length is divided into four segmented regions where their fabrication energy densities change. Part V+ begins at 33 J/mm^3 and increases in energy density by 3 J/mm^3 with each segment. Part C remains at a constant 33 J/mm^3 throughout the part. Part V- begins at 33 J/mm^3 and decreases in energy density by 3 J/mm^3 with each segment. We have found that the ultrasonic waves in regions with higher energy densities will traverse more quickly, as shown by the shorter time of flight, and vice versa. Similarly, nanoindentation measurements of reduced modulus and hardness follow this trend. We have measured the direction in which the parts were fabricated (build direction) and the segmented regions that are perpendicular to the build direction (transverse directions). The build direction longitudinal wave velocity measurements are considerably lower than their transverse direction counterparts, indicating anisotropy between the two directions. In Part C, we have removed material from one of its surfaces and recorded its ultrasonic and nanoindentation measurements with each material removal iteration. Changes in the material properties are more prominent by nanoindentation, suggesting material heterogeneity between the part's surface and interior. The README.txt file has information on how to navigate the files and process the data.
Traveling Wave Protection based on the Internet of Things (TWP-IoT) is developed to enable ultra-resilient microgrids. TWP-IoT adopts a directional traveling wave approach that utilizes the voltage and current modal components obtained from Discrete Hilbert Transform to identify the fault features via ultra-lightweight IoT hardware. Furthermore, the main contributions of this work include (1) a novel TWP-IoT for microgrid and distribution networks that can achieve the time-of-arrival fault detection and localization with composite wave impedance and conductance using Discrete Hilbert Transform, (2) a cost-effective smart IoT prototype with functionalities outperforming commercial relays, and (3) a real-time hardware-in-the-loop prototype based on RTDS that verifies the robustness and efficacy of TWP-IoT subject to a diverse set of working conditions. TWP-IoT is found to have excellent performance under a wide spectrum of internal and external fault conditions across a myriad of microgrid locations, which is unattainable by today’s TWP products.
Here, the design, fabrication, and validation of a thin-film lithium niobate on insulator (LNOI) electro-optic (EO) time-domain terahertz (THz) wave detector is reported. LNOI offers unprecedented properties for the EO detection of freely propagating THz wave radiation pulses and transient electric fields because of the large EO coefficient of the material, engineering of the velocity matching of the THz wave and optical wave, and much reduced detector size. The proof-of-concept device is realized using thin-film lithium niobate optical waveguides forming a Mach–Zehnder interferometer with interferometer arms electrically poled in opposite directions. THz waves are coupled effectively to the fully dielectric device from free space without using antennas or plasmonics. The detection of THz waves with frequencies up to 800 GHz is successfully demonstrated. The detector allows for the detection of THz frequency electric fields up to 4.6 MV/m. The observed frequency response of the device agrees well with theoretical predictions.
Gravitational memory, which describes the permanent shift in the strain after the passage of gravitational waves, is directly related to Weinberg’s soft graviton theorems and the Bondi-Metzner-Sachs (BMS) symmetry group of asymptotically flat space-times. In this work, we provide an equivalent description of the phenomenon in local coordinates around gravitational wave detectors, such as transverse-traceless (TT) gauge. We show that gravitational memory is encoded in large residual diffeomorphisms in this gauge, which include time-dependent anisotropic spatial rescalings, and prove their equivalence to BMS transformations when translated to TT gauge. We then derive the associated Ward identities and associated soft theorems, for both scattering amplitudes and equal-time (in-in) correlation functions, and explicitly check their validity for planar gravitational waves. Furthermore, the in-in identities are recognized as the flat-space analog of the well-known inflationary consistency relations.
In this work, we present a simple and automated approach to estimate primary site-response resonance, layer thickness, and shear-wave velocity directly from a dispersion image for a layer over half-space problem. We demonstrate this for high-impedance boundary conditions that lie in the upper tens of meters. Our approach eliminates the need for time-consuming dispersion curve picking and 1D shear-wave velocity inversion for large data volumes that can capture velocity structure in profile. We highlight important relationships between dispersion characteristics and resonance parameters through synthetic modeling and field data acquired over Atlantic Coastal Plain sediments. In this environment, shallow soil conditions are critical to accurately estimate earthquake site response. We suggest that this image processing approach can be applied to a range of high-impedance conditions, at a range of scales, or can provide model constraints for more complex velocity structures.
Dynamically compressed materials in longitudinal waves are described by two physical models: hydrostatic pressure, with equal, normal, principal stresses or material uniaxially strained in the wave propagation direction. These models are disparate, so experimental comparisons and evaluations are important. Polycrystalline material in a state of hydrostatic pressure will have no eccentricity of x-ray-diffracted Debye–Scherrer rings. A general three-dimensional solution of Bragg diffracted x rays based on principal crystallographic strains in the compression wave was found. The distortion of x-ray diffraction beams has been used for strain measurements; the analysis developed incorporates a strained reciprocal lattice and the incident x-ray beam. Strain-distorted Polanyi surfaces form an annulus of compression with an ellipsoid of revolution in reciprocal space, which is intersected by Ewald's sphere for Bragg diffraction. The in situ measurements for strain describe nanosecond diffraction evaluated using two planes, (hkl) and (h′k′l′), both in the same crystallographic phase. Diffraction from Al, Ni, Na, and Invar quantifies the compression axial strains in these materials: The compression axial ratios are 0.65, 1.05, 0.88, and 1.58 at pressures of 291, 402, 409, and 367 GPa for the respective materials. Crystal structure transformations with homogeneous pressurized stresses, mandating equal normal strains, should not be anticipated to agree with heterogeneous, uniaxially strained, and sheared crystalline phases. Measurements support in-plane strains increasing with pressure, p , in fcc and hcp aluminum as ε 22 = −0.21 p /TPa.
Density matrix electronic structure theory is used in many quantum chemistry methods to “alleviate” the computational cost that arises from directly using wave functions. Although density matrix based methods are computationally more efficient than wave function based methods, significant computational effort is involved. Because the Schrödinger equation needs to be solved as an eigenvalue problem, the time-to-solution scales cubically with the system size in mean-field type approaches such as Hartree–Fock and density functional theory and is solved as many times in order to reach charge or field self-consistency. We hereby propose and study a method to compute the density matrix by using a quadratic unconstrained binary optimization (QUBO) solver. This method could be useful to solve the problem with quantum computers and, more specifically, quantum annealers. Our proposed approach is based on a direct construction of the density matrix using a QUBO eigensolver. We explore the main parameters of the algorithm focusing on precision and efficiency. We show that, while direct construction of the density matrix using a QUBO formulation is possible, the efficiency and precision have room for improvement. Moreover, calculations performed with quantum annealing on D-Wave’s new Advantage quantum computer are compared with results obtained with classical simulated annealing, further highlighting some problems of the proposed method. Finally, we also suggest alternative methods that could lead to a more efficient QUBO-based density matrix construction.
A photonic Bragg grating is a fundamental building block that reflects the direction of wave propagation through spatial phase modulation and can be implemented using sidewall corrugation. However, due to the asymmetric aspect ratio of a waveguide cross section, typical Bragg gratings exhibit a strong polarization sensitivity. Here, we show that photonic Bragg gratings with cladding asymmetry can enable polarization-independent notch filters by rotating input polarizations. Such Bragg gratings strongly couple transverse electric (TE) and transverse magnetic (TM) modes propagating in opposite directions, filtering the input signal and reflecting the rotated mode. We analyzed this polarization-rotating Bragg grating using the coupled-mode theory and experimentally demonstrated it on a silicon-on-insulator platform. Our device concept is simple to implement and compatible with other platforms, readily available as polarization transparent Bragg components.
This award supported an experimental study of the properties of the dust acoustic wave, low-frequency, l ongitudinal wave mode that propagates through the dust component of a dusty plasma system and is self-excited by the free energy from the ion streaming p ast the dust, in the presence of a high magnetic fiel d. In particular, this project measured the thermal state of dusty plasma system by measuring the dispersi on relation of the dust acoustic wave while varying (1) the strength and direction of the magnetic field relative to the propagation direction of the wave and (2) the neutral gas pressure. By exploiting the unique ability of the MDPX facility to generate a dusty plasma in a high magnetic field environment, it is possible to create a preferential ion flow along the magnetic field lines as the ions become magnetized. In this way, the magnetic field strength can be used to control t he strength of the ion flow through the dust cloud an d provide insight into the contribution that one of the physical mechanisms, a two-stream ion-dust instability, thought to be responsible for high temperatures that have been previously observed might have. This work was completed over two, one-week long runs o n the Magnetized Dusty Plasma eXperiment (MDPX) in the Magnetized Plasma Research Laboratory at Auburn University and supported by work at Wittenberg University (development of an electrode system that would allow for the wave mode to propagate in th e direction of gravity as the orientation of the magnetic field is changed by rotating the magnets on the M DPX device, analysis of data acquired during the experimental campaigns, testing alternative strategies to form a dust cloud based on what was learned the first experimental campaign, etc.).
Current commercial long-range guided wave inspection and monitoring technologies for pipes usually require establishment of certain specific acoustic modes for wave propagation and use a large number of sensors in a certain geometric arrangement (for example, uniformly-spaced arrays of sensors on a collar, placed around a pipe). This paper, in contrast, presents a long-range corrosion detection method in pipes based on multi-mode guided acoustic waves, generated and received using simple, single-element sensors. In this technique, available guided wave modes in any interrogated zone of the structure are automatically down-selected by the medium, and changes in energy distribution over the selected modes are analyzed through signal processing to monitor material loss or addition. Here, to investigate the feasibility of the proposed method, we have performed both experimental and numerical simulation studies on 50.8-mm (2-in)-diameter schedule-40 carbon steel pipes. In these studies, a linear chirp signal is injected in compression mode, perpendicular to the surface of the pipe, generating multi-mode, multi-directional guided waves that travel within the pipe wall. Corrosion effects in the pipes are simulated numerically and experimentally. To reduce errors in the measurement technique due to changes in the local or global temperature of the measured media, a universal temperature compensation algorithm is introduced. The experimental results show good agreement with simulation results, indicating that the proposed multi-mode wave generation and detection in a pipe could be used for long-range detection and monitoring of corrosion effects, such as pipe wall losses or solid material deposits.
Efforts to streamline and codify wave energy resource characterization and assessment for regional energy planning and wave energy converter (WEC) project development have motivated the recent development of resource classification systems. Given the unique interplay between WEC absorption and resource attributes, viz, available wave power frequency, directionality, and seasonality, various consensus resource classification metrics have been introduced. However, the main international standards body for the wave energy industry has not reached consensus on a wave energy resource classification system designed with clear goals to facilitate resource assessment, regional energy planning, project site selection, project feasibility studies, and selection of WEC concepts or archetypes that are most suitable for a given wave energy climate. In this work, a primary consideration of wave energy generation is the available energy that can be captured by WECs with different resonant frequency and directional bandwidths. Therefore, the proposed classification system considers combinations of three different wave power classifications: the total wave power, the frequency-constrained wave power, and the frequency-directionally constrained wave power. The dominant wave period bands containing the most wave power are sub-classification parameters that provide useful information for designing frequency and directionally constrained WECs. The bulk of the global wave energy resource is divided into just 22 resource classes representing distinct wave energy climates that could serve as a common language and reference framework for wave energy resource assessment if codified within international standards.
This numerical study compares the wave field generated by the spectral wave action balance code, SNL-SWAN, to the linear-wave boundary-element method (BEM) code, WAMIT. The objective of this study is to assess the performance of SNL-SWAN for modeling wave field effects produced by individual wave energy converters (WECs) and wave farms comprising multiple WECs by comparing results from SNL-SWAN with those produced by the BEM code WAMIT. BEM codes better model the physics of wave-body interactions and thus simulate a more accurate near-field wave field than spectral codes. In SNL-SWAN, the wave field’s energy extraction is modeled parametrically based on the WEC’s power curve. The comparison between SNL-SWAN and WAMIT is made over a range of incident wave conditions, including short-, medium-, and long-wavelength waves with various amounts of directional spreading, and for three WEC archetypes: a point absorber (PA), a pitching flap (PF) terminator, and a hinged raft (HR) attenuator. Individual WECs and wave farms of five WECs in various configuration were studied with qualitative comparisons made of wave height and spectra at specific locations, and quantitative comparisons of the wave fields over circular arcs around the WECs as a function of radial distance. Results from this numerical study demonstrate that in the near-field, the difference between SNL-SWAN and WAMIT is relatively large (between 20% and 50%), but in the far-field from the array the differences are minimal (between 1% and 5%). The resultant wave field generated by the two different numerical approaches is highly dependent on parameters such as: directional wave spreading, wave reflection or scattering, and the WEC’s power curve.
Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply-scattered or reflected waves. However, a challenge arises when multiply-scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The twin neural network with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. These results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.
Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply scattered or reflected waves. However, a challenge arises when multiply scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural or induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (TNN also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The TNN with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. Furthermore, these results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.
With ground-based interferometers detecting hundreds of gravitational-wave (GW) events, GW astronomy has continued to blossom. U.S. and European scientists are developing plans to construct third-generation ground-based interferometers. ESA has proceeded through the mission formulation phase of space-based interferometer in a lower-frequency band, 10 –4 –0.1 Hz. Despite all these exciting developments, there is still a missing frequency band, 0.1–10 Hz. This mid-frequency band is rich with interesting astrophysical events. Coalescence and merger of intermediate-mass black holes (IMBHs) will occur in this frequency band. Coalescing stellar-mass BHs will pass through this frequency band days before they reach the frequency band of LIGO and Virgo. Detection of such signals would enable a mid-frequency detector to issue an advance notice to the high-frequency GW detectors, as well as to optical, x-ray and γ-ray telescopes. We propose Lunar Accelerometer Network Gravitational Observatory (LANGO) to detect GWs in this frequency band. In the first phase (LANGO 1), we propose to deploy four ambient-temperature (250 K) accelerometers in the tetrahedral or in a square configuration on the hemisphere facing the Earth. After successful operation at 250 K, LANGO would be upgraded to a cryogenic (4 K) version with over two orders of magnitude increased sensitivity (LANGO 2), with coherent rejection of seismic noise implemented. LANGO is a full-tensor detector, capable of determining the source direction and wave polarization. Each test mass (TM) is suspended as a pendulum with resonance frequency ∼ 0.01 Hz, thus is only weakly coupled to the lunar surface horizontally. LANGO is designed to detect the relative motion of globally separated, nearly free, TMs by using the Moon as a large quiet platform. LANGO 1 and 2 accelerometers aim at sensitivities ⩽ 10 –11 m s –2 Hz –1/2 and ⩽ 10 –13 m s –2 Hz –1/2 in the horizontal axes over the frequency band of 1 mHz–10 Hz, which yield GW sensitivities 1.1 x 10 -21 Hz -1/2 and 3.8 x 10 -24 Hz -1/2 at 1 Hz, respectively. The LANGO accelerometers will be 10 3 –10 5 times more sensitive than Apollo seismometers. With such sensitivity, LANGO will also make great contribution to the advancement of lunar geophysics.
Abstract The northwestern part of the Gulf of Mexico has undergone two episodes of continental rifting and collision and produced structural artifacts that are now buried under many kilometers of sediments, complicating investigations of the region. The deep sedimentary package precludes outcrops and points to a need for the application of seismic techniques, but low rates of seismicity in the region and sparse seismic monitoring limit the utility of traditional seismic methods. We therefore use diverse data to perform two‐dimensional seismic tomography across the dry land portion of the margin. Data are gleaned from teleseismic and regional earthquakes, postcritical SsPmp arrivals, and direct P wave energy identified with seismic interferometry that were recorded by a broadband, three‐component array and partially overlapping short‐period, vertical‐component array. The Pn and postcritical SsPmp phase help constrain the Moho discontinuity, which a previous receiver function study suggested was absent beneath the seaward portion of this transect. A high‐velocity body is observed in the crust at the same location as the Houston Magnetic Anomaly, possibly marking rocks from the Alleghenian continental assembly. The crust thins from NW to SE, indicating that extension occurred mostly to the south of the Ouachita orogeny. Our model indicates that the margin's sediment package reaches a maximum thickness of ~15 km at the coast and becomes unresolvably thin near the Llano Uplift.
Abstract On September 12, 2018 and April 14, 2019 shallow (<2.5 km from the surface, <1 km relative to sea level) earthquakes ( M w 4.0 and 4.1, respectively) struck a relatively aseismic area in south‐central Utah, within the Black Rock Desert volcanic field. Seismic waveforms are dominated by long‐period energy, no clear direct S‐wave arrivals, and signal durations of ∼100 s. Using co‐seismic interferograms, we detect clear and distinguishable surface deformation confirming the shallowness of the two mainshocks, while full moment tensor solutions rule out an isotropic source, such as an explosion or collapse. For the April 2019 sequence, we build an earthquake catalog using a dense temporary geophone array. The enhanced earthquake catalog shows a rapidly decaying aftershock sequence, with 85% of aftershocks occurring during the first three hours. Combining these observations, we argue that these earthquake sequences are related to the Black Rock Desert volcanic field.