Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PHASE VELOCITY”

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 235 records · Page 13

A neoclassical drift-magnetohydrodynamical fluid model of the interaction of a magnetic island chain with a resonant error-field in a high temperature tokamak plasma

A two-fluid, neoclassical theory of the interaction of a single magnetic island chain with a resonant error-field in a quasi-cylindrical, low-β, tokamak plasma is presented. The plasmas typically found in large hot tokamaks lie in the so-called weak neoclassical flow-damping regime in which the neoclassical ion stress tensor is not the dominant term in the ion parallel equation of motion. Nevertheless, flow-damping in such plasmas dominates ion perpendicular viscosity, and is largely responsible for determining the phase velocity of a freely rotating island chain (which is in the ion diamagnetic direction relative to the local E × B frame at the rational surface). The critical vacuum island width required to lock the island chain is mostly determined by the ion neoclassical poloidal flow damping rate at the rational surface. The stabilizing effect of the average field-line curvature, as well as the destabilizing effect of the perturbed bootstrap current, is the same for a freely rotating, a non-uniformly rotating, and a locked island chain. The destabilizing effect of the error-field averages to zero when the chain is rotating and only manifests itself when the chain locks. The perturbed ion polarization current has a small destabilizing effect on a freely rotating island chain, but a large destabilizing effect on both a non-uniformly rotating and a locked island chain. Finally, this behavior may account for the experimentally observed fact that locked island chains are much more unstable than corresponding freely rotating chains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Harnessing the relativistic Buneman instability for laser-ion acceleration in the transparency regime

Here we examine the relativistic Buneman instability in systems relevant to high-intensity laser-plasma interactions under conditions of relativistically-induced transparency, as this instability can generate large-amplitude electrostatic waves at low frequencies that are pertinent to ion dynamics in these systems. Ion flows are shown to significantly alter the range of unstable wave numbers and to increase the phase velocities of the unstable modes; we particularly highlight the relativistic effects from both the ion and electron (with transverse motion) populations. These findings are related to the mode structure seen in particle-in-cell simulation results of a short-pulse laser breaking through an initially opaque target with the onset of relativistic transparency. Additionally, driving mechanisms from free energy present in density and velocity gradients are shown to be capable of significantly enhancing the growth rates, and these instabilities furthermore extend the breadth of the unstable wave number range. We discuss how the transverse self-generated magnetic fields characteristic of short-pulse interactions can potentially constrain the unstable wave numbers in a non-trivial manner.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A detailed examination of laser-ion acceleration mechanisms in the relativistic transparency regime using tracers

Here we present a particle-in-cell study of linearly polarized laser-ion acceleration systems, in which we use both two-dimensional (2D) and three-dimensional (3D) simulations to characterize the ion acceleration mechanisms in targets which become transparent to the laser pulse during irradiation. First, we perform a target length scan to optimize the peak ion energies in both 2D and 3D, and the predictive capabilities of 2D simulations are discussed. Tracer analysis allows us to isolate the acceleration into stages of target normal sheath acceleration (TNSA), hole boring (HB), and break-out afterburner (BOA) acceleration, which vary in effectiveness based on the simulation parameters. The thinnest targets reveal that enhanced TNSA is responsible for accelerating the most energetic ions, whereas the thickest targets have ions undergoing successive phases of HB and TNSA (in 2D) or BOA and TNSA (in 3D); HB is not observed to be a dominant acceleration mechanism in the 3D simulations. It is in the intermediate optimal regime, both when the laser breaks through the target with appreciable amplitude and when there is enough plasma to form a sustained high density flow, that BOA is most effective and is responsible for the most energetic ions. Eliminating the transverse laser spot size effects by performing a plane wave simulation, we can isolate with greater confidence the underlying physics behind the ion dynamics we observe. Specifically, supplemented by wavelet and FFT analyses, we match the post-transparency BOA acceleration with a wave-particle resonance with a high-amplitude low-frequency electrostatic wave of increasing phase velocity, consistent with that predicted by the Buneman instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extracting the turbulent flow-field from beam emission spectroscopy images using velocimetry

The 2D turbulent E × B flow-field is inferred from density fluctuation images obtained with the beam emission spectroscopy diagnostic on DIII-D using the orthogonal dynamic programming velocimetry algorithm. A synthetic turbulence model is used to test the algorithm and optimize it for measuring zonal flows. Zonal flow measurements are found to require a signal-to-noise ratio above ~10 and a zonal flow wavelength longer than ~2 cm. Comparison between the velocimetry-estimated flow-field and the E × B flow-field using a nonlinear gyrokinetic GENE simulation finds that the flow-fields have identical spatial structure and differ only by the mean turbulence phase velocity, which is spatially uniform in this flux tube simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Angular streaking of betatron X-rays in a transverse density gradient laser-wakefield accelerator

In a plasma with a transverse density gradient, laser wavefront tilt develops gradually due to phase velocity differences in different plasma densities. The wavefront tilt leads to a parabolic trajectory of the plasma wakefield and hence the accelerated electron beam, which leads to an angular streaking of the emitted betatron radiation. In this way, the temporal evolution of the betatron X-ray spectra will be converted into angular “streak,” i.e., having a critical energy-angle correlation. An analytical model for the curved trajectory of a laser pulse in a transverse density gradient is presented. This gives the deflection angle of the electron beam and the betatron X-rays as a function of the plasma and laser parameters, and it was verified by particle-in-cell simulations. This angular streaking could be used as a single-shot diagnostic technique to reveal the temporal evolution of betatron X-ray spectra and hence the electron acceleration itself.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measuring heat flux from collective Thomson scattering with non-Maxwellian distribution functions

Here, heat flux was measured in coronal plasmas using collective Thomson scattering from electron-plasma waves. A laser-produced plasma from a planar aluminum target created a temperature gradient along the target normal. Thomson scattering probed electron-plasma waves in the direction of the temperature gradient with phase velocities relevant to heat flux. The heat-flux measurements were reduced from classical values inferred from the measured plasma conditions in regions with large temperature gradients and agreed with classical values for weak gradients. In regions where classical theory was invalid, the heat flux was determined by reproducing the measured Thomson-scattering spectra using electron distribution functions consistent with nonlocal thermal transport. Lastly, full-scale hydrodynamic simulations using flux-limited thermal transport ( FLASH ) and the multigroup nonlocal Schurtz, Nicolai, and Busquet (SNB) models both underestimated the heat flux at all locations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear response of a Hall magnetic drift wave for verification of Hall MHD algorithms

Numerical implementations of Hall magnetohydrodynamics (Hall MHD) can be challenging due to the nonlinear multidimensional nature of the Hall term. Here, a model problem is presented that couples the hydrodynamic motion of the plasma to Hall MHD evolution of the magnetic field. The Hall MHD equations are linearized about unperturbed solutions in both cylindrical and Cartesian coordinates in two dimensions. The magnetic field is assumed to lie in the ignorable direction, and the linear response about the unperturbed solution is considered. The resulting ordinary differential equation is used to numerically compute the eigenfunctions and eigenfrequencies of the mode. The resulting eigenfunctions do not make the local wave approximation but are instead global solutions that depend on the spatial dependence of the unperturbed Alfvén speed. Hall MHD simulations are then performed in the Ares multiphysics code and shown to agree with the predicted phase velocities of the wave, and the simulated solutions are shown to numerically converge to the semianalytic modes. By varying the background density of the plasma (and correspondingly, the ion inertial length), the importance of Hall physics can be varied. This allows the test problem to transition from the classical MHD limit to the extreme Hall MHD limit. Furthermore, this problem is a useful tool for the verification of Hall MHD routines implemented in various codes, and the robustness of a routine can be tested in regimes in which Hall physics is dominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Instability of nonlinear Trivelpiece-Gould waves. II. Weakly trapped particles

This paper discusses a novel parametric instability mechanism caused by particles that are weakly trapped in the potential wells of a nonlinear “pump” wave. The pump wave is unstable to the growth of daughter waves with longer wavelength and nearly the same phase velocity as the pump. This induces adjacent potential peaks in the wave to slowly approach one-another, receding from other pairs of peaks. Particles that are weakly trapped between approaching peaks, with kinetic energies just below the potential maxima, are heated by compression and escape the well, and then become retrapped on the other side of the approaching peaks, where they amplify the compression by pushing the peaks together. The mechanism applies to low-collisionality plasmas supporting waves with near-acoustic dispersion relations such as ion sound waves, magnetized Langmuir waves, or Alfven waves. The theory is compared to particle in cell simulations of TrivelpieceGould (TG) traveling waves, as well as to experiments on pure ion plasmas that observe parametric instability in TG standing waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Laser-plasma acceleration beyond wave breaking

Laser wakefield accelerators rely on the extremely high electric fields of nonlinear plasma waves to trap and accelerate electrons to relativistic energies over short distances. When driven strongly enough, plasma waves break, trapping a large population of the background electrons that support their motion. Aside from limiting the maximum electric field, this trapping can lead to accelerated electron bunches with large energy spreads. Here we introduce a novel regime of plasma wave excitation and wakefield acceleration that allows for arbitrarily high electric fields while avoiding the deleterious effects of unwanted trapping. Here, the regime, enabled by spatiotemporal shaping of laser pulses, exploits the property that nonlinear plasma waves with superluminal phase velocities cannot trap charged particles and are therefore immune to wave breaking.

36 MATERIALS SCIENCE↗

Filamentation and focusing of electron beams due to interactions with plasma waves

Results of numerical modeling of the interaction of an electron beam propagating across relativistic plasma waves indicate that electron beam filamentation and focusing may occur under certain conditions. The model is based on solving the relativistic equation of motion in three dimensions for the individual electrons in a tenuous Gaussian beam, as they pass through a relativistic plasma wave. Several electron beam and plasma wave parameters were varied, and the results are summarized. One of the results is that the spacing of the electron beam filaments correlates with the wavelength of the plasma wave. The electron beam filaments appear as vertical slabs after the beam exits the plasma waves. The electron beam also compresses to a focus after it exits the plasma, and the focal distance depends on several parameters including the electron beam energy and phase velocity of the relativistic plasma wave. It is suggested that these focusing and filamentation phenomena may be the basis for diagnostics schemes for laser plasma interactions. The parameters used in the model were electron beam energies in the 5–50 keV range and plasma wave properties typical for the beat-wave produced by CO2 lasers, which correspond to the facilities available in our laboratory. The limitations of these results to lower energy density beam and plasma regimes and to higher energy density regimes will be discussed.

47 OTHER INSTRUMENTATION↗

A multi-sheath model for highly nonlinear plasma wakefields

An improved description for nonlinear plasma wakefields with phase velocities near the speed of light is presented and compared against fully kinetic particle-in-cell simulations. These wakefields are excited by intense particle beams or lasers pushing plasma electrons radially outward, creating an ion bubble surrounded by a sheath of electrons characterized by the source term S≡−1enp(ρ−Jz/c), where ρ and Jz are the charge and axial current densities, respectively. Previously, the sheath source term was described phenomenologically with a positive-definite function, resulting in a positive definite wake potential. In reality, the wake potential is negative at the rear of the ion column which is important for self-injection and accurate beam loading models. To account for this, we introduce a multi-sheath model in which the source term, S, of the plasma wake can be negative in regions outside the ion bubble. Using this model, we obtain a new expression for the wake potential and a modified differential equation for the bubble radius. Numerical results obtained from these equations are validated against particle-in-cell simulations for unloaded and loaded wakes. The new model provides accurate predictions of the shape and duration of trailing bunch current profiles that flatten plasma wakefields. It is also used to design a trailing bunch for a desired longitudinally varying loaded wakefield. We present beam loading results for laser wakefields and discuss how the model can be improved for laser drivers in future work. Finally, we discuss differences between the predictions of the multi- and single-sheath models for beam loading.

Dalichaouch, T. N. (ORCID:0000000247350150)↗

A statistical study of three-second foreshock ULF waves observed by the Magnetospheric Multiscale mission

We perform a statistical study of 3-s ultra-low frequency (ULF) waves using Magnetospheric Multiscale observations in the Earth's foreshock region. The average phase velocity in the plasma rest frame is determined to be anti-sunward, and the intrinsic polarization is right-handed. We further examine the linear instability conditions based on the drift-bi-Maxwellian distribution functions according to the observed plasma conditions. The resulting instability is a solution to the common dispersion equation of the ion/ion right-hand non-resonant and left-hand resonant instabilities. The predicted wave propagation is also predominantly anti-sunward. The cyclotron resonant conditions of the solar wind and backstreaming beam ions are evaluated, and we find that, in some cases, the anti-sunward propagating waves can be resonant with beam ions, which was overlooked in previous studies. The study suggests that the dispersion equation provides the 3-s ULF waves a fundamental explanation that unifies a rich variety of resonant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shivalik Plasma Device-I, a glow discharge device to study the collective dynamics of dusty plasma

Using a uniquely configured glow discharge-based Shivalik Plasma Device-I, we demonstrate a variety of collective phenomena in dusty plasma away from the glow discharge region. The cylindrical glass device produces plasma using parallel disc-shaped electrodes with a smaller anode size than the cathode. The dust microparticles are initially sprinkled over the grounded cathode. These particles acquire a significant negative charge upon plasma formation, resulting in their levitation due to the balance between the Coulomb force and gravity. The new device supports the levitation of a big-sized (10 × 8 × 5 cm 3 ) three-dimensional dust cloud over the glass surface. It contrasts the dusty plasma formations in-between electrodes reported earlier. As the discharge voltage varies from high to low, the dust cloud travels from over the glass surface to between the electrodes. A complex interplay of dust void over the cathode, a sharp density gradient, and gravity lead to self-excitation of collective dust phenomena. It includes dust density waves (phase velocity, v ph ~4 cm/s), dust cloud oscillation (frequency, f = 5 Hz), sheared flow (flow velocity, v f ~ 1 cm/s), and multiple-sized dust vortices. These dust vortices provided an excellent platform for studying turbulent mixing phenomena. The power spectrum analysis agreed with two-dimensional Kolmogorov power-law scaling. This is an ideal dusty plasma apparatus where we can create or move the dust cloud to a location of choice from the glass surface to in-between the electrodes and excite one among many collective dust dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ELM-free enhanced D α H-mode with near zero NBI torque injection in DIII-D tokamak

Enhanced D α H-mode (EDA H-mode), an ELM-free H-mode regime, and the concomitant quasi-coherent mode (QCM) are explored in neutral beam heated, lower single null plasmas with near zero external torque injection. This regime exhibits good energy confinement (H 98y2 ~ 1) with β N ~ 2, high density, regime access at low input power, and no ELMs. New time–resolved measurements of electron and ion density, temperature, plasma rotation, and radial electric field as well as turbulence, flows, thermal and particle transport, and linear instability calculations examining the dynamics of the QCM are presented. Measurements using multiple fluctuation diagnostics reveal the QCM to be a separatrix spanning mode, peaking just inside the separatrix, existing in a wide range of k ⊥ ρ s ~ 0.1–1.2 with multiple harmonics, and propagating with a very small phase velocity in the plasma frame, where k ⊥ is the binormal wavenumber and ρ s is the ion sound radius. Linear gyrokinetic simulations of this EDA H-mode discharge with CGYRO indicate that the trapped electron mode (TEM) and electron temperature gradient mode are the dominant instabilities in the radial region where QCM is also unstable. Qualitative analysis indicates that the properties of TEM are consistent with the experimentally observed characteristics of the QCM. These similarities suggest that the QCM is consistent with a TEM instability existing in the edge region of the EDA H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine-Learning Analysis of Moisture Carryover in Boiling Water Reactors

Moisture carryover (MCO) is modeled in the General Electric Type-4 boiling water reactor (BWR) using machine-learning methods and data from operating plants. Understanding MCO and the conditions that give rise to an elevated value is important since excessive MCO can damage critical turbine components, can result in elevated dose levels to on-site personnel, and can interfere with late-cycle power management. The analysis of MCO takes into account simplifying reactor symmetries and important geometric dependencies. The plant data are taken from several reactors and were collected over multiple years and multiple fuel cycles. A brief description of the origin of MCO in U.S. BWR plants is given. A machine-learning model is constructed from the data using applicable algorithms and data-reduction techniques. Matching model complexity with available data is one of the more challenging machine-learning tasks. Too many features and too little data will lead to overfitting. The data for each fuel cycle included over 6876 original features, 9 for each fuel bundle. Two approaches are used to reduce the data set into a manageable number of features. The first was an engineering analysis that resulted in the selection of steam quality Q and steam liquid phase velocity V L as the main features driving MCO. Using a Q and a V L for each fuel bundle gives 1528 Q and a V L feature describing the reactor behavior. An analysis of different functional forms of these two variables led to the actual inputs to the neural network model. The second approach involved the use of statistical techniques such as Pearson’s correlation and k-means analysis. The identified groupings of bundles behaved similarly. Treating each grouping as a single feature further reduced the input variable set to a manageable number. Here, a model selection criterion is proposed, and results are presented along with a discussion of related issues.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Characterization of density fluctuations during the search for an I-mode regime on the DIII-D tokamak

The I-mode regime, routinely observed on the Alcator C-Mod tokamak, is characterized by an edge energy transport barrier without an accompanying particle barrier and with broadband instabilities, known as weakly coherent modes (WCM), believed to regulate particle transport at the edge. Recent experiments on the DIII-D tokamak exhibit I-mode characteristics in various physical quantities. These DIII-D plasmas evolve over long periods, lasting several energy confinement times, during which the edge electron temperature slowly evolves towards an H-mode-like profile, while maintaining a typical L-mode edge density profile. During these periods, referred to as I-mode phases, the radial electric field at the edge also gradually reaches values typically observed in H-mode. Density fluctuations measured with the phase contrast imaging diagnostic during I-mode phases exhibit three features typically observed in H-mode on DIII-D, although they develop progressively with time and without a sharp transition: the intensity of the fluctuations is reduced; the frequency spectrum is broadened and becomes non-monotonic; two dimensional space-time spectra appear to approach those in H-mode, showing phase velocities of density fluctuations at the edge increasing to about 10 km s -1 . However, in DIII-D there is no clear evidence of the WCM. Preliminary linear gyro-kinetic simulations are performed in the pedestal region with the GS2 code and its recently upgraded model collision operator that conserves particles, energy and momentum. The increased bootstrap current and flow shear generated by the temperature pedestal are shown to decrease growth rates, thus possibly generating a feedback mechanism that progressively stabilizes fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic particle simulation of the effect of toroidicity on linear mode conversion and absorption of lower hybrid waves

Effects of toroidicity on linear mode conversion and absorption of lower hybrid (LH) waves in fusion plasmas have been studied using electromagnetic particle simulation. The simulation confirms that the toroidicity induces an upshift of parallel refractive index when LH waves propagate from the tokamak edge toward the core, which affects the radial position for the mode conversion between slow and fast LH waves. Furthermore, moving LH antenna launch position from low field side toward high field side leads to a larger upshift of the parallel refractive index, which helps the slow LH wave penetration into the tokamak core. The broadening of the poloidal spectrum of the wave-packet due to wave diffraction is also verified in the simulation. Both the upshift and broadening effects of the parallel spectrum of the wave-packet modify the parallel phase velocity and thus the linear absorption of LH waves by electron Landau resonance.

Physics↗

Evolution of plasma wakes in density up- and down-ramps

The time evolution of plasma wakes in density up- and down-ramps is examined through theory and particle-in-cell simulations. Motivated by observation of the reversal of a linear plasma wake in a plasma density upramp in a recent experiment we have examined the behaviour of wakes in plasma ramps that always accompany any plasma source used for plasma-based acceleration. In the up-ramp case it is found that, after the passage of the drive pulse, the wavnumber/wavelength of the wake starts to decrease/increase with time until it eventually tends to zero/infinity, then the wake reverses its propagation direction and the wavenunber/wavelength of the wake begins to increase/shrink. Here, the evolutions of the wavenumber and the phase velocity of the wake as functions of time are shown to be significantly different in the up-ramp and the down-ramp cases. In the latter case the wavenumber of the wake at a particular position in the ramp increases until the wake is eventually damped. It is also shown that the waveform of the wake at a particular time after being excited can be precisely controlled by tuning the initial plasma density profile, which may enable a new type of plasma-based ultrafast optics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗