Modeling Hydrodynamic Instabilities, Shocks, and Radiation Waves in High Energy Density Experiments [Slides]
Abstract not provided.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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.
Practical approaches to modeling and control design for a wave energy converter.
Abstract not provided.
Coe, Ryan↗
Ground Motion Models (GMMs) Improvements Using Earthquake Simulations on High Performance Computers
A computationally efficient simulation platform was developed that can provide representative synthetic ground motions from crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using 3D modeling and high-performance computing. The main objective was to use synthetic ground motion to provide constrains to refinements of exiting ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-5Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. As part of a strategy for selecting a reginal velocity model and validation of developed rupture modeling technique, ground motions from the moment magnitude Mw5.0 November 7, 2016, Cushing Oklahoma, and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes were simulated. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model. Saikia’s model demonstrated better performance in modelling high frequency regional wave propagation for CEUS region. The proposed 3D model includes lateral variations added to the 1D background model using the stochastic scheme of Pitarka and Mellors. Comparisons of the simulations with recordings of both earthquakes demonstrated the reliability of our deterministic simulation approach while emphasizing the importance of including small-scale variability in the regional velocity model needed to reproduce the observed high-frequency wave scattering effects. As part of validation analysis, comparisons with different GMMs for a Mw6.5 earthquake in the CESUS region resulted in a very good match between the simulated and empirical ground motion models. Initial investigations of within-event and between-event ground motion variabilities for Mw6.5 scenario earthquakes on a strike-slip fault, suggest that they are strongly related to spatial slip and slip rate variations, average rupture velocity, rupture area and rupture initiation location. For certain scenarios we found that the ground motion variability observed at near-fault distances (< 5 km) also persists at longer distances. Regardless of the rupture scenario, the simulated ground motion tends to fully saturate at short distances and for all periods. The near-fault saturation has to do with the attenuation of waves propagating along the fault and local rupture radiation pattern that also contribute to stronger ground motion variation at such distances. Analysis of effects of rupture initiation location suggest that the peak ground motion (PGV) and spectral acceleration (SA) can be quite variable due to rupture directivity effects. Such effects are stronger at periods longer than 1s. The effect of the 1D velocity models and surface topography on simulated ground motion were investigated by comparing three component synthetic seismograms computed at selected sites. Effect of surface topography was considered using the ratio between spectral accelerations simulated for two 1D models with flat surface topography and realistic model with surface topography. Overall, the topography slightly amplifies (by ~30%) the ground motion amplitude in the frequency range 1-3Hz. The effect of topography is more visible in the surface and coda waves portion of the seismograms.
58 GEOSCIENCES↗
Trends of Sediment Resuspension and Budget in Southern Lake Michigan Under Changing Wave Climate and Hydrodynamic Environment
Sediment suspension and transport driven by waves and currents play a significant role in both the ecological and physical environments of large lakes. Lake Michigan has faced a rapidly increasing water level associated with intensified wind waves in the past decade. To investigate the spatiotemporal characteristics of suspended sediment concentration (SSC) and associated coastal sediment budgets in southern Lake Michigan, a 30-year (1991–2020) hindcast was performed using a coupled wave-current-sediment model (SWAN-FVCOM-CSTMS). We found that in southern Lake Michigan, the basin-wide mean SSC increased, and the coastal sediment loss accelerated dramatically, corresponding with intensified waves, currents and lake water level rises over the past decade. The basin-wide mean SSC, coastal sediment loss, wave height, wind speed, current speed, and water level in southern Lake Michigan are highly correlated. Spatially, the results reveal decreases in coastal SSC and sediment loss in the western portion of the southern basin, while the eastern sectors show an increase in both metrics. This reflects a clear shift in the wave climate and hydrodynamic environment. The alterations in long-term coastal sediment budgets imply that considerable shoreline transformations are being influenced by modifications in the wave climate. Understanding the spatiotemporal characteristics of SSC and coastal sediment budgets is crucial for strategic water resource management and coastal infrastructure planning.
54 ENVIRONMENTAL SCIENCES↗
Data Files for P. Mai et al., "Fluctuating charge-density-wave correlations in the three-band Hubbard model" (2024)
The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model’s fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.
cuprates↗
Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code
A passive yaw implementation is developed, validated, and explored for the WEC-Sim, an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01° was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than 1°, for which computation time is only increased by 2.2x.
50 EE - Wind and Water Power Program - Water (EE-4↗
Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code: Preprint
A passive yaw implementation is developed, validated, and explored for the Wave Energy Converter Simulator (WEC-Sim), an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01 was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than or equal to 1, for which computation time is only increased by 2.2x.
50 EE - Wind and Water Power Program - Water (EE-4↗
Anomalous compressible mode generation by global frame projections of pure Alfven mode
Alfven wave is the single most important physical phenomenon of magneto-hydrodynamic turbulence and has far-reaching impact to almost all studies related to astrophysical magnetic field. Yet the restoration of the Alfven wave fluctuations from a given magnetic field, aka the local Alfven wave problem, is never properly addressed in literature albeit its importance. Previous works model the Alfven wave fluctuation as the perturbation along a straight-line, constant magnetic field. However, Lazarian & Pogosyan (2012) suggested that the decomposition of Alfven wave along a straight line, aka. the global frame decomposition, has a factor of discrepancy to the true local Alfven wave fluctuation. Here, in this paper, we provide a geometric interpretation on how the local Alfven wave is related to the global frame through the use of vector frame formulation. We prove both analytically and numerically that the local frame Alfven wave is an orthogonal transformation of that of the global frame and related by the local Alfvenic Mach number. In other words, when we observe Alfven wave in the global frame of reference, some of the Alfven wave will be mistaken as compressible waves. The importance of frame choices have a far-reaching impact to the analytical studies of MHD turbulence. Combining the frame formalism and the new techniques we can have accurate measurement to some of the fundamental turbulence properties like the inclination angle of mean magnetic field relative to the line of sight.
79 ASTRONOMY AND ASTROPHYSICS↗
Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F
The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.
6 ≤ A ≤ 19↗
Computational modelling and experimental tank testing of the multi float WaveSub under regular wave forcing
A submerged wave device generates energy from the relative motion of floating bodies. In WaveSub, three floats are joined to a reactor; each connected to a spring and generator. Electricity generated damps the orbital movements of the floats. The forces are non-linear and each float interacts with the others. Tuning to the wave climate is achieved by changing the line lengths, so there is a need to understand the performance trade-offs for a large number of configurations. This requires an efficient, large displacement, multidirectional, multi-body numerical scheme. Results from a 1/25 scale wave basin experiment are described. Here, we show that a time domain linear potential flow formulation (Nemoh, WEC-Sim) can match the tank testing provided that suitably tuned drag coefficients are employed. Inviscid linear potential models can match some wave device experiments; however, additional viscous terms generally provide better accuracy. Scale experiments are also prone to mechanical friction, and we estimate friction terms to improve the correlation further. The resulting error in mean power between numerical and physical models is approximately 10%. Predicted device movement shows a good match. Overall, drag terms in time domain wave energy modelling will improve simulation accuracy in wave renewable energy device design. to understand the performance trade-offs for a large number of configurations. This requires an efficient, large displacement, multidirectional, multi-body numerical scheme. Results from a 1/25 scale wave basin experiment are described. Here, we show that a time domain linear potential flow formulation (Nemoh, WEC-Sim) can match the tank testing provided that suitably tuned drag coefficients are employed. Inviscid linear potential models can match some wave device experiments; however, additional viscous terms generally provide better accuracy. Scale experiments are also prone to mechanical friction, and we estimate friction terms to improve the correlation further. The resulting error in mean power between numerical and physical models is approximately 10%. Predicted device movement shows a good match. Overall, drag terms in time domain wave energy modelling will improve simulation accuracy in wave renewable energy device design.
16 TIDAL AND WAVE POWER↗
Charge density wave and superconductivity in the disordered Holstein model
The interplay between electron-electron correlations and disorder has been a central theme of condensed matter physics over the past several decades, with particular interest in the possibility that interactions might cause delocalization of an Anderson insulator into a metallic state, and the disrupting effects of randomness on magnetic order and the Mott phase. In this work, we extend this physics to explore electron-phonon interactions and show, via exact quantum Monte Carlo simulations, that the suppression of the charge density wave correlations in the half-filled Holstein model by disorder can stabilize a superconducting phase. Our simulations thus capture qualitatively the suppression of charge ordered phases and emergent superconductivity recently seen experimentally.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Verifying Ray Tracing Amplitude Methods for Global Magnetospheric Modeling
Abstract Ray tracing is a commonly used method for modeling the propagation of electromagnetic waves in Earth's magnetosphere. To apply ray tracing results to global models of wave‐particle interaction such as energetic electron scattering, it is useful to map the discrete rays to a volume filling mesh. However, some methods have inherent losses of energy from the wave source, or do not account for the full range of wave properties within a sample volume. We have developed and tested a 3D magnetospheric ray tracing code “MESHRAY” which resolves these issues. MESHRAY uses the conservation of Poynting flux through ray triplets with finite volume to determine the local field amplitudes. Electromagnetic wave energy density from all ray data points is mapped to a mesh and verified against the wave source power for energy conservation varying time step length, number of rays, and total time steps. We find that the method is self‐consistent and numerically robust. We further investigate whether the neglect of phase information and superposition has a significant impact on the accuracy of mapping wave intensity to a mesh. We find excellent agreement between the analytic solution for waves emitted by a line source in a plane‐stratified medium and an equivalent ray tracing solution. When phase information is excluded, ray tracing reproduces an average amplitude spread over regions of coherent constructive and destructive interference. This may be an important consideration for interpolating ray tracing results of longer wavelength waves such as magnetosonic, electromagnetic ion cyclotron, or ULF waves.
79 ASTRONOMY AND ASTROPHYSICS↗
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↗
High-resolution lithospheric structure of continental China from joint inversion of surface wave and gravity data
The lithospheric structure of continental China has been previously determined by seismic travel time tomography, surface wave tomography, and joint inversion of body wave and surface wave data. However, due to the inherent limitations of seismic data, the lithospheric structure of continental China is still not well resolved in the shallow part and in some regions where the station coverage is relatively sparse. In this study, we aim at improving the lithospheric structure by joint inversion of seismic surface wave data and satellite gravity data to take advantage of the uniform distribution and complementary strength of the gravity data. Here, the empirical relationship between velocity and density is used as a bridge for joint inversion of surface wave and gravity data. The joint inversion shear-wave velocity Vs (density) model, named as USTClitho1.0g, can fit both surface wave and gravity data well. This high-resolution Vs model can better fit the active airgun source seismic arrival times and better delineate some features in continental China, such as the velocity contrast across the north–south gravity lineament (NSGL), the lithosphere thinning in eastern China, the crustal footprint of the Hainan mantle plume, the likely magma chamber beneath volcanos in northeast China, the middle-lower crust low velocity layer beneath the Tibetan plateau, and the tearing of subducted Indian mantle lithosphere. Our joint inversion Vs model can provide a reference model for geosciences in continental China and surrounding areas.
58 GEOSCIENCES↗
Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses
SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.
Liu, Xuejian↗
The Collaborative Seismic Earth Model: Generation 2
Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.
58 GEOSCIENCES↗
ENSO-Induced Teleconnection: Process-Oriented Diagnostics to Assess Rossby Wave Sources and Ambient Flow Properties in Climate Models
Climate model fidelity in representing ENSO-induced teleconnection is assessed with process-oriented diagnostics that examine a chain of processes, from equatorial Pacific precipitation to the midlatitude circulation pattern over the Pacific–North American regions. Such processes are rarely addressed during model development. Using an upper-tropospheric divergent level, local vorticity gradient of the ambient zonal flow (∂ 2 U¯/∂y 2 ) and a restoring force for Rossby waves (β * ) are estimated, the equivalent barotropic vorticity equation is solved, and an anomalous Rossby wave source (RWS') quantified. The analysis is applied to AMIP5 and AMIP6 simulations. For a realistic circulation response representation, the hypothesis that models accurately represent the strength and location of RWS', and spatial variations in β * is tested. Compared to AMIP5, in AMIP6 there are clear improvements in representing RWS' and β * . To validate the hypothesis, the analysis identifies two metrics: spatially coherent RWS' in the subtropical North Pacific, and longitudes of negative β * over the western-central North Pacific. By projecting these metrics in two and three-dimensional views, improvements or degradations in model versions are apparent. If a model’s fidelity in representing ∂ 2 U¯/∂y 2 and RWS' are compromised, then radiated Rossby waves are reflected more equatorward, resulting in zonally elongated circulation anomalies over the central North Pacific. Thus, during climate model development, applying this analysis frequently will keep a regular check on the fidelity of the modeled response to anomalous El Niño convection in conjunction with changing model ambient flow dependencies. As a result, this analysis is intended to form a process-oriented diagnostics package, a community contribution to the NOAA Model Diagnostics Task Force.
54 ENVIRONMENTAL SCIENCES↗
Quantum simulations of hadron dynamics in the Schwinger model using 112 qubits
Hadron wave packets are prepared and time evolved in the Schwinger model using 112 qubits of IBM’s 133-qubit Heron quantum computer ibm_torino. The initialization of the hadron wave packet is performed in two steps. First, the vacuum is prepared across the whole lattice using the recently developed SC-ADAPT-VQE algorithm and workflow. SC-ADAPT-VQE is then extended to the preparation of localized states, and used to establish a hadron wave packet on top of the vacuum. This is done by adaptively constructing low-depth circuits that maximize the overlap with an adiabatically prepared hadron wave packet. Due to the localized nature of the wavepacket, these circuits can be determined on a sequence of small lattices using classical computers, and then robustly scaled to prepare wave packets on large lattices for simulations using quantum computers. Time evolution is implemented with a second-order Trotterization. To reduce both the required qubit connectivity and circuit depth, an approximate quasilocal interaction is introduced. This approximation is made possible by the emergence of confinement at long distances, and converges exponentially with increasing distance of the interactions. Using multiple error-mitigation strategies, up to 14 Trotter steps of time evolution are performed, employing 13,858 two-qubit gates (with a CNOT depth of 370). The propagation of hadrons is clearly identified, with results that compare favorably with Matrix Product State simulations. Finally, prospects for a near-term quantum advantage in simulations of hadron scattering are discussed.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗