Magnetohydrodynamic turbulence
Biskamp, Dieter
2003-01-01
This book presents an introduction to, and modern account of, magnetohydrodynamic (MHD) turbulence, an active field both in general turbulence theory and in various areas of astrophysics. The book starts by introducing the MHD equations, certain useful approximations and the transition to turbulence. The second part of the book covers incompressible MHD turbulence, the macroscopic aspects connected with the different self-organization processes, the phenomenology of the turbulence spectra, two-point closure theory, and intermittency. The third considers two-dimensional turbulence and compressi
Magnetohydrodynamic Turbulence
Montgomery, David C.
2004-01-01
Magnetohydrodynamic (MHD) turbulence theory is modeled on neutral fluid (Navier-Stokes) turbulence theory, but with some important differences. There have been essentially no repeatable laboratory MHD experiments wherein the boundary conditions could be controlled or varied and a full set of diagnostics implemented. The equations of MHD are convincingly derivable only in the limit of small ratio of collision mean-free-paths to macroscopic length scales, an inequality that often goes the other way for magnetofluids of interest. Finally, accurate information on the MHD transport coefficients-and thus, the Reynolds-like numbers that order magnetofluid behavior-is largely lacking; indeed, the algebraic expressions used for such ingredients as the viscous stress tensor are often little more than wishful borrowing from fluid mechanics. The one accurate thing that has been done extensively and well is to solve the (strongly nonlinear) MHD equations numerically, usually in the presence of rectangular periodic boundary conditions, and then hope for the best when drawing inferences from the computations for those astrophysical and geophysical MHD systems for which some indisputably turbulent detailed data are available, such as the solar wind or solar prominences. This has led to what is perhaps the first field of physics for which computer simulations are regarded as more central to validating conclusions than is any kind of measurement. Things have evolved in this way due to a mixture of the inevitable and the bureaucratic, but that is the way it is, and those of us who want to work on the subject have to live with it. It is the only game in town, and theories that have promised more-often on the basis of some alleged ``instability''-have turned out to be illusory.
Turbulence in magnetohydrodynamics
Beresnyak, Andrey
2016-01-01
Magnetohydrodynamics describes dynamics in electrically conductive fluids. These occur in our environment as well as in our atmosphere and magnetosphere, and play a role in the sun's interaction with our planet. This work gives the basic information on turbulence in nature, comprising the needed equations, notions and numerical simulations. The current state of our knowledge and future implications of MHD turbulence are outlined systematically. It is indispensable for all scientists engaged in research of our atmosphere and in space science.
Scalings of intermittent structures in magnetohydrodynamic turbulence
Zhdankin, Vladimir; Uzdensky, Dmitri A
2016-01-01
Turbulence is ubiquitous in plasmas, leading to rich dynamics characterized by irregularity, irreversibility, energy fluctuations across many scales, and energy transfer across many scales. Another fundamental and generic feature of turbulence, although sometimes overlooked, is the inhomogeneous dissipation of energy in space and in time. This is a consequence of intermittency, the scale-dependent inhomogeneity of dynamics caused by fluctuations in the turbulent cascade. Intermittency causes turbulent plasmas to self-organize into coherent dissipative structures, which may govern heating, diffusion, particle acceleration, and radiation emissions. In this paper, we present recent progress on understanding intermittency in incompressible magnetohydrodynamic turbulence with a strong guide field. We focus on the statistical analysis of intermittent dissipative structures, which occupy a small fraction of the volume but arguably account for the majority of energy dissipation. We show that, in our numerical simulat...
Magnetohydrodynamics turbulence: An astronomical perspective
S Sridhar
2011-07-01
Early work on magnetohydrodynamic (MHD) turbulence in the 1960s due, independently, to Iroshnikov and Kraichnan (IK) considered isotropic inertial-range spectra. Whereas laboratory experiments were not in a position to measure the spectral index, they showed that the turbulence was strongly anisotropic. Theoretical horizons correspondingly expanded in the 1980s, to accommodate both the isotropy of the IK theory and the anisotropy suggested by the experiments. Since the discovery of pulsars in 1967, many years of work on interstellar scintillation suggested that small-scale interstellar turbulence must have a hydromagnetic origin; but the IK spectrum was too ﬂat and the ideas on anisotropic spectra too qualitative to explain the observations. In response, new theories of balanced MHD turbulence were proposed in the 1990s, which argued that the IK theory was incorrect, and made quantitative predictions of anisotropic inertial-range spectra; these theories have since found applications in many areas of astrophysics. Spacecraft measurements of solar-wind turbulence show that there is more power in Alfvén waves that travel away from the Sun than towards it. Theories of imbalanced MHD turbulence have now been proposed to address interplanetary turbulence. This very active area of research continues to be driven by astronomy.
Magnetohydrodynamic turbulence: Observation and experiment
Brown, M. R.; Schaffner, D. A.; Weck, P. J. [Department of Physics and Astronomy, Swarthmore College, 500 College Avenue, Swarthmore, Pennsylvania 19081 (United States)
2015-05-15
We provide a tutorial on the paradigms and tools of magnetohydrodynamic (MHD) turbulence. The principal paradigm is that of a turbulent cascade from large scales to small, resulting in power law behavior for the frequency power spectrum for magnetic fluctuations E{sub B}(f). We will describe five useful statistical tools for MHD turbulence in the time domain: the temporal autocorrelation function, the frequency power spectrum, the probability distribution function of temporal increments, the temporal structure function, and the permutation entropy. Each of these tools will be illustrated with an example taken from MHD fluctuations in the solar wind. A single dataset from the Wind satellite will be used to illustrate all five temporal statistical tools.
Anisotropic Intermittency of Magnetohydrodynamic Turbulence
Osman, K T; Chapman, S C; Hnat, B
2013-01-01
A higher-order multiscale analysis of spatial anisotropy in inertial range magnetohydrodynamic turbulence is presented using measurements from the STEREO spacecraft in fast ambient solar wind. We show for the first time that, when measuring parallel to the local magnetic field direction, the full statistical signature of the magnetic and Els\\"asser field fluctuations is that of a non-Gaussian globally scale-invariant process. This is distinct from the classic multi-exponent statistics observed when the local magnetic field is perpendicular to the flow direction. These observations are interpreted as evidence for the weakness, or absence, of a parallel magnetofluid turbulence energy cascade. As such, these results present strong observational constraints on the statistical nature of intermittency in turbulent plasmas.
Studying Magnetohydrodynamic Turbulence with Synchrotron Polarization Dispersion
Zhang, Jian-Fu; Lazarian, Alex; Lee, Hyeseung; Cho, Jungyeon
2016-07-01
We test a new technique for studying magnetohydrodynamic turbulence suggested by Lazarian & Pogosyan, using synthetic observations of synchrotron polarization. This paper focuses on a one-point statistics, which is termed polarization frequency analysis, that is characterized by the variance of polarized emission as a function of the square of the wavelength along a single line of sight. We adopt the ratio η of the standard deviation of the line-of-sight turbulent magnetic field to the line-of-sight mean magnetic field to depict the level of turbulence. When this ratio is large (η \\gg 1), which characterizes a region dominated by turbulent field, or small (η ≲ 0.2), which characterizes a region dominated by the mean field, we obtain the polarization variance \\propto {λ }-2 or \\propto {λ }-2-2m, respectively. At small η, i.e., in the region dominated by the mean field, we successfully recover the turbulent spectral index from the polarization variance. We find that our simulations agree well with the theoretical prediction of Lazarian & Pogosyan. With existing and upcoming data cubes from the Low-Frequency Array for Radio Astronomy (LOFAR) and the Square Kilometer Array (SKA), this new technique can be applied to study the magnetic turbulence in the Milky Way and other galaxies.
Magnetohydrodynamic turbulent model for LMFRS
Ohira, H.; Ara, K. [O-arai Engineering Center, Japan Nuclear Cycle Development Institute, O-arai (Japan)
2004-07-01
In order to study magnetohydrodynamic behavior in electromagnetic pumps, electromagnetic flow meters, etc. for Liquid Metal Fast Reactors (LMFR), a large eddy simulation method using an artificial wall boundary condition was developed. In this study, Spalding's law of the wall and the eddy viscosity for uniform magnetic fields, which was proposed by Shimomura, was applied to Finite Element Method of Generalized Simplified Marker and Cell (GSMAC-FEM). We calculated MHD channel flow in various element sizes on the conditions of Hartmann numbers of 0, 52.5 and 125, whose Reynolds numbers based on the average velocity were all about 29000. These results showed the average velocity profiles were in good agreement with both the experimental results by Brouillette-Lykoudis and the detail calculation results by Shimomura, although farther calculations were needed to verify the turbulence intensities. (authors)
Magnetohydrodynamic turbulent model for LMFRS
In order to study magnetohydrodynamic behavior in electromagnetic pumps, electromagnetic flow meters, etc. for Liquid Metal Fast Reactors (LMFR), a large eddy simulation method using an artificial wall boundary condition was developed. In this study, Spalding's law of the wall and the eddy viscosity for uniform magnetic fields, which was proposed by Shimomura, was applied to Finite Element Method of Generalized Simplified Marker and Cell (GSMAC-FEM). We calculated MHD channel flow in various element sizes on the conditions of Hartmann numbers of 0, 52.5 and 125, whose Reynolds numbers based on the average velocity were all about 29000. These results showed the average velocity profiles were in good agreement with both the experimental results by Brouillette-Lykoudis and the detail calculation results by Shimomura, although farther calculations were needed to verify the turbulence intensities. (authors)
Magnetohydrodynamic Turbulence and the Geodynamo
Shebalin, John V.
2014-01-01
The ARES Directorate at JSC has researched the physical processes that create planetary magnetic fields through dynamo action since 2007. The "dynamo problem" has existed since 1600, when William Gilbert, physician to Queen Elizabeth I, recognized that the Earth was a giant magnet. In 1919, Joseph Larmor proposed that solar (and by implication, planetary) magnetism was due to magnetohydrodynamics (MHD), but full acceptance did not occur until Glatzmaier and Roberts solved the MHD equations numerically and simulated a geomagnetic reversal in 1995. JSC research produced a unique theoretical model in 2012 that provided a novel explanation of these physical observations and computational results as an essential manifestation of broken ergodicity in MHD turbulence. Research is ongoing, and future work is aimed at understanding quantitative details of magnetic dipole alignment in the Earth as well as in Mercury, Jupiter and its moon Ganymede, Saturn, Uranus, Neptune, and the Sun and other stars.
Studying Magnetohydrodynamic Turbulence with Synchrotron Polarization Dispersion
Zhang, Jian-Fu; Lee, Hyeseung; Cho, Jungyeon
2016-01-01
We test a new technique of studying magnetohydrodynamic (MHD) turbulence suggested by Lazarian \\& Pogosyan, using synthetic synchrotron polarization observations. This paper focuses on a one-point statistics, which is termed the polarization frequency analysis, that is characterized by the variance of polarized emission as a function of the square of wavelengths along a single line of sight. We adopt a ratio $\\eta$ of the standard deviation of the line-of-sight turbulent magnetic field to the line-of-sight mean magnetic field to depict the level of turbulence. When this ratio is either large ($\\eta\\gg1$), which characterizes a turbulent field dominated region, or small ($\\eta\\lesssim0.2$), which characterizes a mean field dominated region, we obtain the polarization variance $\\left\\propto\\lambda^{-2}$ and $\\left\\propto\\lambda^{-2-2m}$, respectively. At small $\\eta$, i.e., the mean field dominated region, we successfully recover the turbulent spectral index by the polarization variance. We find that our si...
Electron magnetohydrodynamics: dynamics and turbulence.
Lyutikov, Maxim
2013-11-01
We consider dynamics and turbulent interaction of whistler modes within the framework of inertialess electron magnetohydrodynamics (EMHD). We argue that there is no energy principle in EMHD: any stationary closed configuration is neutrally stable. On the other hand, the relaxation principle, the long term evolution of a weakly dissipative system towards Taylor-Beltrami state, remains valid in EMHD. We consider the turbulent cascade of whistler modes. We show that (i) harmonic whistlers are exact nonlinear solutions; (ii) collinear whistlers do not interact (including counterpropagating); (iii) waves with the same value of the wave vector k(1)=k(2) do not interact; (iv) whistler modes have a dispersion that allows a three-wave decay, including into a zero frequency mode; (v) the three-wave interaction effectively couples modes with highly different wave numbers and propagation angles. In addition, linear interaction of a whistler with a single zero mode can lead to spatially divergent structures via parametric instability. All these properties are drastically different from MHD, so that the qualitative properties of the Alfvén turbulence can not be transferred to the EMHD turbulence. We derive the Hamiltonian formulation of EMHD, and using Bogoliubov transformation reduce it to the canonical form; we calculate the matrix elements for the three-wave interaction of whistlers. We solve numerically the kinetic equation and show that, generally, the EMHD cascade develops within a broad range of angles, while transiently it may show anisotropic, nearly two-dimensional structures. Development of a cascade depends on the forcing (nonuniversal) and often fails to reach a steady state. Analytical estimates predict the spectrum of magnetic fluctuations for the quasi-isotropic cascade [proportionality]k(-2). The cascade remains weak (not critically balanced). The cascade is UV local, while the infrared locality is weakly (logarithmically) violated. PMID:24329368
Scalings of intermittent structures in magnetohydrodynamic turbulence
Zhdankin, Vladimir; Boldyrev, Stanislav; Uzdensky, Dmitri A.
2016-05-01
Turbulence is ubiquitous in plasmas, leading to rich dynamics characterized by irregularity, irreversibility, energy fluctuations across many scales, and energy transfer across many scales. Another fundamental and generic feature of turbulence, although sometimes overlooked, is the inhomogeneous dissipation of energy in space and in time. This is a consequence of intermittency, the scale-dependent inhomogeneity of dynamics caused by fluctuations in the turbulent cascade. Intermittency causes turbulent plasmas to self-organize into coherent dissipative structures, which may govern heating, diffusion, particle acceleration, and radiation emissions. In this paper, we present recent progress on understanding intermittency in incompressible magnetohydrodynamic turbulence with a strong guide field. We focus on the statistical analysis of intermittent dissipative structures, which occupy a small fraction of the volume but arguably account for the majority of energy dissipation. We show that, in our numerical simulations, intermittent structures in the current density, vorticity, and Elsässer vorticities all have nearly identical statistical properties. We propose phenomenological explanations for the scalings based on general considerations of Elsässer vorticity structures. Finally, we examine the broader implications of intermittency for astrophysical systems.
Turbulent Magnetohydrodynamic Jet Collimation and Thermal Driving
Williams, Peter T.
2003-01-01
We have argued that magnetohydrodynamic (MHD) turbulence in an accretion disk naturally produces hoop-stresses, and that in a geometrically-thick flow these stresses could both drive and collimate an outflow. We based this argument on an analogy of turbulent MHD fluids to viscoelastic fluids, in which azimuthal shear flow creates hoop-stresses that cause a variety of flow phenomena, including the Weissenberg effect in which a fluid climbs a spinning rod. One of the more important differences ...
Two-dimensional electron magnetohydrodynamic turbulence
A novel type of turbulence, which arises in 2D electron magnetohydrodynamics, is studied by numerical simulation. Energy dissipation rates are found to be independent of the dissipation coefficients. The energy spectrum Ek follows the basic Kolmogorov-type predictions, k-5/3 for kde > 1 and k-7/3 for kde e = electron inertial length) and is hence independent of the linear wave properties. Results are compared with other 2D turbulent systems. (author)
Magnetohydrodynamic turbulence in the solar wind
Goldstein, Melvyn L.
1995-01-01
The fluctuations in magnetic field and plasma velocity in solar wind, which possess many features of fully developed magnetohydrodynamic (MHD) turbulence, are discussed. Direct spacecraft observations from 0.3 to over 20 AU, remote sensing radio scintillation observations, numerical simulations, and various models provide complementary methods that show that the fluctuations in the wind parameters undergo significant dynamical evolution independent of whatever turbulence might exist in the solar photosphere and corona. The Cluster mission, with high time resolution particle and field measurements and its variable separation strategies, should be able to provide data for answering many questions on MHD turbulence.
Magnetohydrodynamic turbulence in warped accretion discs
Torkelsson, U; Brandenburg, A; Pringle, J E; Nordlund, A A; Stein, R F; Nordlund, AA.
2001-01-01
Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.
Weak turbulence in two-dimensional magnetohydrodynamics
Tronko, Natalia; Nazarenko, Sergey; Galtier, Sebastien
2012-01-01
A weak wave turbulence theory is developed for two-dimensional (2D) magnetohydrodynamics (MHD). We derive and analyze the kinetic equation describing the three-wave interactions of pseudo-Alfv\\'en waves. Our analysis is greatly helped by the fortunate fact that in 2D the wave-kinetic equation is integrable. In contrast with the 3D case, in 2D the wave interactions are nonlocal. Another distinct feature is that strong derivatives of spectra tend to appear in the region of small parallel (i.e. ...
Imbalanced relativistic force-free magnetohydrodynamic turbulence
When magnetic energy density is much larger than that of matter, as in pulsar/black hole magnetospheres, the medium becomes force-free and we need relativity to describe it. As in non-relativistic magnetohydrodynamics (MHD), Alfvénic MHD turbulence in the relativistic limit can be described by interactions of counter-traveling wave packets. In this paper, we numerically study strong imbalanced MHD turbulence in such environments. Here, imbalanced turbulence means the waves traveling in one direction (dominant waves) have higher amplitudes than the opposite-traveling waves (sub-dominant waves). We find that (1) spectrum of the dominant waves is steeper than that of sub-dominant waves, (2) the anisotropy of the dominant waves is weaker than that of sub-dominant waves, and (3) the dependence of the ratio of magnetic energy densities of dominant and sub-dominant waves on the ratio of energy injection rates is steeper than quadratic (i.e., b+2/b−2∝(ϵ+/ϵ−)n with n > 2). These results are consistent with those obtained for imbalanced non-relativistic Alfvénic turbulence. This corresponds well to the earlier reported similarity of the relativistic and non-relativistic balanced magnetic turbulence.
Global invariants in ideal magnetohydrodynamic turbulence
Magnetohydrodynamic (MHD) turbulence is an important though incompletely understood factor affecting the dynamics of many astrophysical, geophysical, and technological plasmas. As an approximation, viscosity and resistivity may be ignored, and ideal MHD turbulence may be investigated by statistical methods. Incompressibility is also assumed and finite Fourier series are used to represent the turbulent velocity and magnetic field. The resulting model dynamical system consists of a set of independent Fourier coefficients that form a canonical ensemble described by a Gaussian probability density function (PDF). This PDF is similar in form to that of Boltzmann, except that its argument may contain not just the energy multiplied by an inverse temperature, but also two other invariant integrals, the cross helicity and magnetic helicity, each multiplied by its own inverse temperature. However, the cross and magnetic helicities, as usually defined, are not invariant in the presence of overall rotation or a mean magnetic field, respectively. Although the generalized form of the magnetic helicity is known, a generalized cross helicity may also be found, by adding terms that are linear in the mean magnetic field and angular rotation vectors, respectively. These general forms are invariant even in the presence of overall rotation and a mean magnetic field. We derive these general forms, explore their properties, examine how they extend the statistical theory of ideal MHD turbulence, and discuss how our results may be affected by dissipation and forcing
Structures in magnetohydrodynamic turbulence: detection and scaling.
Uritsky, V M; Pouquet, A; Rosenberg, D; Mininni, P D; Donovan, E F
2010-11-01
We present a systematic analysis of statistical properties of turbulent current and vorticity structures at a given time using cluster analysis. The data stem from numerical simulations of decaying three-dimensional magnetohydrodynamic turbulence in the absence of an imposed uniform magnetic field; the magnetic Prandtl number is taken equal to unity, and we use a periodic box with grids of up to 1536³ points and with Taylor Reynolds numbers up to 1100. The initial conditions are either an X -point configuration embedded in three dimensions, the so-called Orszag-Tang vortex, or an Arn'old-Beltrami-Childress configuration with a fully helical velocity and magnetic field. In each case two snapshots are analyzed, separated by one turn-over time, starting just after the peak of dissipation. We show that the algorithm is able to select a large number of structures (in excess of 8000) for each snapshot and that the statistical properties of these clusters are remarkably similar for the two snapshots as well as for the two flows under study in terms of scaling laws for the cluster characteristics, with the structures in the vorticity and in the current behaving in the same way. We also study the effect of Reynolds number on cluster statistics, and we finally analyze the properties of these clusters in terms of their velocity-magnetic-field correlation. Self-organized criticality features have been identified in the dissipative range of scales. A different scaling arises in the inertial range, which cannot be identified for the moment with a known self-organized criticality class consistent with magnetohydrodynamics. We suggest that this range can be governed by turbulence dynamics as opposed to criticality and propose an interpretation of intermittency in terms of propagation of local instabilities. PMID:21230595
Extended inertial range phenomenology of magnetohydrodynamic turbulence
Matthaeus, William H.; Zhou, YE
1989-01-01
A phenomenological treatment of the inertial range of isotropic statistically steady magnetohydrodynamic turbulence is presented, extending the theory of Kraichnan (1965). The role of Alfven wave propagation is treated on equal footing with nonlinear convection, leading to a simple generalization of the relations between the times characteristic of wave propagation, convection, energy transfer, and decay of triple correlations. The theory leads to a closed-form steady inertial range spectral law that reduces to the Kraichnan and Kolmogorov laws in appropriate limits. The Kraichnan constant is found to be related in a simple way to the Kolmogorov constant; for typical values of the latter constant, the former has values in the range 1.22-1.87. Estimates of the time scale associated with spectral transfer of energy also emerge from the new approach, generalizing previously presented 'golden rules' for relating the spectral transfer time scale to the Alfven and eddy-turnover time scales.
Classes of hydrodynamic and magnetohydrodynamic turbulent decay
Brandenburg, Axel
2016-01-01
We perform numerical simulations of decaying hydrodynamic and magnetohydrodynamic turbulence. We classify our time-dependent solutions by their evolutionary tracks in parametric plots between instantaneous scaling exponents. We find distinct classes of solutions evolving along specific trajectories toward points on a line of self-similar solutions. These trajectories are determined by the underlying physics governing individual cases, and not by the initial conditions, as is widely assumed. In the helical case, even for a scale-invariant initial spectrum (inversely proportional to wavenumber k), the solution evolves along the same trajectory as for a Batchelor spectrum (proportional to k^4). All of our self-similar solutions have an intrinsic subinertial range close to k^4$.
Structure of homogeneous nonhelical magnetohydrodynamic turbulence
Miller, R.S.; Mashayek, F.; Adumitroaie, V.; Givi, P. [Department of Mechanical and Aerospace Engineering, State University of New York at Buffalo, Buffalo, New York 14260-4400 (United States)
1996-09-01
Results are presented for three-dimensional direct numerical simulations of nonhelical magnetohydrodynamic (MHD) turbulence for both stationary isotropic and homogeneous shear flow configurations with zero mean induction and unity magnetic Prandtl number. Small scale dynamo action is observed in both flows, and stationary values for the ratio of magnetic to kinetic energy are shown to scale nearly linearly with the Taylor microscale Reynolds numbers above a critical value of Re{sub {lambda}}{approx_equal}30. The presence of the magnetic field has the effect of decreasing the kinetic energy of the flow, while simultaneously increasing the Taylor microscale Reynolds number due to enlargement of the hydrodynamic length scales. For shear flows, both the velocity and the magnetic fields become increasingly anisotropic with increasing initial magnetic field strength. The kinetic energy spectra show a relative increase in high wave-number energy in the presence of a magnetic field. The magnetic field is found to portray an intermittent behavior, with peak values of the flatness near the critical Reynolds number. The magnetic field of both flows is organized in the form of {open_quote}{open_quote}flux tubes{close_quote}{close_quote} and magnetic {open_quote}{open_quote}sheets.{close_quote}{close_quote} These regions of large magnetic field strength show a small correlation with moderate vorticity regions, while the electric current structures are correlated with large amplitude strain regions of the turbulence. Some of the characteristics of small scale MHD turbulence are explained via the {open_quote}{open_quote}structural{close_quote}{close_quote} description of turbulence. {copyright} {ital 1996 American Institute of Physics.}
Structure of homogeneous nonhelical magnetohydrodynamic turbulence
Results are presented for three-dimensional direct numerical simulations of nonhelical magnetohydrodynamic (MHD) turbulence for both stationary isotropic and homogeneous shear flow configurations with zero mean induction and unity magnetic Prandtl number. Small scale dynamo action is observed in both flows, and stationary values for the ratio of magnetic to kinetic energy are shown to scale nearly linearly with the Taylor microscale Reynolds numbers above a critical value of Reλ≅30. The presence of the magnetic field has the effect of decreasing the kinetic energy of the flow, while simultaneously increasing the Taylor microscale Reynolds number due to enlargement of the hydrodynamic length scales. For shear flows, both the velocity and the magnetic fields become increasingly anisotropic with increasing initial magnetic field strength. The kinetic energy spectra show a relative increase in high wave-number energy in the presence of a magnetic field. The magnetic field is found to portray an intermittent behavior, with peak values of the flatness near the critical Reynolds number. The magnetic field of both flows is organized in the form of open-quote open-quote flux tubes close-quote close-quote and magnetic open-quote open-quote sheets.close-quote close-quote These regions of large magnetic field strength show a small correlation with moderate vorticity regions, while the electric current structures are correlated with large amplitude strain regions of the turbulence. Some of the characteristics of small scale MHD turbulence are explained via the open-quote open-quote structural close-quote close-quote description of turbulence. copyright 1996 American Institute of Physics
Structures in magnetohydrodynamic turbulence: detection and scaling
Uritsky, Vadim M; Rosenberg, Duane; Mininni, Pablo D; Donovan, Eric
2010-01-01
We present a systematic analysis of statistical properties of turbulent current and vorticity structures at a given time using cluster analysis. The data stems from numerical simulations of decaying three-dimensional (3D) magnetohydrodynamic turbulence in the absence of an imposed uniform magnetic field; the magnetic Prandtl number is taken equal to unity, and we use a periodic box with grids of up to 1536^3 points, and with Taylor Reynolds numbers up to 1100. The initial conditions are either an X-point configuration embedded in 3D, the so-called Orszag-Tang vortex, or an Arn'old-Beltrami-Childress configuration with a fully helical velocity and magnetic field. In each case two snapshots are analyzed, separated by one turn-over time, starting just after the peak of dissipation. We show that the algorithm is able to select a large number of structures (in excess of 8,000) for each snapshot and that the statistical properties of these clusters are remarkably similar for the two snapshots as well as for the two...
Magnetic Discontinuities in Magnetohydrodynamic Turbulence and in the Solar Wind
Zhdankin, Vladimir; Boldyrev, Stanislav; Mason, Joanne; Perez, Jean Carlos
2012-01-01
Recent measurements of solar wind turbulence report the presence of intermittent, exponentially distributed angular discontinuities in the magnetic field. In this Letter, we study whether such discontinuities can be produced by magnetohydrodynamic (MHD) turbulence. We detect the discontinuities by measuring the fluctuations of the magnetic field direction, Delta theta, across fixed spatial increments Delta x in direct numerical simulations of MHD turbulence with an imposed uniform guide field...
Generalized similarity in finite range solar wind magnetohydrodynamic turbulence
Chapman, S. C.; Nicol, R. M.
2009-01-01
Extended or generalized similarity is a ubiquitous but not well understood feature of turbulence that is realized over a finite range of scales. ULYSSES spacecraft solar polar passes at solar minimum provide \\textit{in situ} observations of evolving anisotropic magnetohydrodynamic turbulence in the solar wind under ideal conditions of fast quiet flow. We find a single generalized scaling function characterises this finite range turbulence and is insensitive to plasma conditions. The recent un...
The small-scale turbulent dynamo in smoothed particle magnetohydrodynamics
Tricco, Terrence S.; Price, Daniel J.; Federrath, Christoph
2016-01-01
Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small-scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of star-forming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few ...
The inverse cascade of magnetic helicity in magnetohydrodynamic turbulence
Müller, W.; Malapaka, S.; Busse, A.
2012-01-01
The nonlinear dynamics of magnetic helicity, $H^M$, which is responsible for large-scale magnetic structure formation in electrically conducting turbulent media is investigated in forced and decaying three-dimensional magnetohydrodynamic turbulence. This is done with the help of high resolution direct numerical simulations and statistical closure theory. The numerically observed spectral scaling of $H^M$ is at variance with earlier work using a statistical closure model [Pouquet et al., J. Fl...
Intermittency of quasi-static magnetohydrodynamic turbulence: A wavelet viewpoint
Intermittency of quasi-static magnetohydrodynamic (MHD) turbulence in an imposed magnetic field is examined, using three-dimensional orthonormal wavelets. The wavelet analysis is applied to two turbulent MHD flows computed by direct numerical simulation with 5123 grid points and with different intensities of the imposed magnetic field. It is found that the imposed magnetic field leads to a substantial amplification of intermittency of the flow, especially in the direction of the imposed magnetic field.
The residual energy in freely decaying magnetohydrodynamic turbulence
Within the framework of the eddy-damped quasi-normal Markovian approximation for incompressible isotropic magnetohydrodynamic (MHD) turbulence a prediction for the inertial range scaling of the residual energy spectrum, ERk= vertical bar Ekkinetic - Emagnetick vertical bar ∼ k-7/3, is obtained. This scaling, while in contradiction to earlier theoretical results, is shown to be in agreement with high-resolution direct numerical simulations of nonhelical decaying MHD turbulence. The underlying phenomenology states a dynamic quasi-equilibrium of the small-scale turbulent dynamo and the Alfven effect
Perpendicular diffusion of energetic particles in noisy reduced magnetohydrodynamic turbulence
Shalchi, Andreas; Hussein, Mohammad
2014-01-01
Recently a model for noisy reduced magnetohydrodynamic turbulence was proposed. The latter model was already used to study the random walk of magnetic field lines. In the current article we use the same model to investigate the diffusion of energetic particles across the mean magnetic field. To compute the perpendicular diffusion coefficient two analytical theories are used, namely the Non- Linear Guiding Center (NLGC) theory and the Unified Non-Linear Transport (UNLT) theory. It is shown tha...
Current and vorticity dynamics in three-dimensional magnetohydrodynamic turbulence
Spectral numerical simulations of homogeneous incompressible magnetohydrodynamic turbulence at Reynolds mumbers up to about 500, are performed using a uniform grid of 1803 collocation points. Strong vorticity and current sheets obtain both in the presence and in the absence of magnetic nulls. Contrary to vortex sheets in hydrodynamics, these structures do not destabilize into filaments, but are locally disrupted. They are the main loci of kinetic and magnetic dissipations. copyright 1995 American Institute of Physics
High Reynolds number magnetohydrodynamic turbulence using a Lagrangian model
Graham, J. Pietarila; P. D. Mininni; Pouquet, A.
2011-01-01
With the help of a model of magnetohydrodynamic (MHD) turbulence tested previously, we explore high Reynolds number regimes up to equivalent resolutions of 6000^3 grid points in the absence of forcing and with no imposed uniform magnetic field. For the given initial condition chosen here, with equal kinetic and magnetic energy, the flow ends up being dominated by the magnetic field, and the dynamics leads to an isotropic Iroshnikov-Kraichnan energy spectrum. However, the locally anisotropic m...
Energy Decay Laws in Strongly Anisotropic Magnetohydrodynamic Turbulence
We investigate the influence of a uniform magnetic field B0=B0eparallel on energy decay laws in incompressible magnetohydrodynamic (MHD) turbulence. The nonlinear transfer reduction along B0 is included in a model that distinguishes parallel and perpendicular directions, following a phenomenology of Kraichnan. We predict a slowing down of the energy decay due to anisotropy in the limit of strong B0, with distinct power laws for energy decay of shear- and pseudo-Alfven waves. Numerical results from the kinetic equations of Alfven wave turbulence recover these predictions, and MHD numerical results clearly tend to follow them in the lowest perpendicular planes
Magnetohydrodynamic turbulence and enhanced atomic processes in astrophysical plasmas
Spangler, Steven R.
1998-08-01
This article discusses a way in which enhanced atomic physics processes, including radiative energy losses, may occur in an astrophysical plasma containing magnetohydrodynamic turbulence. Two-dimensional (2D) magnetohydrodynamics (MHD) is adopted as a model. A major characteristic feature of 2D MHD turbulence is the development of strong current sheets on a dynamical time scale L/V0 where L is the spatial scale of the turbulent fluid and V0 is the scale of the velocity fluctuations. The current contained in the sheets will be carried by an electron drift relative to the ions. The case of a plasma containing minority atoms or ions with an excited state accessible to collisions from the tail of the electron distribution is considered. In the current carrying sheets or filaments, the electron distribution function will be perturbed such that collisional excitations will be enhanced relative to the current-free plasma. Subsequent radiative de-excitation of the atoms or ions removes energy from the turbulence. Expressions are presented for the electron drift velocity arising in 2D turbulence, the enhancement of collisional excitations of a trace atom or ion, and the energy lost to the plasma turbulence by radiative de-excitation of these atoms or ions. The mechanism would be most pronounced in plasmas for which the magnitude of the magnetic field is large, the outer scale of the turbulence is small, and the electron density and temperature are low. A brief discussion of the relevance of this mechanism to some specific astrophysical plasmas is given.
Minimal flow units for magnetohydrodynamic turbulence
Orlandi, P.
2016-08-01
We present direct numerical simulations of two minimal flow units (MFUs) to investigate the differences between inviscid and viscous simulations, and the different behavior of the evolution for conducting fluids. In these circumstances the introduction of the Lorentz force in the momentum equation produces different scenarios. The Taylor–Green vortex, in the past, was an MFU widely considered for both conducting and non-conducting fluids. The simulations were performed by pseudo-spectral numerical methods; these are repeated here by using a finite difference second-order accurate, energy-conserving scheme for ν =0. Having observed that this initial condition could be inefficient for capturing the eventual occurrence of a finite time singularity a potentially more efficient MFU consisting of two interacting Lamb dipoles was considered. It was found that the two flows have a different time evolution in the vortical dominated stage. In this stage, turbulent structures of different size are generated leading to spectra, in the inviscid conditions, with a {k}-3 range. In real conditions the viscosity produces smaller scales characteristic of fully developed turbulence with energy spectra with well defined exponential and inertial ranges. In the presence of non-conducting conditions the passive vector behaves as the vorticity. The evolution is different in the presence of conducting conditions. Although the time evolution is different, both flows lead to spectra in Kolmogorov units with the same shape at high and intermediate wave numbers.
Magnetohydrodynamic Turbulence and Reconnection in the Magnetotail
El-Alaoui, Mostafa; Richard, Robert L.; Ashour-Abdalla, Maha; Goldstein, Melvyn L.; Walker, Raymond J.
2011-01-01
We have used a global MHD simulation with high spatial resolution to investigate the origin and properties of turbulence in the plasma sheet. In this simulation we imposed a steady southward IMF with a magnitude of 5 nT at the upstream simulation boundary for more than three hours followed by ninety minutes of northward IMF of the same magnitude. The solar wind number density was 20 cm-3, the thermal pressure was 20 pPa, and the velocity was 500 km/s in the x direction. The moderately high dynamic pressure confined the magnetotail to the high-grid resolution region. Even for these nominal solar wind parameters and steady driving the plasma sheet became turbulent. The power spectral densities and probability distribution functions computed from the simulations were comparable to those obtained from spacecraft observations. The largest scale vortices were associated with reconnection outflows and, in the southward IMF case, with the diversion of high speed flows in the near-Earth region. Both time and space domain analyses revealed that there were three scales present, the large scale of the driving processes, the intermediate inertial scale and the dissipative scale.
Temporal Intermittency of Energy Dissipation in Magnetohydrodynamic Turbulence
Zhdankin, Vladimir; Boldyrev, Stanislav
2015-01-01
Energy dissipation in magnetohydrodynamic (MHD) turbulence is known to be highly intermittent in space, being concentrated in sheet-like coherent structures. Much less is known about intermittency in time, another fundamental aspect of turbulence which has great importance for observations of solar flares and other space/astrophysical phenomena. In this Letter, we investigate the temporal intermittency of energy dissipation in numerical simulations of MHD turbulence. We consider four-dimensional spatiotemporal structures, "flare events", responsible for a large fraction of the energy dissipation. We find that although the flare events are often highly complex, they exhibit robust power-law distributions and scaling relations. We find that the probability distribution of dissipated energy has a power law index close to -1.75, similar to observations of solar flares, indicating that intense dissipative events dominate the heating of the system. We also discuss the temporal asymmetry of flare events as a signatu...
Electromotive force due to magnetohydrodynamic fluctuations in sheared rotating turbulence
Squire, J.; Bhattacharjee, A.
2015-11-01
This article presents a calculation of the mean electromotive force arising from general small-scale magnetohydrodynamical turbulence, within the framework of the second-order correlation approximation. With the goal of improving understanding of the accretion disk dynamo, effects arising through small-scale magnetic fluctuations, velocity gradients, density and turbulence stratification, and rotation, are included. The primary result, which supplements numerical findings, is that an off-diagonal turbulent resistivity due to magnetic fluctuations can produce large-scale dynamo action—the magnetic analog of the "shear-current" effect. In addition, consideration of α effects in the stratified regions of disks gives the puzzling result that there is no strong prediction for a sign of α , since the effects due to kinetic and magnetic fluctuations, as well as those due to shear and rotation, are each of opposing signs and tend to cancel each other.
Scaling of compressible magnetohydrodynamic turbulence in the fast solar wind
Banerjee, Supratik; Sahraoui, Fouad; Galtier, Sebastien
2016-01-01
The role of compressible fluctuations in the energy cascade of fast solar wind turbulence is studied using a reduced form of an exact law derived recently (Banerjee and Galtier, PRE, 2013) for compressible isothermal magnetohydrodynamics and in-situ observations from the THEMIS B/ARTEMIS P1 spacecraft. A statistical survey of the data revealed a turbulent energy cascade over two decades of scales, which is broader than the previous estimates made from an exact incompressible law. A term-by-term analysis of the compressible model reveals new insight into the role played by the compressible fluctuations in the energy cascade. The compressible fluctuations are shown to amplify (2 to 4 times) the turbulent cascade rate with respect to the incompressible model in 10 % of the analyzed samples. This new estimated cascade rate is shown to provide the adequate energy dissipation required.
Electromotive force due to magnetohydrodynamic fluctuations in sheared rotating turbulence
Squire, J. [Dept Astrophys Sci, Max Planck Princeton Ctr Plasma Phys, Princeton, NJ 08543 USA; Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States); Bhattacharjee, A. [Dept Astrophys Sci, Max Planck Princeton Ctr Plasma Phys, Princeton, NJ 08543 USA; Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
2015-11-01
This article presents a calculation of the mean electromotive force arising from general small-scale magnetohydrodynamical turbulence, within the framework of the second-order correlation approximation. With the goal of improving understanding of the accretion disk dynamo, effects arising through small-scale magnetic fluctuations, velocity gradients, density and turbulence stratification, and rotation, are included. The primary result, which supplements numerical findings, is that an off-diagonal turbulent resistivity due to magnetic fluctuations can produce large-scale dynamo action-the magnetic analog of the "shear-current" effect. In addition, consideration of alpha effects in the stratified regions of disks gives the puzzling result that there is no strong prediction for a sign of alpha, since the effects due to kinetic and magnetic fluctuations, as well as those due to shear and rotation, are each of opposing signs and tend to cancel each other.
Hierarchical Structure of Magnetohydrodynamic Turbulence In Position-Position-Velocity Space
Burkhart, Blakesley; Goodman, Alyssa; Rosolowsky, Erik
2012-01-01
Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic Position-Position-Velocity (PPV) emission cubes of optically thin isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the physics of the gas, i.e. self-gravity and the global sonic and Alfvenic Mach number. Simulations with higher Alfvenic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observed a strong sonic and Alfvenic dependency when we apply the the statistical moments (i.e. mean, variance, skewness, kurtosis) to the dendrogram distribution. Larger magnetic field and sonic Mach number correspond to larger values of the moments. Application of the dendrogram to 3D density cubes, also known as Position-Position-Position cube...
Leprovost, Nicolas; Kim, Eun-Jin
2009-08-01
We investigate three-dimensional magnetohydrodynamics turbulence in the presence of velocity and magnetic shear (i.e., with both a large-scale shear flow and a nonuniform magnetic field). By assuming a turbulence driven by an external forcing with both helical and nonhelical spectra, we investigate the combined effect of these two shears on turbulence intensity and turbulent transport represented by turbulent diffusivities (turbulent viscosity, alpha and beta effect) in Reynolds-averaged equations. We show that turbulent transport (turbulent viscosity and diffusivity) is quenched by a strong flow shear and a strong magnetic field. For a weak flow shear, we further show that the magnetic shear increases the turbulence intensity while decreasing the turbulent transport. In the presence of a strong flow shear, the effect of the magnetic shear is found to oppose the effect of flow shear (which reduces turbulence due to shear stabilization) by enhancing turbulence and transport, thereby weakening the strong quenching by flow shear stabilization. In the case of a strong magnetic field (compared to flow shear), magnetic shear increases turbulence intensity and quenches turbulent transport. PMID:19792244
Broken Ergodicity in Two-Dimensional Homogeneous Magnetohydrodynamic Turbulence
Shebalin, John V.
2010-01-01
Two-dimensional (2-D) homogeneous magnetohydrodynamic (MHD) turbulence has many of the same qualitative features as three-dimensional (3-D) homogeneous MHD turbulence.The se features include several ideal invariants, along with the phenomenon of broken ergodicity. Broken ergodicity appears when certain modes act like random variables with mean values that are large compared to their standard deviations, indicating a coherent structure or dynamo.Recently, the origin of broken ergodicity in 3-D MHD turbulence that is manifest in the lowest wavenumbers was explained. Here, a detailed description of the origins of broken ergodicity in 2-D MHD turbulence is presented. It will be seen that broken ergodicity in ideal 2-D MHD turbulence can be manifest in the lowest wavenumbers of a finite numerical model for certain initial conditions or in the highest wavenumbers for another set of initial conditions.T he origins of broken ergodicity in ideal 2-D homogeneous MHD turbulence are found through an eigen analysis of the covariance matrices of the modal probability density functions.It will also be shown that when the lowest wavenumber magnetic field becomes quasi-stationary, the higher wavenumber modes can propagate as Alfven waves on these almost static large-scale magnetic structures
Spectral properties of decaying turbulence in electron magnetohydrodynamics
The spectral properties of decaying turbulence in 2(1/2)-dimensional electron magnetohydrodynamics are studied numerically. In the range kde-7/3 and k-13/3, respectively. The self-similar decay state of the turbulence is reached after an initial phase of fast exchange between the axial and poloidal magnetic energies. The time behavior t-2/3 of the total energy is found to be consistent with that obtained from selective decay. The maximum of the energy spectrum shifts towards low mode numbers and decays in time as t-1, in agreement with the infrared scaling of the turbulence. In the large de limit, both energy and mean square generalized momentum exhibit direct cascades. No stationary turbulent state could be found as long as the axial kinetic energy is large as compared to the poloidal kinetic energy initially. The global physical quantities decay well before turbulent macroscopic quantities have established similar space-time behavior, and the turbulence is infected by the lack of stationarity. The system decouples into a Navier-Stokes equation and a passive scalar equation only if the poloidal kinetic energy is larger than or equal to the axial kinetic energy. In this limit the k-5/3 and k-3 spectra of the poloidal kinetic energy are recovered
Cosmic Ray propagation in sub-Alfvenic magnetohydrodynamic turbulence
Cohet, Romain
2016-01-01
This work has the main objective to provide a detailed investigation of cosmic ray propagation in magnetohydrodynamic turbulent fields generated by forcing the fluid velocity field at large scales. It provides a derivation of the particle mean free path dependences in terms of the turbulence level described by the Alfv\\'enic Mach number and in terms of the particle rigidity. We use an upgrade version of the magnetohydrodynamic code {\\tt RAMSES} which includes a forcing module and a kinetic module and solve the Lorentz equation for each particle. The simulations are performed using a 3 dimension periodical box in the test-particle and magnetostatic limits. The forcing module is implemented using an Ornstein-Uhlenbeck process. An ensemble average over a large number of particle trajectories is applied to reconstruct the particle mean free paths. We derive the cosmic ray mean free paths in terms of the Alfv\\'enic Mach numbers and particle reduced rigidities in different turbulence forcing geometries. The reduced...
Energy transfers in shell models for magnetohydrodynamics turbulence.
Lessinnes, Thomas; Carati, Daniele; Verma, Mahendra K
2009-06-01
A systematic procedure to derive shell models for magnetohydrodynamic turbulence is proposed. It takes into account the conservation of ideal quadratic invariants such as the total energy, the cross helicity, and the magnetic helicity, as well as the conservation of the magnetic energy by the advection term in the induction equation. This approach also leads to simple expressions for the energy exchanges as well as to unambiguous definitions for the energy fluxes. When applied to the existing shell models with nonlinear interactions limited to the nearest-neighbor shells, this procedure reproduces well-known models but suggests a reinterpretation of the energy fluxes. PMID:19658594
Nonlinear closures for scale separation in supersonic magnetohydrodynamic turbulence
Grete, Philipp; Schmidt, Wolfram; Schleicher, Dominik R G; Federrath, Christoph
2015-01-01
Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment -- magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures -- models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The subgrid-scale energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pre...
The small-scale turbulent dynamo in smoothed particle magnetohydrodynamics
Tricco, Terrence S; Federrath, Christoph
2016-01-01
Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small-scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of star-forming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning met...
Dynamics of decaying two-dimensional magnetohydrodynamic turbulence
High-resolution numerical studies of decaying two-dimensional magnetohydrodynamic turbulence using up to 7682 modes in general periodic systems reveal the following properties: I) the evolution proceeds in a quasi-selfsimilar way with constant kinetic to magnetic energy ratio and constant micro- and macro-scale Reynolds numbers; II) the energy dissipation rate is independent of the values of the dissipation coefficients η, μ; III) the inertial-range energy spectra follow a Kolmogorov law, Ek=Cε2/3k-5/3 with C=3.7±0.3. Small-scale fluctuations are concentrated in the region of weak large-scale magnetic fields. The resulting strong intermittency is analogous to the behavior recently observed in two-dimensional hydrodynamic turbulence (see, for instance, J. Fluid Mech. (1988) 194, 333), with the magnetic field intensity taking the role vorticity plays in hydrodynamic systems. (orig.)
Non-elliptic wavevector anisotropy for magnetohydrodynamic turbulence
Narita, Y.
2015-11-01
A model of non-elliptic wavevector anisotropy is developed for the inertial-range spectrum of magnetohydrodynamic turbulence and is presented in the two-dimensional wavevector domain spanning the directions parallel and perpendicular to the mean magnetic field. The non-elliptic model is a variation of the elliptic model with different scalings along the parallel and the perpendicular components of the wavevectors to the mean magnetic field. The non-elliptic anisotropy model reproduces the smooth transition of the power-law spectra from an index of -2 in the parallel projection with respect to the mean magnetic field to an index of -5/3 in the perpendicular projection observed in solar wind turbulence, and is as competitive as the critical balance model to explain the measured frequency spectra in the solar wind. The parameters in the non-elliptic spectrum model are compared with the solar wind observations.
Magnetic reversals in a modified shell model for magnetohydrodynamics turbulence.
Nigro, Giuseppina; Carbone, Vincenzo
2010-07-01
The aim of the paper is the study of dynamo action using a simple nonlinear model in the framework of magnetohydrodynamic turbulence. The nonlinear behavior of the system is described by using a shell model for velocity field and magnetic field fluctuations, modified for the magnetic field at the largest scale by a term describing a supercritical pitchfork bifurcation. Turbulent fluctuations generate a dynamical situation where the large-scale magnetic field jumps between two states which represent the opposite polarities of the magnetic field. Despite its simplicity, the model has the capability to describe a long time series of reversals from which we infer results about the statistics of persistence times and scaling laws of cancellations between opposite polarities for different magnetic diffusivity coefficients. These properties of the model are compared with real paleomagnetic data, thus revealing the origin of long-range correlations in the process. PMID:20866731
Spontaneous chiral symmetry breaking of Hall magnetohydrodynamic turbulence.
Meyrand, Romain; Galtier, Sébastien
2012-11-01
Hall magnetohydrodynamics (MHD) is investigated through three-dimensional direct numerical simulations. We show that the Hall effect induces a spontaneous chiral symmetry breaking of the turbulent dynamics. The normalized magnetic polarization is introduced to separate the right- (R) and left-handed (L) fluctuations. A classical k(-7/3) spectrum is found at small scales for R magnetic fluctuations which corresponds to the electron MHD prediction. A spectrum compatible with k(-11/3) is obtained at large-scales for the L magnetic fluctuations; we call this regime the ion MHD. These results are explained heuristically by rewriting the Hall MHD equations in a succinct vortex dynamical form. Applications to solar wind turbulence are discussed. PMID:23215387
Magnetic reversals in a modified shell model for magnetohydrodynamics turbulence
Nigro, Giuseppina; Carbone, Vincenzo
2010-07-01
The aim of the paper is the study of dynamo action using a simple nonlinear model in the framework of magnetohydrodynamic turbulence. The nonlinear behavior of the system is described by using a shell model for velocity field and magnetic field fluctuations, modified for the magnetic field at the largest scale by a term describing a supercritical pitchfork bifurcation. Turbulent fluctuations generate a dynamical situation where the large-scale magnetic field jumps between two states which represent the opposite polarities of the magnetic field. Despite its simplicity, the model has the capability to describe a long time series of reversals from which we infer results about the statistics of persistence times and scaling laws of cancellations between opposite polarities for different magnetic diffusivity coefficients. These properties of the model are compared with real paleomagnetic data, thus revealing the origin of long-range correlations in the process.
Helical mode interactions and spectral transfer processes in magnetohydrodynamic turbulence
Linkmann, Moritz F; McKay, Mairi E; Jäger, Julia
2015-01-01
Spectral transfer processes in magnetohydrodynamic (MHD) turbulence are investigated analytically by decomposition of the velocity and magnetic fields in Fourier space into helical modes. Steady solutions of the dynamical system which governs the evolution of the helical modes are determined, and a stability analysis of these solutions is carried out. The interpretation of the analysis is that unstable solutions lead to energy transfer between the interacting modes while stable solutions do not. From this, a dependence of possible interscale energy and helicity transfers on the helicities of the interacting modes is derived. As expected from the inverse cascade of magnetic helicity in 3D MHD turbulence, mode interactions with like helicities lead to transfer of energy and magnetic helicity to smaller wavenumbers. However, some interactions of modes with unlike helicities also contribute to an inverse energy transfer. As such, an inverse energy cascade for nonhelical magnetic fields is shown to be possible. Fu...
Direct Evidence of the Transition from Weak to Strong Magnetohydrodynamic Turbulence.
Meyrand, Romain; Galtier, Sébastien; Kiyani, Khurom H
2016-03-11
One of the most important predictions in magnetohydrodynamics is that in the presence of a uniform magnetic field b_{0}e[over ^]_{∥} a transition from weak to strong wave turbulence should occur when going from large to small perpendicular scales. This transition is believed to be a universal property of several anisotropic turbulent systems. We present, for the first time, direct evidence of such a transition using a decaying three-dimensional direct numerical simulation of incompressible balanced magnetohydrodynamic turbulence with a grid resolution of 3072^{2}×256. From large to small scales, the change of regime is characterized by (i) a change of slope in the energy spectrum going from approximately -2 to -3/2, (ii) an increase of the ratio between the wave and nonlinear times, with a critical ratio of χ_{c}∼1/3, (iii) a modification of the isocontours of energy revealing a transition from a purely perpendicular cascade to a cascade compatible with the critical-balance-type phenomenology, and (iv) an absence followed by a dramatic increase of the communication between Alfvén modes. The changes happen at approximately the same transition scale and can be seen as manifest signatures of the transition from weak to strong wave turbulence. Furthermore, we observe a significant nonlocal three-wave coupling between strongly and weakly nonlinear modes resulting in an inverse transfer of energy from small to large scales. PMID:27015486
Direct Evidence of the Transition from Weak to Strong Magnetohydrodynamic Turbulence
Meyrand, Romain; Galtier, Sébastien; Kiyani, Khurom H.
2016-03-01
One of the most important predictions in magnetohydrodynamics is that in the presence of a uniform magnetic field b0e^∥ a transition from weak to strong wave turbulence should occur when going from large to small perpendicular scales. This transition is believed to be a universal property of several anisotropic turbulent systems. We present, for the first time, direct evidence of such a transition using a decaying three-dimensional direct numerical simulation of incompressible balanced magnetohydrodynamic turbulence with a grid resolution of 30722×256 . From large to small scales, the change of regime is characterized by (i) a change of slope in the energy spectrum going from approximately -2 to -3 /2 , (ii) an increase of the ratio between the wave and nonlinear times, with a critical ratio of χc˜1 /3 , (iii) a modification of the isocontours of energy revealing a transition from a purely perpendicular cascade to a cascade compatible with the critical-balance-type phenomenology, and (iv) an absence followed by a dramatic increase of the communication between Alfvén modes. The changes happen at approximately the same transition scale and can be seen as manifest signatures of the transition from weak to strong wave turbulence. Furthermore, we observe a significant nonlocal three-wave coupling between strongly and weakly nonlinear modes resulting in an inverse transfer of energy from small to large scales.
Long-term evolution of decaying magnetohydrodynamic turbulence in the multiphase interstellar medium
Supersonic turbulence in the interstellar medium (ISM) is believed to decay rapidly within a flow crossing time irrespective of the degree of magnetization. However, this general consensus of decaying magnetohydrodynamic (MHD) turbulence relies on local isothermal simulations, which are unable to take into account the roles of the global structures of magnetic fields and the ISM. Utilizing three-dimensional MHD simulations including interstellar cooling and heating, we investigate decaying MHD turbulence within cold neutral medium sheets embedded in a warm neutral medium. The early evolution of turbulent kinetic energy is consistent with previous results for decaying compressible MHD turbulence characterized by rapid energy decay with a power-law form of E∝t –1 and by a short decay time compared with the flow crossing time. If initial magnetic fields are strong and perpendicular to the sheet, however, long-term evolution of the kinetic energy shows that a significant amount of turbulent energy (∼0.2E 0) still remains even after 10 flow crossing times for models with periodic boundary conditions. The decay rate is also greatly reduced as the field strength increases for such initial and boundary conditions, but not if the boundary conditions are those for a completely isolated sheet. We analyze velocity power spectra of the remaining turbulence to show that in-plane, incompressible motions parallel to the sheet dominate at later times.
RESONANCE BROADENING AND HEATING OF CHARGED PARTICLES IN MAGNETOHYDRODYNAMIC TURBULENCE
The heating, acceleration, and pitch-angle scattering of charged particles by magnetohydrodynamic (MHD) turbulence are important in a wide range of astrophysical environments, including the solar wind, accreting black holes, and galaxy clusters. We simulate the interaction of high-gyrofrequency test particles with fully dynamical simulations of subsonic MHD turbulence, focusing on the parameter regime with β ∼ 1, where β is the ratio of gas to magnetic pressure. We use the simulation results to calibrate analytical expressions for test particle velocity-space diffusion coefficients and provide simple fits that can be used in other work. The test particle velocity diffusion in our simulations is due to a combination of two processes: interactions between particles and magnetic compressions in the turbulence (as in linear transit-time damping; TTD) and what we refer to as Fermi Type-B (FTB) interactions, in which charged particles moving on field lines may be thought of as beads sliding along moving wires. We show that test particle heating rates are consistent with a TTD resonance that is broadened according to a decorrelation prescription that is Gaussian in time (but inconsistent with Lorentzian broadening due to an exponential decorrelation function, a prescription widely used in the literature). TTD dominates the heating for vs >> vA (e.g., electrons), where vs is the thermal speed of species s and vA is the Alfvén speed, while FTB dominates for vs A (e.g., minor ions). Proton heating rates for β ∼ 1 are comparable to the turbulent cascade rate. Finally, we show that velocity diffusion of collisionless, large gyrofrequency particles due to large-scale MHD turbulence does not produce a power-law distribution function.
Small-scale behavior of Hall magnetohydrodynamic turbulence
Stawarz, Julia E.; Pouquet, Annick
2015-12-01
Decaying Hall magnetohydrodynamic (HMHD) turbulence is studied using three-dimensional (3D) direct numerical simulations with grids up to 7683 points and two different types of initial conditions. Results are compared to analogous magnetohydrodynamic (MHD) runs and both Laplacian and Laplacian-squared dissipative operators are examined. At scales below the ion inertial length, the ratio of magnetic to kinetic energy as a function of wave number transitions to a magnetically dominated state. The transition in behavior is associated with the advection term in the momentum equation becoming subdominant to dissipation. Examination of autocorrelation functions reveals that, while current and vorticity structures are similarly sized in MHD, HMHD current structures are narrower and vorticity structures are wider. The electric field autocorrelation function is significantly narrower in HMHD than in MHD and is similar to the HMHD current autocorrelation function at small separations. HMHD current structures are found to be significantly more intense than in MHD and appear to have an enhanced association with strong alignment between the current and magnetic field, which may be important in collisionless plasmas where field-aligned currents can be unstable. When hyperdiffusivity is used, a longer region consistent with a k-7 /3 scaling is present for right-polarized fluctuations when compared to Laplacian dissipation runs.
Protostellar jets and magnetised turbulence with smoothed particle magnetohydrodynamics
Tricco, Terrence
2016-01-01
Magnetic fields are an integral component of the formation of stars. During my thesis work, I built new methods to model magnetic fields in smoothed particle magnetohydrodynamics which enforce the divergence-free constraint on the magnetic field and reduce numerical dissipation of the magnetic field. Using these methods, we have performed simulations of isolated protostar formation, studying the production of jets and outflows of material and their effect on transporting angular momentum away from the protostar and reducing the efficiency of star formation. A major code comparison project on the small-scale turbulent dynamo amplification of magnetic fields was performed, using conditions representative of molecular clouds, the formation site of stars. The results were compared against results from grid-based methods, finding excellent agreement on their statistics and qualitative behaviour. I will outline the numerical methods developed, and present the results from our protostar and molecular cloud simulations.
Turbulent energy dissipation and intermittency in ambipolar diffusion magnetohydrodynamics
Momferratos, Georgios; Falgarone, Edith; Forêts, Guillaume Pineau des
2015-01-01
The dissipation of kinetic and magnetic energy in the interstellar medium (ISM) can proceed through viscous, Ohmic or ambipolar diffusion (AD). It occurs at very small scales compared to the scales at which energy is presumed to be injected. This localized heating may impact the ISM evolution but also its chemistry, thus providing observable features. Here, we perform 3D spectral simulations of decaying magnetohydrodynamic turbulence including the effects of AD. We find that the AD heating power spectrum peaks at scales in the inertial range, due to a strong alignment of the magnetic and current vectors in the dissipative range. AD affects much greater scales than the AD scale predicted by dimensional analysis. We find that energy dissipation is highly concentrated on thin sheets. Its probability density function follows a lognormal law with a power-law tail which hints at intermittency, a property which we quantify by use of structure function exponents. Finally, we extract structures of high dissipation, de...
Magnetic moment non-conservation in magnetohydrodynamic turbulence models
Dalena, S; Rappazzo, A F; Mace, R L; Matthaeus, W H
2012-01-01
The fundamental assumptions of the adiabatic theory do not apply in presence of sharp field gradients as well as in presence of well developed magnetohydrodynamic turbulence. For this reason in such conditions the magnetic moment $\\mu$ is no longer expected to be constant. This can influence particle acceleration and have considerable implications in many astrophysical problems. Starting with the resonant interaction between ions and a single parallel propagating electromagnetic wave, we derive expressions for the magnetic moment trapping width $\\Delta \\mu$ (defined as the half peak-to-peak difference in the particle magnetic moment) and the bounce frequency $\\omega_b$. We perform test-particle simulations to investigate magnetic moment behavior when resonances overlapping occurs and during the interaction of a ring-beam particle distribution with a broad-band slab spectrum. We find that magnetic moment dynamics is strictly related to pitch angle $\\alpha$ for a low level of magnetic fluctuation, $\\delta B/B_0...
Chiral Exact Relations for Helicities in Hall Magnetohydrodynamic Turbulence
Banerjee, Supratik
2016-01-01
Besides total energy, three-dimensional incompressible Hall magnetohydrodynamics (MHD) possesses two inviscid invariants which are the magnetic helicity and the generalized helicity. New exact relations are derived for homogeneous (non-isotropic) stationary Hall MHD turbulence (and also for its inertialess electron MHD limit) with non-zero helicities and in the asymptotic limit of large Reynolds numbers. The universal laws are written only in terms of mixed second-order structure functions, i.e. the scalar product of two different increments. It provides, therefore, a direct measurement of the dissipation rates for the corresponding invariant flux. This study shows that the generalized helicity cascade is strongly linked to the left polarized fluctuations while the magnetic helicity cascade is linked to the right polarized fluctuations.
Energy cascade and its locality in compressible magnetohydrodynamic turbulence
Yang, Yan; Shi, Yipeng; Wan, Minping; Matthaeus, William H.; Chen, Shiyi
2016-06-01
We investigate energy transfer across scales in three-dimensional compressible magnetohydrodynamic (MHD) turbulence, a model often used to study space and astrophysical plasmas. Analysis shows that kinetic and magnetic energies cascade conservatively from large to small scales in cases with varying degrees of compression. With more compression, energy fluxes due to pressure dilation and subscale mass flux are greater, but conversion between kinetic and magnetic energy by magnetic line stretching is less efficient. Energy transfer between the same fields is dominated by local contributions regardless of compressive effects. In contrast, the conversion between kinetic and internal energy by pressure dilation is dominated by the largest scale contributions. Energy conversion between the velocity and magnetic fields is weakly local.
Sub-grid-scale description of turbulent magnetic reconnection in magnetohydrodynamics
Widmer, F.; Büchner, J.; Yokoi, N.
2016-04-01
Magnetic reconnection requires, at least locally, a non-ideal plasma response. In collisionless space and astrophysical plasmas, turbulence could transport energy from large to small scales where binary particle collisions are rare. We have investigated the influence of small scale magnetohydrodynamics (MHD) turbulence on the reconnection rate in the framework of a compressible MHD approach including sub-grid-scale (SGS) turbulence. For this sake, we considered Harris-type and force-free current sheets with finite guide magnetic fields directed out of the reconnection plane. The goal is to find out whether unresolved by conventional simulations MHD turbulence can enhance the reconnection process in high-Reynolds-number astrophysical plasmas. Together with the MHD equations, we solve evolution equations for the SGS energy and cross-helicity due to turbulence according to a Reynolds-averaged turbulence model. The SGS turbulence is self-generated and -sustained through the inhomogeneities of the mean fields. By this way, the feedback of the unresolved turbulence into the MHD reconnection process is taken into account. It is shown that the turbulence controls the regimes of reconnection by its characteristic timescale τt. The dependence on resistivity was investigated for large-Reynolds-number plasmas for Harris-type as well as force-free current sheets with guide field. We found that magnetic reconnection depends on the relation between the molecular and apparent effective turbulent resistivity. We found that the turbulence timescale τt decides whether fast reconnection takes place or whether the stored energy is just diffused away to small scale turbulence. If the amount of energy transferred from large to small scales is enhanced, fast reconnection can take place. Energy spectra allowed us to characterize the different regimes of reconnection. It was found that reconnection is even faster for larger Reynolds numbers controlled by the molecular resistivity η, as
Turbulent Magnetohydrodynamic Reconnection Mediated by the Plasmoid Instability
Huang, Yi-Min
2015-01-01
It has been established that the Sweet-Parker current layer in high Lundquist number reconnection is unstable to the super-Alfv\\'enic plasmoid instability. Past two-dimensional magnetohydrodynamic simulations have demonstrated that the plasmoid instability leads to a new regime where the Sweet-Parker current layer changes into a chain of plasmoids connected by secondary current sheets, and the averaged reconnection rate becomes nearly independent of the Lundquist number. In this work, three-dimensional simulation with a guide field shows that the additional degree of freedom allows plasmoid instabilities to grow at oblique angles, which interact and lead to self-generated turbulent reconnection. The averaged reconnection rate in the self-generated turbulent state is of the order of a hundredth of the characteristic Alfv\\'en speed, which is similar to the two-dimensional result but is an order of magnitude lower than the fastest reconnection rate reported in recent studies of externally driven three-dimensiona...
The small-scale turbulent dynamo in smoothed particle magnetohydrodynamics
Tricco, T. S.; Price, D. J.; Federrath, C.
2016-05-01
Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small- scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of starforming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning method may be used to achieve ∇ • B = 0 to machine precision, albeit at significant computational expense.
Particle acceleration by turbulent magnetohydro-dynamic reconnection
Matthaeus, W. H.; Ambrosiano, J. J.; Goldstein, M. L.
1984-01-01
Test particles in a two dimensional, turbulent MHD simulation are found to undergo significant acceleration. The magnetic field configuration is a periodic sheet pinch which undergoes reconnection. The test particles are trapped in the reconnection region for times of order an Alfven transit time in the large electric fields that characterize the turbulent reconnection process at the relatively large magnetic Reynolds number used in the simulation. The maximum speed attained by these particles is consistent with an analytic estimate which depends on the reconnection electric field, the Alfven speed, and the ratio of Larmor period to the Alfven transit time.
High Reynolds number magnetohydrodynamic turbulence using a Lagrangian model.
Graham, J Pietarila; Mininni, P D; Pouquet, A
2011-07-01
With the help of a model of magnetohydrodynamic (MHD) turbulence tested previously, we explore high Reynolds number regimes up to equivalent resolutions of 6000(3) grid points in the absence of forcing and with no imposed uniform magnetic field. For the given initial condition chosen here, with equal kinetic and magnetic energy, the flow ends up being dominated by the magnetic field, and the dynamics leads to an isotropic Iroshnikov-Kraichnan energy spectrum. However, the locally anisotropic magnetic field fluctuations perpendicular to the local mean field follow a Kolmogorov law. We find that the ratio of the eddy turnover time to the Alfvén time increases with wave number, contrary to the so-called critical balance hypothesis. Residual energy and helicity spectra are also considered; the role played by the conservation of magnetic helicity is studied, and scaling laws are found for the magnetic helicity and residual helicity spectra. We put these results in the context of the dynamics of a globally isotropic MHD flow that is locally anisotropic because of the influence of the strong large-scale magnetic field, leading to a partial equilibration between kinetic and magnetic modes for the energy and the helicity. PMID:21867311
HIERARCHICAL STRUCTURE OF MAGNETOHYDRODYNAMIC TURBULENCE IN POSITION-POSITION-VELOCITY SPACE
Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium (ISM) that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic position-position-velocity (PPV) emission cubes of isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the presence of self-gravity and the global sonic and Alfvénic Mach numbers. Simulations with higher Alfvénic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observe a strong dependency on the sonic and Alfvénic Mach numbers and self-gravity when we apply the statistical moments (i.e., mean, variance, skewness, kurtosis) to the leaf and node distribution of the dendrogram. Simulations with self-gravity, larger magnetic field and higher sonic Mach number have dendrogram distributions with higher statistical moments. Application of the dendrogram to three-dimensional density cubes, also known as position-position-position (PPP) cubes, reveals that the dominant emission contours in PPP and PPV are related for supersonic gas but not for subsonic. We also explore the effects of smoothing, thermal broadening, and velocity resolution on the dendrograms in order to make our study more applicable to observational data. These results all point to hierarchical tree diagrams as being a promising additional tool for studying ISM turbulence and star forming regions for obtaining information on the degree of self-gravity, the Mach numbers and the complicated relationship between PPV and PPP data.
Depletion of Nonlinearity in Magnetohydrodynamic Turbulence: Insights from Analysis and Simulations
Gibbon, J; Krstulovic, G; Pandit, R; Politano, H; Ponty, Y; Pouquet, A; Sahoo, G; Stawarz, J
2015-01-01
We build on recent developments in the study of fluid turbulence [Gibbon \\textit{et al.} Nonlinearity 27, 2605 (2014)] to define suitably scaled, order-$m$ moments, $D_m^{\\pm}$, of $\\omega^\\pm= \\omega \\pm j$, where $\\omega$ and $j$ are, respectively, the vorticity and current density in three-dimensional magnetohydrodynamics (MHD). We show by mathematical analysis, for unit magnetic Prandtl number $P_M$, how these moments can be used to identify three possible regimes for solutions of the MHD equations; these regimes are specified by inequalities for $D_m^{\\pm}$ and $D_1^{\\pm}$. We then compare our mathematical results with those from our direct numerical simulations (DNSs) and thus demonstrate that 3D MHD turbulence is like its fluid-turbulence counterpart insofar as all solutions, which we have investigated, remain in \\textit{only one of these regimes}; this regime has depleted nonlinearity. We examine the implications of our results for the exponents $q^{\\pm}$ that characterize the power-law dependences of...
Zhdankin, Vladimir
2015-11-01
Energy dissipation is highly intermittent in large-scale turbulent plasmas, being localized in space and in time. This intermittency is manifest by the presence of coherent structures such as current (and vorticity) sheets, which account for a large fraction of the overall energy dissipation and may serve as sites for magnetic reconnection and particle acceleration. The statistical analysis of these dissipative structures is a robust and informative methodology for probing the underlying dynamics, both in numerical simulations and in observations. In this talk, the statistical properties of current sheets in numerical simulations of driven magnetohydrodynamic (MHD) turbulence are described, including recent results obtained from applying new methods for characterizing their morphology. Instantaneously, the overall energy dissipation is found to be evenly spread among current sheets spanning a continuum of energy dissipation rates and inertial-range sizes, while their thicknesses are localized deep inside the dissipation range. The temporal dynamics are then investigated by tracking the current sheets in time and considering the statistics of the resulting four-dimensional spatiotemporal structures, which correspond to dissipative events or flares in astrophysical systems. These dissipative events are found to exhibit robust power-law distributions and scaling relations, and are often highly complex, long-lived, and weakly asymmetric in time. Based on the distribution for their dissipated energies, the strongest dissipative events are found to dominate the overall energy dissipation in the system. These results are compared to the observed statistics of solar flares, and some possible implications for the solar wind are also described.
Plasmoid Instabilities Mediated Three-Dimensional Magnetohydrodynamic Turbulent Reconnection
Huang, Yi-min [Princeton University; Guo, Fan [Los Alamos National Laboratory
2015-07-21
After some introductory remarks on fast reconnection in resistive MHD due to plasmoid instability, oblique tearing modes in 3D, and previous studies on 3D turbulent reconnection, the subject is presented under the following topics: 3D simulation setup, time evolution of the 3D simulation, comparison with Sweet-Parker and 2D plasmoid reconnection, and diagnostics of the turbulent state (decomposition of mean fields and fluctuations, power spectra of energy fluctuations, structure function and eddy anisotropy with respect to local magnetic field). Three primary conclusions were reached: (1) The results suggest that 3D plasmoid instabilities can lead to self-generated turbulent reconnection (evidence of energy cascade and development of inertial range, energy fluctuations preferentially align with the local magnetic field, which is one of the characteristics of MHD turbulence); (2) The turbulence is highly inhomogeneous, due to the presence of magnetic shear and outflow jets (conventional MHD turbulence theories or phenomenologies may not be applicable – e.g. scale-dependent anisotropy as predicted by Goldreich & Sridhar is not found); (3) 3D turbulent reconnection is different from 2D plasmoid-dominated reconnection in many aspects. However, in fully developed state, reconnection rates in 2D and 3D are comparable — this result needs to be further checked in higher S.
Cluster observation of magnetohydrodynamic turbulence in the plasma sheet boundary layer
Narita, Y.
2016-04-01
Measurement of turbulent magnetic field is presented from the Earth magnetotail crossing of the Cluster spacecraft on August 25, 2006, as an ideal case study of magnetohydrodynamic turbulence in the plasma sheet boundary layer on a spatial scale of about 10,000 km. The fluctuation energy of the magnetic field is evaluated in both the frequency and wavevector domains. The observed plasma sheet turbulence event shows anisotropy in the wavevector domain with a spectral extension perpendicular to the mean magnetic field. The analyses of the dispersion relation and phase speed diagrams indicate that the coherent wave components should be regarded as a set of the linear-mode waves and the other fluctuation components in magnetohydrodynamics. Although the magnetic field fluctuation amplitudes are sufficiently small compared to the large-scale field strength, there is no clear indication of the linear-mode dominance in the plasma sheet. As a lesson, magnetohydrodynamic turbulence must be modeled by including both linear-mode waves and nonlinear wave components such as sideband waves.
We have developed a three-fluid, three-dimensional magnetohydrodynamic solar wind model that incorporates turbulence transport, eddy viscosity, turbulent resistivity, and turbulent heating. The solar wind plasma is described as a system of co-moving solar wind protons, electrons, and interstellar pickup protons, with separate energy equations for each species. Numerical steady-state solutions of Reynolds-averaged solar wind equations coupled with turbulence transport equations for turbulence energy, cross helicity, and correlation length are obtained by the time relaxation method in the corotating with the Sun frame of reference in the region from 0.3 to 100 AU (but still inside the termination shock). The model equations include the effects of electron heat conduction, Coulomb collisions, photoionization of interstellar hydrogen atoms and their charge exchange with the solar wind protons, turbulence energy generation by pickup protons, and turbulent heating of solar wind protons and electrons. The turbulence transport model is based on the Reynolds decomposition and turbulence phenomenologies that describe the conversion of fluctuation energy into heat due to a turbulent cascade. In addition to using separate energy equations for the solar wind protons and electrons, a significant improvement over our previous work is that the turbulence model now uses an eddy viscosity approximation for the Reynolds stress tensor and the mean turbulent electric field. The approximation allows the turbulence model to account for driving of turbulence by large-scale velocity gradients. Using either a dipole approximation for the solar magnetic field or synoptic solar magnetograms from the Wilcox Solar Observatory for assigning boundary conditions at the coronal base, we apply the model to study the global structure of the solar wind and its three-dimensional properties, including embedded turbulence, heating, and acceleration throughout the heliosphere. The model results are
Nearly incompressible fluids. II - Magnetohydrodynamics, turbulence, and waves
Zank, G. P.; Matthaeus, W. H.
1993-01-01
The theory of nearly incompressible (NI) fluid dynamics developed previously for hydrodynamics is extended to magnetohydrodynamics (MHD). Based on a singular expansion technique, modified systems of fluid equations are obtained for which the effects of compressibility are admitted only weakly in terms of the different possible incompressible solutions. NI MHD represents the interface between the compressible and incompressible magnetofluid descriptions in the subsonic regime. It is shown that three distinct NI descriptions exist corresponding to each of the three possible plasma beta regimes. The detailed theory of weakly compressible corrections to the various incompressible MHD descriptions is presented, and the implications for the solar wind are discussed.
Gomez, T; Sagaut, P; Schilling, O; Zhou, Y
2006-07-05
A spectral subggrid-scale eddy viscosity and magnetic resisitivity model based on the eddy-damped quasi-normal Markovian (EDQNM) spectral kinetic and magnetic energy transfer presented in [12] is used in large-eddy simulation (LES) of large kinetic and magnetic Reynold number magneto-hydrodynamic (MHD) turbulence. The proposed model is assessed via a posteri tests on three-dimensional, incompressible, isotropic, non-helical, freely-decaying MHD turbulence at asymptotically large Reynolds numbers. Using LES with an initial condition characterized by an Alfv{acute e}n ratio of kinetic to magnetic energy {tau}{sub A} equal to unity, it is shown that at the kinetic energy spectrum E{sub K}(k) and magnetic energy spectrum E{sub M}(k) exhibit Kolmogorov -5/3 inertial subrange scalings in the LES, consistent with the EDQNM model.
Residual energy in magnetohydrodynamic turbulence and in the solar wind
Boldyrev, Stanislav; Zhdankin, Vladimir
2011-01-01
Recent observations indicate that kinetic and magnetic energies are not in equipartition in the solar wind turbulence. Rather, magnetic fluctuations are more energetic and have somewhat steeper energy spectrum compared to the velocity fluctuations. This leads to the presence of the so-called residual energy E_r=E_v-E_b in the inertial interval of turbulence. This puzzling effect is addressed in the present paper in the framework of weak turbulence theory. Using a simple model of weakly colliding Alfv\\'en waves, we demonstrate that the kinetic-magnetic equipartition indeed gets broken as a result of nonlinear interaction of Alfv\\'en waves. We establish that magnetic energy is indeed generated more efficiently as a result of these interactions, which proposes an explanation for the solar wind observations.
On the compressibility effect in test particle acceleration by magnetohydrodynamic turbulence
C.A. González; Dmitruk, P.; Mininni, P.D.; Matthaeus, W. H.
2016-01-01
The effect of compressibility in charged particle energization by magnetohydrodynamic (MHD) fields is studied in the context of test particle simulations. This problem is relevant to the solar wind and the solar corona due to the compressible nature of the flow in those astrophysical scenarios. We consider turbulent electromagnetic fields obtained from direct numerical simulations of the MHD equations with a strong background magnetic field. In order to explore the compressibilty effect over ...
Pratt, J.; Busse, A.; Müller, W. -C.
2013-01-01
Intermittent large-scale high-shear flows are found to occur frequently and spontaneously in direct numerical simulations of statistically stationary turbulent Boussinesq magnetohydrodynamic (MHD) convection. The energetic steady state of the system is sustained by convective driving of the velocity field and small-scale dynamo action. The intermittent emergence of flow structures with strong velocity and magnetic shearing generates magnetic energy at an elevated rate on time scales ...
Large-eddy simulations of fluid and magnetohydrodynamic turbulence using renormalized parameters
Mahendra K Verma; Shishir Kumar
2004-09-01
In this paper a procedure for large-eddy simulation (LES) has been devised for fluid and magnetohydrodynamic turbulence in Fourier space using the renormalized parameters; The parameters calculated using field theory have been taken from recent papers by Verma [1, 2]. We have carried out LES on 643 grid. These results match quite well with direct numerical simulations of 1283. We show that proper choice of parameter is necessary in LES.
Fluctuation dynamo driven by shear-bursts in convectively-driven magnetohydrodynamic turbulence
Pratt, J; Mueller, W -C
2013-01-01
Intermittent large-scale high-shear flows are found to occur frequently and spontaneously in direct numerical simulations of statistically stationary turbulent Boussinesq magnetohydrodynamic (MHD) convection. The energetic steady-state of the system is sustained by convective driving of the velocity field and small-scale dynamo action. The intermittent emergence of flow structures with strong velocity and magnetic shearing generates magnetic energy at an elevated rate over time-scales longer than the characteristic time of the large-scale convective motion. The resilience of magnetic energy amplification suggests that intermittent shear-bursts are a significant driver of dynamo action in turbulent magnetoconvection.
A numerical study of the alpha model for two-dimensional magnetohydrodynamic turbulent flows
Mininni, P D; Pouquet, A G
2004-01-01
We explore some consequences of the ``alpha model,'' also called the ``Lagrangian-averaged'' model, for two-dimensional incompressible magnetohydrodynamic (MHD) turbulence. This model is an extension of the smoothing procedure in fluid dynamics which filters velocity fields locally while leaving their associated vorticities unsmoothed, and has proved useful for high Reynolds number turbulence computations. We consider several known effects (selective decay, dynamic alignment, inverse cascades, and the probability distribution functions of fluctuating turbulent quantities) in magnetofluid turbulence and compare the results of numerical solutions of the primitive MHD equations with their alpha-model counterparts' performance for the same flows, in regimes where available resolution is adequate to explore both. The hope is to justify the use of the alpha model in regimes that lie outside currently available resolution, as will be the case in particular in three-dimensional geometry or for magnetic Prandtl number...
Arthur, S J; Mellema, G; De Colle, F; Vázquez-Semadeni, E
2011-01-01
We present the results of radiation-magnetohydrodynamic simulations of the formation and expansion of HII regions and their surrounding photodissociation regions in turbulent, magnetised, molecular clouds on scales of up to 4 parsecs. We include the effects of ionising and non-ionising ultraviolet radiation and x rays from young star clusters. We find that the HII region expansion reduces the disordered component of the B field, imposing a large-scale order on the field around its border, with the field in the neutral gas tending to lie along the ionisation front, while the field in the ionised gas tends to be perpendicular to the front. The highest pressure compressed neutral and molecular gas is driven towards approximate equipartition between thermal/magnetic/turbulent energy densities, whereas lower pressure neutral/molecular gas divides into quiescent, magnetically dominated regions, and, on the other hand, turbulent, demagnetised regions. The ionised gas shows approximate thermal/turbulent equipartition...
Isichenko, M B
1994-01-01
The long-time relaxation of ideal two dimensional magnetohydrodynamic turbulence subject to the conservation of two infinite families of constants of motion---the magnetic and the "cross" topology invariants--is examined. The analysis of the Gibbs ensemble, where all integrals of motion are respected, predicts the initial state to evolve into an equilibrium, stable coherent structure (the most probable state) and decaying Gaussian turbulence (fluctuations) with a vanishing, but always positive temperature. The non-dissipative turbulence decay is accompanied by decrease in both the amplitude and the length scale of the fluctuations, so that the fluctuation energy remains finite. The coherent structure represents a set of singular magnetic islands with plasma flow whose magnetic topology is identical to that of the initial state, while the energy and the cross topology invariants are shared between the coherent structure and the Gaussian turbulence. These conservation laws suggest the variational principle of i...
Ghosh, Sanjoy [Johns Hopkins University Applied Physics Laboratory, Laurel, Maryland 20723 (United States); Parashar, Tulasi N. [University of Delaware, Newark, Delaware 19716 (United States)
2015-04-15
The local k-space ratio of linear and nonlinear accelerations associated with a variety of initial conditions undergoing steady relaxation is investigated for the Hall–finite-Larmor-radius magnetohydrodynamics (MHD) system in the presence of a mean magnetic field. Building on a related study (Paper I) where it was shown that discrepancies exist between describing the global and local characterizations of the pure MHD system with mean magnetic field, we find regions of the Fourier space that are consistently dominated by linear acceleration and other regions that are consistently dominated by nonlinear acceleration, independent of the overall system's description as linear, weakly nonlinear, or turbulent. In general, dynamics within a certain angular range of the mean magnetic field direction are predominantly linear, while dynamics adjacent the Hall scales along the field-parallel direction and dynamics adjacent the finite Larmor radius scales in the field-perpendicular direction can become strongly nonlinear. The nonlinear influences are particularly significant as the plasma beta increases from unity to higher values.
The local k-space ratio of linear and nonlinear accelerations associated with a variety of initial conditions undergoing steady relaxation is investigated for the Hall–finite-Larmor-radius magnetohydrodynamics (MHD) system in the presence of a mean magnetic field. Building on a related study (Paper I) where it was shown that discrepancies exist between describing the global and local characterizations of the pure MHD system with mean magnetic field, we find regions of the Fourier space that are consistently dominated by linear acceleration and other regions that are consistently dominated by nonlinear acceleration, independent of the overall system's description as linear, weakly nonlinear, or turbulent. In general, dynamics within a certain angular range of the mean magnetic field direction are predominantly linear, while dynamics adjacent the Hall scales along the field-parallel direction and dynamics adjacent the finite Larmor radius scales in the field-perpendicular direction can become strongly nonlinear. The nonlinear influences are particularly significant as the plasma beta increases from unity to higher values
Quasi-static magnetohydrodynamic turbulence at high Reynolds number
We analyse the anisotropy of turbulence in an electrically conducting fluid submitted to a uniform magnetic field, for low magnetic Reynolds number, using the quasi-static approximation. In the linear limit, the kinetic energy of velocity components normal to the magnetic field decays faster than the kinetic energy of the component along the magnetic field (Moffatt, 1967). However, numerous numerical studies predict a different behaviour, wherein the final state is characterised by dominant horizontal energy. We investigate the corresponding nonlinear phenomenon using Direct Numerical Simulations (DNS) and spectral closures based on Eddy Damping Quasi-Normal Markovian (EDQNM) models. The initial temporal evolution of the decaying flow indicates that the turbulence is very similar to the so-called 'two-and-a-half-dimensional' flow (Montgomery and Turner, 1982) which explains the observations in numerical studies. EDQNM models confirm this statement at higher Reynolds number.
Kolmogorov's law for two-dimensional electron-magnetohydrodynamic turbulence
Celani, A.; Prandi, R.; Boffetta, G
1997-01-01
The analogue of the Kolmogorov's four-fifths law is derived for two-dimensional, homogeneous, isotropic EMHD turbulence in the energy cascade inertial range. Direct numerical simulations for the freely decaying case show that this relation holds true for different values of the adimensional electron inertial length scale, $d_e$. The energy spectrum is found to be close to the expected Kolmogorov spectrum.
Usmanov, Arcadi V.; Goldstein, Melvyn L.; Matthaeus, William H.
2012-01-01
To study the effects of interstellar pickup protons and turbulence on the structure and dynamics of the solar wind, we have developed a fully three-dimensional magnetohydrodynamic solar wind model that treats interstellar pickup protons as a separate fluid and incorporates the transport of turbulence and turbulent heating. The governing system of equations combines the mean-field equations for the solar wind plasma, magnetic field, and pickup protons and the turbulence transport equations for the turbulent energy, normalized cross-helicity, and correlation length. The model equations account for photoionization of interstellar hydrogen atoms and their charge exchange with solar wind protons, energy transfer from pickup protons to solar wind protons, and plasma heating by turbulent dissipation. Separate mass and energy equations are used for the solar wind and pickup protons, though a single momentum equation is employed under the assumption that the pickup protons are comoving with the solar wind protons.We compute the global structure of the solar wind plasma, magnetic field, and turbulence in the region from 0.3 to 100 AU for a source magnetic dipole on the Sun tilted by 0 deg - .90 deg and compare our results with Voyager 2 observations. The results computed with and without pickup protons are superposed to evaluate quantitatively the deceleration and heating effects of pickup protons, the overall compression of the magnetic field in the outer heliosphere caused by deceleration, and the weakening of corotating interaction regions by the thermal pressure of pickup protons.
Mininni, Pablo; Odier, Philippe; Pinton, Jean-François; Plihon, Nicolas; Verhille, Gautier; Volk, Romain; Bourgoin, Mickael
2014-01-01
We analyze time series stemming from experiments and direct numerical simulations of hydrodynamic and magnetohydrodynamic turbulence. Simulations are done in periodic boxes, but with a volumetric forcing chosen to mimic the geometry of the flow in the experiments, the von K\\'arm\\'an swirling flow between two counter-rotating impellers. Parameters in the simulations are chosen to (within computational limitations) allow comparisons between the experiments and the numerical results. Conducting fluids are considered in all cases. Two different configurations are considered: a case with a weak externally imposed magnetic field, and a case with self-sustained magnetic fields. Evidence of long-term memory and $1/f$ noise is observed in experiments and simulations, in the case with weak magnetic field associated with the hydrodynamic behavior of the shear layer in the von K\\'arm\\'an flow, and in the dynamo case associated with slow magnetohydrodynamic behavior of the large scale magnetic field.
Cosmic ray propagation in sub-Alfvénic magnetohydrodynamic turbulence
Cohet, R.; Marcowith, A.
2016-04-01
Context. The propagation of cosmic rays or energetic charged particles in magnetized turbulence is a complex problem which involves non-linear wave-particle interactions and chaotic magnetic field lines transport. This problem has been addressed until recently using either analytical calculations or simulations using prescribed turbulence models. With the advent of super computers it is now possible to investigate energetic charged particle propagation using direct computation of electromagnetic fields. This is in particular the case for high-energy particles propagation in magnetohydrodynamic turbulence. Aims: This work has the main objective to provide a detailed investigation of cosmic ray propagation in magnetohydrodynamic turbulent fields generated by forcing the fluid velocity field at large scales. It provides a derivation of the particle mean free path dependences in terms of the turbulence level described by the Alfvénic Mach number and in terms of the particle rigidity. Methods: We use an upgrade version of the magnetohydrodynamic code RAMSES which includes a forcing module and a kinetic module and solve the Lorentz equation for each particle. The simulations are performed using a 3 dimension periodical box in the test-particle and magnetostatic limits. The forcing module is implemented using an Ornstein-Uhlenbeck process. An ensemble average over a large number of particle trajectories is applied to reconstruct the particle mean free paths. Results: We derive the cosmic ray mean free paths in terms of the Alfvénic Mach numbers and particle reduced rigidities in different turbulence forcing geometries. The reduced particle rigidity is ρ = rL/L where rL is the particle Larmor radius and L is the simulation box length related to the turbulence coherence or injection scale Linj by L ~ 5 Linj. We have investigated with a special attention compressible and solenoidal forcing geometries. Conclusions: We find that compressible forcing solutions are compatible
Banerjee, Supratik
2013-01-01
Compressible isothermal magnetohydrodynamic turbulence is analyzed under the assumption of statistical homogeneity and in the asymptotic limit of large kinetic and magnetic Reynolds numbers. Following Kolmogorov we derive an exact relation for some two-point correlation functions which generalizes the expression recently found for hydrodynamics. We show that the magnetic field brings new source and flux terms into the dynamics which may act on the inertial range similarly as a source or a sink for the mean energy transfer rate. The introduction of a uniform magnetic field simplifies significantly the exact relation for which a simple phenomenology may be given. A prediction for axisymmetric energy spectra is eventually proposed.
Chandran
2000-11-27
Scattering rates for a Goldreich-Sridhar (GS) spectrum of anisotropic, incompressible, magnetohydrodynamic turbulence are calculated in the quasilinear approximation. Because the small-scale fluctuations are constrained to have wave vectors nearly perpendicular to the background magnetic field, scattering is too weak to provide either the mean-free paths commonly used in Galactic cosmic-ray propagation models or the mean-free paths required for acceleration of cosmic rays at quasiparallel shocks. Where strong pitch-angle scattering occurs, it is due to fluctuations not described by the GS spectrum, such as fluctuations generated by streaming cosmic rays. PMID:11082620
Direct Numerical Simulation of turbulent magnetohydrodynamic flows in an open-channel
Yamamoto, Y.; Kunugi, T.; Serizawa, A. [Kyoto Univ. (Japan). Dept. of Nuclear Engineering
2000-10-01
In this study, Direct Numerical Simulation (DNS) of magnetohydrodynamic (MHD) flows simulated the liquid wall concept in nuclear fusion applications, i.e., turbulent open-channel flows heating under the isoflux condition at the free-surface and insulting condition at the bottom wall, was employed. The magnetic field was imposed in the spanwise direction and the erect of Lorenz force was expressed in terms of the electrostatic potential with an assumption of a low magnetic Reynolds number. As the result, new method for the thermal boundary condition at the free-surface is advanced and numerical investigation of MHD flow was conducted. (author)
Reconnection-Driven Magnetohydrodynamic Turbulence in a Simulated Coronal-Hole Jet
Uritsky, Vadim M; DeVore, C Richard; Karpen, Judith T
2016-01-01
Extreme-ultraviolet and X-ray jets occur frequently in magnetically open coronal holes on the Sun, especially at high solar latitudes. Some of these jets are observed by white-light coronagraphs as they propagate through the outer corona toward the inner heliosphere, and it has been proposed that they give rise to microstreams and torsional Alfv\\'{e}n waves detected in situ in the solar wind. To predict and understand the signatures of coronal-hole jets, we have performed a detailed statistical analysis of such a jet simulated with an adaptively refined magnetohydrodynamics model. The results confirm the generation and persistence of three-dimensional, reconnectiondriven magnetic turbulence in the simulation. We calculate the spatial correlations of magnetic fluctuations within the jet and find that they agree best with the M\\"{u}ller - Biskamp scaling model including intermittent current sheets of various sizes coupled via hydrodynamic turbulent cascade. The anisotropy of the magnetic fluctuations and the sp...
On the compressibility effect in test particle acceleration by magnetohydrodynamic turbulence
González, C A; Mininni, P D; Matthaeus, W H
2016-01-01
The effect of compressibility in charged particle energization by magnetohydrodynamic (MHD) fields is studied in the context of test particle simulations. This problem is relevant to the solar wind and the solar corona due to the compressible nature of the flow in those astrophysical scenarios. We consider turbulent electromagnetic fields obtained from direct numerical simulations of the MHD equations with a strong background magnetic field. In order to explore the compressibilty effect over the particle dynamics we performed different numerical experiments: an incompressible case, and two weak compressible cases with Mach number M = 0.1 and M = 0.25. We analyze the behavior of protons and electrons in those turbulent fields, which are well known to form aligned current sheets in the direction of the guide magnetic field. We show that compressibility enhances the efficiency of proton acceleration, and that the energization is caused by perpendicular electric fields generated between currents sheets. On the ot...
Self-organisation and non-linear dynamics in driven magnetohydrodynamic turbulent flows
Dallas, Vassilios
2014-01-01
Magnetohydrodynamic turbulent flows driven by random mechanical and electromagnetic external forces of zero helicities are investigated by means of direct numerical simulations. It is shown that despite the absence of helicities in the forcing, the system is attracted to self-organized helical states that exhibit laminar behaviour despite the large value of the Reynolds numbers examined. We demonstrate that the correlation time of the external forces is controlling the time spent on these states, i.e. for short correlation times the system remains in the turbulent state while as the correlation time is increased the system spends more and more time in the self-organised states. As a result, time averaged statistics can significantly be affected by the time spent on these states. These results have important theoretical implications for the understanding of the suppression of non-linearities in plasma fusion devises as well as in astrophysical observations.
Discrete filters for large-eddy simulation of forced compressible magnetohydrodynamic turbulence
Chernyshov, A. A.; Petrosyan, A. S.
2016-06-01
We discuss results of the applicability of discrete filters for the large-eddy simulation (LES) method of forced compressible magnetohydrodynamic (MHD) turbulent flows with the scale-similarity model. New results are obtained for cross-helicity and residual energy. Cross-helicity and residual energy are important quantities in magnetohydrodynamic turbulence and have no hydrodynamic counterpart. The influences and effects of discrete filter shapes on the scale-similarity model are examined in physical space using finite-difference numerical schemes. We restrict ourselves to the Gaussian filter and the top-hat filter. Representations of this subgrid-scale model, which correspond to various 3- and 5-point approximations of both Gaussian and top-hat filters for different values of parameter ε (the ratio of the cut-off length-scale of the filter to the mesh size), are investigated. Discrete filters produce more discrepancies for the magnetic field. It is shown that the Gaussian filter is more sensitive to the parameter ɛ than the top-hat filter in compressible forced MHD turbulence. The 3-point filters at ε =2 and ε =3 give the least accurate results whereas the 5-point Gaussian filter shows the best results at ε =2 and ε =3. There are only very small differences deep into the dissipation region in favor of ε =2. For cross-helicity, the 5-point discrete filters are in good agreement with the results of direct numerical simulation (DNS), while the 3-point filter produces the largest discrepancies with DNS results. There is no strong dependence on the choice of the parameter ε and order approximation is a much more important factor for the cross-helicity. The difference between the filters is less for the residual energy compared with total energy. Thus, the total energy is more sensitive to the choice of discrete filter in the modeling of compressible MHD turbulence using the LES method.
Field theoretic calculation of energy cascade rates in non-helical magnetohydrodynamic turbulence
Mahendra K Verma
2003-09-01
Energy cascade rates and Kolmogorov’s constant for non-helical steady magnetohydrodynamic turbulence have been calculated by solving the ﬂux equations to the ﬁrst order in perturbation. For zero cross helicity and space dimension = 3, magnetic energy cascades from large length-scales to small length-scales (forward cascade). In addition, there are energy ﬂuxes from large-scale magnetic ﬁeld to small-scale velocity ﬁeld, large-scale velocity ﬁeld to small-scale magnetic ﬁeld, and large-scale velocity ﬁeld to large-scale magnetic ﬁeld. Kolmogorov’s constant for magnetohydrodynamics is approximately equal to that for ﬂuid turbulence (≈ 1.6) for Alfvén ratio 0.5 ≤ A ≤ ∞. For higher space-dimensions, the energy ﬂuxes are qualitatively similar, and Kolmogorov’s constant varies as 1/3. For the normalized cross helicity c → 1, the cascade rates are proportional to (1-c)/(1+c), and the Kolmogorov’s constants vary signiﬁcantly with c.
Hall effects and sub-grid-scale modeling in magnetohydrodynamic turbulence simulations
Miura, Hideaki; Araki, Keisuke; Hamba, Fujihiro
2016-07-01
Effects of the Hall term on short-wave components of magnetohydrodynamic turbulence and sub-grid-scale modeling of the effects are studied. Direct numerical simulations of homogeneous magnetohydrodynamic turbulence with and without the Hall term are carried out. The Hall term excites short-wave components in the magnetic field, demanding a high numerical resolution to resolve the scales smaller than the ion skin depth. A k 7 / 3-like scaling-law in the magnetic energy spectrum associated with the excitation of the short-wave components is clearly shown by the use of both an isotropic spectrum and a one-dimensional spectrum. It is also shown that the introduction of the Hall term can cause a structural transition in the vorticity field from tubes to sheets. In order to overcome a strong demand on high-resolution in space and time and to enable quicker computations, large eddy simulations with a Smagorinsky-type sub-grid-scale model are carried out. It is shown that our large eddy simulations successfully reproduce not only the energy spectrum but also tubular vortex structures, reducing the computational cost considerably.
Field theoretic calculation of energy cascade rates in non-helical magnetohydrodynamic turbulence
Mahendra K Verma
2004-06-01
Energy cascade rates and Kolmogorov’s constant for non-helical steady magnetohydrodynamic turbulence have been calculated by solving the ﬂux equations to the ﬁrst order in perturbation. For zero cross helicity and space dimension $d = 3$, magnetic energy cascades from large length-scales to small length-scales (forward cascade). In addition, there are energy ﬂuxes from large-scale magnetic ﬁeld to small-scale velocity ﬁeld, large-scale velocity ﬁeld to small-scale magnetic ﬁeld, and large-scale velocity ﬁeld to large-scale magnetic ﬁeld. Kolmogorov’s constant for magnetohydrodynamics is approximately equal to that for ﬂuid turbulence $(≈ 1.6)$ for Alfvén ratio $0.5≤ r_{A}≤ ∞$. For higher space-dimensions, the energy ﬂuxes are qualitatively similar, and Kolmogorov’s constant varies as $d^{1/3}$. For the normalized cross helicity $_{c}→ 1$, the cascade rates are proportional to $(1-_{c})/(1+_{c})$, and the Kolmogorov’s constants vary signiﬁcantly with $_{c}$.
Solar wind low-frequency magnetohydrodynamic turbulence: extended self-similarity and scaling laws
P. Veltri
1996-01-01
Full Text Available In this paper we review some of the work done in investigating the scaling properties of Magnetohydrodynamic turbulence, by using velocity fluctuations measurements performed in the interplanetary space plasma by the Helios spacecraft. The set of scaling exponents ξq for the q-th order velocity structure functions, have been determined by using the Extended Self-Similarity hypothesis. We have found that the q-th order velocity structure function, when plotted vs. the 4-th order structure function, displays a range of self-similarity which extends over all the lengths covered by measurements, thus allowing for a very good determination of ξq. Moreover the results seem to show that the scaling exponents are the same regardless the various observation periods considered. The obtained scaling exponents have been compared with the results of some intermittency models for Kraichnan's turbulence, derived in the framework of infinitely divisible fragmentation processes, showing the good agreement between these models and our observations. Finally, on the basis of the actually available data sets, we show that scaling laws in Solar Wind turbulence seem to be different from turbulent scaling laws in the ordinary fluid flows. This is true for high-order velocity structure functions, while low-order velocity structure functions show the same scaling laws. Since our measurements involve length scales which extend over many order of magnitude where dissipation is practically absent, our results show that Solar Wind turbulence can be regarded as a testing bench for the investigation of general scaling behaviour in turbulent flows. In particular our results strongly support the point of view which attributes a key role to the inertial range dynamics in determining the intermittency characteristics in fluid flows, in contrast with the point of view which attributes intermittency to a finite Reynolds number effect.
Weakly Compressible Magnetohydrodynamic Turbulence in the Solar Wind and the Interstellar Medium
Bhattacharjee, A.; Ng, C. S.; Spangler, S. R.
1998-02-01
A new four-field system of equations is derived from the compressible magnetohydrodynamic (MHD) equations for low Mach number turbulence in the solar wind and the interstellar medium, permeated by a spatially varying magnetic field. The plasma beta is assumed to be of order unity or less. It is shown that the full MHD equations can be reduced rigorously to a closed system for four fluctuating field variables: magnetic flux, vorticity, pressure, and parallel flow. Although the velocity perpendicular to the magnetic field is shown to obey a two-dimensional incompressibility condition (analogous to the Proudman-Taylor theorem in hydrodynamics), the three-dimensional dynamics exhibit the effects of compressibility. In the presence of spatial inhomogeneities, the four dynamical equations are coupled to each other, and pressure fluctuations enter the weakly compressible dynamics at leading order. If there are no spatial inhomogeneities and or if the plasma beta is low, the four-field equations reduce to the well-known equations of reduced magnetohydrodynamics (RMHD). For pressure-balanced structures, the four-field equations undergo a remarkable simplification which provides insight on the special nature of the fluctuations driven by these structures. The important role of spatial inhomogeneities is elucidated by 2.5-dimensional numerical simulations. In the presence of inhomogeneities, the saturated pressure and density fluctuations scale with the Mach number of the turbulence, and the system attains equipartition with respect to the kinetic, magnetic, and thermal energy of the fluctuations. The present work suggests that if heliospheric and interstellar turbulence exists in a plasma with large-scale, nonturbulent spatial gradients, one expects the pressure and density fluctuations to be of significantly larger magnitude than suggested in nearly incompressible models such as pseudosound.
Scaling and anisotropy in magnetohydrodynamic turbulence in a strong mean magnetic field
We present an analysis of the anisotropic spectral energy distribution in incompressible magnetohydrodynamic turbulence permeated by a strong mean magnetic field. The turbulent flow is generated by high-resolution pseudospectral direct numerical simulations with large-scale isotropic forcing. Examining the radial energy distribution for various angles θ with respect to B0 reveals a specific structure which remains hidden when not taking axial symmetry with respect to B0 into account. For each direction, starting at the forced large scales, the spectrum first exhibits an amplitude drop around a wave number k0 which marks the start of a scaling range and goes on up to a dissipative wave number kd(θ). The three-dimensional spectrum for k≥k0 is described by a single θ-independent functional form F(k/kd), with the scaling law being the same in every direction. The previous properties still hold when increasing the mean field from B0=5 up to B0=10brms, as well as when passing from resistive to ideal flows. We conjecture that at fixed B0 the direction-independent scaling regime is reached when increasing the Reynolds number above a threshold which raises with increasing B0. Below that threshold critically balanced turbulence is expected.
Lagrangian-averaged model for magnetohydrodynamic turbulence and the absence of bottlenecks.
Pietarila Graham, Jonathan; Mininni, Pablo D; Pouquet, Annick
2009-07-01
We demonstrate that, for the case of quasiequipartition between the velocity and the magnetic field, the Lagrangian-averaged magnetohydrodynamics (LAMHD) alpha model reproduces well both the large-scale and the small-scale properties of turbulent flows; in particular, it displays no increased (superfilter) bottleneck effect with its ensuing enhanced energy spectrum at the onset of the subfilter scales. This is in contrast to the case of the neutral fluid in which the Lagrangian-averaged Navier-Stokes alpha model is somewhat limited in its applications because of the formation of spatial regions with no internal degrees of freedom and subsequent contamination of superfilter-scale spectral properties. We argue that, as the Lorentz force breaks the conservation of circulation and enables spectrally nonlocal energy transfer (associated with Alfvén waves), it is responsible for the absence of a viscous bottleneck in magnetohydrodynamics (MHD), as compared to the fluid case. As LAMHD preserves Alfvén waves and the circulation properties of MHD, there is also no (superfilter) bottleneck found in LAMHD, making this method capable of large reductions in required numerical degrees of freedom; specifically, we find a reduction factor of approximately 200 when compared to a direct numerical simulation on a large grid of 1536;{3} points at the same Reynolds number. PMID:19658812
Driving Solar Spicules and Jets with Magnetohydrodynamic Turbulence: Testing a Persistent Idea
Cranmer, Steven R
2015-01-01
The solar chromosphere contains thin, highly dynamic strands of plasma known as spicules. Recently, it has been suggested that the smallest and fastest (Type II) spicules are identical to intermittent jets observed by the Interface Region Imaging Spectrograph. These jets appear to expand out along open magnetic field lines rooted in unipolar network regions of coronal holes. In this paper we revisit a thirty-year-old idea that spicules may be caused by upward forces associated with Alfven waves. These forces involve the conversion of transverse Alfven waves into compressive acoustic-like waves that steepen into shocks. The repeated buffeting due to upward shock propagation causes nonthermal expansion of the chromosphere and a transient levitation of the transition region. Some older models of wave-driven spicules assumed sinusoidal wave inputs, but the solar atmosphere is highly turbulent and stochastic. Thus, we model this process using the output of a time-dependent simulation of reduced magnetohydrodynamic...
On the spatio-temporal behavior of magnetohydrodynamic turbulence in a magnetized plasma
Lugones, R; Mininni, P D; Wan, M; Matthaeus, W H
2016-01-01
Using direct numerical simulations of three-dimensional magnetohydrodynamic (MHD) turbulence the spatio-temporal behavior of magnetic field fluctuations is analyzed. Cases with relatively small, medium and large values of a mean background magnetic field are considered. The (wavenumber) scale dependent time correlation function is directly computed for different simulations, varying the mean magnetic field value. From this correlation function the time decorrelation is computed and compared with different theoretical times, namely, the local non-linear time, the random sweeping time, and the Alfv\\'enic time, the latter being a wave effect. It is observed that time decorrelations are dominated by sweeping effects, and only at large values of the mean magnetic field and for wave vectors mainly aligned with this field time decorrelations are controlled by Alfv\\'enic effects.
Basu, Abhik; Naji, Ali; Pandit, Rahul
2014-01-01
We generalize the method of A. M. Polyakov, [ Phys. Rev. E 52 6183 (1995)] for obtaining structure-function relations in turbulence in the stochastically forced Burgers equation, to develop structure-function hierarchies for turbulence in three models for magnetohydrodynamics (MHD). These are the Burgers analogs of MHD in one dimension [ Eur. Phys. J. B 9 725 (1999)], and in three dimensions (3DMHD and 3D Hall MHD). Our study provides a convenient and unified scheme for the development of structure-function hierarchies for turbulence in a variety of coupled hydrodynamical equations. For turbulence in the three sets of MHD equations mentioned above, we obtain exact relations for third-order structure functions and their derivatives; these expressions are the analogs of the von Kármán-Howarth relations for fluid turbulence. We compare our work with earlier studies of such relations in 3DMHD and 3D Hall MHD. PMID:24580182
Magnetic field amplification and saturation in two-dimensional magnetohydrodynamic turbulence
Two-dimensional (2-D) magnetohydrodynamic turbulence is investigated for weak initial magnetic fields using numerical simulation. It is found that the magnetic field is amplified owing to the formation of flux sheets with saturation due either to resistive diffusion (kinematic regime) or to nonlinear effects (dynamic regime). In the kinematic regime, which corresponds to the problem of passive scalar convection by 2-D Navier--Stokes turbulence, the saturation value of the magnetic energy is observed to scale as EMmax∝η-0.8 in approximate agreement with a simple theoretical estimate, EMmax/EM(0)congruent Rm, where Rm is the magnetic Reynolds number. Because of the strongly disparate kinetic and magnetic energy spectra in the kinematic regime, roughly EVk∼k-3, EMk∼k, dynamic interaction on small scales already occurs at very small global energy ratios EM/EV, giving rise to strongly enhanced kinetic energy dissipation. In the fully dynamic regime (reached for EM/EV|t=0>R-1m) global magnetic and kinetic energies become tightly coupled, with EM/EV being approximately constant in time and the energy dissipation rates being independent of the collisional diffusion coefficients. Finally, the effect of the magnetic Prandtl number Pr=μ/η is discussed
Driving Solar Spicules and Jets with Magnetohydrodynamic Turbulence: Testing a Persistent Idea
Cranmer, Steven R.; Woolsey, Lauren N.
2015-10-01
The solar chromosphere contains thin, highly dynamic strands of plasma known as spicules. Recently, it has been suggested that the smallest and fastest (Type II) spicules are identical to intermittent jets observed by the Interface Region Imaging Spectrograph. These jets appear to expand out along open magnetic field lines rooted in unipolar network regions of coronal holes. In this paper we revisit a thirty-year-old idea that spicules may be caused by upward forces associated with Alfvén waves. These forces involve the conversion of transverse Alfvén waves into compressive acoustic-like waves that steepen into shocks. The repeated buffeting due to upward shock propagation causes nonthermal expansion of the chromosphere and a transient levitation of the transition region (TR). Some older models of wave-driven spicules assumed sinusoidal wave inputs, but the solar atmosphere is highly turbulent and stochastic. Thus, we model this process using the output of a time-dependent simulation of reduced magnetohydrodynamic turbulence. The resulting mode-converted compressive waves are strongly variable in time, with a higher TR occurring when the amplitudes are large and a lower TR when the amplitudes are small. In this picture, the TR bobs up and down by several Mm on timescales less than a minute. These motions produce narrow, intermittent extensions of the chromosphere that have similar properties as the observed jets and Type II spicules.
Antonov, N. V.; Kostenko, M. M.
2015-11-01
The field-theoretic renormalization group and the operator product expansion are applied to the model of passive vector (magnetic) field advected by a random turbulent velocity field. The latter is governed by the Navier-Stokes equation for compressible fluid, subject to external random force with the covariance ∝δ (t -t') k4 -d -y , where d is the dimension of space and y is an arbitrary exponent. From physics viewpoints, the model describes magnetohydrodynamic turbulence in the so-called kinematic approximation, where the effects of the magnetic field on the dynamics of the fluid are neglected. The original stochastic problem is reformulated as a multiplicatively renormalizable field-theoretic model; the corresponding renormalization group equations possess an infrared attractive fixed point. It is shown that various correlation functions of the magnetic field and its powers demonstrate anomalous scaling behavior in the inertial-convective range already for small values of y . The corresponding anomalous exponents, identified with scaling (critical) dimensions of certain composite fields ("operators" in the quantum-field terminology), can be systematically calculated as series in y . The practical calculation is performed in the leading one-loop approximation, including exponents in anisotropic contributions. It should be emphasized that, in contrast to Gaussian ensembles with finite correlation time, the model and the perturbation theory presented here are manifestly Galilean covariant.
Scaling and anisotropy of magnetohydrodynamic turbulence in a strong mean magnetic field
Grappin, Roland
2010-01-01
We present a new analysis of the anisotropic spectral energy distribution in incompressible magnetohydrodynamic (MHD) turbulence permeated by a strong mean magnetic field. The turbulent flow is generated by high-resolution pseudo-spectral direct numerical simulations with large-scale isotropic forcing. Examining the radial energy distribution for various angles $\\theta$ with respect to $\\mathbf{B}_0$ reveals a specific structure which remains hidden when not taking axial symmetry with respect to $B_0$ into account. For each direction, starting at the forced large-scales, the spectrum first exhibits an amplitude drop around a wavenumber $k_0$ which marks the start of a scaling range and goes on up to a dissipative wavenumber $k_d(\\theta)$. The 3D spectrum for $k \\ge k_0$ is described by a single $\\theta$-independent functional form $F(k/k_d)$, the scaling law being the same in every direction. The previous properties still hold when increasing the mean field from $B_0=5$ up to $B_0=10 \\ b_{rms}$, as well as wh...
Antonov, N V; Kostenko, M M
2015-11-01
The field-theoretic renormalization group and the operator product expansion are applied to the model of passive vector (magnetic) field advected by a random turbulent velocity field. The latter is governed by the Navier-Stokes equation for compressible fluid, subject to external random force with the covariance ∝ δ(t-t')k(4-d-y), where d is the dimension of space and y is an arbitrary exponent. From physics viewpoints, the model describes magnetohydrodynamic turbulence in the so-called kinematic approximation, where the effects of the magnetic field on the dynamics of the fluid are neglected. The original stochastic problem is reformulated as a multiplicatively renormalizable field-theoretic model; the corresponding renormalization group equations possess an infrared attractive fixed point. It is shown that various correlation functions of the magnetic field and its powers demonstrate anomalous scaling behavior in the inertial-convective range already for small values of y. The corresponding anomalous exponents, identified with scaling (critical) dimensions of certain composite fields ("operators" in the quantum-field terminology), can be systematically calculated as series in y. The practical calculation is performed in the leading one-loop approximation, including exponents in anisotropic contributions. It should be emphasized that, in contrast to Gaussian ensembles with finite correlation time, the model and the perturbation theory presented here are manifestly Galilean covariant. PMID:26651785
V. Carbone
Full Text Available Some signed measures in turbulence are found to be sign-singular, that is their sign reverses continuously on arbitrary finer scales with a reduction of the cancellation between positive and negative contributions. The strength of the singularity is characterized by a scaling exponent κ, the cancellation exponent. In the present study by using some turbulent samples of the velocity field obtained from spacecraft measurements in the interplanetary medium, we show that sign-singularity is present everywhere in low-frequency turbulent samples. The cancellation exponent can be related to the characteristic scaling laws of turbulence. Differences in the values of κ, calculated in both high- and low-speed streams, allow us to outline some physical differences in the samples with different velocities.
On the compressibility effect in test particle acceleration by magnetohydrodynamic turbulence
González, C. A.; Dmitruk, P.; Mininni, P. D.; Matthaeus, W. H.
2016-08-01
The effect of compressibility in a charged particle energization by magnetohydrodynamic (MHD) fields is studied in the context of test particle simulations. This problem is relevant to the solar wind and the solar corona due to the compressible nature of the flow in those astrophysical scenarios. We consider turbulent electromagnetic fields obtained from direct numerical simulations of the MHD equations with a strong background magnetic field. In order to explore the flow compressibility effect over the particle dynamics, we performed different numerical experiments: an incompressible case and two weak compressible cases with Mach number M = 0.1 and M = 0.25. We analyze the behavior of protons and electrons in those turbulent fields, which are well known to form aligned current sheets in the direction of the guide magnetic field. What we call protons and electrons are test particles with scales comparable to (for protons) and much smaller than (for electrons) the dissipative scale of MHD turbulence, maintaining the correct mass ratio m e / m i . For these test particles, we show that compressibility enhances the efficiency of proton acceleration, and that the energization is caused by perpendicular electric fields generated between currents sheets. On the other hand, electrons remain magnetized and display an almost adiabatic motion, with no effect of compressibility observed. Another set of numerical experiments takes into account two fluid modifications, namely, electric field due to Hall effect and electron pressure gradient. We show that the electron pressure has an important contribution to electron acceleration allowing highly parallel energization. In contrast, no significant effect of these additional terms is observed for the protons.
Arthur, S. J.; Henney, W. J.; Mellema, G.; de Colle, F.; Vázquez-Semadeni, E.
2011-06-01
We present the results of radiation-magnetohydrodynamic simulations of the formation and expansion of H II regions and their surrounding photodissociation regions (PDRs) in turbulent, magnetized, molecular clouds on scales of up to 4 pc. We include the effects of ionizing and non-ionizing ultraviolet radiation and X-rays from population synthesis models of young star clusters. For all our simulations we find that the H II region expansion reduces the disordered component of the magnetic field, imposing a large-scale order on the field around its border, with the field in the neutral gas tending to lie along the ionization front, while the field in the ionized gas tends to be perpendicular to the front. The highest pressure-compressed neutral and molecular gas is driven towards approximate equipartition between thermal, magnetic and turbulent energy densities, whereas lower pressure neutral/molecular gas bifurcates into, on the one hand, quiescent, magnetically dominated regions and, on the other hand, turbulent, demagnetized regions. The ionized gas shows approximate equipartition between thermal and turbulent energy densities, but with magnetic energy densities that are 1-3 orders of magnitude lower. A high velocity dispersion (˜8 km s-1) is maintained in the ionized gas throughout our simulations, despite the mean expansion velocity being significantly lower. The magnetic field does not significantly brake the large-scale H II region expansion on the length and time-scales accessible to our simulations, but it does tend to suppress the smallest scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the H II region. However, the relative luminosity of ionizing and non-ionizing radiation has a much larger influence than the presence or absence of the magnetic field. When the star cluster radiation field is relatively soft (as in the case of a lower mass cluster, containing an earliest spectral
Burkhart, Blakesley; Ossenkopf, V.; Lazarian, A.; Stutzki, J.
2013-07-01
We study the effects of radiative transfer on the probability distribution functions (PDFs) of simulations of magnetohydrodynamic turbulence in the widely studied 13CO 2-1 transition. We find that the integrated intensity maps generally follow a log-normal distribution, with the cases that have τ ≈ 1 best matching the PDF of the column density. We fit a two-dimensional variance-sonic Mach number relationship to our logarithmic PDFs of the form \\sigma _{\\ln (\\Sigma /\\Sigma _0)}^2=A\\times \\ln (1+b^2{\\cal M}_s^2) and find that, for parameter b = 1/3, parameter A depends on the radiative transfer environment. We also explore the variance, skewness, and kurtosis of the linear PDFs finding that higher moments reflect both higher sonic Mach number and lower optical depth. Finally, we apply the Tsallis incremental PDF function and find that the fit parameters depend on both Mach numbers, but also are sensitive to the radiative transfer parameter space, with the τ ≈ 1 case best fitting the incremental PDF of the true column density. We conclude that, for PDFs of low optical depth cases, part of the gas is always subthermally excited so that the spread of the line intensities exceeds the spread of the underlying column densities and hence the PDFs do not reflect the true column density. Similarly, PDFs of optically thick cases are dominated by the velocity dispersion and therefore do not represent the true column density PDF. Thus, in the case of molecules like carbon monoxide, the dynamic range of intensities, structures observed, and, consequently, the observable PDFs are less determined by turbulence and more often determined by radiative transfer effects.
We study the effects of radiative transfer on the probability distribution functions (PDFs) of simulations of magnetohydrodynamic turbulence in the widely studied 13CO 2-1 transition. We find that the integrated intensity maps generally follow a log-normal distribution, with the cases that have τ ≈ 1 best matching the PDF of the column density. We fit a two-dimensional variance-sonic Mach number relationship to our logarithmic PDFs of the form σln2(Σ/Σ0) = A x ln(1+b2Ms2) and find that, for parameter b = 1/3, parameter A depends on the radiative transfer environment. We also explore the variance, skewness, and kurtosis of the linear PDFs finding that higher moments reflect both higher sonic Mach number and lower optical depth. Finally, we apply the Tsallis incremental PDF function and find that the fit parameters depend on both Mach numbers, but also are sensitive to the radiative transfer parameter space, with the τ ≈ 1 case best fitting the incremental PDF of the true column density. We conclude that, for PDFs of low optical depth cases, part of the gas is always subthermally excited so that the spread of the line intensities exceeds the spread of the underlying column densities and hence the PDFs do not reflect the true column density. Similarly, PDFs of optically thick cases are dominated by the velocity dispersion and therefore do not represent the true column density PDF. Thus, in the case of molecules like carbon monoxide, the dynamic range of intensities, structures observed, and, consequently, the observable PDFs are less determined by turbulence and more often determined by radiative transfer effects.
Hall effects on local structures in homogeneous, isotropic, and incompressible magnetohydrodynamic turbulence are studied numerically. The transition of vortices from sheet-like to tubular structures induced by the Hall term is found, while the kinetic energy spectrum does not distinguish the two types of structures. It is shown by the use of the sharp low-pass filter that the transition occurs not only in the scales smaller than the ion skin depth but also in a larger scale. The transition is related with the forward energy transfer in the spectral space. Analyses by the use of the sharp low-pass filter show that the nonlinear energy transfer associated with the Hall term is dominated by the forward transfer and relatively local in the wave number space. A projection of the simulation data to a Smagorinsky-type sub-grid-scale model shows that the high wave number component of the Hall term may possibly be replaced by the model effectively
Miura, H., E-mail: miura.hideaki@nifs.ac.jp [Department of Helical Plasma Research, National Institute for Fusion Science, 322-6 Oroshi, Toki, Gifu 509-5292 (Japan); Araki, K. [Faculty of Engineering, Okayama University of Science, 1-1 Ridai-cho, Okayama 700-0005 (Japan)
2014-07-15
Hall effects on local structures in homogeneous, isotropic, and incompressible magnetohydrodynamic turbulence are studied numerically. The transition of vortices from sheet-like to tubular structures induced by the Hall term is found, while the kinetic energy spectrum does not distinguish the two types of structures. It is shown by the use of the sharp low-pass filter that the transition occurs not only in the scales smaller than the ion skin depth but also in a larger scale. The transition is related with the forward energy transfer in the spectral space. Analyses by the use of the sharp low-pass filter show that the nonlinear energy transfer associated with the Hall term is dominated by the forward transfer and relatively local in the wave number space. A projection of the simulation data to a Smagorinsky-type sub-grid-scale model shows that the high wave number component of the Hall term may possibly be replaced by the model effectively.
Alfvén-dynamo balance and magnetic excess in magnetohydrodynamic turbulence
Grappin, Roland; Müller, Wolf-Christian; Verdini, Andrea
2016-05-01
Context. Three-dimensional magnetohydrodynamic (3D MHD) turbulent flows with initially magnetic and kinetic energies at equipartition spontaneously develop a magnetic excess (or residual energy) in both numerical simulations and the solar wind. Closure equations obtained in 1983 describe the residual spectrum as resulting from a balance between a dynamo source proportional to the total energy spectrum and a linear Alfvén damping term. A good agreement was found in 2005 with incompressible simulations; however, recent solar wind measurements disagree with these results. Aims: The previous dynamo-Alfvén theory is generalized to a family of models, leading to simple relations between residual and total energy spectra. We want to assess these models in detail against MHD simulations and solar wind data. Methods: We tested the family of models against compressible decaying MHD simulations with a low Mach number, low cross-helicity, and zero-mean magnetic field with or without expansion terms (EBM; expanding box model). Results: A single dynamo-Alfvén model is found to describe correctly both solar wind scalings and compressible simulations without or with expansion. This model is equivalent to the 1983-2005 closure equation, but it incorporates the critical balance of nonlinear turnover and linear Alfvén times, while the dynamo source term remains unchanged. We elucidate the discrepancy with previous incompressible simulations. The model predicts a linear relation between the spectral slopes of total and residual energies mR = -1/2 + 3/2mT. By examining previous solar wind data, our relation is found to be valid for any cross-helicity, and is even better at high cross-helicity with the total energy slope varying from 1.7 to 1.55.
Podesta, J. J.
It is known that Kolmogorov's four-fifths law for statistically homogeneous and isotropic turbulence can be generalized to anisotropic turbulence. This fundamental result for homogeneous anisotropic turbulence says that in the inertial range the divergence of the vector third-order moment |v(r) is constant and is equal to -4, where is the dissipation rate of the turbulence. This law can be extended to incompressible magnetohydrodyamic (MHD) turbulence where statistical isotropy is often not valid due, for example, to the presence of a large-scale magnetic field. Laws for anisotropic incompressible MHD turbulence were first derived by Politano and Pouquet. In this paper, the laws for vector third-order moments in homogeneous non-isotropic incompressible MHD turbulence are derived by a technique due to Frisch that clarifies the relationship between the energy flux in Fourier space and the vector third-order moments in physical space. This derivation is different from the original derivation of Politano and Pouquet which is based on the Kn-Howarth equation, and provides some new physical insights. Separate laws are derived for the cascades of energy, cross-helicity and magnetic-helicity, the three ideal invariants of incompressible MHD for flows in three dimensions. These laws are of fundamental importance in the theory of MHD turbulence where non-isotropic turbulence is much more prevalent than isotropic turbulence.
SHOCKFIND - An algorithm to identify magnetohydrodynamic shock waves in turbulent clouds
Lehmann, Andrew; Wardle, Mark
2016-01-01
The formation of stars occurs in the dense molecular cloud phase of the interstellar medium. Observations and numerical simulations of molecular clouds have shown that supersonic magnetised turbulence plays a key role for the formation of stars. Simulations have also shown that a large fraction of the turbulent energy dissipates in shock waves. The three families of MHD shocks --- fast, intermediate and slow --- distinctly compress and heat up the molecular gas, and so provide an important probe of the physical conditions within a turbulent cloud. Here we introduce the publicly available algorithm, SHOCKFIND, to extract and characterise the mixture of shock families in MHD turbulence. The algorithm is applied to a 3-dimensional simulation of a magnetised turbulent molecular cloud, and we find that both fast and slow MHD shocks are present in the simulation. We give the first prediction of the mixture of turbulence-driven MHD shock families in this molecular cloud, and present their distinct distributions of s...
Sujovolsky, N E
2016-01-01
We study the transition in dimensionality of a three-dimensional magnetohydrodynamic flow forced only mechanically, when the strength of a magnetic guiding field is gradually increased. We use numerical simulations to consider cases in which the mechanical forcing injects (or not) helicity in the flow. As the guiding field is increased, the strength of the magnetic field fluctuations decrease as a power law of the guiding field intensity. We show that for strong enough guiding fields, the helical magnetohydrodynamic flow can become almost two-dimensional. In this case, the mechanical energy can undergo a process compatible with an inverse cascade, being transferred preferentially towards scales larger than the forcing scale. The presence of helicity changes the spectral scaling of the small magnetic field fluctuations, and affects the statistics of the velocity field and of the velocity gradients. Moreover, at small scales the dynamics of the flow becomes dominated by a direct cascade of helicity, which can b...
Inertial ranges and resistive instabilities in two-dimensional magnetohydrodynamic turbulence
Direct numerical simulations of decaying two-dimensional magnetohydrodynamic flows at Reynolds numbers of several thousand are performed, using resolutions of 10242 collocation points. An inertial range extending to about one decade is observed, with spectral properties depending on the velocity--magnetic field correlation. At very small scales, resistive tearing destabilizes current sheets generated by the inertial dynamics and leads to the formation of small-scale magnetic islands, which may then grow and reach the size of inertial scales
Vetcha, N.; Smolentsev, S.; Abdou, M.; Moreau, R.
2013-02-01
We consider magnetohydrodynamic (MHD) rectangular duct flows with volumetric heating. The flows are upward, subject to a strong transverse magnetic field perpendicular to the temperature gradient, such that the flow dynamics is quasi-two-dimensional. The internal volumetric heating imitates conditions of a blanket of a fusion power reactor, where a buoyancy-driven flow is imposed on the forced flow. Studies of this mixed-convection flow include analysis for the basic flow, linear stability analysis and Direct Numerical Simulation (DNS)-type computations. The parameter range covers the Hartmann number (Ha) up to 500, the Reynolds number (Re) from 1000 to 10 000, and the Grashof number (Gr) from 105 to 5 × 108. The linear stability analysis predicts two primary instability modes: (i) bulk instability associated with the inflection point in the velocity profile near the "hot" wall and (ii) side-wall boundary layer instability. A mixed instability mode is also possible. An equation for the critical Hartmann number has been obtained as a function of Re and Gr. Effects of Ha, Re, and Gr on turbulent flows are addressed via nonlinear computations that demonstrate two characteristic turbulence regimes. In the "weak" turbulence regime, the induced vortices are localized near the inflection point of the basic velocity profile, while the boundary layer at the wall parallel to the magnetic field is slightly disturbed. In the "strong" turbulence regime, the bulk vortices interact with the boundary layer causing its destabilization and formation of secondary vortices that may travel across the flow, even reaching the opposite wall. In this regime, the key phenomena are vortex-wall and various vortex-vortex interactions. Flow and magnetic field effects on heat transfer are also analyzed.
Sub-Grid-Scale Description of Turbulent Magnetic Reconnection in Magnetohydrodynamics
Widmer, Fabien; Yokoi, Nobumitsu
2015-01-01
Magnetic reconnection requires, at least locally, a non-ideal plasma response. In collisionless space and astrophysical plasmas, turbulence could permit this instead of the too rare binary collisions. We investigated the influence of turbulence on the reconnection rate in the framework of a single fluid compressible MHD approach. The goal is to find out, whether unresolved, sub-grid for MHD simulations, turbulence can enhance the reconnection process in high Reynolds number astrophysical plasma. We solve, simultaneously with the grid-scale MHD equations, evolution equations for the sub-grid turbulent energy and cross helicity according to Yokoi's model (Yokoi (2013)) where turbulence is self-generated and -sustained through the inhomogeneities of the mean fields. Simulations of Harris and force free sheets confirm the results of Higashimori et al. (2013) and new results are obtained about the dependence on resistivity for large Reynolds number as well as guide field effects. The amount of energy transferred f...
The Energy Dissipation Rate of Supersonic, Magnetohydrodynamic Turbulence in Molecular Clouds
MacLow, M M
1998-01-01
Molecular clouds have broad linewidths suggesting turbulent supersonic motions in the clouds. These motions are usually invoked to explain why molecular clouds take much longer than a free-fall time to form stars. It has classically been thought that supersonic hydrodynamical turbulence would dissipate its energy quickly, but that the introduction of strong magnetic fields could maintain these motions. In a previous paper it has been shown, however, that isothermal, compressible, MHD and hydrodynamical turbulence decay at virtually the same rate, requiring that constant driving occur to maintain the observed turbulence. In this paper direct numerical computations of uniformly driven turbulence with the ZEUS astrophysical MHD code are used to derive the absolute value of energy dissipation as a function of the driving wavelength and amplitude. The ratio of the formal decay time of turbulence E_{kin}/\\dot{E}_{kin} to the free-fall time of the gas can then be derived as a function of the ratio of driving wavelen...
SHOCKFIND - An algorithm to identify magnetohydrodynamic shock waves in turbulent clouds
Lehmann, Andrew; Federrath, Christoph; Wardle, Mark
2016-08-01
The formation of stars occurs in the dense molecular cloud phase of the interstellar medium. Observations and numerical simulations of molecular clouds have shown that supersonic magnetised turbulence plays a key role for the formation of stars. Simulations have also shown that a large fraction of the turbulent energy dissipates in shock waves. The three families of MHD shocks - fast, intermediate and slow - distinctly compress and heat up the molecular gas, and so provide an important probe of the physical conditions within a turbulent cloud. Here we introduce the publicly available algorithm, SHOCKFIND, to extract and characterise the mixture of shock families in MHD turbulence. The algorithm is applied to a 3-dimensional simulation of a magnetised turbulent molecular cloud, and we find that both fast and slow MHD shocks are present in the simulation. We give the first prediction of the mixture of turbulence-driven MHD shock families in this molecular cloud, and present their distinct distributions of sonic and Alfvénic Mach numbers. Using subgrid one-dimensional models of MHD shocks we estimate that ˜0.03 % of the volume of a typical molecular cloud in the Milky Way will be shock heated above 50 K, at any time during the lifetime of the cloud. We discuss the impact of this shock heating on the dynamical evolution of molecular clouds.
Hydrodynamical and magnetohydrodynamic global bifurcations in a highly turbulent von Karman flow
We report experimental studies of the turbulent von Karman flow, inertially stirred between counter-rotating impellers. We first study the flow and its transition from laminar to turbulent regime. We highlight the role of slowly varying large scales, due to the presence of an azimuthal mixing layer. The large scales of this flow can be unstable in turbulent regime. We study the statistics of the transitions between the different mean states. The second part is dedicated to an experiment in liquid sodium, called VKS2. We optimize the time-averaged flow in order to allow kinematic dynamo action. We report the very first results of the experiment, and discuss the role of the large scales temporal non-stationariness. (author)
Park, Kiwan
2015-01-01
The inverse cascade of magnetic energy occurs when helicity or rotational instability exists in the magnetohydrodynamic (MHD) system. This well known phenomenon provides a basis for the large scale magnetic field in space. However even the decaying nonhelical magnetic energy can evolve to expand its scale. This phenomenon, inverse transfer of decaying nonhelical magnetic field may hold some vital clues to the origin of large scale magnetic field in the astrophysical system without helicity nor any significant driving source. Zeldovich's rope model has been considered as the basic principle with regard to the amplification of magnetic field. However, since the rope model assuming a driving force is not appropriate to the decaying system, we suggest a supplementary dynamo model based on the magnetic induction equation. The model explicitly shows the basic principle of migration and amplification of magnetic field. The expansion of scale and intensity of magnetic field is basically the consequent result of the r...
Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Astrophysics and Space Physics
Miesch, Mark S; Brandenburg, Axel; Petrosyan, Arakel; Pouquet, Annick; Cambon, Claude; Jenko, Frank; Uzdensky, Dmitri; Stone, James; Tobias, Steve; Toomre, Juri; Velli, Marco
2015-01-01
We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in astrophysics and space physics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, astrophysical and heliophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and...
Dynamics and statistics of inverse cascade processes in 2D magnetohydrodynamic turbulence
The inverse cascade of the mean square potential A in a 2D magnetofluid randomly forced at small scales is studied by numerical simulation. One finds the spectrum Ak ≅ 2.6 element of A2/3k-7/3. The cascade proceeds by coalescence of current filaments, which is a fast reconnection process owing to high turbulent resistivity. Statistics of δvl and δBl are strictly Gaussian, also in the condensation phase of Ak at k 1. Only when the coherent magnetic field intensity exceeds that of the fluctuations, non-Gaussian statistics in δB6 occur, which are, however, entirely due to the static magnetic structure and not associated with intermittency of the small-scale turbulence, the latter remaining Gaussian. (orig.)
On the effect of rotation on magnetohydrodynamic turbulence at high magnetic Reynolds number
Favier, Benjamin F N; Cambon, Claude; 10.1080/03091929.2010.544655
2011-01-01
This article is focused on the dynamics of a rotating electrically conducting fluid in a turbulent state. As inside the Earth's core or in various industrial processes, a flow is altered by the presence of both background rotation and a large scale magnetic field. In this context, we present a set of 3D direct numerical simulations of incompressible decaying turbulence. We focus on parameters similar to the ones encountered in geophysical and astrophysical flows, so that the Rossby number is small, the interaction parameter is large, but the Elsasser number, defining the ratio between Coriolis and Lorentz forces, is about unity. These simulations allow to quantify the effect of rotation and thus inertial waves on the growth of magnetic fluctuations due to Alfv\\'en waves. Rotation prevents the occurrence of equipartition between kinetic and magnetic energies, with a reduction of magnetic energy at decreasing Elsasser number {\\Lambda}. It also causes a decrease of energy transfer mediated by cubic correlations....
The amplification of magnetic fields (MFs) in the intracluster medium (ICM) is attributed to turbulent dynamo (TD) action, which is generally derived in the collisional-MHD framework. However, this assumption is poorly justified a priori, since in the ICM the ion mean free path between collisions is of the order of the dynamical scales, thus requiring a collisionless MHD description. The present study uses an anisotropic plasma pressure that brings the plasma within a parametric space where collisionless instabilities take place. In this model, a relaxation term of the pressure anisotropy simulates the feedback of the mirror and firehose instabilities, in consistency with empirical studies. Our three-dimensional numerical simulations of forced transonic turbulence, aiming the modeling of the turbulent ICM, were performed for different initial values of the MF intensity and different relaxation rates of the pressure anisotropy. We found that in the high-β plasma regime corresponding to the ICM conditions, a fast anisotropy relaxation rate gives results that are similar to the collisional-MHD model, as far as the statistical properties of the turbulence are concerned. Also, the TD amplification of seed MFs was found to be similar to the collisional-MHD model. The simulations that do not employ the anisotropy relaxation deviate significantly from the collisional-MHD results and show more power at the small-scale fluctuations of both density and velocity as a result of the action of the instabilities. For these simulations, the large-scale fluctuations in the MF are mostly suppressed and the TD fails in amplifying seed MFs.
Ptuskin, V.S.; /Troitsk, IZMIRAN /Maryland U.; Moskalenko, Igor V.; /Stanford U., HEPL; Jones, F.C.; /NASA, Goddard; Strong, A.W.; /Garching, Max Planck Inst., MPE; Zirakashvili, V.N.; /Troitsk, IZMIRAN /Heidelberg, Max Planck Inst. Astron.
2006-01-17
The physical processes involved in diffusion of Galactic cosmic rays in the interstellar medium are addressed. We study the possibility that the nonlinear MHD cascade sets the power-law spectrum of turbulence which scatters charged energetic particles. We find that the dissipation of waves due to the resonant interaction with cosmic ray particles may terminate the Kraichnan-type cascade below wavelengths 10{sup 13} cm. The effect of this wave dissipation has been incorporated in the GALPROP numerical propagation code in order to asses the impact on measurable astrophysical data. The energy-dependence of the cosmic-ray diffusion coefficient found in the resulting self-consistent model may explain the peaks in the secondary to primary nuclei ratios observed at about 1 GeV/nucleon.
Optical Turbulence Characterization at LAMOST Site: Observations and Models
Liu, L -Y; Yao, Y -Q; Vernin, J; Chadid, M; Wang, H -S; Yin, J; Wang, Y -P
2015-01-01
Atmospheric optical turbulence seriously limits the performance of high angular resolution instruments. An 8-night campaign of measurements was carried out at the LAMOST site in 2011, to characterize the optical turbulence. Two instruments were set up during the campaign: a Differential Image Motion Monitor (DIMM) used to measure the total atmospheric seeing, and a Single Star Scidar (SSS) to measure the vertical profiles of the turbulence C_n^2(h) and the horizontal wind velocity V(h). The optical turbulence parameters are also calculated with the Weather Research and Forecasting (WRF) model coupled with the Trinquet-Vernin model, which describes optical effects of atmospheric turbulence by using the local meteorological parameters. This paper presents assessment of the optical parameters involved in high angular resolution astronomy. Its includes seeing, isoplanatic angle, coherence time, coherence etendue, vertical profiles of optical turbulence intensity _n^2(h)$ and horizontal wind speed V(h). The median...
Generation and characterization of quantum turbulence
Thompson, Kyle J.
Turbulence is a specific type of fluid flow. We all know what it is but it is difficult to define. Turbulence is the rule, rather than the exception and is ubiquitous in nature. It is prevalent in the interstellar dust clouds reaching between the stars and crucial to the plasma dynamics inside them. It can't be ignored in the dynamics of planetary atmospheres, oceans, or cores. It is ubiquitous, from the observable environment down to the internal structure of the human body and every living organism known and unknown. We are composed of and encompassed in an ever changing mess of fluid motion. Nearly the entire universe is in a fluid state, and the majority of fluids are turbulent. Yet there exists no analytic solution to the equations that govern the chaos. This is incredible when pondered; almost everything that can be experienced is governed by the same rules, independent of size, shape or constituents. However, to describe our surroundings we are confined to brute force calculation and experiments. The theories of fluid dynamics are truly theories of everything, and yet, so many properties are still not understood. Presented in this thesis are: an experimental implementation of an existing linear motor to create the first studies of quantum turbulence created by this means in the low temperature limit; identification of the possible shortfalls of this preliminary motor system; the development of a new, more sophisticated, drive apparatus, and the creation, exploration, and development of new and existing sensors for measuring the properties of quantum turbulence. Each of the above parts in this thesis are different prongs of a multifaceted approach underway at the University of Florida to understand quantum turbulence across its entire temperature range. The work here builds on the previous projects in the group, to develop millikelvin turbulence machinery and techniques, but also goes forward in a new direction to explore the higher temperature regime with
The aim of this study was to analyze the feedback process between the magnetic turbulence and the pressure gradients in Large Helical Device (LHD) inward-shifted configurations as well as its role in the transition between the soft-hard magnetohydrodynamic (MHD) regimes for instabilities driven by the mode 1/2 in the middle plasma. In the present paper, we summarize the results of two simulations with different Lundquist numbers, S=2.5×105 and 106, assuming a plasma in the slow reconnection regime. The results for the high Lundquist number simulation show that the magnetic turbulence and the pressure gradient in the middle plasma region of LHD are below the critical value to drive the transition to the hard MHD regime, therefore only relaxations in the soft MHD limit are triggered (1/2 sawtooth-like events) [Phys. Plasmas 19, 082512 (2012)]. In the case of the simulation with low Lundquist number, the system reaches the hard MHD limit and a plasma collapse is observed
M. El-Alaoui; R. L. Richard; Ashour-Abdalla, M.; Walker, R. J.; Goldstein, M. L.
2012-01-01
We report the results of MHD simulations of Earth's magnetosphere for idealized steady solar wind plasma and interplanetary magnetic field (IMF) conditions. The simulations feature purely northward and southward magnetic fields and were designed to study turbulence in the magnetotail plasma sheet. We found that the power spectral densities (PSDs) for both northward and southward IMF had the characteristics of turbulent flow. In both cases, the PSDs showed the three scale ranges expected from ...
Characterizing the mean-field dynamo in turbulent accretion disks
Gressel, Oliver
2015-01-01
The formation and evolution of a wide class of astrophysical objects is governed by turbulent, magnetized accretion disks. Understanding their secular dynamics is of primary importance. Apart from enabling mass accretion via the transport of angular momentum, the turbulence affects the long-term evolution of the embedded magnetic flux, which in turn regulates the efficiency of the transport. In this paper, we take a comprehensive next step towards an effective mean-field model for turbulent astrophysical disks by systematically studying the key properties of magnetorotational turbulence in vertically-stratified, isothermal shearing boxes. This allows us to infer emergent properties of the ensuing chaotic flow as a function of the shear parameter as well as the amount of net-vertical flux. Using the test-field method, we furthermore characterize the mean-field dynamo coefficients that describe the long-term evolution of large-scale fields. We simultaneously infer the vertical shape and the spectral scale depen...
Ravelet, F
2005-09-15
We report experimental studies of the turbulent von Karman flow, inertially stirred between counter-rotating impellers. We first study the flow and its transition from laminar to turbulent regime. We highlight the role of slowly varying large scales, due to the presence of an azimuthal mixing layer. The large scales of this flow can be unstable in turbulent regime. We study the statistics of the transitions between the different mean states. The second part is dedicated to an experiment in liquid sodium, called VKS2. We optimize the time-averaged flow in order to allow kinematic dynamo action. We report the very first results of the experiment, and discuss the role of the large scales temporal non-stationariness. (author)
Characterizing Ocean Turbulence from Argo, Acoustic Doppler, and Simulation Data
McCaffrey, Katherine
Turbulence is inherently chaotic and unsteady, so observing it and modeling it are no easy tasks. The ocean's sheer size makes it even more difficult to observe, and its unpredictable and ever-changing forcings introduce additional complexities. Turbulence in the oceans ranges from basin scale to the scale of the molecular viscosity. The method of energy transfer between scales is, however, an area of active research, so observations of the ocean at all scales are crucial to understanding the basic dynamics of its motions. In this collection of work, I use a variety of datasets to characterize a wide range of scales of turbulence, including observations from multiple instruments and from models with different governing equations. I analyzed the largest scales of the turbulent range using the global salinity data of the Argo profiling float network. Taking advantage of the scattered and discontinuous nature of this dataset, the second-order structure function was calculated down to 2000m depth, and shown to be useful for predicting spectral slopes. Results showed structure function slopes of 2/3 at small scales, and 0 at large scales, which corresponds with spectral slopes of -5/3 at small scales, and -1 at large scales. Using acoustic Doppler velocity measurements, I characterized the meter- to kilometer-scale turbulence at a potential tidal energy site in the Puget Sound, WA. Acoustic Doppler current profiler (ADCP) and acoustic Doppler velocimeter (ADV) observations provided the data for an analysis that includes coherence, anisotropy, and intermittency. In order to more simply describe these features, a parameterization was done with four turbulence metrics, and the anisotropy magnitude, introduced here, was shown to most closely capture the coherent events. Then, using both the NREL TurbSim stochastic turbulence generator and the NCAR large-eddy simulation (LES) model, I calculated turbulence statistics to validate the accuracy of these methods in reproducing
Mean-field magnetohydrodynamics and dynamo theory
Krause, F
2013-01-01
Mean-Field Magnetohydrodynamics and Dynamo Theory provides a systematic introduction to mean-field magnetohydrodynamics and the dynamo theory, along with the results achieved. Topics covered include turbulence and large-scale structures; general properties of the turbulent electromotive force; homogeneity, isotropy, and mirror symmetry of turbulent fields; and turbulent electromotive force in the case of non-vanishing mean flow. The turbulent electromotive force in the case of rotational mean motion is also considered. This book is comprised of 17 chapters and opens with an overview of the gen
Haas, Fernando
2005-01-01
The quantum hydrodynamic model for charged particle systems is extended to the cases of non zero magnetic fields. In this way, quantum corrections to magnetohydrodynamics are obtained starting from the quantum hydrodynamical model with magnetic fields. The quantum magnetohydrodynamics model is analyzed in the infinite conductivity limit. The conditions for equilibrium in ideal quantum magnetohydrodynamics are established. Translationally invariant exact equilibrium solutions are obtained in t...
El-Alaoui, M.; Richard, R. L.; Ashour-Abdalla, M.; Walker, R. J.; Goldstein, M. L.
2012-01-01
We report the results of MHD simulations of Earth's magnetosphere for idealized steady solar wind plasma and interplanetary magnetic field (IMF) conditions. The simulations feature purely northward and southward magnetic fields and were designed to study turbulence in the magnetotail plasma sheet. We found that the power spectral densities (PSDs) for both northward and southward IMF had the characteristics of turbulent flow. In both cases, the PSDs showed the three scale ranges expected from theory: the energy-containing scale, the inertial range, and the dissipative range. The results were generally consistent with in-situ observations and theoretical predictions. While the two cases studied, northward and southward IMF, had some similar characteristics, there were significant differences as well. For southward IMF, localized reconnection was the main energy source for the turbulence. For northward IMF, remnant reconnection contributed to driving the turbulence. Boundary waves may also have contributed. In both cases, the PSD slopes had spatial distributions in the dissipative range that reflected the pattern of resistive dissipation. For southward IMF there was a trend toward steeper slopes in the dissipative range with distance down the tail. For northward IMF there was a marked dusk-dawn asymmetry with steeper slopes on the dusk side of the tail. The inertial scale PSDs had a dusk-dawn symmetry during the northward IMF interval with steeper slopes on the dawn side. This asymmetry was not found in the distribution of inertial range slopes for southward IMF. The inertial range PSD slopes were clustered around values close to the theoretical expectation for both northward and southward IMF. In the dissipative range, however, the slopes were broadly distributed and the median values were significantly different, consistent with a different distribution of resistivity.
M. El-Alaoui
2012-03-01
Full Text Available We report the results of MHD simulations of Earth's magnetosphere for idealized steady solar wind plasma and interplanetary magnetic field (IMF conditions. The simulations feature purely northward and southward magnetic fields and were designed to study turbulence in the magnetotail plasma sheet. We found that the power spectral densities (PSDs for both northward and southward IMF had the characteristics of turbulent flow. In both cases, the PSDs showed the three scale ranges expected from theory: the energy-containing scale, the inertial range, and the dissipative range. The results were generally consistent with in-situ observations and theoretical predictions. While the two cases studied, northward and southward IMF, had some similar characteristics, there were significant differences as well. For southward IMF, localized reconnection was the main energy source for the turbulence. For northward IMF, remnant reconnection contributed to driving the turbulence. Boundary waves may also have contributed. In both cases, the PSD slopes had spatial distributions in the dissipative range that reflected the pattern of resistive dissipation. For southward IMF there was a trend toward steeper slopes in the dissipative range with distance down the tail. For northward IMF there was a marked dusk-dawn asymmetry with steeper slopes on the dusk side of the tail. The inertial scale PSDs had a dusk-dawn symmetry during the northward IMF interval with steeper slopes on the dawn side. This asymmetry was not found in the distribution of inertial range slopes for southward IMF. The inertial range PSD slopes were clustered around values close to the theoretical expectation for both northward and southward IMF. In the dissipative range, however, the slopes were broadly distributed and the median values were significantly different, consistent with a different distribution of resistivity.
Magnetohydrodynamics of Fractal Media
Vasily E. Tarasov
2007-01-01
The fractal distribution of charged particles is considered. An example of this distribution is the charged particles that are distributed over fractal. The fractional integrals are used to describe fractal distribution. These integrals are considered as approximations of integrals on fractals. Typical turbulent media could be of a fractal structure and the corresponding equations should be changed to include the fractal features of the media. The magnetohydrodynamics equations for fractal me...
Hydrocarbon characterization experiments in fully turbulent fires.
Ricks, Allen; Blanchat, Thomas K.
2007-05-01
As the capabilities of numerical simulations increase, decision makers are increasingly relying upon simulations rather than experiments to assess risks across a wide variety of accident scenarios including fires. There are still, however, many aspects of fires that are either not well understood or are difficult to treat from first principles due to the computational expense. For a simulation to be truly predictive and to provide decision makers with information which can be reliably used for risk assessment the remaining physical processes must be studied and suitable models developed for the effects of the physics. The model for the fuel evaporation rate in a liquid fuel pool fire is significant because in well-ventilated fires the evaporation rate largely controls the total heat release rate from the fire. A set of experiments are outlined in this report which will provide data for the development and validation of models for the fuel regression rates in liquid hydrocarbon fuel fires. The experiments will be performed on fires in the fully turbulent scale range (> 1 m diameter) and with a number of hydrocarbon fuels ranging from lightly sooting to heavily sooting. The importance of spectral absorption in the liquid fuels and the vapor dome above the pool will be investigated and the total heat flux to the pool surface will be measured. The importance of convection within the liquid fuel will be assessed by restricting large scale liquid motion in some tests. These data sets will provide a sound, experimentally proven basis for assessing how much of the liquid fuel needs to be modeled to enable a predictive simulation of a fuel fire given the couplings between evaporation of fuel from the pool and the heat release from the fire which drives the evaporation.
Characterizing Electron Temperature Gradient Turbulence Via Numerical Simulation
Nevins, W M; Candy, J; Cowley, S; Dannert, T; Dimits, A; Dorland, W; Estrada-Mila, C; Hammett, G W; Jenko, F; Pueschel, M J; Shumaker, D E
2006-05-22
Numerical simulations of electron temperature gradient (ETG) turbulence are presented which characterize the ETG fluctuation spectrum, establish limits to the validity of the adiabatic ion model often employed in studying ETG turbulence, and support the tentative conclusion that plasmaoperating regimes exist in which ETG turbulence produces sufficient electron heat transport to be experimentally relevant. We resolve prior controversies regarding simulation techniques and convergence by benchmarking simulations of ETG turbulence from four microturbulence codes, demonstrating agreement on the electron heat flux, correlation functions, fluctuation intensity, and rms flow shear at fixed simulation cross section and resolution in the plane perpendicular to the magnetic field. Excellent convergence of both continuum and particle-in-cell codes with time step and velocity-space resolution is demonstrated, while numerical issues relating to perpendicular (to the magnetic field) simulation dimensions and resolution are discussed. A parameter scan in the magnetic shear, s, demonstrates that the adiabatic ion model is valid at small values of s (s < 0.4 for the parameters used in this scan) but breaks down at higher magnetic shear. A proper treatment employing gyrokinetic ions reveals a steady increase in the electron heat transport with increasing magnetic shear, reaching electron heat transport rates consistent with analyses of experimental tokamak discharges.
Characterization of Fibre Channel over Highly Turbulent Optical Wireless Links
Johnson, G W; Henderer, B D; Wilburn, J W; Ruggiero, A J
2003-07-28
We report on the performance characterization and issues associated with using Fibre Channel (FC) over a highly turbulent free-space optical (FSO) link. Fibre Channel is a storage area network standard that provides high throughput with low overhead. Extending FC to FSO links would simplify data transfer from existing high-bandwidth sensors such as synthetic aperture radars and hyperspectral imagers. We measured the behavior of FC protocol at 1 Gbps in the presence of synthetic link dropouts that are typical of turbulent FSO links. Results show that an average bit error rate of less than 2 x 10{sup -8} is mandatory for adequate throughput. More importantly, 10 ns dropouts at a 2 Hz rate were sufficient to cause long (25 s) timeouts in the data transfer. Although no data was lost, this behavior is likely to be objectionable for most applications. Prospects for improvements in hardware and software will be discussed.
Optical measurements for turbulence characterization in RFX-mod edge
Agostini, M.; Scarin, P.; Cavazzana, R.; Fassina, A.; Alfier, A.; Cervaro, V. [Consorzio RFX, Associazione EURATOM-ENEA sulla Fusione, 35127 Padova (Italy)
2010-10-15
The edge of the reversed field pinch experiment RFX-mod (R=2 m, a=0.46 m) is characterized by a complex magnetic topology and strong pressure gradients. A set of fast spectroscopic diagnostics has been designed to characterize both main plasma edge parameters and turbulence, including coherent structures (blobs) and radial profiles of the thermodynamic variables. Measurements of the edge radial profiles of electron density (n{sub e}), temperature (T{sub e}), and pressure (P{sub e}) are obtained by a thermal helium beam diagnostic, which locally measures the ratios of the spectral lines emitted by a neutral He cloud puffed into hydrogen plasma discharges. A gas puff imaging system allows the high frequency (1 MHz bandwidth) measurement of the fluctuations of He I emissivity in the same positions. The system measures both the intermittent behavior of the edge turbulence and n{sub e}, T{sub e}, and P{sub e} profiles of the corresponding coherent structures, with a resolution down to 6 {mu}s. The optical diagnostics are complemented by a set of triaxial magnetic coils to simultaneously characterize the edge magnetic field.
Bailly, Christophe
2015-01-01
This book covers the major problems of turbulence and turbulent processes, including physical phenomena, their modeling and their simulation. After a general introduction in Chapter 1 illustrating many aspects dealing with turbulent flows, averaged equations and kinetic energy budgets are provided in Chapter 2. The concept of turbulent viscosity as a closure of the Reynolds stress is also introduced. Wall-bounded flows are presented in Chapter 3, and aspects specific to boundary layers and channel or pipe flows are also pointed out. Free shear flows, namely free jets and wakes, are considered in Chapter 4. Chapter 5 deals with vortex dynamics. Homogeneous turbulence, isotropy, and dynamics of isotropic turbulence are presented in Chapters 6 and 7. Turbulence is then described both in the physical space and in the wave number space. Time dependent numerical simulations are presented in Chapter 8, where an introduction to large eddy simulation is offered. The last three chapters of the book summarize remarka...
Magnetohydrodynamic cosmologies
We analyse a class of cosmological models in magnetohydrodynamic regime extending and completing the results of a previous paper. The material content of the models is a perfect fluid plus electromagnetic fields. The fluid is neutral in average but admits an electrical current which satisfies Ohm's law. All models fulfil the physical requirements of near equilibrium thermodynamics and can be favourably used as a more realistic description of the interior of a collapsing star in a magnetohydrodynamic regime with or without a magnetic field. (author)
Z. Lin; R.E. Waltz
2007-01-01
@@ Turbulent transport driven by plasma pressure gradients [Tangl978] is one of the most important scientific challenges in burning plasma experiments since the balance between turbulent transport and the self-heating by the fusion products (a-particles) determines the performance of a fusion reactor like ITER.
Lingam, M; Morrison, Philip; E. Tassi
2015-01-01
A version of extended magnetohydrodynamics (MHD) that incorporates electron inertia is obtained by constructing an action principle. Unlike MHD which freezes in magnetic flux, the present theory freezes in an alternative flux related to the electron canonical momentum. The associated Hamiltonian formulation is derived and reduced models that have previously been used to describe collisionless reconnection are obtained.
A new neutrino magnetohydrodynamics (NMHD) model is formulated, where the effects of the charged weak current on the electron-ion magnetohydrodynamic fluid are taken into account. The model incorporates in a systematic way the role of the Fermi neutrino weak force in magnetized plasmas. A fast neutrino-driven short wavelengths instability associated with the magnetosonic wave is derived. Such an instability should play a central role in strongly magnetized plasma as occurs in supernovae, where dense neutrino beams also exist. In addition, in the case of nonlinear or high frequency waves, the neutrino coupling is shown to be responsible for breaking the frozen-in magnetic field lines condition even in infinite conductivity plasmas. Simplified and ideal NMHD assumptions were adopted and analyzed in detail
Haas, Fernando; Pascoal, Kellen Alves [Instituto de Física, Universidade Federal do Rio Grande do Sul, Av. Bento Gonçalves 9500, 91501-970 Porto Alegre, RS (Brazil); Mendonça, José Tito [IPFN, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal and Instituto de Física, Universidade de São Paulo, 05508-090 São Paulo, SP (Brazil)
2016-01-15
A new neutrino magnetohydrodynamics (NMHD) model is formulated, where the effects of the charged weak current on the electron-ion magnetohydrodynamic fluid are taken into account. The model incorporates in a systematic way the role of the Fermi neutrino weak force in magnetized plasmas. A fast neutrino-driven short wavelengths instability associated with the magnetosonic wave is derived. Such an instability should play a central role in strongly magnetized plasma as occurs in supernovae, where dense neutrino beams also exist. In addition, in the case of nonlinear or high frequency waves, the neutrino coupling is shown to be responsible for breaking the frozen-in magnetic field lines condition even in infinite conductivity plasmas. Simplified and ideal NMHD assumptions were adopted and analyzed in detail.
Over the past 2 decades, ideal magnetohydrodynamics (MHD), has developed into a relatively mature theory within the field of plasma physics. MHD represents the simplest, self-consistent model describing the macroscopic equilibrium and stability properties of plasma. Thus, the role of ideal MHD in magnetic fusion is the discovery of magnetic geometries that possess attractive equilibrium and stability properties for fusion reactors. This book is based on the author's two-semester sequence in MHD theory offered at MIT. The first semester is devoted to ideal magnetohydrodynamics and is largely involved in studying equilibrium and fast, often catastrophic, instabilities. The second semester treats slower, nonideal effects (such as resistivity, pressure anisotropy, and finite Larmor radius effects) and their influence on plasmas that are stable or at most, weakly unstable against the ideal modes. This textbook is concerned with the first semester's material. The author's goal is to provide an in-depth introduction to the subject of ideal magnetohydrodynamics. The style has a theoretical emphasis with considerable discussion of the application of theory to experiments and generous examples of the experiments themselves. The book has been written to satisfy the needs of the pedagogical requirements of the student. Accordingly, there are many descriptive explanations and virtually every equation is derived - either in detail, in outline form, or in the Appendices. Chapters are basically divided between equilibrium: general considerations, one-dimensional configurations, two-dimensional configurations, three-dimensional configurations; and stability; general considerations, one-dimensional configurations, and multidimensional configurations
Topological characterization of flow structures in resistive pressure-gradient-driven turbulence
Carreras, B. A.; Llerena, I.; Garcia, L.; Calvo, I.
2008-12-01
Visualization of turbulent flows is a powerful tool to help understand the turbulence dynamics and induced transport. However, it does not provide a quantitative description of the observed structures. In this paper, an approach to characterize quantitatively the topology of the flows is given. The technique, which can be applied to any type of turbulence dynamics, is illustrated through the example of resistive ballooning instabilities.
Topological characterization of flow structures in resistive pressure-gradient-driven turbulence
Carreras, Benjamin A; Garcia, Luis; Calvo, Ivan
2008-01-01
Visualization of turbulent flows is a powerful tool to help understand the turbulence dynamics and induced transport. However, it does not provide a quantitative description of the observed structures. In this paper, an approach to characterize quantitatively the topology of the flows is given. The technique, which can be applied to any type of turbulence dynamics, is illustrated through the example of resistive ballooning instabilities.
Ben-Nun, R.; Sheintuch, M.; Kysela, B.; Konfršt, J. (Jiří); I. Fořt
2015-01-01
We characterize turbulent flow from radial impellers in stirred tanks, formulating semi-analytical predictions based on conventional turbulent-jet theory and the general framework of scalar dispersion in turbulent shear flows. This work introduces a novel formulation of the radial-discharge flow as two separate spatial regions along the radial axis: the zone of flow establishment (ZFE) and the zone of established flow (ZEF). We discover and analyze resulting critical improvements in the predi...
Dissipation and reconnection in boundary-driven reduced magnetohydrodynamics
We study the statistics of coherent current sheets, the population of X-type critical points, and reconnection rates in a coronal loop geometry, via numerical simulations of reduced magnetohydrodynamic turbulence. Current sheets and sites of reconnection (magnetic X-points) are identified in two-dimensional planes of the three-dimensional simulation domain. The geometry of the identified current sheets—including area, length, and width—and the magnetic dissipation occurring in the current sheets are statistically characterized. We also examine the role of magnetic reconnection, by computing the reconnection rates at the identified X-points and investigating their association with current sheets.
Magnetohydrodynamics of accretion disks
The thesis consists of an introduction and summary, and five research papers. The introduction and summary provides the background in accretion disk physics and magnetohydrodynamics. The research papers describe numerical studies of magnetohydrodynamical processes in accretion disks. Paper 1 is a one-dimensional study of the effect of magnetic buoyancy on a flux tube in an accretion disk. The stabilizing influence of an accretion disk corona on the flux tube is demonstrated. Paper 2-4 present numerical simulations of mean-field dynamos in accretion disks. Paper 11 verifies the correctness of the numerical code by comparing linear models to previous work by other groups. The results are also extended to somewhat modified disk models. A transition from an oscillatory mode of negative parity for thick disks to a steady mode of even parity for thin disks is found. Preliminary results for nonlinear dynamos at very high dynamo numbers are also presented. Paper 3 describes the bifurcation behaviour of the nonlinear dynamos. For positive dynamo numbers it is found that the initial steady solution is replaced by an oscillatory solution of odd parity. For negative dynamo numbers the solution becomes chaotic at sufficiently high dynamo numbers. Paper 4 continues the studies of nonlinear dynamos, and it is demonstrated that a chaotic solution appears even for positive dynamo numbers, but that it returns to a steady solution of mixed parity at very high dynamo numbers. Paper 5 describes a first attempt at simulating the small-scale turbulence of an accretion disk in three dimensions. There is only find cases of decaying turbulence, but this is rather due to limitations of the simulations than that turbulence is really absent in accretion disks
A version of extended magnetohydrodynamics (MHD) that incorporates electron inertia is obtained by constructing an action principle. Unlike MHD which freezes in magnetic flux, the present theory freezes in an alternative flux related to the electron canonical momentum. The associated Hamiltonian formulation is derived and reduced models that have previously been used to describe collisionless reconnection are obtained. - Highlights: • A rare, and possibly novel, version of MHD presented, and similarities to MHD highlighted via an appropriate frozen flux constraint. • Hamiltonian and action principle formalisms for this model, inertial MHD, are presented. • Shown to reproduce the Ottaviani–Porcelli model of reconnection as a special case
Characterization of Gigabit Ethernet Over Highly Turbulent Optical Wireless Links
Johnson, G W; Cornish, J P; Wilburn, J W; Young, R A; Ruggiero, A J
2002-07-01
We report on the performance characterization and issues associated with using Gigabit Ethernet (GigE) over a highly turbulent (C{sub n}{sup 2} > 10{sup -12}) 1.3 km air-optic lasercom links. Commercial GigE hardware is a cost-effective and scalable physical layer standard that can be applied to air-optic communications. We demonstrate a simple GigE hardware interface to a singlemode fiber-coupled, 1550 nm, WDM air-optic transceiver. TCPAP serves as a robust and universal foundation protocol that has some tolerance of data loss due to atmospheric fading. Challenges include establishing and maintaining a connection with acceptable throughput under poor propagation conditions. The most useful link performance diagnostic is shown to be scintillation index, where a value of 0.2 is the maximum permissible for adequate GigE throughput. Maximum GigE throughput observed was 49.7% of that obtained with a fiber jumper when scintillation index is 0.1. Shortcomings in conventional measurements such as bit error rate are apparent. Prospects for forward mor correction and other link enhancements will be discussed.
Logarithmic Conformal Field Theory Solutions of Two Dimensional Magnetohydrodynamics
Skoulakis, Spyros; Thomas, Steven
1998-01-01
We consider the application of logarithmic conformal field theory in finding solutions to the turbulent phases of 2-dimensional models of magnetohydrodynamics. These arise upon dimensional reduction of standard (infinite conductivity) 3-dimensional magnetohydrodynamics, after taking various simplifying limits. We show that solutions of the corresponding Hopf equations and higher order integrals of motion can be found within the solutions of ordinary turbulence proposed by Flohr, based on the ...
Statistical characterization of turbulence in the boundary plasma of EAST
Yan, Ning; Nielsen, Anders Henry; Xu, G.S.;
2013-01-01
In Ohmic heated low confinement mode (L-mode) discharges, the intermittent statistical characteristics of turbulent fluctuations have been investigated in the edge and the scrape-off layer (SOL) plasma on EAST (the experimental advanced superconducting tokamak) by fast reciprocating Langmuir probe...
Maron, Jason L.; Chandran, Benjamin D. G.; Blackman, Eric G.
2003-01-01
We investigate field-line separation in strong MHD turbulence using direct numerical simulations. We find that in the static-magnetic-field approximation the thermal conductivity in galaxy clusters is reduced by a factor of about 50 relative to the Spitzer thermal conductivity of a non-magnetized plasma. This value is too small for heat conduction to balance radiative cooling in clusters.
Experimental characterization of turbulent subsonic transitional-open cavity flow
Rokita, T.; Elimelech, Y.; Arieli, R.; Levy, Y.; Greenberg, J. B.
2016-04-01
Turbulent subsonic "transitional-open" cavity flow was investigated by wind-tunnel tests. The investigated cavity configuration had a length-to-depth ratio of 6.25 and a width-to-depth ratio of 2. The cavity was exposed to a free-stream Mach number of 0.40 and a Reynolds number (based on cavity depth) of 1.6× 10^6, with a turbulent incoming boundary layer. Measurements of velocity and wall pressures were taken simultaneously, which enabled the analysis of velocity-pressure cross-correlations. Special attention is paid to the shear layer that develops over the cavity and an emphasis is placed on the analysis of its characteristics and its stability. Application of linear hydrodynamic stability theory, together with examining velocity-pressure cross correlations, revealed that the behavior of the cavity shear layer is analogous to a free shear layer, approximately up to mid-length of the cavity, where further downstream nonlinear interactions occur.
Edge turbulence characterization in RFX-mod with optical diagnostics
Agostini, M; Scarin, P; Cavazzana, R; Sattin, F; Serianni, G; Spolaore, M; Vianello, N [Consorzio RFX, Associazione EURATOM-ENEA sulla fusione, C.so Stati Uniti 4 35127, Padova (Italy)
2009-10-15
In the edge region of fusion devices coherent structures (blobs) are commonly observed. In order to investigate their evolution and the possible driving mechanisms, we studied the relationship between the size of the coherent structures and the local characteristic length L{sub p} of the edge electron pressure profile P{sub e}(r), together with the effect of magnetic fluctuations on nablaP{sub e}. The results show an increase in the blob size with both L{sub p} and magnetic fluctuations, suggesting an interplay between electron pressure gradient, magnetic fluctuations and edge turbulence generation. Measurements of the spatial scales of the edge turbulence structures and the corresponding toroidal phase velocity are obtained in the reversed field pinch RFX-mod by means of a gas puffing imaging diagnostic, while the local edge profiles of electron density n{sub e}(r) and temperature T{sub e}(r) and therefore P{sub e}(r) have been measured via the line intensity ratio technique of helium emission lines.
Filamentary magnetohydrodynamic plasmas
A filamentary construct of magnetohydrodynamical plasma dynamics, based on the Elsasser variables was developed. This approach is modeled after discrete vortex models of hydrodynamical turbulence, which cannot be expected in general to produce results identical to ones based on a Fourier decomposition of the fields. In a highly intermittent plasma, the induction force is small compared to the convective motion, and when this force is neglected. the plasma vortex system is described by a Hamiltonian. For a system with many such vortices we present a statistical treatment of a collection of discrete current-vorticity concentrations. Canonical and microcanonical statistical calculations show that both the vorticity and the current spectra are peaked at long wavelengths, and the expected states revert to known hydrodynamical states as the magnetic field vanishes. These results differ from previous Fourier-based statistical theories. but it is found that when the filament calculation is expanded to include the inductive force, the results approach the Fourier equilibria in the low-temperature limit, and the previous Hamiltonian plasma vortex results in the high-temperature limit. Numerical simulations of a large number of filaments are carried out and support the theory. A three-dimensional vortex model is outlined as well, which is also Hamiltonian when the inductive force is neglected
Attractors of magnetohydrodynamic flows in an Alfvenic state
Nunez, Manuel; Sanz, Javier [Departamento de Analisis Matematico, Universidad de Valladolid, Valladolid (Spain)
1999-08-13
We present a simplified form of the magnetohydrodynamic system which describes the evolution of a plasma where the small-scale velocity and magnetic field are aligned in the form of Alfven waves, such as happens in several turbulent situations. Bounds on the dimension of the global attractor are found, and are shown to be an improvement of the standard ones for the full magnetohydrodynamic equations. (author)
Characterization of an incipiently separated shock wave/turbulent boundary layer interaction
Schreyer, A.-M.; Dussauge, J.-P.; Krämer, E.
2016-05-01
The turbulence structure in a shock wave/turbulent boundary layer interaction at incipient separation was investigated in order to get insight into turbulence generation and amplification mechanisms in such flow fields. The flow along a two-dimensional 11.5° compression corner was studied experimentally at a Mach number of M=2.53 and with a momentum-thickness Reynolds number of Re_{θ }=5370 . From hot-wire boundary layer traverses and surface heat-flux density fluctuation measurements with the fast-response atomic layer thermopile, the turbulence structure and amplification was described. Space-time correlations of the mass-flux fluctuations across the boundary layer and the surface heat-flux density fluctuations were measured to further characterize the development of the turbulence structure across the interaction. The large-scale boundary layer structures are concealed by shock-related effects in the strongly disturbed shock-foot region. Shortly downstream, however, large-scale structures dominate the signal again, just as in the incoming flow. A mechanism explaining this behavior is suggested.
Turbulent transport and dynamo in sheared MHD turbulence with a non-uniform magnetic field
Leprovost, Nicolas; Kim, Eun-Jin
2009-01-01
We investigate three-dimensional magnetohydrodynamics turbulence in the presence of velocity and magnetic shear (i.e., with both a large-scale shear flow and a nonuniform magnetic field). By assuming a turbulence driven by an external forcing with both helical and nonhelical spectra, we investigate the combined effect of these two shears on turbulence intensity and turbulent transport represented by turbulent diffusivities (turbulent viscosity, α and β effect) in Reynolds-averaged equations. ...
For the resistive pressure-gradient-driven turbulence model, the transition from laminar regime to fully developed turbulence is not simple and goes through several phases. For low values of the plasma parameter β, a single quasicoherent structure forms. As β increases, several of these structures may emerge and in turn take the dominant role. Finally, at high β, fully developed turbulence with a broad spectrum is established. A suitable characterization of this transition can be given in terms of topological properties of the flow. Here, we analyze these properties that provide an understanding of the turbulence-induced transport and give a measure of the breaking of the homogeneity of the turbulence. To this end, an approach is developed that allows discriminating between topological properties of plasma turbulence flows that are relevant to the transport dynamics and the ones that are not. This is done using computational homology tools and leads to a faster convergence of numerical results for a fixed level of resolution than previously presented in Phys. Rev. E 78, 066402 (2008).
Garcia, L.; Carreras, B. A.; Llerena, I.; Calvo, I.
2009-10-01
For the resistive pressure-gradient-driven turbulence model, the transition from laminar regime to fully developed turbulence is not simple and goes through several phases. For low values of the plasma parameter β , a single quasicoherent structure forms. As β increases, several of these structures may emerge and in turn take the dominant role. Finally, at high β , fully developed turbulence with a broad spectrum is established. A suitable characterization of this transition can be given in terms of topological properties of the flow. Here, we analyze these properties that provide an understanding of the turbulence-induced transport and give a measure of the breaking of the homogeneity of the turbulence. To this end, an approach is developed that allows discriminating between topological properties of plasma turbulence flows that are relevant to the transport dynamics and the ones that are not. This is done using computational homology tools and leads to a faster convergence of numerical results for a fixed level of resolution than previously presented in Phys. Rev. E 78, 066402 (2008).
Structural characterization of wind-sheared turbulent flow using self-organized mapping
Scott, Nicholas V.; Handler, Robert A.
2016-05-01
A nonlinear cluster analysis algorithm is used to characterize the spatial structure of a wind-sheared turbulent flow obtained from the direct numerical simulation (DNS) of the three-dimensional temperature and momentum fields. The application of self-organizing mapping to DNS data for data reduction is utilized because of the dimensional similitude in structure between DNS data and remotely sensed hyperspectral and multispectral data where the technique has been used extensively. For the three Reynolds numbers of 150, 180, and 220 used in the DNS, self-organized mapping is successful in the extraction of boundary layer streaky structures from the turbulent temperature and momentum fields. In addition, it preserves the cross-wind scale structure of the streaks exhibited in both fields which loosely scale with the inverse of the Reynolds number. Self-organizing mapping of the along wind component of the helicity density shows a layer of the turbulence field which is spotty suggesting significant direct coupling between the large and small-scale turbulent structures. The spatial correlation of the temperature and momentum fields allows for the possibility of the remote extrapolation of the momentum structure from thermal structure.
Collisionless magnetohydrodynamics with gyrokinetic effects
Anisotropic magnetohydrodynamics equations, which also capture the dynamics of quasi-transverse small scales obeying the gyrokinetic ordering, are derived using fourth-rank moment closures, based on a refined description of linear Landau damping and finite Larmor radius (FLR) corrections. This 'FLR-Landau fluid model' reproduces the dispersion relation of low-frequency waves, up to scales that, in the case of quasi-transverse kinetic Alfven waves, can be much smaller than the ion gyroradius. The mirror instability, which requires temperature anisotropy, is also captured, together with its quenching at small scales. This model that accurately reproduces the collisionless dissipation of low-frequency modes, should provide an efficient tool to simulate mesoscale turbulence in a magnetized collisionless plasma
Characterization of Turbulent Open Channel Flow in a Full-Scale Spiral Corrugated Culvert
Guensch, Greg R.; Richmond, Marshall C.; Tritico, Hans; Pearson, Walter H.
2004-02-26
A micro-Acoustic Doppler Velocimeter (ADV) was used to characterize the three-dimensional velocity and turbulence characteristics in a full-scale culvert with spiral corrugations. The culvert was set up in a test bed constructed to examine juvenile salmon passage success in various culvert types. The test culvert was 12.2 m long and 1.83m in diameter, and set at a 1.14% slope. The corrugations were 2.54 cm deep by 7.62 cm peak to peak with a 5° right-handed pitch. Tailwater elevation was adjustable with a stop-log system and was set slightly above the water surface level at the culvert exit. Cross-sectional grids of ADV measurements were taken at discharges of 28, 57, 113, 227, and 453 lps at 9 locations within the culvert barrel and just inside the headwater and tailwater tanks. Results indicated that the spiral corrugations induced asymmetries in the velocity and turbulence distributions. These asymmetries caused a Reduced Velocity Zone (RVZ) on the right side of the culvert as seen looking upstream, which could aid small fish during upstream passage. Velocity and turbulence levels in the RVZ were found to be less than in mid channel or on the left side of the culvert, and the difference became greater at increased flow rates. Lateral and vertical velocity components were very small relative to the axial component, while lateral and vertical turbulence intensities were comparable to the axial component. Inlet loss coefficients were calculated as well and ranged from 0.32 to 0.42. Relationships between the average velocity and the velocity and turbulence intensity in the RVZ were developed, which may be useful for evaluating whether the barrel of a culvert is passable for juvenile fish.
Richmond, M. C.; Thomson, J. M.; Durgesh, V.; Polagye, B. L.
2011-12-01
Field measurements are essential for developing an improved understanding of turbulent inflow conditions that affect the design and operation of marine and hydrokinetic (MHK) devices. The Marrowstone Island site in Puget Sound, Washington State is a potential location for installing MHK devices, as it experiences strong tides and associated currents. Here, field measurements from Nodule Point on the eastern side of Marrowstone Island are used to characterize the turbulence in terms of velocity variance as a function of length and time scales. The field measurements were performed using Acoustic Doppler Velocimetry (ADV) and Acoustic Doppler Current Profiler (ADCP) instruments. Both were deployed on a bottom-mounted tripod at the site by the Applied Physics Lab at the University of Washington (APL-UW). The ADV acquired single point, temporally resolved velocity data from 17-21 Feb 2011, at a height of 4.6 m above the seabed at a sampling frequency of 32 Hz. The ADCP measured the velocity profile over the water column from a height of 2.6 m above the seabed up to the sea-surface in 36 bins, with each bin of 0.5 m size. The ADCP acquired data from 11-27 Feb 2011 at a sampling frequency of 2 Hz. Analysis of the ADV measurements shows distinct dynamic regions by scale: anisotropic eddies at large scales, an isotropic turbulent cascade (-5/3 slope in frequency spectra) at mesoscales, and contamination by Doppler noise at small scales. While Doppler noise is an order of magnitude greater for the ADCP measurements, the turbulence bulk statistics are consistent between the two instruments. There are significant variations in turbulence statistics with stage of the tidal currents (i.e., from slack to non-slack tidal conditions), however an average turbulent intensity of 10% is a robust, canonical value for this site. The ADCP velocity profiles are useful in quantifying the variability in velocity along the water column, and the ensemble averaged velocity profiles may be
Experiments in Magnetohydrodynamics
Rayner, J. P.
1970-01-01
Describes three student experiments in magnetohydrodynamics (MHD). In these experiments, it was found that the electrical conductivity of the local water supply was sufficient to demonstrate effectively some of the features of MHD flowmeters, generators, and pumps. (LC)
Parallel Simulations in Turbulent MHD
The large-scale dynamics of plasma flows can often be described within a fluidistic approximation known as one-fluid magnetohydrodynamics. Complex flows such as those corresponding to turbulent regimes are ubiquitous in laboratory plasmas and in astrophysics, because of their typically very large Reynolds numbers. Numerical simulations have become a powerful tool for the study of complex plasma flows in recent years. The aim of the present paper is to introduce the reader to some of the standard numerical approximations used for the integration of the magnetohydrodynamic equations. In particular, we focus on pseudo-spectral methods and on how to develop parallel codes to speed up large Reynolds number simulations. We show the results arising from numerical simulations of astrophysical interest such as the development of turbulent flows in reduced magnetohydrodynamics and the generation of magnetic fields by dynamo mechanisms in three dimensional magnetohydrodynamics
Parallel Simulations in Turbulent MHD
Gomez, Daniel O. [C. Universitaria, Buenos Aires (Argentina). Dept. of Physics, Pabellon I; Mininni, Pablo D. [National Center for Atmospheric Research, Boulder, CO (United States). Advanced Study Program; Dmitruk, Pablo [Univ. of Delaware, Newark (United States). Bartol Research Inst.
2005-04-01
The large-scale dynamics of plasma flows can often be described within a fluidistic approximation known as one-fluid magnetohydrodynamics. Complex flows such as those corresponding to turbulent regimes are ubiquitous in laboratory plasmas and in astrophysics, because of their typically very large Reynolds numbers. Numerical simulations have become a powerful tool for the study of complex plasma flows in recent years. The aim of the present paper is to introduce the reader to some of the standard numerical approximations used for the integration of the magnetohydrodynamic equations. In particular, we focus on pseudo-spectral methods and on how to develop parallel codes to speed up large Reynolds number simulations. We show the results arising from numerical simulations of astrophysical interest such as the development of turbulent flows in reduced magnetohydrodynamics and the generation of magnetic fields by dynamo mechanisms in three dimensional magnetohydrodynamics.
Gradient Particle Magnetohydrodynamics
Maron, Jason L.; Howes, Gregory G.
2001-01-01
We introduce Gradient Particle Magnetohydrodynamics (GPM), a new Lagrangian method for magnetohydrodynamics based on gradients corrected for the locally disordered particle distribution. The development of a numerical code for MHD simulation using the GPM algorithm is outlined. Validation tests simulating linear and nonlinear sound waves, linear MHD waves, advection of magnetic fields in a magnetized vortex, hydrodynamical shocks, and three-dimensional collapse are presented, demonstrating th...
Roberts, J. Brent; Robertson, Franklin R.; Clayson, Carol Anne; Bosilovich, Michael G.
2012-01-01
Turbulent fluxes of heat and moisture across the atmosphere-ocean interface are fundamental components of the Earth s energy and water balance. Characterizing both the spatiotemporal variability and the fidelity of these exchanges of heat and moisture is critical to understanding the global water and energy cycle variations, quantifying atmosphere-ocean feedbacks, and improving model predictability. This study examines the veracity of the recently completed NASA Modern-Era Retrospective analysis for Research and Applications (MERRA) product with respect to its representation of the surface turbulent heat fluxes. A validation of MERRA turbulent heat fluxes and near-surface bulk variables at local, high-resolution space and time scales is achieved by making comparisons to a large suite of direct observations. Both in situ and satellite-observed gridded surface heat flux estimates are employed to investigate the spatial and temporal variability of the surface fluxes with respect to their annual mean climatologies, their seasonal covariability of near-surface bulk parameters, and their representation of extremes. The impact of data assimilation on the near-surface parameters is assessed through evaluation of incremental analysis update tendencies produced by the assimilation procedure. It is found that MERRA turbulent surface heat fluxes are relatively accurate for typical conditions but have systematically weak vertical gradients in moisture and temperature and have a weaker covariability between the near-surface gradients and wind speed than found in observations. This results in an underestimate of the surface latent and sensible heat fluxes over the western boundary current and storm track regions. The assimilation of observations mostly acts to bring MERRA closer to observational products by increasing moisture and temperature near the surface and decreasing the near-surface wind speeds. The major patterns of spatial and temporal variability of the turbulent heat
Hamiltonian formulation of reduced magnetohydrodynamics
Reduced magnetohydrodynamics (RMHD) has become a principal tool for understanding nonlinear processes, including disruptions, in tokamak plasmas. Although analytical studies of RMHD turbulence have been useful, the model's impressive ability to simulate tokamak fluid behavior has been revealed primarily by numerical solution. The present work describes a new analytical approach, not restricted to turbulent regimes, based on Hamiltonian field theory. It is shown that the nonlinear (ideal) RMHD system, in both its high-beta and low-beta versions, can be expressed in Hanmiltonian form. Thus a Poisson bracket, [ , ], is constructed such that each RMHD field quantitity, xi/sub i/, evolves according to xi/sub i/ = [xi/sub i/,H], where H is the total field energy. The new formulation makes RMHD accessible to the methodology of Hamiltonian mechanics; it has lead, in particular, to the recognition of new RMHD invariants and even exact, nonlinear RMHD solutions. A canonical version of the Poisson bracket, which requires the introduction of additional fields, leads to a nonlinear variational principle for time-dependent RMHD
Lattice Boltzmann approaches to magnetohydrodynamics and electromagnetism
Dellar, Paul
2010-03-01
J u B E g We present a lattice Boltzmann approach for magnetohydrodynamics and electromagnetism that expresses the magnetic field using a discrete set of vector distribution functions i. The i were first postulated to evolve according to a vector Boltzmann equation of the form ti+ ξi.∇i= - 1τ ( i- i^(0) ), where the ξi are a discrete set of velocities. The right hand side relaxes the i towards some specified functions i^(0) of the fluid velocity , and of the macroscopic magnetic field given by = ∑ii. Slowly varying solutions obey the equations of resistive magnetohydrodynamics. This lattice Boltzmann formulation has been used in large-scale (up to 1800^3 resolution) simulations of magnetohydrodynamic turbulence. However, this is only the simplest form of Ohm's law. We may simulate more realistic extended forms of Ohm's law using more complex collision operators. A current-dependent relaxation time yields a current-dependent resistivity η(|∇x|), as used to model ``anomalous'' resistivity created by small-scale plasma processes. Using a hydrodynamic matrix collision operator that depends upon the magnetic field , we may simulate Braginskii's magnetohydrodynamics, in which the viscosity for strains parallel to the magnetic field lines is much larger than the viscosity for strains in perpendicular directions. Changing the collision operator again, from the above vector Boltzmann equation we may derive the full set of Maxwell's equations, including the displacement current, and Ohm's law, - 1c^2 tE+ ∇x= μo,= σ( E + x). The original lattice Boltzmann scheme was designed to reproduce resistive magnetohydrodynamics in the non-relativistic limit. However, the kinetic formulation requires a system of first order partial differential equations with collision terms. This system coincides with the full set of Maxwell's equations and Ohm's law, so we capture a much wider range of electromagnetic phenomena, including electromagnetic waves.
BOOK REVIEW: Nonlinear Magnetohydrodynamics
Shafranov, V.
1998-08-01
Nonlinear magnetohydrodynamics by Dieter Biskamp is a thorough introduction to the physics of the most impressive non-linear phenomena that occur in conducting magnetoplasmas. The basic systems, in which non-trivial dynamic processes are observed, accompanied by changes of geometry of the magnetic field and the effects of energy transformation (magnetic energy into kinetic energy or the opposite effect in magnetic dynamos), are the plasma magnetic confinement systems for nuclear fusion and space plasmas, mainly the solar plasma. A significant number of the examples of the dynamic processes considered are taken from laboratory plasmas, for which an experimental check of the theory is possible. Therefore, though the book is intended for researchers and students interested in both laboratory, including nuclear fusion, and astrophysical plasmas, it is most probably closer to the first category of reader. In the Introduction the author notes that unlike the hydrodynamics of non-conducting fluids, where the phenomena caused by rapid fluid motions are the most interesting, for plasmas in a strong magnetic field the quasi-static configurations inside which the local dynamic processes occur are often the most important. Therefore, the reader will also find in this book rather traditional material on the theory of plasma equilibrium and stability in magnetic fields. In addition, it is notable that, as opposed to a linear theory, the non-linear theory, as a rule, cannot give quite definite explanations or predictions of phenomena, and consequently there are in the book many results obtained by consideration of numerical models with the use of supercomputers. The treatment of non-linear dynamics is preceded by Chapters 2 to 4, in which the basics of MHD theory are presented with an emphasis on the role of integral invariants of the magnetic helicity type, a derivation of the reduced MHD equations is given, together with examples of the exact solutions of the equilibrium
K. C. Lee
2013-02-01
Full Text Available Multifractional Brownian motions have become popular as flexible models in describing real-life signals of high-frequency features in geoscience, microeconomics, and turbulence, to name a few. The time-changing Hurst exponent, which describes regularity levels depending on time measurements, and variance, which relates to an energy level, are two parameters that characterize multifractional Brownian motions. This research suggests a combined method of estimating the time-changing Hurst exponent and variance using the local variation of sampled paths of signals. The method consists of two phases: initially estimating global variance and then accurately estimating the time-changing Hurst exponent. A simulation study shows its performance in estimation of the parameters. The proposed method is applied to characterization of atmospheric stability in which descriptive statistics from the estimated time-changing Hurst exponent and variance classify stable atmosphere flows from unstable ones.
Lee, K. C.
2013-02-01
Multifractional Brownian motions have become popular as flexible models in describing real-life signals of high-frequency features in geoscience, microeconomics, and turbulence, to name a few. The time-changing Hurst exponent, which describes regularity levels depending on time measurements, and variance, which relates to an energy level, are two parameters that characterize multifractional Brownian motions. This research suggests a combined method of estimating the time-changing Hurst exponent and variance using the local variation of sampled paths of signals. The method consists of two phases: initially estimating global variance and then accurately estimating the time-changing Hurst exponent. A simulation study shows its performance in estimation of the parameters. The proposed method is applied to characterization of atmospheric stability in which descriptive statistics from the estimated time-changing Hurst exponent and variance classify stable atmosphere flows from unstable ones.
Magnetohydrodynamics in rectangular ducts
Magnetohydrodynamic flow in straight ducts or bends is a key issue, which has to be investigated for developing self-cooled liquid metal blankets of fusion reactors. The code presented solves the full set of governing equations and simulates all phenomena of such flows, including inertial effects. The range of application is limited by computer storage only. (orig./WL)
Topological soliton in magnetohydrodynamics
Kamchatnov, A. M.
2004-01-01
We use the Hopf mapping to construct a magnetic configuration consisting of closed field lines, each of which is linked with all the other ones. We obtain in this way a solution of the equations of magnetohydrodynamics of an ideal incompressible fluid with infinite conductivity, which describes a localized topological soliton.
The spherical torus edge region is among the most challenging regimes for plasma turbulence simulations. Here, we measure the spatial and temporal properties of ion-scale turbulence in the steep gradient region of H-mode pedestals during edge localized mode-free, MHD quiescent periods in the National Spherical Torus Experiment. Poloidal correlation lengths are about 10 ρi, and decorrelation times are about 5 a/cs. Next, we introduce a model aggregation technique to identify parametric dependencies among turbulence quantities and transport-relevant plasma parameters. The parametric dependencies show the most agreement with transport driven by trapped-electron mode, kinetic ballooning mode, and microtearing mode turbulence, and the least agreement with ion temperature gradient turbulence. In addition, the parametric dependencies are consistent with turbulence regulation by flow shear and the empirical relationship between wider pedestals and larger turbulent structures
Introduction to quantum turbulence
Barenghi, Carlo F.; Skrbek, Ladislav; Sreenivasan, Katepalli R.
2014-01-01
The term quantum turbulence denotes the turbulent motion of quantum fluids, systems such as superfluid helium and atomic Bose–Einstein condensates, which are characterized by quantized vorticity, superfluidity, and, at finite temperatures, two-fluid behavior. This article introduces their basic properties, describes types and regimes of turbulence that have been observed, and highlights similarities and differences between quantum turbulence and classical turbulence in ordinary fluids. Our ai...
D. Falceta-Gonçalves
2011-01-01
Full Text Available The Interstellar Medium (ISM is a complex, multi-phase system, where the history of the stars occurs. The processes of birth and death of stars are strongly coupled to the dynamics of the ISM. The observed chaotic and diffusive motions of the gas characterize its turbulent nature. Understanding turbulence is crucial for understanding the star-formation process and the energy-mass feedback from evolved stars. Magnetic fields, threading the ISM, are also observed, making this effort even more difficult. In this work, I briefly review the main observations and the characterization of turbulence from these observable quantities. Following on, I provide a review of the physics of magnetized turbulence. Finally, I will show the main results from theoretical and numerical simulations, which can be used to reconstruct observable quantities, and compare these predictions to the observations.
Magnetohydrodynamics and Plasma Cosmology
Kleidis, K; Papadopoulos, D B; Vlahos, L
2005-01-01
We study the linear magnetohydrodynamic (MHD) equations, both in the Newtonian and the general-relativistic limit, as regards a viscous magnetized fluid of finite conductivity and discuss instability criteria. In addition, we explore the excitation of cosmological perturbations in anisotropic spacetimes, in the presence of an ambient magnetic field. Acoustic, electromagnetic (e/m) and fast-magnetosonic modes, propagating normal to the magnetic field, can be excited, resulting in several implications of cosmological significance.
Solar Flares: Magnetohydrodynamic Processes
Kazunari Shibata; Tetsuya Magara
2011-01-01
This paper outlines the current understanding of solar flares, mainly focused on magnetohydrodynamic (MHD) processes responsible for producing a flare. Observations show that flares are one of the most explosive phenomena in the atmosphere of the Sun, releasing a huge amount of energy up to about 10^32 erg on the timescale of hours. Flares involve the heating of plasma, mass ejection, and particle acceleration that generates high-energy particles. The key physical processes for producing a fl...
Dual strings and magnetohydrodynamics
Olesen, P.
1995-01-01
We investigate whether dual strings could be solutions of the magnetohydrodynamics equations in the limit of infinite conductivity. We find that the induction equation is satisfied, and we discuss the Navier-Stokes equation (without viscosity) with the Lorentz force included. We argue that the dual string equations (with a non-universal maximum velocity) should describe the large scale motion of narrow magnetic flux tubes, because of a large reparametrization (gauge) invariance of the magneti...
Zunino, Luciano; Perez, Dario G.; Rosso, Osvaldo A.; Garavaglia, Mario
2003-01-01
The propagation of a laser beam through turbulent media is modeled as a fractional Brownian motion (fBm). Time series corresponding to the center position of the laser spot (coordinates x and y) after traveling across air in turbulent motion, with different strength, are analyzed by the wavelet theory. Two quantifiers are calculated, the Hurst exponent and the mean Normalized Total Wavelet Entropy. It is verified that both quantifiers gives complementary information about the turbulence state.
Nonresonant Grain Acceleration in MHD Turbulence
Yan, Huirong
2009-01-01
We discuss a new type of dust acceleration mechanism that acts in a turbulent magnetized medium. The magnetohydrodynamic (MHD) turbulence can accelerate grains through resonant as well as nonresonant interactions. We show that the magnetic compression provides higher velocities for super-Alfvenic turbulence and can accelerate an extended range of grains in warm media compared to gyroresonance. While fast modes dominate the acceleration for the large grains, slow modes can be important for sub...
Magnetohydrodynamic study of electromagnetic separation of nonmetallic inclusions from aluminum melt
SHU; Da(疏达); SUN; Baode(孙宝德); WANG; Jun(王俊); ZHANG; Xueping(张雪萍); ZHOU; Yaohe(周尧和)
2002-01-01
Magnetohydrodynamic flow around the nonmetallic inclusions in aluminum melt and the force exerted on the inclusions were explored by dimensional analysis and numerical calculations. Dimensional analysis shows that the invariant characterizes the force exerted on the inclusions and the flow intensity of the melt. The physical significance of A is represented as a modified particle Reynolds number that reflects the effects of electromagnetic force. The fluid flow around the particle becomes unstable when A＞2×103. It is shown that the neglect of the inertial terms has little effect on the force exerted on the inclusions in the range of A≤1×106. However, the analytical solution of the maximum velocity inside the melt does not apply due to the appearance of turbulent flow in the case of A＞2×103.
Cancellation properties in Hall magnetohydrodynamics with a strong guide magnetic field.
Martin, L N; De Vita, G; Sorriso-Valvo, L; Dmitruk, P; Nigro, G; Primavera, L; Carbone, V
2013-12-01
We present a signed measure analysis of compressible Hall magnetohydrodynamic turbulence with an external guide field. Signed measure analysis allows us to characterize the scaling behavior of the sign-oscillating flow structures and their geometrical properties (fractal dimensions of structures). A reduced numerical model, valid when a strong guide magnetic field is present, is used here. In order to discuss the effect of the Hall term, different values for the ion skin depth are considered in the simulations. Results show that as the Hall term is increased, the fractal dimension of the current and vorticity sheets decreases. This observation, together with previous analysis of the same fields, provides a comprehensive description of the effect of the Hall force on the formation of structures. Two main processes are identified, namely, the widening and unraveling of the sheets. PMID:24483577
Magnetic reversals in a simple model of magnetohydrodynamics.
Benzi, Roberto; Pinton, Jean-François
2010-07-01
We study a simple magnetohydrodynamical approach in which hydrodynamics and MHD turbulence are coupled in a shell model, with given dynamo constraints in the large scales. We consider the case of a low Prandtl number fluid for which the inertial range of the velocity field is much wider than that of the magnetic field. Random reversals of the magnetic field are observed and it shown that the magnetic field has a nontrivial evolution--linked to the nature of the hydrodynamics turbulence. PMID:20867710
Radio Wave Scintillations and Models of Interstellar Turbulence
Spangler, Steven R.
1998-05-01
There are a number of well-established observational results from radio scintillations which have implications for the nature of interstellar turbulence. Among such results are evidence for anisotropy and a Kolmogorov spectrum for the density irregularities. It is probable the galactic magnetic field organizes these irregularities so that spatial gradients along the field are much less than those perpendicular to the field. Such a behavior for turbulence is predicted by theories of magnetohydrodynamic turbulence in which the amplitude is small. The turbulence is then described by a theory termed reduced magnetohydrodynamics. A limiting case of reduced magnetohydrodynamics is two dimensional magnetohydrodynamics, in which the direction of the large scale magnetic field z defines the ignorable coordinate. Two dimensional magnetohydrodynamics consists of a pair of coupled nonlinear partial differential equations for the velocity stream function psi and the z component of the magnetic vector potential A_z. A number of observed features of interstellar turbulence can be identified with solutions to the equations of two dimensional magnetohydrodynamics. Examples are the development of Kolmogorov-like spectra for the velocity and magnetic field from a wide class (although not totally general) initial conditions, a natural explanation for the formation of intermittancy in turbulence, and the rapid development of small scale, large spatial wavenumber fluctuations, in contrast to the eddy cascade of hydrodynamic turbulence. The equations of two dimensional magnetohydrodynamics may serve as a simple but tractable model of interstellar plasma turbulence that may complement and be superior to the traditional model of an ensemble of magnetohydrodynamic waves.
Intermittent dissipation and heating in 3D kinetic plasma turbulence
Wan, M.; Matthaeus, W. H.; Roytershteyn, V.; Karimabadi, H.; Parashar, T.; Wu, P.; Shay, M. A.
2014-12-01
The nature of collisionless dissipation has been hotlydebated in recent years, with alternative ideas posed interms of various wave modes, such as kinetic Alfven waves,whistlers, linear Vlasov instabilities, cyclotron resonance,and Landau damping. Here we use large scale, fully kinetic3D simulations of collisionless plasma turbulence which showthe development of turbulence characterized by sheet-likecurrent density structures spanning a range of scales.We present evidence that these structures are sites for heatingand dissipation, and that stronger current structures signifyhigher dissipation rates. The analyses focus on quantities such as J.E, electron and proton temperatures, and PVI of the magnetic field. Evidently, kinetic scale plasma,like magnetohydrodynamics, becomes intermittent due tocurrent sheet formation, leading to the expectationthat heating and dissipation in astrophysical and space plasmasmay be highly nonuniform. Comparison with previousresults from 2D kinetic simulations, as well as high frequencysolar wind observational data will also be discussed.
Conservation of Circulation in Magnetohydrodynamics
Bekenstein, J D; Bekenstein, Jacob D.; Oron, Asaf
2000-01-01
We demonstrate, both at the Newtonian and (general) relativistic levels, theexistence of a generalization of Kelvin's circulation theorem (for pure fluids)which is applicable to perfect magnetohydrodynamics. The argument is based onthe least action principle for magnetohydrodynamic flow. Examples of the newconservation law are furnished. The new theorem should be helpful inidentifying new kinds of vortex phenomena distinct from magnetic ropes or fluidvortices.
The Modified Magnetohydrodynamical Equations
EvangelosChaliasos
2003-01-01
After finding the really self-consistent electromagnetic equations for a plasma, we proceed in a similar fashion to find how the magnetohydrodynamical equations have to be modified accordingly. Substantially this is done by replacing the "Lorentz" force equation by the correct (in our case) force equation. Formally we have to use the vector potential instead of the magnetic field intensity. The appearance of the formulae presented is the one of classical vector analysis. We thus find a set of eight equations in eight unknowns, as previously known concerning the traditional MHD equations.
The Modified Magnetohydrodynamical Equations
Evangelos Chaliasos
2003-01-01
After finding the really self-consistent electromagnetic equations for a plasma, we proceed in a similarfashion to find how the magnetohydrodynamical equations have to be modified accordingly. Substantially this is doneby replacing the "Lorentz" force equation by the correct (in our case) force equation. Formally we have to use the vectorpotential instead of the magnetic field intensity. The appearance of the formulae presented is the one of classical vectoranalysis. We thus find a set of eight equations in eight unknowns, as previously known concerning the traditional MHDequations.
Accurate, meshless methods for magnetohydrodynamics
Hopkins, Philip F.; Raives, Matthias J.
2016-01-01
Recently, we explored new meshless finite-volume Lagrangian methods for hydrodynamics: the `meshless finite mass' (MFM) and `meshless finite volume' (MFV) methods; these capture advantages of both smoothed particle hydrodynamics (SPH) and adaptive mesh refinement (AMR) schemes. We extend these to include ideal magnetohydrodynamics (MHD). The MHD equations are second-order consistent and conservative. We augment these with a divergence-cleaning scheme, which maintains nabla \\cdot B≈ 0. We implement these in the code GIZMO, together with state-of-the-art SPH MHD. We consider a large test suite, and show that on all problems the new methods are competitive with AMR using constrained transport (CT) to ensure nabla \\cdot B=0. They correctly capture the growth/structure of the magnetorotational instability, MHD turbulence, and launching of magnetic jets, in some cases converging more rapidly than state-of-the-art AMR. Compared to SPH, the MFM/MFV methods exhibit convergence at fixed neighbour number, sharp shock-capturing, and dramatically reduced noise, divergence errors, and diffusion. Still, `modern' SPH can handle most test problems, at the cost of larger kernels and `by hand' adjustment of artificial diffusion. Compared to non-moving meshes, the new methods exhibit enhanced `grid noise' but reduced advection errors and diffusion, easily include self-gravity, and feature velocity-independent errors and superior angular momentum conservation. They converge more slowly on some problems (smooth, slow-moving flows), but more rapidly on others (involving advection/rotation). In all cases, we show divergence control beyond the Powell 8-wave approach is necessary, or all methods can converge to unphysical answers even at high resolution.
Magnetohydrodynamic process in solar activity
Jingxiu Wang
2014-01-01
Full Text Available Magnetohydrodynamics is one of the major disciplines in solar physics. Vigorous magnetohydrodynamic process is taking place in the solar convection zone and atmosphere. It controls the generating and structuring of the solar magnetic fields, causes the accumulation of magnetic non-potential energy in the solar atmosphere and triggers the explosive magnetic energy release, manifested as violent solar flares and coronal mass ejections. Nowadays detailed observations in solar astrophysics from space and on the ground urge a great need for the studies of magnetohydrodynamics and plasma physics to achieve better understanding of the mechanism or mechanisms of solar activity. On the other hand, the spectacular solar activity always serves as a great laboratory of magnetohydrodynamics. In this article, we reviewed a few key unresolved problems in solar activity studies and discussed the relevant issues in solar magnetohydrodynamics.
Ultrasonic Technique in Characterization of the Grid-Generated Turbulent Flow
Andreeva, Tatiana; Durgin, William
2002-11-01
Ultrasonic time-of-flight method using dual transducers is utilized to determine some characteristics of grid-generated turbulent flow produced in a wind tunnel. The ultrasonic flowmeter equation is reconsidered, where the effects of turbulent velocity and sound speed fluctuations are included. The result is an integral equation for the corresponding correlation functions. The influence of temperature inhomogeneous on ultrasonic wave propagation is investigated using a set of experiments with a heated grid. In this paper experimentally measured travel time data are used to solve integral equation analytically in terms of correlation functions of turbulent velocity and sound speed fluctuations, and demonstrate qualitatively and quantitatively effect of turbulence on ultrasound wave propagation. First, the auto-correlation function of the travel time is constructed experimentally and is an integral of the unknown auto-correlation function of turbulent velocity. Such a relationship is known as the 2D Volterra integral equation and can be solved numerically to find the unknown auto-correlation functions of turbulent velocity and sound speed fluctuations. This leads to a new method of finding a spectral density of turbulent velocity from the flowmeter equation.
Turbulent Amplification and Structure of the Intracluster Magnetic Field
Beresnyak, Andrey; Miniati, Francesco
2015-01-01
We compare DNS calculations of homogeneous isotropic turbulence with the statistical properties of intra-cluster turbulence from the Matryoshka Run (Miniati 2014) and find remarkable similarities between their inertial ranges. This allowed us to use the time dependent statistical properties of intra-cluster turbulence to evaluate dynamo action in the intra-cluster medium, based on earlier results from numerically resolved nonlinear magneto-hydrodynamic turbulent dynamo (Beresnyak 2012). We ar...
Impact of Winds from Intermediate-Mass Stars on Molecular Cloud Structure and Turbulence
Offner, S S R
2015-01-01
Observations of nearby molecular clouds detect "shells", which are likely caused by winds from young main sequence stars. However, the progenitors of these observed features are not well characterized and the mass-loss rates inferred from the gas kinematics are several orders of magnitude greater than those predicted by atomic line-driven stellar wind models. We use magnetohydrodynamic simulations to model winds launching within turbulent molecular clouds and explore the impact of wind properties on cloud morphology and turbulence. We find that winds do not produce clear features in turbulent statistics such as the Fourier spectra of density and momentum but do impact the Fourier velocity spectrum. The density and velocity distribution functions, especially as probed by CO spectral lines, strongly indicate the presence and influence of winds. We show that stellar mass-loss rates for individual stars must be $\\dot m_w \\gtrsim 10^{-7}$ Msun yr$^{-1}$, similar to those estimated from observations, to reproduce s...
Finite dissipation and intermittency in magnetohydrodynamics
Mininni, P D
2009-01-01
We present an analysis of data stemming from numerical simulations of decaying magnetohydrodynamic (MHD) turbulence up to grid resolution of 1536^3 points and up to Taylor Reynolds number of 1200. The initial conditions are such that the initial velocity and magnetic fields are helical and in equipartition, while their correlation is negligible. Analyzing the data at the peak of dissipation, we show that the dissipation in MHD seems to asymptote to a constant as the Reynolds number increases, thereby strengthening the possibility of fast reconnection events in the solar environment for very large Reynolds numbers. Furthermore, intermittency of MHD flows, as determined by the spectrum of anomalous exponents of structure functions of the velocity and the magnetic field, is stronger than for fluids, confirming earlier results; however, we also find that there is a measurable difference between the exponents of the velocity and those of the magnetic field, as observed recently in the solar wind. Finally, we discu...
Peterson, M A
1997-01-01
The possibility is considered that turbulence is described by differential equations for which uniqueness fails maximally, at least in some limit. The inviscid Burgers equation, in the context of Onsager's suggestion that turbulence should be described by a negative absolute temperature, is such a limit. In this picture, the onset of turbulence coincides with the proliferation of singularities which characterizes the failure of uniqueness.
L. Souinida
2015-01-01
Full Text Available The aim of this work is a study of the dynamics of the diesel spray in a combustion chamber of a Diesel engine of direct injection during the injection spray. To do it, the computer code Fluent (simulator of turbulent multiphase, multi-dimensional and unsteadies flows is used to model the behavior of the spray in two dimensions. The spray evolution is simulated by using the Reynolds Averaged Navier-Stokes equations with a many models of closure of these equations, such as: Spalart Allmaras, k-ε, k-ω, k-kL-ω, SAS, RSM. The twophase flow is modeled by using the Volume of Fluid model and its coupling with Level-Set model, where the two phases behave as a pseudo-fluid with an indicator function determining the volume fraction of each phase. In this study, we demonstrate the evolution of the volume fraction field of the liquid and the averaged velocity field of the mixture, characterizing the behavior of the dynamic spray. Also, we drew the temporal evolution of the penetration length of our calculation (for all turbulence models cited above with the experiment curve and that of computer code AVBP. Finally, we deduced the appropriate model of turbulence and two-phase flow for better characterizing the dynamic of the diesel spray.
Global simulations of magnetorotational turbulence II: turbulent energetics
Parkin, E R
2013-01-01
Magnetorotational turbulence draws its energy from gravity and ultimately releases it via dissipation. However, the quantitative details of this energy flow have not been assessed for global disk models. In this work we examine the energetics of a well-resolved, three-dimensional, global magnetohydrodynamic accretion disk simulation by evaluating statistically-averaged mean-field equations for magnetic, kinetic, and internal energy using simulation data. The results reveal that turbulent magnetic (kinetic) energy is primarily injected by the correlation between Maxwell (Reynolds) stresses and shear in the (almost Keplerian) mean flow, and removed by dissipation. This finding differs from previous work using local (shearing-box) models, which indicated that turbulent kinetic energy was primarily sourced from the magnetic energy reservoir. Lorentz forces provide the bridge between the magnetic and kinetic energy reservoirs, converting ~ 1/5 of the total turbulent magnetic power input into turbulent kinetic ener...
New trends in turbulence; Turbulence: nouveaux aspects
Lesieur, M. [Institut National Polytechnique, LEGI/INPG, Institut de Mecanique, UMR 101, 38 - Grenoble (France); Yaglom, A. [Institut of Atmospheric Physics, Russian Academy of Sciences, Moscow (Russian Federation)]|[MIT, Dept. of Aeronautics and Astronautics, Cambridge, MA (United States); David, F. [CEA Saclay, SPhT, 91 - Gif-sur-Yvette (France)
2001-07-01
According to a Russian scientist, the flow of fluids actually met both in nature and engineering practice are turbulent in the overwhelmingly majority of cases. This document that reviews all the progress made recently in the understanding of turbulence, is made up of 10 courses. Course 1 ''a century of turbulence'' deals with the linear and non-linear points of views. In course 2 ''measures of anisotropy and the universal properties of turbulence'' the author gives a very complete account of fully developed turbulence experimental data both in the laboratory and in the atmosphere. Course 3 ''large-eddy simulations of turbulence (LES)'', LES are powerful tools to simulate the coherent vortices formation and evolution in a deterministic way. In Course 4 ''statistical turbulence modelling for the computation of physically complex flows'' the author describes methods used for predicting statistical industrial flows, where the geometry is right now too complex to allow the use of LES. In course 5 ''computational aero-acoustics'' an informative review of computational aero-acoustics with many applications to aircraft noise, is made. In course 6 ''the topology of turbulence'' the author presents the basis of topological fluid dynamics and stresses the importance of helicity in neutral and in magnetohydrodynamics (MHD) flows. In course 7 ''burgulence'' the authors deal with finite-time singularities, but mostly on the basis of Burger equations in one or several dimensions with the formation of multiple shocks. In course 8 ''2-dimensional turbulence'' the author presents numerous examples of 2D turbulence in the laboratory (rotating or MHD flows, plasmas), in the ocean and in the planetary atmosphere. Course 9 ''analysing and computing turbulent flows using wavelets'' is a useful presentation of
Magnetohydrodynamic inertial reference system
Eckelkamp-Baker, Dan; Sebesta, Henry R.; Burkhard, Kevin
2000-07-01
Optical platforms increasingly require attitude knowledge and optical instrument pointing at sub-microradian accuracy. No low-cost commercial system exists to provide this level of accuracy for guidance, navigation, and control. The need for small, inexpensive inertial sensors, which may be employed in pointing control systems that are required to satisfy angular line-of-sight stabilization jitter error budgets to levels of 1-3 microradian rms and less, has existed for at least two decades. Innovations and evolutions in small, low-noise inertial angular motion sensor technology and advances in the applications of the global positioning system have converged to allow improvement in acquisition, tracking and pointing solutions for a wide variety of payloads. We are developing a small, inexpensive, and high-performance inertial attitude reference system that uses our innovative magnetohydrodynamic angular rate sensor technology.
Spectrum of anomalous magnetohydrodynamics
Giovannini, Massimo
2016-05-01
The equations of anomalous magnetohydrodynamics describe an Abelian plasma where conduction and chiral currents are simultaneously present and constrained by the second law of thermodynamics. At high frequencies the magnetic currents play the leading role, and the spectrum is dominated by two-fluid effects. The system behaves instead as a single fluid in the low-frequency regime where the vortical currents induce potentially large hypermagnetic fields. After deriving the physical solutions of the generalized Appleton-Hartree equation, the corresponding dispersion relations are scrutinized and compared with the results valid for cold plasmas. Hypermagnetic knots and fluid vortices can be concurrently present at very low frequencies and suggest a qualitatively different dynamics of the hydromagnetic nonlinearities.
Essential Magnetohydrodynamics for Astrophysics
Spruit, H C
2013-01-01
This text is intended as an introduction to magnetohydrodynamics in astrophysics, emphasizing a fast path to the elements essential for physical understanding. It assumes experience with concepts from fluid mechanics: the fluid equation of motion and the Lagrangian and Eulerian descriptions of fluid flow. In addition, the basics of vector calculus and elementary special relativity are needed. Not much knowledge of electromagnetic theory is required. In fact, since MHD is much closer in spirit to fluid mechanics than to electromagnetism, an important part of the learning curve is to overcome intuitions based on the vacuum electrodynamics of one's high school days. The first chapter (only 36 pp) is meant as a practical introduction including exercises. This is the `essential' part. The exercises are important as illustrations of the points made in the text (especially the less intuitive ones). Almost all are mathematically unchallenging. The supplement in chapter 2 contains further explanations, more specialize...
Computational fusion magnetohydrodynamics
Simple magnetohydrodynamic models provide the framework for much of our understanding of the macroscopic behavior of magnetically confined laboratory plasmas. In even the simplest of models, however, the many different time and spatial scales, the multidimensionality, and the nonlinearity of the equations make finding solutions difficult. In realistic geometries obtaining quantitative results to aid our understanding, to interpret experiment, and to design new devices, involves the development of large scale numerical codes. During the past decade considerable effort has been extended in the fusion community to develop equilibrium, linear stability, and nonlinear time evolution codes in two and three dimensions, some of which have had a considerable impact on the fusion program. An overview of the various types of codes and numerical methods is given. Emphasis is on the spectrum of linear perturbations and ideal MHD stability, boundary layer methods and resistive MHD stability, and modeling of nonlinear, time evolution resistive MHD phenomena in tokamak configurations
Potential vorticity in magnetohydrodynamics
Webb, G. M.; Mace, R. L.
2015-01-01
A version of Noether's second theorem using Lagrange multipliers is used to investigate fluid relabelling symmetries conservation laws in magnetohydrodynamics (MHD). We obtain a new generalized potential vorticity type conservation equation for MHD which takes into account entropy gradients and the J × B force on the plasma due to the current J and magnetic induction B. This new conservation law for MHD is derived by using Noether's second theorem in conjunction with a class of fluid relabelling symmetries in which the symmetry generator for the Lagrange label transformations is non-parallel to the magnetic field induction in Lagrange label space. This is associated with an Abelian Lie pseudo algebra and a foliated phase space in Lagrange label space. It contains as a special case Ertel's theorem in ideal fluid mechanics. An independent derivation shows that the new conservation law is also valid for more general physical situations.
Conservation of circulation in magnetohydrodynamics
Bekenstein; Oron
2000-10-01
We demonstrate at both the Newtonian and (general) relativistic levels the existence of a generalization of Kelvin's circulation theorem (for pure fluids) that is applicable to perfect magnetohydrodynamics. The argument is based on the least action principle for magnetohydrodynamic flow. Examples of the new conservation law are furnished. The new theorem should be helpful in identifying new kinds of vortex phenomena distinct from magnetic ropes or fluid vortices. PMID:11089118
Magnetohydrodynamic process in solar activity
Jingxiu Wang; Jie Jiang
2014-01-01
Magnetohydrodynamics is one of the major disciplines in solar physics. Vigorous magnetohydrodynamic process is taking place in the solar convection zone and atmosphere. It controls the generating and structuring of the solar magnetic fields, causes the accumulation of magnetic non-potential energy in the solar atmosphere and triggers the explosive magnetic energy release, manifested as violent solar flares and coronal mass ejections. Nowadays detailed observations in solar astrophysics from s...
Ben-Nun, R.; Sheintuch, M.; Kysela, Bohuš; Konfršt, Jiří; Fořt, I.
2015-01-01
Roč. 61, č. 4 (2015), s. 1413-1426. ISSN 0001-1541 R&D Projects: GA ČR GAP101/12/2274 Institutional support: RVO:67985874 Keywords : turbulent jet * mixing * radial impeller * stirred tank * Rushton turbine Subject RIV: BK - Fluid Dynamics Impact factor: 2.748, year: 2014
Wavelet characterization of 2D turbulence and intermittency in magnetized electron plasmas
Romé, M.; Chen, S.; Maero, G.
2016-06-01
A study of the free relaxation of turbulence in a two-dimensional (2D) flow is presented, with a focus on the role of the initial vorticity conditions. Exploiting a well-known analogy with 2D inviscid incompressible fluids, the system investigated here is a magnetized pure electron plasma. The dynamics of this system are simulated by means of a 2D particle-in-cell code, starting from different spiral density (vorticity) distributions. A wavelet multiresolution analysis is adopted, which allows the coherent and incoherent parts of the flow to be separated. Comparison of the turbulent evolution in the different cases is based on the investigation of the time evolution of statistical properties, including the probability distribution functions and structure functions of the vorticity increments. It is also based on an analysis of the enstrophy evolution and its spectrum for the two components. In particular, while the statistical features assess the degree of flow intermittency, spectral analysis allows us not only to estimate the time required to reach a state of fully developed turbulence, but also estimate its dependence on the thickness of the initial spiral density distribution, accurately tracking the dynamics of both the coherent structures and the turbulent background. The results are compared with those relevant to annular initial vorticity distributions (Chen et al 2015 J. Plasma Phys. 81 495810511).
Magnetohydrodynamic Augmented Propulsion Experiment
Litchford, Ron J.; Cole, John; Lineberry, John; Chapman, Jim; Schmidt, Harold; Cook, Stephen (Technical Monitor)
2002-01-01
A fundamental obstacle to routine space access is the specific energy limitations associated with chemical fuels. In the case of vertical take-off, the high thrust needed for vertical liftoff and acceleration to orbit translates into power levels in the 10 GW range. Furthermore, useful payload mass fractions are possible only if the exhaust particle energy (i.e., exhaust velocity) is much greater than that available with traditional chemical propulsion. The electronic binding energy released by the best chemical reactions (e.g., LOX/LH2 for example, is less than 2 eV per product molecule (approx. 1.8 eV per H2O molecule), which translates into particle velocities less than 5 km/s. Useful payload fractions, however, will require exhaust velocities exceeding 15 km/s (i.e., particle energies greater than 20 eV). As an added challenge, the envisioned hypothetical RLV (reusable launch vehicle) should accomplish these amazing performance feats while providing relatively low acceleration levels to orbit (2-3g maximum). From such fundamental considerations, it is painfully obvious that planned and current RLV solutions based on chemical fuels alone represent only a temporary solution and can only result in minor gains, at best. What is truly needed is a revolutionary approach that will dramatically reduce the amount of fuel and size of the launch vehicle. This implies the need for new compact high-power energy sources as well as advanced accelerator technologies for increasing engine exhaust velocity. Electromagnetic acceleration techniques are of immense interest since they can be used to circumvent the thermal limits associated with conventional propulsion systems. This paper describes the Magnetohydrodynamic Augmented Propulsion Experiment (MAPX) being undertaken at NASA Marshall Space Flight Center (MSFC). In this experiment, a 1-MW arc heater is being used as a feeder for a 1-MW magnetohydrodynamic (MHD) accelerator. The purpose of the experiment is to demonstrate
Influences of short-wave truncations to spectral energy budget in hall MHD turbulence
The effect of a sharp short-wave truncation on Hall magnetohydrodynamic (MHD) turbulence is studied numerically to obtain basic information for constructing sub-grid-scale models of the Hall MHD equations. Hall MHD turbulence is found to be less sensitive to truncation than MHD turbulence, because the Hall term suppresses energy transfer in the magnetic field at relatively low wave numbers. (author)
Intermittency, coherent structures and dissipation in plasma turbulence
Wan, M.; Matthaeus, W. H.; Roytershteyn, V.; Parashar, T. N.; Wu, P.; Karimabadi, H.
2016-04-01
Collisionless dissipation in turbulent plasmas such as the solar wind and the solar corona has been an intensively studied subject recently, with new insights often emerging from numerical simulation. Here we report results from high resolution, fully kinetic simulations of plasma turbulence in both two (2D) and three (3D) dimensions, studying the relationship between intermittency and dissipation. The simulations show development of turbulent coherent structures, characterized by sheet-like current density structures spanning a range of scales. An approximate dissipation measure is employed, based on work done by the electromagnetic field in the local electron fluid frame. This surrogate dissipation measure is highly concentrated in small subvolumes in both 2D and 3D simulations. Fully kinetic simulations are also compared with magnetohydrodynamics (MHD) simulations in terms of coherent structures and dissipation. The interesting result emerges that the conditional averages of dissipation measure scale very similarly with normalized current density J in 2D and 3D particle-in-cell and in MHD. To the extent that the surrogate dissipation measure is accurate, this result implies that on average dissipation scales as ˜J2 in turbulent kinetic plasma. Multifractal intermittency is seen in the inertial range in both 2D and 3D, but at scales ˜ion inertial length, the scaling is closer to monofractal.
Magnetohydrodynamic power generation
The paper describes research and development in the field of magnetohydrodynamic power generation technology, based on discussions held in the Joint IAEA/UNESCO International Liaison Group on MHD electrical power generation. Research and development programmes on open cycle, closed cycle plasma and liquid-metal MHD are described. Open cycle MHD has now entered the engineering development stage. The paper reviews the results of cycle analyses and economic and environmental evaluations: substantial agreement has been reached on the expected overall performance and necessary component specifications. The achievement in the Soviet Union on the U-25 MHD pilot plant in obtaining full rated electrical power of 20.4 MW is described, as well as long duration testing of the integrated operation of MHD components. Work in the United States on coal-fired MHD generators has shown that, with slagging of the walls, a run time of about one hundred hours at the current density and electric field of a commercial MHD generator has been achieved. Progress obtained in closed cycle plasma and liquid metal MHD is reviewed. Electrical power densities of up to 140 MWe/m3 and an enthalpy extraction as high as 24 per cent have been achieved in noble gas MHD generator experiments. (Auth.)
JET snake magnetohydrodynamic equilibria
Magnetohydrodynamic (MHD) equilibrium states with a three-dimensional helical core that display the characteristics of a saturated ideal internal kink mode are computed to model snake structures that have been observed in the JET tokamak (Weller et al 1987 Phys. Rev. Lett. 59 2303). The equilibrium states are calculated with a peaked pressure profile and a weak to moderate reversed core magnetic shear with a minimum safety factor qmin near unity in the neighbourhood of the mid-radius of the plasma. Snake equilibrium states are computed in the range 0.94 min < 1.03. This range aligns with linearly unstable ideal MHD internal kink solutions of the purely axisymmetric branch of the equilibrium states. The energy difference between the bifurcated axisymmetric and helical snake equilibrium solutions is minimal. One very important novelty is that the helical structures are computed with an equilibrium code developed for three-dimensional (3D) stellarator applications in a tokamak context and cannot be obtained with standard Grad-Shafranov equation solvers. (letter)
Nonlinear ideal magnetohydrodynamics instabilities
Explosive phenomena such as internal disruptions in toroidal discharges and solar flares are difficult to explain in terms of linear instabilities. A plasma approaching a linear stability limit can, however, become nonlinearly and explosively unstable, with noninfinitesimal perturbations even before the marginal state is reached. For such investigations, a nonlinear extension of the usual MHD (magnetohydrodynamic) energy principle is helpful. (This was obtained by Merkel and Schlueter, Sitzungsberichted. Bayer. Akad. Wiss., Munich, 1976, No. 7, for Cartesian coordinate systems.) A coordinate system independent Eulerian formulation for the Lagrangian allowing for equilibria with flow and with built-in conservation laws for mass, magnetic flux, and entropy is developed in this paper which is similar to Newcomb's Lagrangian method of 1962 [Nucl. Fusion, Suppl., Pt. II, 452 (1962)]. For static equilibria nonlinear stability is completely determined by the potential energy. For a potential energy which contains second- and nth order or some more general contributions only, it is shown in full generality that linearly unstable and marginally stable systems are explosively unstable even for infinitesimal perturbations; linearly absolutely stable systems require finite initial perturbations. For equilibria with Abelian symmetries symmetry breaking initial perturbations are needed, which should be observed in numerical simulations. Nonlinear stability is proved for two simple examples, m=0 perturbations of a Bennet Z-pinch and z-independent perturbations of a θ pinch. The algebra for treating these cases reduces considerably if symmetries are taken into account from the outset, as suggested by M. N. Rosenbluth (private communication, 1992)
Generalized reduced magnetohydrodynamic equations
A new derivation of reduced magnetohydrodynamic (MHD) equations is presented. A multiple-time-scale expansion is employed. It has the advantage of clearly separating the three time scales of the problem associated with (1) MHD equilibrium, (2) fluctuations whose wave vector is aligned perpendicular to the magnetic field, and (3) those aligned parallel to the magnetic field. The derivation is carried out without relying on a large aspect ratio assumption; therefore this model can be applied to any general configuration. By accounting for the MHD equilibrium and constraints to eliminate the fast perpendicular waves, equations are derived to evolve scalar potential quantities on a time scale associated with the parallel wave vector (shear-Alfven wave time scale), which is the time scale of interest for MHD instability studies. Careful attention is given in the derivation to satisfy energy conservation and to have manifestly divergence-free magnetic fields to all orders in the expansion parameter. Additionally, neoclassical closures and equilibrium shear flow effects are easily accounted for in this model. Equations for the inner resistive layer are derived which reproduce the linear ideal and resistive stability criterion of Glasser, Greene, and Johnson. The equations have been programmed into a spectral initial value code and run with shear flow that is consistent with the equilibrium input into the code. Linear results of tearing modes with shear flow are presented which differentiate the effects of shear flow gradients in the layer with the effects of the shear flow decoupling multiple harmonics
Environmental Development Plan (EDP): magnetohydrodynamics program, FY 1977
This magnetohydrodynamics (MHD) EDP identifies and examines the environmental, health, and safety issues concerning the development of the ERDA Magnetohydrodynamics Program, the environmental activities needed to resolve these issues, applicable ongoing and completed research, and a time-phased action plan for the evaluation and mitigation of environmental impacts. A schedule for environmental research, assessment, and other activities is laid out. The purpose of the EDP is to identify environmental issues and to specify actions to ensure the environmental acceptability of commercial energy technologies being developed by ERDA. The EDP also will assist in coordinating ERDA's environmental activities with those of other government agencies. This document addresses the following technologies associated with ERDA's MHD program: (1) open-cycle magnetohydrodynamics; (2) closed-cycle plasma magnetohydrodynamics; and (3) closed-cycle liquid metal magnetohydrodynamics. The proposed environmental action plan is designed to meet the following objectives: (1) develop methods for monitoring and measuring emissions; (2) characterize air emissions, water effluents, and solid wastes from MHD; (3) determine potential environmental impacts and health hazards associated with MHD; (4) model pollutant transport and transformation; (5) ensure adequate control of pollutant emissions; (6) identify and minimize occupational health and safety hazards; (7) prepare NEPA compliance documents; and (8) assess the environmental, health, and safety impacts of the commercialized industry. This EDP will be updated and revised annually to take into account the progress of technologies toward commercialization, the environmental work accomplished, and the resolution of outstanding environmental issues concerning the technologies
Introduction to quantum turbulence.
Barenghi, Carlo F; Skrbek, Ladislav; Sreenivasan, Katepalli R
2014-03-25
The term quantum turbulence denotes the turbulent motion of quantum fluids, systems such as superfluid helium and atomic Bose-Einstein condensates, which are characterized by quantized vorticity, superfluidity, and, at finite temperatures, two-fluid behavior. This article introduces their basic properties, describes types and regimes of turbulence that have been observed, and highlights similarities and differences between quantum turbulence and classical turbulence in ordinary fluids. Our aim is also to link together the articles of this special issue and to provide a perspective of the future development of a subject that contains aspects of fluid mechanics, atomic physics, condensed matter, and low-temperature physics. PMID:24704870
Characterization of radial turbulent fluxes in the Santander linear plasma machine
Mier, J. A., E-mail: mierja@unican.es; Anabitarte, E.; Sentíes, J. M. [Departamento de Física Aplicada, Universidad de Cantabria, 39005 Santander (Spain); Sánchez, R. [Departamento de Física, Universidad Carlos III de Madrid, 28911 Leganés, Madrid (Spain); Newman, D. E. [Department of Physics, University of Alaska, Fairbanks, Alaska 99775-5920 (United States); Castellanos, O. F. [Instituto de Hidráulica Ambiental, Universidad de Cantabria, 39005 Santander (Spain); Milligen, B. Ph. van [Laboratorio Nacional de Fusión, Asociación EURATOM-CIEMAT, 28040 Madrid (Spain)
2014-05-15
It is shown that the statistical and correlation properties of the local turbulent flux measured at different radial locations of the cold, weakly ionized plasmas inside the Santander Linear Plasma Machine [Castellanos et al., Plasma Phys. Control. Fusion 47, 2067 (2005)] are consistent with diffusive-like transport dynamics. This is in contrast to the dynamical behavior inferred from similar measurements taken in hotter, fully ionized tokamak and stellarator edge plasmas, in which long-term correlations and other features characteristic of complex, non-diffusive transport dynamics have been reported in the past. These results may shed some light on a recent controversy regarding the possible universality of the dynamics of turbulent transport in magnetized plasmas.
In this work we numerically test a model of Hall magnetohydrodynamics in the presence of a strong mean magnetic field: the reduced Hall magnetohydrodynamic model (RHMHD) derived by [Gomez et al., Phys. Plasmas 15, 102303 (2008)] with the addition of weak compressible effects. The main advantage of this model lies in the reduction of computational cost. Nevertheless, up until now the degree of agreement with the original Hall MHD system and the range of validity in a regime of turbulence were not established. In this work direct numerical simulations of three-dimensional Hall MHD turbulence in the presence of a strong mean magnetic field are compared with simulations of the weak compressible RHMHD model. The results show that the degree of agreement is very high (when the different assumptions of RHMHD, such as spectral anisotropy, are satisfied). Nevertheless, when the initial conditions are isotropic but the mean magnetic field is maintained strong, the results differ at the beginning but asymptotically reach a good agreement at relatively short times. We also found evidence that the compressibility still plays a role in the dynamics of these systems, and the weak compressible RHMHD model is able to capture these effects. In conclusion the weak compressible RHMHD model is a valid approximation of the Hall MHD turbulence in the relevant physical context.
A Data-Driven Method to Characterize Turbulence-Caused Uncertainty in Wind Power Generation
Zhang, Jie; Jain, Rishabh; Hodge, Bri-Mathias
2016-10-01
A data-driven methodology is developed to analyze how ambient and wake turbulence affect the power generation of wind turbine(s). Using supervisory control and data acquisition (SCADA) data from a wind plant, we select two sets of wind velocity and power data for turbines on the edge of the plant that resemble (i) an out-of-wake scenario and (ii) an in-wake scenario. For each set of data, two surrogate models are developed to represent the turbine(s) power generation as a function of (i) the wind speed and (ii) the wind speed and turbulence intensity. Three types of uncertainties in turbine(s) power generation are investigated: (i) the uncertainty in power generation with respect to the reported power curve; (ii) the uncertainty in power generation with respect to the estimated power response that accounts for only mean wind speed; and (iii) the uncertainty in power generation with respect to the estimated power response that accounts for both mean wind speed and turbulence intensity. Results show that (i) the turbine(s) generally produce more power under the in-wake scenario than under the out-of-wake scenario with the same wind speed; and (ii) there is relatively more uncertainty in the power generation under the in-wake scenario than under the out-of-wake scenario.
Broken Ergodicity in MHD Turbulence
Shebalin, John V.
2010-01-01
Ideal magnetohydrodynamic (MHD) turbulence may be represented by finite Fourier series, where the inherent periodic box serves as a surrogate for a bounded astrophysical plasma. Independent Fourier coefficients form a canonical ensemble described by a Gaussian probability density function containing a Hermitian covariance matrix with positive eigenvalues. The eigenvalues at lowest wave number can be very small, resulting in a large-scale coherent structure: a turbulent dynamo. This is seen in computations and a theoretical explanation in terms of 'broken ergodicity' contains Taylor s theory of force-free states. An important problem for future work is the case of real, i.e., dissipative flows. In real flows, broken ergodicity and coherent structure are still expected to occur in MHD turbulence at the largest scale, as suggested by low resolution simulations. One challenge is to incorporate coherent structure at the largest scale into the theory of turbulent fluctuations at smaller scales.
Statistical properties of turbulence: An overview
Rahul Pandit; Prasad Perlekar; Samriddhi Sankar Ray
2009-07-01
We present an introductory overview of several challenging problems in the statistical characterization of turbulence. We provide examples from fluid turbulence in three and two dimensions, from the turbulent advection of passive scalars, turbulence in the one-dimensional Burgers equation, and fluid turbulence in the presence of polymer additives.
Solar Flares: Magnetohydrodynamic Processes
Kazunari Shibata
2011-12-01
Full Text Available This paper outlines the current understanding of solar flares, mainly focused on magnetohydrodynamic (MHD processes responsible for producing a flare. Observations show that flares are one of the most explosive phenomena in the atmosphere of the Sun, releasing a huge amount of energy up to about 10^32 erg on the timescale of hours. Flares involve the heating of plasma, mass ejection, and particle acceleration that generates high-energy particles. The key physical processes for producing a flare are: the emergence of magnetic field from the solar interior to the solar atmosphere (flux emergence, local enhancement of electric current in the corona (formation of a current sheet, and rapid dissipation of electric current (magnetic reconnection that causes shock heating, mass ejection, and particle acceleration. The evolution toward the onset of a flare is rather quasi-static when free energy is accumulated in the form of coronal electric current (field-aligned current, more precisely, while the dissipation of coronal current proceeds rapidly, producing various dynamic events that affect lower atmospheres such as the chromosphere and photosphere. Flares manifest such rapid dissipation of coronal current, and their theoretical modeling has been developed in accordance with observations, in which numerical simulations proved to be a strong tool reproducing the time-dependent, nonlinear evolution of a flare. We review the models proposed to explain the physical mechanism of flares, giving an comprehensive explanation of the key processes mentioned above. We start with basic properties of flares, then go into the details of energy build-up, release and transport in flares where magnetic reconnection works as the central engine to produce a flare.
Effect of Externally Driven Magnetic Islands on Resistive Ballooning Turbulence
NISHIMURA, Seiya; Yagi, Masatoshi
2011-01-01
Turbulent transport in the edge region of tokamak plasmas is simulated using a reduced set of magnetohydrodynamic equations. Repetitive and intermittent transport bursts driven by resistive ballooning turbulence with external heating are observed. The effect of a resonant magnetic perturbation (RMP) on turbulent heat transport is examined, where the electromagnetic response of the plasma to the RMP is solved consistently. The penetration of the RMP excites a magnetic island chain and damps th...
Extended Kelvin theorem in relativistic magnetohydrodynamics
Bekenstein, Jacob D.; Oron, Asaf
2000-01-01
We prove the existence of a generalization of Kelvin's circulation theorem in general relativity which is applicable to perfect isentropic magnetohydrodynamic flow. The argument is based on a new version of the Lagrangian for perfect magnetohydrodynamics. We illustrate the new conserved circulation with the example of a relativistic magnetohydrodynamic flow possessing three symmetries.