Spectral Direct Numerical Simulations of low Rm MHD channel flows based on the least dissipative modes

TL;DR

This paper introduces a spectral method using eigenfunctions of the dissipation operator to efficiently simulate low Rm MHD channel flows with boundary layers, reducing computational costs associated with high magnetic fields.

Contribution

The paper develops a novel spectral method based on dissipation eigenfunctions, enabling efficient simulation of MHD flows with boundary layers without increased computational cost at high magnetic fields.

Findings

Accurately simulates wall-bounded MHD turbulence with fewer modes.

Validates the method against finite volume simulations.

Demonstrates potential for high magnetic field MHD turbulence simulations.

Abstract

We put forward a new type of spectral method for the direct numerical simulation of flows where anisotropy or very fine boundary layers are present. The mean idea is to take advantage of the fact that such structures are dissipative and that their presence should reduce the number of degrees of freedom of the flow, when paradoxically, their fine resolution incurs extra computational cost in most current methods. The principle of this method is to use a functional basis with elements that already include these fine structure so as to avoid these extra costs. This leads us to develop an algorithm to implement a spectral method for arbitrary functional bases, and in particular, non-orthogonal ones. We construct a basic implementation of this algorithm to simulate Magnetohydrodynamic (MHD) channel flows with an externally imposed, transverse magnetic field, where very thin boundary layers are known to develop along the channel walls. In this case, the sought functional basis can be built out of the eigenfunctions of the dissipation operator, which incorporate these boundary layers, and it turns out to be non-orthogonal. We validate this new scheme against numerical simulations of freely decaying MHD turbulence based on a Finite Volume code and it is found to provide accurate results. Its ability to fully resolve wall-bounded turbulence with a number of modes close to that required by the dynamics is demonstrated on a simple example. This opens the way to full blown simulations of MHD turbulence under very high magnetic fields, which were until now too computationally expensive, as the computational cost of the proposed method, in contrast to traditional methods, does not depend on the intensity of the magnetic field.