Detonation Structure in Ethylene/Air-Based Non-Premixed Rotating Detonation Engine
Not provided.
Engineering↗
SEARCH · Engineering Papers
Results for “simulated waves”
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.
HERO WEC V1.0 - WEC-Sim Model (July 2024)
**This submission supersedes submission MHKDR-483** This submission file contains the files that are needed to simulate NREL's HERO WEC (hydraulic and electric reverse osmosis wave energy converter). This requires the user to have already installed WEC-Sim. In addition to the standard toolboxes that are required to run WEC-Sim the user will also need the Simscape Fluids and Simscape Driveline packages. The zip file (HERO_V1_WECSim_2024.zip) contains the following: - HERO_HPTO_2024.slx: Simulink-based WEC Sim model of the first gen (V1.0) Hydraulic PTO (power take-off) that was designed for the HERO WEC. This model has been updated since submission #483 based on in-laboratory experimental results. - wecSimInputFile.m: Input file needed to run the model - userDefinedFunctionsMCR.m: MCR (multi condition run) script that is needed if a use wants to simulate multiple wave conditions. - geometry (folder): Includes the geometry file that is needed for visualization - hydroData (folder): Includes the required WAMIT data to run WEC-Sim -HydVisualization.mlx: Visualization script to plot simulation results (not needed to run)
16 TIDAL AND WAVE POWER↗
Influence of Time and Frequency Domain Wave Forcing on the Power Estimation of a Wave Energy Converter Array
Industry-specific tools for analyzing and optimizing the design of wave energy converters (WECs) and associated power systems are essential to advancing marine renewable energy. This study aims to quantify the influence of phase information on the device power output of a virtual WEC array. We run the phase-resolving wave model FUNWAVE-TVD (Total Variation Diminishing) to generate directional waves at the PacWave South site offshore from Newport, Oregon, where future WECs are expected to be installed for testing. The two broad cases presented correspond to mean wave climates during warm months (March–August) and cold months (September–February). FUNWAVE-TVD time series of sea-surface elevation are then used in WEC-Sim, a time domain numerical model, to simulate the hydrodynamic response of each device in the array and estimate their power output. For comparison, WEC-Sim is also run with wave energy spectra calculated from the FUNWAVE-TVD simulations, which do not retain phase information, and with wave spectra computed using the phase-averaged model Simulating WAves Nearshore (SWAN). The use of spectral data in WEC-Sim requires a conversion from frequency to time domain by means of random superposition of wave components, which are not necessarily consistent because of the linear assumption implicit in this method. Thus, power response is characterized by multiple realizations of the wave climates.
13 HYDRO ENERGY↗
Shock compression/release of magnesium single crystals along a low-symmetry orientation: Role of basal slip
To gain insight into the relative contributions of different plastic deformation mechanisms, particularly basal slip, for shocked hexagonal close-packed (hcp) metals, magnesium (Mg) single crystals were subjected to shock compression and release along a low-symmetry (LS) orientation to 1.9 and 4.8 GPa elastic impact stresses. LS-axis is a “non-specific” direction resulting in propagation of quasi-longitudinal and quasi-shear waves. Wave profiles, measured using laser interferometry, show a small elastic wave followed by two plastic waves in compression; release wave profiles exhibited a structured response for the higher stress and a smooth response for the lower stress. The LS-axis wave profiles are significantly different than profiles published previously for c- and a -axes, demonstrating that Mg single crystals exhibit strong anisotropy under shock compression/release. Numerical simulations, using a time-dependent anisotropic modeling framework, show that shock wave loading along the LS-axis involves simultaneous operation of multiple deformation mechanisms. Shock compression along LS-axis is dominated by basal slip while prismatic and pyramidal I {$10\bar{1}1$} $\langle$ $11\bar{2}3$ $\rangle$ slip play a smaller role; coupling between longitudinal and shear deformations was observed. The unloading response is dominated by basal slip with some contribution from prismatic slip; pyramidal I slip is not activated. The present results, unlike results obtained for c - and a -axes, show that the deformation mechanism observed under quasi-static loading conditions along LS-axis is not sufficient to determine shock response along this orientation. Finally, although requiring numerical simulations for wave analysis, shock propagation along a LS orientation provides new insights into the plastic deformation response of hcp metal single crystals.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Propagation of ion cyclotron emission in the DIII-D tokamak
Experiments in L- and H-mode plasmas performed on the DIII-D tokamak explored ion cyclotron emission (ICE) propagation via the recently upgraded ICE diagnostic. The distance between the plasma and the outer wall was scanned to alter the evanescent region, which was seen to impact the detection of runaway electron-driven whistler modes in a comparable frequency range to typical ICE harmonics (f≳10fci) [Heidbrink et al., Plasma Phys. Controlled Fusion 61, 014007 (2018)]. In L-mode plasmas, core ICE did not exhibit clear trends as the plasma-wall distance was altered. Instead, inadvertent changes to the fast ion density correlated with different ICE fine structures; the most spectral bands (spaced ∼150–200 kHz apart) were observed at the highest fast ion density, and a just single band when nfast decreased by ∼ 40%. The full-wave Petra-M code simulated core ICE propagation through the plasma and to ICE probes, agreeing with experimental findings that these are likely compressional fast waves, which propagate to probe locations with no dependence on plasma-wall gap. Edge ICE in H-mode plasmas was similarly unaffected by plasma-wall distance. Unlike in core ICE cases, there were negligible changes to the edge fast ion distribution and edge bulk plasma profiles. However, changes in ELM frequency seemed most closely tied to ICE behavior. Ultimately, core and edge ICE did not exhibit clear dependencies on the plasma-wall distance, and ICE harmonics were observed with even the largest plasma-wall gaps. This independence bodes well for similar diagnostics in future tokamaks, which might be placed further from high radiation areas without adversely impacting measurement capabilities.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Study of spatiotemporal development of global distribution of magnetospheric ELF/VLF waves using ground-based and satellite observations, and RAM-SCB simulations, for the March and November 2017 storms
Explore the source record for details and available documents.
54 ENVIRONMENTAL SCIENCES↗
MBARI-WEC September and October 2022 Field Data
This data is needed to simulate a model of the MBARI-WEC (Monterey Bay Aquarium Research Institute, Wave Energy Converter device) in a simulation environment (e.g. Gazebo) for 56 observation dates in the time between September and October 2022, and to compare the simulation outputs to the corresponding field data of the physical MBARI-WEC. There were 50 observations chosen in Sept and 6 observations in Oct. To help understand terms below, a summary of the system can be found at the github link in the downloads section below. The Gazebo MBARI-WEC model is also provided, should users wish to simulate using this platform. There are 4 *mat files included. ................................................................................................................................................................................................................................... Spectrum and Simulation Inputs: September2022_spectrum_siminputs.mat and October2022_spectrum_siminputs.mat has data needed for simulation inputs in table format. These include the ocean spectrum for an observation and operating parameters of the MBARI-WEC during that observation. They are organized as rows representing an observation and columns representing data. For example, for the September *mat there are 50 rows. The first 7 columns are Datetime, sig_waveheight, peak_period, mean_period, heaveconedoor_status, pistonpos_mean, and scale_factor: - Datetime is the date and time the observation occurred in PST - sig_waveheight is the significant wave height of the ocean spectrum during that observation in meters - peak_period is the peak period of the ocean spectrum during that observation in seconds - mean_period is the mean period of the ocean spectrum during that observation in seconds - heaveconedoor_status is the status of the heave cone doors where 0 represents the doors are open and 1 represents they are closed - pistonpos_mean is the mean position of the PTO ram (piston) in meters - scale_factor is an additional factor of 0.5 --1.4 applied to a default damping relationship The next columns are data needed to represent the ocean spectrum. First are the frequencies [Hz] labeled as "f0-f38", then the variance density [m2/Hz] labeled as "vardens0-vardens38". October2022_spectrum_siminputs.mat follows as a similar format as above, but includes a larger amount of ocean spectrum frequencies and variance density elements. ................................................................................................................................................................................................................................... Field data: The field data is found in MBARIWEC_septdata.mat and MBARIWEC_octdata.mat for the observations of September and October, respectively. These contain data in a struct format. The struct contains the following fields for each observation: PC_BattCurr, PC_LoadCurr, PC_RPM, PC_Voltage, SC_Range, SC_Velocity, DateTime, where: - PC_BattCurr is the current flowing to or from the onboard batteries in Amps - PC_LoadCurr is the current flowing to the load dump in Amps - PC_RPM is the electric/hydraulic motor shaft speed (directly coupled) in RPM - PC_Voltage is the bus voltage at the power converter in Volts - SC_Range is the PTO ram (piston) position in meters where 0 is fully retracted and 2.03 is fully extended - SC_Velocity is the PTO ram (piston) velocity in meters/sec - DateTime is the date and time of the sampled field data in each observation in PST - Electric Power is equal to: PC_Voltage*(PC_BattCurr + PC_LoadCurr) in Watts For example, upon loading MBARIWEC_octdata.mat, the aforementioned fields would be loaded, each with {6x1} cells for the 6 observations chosen in October. Within the first cell of e.g. SC_Range would be sampled data representing the field data of the MBARIWEC PTO piston position for say, one hour, of the first October observation. The corresponding field DateTime would...
16 TIDAL AND WAVE POWER↗
Implementation and Verification of Cable Bending Stiffness in MoorDyn: Preprint
Cable bending stiffness capabilities have been added to MoorDyn to enable modeling of dynamic power cables. The relatively large motions experienced by floating wind turbines and wave energy converters pose a challenge for power cables, whose internal components provide significant bending resistance and are sensitive to deformation. The behavior and associated design considerations of power cables in these highly dynamic applications make coupled analysis relevant for design. MoorDyn is a lumped-mass mooring dynamics model that is used for floating wind energy simulation (often coupled with OpenFAST) and floating wave energy converter simulation (often coupled with WEC-Sim). Its existing line elasticity formulation considers axial stiffness only. To properly capture the dynamics of power cables, a cable bending stiffness model has been added that approximates cable curvature based on the difference in tangent vectors of adjacent line elements. The resulting bending moment is realized by applying forces on adjacent nodes, leaving the underlying lumped-mass formulation unchanged. The bending stiffness implementation is verified in static conditions against analytical solutions and then in a dynamic power cable scenario in comparison with OrcaFlex. The dynamic scenario uses prescribed motions and includes wave loadings on the cable. Results indicate correct implementation of bending stiffness and show close agreement with OrcaFlex.
50 EE - Wind and Water Power Program - Water (EE-4↗
Machine learning the dynamics of quantum kicked rotor
Highlights: • Quantum kicked rotor. • Quantum phase transition. • Machine learning. • LSTM network. • Convolutional neural network. Using the multilayer convolutional neural network (CNN), we can detect the quantum phases in random electron systems, and phase diagrams of two and higher dimensional Anderson transitions and quantum percolations as well as disordered topological systems have been obtained. Here, instead of using CNN to analyze the wave functions, we analyze the dynamics of wave packets via long short-term memory network (LSTM). We adopt the quasi-periodic quantum kicked rotors, which simulate the three and four dimensional Anderson transitions. By supervised training, we let LSTM extract the features of the time series of wave packet displacements in localized and delocalized phases. We then simulate the wave packets in unknown phases and let LSTM classify the time series to localized and delocalized phases. We compare the phase diagrams obtained by LSTM and those obtained by CNN.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Focused wave interactions with floating structures: a blind comparative study
The paper presents results from the Collaborative Computational Project in Wave Structure Interaction (CCP-WSI) Blind Test Series 2. Without prior access to the physical data, participants, with numerical methods ranging from low-fidelity linear models to fully non-linear Navier-Stokes (NS) solvers, simulate the interaction between focused wave events and two separate, taut-moored, floating structures: a hemispherical-bottomed cylinder and a cylinder with a moonpool. The ‘blind’ numerical predictions for heave, surge, pitch and mooring load, are compared against physical measurements. Dynamic time warping is used to quantify the predictive capability of participating methods. In general, NS solvers and hybrid methods give more accurate predictions; however, heave amplitude is predicted reasonably well by all methods; and a WEC-Sim implementation, with CFD-informed viscous terms, demonstrates comparable predictive capability to even the stronger NS solvers. Large variations in the solutions are observed (even among similar methods), highlighting a need for standardisation in the numerical modelling of WSI problems.
42 ENGINEERING↗
Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations
One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.
1D assumption↗
Source Discretization Near Curvilinear Mesh Refinements in One Dimeions (Summer Internship Final Report)
Lawrence Livermore National Lab developed the seismic wave code SW4 for earthquake simulations, and has recently added the ability to specify mesh-refinement interfaces in the curvilinear mesh. The primary purpose of this is to allow for high resolution near the surface, with relatively fewer total points in the mesh leading to shorter computation time. Presently, the code does not allow for point-source forcing terms (commonly used in seismic simulations) to be located near these curvilinear mesh refinement interfaces. The goal of this internship was to investigate the discretization of point-source forcing terms in the seismic wave equation when the sources are located near the interfaces between two curvilinear meshes. Ultimately this will be incorporated into SW4 to allow for more realistic simulations to be carried out in this faster implementation.
97 MATHEMATICS AND COMPUTING↗
Parametric Instability of Alfvén Waves and Wave Packets in Periodic and Open Systems
The parametric decay instability of Alfvén waves has been widely studied, but few investigations have examined wave packets of finite size and the effect of different boundary conditions on the growth rate. In this paper, we perform a linear analysis of circular and arc-polarized wave trains and wave packets in periodic and open boundary systems in a low- β plasma. We find that both types of wave are 3–5 times more stable in open boundary conditions compared to periodic. Additionally, once the wave packet width ℓ becomes smaller than the system size L , the growth rate decreases nearly with a power law γ ∝ ℓ / L . This study demonstrates that the stability of a pump wave cannot be separated from the laboratory settings, and that the growth rate of daughter waves depends on the conditions downstream and upstream of the pump wave and on the fraction of volume it fills. Our results can explain simulations and experiments of localized Alfvén waves. They also suggest that Alfvénic fluctuations in the solar wind, including sharp impulses known as switchbacks, can be more stable than traditional theory suggests depending on wind conditions.
Alfven waves↗
Weak shock compaction on granular salt
This study conducted integrated experiments and computational modeling to investigate the speeds of a developing shock within granular salt and analyzed the effect of various impact velocities up to 245 m/s. Experiments were conducted on table salt utilizing a novel setup with a considerable bore length for the sample, enabling visualization of a moving shock wave. Experimental analysis using particle image velocimetry enabled the characterization of shock velocity and particle velocity histories. Mesoscale simulations further enabled advanced analysis of the shock wave’s substructure. In simulations, the shock front’s precursor was shown to have a heterogeneous nature, which is usually modeled as uniform in continuum analyses. The presence of force chains results in a spread out of the shock precursor over a greater ramp distance. With increasing impact velocity, the shock front thickness reduces, and the precursor of the shock front becomes less heterogeneous. Furthermore, mesoscale modeling suggests the formation of force chains behind the shock front, even under the conditions of weak shock. This study presents novel mesoscale simulation results on salt corroborated with data from experiments, thereby characterizing the compaction front speeds in the weak shock regime.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Effects of sub‐mm cylindrical voids on detonation performance in PBX 9501
Abstract Internal features of varying scale and geometry are always present in explosives systems. Below a critical length, dependent on the explosive, these features can operate as a driving force for energy concentration and reaction, known as hotspots. At larger length scales, internal features result in jetting and detonation wave shaping, allowing for bulk work to be done by the explosive as is seen in shape charges. To date, a large volume of work has been performed to simulate hotspot ignition and large‐scale wave shaping. However, little experimental data exists on the effects of features in intermediate length scales, 0.1 to 1 millimeter. It has been observed in many tests that these small‐scale features can influence the high explosive (HE) performance and in some cases cause substantial damage to adjacent systems. This work provides quantitative data on the effects of machined voids moderately above more typical hotspot lengths, 0.3–0.8 mm in diameter, in PBX 9501 pressed to 1.785 g cm −3 ±2.5 mg cm −3 . Streak imaging was used to visualize void collapse, jet velocity, re‐initiation time, and wave shape evolution. Delay in detonation front propagation time was found to be linearly dependent on the void diameter and jet velocity was found to be independent within the tested range and resolution. Cut‐back experiments were used to investigate wave shape distortion and evolution downstream of the void, showing consistent growth and decay shapes across all void sizes. Simulations using CTH, a hydrocode by Sandia National Laboratory, were used to investigate void collapse showing agreement with the trends of experimental results but yielded inaccurate wave shape development delay values. Both simulation and experimental results identified several re‐initiation mechanisms with jetting and subsequent double shocking of localized HE being the dominant mechanisms for the void sizes that were studied.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
The Relationship Between African Easterly Waves and Tropical Cyclones in Historical and Future Climates in the HighResMIP‐PRIMAVERA Simulations
Abstract The deadly and destructive nature of tropical cyclones (TCs) makes understanding their response to future climate change of the utmost importance. TC genesis hinges on multiple factors, including an initial disturbance. African easterly waves (AEWs) have been shown to serve as such disturbances for TCs developing in the North Atlantic. It is therefore crucial to understand the relationship between AEWs and TCs and how this relationship may be affected by climate change. In this study, we examine the AEW‐TC relationship in historical and future climates using three models from the HighResMIP PRIMAVERA simulations. The AEWs and TCs were tracked in the model data using objective tracking algorithms, and AEW and TC tracks were then matched together if they were close to each other in space and time. The strength of the AEWs was measured using the eddy kinetic energy and the curvature vorticity of the waves. TC strength and intensity were measured using potential intensity and lifetime maximum 10 m windspeed. We found that future changes in the frequency of AEWs are not a good indicator of future TC activity. However, AEW strength, as well as environmental conditions conducive to strong TCs, are good indicators of AEWs that develop into TCs in both historical and future climates.
54 ENVIRONMENTAL SCIENCES↗
TEAMER Technical Support of Ramboll's for Numerical Modeling of WECs to Support OES Task 10: Cooperative Research and Development (Final Report)
National Technology & Engineering Solutions of Sandia, LLC (NTESS) in collaboration with the National Renewable Energy Laboratory (NREL), and with guidance from Ramboll, will perform fluid dynamics simulations to support of The Ocean Energy Systems (OES) Energy Technology Collaboration Program Task 10 Wave Energy Converters (WEC) Modelling Verification and Validation effort. Specific numerical simulations include the fixed device wave impingement studies to benchmark the performance of simulation techniques against physical testing results.
16 TIDAL AND WAVE POWER↗
Langmuir turbulence and filament frontogenesis in the oceanic surface boundary layer
Submesoscale currents, small-scale turbulence and surface gravity waves co-exist in the upper ocean and interact in complex ways. To expose the couplings, the frontogenetic life cycle of an idealized cold dense submesoscale filament interacting with upper ocean Langmuir turbulence is investigated in large-eddy simulations (LESs) based on the incompressible wave-averaged equations. The simulations utilize large domains and fine meshes with $6.4\times 10^{9}$ grid points. Case studies are made with surface winds or surface cooling with waves oriented in across-filament (perpendicular) or down-filament (parallel) directions relative to the two-dimensional filament axis. The currents $u$ , $v$ and $w$ are aligned with the across-filament, down-filament and vertical directions, respectively. Frontogenesis is induced by across-filament Lagrangian secondary circulations in the boundary layer, and it is shown to be strongly impacted by surface waves, in particular the propagation direction relative to the filament axis. In a horizontally heterogeneous boundary layer, surface waves induce both mean and fluctuating Stokes-drift vortex forces that modify a linear, hydrostatic turbulent thermal wind (TTW) approximation for momentum. Down-filament winds and waves are found to be especially impactful, they significantly reduce the peak level of frontogenesis by fragmenting the filament into primary and secondary down-welling sites in a broad frontal zone over a width ${\sim}500~\text{m}$ . At peak frontogenesis, opposing down-filament jets $\langle v\rangle$ overlie each other resulting in a vigorous vertical shear layer $\unicode[STIX]{x2202}_{z}\langle v\rangle$ with large vertical momentum flux $\langle v^{\prime }w^{\prime }\rangle$ . Filament arrest is induced by a lateral shear instability that generates horizontal momentum flux $\langle u^{\prime }v^{\prime }\rangle$ at low wavenumbers. The turbulent vertical velocity patterns, indicative of coherent Langmuir cells, change markedly across the horizontal domain with both across-filament and down-filament winds under the action of submesoscale currents.
Mechanics↗