No Arabic abstract
We propose two sets of initial conditions for magnetohydrodynamics (MHD) in which both the velocity and the magnetic fields have spatial symmetries that are preserved by the dynamical equations as the system evolves. When implemented numerically they allow for substantial savings in CPU time and memory storage requirements for a given resolved scale separation. Basic properties of these Taylor-Green flows generalized to MHD are given, and the ideal non-dissipative case is studied up to the equivalent of 2048^3 grid points for one of these flows. The temporal evolution of the logarithmic decrements, delta, of the energy spectrum remains exponential at the highest spatial resolution considered, for which an acceleration is observed briefly before the grid resolution is reached. Up to the end of the exponential decay of delta, the behavior is consistent with a regular flow with no appearance of a singularity. The subsequent short acceleration in the formation of small magnetic scales can be associated with a near collision of two current sheets driven together by magnetic pressure. It leads to strong gradients with a fast rotation of the direction of the magnetic field, a feature also observed in the solar wind.
Evolution of stochastically homogeneous magnetic field advected by incompressible turbulent flow with large magnetic Prandtl numbers is considered at the scales less than Kolmogorov viscous scale. It is shown that, despite unlimited growth of the magnetic field, its feedback on the fluids dynamics remains negligibly small.
In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We include these non-ideal terms for two MHD techniques: the Powell 8-wave formalism and a constrained transport scheme, which evolves the cell-centred magnetic vector potential. We test our implementation against problems of increasing complexity, such as one- and two-dimensional diffusion problems, and the evolution of progressive and stationary Alfven waves. On these test problems, our implementation recovers the analytic solutions to second-order accuracy. As first applications, we investigate the tearing instability in magnetized plasmas and the gravitational collapse of a rotating magnetized gas cloud. In both systems, resistivity plays a key role. In the former case, it allows for the development of the tearing instability through reconnection of the magnetic field lines. In the latter, the adopted (constant) value of ohmic resistivity has an impact on both the gas distribution around the emerging protostar and the mass loading of magnetically driven outflows. Our new non-ideal MHD implementation opens up the possibility to study magneto-hydrodynamical systems on a moving mesh beyond the ideal MHD approximation.
Dielectric barrier discharge (DBD) plasma actuators are an attractive option for separation control, lift enhancement, and drag reduction. Some plasma actuators feature optimized electrode shapes, electrical waveforms to maximize the aerodynamic forces at higher angles of attack. Here, we analyze the performance of a direct current augmented DBD (DBD - DCA) actuator with a sawtooth shape exposed electrode. The active electrode was positioned at 18% chord and the electrode at 48% chord of NACA 0012 airfoil. Wind tunnel experiments were conducted at wind speeds of 15 - 25 m/s, corresponding to Reynolds number Re = 201k - 335k. Lift coefficient (C$_L$), drag coefficient (C$_D$), and pitching moment coefficients (C$_M$), were measured with and without plasma actuation for angles of attack $alpha$ = 0$^o$ - 8$^o$ and the DCA electrode potential ($varphi_{DC}$) was varied from 0 kV to -15 kV. With energized DCA electrode, the C$_L$ increases up to 0.03 and the C$_D$ decreases by 50% at 15 m/s flow speeds and 0$^o$ angle of attack, the results are similar throughout the range of $alpha$. The effect of the actuator at higher Re diminishes, suggesting that the maximum control authority could be achieved at lower wind speeds.
We present a model describing evolution of the small-scale Navier-Stokes turbulence due to its stochastic distortions by much larger turbulent scales. This study is motivated by numerical findings (laval, 2001) that such interactions of separated scales play important role in turbulence intermittency. We introduce description of turbulence in terms of the moments of the k-space quantities using a method previously developed for the kinematic dynamo problem (Nazarenko, 2003). Working with the $k$-space moments allows to introduce new useful measures of intermittency such as the mean polarization and the spectral flatness. Our study of the 2D turbulence shows that the energy cascade is scale invariant and Gaussian whereas the enstrophy cascade is intermittent. In 3D, we show that the statistics of turbulence wavepackets deviates from gaussianity toward dominance of the plane polarizations. Such turbulence is formed by ellipsoids in the $k$-space centered at its origin and having one large, one neutral and one small axes with the velocity field pointing parallel to the smallest axis.
The self-similar Richardson cascade admits two logically possible scenarios of small-scale turbulence at high Reynolds numbers. In the first scenario, eddies population densities vary as a function of eddies scales. As a result, one or a few eddy types dominate at small scales, and small-scale turbulence lacks diversity. In the second scenario, eddies population densities are scale-invariant across the inertial range, resulting in small-scale diversity. That is, there are as many types of eddies at the small scales as at the large scales. In this letter, we measure eddies population densities in three-dimensional isotropic turbulence and determine the nature of small-scale turbulence. The result shows that eddies population densities are scale-invariant.