Straining of magnetic fields by large-scale shear flow, which is generally assumed to lead to intensification and generation of small scales, is reexamined in light of the persistent observation of large-scale magnetic fields in astrophysics. It is shown that, in magnetohydrodynamic turbulence, unstable shear flows have the unexpected effect of sequestering magnetic energy at large scales due to counteracting straining motion of nonlinearly excited large-scale stable eigenmodes. This effect is quantified via dissipation rates, energy transfer rates, and visualizations of magnetic field evolution by artificially removing the stable modes. These analyses show that predictions based upon physics of the linear instability alone miss substantial dynamics, including those of magnetic fluctuations.

Turbulence is a fundamentally multiscale process in which nonlinear dynamical motions carry energy across scales.^{1,2} One of the most robust and extensively studied mechanisms for cross-scale energy transfer is straining by shear flow,^{2,3} long recognized for its role in cascades to small scale.^{4} This process is especially active in magnetohydrodynamic (MHD) turbulence.^{5} The stretching, squeezing, and folding of magnetic field lines by shear flow are readily discernible in visualizations and lead to, for example, the enhancement of small-scale resistive dissipation and the intensification of magnetic energy at small scales beyond the viscous cutoff in turbulence with the magnetic Prandtl number greater than unity.^{5,6} The ubiquity of turbulence and turbulent straining is at odds with the well-established observations of magnetic fields at large scales in the universe (stellar, galactic, and beyond)^{7,8} and part of a conundrum as to how magnetic fields exist at such large scales.^{9–11}

In this Letter, we report the first observations that magnetic energy is, to a considerable degree, counterintuitively sequestered at large scales in straining by a paradigmatic turbulent shear flow. The result applies to a Kelvin–Helmholtz (KH)-unstable flow maintained by an external force with a uniform magnetic field that is initially flow-aligned and, therefore, optimally configured to promote small-scale generation.

The process that counteracts the small scale generation of magnetic energy, which is quantified in this Letter, can be traced to the nonlinear excitation of large-scale stable linear eigenmodes.^{12} The dispersion relation that is associated with the shear-flow instability present in this study has a stable (damped) root, which acts to return energy and momentum to the mean profile and which is excited to a high level by the nonlinearity.^{13,14} The nonlinear excitation of stable modes and its effect on turbulence levels and transport have been studied previously, particularly in the context of fusion-relevant gyroradius-scale turbulence.^{12,15–23} Recent studies have shown that stable modes are also excited in macroscopic shear-flow driven turbulence and transiently affect momentum transport.^{13,14} However, critical features of counteracting straining motion of the stable modes, essential for large-scale sequestering of magnetic energy, have only come to light in the study described here.

Moreover, stable-mode effects have been linked to small-scale dissipation,^{20} suggesting that they affect turbulence only if there is an increase in entropy. However, we show here that the reversible process by which stable modes return energy to the mean flow, in opposition to extraction by the instability, is quite effective at blunting the energy cascade to small scales. Stable modes have been connected to the shear-layer contraction (i.e., build-up of mean flow)^{14} but only as a transient process. Here, we demonstrate that for a driven flow profile, the return of fluctuation energy to the mean profile occurs continuously at a rate that nearly matches the rate at which unstable modes attempt to flatten the flow profile.

The principal result of this Letter is that, in driven shear-flow MHD turbulence, the nonlinearly excited large-scale stable eigenmodes efficiently strain the magnetic fields in a way that counteracts the effect of unstable modes and, thus, greatly weakens the magnetic energy cascade to small scales.

We study an incompressible, two-dimensional (2D) system with the initial fluid velocity given as $v(x,z,t=0)=Uref(z)x\u0302=U0\u2009tanh(z/a)x\u0302$ and an initially uniform flow-aligned magnetic field as $B(x,z,t=0)=B0x\u0302$. The parameters $a,U0,and\u2009B0$ represent the half-width of the flow-shear, the maximum fluid velocity, and the magnetic field strength, respectively, which are used to non-dimensionalize all the variables henceforth. Consequently, time and energy (per unit mass) have units of $a/U0$ and $U02$. We describe the flow and the magnetic field by a stream function $\varphi $ and a flux function *ψ*, where $v=y\u0302\xd7\u2207\varphi $ and $B=y\u0302\xd7\u2207\psi $. We study their evolution using the momentum and induction equations of the standard MHD (see Ref. 24 and the supplementary material)

where $\u22072\varphi $ and $\u22072\psi $ are the vorticity and the current density, respectively; the Poisson bracket is ${P,Q}\u2261\u2202xP\xb7\u2202zQ\u2212\u2202xQ\xb7\u2202zP$; the Alfvénic Mach number is $MA\u221dU0/B0$; the fluid and magnetic Reynolds numbers are defined as $Re=U0a/\nu $ and $Rm=U0a/\eta $. In this study, we take $Re=500$. Except for simulations where $Rm=50$ and $Rm=5\u2009000$ (which will be stated explicitly), all others have $Rm=500$. In Eq. (1a), *f* represents a body force that continuously regenerates the mean-flow shear. The expression $kx=0$ in it indicates that the force acts only on the (*x*-averaged) mean flow, where *k _{x}* is the Fourier wavenumber. The forcing, thus, prevents gradual flattening of the mean flow as the instability extracts energy from it. A similar process occurs in astrophysical flows, where forces like gravitation replenish continually the shear profile, e.g., shear layers in accretion disks, stars, and planetary atmospheres. We represent this force by a Krook-like operator,

^{25–27}sometimes referred to as a profile relaxation term,

^{28}

where $DKrook$ represents the rate at which the mean flow is forced toward the reference flow profile. We choose $DKrook=2$. Detailed analyses of a complete scan of different parameters will be reported elsewhere. The force *F*_{0} balances the viscous diffusion of the mean shear layer, $Re\u22121\u22072Uref(z)+F0=0$, in Eq. (1a), to ensure an initial equilibrium state.

We add small-amplitude perturbations to the flow and the magnetic field (see Ref. 14 and the supplementary material) and perform time-integration of the above set of equations using the pseudospectral code Dedalus^{29} for long times (*t *>* *1000; the $e$-folding growth time of the instability is $\u22486$). Apart from the simulation of $Pm=10$ that uses a spectral resolution of 4096 Fourier (Chebyshev) modes along the *x*-(*z*-) axis for the box size of $(Lx,Lz)=(6\pi ,8\pi )$, all other simulations are performed at 2048 × 2048 spectral resolution for the box size of $(Lx,Lz)=(10\pi ,20\pi )$, while convergence checks utilize resolutions as high as 8192 × 8192 with no substantial impact on dissipation rates and spectral energy densities. (All simulations additionally use a 3/2 dealiasing rule.^{29})

We also perform a corresponding eigenvalue calculation of the non-dissipative equations, derived from Eqs. (1a) and (1b), by linearizing about the initial flow and magnetic field profiles.^{14} *A priori*, one does not know whether such an eigenmode basis is useful to understand turbulent features. Nevertheless, since the eigenmodes of this linear operator form a complete (albeit non-orthogonal) basis,^{30,31} we can expand the time-evolving state space of nonlinear simulations in this basis after determining amplitude $\beta j(kx,t)$ of each eigenmode *j* at each *k _{x}* as a function of time, $\chi \u0303turbulent(x,z,t)=\u2211kx:kx\u22600eikxx[\u2211j\beta j(kx,t)\chi j(kx,z)],$ where $\chi \u0303turbulent=(\varphi \u0303turbulent,\psi \u0303turbulent)$ captures the turbulent state space and

*χ*is the $jth$ eigenmode.

_{j}^{14,32–34}This procedure is illustrated in Fig. 1(a), where the eigenmodes of the flow fluctuations are shown. (The magnetic fluctuations have eigenmodes

^{14}and parity

^{37,38}related to those of the flow fluctuations.)

At each wavenumber *k _{x}*, linearly unstable to perturbations (i.e., $0<kx<1$ for the present shear-flow instability), there is an unstable eigenmode $\varphi 1$ and its conjugate stable eigenmode $\varphi 2$ with eigenvalues complex conjugate to each other.

^{13,39}All the remaining eigenmodes have purely real eigenfrequencies,

^{32}belong to the continuum

^{40}and have narrow structures along the

*z*-axis.

Tracking the time evolution of each eigenmode amplitude $\beta j(kx,t)$ in the nonlinear simulations, the initially exponentially decaying mode has been found to be nonlinearly excited to almost the same level as the unstable mode.^{14} The driving mechanism for this excitation is, at early times, the nonlinear interaction between the unstable modes, see, e.g., Refs. 14, 22, and 23. In the fully nonlinear phase, however, all modes that have large energy interact dominantly.

For a given KH-unstable wavenumber, it is found, in Fig. 1(b), that a combination of two eigenmodes from the above basis reconstructs the large-scale turbulent flow at that wavenumber. Figure 1(b) shows a reconstruction of the turbulent flow features using unstable eigenmodes at each unstable wavenumber (left column) and with their conjugate stable eigenmodes added (right column). Simply adding one additional mode drastically decreases the difference between the true stream function and its reduced representation compared to what is achieved with a reduced representation based on the unstable modes alone (such as those predicted by quasilinear models).^{41}

Note that the stable and unstable eigenmodes of the KH linear operator, employed above, are related to one another by a space-time-reversal symmetry operation.^{31,42} Thus, the unstable and stable eigenmodes transfer energy in opposite directions—the former from the mean flow to the fluctuations, while the latter in the reverse direction. The energy that is available to cascade to small scales, thus, depends on this competition. The magnitude of the energy transfer rates between the background shear-flow $Uref$ and the fluctuation scale $kx\u22600$ by the unstable and stable eigenmodes can be derived from the MHD equations (see the supplementary material): $Q1(kx)=\Gamma (kx)|\beta 1(kx)|2$ and $Q2(kx)=\Gamma (kx)|\beta 2(kx)|2$, where $\Gamma (kx)$ is the linear growth rate of the unstable mode at that wavenumber *k _{x}*.

The time traces of *Q*_{1} and *Q*_{2} for different $Pm$ (resistivity) over a challenging two orders of magnitude are shown in Fig. 2. Summation over the KH-unstable wavenumbers, $0<kx<1$, is used in determining *Q*_{1} and *Q*_{2}. Within each subplot, it is observed that *Q*_{1} and *Q*_{2} are nearly in equipartition in the nonlinear phase.

A stable-mode-excitation process may also be impacted by magnetic fields. To analyze such, we first observe the visco-resistive dissipation with varying field strengths in Figs. 3(a) and 3(b). The small-scale dissipation—a proxy of small-scale cascade—enhances with stronger fields, despite the fields reducing the linear growth rate of the instability, which sometimes is argued to lead to subdued turbulent fluctuations as the energy extraction rate by the unstable mode becomes lower.

To understand such enhanced dissipation rate, it may be instructive to compute the energy transfer rates between the background flow and the fluctuations due to both the unstable and stable eigenmodes. As seen in Fig. 3(c), the (time-averaged) energy injection rate by the instability, *Q*_{1}, decreases with stronger magnetic fields, as anticipated. However, the energy return rate by the stable modes, *Q*_{2}, is impacted more by the stronger fields, allowing larger energy cascade to small scales ($\u221dQ1\u2212Q2$) that manifests as enhanced dissipation rate at such scales.^{43} Thus, a stronger magnetic field, until it nearly stabilizes the system, allows more net energy extraction from the free-energy source in the *non*linear state. This questions models that assume the linear growth rate of an instability as a surrogate for the cascade rate or a proxy for nonlinear time, see, e.g., Refs. 44–47. Predictions based upon conventional wisdom of instability-saturation are, thus, challenged.

The ratio $Q2/Q1$ of the stable-eigenmode efficiency to the unstable-eigenmode efficiency in transferring energy between the fluctuations and the background flow is more than 60% even for the strongest magnetic fields in Fig. 3(c), see the red dashed line. Such a strong Lorentz back-reaction is most potent in disrupting the vortex^{48} and creating more small scales, leading to more reduced *Q*_{2} than the corresponding reduction in *Q*_{1}.

We now show visualizations of a highly effective role of stable modes in structuring the magnetic fields. After carrying out an initial-value simulation up to a time when the dynamics are nonlinear and turbulent, we perform two distinct simulation continuations—one unchanged with the stable modes untouched and another with the stable modes removed. For the latter, we set only the large-scale ($0<|kx|<1$) stable-eigenmode amplitudes to zero at the instant when the simulation is restarted. The differences immediately afterward are significant: Removing stable eigenmodes rapidly enhances dissipation rates, see Figs. 4(a) and 4(b) (Multimedia view). To understand this increase in the dissipation, we plot the squared current density $|\u22072\psi |2$ at small scales ($|kx|>1$) in Figs. 4(c)–4(f). Observe in Figs. 4(d) and 4(e) that the magnetic fields get rapidly distorted by the unstable modes of the flow, creating spirals as the counter-straining motion of the stable modes is removed. Note that the (un)stable modes have the largest strain around $(x,z)\u2248(10,0)$ for the shown time, as can also be seen in the figure in the supplementary material. The stable modes here are deleted only at the instant when the simulation is paused; they are quickly nonlinearly driven back. In another simulation where the stable modes are untouched (retained), the spirals do not appear, however, and the time evolution immediately afterward remains almost identical to the initial stage [cf. Figs. 4(f) and 4(c)]. This exercise of instantaneous deletion of stable modes exposes the true nonlinear impact of the stable modes on magnetic field evolution, which remained hidden in Fig. 1(a).

This numerical experiment illustrates that the energy injected by the unstable modes into the large-scale fluctuations ($|kx|<1$) via mean-fluctuation interactions, *in the absence of stable modes*, cascades in its entirety to smaller scales nonlinearly where it dissipates visco-resistively. This is a statement of the classical small-scale cascade by the advective nonlinearity, $(v\xb7\u2207)\psi $, as commonly understood. When **v** is composed of comparable unstable and stable mode-amplitudes, though, the cascade is weakened by a factor of up to ten [Figs. 2, 3(c), and 4(b)]. Such a finding informs and can be used to improve reduced models of geo- and astrophysical instability-driven turbulence,^{44,50} where all the energy injected by the instability is assumed to pass through an inertial cascade, e.g., Refs. 44–47.

It may be noted that, although 2D and 3D turbulence have different conserved quantities and corresponding cascade processes, the stable modes of this study are the conjugate roots of the KH-instability, which exist in both 2D and 3D.^{39} Our preliminary studies of 3D turbulence have indeed shown that virtually identical stable mode excitation and its impact in turbulence occur in 3D as well. Details of such an investigation will be left for a future publication.^{49}

This Letter demonstrates, for the first time, that stable modes via their counteracting straining motion of the flow, sequester magnetic energy at large scales despite the conventional association of shear-flow with intensification and generation of small scales in MHD turbulence. Our analyses show that the stable modes greatly weaken the magnetic energy cascade to small scales and, hence, can thwart the magnetic-energy accumulation at sub-viscous scales in the high-magnetic-Prandtl-number regime, typically encountered in diffuse astrophysical plasmas, thus possibly allowing operation of a dynamo with magnetic fields concentrated at large scales, in accord with astrophysical observations.^{6,10} This work opens a new direction toward such studies.

See the supplementary material for the details of (I) simulation set-up, (II) modal energy evolution equation, (III) energy transfer rates between the background flow and the fluctuations, and (IV) residuals in low-order representation of the turbulent flow in Fig. 1(b). A simulation movie supplements Fig. 4 of this Letter.

The authors thank K. Burns and Dedalus developers for assistance. The authors also thank N. Hurst for offering helpful discussions and a careful reading of the manuscript. This material is based upon the work funded by the U.S. Department of Energy (No. DE-SC0022257) through the National Science Foundation (NSF)/DOE Partnership in Basic Plasma Science and Engineering. We also gratefully acknowledge the support by the U.S. DOE under Award No. DE-FG02-04ER54742 and by the U.S. NSF under Grant Nos. AST-1814327 and AST-1908338. All simulations reported herein were performed using the supercomputing resources provided by the Extreme Science and Engineering Discovery Environment via the Allocation No. TG-PHY130027.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

### Author Contributions

**B. Tripathi:** Conceptualization (equal); Data curation (lead); Formal analysis (lead); Investigation (lead); Methodology (equal); Visualization (lead); Writing - original draft (lead); Writing - review and editing (equal). **A. E. Fraser:** Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Validation (equal); Writing - review and editing (equal). **P. W. Terry:** Conceptualization (equal); Formal analysis (equal); Funding acquisition (lead); Investigation (equal); Methodology (equal); Resources (equal); Supervision (lead); Validation (equal); Writing - review and editing (equal). **E. G. Zweibel:** Conceptualization (equal); Formal analysis (equal); Funding acquisition (lead); Investigation (equal); Methodology (equal); Resources (equal); Supervision (lead); Validation (equal); Writing - review and editing (equal). **M. J. Pueschel:** Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Resources (equal); Supervision (equal); Validation (equal); Writing - review and editing (equal).

## DATA AVAILABILITY

The data that support the findings of this study are available from the corresponding author upon reasonable request.

## References

*operator*respects the parity-time-reversal symmetry operation, and hence, its eigenmodes can be proven to form a complete basis, see Ref. 31.

*Q*

_{1}and

*Q*

_{2}where they refer to the transfer rates between the fluctuations and the

*instantaneous*mean profiles.