Engineering Papers⌕ Search

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.

At least 109 records · Page 6

Wave Dark Matter

Here we review the physics and phenomenology of wave dark matter: a bosonic dark matter candidate lighter than about 30 eV. Such particles have a de Broglie wavelength exceeding the average interparticle separation in a galaxy like the Milky Way and are, thus, well described as a set of classical waves. We outline the particle physics motivations for such particles, including the quantum chromodynamics axion as well as ultralight axion-like particles such as fuzzy dark matter. The wave nature of the dark matter implies a rich phenomenology: Wave interference gives rise to order unity density fluctuations on de Broglie scale in halos. One manifestation is vortices where the density vanishes and around which the velocity circulates. There is one vortex ring per de Broglie volume on average. For sufficiently low masses, soliton condensation occurs at centers of halos. The soliton oscillates and undergoes random walks, which is another manifestation of wave interference. The halo and subhalo abundance is expected to be suppressed at small masses, but the precise prediction from numerical wave simulations remains to be determined. For ultralight ∼10 −22 eV dark matter, the wave interference substructures can be probed by tidal streams or gravitational lensing. The signal can be distinguished from that due to subhalos by the dependence on stream orbital radius or image separation. Axion detection experiments are sensitive to interference substructures for wave dark matter that is moderately light. The stochastic nature of the waves affects the interpretation of experimental constraints and motivates the measurement of correlation functions. Current constraints and open questions, covering detection experiments and cosmological, galactic, and black hole observations, are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Three-dimensional simulations of large eddies in the compressible mixing layer

Consideration is given to the effect of compressibility on a plane mixing layer that is a prototype free shear layer, amenable to study by numerical simulation and experiment. The full time-dependent compressible Navier-Stokes equations are solved numerically for a temporally evolving mixing layer employing a mixed spectral and high-order finite difference method. Simulations with random initial conditions confirm the prediction of linear stability theory that at high Mach numbers oblique waves grow faster than two-dimensional waves. Simulations are presented of the nonlinear temporal evolution of the most rapidly amplified linear instability waves.

Sandham, N. D.↗

Magnetic Tilt Effect on Externally Driven Electromagnetic Ion Cyclotron (EMIC) Waves

Abstract We examine coupling of fluctuations in the solar wind with electromagnetic ion cyclotron (EMIC) waves in the magnetosphere using an advanced full‐wave simulation code, Petra‐M. Dipole tilt dramatically affects the coupling process. While very little wave power can reach the inner magnetosphere without tilt effects, a tilted dipole field dramatically increases the efficiency of the coupling process. Solar wind fluctuations incident at high magnetic latitude effectively reaches the ground along the field line and mode‐convert to linearly polarized field‐aligned propagating waves at the Alfvén and IIH resonances. Therefore, solar wind compressions efficiently drive linearly polarized EMIC waves when the dipole angle is tilted toward or away from the Sun‐Earth direction.

58 GEOSCIENCES↗

A simulation of synthetic aperture radar imaging of ocean waves

A simulation of radar imaging of ocean waves with synthetic aperture techniques is presented. The modelling is simplistic from the oceanographic and electromagnetic viewpoint in order to minimize the computational problems, yet reveal some of the physical problems associated with the imaging of moving ocean waves. The model assumes: (1) The radar illuminates a one-dimensional, one harmonic ocean wave. (2) The scattering is assumed to be governed by geometrical optics. (3) The radar is assumed to be down-looking, with Doppler processing (range processing is suppressed due to the one-dimensional nature of the problem). (4) The beamwidth of the antenna (or integration time) is assumed to be sufficiently narrow to restrict the specular points of the peaks and troughs of the wave. The results show that conventional processing of the image gives familiar results if the ocean waves are stationary. When the ocean wave dispersion relationship is satisfied, the image is smeared due to the motion of the specular points over the integration time. In effect, the image of the ocean is transferred to the near field of the synthetic aperture.

Swift, C. T.↗

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scoping study of detecting high harmonic fast waves in NSTX-U hot core plasma directly using beam emission spectroscopy

In this paper, we used radio frequency (RF) full wave simulation to explore the feasibility of generating a density beat oscillation formed by high harmonic fast waves (HHFWs) in the National Spherical Tokamak Experiment-Upgrade (NSTX-U) and of measuring it numerically with 2D beam emission spectroscopy (BES). This oscillation is of interest as it can be used to determine the injected RF fields. We began by computing HHFW fields in NSTX-U under a variety of experimental conditions and antenna configurations with Petra-M, a 3D RF wavefield solver. These fields were then used to evaluate the amplitude of δne,beat/ne0 in the volume of measurement, which we predict to be between 5 × 10−10 and 10−9, although it is highly dependent on plasma parameters. Trends in magnetic field, core temperature, core density, and antenna phasing are presented, and the parameters most suitable to this measurement are identified. Finally, we developed a synthetic BES diagnostic to compare the experimental BES signal against.

Rutherford, Grant (ORCID:0000000163620037)↗

High order hybrid numerical simulations of two dimensional detonation waves

In order to study multi-dimensional unstable detonation waves, a high order numerical scheme suitable for calculating the detailed transverse wave structures of multidimensional detonation waves was developed. The numerical algorithm uses a multi-domain approach so different numerical techniques can be applied for different components of detonation waves. The detonation waves are assumed to undergo an irreversible, unimolecular reaction A yields B. Several cases of unstable two dimensional detonation waves are simulated and detailed transverse wave interactions are documented. The numerical results show the importance of resolving the detonation front without excessive numerical viscosity in order to obtain the correct cellular patterns.

Cai, Wei↗

Computer simulation of Alfven waves and double layers along auroral magnetic field lines

A plasma simulation has been developed to model interactions between inertial Alfven waves and double layers and to investigate their relative contributions at auroral particle acceleration. We use a novel one-dimensional particle-in-cell code, with periodic boundary conditions, to model the nonlinear excitation of current-driven weak double layers via free energy supplied by an inertial Alfven wave. Analysis of the simulation output shows that double layers are not the agent primarily responsible for electron acceleration. Rather, the inertial Alfven wave accelerates groups of electrons into a steepening beam as it encounters them. As the beam electrons reenter the main distribution, decelerated by anomalous resistive effects, they are replaced by electrons farther downstream. Hence, the particles do not free-stream over the length of the channel. Furthermore, this wave action persists even when the system is linearly stable to ion-acoustic modes, precluding the possibility that this behavior is brought about by the formation of ion-acoustic double layers.

Silberstein, M.↗

Solitary waves in the resonant phenomenon between a surface gravity wave packet and an internal gravity wave

A two-layer inviscid incompressible fluid system of intermediate depth is considered. A multiple-scales perturbation technique is applied to the basic equations and boundary conditions for a two-layer fluid system to derive a system of weakly nonlinear partial integrodifferential equations governing the resonant interaction between a surface gravity wave packet and an internal gravity wave at an intermediate depth, providing a bridge between the existing shallow and deep fluid theories. The convolution integral term in these equations accounts for the dispersion in the lower-layer fluid. An iterative fast Fourier transform scheme is developed to find solitary wave solutions to this system of equations. The overtaking collision of two pairs of solitary waves, simulated using a spectral method, is found to be inelastic. It is found that the amplitude of the solitary waves changes slightly after the collision. The phase shifts these solitary waves undergo was calculated numerically.

Sepulveda, Nicasio↗

The Application of a Massively Parallel Computer to the Simulation of Electrical Wave Propagation Phenomena in the Heart Muscle Using Simplified Models

The simulation of heart arrhythmia and fibrillation are very important and challenging tasks. The solution of these problems using sophisticated mathematical models is beyond the capabilities of modern super computers. To overcome these difficulties it is proposed to break the whole simulation problem into two tightly coupled stages: generation of the action potential using sophisticated models. and propagation of the action potential using simplified models. The well known simplified models are compared and modified to bring the rate of depolarization and action potential duration restitution closer to reality. The modified method of lines is used to parallelize the computational process. The conditions for the appearance of 2D spiral waves after the application of a premature beat and the subsequent traveling of the spiral wave inside the simulated tissue are studied.

Karpoukhin, Mikhii G.↗

Characterization of Extreme Wave Conditions for Wave Energy Converter Design and Project Risk Assessment

Best practices and international standards for determining n-year return period extreme wave (sea states) conditions allow wave energy converter designers and project developers the option to apply simple univariate or more complex bivariate extreme value analysis methods. The present study compares extreme sea state estimates derived from univariate and bivariate methods and investigates the performance of spectral wave models for predicting extreme sea states at buoy locations within several regional wave climates along the US East and West Coasts. Two common third-generation spectral wave models are evaluated, a WAVEWATCH III® model with a grid resolution of 4 arc-minutes (6–7 km), and a Simulating WAves Nearshore model, with a coastal resolution of 200–300 m. Both models are used to generate multi-year hindcasts, from which extreme sea state statistics used for wave conditions characterization can be derived and compared to those based on in-situ observations at National Data Buoy Center stations. Comparison of results using different univariate and bivariate methods from the same data source indicates reasonable agreement on average. Discrepancies are predominantly random. Large discrepancies are common and increase with return period. There is a systematic underbias for extreme significant wave heights derived from model hindcasts compared to those derived from buoy measurements. This underbias is dependent on model spatial resolution. However, simple linear corrections can effectively compensate for this bias. A similar approach is not possible for correcting model-derived environmental contours, but other methods, e.g., machine learning, should be explored.

16 TIDAL AND WAVE POWER↗

Identifying Ultra Low Frequency Waves in the Lunar Plasma Environment Using Trajectory Analysis and Resonance Conditions

Abstract Recent studies show that localized crustal magnetic fields on the lunar surface can reflect a significant portion of the incoming solar wind protons. These reflected ions can drive a wide range of plasma waves. It is difficult to determine the intrinsic properties of low-frequency waves with single-spacecraft observations, which can be heavily Doppler shifted. We describe a technique to combine trajectory analysis of reflected protons with the Doppler shift and resonance conditions to identify ultralow-frequency waves at the Moon. On 31 January 2014 plasma waves were detected by one of the Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun (ARTEMIS) probes as it approached the lunar wake; these waves were not detected by the second ARTEMIS probe located upstream in the undisturbed solar wind. The observed waves had a frequency below the local ion cyclotron frequency and had right-hand circular polarization in the reference frame of the Moon. By solving the Doppler shift and the cyclotron resonance equations, we determined the conditions for reflected ions to excite the observed waves. Simulated trajectories of reflected ions correspond to ARTEMIS ion observations and support the hypothesis that reflected ions are the primary driver of the waves. By combining trajectory analysis with the resonance conditions, we identify scenarios where ions that satisfy the resonance conditions are present in the right location to generate the observed waves. Using this method, we can uniquely identify the observed waves as upstream propagating right-hand polarized waves, subject to the assumption that they are generated by cyclotron resonance with ions.

Magnetic↗

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

X-ray free-electron laser based dark-field X-ray microscopy: a simulation-based study

Dark-field X-ray microscopy (DFXM) is a nondestructive full-field imaging technique providing three-dimensional mapping of microstructure and local strain fields in deeply embedded crystalline elements. This is achieved by placing an objective lens in the diffracted beam, giving a magnified projection image. So far, the method has been applied with a time resolution of milliseconds to hours. In this work, the feasibility of DFXM at the picosecond time scale using an X-ray free-electron laser source and a pump–probe scheme is considered. Thermomechanical strain-wave simulations are combined with geometrical optics and wavefront propagation optics to simulate DFXM images of phonon dynamics in a diamond single crystal. Using the specifications of the XCS instrument at the Linac Coherent Light Source as an example results in simulated DFXM images clearly showing the propagation of a strain wave.

47 OTHER INSTRUMENTATION↗

Equatorial waves in a stratospheric GCM - Effects of vertical resolution

It is noted that equatorially trapped wave modes, e.g., Kelvin and Rossby-gravity waves, are thought to play a critical role in forcing the quasi-biennial oscillation of the lower tropical stratosphere. The ability of a GCM to simulate these waves is examined and the changes in the wave properties as a function of the vertical resolution of the model are investigated. The dependence of the equatorial wave simulation of vertical resolution is examined utilizing three experiments with vertical grid spacings of about 2.8, 1.4, and 0.7 km. At high vertical resolution, the simulated waves are shown to agree with the available observations.

Boville, Byron A.↗

Coupling between Alfven wave and Kelvin-Helmholtz waves in the low latitude boundary layer

The Kelvin-Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfven velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfv\'en Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfven wave. There are remarkable differences between the primary and the secondary KH waves including wave frequency, the growth rate, and the ratio between transverse and the compressional component. The secondary KH wave energy is concentrated near the shear Alfven wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at lower Mach number. Because the transverse component of the secondary KH waves is stronger than the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

Alfven wave↗

Coupling Between Alfvén Wave and Kelvin–Helmholtz Waves in the Low Latitude Boundary Layer

The Kelvin–Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfvén velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfvén Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfvén wave. There are remarkable differences between the primary and the secondary KH waves, including wave frequency, the growth rate, and the ratio between the transverse and compressional components. The secondary KH wave energy is concentrated near the shear Alfvén wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at a lower Mach number. Because the transverse component of the secondary KH waves is stronger than that of the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

79 ASTRONOMY AND ASTROPHYSICS↗