激光驱动亥姆霍兹电容线圈靶磁重联实验的数值模拟

激光驱动亥姆霍兹电容线圈靶的磁重联实验已经提出并进行了多年.当实验中的金属板被强激光照射时产生自由电子,这些自由电子的运动在连接两金属板的两个平行线圈中产生电流,由两个平行线圈内部电流产生的磁场之间随即发生重联.该实验不同于其他直接由Biermann电池效应所产生高β(等离子体热压与磁压的比值)环境下的磁重联实验.对该类实验进行了3维磁流体动力学数值模拟,首次展示了亥姆霍兹电容器线圈靶如何驱动磁重联的过程.数值模拟结果清楚地表明,磁重联的出流等离子体在线圈周围发生与实验结果相一致的堆积现象.线圈电流产生的磁场可高达100 T,使得磁重联区域周围的等离子体β值达到10^-2.与实验室结果进行比较,数值模拟重复了实验展示的大多数特征,可有助于深入认识和理解实验结果背后的物理学原理.     

Numerical Experiments Base on the Catastrophe Model of Solar Eruptionstwo

On the basis of the catastrophe model developed by Isenberg et al., we have used the NIRVANA code to perform the magnetohydrodynamics (MHD) numerical experiments to look into the various behaviors of the coronal magnetic configuration that includes a current-carrying flux rope for modelling the prominence levitation in the corona. These behaviors include the evolution of the equilibrium height of magnetic flux rope with the background magnetic field, the corresponding interior equilibrium of magnetic flux rope, the dynamic properties of magnetic flux rope after the system loses equilibrium, as well as the impact of the reference radius on the equilibrium height of magnetic flux rope. In our calculations, an empirical model of the coronal density distribution given by Sittler & Guhathakurta is used, and the physical dissipation is included. Our experiments show that a deviation between the simulated equilibrium height of magnetic flux rope and the theoretical result of Isenberg et al. exists, but it is not apparent, and the evolutionary features of the two results are similar. If the magnetic flux rope is initially located at the stable branch of the theoretical equilibrium curve, the magnetic flux rope will quickly reach the equilibrium position after several rounds of oscillations as a result of the self-adjustment of the system; when the system is located at the critical point it will quickly lose equilibrium and evolve to the eruptive state; the impact of the variation of reference radius on the equilibrium height of magnetic flux rope is consistent with the prediction of the theory; in the eruptive state, the kinetic properties of magnetic flux rope are consistent with the results given by the Lin-Forbes model and observation, and the fast-mode shock in front of the magnetic flux rope is observed in our experiments; furthermore, because that the dissipation is included in our numerical experiments, the energy conversion from the magnetic energy to other forms of energy is very apparent in the eruptive process.     

The screw dynamo in a thick torus

The problem of magnetic field generation under screw motion in a toroidal channel is studied numerically. The screw dynamo in the cylinder with periodical boundary conditions was found to be a suitable approximation for generation of the magnetic field by a screw flow in a thin torus. For the thick torus, a principally new solution of the screw dynamo problem was obtained. In this case the growing global magnetic field mode has the scale of a maximal geometrical size of the torus and does not vanish on the axis of the torus (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modeling spatio‐temporal nonlocality in mean‐field dynamos

When scale separation in space or time is poor, the mean‐field α effect and turbulent diffusivity have to be replaced by integral kernels by which the dependence of the mean electromotive force on the mean magnetic field becomes nonlocal. Earlier work in computing these kernels using the test‐field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel for isotropic stationary turbulence is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean‐field equations are solved for oscillatory α –shear dynamos as well as α² dynamos with α linearly depending on position, which makes this dynamo oscillatory, too. In both cases, the critical values of the dynamo number is lowered due to spatio‐temporal nonlocality.When scale separation in space or time is poor, the mean‐field α effect and turbulent diffusivity have to be replaced by integral kernels by which the dependence of the mean electromotive force on the mean magnetic field becomes nonlocal. Earlier work in computing these kernels using the test‐field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel for isotropic stationary turbulence is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean‐field equations are solved for oscillatory α –shear dynamos as well as α² dynamos     

实验室模拟磁力线重联的初步结果

在中国科学技术大学的线性磁化等离子体装置上,通过对两个平行电流板施加同向电流,实现重联磁场位型的构造,进而开展实验室等离子体中磁力线重联过程的研究.利用发射探针测量了重联过程中的平行(轴向)电场,实验验证了重联电流与通行粒子的依赖关系.利用磁探针测量了磁场通量的演化,未发现通量堆积现象,与数值预言相符.     

Reduction of boundary effects in the spiralMRI experiment PROMISE

Magnetorotational instability (MRI) is one of the most important and most common instabilities in astrophysics. It is widely accepted that it serves as a source of turbulent viscosity in accretion disks – the most energy efficient objects in the Universe. However it is very difficult to bring this process down on earth and model it in a laboratory experiment. Several different approaches have been proposed, one of the most recent is PROMISE (Potsdam‐ROssendorf Magnetorotational InStability Experiment). It consists of a flow of a liquid metal between two rotating cylinders under applied current‐free spiral magnetic field. The cylinders must be covered with plates which introduce additional end‐effects which alter the flow and make it more difficult to clearly distinguish between MRI stable and unstable state. In this paper we propose simple and inexpensive improvement to the PROMISE experiment which would reduce those undesirable effects. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The SHASTA Code Modified by Self-adaptive Mesh and Numerical Experiment of Magnetic Reconnections

SHASTA (Sharp and Smooth Transport Algorithm) is a code with single mesh to solve the 2-dimensional magnetohydrodynamic (MHD) equations. When SHASTA is used to the numerical simulation of magnetic reconnection problem, it is modified to be the code which adopts the method of the selfadaptive mesh. The modified code can carry out refined calculations in diffusion regions. In the process of the self-adaptive calculations with SHASTA, a “plugand-play” strategy is adopted and the original algorithm to solve 2-dimensional MHD partial differential equations is treated as an independent cell. In addition, the hierarchical data structure is used in this modification and parameters in each refined level are described by a 2-dimensional variable array. The regions where the distributions of magnetic field and pressure exhibit steep variations are marked as the refined regions. Then, the distributions of physical quantities and the boundary conditions in the grid points of refined levels are deduced via interpolation method. Finally, the refined calculated results of refined regions are assigned to the previous level of mesh and the existing results are updated. The numerical experiment of magnetic reconnections which adopts refined calculations indicates that compared with the code with single mesh, the resolution of details is improved and the corresponding increment of computing time is related to the selection of parameters in the simulation. The calculation accuracy and effect on instability, which are caused by a part of the self-adaptive code, depend on the boundary settings, push strategy over each single step as well as the interpolation algorithm.     

耀斑大尺度电流片的湍流耗散和能谱分析

磁重联在宇宙的许多动力学现象中都是非常核心的过程.磁流体动力学(MHD)数值模拟是研究磁重联过程以及相应物理图像的一种很有效的手段.通过不同的参数组合,来研究MHD数值模拟中磁雷诺数和空间分辨率对磁重联率、数值耗散和能谱分布的影响.对得到的数据进行分析后,发现磁雷诺数对磁重联率和能谱分布有一定的影响.磁雷诺数越大,磁重联过程进入非线性阶段所需的特征时间越短,磁重联率就越早发生跃升.磁雷诺数R_m对耗散开始发挥作用的Kolmogorov微观尺度l_ko有明显影响:R_m越大,l_ko就越小.研究了磁重联过程中包括数值耗散在内的额外耗散对重联过程的影响.结果表明,撕裂模不稳定性开始之前的额外耗散以纯数值耗散为主,撕裂模不稳定性出现之后,额外耗散出现同步跃升,说明不稳定性导致的湍流明显增强了耗散的效果,相当于在局部湍流区引入了超电阻.能谱分析进一步表明,大尺度电流片的l_ko完全可能出现在宏观的MHD尺度上.     

3-D non-linear evolution of a magnetic flux tube in a spherical shell: The isentropic case

We present recent 3-D MHD numerical simulations of the non-linear dynamical evolution of magnetic flux tubes in an adiabatically stratified convection zone in spherical geometry, using the anelastic spherical harmonic (ASH) code.We seek to understand the mechanism of emergence of strong toroidal fields from the base of the solar convection zone to the solar surface as active regions. We confirm the results obtained in cartesian geometry that flux tubes that are not twisted split into two counter vortices before reaching the top of the convection zone. Moreover, we find that twisted tubes undergo the poleward-slip instability due to an unbalanced magnetic curvature force which gives the tube a poleward motion both in the non-rotating and in the rotating case. This poleward drift is found to be more pronounced on tubes originally located at high latitudes. Finally, rotation is found to decrease the rise velocity of the flux tubes through the convection zone, especially when the tube is introduced at low latitudes. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonlinear simulations of magnetic Taylor‐Couette flow with currentfree helical magnetic fields

Themagnetorotational instability (MRI) in cylindrical Taylor‐Couette flow with external helical magnetic field is simulated for infinite and finite aspect ratios. We solve the MHD equations in their small Prandtl number limit and confirm with timedependent nonlinear simulations that the additional toroidal component of the magnetic field reduces the critical Reynolds number from O (10⁶) (axial field only) to O (10³) for liquid metals with their small magnetic Prandtl number. Computing the saturated state we obtain velocity amplitudes which help designing proper experimental setups. Experiments with liquid gallium require axial field ∼50 Gauss and axial current ∼4 kA for the toroidal field. It is sufficient that the vertical velocity u_z of the flow can be measured with a precision of 0.1 mm/s.We also show that the endplates enclosing the cylinders do not destroy the traveling wave instability which can be observed as presented in earlier studies. For TC containers without and with endplates the angular momentum transport of the MRI instability is shown as to be outwards. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Transition-region explosive events produced by plasmoid instability

Magnetic reconnection is thought to be a key process in most solar eruptions. Thanks to highresolution observations and simulations, the studied scale of the reconnection process has become smaller and smaller. Spectroscopic observations show that the reconnection site can be very small, which always exhibits a bright core and two extended wings with fast speeds, i.e., transition-region explosive events.In this paper, using the PLUTO code, we perform a 2-D magnetohydrodynamic simulation to investigate small-scale reconnection in double current sheets. Based on our simulation results, such as the line-of-sight velocity, number density and plasma temperature, we can synthesize the line profile of SiIV 1402.77? which is a well known emission line used to study transition-region explosive events on the Sun. The synthetic line profile of Si IV 1402.77? is complex with a bright core and two broad wings which can extend to nearly 200 km s^-1. Our simulation results suggest that the transition-region explosive events on the Sun are produced by plasmoid instability during small-scale magnetic reconnection.     

Dynamo saturation in direct simulations of the multi‐phase turbulent interstellar medium

The ordered magnetic field observed via polarised synchrotron emission in nearby disc galaxies can be explained by a mean‐field dynamo operating in the diffuse interstellar medium (ISM). Additionally, vertical‐flux initial conditions are potentially able to influence this dynamo via the occurrence of the magnetorotational instability (MRI). We aim to study the influence of various initial field configurations on the saturated state of the mean‐field dynamo. This is motivated by the observation that different saturation behaviour was previously obtained for different supernova rates. We perform direct numerical simulations (DNS) of three‐dimensional local boxes of the vertically stratified, turbulent interstellar medium, employing shearing‐periodic boundary conditions horizontally. Unlike in our previous work, we also impose a vertical seed magnetic field. We run the simulations until the growth of the magnetic energy becomes negligible. We furthermore perform simulations of equivalent 1D dynamo models, with an algebraic quenching mechanism for the dynamo coefficients. We compare the saturation of the magnetic field in the DNS with the algebraic quenching of a mean‐field dynamo. The final magnetic field strength found in the direct simulation is in excellent agreement with a quenched α) dynamo. For supernova rates representative of the Milky Way, field losses via a Galactic wind are likely responsible for saturation. We conclude that the relative strength of the turbulent and regular magnetic fields in spiral galaxies may depend on the galaxy's star formation rate. We propose that a mean field approach with algebraic quenching may serve as a simple sub‐grid scale model for galaxy evolution simulations including a prescribed feedback from magnetic fields. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Catastrophe of coronal magnetic rope in partly open multipolar magnetic field

Catastrophe of coronal magnetic rope embedded in a partly open multipolar background magnetic field is studied by using a 2-dimensional, 3-component ideal MHD model in spherical coordinates. The background field is composed of three closed bipolar fields of a coronal streamer and an open field with an equatorial current sheet. The magnetic rope lies below the central bipolar field, and it is characterized by its annular and axial magnetic fluxes. For a given annual flux, there is a critical value of the axial flux, and for a given axial flux, there is a critical value of annual flux such that, below the critical value, the magnetic rope is attached to the solar surface and the system stays in equilibrium, but when the critical value is exceeded, the magnetic rope breaks free and erupts upward. This implies that catastrophe can occur in a coronal magnetic rope embedded in a partly open multipolar background magnetic field. Our computation gives a threshold value of magnetic energy that is about 15% greater than the energy of the partly open magnetic field (the central bipolar field open and the fields on either side closed). The excess energy may serve as source for solar explosions such as coronal mass ejections.     

无碰撞重联区中的混沌感应电阻

无碰撞磁场重联作为一种将磁能有效转化为等离子体动能和热能的机制,已经被广泛应用于解释太阳耀斑、地球磁暴等各类等离子体的爆发活动.然而,在无碰撞重联区中反常电阻的微观物理机制仍然是尚未解决的基本问题.在众多反常电阻的形成机制中,基于磁零点附近粒子轨道混沌性产生的混沌感应电阻,虽然不是最普遍流行的形成机制,但它的微观物理图像却是最为清晰的.回顾了无碰撞重联区中混沌感应电阻的早期研究和基本理论模型,介绍了关于混沌感应电阻研究的新进展并阐述了混沌感应电阻未来的研究方向.     

Self‐similar evolutionary solutions for an accreting magneto‐fluid around a compact object with finite electrical conductivity

In this paper, we investigate the time evolution of an accreting magneto‐fluid with finite conductivity. For the case of a thin disk, the fluid equations along with Maxwell's equations are derived in a simplified, one‐dimensional model that neglects the latitudinal dependence of the flow. The finite electrical conductivity of the plasma is taken into account by Ohm's law; however, the shear viscous stress is neglected, as well as the self‐gravity of the disk. In order to solve the integrated equations that govern the dynamical behaviour of the magneto‐fluid, we have used a self‐similar solution. We introduce two dimensionless variables, S₀ and ε_ϱ, which represent the size of the electrical conductivity and the density behaviour with time, respectively. The effect of each of these on the structure of the disk is studied. While the pressure is obtained simply by solving an ordinary differential equation, the density, the magnetic field, the radial velocity, and the rotational velocity are presented analytically. The solutions show that the S0 and ε_ϱ parameters affect the radial thickness of the disk. Also, radial velocity and gas pressure are more sensitive to the electrical conductivity in the inner regions of disk. Moreover, the parameter ε_ϱ has a more significant effect on the physical quantities for small radii. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulence Diffusion and Energy Spectrum Analysis of Large-scale Current Sheets in Flares

Magnetic reconnection (MR) is one of the most important physical processes for many dynamical phenomena in the universe. Magnetohydrodynamical (MHD) simulation is an effective way to study the MR process and the physical pictures related to the MR. With different parameter setups, we investigate the influences of the Magnetic Reynolds number and spatial resolution on the reconnection rate, numerical dissipation, and energy spectrum distribution in the MHD simulation. We have found that the magnetic Reynolds number $R_m$ has definite impact on the reconnection rate and energy spectrum distribution. The characteristic time for entering into the non-linear phase will be earlier as the Reynolds number increases. When it comes to the tearing phase, the reconnection rate will increase rapidly. On the other hand, the magnetic Reynolds number affects significantly the Kolmogorov microscopic scale $l_{\mathit{ko}}$, which becomes smaller as $R_m$ increases. An extra dissipation is defined as the combined effect of the numerical diffusion and turbulence dissipation. It is shown that the extra dissipation is dominated by the numerical diffusion before the tearing mode instability takes place. After the instability develops, the extra dissipation rises vastly, which indicates that turbulence caused by the instability can enhance the diffusion obviously. Furthermore, the energy spectrum analysis indicates that $l_{\mathit{ko}}$ of the large-scale current sheet may appear at a macroscopic MHD scale very possibly.     

The problem of phase mixed shear Alfvén waves in the solar corona revisited

The problem of phase mixing of shear Alfvén waves is revisited taking into account dissipative phenomena specific for the solar corona. In regions of space plasmas where the dynamics is controlled by the magnetic field, transport coefficients become anisotropic with transport mechanism having different behavior and magnitude depending on the orientation with respect to the ambient magnetic field. Taking into account realistic values for dissipative coefficients we obtain that the previous results derived in context of torsional Alfvén wave phase mixing are actually heavily underestimated so phase mixing cannot be used to explain the damping of torsional Alfvén waves and heating of open coronal structures. The presented results indicate that in order for phase mixing to still be a viable mechanism to explain heating or wave damping unrealistic assumptions have to be made. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Suppression of the large‐scale Lorentz force by turbulence

The components of the total stress tensor (Reynolds stress plus Maxwell stress) are computed within the quasilinear approximation for a driven turbulence influenced by a large‐scale magnetic background field. The conducting fluid has an arbitrary magnetic Prandtl number and the turbulence without the background field is assumed as homogeneous and isotropic with a free Strouhal number St. The total large‐scale magnetic tension is always reduced by the turbulence with the possibility of a ‘catastrophic quenching’ for large magnetic Reynolds number Rm so that even its sign is reversed. The total magnetic pressure is enhanced by turbulence in the high‐conductivity limit but it is reduced in the low‐conductivity limit. Also in this case the sign of the total pressure may reverse but only for special turbulences with sufficiently large St > 1. The turbulence‐induced terms of the stress tensor are suppressed by strong magnetic fields. For the tension term this quenching grows with the square of the Hartmann number of the magnetic field. For microscopic (i.e. small) diffusivity values the magnetic tension term becomes thus highly quenched even for field amplitudes much smaller than their equipartition value. In the opposite case of large‐eddy simulations the magnetic quenching is only mild but then also the turbulence‐induced Maxwell tensor components for weak fields remain rather small (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Simulation Study of the Magnetic Reconnections Between Closed Loops and Open Magnetic Fields Driven by Horizontal Flows on Solar Surface

Numerous observational events in the solar atmosphere (e.g., solar ?ares and jets) are attributed to the energy conversion due to magnetic reconnec- tions. Magnetic reconnections are also involved in a new scenario of solar wind origin to play a crucial role in opening the closed magnetic loop and releasing its mass into the open magnetic funnel. In this scenario, the closed magnetic loop moves towards the supergranular boundary by the supergranular convection, and collides with the open magnetic funnel there to trigger the magnetic reconnec- tion between each other. This work aims at studying the occurrence and effect of magnetic reconnection in this scenario in detail. The magnetohydrodynamic (MHD) numerical simulation is an important approach to investigate the mag- netic reconnection process in the solar atmosphere. A two-dimensional MHD numerical model has been established, and in combination with the strati?ed temperature and density distributions in the solar atmosphere, the numerical simulation on the process of magnetic reconnection of the closed magnetic loops driven by the horizontal ?ows with the open magnetic ?elds has been performed on the scale of supergranulation. Based on a quantitative analysis of the simula- tion result, it is suggested that the process of magnetic reconnection can really realize the mass release of closed magnetic loops, and further supply to the new open magnetic structures to produce upward mass ?ows. Our results provide a basis for the further modeling of solar wind origin.