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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Global MHD Models of the Solar Corona

Global magnetohydrodynamic (MHD) models of the solar corona are computationally intensive, numerically complex simulations that have produced important new results over the past few years. After a brief overview of how these models usually work, I will address three topics: (1) How these models are now routinely used to predict the morphology of the corona and analyze Earth and space-based remote observations of the Sun; (2) The direct application of these models to the analysis of physical processes in the corona and chromosphere and to the interpretation of in situ solar wind observations; and (3) The use of results from global models to validate the approximations used to make detailed studies of physical processes in the corona that are not otherwise possible using the global models themselves.

Suess, S. T.↗
NASA NTRS2012

Polarization and Compressibility of Oblique Kinetic Alfven Waves

Even though solar wind, as a collisionless plasma, is properly described by the kineticMaxwell-Vlasov description, it can be argued that much of our understanding of solar wind observational data comes from an interpretation and numerical modeling which is based on a fluid description of magnetohydrodynamics. In recent years, there has been a significant interest in better understanding the importance of kinetic effects, i.e. the differences between the kinetic and usual fluid descriptions. Here we concentrate on physical properties of oblique kinetic Alfvn waves (KAWs), which are often recognized as one of the key ingredients in the solar wind turbulence cascade. We use three different fluid models with various degrees of complexity and calculate polarization and magnetic compressibility of oblique KAWs (propagation angle q = 88), which we compare to solutions derived from linear kinetic theory. We explore a wide range of possible proton plasma b = [0.1,10.0] and a wide range of length scales krL = [0.001,10.0]. It is shown that the classical isotropic two-fluid model is very compressible in comparison with kinetic theory and that the largest discrepancy occurs at scales larger than the proton gyroscale. We also show that the two-fluid model contains a large error in the polarization of electric field, even at scales krL 1. Furthermore, to understand these discrepancies between the two-fluid model and the kinetic theory, we employ two versions of the Landau fluid model that incorporate linear low-frequency kinetic effects such as Landau damping and finite Larmor radius (FLR) corrections into the fluid description. It is shown that Landau damping significantly reduces the magnetic compressibility and that FLR corrections (i.e. nongyrotropic contributions) are required to correctly capture the polarization.We also show that, in addition to Landau damping, FLR corrections are necessary to accurately describe the damping rate of KAWs. We conclude that kinetic effects are important even at scales which are significantly larger than the proton gyroscale krL 1.

Turbulence↗
NASA NTRS2008

Gas Core Reactor Numerical Simulation Using a Coupled MHD-MCNP Model

Analysis is provided in this report of using two head-on magnetohydrodynamic (MHD) shocks to achieve supercritical nuclear fission in an axially elongated cylinder filled with UF4 gas as an energy source for deep space missions. The motivation for each aspect of the design is explained and supported by theory and numerical simulations. A subsequent report will provide detail on relevant experimental work to validate the concept. Here the focus is on the theory of and simulations for the proposed gas core reactor conceptual design from the onset of shock generations to the supercritical state achieved when the shocks collide. The MHD model is coupled to a standard nuclear code (MCNP) to observe the neutron flux and fission power attributed to the supercritical state brought about by the shock collisions. Throughout the modeling, realistic parameters are used for the initial ambient gaseous state and currents to ensure a resulting supercritical state upon shock collisions.

Kazeminezhad, F.↗
NASA NTRS2004

Nonlinear filtering and limiting in high order methods for ideal and non-ideal MHD

The various filtering mechanisms and base scheme options of the newly developed adaptive numerical dissipation control in spatially high order filter schemes for the ideal and non-ideal magnetohydrodynamics (MHD) equations are investigated. These filter schemes are applicable to complex unsteady MHD high-speed shock/shear/turbulence problems. They also provide a natural and efficient way for the minimization of Div(B) numerical error. The type of spatial base scheme to be used in conjunction with our filter idea is very general. For example, spectral, compact and non-compact spatially central finite difference schemes are possible candidates. The adaptive numerical dissipation mechanism consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and to leave the rest of the region free from numerical dissipation contamination. The numerical dissipation considered consists of high order linear dissipation for the suppression of high frequency oscillation and the nonlinear dissipative portion of high-resolution shock-capturing methods for discontinuity capturing. The applicable nonlinear dissipative portion of high-resolution shock-capturing methods is also very general. The objective of this paper is to investigate the performance of using compact and non-compact central base schemes in conjunction with three commonly used types of nonlinear numerical dissipation for both the ideal and non-ideal MHD. This extended abstract shows the performance of three nonlinear filters in conjunction with a sixth-order non-compact spatial central base scheme. In the final paper, the high order compact spatial central base scheme will be illustrated and compared with the non-compact base scheme. The reason for the investigation of the high order compact spatial central base scheme over the non-compact base scheme is to evaluate if additional accuracy can be gained in regions of fine scale turbulence that are away from shocks/shears.

Yee,H. C.↗
NASA NTRS2001

Computational Relativistic Astrophysics Using the Flowfield-Dependent Variation Theory

Theoretical models, observations and measurements have preoccupied astrophysicists for many centuries. Only in recent years, has the theory of relativity as applied to astrophysical flows met the challenges of how the governing equations can be solved numerically with accuracy and efficiency. Even without the effects of relativity, the physics of magnetohydrodynamic flow instability, turbulence, radiation, and enhanced transport in accretion disks has not been completely resolved. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks and also in the study of Gamma-Ray bursts (GRB). Thus, our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flowfield-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for Computational Relativistic Astrophysics (CRA) are demonstrated.

Richardson, G. A.↗
NASA NTRS1989

Boundary value problems in magnetohydrodynamics (and fluid dynamics). I - Radiation boundary condition

In view of the difficulties that emerge in numerical simulations of time-dependent phenomena involving MHD waves due to waves' artificial reflection at numerical boundaries, where time-dependent magnetic fields outside those boundaries are unknown, a method previously developed for the determination of radiative boundary conditions is generalized for the treatment of isothermal MHD waves. The illustrative problem thus treated is that of the magnetic braking of a rigidly rotating disk which is embedded in a differentially-rotating fluid. The analytic results obtained are compared with those from numerical simulations with and without the present boundary treatment.

Vanajakshi, T. C.↗
NASA NTRSDate not supplied

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

Missions to the Ice Giants are a top priority for flagship missions this coming decade. However, a fully propulsive orbital insertion into these planets requires an immense amount of fuel, taking a significant portion of the spacecraft mass and restricting the scientific payload. To mitigate this, aerocapture has been heavily investigated. Although simulations have shown that aerodynamically controlled aerocapture can successfully insert into an orbit around both Ice Giants, the deep atmospheric pass required necessitates a complex, mass expensive, and sometimes prohibitive thermal protection system. Magnetohydrodynamic drag modulation serves as a potential alternative control method for aerocapture which could not only save propellant mass compared to fully propulsive orbital insertion, but also save thermal protection system mass compared to conventional aerocapture methods. Both aerodynamically controlled and magnetohydrodynamically controlled aerocapture methodologies were simulated in NASA Langley’s high-fidelity six degree-of-freedom flight dynamics code, the Program to Optimize Simulated Trajectories II. Each method was simulated to identical missions to Neptune using a numerical predictor-corrector algorithm to optimize the control towards the target orbit. The results concluded that magnetohydrodynamic drag modulation can successfully capture and performs on par with aerodynamic drag modulation while significantly reducing the heatflux and aeroshell complexity.

Danny N Nguyen↗
NASA NTRS1996

Numerical Simulations of Mass Loading in the Solar Wind Interaction with Venus

Numerical simulations are performed in the framework of nonlinear two-dimensional magnetohydrodynamics to investigate the influence of mass loading on the solar wind interaction with Venus. The principal physical features of the interaction of the solar wind with the atmosphere of Venus are presented. The formation of the bow shock, the magnetic barrier, and the magnetotail are some typical features of the interaction. The deceleration of the solar wind due to the mass loading near Venus is an additional feature. The effect of the mass loading is to push the shock farther outward from the planet. The influence of different values of the magnetic field strength on plasma evolution is considered.

Murawski, K.↗
NASA NTRS1980

Numerical studies of new stellarator concepts

A three-dimensional computer code has been developed to study the magnetohydrodynamic equilibrium and stability of a diffuse or sharp boundary plasma in toroidal geometry. It is shown how equilibria with net toroidal current identically zero can be determined and how growth rates of instabilities can be calculated. Applications are made to an l = 2, 3 stellarator configuration that offers the possibility of achieving a critical value as high as 10% for the plasma parameter

Bauer, F.↗
NASA NTRS1991

Flux-corrected transport techniques for multidimensional compressible magnetohydrodynamics

A prescription is given for conservatively integrating generalized hydromagnetic equations using flux-corrected transport (FCT) techniques. By placing the magnetic-field components at the interface locations of the finite-difference grid, the field is kept divergence-free to within machine roundoff error. The use of FCT techniques allows an integration scheme of high accuracy to be employed, while the numerical ripples associated with large dispersion errors are avoided. The method is particularly well suited for problems involving magnetohydrodynamic shocks and other discontinuities.

Devore, C. Richard↗
NASA NTRS1972

The responses of the thermosphere due to a geomagnetic storm: A MHD model

A magnetohydrodynamics theory was used to study the dynamic response of the neutral atmosphere to a geomagnetic storm. A full set of magnetohydrodynamic equations appropriate for the present problem is derived and their various orders of approximation are discussed in some detail. In order to demonstrate the usefulness of this theoretical model, the May 1967 geomagnetic storm data were used in the resulting set of nonlinear, time dependent, partial differential magnetohydrodynamic equations to calculate variations of the thermosphere due to the storm. The numerical results are presented for wind speeds, electric field strength, and amount of joule heating at a constant altitude for the data recorded. Data show that the strongest thermospheric responses are at the polar region becoming weaker in the equatorial region. This may lead to the speculation that a thermospheric wave is generated in the polar region due to the geomagnetic storm which propagates towards the equator.

Wu, S. T.↗
NASA NTRS2004

Annual Research Briefs, 2004: Center for Turbulence Research

This report contains the 2004 annual progress reports of the Research Fellows and students of the Center for Turbulence Research in its eighteenth year of operation. Since its inception in 1987, the objective of the CTR has been to advance the physical understanding of turbulent flows and development of physics based predictive tools for engineering analysis and turbulence control. Turbulence is ubiquitous in nature and in engineering devices. The studies at CTR have been motivated by applications where turbulence effects are significant; these include a broad range of technical areas such as planetary boundary layers, formation of planets, solar convection, magnetohydrodynamics, environmental and eco systems, aerodynamic noise, propulsion systems and high speed transportation. Numerical simulation has been the predominant research tool at CTR which has required a critical mass of researchers in numerical analysis and computer science in addition to core disciplines such as applied mathematics, chemical kinetics and fluid mechanics. Maintaining and promoting this interdisciplinary culture has been a hallmark of CTR and has been responsible for the realization of the results of its basic research in applications. The first group of reports in this volume are directed towards development, analysis and application of novel numerical methods for ow simulations. Development of methods for large eddy simulation of complex flows has been a central theme in this group. The second group is concerned with turbulent combustion, scalar transport and multi-phase ows. The nal group is devoted to geophysical turbulence where the problem of solar convection has been a new focus of considerable attention recently at CTR.

Moin, Parviz↗
NASA NTRS2006

GRMHD Simulations of Jet Formation with RAISHIN

We have developed a new three dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high-resolution shock capturing scheme. Numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. We describe code performance on some test problems in both special and general relativity. Our new GRMHD code has proven to be accurate to second order and has successfully passed several numerical test problems including highly relativistic and magnetized tests in both special and general relativity. We have performed several simulations of non-rotating and rotating black hole systems with a geometrically thin accretion disk. The simulations show the formation of jets driven by the Lorentz force and the gas pressure. It appears that the rotating black hole creates an additional faster, and more collimated outflow inside a broader, slower outflow that is also generated by the rotating accretion disk around a non-rotating black hole. The kinematic jet structure could thus be a sensitive function of black hole rotation.

Mizuno, Y.↗
NASA NTRS1989

Magnetic dynamo activity in mechanically driven compressible magnetohydrodynamic turbulence

Magnetic dynamo activity in a homogeneous, dissipative, polytropic, two-dimensional, turbulent magneto-fluid is simulated numerically. The magneto-fluid is simulated numerically. The magneto-fluid is, in a number of cases, mechanically forced so that energy input balances dissipation, thereby maintaining constant energy. In the presence of a mean magnetic field, a magneto-fluid whose initial turbulent magnetic energy is zero quickly arrives at a state of non-zero turbulent magnetic energy. If the mean magnetic field energy density is small, the turbulent magnetic field can achieve a local energy density more than four hundred times larger; if the mean magnetic field energy density is large, then equipartition between the turbulent magnetic and kinetic energy is achieved. Compared to the presence of a mean magnetic field, compressibility appears to have only a marginal effect in mediating the transfer of turbulent kinetic energy into magnetic energy.

Shebalin, John V.↗
NASA NTRS2010

Interchange Reconnection and Coronal Hole Dynamics

We investigate the effect of magnetic reconnection between open and closed fields, often referred to as interchange reconnection, on the dynamics and topology of coronal hole boundaries. The most important and most prevalent three-dimensional topology of the interchange process is that of a small-scale bipolar magnetic field interacting with a large-scale background field. We determine the evolution of such a magnetic topology by numerical solution of the fully three-dimensional MHD equations in spherical coordinates. First, we calculate the evolution of a small-scale bipole that initially is completely inside an open field region and then is driven across a coronal hole boundary by photospheric motions. Next the reverse situation is calculated in which the bipole is initially inside the closed region and driven toward the coronal hole boundary. In both cases, we find that the stress imparted by the photospheric motions results in deformation of the separatrix surface between the closed field of the bipole and the background field, leading to rapid current sheet formation and to efficient reconnection. When the bipole is inside the open field region, the reconnection is of the interchange type in that it exchanges open and closed fields. We examine, in detail, the topology of the field as the bipole moves across the coronal hole boundary and find that the field remains well connected throughout this process. Our results, therefore, provide essential support for the quasi-steady models of the open field, because in these models the open and closed flux are assumed to remain topologically distinct as the photosphere evolves. Our results also support the uniqueness hypothesis for open field regions as postulated by Antiochos et al. On the other hand, the results argue against models in which open flux is assumed to diffusively penetrate deeply inside the closed field region under a helmet streamer. We discuss the implications of this work for coronal observations.

MHD (MAGNETOHYDRODYNAMICS)↗
NASA NTRS1985

Particle behavior in the magnetosphere

The Rice Convection Model deals with large-scale processes in the earth's inner and middle magnetosphere, including coupling to the ionosphere. Starting from appropriate initial and boundary conditions, the model computes the following physical parameters: ionospheric electric fields and currents; magnetospheric particle distributions, electric fields, and electric currents; and magnetic-field-aligned (Birkeland) currents connecting the two regions. This paper evaluates work on the model, with emphasis on the assumptions made, the basic equations, and the numerical methods. The theoretical basis of the model is compared and contrasted with standard magnetohydrodynamics. The limitations imposed by the major assumptions are discussed. Model inputs and boundary conditions are listed, and the methods of specifying them discussed. Some physical conclusions and insights that have been gained from the model are listed and described very briefly. References are given to published discussions of the major points of physics.

Wolf, R. A.↗
NASA NTRS1994

Coronal heating by the resonant absorption of Alfven waves: The effect of viscous stress tensor

The time-dependent linearized magnetohydrodynamics (MHD) equations for a fully compressible, low-beta, viscoresistive plasma are solved numerically using an implicit integration scheme. The full viscosity stress tensor (Braginskii 1965) is included with the five parameters eta(sub i) i = 0 to 4. In agreement with previous studies, the numerical simulations demonstrate that the dissipation on inhomogeneities in the background Alfven speed occurs in a narrow resonant layer. For an active region in the solar corona the values of eta(sub i) are eta(sub o) = 0.65 g/cm/s, eta(sub 1) = 3.7 x 10(exp -12) g/cm/s, eta(sub 2) = 4 eta(sub 1), eta(sub 3) = 1.4 x 10(exp -6) g/cm/s, eta(sub 4) = 2 eta(sub 3), with n = 10(exp 10)/cu cm, T = 2 x 10(exp 6) K, and B = 100 G. When the Lundquist number S = 10(exp 4) and R(sub 1) much greater than S (where R(sub 1) is the dimensionless shear viscous number) the width of the resistive dissipation layer d(sub r) is 0.22a (where a is the density gradient length scale) and d(sub r) approximately S(exp -1/3). When S much greater than R(sub 1) the shear viscous dissipation layer width d(sub r) scales as R(sub 1)(exp -1/3). The shear viscous and the resistive dissipation occurs in an overlapping narrow region, and the total heating rate is independent of the value of the dissipation parameters in agreement with previous studies. Consequently, the maximum values of the perpendicular velocity and perpendicular magnetic field scale as R(sub 1)(exp -1/3). It is evident from the simulations that for solar parameters the heating due to the compressive viscosity (R(sub 0) = 560) is negligible compared to the resistive and the shear viscous (R(sub 1)) dissipation and it occurs in a broad layer of order a in width. In the solar corona with S approximately equals 10(exp 4) and R(sub 1) approximately equals 10(exp 14) (as calculated from the Braginskii expressions), the shear viscous resonant heating is of comparable magnitude to the resistive resonant heating.

Ofman, L.↗