Search NASA⌕ Search

SEARCH · Search NASA

Results for “Anisotropic plasmas”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2021

Modification of favorable average curvature effect by changing parallel sound wave behavior in tokamak plasmas

Abstract The favorable average curvature effect, also known as the GGJ effect (Glasser et al 1975 Phys. Fluids 18 875), is intrinsically associated with parallel sound wave propagation in a tokamak plasma. This work investigates how the GGJ effect is modified by changing the parallel sound wave behavior. Two physics models beyond the standard single fluid theory, i.e. an anisotropic thermal transport model and a parallel sound wave damping model, are employed to change parallel sound waves in a toroidal plasma, and the consequence on the GGJ effect is demonstrated for two important classes of problems, i.e. the resistive plasma response to the applied resonant magnetic perturbation and the stability of the tearing mode. Toroidal modeling reveals that the GGJ effect is significantly altered by both of the aforementioned physics effects. Compared to the thermal transport physics, which completely removes the GGJ effect, the sound wave damping effect only offers partial mitigation. The differences between these two models are further illustrated in terms of the radial structure of the shielding current and the eigenfunction of the tearing instability. In particular, a fundamental reason for complete suppression of the GGJ effect by the thermal transport is identified as an extra toroidal coupling of the poloidal harmonics.

Physics↗
NASA NTRS1976

Filamentation instability of a large amplitude whistler wave

Experiments performed in a pulsed afterglow plasma column with specified parameters revealed a filamentation instability of a large-amplitude whistler wave launched from antennas which produce a diverging energy flow in the linear regime. The difficult problem of diagnosing local density perturbations in a magnetized plasma in the presence of large-amplitude RF signals and nonuniform anisotropic electron distributions is discussed. Since Langmuir probes are too unreliable under these conditions, a microwave probe is developed based on the principle of the cavity shift method. The temporal and spatial evolution of the duct formation and wave propagation are shown. The role of the observed electron heating in the filamentation process is examined. The observed self-focusing process appears to be the result of the reinforcing interaction between the wave-induced density depression and the density-induced wave refraction. Interesting applications are noted.

Stenzel, R. L.↗
NASA NTRS1994

Gamma-ray and millimeter-wave emissions from the 1991 June X-class solar flares

We have studied the spectacular 1991 June X-class flares using gamma-ray data from the Charged Particle Detectors (CPDs) of the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory (CGRO) and 80 GHz millimeter data from Nobeyama, Japan. The CPDs were the only CGRO instrument that did not saturate during the extremely intense 1991 June 4 flare. We have shown that for this flare the CPDs respond to MeV photons, most of which are due to bremsstrahlung produced by relativistic electrons at the Sun. We have further shown that the gamma-ray and millimeter observations agree numerically if the 80 GHz radiation is gyrosynchrotron radiation produced by trapped electrons and the gamma rays are thick-target bremsstrahlung due to electrons precipitating out of the trap. The requirement that the trapping time obtained from the numerical comparison be consistent with the observed time profiles implies a magnetic field between about 200 and 300 G and an electron spectral index between about 3 to 5. By comparing the CPD observations with both the 80 GHz data and nuclear line data from the Energetic Gamma Ray Experiment Telescope (EGRET) and the Oriented Scintillation Spectroscopy Experiment (OSSE) on CGRO for the flares of June 4, 6, 9, and 11, we found that the ratio of the CPD counts to both the millimeter flux densities and the nuclear line fluences decreases with decreasing flare heliocentric angle. All of these flares were produced in the same active region. We interpreted this result in terms of a loop model in which the gyrosynchrotron emission is produced in the coronal portion of the loop where the electrons are kept isotropic by pitch angle scattering due to plasma turbulence, while the bremsstrahlung is produced by precipitating electrons that interact anisotropically. We found that the trapping time in the coronal portion is time dependent, reaching a minimum of about 10 s at the peak of the CPD count rate. We suggested the damping of the turbulence as a possible reason for the variation of the trapping time. turbulence as a possible reason for the variation of the trapping time.

Ramaty, Reuven↗
DOE OSTI2024

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Hybrid Model of Inhomogeneous Solar Wind Plasma Heating by Alfven Wave Spectrum: Parametric Studies

Observations of the solar wind plasma at 0.3 AU and beyond show that a turbulent spectrum of magnetic fluctuations is present. Remote sensing observations of the corona indicate that heavy ions are hotter than protons and their temperature is anisotropic (T(sub perpindicular / T(sub parallel) >> 1). We study the heating and the acceleration of multi-ion plasma in the solar wind by a turbulent spectrum of Alfvenic fluctuations using a 2-D hybrid numerical model. In the hybrid model the protons and heavy ions are treated kinetically as particles, while the electrons are included as neutralizing background fluid. This is the first two-dimensional hybrid parametric study of the solar wind plasma that includes an input turbulent wave spectrum guided by observation with inhomogeneous background density. We also investigate the effects of He++ ion beams in the inhomogeneous background plasma density on the heating of the solar wind plasma. The 2-D hybrid model treats parallel and oblique waves, together with cross-field inhomogeneity, self-consistently. We investigate the parametric dependence of the perpendicular heating, and the temperature anisotropy in the H+-He++ solar wind plasma. It was found that the scaling of the magnetic fluctuations power spectrum steepens in the higher-density regions, and the heating is channeled to these regions from the surrounding lower-density plasma due to wave refraction. The model parameters are applicable to the expected solar wind conditions at about 10 solar radii.

Ofman, L.↗
DOE OSTI2021

Erosion of focus rings in capacitively coupled plasma etching reactors

In plasma etching reactors, the structure surrounding the wafer, often called a focus ring (FR), plays an important role in maintaining uniform fluxes of reactants across the wafer. The FR is typically made of dielectric materials. During the plasma etching process, the sheath that forms over the wafer to accelerate ions anisotropically into the surface extends over the FR. Electrical charging of the FR modifies the sheath relative to that over the wafer. On the one hand, one wants the sheath to be uniform across the wafer-FR boundary to enable uniform fluxes to the edge of the wafer. On the other hand, maintaining a high voltage sheath over the FR will erode the FR, which is undesirable as the FR is a consumable component that must be periodically replaced in high volume manufacturing. In this work, we computationally investigated the consequences of dielectric constant ε of the FR materials on erosion of the FR. The series capacitance of the FR and its underlying structure is typically smaller than that of the wafer and its underlying structure. As a result, the FR charges quickly relative to the wafer, which then reduces the voltage across the sheath on top of the FR. The ion energy and angular distributions (IEADs) striking the FR are, therefore, generally lower in energy with a broader angular distribution. With ε = 2, the ion energies striking the middle to the outer edge of the FR are 30–180 eV, whereas for ε = 100, the ion energies are 120–380 eV. At the transition between the wafer and the FR, there is a skew in the IEAD as large as 15° that results from the difference in sheath thickness above the wafer and the FR. This skew and the erosion rate across the FR are functions of the dielectric constant of the FR material. With low ε, the FR charges quickly, less plasma is produced above the FR, and there is less voltage across the sheath that results in less FR erosion. Increasing ε of the FR produces a higher sheath voltage as well as higher ion fluxes over the FR, which increases erosion, while the skew at the edge of the wafer is less severe. Further, the material of the subsurface portion of the FR, which dominates its capacitance, is an important consideration in the design of the substrate assembly.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles in $\mathrm {PbPb}$ collisions at $\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 2.76$ and 5.02$\,\text {TeV}$

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective anisotropic flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in $\mathrm {PbPb}$ collisions at nucleon-nucleon center-of-mass energies $\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 2.76$ and 5.02$\,\text {TeV}$ with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Advanced low blaze angle x-ray gratings via nanoimprint replication and plasma etch

We developed a new method of making ultra-low blaze angle diffraction gratings for x-ray applications. The method is based on reduction of the blaze angle of a master grating by nanoimprint replication followed by a plasma etch. A master blazed grating with a relatively large blaze angle is fabricated by anisotropic wet etching of a Si single crystal substrate. The surface of the master grating is replicated by a polymer material on top of a quartz substrate by nanoimprinting. Then a second nanoimprinting is performed using the 1st replica as a mold to replicate the saw-tooth surface into a resist layer on top of a Si grating substrate. A reactive ion etch is used to transfer the grating grooves into the Si substrate. The plasma etch provides reduction of the groove depth by a factor defined by the ratio of the etch rates for the resist and Si. We demonstrate reduction of the blaze angle of a master grating by a factor of 5 during fabrication of a 200 lines/mm blazed grating with a blaze angle of 0.2°. We investigated the quality and performance of the fabricated low blaze angle gratings and evaluate process accuracy and reproducibility. The new blaze angle reduction method preserves the planarity of the optical surface of the grating substrate and at the same time provides improvement in the grating groove quality during the reduction process.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Multipoint Observations of Quasiperiodic Emission Intensification and Effects on Energetic Electron Precipitation

The two Van Allen Probes simultaneously recorded a coherently modulated quasiperiodic (QP) emission that persisted for 3 h. The magnetic field pulsation at the locations of the two satellites showed a substantial difference, and their frequencies were close to but did not exactly match the repetition frequency of QP emissions for most of the time, suggesting that those coherent QP emissions probably originated from a common source, which then propagated over a broad area in the magnetosphere. The QP emissions were amplified by local anisotropic electron distributions, and their large-scale amplitudes were modulated by the plasma density. A novel observation of this event is that chorus waves at frequencies above QP emissions exhibit a strong correlation with QP emissions. Those chorus waves intensified when the QP emissions reach their peak frequency. This indicates that embryonic QP emissions may be critical for its own intensification as well as chorus waves under certain circumstances. A Low-Earth-Orbit POES satellite observed enhanced energetic electron precipitation in conjunction with the Van Allen Probes, providing direct evidence that QP emissions precipitate energetic electrons into the atmosphere. Finally, this scenario is quantitatively confirmed by our quasilinear diffusion simulation results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗
NASA NTRS1976

Excitation of the ordinary electromagnetic mode in low-beta plasmas

The excitation of the ordinary mode propagating perpendicular to an external magnetic field in a plasma is discussed for the case where the ratio of the kinetic energy associated with electron motion parallel to the external magnetic field to the energy density of the magnetic field is low. The analysis is performed for plasma in which in addition to the thermal electrons there is a hot anisotropic electron component. Relativistic effects are taken into account in the stability analysis. Both a ring distribution and a loss cone distribution are used to model the hot electron species. For the ring distribution typical growth rates are larger than that for the case of the loss cone distribution. Excitation of the ordinary mode is found at numerous harmonics of the cyclotron frequency for the ring distribution, however, the loss cone distribution is found to give instability primarily at the fundamental harmonic.

Freund, H. P.↗
DOE OSTI2025

Imaging nuclear shape through anisotropic and radial flow in high-energy heavy-ion collisions

Most atomic nuclei exhibit ellipsoidal shapes characterized by quadrupole deformation β 2 and triaxiality γ, and sometimes even a pear-like octupole deformation β 3 . The STAR experiment introduced a new ‘imaging-by-smashing’ technique to image the nuclear global shape by colliding nuclei at ultra-relativistic speeds and analyzing outgoing debris. Features of nuclear shape manifest in collective observables like anisotropic flow vn and radial flow via mean transverse momentum [p T ]. We present new measurements of the variances of v n (n = 2, 3, and 4) and [p T ], and the covariance of $v$$^{2}_{n}$ with [p T ], in collisions of highly deformed 238 U and nearly spherical 197 Au. Ratios of these observables between the two systems effectively suppress common final-state effects, isolating the strong impact of uranium’s deformation. By comparing results with state-of-the-art hydrodynamic model calculations, we extract β 2U and γ U values consistent with those deduced from low-energy nuclear structure measurements. Measurements of v 3 and its correlation with [p T ] also provide the first experimental suggestion of a possible octupole deformation for 238 U. These findings provide significant support for using high-energy collisions to explore nuclear shapes on femtosecond timescales, with implications for both nuclear structure and quark-gluon plasma studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Modelling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

ABSTRACT Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modelling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron–ion or electron–positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann–Vlasov equation. This new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes, not present in previous formulations of two-fluid magnetohydrodynamics. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. We keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The out-of-equilibrium corrections take the form of a collisional 14-moment closure previously described in the context of viscous single fluids. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines

We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Using Direct Laboratory Measurements of Electron Temperature Anisotropy to Identify the Heating Mechanism in Electron-Only Guide Field Magnetic Reconnection

Anisotropic electron heating during electron-only magnetic reconnection with a large guide magnetic field is directly measured in a laboratory plasma through in situ measurements of electron velocity distribution functions. Electron heating preferentially parallel to the magnetic field is localized to one separatrix, and anisotropies of 1.5 are measured. The mechanism for electron energization is identified as the parallel reconnection electric field because of the anisotropic nature of the heating and spatial localization. These characteristics are reproduced in a 2D particle-in-cell simulation and are also consistent with numerous magnetosheath observations. A measured increase in the perpendicular temperature along both separatrices is not reproduced by our 2D simulations. Here, this work has implications for energy partition studies in magnetosheath and laboratory reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Anisotropic shock jump conditions: Theory and observations

The MHD Rankine-Hugoniot (RH) relations for shock waves in a collisionless plasma with bi-Maxwellian distribution functions are considered. While by introducing the pressure anisotropy parameter xi in the RH relations, the number of unknowns -- B, V, n, p and xi (a total of 9) -- becomes one more than the total number of the conservation equations, it is possible to use the observed quantities on both sides of the shock to study the anisotropy changes across the shock. A simple relation for the anisotropy change across the shock is derived as a function of the ratio of magnetic fields m(= B'/B), the shock normal angle theta(sub Bn) and the plasma beta and beta' (primes are downstream values). Since m and theta(sub Bn) can be determined accurately in observation, the reliability of the anisotropy change deduced is mostly dependent on the accuracy of the measurements beta and beta'. We have applied the results to six low-beta quasi-perpendicular (Q perpendicular) laminar bow shock crossings with temperature anisotropy measured in the magnetosheath. In the six test cases, it is found that the predicted pressure anisotropies agree well with those observed in the magnetosheath.

Chao, J. K.↗
NASA NTRS2008

Mixed Mode Fracture of Plasma Sprayed Thermal Barrier Coatings: Effects of Anisotropy and Heterogeneity

The combined mode I-mode II fracture behavior of anisotropic ZrO2-8wt%Y2O3 thermal barrier coatings was determined in asymmetric flexure loading at both ambient and elevated temperatures. A fracture envelope of KI versus KII was determined for the coating material at ambient and elevated temperatures. Propagation angles of fracture as a function of KI/KII were also determined. The mixed-mode fracture behavior of the microsplat coating material was modeled using Finite Element approach to account for anisotropy and micro cracked structures, and predicted in terms of fracture envelope and propagation angle using mixed-mode fracture theories.

Zhu, Dongming↗