A finite-volume shallow layer method for the MHD instabilities in an aluminium production cell

An electrolytic cell for Aluminium production contains molten metal and molten electrolyte, which are subject to high dc-currents and magnetic fields. Lorentz forces arising from the cross product of current and magnetic field may amplify natural gravity waves at the interface between the two fluids, leading to short circuits in extreme cases. The external magnetic field and current distribution in the production cell is computed through a detailed finite element analysis at Torino Polytechnic. The results are then used to compute the magnetohydrodynamic and thermal effects in the aluminium/electrolyte bath. Each cell has lateral dimensions of 6m x 2m, whilst the bath depth is only 30cm. the electrically resistive electrolyte path, which is critical in the operation of the cell, has layer depth of only a few centimetres below each carbon anode. Because the shallow dimensions of the liquid layer a finite-volume shallow-layer technique has been used at Greenwich to compute the resulting flow-field and interface perturbations. The information obtained from this method, i.e. depth averaged velocities and aluminium/electrolyte interface position is then embedded in the three-dimensional finite volume code PHYSICA and will be used to compute the heat transfer and phase change in the cell.

Les instabilités magnétohydrodynamiques dans EULAG-MHD.

Les mécanismes qui entretiennent le cycle magnétique solaire sont encore aujourd’hui relativement mal compris. Entre autres, certains travaux suggèrent la présence d’insta- bilités magnétohydrodynamiques qui pourraient avoir une influence significative sur la période du cycle par leur capacité à accélérer la destruction des structures magnétiques à grandes échelles. Nous analysons la présence des instabilités au sein des simulations effectuées à l’aide du modèle EULAG-MHD en utilisant premièrement une méthodologie basée sur un proxy spécifique associé à l’instabilité et en le comparant à un proxy similaire, mais pour le cycle magnétique solaire observé dans notre modèle. Cette méthodologie fait ressortir une évolution temporellement cyclique du proxy de l’instabilité avec des amplitudes similaires au proxy du cycle magnétique, mais présentant un léger déphasage. Nous poursuivons cette analyse en appliquant une méthode se basant sur le découpage de “zones instables” via le critère de Tayler dans la zone stable de notre modèle. L’application expose une migration équatoriale d’une zone instable qui débute à très hautes latitudes aux pôles, passe par le champ toroïdal classique, accélère et atteint l’équateur. Cette zone instable semble accélérer la destruction du champ magnétique présent, laissant place au nouveau champ pour la prochaine demie-période du cycle. La combinaison de ces deux analyses permet d’énoncer un scénario plausible qui inclut les effets d’une telle instabilité sur le cycle magnétique ainsi que sur la stabilité globale de notre simulation. Dans ce scénario, il est important de noter que les inversions de polarités semblent indépendantes de cette instabilité, qui ne ferait qu’accélérer le processus de destruction du champ déjà en place.

Instabilidades MHD no Tokamak TCABR

Este trabalho descreve o estudo das instabilidades magneto-hidro-dinâmicas (MHD) comumente observadas nas descargas elétricas de plasma no tokamak TCABR, do Instituto de Física da USP. Dois diagnósticos principais foram empregados para observar essas instabilidades: um conjunto poloidal de 24 bobinas magnéticas (bobinas de Mirnov) colocadas próximas à borda do plasma e um medidor de emissões na faixa do Ultra Violeta e de raios X moles com 20 canais (sistema SXR), cujo circuito de condicionamento de sinais foi aprimorado como parte deste trabalho. Esses diagnósticos foram escolhidos porque fornecem informações complementares, uma vez que o sistema SXR observa a parte central da coluna de plasma, enquanto as bobinas de Mirnov detectam as instabilidades MHD na região mais externa da coluna. As informações coletadas por esses diagnósticos foram submetidas à análise espectral com resolução temporal e espacial, possibilitando determinar a evolução das características espectrais e espaciais das instabilidades MHD observadas. Essas análises revelaram que durante a etapa inicial da formação do plasma (quando a corrente de plasma ainda está aumentando) ilhas magnéticas com números de onda decrescente, identificadas como sendo modos kink de borda, são detectadas nas bobinas de Mirnov. Após a formação do plasma, quando os parâmetros de equilíbrio estão relativamente estáveis (platô), oscilações são detectadas tanto nas bobinas de Mirnov quanto no sistema de SXR, indicando a presença de instabilidades MHD em toda a coluna de plasma. Em geral as oscilações medidas nas bobinas de Mirnov tem baixa amplitude e correspondem a pequenas ilhas magnéticas que foram identificadas como sendo modos de ruptura (modos tearing). Por outro lado, as instabilidades na região central foram identificadas como dentes de serra, que correspondem a relaxações periódicas da região interna à superfície magnética com fator de segurança q=1 e que são acompanhadas de oscilações precursoras, cuja amplitude depende da fase do ciclo de relaxação. Devido à essa modulação de amplitude, aparecem picos de frequência satélite nos espectrogramas dos sinais do SXR. Além disso, devido ao fato dos ciclos de relaxação não serem sinusoidais, os harmônicos da frequência de relaxação também aparecem nesses espectrogramas. No entanto, em muitas descargas do TCABR, a intensidade das oscilações medidas nas bobinas de Mirnov aumentam significativamente durante o platô, com efeitos sobre a frequência de todas as instabilidades MHD, até mesmo sobre os dentes de serra localizados na região central da coluna. Em todos os casos, observou-se que durante o platô a frequência das ilhas magnéticas coincide com a frequência das oscilações precursoras do dente de serra, apesar de serem duas instabilidades distintas, localizadas em posições radiais muito diferentes. Essa coincidência de frequências possibilitou descrever a evolução em frequência de todas as oscilações detectadas em diversos diagnósticos com base em apenas duas frequências básicas: a dos ciclos de relaxação dente de serra e a das ilhas magnéticas.

Experimental, Numerical and Analytical Studies of the MHD-driven plasma jet, instabilities and waves

This thesis describes a series of experimental, numerical, and analytical studies involving the Caltech magnetohydrodynamically (MHD)-driven plasma jet experiment. The plasma jet is created via a capacitor discharge that powers a magnetized coaxial planar electrodes system. The jet is collimated and accelerated by the MHD forces.

We present three-dimensional ideal MHD finite-volume simulations of the plasma jet experiment using an astrophysical magnetic tower as the baseline model. A compact magnetic energy/helicity injection is exploited in the simulation analogous to both the experiment and to astrophysical situations. Detailed analysis provides a comprehensive description of the interplay of magnetic force, pressure, and flow effects. We delineate both the jet structure and the transition process that converts the injected magnetic energy to other forms.

When the experimental jet is sufficiently long, it undergoes a global kink instability and then a secondary local Rayleigh-Taylor instability caused by lateral acceleration of the kink instability. We present an MHD theory of the Rayleigh-Taylor instability on the cylindrical surface of a plasma flux rope in the presence of a lateral external gravity. The Rayleigh-Taylor instability is found to couple to the classic current-driven instability, resulting in a new type of hybrid instability. The coupled instability, produced by combination of helical magnetic field, curvature of the cylindrical geometry, and lateral gravity, is fundamentally different from the classic magnetic Rayleigh-Taylor instability occurring at a two-dimensional planar interface.

In the experiment, this instability cascade from macro-scale to micro-scale eventually leads to the failure of MHD. When the Rayleigh-Taylor instability becomes nonlinear, it compresses and pinches the plasma jet to a scale smaller than the ion skin depth and triggers a fast magnetic reconnection. We built a specially designed high-speed 3D magnetic probe and successfully detected the high frequency magnetic fluctuations of broadband whistler waves associated with the fast reconnection. The magnetic fluctuations exhibit power-law spectra. The magnetic components of single-frequency whistler waves are found to be circularly polarized regardless of the angle between the wave propagation direction and the background magnetic field.

The effect of temperature dependent viscosity on MHD natural convection flow from an isothermal sphere

Laminar magnetohydrodynamic (MHD) natural convection flow from an isothermal sphere immersed in a fluid with viscosity proportional to linear function of temperature has been studied. The governing boundary layer equations are transformed into a non-dimensional form and the resulting nonlinear system of partial differential equations are reduced to convenient form which are solved numerically by two very efficient methods, namely, (i) Implicit finite difference method together with Keller box scheme and (ii) Direct numerical scheme. Numerical results are presented by velocity and temperature distribution, streamlines and isotherms of the fluid as well as heat transfer characteristics, namely the local skin-friction coefficients and the local heat transfer rate for a wide range of magnetohydrodynamic paramagnet and viscosity-variation parameter.

Effect of MHD and heat generation on natural convection flow in an open square cavity under microgravity condition

Purpose - Thermo-magnetic convection and heat transfer of paramagnetic fluid placed in a micro-gravity condition (g = 0) and under a uniform vertical gradient magnetic field in an open square cavity with three cold sidewalls have been studied numerically. Design/methodology/approach - This magnetic force is proportional to the magnetic susceptibility and the gradient of the square of the magnetic induction. The magnetic susceptibility is inversely proportional to the absolute temperature based on Curie’s law. Thermal convection of a paramagnetic fluid can therefore take place even in zero-gravity environment as a direct consequence of temperature differences occurring within the fluid due to a constant internal heat generation placed within a magnetic field gradient. Findings - Effects of magnetic Rayleigh number, Ra, Prandtl number, Pr, and paramagnetic fluid parameter, m, on the flow pattern and isotherms as well as on the heat absorption are presented graphically. It is found that the heat transfer rate is suppressed in increased of the magnetic Rayleigh number and the paramagnetic fluid parameter for the present investigation. Originality/value - It is possible to control the buoyancy force by using the super conducting magnet. To the best knowledge of the author no literature related to magnetic convection for this configuration is available.

MHD natural convection flow from an isothermal horizontal circular cylinder under consideration of temperature dependent viscosity

Magnetohydrodynamic (MHD) natural convection laminar flow from an iso-thermal horizontal circular cylinder immersed in a fluid with viscosity proportional to a linear function of temperature will be discussed with numerical simulations. The governing boundary layer equations are transformed into a non-dimensional form and the resulting nonlinear system of partial differential equa-tions are reduced to convenient form, which are solved numerically by two very efficient methods, namely, (i) Implicit finite difference method together with Keller box scheme and (ii) Direct numerical scheme. Numerical results are presented by velocity and temperature distributions of the fluid as well as heat transfer characteristics, namely the shearing stress and the local heat transfer rate in terms of the local skin-friction coefficient and the local Nusselt number for a wide range of magnetohydrodynamic parameter, viscosity-variation parameter and viscous dissipation parameter. MHD flow in this geometry with temperature dependent viscosity is absent in the literature. The results obtained from the numerical simulations have been veri-fied by two methodologies.

A phase field study of morphological instabilities in multilayer thin films

We studied the microstructural evolution of multiple layers of elastically stiff films embedded in an elastically soft matrix using a phase field model. The coherent and planar film/matrix interfaces are rendered unstable by the elastic stresses due to a lattice parameter mismatch between the film and matrix phases, resulting in the break-up of the films into particles. With an increasing volume fraction of the stiff phase, the elastic interactions between neighbouring layers lead to: (i) interlayer correlations from an early stage; (ii) a longer wavelength for the maximally growing wave; and therefore (iii) a delayed break-LIP. Further, they promote a crossover in the mode of instability from a predominantly anti-symmetric (in phase) one to a symmetric (out of phase) one. We have computed a stability diagram for the most probable mode of break-up in terms of elastic modulus Mismatch and Volume fraction. We rationalize our results in terms of the initial driving force for destabilization, and corroborate our conclusions using simulations in elastically anisotropic systems.

Instabilities and Waves in Thin Films of Living Fluids

We formulate the thin-film hydrodynamics of a suspension of polar self-driven particles and show that it is prone to several instabilities through the interplay of activity, polarity, and the existence of a free surface. Our approach extends, to self-propelling systems, the work of Ben Amar and Cummings [Phys. Fluids 13 1160 (2001)] on thin-film nematics. Based on our estimates the instabilities should be seen in bacterial suspensions and the lamellipodium, and are potentially relevant to the morphology of biofilms. We suggest several experimental tests of our theory.

Analytical solution of unsteady three-dimensional MHD boundary layer flow and heat transfer due to impulsively stretched plane surface

The unsteady magnetohydrodynamic viscous flow and heat transfer of Newtonian fluids induced by an impulsively stretched plane surface in two lateral directions are studied by using an analytic technique, namely, the homotopy method. The analytic series solution presented here is highly accurate and uniformly valid for all time in the entire region. The effects of the stretching ratio and the magnetic field on the surface shear stresses and heat transfer are studied. The surface shear stresses in x- and y-directions and the surface heat transfer are enchanced by increasing stretching ratio for a fixed value of the magnetic parameter. For a fixed stretching ratio, the surface shear stresses increase with the magnetic parameter, but the heat transfer decreases. The Nusselt number takes longer time to reach the steady state than the skin friction coefficients. There is a smooth transition from the initial unsteady state to the steady state.

MHD Flow of a 2nd Grade Fluid in an Orthogonal Rheometer

The magnetohydrodynamics (MHD) flow of a conducting, homogeneous incompressible Rivlin-Ericksen fluid of second grade contained between two infinite, parallel, insulated disks rotating with the same angular velocity about two noncoincident axes, under the application of a uniform transverse magnetic field, is investigated. This model represents the MHD flow of the fluid in the instrument called an orthogonal rheometer, except for the fact that in the rheometer the rotating plates are necessarily finite. An exact solution of the governing equations of motion is presented. The force components in the x and y directions on the disks are calculated. The effects of magnetic field and the viscoelastic parameter on the forces are discussed in detail.

Heat transfer in a MHD flow between eccentric disks rotating at different speeds

Heat transfer in a MHD flow between two infinite eccentric disks rotating with different speeds is considered when the plates are maintained at different temperatures. The results for the corresponding nonmagnetic case presented wrongly by Banerjee and Borkakati [7] are corrected. It is observed that the eccentric rotation reduces the heat transfer on the disks.

Instabilities of a compressible stratified fluid in horizontal sheared motion

The stability characteristics of a Helmholtz velocity profile in a stably stratified, compressible atmosphere in the presence of a lower boundary are studied. A jump in the Brunt–Väisälä frequency is introduced and the level at which this jump occurs is assumed to be different from the shear zone, to simulate sharp temperature discontinuities in the atmosphere. The results are compared with those of Pellacani, Tebaldi, and Tosi and Lindzen and Rosenthal. In the present configuration, new unstable modes with larger growth rates are found. The wavelengths of the most unstable gravity waves for the parameters pertaining to observed cases of clear air turbulence agree quite closely with the experimental values. Physics of Fluids is copyrighted by The American Institute of Physics

MHD Flow with Heat and Mass Transfer Due to a Point Sink

The magnetofluid dynamic steady incompressible laminar boundary layer flow for a point sink with an applied magnetic field and mass transfer has been studied. The two-point boundary-value problem governed by self-similar equations has been solved numerically. It is observed that the magnetic field increases the skin friction, but reduces the heat transfer and mass flux diffusion. However, the skin friction, heat transfer and mass flux diffusion increase due to suction and the effect of injection is just opposite. Prandtl and Schmidt numbers affect the temperature and concentration, respectively.