Search NASA⌕ Search

SEARCH · Search NASA

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.

At least 307 records · Page 17

Microscopic optical potentials from a Green's function approach

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.

Direct reactions↗

Charge singlets and orbital-selective charge density wave transitions

The possibility of “orbitally selective Mott transitions” within a multiband Hubbard model, in which one orbital with large on-site electron-electron repulsion U 1 is insulating and another orbital, to which it is hybridized, with small U –1 , is metallic, is a problem of long-standing debate and investigation. In this paper we study an analogous phenomenon, the coexistence of metallic and insulating bands in a system of orbitals with different electron-phonon coupling. To this end, we examine two variants of the bilayer Holstein model: a uniform bilayer and a “Holstein-metal interface” where the electron-phonon coupling, λ, is zero in the “metallic” layer. In the uniform bilayer Holstein model, charge density wave (CDW) order dominates at small interlayer hybridization t 3 , but decreases and eventually vanishes as t 3 grows, providing a charge analog of singlet (spin liquid) physics. In the interface case, we show that CDW order penetrates into the metal layer and forms long-range CDW order at an intermediate ratio of inter- to intralayer hopping strengths, 1.4 ≲ t 3 /t ≲ 3.4. Furthermore, this is consistent with the occurrence of an “orbitally selective CDW” regime at weak t3 in which the layer with λ 1 ≠ 0 exhibits long-range charge order, but the “metallic layer” with λ –1 = 0, to which it is hybridized, does not.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Frequency Domain Methodology for Quantitative Evaluation of Diffuse Wavefield With Applications to Seismic Imaging

Abstract Ambient Noise Imaging (ANI) of subsurface structures relies on seismic interferometry of diffuse seismic wavefields. However, the lack of effective methods to quantify and identify highly diffuse waves hampers applications of ANI, particularly in evaluating seismic attenuation and monitoring structural changes with high temporal resolution. Conventional ANI approaches require data normalization, which effectively suppresses the non‐diffuse component with large amplitude but also results in significant loss of amplitude and phase information in the continuous seismic records. In this study, we propose a frequency domain method to quantitatively evaluate the degree of diffuseness of seismic wavefields by analyzing their statistical characteristics of modal amplitudes for stationarity and randomness. Tests on synthetic waveform and field nodal records show that the proposed method can effectively distinguish between diffuse and non‐diffuse waveforms for either single‐ or three‐component data. As an application, we identify a 60‐s‐long diffuse coda of a local M 2.2 earthquake recorded by a dense nodal array on the San Jacinto Fault Zone, and successfully extract high‐quality dispersion curve andQ‐value without performing data normalization. These results are consistent with those obtained by conventional methods that assess the correlation between coherency and the Green's function, and by modeling ballistic waves generated by road traffic. Our proposed method can advance the imaging of subsurface velocity and attenuation structures as well as monitoring temporal changes for scientific studies and engineering applications.

Geochemistry & Geophysics↗

Uncertainty quantification of optical models in fission fragment deexcitation

Here, we take the first step towards incorporating compound nuclear observables at astrophysically relevant energies into the experimental evidence used to constrain optical models, by propagating the uncertainty in two global optical potentials, one phenomenological and one microscopic, to correlated fission observables using the Monte Carlo Hauser-Feshbach formalism. We compare to a wide range of historic and recent experimental fission measurements, and discuss in detail regions of disagreement. We find that the parametric optical model uncertainty in neutron-fragment correlated observables involving neutron energy is significant. On the other hand, we observe that other experimental features, particularly neutron-fragment correlations near the 132 Sn shell closure and the high energy component of neutron spectra, are unlikely to be explained by the optical potential, and will require further experimental and theoretical effort to explain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Elastic scattering in the 12 N + 197 Au system at 𝐸 lab = 70 MeV

Halo structure is an interesting exotic configuration developed in some light weakly bound nuclei, where a valence particle orbits a nuclear core. Signatures of halo structure can be observed in the angular distributions of the elastic scattering induced by these nuclei at energies around the Coulomb barrier. There are some well-studied reactions with neutron-rich halo nuclei, such as 6 He and 11 Li . However, the information is scarce on the proton-rich side. Recent works confirm the halo structure in the 8B nuclei but still lack more experimental studies for other proton-halo candidates, such as 12 N and 17 F. In this work, we report experimental data for the elastic scattering of 12 N on 197 Au target at 𝐸 lab = 70 MeV. The 12 N is a proton-rich nucleus with proton separation energy 𝑆 𝑝 = 600 keV, which is higher than the 8 B (𝑆 𝑝 = 137 keV) and almost the same as the 17 F (𝑆 𝑝 = 601 keV). Data were obtained at the Cyclotron Institute of Texas A&M University where the 12 N radioactive beam was produced by the momentum achromatic recoil spectrometer. The calculation of the optical model was used to fit the measured angular distribution for the elastic scattering and to obtain the reaction cross section 𝜎 𝑅 . We also performed continuum discretized coupled-channel calculations to compare with the experimental data. The angular distribution of the elastic cross sections exhibits a suppression of the Fresnel peak. From the fitting of the optical model, we obtain the total reaction cross section, 𝜎 𝑅 = 1269 ± 41 mb. The agreement between the CDCC calculation and the experimental elastic scattering data is limited and the breakdown does not exhaust the measured 𝜎 𝑅 . The resulting reduced reaction cross section 𝜎 Red for 12 N + 197 Au is large and comparable to the one obtained for the 8 B + 208 Pb system. This suggests a strong decoupling of the valence proton from the core because of the low binding energy and a dynamic polarization effect. Here, more research is required to estimate the contribution of core excitation.

Elastic scattering reactions↗

A High-resolution Regional Wave Resource Characterization For The U.S. West Coast

Objectives/Scope: Wave resource characterization is a critical step for wave energy converter deployment in the coastal ocean and relies on long-term, high-resolution wave datasets. This study presents a detailed modeling study of the wave resource along the U.S. West Coast (Washington, Oregon, and California), a coastal region that was identified with high wave energy potential in earlier studies. Methods, Procedures, Process: The wave hindcast covers a 32-year period from 1979 to 2010 and is based on a multi-resolution, unstructured-grid SWAN model framework. Model configuration closely follows and meets the requirements recommended by the International Electrotechnical Commission Technical Specification (IEC TS) for wave energy resource assessment and characterization (Class 2 - feasibility study). The model domain covers the entire U.S. Exclusive Economic Zone (EEZ) in the West Coast and has a spatial resolution varying from ~300 m in the nearshore region (20 km from the shoreline) to ~2500 m within the EEZ and ~5000 m at the open boundary, which extends beyond the EEZ. The model was forced by hourly 2-D wave spectra produced by a two-way nested WaveWatch III model, which covers the global ocean domain and the broader U.S. West Coast region domain with spatial resolutions of 0.5 degree and 10 arc-minutes, respectively. Both wave models are forced by hourly, 0.5-degree wind forcing obtained from NCEP’s Climate Forecast System Reanalysis (CFSR) product. Results, Observations, Conclusions: The standard model output for the SWAN model includes 3-hourly output for the six IEC wave resource parameters (e.g., omnidirectional wave power) at each grid point and hourly 2-D spectra at more than 50 NDBC buoys. Extensive model validation was achieved by comparing the six model-predicted IEC parameters with those derived from field observations at representative NDBC buoys. The error statistics indicated the model’s satisfactory performance. Further analyses were conducted to systematically evaluate the temporal and spatial distributions of wave energy potential and wave climate along the U.S. West Coast. Results suggest that Washington and Oregon coasts have similar nearshore wave resource, which is significantly higher than resources in Southern California. Strong seasonal variations are also observed, e.g., high wave energy tends to occur in the winter months. In summary, this study produced the first high-resolution, comprehensive dataset on wave energy distribution along the U.S. West Coast. Novel/Additive Information: The results are being used by the National Renewable Energy Laboratory to update the MHK Atlas, which was originally derived from NOAA’s 4-arc-minute WaveWatch III model output. In addition, the monthly averaged wave energy climatology dataset can be readily shared to support a variety of research and application efforts within the EEZ of the U.S. West Coast.

Wang, Taiping↗

TEAMER - AquaHarmonics High Fidelity WEC Sim PTO and Control Model Validation, Test Logs and Results

Collaborative effort between AquaHarmonics, Sandia National Laboratories (SNL), and the National Renewable Energy Laboratory (NREL) to revise and validate Aquaharmonics' full wave to wire model, allowing for reduced uncertainty and increased understanding of design requirements of a utility scale wave energy converter (WEC). SNL and NREL in collaboration with AquaHarmonics, will set up and run WEC Simulator (WEC-Sim) models of the AquaHarmonics WEC, building off past model developments for inclusion of custom PTO (power take-off) dynamics. The intent is to review, update, and verify or validate a new WEC-Sim model against wave tank experimental data. Furthermore, the WEC-Sim model will be coupled to an energy storage system model to better understand the wave-to-wire functionality. This data set is described in the "Test Log" excel file. Please refer to that document for details on each specific test date/time, constraint parameters and model hardware setup details. Sim model can be found in the associated MHKDR link below.

16 TIDAL AND WAVE POWER↗

Validation of (not-historical) large-event near-fault ground-motion simulations for use in civil engineering applications

Ground-motion simulations generated from physics-based wave propagation models are gaining increasing interest in the engineering community for their potential to inform the performance-based design and assessment of infrastructure residing in active seismic areas. A key prerequisite before the ground-motion simulations can be used with confidence for application in engineering domains is their comprehensive and rigorous investigation and validation. This article provides a four-step methodology and acceptance criteria to assess the reliability of simulated ground motions of not historical events, which includes (1) the selection of a population of real records consistent with the simulated scenarios, (2) the comparison of the distribution of Intensity Measures (IMs) from the simulated records, real records, and Ground-Motion Prediction Equations (GMPEs), (3) the comparison of the distribution of simple proxies for building response, and (4) the comparison of the distribution of Engineering Demand Parameters (EDPs) for a realistic model of a structure. Specific focus is laid on near-field ground motions (<10km) from large earthquakes (M w 7), for which the database of real records for potential use in engineering applications is severely limited. The methodology is demonstrated through comparison of (2490) near-field synthetic records with 5 Hz resolution generated from the Pitarka et al (2019) kinematic rupture model with a population of (38) pulse-like near-field real records from multiple events and, when applicable, with NGA-W2 GMPEs. Lastly, the proposed procedure provides an effective method for informing and advancing the science needed to generate realistic ground-motion simulations, and for building confidence in their use in engineering domains.

42 ENGINEERING↗

Controls on the frequency content of near-source infrasound at open-vent volcanoes: a case study from Volcán Villarrica, Chile

Infrasound signals from open-vent volcanoes can contain information related to that volcano’s eruption dynamics and future activity. Studying a specific volcano’s acoustic signature during periods of relative quiescence is thus important for understanding potential warning signs of impending eruptions. Volcán Villarrica, located in southern Chile, has a long-lived active lava lake that produces continuous infrasound with spectral peaks centered around 1 Hz that vary by + / - ~ 0.2 Hz over day-to-week-long time scales. The infrasound frequency content has been shown to illuminate key volcanic properties such as eruption style and crater shape. Leading up to Villarrica’s most recent paroxysm in 2015, for instance, infrasound spectral changes coincided with a rise in Villarrica’s lava lake level. Quantifying and understanding the regular fluctuation in recorded infrasonic frequencies from Villarrica and other open-vent volcanoes is thus imperative during periods of relative calm. Our experiment entailed a week-long period of crater rim infrasound observations associated with stable, open-vent activity and revealed two independent source processes: spatter bursting events (relatively broadband infrasound) and lava lake-induced tremor (sustained signals peaked at 1 Hz). A comparison of these acoustic signals with results from 3-D finite-difference time-domain wave propagation models shows that while sound speed variability can influence Villarrica’s spectrum, it cannot explain the full extent of the observed frequency excursions. We conclude that source spectrum variability is primarily responsible for the frequency excursions. In conclusion, our work highlights the utility of data collected from a dense distribution of twenty infrasound sensors operating at the summit coupled with numerical modeling of sound radiation.

58 GEOSCIENCES↗

Barrier-free predictions of short-range ordering/clustering kinetics in binary FCC solid solutions

We present comparisons of kinetic Monte Carlo (kMC) simulations of isothermal short-range ordering (SRO) and clustering (SRC) kinetics in binary FCC alloys with a mean-field concentration wave (CW) model. We find that the CW model is able to give order-of-magnitude agreement with kMC simulations for ordering/clustering relaxation times over a wide range of temperatures and compositions. The advantage of the CW model is that it does not require parameterization of vacancy hopping energy barriers, which, for a concentrated alloy, becomes prohibitive. We assess limits in the accuracy of the model and discuss the effect of cooling rates as well as the extension to multi-component systems. Ultimately, the simplicity and performance of the CW model compared to kMC simulations suggests that it is a useful tool to connect with models of properties dependent on SRO/SRC as well as for designing thermal treatments to control formation of SRO/SRC.

36 MATERIALS SCIENCE↗

Simulating nonlinear optical processes on a superconducting quantum device

Simulating plasma physics on quantum computers is difficult because most problems of interest are nonlinear, but quantum computers are not naturally suitable for nonlinear operations. In weakly nonlinear regimes, plasma problems can be modelled as wave–wave interactions. In this paper, we develop a quantization approach to convert nonlinear wave–wave interaction problems to Hamiltonian simulation problems. We demonstrate our approach using two qubits on a superconducting device. Unlike a photonic device, a superconducting device does not naturally have the desired interactions in its native Hamiltonian. Nevertheless, Hamiltonian simulations can still be performed by decomposing required unitary operations into native gates. To improve experimental results, we employ a range of error-mitigation techniques. Apart from readout error mitigation, we use randomized compilation to transform undiagnosed coherent errors into well-behaved stochastic Pauli channels. Moreover, to compensate for stochastic noise, we rescale exponentially decaying probability amplitudes using rates measured from cycle benchmarking. We carefully consider how different choices of product-formula algorithms affect the overall error and show how a trade-off can be made to best utilize limited quantum resources. This study provides an example of how plasma problems may be solved on near-term quantum computing platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha'apai

Understanding the forces and magma system dynamics on timescales of seconds to minutes remains challenging. In the January 2022 phreatoplinian Hunga Tonga-Hunga Ha'apai eruption, four remarkably similar seismic subevents within a 5-min interval occurred during the intensifying early eruptive phase. The subevents are similar in waveforms and durations (~25 s each). Each subevent begins with an unusual negative P-wave polarity which is inferred, using full-wave seismic modeling, to be caused by an upward single-force mechanism at the volcano created by a magma hammer likely in response to magma flow blockage/constriction during the early part of the eruption as discharge rapidly increased over orders of magnitude with concomitant conduit geometry evolution and instability. Our proposed episodic magma hammer model is consistent with thermodynamic and phase properties of the magmatic mixture, and yields an estimate of conduit mass flow in agreement with vent discharge rates derived from satellite imagery of plume heights.

58 GEOSCIENCES↗

Microseismicity Modulation Due To Changes in Geothermal Production at San Emidio, Nevada, USA

Brief cessations of geothermal production can induce seismicity, a phenomenon that has drawn increasing attention in recent years. Such observations are rare, and the underlying mechanism requires careful analysis. In April 2022, a dense seismic and hydrologic monitoring system was deployed at the San Emidio geothermal field, Nevada, to accompany a planned power plant shutdown. Using the dense seismic array data, we detected and located ∼1,800 microseismic events (MSEs) and developed a high-resolution tomographic P-wave velocity model. We observed substantially increased microseismicity during shutdown. Most MSEs occurred on pre-existing normal faults, which are contained within extremely low-velocity zones that are likely damaged, fluid-filled, and hydraulically connected to nearby production wells. Hydrologic data show rapid fluid pressure increases of <60 kPa following the shutdown. We suggest that the cessation of production rapidly increased fluid pressures along pre-existing fault zones, activating critically stressed fault patches and fractures and producing microseismicity.

Guo, Hao [University of Wisconsin-Madison, WI (Uni↗

Superconducting gap of H 3 S measured by tunnelling spectroscopy

Several hydrogen-rich superconductors have been found to show unprecedentedly high critical temperatures, stimulating investigations into the nature of the superconductivity in these materials. Although their macroscopic superconducting properties are established, microscopic insights into the pairing mechanism remains unclear. Here we characterize the superconducting gap structure in the high-temperature superconductor H 3 S and its deuterium counterpart D 3 S by performing tunnelling spectroscopy measurements. The tunnelling spectra reveal that H 3 S and D 3 S both have a fully gapped structure, which could be well described by a single s-wave Dynes model, with gap values 2Δ of approximately 60 meV and 44 meV, respectively. Furthermore, we observed gap features of another likely H-depleted H x S superconducting phase in a poorly synthesized hydrogen sulfide sample. Our work offers direct experimental evidence for superconductivity in the hydrogen-rich superconductor H 3 S from a microscopic perspective. It validates the phonon-mediated mechanism of superconducting pairing and provides a foundation for further understanding the origins of high-temperature superconductivity in hydrogen-rich compounds.

36 MATERIALS SCIENCE↗

Elucidating plasma dynamics in Hasegawa–Wakatani turbulence by information geometry

In this work, the impact of adiabatic electrons on drift-wave turbulence, modeled by the Hasegawa–Wakatani equations, is studied using information length. Information length is a novel theoretical method for measuring distances between statistical states represented by different probability distribution functions (PDFs) along the path of a system and represents the total number of statistically different states that a system evolves through in time. Specifically, the time-dependent PDFs of turbulent fluctuations for a given adiabatic index $A$ are computed. The changes in fluctuation statistics are then quantified in time by using information length. The numerical results provide time traces exhibiting intermittent plasma dynamics, and such behavior is identified by a rapid change in the information length. The effects of $A$ are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A super-resolution technique to analyze single-crystal inelastic neutron scattering measurements using direct-geometry chopper spectrometers

Direct-geometry time-of-flight chopper neutron spectroscopy is instrumental in studying dynamics in liquid, powder, and single crystal systems. We report here that real-space techniques in optical imagery can be adapted to obtain reciprocal-space super resolution dispersion for phonon or magnetic excitations from single-crystal neutron spectroscopy measurements. The procedure to reconstruct super-resolution energy dispersion of excitations relies on an accurate determination of the momentum and energy-dependent point spread function and a dispersion correction technique inspired by an image disparity calculation technique commonly used in stereo imaging. Here by, applying these methods to spinwave dispersion data from a virtual neutron experiment demonstrates ~5-fold improvement over nominal energy resolution.

47 OTHER INSTRUMENTATION↗

Quantum computing for fusion energy science applications

This is a review of recent research exploring and extending present-day quantum computing capabilities for fusion energy science applications. We begin with a brief tutorial on both ideal and open quantum dynamics, universal quantum computation, and quantum algorithms. Then, we explore the topic of using quantum computers to simulate both linear and nonlinear dynamics in greater detail. Because quantum computers can only efficiently perform linear operations on the quantum state, it is challenging to perform nonlinear operations that are generically required to describe the nonlinear differential equations of interest. In this work, we extend previous results on embedding nonlinear systems within linear systems by explicitly deriving the connection between the Koopman evolution operator, the Perron–Frobenius evolution operator, and the Koopman–von Neumann evolution (KvN) operator. We also explicitly derive the connection between the Koopman and Carleman approaches to embedding. Extension of the KvN framework to the complex-analytic setting relevant to Carleman embedding, and the proof that different choices of complex analytic reproducing kernel Hilbert spaces depend on the choice of Hilbert space metric are covered in the appendixes. Finally, we conclude with a review of recent quantum hardware implementations of algorithms on present-day quantum hardware platforms that may one day be accelerated through Hamiltonian simulation. We discuss the simulation of toy models of wave–particle interactions through the simulation of quantum maps and of wave–wave interactions important in nonlinear plasma dynamics.

Joseph, I. (ORCID:0000000255400840)↗