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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simultaneous Observation of Ion-scale Wave Packets with Opposite Polarizations and Their Implications on the Generation Region in the Inner Heliosphere

This paper reports a dispersion analysis of two wave packets simultaneously observed near the local proton gyrofrequency by the Parker Solar Probe. The observed wave event exhibits clear two-banded wave packets both propagating along the magnetic field, characterized by left-handed (L-mode) and right-handed (R-mode) polarizations simultaneously. By incorporating the Doppler shift effect into a linear dispersion analysis, we find two possible scenarios that explain these simultaneous opposite polarizations: (1) Two inherently L-mode waves in the plasma frame, propagate parallel and antiparallel to the solar wind velocity, with similar wave frequencies and wave numbers. The polarization of the antiparallel propagating wave reverses as it moves sunward in the plasma frame while still comoving with the solar wind in the stationary frame. This reversal manifests the polarization of the wave as an R-mode in the spacecraft frame. (2) Simultaneous L-mode and R-mode waves propagate parallel to the solar wind velocity, with different wave frequencies and wave numbers. Concurrent proton observations during the wave event reveal a dominant anisotropic ($T$⟂/$T$ ∥ > 1) core distribution with a drifting beam population. Estimation of the linear growth rate for both L-mode and R-mode waves suggests that both scenarios are plausible, indicating that the observation is near the wave-generation region. We explore the potential impact of these simultaneous waves on solar wind heating and scattering effects, hypothesizing that such waves might enhance efficiency compared to waves with a single wave packet, contingent upon the statistical significance of such waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Weak magnetohydrodynamic turbulence theory revisited

Two recent papers, P. H. Yoon and G. Choe, Phys. Plasmas 28, 082306 (2021) and Yoon et al., Phys. Plasmas 29, 112303 (2022), utilized in the derivation of the kinetic equation for the intensity of turbulent fluctuations the assumption that the wave spectra are isotropic, that is, the ensemble-averaged magnetic field tensorial fluctuation intensity is given by the isotropic diagonal form, ⟨δBiδBj⟩k=⟨δB2⟩kδij. However, it is more appropriate to describe the incompressible magnetohydrodynamic turbulence involving shear Alfvénic waves by modeling the turbulence spectrum as being anisotropic. That is, the tensorial fluctuation intensity should be different in diagonal elements across and along the direction of the wave vector, ⟨δBiδBj⟩k=12 ⟨δB⊥2⟩k(δij−kikj/k2)+⟨δB∥2⟩k(kikj/k2). In the present paper, we thus reformulate the weak magnetohydrodynamic turbulence theory under the assumption of anisotropy and work out the form of nonlinear wave kinetic equation.

Physics↗
DOE OSTI2022

Dynamics of pectic homogalacturonan in cellular morphogenesis and adhesion, wall integrity sensing and plant development

Homogalacturonan (HG) is the most abundant pectin subtype in plant cell walls. Although it is a linear homopolymer, its modification states allow for complex molecular encoding. HG metabolism affects its structure, chemical properties, mobility and binding capacity, allowing it to interact dynamically with other polymers during wall assembly and remodelling and to facilitate anisotropic cell growth, cell adhesion and separation, and organ morphogenesis. HGs have also recently been found to function as signalling molecules that transmit information about wall integrity to the cell. Here we highlight recent advances in our understanding of the dual functions of HG as a dynamic structural component of the cell wall and an initiator of intrinsic and environmental signalling. Here, we also predict how HG might interconnect the cell wall, plasma membrane and intracellular components with transcriptional networks to regulate plant growth and development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Nature of Non-Adiabatic Electron-Ion Forces in Liquid Metals

We report that an accurate description of electron-ion interactions in materials is crucial for our understanding of their equilibrium and non-equilibrium properties. Here, we assess the properties of frictional forces experienced by ions in non-crystalline metallic systems, including liquid metals and warm dense plasmas, that arise from electronic excitations driven by the nuclear motion due to the presence of a continuum of low-lying electronic states. To this end, we perform detailed ab-initio calculations of the full friction tensor that characterizes the set of friction forces. The nonadiabatic electron-ion interactions introduce hydrodynamic couplings between the ionic degrees of freedom, which are sizeable between nearest neighbors. The friction tensor is generally inhomogeneous, anisotropic and non-diagonal, especially at lower densities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Solving the problem of overdetermination of quasisymmetric equilibrium solutions by near-axis expansions. I. Generalized force balance

It is well known that the process of construction of quasisymmetric magnetic fields in magnetostatic equilibrium with isotropic pressure suffers from the problem of overdetermination. This has led to the widespread belief that global quasisymmetric solutions are likely not to exist. Here we develop a general near-axis expansion procedure that does not rely on the assumption of magnetostatic equilibria with isotropic pressure. We then demonstrate that in equilibria with anisotropic pressure, it is possible to circumvent the problem of overdetermination and carry out the power-series solutions to higher order. This suggests, contrary to current belief, that the existence of globally quasisymmetric fields is likely if one relaxes the assumption of magnetostatic equilibria with isotropic pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Bayesian calibration of viscous anisotropic hydrodynamic simulations of heavy-ion collisions

Owing to large pressure gradients at early times, standard hydrodynamic model simulations of relativistic heavy-ion collisions do not become reliable until O(1) fm/c after the collision. To address this one often introduces a prehydrodynamic stage that models the early evolution microscopically, typically as a conformal, weakly interacting gas. In such an approach the transition from the prehydrodynamic to the hydrodynamic stage is discontinuous, introducing considerable theoretical model ambiguity. Alternatively, fluids with large anisotropic pressure gradients can be handled macroscopically using the recently developed viscous anisotropic hydrodynamics (VAH). In high-energy heavy-ion collisions VAH is applicable already at very early times, and at later times transitions smoothly into conventional second-order viscous hydrodynamics. Here we present a Bayesian calibration of the VAH model with experimental data for Pb-Pb collisions at the LHC at $\sqrt{^{S}NN}$ = 2.76 A TeV. We find that the VAH model has the unique capability of constraining the specific viscosities of the quark-gluon plasma at higher temperatures than other previously used models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Electrical, thermal, and mechanical properties of spatially tailored fiber orientations in 3D printed carbon-carbon composites for EFAS/SPS

Using C-C synthesized at INL from continuous fiber 3D printed preforms, tooling was fabricated, tested, and compared to graphite. The anisotropic properties of the C-C material were leveraged to make stronger and more energy efficient tooling. Future work involves tailoring the fiber orientations to create targeted heating zones and thermal gradients.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Linear gyrokinetic simulations of toroidal Alfvén eigenmodes in the Mega-Amp Spherical Tokamak

Linear gyrokinetic (GK) simulations using the Gyrokinetic Toroidal Code (GTC) have been performed to investigate Toroidicity-driven Alfvén Eigenmodes (TAEs) driven by the neutral beam injection (NBI) induced fast ions in the Mega-Amp Spherical Tokamak (MAST) to identify the non-perturbative and kinetic effects of thermal plasma. A specific TAE in MAST discharge 26887, with an on-axis NBI power of approximately 1.5 MW and plasma current around 800 kA, exhibited frequency chirping, and the tangential soft x-ray camera array resolved the radial mode structure peaked near |q|=1.5. Various excitation methods were used in the GTC linear simulations, illustrating this code's capability to realistically represent the mechanisms and behaviors of fast ion-driven TAEs in spherical tokamaks. The radial structures from these GK simulations closely match measurements and calculations performed using the NOVA ideal MHD code, though with the frequencies approximately 10 kHz lower, likely due to various kinetic and non-perturbative effects. The simulations measured the damping rates due to continuum damping, radiative damping, and ion Landau damping, revealing that ion Landau damping has the most significant contribution to the total damping rate of the TAE. A comparison of growth rates of TAEs excited by fast ion Maxwellian and slowing-down distributions shows that the TAEs excited by a fast ion anisotropic pitch distribution (as part of the slowing-down distributions) are more unstable compared to those excited by a Maxwellian distribution with an equivalent fast ion beta. This shows that the use of fast ion anisotropy alters the number of fast ions to be in shear Alfvén resonance, and hence, it can greatly affect the stability of TAEs. These tests can be performed with the GTC but impossible with ideal MHD simulations, highlighting the necessity of kinetic simulations such as the GTC for a precise prediction of the TAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

3D orthorhombic earth model effects on seismic source characterization

Most earth materials are anisotropic with regard to seismic wave-speeds, especially materials such as shales, or where oriented fractures are present. However, the base assumption for many numerical simulations is to treat earth materials as isotropic media. This is done for simplicity, the apparent weakness of anisotropy in the far field, and the lack of well-characterized anisotropic material properties for input into numerical simulations. One approach for addressing the higher complexity of actual geologic regions is to model the material as an orthorhombic medium. We have developed an explicit time-domain, finite-difference (FD) algorithm for simulating three-dimensional (3D) elastic wave propagation in a heterogeneous orthorhombic medium. The objective of this research is to investigate the errors and biases that result from modeling a non-isotropic medium as an isotropic medium. This is done by computing “observed data” by using synthetic, anisotropic simulations with the assumption of an orthorhombic, anisotropic earth model. Green’s functions for an assumed isotropic earth model are computed and then used an inversion designed to estimate moment tensors with the “observed” data. One specific area of interest is how shear waves, which are introduced in an anisotropic model even for an isotropic explosion, affect the characterization of seismic sources when isotropic earth assumptions are made. This work is done in support of the modeling component of the Source Physics Experiment (SPE), a series of underground chemical explosions at the Nevada National Security Site (NNSS).

58 GEOSCIENCES↗
DOE OSTI2024

Effects of Nonzero-frequency Fluctuations on Turbulence Spectral Observations

In situ observations of turbulence spectra in space plasmas are usually interpreted as wavenumber spectra, assuming that the fluctuation frequency is negligible in the plasma flow frame. We explore the effects of nonzero frequency in the plasma flow frame on turbulence spectral observations. The finite frequency can be caused by either propagating waves or nonlinear broadening of nonpropagating structures. We show that the observed frequency spectrum can be modified by the nonzero frequency of turbulent fluctuations in several ways. Specifically, (i) frequency broadening results in a minor modification to the observed spectrum, primarily acting as a smoothing kernel of the spectrum near the spectral break, while the asymptotic spectral index remains unchanged; (ii) wave propagation can affect the observed spectral index for anisotropic turbulence. The effect is significant at low frequencies and weaker at high frequencies, leading to a "concave" shape of the observed perpendicular spectrum; (iii) the Doppler shift for forward- and backward-propagating Elsasser modes can result in a nonzero cross helicity for critical-balanced turbulence since the effect of the Doppler shift favors outward-propagating waves systematically, resulting in an observed imbalance. These results may have important implications for the interpretation of solar wind flows observed by Parker Solar Probe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Constraining the phase shift of relativistic species in DESI BAOs

In the early Universe, neutrinos decouple quickly from the primordial plasma and propagate without further interactions. The impact of free-streaming neutrinos is to create a temporal shift in the gravitational potential that impacts the acoustic waves known as baryon acoustic oscillations (BAOs), resulting in a non-linear spatial shift in the Fourier-space BAO signal. In this work, we make use of and extend upon an existing methodology to measure the phase shift amplitude $\beta _{\phi }$ and apply it to the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) BAOs with an anisotropic BAO fitting pipeline. We validate the fitting methodology by testing the pipeline with two publicly available fitting codes applied to highly precise cubic box simulations and realistic simulations representative of the DESI DR1 data. We find further study towards the methods used in fitting the BAO signal will be necessary to ensure accurate constraints on $\beta _{\phi }$ in future DESI data releases. Using DESI DR1, we present individual measurements of the anisotropic BAO distortion parameters and the $\beta _{\phi }$ for the different tracers, and additionally a combined fit to $\beta _{\phi }$ resulting in $\beta _{\phi } = 2.7 \pm 1.7$. After including a prior on the distortion parameters from constraints using Planck we find $\beta _{\phi } = 2.7^{+0.60}_{-0.67}$ suggesting $\beta _{\phi } > 0$ at 4.3$\sigma$ significance. This result may hint at a phase shift that is not purely sourced from the standard model expectation for $N_{\rm {eff}}$ or could be a upwards statistical fluctuation in the measured $\beta _{\phi }$; this result relaxes in models with additional freedom beyond Lambda-cold dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Solving the problem of overdetermination of quasisymmetric equilibrium solutions by near-axis expansions. II. Circular axis stellarator solutions

Here we apply the near-axis expansion method for quasisymmetric magnetic fields with anisotropic pressure to construct numerical solutions to circular axis stellarators. The solutions are found to second order in the distance from the axis, not possible in the standard Garren–Boozer construction, which assumes magnetostatic equilibria with isotropic pressure. In the limit of zero anisotropy, it is shown that a subset of coefficients can be chosen to avoid the overdetermination problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Pressure-driven tearing and thermal transport in finite-beta reversed field pinch computations

In this work, nonlinear resistive-magnetohydrodynamics (MHD) computation with heating and anisotropic transport is applied to examine the interaction between thermal energy and magnetic fluctuations in inductively driven reversed-field pinches (RFPs). The magnetic fluctuations underlie magnetic field reversal through dynamo-like correlations, and they enhance thermal energy transport through fluctuations of parallel heat flux density. With the unfavorable magnetic curvature that exists across the RFP profile, thermal energy also affects the magnetic fluctuations. Computations with the NIMROD code [Sovinec et al., J. Comput. Phys. 195, 355–386 (2004)] integrate nonlinear MHD dynamics with energy transport and reproduce an RFP state with experimentally relevant values of plasma-β. Equilibria constructed from results of the 3D computations are analyzed to assess the sources of free energy in the saturated nonlinear state. Linear computations for these profiles show unstable modes of tearing parity. Their eigenfunctions are used to evaluate and compare stabilizing and destabilizing contributions to the kinetic energy integral. An assessment of the drives in the integral reveals that the pressure gradient drive is of comparable magnitude to the parallel current drive, and only the sum of the two surpasses the stabilizing contributions. Correlation of magnetic and parallel heat flux density fluctuations in the nonlinear computations shows that fluctuation-induced thermal conduction is the dominant mode of energy loss, as expected from experimental evidence. Decomposition of the fluctuating heat flux density shows that second-order correlations, alone, do not explain the total energy transport. Higher-order correlations are also important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Extended-MHD simulations of disruption mitigation via massive gas injection in SPARC

Recent developments to the M3D-C1 code enable higher fidelity modeling of disruptions, and can be applied in the design verification of reactor-scale tokamaks. Among these new capabilities is a method to mesh conducting vessel structures such as coils and passive plates, packing of the toroidal mesh around gas injectors, as well as anisotropic resistivity inside the vessel structures. We present extended-magnetohydrodynamic (MHD) simulations of disruption mitigation via massive gas injection (MGI) in SPARC. The goal of this study is to inform the disruption mitigation layout of SPARC and aid in the design of an effective gas injector configuration. Fully three-dimensional simulations with M3D-C1 are carried out for various injector configurations with the primary goal of determining the effect of different MGI parameters on heat loads and vessel forces. The simulations include a model for impurity ionization, recombination, advection and radiation, as well as spatially resolved conducting structures around the plasma. A localized mixture of deuterium and neon with a small toroidal and poloidal width is injected in up to six locations. We demonstrate that M3D-C1 can model a rapid shutdown via MGI using narrow and more realistic gas plumes than in previous simulations. As a result of the q = 1 surface in the SPARC baseline case a sawtooth is observed early in the simulations. Despite the sawtooth and the onset of edge MHD instabilities, the impurity distribution remains localized around the injector locations, but enables a radiative shutdown of the plasma. We find that using the maximum of six gas injectors results in a lower peaking factor and leads to a more even distribution of radiation toroidally than using two injectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An efficient high-order numerical solver for diffusion equations with strong anisotropy

In this paper, we present an interior penalty discontinuous Galerkin finite element scheme for solving diffusion problems with strong anisotropy arising in magnetized plasmas for fusion applications. Additionally, we demonstrate the accuracy produced by the high-order scheme and develop an efficient preconditioning technique to solve the corresponding linear system, which is robust to the mesh size and anisotropy of the problem. Several numerical tests are provided to validate the accuracy and efficiency of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗
DOE OSTI2024

Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy

Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.

Astronomy & Astrophysics↗