scholarly journals The influence of a mean magnetic field on three-dimensional magnetohydrodynamic turbulence

1994 ◽  
Vol 280 ◽  
pp. 95-117 ◽  
Author(s):  
Sean Oughton ◽  
Eric R. Priest ◽  
William H. Matthaeus

Building on results from two-dimensional magnetohydrodynamic (MHD) turbulence (Shebalin, Matthaeus & Montgomery 1983), the development of anisotropic states from initially isotropic ones is investigated numerically for fully three-dimensional incompressible MHD turbulence. It is found that when an external d.c. magnetic field (B0) is imposed on viscous and resistive MHD systems, excitations are preferentially transferred to modes with wavevectors perpendicular to B0). The anisotropy increases with increasing mechanical and magnetic Reynolds numbers, and also with increasing wavenumber. The tendency of B0 to inhibit development of turbulence is also examined.

1983 ◽  
Vol 29 (3) ◽  
pp. 525-547 ◽  
Author(s):  
John V. Shebalin ◽  
William H. Matthaeus ◽  
David Montgomery

The development of anisotropy in an initially isotropie spectrum is studied numerically for two-dimensional magnetohydrodynamic turbulence. The anisotropy develops through the combined effects of an externally imposed d.c. magnetic field and viscous and resistive dissipation at high wavenumbers. The effect is most pronounced at high mechanical and magnetic Reynolds numbers. The anisotropy is greater at the higher wavenumbers.


2019 ◽  
Vol 876 ◽  
pp. 5-18 ◽  
Author(s):  
Riddhi Bandyopadhyay ◽  
William H. Matthaeus ◽  
Sean Oughton ◽  
Minping Wan

In an earlier paper (Wan et al., J. Fluid Mech., vol. 697, 2012, pp. 296–315), the authors showed that a similarity solution for anisotropic incompressible three-dimensional magnetohydrodynamic (MHD) turbulence, in the presence of a uniform mean magnetic field $\boldsymbol{B}_{0}$, exists if the ratio of parallel to perpendicular (with respect to $\boldsymbol{B}_{0}$) similarity length scales remains constant in time. This conjecture appears to be a rather stringent constraint on the dynamics of decay of the energy-containing eddies in MHD turbulence. However, we show here, using direct numerical simulations, that this hypothesis is indeed satisfied in incompressible MHD turbulence. After an initial transient period, the ratio of parallel to perpendicular length scales fluctuates around a steady value during the decay of the eddies. We show further that a Taylor–Kármán-like similarity decay holds for MHD turbulence in the presence of a mean magnetic field. The effect of different parameters, including Reynolds number, mean field strength, and cross-helicity, on the nature of similarity decay is discussed.


1998 ◽  
Vol 358 ◽  
pp. 299-333 ◽  
Author(s):  
OLEG ZIKANOV ◽  
ANDRE THESS

The transformation of initially isotropic turbulent flow of electrically conducting incompressible viscous fluid under the influence of an imposed homogeneous magnetic field is investigated using direct numerical simulation. Under the assumption of large kinetic and small magnetic Reynolds numbers (magnetic Prandtl number Pm[Lt ]1) the quasi-static approximation is applied for the computation of the magnetic field fluctuations. The flow is assumed to be homogeneous and contained in a three-dimensional cubic box with periodic boundary conditions. Large-scale forcing is applied to maintain a statistically steady level of the flow energy. It is found that the pathway traversed by the flow transformation depends decisively on the magnetic interaction parameter (Stuart number). If the magnetic interaction number is small the flow remains three-dimensional and turbulent and no detectable deviation from isotropy is observed. In the case of a strong magnetic field (large magnetic interaction parameter) a rapid transformation to a purely two-dimensional steady state is obtained in agreement with earlier analytical and numerical results for decaying MHD turbulence. At intermediate values of the magnetic interaction parameter the system exhibits intermittent behaviour, characterized by organized quasi-two-dimensional evolution lasting several eddy-turnover times, which is interrupted by strong three-dimensional turbulent bursts. This result implies that the conventional picture of steady angular energy transfer in MHD turbulence must be refined. The spatial structure of the steady two-dimensional final flow obtained in the case of large magnetic interaction parameter is examined. It is found that due to the type of forcing and boundary conditions applied, this state always occurs in the form of a square periodic lattice of alternating vortices occupying the largest possible scale. The stability of this flow to three-dimensional perturbations is analysed using the energy stability method.


2014 ◽  
Vol 761 ◽  
pp. 168-205 ◽  
Author(s):  
Alban Pothérat ◽  
Rico Klein

AbstractMagnetohydrodynamic (MHD) turbulence at low magnetic Reynolds number is experimentally investigated by studying a liquid metal flow in a cubic domain. We focus on the mechanisms that determine whether the flow is quasi-two-dimensional, three-dimensional or in any intermediate state. To this end, forcing is applied by injecting a DC current $I$ through one wall of the cube only, to drive vortices spinning along the magnetic field. Depending on the intensity of the externally applied magnetic field, these vortices extend part or all of the way through the cube. Driving the flow in this way allows us to precisely control not only the forcing intensity but also its dimensionality. A comparison with the theoretical analysis of this configuration singles out the influences of the walls and of the forcing on the flow dimensionality. Flow dimensionality is characterised in several ways. First, we show that when inertia drives three-dimensionality, the velocity near the wall where current is injected scales as $U_{b}\sim I^{2/3}$. Second, we show that when the distance $l_{z}$ over which momentum diffuses under the action of the Lorentz force (Sommeria & Moreau, J. Fluid Mech., vol. 118, 1982, pp. 507–518) reaches the channel width $h$, the velocity near the opposite wall $U_{t}$ follows a similar law with a correction factor $(1-h/l_{z})$ that measures three-dimensionality. When $l_{z}<h$, by contrast, the opposite wall has less influence on the flow and $U_{t}\sim I^{1/2}$. The central role played by the ratio $l_{z}/h$ is confirmed by experimentally verifying the scaling $l_{z}\sim N^{1/2}$ put forward by Sommeria & Moreau ($N$ is the interaction parameter) and, finally, the nature of the three-dimensionality involved is further clarified by distinguishing weak and strong three-dimensionalities previously introduced by Klein & Pothérat (Phys. Rev. Lett., vol. 104 (3), 2010, 034502). It is found that both types vanish only asymptotically in the limit $N\rightarrow \infty$. This provides evidence that because of the no-slip walls, (i) the transition between quasi-two-dimensional and three-dimensional turbulence does not result from a global instability of the flow, unlike in domains with non-dissipative boundaries (Boeck et al. Phys. Rev. Lett., vol. 101, 2008, 244501), and (ii) it does not occur simultaneously at all scales.


2015 ◽  
Vol 783 ◽  
pp. 605-636 ◽  
Author(s):  
Alban Pothérat ◽  
Kacper Kornet

We present direct numerical simulations of decaying magnetohydrodynamic (MHD) turbulence at low magnetic Reynolds number. The domain considered is bounded by periodic boundary conditions in the two directions perpendicular to the magnetic field and by two plane Hartmann walls in the third direction. Regimes of high magnetic fields (Hartmann number of up to 896) are reached thanks to a new spectral method using the eigenvectors of the dissipation operator. The decay is found to proceed through two phases: first, energy and integral length scales vary rapidly during a two-dimensionalisation phase extending over approximately a Hartmann friction time. During this phase, the evolution of the former appears significantly more impeded by the presence of walls than that of the latter. Once the large scales are nearly quasi-two-dimensional, the decay results from the competition of a two-dimensional dynamics driven by dissipation in the Hartmann boundary layers and the three-dimensional dynamics of smaller scales. In the later stages of the decay, three-dimensionality subsists under the form of barrel-shaped structures. A purely quasi-two dimensional decay entirely dominated by friction in the Hartmann layers is not reached because of residual dissipation in the bulk. However, this dissipation is not generated by the three-dimensionality that subsists, but by residual viscous friction due to horizontal velocity gradients. In the process, the energy in the velocity component aligned with the magnetic field is found to be strongly suppressed, as is skewness. This result reproduces the experimental findings of Kolesnikov & Tsinober (Fluid Dyn., vol. 9, 1974, pp. 621–624), where, as in the present simulations, Hartmann walls were present.


2000 ◽  
Vol 64 (2) ◽  
pp. 179-193 ◽  
Author(s):  
SEAN OUGHTON ◽  
ROSSELLA PRANDI

The issue of dynamical anisotropy in helical three-dimensional magnetohydrodynamic turbulence with a mean magnetic field B0 is investigated. Using high-resolution direct numerical simulations, we follow the evolution of various isotropic initial states characterized by their different values of the kinetic helicity. The cross helicity and magnetic helicity of the initial conditions are also varied. In agreement with earlier work, we find that such initial states become anisotropic in of order an eddy-turnover time, with correlation lengths parallel to B0 remaining largely unchanged while finer scales are excited in the perpendicular directions. Moreover, it is found that the development of both the anisotropy and the energy are essentially independent of the initial level of kinetic helicity. The physics associated with this latter feature is discussed.


Author(s):  
Francine Battaglia ◽  
George Papadopoulos

The effect of three-dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio of in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations fail to capture completely the total expansion effect on the flow, which couples both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When the two-dimensional geometry was consistent with the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations was compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, helping to sustain laminar flow symmetry to higher Reynolds numbers in comparison to nominally two-dimensional double-expansion geometries. Lastly, and most importantly, when the logarithm of the critical Reynolds number from all these studies was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.


1969 ◽  
Vol 3 (2) ◽  
pp. 255-267 ◽  
Author(s):  
M. P. Srivastava ◽  
P. K. Bhat

We have studied the behaviour of a charged particle in an axially symmetric magnetic field having a neutral point, so as to find a possibility of confining a charged particle in a thermonuclear device. In order to study the motion we have reduced a three-dimensional motion to a two-dimensional one by introducing a fictitious potential. Following Schmidt we have classified the motion, as an ‘off-axis motion’ and ‘encircling motion’ depending on the behaviour of this potential. We see that the particle performs a hybrid type of motion in the negative z-axis, i.e. at some instant it is in ‘off-axis motion’ while at another instant it is in ‘encircling motion’. We have also solved the equation of motion numerically and the graphs of the particle trajectory verify our analysis. We find that in most of the cases the particle is contained. The magnetic moment is found to be moderately adiabatic.


2017 ◽  
Vol 825 ◽  
pp. 631-650 ◽  
Author(s):  
Francesco Romanò ◽  
Arash Hajisharifi ◽  
Hendrik C. Kuhlmann

The topology of the incompressible steady three-dimensional flow in a partially filled cylindrical rotating drum, infinitely extended along its axis, is investigated numerically for a ratio of pool depth to radius of 0.2. In the limit of vanishing Froude and capillary numbers, the liquid–gas interface remains flat and the two-dimensional flow becomes unstable to steady three-dimensional convection cells. The Lagrangian transport in the cellular flow is organised by periodic spiralling-in and spiralling-out saddle foci, and by saddle limit cycles. Chaotic advection is caused by a breakup of a degenerate heteroclinic connection between the two saddle foci when the flow becomes three-dimensional. On increasing the Reynolds number, chaotic streamlines invade the cells from the cell boundary and from the interior along the broken heteroclinic connection. This trend is made evident by computing the Kolmogorov–Arnold–Moser tori for five supercritical Reynolds numbers.


Author(s):  
Jasim Mohmed Jasim Jasim ◽  
Iryna Shvedchykova ◽  
Igor Panasiuk ◽  
Julia Romanchenko ◽  
Inna Melkonova

An approach is proposed to carry out multivariate calculations of the magnetic field distribution in the working gaps of a plate polygradient matrix of an electromagnetic separator, based on a combination of the advantages of two- and three-dimensional computer modeling. Two-dimensional geometric models of computational domains are developed, which differ in the geometric dimensions of the plate matrix elements and working air gaps. To determine the vector magnetic potential at the boundaries of two-dimensional computational domains, a computational 3D experiment is carried out. For this, three variants of the electromagnetic separator are selected, which differ in the size of the working air gaps of the polygradient matrices. For them, three-dimensional computer models are built, the spatial distribution of the magnetic field in the working intervals of the electromagnetic separator matrix and the obtained numerical values of the vector magnetic potential at the boundaries of the computational domains are investigated. The determination of the values of the vector magnetic potential for all other models is carried out by interpolation. The obtained values of the vector magnetic potential are used to set the boundary conditions in a computational 2D experiment. An approach to the choice of a rational version of a lamellar matrix is substantiated, which provides a solution to the problem according to the criterion of the effective area of the working area. Using the method of simple enumeration, a variant of the structure of a polygradient matrix with rational geometric parameters is selected. The productivity of the electromagnetic separator with rational geometric parameters of the matrix increased by 3–5 % with the same efficiency of extraction of ferromagnetic inclusions in comparison with the basic version of the device


Sign in / Sign up

Export Citation Format

Share Document