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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Gravitational waves from current-carrying cosmic strings

Cosmic strings are predicted by many Standard Model extensions involving the cosmological breaking of a symmetry with nontrivial first homotopy group and represent a potential source of primordial gravitational waves (GWs). Present efforts to model the GW signal from cosmic strings are often based on minimal models, such as, e.g., the Nambu-Goto action that describes cosmic strings as exactly one-dimensional objects without any internal structure. In order to arrive at more realistic predictions, it is therefore necessary to consider nonminimal models that make an attempt at accounting for the microscopic properties of cosmic strings. With this goal in mind, we derive in this paper the GW spectrum emitted by current-carrying cosmic strings (CCCSs), which may form in a variety of cosmological scenarios. Our analysis is based on a generalized version of the velocity-dependent one-scale (VOS) model, which, in addition to the mean velocity and correlation length of the string network, also describes the evolution of a chiral (light-like) current. Here we are able to show, the solutions of the VOS equations imply a temporarily growing fractional cosmic-string energy density, Ω cs . This results in an enhanced GW signal across a broad frequency interval, whose boundaries are determined by the times of generation and decay of cosmic-string currents. Our findings have important implications for GW experiments in the Hz to MHz band and motivate the construction of realistic particle physics models that give rise to large currents on cosmic strings.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Laboratory Upgrade Point Absorber WEC-Sim Model with MoorDyn Moorings

This dataset includes a WEC-Sim and MoorDyn model of the Laboratory Upgrade Point Absorber (LUPA). LUPA is an open-source wave energy converter designed and tested by Oregon State University. The files provided here constitute a stable LUPA configuration with three mooring lines. This model is 1/20 scale, optimized for the O.H. Hinsdale Wave Lab at Oregon State University. This model of LUPA adds MoorDyn functionality for more accurate mooring predictions and uses a more stable, updated version of LUPA's current physical configuration. This model is for WEC-Sim Version 6.0. A recent update of WEC-Sim has changed some functionality of MoorDyn such that this model will not work with WEC-Sim Version 6.1.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Effects of polarization-reversed EMIC waves on the ring current dynamics

Electromagnetic ion cyclotron (EMIC) waves are widely believed to play an important role in influencing the radiation belt and ring current dynamics. Most studies have investigated the effects or characteristics of EMIC waves by assuming their left-handed polarization. However, recent studies have found that the reversal of polarization, which occurs at higher latitudes along the wave propagation path, can change the wave-induced pitch angle diffusion coefficients. Whether such a polarization reversal can influence the global ring current dynamics remains unknown. In this study, we investigate the ring current dynamics and proton precipitation loss in association with polarization-reversed EMIC waves by using the ring current–atmosphere interactions model (RAM). The results indicate that the polarization reversal of H-band EMIC waves can truly decrease the scattering rates of protons of 10 to 50 keV or >100 keV in comparison with the scenario in which the EMIC waves are considered purely left-handed polarized. Additionally, the global ring current intensity and proton precipitation may be slightly affected by the polarization reversal, especially during prestorm time and the recovery phase, but the effects are not large during the main phase. This is probably because the H-band EMIC waves contribute to the proton scattering loss primarily at E < 10 keV, an energy range that is not strongly affected by the polarization reversal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Short-depth QAOA circuits and quantum annealing on higher-order ising models

Abstract We present a direct comparison between QAOA (Quantum Alternating Operator Ansatz), and QA (Quantum Annealing) on 127 qubit problem instances. QAOA with p = 1, 2 rounds is executed on the 127 qubit heavy-hex graph gate-model quantum computer ibm_washington, using on-device grid-searches for angle finding, and QA is executed on two Pegasus-chip D-Wave quantum annealers. The problems are random Ising models whose connectivity matches heavy-hex graphs and the Pegasus graph connectivity, and optionally include hardware-compatible cubic terms ( Z Z Z terms). The QAOA circuits are heavily optimized and of extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of the heavy-hex lattice. QAOA and QA are both compared against simulated annealing and the optimal solutions are computed exactly using CPLEX. The noiseless mean QAOA expectation values for p = 1, 2 are computed using classical light-cone based simulations. We find QA outperforms QAOA on the evaluated devices.

127 qubits↗
DOE OSTI2023

Direct imaging of shock wave splitting in diamond at Mbar pressure

Understanding the behavior of matter at extreme pressures of the order of a megabar (Mbar) is essential to gain insight into various physical phenomena at macroscales—the formation of planets, young stars, and the cores of super-Earths, and at microscales—damage to ceramic materials and high-pressure plastic transformation and phase transitions in solids. Under dynamic compression of solids up to Mbar pressures, even a solid with high strength exhibits plastic properties, causing the induced shock wave to split in two: an elastic precursor and a plastic shock wave. This phenomenon is described by theoretical models based on indirect measurements of material response. The advent of x-ray free-electron lasers (XFELs) has made it possible to use their ultrashort pulses for direct observations of the propagation of shock waves in solid materials by the method of phase-contrast radiography. However, there is still a lack of comprehensive data for verification of theoretical models of different solids. Here, we present the results of an experiment in which the evolution of the coupled elastic–plastic wave structure in diamond was directly observed and studied with submicrometer spatial resolution, using the unique capabilities of the x-ray free-electron laser (XFEL). The direct measurements allowed, for the first time, the fitting and validation of the 2D failure model for diamond in the range of several Mbar. Our experimental approach opens new possibilities for the direct verification and construction of equations of state of matter in the ultra-high-stress range, which are relevant to solving a variety of problems in high-energy-density physics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Frequency-domain modeling of floating wind arrays with shared mooring lines

A frequency-domain model for floating wind turbine dynamics has been extended to model floating wind farms with couplings from shared mooring systems. The model, called RAFT, could previously calculate the mean offsets and wave-induced response spectra for single floating wind turbines. Now, the model supports multiple floating wind turbines, each with their own properties and responses, along with mooring lines that run directly between floating wind turbines in the array, meaning that shared mooring lines or fully suspended dynamic power cables can be included. This capability is achieved by setting up an array-level solution of the system mean offsets and assembling the full system matrices for solving the dynamic response. The quasi-static mooring model MoorPy is used to linearize the mooring system properties. To compute the floating wind turbine relative motions, phase offsets are applied to each turbine's response as a function of wave frequency based on the wavelength and relative positions in the array. These differential motions are then applied to mooring system tension Jacobians to compute the tension loads in the shared mooring lines. Overall, the capability provides a frequency-domain analog to the modeling capabilities of the floating support structure in FAST.Farm. Mean offsets and power spectral density plots of responses are compared between RAFT and FAST.Farm to verify the implementation. The results indicate good agreement within the expectations of a frequency-domain modeling approach and suggest correct implementation of the shared mooring aspects. Additionally, a unique comb-like frequency response in the shared mooring line tensions has been observed. This phenomenon has a clear physical basis and may be an important design consideration for future shared mooring systems.

17 WIND ENERGY↗
DOE OSTI2020

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions

In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25 M ⊙

Here, we present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25M ⊙ , are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity 9.6M ⊙ progenitor model is distinct from the temporal evolution of our solar metallicity 15M ⊙ progenitor model and our zero metallicity 25M ⊙ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period ~75 ms after bounce, whereas in the latter two cases high-frequency emission does not commence until ~125 ms after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below ~250 Hz, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at ~125, ~500, and ~250 ms after bounce in the 9.6, 15, and 25M ⊙ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below ~10 Hz. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our 9.6M ⊙ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5–10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today’s detectors. Finally, in our 9.6M ⊙ progenitor model, we see very high-frequency gravitational radiation, extending up to ~2000 Hz. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the 9.6M ⊙ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The physics of desiccation cracks 2: Modeling and prediction of the crack patterns

Here, this paper extends the experimental results of our companion paper by modeling and predicting the onset and pattern formation of desiccation cracks in geomaterials (Ruoyu et al., 2023). Thin-layer silt samples in controlled atmospheric conditions were tested to obtain the surface strain maps with the digital image correlation (DIC) method during dehydration. Support experiments, including consolidation and displacement-controlled triaxial tests, were conducted for the properties of geomaterials. These experimental results were used to validate a viscoplastic theoretical model by comparing cracks and singularities locations that distribute following the Cnoidal wave pattern. A critical value exists in the viscoplastic model that determines the number of singularities, which accurately predicts the number of cracks in the experiment. The critical values contain two crucial parameters: rate sensitivity and λ. Rate sensitivity describes the rate-dependent stress–strain relation, while λ shows the ability of pore pressure redistribution under the external mechanical loading rate. These results provide a new view to analyze the desiccation cracks considering the rate-dependent viscoplasticity.

42 ENGINEERING↗
DOE OSTI2021

Raman Spectroscopy Based On-Line, Real-Time Monitoring to Reduce Composition Uncertainties: Enhanced sensitivity through optimization of Raman Parameters

Optical spectroscopy-based on-line monitoring of Hanford processing streams can enable real-time characterization of chemical composition of process streams and batches, ultimately enabling and enhancing process control. It can provide immediate feedback on process conditions and has the potential to reduce the needed number of grab sample collections, thereby reducing times and costs associated with laboratory processing. Here we discuss the utilization of Raman spectroscopy to quantify multiple target analytes that are common within Hanford tanks and waste processing streams. Analytes include: nitrate, nitrite, carbonate, chromate, sulfate, phosphate, hydroxide, oxalate, ammonia, and aluminate. Most notably in this work, Raman applications to low-concentration streams are explored and optimized. Raman instrument specifications are compared; specifically, the impact of utilizing three different Raman excitation wavelengths, 405, 532, and 671 nm, is discussed. Also, Raman data collection parameters such as collection time and spectral averaging are measured and discussed. Finally, optical libraries of chemical targets were collected using optimized collection parameters and chemometric models were built to automate quantification of chemical targets. These models were validated through application to simulants and real Hanford process samples. Chemometric models performed well on both training and validation sets, suggesting these approaches can be successfully applied to on-line and real-time monitoring of low concentration Hanford processing streams. The use of the combine Raman wave lengths with the enhanced chemometric models for the low concentration streams significantly improved the chemical detection levels and significantly reduced uncertainties of those measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Linear Time-Invariant Models of a Large Cumulus Ensemble

Abstract Methods in system identification are used to obtain linear time-invariant state-space models that describe how horizontal averages of temperature and humidity of a large cumulus ensemble evolve with time under small forcing. The cumulus ensemble studied here is simulated with cloud-system-resolving models in radiative–convective equilibrium. The identified models extend steady-state linear response functions used in past studies and provide accurate descriptions of the transfer function, the noise model, and the behavior of cumulus convection when coupled with two-dimensional gravity waves. A novel procedure is developed to convert the state-space models into an interpretable form, which is used to elucidate and quantify memory in cumulus convection. The linear problem studied here serves as a useful reference point for more general efforts to obtain data-driven and interpretable parameterizations of cumulus convection.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Graph Neural Networks for Surrogate Modeling of Offshore Floating Platforms

Floating offshore wind turbines (FOWTs) present an significant opportunity to increase renewable energy generation. However, significant challenges remain before FOWTs can be widely commercialized and deployed. In particular, hydrodynamic loading on the platforms can stress the overall structure, damage the mooring systems, and impact power generation. Studying these loads is difficult and often relies on computationally expensive models or experiments. In this work, we explore the use of graph neural networks (GNNs) to construct flexible, data-driven surrogates for hydrodynamic loads on platforms. We leverage the natural graph-like structure of offshore wind platform designs to enable the GNN model to learn to approximate the loads for different wave conditions and structural designs. We demonstrate potential uses for the surrogate by performing parameter sweeps and ridge analysis on the trained model to identify the impacts of different wave and structural features on the loads.

floating offshore wind turbines↗
DOE OSTI2022

Review of WEC-Sim development and applications

WEC-Sim (Wave Energy Converter Simulator) is an open-source code for simulating wave energy converters, which has been actively developed and applied to simulate a wide variety of device archetypes, and has become a popular tool since its release. This paper reviewed the development efforts and the usage of WEC-Sim. The publications considered in this study have been broken down into six topic areas, namely feature development, experimental validation, device modeling, control modeling, PTO and grid modeling, and novel applications, which even includes some non-wave energy applications. This review paper has also attempted to recognize the contributions of the broader WEC-Sim development effort, meaning not only the internal WEC-Sim development team but also the external efforts from the academia researchers and technology developers around the world. The growing trend of external applications of WEC-Sim has demonstrated the broader acceptance of the open-source code, and how WEC-Sim has been used in a certain topic area also highlights the potential future development needs.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Improving Common PV Module Temperature Models by Incorporating Radiative Losses to the Sky

PV module operating temperature is the second-most important factor influencing PV system yield–after irradiance–and a substantial contributor to uncertainty in energy system yield predictions. Models commonly used to predict operating temperature in system simulations are based on a simplified energy balance that lumps together different heat loss mechanisms–including radiation–and assumes an overall linear behavior. Radiative heat loss to the sky is usually substantial, but modeling it accurately requires additional information about down-welling long-wave radiation or sky temperature and increases the complexity of temperature model equations. In this work we show how radiative losses to the sky can be separated into two parts to improve the accuracy of modeling without additional complexity. We also predict and demonstrate the variation of these losses at different tilt angles and show that the effective view factor is reduced by the non- isotropic distribution of down-welling long-wave radiation. Finally, we demonstrate substantial reduction in bias (MBE) and scatter (RMSE) when the new radiative loss term is added to the Faiman model using one year of measurements at Sandia National Labs.

14 SOLAR ENERGY↗
DOE OSTI2020

Langevin simulations of the half-filled cubic Holstein model

Over the past several years, reliable quantum Monte Carlo results for the charge density wave transition temperature T cdw of the half-filled two-dimensional Holstein model in square and honeycomb lattices have become available for the first time. Here, exploiting the further development of numerical methodology, here we present results in three dimensions, which are made possible through the use of Langevin evolution of the quantum phonon degrees of freedom. In addition to determining T cdw from the scaling of the charge correlations, we also examine the nature of charge order at general wave vectors for different temperatures, couplings, and phonon frequencies, and the behavior of the spectral function and specific heat.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M ⊙ neutron star to be $R$ $=$ $11.98^{+0.35}_{-0.04}$ km.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration

SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.

Feng, Zongcai↗
DOE OSTI2022

Influence on Structural Loading of a Wave Energy Converter by Controlling Variable-Geometry Components and the Power Take-Off: Preprint

Oceans are a harsh environment and can impose significant loads on deployed structures. Deployment of wave energy converters (WECs) faces a design challenge with apparently contradictory goals. A WEC should be designed to maximize the energy absorbed while ensuring the operating wave condition does not exceed the failure limits of the device itself. Therefore, the loads endured by the support structure are a design constraint for the system. Adaptability to different sea states is, therefore, highly desirable. This work uses a WEC-Sim model of a Variable Geometry Oscillating Wave Energy Converter (VGOSWEC) mounted on a support structure simulated under different wave scenarios. A VGOSWEC resembles a paddle pitching about a fixed hinge perpendicular to the incoming waves-fronts. Therefore, the hinge experiences loads perpendicular to its axis as it maintains its position. The geometry of the VGOSWEC was varied by opening a series of controllable flaps on the pitching paddle when the structure experiences threshold loads. Since opening the flaps lets the waves transmit through the paddle, it is hypothesized that opening the flaps should result in load shedding at the base of the support structure. This load shedding is achieved by reducing the moments about the hinge axis. This work compares the hydrodynamic coefficients, natural periods, and Response Amplitude Operators (RAOs) from completely closed to completely open configurations of the controllable flaps. These comparisons quantify the effects of letting the waves transmit through the VGOSWEC. This work shows that the completely open configuration can reduce the pitch and surge loads on the base of the support structure by as much as 80%. It was observed that at the paddle's resonance frequency, the loads on the structure increased substantially. This increase in loads can be mitigated by a rotational Power Take-off (PTO) damping about the hinge axis. Changing the rotational PTO damping was identified as an additional design parameter that can be used to control the loads experienced by the WEC's support structure.

fixed and floating structures↗