Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MAGNETOHYDRODYNAMIC STABILITY”

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

The stability of a compressible stratified shear layer

The stability of a shear layer under the effect of gravity is investigated using the compressible magnetohydrodynamic (MHD) equations, including an effective gravity term to represent the curvature effects of the flow and magnetic field line geometry. A general eigenmode equation is derived for a two-dimensional MHD fluid, and an energy-principle analysis to explain the effect of compressibility on the critical Richardson number is presented. For the case of a hyperbolic tangent shear flow and exponential density profile, it was found that, in the Boussinesq approximation, the compressibility raises the critical Richardson number from 1/4 to as much as 1/2, with the exact value depending on the value of the magnetic field at infinity. Under approximation of a strong asymptotic magnetic field, without invoking the Boussinesq approximation, it is shown both analytically and numerically that the density gradient terms cause the shear instability to be dispersive. The long-wavelength stability boundary for the Richardson number J = 0 is characterized by a normalized phase velocity c =

Wang, Z.↗
NASA NTRS1981

Magnetic flux tubes

The magnetohydrodynamics of flux tubes are considered. The sections on equilibrium of flux tubes, and stability and waves deal with sunspots, the largest members of the general class of photospheric flux tubes.

Spruit, H. C.↗
NASA NTRS1992

The stability of the near-Earth magnetotail against the ballooning instability

In an extension to a previous study, conditions for the ballooning instability in the magnetohydrodynamic approximation are derived. It is necessary for exciting the ballooning instability that the plasma pressure is higher where plasma is convected outward. However, the expansion of plasma, which results from the divergence of outward flows in the curved geometry of field lines, may overcompensate for the convective change in the plasma pressure. Conditions for the instability are derived from the equation describing the linear coupling between the Alfven and slow magnetosonic waves. The result indicates that excessive curvature stabilizes, rather than destabilizes, perturbations, as well as the field aligned current and the field aligned flow. That is the field line curvature should be smaller than the pressure gradient for exciting the ballooning instability. It is inferred that this condition is not satisfied during quiet times in the near Earth tail. The plasma sheet becomes more stable to the instability in accordance with the thinning of the plasma sheets during the growth phase of substorms, while in the region nearer to the Earth (the inner edge of the plasma sheet or the ring current region) the plasma pressure gradient, which is possibly localized during the growth phase, may drive the instability.

Ohtani, Shin-Ichi↗
NASA NTRS1994

Anisotropic Alfven-ballooning modes in Earth's magnetosphere

We have carried out a theoretical analysis of the stability and parallel structure of coupled shear Alfven and slow magnetosonic waves in Earth's inner magnetopause (i.e., at equatorial distances between about five and ten Earth radii) including effects of finite anisotropic Grad-Shafranov equation yields an approximate self-consistent magnetohydrodynamic (MHD) equilibrium. This MHD equilibrium is used in the numerical solution of a set of eigenmode equations which describe the field line eigenfrequency, linear stability, and parallel eigenmode structure. We call these modes anisotropic Alfven-ballooning modes. The main results are: (1) The field line eigenfrequency can be significantly lowered by finite pressure effects. (2) The parallel mode structure of the transverse wave components is fairly insensitive to changes in the plasma pressure, but the compressional magnetic component can become highly peaked near the magnetic equator as a result of increased pressure, especially when P(sub perpendicular to) is greater than P(sub parallel) (here P(sub perpendicular to) and P(sub parallel) are the perpendicular and parallel plasma pressure). (3) For the isotropic (P(sub parallel) = P(sub perpendicular to) = P) case ballooning instability can occur when the ratio of the plasma presure to the magnetic pressure, beta = P/(B squared/8 pi), exceeds a critical value beta(sup B)(sub 0) is approximately equal to 3.5 at the equator. (4) Compared to the isotropic case the critical beta value is lowered by anisotropy, either due to decreased field line bending stabilization when P(sub parallel) is greater than P(sub perpendicular to) or due to increased ballooning-mirror destabilization when P(sub perpendicular to) is greater than P(sub parallel). (5) We use a beta-delta stability diagram to display the regions of instability with respect to the equatorial values of the parameters bar beta and delta, where bar beta = (1/3)(beta(sub parallel) + 2 beta(sub perpendicular to)) is an average beta value and delta = 1 - P(sub parallel)/P(sub perpendicular to) is a measure of the plasma anisotropy. The diagram is divided into regions corresponding to the firehose, mirror and ballooning instabilities. It appears that observed values of the plasma pressure are below the critical value for the isotropic ballooning instability but it may be possible to approach a ballooning-mirror instability when P(sub perpendicular to)/P(sub parallel) is greater than or approximately 2.

Chan, Anthony A.↗
NASA NTRS2011

Sloshing of the Magnetized Cool Gas in the Cores of Galaxy Clusters

X-ray observations of many clusters of galaxies reveal the presence of edges in surface brightness and temperature, known as "cold fronts". In relaxed clusters with cool cores, these edges have been interpreted as evidence for the "sloshing" of the core gas in the cluster's gravitational potential. The smoothness of these edges has been interpreted as evidence for the stabilizing effect of magnetic fields "draped" around the front surfaces. To check this hypothesis, we perform high-resolution magnetohydrodynamics simulations of magnetized gas sloshing in galaxy clusters initiated by encounters with subclusters. We go beyond previous works on the simulation of cold fronts in a magnetized intracluster medium by simulating their formation in realistic, idealized mergers with high resolution ((Delta)x approx. 2 kpc). Our simulations sample a parameter space of plausible initial magnetic field strengths and field configurations. In the simulations, we observe strong velocity shears associated with the cold fronts amplifying the magnetic field along the cold front surfaces, increasing the magnetic field strength in these layers by up to an order of magnitude, and boosting the magnetic pressure up to near-equipartition with thermal pressure in some cases. In these layers, the magnetic field becomes strong enough to stabilize the cold fronts against Kelvin-Helmholtz instabilities, resulting in sharp, smooth fronts as those seen in observations of real clusters. These magnetic fields also result in strong suppression of mixing of high and low-entropy gas in the cluster, seen in our simulations of mergers in the absence of a magnetic field. As a result, the heating of the core due to sloshing is very modest and is unable to stave off a cooling catastrophe.

ZuHone, J. A.↗
DOE OSTI2026

Verification of global gyrokinetic simulation of low frequency mode excited by thermal plasma in spherical tokamak

A low frequency Alfven eigenmode in the ST40 spherical tokamak (ST) is found to be resonantly excited by thermal ions from global gyrokinetic GTC simulations. The analyses of phase-space resonances and wave-particle energy exchanges show that the Beta-induced Alfvén-Acoustic Eigenmode (BAAE) is primarily driven unstable by thermal ions, but is stabilized by energetic particles (EPs). When all kinetic effects are artificially suppressed, GTC simulations in the ideal magnetohydrodynamic (MHD) limit confirm the existence of the MHD eigenmode with a frequency in the BAAEs gap, as previously found by the ideal MHD code NOVA. The kinetic effects of thermal ions induce a frequency downshift compared to the stable BAAE in the MHD simulations. These results support the verification for global gyrokinetic simulations of the BAAE in the ST. However, the simulations do not explain the experimental observations of the BAAE in the ST40, possibly owing to the limitation of the collisionless gyrokinetic Monte Carlo model for the EPs.

Alfvén eigenmode↗
DOE OSTI2025

Omnigenous stellarators with improved ideal and kinetic ballooning stability

Omnigenity is a property of a magnetic field which ensures confinement of trapped particles. It is a necessary requirement for any high-performance stellarator. After creating an omnigenous equilibrium, one must also ensure reduced transport resulting from kinetic and magnetohydrodynamic (MHD) instabilities. To this end, we leverage the GPU-accelerated DESC optimization suite, which is used to design stable, finite-β omnigenous equilibria with poloidal, toroidal, and helical symmetry, achieving Mercier, ideal ballooning, and as a consequence, improved kinetic ballooning stability. We discover stellarators with second stability, a regime of large pressure gradient where an equilibrium becomes ideal ballooning stable, and demonstrate and explore both using theory and gyrokinetic simulations the connection between ideal and kinetic ballooning stability.

optimization↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An experimental platform for investigating astrophysically relevant magnetized plasma jets on the COBRA facility

A new platform has been developed for the 1-MA COBRA generator to investigate the physical processes affecting the formation, collimation, and stability of high-speed outflows in magnetically driven laboratory plasma jets. Such experiments serve as diagnostically accessible surrogates for astrophysical jets under the assumption that the underlying dynamics are scale invariant. In contrast to previous current driven high energy density laboratory jet experiments that use radial/conical wire arrays or foils, the platform described here uses azimuthally symmetric gas-puff injection. This avoids the ablation phase from a solid target, allowing the jets to develop earlier and be driven longer without depleting their mass source and disrupting. A permanent magnet provides an initial poloidal magnetic field, which links the two concentric electrodes and mimics the boundary conditions of a star-accretion disk system. Extended magnetohydrodynamic effects can be assessed using a polarity convolute, which allows for reversal of the electrode bias. The resulting plasma jets exhibit remarkable stability, persisting for hundreds of nanoseconds and achieving aspect ratios ≳30:1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Hydromagnetic stability of coronal arcade structures The effects of photospheric line tying

A model is given of the magnetic-field equilibrium and possible dynamic excitations of a solar coronal arcade. Such structures are well observed in the spectral range from H-alpha to X-rays and often give rise to two-ribbon flares. However, the preflare state must be stable to ideal magnetohydrodynamic disturbances, and this problem is treated with particular attention to the necessary foot-point boundary conditions. With reasonably general perturbation set, an energy-principle analysis is used to show the strong stabilizing influence of inertial field-line tying at the photosphere.

Ray, A.↗
NASA NTRS1974

External field effects on solidification - Macroscopic and microscopic models

Comparison of theoretical effects of magnetic and gravitational fields on melts of metals and semiconductors provides two kinds of insight into space processing materials problems. First, the task of comparing effects of the two fields yields mathematical models which may be used to estimate material properties and responses to processing under various field conditions. Then, by using these models, microgravity effects can be inferred from magnetic field data already obtained in laboratories on earth. Magnetohydrodynamics (static and oscillating fields) and the Free Volume Model are used to study macroscopic effects, and Sekerka's interface stability theory is used to show how external fields can affect microsegregation. It is demonstrated that the mathematical formalism is essentially the same for describing macroscopic effects of both microgravity and magnetic fields.

Miller, R. I.↗
NASA NTRS2004

Studying Turbulence Using Numerical Simulation Databases - X Proceedings of the 2004 Summer Program

This Proceedings volume contains 32 papers that span a wide range of topics that reflect the ubiquity of turbulence. The papers have been divided into six groups: 1) Solar Simulations; 2) Magnetohydrodynamics (MHD); 3) Large Eddy Simulation (LES) and Numerical Simulations; 4) Reynolds Averaged Navier Stokes (RANS) Modeling and Simulations; 5) Stability and Acoustics; 6) Combustion and Multi-Phase Flow.

Moin, Parviz↗
DOE OSTI2026

Spherical tokamak physics research in preparation for the operation of NSTX-U

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high-β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations (δB ∥ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

NSTX-U↗
NASA NTRS1982

Magnetohydrodynamic equilibrium. III - Helically symmetric fields. IV - Nonequilibrium of nonsymmetric hydrodynamic topologies

It is pointed out that plasma confinement in stable equilibrium states constitutes a fundamental and still unresolved question in plasma astrophysics and thermonuclear fusion research. The problem has two parts related to the equilibrium states themselves and their mechanical stability. The question of the existence of general solutions of the field and fluid equations for the steady dynamical interaction of inviscid compressible fluids of high electrical conductivity with magnetic and gravity fields is considered. In the absence of fluid motions, the presented equations become the familiar equations of magnetostatics supplemented by an equation of state. Starting from this simplest case of magnetostatic equilibrium, the investigation proceeds to the more complex case of magnetohydrodynamic equilibrium. Examples of helically symmetric fields are presented to illustrate the use of the formulation for treating the dynamics of helically symmetric hydromagnetic flows.

Tsinganos, K. C.↗
NASA NTRS1994

A note on two-dimensional asymptotic magnetotail equilibria

In order to understand, on the fluid level, the structure, the time evolution, and the stability of current sheets, such as the magnetotail plasma sheet in Earth's magnetosphere, one has to consider magnetic field configurations that are in magnetohydrodynamic (MHD) force equilibrium. Any reasonable MHD current sheet model has to be two-dimensional, at least in an asymptotic sense (B(sub z)/B (sub x)) = epsilon much less than 1. The necessary two-dimensionality is described by a rather arbitrary function f(x). We utilize the free function f(x) to construct two-dimensional magnetotail equilibria are 'equivalent' to current sheets in empirical three-dimensional models. We obtain a class of asymptotic magnetotail equilibria ordered with respect to the magnetic disturbance index Kp. For low Kp values the two-dimensional MHD equilibria reflect some of the realistic, observation-based, aspects of three-dimensional models. For high Kp values the three-dimensional models do not fit the asymptotic MHD equlibria, which is indicative of their inconsistency with the assumed pressure function. This, in turn, implies that high magnetic activity levels of the real magnetosphere might be ruled by thermodynamic conditions different from local thermodynamic equilibrium.

Voigt, Gerd-Hannes↗
DOE OSTI2025

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗
DOE OSTI2024

Validation study of RWM stability in DIII-D high- β N plasmas

The n = 1 (n is the toroidal mode number) resistive wall mode (RWM) stability is numerically investigated for two DIII-D high-β N discharges 176440 and 172461, utilizing the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) codes. Systematic validation efforts are attempted, for the first time, for discharges with very slow or vanishing toroidal flow for a large fraction of the plasma volume. While gaining physics insights in accessing stable operation regime at β N exceeding the Troyon no-wall limit in these slow-rotation experiments, the predictive capability of fluid and non-perturbative magnetohydrodynamic-kinetic hybrid models for the RWM is further confirmed. The MARS-F fluid model, with a strong but numerically tunable viscosity mimicking ion Landau damping of parallel sound waves, finds complete stabilization of the n = 1 RWM in the considered DIII-D plasmas under the experimental flow conditions. Similarly, either full stabilization (for discharge 176440) or marginal stability (for discharge 172461) of the mode is computed by the MARS-K hybrid model, which is first-principle based without free model parameters. In particular, all drift kinetic resonances, including those of thermal and energetic particles, are found to synergistically act to marginally stabilize the RWM in discharge 172461. These MARS-F/K modeling results explain the experimentally observed stable operational regime in DIII-D, as far as the RWM stability is concerned. Extensive numerical sensitivity studies, with respect to the plasma toroidal flow speed as well as the radial location of the resistive wall, are also carried out to further support the validation study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗