Search NASA⌕ Search

SEARCH · Search NASA

Results for “MHD modeling”

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

A numerical code for a three-dimensional magnetospheric MHD equilibrium model

Two dimensional and three dimensional MHD equilibrium models were begun for Earth's magnetosphere. The original proposal was motivated by realizing that global, purely data based models of Earth's magnetosphere are inadequate for studying the underlying plasma physical principles according to which the magnetosphere evolves on the quasi-static convection time scale. Complex numerical grid generation schemes were established for a 3-D Poisson solver, and a robust Grad-Shafranov solver was coded for high beta MHD equilibria. Thus, the effects were calculated of both the magnetopause geometry and boundary conditions on the magnetotail current distribution.

Voigt, G.-H.↗

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

Magnetospheric equilibrium configurations and slow adiabatic convection

This review paper demonstrates how the magnetohydrostatic equilibrium (MHE) theory can be used to describe the large-scale magnetic field configuration of the magnetosphere and its time evolution under the influence of magnetospheric convection. The equilibrium problem is reviewed, and levels of B-field modelling are examined for vacuum models, quasi-static equilibrium models, and MHD models. Results from two-dimensional MHE theory as they apply to the Grad-Shafranov equation, linear equilibria, the asymptotic theory, magnetospheric convection and the substorm mechanism, and plasma anisotropies are addressed. Results from three-dimensional MHE theory are considered as they apply to an intermediate analytical magnetospheric model, magnetotail configurations, and magnetopause boundary conditions and the influence of the IMF.

Voigt, Gerd-Hannes↗

Models of Plumes: Their Flow, Their Geometric Spreading, and Their Mixing with Interplume Flow

There are two types of plume flow models: (1) 1D models using ad hoc spreading functions, f(r); (2) MagnetoHydroDynamics (MHD) models. 1D models can be multifluid, time dependent, and incorporate very general descriptions of the energetics. They confirm empirical results that plume flow is slow relative to requirements for high speed wind. But, no published 1 D model incorporates the rapid local spreading at the base (fl(r)) which has an important effect on mass flux. The one published MHD model is isothermal, but confirms that if b=8*pi*p/absolute value(B)2<<l then the field is nearly potential below -70,000 km. Building on the MHD result, we apply a two scale approximation to calculate fl(r). We also compute the global spreading (fg(r)) out to 5.0 RSUN imposed by coronal hole geometry. Global MHD models provide a potent method of calculating fg(r). Unambiguous plume signatures have not yet been found in the solar wind. This is probably due to strong mixing of plume and interplume flows near the Sun. We describe a physical source for strong mixing due to the observed flows being unstable to shear instabilities that lead to rapid disruption.

Suess, Steven T.↗

The effects of dipole tilt on the structure of the magnetosphere

The structure of the magnetosphere as a function of dipole tilt was studied by using a global magnetohydrodynamics (MHD) model. We show here that the cusp current sheet in the MHD model is related to the tail current sheet and that the calculated shape and positions of the tail current sheet are consistent with observations. In addition, we discuss the difference between the MHD model and the Chapman-Ferraro type model which manifests as a function of tilt angle.

Wu, C. C.↗

Shape of the magnetosphere

The shape of the magnetosphere has been studied by using a global MHD model of the magnetosphere. It is found that MHD and the Chapman-Ferraro models agree in total magnetopause current distributions, but that the MHD model presents a very different physical picture from that of the Chapman-Ferraro model. In the MHD model with B(IMF) = 0 or northward, the cusp is not exposed directly to the solar wind because the dayside magnetic field lines drape around the nightside magnetosphere. Furthermore, in the MHD model a current sheet exists above the cusp region.

Wu, C. C.↗

Observations of plasma dynamics in the coma of P/Halley by the Giotto Ion Mass Spectrometer

Observations in the coma of P/Halley by the Giotto Ion Mass Spectrometer (IMS) are reported. The High Energy Range Spectrometer (HERS) of the IMS obtained measurements of protons and alpha particles from the far upstream region to the near ionopause region and of ions from mass 12 to 32 at distances of about 250,000 to 40,000 km from the nucleus. Plasma parameters from the High Intensity Spectrometer (HIS) of the IMS obtained between 150,000 to 5000 km from the nucleus are also discussed. The distribution functions of water group ions (water group will be used to refer to ions of 16 to 18 m/q, where m is in AMU and q is in unit charges) are observed to be spherically symmetric in velocity space, indicating strong pitch angle scattering. The discontinuity known as the magnetic pile-up boundary (MPB) is apparent only in proton, alpha, and magnetometer data, indicating that it is a tangential discontinuity of solar wind origin. HERS observations show no significant change in the properties of the heavy ions across the MPB. A comparison of the observations to an MHD model is made. The plasma flow directions at all distances greater than 30,000 km from the nucleus are in agreement with MHD calculations. However, despite the agreement in flow direction, within 200,000 km of the nucleus the magnitude of the velocity is lower than predicted by the MHD model and the density is much larger (a factor of 4). Within 30,000 km of the nucleus there are large theoretical differences between the MHD model flow calculations for the plane containing the magnetic field and for the plane perpendicular to the magnetic field. The observations agreed much better with the pattern calculated for the plane perpendicular to the magnetic field. The data obtained by the High Energy Range Spectrometer (HERS) of the IMS that are published herein were provided to the International Halley Watch archive.

Goldstein, B. E.↗

Development of Solar Wind Model Driven by Empirical Heat Flux and Pressure Terms

We are developing a time stationary self-consistent 2D MHD model of the solar corona and solar wind as suggested by Sittler et al. (2003). Sittler & Guhathakurta (1999) developed a semiempirical steady state model (SG model) of the solar wind in a multipole 3-streamer structure, with the model constrained by Skylab observations. Guhathakurta et al. (2006) presented a more recent version of their initial work. Sittler et al. (2003) modified the SG model by investigating time dependent MHD, ad hoc heating term with heat conduction and empirical heating solutions. Next step of development of 2D MHD models was performed by Sittler & Ofman (2006). They derived effective temperature and effective heat flux from the data-driven SG model and fit smooth analytical functions to be used in MHD calculations. Improvements of the Sittler & Ofman (2006) results now show a convergence of the 3-streamer topology into a single equatorial streamer at altitudes > 2 R(sub S). This is a new result and shows we are now able to reproduce observations of an equatorially confined streamer belt. In order to allow our solutions to be applied to more general applications, we extend that model by using magnetogram data and PFSS model as a boundary condition. Initial results were presented by Selwa et al. (2008). We choose solar minimum magnetogram data since during solar maximum the boundary conditions are more complex and the coronal magnetic field may not be described correctly by PFSS model. As the first step we studied the simplest 2D MHD case with variable heat conduction, and with empirical heat input combined with empirical momentum addition for the fast solar wind. We use realistic magnetic field data based on NSO/GONG data, and plan to extend the study to 3D. This study represents the first attempt of fully self-consistent realistic model based on real data and including semi-empirical heat flux and semi-empirical effective pressure terms.

Sittler, Edward C., Jr.↗

Predicting the Structure of the Solar Corona During the December 4, 2002 Total Solar Eclipse

The solar magnetic field plays a key role in determining coronal. The principal input to MHD models is the observed solar magnetic field. 3D MHD models can be used to compare with eclipse and coronograph images, SOHO images (LOSCO, EIT), Ulysses and WIND spacecraft data, and interplanetary scintillation (IPS) measurements. MHD computations can tell us about the structure of the corona. Eclipses can help us to verify the accuracy of the models. 4 December, 2002 total eclipce: visible in the southern hemisphere (South Atlantic, southern Africa, Indian Ocean, and Australia). Total in center Angola is at 06:00 UT.

Mikic, Zoran↗

A Leaky Waveguide Model for MHD Wave Driven Winds from Coronal Holes

Magnetohydrodynamic (MHD) waves, driven by the large scale convective motions of the photosphere are suggested as a possible source of additional acceleration for the stellar wind. Most of the turbulent power in a coronal hole is carried by MHD waves with periods of a few hundred seconds or longer. This is evident from direct observations of turbulence in the solar photosphere, as well as in situ observations of turbulence in the solar wind. But waves with periods this long have wavelengths which are typically as large as the transverse scale of the coronal hole flux tube itself. For these waves boundary effects are important and the coronal hole must be treated as a waveguide. The propagation of MHD waves using this waveguide approach is discussed. The simple model presented demonstrates that coronal holes can act as waveguides for MHD waves. For typical solar parameters the waves are compressible and can generate a wave tensile force which tends to cancel at least part of the wave pressure force. This effect tends to decrease the efficiency of MHD wave acceleration.

Davila, J. M.↗

Space Weather Modeling and Prediction for Intermediate Time-scales

Much of solar activity within a sunspot cycle occurs as bursts, or 'seasons' of strong activity over several months, separated by periods of much less activity. The most important space weather effects occur during these bursts. Previous modeling and forecasting efforts have focused on time-scales of hours-to-days and decades-to-centuries. The recent discovery of Rossby waves in the Sun, together with recently developed global models of solar MHD Rossby waves and their interactions with differential rotation and spot-producing magnetic fields, reveal the opportunity to simulate and predict the occurrence, strength and location of enhanced activity bursts a few weeks up to several months in advance. We now have a golden opportunity to fill in this gap in time-scales of forecasting space weather. This requires a) continuous observations of solar Rossby waves by various techniques; b) development of coupled nonlinear MHD models that simulate both global Rossby waves and the much smaller spatial scale emergence of new active regions; c) application of advanced data assimilation techniques to couple surface observations to update the model-system to integrate forward in time for creating forecasts months ahead. Then it will be possible to build operational prediction models to meet the needs of customers and stakeholders, including support of future NASA missions, regarding what kind and level of space weather to expect a few weeks to several months ahead.

SMD↗

Magnetotail Current Sheet Thinning and Magnetic Reconnection Dynamics in Global Modeling of Substorms

Magnetotail current sheet thinning and magnetic reconnection are key elements of magnetospheric substorms. We utilized the global MHD model BATS-R-US with Adaptive Mesh Refinement developed at the University of Michigan to investigate the formation and dynamic evolution of the magnetotail thin current sheet. The BATSRUS adaptive grid structure allows resolving magnetotail regions with increased current density up to ion kinetic scales. We investigated dynamics of magnetotail current sheet thinning in response to southwards IMF turning. Gradual slow current sheet thinning during the early growth phase become exponentially fast during the last few minutes prior to nightside reconnection onset. The later stage of current sheet thinning is accompanied by earthward flows and rapid suppression of normal magnetic field component $B-z$. Current sheet thinning set the stage for near-earth magnetic reconnection. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. One of the major challenges in global MHD modeling of the magnetotail magnetic reconnection is to reproduce fast reconnection rates typically observed in smallscale kinetic simulations. Bursts of fast reconnection cause fast magnetic field reconfiguration typical for magnetospheric substorms. To incorporate nongyritropic effects in diffusion regions we developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated self-consistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites and spatial scales of the diffusion region are updated during the simulations. Such an approach allows quantifying the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions localized near reconnection sites. To clarify the role of smallscale non-MHD effects in diffusion region on the global magnetospheric dynamic and to test different models of dissipation we perform simulations with steady southward IMF driving of the magnetosphere.

Kuznetsova, M. M.↗

Nonlinear Fast Magnetosonic Waves in Solar Prominence Pillars

We investigate the properties of nonlinear fast magnetosonic (NFM) waves in a solar prominence, motivated by recent high-resolution and high-cadence Hinode/SOT observations of small-scale oscillations in a prominence pillar. As an example, we analyze the details of the 2012 February 14 Hinode/SOT observations of quasi-periodic propagating features consistent with NFM waves, imaged in emission in Ca II and in the far blue wing of Hα. We perform wavelet analysis and find oscillations in the 1-3 min period range. Guided by these observations, we model the NFM waves with a three-dimensional magnetohydrodynamics (3D MHD) model, extending previous 2.5D MHD studies. The new model includes the structure of the high-density, low-temperature material of the prominence pillar embedded in the hot corona, in both potential and non-force-free sheared magnetic field configurations. The nonlinear model demonstrates the effects of mode coupling and the propagating density compressions associated with linear and NFM waves. The guided fast magnetosonic waves, together with density compressions and currents, are reproduced in the 3D pillar structure. We demonstrate or the first time the dynamic effects of the Lorentz force due to the magnetic shear in the non-force-free field on the pillar structure and on the propagation of the waves. The insights gained from the 3D MHD modeling are useful for improving coronal seismology of prominence structures that exhibit fast MHD wave activity.

Leon Ofman↗

A simulation of the IPS variations from a magnetohydrodynamical simulation

Calculations of the variations of interplanetary scintillation (IPS) from a disturbance simulated by a 3-D magnetohydrodynamical (MHD) model of the solar wind are presented. The simulated maps are compared with observations and it is found that the MHD model reproduces the qualitative features of observed disturbances. The disturbance produced by the MHD simulation is found to correspond in strength with the weakest disturbance which can be reliably detected by existing single station IPS observations.

Tappin, S. J.↗

Model development for Ulysses and SOHO

The purpose of this research is to provide scientific expertise in solar physics and in the development and use of magnetohydrodynamic (MHD) models of coronal structures for the computation of Lyman alpha scattered radiation in these structures. The specific objectives will be to run MHD models with new boundary conditions and compute resulting scattered solar Lyman alpha intensities, guided by results from the first series of boundary conditions.

Wu, S. T.↗

The solar/interplanetary/magnetosphere/ionosphere connection - A strategy for prediction of geomagnetic storms

A physically-based strategy for the prediction of geomagnetic/ionospheric disturbances, the Solar-Terrestrial/Environmental Model (STEM 2000), is proposed, with application to the prediction of periods of spacecraft charging. Synoptic solar observations provide input for MHD models for flare occurrence, propagation of coronal disturbances, and high speed solar wind streams. A three-dimensional interplanetary MHD model determines solar wind parameters including the energy flux and the cross-magnetosphere tail electric field. Observational earth data are used to predict local time, high latitude ionospheric disturbances which have an impact on the ionospheric structure.

Dryer, M.↗

Kinematic Treatment of CME Evolution in the Solar Wind

We present a kinematic study of the evolution of coronal mass ejections (CMEs) in the solar wind. Specifically, we consider the effects of: (1) spherical expansion; and (2) uniform expansion due to pressure gradients between the Interplanetary CME (ICME) and the ambient solar wind. We compare these results with an MHD model, which allows us to isolate these effects from the combined kinematic and dynamical effects, which are included in MHD models. They also provide compelling evidence that the fundamental cross section of so-called "force-free" flux ropes (or magnetic clouds) is neither circular or elliptical, but rather a convex-outward, "pancake" shape. We apply a force-free fitting to the magnetic vectors from the MHD simulation to assess how the distortion of the flux rope affects the fitting. In spite of these limitations, force-free fittings, which are straightforward to apply, do provide an important description of a number of parameters, including the radial dimension, orientation and chirality of the ICME.

Riley, Pete↗

KINEMATIC TREATMENT OF CORONAL MASS EJECTION EVOLUTION IN THE SOLAR WIND

We present a kinematic study of the evolution of coronal mass ejections (CMEs) in the solar wind. Specifically, we consider the effects of (1) spherical expansion and (2) uniform expansion due to pressure gradients between the interplanetary CME (ICME) and the ambient solar wind. We compare these results with an MHD model that allows us to isolate these effects h m the combined kinematic and dynamical effects, which are included in MHD models. They also provide compelling evidence that the fundamental cross section of so-called "force-free" flux ropes (or magnetic clouds) is neither circular or elliptical, but rather a convex-outward, "pancake" shape. We apply a force-free fit to the magnetic vectors from the MHD simulation to assess how the distortion of the flux rope affects the fit. In spite of these limitations, force-free fits, which are straightforward to apply, do provide an important description of a number of parameters, including the radial dimension, orientation, and chirality of the ICME. Subject headings: MHD - solar wind - Sun: activity - Sun: corona - Sun: coronal mass ejections (CMEs) - On-line material color figures Sun: magnetic fields

Riley, Pete↗