Resonant MHD Waves in the Solar Atmosphere

The linear theory of MHD resonant waves in inhomogeneous plasmas is reviewed. The review starts from discussing the properties of driven resonant MHD waves. The dissipative solutions in Alfvén and slow dissipative layers are presented. The important concept of connection formulae is introduced. Next, we proceed on to non-stationary resonant MHD waves. The relation between quasi-modes of ideal MHD and eigenmodes of dissipative MHD are discussed. The solution describing the wave motion in non-stationary dissipative layers is given. It is shown that the connection formulae remain valid for non-stationary resonant MHD waves. The initial-value problem for resonant MHD waves is considered. The application of theory of resonant MHD waves to solar physics is discussed.     

Why space physics needs to go beyond the MHD box

Magnetohydrodynamic (MHD) theory has been used in space physics for more than forty years, yet many important questions about space plasmas remain unanswered. We still do not understand how the solar wind is accelerated, how mass, momentum and energy are transported into the magnetosphere and what mechanisms initiate substorms. Questions have been raised from the beginning of the space era whether MHD theory can describe correctly space plasmas that are collisionless and rarely in thermal equilibrium. Ideal MHD fluids do not induce electromotive force, hence they lose the capability to interact electromagnetically. No currents and magnetic fields are generated, rendering ideal MHD theory not very useful for space plasmas. Observations from the plasma sheet are used as examples to show how collisionless plasmas behave. Interpreting these observations using MHD and ideal MHD concepts can lead to misleading conclusions. Notably, the bursty bulk flows (BBF) with large mean velocities left(〈 v 〉≥400 km s right) that have been interpreted previously as E×B flows are shown to involve much more complicated physics. The sources of these nonvanishing 〈 v 〉 events, while still not known, are intimately related to mechanisms that create large phase space gradients that include beams and acceleration of ions to MeV energies. The distributions of these nonvanishing 〈 v 〉 events are associated with large amplitude variations of the magnetic field at frequencies up to and exceeding the local Larmor frequency where MHD theory is not valid. Understanding collisionless plasma dynamics such as substorms in the plasma sheet requires the self-consistency that only kinetic theory can provide. Kinetic modeling is still undergoing continual development with many studies limited to one and two dimensions, but there is urgent need to improve these models as more and more data show kinetic physics is fundamentally important. Only then will we be able to make progress and obtain a correct picture of how collisionless plasmas work in space.     

Streamer Evaporation

Streamer evaporation is the consequence of heating in ideal MHD models because plasma is weakly contained by the magnetic field. Heating causes inflation, opening of field lines, and release of solar wind. It was discovered in simulations and, due to the absence of loss mechanisms, the ultimate end point is the complete evaporation of the streamer. Of course streamers do not behave in this way because of losses by thermal conduction and radiation. Heating is also expected to depend on ambient conditions. We use a global MHD model with thermal conduction to examine the effect of changing the heating scale height. We also extend an analytic model of streamers developed by Pneuman (1968) to show that steady streamers are unable to contain plasma for temperatures near the cusp greater than ∼ 2 × 10⁶ K. This revised version was published online in June 2006 with corrections to the Cover Date.     

三维磁流体强化超燃冲压发动机数值模拟

建立了三维磁流体强化超燃冲压发动机内部黏性流场的求解模型.针对马赫数为6设计了联合应用磁控进气道和磁流体能量旁路的磁流体强化超燃冲压发动机模型.针对该模型进行了数值模拟研究,分析其中的三维流场结构、电参数分布规律以及能量转换特性.结果表明:当飞行马赫数为8时,磁控进气道的应用能够使头部压缩激波回到唇口,使分离区消失,内进气道中的流动恢复到设计状态.磁流体能量旁路可有效降低燃烧室入口处的马赫数,从而改善发动机性能.其中发生器中的流动参数和电参数的分布比较理想,效果显著;而加速器要取得显著的加速效果则需要大量的能量输入.在加速器中,电极附近焦耳耗散严重,导致局部高温以及流动的复杂性,性能不够理想.      

基于HLLC格式的理想MHD方程数值模拟

针对理想MHD(magnetohydrodynamics)方程数值求解困难的问题,基于原始的HLLC (Harten-Lax-Van Leer Contact wave)近似黎曼解方法,发展出一种新的适用于理想MHD问题的通量计算格式。控制方程采用有限体积法离散,时间推进采用隐式的LU-SGS（lower-upper symmetric Gauss-Seidel）格式,并且采用双曲型散度清除技术来抑制磁场散度的累积。通过一维激波管问题的数值模拟表明,HLLC-MHD格式能准确的分辨并捕捉复杂的磁流体力学波系结构,保证数值计算的精度,最大的数值计算误差不超过10%;通过二维的Rotor问题的数值模拟研究表明,HLLC-MHD格式能够应用于多维理想MHD问题的数值模拟,并且能够准确捕捉多维情况下磁流场中的阿尔文波;通过对比有、无散度清除的Rotor问题计算结果表明,双曲型散度清除技术可以将磁场散度峰值从50降低到2,有效抑制高磁场散度区域的散度累积,同时也会将误差传播到低散度区域,并且引起边界处散度的累积,影响计算的稳定性。      

超声速磁流体加速实验及一维模型分析

为了获得负载系数、电导率等参数变化对超声速磁流体加速效果的影响规律,利用激波风洞,采用氩气与碳酸钾作为工质,电容提供电能的方式,在磁感应强度为0.5T的条件下,进行了不同电容充电电压下的超声速磁流体加速实验研究,并对一维定常理想分段法拉第型磁流体加速模型进行了分析.通过实验获得了不同电容充电电压下#10电极间的电压、电流、负载系数、电导率及#20电极开路电压等数据,在300, 400V电容充电电压下,气流速度分别增加11.4%和24.0%,在500V电容充电电压下气流速度减小11.1%.实验及模型分析得出不同的负载系数会使超声速磁流体处于加速或减速的不同状态,而电导率会影响注入总能量的大小,使磁流体流动的速度梯度大小发生改变.      

三维磁流体动力学管道流动加减速控制数值研究

阐述了磁流体动力学(MHD)控制流场作用机理.在小磁雷诺数条件下,运用数值模拟方法,对不同外加电磁场条件下三维MHD管道流动的流场情况进行研究,得出不同磁场、电场等MHD参数作用下MHD加速器的性能.对MHD加速器的应用前景进行了展望.计算结果显示,电场取6000V/m,磁场为0.92T 情况下,流场的加速性能可达13.46%,并且可以通过减小磁场或增大电场进一步提高加速性能;在仅添加0.92T磁场条件下,速度减速可达16.35%.     

MHD ANALYSIS OF PETSCHEK-TYPE RECONNECTION IN NON-UNIFORM FIELD AND FLOW GEOMETRIES

Analytical studies of reconnection have, for the most part, been confined to steady and uniform current sheet geometries. In contrast to these implifications, natural phenomena associated with the presence of current sheets indicate highly non-uniform structure and time-varying behaviour. Examples include the violent outbursts of energy on the Sun known as solar flares, and magnetospheric phenomena such as flux transfer events, plasmoids, and auroral activity. Unlike the theoretical models, reconnection therefore occurs in a highly dynamic and structured plasma environment. In this article we review the mathematical tools and techniques which are available to formulate models capable of describing the effects of reconnection in such situations. We confine attention to variants of the reconnection model first discussed by Petschek in the 1960s, in view of its successful application in predicting and interpreting phenomena in the terrestrial magnetosphere. The analysis of Petschek-type reconnection is based on the equations of ideal magnetohydrodynamics (MHD), which describe the large-scale behaviour of the magnetic field and plasma flow outside the diffusion region, which we determine as a localised part of the current sheet in which reconnection is initiated. The approach we adopt here is to transform the MHD equations into a Lagrangian or so-called 'frozen-in' coordinate system. In this coordinate system, the equation of motion transforms into a set of coupled nonlinear equations, in which the presence of inhomogeneous magnetic fields and/or plasma flows gives rise to a term similar to that which appears in the study of the ordinary string equation in a non-homogeneous medium. As demonstrated here, this approach not only clarifies and highlights the effects of such non-uniformities, it also simplifies the solution of the original set of MHD equations. In particular, this is true for those types of problem in which the total pressure can be considered as a known quantity from the outset. To illustrate the method, we solve several 2D problems involving magnetic field and flow non-uniformities: reconnection in a stagnation-point flow geometry with antiparallel magnetic fields; reconnection in a Y-type magnetic field geometry with and without velocity shear across the current sheet; and reconnection in a force-free magnetic field geometry with field lines of the form xy = const. These case examples, chosen for their tractability, each incorporate some aspects of the field and flow geomtries encountered in solar-terrestrial applications, and they provide a starting point for further analytical as well as numerical studies of reconnection.     

Stability and waves of transonic laboratory and space plasmas

The properties of magnetohydrodynamic waves and instabilities of laboratory and space plasmas are determined by the overall magnetic confinement geometry and by the detailed distributions of the density, pressure, magnetic field, and background velocity of the plasma. Consequently, measurement of the spectrum of MHD waves (MHD spectroscopy) gives direct information on the internal state of the plasma, provided a theoretical model is available to solve the forward as well as the inverse spectral problems. This terminology entails a program, viz. to improve the accuracy of our knowledge of plasmas, both in the laboratory and in space. Here, helioseismology (which could be considered as one of the forms of MHD spectroscopy) may serve as a luminous example. The required study of magnetohydrodynamic waves and instabilities of both laboratory and space plasmas has been conducted for many years starting from the assumption of static equilibrium. Recently, there is a outburst of interest for plasma states where this assumption is violated. In fusion research, this interest is due to the importance of neutral beam heating and pumped divertor action for the extraction of heat and exhaust needed in future tokamak reactors. Both result in rotation of the plasma with speeds that do not permit the assumption of static equilibrium anymore. In astrophysics, observations in the full range of electromagnetic radiation has revealed the primary importance of plasma flows in such diverse situations as coronal flux tubes, stellar winds, rotating accretion disks, and jets emitted from radio galaxies. These flows have speeds which substantially influence the background stationary equilibrium state, if such a state exists at all. Consequently, it is important to study both the stationary states of magnetized plasmas with flow and the waves and instabilities they exhibit. We will present new results along these lines, extending from the discovery of gaps in the continuous spectrum and low-frequency Alfvén waves driven by rotation to the nonlinear flow patterns that occur when the background speed traverses the full range from sub-slow to super-fast. This revised version was published online in August 2006 with corrections to the Cover Date.     

Modeling coronal streamers and their eruption

Numerical solutions of the time-dependent MHD equations are used to generate ambient coronal streamer structures in a corona characteristic of that near solar minimum. The streamers are then disrupted by slow photospheric shear motion at the base of magnetic field lines within the closed field region, which is currently believed to be responsible for producing at least some CMEs. In contrast to several other simulations of this phenomena, the polytropic index is maintained at a value of 5/3 through the addition of coronal heating. Observations are used as a guide in determining the thermodynamic structure and plasma beta in the ambient corona. For a shear speed of 2.5 km/sec, the streamer configuration evolves slowly for about 65 hours before erupting outward with the formation of a CME. The bright CME leading edge travels outward at a speed of about 240 km/sec, and the sheared field lines follow at a somewhat slower speed. A closed magnetic field region is ejected as the magnetic field lines that were opened by the CME reconnect and reform the streamer.     

Alfvén wave heating

This review discusses Alfvén wave heating in non-uniform plasmas as a possible means for explaining the heating of the solar corona. It focusses on recent analytical results that enable us to understand the basic physics of Alfvén wave heating and help us with the interpretation of results of numerical simulations. First we consider the singular wave solutions that are found in linear ideal MHD at the resonant magnetic surface where the frequency of the wave equals the local Alfvén frequency. Next, we use linear resistive MHD for describing the waves in the dissipative region and explain how dissipation modifies the singular solutions found in linear ideal MHD.     

Electrodynamic coupling of the magnetosphere and ionosphere

The auroral zone ionosphere is coupled to the outer magnetosphere by means of field-aligned currents. Parallel electric fields associated with these currents are now widely accepted to be responsible for the acceleration of auroral particles. This paper will review the theoretical concepts and models describing this coupling. The dynamics of auroral zone particles will be described, beginning with the adiabatic motions of particles in the converging geomagnetic field in the presence of parallel potential drops and then considering the modifications to these adiabatic trajectories due to wave-particle interactions. The formation of parallel electric fields can be viewed both from microscopic and macroscopic viewpoints. The presence of a current carrying plasma can give rise to plasma instabilities which in a weakly turbulent situation can affect the particle motions, giving rise to an effective resistivity in the plasma. Recent satellite observations, however, indicate that the parallel electric field is organized into discrete potential jumps, known as double layers. From a macroscopic viewpoint, the response of the particles to a parallel potential drop leads to an approximately linear relationship between the current density and the potential drop.The currents flowing in the auroral circuit must close in the ionosphere. To a first approximation, the ionospheric conductivity can be considered to be constant, and in this case combining the ionospheric Ohm's Law with the linear current-voltage relation for parallel currents leads to an outer scale length, above which electric fields can map down to the ionosphere and below which parallel electric fields become important. The effects of particle precipitation make the picture more complex, leading to enhanced ionization in upward current regions and to the possibility of feedback interactions with the magnetosphere.Determining adiabatic particle orbits in steady-state electric and magnetic fields can be used to determine the self-consistent particle and field distributions on auroral field lines. However, it is difficult to pursue this approach when the fields are varying with time. Magnetohydrodynamic (MHD) models deal with these time-dependent situations by treating the particles as a fluid. This class of model, however, cannot treat kinetic effects in detail. Such effects can in some cases be modeled by effective transport coefficients inserted into the MHD equations. Intrinsically time-dependent processes such as the development of magnetic micropulsations and the response of the magnetosphere to ionospheric fluctuations can be readily treated in this framework.The response of the lower altitude auroral zone depends in part on how the system is driven. Currents are generated in the outer parts of the magnetosphere as a result of the plasma convection. The dynamics of this region is in turn affected by the coupling to the ionosphere. Since dissipation rates are very low in the outer magnetosphere, the convection may become turbulent, implying that nonlinear effects such as spectral transfer of energy to different scales become important. MHD turbulence theory, modified by the ionospheric coupling, can describe the dynamics of the boundary-layer region. Turbulent MHD fluids can give rise to the generation of field-aligned currents through the so-called -effect, which is utilized in the theory of the generation of the Earth's magnetic field. It is suggested that similar processes acting in the boundary-layer plasma may be ultimately responsible for the generation of auroral currents.     

Numerical simulation of magnetohydrodynamic shock propagation in the corona

Recent developments in the field of numerical simulation models for the study of shock wave propagation in the corona are presented. These models are based on gasdynamic (GD) and ideal (that is, dissipationless, except at shocks) magnetohydrodynamic (MHD) theories. The characteristics and physical interpretations of the results derived from these models are discussed in some detail.The most significant physical results obtained to date are provided by the two-dimensional non-planar, time-dependent, MHD numerical simulation model. In this model, the non-linear interaction among the three essential MHD waves, i.e., fast-, slow-, and Alfvén waves are demonstrated. Finally, the physical relevance of these numerical simulation models in relation to observed solar activity is presented.An invited paper presented at STIP Workshop on Shock Waves in the Solar Corona and Interplanetary Space, 15–19 June, 1980, Smolenice, Czechoslovakia.     

Progress, Challenges, and Perspectives of the 3D MHD Numerical Modeling of Oscillations in the Solar Corona

Recent high temporal and spatial resolution satellite observations of the solar corona provide ample evidence of oscillations in coronal structures. The observed waves and oscillations can be used as a diagnostic tool of the poorly known coronal parameters, such as magnetic field, density, and temperature. The emerging field of coronal seismology relies on the interpretation of the various coronal oscillations in terms of theoretically known wave modes, and the comparison of observed and theoretical wave mode properties for the determination of the coronal parameters. However, due to complexity of coronal structures the various modes are coupled, and the application of linear theory of idealized structures to coronal loops and active regions limits the usefulness of such methods. Improved coronal seismology can be achieved by the development of full 3D MHD dynamical model of relevant coronal structures and the oscillation phenomena. In addition to improved accuracy compared to linear analysis, 3D MHD models allow the diagnostic method to include nonlinearity, compressibility, and dissipation. The current progress made with 3D MHD models of waves in the corona is reviewed, and the challenges facing further development of this method are discussed in the perspective of future improvement that will be driven by new high resolution and high cadence satellite data, such as received from Hinode and STEREO, and expected from SDO.     

Investigation of MHD power generation with supersonic non-equilibrium RF discharge

Magnetohydrodynamic (MHD) power generation with supersonic non-equilibrium plasma is demonstrated. Capacitively coupled radio frequency (RF) discharge (6 MHz, maximum continual power output of 200 W) was adopted to ionize the Mach number 3.5 (650 m/s), 0.023 kg/m3 airflow. In a MHD channel of 16 mm × 10 mm × 20 mm, MHD open voltage of 10 V is realized in the magnetic field of 1.25 T, and power of 0.12 mW is extracted steadily and con-tinuously in the magnetic field of 1 T. The reasons for limited power generation are proposed as:low conductivity of RF discharge; large touch resistance between MHD electrode and plasma;strong current eddies due to flow boundary layer. In addition, the cathode voltage fall is too low to have obvious effects on MHD power generation.     

磁流体湍流中的结构和能量传输

磁流体湍流中存在多种相干结构,包括电流和涡旋结构。文章主要对近几年有关磁流体湍流相干结构及其相关的能量传输的一些工作进行了介绍。不同于中性流体湍流中的管状结构,磁流体中的相干结构多呈现为片状,强电流片附近也存在强涡片结构。磁流体湍流中的能量级串使能量跨尺度从大尺度传递至小尺度;与中性流体湍流能量传输主要集中在相近尺度上不同,磁流体湍流除了动能传输,还有磁能传输,以及动能与磁能之间的跨尺度互相转换。研究表明,磁流体湍流中的能量传输及耗散具有强间歇性,集中在较小区域范围,且与相干结构存在关联。在存在背景磁场的情况下,磁流体湍流中的相干结构在背景磁场方向上被拉长,而同样方向上的能量传输受到抑制。在高马赫数的情况下,激波的产生会使能量传输明显增强,同时带来很强的间歇性,结构函数标度律远远偏离线性标度关系,并且随着阶数的增加,标度指数会达到饱和值。     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

Reconnection Diffusion in Turbulent Fluids and Its Implications for Star Formation

Astrophysical fluids are turbulent a fact which changes the dynamics of many key processes, including magnetic reconnection. Fast reconnection of magnetic field in turbulent fluids allows the field to change its topology and connections. As a result, the traditional concept of magnetic fields being frozen into the plasma is no longer applicable. Plasma associated with a given magnetic field line at one instant is distributed along a different set of magnetic field lines at the next instant. This diffusion of plasmas and magnetic field is enabled by reconnection and therefore is termed “reconnection diffusion”. The astrophysical implications of this concept include heat transfer in plasmas, advection of heavy elements in interstellar medium, magnetic field generation etc. However, the most dramatic implications of the concept are related to the star formation process. The reason is that magnetic fields are dynamically important for most of the stages of star formation. The existing theory of star formation has been developed ignoring the possibility of reconnection diffusion. Instead, it appeals to the decoupling of mass and magnetic field arising from neutrals drifting in respect to ions entrained on magnetic field lines, i.e. through the process that is termed “ambipolar diffusion”. The predictions of ambipolar diffusion and reconnection diffusion are very different. For instance, if the ionization of media is high, ambipolar diffusion predicts that the coupling of mass and magnetic field is nearly perfect. At the same time, reconnection diffusion is independent of the ionization but depends on the scale of the turbulent eddies and on the turbulent velocities. In the paper we explain the physics of reconnection diffusion both from macroscopic and microscopic points of view, i.e. appealing to the reconnection of flux tubes and to the diffusion of magnetic field lines. We make use of the Lazarian and Vishniac (Astrophys. J. 517:700, 1999) theory of magnetic reconnection and show that this theory is applicable to the partially ionized gas. We quantify the reconnection diffusion rate both for weak and strong MHD turbulence and address the problem of reconnection diffusion acting together with ambipolar diffusion. In addition, we provide a criterion for correctly representing the magnetic diffusivity in simulations of star formation. We discuss the intimate relation between the processes of reconnection diffusion, field wandering and turbulent mixing of a magnetized media and show that the role of the plasma effects is limited to “breaking up lines” on small scales and does not affect the rate of reconnection diffusion. We address the existing observational results and demonstrate how reconnection diffusion can explain the puzzles presented by observations, in particular, the observed higher magnetization of cloud cores in comparison with the magnetization of envelopes. We also outline a possible set of observational tests of the reconnection diffusion concept and discuss how the application of the new concept changes our understanding of star formation and its numerical modeling. Finally, we outline the differences of the process of reconnection diffusion and the process of accumulation of matter along magnetic field lines that is frequently invoked to explain the results of numerical simulations.