SEARCH · Engineering Papers

Results for “constrained wave”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Physics-constrained Gaussian process model for prediction of hydrodynamic interactions between wave energy converters in an array

To improve the efficiency of wave farms and achieve maximum power generation, the layout of wave energy converters (WECs) in an array needs to be carefully designed so that the hydrodynamic interactions can be positively exploited. For this, the hydrodynamic characteristics of the WEC array in different layouts need to be calculated. However, such calculations using numerical models usually entail significant computational cost, especially for large arrays of WECs. To address the computational challenge, a physics-constrained Gaussian process (GP) model is proposed to replace the original expensive numerical model and predict the hydrodynamic characteristics of the WECs for any array layout. By exploring the relationship between the WEC array (i.e., the input) and different hydrodynamic characteristics (i.e., the output), here we summarize a set of physical constraints/features, including invariance, symmetry, and additivity. This prior knowledge about the input-output relationship is then directly embedded in the constructed GP model through the design of physics-constrained kernels. In particular, a double-sum invariant kernel is first developed to incorporate the invariance and symmetry features, and then an additive kernel is developed to incorporate the additive feature of the problem. The invariant kernel and the additive kernel are then integrated to construct the physics-constrained GP model. Compared to the standard GP model, the proposed physics-constrained GP models require less training data to achieve the desired accuracy in predicting the hydrodynamic characteristics and are also less vulnerable to the curse of dimensionality (i.e., good scalability for large arrays) due to the use of an additive kernel. The efficiency, accuracy, and scalability of the proposed approach are demonstrated through an application to predict the hydrodynamic characteristics for WEC arrays of different sizes and layouts.

16 TIDAL AND WAVE POWER↗

DOE OSTI2023

Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations

SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.

Pinzón, Juan I. (ORCID:0000000338078068)↗

NASA NTRS2011

Deep Moonquake Focal Mechanisms: Recovery and Implications

A defining characteristic of deep moonquakes is their tendency to occur with tidal periodicity, prompting previous studies to infer that they are related to the buildup and release of tidal stress within the Moon. In studies of tidal forcing, a key constraint is the focal mechanism: the fault parameters describing the type of failure moonquakes represent. The quality of the lunar seismic data and the limited source/receiver geometries of the Apollo seismic network prohibit the determination of deep moonquake fault parameters using first-motion polarities, as is typically done in terrestrial seismology. Without being able to resolve tidal stress onto a known failure plane, we can examine only gross qualities of the tidal stress tensor with respect to moonquake occurrence, so we cannot fully address the role of tidal stress in moonquake generation. We will examine the extent to which shear (S) and compression (P) wave amplitude ratios can constrain moonquake fault geometry by determining whether, for a given cluster, there exists a focal mechanism that can produce a radiation pattern consistent with the amplitudes measured by the Apollo instruments. Amplitudes are read in the ray coordinate frame, directly from seismograms for which the P and S arrivals are clearly identifiable on all long-period channels of the four Apollo stations. We apply an empirical station correction to account for site effects and the differences between P- and S-wave attenuation. Instead of focusing on the best fitting solution only, we formulate the inverse problem using a falsification criterion: all source orientations that do not reproduce the observed SV/P ratios within an error margin derived from the uncertainty of amplitude readings are rejected. All others are accepted as possible solutions. The inversion is carried out using an exhaustive grid search on a regular grid with predefined step size, encompassing all possible combinations of strike, dip and slip. To assess the sensitivity of the inversion for the uncertainty of the lunar interior structure, we carry out repeated inversions with different velocity structures. Our data set consist of a total of 106 events from 25 deep moonquake clusters. The largest contribution of 37 events originates from the most active cluster, A001, while other clusters are represented by 1 to 9 events. Since the definition of a cluster implies that all events share the same source orientation, a comparison of the inversion results of all events from one cluster will reduce ambiguities of the inversion. Once we obtain a suite of fault parameters for a given source, we can attempt to further constrain the focal mechanism with refined analyses of tidal stresses and predictions based on synthetic seismograms.

Knapmeyer, Martin↗

DOE OSTI2025

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES↗

DOE OSTI2023

Constraining properties of asymmetric dark matter candidates from gravitational-wave observations

The accumulation of certain types of dark matter particles in neutron star cores due to accretion over long timescales can lead to the formation of a mini black hole. In this scenario, the neutron star is destabilized and implodes to form a black hole without significantly increasing its mass. When this process occurs in neutron stars in coalescing binaries, one or both stars might be converted to a black hole before they merge. Thus, in the mass range of ~1–2M ⊙ , the Universe might contain three distinct populations of compact binaries: one containing only neutron stars, the second population of only black holes, and a third, mixed population consisting of a neutron star and a black hole. However, it is unlikely to have a mixed population as the various timescales allow for both neutron stars to remain or collapse within a short timescale. In this paper, we explore the capability of future gravitational-wave detector networks, including upgrades of Advanced LIGO and Virgo, and new facilities such as the Cosmic Explorer and Einstein Telescope (XG network), to discriminate between different populations by measuring the effective tidal deformability of the binary, which is zero for binary black holes but nonzero for binary neutron stars. Furthermore, we show that observing the relative abundances of the different populations can be used to infer the timescale for neutron stars to implode into black holes, and in turn, provide constraints on the particle nature of dark matter. Furthermore, the XG network will infer the implosion timescale to within an accuracy of 0.01 Gyr at 90% credible interval and determine the dark matter mass and interaction cross section to within a factor of 2 GeV and 10 cm –2 , respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1989

Stochastic acceleration of cometary pickup ions - The classic leaky box model

The acceleration of cometary pickup ions by magnetohydrodynamic waves at P/Giacobini-Zinner is examined in a model where acceleration predominantly occurs downstream of the bow shock throughout the cometosheath where intense magnetic turbulence exists. The mean free path for scattering by the magnetic fluctuations in this region is less than the characteristic dimension of the cometosheath so that pickup ions are rendered isotropic and energized by a modest amount in the process. This principal loss mechanism for the ions is spatial diffusion out of the acceleration region moderated by the self-same accelerating waves. This particular feature constrains the model in a way that the predicted ion spectrum is uniquely determined by the power spectrum of the magnetic turbulence. At both P/Giacobini-Zinner and P/Halley, the turbulence is non-Kolmogoroff with a spectral index of 2 resulting in an ion spectral behavior that is approximately an exponential in ion speed, consistent with a recent analysis of the Giacobini-Zinner data (see Richardson et al.).

Barbosa, D. D.↗

NASA NTRS2020

Differences in Tropical High Clouds Among Reanalyses: Origins and Radiative Impacts

We examine differences among reanalysis highcloud products in the tropics, assess the impacts of these differences on radiation budgets at the top of the atmosphere and within the tropical upper troposphere and lower stratosphere (UTLS), and discuss their possible origins in the context of the reanalysis models. We focus on the ERA5 (fifthgeneration European Centre for Medium-range Weather Forecasts – ECMWF – reanalysis), ERA-Interim (ECMWF Interim Reanalysis), JRA-55 (Japanese 55-year Reanalysis), MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2), and CFSR/CFSv2 (Climate Forecast System Reanalysis/Climate Forecast System Version 2) reanalyses. As a general rule, JRA-55 produces the smallest tropical high-cloud fractions and cloud water contents among the reanalyses, while MERRA-2 produces the largest. Accordingly, long-wave cloud radiative effects are relatively weak in JRA-55 and relatively strong in MERRA-2. Only MERRA-2 and ERA5 among the reanalyses produce tropical-mean values of outgoing long-wave radiation (OLR) close to those observed, but ERA5 tends to underestimate cloud effects, while MERRA-2 tends to overestimate variability. ERA5 also produces distributions of longwave, short-wave, and total cloud radiative effects at the top of the atmosphere that are very consistent with those observed. The other reanalyses all exhibit substantial biases in at least one of these metrics, although compensation between the long-wave and short-wave effects helps to constrain biases in the total cloud radiative effect for most reanalyses. The vertical distribution of cloud water content emerges as a key difference between ERA-Interim and other reanalyses. Whereas ERA-Interim shows a monotonic decrease of cloud water content with increasing height, the other reanalyses all produce distinct anvil layers. The latter is in better agreement with observations and yields very different profiles of radiative heating in the UTLS. For example, whereas the altitude of the level of zero net radiative heating tends to be lower in convective regions than in the rest of the tropics in ERAInterim, the opposite is true for the other four reanalyses. Differences in cloud water content also help to explain systematic differences in radiative heating in the tropical lower stratosphere among the reanalyses. We discuss several ways in which aspects of the cloud and convection schemes impact the tropical environment. Discrepancies in the vertical profiles of temperature and specific humidity in convective regions are particularly noteworthy, as these variables are directly constrained by data assimilation, are widely used, and feed back to convective behaviour through their relationships with thermodynamic stability.

High Clouds↗

DOE OSTI2024

Characterizing Interference in a Wave-Powered Buoy

Oceanographic measurement buoys are power constrained by the capacity of their batteries and the available solar power. Many buoys are in locations with abundant wave energy, but wave energy converters are not being used to power them. One reason wave energy converters may not have not been developed for wave measurement buoys is that many buoys measure waves by precisely following the wave motion. Wave energy converters may inhibit this wave-following ability by altering a buoy's response to incoming waves. In this report, we present the results of field tests that measured the extent to which buoy motion was changed by modifications to the buoy for wave energy harvesting. The motion of the buoy was significantly modified by the addition of the hardware components necessary for harvesting wave energy. The ratio of the mass of the wave energy converter to the mass of the buoy was high for the system tested, which strongly affected buoy motion. A lower value for this ratio would decrease the effect of the wave energy converter on the motion of the buoy, thereby decreasing the impact on the buoy's ability to measure waves but also decreasing the power production capabilities. Therefore, the use of wave energy converters to power wave measurement buoys may be feasible for large buoys that measure low-frequency waves and have low power requirements.

16 TIDAL AND WAVE POWER↗

NASA NTRS1984

Tropical response to lateral forcing with a latitudinally and zonally nonuniform basic state

The dynamics of lateral forcing of a stable tropical atmosphere by incident waves from the midlatitudes are examined. Attention is focused on conditions where the WKB method is invalid, and only given to barotropic flows. The tropical atmosphere is treated as a steady planetary wave and a zonal shear flow. The 200 mb level is investigated in terms of a midlatitude wave with self-compatible parametric values that is constrained to satisfy the radiation conditions at the boundary. Wave modes are found to be susceptible to trapping between two turning points in the equatorial zone in the presence of a zonal shear flow. Quasi-resonant waves can then appear over a wide range of parametric conditions of the forcing wave and the basic state. The zonal wind structure determines the resonance conditions.

Wilson, J. D.↗

NASA NTRS2012

Discovery Of A Rossby Wave In Jupiter's South Equatorial Region

A detailed study of the chevron-shaped dark spots on the strong southern equatorial wind jet near 7.5 deg S planetographic latitude shows variations in velocity with longitude and time. The chevrons move with velocities near the maximum wind jet velocity of approx.140 m/s, as deduced by the history of velocities at this latitude and the magnitude of the symmetric wind jet near 7 deg N latitude. Their repetitive nature is consistent with an inertia-gravity wave (n = 75-100) with phase speed up to 25 m/s, relative to the local flow, but the identity of this wave mode is not well constrained. However, high spatial resolution movies from Cassini images show that the chevrons oscillate in latitude with a approx.7-day period. This oscillating motion has a wavelength of approx.20 deg and a speed of approx.100 m/s, following a pattern similar to that seen in the Rossby wave plumes of the North Equatorial Zone, and possibly reinforced by it, though they are not perfectly in phase. The transient anticyclonic South Equatorial Disturbance (SED) may be a similar wave feature, but moves at slower velocity. All data show chevron latitude variability, but it is unclear if this Rossby wave is present during other epochs, without time series movies that fully delineate it. In the presence of multiple wave modes, the difference in dominant cloud appearance between 7 deg N and 7.5 deg S may be due to the presence of the Great Red Spot, either through changes in stratification and stability or by acting as a wave boundary.

Simon-Miller, Amy A.↗

NASA NTRS2018

Constraining the Models' Response of Tropical Clouds to SST Forcings Using CALIPSO Observations

Here we present preliminary results from the analysis of the low cloud cover (LCC) and cloud radiative effect (CRE) interannual changes in response to sea surface temperature (SST) forcings in two GISS climate models, and 12 other climate models. We further classify them as a function of their ability to reproduce the vertical structure of the cloud response to SST change against 10 years of CALIPSO observations: “the constrained models, which match the observation constraint, and the unconstrained models”. The constrained models replicate the observed interannual LCC change particularly well (ΔLCC(sub con)=-3.49 ±1.01 %/K vs. ΔLCC(sub obs)=-3.59 ±0.28 %/K) as opposed to the unconstrained models, which largely underestimate it (ΔLCC(sub unc) = -1.32 ± 1.28 %/K). As a result, the amount of short-wave warming simulated by the constrained models (ΔCRE(sub con)=2.60 ±1.13 W/m2/K) is in better agreement with the observations (ΔCRE(sub obs)=3.05 ±0.28 W/m2/K) than the unconstrained models (ΔCRE(sub con)=0.87 ±2.63 W/m2/K). Depending on the type of low cloud, the observed relationship between cloud/radiation and surface temperature varies. Over the stratocumulus regions, increasing SSTs generate higher cloud top height along with a large decrease of the cloud fraction below as opposed to a slight decrease of the cloud fraction at each level over the trade cumulus regions. Our results suggest that the models must generate sustainable stratocumulus decks and moist processes in the planetary boundary layer to reproduce these observed features. Future work will focus on defining a method to objectively discriminate these cloud types that can be applied consistently in both the observations and the models.

Cesana, Gregory↗

NASA NTRS2012

Longitudinal Variation and Waves in Jupiter's South Equatorial Wind Jet

A detailed study of the chevron-shaped dark spots on the strong southern equatorial wind jet near 7.5 S planetographic latitude shows variations in velocity with longitude and time. The presence of the large anticyclonic South Equatorial Disturbance (SED) has a profound effect on the chevron velocity, causing slower velocities to its east and accelerations over distance from the disturbance. The chevrons move with velocities near the maximum wind jet velocity of approx 140 m/s, as deduced by the history of velocities at this latitude and the magnitude of the symmetric wind jet near 7 N latitude. Their repetitive nature is consistent with a gravity-inertia wave (n = 75 to 100) with phase speed up to 25 m/s, relative to the local flow, but the identity of this wave mode is not well constrained. However, for the first time, high spatial resolution movies from Cassini images show that the chevrons oscillate in latitude with a 6.7 +/- 0.7-day period. This oscillating motion has a wavelength of approx 20 and a speed of 101 +/- 3 m/s, following a pattern similar to that seen in the Rossby wave plumes of the North Equatorial Zone, and possibly reinforced by it. All dates show chevron latitude variability, but it is unclear if this larger wave is present during other epochs, as there are no other suitable time series movies that fully delineate it. In the presence of mUltiple wave modes, the difference in dominant cloud appearance between 7 deg N and 7.5 deg S is likely due to the presence of the Great Red Spot, either through changes in stratification and stability or by acting as a wave boundary.

Simon-Miller, Amy A.↗

DOE OSTI2025

Mathematical Foundation for Quantum Computing of Electromagnetic Wave Propagation in Dielectric Media

Can quantum computers effectively simulate the propagation and scattering of electromagnetic waves in a classical plasma? This chapter introduces some of the basic concepts in mathematics and physics essential to answering that question. The numerical simulations of Maxwell equations for wave propagation in dielectrics are constrained by technological limitations of the present-day computers. In contrast, there has been ample fanfare around quantum computers and their potential to far exceed the performance of traditional computers. Whether the enhanced capabilities of a quantum computer can be put to use for simulating topics in classical physics is a source of intrigue and curiosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Monitoring Seasonal Shear Wave Velocity Changes in the Top 6 m at Garner Valley in Southern California With Borehole Data

Abstract Subsurface structures play important roles in seismic ground motion, crustal hydrology, stability of the built environment, and more. Constraining temporal changes of subsurface shear wave velocity ( V S ) can provide useful information to all these topics and the growing field of hydrological monitoring with seismic velocity. Using borehole records at Garner Valley, CA, we estimate seasonal subsurface V S variations from impulse response functions (IRFs) of earthquake data (2005–2018) along with IRFs and cross‐correlation of cross‐hole experiment data (2015–2018). The inferred V S variations are up to ∼25% in the top 6 m and ∼10% at 2–5 m in depth. The V S variations correlate strongly with the water table depth changes, suggesting that the changes are mostly due to fluctuations of pore pressure in the shallow material. The shallow velocity changes alter the near‐surface conditions, can affect seismic hazard estimation, and may be improperly attributed to deeper processes without careful analysis.

58 GEOSCIENCES↗

NASA NTRS2000

Modeling the Diurnal Tides in the MLT Region with the Doppler Spread Parameterization of Gravity Waves

Special Session: SA03 The mesosphere/lower thermosphere region: Structure, dynamics, composition, and emission. Ground based and satellite observations in the upper mesosphere and lower thermosphere (MLT) reveal large seasonal variations in the horizontal wind fields of the diurnal and semidiurnal tides. To provide an understanding of the observations, we discuss results obtained with our Numerical Spectral Model (NMS) that incorporates the gravity wave Doppler Spread Parameterization (DSP) of Hines. Our model reproduces many of the salient features observed, and we discuss numerical experiments that delineate the important processes involved. Gravity wave momentum deposition and the seasonal variations in the tidal excitation contribute primarily to produce the large equinoctial amplitude maxima in the diurnal tide. Gravity wave induced variations in eddy viscosity, not accounted for in the model, have been shown by Akmaev to be important too. For the semidiurnal tide, with amplitude maximum observed during winter solstice, these processes also contribute, but filtering by the mean zonal circulation is more important. A deficiency of our model is that it cannot reproduce the observed seasonal variations in the phase of the semidiurnal tide, and numerical experiments are being carried out to diagnose the cause and to alleviate this problem. The dynamical components of the upper mesosphere are tightly coupled through non-linear processes and wave filtering, and this may constrain the model and require it to reproduce in detail the observed phenomenology.

Mayr, H. G.↗

NASA NTRS1998

Shock Temperatures in Calcite at 115 to 185 GPa

We present optical pyrometer measurements of shock temperatures in crystalline calcite (CaC03) which constrain its behavior in shock waves with amplitudes of 115 to 185 GPa.

shock temperatures optical pyrometer crystalline c↗

NASA NTRS2017

The Nanograv Nine-Year Data Set: Measurement and Analysis of Variations in Dispersion Measures

We analyze dispersion measure(DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause.

pulsars: general↗

NASA NTRS2017

North American Drought and Links to Northern Eurasia: The Role of Stationary Rossby Waves

This chapter provides an overview of the role of stationary Rossby waves in the sub-seasonal development of warm season drought over North America and subsequent downstream development of climate anomalies over northern Eurasia. The results are based on a case study of a stationary Rossby wave event that developed during 20 May 15 June 1988. Simulations with the NASA Goddard Earth Observing System (GEOS-5) atmospheric general circulation model highlight the importance of the mean jet streams in guiding and constraining the path and speed of wave energy propagation. In particular, convective anomalies that developed over the western Pacific in late May (in the presence of the strong North Pacific jet) produce a predilection for persistent upper-level high anomalies over central North America about ten days later, leading to the rapid development of severe dry conditions there. There are indications of continued downstream wave energy propagation that reaches northern Eurasia about two weeks later, leading to the development of dry conditions over eastern Europe and western Russia, and cool and wet conditions over western Europe and central northern Eurasia. The results suggest that stationary Rossby waves can serve as a source of predictability for sub-seasonal development of droughts over North America and northern Eurasia.

GEOS-5↗