Free oscillations of new lunar models
Free oscillations and surface wave dispersions from lunar models based on Orbiter data and corresponding pressure region of earth crust
Derr, J. S.↗
SEARCH · Engineering Papers
Results for “wave modeling”
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.
Physics of the Seasonal Sea Ice Zone
The seasonal sea ice zone encompasses the region between the winter maximum and summer minimum sea ice extent. In both the Arctic and Antarctic, the majority of the ice cover can now be classified as seasonal. Here, we review the sea ice physics that governs the evolution of seasonal sea ice in the Arctic and Antarctic, spanning sea ice growth, melt, and dynamics and including interactions with ocean surface waves as well as other coupled processes. The advent of coupled wave–ice modeling and discrete-element modeling, together with improved and expanded satellite observations and field campaigns, has yielded advances in process understanding. Many topics remain in need of further investigation, including rheologies appropriate for seasonal sea ice, wave-induced sea ice fracture, welding for sea ice freeze-up, and the distribution of snow on seasonal sea ice. Future research should aim to redress biases (such as disparities in focus between the Arctic and Antarctic and between summer and winter processes) and connect observations to modeling across spatial scales.
sea ice↗
Photochemical damping of inertio-gravity waves
A simple model is developed for the coupling effect of radiative heating and ozone photochemistry on inertio-gravity waves in the region 16-70 km. It is found that the photochemical damping rate is the weighted sum of the rates caused by two individual processes. One, suggested by Craig and Ohring, involves the temperature dependence of ozone reaction rate coefficients; the other, suggested by Leovy, involves the vertical gradient of the ozone mixing ratio. The relative importance of these two processes depends on the ratio of dynamical and chemical relaxation times. Calculations show that the most important region for coupling is near 48 km where the total relaxation time is 8 days for equinox and 6 days for summer midlatitude.
Zhu, Xun↗
Numerical simulations of convectively generated stratospheric gravity waves
A 2D model of a mesoscale convective flow is used to simulate the excitation and vertical propagation of gravity waves. Data obtained show that, in the absence of storm-relative mean winds in the stratosphere, the primary mode of excitation of gravity waves is by mechanical forcing owing to oscillatory updrafts. The stratospheric response consists of waves whose periods match the primary periods of the forcing. Due to the tendency of the oscillating updrafts to propagate toward the rear of the storm, gravity wave propagation is limited primarily to the rearward direction. Results suggest that squall-line-generated gravity waves arise from mechanical forcing rather than thermal effects.
Fovell, R.↗
Measurements of wave height statistics and radar, cross-section in a wind wave tank
There is currently wide interest among oceanographers and meteorologists in remote sensing of ocean surface characteristics. A wind wave tank developed at Langley Research Center is used to evaluate various remote sensing techniques based on electromagnetic scattering phenomena, and in the development and evaluation of theoretical scattering models. The wave tank is described, the statistics of the rough water surface are documented, and microwave radar cross-section measurement results are presented. The water surface statistics are similar in key respects to the open ocean, and the microwave scattering measurements show, qualitatively, theoretically predicted large and small scale scattering effects.
Johnson, J. W.↗
Numerical simulation of synthetic aperture radar image spectra for ocean waves
A numerical model for predicting the synthetic aperture radar (SAR) image of a moving ocean surface is described, and results are presented for two SIR-B data sets collected off the coast of Chile. Wave height spectra measured by the NASA radar ocean wave spectrometer (ROWS) were used as inputs to this model, and results are compared with actual SIR-B image spectra from orbits 91 and 106. Additional parametric variations are presented to illustrate the effects of nonlinearities in the imaging process.
Lyzenga, D. R.↗
Quantifying uncertainties due to optical potentials in one-neutron knockout reactions
One-neutron knockout reactions have been widely used to extract information about the single-particle structure of nuclei from the valley of stability to the drip lines. The interpretation of knockout data relies on reaction models, where the uncertainties are typically not accounted for. In this work, we quantify uncertainties of optical potentials used in these reaction models and propagate them, for the first time, to knockout observables using a Bayesian analysis. We study two reactions in the present paper, the first of which involves a loosely bound halo projectile, 11 Be, and the second a tightly bound projectile, 12 C. We first quantify the parametric uncertainties associated with phenomenological optical potentials. Complementing this approach, we also quantify the model uncertainties associated with the chiral forces that can be used to construct microscopic optical potentials. For the phenomenological study, we investigate the impact of the imaginary terms of the optical potential on the breakup and stripping components of the knockout cross sections as well as the impact of the angular range. For the 11 Be case, the theoretical uncertainty from the phenomenological method is on the order of the experimental uncertainty of the knockout observables; however, for the 12 C case, the theoretical uncertainty is significantly larger. Additionally, the widths of the uncertainty bands for the knockout observables obtained for the microscopic study and the phenomenological approach are of similar orders of magnitude. Based on this work we conclude that structure information inferred from the ratio of the knockout cross sections will carry a theoretical uncertainty of at least 20% for halo nuclei and at least 40% for tightly bound nuclei.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Computational Modeling of Semiconductor Dynamics at Femtosecond Time Scales
The Interchange No. NCC2-5149 deals with the emerging technology of photonic (or optoelectronic) integrated circuits (PICs or OEICs). In PICs, optical and electronic components are grown together on the same chip. To build such devices and subsystems, one needs to model the entire chip. PICs are useful for building components for integrated optical transmitters, integrated optical receivers, optical data storage systems, optical interconnects, and optical computers. For example, the current commercial rate for optical data transmission is 2.5 gigabits per second, whereas the use of shorter pulses to improve optical transmission rates would yield an increase of 400 to 1000 times. The improved optical data transmitters would be used in telecommunications networks and computer local-area networks. Also, these components can be applied to activities in space, such as satellite to satellite communications, when the data transmissions are made at optical frequencies. The research project consisted of developing accurate computer modeling of electromagnetic wave propagation in semiconductors. Such modeling is necessary for the successful development of PICs. More specifically, these computer codes would enable the modeling of such devices, including their subsystems, such as semiconductor lasers and semiconductor amplifiers in which there is femtosecond pulse propagation. Presently, there are no computer codes that could provide this modeling. Current codes do not solve the full vector, nonlinear, Maxwell's equations, which are required for these short pulses and also current codes do not solve the semiconductor Bloch equations, which are required to accurately describe the material's interaction with femtosecond pulses. The research performed under NCC2-5149 solves the combined Maxwell's and Bloch's equations.
Agrawal, Govind P.↗
Influence of transient baroclinic eddies on planetary-scale waves
Low-order models with a mean zonal flow forced by zonally symmetric heating and one planetary-scale wave forced by topography have multiple flow equilibria. Two stable equilibria are characterized by strong zonal flow with low wave amplitude (high index circulation) and weak zonal flow with a high wave amplitude (low index circulation) fixed in phase with the topography. These two states presumably represent normal zonal circulation and a blocking configuration, respectively. When a shorter, baroclinically unstable wave is dded to this low-order model, the planetary-scale wave no longer stays in stable equilibrium states. Instead, the long wave remains in weather regimes or preferred regions in phase space. These regimes are, in general, different from the equilibria of the model with the planetary-scale wave only. Thus, the short unstable wave adds some randomness to planetary-scale circulation and changes its position with respect to the planetary-scale topography. The Reinhold and Pierrehumbert model is further explored here by adding a long-wave in the thermal forcing and exploring wider parameter space. When the symmetric thermal forcing is weak so that the short wave is stable, the planetary scale wave is fixed in amplitude and phase by the asymmetric forcing. However, when the symmetric forcing is increased to an unstable level, the planetary-scale wave becomes less organized. With asymmetric thermal forcing only (i.e., no topography) the long-wave is randomly distributed through phase space. Thermal forcing seems to be less effective in organizing weather regimes than topographic forcing.
Branscome, L. E.↗
Sequence modeling of higher-order wave modes of quasi-circular, spinning, non-precessing binary black hole mergers
Higher-order gravitational wave modes from quasi-circular, spinning, non-precessing binary-black-hole (BBH) mergers encode rich information about the nonlinear dynamics of strong-field gravity. We present a transformer-based sequence-completion surrogate that, given an early-inspiral segment, forecasts the subsequent late inspiral, merger, and ringdown. The intended applications are (i) patching or completing expensive or interrupted numerical-relativity (NR) simulations and (ii) providing late-time cross-checks and rapid hybridization studies. The training set is built from the NRHybSur3dq8 surrogate, which provides spherical-harmonic modes up to $\ell$ ≤ 4 (excluding (4, 0) and (4,±1), and including (5, 5)) for mass ratios q ≤ 8, dimensionless spin components s$^{z}_{1,2}$ ϵ[–0.8, 0.8], and inclination angles θ ϵ [0, π]. Waveforms are supplied on the interval t ϵ [–5000M, –100 M) and the model autoregressively generates the plus and cross polarizations (h + , h x ) on t ϵ [–100 M, 130M]. Training on the Delta supercomputer with 16 NVIDIA A100 GPUs required ~15 h on more than 14 million hybrid waveforms. Evaluation on a held-out test set of 840,000 samples yields mean and median overlaps of 0.996 and 0.997, respectively, with respect to the surrogate ground truth.
black-hole merger↗
Revised NASA axially symmetric ring model for coupled-cavity traveling-wave tubes
A versatile large-signal, two-dimensional computer program is used by NASA to model coupled-cavity travelling-wave tubes (TWTs). In this model, the electron beam is divided into a series of disks, each of which is further divided into axially symmetric rings which can expand and contract. The trajectories of the electron rings and the radiofrequency (RF) fields are determined from the calculated axial and radial space-charge, RF, and magnetic forces as the rings pass through a sequence of cavities. By varying electrical and geometric properties of individual cavities, the model is capable of simulating severs, velocity tapers, and voltage jumps. The calculated electron ring trajectories can be used in designing magnetic focusing and multidepressed collectors. The details of using the program are presented, and results are compared with experimental data.
Wilson, Jeffrey D.↗
Graph-learning approach to combine multiresolution seismic velocity models
SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.
Geochemistry & Geophysics↗
Seismic properties of the crust and uppermost mantle of North America
Seismic refraction profiles for the North American continent were compiled. The crustal models compiled data on the upper mantle seismic velocity (P sub n), the crustal thickness (H sub c) and the average seismic velocity of the crystalline crust (V sub p). Compressional wave parameters were compared with shear wave data derived from surface wave dispersion models and indicate an average value for Poisson's ratio of 0.252 for the crust and of 0.273 for the uppermost mantle. Contour maps illustrate lateral variations in crustal thickness, upper mantle velocity and average seismic velocity of the crystalline crust. The distribution of seismic parameters are compared with a smoothed free air anomaly map of North America and indicate that a complidated mechanism of isostatic compensation exists for the North American continent. Several features on the seismic contour maps also correlate with regional magnetic anomalies.
Braile, L. W.↗
Weakly nonlinear models for turbulent mixing in a plane mixing layer
New closure models for turbulent free shear flows are presented in this paper. They are based on a weakly nonlinear theory with a description of the dominant large-scale structures as instability waves. Two models are presented that describe the evolution of the free shear flows in terms of the time-averaged mean flow and the dominant large-scale turbulent structure. The local characteristics of the large-scale motions are described using linear theory. Their amplitude is determined from an energy integral analysis. The models have been applied to the study of an incompressible mixing layer. For both models, predictions of the mean flow developed are made. In the second model, predictions of the time-dependent motion of the large-scale structures in the mixing layer are made. The predictions show good agreement with experimental observations.
Liou, William W.↗
Two-component model of strong Langmuir turbulence - Scalings, spectra, and statistics of Langmuir waves
A simple two-component model of strong turbulence that makes clear predictions for the scalings, spectra, and statistics of Langmuir waves is developed. Scalings of quantities such as energy density, power input, dissipation power wave collapse, and number density of collapsing objects are investigated in detail and found to agree well with model predictions. The nucleation model of wave-packet formation is strongly supported by the results. Nucleation proceeds with energy flowing from background to localized states even in the absence of a driver. Modulational instabilities play little or no role in maintaining the turbulent state when significant density nonuniformities are present.
Robinson, P. A.↗
Limits on the Range and Rate of Change in Power Take-Off Load in Ocean Wave Energy Conversion: A Study Using Model Predictive Control
Previous work comparing power take-off (PTO) architectures for ocean wave-powered reverse osmosis suggests that variable displacement in the wave energy converter (WEC)-driven pump does not offer a significant performance advantage. A limitation of that study is that the WEC was subject to a constant load within a given sea state (“Coulomb damping”) and did not account for controlled, moment-to-moment variation of the PTO load enabled by a variable displacement pump. This study explores the potential performance advantage of a variable PTO load over Coulomb damping. Model predictive control is used to provide optimal load control with constraints on the PTO load. The constraints include minimum and maximum loads and a limit on the rate of load adjustment. Parameter studies on these constraints enable conclusions about PTO design requirements in addition to providing an estimated performance advantage over Coulomb damping. Numerical simulation of the Oyster 1 WEC is carried out with performance weighted by historical sea state data from Humboldt Bay, CA. The results show a performance advantage of up to 20% higher yearly-average power absorption over Coulomb damping. Additionally, the parameter studies suggest that the PTO load should be adjustable down to at least 25% of the maximum load and should be adjustable between the minimum and maximum loads within a few seconds.
16 TIDAL AND WAVE POWER↗
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)↗
Imaging the crust and uppermost mantle structure of Portugal (West Iberia) with seismic ambient noise
SUMMARY We present a new high-resolution 3-D shear wave velocity (Vs) model of the crust and uppermost mantle beneath Portugal, inferred from ambient seismic noise tomography. We use broad-band seismic data from a dense temporary deployment covering the entire Portuguese mainland between 2010 and 2012 in the scope of the WILAS project. Vertical component data are processed using phase correlation and phase weighted stack to obtain empirical Green functions (EGFs) for 2016 station pairs. Further, we use a random sampling and subset stacking strategy to measure robust Rayleigh-wave group velocities in the period range 7–30 s and associated uncertainties. The tomographic inversion is performed in two steps: First, we determine group-velocity lateral variations for each period. Next, we invert them at each grid point using a new trans-dimensional inversion scheme to obtain the 3-D shear wave velocity model. The final 3-D model extends from the upper crust (5 km) down to the uppermost mantle (60 km) and has a lateral resolution of ∼50 km. In the upper and middle crusts, the Vs anomaly pattern matches the tectonic units of the Variscan Massif and Alpine basins. The transition between the Lusitanian Basin and the Ossa Morena Zone is marked by a contrast between moderate- and high-velocity anomalies, in addition to two arched earthquake lineations. Some faults, namely, the Manteigas–Vilariça–Bragança fault and the Porto–Tomar–Ferreira do Alentejo fault, have a clear signature from the upper crust down to the uppermost mantle (60 km). Our 3-D shear wave velocity model offers new insights into the continuation of the main tectonic units at depth and contributes to better understanding the seismicity of Portugal.
Silveira, Graça (ORCID:0000000221102554)↗