No Arabic abstract
In this study, the evolution of a highly unstable m = 1 resistive tearing mode, leading to plasmoid formation in a Harris sheet is studied in the framework of full MHD model using the NIMROD simulation. Following the initial nonlinear growth of the primary m = 1 island, the X-point develops into a secondary elongated current sheet that eventually breaks into plasmoids. Two distinctive viscous regimes are found for the plasmoid formation and saturation. In the low viscosity regime (i.e. P r . 1), the plasmoid width increases sharply with viscosity, whereas in the viscosity dominant regime (i.e. P r & 1 ), the plasmoid size gradually decreases with viscosity. Such a finding quantifies the role of viscosity in modulating the plasmoid formation process through its effects on the plasma flow and the reconnection itself.
Axisymmetric current-carrying plasmoids are formed in the presence of nonaxisymmetric fluctuations during nonlinear three-dimensional resistive MHD simulations in a global toroidal geometry. We utilize the helicity injection technique to form an initial poloidal flux in the presence of a toroidal guide field. As helicity is injected, two types of current sheets are formed from 1) the oppositely directed field lines in the injector region (primary reconnecting current sheet), and 2) the poloidal flux compression near the plasma edge (edge current sheet). We first find that nonaxisymmetic fluctuations arising from the current-sheet instability isolated near the plasma edge have tearing parity but can nevertheless grow fast (on the poloidal Alfven time scale). These modes saturate by breaking up the current sheet. Second, for the first time a dynamo poloidal flux amplification is observed at the reconnetion site (in the region of the oppositely directed magnetic field). This fluctuation-induced flux amplification increases the local Lundquist number, which then triggers a plasmoid instability and breaks the primary current sheet at the reconnection site. The plasmoids formation driven by large-scale flux amplification, i.e. a large-scale dynamo, observed here has strong implications for astrophysical reconnection as well as fast reconnection events in laboratory plasmas.
The effects of line-tying on resistive tearing instability in slab geometry is studied within the framework of reduced magnetohydrodynamics (RMHD).citep{KadomtsevP1974,Strauss1976} It is found that line-tying has a stabilizing effect. The tearing mode is stabilized when the system length $L$ is shorter than a critical length $L_{c}$, which is independent of the resistivity $eta$. When $L$ is not too much longer than $L_{c}$, the growthrate $gamma$ is proportional to $eta$ . When $L$ is sufficiently long, the tearing mode scaling $gammasimeta^{3/5}$ is recovered. The transition from $gammasimeta$ to $gammasimeta^{3/5}$ occurs at a transition length $L_{t}simeta^{-2/5}$.
New non-linear, spatially periodic, long wavelength electrostatic modes of an electron fluid oscillating against a motionless ion fluid (Langmuir waves) are given, with viscous and resistive effects included. The cold plasma approximation is adopted, which requires the wavelength to be sufficiently large. The pertinent requirement valid for large amplitude waves is determined. The general non-linear solution of the continuity and momentum transfer equations for the electron fluid along with Poissons equation is obtained in simple parametric form. It is shown that in all typical hydrogen plasmas, the influence of plasma resistivity on the modes in question is negligible. Within the limitations of the solution found, the non-linear time evolution of any (periodic) initial electron number density profile n_e(x, t=0) can be determined (examples). For the modes in question, an idealized model of a strictly cold and collisionless plasma is shown to be applicable to any real plasma, provided that the wavelength lambda >> lambda_{min}(n_0,T_e), where n_0 = const and T_e are the equilibrium values of the electron number density and electron temperature. Within this idealized model, the minimum of the initial electron density n_e(x_{min}, t=0) must be larger than half its equilibrium value, n_0/2. Otherwise, the corresponding maximum n_e(x_{max},t=tau_p/2), obtained after half a period of the plasma oscillation blows up. Relaxation of this restriction on n_e(x, t=0) as one decreases lambda, due to the increase of the electron viscosity effects, is examined in detail. Strong plasma viscosity is shown to change considerably the density profile during the time evolution, e.g., by splitting the largest maximum in two.
A method for solving model nonlinear equations describing plasma oscillations in the presence of viscosity and resistivity is given. By first going to the Lagrangian variables and then transforming the space variable conveniently, the solution in parametric form is obtained. It involves simple elementary functions. Our solution includes all known exact solutions for an ideal cold plasma and a large class of new ones for a more realistic plasma. A new nonlinear effect is found of splitting of the largest density maximum, with a saddle point between the peaks so obtained. The method may sometimes be useful where Inverse Scattering fails.
Magnetic reconnection may be the fundamental process allowing energy stored in magnetic fields to be released abruptly, solar flares and coronal mass ejection (CME) being archetypal natural plasma examples. Magnetic reconnection is much too slow a process to be efficient on the large scales, but accelerates once small enough scales are formed in the system. For this reason, the fractal reconnection scenario was introduced (Shibata and Tanuma 2001) to explain explosive events in the solar atmosphere: it was based on the recursive triggering and collapse via tearing instability of a current sheet originally thinned during the rise of a filament in the solar corona. Here we compare the different fractal reconnection scenarios that have been proposed, and derive generalized scaling relations for the recursive triggering of fast, `ideal - i.e. Lundquist number independent - tearing in collapsing current sheet configurations with arbitrary current profile shapes. An important result is that the Sweet-Parker scaling with Lundquist number, if interpreted as the aspect ratio of the singular layer in an ideally unstable sheet, is universal and does not depend on the details of the current profile in the sheet. Such a scaling however must not be interpreted in terms of stationary reconnection, rather it defines a step in the accelerating sequence of events of the ideal tearing mediated fractal cascade. We calculate scalings for the expected number of plasmoids for such generic profiles and realistic Lundquist numbers.