Search NASA⌕ Search

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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Numerical Modeling of Three-Dimensional Confined Flows

A three dimensional confined flow model is presented. The flow field is computed by calculating velocity and enthalpy along a set of streamlines. The finite difference equations are obtained by applying conservation principles to streamtubes constructed around the chosen streamlines. With appropriate substitutions for the body force terms, the approach computes three dimensional magnetohydrodynamic channel flows. A listing of a computer code, based on this approach is presented in FORTRAN IV language. The code computes three dimensional compressible viscous flow through a rectangular duct, with the duct cross section specified along the axis.

Greywall, M. S.↗
NASA NTRS1984

Numerical Simulations of Convectively Excited Gravity Waves

Magneto-convection and gravity waves are numerically simulated with a nonlinear, three-dimensional, time-dependent model of a stratified, rotating, spherical fluid shell heated from below. A Solar-like reference state is specified while global velocity, magnetic field, and thermodynamic perturbations are computed from the anelastic magnetohydrodynamic equations. Convective overshooting from the upper (superadiabatic) part of the shell excites gravity waves in the lower (subadiabatic) part. Due to differential rotation and Coriolis forces, convective cell patterns propagate eastward with a latitudinally dependent phase velocity. The structure of the excited wave motions in the stable region is more time-dependent than of the convective motions above. The magnetic field tends to be concentrated over giant-cell downdrafts in the convective zone but is affected very little by the wave motion in the stable region.

Glatzmaier, G. A.↗
NASA NTRS2011

CUDA Simulation of Homogeneous, Incompressible Turbulence

We discuss very fast Compute Unified Device Architecture (CUDA) simulations of ideal homogeneous incompressible turbulence based on Fourier models. These models have associated statistical theories that predict that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. Prior numerical simulations have shown that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We review the theoretical basis of this "broken ergodicity" as applied to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence. Our new simulations examine the phenomenon of broken ergodicity through very long time and large grid size runs performed on a state-of-the-art CUDA platform. Results comparing various CUDA hardware configurations and grid sizes are discussed. NS and MHD results are compared.

Morin, Lee↗
NASA NTRS1989

Magnetohydrodynamic waves and particle acceleration in a coronal loop

A model is presented for the acceleration of electrons in a flaring solar coronal loop. It is suggested that the primary energy release event in the flare establishes a spectrum of MHD waves within the loop which accelerate electrons by wave-particle resonant interactions. The nature of the MHD oscillations are described and numerical and analytic solution of the dispersion relation are examined. It is found that the behavior of a particle depends on whether or not it is in resonance with the wave spectrum. Resonant particles follow chaotic trajectories and can diffuse in velocity space between wave effective potentials, thereby gaining energy.

De La Beaujardiere, Jean-Francois↗
NASA NTRS1974

Solar Physics - Plasma Physics Workshop

A summary of the proceedings of a conference whose purpose was to explore plasma physics problems which arise in the study of solar physics is provided. Sessions were concerned with specific questions including the following: (1) whether the solar plasma is thermal or non-themal; (2) what spectroscopic data is required; (3) what types of magnetic field structures exist; (4) whether magnetohydrodynamic instabilities occur; (5) whether resistive or non-magnetohydrodynamic instabilities occur; (6) what mechanisms of particle acceleration have been proposed; and (7) what information is available concerning shock waves. Very few questions were answered categorically but, for each question, there was discussion concerning the observational evidence, theoretical analyses, and existing or potential laboratory and numerical experiments.

Baum, P. J.↗
NASA NTRS1994

Stochastic acceleration of energetic ions in Jupiter's magnetosphere

An equation governing the combined radial diffusion and stochastic acceleration of super-Alfvenic ions by magnetohydrodynamic waves in Jupiter's outer magnetosphere is derived. The formulation is based upon a total energy invariant of the adiabatic transport which applies to an isotropic distribution undergoing rapid pitch angle scattering by waves. An analytic solution to the double diffusion equation is obtained and numerical results are presented for two models of ion injection. The first model assumes S(+) and O(+) are injected throughout a broad region of space through photoionization of Jupiter's magnetospheric neutral wind and obtain an initial energy corresponding to the local corotation energy after pickup by the planetary magnetic field. The second model assumes a monoenergetic distribution of energetic protons is implanted in the middle magnetosphere by the action of field-aligned potential drops in Jupiter's auroral ionosphere. For both light and heavy ions the injection process creates a seed population of particles which are further accelerated nonadiabatically by the MHD waves and adiabatically through radial diffusion. A comparison of the theoretical results with a recent data analysis of Voyager low-energy charged particle measurements is made with very good agreement, thus providing a rigorous quantitative account of and definitive explanation for the high-energy ion component of Jupiter's magnetosphere.

Barbosa, D. D.↗
NASA NTRS1982

Interaction of small-amplitude fluctuations with a strong magnetohydrodynamic shock

The interaction of small fluctuations with a strong magnetohydrodynamic shock is considered theoretically in the most general situation using an ideal model. An analytical treatment of the cutoff conditions for penetration of the upstream fluctuation through the bow shock is presented. The situation in which there exist unique solutions of the interaction is investigated quantitatively. Numerical examples of the interaction of waves in the solar wind with the terrestrial bow shock are given to illustrate the refracted waves to be found in the earth's magnetosheath.

Zhuang, H.-C.↗
NASA NTRS1992

The theory of sunspots

This review covers the present state of our theoretical understanding of the physics of sunspots, along with the principal observational results that need to be explained. The topics covered range from the detailed structure of an individual sunspot to the broad connection between sunspots and the global solar magnetic field and the solar cycle. Our aim is to give a critical discussion of the theoretical ideas and models without presenting mathematical details. After outlining the historical development of the basic concepts associated with the magnetohydrodynamic theory of sunspots, we discuss recent treatments of their properties and structure, placing special emphasis on developments that have occurred within the last ten years. There have been remarkable improvements in the theoretical modelling of sunspots, led by new ideas and by more elaborate and realistic numerical simulations. At the same time, new observations have raised new theoretical questions or caused old ones to be reconsidered. In particular, measurements of oscillations in and around sunspots have opened up the new field of sunspot seismology, while recent high-resolution observations have forced us to rethink the structure of a sunspot penumbra.

Thomas, John H.↗
NASA NTRS1992

Model of a rotating magnetic cloud

The possibility that magnetic clouds rotate while they propagate antisunward was investigated. Magnetic clouds are modeled as magnetic flux ropes which rotate rigidly about the axis of symmetry. An ideal magnetohydrodynamic model, in which the evolution of the magnetic structure is related to the time evolution of the angular frequency, is developed. A class of 'separable' magnetic fields is employed to reduce the problem to a nonlinear ordinary differential equation for the evolution function, and it is solved numerically. The corresponding effective potential gives rise to two modes of evolution--expansion and oscillation--depending on the energy and on the value of a dimensionless parameter, k. Parameter k depends on the gas pressure, the ratio of the magnetic field components, and the frequency of rotation. There is a critical value of k, k(sub c), above which the oscillatory regime disappears and the flux rope invariably expands, regardless of the energy. Below k(sub c) the energy determines whether the configuration is confined or unbounded. Rotation always helps expansion by lowering the potential barrier. A data example was studied and features which are interpreted as signatures of rotation are presented. The angular speed is comparable to the Alfven speed, and the core of the rotating cloud completes on average one full revolution every three days at 1 AU. The parameter k is calculated from observations, and it is found to be close to, but below, critical. Only three out of the nine clouds examined showed signatures of rotation. Theoretical analysis suggests that close to the Sun rotation effects may play a more important role in the evolution of magnetic clouds than 1 AU.

Farrugia, C. J.↗
NASA NTRS1988

The role of magnetohydrodynamics in heliospheric space plasma physics research

Magnetohydrodynamics (MHD) is a fairly recent extension of the field of fluid mechanics. While much remains to be done, it has successfully been applied to the contemporary field of heliospheric space plasma research to evaluate the 'macroscopic picture' of some vital topics via the use of conducting fluid equations and numerical modeling and simulations. Some representative examples from solar and interplanetary physics are described to demonstrate that the continuum approach to global problems (while keeping in mind the assumptions and limitations therein) can be very successful in providing insight and large scale interpretations of otherwise intractable problems in space physics.

Dryer, Murray↗
NASA NTRS2007

GRMHD/RMHD Simulations and Stability of Magnetized Spine-Sheath Relativistic Jets

A new general relativistic magnetohydrodynamics (GRMHD ) code "RAISHIN" used to simulate jet generation by rotating and non-rotating black holes with a geometrically thin Keplarian accretion disk finds that the jet develops a spine-sheath structure in the rotating black hole case. Spine-sheath structure and strong magnetic fields significantly modify the Kelvin-Helmholtz (KH) velocity shear driven instability. The RAISHIN code has been used in its relativistic magnetohydrodynamic (RMHD) configuration to study the effects of strong magnetic fields and weakly relativistic sheath motion, cl2, on the KH instability associated with a relativistic, Y = 2.5, jet spine-sheath interaction. In the simulations sound speeds up to ? c/3 and Alfven wave speeds up to ? 0.56 c are considered. Numerical simulation results are compared to theoretical predictions from a new normal mode analysis of the RMHD equations. Increased stability of a weakly magnetized system resulting from c/2 sheath speeds and stabilization of a strongly magnetized system resulting from d 2 sheath speeds is found.

Hardee, Philip↗
NASA NTRS1990

A new class of solutions for interstellar magnetohydrodynamic shock waves

An analysis is presented of the equations of motion for steady MHD shock waves proopagating in interstellar clouds, for boundary conditions that preclude C shocks. In addition to J shocks, in which the neutral fluid component becomes subsonic at an adiabatic jump front, the equations admit a new class of solutions, called C-asterisk shocks, in which the transition to subsonic flow occurs continuously at a sonic point. Numerical methods are developed for computing the structure of J and C-asterisk shocks propagating in diffuse interstellar clouds. The effects of chemical, ionization, and recombination processes are included in this treatment. An alternative numerical method, which uses artificial viscosity to facilitate integration through sonic points, is analyzed and shown to be invalid. A set of exemplary solutions, computed for realistic shock parameters, shows that C-asterisk shocks occur for a broad range of conditions relevant to diffuse interstellar clouds.

Roberge, W. G.↗
NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗
NASA NTRS1988

Cosmic-ray-modified stellar winds. III - A numerical iterative approach

A numerical iterative method is used to determine the modification of a stellar wind flow with a termination shock by the galactic cosmic rays. A two-fluid model consisting of cosmic rays and thermal stellar wind gas is used in which the cosmic rays are coupled to the background flow via scattering with magnetohydrodynamic waves or irregularities. A polytropic model is used to describe the thermal stellar wind gas, and the cosmic-rays are modeled as a hot, low-density gas with negligible mass flux. The positive galactic cosmic-ray pressure gradient serves to brake the outflowing stellar wind gas, and the cosmic rays modify the location of the critical point of the wind, the location of the shock, the wind fluid velocity profile, and the thermal gas entropy constants on both sides of the shock. The transfer of energy to the cosmic rays results in an outward radial flux of cosmic-ray energy.

Ko, C. M.↗
NASA NTRS1977

Dynamic modeling of coronal and interplanetary responses to solar events

Recent progress in the dynamic modeling of responses of the corona and interplanetary medium to solar events (such as surges, eruptive prominences, flares, etc.) is reviewed. In particular, coronal transients and wave phenomena are discussed in some detail including pertinent mathematical requirements. Within the context of hydrodynamics and magnetohydrodynamics, a summary of both one- and two-dimensional time-dependent models is presented. A comparison of theoretical results with ground-based optical, radio, ATM Skylab observations, Pioneer 9 and Pioneer 10 solar wind observations is also presented. It is illustrated that: (1) substantial progress has been made in the theoretical (i.e., numerical) modeling of coronal and interplanetary responses within the last few years; (2) two-dimensional, time-dependent modeling is needed to examine the details of nonlinear wave coupling; and (3) theoretical results of modeling appear to reproduce physical consequences successfully.

Wu, S. T.↗
NASA NTRS1990

Investigation of turbulent processes in magnetospheric boundary layers

A self-consistent non-evolving two dimensional slab model of a viscous low-latitude boundary layer (LLBL) coupled to the ionosphere was developed by Phan, et al., (1989). Numerical results from the model and possible use of observations to determine the model parameters are discussed. The dynamical model developed by Lotko, et al., (1987) was used by Lotko and Shen (1991) to examine dynamical processes relevant to the LLBL with particular application to post-noon auroral shear layers. Initial results from a magnetohydrodynamic study of flank-side mangetopause boundary configuration are described. Effects of compressibility, scalar viscosity, and electrical resistivity are included in the MHD equations.

Lotko, William↗