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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2005

A Mechanism for the Loading-Unloading Substorm Cycle Missing in MHD Global Magnetospheric Simulation Models

Loading and consequent unloading of magnetic flux is an essential element of the substorm cycle in Earth's magnetotail. We are unaware of an available global MHD magnetospheric simulation model that includes a loading- unloading cycle in its behavior. Given the central role that MHD models presently play in the development of our understanding of magnetospheric dynamics, and given the present plans for the central role that these models will play in ongoing space weather prediction programs, it is clear that this failure must be corrected. A 2-dimensional numerical driven current-sheet model has been developed that incorporates an idealized current- driven instability with a resistive MHD system. Under steady loading, the model exhibits a global loading- unloading cycle. The specific mechanism for producing the loading-unloading cycle will be discussed. It will be shown that scale-free avalanching of electromagnetic energy through the model, from loading to unloading, is carried by repetitive bursts of localized reconnection. Each burst leads, somewhat later, to a field configuration that is capable of exciting a reconnection burst again. This process repeats itself in an intermittent manner while the total field energy in the system falls. At the end of an unloading interval, the total field energy is reduced to well below that necessary to initiate the next unloading event and, thus, a loading-unloading cycle results. It will be shown that, in this model, it is the topology of bursty localized reconnection that is responsible for the appearance of the loading-unloading cycle.

Klimas, A. J.↗

NASA NTRS1985

Current evolution in a numerical emerging-magnetic-flux model

The resistive-MHD equations are numerically solved in two-dimensions for an initial-boundary-value problem which models the emergence of magnetic flux from the photosphere into the corona. As the emergence begins a current sheet forms around the emerging region which separates the emerging region from the overlying coronal magnetic field. This current sheet is the source of the free-magnetic energy in the system, and in the limit of zero resistivity it is a simple tangential discontinuity. However, when the resistivity is finite, reconnection between the magnetic field in the emerging region and the overlying coronal magnetic field ensures, and the subsequent evolution of the enveloping current sheet becomes complex. The overall time history of the current evolution is suggestive of the expected current evolution for the pre-flare, impulsive, and main phases of flares.

Forbes, T.↗

NASA NTRS1994

Influence of plasma sheet thickening toward the tail flanks on magnetotail stability and dynamics

On the basis of resistive MHD simulations we investigate the stability and dynamic evolution of the magnetotail for configurations with various degrees of plasma sheet thickening from midnight toward the tail flanks, associated with a similar increase of the B(sub z) component. This increase is varied by factors between 2 and 3.5, while the magnitude of B(sub z) at midnight is left unchanged. The increase of B(sub z) has a strong effect on the spatial extend of resistive tearing and reconnection across the tail and on the unstable evolution. For stronger thickening the cross-tail extent of the near-Earth neutral line, formed by reconnection, and correspondingly the amount of reconnected magnetic flux get reduced and the dipolarization effects earthward of the neutral line are more concentrated near midnight. An increase of B(sub z) from midnight toward the flanks by more than a factor of about 3 can possibly even suppress the resistive tearing instability. This indicates the possibility of a stability transition due to a reduction of B(sub z) in the flank regions in response to changes in the solar wind, even when the magnitude of B(sub z) near midnight and the fluctuation level in that region are unchanged (provided that a sufficient fluctuation level exists).

Birn, Joachim↗

NASA NTRS2016

Separator Reconnection at the Magnetopause for Predominantly Northward and Southward IMF: Techniques and Results

In this work, we demonstrate how to track magnetic separators in three-dimensional simulated magnetic fields with or without magnetic nulls, apply these techniques to enhance our understanding of reconnection at the magnetopause. We present three methods for locating magnetic separators and apply them to 3-D resistive MHD simulations of the Earth's magnetosphere using the Block-Adaptive-Tree Solar-wind Roe-type Upwind Scheme code. The techniques for finding separators and determining the reconnection rate are insensitive to interplanetary magnetic field (IMF) clock angle and can in principle be applied to any magnetospheric model. Moreover, the techniques have a number of advantages over prior separator finding techniques applied to the magnetosphere. The present work examines cases of high and low resistivity for two clock angles. We go beyond previous work examine the separator during Flux Transfer Events (FTEs). Our analysis of reconnection on the magnetopause yields a number of interesting conclusions: Reconnection occurs all along the separator even during predominately northward IMF cases. Multiple separators form in low-resistivity conditions, and in the region of an FTE the separator splits into distinct branches. Moreover, the local contribution to the reconnection rate, as determined by the local parallel electric field, drops in the vicinity of the FTE with respect to the value when there are none.

modeling↗

NASA NTRS1984

Numerical simulation of reconnection in an emerging magnetic flux region

Numerical solutions in two dimensions are presented for the resistive MHD equations of an initial boundary value problem, simulating reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. This numerical model displays four main phases, which are interpreted in terms of (1) a slowly evolving quasi-steady phase during which most of the magnetic flux emerges, with approximate equilibrium between magnetic and pressure forces; (2) an impulsive phase in which either the reconnection or continuing emergence of the first phase comes to disrupt the equilibrium, and extensive acceleration occurs as the high pressure region within the emerging region drives the fluid upwards and outwards; (3) a second quasi-steady phase; and (4) a potential-static phase in which continuing reconnection, ohmic dissipation, and fluid transport through boundaries depletes the system of all currents and flows.

Forbes, T. G.↗

NASA NTRS1986

A shock condensation mechanism for loop prominences

Self-consistent, numerical solutions of the resistive MHD equations in two dimensions show that a quasi-stationary, fast-mode shock is a characteristic feature of the reconnection dynamics of the Kopp-Pneuman model of two-ribbon flares. A preliminary analysis of the effects of radiative cooling and thermal conduction suggests that the fast shock can help trigger a thermal condensation (i.e., a loop prominence)if the reconnecting magnetic fields are sufficiently strong.

Forbes, T. G.↗

NASA NTRS1986

Can prominences form in current sheets

Two-dimensional numerical simulations of the formation of cold condensations in a vertical current sheet have been performed using the radiative, resistive MHD equations with line-tied boundary conditions at one end of the sheet. Prominence-like condensations are observed to appear above and below an X-line produced by the onset of the tearing-mode instability. Cooling in the sheet is initiated by Ohmic decay, with the densest condensations occurring in the region downstream of a fast-mode shock. This shock, which is due to the line-tied boundary conditions, terminates one of the two supermagnetosonic reconnection jets that develop when the tearing is fully developed. The condensation properties of shock waves, which may trigger or considerably enhance the conditions for thermal condensation are emphasized.

Malherbe, J. M.↗

NASA NTRS1986

On the influence of an external electric field on magnetotail reconnection

Magnetotail reconnection is studied under the influence of a dawn to dusk electric field applied at the high-latitude tail boundary using a two-dimensional resistive MHD code. It is found that the basic unstable sequence of tearing and plasmoid formation and ejection was not changed if the external electric field was fairly uniform. The externally forced compression, however, could shorten the onset time and the growth time of the evolution. The faster growth can be understood by the fact that the forced compression gradually alters the lobe field strength and the plasma sheet width and reduces thereby the relevant Alfven time. In fact, if the evolution of the cases with imposed electric field is rescaled on the basis of these characteristic parameters, a good quantitative agreement with the case without external electric field is found.

Birn, J.↗

NASA NTRS1990

Buoyancy-driven instabilities and the nonlinear breakup of a sheared magnetic layer

Motivated by problems concerning the storage and subsequent escape of the solar magnetic field, a study has been made of how a magnetic layer embedded in a convectively stable atmosphere evolves due to axisymmetric instabilities driven by magnetic buoyancy. The initial equilibrium consists of a toroidal field sheared by a weaker poloidal component. The linear stability problem is investigated for both ideal and resistive MHD, and the nonlinear evolution is followed by numerical integration of the equations of motion. In all cases, the instability is greatly affected by the distribution and strength of the poloidal field. In particular, both the horizontal and vertical scales of the motions are controlled by the location of the surface on which the poloidal field vanishes: the resonant surface. In the nonlinear regime, a resonant surface close to the interface between the magnetized and field-free fluid leads to the localization of the instability, so that only a fraction of the magnetic region is disrupted by the motions. By contrast, a deeply seated resonant surface leads to the complete disruption of the layer and to the formation of large, helical magnetic fragments whose identity is preserved for the entire simulation.

Cattaneo, Fausto↗

NASA NTRS1995

Physics of intermediate shocks: A review

Intermediate shocks (ISs) lead to a transition from super-Alfvenic to sub-Alfvenic flow and are different from slow and fast shocks in that an IS rotates the component of the magnetic field tangent to the shock plane by 180 deg. Another peculiarity of ISs is that for the same upstream conditions an IS can have two different downstream states. There also exist a second class of ISs which rotate the magnetic field by an angle other than 180 deg. Due to their noncoplanar nature they cannot be time-stationary and are referred to as time-dependent intermediate shocks (TDIS). The existence of ISs has been the subject of much controversy over the years. Early studies questioned the physical reality of ISs. However, the studies of ISs found a new impetus when C.C. Wu showed that ISs do exist and are stable within the resistive MHD framework. In this paper, after a brief historical overview of the subject, we will review the latest developments in the study of ISs. In particular, we will address the questions of stability and structure of ISs and the relationship between ISs and other discontinuities. One of the recent developments has been the finding that ISs can be unsteady, reforming in time. Details of this process will be discussed. Finally, we examine the effect of anisotropy on the resolutions and discuss the relevance of ISs to the observed field rotations at the Earth's magnetopause.

Karimabadi, H.↗

NASA NTRS1995

Heating of coronal holes by the resonant absorption and dissipation of Alfven waves and its relation to solar wind acceleration

Coronal hole regions are well known sources of high-speed solar wind, however to account for the observed properties of the solar wind a source of momentum and heat must be included. Alfven waves were suggested as the possible source of heating that accelerates the solar wind. We investigate the propagation of the Alfven waves in coronal holes via numerical solution of the linearized 2-D resistive MHD equations in slab geometry. The Alfven waves are driven at the lower boundary of the coronal hole and propagate into the corona. The waves are reflected at the coronal hole boundary and part of the wave energy leaks out of the coronal hole. We compare the calculated wavelengths and the attenuation rate of the fast mode Alfven waves in the leaky waveguide formed by the coronal hole with the analytical ideal MHD solutions. The formation of resonance heating layers is found to occur when shear Alfven waves propagate in an inhomogeneous coronal hole. The heating is enhanced when fast mode waves couple to the shear Alfven waves. The narrow heating layers are formed near the location of the ideal resonance, which might occur near the coronal hole boundary for a nearly constant density coronal hole, surrounded by a higher density plasma. We investigate the dependence of the heating on the driver frequency, the Lundquist number, and on the heliocentric distance. and find that the low frequency Alfven waves can be an efficient source of heating at large distances from the Sun. We discuss the relation of our results to the observed properties of high-speed solar wind and coronal holes.

Ofman, L.↗

NASA NTRS2005

Intermittent Turbulence and SOC Dynamics in a 2-D Driven Current-Sheet Model

Borovsky et al. have shown that Earth's magnetotail plasma sheet is strongly turbulent. More recently, Borovsky and Funsten have shown that eddy turbulence dominates and have suggested that the eddy turbulence is driven by fast flows that act as jets in the plasma. Through basic considerations of energy and magnetic flux conservation, these fast flows are thought to be localized to small portions of the total plasma sheet and to be generated by magnetic flux reconnection that is similarly localized. Angelopoulos et al., using single spacecraft Geotail data, have shown that the plasma sheet turbulence exhibits signs of intermittence and Weygand et al., using four spacecraft Cluster data, have confirmed and expanded on this conclusion. Uritsky et al., using Polar UVI image data, have shown that the evolution of bright, nightside, UV auroral emission regions is consistent with many of the properties of systems in self-organized criticality (SOC). Klimas et al. have suggested that the auroral dynamics is a reflection of the dynamics of the fast flows in the plasma. sheet. Their hypothesis is that the transport of magnetic fludenergy through the magnetotail is enabled by scale-free avalanches of localized reconnection whose SOC dynamics are reflected in the auroral UV emission dynamics. A corollary of this hypothesis is that the strong, intermittent, eddy turbulence of the plasma sheet is closely related to its critical dynamics. The question then arises: Can in situ evidence for the SOC dynamics be found in the properties of the plasma sheet turbulence? A 2-dimensional numerical driven current-sheet model of the central plasma sheet has been developed that incorporates an idealized current-driven instability with a resistive MHD system. It has been shown that the model can evolve into SOC in a physically relevant parameter regime. Initial results from a study of intermittent turbulence in this model and the relationship of this turbulence to the model's known SOC dynamics will be discussed.

Klimas, A. J.↗

NASA NTRS2010

Impact of Magnetic Draping, Convection, and Field Line Tying on Magnetopause Reconnection Under Northward IMF

We simulate a northward IMF cusp reconnection event at the magnetopause using the OpenGGCM resistive MHD code. The ACE input data, solar wind parameters, and dipole tilt belong to a 2002 reconnection event observed by IMAGE and Cluster. Based on a fully three-dimensional skeleton separators, nulls, and parallel electric fields, we show magnetic draping, convection, ionospheric field line tying play a role in producing a series of locally reconnecting nulls with flux ropes. The flux ropes in the cusp along the global separator line of symmetry. In 2D projection, the flux ropes the appearance of a tearing mode with a series of 'x's' and 'o's' but bearing a kind of 'guide field' that exists only within the magnetopause. The reconnecting field lines in the string of ropes involve IMF and both open and closed Earth magnetic field lines. The observed magnetic geometry reproduces the findings of a superposed epoch impact parameter study derived from the Cluster magnetometer data for the same event. The observed geometry has repercussions for spacecraft observations of cusp reconnection and for the imposed boundary conditions reconnection simulations.

Wendel, Deirdre E.↗

NASA NTRS2013

Computational Studies of Magnetic Nozzle Performance

An extensive literature review of magnetic nozzle research has been performed, examining previous work, as well as a review of fundamental principles. This has allow us to catalog all basic physical mechanisms which we believe underlie the thrust generation process. Energy conversion mechanisms include the approximate conservation of the magnetic moment adiabatic invariant, generalized hall and thermoelectric acceleration, swirl acceleration, thermal energy transformation into directed kinetic energy, and Joule heating. Momentum transfer results from the interaction of the applied magnetic field with currents induced in the plasma plume., while plasma detachment mechanisms include resistive diffusion, recombination and charge exchange collisions, magnetic reconnection, loss of adiabaticity, inertial forces, current closure, and self-field detachment. We have performed a preliminary study of Hall effects on magnetic nozzle jets with weak guiding magnetic fields and weak expansions (p(sub jet) approx. = P(sub background)). The conclusion from this study is that the Hall effect creates an azimuthal rotation of the plasma jet and, more generally, creates helical structures in the induced current, velocity field, and magnetic fields. We have studied plasma jet expansion to near vacuum without a guiding magnetic field, and are presently including a guiding magnetic field using a resistive MHD solver. This research is progressing toward the implementation of a full generalized Ohm's law solver. In our paper, we will summarize the basic principle, as well as the literature survey and briefly review our previous results. Our most recent results at the time of submittal will also be included. Efforts are currently underway to construct an experiment at the University of Michigan Plasmadynamics and Electric Propulsion Laboratory (PEPL) to study magnetic nozzle physics for a RF-thruster. Our computational study will work directly with this experiment to validate the numerical model, in order to study magnetic nozzle physics and optimize magnetic nozzle design. Preliminary results from the PEPL experiment will also be presented.

Ebersohn, Frans H.↗

NASA NTRS2022

Impact of 3D Structure on Magnetic Reconnection

Results from 2.5D and 3D studies of the onset and development of the tearing instability are presented, using high fidelity resistive MHD simulations. A limited parameter study of the strength of the reconnecting field (or shear angle) was performed. An initially simple 1D equilibrium was used, consisting of a modified force-free current sheet, with periodic boundary conditions in all directions. In all cases, the linear and nonlinear evolution led to a primary current sheet between two large flux ropes. The global reconnection rate during this later stage was analyzed in all simulations. It was found that in 2.5D the primary current sheet fragmented owing to plasmoids, and as expected, the global reconnection rate, calculated using multiple methods, increases with the strength of the reconnecting field owing to a stronger Alfvén speed. In 3D, the presence of interacting oblique modes of the tearing instability complicates the simple 2.5D picture, entangling the magnetic field of the inflow and introducing a negative effect on the reconnection rate. The two competing effects of stronger Alfvén speed and entangling, which both increase with the strength of the reconnecting field, resulted in a decrease in the reconnection rate with increasing reconnecting field. For all simulations, the 3D rates were less than in 2.5D but suggest that as one goes to weak reconnecting field (or strong guide field) the system becomes more 2.5D-like and the 2.5D and 3D rates converge. These results have relevance to situations like nanoflare heating and flare current sheets in the corona.

Lars K S Daldorff↗

NASA NTRS2005

A Numerical Study of Resistivity and Hall Effects for a Compressible MHD Model

The effect of resistive, Hall, and viscous terms on the flow structure compared with compressible ideal MHD is studied numerically for a one-fluid non-ideal MHD model. The goal of the present study is to shed some light on the emerging area of non-ideal MHD modeling and simulation. Numerical experiments are performed on a hypersonic blunt body flow with future application to plasma aerodynamics flow control in reentry vehicles. Numerical experiments are also performed on a magnetized time-developing mixing layer with possible application to magnetic/turbulence mixing.

Yee, H. C.↗

NASA NTRS1995

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗

NASA NTRS1994

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗