Search NASA⌕ Search

SEARCH · Search NASA

Results for “magnetic fluids”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Kinetic effects of Alfven wave nonlinearity. I - Ponderomotive density fluctuations

The Vlasov theory is used to study kinetic corrections to fluid descriptions of Alfven wave nonlinearity. The method is to obtain an expression for the second-order perturbed distribution function produced by a nonlinear Alfven wave. From this distribution function a kinetically correct expression is obtained for the plasma density perturbation associated with an envelope-modulated Alfven wave. This kinetic theory result differs substantially from the fluid expression when the plasma beta is greater than about 1, and the electron and ion temperatures are approximately equal. This result is of interest because density fluctuations are an observationally accessible indicator of wave nonlinearity in solar system Alfven waves. It also will assist in the determination of properties of Alfven waves in the interstellar medium. Finally, this analysis also yields a kinetically correct expression for u, the magnetic field-aligned component of the plasma fluid velocity.

Spangler, Steven R.↗
NASA NTRS1990

Self-similar evolution of the nonlinear magnetic buoyancy instability

A new type of self-similar solution of ideal magnetohydrodynamics (MHD) in the nonlinear stage of the undular model (k parallel to B) of the magnetic buoyancy instability (the ballooning instability in fusion plasma physics or the Parker instability in astrophysics) is found through MHD simulation and theory. The linear theory developed agrees well with the simulation in the early (linear) stage. The nonlinear stages of the instability in the simulation show the self-similar evolution. One of the solutions obtained from the nonlinear analysis has the characteristics of nonlinear instability in Lagrangian coordinates; the fluid velocity and the Alfven speed on each magnetic loop increase exponentially with time, because the loop is evacuated by the field-aligned motion of matter resulting from gravitational acceleration. In the later stage of the nonlinear evolution, the solution property changes from exponential to power-law time dependence. The latter corresponds to a force-free expansion solution. The later saturation of the velocity increment is also discussed.

Shibata, K.↗
DOE OSTI2020

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Mapping the space of quasisymmetric stellarators using optimized near-axis expansion

A method is demonstrated to rapidly calculate the shapes and properties of quasi-axisymmetric and quasi-helically symmetric stellarators. In this approach, optimization is applied to the equations of magnetohydrodynamic equilibrium and quasisymmetry, expanded in the small distance from the magnetic axis, as formulated by Garren & Boozer [ Phys. Fluids B, vol. 3, 1991, p. 2805]. Due to the reduction of the equations by the expansion, the computational cost is significantly reduced, to times of the order of 1 cpu second, enabling wide and high-resolution scans over parameter space. In contrast to traditional stellarator optimization, here, the cost function serves to maximize the volume in which the expansion is accurate. A key term in the cost function is $\| \boldsymbol {\nabla }\boldsymbol B \|$ , the norm of the magnetic field gradient, to maximize scale lengths in the field. Using this method, a database of $5\times 10^5$ optimized configurations is calculated and presented. Quasisymmetric configurations are observed to exist in continuous bands, varying in the ratio of the magnetic axis length to average major radius. Several qualitatively new types of configuration are found, including quasi-helically symmetric fields in which the number of field periods is two or more than six.

Landreman, Matt (ORCID:000000027233577X)↗
DOE OSTI2022

Design of a high energy density experiment to measure the suppression of hydrodynamic instability in an applied magnetic field

A transverse magnetic field can suppress the hydrodynamic instability growth of an unstable plasma interface. This effect, of interest to inertial confinement fusion and astrophysics, has mostly been studied via simulation. Here, we present the design of an experiment at the National Ignition Facility (NIF) to demonstrate this effect in a laboratory. Simulations indicate that the timescale for the diffusion of the magnetic field across the mixing region should be at least comparable to the timescale of the instability growth in order to have a measurable suppression effect. This motivates the use of lower density target materials than usual high energy density (HED) hydrodynamics experiments to permit faster hydrodynamics and higher plasma conductivities (through higher temperature), for a given laser drive and magnetic field. We discuss a target design for creating a Rayleigh–Taylor unstable HED plasma interface that uses 320 mg/cc iodine-doped carbon foam as the heavy material, 20 mg/cc carbon foam as the light material, and a 6 lm amplitude, a 120 lm wavelength ripple machined at the interface, which shows a measurable suppression effect with a nominal NIF drive and a 30 T magnetic field (the present facility limit). Here, models indicate lower density foams that may display even larger suppression effects, as the Hohlraum drive also radiatively preheats the foam to permit even higher temperatures and, hence, higher conductivities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Coupled Electro-Thermal Analysis of Permanent Magnet Synchronous Motor for Electric Vehicles

Permanent magnet synchronous motors (PMSM) are extensively used in electric vehicles. However, high internal heat generation and inefficient heat dissipation often limit the operational reliability, and longevity of the PMSM. Therefore, proper quantification of heat generation in electric motor and advanced embedded motor cooling techniques remain topics of immense interest. In order to accurately predict electro-magnetic performance, i.e., efficiency, component-wise heat losses and the corresponding temperature distribution of a jacket cooled machine, this paper presents a two-way iteratively coupled electro-thermal modeling framework. Finite element based software Motor-CAD has been utilized for electro-magnetic performance calculation of BMW i3 PMSM. Finite volume based computational fluid dynamics/heat transfer (CFD/HT) software ANSYS® FLUENT® has been employed to simulate the temperature distribution of the PMSM, using the electro-magnetic losses as heat input. Computed heat losses, stator, winding, rotor, and magnet temperatures are utilized as coupling parameters between the electro-magnetic and thermal models. A conventional one-way coupling algorithm has also been developed and compared to the newly proposed two-way coupling algorithm. Numerical results confirm that at high current density, one-way coupling algorithm significantly over-predicts the motor temperature compared to the two-way algorithm. A comprehensive analysis has been carried out to characterize the influences of current density, speed, and forced convection heat transfer coefficient on the heat losses, overall efficiency, and maximum temperature of the PMSM. Lastly, an efficiency map has been interpreted from the coupled electro-magnetic simulation.

42 ENGINEERING↗
NASA NTRS1982

Investigation of geomagnetic field forecasting and fluid dynamics of the core

It was established that the total absolute magnetic flux crossing the core- mantle boundary has been a constant of the core motion for the last 50 years. This provides a scalar constraint that could be added to the geometric modelling procedure. The GSFC 8 8/80 model is being evaluated. The absolute magnetic flux linking the CMB to that model was plotted as a function of time during the span covered by the data, and increasing truncation level. The inclusion of the standard error of each Gauss coefficient derived from the statistics of fit in the GSFC 9/80 model is useful. The magnitude and sense (upwelling or downe welling) of vertical fluid motion adjacent to the core-mantle boundary was calculated using the model. Standard errors were found to be sufficiently small at all but one or two of the 40 or more critical points of B sub r. They do not nearly overlap the value gamma u/gamma r = 0. It is concluded that the core is upwelling and downwelling at an observationally detectable level.

Benton, E. R.↗
DOE OSTI2020

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave–particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel & Engelmann ( Phys. Fluids , vol. 9, 1966, pp. 2377–2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that scatter them out of and into resonance, as for Landau damping as shown by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). As a result, the resonant particle–wave interactions occur in the centre of a narrow collisional boundary when the collision frequency $\unicode[STIX]{x1D708}$ is very small compared to the wave frequency $\unicode[STIX]{x1D714}$ . The diffusive nature of the pitch angle scattering combined with the wave–particle resonance condition enhances the collision frequency by $(\unicode[STIX]{x1D714}/\unicode[STIX]{x1D708})^{2/3}\gg 1$ , resulting in an effective resonant particle collisional interaction time of $\unicode[STIX]{x1D70F}_{\text{int}}\sim (\unicode[STIX]{x1D708}/\unicode[STIX]{x1D714})^{2/3}/\unicode[STIX]{x1D708}\ll 1/\unicode[STIX]{x1D708}$ . A collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel–Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

Physics↗
NASA NTRS1978

Small amplitude waves in a hot relativistic two-fluid plasma

The dispersion relation for small amplitude waves in a hot relativistic plasma embedded in a uniform magnetic field is derived. Four plasma modes exist for propagation oblique to the magnetic field, and the mode properties are summarized in a hot relativistic generalization of the cold plasma pond diagram of Stix (1962). The two-fluid approximation is considered, and two-fluid equations are applied in this procedure which summarizes properties of linear waves in an unbounded magnetized relativistic plasma by means of a parameter-space diagram of wave phase velocity.

Hyun, S.↗
DOE OSTI2020

Final Report for Tech-X Corporations’s contribution to Enhancing Understanding of High Energy Density Plasmas from Wire Array and Solid Liner Implosions Using Fluid Modeling with Kinetic Closures

Recent results from experiments and simulation of magnetically driven pulsed power liners have explored the role of the early-time electrothermal instability in the evolution of the magneto-Rayleigh-Taylor instability. Our focus will be on understanding the development of such instabilities and the potential stabilization mechanisms via electron thermal conduction, viscosity, and large magnetic fields which we expect could play a significant role in supporting the success of the MagLIF program. Experiments have shown that studies of high-energy density plasmas from wire-array implosions require physics modeling that goes well beyond simple models such as ideal magnetohydrodynamics. The goal of this work is to provide increased understanding of these experiments by employing simulations with a multi-fluid extended-MHD model which uses kinetic closures for thermal conduction, resistivity and viscosity. We will use codes easily available to the wider research community, including university students, with a secondary goal of providing the community with well-benchmarked tools capable of advanced modeling of high-energy-density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Laser-induced-fluorescence imaging of a spin-polarized ultracold neutral plasma in a magnetic field

In this work, we report Doppler-sensitive laser-induced-fluorescence (LIF) imaging of an ultracold neutral plasma in a magnetic field. Local values of ion density, hydrodynamic fluid velocity, temperature, and spin polarization are obtained using a fluorescence model based on velocity-resolved rate equations (REs) including the transfer of ions between states due to laser coupling and spontaneous emission. The RE approach captures optical pumping of ions into states that are not driven by the LIF excitation laser, and this is validated with experimental data. Combined molecular-dynamics and quantum-trajectories simulations verify that velocity-changing collisions have a negligible impact on the state population evolution for typical experimental conditions. Relative intensities of Zeeman components of the LIF spectra provide clear evidence that the ions are electron-spin-polarized when created by photoionization of magnetically trapped 88 Sr atoms. This probe opens many possibilities for studying thermal transport and the equilibration of neutral plasmas in overlapping regimes of strong coupling and magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

Polarization and Compressibility of Oblique Kinetic Alfven Waves

Even though solar wind, as a collisionless plasma, is properly described by the kineticMaxwell-Vlasov description, it can be argued that much of our understanding of solar wind observational data comes from an interpretation and numerical modeling which is based on a fluid description of magnetohydrodynamics. In recent years, there has been a significant interest in better understanding the importance of kinetic effects, i.e. the differences between the kinetic and usual fluid descriptions. Here we concentrate on physical properties of oblique kinetic Alfvn waves (KAWs), which are often recognized as one of the key ingredients in the solar wind turbulence cascade. We use three different fluid models with various degrees of complexity and calculate polarization and magnetic compressibility of oblique KAWs (propagation angle q = 88), which we compare to solutions derived from linear kinetic theory. We explore a wide range of possible proton plasma b = [0.1,10.0] and a wide range of length scales krL = [0.001,10.0]. It is shown that the classical isotropic two-fluid model is very compressible in comparison with kinetic theory and that the largest discrepancy occurs at scales larger than the proton gyroscale. We also show that the two-fluid model contains a large error in the polarization of electric field, even at scales krL 1. Furthermore, to understand these discrepancies between the two-fluid model and the kinetic theory, we employ two versions of the Landau fluid model that incorporate linear low-frequency kinetic effects such as Landau damping and finite Larmor radius (FLR) corrections into the fluid description. It is shown that Landau damping significantly reduces the magnetic compressibility and that FLR corrections (i.e. nongyrotropic contributions) are required to correctly capture the polarization.We also show that, in addition to Landau damping, FLR corrections are necessary to accurately describe the damping rate of KAWs. We conclude that kinetic effects are important even at scales which are significantly larger than the proton gyroscale krL 1.

Turbulence↗
DOE OSTI2022

Electron transport in a collisional plasma with multiple ion species in the presence of a magnetic field

The classical work of Braginskii [Zh. Eksp. Teor. Fiz. 33, 459 (1957)] published almost 65 years ago was the first to provide a complete, closed fluid description of a weakly coupled, fully ionized, collisional plasma immersed in a magnetic field. While this fact is not widely known or appreciated, the Braginskii expressions for the electron transport coefficients can under- or overestimate the said coefficients by up to a factor of two for the electron Hall parameter of order unity (with the Hall parameter being proportional to the electron gyro-frequency over the electron–ion collision frequency) and can provide incorrect Hall-parameter scalings for its large values. Starting with the work of Epperlein and Haines [Phys. Fluids 29, 1029 (1986)], several papers attempted to correct the Braginskii electron results with varying degrees of success. Herein, we present our own effort with a hope to finally put this problem to rest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Magnetic blowout during collisionless reconnection

Analytic calculations and simulations of the steady-state structure of the dissipation region during collisionless reconnection are described. Ions are treated as particles and electrons as a dissipative fluid. A limiting reconnection rate results as the ions 'blowout' the magnetic field in the dissipation area. The implications for comprehending substorm dynamics are discussed.

Drake, J. F.↗
DOE OSTI2021

Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows

The guiding center and gyrokinetic theory of magnetized particle motion is extended to the regime of large electric field gradients perpendicular to the magnetic field. A gradient in the electric field directly modifies the oscillation frequency and causes the Larmor orbits to deform from circular to elliptical trajectories. In order to retain a good adiabatic invariant, there can only be strong dependence on a single coordinate at lowest order, so that resonances do not generate chaotic motion that destroys the invariant. When the gradient across magnetic flux surfaces is dominant, the guiding center drift velocity becomes anisotropic in response to external forces and additional curvature drifts must be included. The electric polarization density remains gyrotropic, but both the polarization and magnetization are modified by the change in gyrofrequency. The theory can be applied to shear flows that are even stronger than those observed in the edge transport barrier of a high-performance tokamak (H-mode) pedestal, even if the toroidal field is as small as or even smaller than the poloidal field. Yet, the theory retains a mathematical form that is similar to the standard case and can readily be implemented within existing simulation tools.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Initial Development of Fusion Magnet Simulation Capabilities for Performance and Safety Evaluation Using the MOOSE Framework

Fusion energy holds the promise of being a transformative technology as a carbon-neutral, sustainable source of energy. Whole device modeling and the development of fusion digital twins will be increasingly important for emerging fusion device concepts at both national laboratories and within the commercial fusion industry. However, meeting the challenge of whole device modeling of fusion energy devices requires robust, multiphysics, multiscale modeling and simulation technologies capable of running on large-scale supercomputers. Detailed analysis of individual systems at-scale is also required to ensure safe and efficient operation as well as provide the safety basis for future device designs and licensing activities. In a tokamak, toroidal and poloidal magnets confine and shape the fusion plasma to promote the fusion reaction. High plasma temperatures and high magnetic field requirements in modern design concepts (leading to high amounts of energy stored within each magnet) impose electrical, thermal, and mechanical loads on the magnet components, which in turn impacts the safety considerations of the magnet and their supporting systems. Idaho National Laboratory (INL) has a history of working in this space, including development and benchmarking of the Magnetic System Circuitry Analysis Program (MSCAP) and Magnet Arcing (MAGARC) codes to study magnet quench events; notably, MAGARC was used to study quenching during the ITER Engineering Design Activity. However, these legacy codes and capabilities are not parallel and scalable, and new tools are required for future advances in this area, which leads to the INL-developed Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Developed originally for fission reactor systems under United States Department of Energy, Office of Nuclear Energy modeling and simulation programs, the MOOSE framework has been well-suited to multiscale, multiphysics modeling and simulation needs for nuclear systems. The framework is open-source, well-tested, under continuous development and deployment, and developed to a Nuclear Quality Assurance, Level 1 software quality standard. MOOSE has also been used in the fusion space previously in several projects: INL’s Tritium Migration Analysis Program, Version 8 (TMAP8) for tritium migration, UK Atomic Energy Authority’s A Unified Resource for OpenMC (fusion) Reactor Applications (AURORA) code for fusion thermo-mechanical and neutronics analysis, and Argonne National Laboratory’s Cardinal for high-fidelity computational fluid dynamics and neutronics. However, to model superconducting magnets, several MOOSE enhancements are required: additions to the current MOOSE electromagnetic capabilities, new material libraries for superconductors of interest (such as YBCO), as well as fusion-specific models for thermo-mechanics. This talk will discuss initial development activities to build these capabilities in MOOSE, focusing on initial validation and benchmarking activities. Proposed coupling workflows and future work to support the simulation of fusion magnets and magnet structural assemblies for performance and safety evaluation in MOOSE will also be discussed.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Sub-Alfvénic Turbulence: Magnetic-to-kinetic Energy Ratio, Modification of Weak Cascade, and Implications for Magnetic Field Strength Measurements

Abstract We study the properties of sub-Alfvénic magnetohydrodynamic (MHD) turbulence, i.e., turbulence with Alfvén mach number M A = V L / V A < 1, where V L is the velocity at the injection scale and V A is the Alfvén velocity. We demonstrate that MHD turbulence can have different properties, depending on whether it is driven by velocity or magnetic fluctuations. If the turbulence is driven by isotropic bulk forces acting upon the fluid, i.e., is velocity driven, in an incompressible conducting fluid we predict that the kinetic energy is M A − 2 times larger than the energy of magnetic fluctuations. This effect arises from the long parallel wavelength tail of the forcing, which excites modes with k ∥ / k ⊥ < M A . We also predict that as the MHD turbulent cascade reaches the strong regime, the energy of slow modes exceeds the energy of Alfvén modes by a factor M A − 1 . These effects are absent if the turbulence is driven through magnetic fluctuations at the injection scale. We confirm our predictions with numerical simulations. Since the assumption of magnetic and kinetic energy equipartition is at the core of the Davis–Chandrasekhar–Fermi (DCF) approach to measuring magnetic field strength in sub-Alfvénic turbulence, we conclude that the DCF technique is not universally applicable. In particular, we suggest that the dynamical excitation of long azimuthal wavelength modes in the galactic disk may compromise the use of the DCF technique. We discuss alternative expressions that can be used to obtain magnetic field strength from observations and consider ways of distinguishing the cases of velocity and magnetically driven turbulence using observational data.

79 ASTRONOMY AND ASTROPHYSICS↗