Publisher’s Note: “Nonadiabatic electron response in the Hasegawa-Wakatani equations”[Phys. Plasmas 20, 082314 (2013)]
No abstract prepared.
SEARCH · Engineering Papers
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.
No abstract prepared.
The linear and nonlinear evolution of magnetic reconnection in collisionless high-temperature plasmas with a strong guide field is analyzed on the basis of a two-dimensional gyrofluid model. The linear growth rate of the reconnecting instability is compared to analytical calculations over the whole spectrum of linearly unstable wave numbers. In the strongly unstable regime (large Δ′), the nonlinear evolution of the reconnecting instability is found to undergo two distinctive acceleration phases separated by a stall phase in which the instantaneous growth rate decreases. The first acceleration phase is caused by the formation of strong electric fields close to the X-point due to ion gyration, while the second acceleration phase is driven by the development of an open Petschek-like configuration due to both ion and electron temperature effects. Furthermore, the maximum instantaneous growth rate is found to increase dramatically over its linear value for decreasing diffusion layers. This is a consequence of the fact that the peak instantaneous growth rate becomes weakly dependent on the microscopic plasma parameters if the diffusion region thickness is sufficiently smaller than the equilibrium magnetic field scale length. When this condition is satisfied, the peak reconnection rate asymptotes to a constant value.
The nonlinear dynamics of filamentation instability and magnetic field in a current-carrying plasma is investigated in the presence of quantum effects using the quantum hydrodynamic model. A new nonlinear partial differential equation is obtained for the spatiotemporal evolution of the magnetic field in the diffusion regime. This equation is solved by applying the Adomian decomposition method, and then the profiles of magnetic field and electron density are plotted. It is shown that the saturation time of filamentation instability increases and, consequently, the instability growth rate and the magnetic field amplitude decrease in the presence of quantum effects.
The obliquely propagating ion acoustic (IA) double-layers are investigated in a magnetized two population electron plasmas. The extended Korteweg–de Vries equation is derived by using the reductive perturbation technique. The effect of obliqueness (l{sub z}) and magnitude of the external magnetic field (ω{sub ci}), as well as the electron number density (β) on the double-layer profile, is studied, and then the ranges of parameters for which the double-layers exist are investigated in detail. We found that the combined effects of l{sub z}, ω{sub ci}, and β significantly modify the basic properties (viz. amplitude and width) of the IA double-layers.
A performance model is presented for magnetic nozzle plasmas driven by electron thermal expansion to investigate how the thrust coefficient and beam divergence efficiency scale with the incoming plasma flow and magnetic field geometry. Using a transformation from cylindrical to magnetic coordinates, an approximate analytical solution is derived to the axisymmetric two-fluid equations for a collisionless plasma flow along an applied magnetic field. This solution yields an expression for the half-width at half-maximum of the plasma density profile in the far-downstream region, from which simple scaling relations for the thrust coefficient and beam divergence efficiency are derived. It is found that the beam divergence efficiency is most sensitive to the density profile of the flow into the nozzle throat, with the highest efficiencies occurring for plasmas concentrated along the nozzle axis. Increasing the expansion ratio of the magnetic field leads to efficiency improvements that are more pronounced for incoming plasmas that are not concentrated along the axis. This implies that the additional magnet required to increase the expansion ratio may be worth the added complexity for plasma sources that exhibit poor confinement.
The dynamics of filaments or blobs in the scrape-off layer of magnetic fusion devices are studied by magnitude estimates of a comprehensive drift-interchange-Alfvén fluid model. The standard blob models are reproduced in the cold ion case. Even though usually neglected, in the scrape-off layer, the ion temperature can exceed the electron temperature by an order of magnitude. The ion pressure affects the dynamics of filaments amongst others by adding up to the interchange drive and the polarisation current. It is shown how both effects modify the scaling laws for filament velocity in dependence of its size. Simplifications for experimentally relevant limit regimes are given. These are the sheath dissipation, collisional, and electromagnetic regime.
A two-dimensional particle-in-cell/Monte Carlo collision model has been developed and used to study low electronegative magnetized hydrogen plasma. A configuration characterized by four electrodes is used: the left electrode is biased at V{sub l} = −100 V, the right electrode is grounded, while the upper and lower transversal electrodes are biased at an intermediate voltage V{sub ud} between 0 and −100 V. A constant and homogeneous magnetic field is applied parallel to the lateral (left/right) electrodes. It is shown that in the magnetized case, the bulk plasma potential is close to the transversal electrodes bias inducing then a reversed sheath in front of the right electrode. The potential drop within the reversed sheath is controlled by the transversal electrodes bias allowing extraction of negative ions with a significant reduction of co-extracted electron current. Furthermore, introducing plasma electrodes, between the transversal electrodes and the right electrode, biased with a voltage just above the plasma bulk potential, increases the negative ion extracted current and decreases significantly the co-extracted electron current. The physical mechanism on basis of this phenomenon has been discussed.
The nonlinear evolution of collisionless plasmas is typically a multi-scale process, where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic (MHD), Hall-MHD, two-fluid, hybrid kinetic, and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligible in the nonlinear regime. This study shows that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the collisionless regime.
New analytic force-free Vlasov-Maxwell equilibria for thin current sheets are presented. The magnetic flux densities are expressed in terms of Jacobian elliptic functions of one Cartesian spatial coordinate. The magnetic flux densities reduce to previously reported hyperbolic functions in one limit and sinusoidal functions in another limit of the modulus k. A much wider class of nonlinear force-free Vlasov-Maxwell equilibria open expanded possibilities for modeling of solar system, astrophysical and laboratory plasmas. Modified Maxwellian distribution functions are determined explicitly in terms of Jacobian elliptic functions. Conditions for double peaked distribution functions that could be unstable are developed.
A previous publication [R. V. Bravenec et al., Phys. Plasmas 18, 122505 (2011)] presented favorable comparisons of linear frequencies and nonlinear fluxes from the Eulerian gyrokinetic codes gyro[J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)] and gs2[W. Dorland et al., Phys. Rev. Lett. 85, 5579 (2000)]. The motivation was to verify the codes, i.e., demonstrate that they correctly solve the gyrokinetic-Maxwell equations. The premise was that it is highly unlikely for both codes to yield the same incorrect results. In this work, we add the Lagrangian particle-in-cell code gem[Y. Chen and S. Parker, J. Comput. Phys. 220, 839 (2007)] to the comparisons, not simply to add another code, but also to demonstrate that the codes' algorithms do not matter. We find good agreement of gem with gyro and gs2 for the plasma conditions considered earlier, thus establishing confidence that the codes are verified and that ongoing validation efforts for these plasma parameters are warranted.
In nature and many laboratory plasmas, a magnetic flux tube threaded by current or a flux rope has a footpoint at a boundary. The current driven kink mode is one of the fundamental ideal magnetohydrodynamic instabilities in plasmas. It has an instability threshold that has been found to strongly depend on boundary conditions (BCs). We provide a theoretical model to explain the transition of this threshold dependence between nonline tied and line tied boundary conditions. We evaluate model parameters using experimentally measured plasma data, explicitly verify several kink eigenfunctions, and validate the model predictions for boundary conditions BCs that span the range between NLT and LT BCs. Based on this model, one could estimate the kink threshold given knowledge of the displacement of a flux rope end, or conversely estimate flux rope end motion based on knowledge of it kink stability threshold.
In a previous contribution to this journal [H. P. Stormberg, J. Appl. Phys. 51(4), 1963 (1980)], Stormberg presented an analytical expression for the convolution of Lorentz and Levy line profiles, which models atomic radiative transitions in high pressure plasmas. Unfortunately, the derivations are flawed with errors and the final expression, while correct, is accompanied by misguiding comments about the meaning of the symbols used therein, in particular the “complex error function.” In this paper, we discuss the broadening mechanisms that give rise to Stormberg's model and present a correct derivation of his final result. We will also provide an alternative expression, based on the Faddeeva function, which has decisive computational advantages and emphasizes the real-valuedness of the result. The MATLAB/Octave scripts of our implementation have been made available on the publisher's website for future reference.
Lower hybrid (LH) ray propagation in toroidal plasma is shown to be controlled by combination of the azimuthal spectrum launched by the antenna, the poloidal variation of the magnetic field, and the scattering of the waves by the drift wave fluctuations. The width of the poloidal and radial radio frequency wave spectrum increases rapidly as the rays penetrate into higher density and scatter from the drift waves. The electron temperature gradient (ETG) spectrum is particularly effective in scattering the LH waves due to its comparable wavelengths and phase velocities. ETG turbulence is also driven by the radial gradient of the electron current profile giving rise to an anomalous viscosity spreading the LH driven plasma currents. The LH wave scattering is derived from a Fokker-Planck equation for the distribution of the ray trajectories with diffusivities derived from the drift wave fluctuations. The condition for chaotic diffusion for the rays is derived. The evolution of the poloidal and radial mode number spectrum of the lower hybrid waves are both on the antenna spectrum and the spectrum of the drift waves. Antennas launching higher poloidal mode number spectra drive off-axis current density profiles producing negative central shear [RS] plasmas with improved thermal confinement from ETG transport. Core plasma current drive requires antennas with low azimuthal mode spectra peaked at m = 0 azimuthal mode numbers.
Isochoric heating of solid-density matter up to a few tens of eV is of interest for investigating astrophysical or inertial fusion scenarios. Such ultra-fast heating can be achieved via the energy deposition of short-pulse laser generated electrons. Here, we report on experimental measurements of this process by means of time- and space-resolved optical interferometry. Our results are found in reasonable agreement with a simple numerical model of fast electron-induced heating.
Results from a 3D self-consistent Particle-In-Cell Monte Carlo Collisions (PIC MCC) model of a high power fusion-type negative ion source are presented for the first time. The model is used to calculate the plasma characteristics of the ITER prototype BATMAN ion source developed in Garching. Special emphasis is put on the production of negative ions on the plasma grid surface. The question of the relative roles of the impact of neutral hydrogen atoms and positive ions on the cesiated grid surface has attracted much attention recently and the 3D PIC MCC model is used to address this question. The results show that the production of negative ions by positive ion impact on the plasma grid is small with respect to the production by atomic hydrogen or deuterium bombardment (less than 10%)
This paper is aimed to contribute to the scientific discussions that have been triggered by the experimental observation of a quadratic relation between the kurtosis and skewness of turbulent fluctuations present in fusion plasmas and other nonlinear physical systems. In this paper, we offer a general statistical model which attributes the observed K=aS{sup 2}+b relation to the varying intermittency of the experimental signals. The model is a two random variable model constructed to catch the essential intermittent feature of the real signal. One of the variables is the amplitude of the underlying intermittent event (e.g., turbulent structure) while the other is connected to the intermittency level of the system. This simple model can attribute physical meaning to the a and b coefficients, as they characterize the spatio-temporal statistics of intermittent events. By constructing a particle-conserving Gaussian model for the underlying coherent structures, the experimentally measured a and b coefficients could be adequately reproduced.
Resistively guiding laser-generated fast electron beams in targets consisting of a resistive wire embedded in lower Z material should allow one to rapidly heat the wire to over 100 eV over a substantial distance without strongly heating the surrounding material. On the multi-ps timescale, this can drive hydrodynamic motion in the surrounding material. Thus, ultra-intense laser solid interactions have the potential as a controlled driver of radiation hydrodynamics in solid density material. In this paper, we assess the laser and target parameters needed to achieve such rapid and controlled heating of the embedded wire.
The counterpropagating shear experiment has previously demonstrated the ability to create regions of shock-driven shear, balanced symmetrically in pressure, and experiencing minimal net drift. This allows for the creation of a high-Mach-number high-energy-density shear environment. New data from the counterpropagating shear campaign is presented, and both hydrocode modeling and theoretical analysis in the context of a Reynolds-averaged-Navier-Stokes model suggest turbulent dissipation of energy from the supersonic flow bounding the layer is a significant driver in its expansion. A theoretical minimum shear flow Mach number threshold is suggested for substantial thermal-turbulence coupling.