Resistive tearing: numerical exploration of nonmodal effects

Elias Pratschke
2024-10-29
Abstract:The fluid dynamics community has found success in explaining both the onset and coherent structure formation in wall-bounded turbulence through examining transient growth and pseudoresonance. Whether similar effects are important in plasmas well-described by magnetohydrodynamics is an open question. In nuclear fusion experiments, the onset of turbulence often enhances undesirable heat transport between the core and the edge, and a vast number of studies have attempted to understand these phenomena via linear and nonlinear numerical codes. If there are plasma experiments where instabilities and turbulence onset are dominated by nonmodal effects, the tools of non-modal stability and resolvent analysis could prove to be essential to practitioners. Towards that goal, in this work we provide a detailed derivation of non-modal stability tools and the resolvent operator for incompressible, resistive magnetohydrodynamics (MHD). Concretely, optimal initial conditions that maximize transient growth are computed, along with bounds on linear transient growth, and are numerically verified with a nonlinear solver to show that these linear effects can be reproduced in nonlinear simulations. The potential of such initial conditions to initiate the tearing mode in a spectrally stable system is explored. We found the approach common in nonmodal hydrodynamic stability of using pseudospectral methods to be infeasible for many parameter combinations due to numerical ill-conditioning. Further, transient growth in the Harris current sheet proved to be highly norm-dependent, with the observed differences being of an order of magnitude. Nonlinear simulations showed that transient growth occurs to some extent in all norms, both in the infinitesimal as well as finite limit. No case was observed of transient growth being strong enough to trigger the full plasmoid instability.
Plasma Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is whether non - modal effects in magnetohydrodynamics (MHD), especially in plasma physics, can explain the rapid onset of certain instabilities and the formation of turbulent structures. Specifically, the paper explores the following points: 1. **Non - modal stability analysis**: Traditional linear stability analysis usually relies on eigenvalue analysis, and this method has limitations in predicting the rapid onset of certain phenomena in plasma (such as tearing modes and magnetic reconnection). By introducing non - modal stability analysis tools, the paper attempts to go beyond the traditional eigenvalue method to better understand these phenomena. 2. **Transient growth**: The paper derives in detail the non - modal stability tools and analytical operators for incompressible, resistive MHD, and calculates the optimal initial conditions that can maximize transient growth. In addition, it also verifies through a nonlinear spectral solver whether these linear effects can be reproduced in nonlinear simulations. 3. **Triggering of tearing modes**: It studies whether these optimal initial conditions can trigger tearing modes in spectrally stable systems. Tearing modes are an important plasma instability and are closely related to the magnetic reconnection process. 4. **Challenges of numerical methods**: The paper points out that under many parameter combinations, it is not feasible to use the pseudospectral method for non - modal hydrodynamic stability analysis because of numerical ill - conditioning problems. In addition, transient growth in Harris current sheets is highly norm - dependent, and the differences under different norms can reach an order of magnitude. 5. **Results of nonlinear simulations**: Nonlinear simulations show that transient growth occurs to some extent in all norms, whether in the case of infinitesimal perturbations or finite perturbations. However, no case has been observed where transient growth is strong enough to trigger a complete plasma instability. In summary, this paper aims to explore the onset mechanisms of certain instabilities and turbulent phenomena in plasma through non - modal stability analysis and the study of transient growth, especially whether these non - modal effects can explain the experimentally observed phenomena. This not only helps to understand the basic principles of plasma physics but also provides theoretical support for the design and control in application fields such as nuclear fusion.