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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Nonparametric solutions to the variational principle of ideal magnetohydrodynamics

In an effort to gain a better understanding of MHD equilibria in three dimensions, the lower dimensional cases are studied. The solution of the three-dimensional problem is based on the classical variational principle of ideal magnetohydrodynamics. The crucial assumption for the numerical method is the existence of a nested set of toroidal flux surfaces, which is then used as a coordinate. This paper studies the nonparametric solutions to this variational problem in those cases when the direct solution is known to have islands. A form of the variational principle for the slab geometry is described; the one-dimensional problem is analyzed; and asymptotic expansions and numerical solutions to the two-dimensional problem are discussed. An example is presented which shows that the assumption of nested flux surfaces need not rule out the occurrence of islands.

Betancourt, O. L.↗
NASA NTRS1994

The evolution of cosmic-ray-mediated magnetohydrodynamic shocks: A two-fluid approach

We study the shock structure and acceleration efficiency of cosmic-ray mediated Magnetohydrodynamic (MHD) shocks both analytically and numerically by using a two-fluid model. Our model includes the dynamical effect of magnetic fields and cosmic rays on a background thermal fluid. The steady state solution is derived by following the technique of Drury & Voelk (1981) and compared to numerical results. We explore the time evolution of plane-perpendicular, piston-driven shocks. From the results of analytical and numerical studies, we conclude that the mean magnetic field plays an important role in the structure and acceleration efficiency of cosmic-ray mediated MHD shocks. The acceleration of cosmic-ray particles becomes less efficient in the presence of strong magnetic pressure since the field makes the shock less compressive. This feature is more prominent at low Mach numbers than at high Mach numbers.

Jun, Byung-Il↗
NASA NTRS1994

Numerical simulation of large-scale field-aligned current generation from finite-amplitude magnetosonic waves

A two-dimensional numerical simulation of finite-amplitude magnetohydrodynamic (MHD) magnetosonic waves is performed under a finite-velocity background convection condition. Isothermal cases are considered for simplicity. External dissipation is introduced by assuming that the field-aligned currents are generated in proportion to the accumulated charges. The simulation results are as follows: Paired field-aligned currents are found from the simulated waves. The flow directions of these field-aligned currents depend on the angle between the background convection and the wave normal, and hence two pairs of field-aligned currents are found from a bowed wave if we look at the overall structure. The majority of these field-aligned currents are closed within each pair rather than between two wings. These features are not observed under slow background convection. The result could be applied to the cusp current system and the substorm current system.

Yamauchi, M.↗
NASA NTRS1982

The steady global corona

The formation of the steady coronal structure which consists of coronal streamers and holes is seen in terms of a model employing numerical solutions of time-dependent, dissipationless, magnetohydrodynamic equations of motion applicable to the meridional plane. A coronal streamer consists of closed magnetic field lines near the solar surface with overlying and adjacent open field lines, and the atmosphere, which is stationary in the closed region, flows outward in the open region, or coronal hole. The steady coronal structure is obtained by starting the numerical calculation with a state comprising a polytropic, hydrodynamical solution to the steady-state radial equation of motion coupled with a dipole magnetic field. Global coronal configurations are calculated for values of the plasma beta which vary from 0.1 to 100. It is found that the height and lateral extension of the closed region are only weakly dependent on the data.

Steinolfson, R. S.↗
NASA NTRS2004

Studying Turbulence Using Numerical Simulation Databases - X Proceedings of the 2004 Summer Program

This Proceedings volume contains 32 papers that span a wide range of topics that reflect the ubiquity of turbulence. The papers have been divided into six groups: 1) Solar Simulations; 2) Magnetohydrodynamics (MHD); 3) Large Eddy Simulation (LES) and Numerical Simulations; 4) Reynolds Averaged Navier Stokes (RANS) Modeling and Simulations; 5) Stability and Acoustics; 6) Combustion and Multi-Phase Flow.

Moin, Parviz↗
NASA NTRS1971

Development of MHD flow in the entrance region of a channel.

A general method has been developed for determining the developing magnetohydrodynamic flow in the entrance region of a parallel-plate channel. Closed form solutions are obtained. Numerical results are presented for the development of velocity distribution, pressure drop, and magnetohydrodynamic entrance length at different Hartmann numbers. At large Hartmann numbers, the magnetic field effect becomes dominant and the magnetohydrodynamic entrance length becomes negligible. Also, the pressure drop is enhanced by the magnetic field. Comparison of results with direct numerical analysis shows excellent agreement.-

Wu, S. T.↗
NASA NTRS1995

Numerical magetohydrodynamics in astronphysics: Algorithm and tests for one-dimensional flow`

We describe a numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference scheme on an Eulerian grid, called the total variation diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. We also describe a non-linear Riemann solver for ideal MHD, which includes rarefractions as well as shocks. The numerical code and the Riemann solver have been used to test each other. Extensive tests encompassing all the possible ideal MHD structures with planar symmetries (i.e., one-dimensional flows) are presented. These include those for which the field structure is two dimensional (i.e., those flows often called '1 + 1/2 dimensional') as well as those for which the magnetic field plane rotates (i.e., those flows often called '1 + 1/2 + 1/2 dimensional'). Results indicate that the code can resolve strong fast, slow, and magnetosonic shocks within two to four cells, but more cells are required if shocks become weak. With proper steepening, we could resolve rotational discontinuities within three to five cells. However, without successful implementation of steepening, contact discontinuities are resolved with approximately 10 cells and tangential discountinuities are resolved with approximately 15 cells. Out tests confirm that slow compound structures with tow-dimensional magnetic fields are composed of intermediate shocks (so-called 2-4 intermediate shocks) followed by slow rarefaction waves. Finally, tests demostrate that in two-dimensional magnetohydrodynamics, fast compound structures, which are composed of intermediate shocks (so-called 1-3 intermediate shocks) preceeded by fast rarefaction waves, are also possible.

Ryu, Dongsu↗
NASA NTRS2012

Dipole Alignment in Rotating MHD Turbulence

We present numerical results from long-term CPU and GPU simulations of rotating, homogeneous, magnetohydrodynamic (MHD) turbulence, and discuss their connection to the spherically bounded case. We compare our numerical results with a statistical theory of geodynamo action that has evolved from the absolute equilibrium ensemble theory of ideal MHD turbulence, which is based on the ideal MHD invariants are energy, cross helicity and magnetic helicity. However, for rotating MHD turbulence, the cross helicity is no longer an exact invariant, although rms cross helicity becomes quasistationary during an ideal MHD simulation. This and the anisotropy imposed by rotation suggests an ansatz in which an effective, nonzero value of cross helicity is assigned to axisymmetric modes and zero cross helicity to non-axisymmetric modes. This hybrid statistics predicts a large-scale quasistationary magnetic field due to broken ergodicity , as well as dipole vector alignment with the rotation axis, both of which are observed numerically. We find that only a relatively small value of effective cross helicity leads to the prediction of a dipole moment vector that is closely aligned (less than 10 degrees) with the rotation axis. We also discuss the effect of initial conditions, dissipation and grid size on the numerical simulations and statistical theory.

Shebalin, John V.↗
NASA NTRS1976

Magnetic dynamo action in two-dimensional turbulent magneto-hydrodynamics

Two-dimensional magnetohydrodynamic turbulence is explored by means of numerical simulation. Previous analytical theory, based on non-dissipative constants of the motion in a truncated Fourier representation, is verified by following the evolution of highly non-equilibrium initial conditions numerically. Dynamo action (conversion of a significant fraction of turbulent kinetic energy into long-wavelength magnetic field energy) is observed. It is conjectured that in the presence of dissipation and external forcing, a dual cascade will be observed for zero-helicity situations. Energy will cascade to higher wave numbers simultaneously with a cascade of mean square vector potential to lower wave numbers, leading to an omni-directional magnetic energy spectrum which varies as 1/k 3 at lower wave numbers, simultaneously with a buildup of magnetic excitation at the lowest wave number of the system. Equipartition of kinetic and magnetic energies is expected at the highest wave numbers in the system.

Fyfe, D.↗
NASA NTRS1976

Depletion of solar wind plasma near a planetary boundary

A mathematical model is presented that describes the squeezing of solar wind plasma out along interplanetary magnetic field lines in the region between the bow shock and the effective planetary boundary (in the case of the earth, the magnetopause). In the absence of local magnetic merging the squeezing process should create a 'depletion layer', a region of very low plasma density just outside the magnetopause. Numerical solutions are obtained for the dimensionless magnetohydrodynamic equations describing this depletion process for the case where the solar wind magnetic field is perpendicular to the solar wind flow direction. For the case of the earth, the theory predicts that the density should be reduced by a factor exceeding 2 in a layer about 700-1300 km thick if the Alfven Mach number in the solar wind, is equal to 8. Scaling of the model calculations to Venus and Mars suggests layer thicknesses about 1/10 and 1/15 those of the earth, respectively, neglecting diffusion and ionospheric effects.

Zwan, B. J.↗
NASA NTRS1977

Magnetic dynamo action in two-dimensional turbulent magneto-hydrodynamics

Two-dimensional magnetohydrodynamic turbulence is explored by means of numerical simulation. Previous analytical theory, based on non-dissipative constants of the motion in a truncated Fourier representation, is verified by following the evolution of highly non-equilibrium initial conditions numerically. Dynamo action (conversion of a significant fraction of turbulent kinetic energy into long-wavelength magnetic field energy) is observed. It is conjectured that in the presence of dissipation and external forcing, a dual cascade will be observed for zero-helicity situations. Energy will cascade to higher wavenumbers simultaneously with a cascade of mean square vector potential to lower wavenumbers, leading to an omni-directional magnetic energy spectrum.

Fyfe, D.↗
NASA NTRS1992

Dynamics and energetics of the solar corona

The primary objective of this research program is to improve our understanding of the dynamics and energetics of the solar corona both in the quiescent dynamic equilibrium state when coronal structure is dominated by the equatorial streamer belt and in the eruptive state when coronal plasma is ejected into the interplanetary medium. Numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations and comparisons of the computed results with observations form the core of the approach to achieving this objective. Some of the specific topics that have been studied are: (1) quiescent coronal streamers in an atmosphere dominated by a dipole magnetic field at large radii, (2) the formation of coronal mass ejections (CMEs) in quiescent streamers due to the emergence of new magnetic flux and due to photospheric shear motion, (3) MHD shock formation near the leading edge of CMEs, (4) coronal magnetic arcade eruption as a result of applied photospheric shear motion, and (5) the three-dimensional structure of CMEs.

Steinolfson, R. S.↗
NASA NTRS2007

Shell Models of RMHD Turbulence and the Heating of Solar Coronal Loops

A simplified nonlinear numerical model for the development of incompressible magnetohydrodynamics in the presence of a strong magnetic field B|| and stratification, nicknamed 'Shell-Atm,' is presented. In planes orthogonal to the mean field, the nonlinear incompressible dynamics is replaced by two-dimensional shell models for the complex variables u and b, allowing one to reach large Reynolds numbers while at the same time carrying out sufficiently long integrations to obtain good statistics at moderate computational cost. The shell models of different planes are coupled by Alfve'n waves propagating along B||. The model may be applied to open or closed magnetic field configurations where the axial field dominates and the plasma pressure is low; here we apply it to the specific case of a magnetic loop of the solar corona heated by means of turbulence driven by photospheric motions, and we use statistics for its analysis. The Alfven waves interact nonlinearly and form turbulent spectra in the directions perpendicular and, through propagation, also parallel to the mean field. A heating function is obtained and shown to be intermittent; the average heating is consistent with values required for sustaining a hot corona and is proportional to the aspect ratio of the loop to the -1.5 power, and haracteristic properties of heating events are distributed as power laws. Crosscorrelations show a delay of dissipation compared with energy content.

corona↗
NASA NTRSDate not supplied

Performance Analysis of Magnetohydrodynamic Drag Modulation for Actively Controlled Aerocapture at Neptune

While the Ice Giants are a top priority for flagship missions in the NASA’s near-term future, current technologies limit the scientific payload and mission timeline for future Uranian and Neptunian missions due to the need for fuel for orbit insertion maneuvers. Thus, to maximize the scientific potential of future missions, atmospheric aerocapture has been heavily researched. While atmospheric aerocapture simulations have proven enabling for capturing around Neptune, its deep atmospheric pass requires an aeroshell with robust thermal protection systems (TPS). Magnetohydrodynamically-controlled aerocapture serves as a potential solution to the limitations of both fully propulsive orbit insertion and aerodynamically controlled aerocapture. Through NASA Langley’s high-fidelity flight dynamics simulation, the Program to Optimize Simulated Trajectories II, both the aerodynamic-only and magnetohydrodynamic aerocapture methods were simulated and compared for identical missions to Neptune. After applying an optimized numerical predictor-corrector guidance algorithm for both methods, the results showed that magnetohydrodynamics has not only the control authority to successfully capture around Neptune, but also the unique advantage of a shallower atmospheric pass, decreasing the maximum heat load and the required TPS mass.

Danny N. Nguyen↗
NASA NTRS1974

Aligned magnetohydrodynamic solution for solar wind flow past the earth's magnetosphere

Exact numerical solutions of the MHD equations for a perfect dissipationless gas with aligned magnetic field are given for conditions representative of steady supersonic solar wind flow past an axisymmetric model of the earth's magnetosphere. The solution is based on use of a transformation that relates without approximation the equations of magnetohydrodynamics to those of gasdynamics of a pseudogas that has an unusual equation of state. The results confirm the applicability of the previously existing gasdynamic solutions for the typically modest intensities of the interplanetary magnetic field that lead to Alfven Mach numbers of about 10 or greater. For smaller values, however, significant effects are indicated with the flanks of the bow wave moving away from the earth and the nose moving toward the earth. The results are consistent with direct observations in space.

Spreiter, J. R.↗
NASA NTRS1979

Numerical experiments on the theta pinch

Numerical calculation of theta pinch problems are presented. Physical processes in theta pinch systems are considered in a one dimensional, two temperature magnetohydrodynamic, approximation with allowance for end losses by longitudinal heat conductivity. The numerical calculations are compared with results of earlier experiments.

Volosevich, P. P.↗
NASA NTRS1985

Streaming tearing mode

The magnetohydrodynamic stability of streaming tearing mode is investigated numerically. A bulk plasma flow parallel to the antiparallel magnetic field lines and localized in the neutral sheet excites a streaming tearing mode more strongly than the usual tearing mode, particularly for the wavelength of the order of the neutral sheet width (or smaller), which is stable for the usual tearing mode. Interestingly, examination of the eigenfunctions of the velocity perturbation and the magnetic field perturbation indicates that the streaming tearing mode carries more energy in terms of the kinetic energy rather than the magnetic energy. This suggests that the streaming tearing mode instability can be a more feasible mechanism of plasma acceleration than the usual tearing mode instability.

Shigeta, M.↗