SEARCH · Search NASA

Results for “Numerical magnetohydrodynamics”

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.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1997

The Ambient and Perturbed Solar Wind: From the Sun to 1 AU

The overall objective of the proposed research was to use numerical solutions of the magnetohydrodynamic (MHD) equations along with comparisons of the computed results with observations to study the following topics: (1) ambient solar wind solutions that extend from the solar surface to 1 astronomical unit (AU), contain closed magnetic structures near the Sun, and are consistent with observed values; (2) magnetic and plasma structures in coronal mass ejections (CMES) as they propagate to the interplanetary medium; (3) relation of MHD shocks to CMEs in the interplanetary medium; (4) interaction of MHD shocks with structures (such as other shocks, corotating interaction regions, current sheets) in the interplanetary plasma; and (5) simulations of observed interplanetary structures. A planned close collaboration with data analysts served to make the model more relevant to the data. The outcome of this research program is an improved understanding of the physical processes occurring in solar-generated disturbances in the interplanetary medium between the Sun and 1 AU.

Steinolfson, R. S.↗

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.↗

NASA NTRS2006

3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

Mizuno, Yosuke↗

NASA NTRS2007

Development of Numerical Tools for the Investigation of Plasma Detachment from Magnetic Nozzles

A multidimensional numerical simulation framework aimed at investigating the process of plasma detachment from a magnetic nozzle is introduced. An existing numerical code based on a magnetohydrodynamic formulation of the plasma flow equations that accounts for various dispersive and dissipative processes in plasmas was significantly enhanced to allow for the modeling of axisymmetric domains containing three.dimensiunai momentum and magnetic flux vectors. A separate magnetostatic solver was used to simulate the applied magnetic field topologies found in various nozzle experiments. Numerical results from a magnetic diffusion test problem in which all three components of the magnetic field were present exhibit excellent quantitative agreement with the analytical solution, and the lack of numerical instabilities due to fluctuations in the value of del(raised dot)B indicate that the conservative MHD framework with dissipative effects is well-suited for multi-dimensional analysis of magnetic nozzles. Further studies will focus on modeling literature experiments both for the purpose of code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗

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 NTRS2002

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Sjoegreen, Bjoern↗

NASA NTRS1988

Turbulent magnetohydrodynamic density fluctuations

A spectral-method numerical code is used to compute mass-density fluctuation spectra in turbulent magnetofluids. The computations are used to test and extend the analytical theory of density variations in slightly compressible magnetofluids given by Montgomery, et al. (1987) and used to infer inertial-range density-fluctuation spectra for the nearby interstellar medium and solar wind. A local equation of state is assumed, relating density to pressure. Constant, scalar resistivities and viscosities are used. In the limit of low Mach numbers and high mechanical-to-magnetic pressure ratios, the fit of the computations to the analytical theory is seen to be close.

Shebalin, John V.↗

NASA NTRS1988

Helical magnetohydrodynamic turbulence and the coronal heating problem

Numerical simulations are used to investigate the relaxation of an unconfined, helically turbulent, fully three-dimensional magnetofluid, with conditions similar to those which are thought to result in the heating of the solar corona. In these simulations, the system evolves through a succession of force free states. After a relatively quiescent period of Ohmic decay, a phase of accelerated magnetic energy dissipation occurs. Some magnetic energy is transformed into kinetic energy, and the magnitude of entrophy created is a nontrival fraction of the mean square electric current. Concentrated vorticity structures are seen to play almost as important a role as electric current sheets in the heating process. Coincident with this accelerated dissipation process, a reorganization of the magnetic fields occurs, with transfer of magnetic energy to both shorter and longer wavelength modes than are initially present. The ratio of the magnetic field to the electric current density, alpha does not in general tend to assume a constant value in the force free regions during the evolution of the magnetofluid.

Dahlburg, R. B.↗

NASA NTRS1976

Nonlinear effects for the Taylor column for a hemisphere

When a flow is forced past an obstacle in a rapidly rotating fluid, a Taylor column forms. This is defined by a set of vertical detached shear layers circumscribing the obstacle which provide the smooth transition from an external inviscid potential flow to a stagnant core above the obstacle. For a hemispherical object, the main adjustment takes place in an external E to the 1/4 power layer and an internal E to the 2/7 power layer; here, the nonlinear flow in these layers is investigated. The problem in the E to the 1/4 power layer is identical to a problem occurring in magnetohydrodynamic flow; in addition, some features of the magnetohydrodynamic problem have been resolved. Numerical solutions are obtained for the steady nonlinear external E to the 1/4 power layer flow up to the point where unsteady flow separation from the Taylor column is imminent. The response of the internal E to the 2/7 power layer to the flow in the E to the 1/4 power layer is calculated, and the results suggest that the internal shear layer is unlikely to play any significant role in the separation process

Crisalli, A. J.↗

NASA NTRS1978

A three-dimensional model of corotating streams in the solar wind. I - Theoretical foundations

The paper is concerned with the development of the theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 Rs where the flow is supersonic) is approximately described by the nonlinear, single-fluid, polytropic magnetohydrodynamic or hydrodynamic equations. Efficient numerical techniques are outlined for solving this complex system of coupled, hyperbolic partial differential equations. The present formulation is inviscid and nonmagnetic, but the methods used allow for the potential inclusion of both features with only modest modifications. A simple, highly idealized hydrodynamic model stream is examined to illustrate the fundamental processes involved in the three-dimensional dynamics of stream evolution. It is found that spatial variations in the rotational stream interaction mechanism produce small nonradial flows on a global scale that lead to the transport of mass, energy, and momentum away from regions of relative compression and into regions of relative rarefaction. Comparison with simpler models demonstrates the essential nonlinear, multidimensional nature of the interplanetary dynamics.

Pizzo, V.↗

NASA NTRS1982

Anisotropy in MHD turbulence due to a mean magnetic field

The development of anisotropy in an initially isotropic spectrum is studied numerically for two-dimensional magnetohydrodynamic turbulence. The anisotropy develops due to the combined effects of an externally imposed dc magnetic field and viscous and resistive dissipation at high wave numbers. The effect is most pronounced at high mechanical and magnetic Reynolds numbers. The anisotropy is greater at the higher wave numbers.

Shebalin, J. V.↗

NASA NTRS1983

Anisotropy in MHD turbulence due to a mean magnetic field

Shebalin, J. V.↗

NASA NTRS1990

A helically distorted MHD flux rope model

A flux rope model is proposed which has a variable degree of helical distortion from axisymmetry. The basis for this suggestion is a series of numerical and analytical investigations of magnetohydrodynamic states which result when an axial electric current is directed down on dc magnetic field. The helically distorted states involve a flow velocity and seem to be favored because of their lower rate of energy dissipation. Emphasis is on the magnetometer and particle energy analyzer traces that might be characteristic of such flux ropes. It is shown that even a fractionally small helical distortion may considerably alter the traces in minimum-variance coordinates. In short, what may be fairly common MHD processes can render a flux rope almost unrecognizable under standard diagnostics, even if the departures from axisymmetry are not great.

Theobald, Michael L.↗

NASA NTRS1993

An MHD study of the interaction between the solar wind and the interstellar medium

Steinolfson, R. S.↗

NASA NTRS1996

An MHD Study of the Interaction Between the Solar Wind and the Interstellar Medium

The overall objective of this research program is to obtain a better understanding of the interaction between the solar wind and the interstellar medium through the use of numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations. The simulated results have been compared with observations where possible and with the results from previous analytic and numerical studies. The primary accomplishment of this project has been the development of codes for 2-D models in both spherical and cylindrical coordinates and the application of the codes to the solar wind/interstellar medium interaction. Computations have been carried out for both a relatively simple gas-dynamic interaction and a flow-aligned interstellar magnetic field. The results have been shown to compare favorably with models that use more approximations and to modify and extend the previous results as would be expected. The simulations have also been used along with a data analysis study to provide a quantitative estimate of the distance to the termination and bow shocks. Some of the specific topics that have been studied are: (1) gas dynamic models of the solar wind/interstellar medium interaction, (2) termination shock response to large-scale solar wind fluctuations, and (3) distances to the termination shock and heliopause. The main results from each of these studies are summarized. The results were published in three papers which are included as attachments.

Steinolfson, R. S.↗

NASA NTRS1987

Magnetohydrodynamic (MHD) analyses of various forms of activity and their propagation through helio spheric space

Theoretical and numerical modeling of solar activity and its effects on the solar atmosphere within the context of magnetohydrodynamics were examined. Specifically, the scientific objectives were concerned with the physical mechanisms for the flare energy build-up and subsequent release. In addition, transport of this energy to the corona and solar wind was also investigated. Well-posed, physically self-consistent, numerical simulation models that are based upon magnetohydrodynamics were sought. A systematic investigation of the basic processes that determine the macroscopic dynamic behavior of solar and heliospheric phenomena was conducted. A total of twenty-three articles were accepted and published in major journals. The major achievements are summarized.

Wu, S. T.↗

NASA NTRS1991

Nonlinear magnetohydrodynamics by Galerkin-method computation

A fully spectral numerical code is used to explore the properties of voltage-driven dissipative magnetofluids inside a periodic cylinder with circular cross section. The trial functions are orthonormal eigenfunctions of the curl (Chandrasekhar-Kendall functions). Transitions are observed from axisymmetric resistive equilibria without flow to helically deformed laminar states with flow, and between pairs of helical laminar states with different pairs of poloidal and toroidal m and n numbers. States of minimum energy dissipation rate seem to be preferred. At high values of the pinch ratio, fully developed magnetohydrodynamic turbulence is observed.

Shan, Xiaowen↗

NASA NTRS1984

Nonlinear evolution of the resistive tearing mode

The nonlinear behavior of the tearing instability is investigated with the aid of numerical solutions of the resistive incompressible magnetohydrodynamic equations. This study is directed toward solar applications, such as flares. Simulations were carried out for values of the Lundquist number S from 100 to 1,000,000 and wavelength parameter alpha from 0.042 to 0.5. The influence of linear conditions, such as the constancy of psi, on nonlinear growth characteristics is investigated. For high S and low alpha, secondary-flow vortices, opposite in direction to the linear vortices, were found to generate a new magnetic island centered at the initial X point. The nonlinear spatial distributions of the physical variables differ greatly from the linear behavior, and growth slows considerably once nonlinear effects become important.

Steinolfson, R. S.↗