Direct numerical simulation is applied to investigate instability and transition to turbulence in the flow of an electrically conducting incompressible fluid between two parallel unbounded insulating walls affected by a wall-normal magnetic field (the Hartmann flow). The linear stability analysis of this flow provided unrealistically high critical Reynolds numbers, about two orders of magnitude higher than those observed in experiments. We propose an explanation based on the streak growth and breakdown mechanism described earlier for other shear flows. The mechanism is investigated using a two-step procedure that includes transient growth of two-dimensional optimal perturbations and the subsequent three-dimensional instability of the modulated streaky flow. In agreement with recent experimental investigations the calculations produce a critical range between 350 and 400 for the Hartmann thickness based Reynolds number, where the transition occurs at realistic amplitudes of two-and three-dimensional perturbations.

Many practical applications exploit an external local magnetic field -magnetic obstacle -as an essential part of their construction. Recently, have demonstrated that the flow of an electrically conducting fluid influenced by an external field can show several kinds of recirculation. The present paper reports a 3D numerical study whose some results are compared with an experiment about such a flow in a rectangular duct. First, we derive equations to compute analytically an external magnetic field and verify these equations by comparing with experimentally measured field intensity. Then, we study flow characteristics for different magnetic field configurations. The flow inside the magnetic gap is dependent mainly on the interaction parameter N , which represents the ratio of the Lorentz force to the inertial force. Depending on the constrainment factor κ = M y /L y , where M y and L y are half-widths of the external magnet and duct, the flow can show different stationary recirculation patterns: two magnetic vortices at small κ, a six-vortex ensemble at moderate κ, and no vortices at large κ. Recirculation appears when N is higher than a critical value N c,m . The driving force for the recirculation is the reverse electromotive force that arises to balance the reverse electrostatic field. The reversion of the electrostatic field is caused by a concurrence of internal and external vorticity correspondingly related to the internal and external slopes in the M -shaped velocity profile. The critical value of N c,m quickly grows as κ increases. For the case of well developed recirculation, the numerical reverse velocity agrees well with that obtained in physical experiments. Two different magnetic systems induce the same electric field and stagnancy region provided these systems have the same power of recirculation given by the N/N c,m ratio. The 3D helical peculiarities of the vortices are elaborated, and an analogy is shown to exist between a helical motion inside the studied recirculation and a secondary motion in the process of the Ekman pumping. Finally, it is shown that a 2D model fails to properly produce stable two and six-vortex structures as found in the 3D system. Interestingly, these recirculation patterns appear only as time dependent and unstable transitional states before the Karman vortex street forms, when one suddenly applies a retarding local magnetic field on a constant flow.

Mean flow properties of turbulent magnetohydrodynamic channel flow with electrically insulating channel walls are studied using high-resolution direct numerical simulations. The Lorentz force due to the homogeneous wall-normal magnetic field is computed in the quasi-static approximation. For strong magnetic fields, the mean velocity profile shows a clear three-layer structure consisting of a viscous region near each wall and a plateau in the middle connected by logarithmic layers. This structure reflects the significance of viscous, turbulent, and electromagnetic stresses in the streamwise momentum balance dominating the viscous, logarithmic, and plateau regions, respectively. The width of the logarithmic layers changes with the ratio of Reynolds- and Hartmann numbers. Turbulent stresses typically decay more rapidly away from the walls than predicted by mixing-length models.

Based on molecular dynamics simulations of a lithium metasilicate glass we study the potential of bond valence sum calculations to identify sites and diffusion pathways of mobile Li ions in a glassy silicate network. We find that the bond valence method is not well suitable to locate the sites, but allows one to estimate the number of sites. Spatial regions of the glass determined as accessible for the Li ions by the bond valence method can capture up to 90% of the diffusion path. These regions however entail a significant fraction that does not belong to the diffusion path. Because of this low specificity, care must be taken to determine the diffusive motion of particles in amorphous systems based on the bond valence method. The best identification of the diffusion path is achieved by using a modified valence mismatch in the BV analysis that takes into account that a Li ion favors equal partial valences to the neighboring oxygen ions. Using this modified valence mismatch it is possible to replace hard geometric constraints formerly applied in the BV method. Further investigations are necessary to better understand the relation between the complex structure of the host network and the ionic diffusion paths.

We use a combination of numerical simulations and experiments to elucidate the structure of the flow of an electrically conducting fluid past a localized magnetic field, called magnetic obstacle. We demonstrate that the stationary flow pattern is considerably more complex than in the wake behind an ordinary body. The steady flow is shown to undergo two bifurcations (rather than one) and to involve up to six (rather than just two) vortices. We find that the first bifurcation leads to the formation of a pair of vortices within the region of magnetic field that we call inner magnetic vortices, whereas a second bifurcation gives rise to a pair of attached vortices that are linked to the inner vortices by connecting vortices.

No abstract

No abstract

An investigation of the decay laws of energy and of higher moments of the Elsässer fields z±=v±b in the self-similar regime of magnetohydrodynamic (MHD) turbulence is presented, using phenomenological models as well as two-dimensional numerical simulations with periodic boundary conditions and up to 20482 grid points. The results are compared with the generalization of the parameter-free model derived by Galtier et al. [Phys. Rev. Lett.79, 2807 (1997)], which takes into account the slowing down of the dynamics due to the propagation of Alfvén waves. The new model developed here allows for a study in terms of one parameter governing the wavenumber dependence of the energy spectrum at scales of the order of (and larger than) the integral scale of the flow. The one-dimensional compressible case is also dealt with in two of its simplest configurations. Computations are performed for a standard Laplacian diffusion as well as with a hyperdiffusive algorithm. The results are sensitive to the amount of correlation between the velocity and the magnetic field, but rather insensitive to all other parameters such as the initial ratio of kinetic to magnetic energy or the presence or absence of a uniform component of the magnetic field. In all cases, the decay is significantly slower than for neutral fluids in a way that favours for MHD flows the phenomenology of Iroshnikov [Soviet Astron.7, 566 (1963)] and Kraichnan [Phys. Fluids8, 1385 (1965)] as opposed to that of Kolmogorov [Dokl. Akad. Nauk. SSSR31, 538 (1941)]. The temporal evolution of q-moments of the generalized vorticities 〈[mid ]ω±[mid ]q〉 =〈[mid ]ω±j[mid ]q〉 up to order q=10 is also given, and is compared with the prediction of the model. Less agreement obtains as q grows – a fact probably due to intermittency and the development of coherent structures in the form of eddies, and of vorticity and current sheets.