Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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.

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

At least 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

NON-EQUILIBRIUM HELIUM IONIZATION IN AN MHD SIMULATION OF THE SOLAR ATMOSPHERE

The ionization state of the gas in the dynamic solar chromosphere can depart strongly from the instantaneous statistical equilibrium commonly assumed in numerical modeling. We improve on earlier simulations of the solar atmosphere that only included non-equilibrium hydrogen ionization by performing a 2D radiation-magnetohydrodynamics simulation featuring non-equilibrium ionization of both hydrogen and helium. The simulation includes the effect of hydrogen Lyα and the EUV radiation from the corona on the ionization and heating of the atmosphere. Details on code implementation are given. We obtain helium ion fractions that are far from their equilibrium values. Comparison with models with local thermodynamic equilibrium (LTE) ionization shows that non-equilibrium helium ionization leads to higher temperatures in wavefronts and lower temperatures in the gas between shocks. Assuming LTE ionization results in a thermostat-like behavior with matter accumulating around the temperatures where the LTE ionization fractions change rapidly. Comparison of DEM curves computed from our models shows that non-equilibrium ionization leads to more radiating material in the temperature range 11–18 kK, compared to models with LTE helium ionization. We conclude that non-equilibrium helium ionization is important for the dynamics and thermal structure of the upper chromosphere and transition region. It might also help resolve the problem that intensities of chromospheric lines computed from current models are smaller than those observed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

s-PROCESSING IN AGB STARS REVISITED. II. ENHANCED {sup 13}C PRODUCTION THROUGH MHD-INDUCED MIXING

Slow neutron captures are responsible for the production of about 50% of elements heavier than iron, mainly occurring during the asymptotic giant branch phase of low-mass stars (1 ≲ M/M{sub ⊙} ≲ 3), where the main neutron source is the {sup 13}C(α, n){sup 16}O reaction. This last reaction is activated from locally produced {sup 13}C, formed by partial mixing of hydrogen into the He-rich layers. We present here the first attempt to describe a physical mechanism for the formation of the {sup 13}C reservoir, studying the mass circulation induced by magnetic buoyancy without adding new free parameters to those already involved in stellar modeling. Our approach represents the application to the stellar layers relevant for s-processing of recent exact analytical 2D and 3D models for magneto-hydrodynamic processes at the base of convective envelopes in evolved stars in order to promote downflows of envelope material for mass conservation during the occurrence of a dredge-up phenomenon. We find that the proton penetration is characterized by small concentrations, but is extended over a large fractional mass of the He-layers, thus producing {sup 13}C reservoirs of several 10{sup −3} M{sub ⊙}. The ensuing {sup 13}C-enriched zone has an almost flat profile, while only a limited production of {sup 14}N occurs. In order to verify the effects of our new findings we show how the abundances of the main s-component nuclei can be accounted for in solar proportions and how our large {sup 13}C-reservoir allows us to solve a few so far unexplained features in the abundance distribution of post-AGB objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A FOUR-FLUID MHD MODEL OF THE SOLAR WIND/INTERSTELLAR MEDIUM INTERACTION WITH TURBULENCE TRANSPORT AND PICKUP PROTONS AS SEPARATE FLUID

We have developed a four-fluid, three-dimensional magnetohydrodynamic model of the solar wind interaction with the local interstellar medium. The unique features of the model are: (a) a three-fluid description for the charged components of the solar wind and interstellar plasmas (thermal protons, electrons, and pickup protons), (b) the built-in turbulence transport equations based on Reynolds decomposition and coupled with the mean-flow Reynolds-averaged equations, and (c) a solar corona/solar wind model that supplies inner boundary conditions at 40 au by computing solar wind and magnetic field parameters outward from the coronal base. The three charged species are described by separate energy equations and are assumed to move with the same velocity. The fourth fluid in the model is the interstellar hydrogen which is treated by separate continuity, momentum, and energy equations and is coupled with the charged components through photoionization and charge exchange. We evaluate the effects of turbulence transport and pickup protons on the global heliospheric structure and compute the distribution of plasma, magnetic field, and turbulence parameters throughout the heliosphere for representative solar minimum and maximum conditions. We compare our results with Voyager 1 observations in the outer heliosheath and show that the relative amplitude of magnetic fluctuations just outside the heliopause is in close agreement with the value inferred from Voyager 1 measurements by Burlaga et al. The simulated profiles of magnetic field parameters in the outer heliosheath are in qualitative agreement with the Voyager 1 observations and with the analytical model of magnetic field draping around the heliopause of Isenberg et al.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SOLAR PROMINENCES EMBEDDED IN FLUX ROPES: MORPHOLOGICAL FEATURES AND DYNAMICS FROM 3D MHD SIMULATIONS

The temporal evolution of a solar prominence inserted in a three-dimensional magnetic flux rope is investigated numerically. Using the model of Titov and Démoulin under the regime of weak twist, the cold and dense prominence counteracts gravity by modifying the initially force-free magnetic configuration. In some cases a quasi-stationary situation is achieved after the relaxation phase, characterized by the excitation of standing vertical oscillations. These oscillations show a strong attenuation with time produced by the mechanism of continuum damping due to the inhomogeneous transition between the prominence and solar corona. The characteristic period of the vertical oscillations does not depend strongly on the twist of the flux rope. Nonlinearity is responsible for triggering the Kelvin–Helmholtz instability associated with the vertical oscillations and that eventually produces horizontal structures. Contrary to other configurations in which the longitudinal axis of the prominence is permeated by a perpendicular magnetic field, like in unsheared arcades, the orientation of the prominence along the flux rope axis prevents the development of Rayleigh–Taylor instabilities and therefore the appearance of vertical structuring along this axis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE COOL SURGE FOLLOWING FLUX EMERGENCE IN A RADIATION-MHD EXPERIMENT

Cool and dense ejections, typically Hα surges, often appear alongside EUV or X-ray coronal jets as a result of the emergence of magnetized plasma from the solar interior. Idealized numerical experiments explain those ejections as being indirectly associated with the magnetic reconnection taking place between the emerging and preexisting systems. However, those experiments miss basic elements that can importantly affect the surge phenomenon. In this paper we study the cool surges using a realistic treatment of the radiation transfer and material plasma properties. To that end, the Bifrost code is used, which has advanced modules for the equation of state of the plasma, photospheric and chromospheric radiation transfer, heat conduction, and optically thin radiative cooling. We carry out a 2.5D experiment of the emergence of magnetized plasma through (meso) granular convection cells and the low atmosphere to the corona. Through detailed Lagrange tracing we study the formation and evolution of the cool ejection and, in particular, the role of the entropy sources; this allows us to discern families of evolutionary patterns for the plasma elements. In the launch phase, many elements suffer accelerations well in excess of gravity; when nearing the apex of their individual trajectories, instead, the plasma elements follow quasi-parabolic trajectories with accelerations close to g{sub ⊙}. We show how the formation of the cool ejection is mediated by a wedge-like structure composed of two shocks, one of which leads to the detachment of the surge from the original emerged plasma dome.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Nonaxisymmetric MHD instabilities of Chandrasekhar states in Taylor-Couette geometry

We consider axially periodic Taylor–Couette geometry with insulating boundary conditions. The imposed basic states are so-called Chandrasekhar states, where the azimuthal flow U {sub ϕ} and magnetic field B {sub ϕ} have the same radial profiles. Mainly three particular profiles are considered: the Rayleigh limit, quasi-Keplerian, and solid-body rotation. In each case we begin by computing linear instability curves and their dependence on the magnetic Prandtl number Pm. For the azimuthal wavenumber m = 1 modes, the instability curves always scale with the Reynolds number and the Hartmann number. For sufficiently small Pm these modes therefore only become unstable for magnetic Mach numbers less than unity, and are thus not relevant for most astrophysical applications. However, modes with m>1 can behave very differently. For sufficiently flat profiles, they scale with the magnetic Reynolds number and the Lundquist number, thereby allowing instability also for the large magnetic Mach numbers of astrophysical objects. We further compute fully nonlinear, three-dimensional equilibration of these instabilities, and investigate how the energy is distributed among the azimuthal (m) and axial (k) wavenumbers. In comparison spectra become steeper for large m, reflecting the smoothing action of shear. On the other hand kinetic and magnetic energy spectra exhibit similar behavior: if several azimuthal modes are already linearly unstable they are relatively flat, but for the rigidly rotating case where m = 1 is the only unstable mode they are so steep that neither Kolmogorov nor Iroshnikov–Kraichnan spectra fit the results. The total magnetic energy exceeds the kinetic energy only for large magnetic Reynolds numbers Rm>100.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Cosmic rays and MHD turbulence generation in interstellar giant molecular clouds

The diffusive propagation of nonrelativistic cosmic ray (CR) protons undergoing energy losses by ionization in a dense homogeneous infinitely extended interstellar molecular cloud (MC) is investigated. The steady-state transport equation for the differential number density of nonrelativistic CR protons is solved with the boundary condition that at the edge of cloud it agrees with the interstellar CR number density. It is shown that giant interstellar MCs with column depths much greater than about 7⋅10{sup 22} cm{sup −2} are an efficient sink of nonrelativistic CRs. At small penetration depths the CRs lose energy by ionizing and heating the molecular gas, whereas at large penetration depths they are collectively dissipated by the streaming instability, which transfers one-half of the energy density of the incoming interstellar nonrelativistic CRs to Alfvénic magnetic field turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

MHD SIMULATIONS OF CORONAL SUPRA-ARCADE DOWNFLOWS INCLUDING ANISOTROPIC THERMAL CONDUCTION

Coronal supra-arcade downflows (SADs) are observed as dark trails descending toward hot turbulent-fan-shaped regions. Due to the large temperature values and gradients in these fan regions, the thermal conduction (TC) should be very efficient. While several models have been proposed to explain the triggering and the evolution of SADs, none of these scenarios address a systematic consideration of TC. Thus, we accomplish this task numerically simulating the evolution of SADs within this framework. That is, SADs are conceived as voided (subdense) cavities formed by nonlinear waves triggered by downflowing bursty localized reconnection events in a perturbed hot fan. We generate a properly turbulent fan, obtained by a stirring force that permits control of the energy and vorticity input in the medium where SADs develop. We include anisotropic TC and consider plasma properties consistent with observations. Our aim is to study whether it is possible to prevent SADs from vanishing by thermal diffusion. We find that this will be the case, depending on the turbulence parameters, in particular if the magnetic field lines are able to envelope the voided cavities, thermally isolating them from the hot environment. Velocity shear perturbations that are able to generate instabilities of the Kelvin–Helmholtz type help to produce magnetic islands, extending the lifetime of SADs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

FARGO3D: A NEW GPU-ORIENTED MHD CODE

We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on the physics of protoplanetary disks and planet–disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite-difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on either graphical processing units (GPUs) or central processing units (CPUs), achieving large speed-up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for CPUs, and have them translated automatically for GPUs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Testing the Accuracy of Data-driven MHD Simulations of Active Region Evolution

Models for the evolution of the solar coronal magnetic field are vital for understanding solar activity, yet the best measurements of the magnetic field lie at the photosphere, necessitating the development of coronal models which are “data-driven” at the photosphere. We present an investigation to determine the feasibility and accuracy of such methods. Our validation framework uses a simulation of active region (AR) formation, modeling the emergence of magnetic flux from the convection zone to the corona, as a ground-truth data set, to supply both the photospheric information and to perform the validation of the data-driven method. We focus our investigation on how the accuracy of the data-driven model depends on the temporal frequency of the driving data. The Helioseismic and Magnetic Imager on NASA’s Solar Dynamics Observatory produces full-disk vector magnetic field measurements at a 12-minute cadence. Using our framework we show that ARs that emerge over 25 hr can be modeled by the data-driving method with only ∼1% error in the free magnetic energy, assuming the photospheric information is specified every 12 minutes. However, for rapidly evolving features, under-sampling of the dynamics at this cadence leads to a strobe effect, generating large electric currents and incorrect coronal morphology and energies. We derive a sampling condition for the driving cadence based on the evolution of these small-scale features, and show that higher-cadence driving can lead to acceptable errors. Future work will investigate the source of errors associated with deriving plasma variables from the photospheric magnetograms as well as other sources of errors, such as reduced resolution, instrument bias, and noise.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Reduced MHD in Astrophysical Applications: Two-dimensional or Three-dimensional?

Originally proposed as an efficient approach to computation of nonlinear dynamics in tokamak fusion research devices, reduced magnetohydrodynamics (RMHD) has subsequently found application in studies of coronal heating, flux tube dynamics, charged particle transport, and, in general, as an approximation to describe plasma turbulence in space physics and astrophysics. Given the diverse set of derivations available in the literature, there has emerged some level of discussion and a lack of consensus regarding the completeness of RMHD as a turbulence model, and its applicability in contexts such as the solar wind. Some of the key issues in this discussion are examined here, emphasizing that RMHD is properly neither 2D nor fully 3D, being rather an incomplete representation that enforces at least one family of extraneous conservation laws.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

MHD-kinetic transition in imbalanced Alfvénic turbulence

Alfvénic turbulence in space is usually imbalanced: amplitudes of waves propagating parallel and anti-parallel to the mean magnetic field B{sub 0} are unequal. It is commonly accepted that the turbulence is driven by (counter-)collisions between these counter-propagating wave fractions. Contrary to this, we found a new ion-scale dynamical range of the turbulence established by (co-)collisions among waves co-propagating in the same direction along B{sub 0}. Co-collisions become stronger than counter-collisions and produce steep non-universal spectra above certain wavenumbers dependent on the imbalance. Spectral indexes of the strong turbulence vary around ≳−3, such that steeper spectra follow larger imbalances. Intermittency steepens the −3 spectra further, up to −3.7. Our theoretical predictions are compatible with steep variable spectra observed in the solar wind at ion kinetic scales, but further verifications are needed by correlating observed spectra with measured imbalances.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

General Relativistic Radiation MHD Simulations of Supercritical Accretion onto a Magnetized Neutron Star: Modeling of Ultraluminous X-Ray Pulsars

By performing 2.5-dimensional general relativistic radiation magnetohydrodynamic simulations, we demonstrate supercritical accretion onto a non-rotating, magnetized neutron star, where the magnetic field strength of dipole fields is 10{sup 10} G on the star surface. We found the supercritical accretion flow consists of two parts: the accretion columns and the truncated accretion disk. The supercritical accretion disk, which appears far from the neutron star, is truncated at around ≃3 R {sub *} ( R {sub *} = 10{sup 6} cm is the neutron star radius), where the magnetic pressure via the dipole magnetic fields balances with the radiation pressure of the disks. The angular momentum of the disk around the truncation radius is effectively transported inward through magnetic torque by dipole fields, inducing the spin up of a neutron star. The evaluated spin-up rate, ∼−10{sup −11} s s{sup −1}, is consistent with the recent observations of the ultraluminous X-ray pulsars. Within the truncation radius, the gas falls onto a neutron star along the dipole fields, which results in a formation of accretion columns onto the northern and southern hemispheres. The net accretion rate and the luminosity of the column are ≃66 L {sub Edd}/ c {sup 2} and ≲10 L {sub Edd}, where L {sub Edd} is the Eddington luminosity and c is the light speed. Our simulations support a hypothesis whereby the ultraluminous X-ray pulsars are powered by the supercritical accretion onto the magnetized neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗