SEARCH · Search NASA

Results for “ideal MHD”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

339 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS2016

FORWARD: A Toolset for Multiwavelength Coronal Magnetometry

Determining the 3D coronal magnetic field is a critical, but extremely difficult problem to solve. Since different types of multiwavelength coronal data probe different aspects of the coronal magnetic field, ideally these data should be used together to validate and constrain specifications of that field. Such a task requires the ability to create observable quantities at a range of wavelengths from a distribution of magnetic field and associated plasma i.e., to perform forward calculations. In this paper we describe the capabilities of the FORWARD SolarSoft IDL package, a uniquely comprehensive toolset for coronal magnetometry. FORWARD is a community resource that may be used both to synthesize a broad range of coronal observables, and to access and compare synthetic observables to existing data. It enables forward fitting of specific observations, and helps to build intuition into how the physical properties of coronal magnetic structures translate to observable properties. FORWARD can also be used to generate synthetic test beds from MHD simulations in order to facilitate the development of coronal magnetometric inversion methods, and to prepare for the analysis of future large solar telescope data.

Sun: magnetic fields↗

NASA NTRS2015

The Substorm Current Wedge: Further Insights from MHD Simulations

Using a recent magnetohydrodynamic simulation of magnetotail dynamics, we further investigate the buildup and evolution of the substorm current wedge (SCW), resulting from flow bursts generated by near-tail reconnection. Each flow burst generates an individual current wedge, which includes the reduction of cross-tail current and the diversion to region 1 (R1)-type field-aligned currents (earthward on the dawn and tailward on the duskside), connecting the tail with the ionosphere. Multiple flow bursts generate initially multiple SCW patterns, which at later times combine to a wider single SCW pattern. The standard SCWmodel is modified by the addition of several current loops, related to particular magnetic field changes: the increase of Bz in a local equatorial region (dipolarization), the decrease of |Bx| away from the equator (current disruption), and increases in |By| resulting from azimuthally deflected flows. The associated loop currents are found to be of similar magnitude, 0.1-0.3 MA. The combined effect requires the addition of region 2 (R2)-type currents closing in the near tail through dawnward currents but also connecting radially with the R1 currents. The current closure at the inner boundary, taken as a crude proxy of an idealized ionosphere, demonstrates westward currents as postulated in the original SCW picture as well as North-South currents connecting R1- and R2-type currents, which were larger than the westward currents by a factor of almost 2. However, this result should be applied with caution to the ionosphere because of our neglect of finite resistance and Hall effects.

reconnection↗

DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

On large-scale dynamos with stable stratification and the application to stellar radiative zones

Our understanding of large-scale magnetic fields in stellar radiative zones remains fragmented and incomplete. Such magnetic fields, which must be produced by some form of dynamo mechanism, are thought to dominate angular-momentum transport, making them crucial to stellar evolution. A major difficulty is the effect of stable stratification, which generally suppresses dynamo action. We explore the effects of stable stratification on mean-field dynamo theory with a particular focus on a non-helical large-scale dynamo (LSD) mechanism known as the magnetic shear-current effect. We find that the mechanism is robust to increasing stable stratification as long as the original requirements for its operation are met: a source of shear and non-helical magnetic fluctuations (e.g. from a small-scale dynamo). Both are plausibly sourced in the presence of differential rotation. Our idealized direct numerical simulations, supported by mean-field theory, demonstrate the generation of near equipartition large-scale toroidal fields. Additionally, a scan over magnetic Reynolds number shows no change in the growth or saturation of the LSD, providing good numerical evidence of a dynamo mechanism resilient to catastrophic quenching, which has been an issue for helical dynamos. These properties – the absence of catastrophic quenching and robustness to stable stratification – make the mechanism a plausible candidate for generating in situ large-scale magnetic fields in stellar radiative zones.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1976

Closed cycle MHD power generation experiments using a helium-cesium working fluid in the NASA Lewis Facility

The MHD channel in the NASA Lewis Research Center was redesigned and used in closed cycle power generation experiments with a helium-cesium working fluid. The cross sectional dimensions of the channel were reduced to 5 by 16.5 cm to allow operation over a variety of conditions. Experiments have been run at temperatures of 1900-2100 K and Mach numbers from 0.3 to 0.55 in argon and 0.2 in helium. Improvements in Hall voltage isolation and seed vaporization techniques have resulted in significant improvements in performance. Typical values obtained with helium are Faraday open circuit voltage 141 V (92% of uBh) at a magnetic field strength of 1.7 T, power outputs of 2.2 kw for tests with 28 electrodes and 2.1 kw for tests with 17 electrodes. Power densities of 0.6 MW/cu m and Hall fields of about 1100 V/m were obtained in the tests with 17 electrodes, representing a factor of 18 improvement over previously reported results. The V-I curves and current distribution data indicate that while near ideal equilibrium performance is obtained under some conditions, no nonequilibrium power has been generated to date.

Sovie, R. J.↗

NASA NTRS2016

The Magnetic Rayleigh-Taylor Instability in Astrophysical Discs

This is our first study of the magnetic Rayleigh-Taylor instability at the inner edge of an astrophysical disc around a central back hole. We derive the equations governing small-amplitude oscillations in general relativistic ideal magnetodydrodynamics and obtain a criterion for the onset of the instability. We suggest that static disc configurations where magnetic field is held by the disc material are unstable around a Schwarzschild black hole. On the other hand, we find that such configurations are stabilized by the space-time rotation around a Kerr black hole. We obtain a crude estimate of the maximum amount of poloidal magnetic flux that can be accumulated around the centre, and suggest that it is proportional to the black hole spin. Finally, we discuss the astrophysical implications of our result for the theoretical and observational estimations of the black hole jet power.

MHD↗

DOE OSTI2021

Formation and evolution of protostellar accretion discs – II. From 3D simulation to a simple semi-analytic model of Class 0/I discs

Here we use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Parallelized real-time physics codes for plasma control on DIII-D

A real-time safe multi-threading library was developed on the DIII-D plasma control system to optimize the real-time TORBEAM and real-time STRIDE physics codes. These physics codes are crucial for future fusion power plant operation as they provide information about electron cyclotron wave propagation and heating as well as inform about ideal plasma stability limits. The real-time TORBEAM code executed consistently in under 20 ms while the real-time STRIDE code computes in 100 ms. The multi-threading library developed in this work can be applied to other real-time physics-based codes that will be crucial for the next generation of fusion devices.

DIII-D↗

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↗

NASA NTRS2019

A Comprehensive Chandra Study of the Disk Wind in the Black Hole Candidate 4U 1630-472

The mechanisms that drive disk winds are a window into the physical processes that underlie the disk. Stellar-mass black holes are an ideal setting in which to explore these mechanisms, in part because their outbursts span a broad range in mass accretion rate. We performed a spectral analysis of the disk wind found in six Chandra/HETG observations of the black hole candidate 4U 1630−472, covering a range of luminosities over two distinct spectral states. We modeled both wind absorption and extended wind re-emission components using PION, a selfconsistent photoionized absorption model. In all but one case, two photoionization zones were required in order to obtain acceptable fits. Two independent constraints on launching radii, obtained via the ionization parameter formalism and the dynamical broadening of the re-emission, helped characterize the geometry of the wind. The innermost wind components (r ≃ 10(sup 2-3)GM/c(exp 2)) tend toward small volume filling factors, high ionization, densities up to n ≃ 10(sup 15-16) cm(exp -3), and outflow velocities of ∼0.003c. These small launching radii and large densities require magnetic driving, as they are inconsistent with numerical and analytical treatments of thermally driven winds. Outer wind components (r ≃ 10(exp 5)GM/c(exp 2)) are significantly less ionized and have filling factors near unity. Their larger launching radii, lower densities (n ≃ 10(exp 12) cm(exp -3)), and outflow velocities (∼0.0007c) are nominally consistent with thermally driven winds. The overall wind structure suggests that these components may also be part of a broader MHD outflow and perhaps better described as magneto-thermal hybrid winds.

Black hole physics↗

DOE OSTI2025

Assessing the numerical stability of physics models to equilibrium variation through database comparisons on DIII-D

High fidelity kinetic equilibria are crucial for tokamak modeling and analysis. Manual workflows for constructing kinetic equilibria are time consuming and subject to user error, motivating development of automated equilibrium reconstruction tools to provide accurate and consistent reconstructions for downstream physics analysis. These automated tools also provide access to kinetic equilibria at large database scales, which enables the quantification of general uncertainties arising from equilibrium reconstruction techniques. In this paper, we compare a large database of DIII-D kinetic equilibria generated manually by physics experts to equilibria from automated kinetic reconstruction tools, assessing the impact of reconstruction method on equilibrium parameters and resulting magnetohydrodynamic stability calculations. We find agreement among scalar parameters, whereas profile quantities, such as the bootstrap current, show larger disagreements. We analyze ideal kink and classical tearing stability with DCON and STRIDE respectively, finding that the kink stability calculation is generally more robust than the tearing index Δ' calculation. We find that in 90% of cases, both kink stability classifications are unchanged between the manual expert and automated kinetic equilibria.

CAKE↗

DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Because flows in stellar convection zones are highly turbulent, multidimensional hydrodynamic simulations are fundamental to accurately capture the physics of mixing processes. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible S EVEN -L EAGUE H YDRO code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of 25 M ⊙ . We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as 0.33, and the average magnetic field strength is ∼10 10 G. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to 60% of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ≈20% as compared to the purely hydrodynamic case.

79 ASTRONOMY AND ASTROPHYSICS↗