Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD simulation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Impact of Substorm Time O+ Outflow on Ring Current Enhancement

Energetic O+ ions (tens of keV) rapidly increase in the inner magnetosphere and contribute significantly to the ring current during substorms. Previously, two source regions of the energetic O+ ions have been proposed. The first one is the dayside polar region. Ions from the dayside polar region are transported to the lobe; then they are injected to the nightside plasma sheet during substorm expansion phase. The second one is the nightside aurora region. After the substorm onset, energetic O+ ions are extracted from the ionosphere with the auroral acceleration processes, and the O+ ions are directly supplied to the nightside plasma sheet. We investigated the relative importance of these two regions on supplying the energetic O+ ions in the inner magnetosphere. We performed a test particle simulation in global MHD electromagnetic fields. We obtained the following results. (1) During the substorm growth phase, O+ ions at tens of eV are extracted from the dayside polar region, resulting in the enhancement of the warm O+ ions (hundreds of eV) in the lobe. After the substorm onset, the warm O+ ions are nonadiabatically accelerated to tens of keV and injected to the inner magnetosphere. These O+ ions contribute to most of the O+ ring current. (2) O+ ions at less than a few keVs are supplied from the nightside aurora region to the plasma sheet. However, their contribution to the O+ ring current remains small. From the results, we concluded that the main source of the energetic O+ ions is the dayside polar region.

outflow↗
NASA NTRS2011

Broken Ergodicity in MHD Turbulence in a Spherical Domain

Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.

Shebalin, John V.↗
NASA NTRS2003

Fundamentals of the Plasma Sail Concept: MHD and Kinetic Studies

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000] predicts that a 15-km standoff distance (or 20-km cross-sectional dimension) of the magnetic bubble will provide for sufficient momentum transfer from the solar wind to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Such velocities will enable travel out of the solar system in period of about seven years-almost an order of magnitude improvement over present chemical based propulsion systems. However, for the parameters of the simulation of Winglee et al. [2000], a fluid model for the interaction of M2P2 with the solar wind is not valid. It is assumed in the MHD fluid model, normally applied to planetary magnetospheres, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the solar wind. In the case of M2P2, the size of the magnetic bubble is actually less than or, comparable to, the scale of these characteristic parameters. Therefore, a kinetic approach, which addresses the small-scale physical mechanisms, must be used. We have adopted a two-component approach to determining a preliminary estimate of the momentum transfer to the plasma sail. The first component is a self-consistent MHD simulation of the small-scale expansion phase of the magnetic bubble. The fluid treatment is valid to roughly 5 km from the source and the steady-state MHD solution at the 5 km boundary was then used as initial conditions for the hybrid simulation. The hybrid simulations showed that the momentum transfer to the innermost regions of the plasma sail is negligible.

Khazanov, G.↗
NASA NTRS1987

Advances in magnetospheric physics by numerical simulations - A critical analysis by an outsider

This is a selective, somewhat editorialized assessment, based on a literature survey, of computer simulation as it exists today in magnetospheric physics. Both large (MHD) and small scale dynamical simulations are described and considered from the perspective of what they are trying to do and with what success. Several specific problem areas where simulations are being carried out are called for commentary: global magnetospheric structure, magnetotail and magnetopause reconnection, Kelvin-Helmholtz instability, hydrodynamic expansion of ionospheric ions, electric double layers, heavy ion heating, and auroral kilometric radiation. It is concluded that simulation is a necessary tool for understanding magnetospheric physics and that significant progress has been made in simulation development. However, results should be evaluated knowing that many factors, some real and physical, others structural, may contribute to such output. A plea is made for greater intercalibration among different simulators working in parallel areas, so that facts can be distinguished from artifacts.

Birmingham, Thomas J.↗
NASA NTRS2012

Statistical Theory of the Ideal MHD Geodynamo

A statistical theory of geodynamo action is developed, using a mathematical model of the geodynamo as a rotating outer core containing an ideal (i.e., no dissipation), incompressible, turbulent, convecting magnetofluid. On the concentric inner and outer spherical bounding surfaces the normal components of the velocity, magnetic field, vorticity and electric current are zero, as is the temperature fluctuation. This allows the use of a set of Galerkin expansion functions that are common to both velocity and magnetic field, as well as vorticity, current and the temperature fluctuation. The resulting dynamical system, based on the Boussinesq form of the magnetohydrodynamic (MHD) equations, represents MHD turbulence in a spherical domain. These basic equations (minus the temperature equation) and boundary conditions have been used previously in numerical simulations of forced, decaying MHD turbulence inside a sphere [1,2]. Here, the ideal case is studied through statistical analysis and leads to a prediction that an ideal coherent structure will be found in the form of a large-scale quasistationary magnetic field that results from broken ergodicity, an effect that has been previously studied both analytically and numerically for homogeneous MHD turbulence [3,4]. The axial dipole component becomes prominent when there is a relatively large magnetic helicity (proportional to the global correlation of magnetic vector potential and magnetic field) and a stationary, nonzero cross helicity (proportional to the global correlation of velocity and magnetic field). The expected angle of the dipole moment vector with respect to the rotation axis is found to decrease to a minimum as the average cross helicity increases for a fixed value of magnetic helicity and then to increase again when average cross helicity approaches its maximum possible value. Only a relatively small value of cross helicity is needed to produce a dipole moment vector that is aligned at approx.10deg with the rotation axis.

Shebalin, J. V.↗
NASA NTRS2014

Magnetohydrodynamic Turbulence and the Geodynamo

The ARES Directorate at JSC has researched the physical processes that create planetary magnetic fields through dynamo action since 2007. The "dynamo problem" has existed since 1600, when William Gilbert, physician to Queen Elizabeth I, recognized that the Earth was a giant magnet. In 1919, Joseph Larmor proposed that solar (and by implication, planetary) magnetism was due to magnetohydrodynamics (MHD), but full acceptance did not occur until Glatzmaier and Roberts solved the MHD equations numerically and simulated a geomagnetic reversal in 1995. JSC research produced a unique theoretical model in 2012 that provided a novel explanation of these physical observations and computational results as an essential manifestation of broken ergodicity in MHD turbulence. Research is ongoing, and future work is aimed at understanding quantitative details of magnetic dipole alignment in the Earth as well as in Mercury, Jupiter and its moon Ganymede, Saturn, Uranus, Neptune, and the Sun and other stars.

Shebalin, John V.↗
NASA NTRS1995

Simulating the interaction of the heliosphere with the local interstellar medium: MHD results from a finite volume approach, first bidimensional results

We have designed a full compressible MHD code working on unstructured meshes in order to be able to compute accurately sharp structures embedded in large scale simulations. The code is based on a finite volume method making use of a kinetic flux splitting. A bidimensional version of the code has been used to simulate the interaction of a moving interstellar medium, magnetized or unmagnetized with a rotating and magnetized heliopspheric plasma source. Being aware that these computations are not realistic due to the restriction to two dimensions, we present it to demonstrate the ability of this new code to handle this problem. An axisymetric version, now under development, will be operational in a few months. Ultimately we plan to run a full 3d version.

Chanteur, G.↗
NASA NTRS1984

Externally driven magnetic reconnection

A model is presented in which externaly driven reconnection is simulated by solving the MHD equations in an initially plane current sheet. Magnetic reconnection in the earth's magnetotail is widely believed to be the direct cause of magnetospheric substorms. Both 2-D and 3-D versions of the model have been developed. It is postulated that connection in the tail is triggered by a local compression of the plasma sheet which results from an invasion of the solar wind into the magnetotail. Thus, the simulation is started by introducing flow from the lobes normal to the plasma sheet. When resistivity is generated in a local region of the neutral sheet, reconnection develops and magnetic energy is converted into plasma bulk flow. Although the driven reconnection model is highly simplified, it can aid in understanding many features of substorms in the tail; in particular, results show that rapid flows both earthward and tailward of the neutral line and the nightside substorm current system are natural consequences of driven magnetic reconnection.

Walker, R. J.↗
NASA NTRS2022

Electromagnetic Time and Wavelength Domain Predictions for Accreting Binary Black Holes

Accreting massive binary black holes in gaseous environments are keymulti-messenger sources for the LISA mission. We will present some ofthe latest predictions of the observational signatures of accreting binary black holes across the electromagnetic spectrum from the fieldand provide some of our own. In particular, we will provide a briefsummary of the progress made in our group to understand these systemstheoretically using high-performance simulations of generalrelativistic magnetohydrodynamics (MHD) and dynamic spacetimes.Simulating these systems realistically is challenging asradiation-coupled MHD must be considered over large dynamic ranges inspace and time. Constraining our view to the circumbinary disk regionfirst, we will show how the evolution duration, accretion disk size,mass ratio, and thermodynamics model affect the structure andvariability of the accretion flow. We particularly emphasize how theseparameters influence the over-density feature, or "lump", which orbitsthe binary near the edge of the cavity, since it is responsible formost of the electromagnetic emission's variability---a key signatureof a system being a binary. Extending to smaller length scales, wewill report on simulations following accretion all the down to theevent horizons so that we may begin to investigate how black hole spinaffects mini-disk dynamics, accretion rate, and jet power. We willalso report on our radiative transfer predictions for how black holespin alters the spectral energy distribution of these systems.Wrapping up, we will itemize the remaining key theoreticaltechnologies and directions needed to be accomplished to betterunderstand these sources, and assess strategies for searching for themin the coming decade.

Scott C Noble↗
NASA NTRS1996

Drag reduction in turbulent MHD pipe flows

This is a preliminary study devoted to verifying whether or not direct simulations of turbulent Magneto-Hydro-Dynamic (MHD) flows in liquid metals reproduce experimental observations of drag reduction. Two different cases have been simulated by a finite difference scheme which is second order accurate in space and time. In the first case, an external azimuthal magnetic field is imposed. In this case, the magnetic field acts on the mean axial velocity and complete laminarization of the flow at N(sub a) = 30 has been achieved. In the second case, an axial magnetic field is imposed which affects only fluctuating velocities, and thus the action is less efficient. This second case is more practical, but comparison between numerical and experimental results is only qualitative.

Orlandi, P.↗
NASA NTRS1988

A Three-Dimensional, Time-Dependent Numerical Modeling of Super-Sonic, Super-Alfvenic MHD Flow

A three-dimensional, time-dependent, MHD disturbance propagating through a steady inhomogeneous MHD flow is analyzed by a numerical model based on the Lax-Wendroff finite difference method. The initial steady MHD flow is assumed to be supersonic and super-Alfvenic. The background steady MHD flow is obtained by dropping the time-dependent terms in the unsteady MHD equations. A transient MHD disturbance corresponding to a MHD shock is generated at the lower radial boundary surface. The resulting development and propagation of 3-D MHD disturbances interacting with the background steady MHD flow are numerically simulated in this first application to the field of interplanetary physics. Detailed descriptions of the mathematical formulation and numerical methodology are presented. Example computations are performed to demonstrate the validity of the formulation and the capabilities of the model.

Han, S. M.↗
NASA NTRS1987

Evolution of polarization in localized nonlinear Alfven waves

A polarization change of localized and finite-amplitude Alfven waves is found by use of both MHD theory and numerical simulation. One of the transverse components of the circularly polarized Alfven wave is stable as it propagates, but the other component is unstable under both self-focusing and diffraction effects. Consequently, the wave changes its polarization from circular to linear.

Hoshino, Masahiro↗
NASA NTRS1989

Waves in low-beta plasmas - Slow shocks

Results from wave theory and numerical simulation of the nonlinear MHD equations are used to study the response of a conducting fluid containing an embedded magnetic field with beta less than 1 to the sudden injection of material along the field lines. It is shown that the injection produces slow shocks with configurations which are concave toward the ejecta driver. Fast-mode waves which have not steepened into the shock precede the slow shock and alter the ambient medium. When beta equals 0.1, the fast mode becomes a transverse wave for parallel propagation, while the slow wave approaches a longitudinal, or sound, wave.

Steinolfson, R. S.↗
NASA NTRS1992

Solar wind-magnetosphere interaction as simulated by a 3-D EM particle code

We present here our first results of simulating the solar wind-magnetosphere interaction with a new three-dimensional electromagnetic particle code. Hitherto such global simulations were done with MHD codes while lower-dimensional particle or hybrid codes served to account for microscopic processes and such transport parameters as have to be introduced ad hoc in MHD. Our kinetic model attempts to combine the macroscopic and microscopic tasks. It relies only on the Maxwell curl equation and the Lorentz equation for particles, which are ideally suited for computers. The preliminary results shown here are for an unmagnetized solar wind plasma streaming past a dipolar magnetic field. The results show the formation of a bow shock and a magnetotail, the penetration of energetic particles into cusp and radiation belt regions, and dawn-dusk asymmetries.

Buneman, Oscar↗
NASA NTRS1992

A two-dimensional MHD global coronal model - Steady-state streamers

A 2D, time-dependent, numerical, MHD model for the simulation of coronal streamers from the solar surface to 15 solar is presented. Three examples are given; for dipole, quadrupole and hexapole (Legendre polynomials P1, P2, and P3) initial field topologies. The computed properties are density, temperature, velocity, and magnetic field. The calculation is set up as an initial-boundary value problem wherein a relaxation in time produces the steady state solution. In addition to the properties of the solutions, their accuracy is discussed. Besides solutions for dipole, quadrupole, and hexapole geometries, the model use of realistic values for the density and Alfven speed while still meeting the requirement that the flow speed be super-Alfvenic at the outer boundary by extending the outer boundary to 15 solar radii.

Wang, A.-H.↗
NASA NTRSDate not supplied

Photometric Properties of Quiet Sun Observed With the DKIST and HINODE

The variability of irradiance is primarily influenced by surface magnetism, such as sunspots, faculae, and network structures. Each of these features contributes differently to the overall variability. However, the impact of small magnetic elements, which are commonly observed in high resolution images of the solar photosphere but remain unresolved in full-disk images used for irradiance modeling, is still a topic of debate. Understanding the influence of these small magnetic elements to the overall brightness, particularly in quiet regions, is fundamental to understanding irradiance fluctuations over decadal and longer temporal scales. We investigate the brightness of small-size magnetic elements using very high spatial resolution spectro-polarimetric observations acquired at different spectral ranges with the National Science Foundations Daniel K. Inouye Solar Telescope (DKIST) and with HINODE/SOT. Our results are compared to 3D MHD Small Solar Dynamo simulations obtained with the MURaM code.

SMD↗
NASA NTRSDate not supplied

An Update on the Storm Time O+ Ring Current Imaging Evolution (STORIE) Mission

The Storm Time O+ Ring Current Imaging Evolution (STORIE) mission launching in 2025 will provide continuous remote sensing of the ring current using Energetic Neutral Atoms (ENAs). The STORIE payload will be hosted on the International Space Station (ISS) where it will image the ring current ENA from the inside out. The mission will address compelling issues of ring current science including issues of buildup and decay, as well as the role of heavy ions. This presentation will overview the mission, and update the community on the current status. We will moreover discuss the science objectives, and elaborate on the advantages of ENA imaging from the unique vantage point of the ISS. We will then present the result of numerical models illustrating the expected ENA signal using coupled global MHD and ring current simulations during a geomagnetic event. Finally, we will discuss the initial approach to inverting the ENA signal to obtain the ring current ion flux.

Alex Glocer↗
NASA NTRS2008

Recovering Photospheric Velocities from Vector Magnetograms by using a Three-Dimensional, Fully Magnetohydrodynamic (MHD) Model

We introduce a numerical simulation method for recovering the photospheric velocity field from the vector magnetograms. The traditional method is local correlation tracking (LCT) which is based on measuring the relative displacements of features in blocks of pixels between successive white-light images or magnetograms. Within this method, there are a variety of implementations. One of recently developed implementations is induction local correlation tracking (ILCT) as described by Welsch et al. (2004). They employ the normal component of magnetic induction equation as a constraint to assure consistent solutions. Our numerical method uses the fully three-dimensional MHD equations to recover the photospheric velocity field with individual vector magnetograms. We compare our method to the ILCT method using NOAA AR8210 as an example. The differences and similarities are discussed in detail.

Wang, A. H.↗