SEARCH · Search NASA

Results for “electromagnetic wave propagation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Nonlinear self-focusing in strongly magnetized pair plasma

An intense radiation field can modify plasma properties and the corresponding refractive index and lead to nonlinear propagation effects such as self-focusing. In this work, we estimate the corresponding effects in pair plasmas for circularly polarized waves, in both unmagnetized and strongly magnetically dominated cases. First, in the unmagnetized pair plasma the ponderomotive force does not lead to charge separation but to density depletion. Second, for astrophysically relevant plasmas of pulsar magnetospheres [and possible loci of fast radio bursts (FRBs)], where the cyclotron frequency $ω_B$ dominates over the plasma frequency $ω_p$ and the frequency of the electromagnetic wave $ω_B \gg ω_p, ω$, we show that (i) there is virtually no nonlinearity due to changing effective mass in the field of the wave; (ii) the ponderomotive force is $F^{(B)}_p = –m_ec^2/4B^2_0∇E^2$, which is reduced by a factor $(ω/ω_B)^2$ if compared to the unmagnetized case ($B_0$ is the external magnetic field and $\textit{E}$ is the electric field of the wave); and (iii) for a radiation beam propagating along a constant magnetic field in the pair plasma with density $n_±$, the ponderomotive force leads to the appearance of circular currents that lead to a decrease of the field within the beam by a factor $ΔB/B_0 = 2πn_±m_ec^2E^2 /B^4_0$. Applications to the physics of FRBs are discussed; we conclude that for the parameters of FRBs, the dominant magnetic field completely suppresses nonlinear self-focusing or filamentation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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 OSTI2025

Prediction of the Wave Normal Angle of Proton‐Band EMIC Waves Near Geosynchronous Orbit

We investigate how the wave normal angle (WNA) and polarization of proton-band electromagnetic ion cyclotron (EMIC) waves change as they travel from their source to Earth. This paper marks a significant milestone as the first full-wave simulation of proton-band EMIC waves reflecting from the ionosphere. Our findings show that the WNA can change rapidly during propagation, primarily due to plasma inhomogeneities, such as variations in the Alfvén speed. The wave polarization is strongly related to the WNA, consistent with theory. Newly generated EMIC waves near the equator propagate with a WNA of 0° , then the WNA gradually shifts to 90° as they move toward Earth. In contrast, reflecting waves having 90° of WNA at Earth maintain a relatively larger WNA even near the magnetic equator. As a result, only the newly generated waves close to the source, where the magnetic latitude is less than approximately 20° , show left-handed polarization, while linear polarization remains dominant throughout the rest of the propagation.

EMIC wave↗

DOE OSTI2020

Real-time time-dependent density functional theory in Quantum Espresso and related codes

A set of subroutines that interact with an existing computer program, Quantum Espresso (QE), that takes a set of electronic orbitals from QE and propagates them forward in time in response to a number of perturbations. These perturbations include the motion of ions and a wide spectrum of applied electromagnetic fields (e.g., static magnetic fields, slowly-varying electric fields, x-rays, etc.). From the time evolution of the electronic orbitals comes time-varying values of physical observables like the electronic charge and current densities and atomic forces. These are further post-processed to yield properties of interest like stopping powers, conductivities, and the dynamic structure factor. The electronic orbitals are described in a plane wave basis and the equations of motion are carried out using a unitary or approximately unitary time propagation algorithm. Related work will include subroutines for interfacing this capability, as well as others that already exist within QE, with a fork of QE that is more extensible and modular.SAND2019-12611 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Baczewski, Andrew↗

DOE OSTI2021

Coherent Electromagnetic Emission from Relativistic Magnetized Shocks

Relativistic magnetized shocks are a natural source of coherent emission, offering a plausible radiative mechanism for fast radio bursts (FRBs). We present first-principles 3D simulations that provide essential information for the FRB models based on shocks: the emission efficiency, spectrum, and polarization. The simulated shock propagates in an e ± plasma with magnetization σ > 1 . The measured fraction of shock energy converted to coherent radiation is ≃ 10 -3 σ -1 , and the energy-carrying wave number of the wave spectrum is ≃ 4ω c /c, where ω c is the upstream gyrofrequency. The ratio of the O-mode and X-mode energy fluxes emitted by the shock is ≃ 0.4σ -1 . The dominance of the X mode at σ >> 1 is particularly strong, approaching 100% in the spectral band around 2ω c . We also provide a detailed description of the emission mechanism for both X and O modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Photon scattering by a 4π -spherically-focused ultrastrong electromagnetic wave

The scattering of high-power probe laser pulses by tightly focused ultrastrong laser pulses has been investigated through a semiclassical approach using a perturbative method based on the Born approximation. Under a 4π-spherically-focused ultrastrong light field, the electric permittivity and magnetic permeability tensors for vacuum are calculated from the Euler-Heisenberg Lagrangian to show the nonlinear birefringent property of vacuum. And, from permittivity and permeability tensors, the scattering potential is derived for the Born approximation. The first-order solution of the Born approximation is taken as an electric field scattered from the 4π-spherically-focused laser pulse when a probe laser pulse propagates through the focused laser field. The differential cross section of the nonlinear birefringent vacuum is derived, and the number of photons scattered from the nonlinear birefringent vacuum is analyzed in the laser power range of 10-1000 PW.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Machined phononic crystals to block high-order Lamb waves and crosstalk in through-metal ultrasonic communication systems

For systems that require complete metallic enclosures (e.g., containment buildings for nuclear reactors), it is impossible to access interior sensors and equipment using standard electromagnetic techniques. A viable way to communicate and supply power through metallic barriers is the use of elastic waves and ultrasonic transducers, introducing several design challenges that must be addressed. Specifically, the use of multiple communication channels on the same enclosure introduces an additional mechanism for signal crosstalk between channels: guided waves propagating in the barrier between channels. This work numerically and experimentally investigates a machined phononic crystal to block MHz Lamb wave propagation between ultrasonic communication channels, greatly reducing wave propagation and the resulting crosstalk voltage. Blind grooves are machined into one or both sides of a metallic barrier to introduce a periodic unit cell, greatly altering the guided wave dispersion in the barrier. Numerical simulations are used to determine a set of groove geometries for testing, and experiments were performed to characterize the wave-blocking performance of each design. The best-performing design was tested using piezoelectric transducers bonded to the barrier, showing a 14.4 dB reduction in crosstalk voltage. Overall, the proposed periodic grooving method is a promising technique for completely isolating ultrasonic power/data transfer systems operating in a narrow frequency range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗

DOE OSTI2025

A mixed-integer PDE-constrained optimization formulation for constructing electromagnetic cloaks with multiple materials

We study the design of an electromagnetic cloak from multiple materials with an additional constraint on the mass of the cloak. Our problem is an example of a topology optimization problem, and we formulate this problem as a mixed-integer partial-differential equation constrained optimization (MIPDECO) problem, where Maxwell’s equation models the propagation of the wave through the cloak and surrounding medium. We use binary variables to model the assignment of the different materials, and their relevant properties (permittivity and density). The mass constraint adds a nontrivial constraint to this problem. We propose a two-phase strategy to solve this problem. In the first phase, we solve a continuous relaxation, and then propose a new variant of the feasibility pump that exploits the structure of the PDE to obtain an initial integral solution candidate. In the second phase, we use a trust-region approach to improve this incumbent. We also consider a continuation or mesh-sequencing approach to find better solutions faster on consecutively finer meshes. We present detailed numerical results to illustrate the effectiveness of our approaches for constructing multi-material cloaks with a mass constraint.

Calculus of Variations and Optimization↗

DOE OSTI2021

Interaction of radio frequency waves with cylindrical density filaments: scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas , vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci. , vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

Physics↗

DOE OSTI2022

Interaction of radio frequency waves with cylindrical density filaments -- scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas, vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci., vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Verification of a fully implicit particle-in-cell method for the  v ∥ -formalism of electromagnetic gyrokinetics in the XGC code

A fully implicit particle-in-cell method for handling the v ∥ -formalism of electromagnetic gyrokinetics has been implemented in XGC. By choosing the v ∥ -formalism, here we avoid introducing the nonphysical skin terms in Ampère's law, which are responsible for the well-known “cancellation problem” in the p ∥ -formalism. The v ∥ -formalism, however, is known to suffer from a numerical instability when explicit time integration schemes are used due to the appearance of a time derivative in the particle equations of motion from the inductive component of the electric field. Here, using the conventional δf scheme, we demonstrate that our implicitly discretized algorithm can provide numerically stable simulation results with accurate dispersive properties. We verify the algorithm using a test case for shear Alfvén wave propagation in addition to a case demonstrating the ion temperature gradient-kinetic ballooning mode (ITG-KBM) transition. The ITG-KBM transition case is compared to results obtained from other δf gyrokinetic codes/schemes, whose verification has already been archived in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Oblique propagation and temperature effects on the resonant right-hand ion beam instability

The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Wind and Temperature Consensus at Horn Point, HU-Beltsville, Piney Run (Maryland) in support of CoURAGE

The Maryland Department of the Environment (MDE) operates a ground-based atmospheric profiling network consisting of collocated radar wind profilers (RWP) and radio acoustic sounding systems (RASS) as part of its Ambient Air Monitoring Program. This network provides continuous observations of wind and temperature structure in the lower troposphere to support air quality forecasting, regulatory analysis, and atmospheric research. The network currently includes three fixed sites across Maryland: Horn Point (HP, lower eastern shore) [38.587525°,-76.141006°], Howard University-Beltsville (HUB, central Maryland) [39.055277°, -76.878632°], and Piney Run (PR, western Maryland) [39.705950°, -79.012000°] The network is designed to capture regional variability in atmospheric transport and boundary-layer processes. These systems measure vertical profiles of horizontal wind speed and direction using Doppler radar techniques, with observations typically spanning from ~100 m above ground level up to approximately 2.5–4 km. Measurements are derived from the Doppler shift of backscattered electromagnetic signals, enabling retrieval of wind vectors at multiple altitudes with high temporal resolution (e.g., 30-minute averages reported every 6 minutes). Each radar wind profiler is paired with a Radio Acoustic Sounding System (RASS) to provide profiles of virtual temperature in the lower atmosphere (~100–200 m AGL) by measuring the propagation speed of acoustic waves. Together, the RWP/RASS system yields a coupled data set of thermodynamic and kinematic atmospheric structure, including additional parameters such as vertical velocity, radial velocity, signal-to-noise ratio, and spectral width for advanced analysis. There are two types of files for each station: wind data (files with a "w" prefix) and virtual temperature RASS data (files with a "t" prefix). The wind data files are in the format wYYDDD.cns, where YY is the 2-digit year and DDD is the day of the year. The RASS virtual temperature data files are in the format tYYDDD.cns. Each record has the following header structure: Line 1 : Station Name RASS files Line 2 : RASS rev DeTect_2.0, WINDS files Line 2 : WINDS rev ATI 5.1 Line 3 : N latitude, W longitude, and site elevation (m) Line 4 : Date and begin time of consensus: yy mm dd hh mn ss plus # minutes to add to get UTC Line 5 : Consensus averaging time (minutes); number of beams; number of range gates Line 6 : Number of records required to make consensus (num) total number of records (tot) and the consensus window size (m/s) in the format: num:tot (window) RASS files Line 7 : no. of coded cells, no. of spec, pulse width (ns), and inter-pulse period (µs), WINDS files Line 7 : No. of coded cells, no. of spectra, pulse width (ns), and inter-pulse period (µs), each with a pair of values: first value is for oblique beams, second for vertical RASS files Line 8 : Full scale Doppler value (m/s) Delay to first gate (ns) Number of gates Spacing of gates (ns), WINDS files Line 8 : Full scale Doppler velocity (m/s), oblique and vertical Vertical correction applied to oblique beams? (0 = no, 1 = yes) Delay to first gate (ns), oblique and vertical Number of gates, oblique and vertical Spacing of gates (ns), oblique and vertical Line 9 : Azimuth and elevation (9s indicate vertical beam not used) RASS files Line 10, values : HT = Height above ground (km), T = Uncorrected virtual temperature consensus (deg C), Tc = Corrected virtual temperature consensus (deg C), W = Vertical wind consensus (9s indicate vertical beam not used, w-component, positive upward, m/s), CNT = Number of records that made consensus (for the 3 values in same order), SNR = Average signal to noise ratio (dB) of records in consensus (same order) WINDS files Line 10, values : HT = Height above ground (km), SPD = Wind speed (m/s), DIR = Wind direction (deg E of N from N), RAD = Radial velocities for each beam (m/s) in order given in azimuth and elevation line (positive toward radar; 9s indicate vertical beam not used, CNT = Number of records that made consensus, SNR = Average signal to noise ratio (dB) of records in consensus

{"wind speed and direction",temperature}↗

DOE OSTI2024

Plasmons in the Kagome metal CsV 3 Sb 5

Plasmon polaritons, or plasmons, are coupled oscillations of electrons and electromagnetic fields that can confine the latter into deeply subwavelength scales, enabling novel polaritonic devices. While plasmons have been extensively studied in normal metals or semimetals, they remain largely unexplored in correlated materials. In this paper, we report infrared (IR) nano-imaging of thin flakes of CsV3Sb5, a prototypical layered Kagome metal. We observe propagating plasmon waves in real-space with wavelengths tunable by the flake thickness. From their frequency-momentum dispersion, we infer the out-of-plane dielectric function that is generally difficult to obtain in conventional far-field optics, and elucidate signatures of electronic correlations when compared to density functional theory (DFT). We propose correlation effects might have switched the real part of ϵ c from negative to positive values over a wide range of middle-IR frequencies, transforming the surface plasmons into hyperbolic bulk plasmons, and have dramatically suppressed their dissipation.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Electromagnetic power of lightning superbolts from Earth to space

Lightning superbolts are the most powerful and rare lightning events with intense optical emission, first identified from space. Superbolt events occurred in 2010-2018 could be localized by extracting the high energy tail of the lightning stroke signals measured by the very low frequency ground stations of the World-Wide Lightning Location Network. Here, we report electromagnetic observations of superbolts from space using Van Allen Probes satellite measurements, and ground measurements, and with two events measured both from ground and space. From burst-triggered measurements, we compute electric and magnetic power spectral density for very low frequency waves driven by superbolts, both on Earth and transmitted into space, demonstrating that superbolts transmit 10-1000 times more powerful very low frequency waves into space than typical strokes and revealing that their extreme nature is observed in space. We find several properties of superbolts that notably differ from most lightning flashes; a more symmetric first ground-wave peak due to a longer rise time, larger peak current, weaker decay of electromagnetic power density in space with distance, and a power mostly confined in the very low frequency range. Their signal is absent in space during day times and is received with a long-time delay on the Van Allen Probes. These results have implications for our understanding of lightning and superbolts, for ionosphere-magnetosphere wave transmission, wave propagation in space, and remote sensing of extreme events.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Net energy gain in direct laser acceleration due to enhanced dephasing induced by an applied magnetic field

Even in the situation where an electron interacts with a single plane wave, the well-known dynamical adiabaticity can be broken when an applied magnetic field is present, which will act to increase the dephasing rate of the electron during the interaction. Here we demonstrate this for the case where there is a uniform static magnetic field which is oriented either parallel or perpendicular to the electric field of the incident plane wave, and perpendicular to the direction of its propagation. The described energy gain phenomenon has direct relevance to laser-plasma interactions that involve external magnetic fields generated by laser-driven capacitor coils.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Synchrotron Polarization Signatures of Surface Waves in Supermassive Black Hole Jets

Abstract Supermassive black holes in active galactic nuclei are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of nonthermal electrons. In this Letter, we present a global general relativistic magnetohydrodynamical (GRMHD) simulation of a magnetically arrested disk (MAD). After the simulation reaches the MAD state, we show that waves are continuously launched from the vicinity of the black hole and propagate along the interface between the jet and the wind. At this interface, a steep gradient in velocity is present between the mildly relativistic wind and the highly relativistic jet. The interface is, therefore, a shear layer, and due to the shear, the waves generate roll-ups that alter the magnetic field configuration and the shear layer geometry. We then perform polarized radiation transfer calculations of our GRMHD simulation and find signatures of the waves in both total intensity and linear polarization, effectively lowering the fully resolved polarization fraction. The telltale polarization signatures of the waves could be observable by future very long baseline interferometric observations, e.g., the next-generation Event Horizon Telescope.

Astronomy & Astrophysics↗