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 127 records · Page 7

A two-dimensional MHD global coronal model - Steady-state streamers

A 2D, time-dependent, numerical, MHD model for the simulation of coronal streamers from the solar surface to 15 solar is presented. Three examples are given; for dipole, quadrupole and hexapole (Legendre polynomials P1, P2, and P3) initial field topologies. The computed properties are density, temperature, velocity, and magnetic field. The calculation is set up as an initial-boundary value problem wherein a relaxation in time produces the steady state solution. In addition to the properties of the solutions, their accuracy is discussed. Besides solutions for dipole, quadrupole, and hexapole geometries, the model use of realistic values for the density and Alfven speed while still meeting the requirement that the flow speed be super-Alfvenic at the outer boundary by extending the outer boundary to 15 solar radii.

Wang, A.-H.↗

Characterization of the Optical and X-ray Properties of the Northwestern Wisps in the Crab Nebula

We have studied the wisps to the northwest of the Crab pulsar as part of a multi-wavelength campaign in the visible and in X-rays. Optical observations were obtained using the Nordic Optical Telescope in La Palma and X-ray observations were made with the Chandra X-ray Observatory. The observing campaign took place from October 2010 until September 2012. About once per year we observe wisps forming and peeling off from (or near) the region commonly associated with the termination shock of the pulsar wind. We find that the exact locations of the northwestern wisps in the optical and in X-rays are similar but not coincident, with X-ray wisps preferentially located closer to the pulsar. This suggests that the optical and X-ray wisps are not produced by the same particle distribution. It is also interesting to note that the optical and radio wisps are also separated from each other (Bietenholz et al. 2004). Our measurements and their implications are interpreted in terms of a Doppler-boosted ring model that has its origin in MHD modeling. While the Doppler boosting factors inferred from the X-ray wisps are consistent with current MHD simulations of PWNe, the optical boosting factors are not, and typically exceed values from MHD simulations by about a factor of 4.

Weisskopf, M. C.↗

The solar-flare induced earth's environment

A composite numerical simulation model developed from a series of MHD models was used to compute the solar-flare-generated disturbances of physical parameters, such as density, temperature, velocity, and magnetic field from the solar surface (i.e., the photospheric level) to the earth's environment. It is shown that the disturbed earth's environment at high latitudes can be approximated by starting with the knowledge of the occurrence and the strength of a solar flare, then simulating the evolutionary consequences of the solar disturbance through interplanetary space up to and through the magnetosphere.

Wu, S. T.↗

Do We Know the Actual Magnetopause Position for Typical Solar Wind Conditions?

We compare predicted magnetopause positions at the subsolar point and four reference points in the terminator plane obtained from several empirical and numerical MHD (magnetohydrodynamics) models. Empirical models using various sets of magnetopause crossings and making different assumptions about the magnetopause shape predict significantly different magnetopause positions (with a scatter greater than 1 Earth radius (R (sub E)) even at the subsolar point. Axisymmetric magnetopause models cannot reproduce the cusp indentations or the changes related to the dipole tilt effect, and most of them predict the magnetopause closer to the Earth than non axisymmetric models for typical solar wind conditions and zero tilt angle. Predictions of two global non axisymmetric models do not match each other, and the models need additional verification. MHD models often predict the magnetopause closer to the Earth than the non axisymmetric empirical models, but the predictions of MHD simulations may need corrections for the ring current effect and decreases of the solar wind pressure that occur in the foreshock. Comparing MHD models in which the ring current magnetic field is taken into account with the empirical Lin et al. model, we find that the differences in the reference point positions predicted by these models are relatively small for B (sub z) equals 0 (note: B (sub z) is when the Earth's magnetic field points north versus Sun's magnetic field pointing south). Therefore, we assume that these predictions indicate the actual magnetopause position, but future investigations are still needed.

Samsonov, A. A.↗

Validation of Model Forecasts of the Ambient Solar Wind

Independent and automated validation is a vital step in the progression of models from the research community into operational forecasting use. In this paper we describe a program in development at the CCMC to provide just such a comprehensive validation for models of the ambient solar wind in the inner heliosphere. We have built upon previous efforts published in the community, sharpened their definitions, and completed a baseline study. We also provide first results from this program of the comparative performance of the MHD models available at the CCMC against that of the Wang-Sheeley-Arge (WSA) model. An important goal of this effort is to provide a consistent validation to all available models. Clearly exposing the relative strengths and weaknesses of the different models will enable forecasters to craft more reliable ensemble forecasting strategies. Models of the ambient solar wind are developing rapidly as a result of improvements in data supply, numerical techniques, and computing resources. It is anticipated that in the next five to ten years, the MHD based models will supplant semi-empirical potential based models such as the WSA model, as the best available forecast models. We anticipate that this validation effort will track this evolution and so assist policy makers in gauging the value of past and future investment in modeling support.

Macneice, P. J.↗

Magnetohydrodynamic modeling of the solar corona

The ideal and resistive magnetohydrodynamic (MHD) model is used to examine the dynamics and structure of the solar corona. When the coronal magnetic field is deformed by photospheric flow it can evolve to states that become unstable to ideal MHD modes. The nonlinear evolution of these instabilities can lead to the generation of current sheets, field line reconnection, and energy release. The disruption of an arcade field and the kinking of coronal loops is described. The braiding of the large-scale coronal field by convective photospheric motions develops fine-scale structure in the magnetic field and leads to the development of intense current filaments. The resistive dissipation of these currents can provide an efficient coronal heating mechanism.

Mikic, Zoran↗

A test suite for magnetohydrodynamical simulations

A collection is presented of MHD problems that will provide researchers who have interests in modeling MHD phenomena with a battery of tests (a 'test suite') for calibrating their numerical algorithms. Use of these tests will provide a common reference for comparison of different numerical MHD algorithms. The test suite includes both 1D and 2D problems. Taken together, these problems test the abilities of a numerical method to propagate accurately all of the MHD wave families in moving and stationary media, to capture MHD shocks and contact discontinuities, and to model accurately Lorentz force terms in multiple dimensions. An example solution of each test problem is presented. Diagnostic convergence-testing procedures that provide a quantitative evaluation of MHD algorithms are demonstrated.

Stone, James M.↗

Solar Wind Charge Exchange During Geomagnetic Storms

On March 31st. 2001, a coronal mass ejection pushed the subsolar magnetopause to the vicinity of geosynchronous orbit at 6.6 RE. The NASA/GSFC Community Coordinated Modeling Center (CCMe) employed a global magnetohydrodynamic (MHD) model to simulate the solar wind-magnetosphere interaction during the peak of this geomagnetic storm. Robertson et aL then modeled the expected 50ft X-ray emission due to solar wind charge exchange with geocoronal neutrals in the dayside cusp and magnetosheath. The locations of the bow shock, magnetopause and cusps were clearly evident in their simulations. Another geomagnetic storm took place on July 14, 2000 (Bastille Day). We again modeled X-ray emission due to solar wind charge exchange, but this time as observed from a moving spacecraft. This paper discusses the impact of spacecraft location on observed X-ray emission and the degree to which the locations of the bow shock and magnetopause can be detected in images.

Robertson, Ina P.↗

Structure of the dayside reconnection layer in resistive MHD and hybrid models

Numerical simulations were performed to investigate the structure of the reconnection layer at the dayside magnetopause. Two typical cases are examined in detail; both are asymmetric in magnetic field and plasma density. In case 1, the guide fields in the magnetosheath and in the magnetosphere are set at zero and thus the tangential magnetic fields on the two sides of the initial current sheet are exactly antiparallel. In case 2, the angle between the tangential magnetic fields on the two sides of the initial current sheet is 145 deg. The results obtained from a resistive MHD model and from a hybrid model are found to be different. In the MHD simulation of case 1, a 2-4 intermediate shock is found to bound the reconnection layer on the magnetosheath side, while an Alfven wave pulse bounds the reconnection layer on the magnetospheric side. In case 2, it is found that a time-dependent intermediate shock (TDIS) bounds the reconnection layer on the magnetosheath side, with a slow expansion wave propagating behind. With the MHD simulations, in the general case in which the tangential magnetic fields on the two sides of the initial current sheet are not exactly antiparallel, a rotational discontinuity across which the tangential magnetic field rotates, a large angle is found to bound the reconnection layer on the magnetosheath side.

Lin, Y.↗

Semiempirical Two-Dimensional Magnetohydrodynamic Model of the Solar Corona and Interplanetary Medium

We have developed a two-dimensional semiempirical MHD model of the solar corona and solar wind. The model uses empirically derived electron density profiles from white-light coronagraph data measured during the Skylub period and an empirically derived model of the magnetic field which is fitted to observed streamer topologies, which also come from the white-light coronagraph data The electron density model comes from that developed by Guhathakurta and coworkers. The electron density model is extended into interplanetary space by using electron densities derived from the Ulysses plasma instrument. The model also requires an estimate of the solar wind velocity as a function of heliographic latitude and radial component of the magnetic field at 1 AU, both of which can be provided by the Ulysses spacecraft. The model makes estimates as a function of radial distance and latitude of various fluid parameters of the plasma such as flow velocity V, effective temperature T(sub eff), and effective heat flux q(sub eff), which are derived from the equations of conservation of mass, momentum, and energy, respectively. The term effective indicates that wave contributions could be present. The model naturally provides the spiral pattern of the magnetic field far from the Sun and an estimate of the large-scale surface magnetic field at the Sun, which we estimate to be approx. 12 - 15 G. The magnetic field model shows that the large-scale surface magnetic field is dominated by an octupole term. The model is a steady state calculation which makes the assumption of azimuthal symmetry and solves the various conservation equations in the rotating frame of the Sun. The conservation equations are integrated along the magnetic field direction in the rotating frame of the Sun, thus providing a nearly self-consistent calculation of the fluid parameters. The model makes a minimum number of assumptions about the physics of the solar corona and solar wind and should provide a very accurate empirical description of the solar corona and solar wind Once estimates of mass density rho, flow velocity V, effective temperature T(sub eff), effective heat flux q(sub eff), and magnetic field B are computed from the model and waves are assumed unimportant, all other plasma parameters such as Mach number, Alfven speed, gyrofrequency, etc. can be derived as a function of radial distance and latitude from the Sun. The model can be used as a planning tool for such missions as Slar Probe and provide an empirical framework for theoretical models of the solar corona and solar wind The model will be used to construct a semiempirical MHD description of the steady state solar corona and solar wind using the SOHO Large Angle Spectrometric Coronagraph (LASCO) polarized brightness white-light coronagraph data, SOHO Extreme Ultraviolet Imaging Telescope data, and Ulysses plasma data.

Sittler, Edward C., Jr.↗

Physics and evolution of constant opening angle jets using a quasi-one-dimensional magnetohydrodynamic model

We discuss the significance of the assumptions of infinite conductivity and time independence in the context of an ideal MHD model for constant opening angle jets. The model is developed by projecting the MHD equations onto the jet axis. We find that for initially sub-Alfvenic flows (i.e., flows emanating from active galactic nuclei and neutron stars) wind-type solutions exist only when the field lines at the origin are wound up in a direction opposite to the direction of rotation of the compact source. We discuss the possibility that the time evolution of these outflows may be a cycle between breeze- and wind-type solutions as a result of continuous changes in the boundary conditions at the origin due to accretion. We propose that such cycles may explain the apparent one-sideness of some jets, especially the ones for which we cannot use arguments of relativistic beaming. We examine the dependence of the wind-type solutions on the following parameters describing the outflow at the origin: the degree of winding of the field lines, the value of the gas pressure, the polytropic index, the strength of the magnetic field, the value of the rotational velocity, the gravitational potential of the compact object, and the injection velocity. We compare results with results obtained previously, and discuss briefly the qualitative features and physical interpretation of the solutions for outflows emanating from neutron stars and protostars.

Koupelis, Theodoros↗

Solar maximum mission: Ground support programs at the Harvard Radio Astronomy Station

Observations of the spectral characteristics of solar radio bursts were made with new dynamic spectrum analyzers of high sensitivity and high reliability, over the frequency range 25-580 MHz. The observations also covered the maximum period of the current solar cycle and the period of international cooperative programs designated as the Solar Maximum Year. Radio data on shock waves generated by solar flares were combined with optical data on coronal transients, taken with equipment on the SMM and other satellites, and then incorporated into computer models for the outward passage of fast-mode MHD shocks through the solar corona. The MHD models are non-linear, time-dependent and for the most recent models, quasi-three-dimensional. They examine the global response of the corona for different types of input pulses (thermal, magnetic, etc.) and for different magnetic topologies (for example, open and closed fields). Data on coronal shocks and high-velocity material ejected from solar flares have been interpreted in terms of a model consisting of three main velocity regimes.

Maxwell, A.↗

Modeling the Magnetic and Thermal Structure of Active Regions: 1st Year 1st Semi-Annual Progress Report

This report covers technical progress during the first six months of the first year of NASA SR&T contract "Modeling the Magnetic and Thermal Structure of Active Regions", NASW-03008, between NASA and Science Applications International Corporation, and covers the period January 14, 2003 to July 13, 2003. Under this contract SAIC has conducted research into theoretical modeling of the properties of active regions using the MHD model.

Mikic, Zoran↗

Dst Index in the 2008 GEM Modeling Challenge - Model Performance for Moderate and Strong Magnetic Storms

The GEM 2008 modeling challenge efforts are expanding beyond comparing in-situ measurements in the magnetosphere and ionosphere to include the computation of indices to be compared. The Dst index measures the largest deviations of the horizontal magnetic field at 4 equatorial magnetometers from the quiet-time background field and is commonly used to track the strength of the magnetic disturbance of the magnetosphere during storms. Models can calculate a proxy Dst index in various ways, including using the Dessler-Parker Sckopke relation and the energy of the ring current and Biot-Savart integration of electric currents in the magnetosphere. The GEM modeling challenge investigates 4 space weather events and we compare models available at CCMC against each other and the observed values of Ost. Models used include SWMF/BATSRUS, OpenGGCM, LFM, GUMICS (3D magnetosphere MHD models), Fok-RC, CRCM, RAM-SCB (kinetic drift models of the ring current), WINDMI (magnetosphere-ionosphere electric circuit model), and predictions based on an impulse response function (IRF) model and analytic coupling functions with inputs of solar wind data. In addition to the analysis of model-observation comparisons we look at the way Dst is computed in global magnetosphere models. The default value of Dst computed by the SWMF model is for Bz the Earth's center. In addition to this, we present results obtained at different locations on the Earth's surface. We choose equatorial locations at local noon, dusk (18:00 hours), midnight and dawn (6:00 hours). The different virtual observatory locations reveal the variation around the earth-centered Dst value resulting from the distribution of electric currents in the magnetosphere during different phases of a storm.

Rastaetter, Lutz↗

Quantifying the Ability of Magnetohydrodynamic Models to Reproduce Observed Birkeland Current and Auroral Electrojet Magnitudes

Although global magnetohydrodynamic (MHD) models have increased in sophistication and are now at the forefront of modeling Space Weather, there is still no clear understanding of how well these models replicate the observed ionospheric current systems. Without a full understanding and treatment of the ionospheric current systems, global models will have significant shortcomings that will limit their use. In this study we focus on reproducing observed seasonal interhemispheric asymmetry in ionospheric currents using the Space Weather Modeling Framework (SWMF). We find that SWMF does reproduce the linear relationship between the electrojets and the FACs, despite the underestimation of the currents’ magnitudes. Quantitatively, we find that at best SWMF is only capturing approximately 60% of the observed current. We also investigate how varying F10.7 effects the ionospheric potential and currents during the summer and winter. We find that simulations ran with higher F10.7 result in lower ionospheric potentials. Additionally, we find that the models do not always replicate the expected behavior of the currents with varying F10.7. This work points to a needed improvement in ionospheric conductance models.

Tre'Shunda James↗

Acceleration of the Fast Solar Wind by Solitary Waves in Coronal Holes

The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear, time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that solitary-like waves are generated in coronal holes nonlinearly by torsional Alfven waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature, and magnetic field geometry) that will become available from the recently launched SOHO spacecraft.

Ofman, Leon↗

Acceleration of the Fast Solar Wind by Solitary Waves in Coronal Holes

The purpose of this investigation is to develop a new model for the acceleration of the fast solar wind by nonlinear. time-dependent multidimensional MHD simulations of waves in solar coronal holes. Preliminary computational studies indicate that nonlinear waves are generated in coronal holes by torsional Alfv\'{e}n waves. These waves in addition to thermal conduction may contribute considerably to the accelerate the solar wind. Specific goals of this proposal are to investigate the generation of nonlinear solitary-like waves and their effect on solar wind acceleration by numerical 2.5D MHD simulation of coronal holes with a broad range of plasma and wave parameters; to study the effect of random disturbances at the base of a solar coronal hole on the fast solar wind acceleration with a more advanced 2.5D MHD model and to compare the results with the available observations; to extend the study to a full 3D MHD simulation of fast solar wind acceleration with a more realistic model of a coronal hole and solar boundary conditions. The ultimate goal of the three year study is to model the, fast solar wind in a coronal hole, based on realistic boundary conditions in a coronal hole near the Sun, and the coronal hole structure (i.e., density, temperature. and magnetic field geometry,) that will become available from the recently launched SOHO spacecraft.

Ofman, Leon↗

Modeling of energy buildup for a flare-productive region

A self-consistent MHD model of shearing magnetic loops is used to investigate magnetic energy buildup in active region AR 2372 (Boulder number), in the period of April 5-7, 1980. The magnetic field and sunspot motions in this region, derived using observational data obtained by the Marshall Space Flight Center Solar Observatory, suggest the initial boundary conditions for the model. It is found that the plasma parameters (i.e., density, temperature, and plasma flow velocity) do not change appreciably during the process of energy buildup as the magnetic loops are sheared. Thus, almost all of the added energy is stored in the magnetic field. Furthermore, it is shown that dynamical processes are not important during a slow buildup (i.e., for a shearing velocity less than 1 km/s). Finally, it is concluded that the amount of magnetic energy stored and the location of this stored magnetic energy depend on the initial magnetic field (whether potential or sheared) and the magnitude of the shearing motion.

Wu, S. T.↗