Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “resistive MHD”

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 19 records

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

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

Effects of flows on viscous and resistive MHD stability

In many solar applications the viscosity appears to be more important than resistivity. In order to discuss the instabilities in solar conditions, an idealized configuration is considered in which the plasma is flowing in the z-direction along the magnetic field B sub 0 with a velocity V sub 0. As far as the velocity is concerned two different velocity profiles, with different hydrodynmaic stability properties are discussed. The results are summarized.

Einaudi, Giorgio↗

Numerical simulation of reconnection in an emerging magnetic flux region

The resistive MHD equations are numerically solved in two dimensions for an initial-boundary-value problem which simulates reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. The solution involves both ideal-MHD and resistive-MHD processes, and the solution shows an evolution which is remarkably suggestive of the preflare, impulsive, and main phases of the flare-cycle.

Forbes, T. C.↗

Magnetic Reconnection in a Compressible MHD Plasma

Using steady-state resistive MHD, magnetic reconnection is reinvestigated for conditions of high resistivity/low magnetic Reynolds number, when the thickness of the diffusion region is no longer small compared to its length. Implicit expressions for the reconnection rate and other reconnection parameters are derived based on the requirements of mass, momentum, and energy conservation. These expressions are solved via simple iterative procedures. Implications specifically for low Reynolds number/high resistivity are being discussed

Hesse, Michael↗

The Role of Self-Organized Criticality in the Substorm Phenomenon and its Relation to Localized Reconnection in the Magnetospheric Plasma Sheet

Evidence is presented that suggests there is a significant self-organized criticality (SOC) component in the dynamics of substorms in the magnetosphere. Observations of BBFs, fast flows, localized dipolarizations, plasma turbulence, etc. are taken to show that multiple localized reconnection sites provide the basic avalanche phenomenon in the establishment of SOC in the plasma sheet. First results are presented from a continuing plasma physical study of this avalanche process. A one-dimensional resistive MHD model of a magnetic field reversal is discussed. Resistivity, in this model, is self-consistently generated in response to the excitation of an idealized current-driven instability. When forced by convection of magnetic flux into the field reversal region, the model yields rapid magnetic field annihilation through a dynamic behavior that is shown to exhibit many of the characteristics of SOC. Over a large range of forcing strengths, the annihilation rate is shown to self-adjust to balance the rate at which flux is convected into the reversal region. Several analogies to magnetotail dynamics are discussed: (1) It is shown that the presence of a localized criticality in the model produces a remarkable stability in the global configuration of the field reversal while simultaneously exciting extraordinarily dynamic internal evolution. (2) Under steady forcing, it is shown that a loading-unloading cycle may arise that, as a consequence of the global stability, is quasi-periodic and, therefore, predictable despite the presence of internal turbulence in the field distribution. Indeed, it is shown that the global loading-unloading cycle is a consequence of the internal turbulence. (3) It is shown that, under steady, strong forcing the loading-unloading cycle vanishes. Instead, a recovery from a single unloading persists indefinitely. The field reversal is globally very steady while internally it is very dynamic as field annihilation goes on at the rate necessary to match the strong forcing. From this result we speculate that steady magnetospheric convection events result when the plasma sheet has been driven close to criticality over an extended spatial domain. During these events, we would expect to find localized reconnection sites distributed over the spatial domain of near criticality and we would expect to find plasma sheet transport in that domain to be closely related to that of BBF and fast flow events.

Klimas, Alex J.↗

Three-dimensional MHD modeling of magnetotail dynamics for different polytropic indices

The paper analyzes the 3D resistive tearing instability in a magnetotail configuration by means of a resistive MHD code for various forms of the energy equation, including ohmic and compressional heating. A complete derivation of the set of MHD equations from a kinetic formalism is presented, with all simplifying assumptions noted. Emphasis is placed on the changes in magnetotail evolution caused by a variation of the assumed value of the polytropic exponent gamma, representing the ratio of specific heat in the energy equation which closes the set of MHD equations. Within the gamma range of 2/3-4, a tearing mode is found to develop in all cases. Certain properties, such as the formation and ejection of a plasmoid, associated with the formation of the reversal region of the north-south component of the magnetic field, and strong earthward and tailward flow fields prevail as qualitative features independent of the value of gamma chosen. Ohmic dissipation is found to be the major source of the heating, while adiabatic or nonadiabatic compression may contribute primarily in the closed field region earthward of the reconnection site.

Hesse, Michael↗

Effect of a localized minimum in equatorial field strength on resistive tearing instability in the geomagnetotail

A two-dimensional, resistive-MHD computer code is used to investigate the spontaneous reconnection of magnetotaillike configurations. The initial conditions adopted in the simulations are of two types: (1) in which the equatorial normal magnetic field component B(ze) declines monotonically down the tail, and (2) in which B(ze) exhibits a deep minimum in the near-earth plasma sheet. To represent the case where the earthward convection stops before the X line forms, zero-flow boundary conditions are imposed at the edges of the computational box. The initial configurations are in equilibrium and table within ideal MHD. The dynamic evolution of the system starts after the resistivity is turned on. The main results of these simulations basically support the neutral-line model of substorms and confirm Birn's (1980) computer studies. Spontaneous formation of an X-type neutral point and a single O-type plasmoid with strong tailward flow on the tailward side of the X point is found. The time interval from the turning on of the resistivity to the formation of a plasmoid is much shorter in the case where there is an initial deep minimum. A simple analytic calculation is also carried out to demonstrate why the configuration with a deep minimum is more susceptible to the development of the neutral point.

Hau, L.-N.↗

Resistive ballooning modes in line-tied coronal arcades

The equations describing the linear evolution of resistive ballooning modes are obtained by using a modified WKB expansion in the short perpendicular wavelength, while variations of the perturbations along the field are described by a slowly varying amplitude, on which the tying boundary conditions are imposed. In general, given an equilibrium, there are certain ranges of magnetic surfaces for which the system predicts instability even without dissipation. The main conclusion is that within the resistive MHD approximation cylindrically symmetric arcades with pressure falling with radius are unstable to resistive localized modes; the growth rates, close to ideal marginal stability, are large, so that it would appear that energy could be released during 10 to 100 Alfven times. The wavelength of the modes is expected to be limited by the ion gyroradius, when stabilizing drift effects must be taken into account. The nonlinear evolution of resistive ballooning modes should be studied to assess their overall relevance to the violent and rapidly evolving phenomena observed on the sun.

Velli, Marco↗

Does the Hall Effect Solve the Flux Pileup Saturation Problem?

It is well known that magnetic flux pileup can significantly speed up the rate of magnetic reconnection in high Lundquist number resistive MHD,allowing reconnection to proceed at a rate which is insensitive to the plasma resistivity over a wide range of Lundquist number. Hence, pileup is a possible solution to the Sweet-Parker time scale problem. Unfortunately, pileup tends to saturate above a critical value of the Lundquist number, S_c, where the value ofS_c depends on initial and boundary conditions, with Sweet-Parker scaling returning above S_c. It has been argued (see Dorelli and Bim [2003] and Dorelli [2003]) that the Hall effect can allow flux pileup to saturate (when the scale of the current sheet approaches ion inertial scale, di) before the reconnection rate begins to stall. However, the resulting saturated reconnection rate, while insensitive to the plasma resistivity, was found to depend strongly on the di. In this presentation, we revisit the problem of magnetic island coalescence (which is a well known example of flux pileup reconnection), addressing the dependence of the maximum coalescence rate on the ratio of di in the "large island" limit in which the following inequality is always satisfied: l_eta di lambda, where I_eta is the resistive diffusion length and lambda is the island wavelength.

Dorelli, John C.↗

Magnetosphere-ionosphere coupling during plasmoid evolution - First results

The influence of magnetosphere-ionosphere coupling on the dynamic evolution of the magnetotail is investigated by a three-dimensional resistive MHD code that includes the effects of the closure of field-aligned currents in a simple resistive model ionosphere. Particular emphasis is on the role of this coupling during substorm evolution and the modification of the latter by the convection driven by the ionospheric electric fields. For comparison, results are presented from a simulation which uses an infinitely conducting ionosphere but is otherwise identical. Comparison of the two simulations shows that the major impact of magnetosphere-ionospheric communication is an acceleration of magnetotail evolution. Otherwise, phenomena in the two models are qualitatively similar. It is concluded that ionospheric effects do not significantly affect substorm associated magnetotail dynamics.

Hesse, Michael↗

MHD simulations of magnetic reconnection in a skewed three-dimensional tail configuration

The dynamic evolution of a nonsymmetric magnetotail configuration initiated by the sudden occurrence of (anomalous) resistivity are examined using a three-dimensional resistive MHD simulation developed by Birn (1990) that includes a net cross-tail field and thus breaks the mirror symmetry around the neutral sheet. Results show that the field evolution is similar to that of a symmetric configuration studied by Birn and Hones (1981), pointing to the formation and ejection of a plasmoid. On the other hand, the topological structure of the magnetic field, defined by the field line connections, was remarkably different from the symmetric case. The plasmoid in this case became 'open', connected initially with the earth but gradually becoming connected with the interplanetary field. The openness of the plasmoid and its magnetic connection with interplanetary field lines suggest the possibility of a heat flux out of the plasmoid on interconnected flux tubes.

Birn, J.↗

Field-aligned plasma flow in MHD simulations of magnetotail reconnection and the formation of boundary layers

The outflow from a reconnection region in a realistic magnetotail geometry is studied using a two-dimensional time-dependent, compressible, resistive MHD code. Two cases are emphasized: (1) the evolution of near-earth reconnection, which grows in the form of an internal unstable mode after gradual externally forced changes have initiated an anomalous dissipation process, and (2) the evolution of more distant reconnection under influence of a nonuniform inflow that forces reconnection to occur at a given location in the distant tail. In both cases, it is demonstrated that plasma flow is primarily parallel to the magnetic field in regions away from the localized area of reconnection and outside a narrow central layer.

Birn, J.↗

The Role of the Hall Effect in Global Structure and Dynamics of Planetary Magnetospheres: Ganymede as a Case Study

We present high-resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall J x B forces accelerate ions to the local Alfven speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian 'double magnetopause' structure, an Alfvenic ion jet extending across the upstream magnetopause, and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves.

Magnetic Reconnection↗

Resistive wave dissipation on magnetic inhomogeneities Normal modes and phase mixing

Numerical solutions of the linearized, resistive MHD equations indicate the presence of two characteristic forms of wave decay on magnetic inhomogeneities. Solutions for relatively long-wavelength disturbances and large values of the resistivity have the behavior of decaying normal modes. In the opposite limits, phase mixing becomes important, with the eventual build-up of large spatial gradients in which the resistive terms dominate. This composite behavior is shown to be conceptually consistent with analytic results, which predict that the complete solution consists of the sum of collective and noncollective contributions. The numerical simulations go beyond analytic theory by defining where each contribution prevails and computing wave decay when phase mixing is effective. The dispersion relation in the normal-mode regime is also determined analytically, with some approximations, and is in good agreement with the numerical solution. The decay time of the normal modes varies with Lundquist number at S1/6, while the phase-mixed decay time scales as S1/3.

Steinolfson, R. S.↗

Advanced MHD Algorithm for Solar and Space Science: lst Year Semi Annual Progress Report

We report progress for the development of MH4D for the first and second quarters of FY2004, December 29, 2002 - June 6, 2003. The present version of MH4D can now solve the full viscous and resistive MHD equations using either an explicit or a semi-implicit time advancement algorithm. In this report we describe progress in the following areas. During the two last quarters we have presented poster at the EGS-AGU-EUG Joint Assembly in Nice, France, April 6-11, 2003, and a poster at the 2003 International Sherwood Theory Conference in Corpus Christi, Texas, April 28-30 2003. In the area of code development, we have implemented the MHD equations and the semi-implicit algorithm. The new features have been tested.

Schnack, Dalton D.↗

Nonlinear evolution of resistive tearing mode instability with shear flow and viscosity

The effect of shear flow on the nonlinear evolution of the tearing mode is investigated via numerical solutions of the resistive MHD equations in slab geometry, using a finite-difference alternative-direction implicit method. It was found that, when the shear flow is small (V less than 0.3), the tearing mode saturates within one resistive time, whereas for larger flows the nonlinear saturation develops on longer time scales. The magnetic energy release decreases and the saturation time increases with increasing values of V for both small and large resistivity. Shear flow was found to decrease the saturated magnetic island width and to generate currents far from the tearing layer. Results suggest that equilibrium shear flow may improve the confinement of tokamak plasma.

Ofman, L.↗

Magnetic Field Effects on Plasma Plumes

Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results

Ebersohn, F.↗