SEARCH · Search NASA
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.
This report covers technical progress during the first quarter of the second year of NASA Sun-Earth Connections Theory Program (SECTP). SAIC and the University of California, Irvine (UCI) have conducted research into theoretical modeling of active regions, the solar corona, and the inner heliosphere, using the MHD model.
We present the results of magnetohydrodynamic (MHD) modeling of winds from luminous late-type stars using a 2.5-dimensional, nonlinear MHD computer code. We assume that the wind is generated within an initially hydrostatic atmosphere and is driven by torsional Alfven waves generated at the stellar surface. Two cases of atmospheric topology are considered: case I has longitudinally uniform density distribution and isotropic radial magnetic field over the stellar surface, and case II has an isotropic, radial magnetic field with a transverse density gradient, which we refer to as an "atmospheric hole." We use the same set of boundary conditions for both models. The calculations are designed to model a cool luminous star, for which we assume an initial hydrostatic pressure scale height of 0.072 Stellar Radius, an Alfven wave speed of 92 km/s at the surface, and a wave period of 76 days, which roughly corresponds with the convective turnover time. For case I the calculations produce a wind with terminal velocity of about 22 km/s and a mass loss rate comparable to the expected value of 10(exp -6) Solar Mass/yr. For case II we predict a two-component wind: a fast (25 km/s) and relatively dense wind outside of the atmospheric hole and a slow (1.5 km/s), rarefied wind inside of the hole.
Using a modification of the one-dimensional multispecies 'one-major-ion' MHD model of Shinagawa et al. (1987), the behaviors of plasma and magnetic field in the dayside ionosphere of Venus was studied for both time-dependent and steady-state conditions. The present model is more complete than the one-major-ion model of Shinagawa et al., although a comparison of the results indicated that the one-major-ion treatment was a fairly good approximation. Two new cases are presented, including steady-state conditions for the magnetized ionosphere, and the inclusion of ion loss due to horizontal transport in the magnetized region. The resulting calculated profiles of the magnetic field and the electron density agree much better with the observations at high altitudes than those without the ion loss terms, indicating the importance of the horizontal transport processes in the ionosphere of Venus at high altitudes.
In recent Ultraviolet Coronagraph Spectrometer/Solar and Heliospheric Observatory (UVCS/SOHO) White Light Channel (WLC) observations we found quasi-periodic variations in the polarized brightness (pB) in the polar coronal holes at heliocentric distances of 1.9-2.45 solar radii. The motivation for the observation is the 2.5D Magnetohydrodynamics (MHD) model of solar wind acceleration by nonlinear waves, that predicts compressive fluctuations in coronal holes. To help identify the waves observed with the UVCS/WLC we model the propagation and dissipation of slow magnetosonic waves in polar plumes using 1D MHD code in spherical geometry, We find that the slow waves nonlinearly steepen in the gravitationally stratified plumes. The nonlinear steepening of the waves leads to enhanced dissipation due to compressive viscosity at the wave-fronts.
The properties of the synchrotron continuum are calculated on the basis of a steady spherically symmetric MHD model of the flow in the Crab nebula. The model accounts for the spectral distribution of the continuum from the optical to gamma-rays, and semiquantitatively for its observed spatial distribution. The model does not account for the radio portion of the Crab nebula's spectrum.
Observations of solar wind proton temperatures indicate that the solar wind is heated as it moves outward toward the orbit of Earth. This heating, which may be the results of electron heat conduction and perhaps MHD waves, has proven difficult to quantify and hence is often neglected in MHD models of the solar wind. An alternate approach to finding explicit heating terms for the MHD energy equation is to use a polytropic approximation. This paper discusses the properties of the polytropic approximation and its application to the solar wind plasma. By using data from the Helios 1 spacecraft, an empirical value for the polytropic index of the free-streaming solar wind is determined. Various corrections to the data are made to account for velocity, nonuniformity in radial sampling, and stream interaction regions. The polytropic index, as derived from proton data, is found to indepedent of speed state, within statistical error, and has an average value of 1.46. If magnetic pressure is included, the polytropic index has an average value of 1.58.
Recent analysis of the magnetic correlation function of solar wind fluctuations at 1 AU suggests the existence of two-component structure near the proton-cyclotron scale. Here we use two-and-one-half dimensional and three-dimensional compressible MHD models to look for two-component structure adjacent the proton-cyclotron scale. Our MHD system incorporates both Hall and Finite Larmor Radius (FLR) terms. We find that strong spectral anisotropies appear adjacent the proton-cyclotron scales depending on selections of initial condition and plasma beta. These anisotropies are enhancements on top of related anisotropies that appear in standard MHD turbulence in the presence of a mean magnetic field and are suggestive of one turbulence component along the inertial scales and another component adjacent the dissipative scales. We compute the relative strengths of linear and nonlinear accelerations on the velocity and magnetic fields to gauge the relative influence of terms that drive the system with wave-like (linear) versus turbulent (nonlinear) dynamics.
Context. Jets are dynamic, impulsive, well-collimated plasma events that develop at many different scales and in different layers of the solar atmosphere. Aims. Jets are believed to be induced by magnetic reconnection, a process central to many astrophysical phenomena. Within the solar atmosphere, jet-like events develop in many different environments, e.g. in the vicinity of active regions as well as in coronal holes, and at various scales, from small photospheric spicules to large coronal jets. In all these events, signatures of helical structure and/or twisting/rotating motions are regularly observed. The present study aims to establish that a single model can generally reproduce the observed properties of these jet-like events. Methods. In this study, using our state-of-the-art numerical solver ARMS, we present a parametric study of a numerical tridimensional magnetohydrodynamic (MHD) model of solar jet-like events. Within the MHD paradigm, we study the impact of varying the atmospheric plasma beta on the generation and properties of solar-like jets. Results. The parametric study validates our model of jets for plasma beta ranging from 10(sup 3) to 1, typical of the different layers and magnetic environments of the solar atmosphere. Our model of jets can robustly explain the generation of helical solar jet-like events at various beta less than or equal to 1. We show that the plasma beta modifies the morphology of the helical jet, explaining the different observed shapes of jets at different scales and in different layers of the solar atmosphere. Conclusions. Our results allow us to understand the energisation, triggering, and driving processes of jet-like events. Our model allows us to make predictions of the impulsiveness and energetics of jets as determined by the surrounding environment, as well as the morphological properties of the resulting jets.
Auto-magnetizing (AutoMag) liners are cylindrical tubes that employ helical current flow to produce strong internal axial magnetic fields prior to radial implosion on ~100 ns timescales. AutoMag liners have demonstrated strong uncompressed axial magnetic field production (>100 T) and remarkable implosion uniformity during experiments on the 20 MA Z accelerator. However, both axial field production and implosion morphology require further optimization to support the use of AutoMag targets in magnetized liner inertial fusion (MagLIF) experiments. Data from experiments studying the initiation and evolution of dielectric flashover in AutoMag targets on the Mykonos accelerator have enabled the advancement of magnetohydrodynamic (MHD) modeling protocols used to simulate AutoMag liner implosions. Implementing these protocols using ALEGRA has improved the comparison of simulations to radiographic data. Specifically, both the liner in-flight aspect ratio and the observed width of the encapsulant-filled helical gaps during implosion in ALEGRA simulations agree more closely with radiography data compared to previous GORGON simulations. Although simulations fail to precisely reproduce the measured internal axial magnetic field production, improved agreement with radiography data inspired the evaluation of potential design improvements with newly developed modeling protocols. Three-dimensional MHD simulation studies focused on improving AutoMag target designs, specifically seeking to optimize the axial magnetic field production and enhance the cylindrical implosion uniformity for MagLIF. Importantly, by eliminating the driver current prepulse and reducing the initial inter-helix gap widths in AutoMag liners, simulations indicate that the optimal 30–50 T range of precompressed axial magnetic field for MagLIF on Z can be accomplished concurrently with improved cylindrical implosion uniformity.
A two-dimensional planar MHD model of solar atmospheric transient events is used to simulate the flare-associated events of August 21, 1973. This event, observed in H-alpha, He II 304 A, soft X-ray, and coronal white light, provided sufficient information (especially in the latter diagnostic) for an assessment to be made of the model's ability to simulate major features of an actual solar event. It was found that a thermodynamic input pulse based on data provided by the NASA Marshall Space Flight Center-Aerospace Corporation X-ray telescope (S-056) on Skylab was sufficient to produce the global geometry, shock and contact-surface velocities, excess mass contours, and energy budget which were, for the most part, observed by the High Altitude Observatory white-light coronagraph (S-052) on Skylab in the form of a forerunner and coronal transient.
The modified gasdynamic convected-magnetic-field MHD model developed by Spreiter and Stahara (1980) to simulate the supersonic flow of the solar wind past planetary magnetoionospheres is extended to account for cases (such as Venus and comets) in which significant numbers of neutral atmospheric atoms become ionized in the surrounding flow and add to its momentum, energy, and mass. The mathematical model and the solution procedures for the nose and tail regions are explained; typical computational grids are shown; and numerical results for a comet and for the Venus ionosheath and bow shock are presented graphically. It is found that the bow shock weakens and moves further upstream of the obstacle as mass loading is increased and the flow upstream of the bow shock becomes more compressed.
Exploding Foil Initiator (EFI) systems function by launching a polymer flyer (often referred to as a slapper) at a high explosive (HE) pellet, which is intended to shock-initiate the HE. Flyer characteristics (e.g., shape, planarity, etc.) are extremely important in this respect, governing their ability to initiate the HE. Polymer models currently used in these systems are largely unvalidated, potentially compromising the integrity of magneto-hydrodynamic (MHD) predictions of EFI function. An improved EFI MHD modelling capability is expected to: (1) expedite optimization of EFI-based initiation system design, (2) enable consideration of EFI-related aging and corresponding lifetime predictions, and (3) potentially minimize the number of costly experiments required to certify EFI designs. To address this knowledge gap, several years ago we began a computational project to develop a higher-fidelity model for Kapton, a Polyimide copolymer most commonly employed as the flyer material in EFI devices.
Non-axisymmetric, flow-driven instabilities in the incompressible Hall-MHD model are studied in a differentially rotating cylindrical plasma. It is found that, in the Hall-MHD regime, both whistler waves and ion-cyclotron waves can extract energy from the flow shear, resulting in two distinct branches of global instability. The non-axisymmetric whistler modes grow significantly faster than non-axisymmetric, ideal MHD modes. A discussion of the global whistler instability mechanism is presented in the large-ion-skin-depth, “electron-MHD” limit. When the magnetic field is azimuthal, a subset of the whistler modes having zero axial wave number are uncovered to be destabilized by the “corotation amplifier” mechanism. It is observed that the effect of the Hall term on the non-axisymmetric modes can be appreciable when d i is on the order of a few percent of the width of the cylindrical annulus. Distinct global modes emerge in the strong Hall-MHD regime at significantly stronger magnetic fields than those required for unstable global MHD modes, as the Hall effect weakens the stabilizing “field-line bending” by decoupling ion motion from the magnetic field. These global non-axisymmetric modes may play an important role in weakly ionized accretion disks.
We present two generalized hybrid kinetic-Hall magnetohydrodynamics (MHD) models describing the interaction of a two-fluid bulk plasma, which consists of thermal ions and electrons, with energetic, suprathermal ion populations described by Vlasov dynamics. The dynamics of the thermal components are governed by standard fluid equations in the Hall MHD limit with the electron momentum equation providing an Ohm's law with Hall and electron pressure terms involving a gyrotropic electron pressure tensor. The coupling of the bulk, low-energy plasma with the energetic particle dynamics is accomplished through the current density (current coupling scheme; CCS) and the ion pressure tensor appearing in the momentum equation (pressure coupling scheme; PCS) in the first and the second model, respectively. The CCS is a generalization of two well-known models, because in the limit of vanishing energetic and thermal ion densities, we recover the standard Hall MHD and the hybrid kinetic-ions/fluid-electron model, respectively. This provides us with the capability to study in a continuous manner, the global impact of the energetic particles in a regime extending from vanishing to dominant energetic particle densities. The noncanonical Hamiltonian structures of the CCS and PCS, which can be exploited to study equilibrium and stability properties through the energy-Casimir variational principle, are identified. As a first application here, we derive a generalized Hall MHD Grad–Shafranov–Bernoulli system for translationally symmetric equilibria with anisotropic electron pressure and kinetic effects owing to the presence of energetic particles using the PCS.
The observations made during one unusual crossing of the magnetopause by the satellite Ogo 5 are compared with the classical reconnection model developed by Levy, Petschek, and Siscoe. The magnetic field observations appear to be generally consistent with this MHD model, although allowance must be made for the fact that the estimated magnetopause thickness was no more than 3.5 ion gyrodiameters. The nature of the finite gyroradius effects in such thin structures is discussed.
On May 11 and 12, 1999, when the solar wind proton density dropped to only about 3% of its nominal value, WIND and IMP 8 (and possibly Lunar Prospector) encountered the Earth's bow shock at unusually large distances. These observations, at 57-63 Re, correspond to the most distant crossings of the bow shock ever recorded. This paper will present observations from Interball-1 and Geotail, showing the bow shock near its nominal position, and from IMP 8 and WIND (and perhaps Lunar Prospector) crossing the bow shock at increasing radial distances. The bow shock is extremely weak by the time of its observation by WIND. However, a clear jump in both the magnetic field strength and solar wind bulk velocity is quite apparent, although with no significant changes in solar wind proton density. The nature of this dissipating shock will be discussed. Also, interball-1 has observed the magnetopause at an unusually large standoff distance. The observed motion of the magnetopause and bow shock will be compared with model predictions. This set of observations provide an excellent test case for the behavior of gas dynamic and MHD models under extreme conditions.
Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.
A disturbance or coronal mass ejection being advected by the solar wind will expand at the fastest local characteristic speed - typically approximately the fast-mode speed. To estimate this characteristic wave speed and the velocity field in the ambient corona, it is necessary to know the magnetic field, temperature, and density. Only the density is known from coronal observations. The temperature, magnetic field, and velocity are not yet directly measured in the outer corona and must be estimated from a model. In this study, it is estimated that the magnetic field, solar wind velocity, and characteristic speeds use the MHD model of coronal expansion between 1 and 5 solar radii (R solar radii) with a dipole magnetic field at the base. This model, for a field strength of about 2 gauss at the base, gives flow speeds at low latitudes (near the heliospheric current sheet) of 250 km/s at 5 R solar radii and, 50 km/s at 2 solar radii, and fast-mode speeds to 400 to 500 km/s everywhere between 2 and 5 solar radii. This suggests that the outer edge of a velocity of mass ejection reported by MacQueen and Fisher (1983) and implies that the acceleration mechanism for coronal mass ejections is other than simple entrainment in the solar wind.