Thomas Boeck

Numerical study of turbulent magnetohydrodynamic channel 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.

Transient Rayleigh-Benard-Marangoni Convection due to Evaporation : a Linear Non-normal Stability Analysis

The convective instability in a plane liquid layer with time-dependent temperature profile is investigated by means of a general method suitable for linear stability analysis of an unsteady basic flow. The method is based on a non-normal approach, and predicts the onset of instability, critical wavenumber and time. The method is applied to transient Rayleigh–Benard–Marangoni convection due to cooling by evaporation. Numerical results as well as theoretical scalings for the critical parameters as function of the Biot number are presented for the limiting cases of purely buoyancy-driven and purely surface-tension-driven convection. Critical parameters from calculations are in good agreement with those from experiments on drying polymer solutions, where the surface cooling is induced by solvent evaporation.

Magnetohydrodynamic turbulence in a channel with spanwise magnetic field

The effect of a uniform spanwise magnetic field on a turbulent channel flow is investigated for the case of a low magnetic Reynolds number. Direct numerical simulation (DNS) and large eddy simulation (LES) computations are performed for two values of the hydrodynamic Reynolds number (104 and 2×104) and with the Hartmann number varying in a wide range. It is shown that the main effect of the magnetic field is the suppression of turbulent velocity fluctuations and momentum transfer in the wall-normal direction. This leads to drag reduction and transformation of the mean flow profile. The centerline velocity grows, the mean velocity gradients near the wall decrease, and the typical horizontal dimensions of the coherent structures enlarge upon increasing the Hartmann number. Comparison between LES and DNS results shows that the dynamic Smagorinsky model accurately reproduces the flow transformation.

Optimal linear growth in magnetohydrodynamic duct flow

Transient linear growth in laminar magnetohydrodynamic duct flow is analysed. The duct is straight with rectangular cross-section and electrically insulating walls. The applied uniform magnetic field is oriented perpendicular to the mean flow direction and parallel to one of the walls. Optimal perturbations and their maximum amplifications over finite time intervals are computed. The optimal perturbations are increasingly damped by the magnetic field, localized in the boundary layers parallel to the magnetic field irrespective of the duct aspect ratio. Typically, the optimal perturbations have non-vanishing streamwise wavenumber as found in magnetohydrodynamic channel flow with spanwise magnetic field. The Hartmann boundary layers perpendicular to the magnetic field do not contribute to the transient growth.

Optimal growth and transition to turbulence in channel flow with spanwise magnetic field

Instability and transition to turbulence in a magnetohydrodynamic channel flow are studied numerically for the case of a uniform magnetic field imposed along the spanwise direction. Optimal perturbations and their maximum amplifications over finite time intervals are computed in the framework of the linear problem using an iterative scheme based on direct and adjoint governing equations. It is shown that, at sufficiently strong magnetic field, the maximum amplification is no longer provided by classical streamwise rolls, but rather by rolls oriented at an oblique angle to the basic flow direction. The angle grows with the Hartmann number Ha and reaches the limit corresponding to purely spanwise rolls at Ha between 50 and 100 depending on the Reynolds number. Direct numerical simulations are applied to investigate the transition to turbulence at a single subcritical Reynolds number Re = 5000 and various Hartmann numbers. The transition is caused by the transient growth and subsequent breakdown of optimal perturbations, which take the form of one or two symmetric optimal modes (streamwise, oblique or spanwise modes depending on Ha) with low-amplitude three-dimensional noise added at the moment of strongest energy amplification. A sufficiently strong magnetic field (Ha larger than approximately 30) is found to completely suppress the instability. At smaller Hartmann numbers, the transition is observed but it is modified in comparison with the pure hydrodynamic case.

Pattern formation and mass transfer under stationary solutal Marangoni instability

According to the seminal theory by Sternling and Scriven, solutal Marangoni convection during mass transfer of surface-active solutes may occur as either oscillatory or stationary instability. With strong support of Manuel G. Velarde, a combined initiative of experimental works, in particular to mention those of Linde, Wierschem and coworkers, and theory has enabled a classification of dominant wave types of the oscillatory mode and their interactions. In this way a rather comprehensive understanding of the nonlinear evolution of the oscillatory instability could be achieved. A comparably advanced state-of-the-art with respect to the stationary counterpart seemed to be out of reach a short time ago. Recent developments on both the numerical and experimental side, in combination with assessing an extensive number of older experiments, now allow one to draw a more unified picture. By reviewing these works, we show that three main building blocks exist during the nonlinear evolution: roll cells, relaxation oscillations and relaxation oscillations waves. What is frequently called interfacial turbulence results from the interaction between these partly coexisting basic patterns which may additionally occur in different hierarchy levels. The second focus of this review lies on the practical importance of such convection patterns concerning their influence on mass transfer characteristics. Particular attention is paid here to the interaction between Marangoni and buoyancy effects which frequently complicates the pattern formation even more. To shed more light on these dependencies, new simulations regarding the limiting case of stabilizing density stratification and vanishing buoyancy are incorporated.

Large-scale intermittency of liquid-metal channel flow in a magnetic field

We predict a novel flow regime in liquid metals under the influence of a magnetic field. It is characterized by long periods of nearly steady, two-dimensional flow interrupted by violent three-dimensional bursts. Our prediction has been obtained from direct numerical simulations in a channel geometry at low magnetic Reynolds number and translates into physical parameters which are amenable to experimental verification under laboratory conditions. The new regime occurs in a wide range of parameters and may have implications for metallurgical applications.

Multiscale structures in solutal Marangoni convection: Three-dimensional simulations and supporting experiments

Transient solutal Marangoni convection in a closed two-layer system is studied by a combination of numerical simulations and supplementary validation experiments. The initially quiescent, equally sized liquid layers are the phases of a cyclohexanol/water mixture. Butanol is additionally dissolved in the upper organic layer. Its diffusion across the interface is sensitive to the Marangoni instability. Complex convective patterns emerge that develop a hierarchical cellular structure in the course of the mass transfer. Our highly resolved simulations based on a pseudospectral method are the first to successfully reproduce the multiscale flow observed in the experiments. We solve the three-dimensional Navier-Stokes-Boussinesq equations with an undeformable interface, which is modeled using the linear Henry relation for the partition of the weakly surface-active butanol. Length scales in the concentration and velocity fields associated with the small and large-scale cells agree well with our experimental data from shadowgraph images. Moreover, the simulations provide detailed information on the local properties of the flow by which the evolution of the patterns and their vertical structure are analyzed. Apart from relatively weak influences due to buoyancy, the evolution of the convective structures is self-similar between different initial butanol concentrations when length and time are appropriately rescaled.

Viscous instability of a sheared liquid-gas interface: Dependence on fluid properties and basic velocity profile

In the framework of linear stability theory, we analyze how a liquid-gas mixing layer is affected by several parameters: viscosity ratio, density ratio, and several length scales. These scales reflect the presence of a velocity defect induced by the wake behind the splitter plate and the presence of boundary layers which develop ahead of the plate trailing edge. Incorporating such effects, we compute the various temporal and spatial instability modes and identify their driving instability mechanism based on their Reynolds number dependence, spatial structure, and energy budget. It is examined how the velocity defect modifies the temporal and the spatial stability properties. In addition, the transition from convective to absolute instability occurs at lower velocity contrast between gas and liquid free streams when a defect is present. This transition is also promoted by surface tension. Compared to inviscid stability computations, our spatial stability analysis displays a better agreement with measured g...

Electromagnetic drag on a magnetic dipole near a translating conducting bar

The electromagnetic drag force and torque acting on a magnetic dipole due to the translatory motion of an electrically conducting bar with square cross section and infinite length is computed by numerical analysis for different orientations and locations of the dipole. The study is motivated by the novel techniques termed Lorentz force velocimetry and Lorentz force eddy current testing for noncontact measurements of the velocity of a conducting liquid and for detection of defects in the interior of solid bodies, respectively. The present, simplified configuration provides and explains important scaling laws and reference results that can be used for verification of future complete numerical simulations of more realistic problems and complex geometries. The results of computations are also compared with existing analytical solutions for an infinite plate and with a newly developed asymptotic theory for large distances between the bar and the magnetic dipole. We finally discuss the optimization problem of f...