Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical Study of Compressible Magnetohydrodynamic Turbulence in Two Dimensions
50
0
0.0
(
0
)
Authors
Hyesook Lee
Ask ChatGPT about the research
No Arabic abstract
We have studied forced turbulence of compressible magnetohydrodynamic (MHD) flows through two-dimensional simulations with different numerical resolutions. First, hydrodynamic turbulence with Mach number $<M_s >_{rm init} equiv < v >_{rm rms}/ c_s = 1$ and density compression ${< deltarho / rho >}_{rm rms} simeq 0.45$ was generated by enforcing a random force. Then, initial, uniform magnetic fields of various strengths were added with Alfvenic Mach number $<M_A >_{rm init} equiv < v >_{rm rms} / c_{A, {rm init}} gg 1$. An isothermal equation of state was employed, and no explicit dissipation was included. After the MHD turbulence is saturated, the resulting flows are categorized as very weak field (VWF), weak field (WF), and strong field (SF) classes, which have $<M_A > equiv < v >_{rm rms} / < c_A >_{rm rms} gg 1$, $<M_A > > 1$, and $<M_A > sim 1$, respectively. Not only in the SF regime but also in the WF regime, turbulent transport is suppressed by the magnetic field. In the SF cases, the energy power spectra in the inertial range, although no longer power-law, exhibit a range with slopes close to $sim1.5$, hinting the Iroshnikov-Kraichnan spectrum. Our simulations were able to produce the SF class behaviors only with high resolution of at least $1024^2$ grid cells. The specific requirements for the simulation of the SF class should depend on the code (and the numerical scheme) as well as the initial setup, but our results do indicate that very high resolution would be required for converged results in simulation studies of MHD turbulence.
rate research
Read More
This paper presents an extension of the hybrid scheme proposed by Wang et al. (J. Comput. Phys. 229 (2010) 169-180) for numerical simulation of compressible isotropic turbulence to flows with higher turbulent Mach numbers. The scheme still utilizes an 8th-order compact scheme with built-in hyperviscosity for smooth regions and a 7th-order WENO scheme for highly compressive regions, but now both in their conservation formulations and for the latter with the Roe type characteristic-wise reconstruction. To enhance the robustness of the WENO scheme without compromising its high-resolution and accuracy, the recursive-order-reduction procedure is adopted, where a new type of reconstruction-failure-detection criterion is constructed. To capture the upwind direction properly in extreme conditions, the global Lax-Friedrichs numerical flux is used. In addition, a new form of cooling function is proposed, which is proved to be positivity-preserving. With these techniques, the new scheme not only inherits the good properties of the original one but also extends largely the computable range of turbulent Mach number, which has been further confirmed by numerical results.
In this work a weak-turbulence closure is used to determine the structure of the two-time power spectrum of weak magnetohydrodynamic (MHD) turbulence from the nonlinear equations describing the dynamics. The two-time energy spectrum is a fundamental quantity in turbulence theory from which most statistical properties of a homogeneous turbulent system can be derived. A closely related quantity, obtained via a spatial Fourier transform, is the two-point two-time correlation function describing the space-time correlations arising from the underlying dynamics of the turbulent fluctuations. Both quantities are central in fundamental turbulence theories as well as in the analysis of turbulence experiments and simulations. However, a first-principles derivation of these quantities has remained elusive due to the statistical closure problem, in which dynamical equations for correlations at order $n$ depend on correlations of order $n+1$. The recent launch of the Parker Solar Probe (PSP), which will explore the near-Sun region where the solar wind is born, has renewed the interest in the scientific community to understand the structure, and possible universal properties of space-time correlations. The weak MHD turbulence regime that we consider in this work allows for a natural asymptotic closure of the two-time spectrum, which may be applicable to other weak turbulence regimes found in fluids and plasmas. An integro-differential equation for the scale-dependent temporal correlation function is derived for incompressible Alfvenic fluctuations whose nonlinear dynamics is described by the reduced MHD equations.
White-light observations of the solar corona show that there are two characteristic types of Coronal Mass Ejections (CMEs) in terms of speed-height profiles: so-called fast CMEs that attain high speeds low in the corona and slow CMEs that gradually accelerate from low initial speeds. Low and Zhang (2002) have recently proposed that fast and slow CMEs result from initial states with magnetic configurations characterized by normal prominences (NPs) and inverse prominences (IPs), respectively. To test their theory, we employed a two-dimensional, time-dependent, resistive magnetohydrodynamic code to simulate the expulsion of CMEs in these two different prominence environments. Our numerical simulations demonstrate that (i) a CME-like expulsion is more readily produced in an NP than in an IP environment, and, (ii) a CME originating from an NP environment tends to have a higher speed early in the event than one originating from an IP environment. Magnetic reconnection plays distinct roles in the two different field topologies of these two environments to produce their characteristic CME speed-height profiles. Our numerical simulations support the proposal of Low and Zhang (2002) although the reconnection development for the NP associated CME is different from the one sketched in their theory. Observational implications of our simulations are discussed.
Under suitable forcing a fluid exhibits turbulence, with characteristics strongly affected by the fluids confining geometry. Here we study two-dimensional quantum turbulence in a highly oblate Bose-Einstein condensate in an annular trap. As a compressible quantum fluid, this system affords a rich phenomenology, allowing coupling between vortex and acoustic energy. Small-scale stirring generates an experimentally observed disordered vortex distribution that evolves into large-scale flow in the form of a persistent current. Numerical simulation of the experiment reveals additional characteristics of two-dimensional quantum turbulence: spontaneous clustering of same-circulation vortices, and an incompressible energy spectrum with $k^{-5/3}$ dependence for low wavenumbers $k$ and $k^{-3}$ dependence for high $k$.
We present a search for conformal invariance in vorticity isolines of two-dimensional compressible turbulence. The vorticity is measured by tracking the motion of particles that float at the surface of a turbulent tank of water. The three-dimensional turbulence in the tank has a Taylor microscale $Re_lambda simeq 160$. The conformal invariance theory being tested here is related to the behavior of equilibrium systems near a critical point. This theory is associated with the work of Lowner, Schramm and others and is usually referred to as Schramm-Lowner Evolution (SLE). The system was exposed to several tests of SLE. The results of these tests suggest that zero-vorticity isolines exhibit noticeable departures from this type of conformal invariance.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
