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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Complex Density Wave Orders and Quantum Phase Transitions in a Model of Square-Lattice Rydberg Atom Arrays

In this work, we describe the zero-temperature phase diagram of a model of a two-dimensional square-lattice array of neutral atoms, excited into Rydberg states and interacting via strong van der Waals interactions. Using the density-matrix renormalization group algorithm, we map out the phase diagram and obtain a rich variety of phases featuring complex density wave orderings, upon varying lattice spacing and laser detuning. While some of these phases result from the classical optimization of the van der Waals energy, we also find intrinsically quantum-ordered phases stabilized by quantum fluctuations. These phases are surrounded by novel quantum phase transitions, which we analyze by finite-size scaling numerics and Landau theories. Our work highlights Rydberg quantum simulators in higher dimensions as promising platforms to realize exotic many-body phenomena.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Probing the evolution of fault properties during the seismic cycle with deep learning

We use seismic waves that pass through the hypocentral region of the 2016 M6.5 Norcia earthquake together with Deep Learning (DL) to distinguish between foreshocks, aftershocks and time-to-failure (TTF). Binary and N-class models defined by TTF correctly identify seismograms in test with > 90% accuracy. We use raw seismic records as input to a 7 layer CNN model to perform the classification. Here we show that DL models successfully distinguish seismic waves pre/post mainshock in accord with lab and theoretical expectations of progressive changes in crack density prior to abrupt change at failure and gradual postseismic recovery. Performance is lower for band-pass filtered seismograms (below 10 Hz) suggesting that DL models learn from the evolution of subtle changes in elastic wave attenuation. Tests to verify that our results indeed provide a proxy for fault properties included DL models trained with the wrong mainshock time and those using seismic waves far from the Norcia mainshock; both show degraded performance. Our results demonstrate that DL models have the potential to track the evolution of fault zone properties during the seismic cycle. If this result is generalizable it could improve earthquake early warning and seismic hazard analysis.

58 GEOSCIENCES↗
DOE OSTI2021

Real-Time Nonlinear Model Predictive Controller for Multiple Degrees of Freedom Wave Energy Converters with Non-Ideal Power Take-Off

An increase in wave energy converter (WEC) efficiency requires not only consideration of the nonlinear effects in the WEC dynamics and the power take-off (PTO) mechanisms, but also more integrated treatment of the whole system, i.e., the buoy dynamics, the PTO system, and the control strategy. It results in an optimization formulation that has a nonquadratic and nonstandard cost functional. This article presents the application of real-time nonlinear model predictive controller (NMPC) to two degrees of freedom point absorber type WEC with highly nonlinear PTO characteristics. The nonlinear effects, such as the fluid viscous drag, are also included in the plant dynamics. The controller is implemented on a real-time target machine, and the WEC device is emulated in real-time using the WECSIM toolbox. The results for the successful performance of the design are presented for irregular waves under linear and nonlinear hydrodynamic conditions.

Haider, Ali S. (ORCID:0000000246513157)↗
DOE OSTI2023

A Moving-Wave Implementation in WRF to Study the Impact of Surface Water Waves on the Atmospheric Boundary Layer

Abstract While numerous modeling studies have focused on the interaction of ocean surface waves with the atmospheric boundary layer, most employ idealized waves that are either monochromatic or synthetically generated from a theoretical wave spectrum, and the atmospheric solvers are typically incompressible. To study wind–wave coupling in real-world scenarios, a model that can simulate both realistic meteorological and wave conditions is necessary. In this paper we describe the implementation of a moving bottom boundary condition into the Weather Research and Forecasting Model for large-eddy simulation applications. We first describe the moving bottom boundary conditions within WRF’s pressure-based vertical coordinate system. We then validate our code with idealized test cases that have analytical solutions, including flow over a monochromatic wave with and without viscosity. Finally, we present results from turbulent flows over a moving monochromatic wave with different wave ages, and demonstrate satisfactory agreement of the wave growth rate with results from the literature. We also compare atmospheric stress and wind parameters from two physically equivalent cases. The first specifies a wind moving in the same direction as a propagating wave, while the second involves a stationary wave with the wind adjusted such that the wind relative to the wave is the same as in the first case. Results indicate that the velocity and Reynolds stress profiles for the two cases match, further validating the moving bottom implementation.

17 WIND ENERGY↗
DOE OSTI2021

Imaging the Shallow Structure of the Yucca Flat at the Source Physics Experiment Phase II Site with Horizontal-to-Vertical Spectral Ratio Inversion and a Large- N Seismic Array

The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.

58 GEOSCIENCES↗
DOE OSTI2025

Simulating nationwide coupled disease and fear spread in an agent-based model

Human cognitive responses, behavioral responses, and disease dynamics co-evolve over the course of any disease outbreak, and can result in complex feedbacks. We present a dynamic agent-based model that explicitly couples the spread of disease with the spread of fear surrounding the disease, implemented within the EpiCast simulation framework. EpiCast models transmission within a realistic synthetic population, capturing individual-level interactions. In our model, fear propagates through both in-person contact and broadcast media, prompting individuals to adopt protective behaviors that reduce disease spread. In order to better understand these coupled dynamics, we create and compare a range of compartmental models to ensure that introducing additional disease states does not prevent the emergence of multiple waves in these simpler models. Additionally, we compare a range of behavioral scenarios within EpiCast, varying the level and intensity of fear and behavior change. Our results show that the addition of asymptomatic, exposed, and pre-symptomatic disease states can impact both the rate at which an outbreak progresses and its overall trajectory in compartmental models. In EpiCast, the combination of non-local fear spread via broadcasters and strong behavioral responses by fearful individuals generally leads to multiple epidemic waves, an outcome that occurs only within a narrow parameter range when fear spreads purely through local contact. Accounting for the coupled spread of fear and disease is critical for understanding disease dynamics and designing timely, targeted responses to emerging infectious threats.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗
DOE OSTI2021

Gravitational waves from first-order phase transitions in Majoron models of neutrino mass

We show how the generation of right-handed neutrino masses in Majoron models may be associated with a first-order phase transition and accompanied by the production of a stochastic background of gravitational waves (GWs). We explore different energy scales with only renormalizable operators in the effective potential. If the phase transition occurs above the electroweak scale, the signal can be tested by future interferometers. We consider two possible energy scales for phase transitions below the electroweak scale. If the phase transition occurs at a GeV, the signal can be tested at LISA and provide a complementary cosmological probe to right-handed neutrino searches at the FASER detector. If the phase transition occurs below 100 keV, we find that the peak of the GW spectrum is two or more orders of magnitude below the putative NANOGrav GW signal at low frequencies, but well within reach of the SKA and THEIA experiments. We show how searches of very low frequency GWs are motivated by solutions to the Hubble tension in which ordinary neutrinos interact with the dark sector. We also present general calculations of the phase transition temperature and Euclidean action that apply beyond Majoron models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

CO 2 rock physics modeling for reliable monitoring of geologic carbon storage

Monitoring, verification, and accounting (MVA) are crucial to ensure safe and long-term geologic carbon storage. Seismic monitoring is a key MVA technique that utilizes seismic data to infer elastic properties of CO 2 -saturated rocks. Reliable accounting of CO 2 in subsurface storage reservoirs and potential leakage zones requires an accurate rock physics model. However, the widely used CO 2 rock physics model based on the conventional Biot-Gassmann equation can substantially underestimate the influence of CO 2 saturation on seismic waves, leading to inaccurate accounting. We develop an accurate CO 2 rock physics model by accounting for both effects of the stress dependence of seismic velocities in porous rocks and CO 2 weakening on the rock framework. We validate our CO 2 rock physics model using the Kimberlina-1.2 model (a previously proposed geologic carbon storage site in California) and create time-lapse elastic property models with our new rock physics method. We compare the results with those obtained using the conventional Biot-Gassmann equation. Our innovative approach produces larger changes in elastic properties than the Biot-Gassmann results. Using our CO 2 rock physics model can replicate shear-wave speed reductions observed in the laboratory. Our rock physics model enhances the accuracy of time-lapse elastic-wave modeling and enables reliable CO 2 accounting using seismic monitoring.

58 GEOSCIENCES↗
DOE OSTI2022

Verifying raytracing/Fokker–Planck lower-hybrid current drive predictions with self-consistent full-wave/Fokker–Planck simulations

Raytracing/Fokker–Planck (FP) simulations used to model lower-hybrid current drive (LHCD) often fail to reproduce experimental results, particularly when LHCD is weakly damped. A proposed reason for this discrepancy is the lack of ‘full-wave’ effects, such as diffraction and interference, in raytracing simulations and the breakdown of the raytracing approximation. Previous studies of LHCD using non-Maxwellian full-wave/FP simulations have been performed, but these simulations were not self-consistent and enforced power conservation between the FP and full-wave code using a numerical rescaling factor. Here, we have created a fully self-consistent full-wave/FP model for LHCD that is automatically power conserving. This was accomplished by coupling an overhauled version of the non-Maxwellian TORLH full-wave solver and the CQL3D FP code using the Integrated Plasma Simulator. We performed converged full-wave/FP simulations of Alcator C-Mod discharges and compared them with raytracing. We found that excellent agreement in the power deposition profiles from raytracing and TORLH could be obtained, however, TORLH had somewhat lower current drive efficiency and broader power deposition profiles in some cases. This discrepancy appears to be a result of numerical limitations present in the TORLH model and a small amount of diffractional broadening of the TORLH wave spectrum. Our results suggest full-wave simulation of LHCD is likely not necessary as diffraction and interference represented only a small correction that could not account for the differences between simulations and experiment.

Physics↗
DOE OSTI2022

Crust and Upper Mantle Structure Beneath the Eastern United States

Abstract The Eastern United States (EUS) has a complex geological history and hosts several seismic active regions. We investigate the subsurface structure beneath the broader EUS. To produce reliable images of the subsurface, we simultaneously invert smoothed P‐wave receiver functions, Rayleigh‐wave phase and group velocity measurements, and Bouguer gravity observations for the 3D shear‐wave speed. Using surface‐wave observations (3–250 s) and spatially smoothed receiver functions, our velocity models are robust, reliable, and rich in detail. The shear‐wave velocity models fit all three types of observations well. The resulting velocity model for the eastern U.S. shows thinner crust beneath New England, the east coast, and the Mississippi Embayment (ME). A relatively thicker crust was found beneath the stable North America craton. A relatively slower upper mantle was imaged beneath New England, the east coast, and western ME. A comparison of crust thickness derived from our model against four recent published models shows first‐order consistency. A relatively small upper mantle low‐speed region correlates with a published P‐wave analysis that has associated the anomaly with a 75 Ma kimberlite volcanic site in Kentucky. We also explored the relationship between the subsurface structure and seismicity in the eastern U.S. We found that earthquakes often locate near regions with seismic velocity variations, but not universally. Not all regions of significant subsurface wave speed changes are loci of seismicity. A weak correlation between upper mantle shear velocity and earthquake focal mechanism has been observed.

58 GEOSCIENCES↗
DOE OSTI2023

Multi-scenario Extreme Weather Simulator application to heat waves: Ko’olauloa community resilience hub

Heat waves are increasing in severity, duration, and frequency. The Multi-Scenario Extreme Weather Simulator (MEWS) models this using historical data, climate model outputs, and heat wave multipliers. In this study, MEWS is applied for planning of a community resilience hub in Hau'ula, Hawaii. The hub will have normal operations and resilience operations modes. Both these modes were modeled using EnergyPlus. The resilience operations mode includes cutting off air conditioning for many spaces to decrease power requirements during emergencies. Results were simulated for 300 future weather files generated by MEWS for 2020, 2040, 2060, and 2080. Shared socioeconomic pathways 2-4.5, 3-7.0 and 5-8.5 were used. The resilience operations mode results show two to six times increase of hours of exceedance beyond 32.2 °C from present conditions, depending on climate scenario and future year. The resulting decrease in thermal resilience enables an average decrease of energy use intensity of 26% with little sensitivity to climate change. The decreased thermal resilience predicted in the future is undesirable, but was not severe enough to require a more energy-intensive resilience mode. Instead, planning is needed to assure vulnerable individuals are given prioritized access to air-conditioned parts of the hub if worst-case heat waves occur.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model

A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Crustal and uppermost mantle structure of Cape Verde from ambient noise tomography

SUMMARY We present a seismic ambient noise tomography of the Cape Verde archipelago, located in the Atlantic Ocean, approximately 600 km west of Senegal. We used 38 seismic broad-band stations that continuously recorded for 10 months, in order to construct the first 3-D model of Sv-wave velocities for the crust and uppermost mantle beneath the Cape Verde region. We started by computing phase cross-correlations for vertical component recordings using all possible inter-island station pairs. Next, a time–frequency phase-weighted stack was applied to obtain robust Rayleigh-wave group-velocity dispersion curves in the period band between 10 and 24 s. Group-velocity maps at different periods are obtained by inverting the dispersion curves. We then inverted the group-velocity maps to obtain the 3-D shear wave velocity structure of the crust and uppermost mantle beneath Cape Verde. The final 3-D model extends from 8 km down to 23 km and has a lateral resolution of about 50 km. The crust in the southwestern sector, encompassing Fogo, presents lower S-wave velocities that may be caused by the presence of melt pockets and/or hydrothermal fluids circulation. The uppermost mantle beneath the northwestern sector is characterized by higher S-wave velocities in agreement with previous results obtained from Ps and Sp receiver functions. Those high-velocity anomalies can reflect non-altered crust or remnants of magma chambers or solidified basaltic intrusions, which fed the volcanism in these islands. Our maps revealed the presence of crustal underplating across the entire archipelago, yet stronger beneath the groups Santo Antão—São Vicente—São Nicolau and Fogo—Santiago—Maio.

Carvalho, J. (ORCID:0000000310778994)↗
DOE OSTI2024

Modeling and Design of the Triton OWC System (Final Report)

Two versions of the Triton Oscillating Water Column type device will be modeled in WEC Sim. First, a model of a wave-tank scale device which can be tuned and validated against tank test data. Second, a model of the deployment-scale device will be constructed following the method of the tank-scale device to ensure that relevant physics are captured. In the latter case, the geometry, PTO architecture, and other design details are not yet finalized, so the model will serve as a platform for design iteration as time and budget allows. A subset of this iteration will be automated using existing WEC-Sim capabilities. A primary focus of this work will be familiarizing Triton personnel with the WEC-Sim workflow and model details so that the model of the deployment device can continue to be enhanced after project end.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Evaluating the effect of meso/submesoscale current–wave interactions on wave energy resource characterization at northeast U.S. coast

Wave energy is a promising renewable resource, but accurate assessment is difficult in regions with strong currents due to wave–current interactions (WCI). Here, this study develops a two-way coupled WCI model within the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) framework at 2 km resolution to improve wave energy characterization along the northeastern U.S. coast, including the Mid-Atlantic Bight and Gulf of Maine. The model integrates WaveWatchIII (WWIII) and the Regional Ocean Modeling System (ROMS) to enhance wave hindcasting by accounting for Doppler-shift, refraction, and nonlinear energy exchanges. Validation against buoy and satellite observations confirms model accuracy. Analysis shows that Doppler-shifting can alter wave power density by over 20%, while strong current gradients and shear distort wave crests via focusing/defocusing and stretching/squeezing, modifying wave direction and frequency. These processes together can induce wave power fluctuations of up to 40% on synoptic scales. Applying a 2.5 MW Ocean Energy Converter power matrix shows that WCI may change harvested energy by up to 100% in shallow-waters and 60% in deep-waters. These results underscore the importance of incorporating nonlinear WCI for reliable wave climate predictions and resource assessments in energetic coastal regions.

doppler-shift↗
DOE OSTI2025

X-ray fluorescence standing wave study of the interaction of the antimicrobial peptide indolicidin with a supported model membrane

Indolicidin, a cationic antimicrobial peptide, interacts with lipid bilayers through electrostatic and hydrophobic interactions, disrupting microbial membranes. We investigated the depth-dependent localization of gold-nanoparticle-labeled indolicidin in a supported model membrane using X-ray fluorescence standing wave (SWXF) analysis. Liposomes composed of DMPC and DMPG were incubated with indolicidin labeled at its C-terminus with a 1.8 nm gold nanoparticle, then deposited onto a Si/Mo multilayer substrate via vesicle bursting. SWXF measurements revealed that at low peptide incubation concentrations (2–5 μM), gold-nanoparticle-labeled indolicidin remains primarily associated with the bilayer's outer leaflet. At higher concentrations (10 μM), the peptide penetrates deeper into the bilayer, with the labeled C-terminal region either localizing near the membrane's hydrophobic core or inducing membrane breakup. In conclusion, these findings suggest a concentration-dependent insertion mechanism.

Basnet, Gobind [Northern Illinois Univ., DeKalb, I↗
DOE OSTI2021

Observations of Gravity Waves in the Middle Atmosphere of Mars

Gravity waves are ubiquitous throughout the atmosphere of Mars. Their propagation and dissipation influence the circulation and thermal structure of the middle and upper atmosphere. Yet there have been few studies of gravity wave characteristics in the middle atmosphere, a region that is critical for their propagation from generation in the lower atmosphere to dissipation and associated exchanges of momentum and energy in the upper atmosphere. There have been few studies because few atmospheric profiles span the middle atmosphere with the kilometer-scale or smaller vertical resolution that is Required to characterize gravity waves. Here we report the characterization of gravity waves in the middle atmosphere using 12 high-resolution atmospheric density profiles. Four of these were acquired from a ground-based stellar occultation from 1976 that yielded immersion and emersion profiles on opposite sides of the planet, and the remaining eight were measured during atmospheric entry by landers and rovers. Predominant wavelengths were 3–14 km, and amplitudes were generally 0.8%–2.5%. Where static stability is large and positive, gravity waves grow efficiently. In other instances, static stability is not large and positive over a wide altitude range, and gravity wave amplitudes do not behave as neatly. These observations of gravity waves in the middle atmosphere of Mars can be used to test gravity wave parameterizations in large-scale general circulation models and to investigate predictions for how gravity wave propagation and dissipation influence the circulation and thermal structure of the middle and upper atmosphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗