The element-free Galerkin method applied to the study of fully developed magnetohydrodynamic duct flows

In this paper, the meshless method is introduced to magnetohydrodynamics. A numerical scheme based on the element-free Galerkin method is used to solve the laminar steady-state two-dimensional fully developed magnetohydrodynamic flow in a rectangular duct. Accurate and convergent solutions are achieved for low to moderately high Hartmann numbers.

Three-dimensional finite element analysis of MHD duct flow by the penalty function formulation

A numerical scheme based on the Finite Element Method (FEM) is presented to calculate the full solution of a three-dimensional steady magnetohydrodynamic (MHD) flow with moderately high Hartmann numbers and interaction parameters. An incompressible, viscous and electrically conducting liquid-metal is considered. Assuming a low magnetic Reynolds number, the solution method solves the coupled Navier-Stokes and Maxwell's equations through the use of a penalty function method. Results are presented for Hartmann numbers in the range 10(2)-10(3).

Application of a bounded upwinding scheme to complex fluid dynamics problems

This paper is concerned with an overview of upwinding schemes, and further nonlinear applications of a recently introduced high resolution upwind differencing scheme, namely the ADBQUICKEST [V.G. Ferreira, F.A. Kurokawa, R.A.B. Queiroz, M.K. Kaibara, C.M. Oishi, J.A.Cuminato, A.F. Castelo, M.F. Tomé, S. McKee, assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems, International Journal for Numerical Methods in Fluids 60 (2009) 1-26]. The ADBQUICKEST scheme is a new TVD version of the QUICKEST [B.P. Leonard, A stable and accurate convective modeling procedure based on quadratic upstream interpolation, Computer Methods in Applied Mechanics and Engineering 19 (1979) 59-98] for solving nonlinear balance laws. The scheme is based on the concept of NV and TVD formalisms and satisfies a convective boundedness criterion. The accuracy of the scheme is compared with other popularly used convective upwinding schemes (see, for example, Roe (1985) [19], Van Leer (1974) [18] and Arora & Roe (1997) [17]) for solving nonlinear conservation laws (for example, Buckley-Leverett, shallow water and Euler equations). The ADBQUICKEST scheme is then used to solve six types of fluid flow problems of increasing complexity: namely, 2D aerosol filtration by fibrous filters; axisymmetric flow in a tubular membrane; 2D two-phase flow in a fluidized bed; 2D compressible Orszag-Tang MHD vortex; axisymmetric jet onto a flat surface at low Reynolds number and full 3D incompressible flows involving moving free surfaces. The numerical simulations indicate that this convective upwinding scheme is a good generic alternative for solving complex fluid dynamics problems. © 2012.

Magnetic flux transport by turbulent reconnection in astrophysical flows

The role of magnetohydrodynamics (MHD) turbulence in astrophysical environments is still highly debated. An important question that permeates this debate is the transport of magnetic flux. This is particularly important, for instance, in the context of star formation. When clouds collapse gravitationally to form stars, there must be some magnetic flux transport. Otherwise, the newborn stars would have magnetic fields several orders of magnitude larger than the observed ones. Also, the magnetic flux that is dragged in the late stages of the formation of a star can remove all the rotational support from the accretion disc that grows around the protostar. The efficiency of the mechanism that is often invoked to allow transport of magnetic fields at different stages of star formation, namely ambipolar diffusion, has recently been put in check. We discuss here an alternative mechanism for magnetic flux transport which is based on turbulent fast magnetic reconnection. We review recent results from three-dimensional MHD numerical simulations that indicate that this mechanism is very efficient in decoupling and transporting magnetic flux from the inner denser regions to the outskirts of collapsing clouds at different stages of star formation. We discuss this mechanism also in the context of dynamo processes and speculate that it can play a role both in solar dynamo and in accretion disc dynamo processes.

Self-organized criticality in MHD driven plasma edge turbulence

We analyze long-range time correlations and self-similar characteristics of the electrostatic turbulence at the plasma edge and scrape-off layer in the Tokamak Chauffage Alfven Bresillien (TCABR), with low and high Magnetohydrodynamics (MHD) activity. We find evidence of self-organized criticality (SOC), mainly in the region near the tokamak limiter. Comparative analyses of data before and during the MHD activity reveals that during the high mHD activity the Hurst parameter decreases. Finally, we present a cellular automaton whose parameters are adjusted to simulate the analyzed turbulence SOC change with the MHD activity variation. (C) 2011 Published by Elsevier B.V.

THE ROLE OF TURBULENT MAGNETIC RECONNECTION IN THE FORMATION OF ROTATIONALLY SUPPORTED PROTOSTELLAR DISKS

The formation of protostellar disks out of molecular cloud cores is still not fully understood. Under ideal MHD conditions, the removal of angular momentum from the disk progenitor by the typically embedded magnetic field may prevent the formation of a rotationally supported disk during the main protostellar accretion phase of low-mass stars. This has been known as the magnetic braking problem and the most investigated mechanism to alleviate this problem and help remove the excess of magnetic flux during the star formation process, the so-called ambipolar diffusion (AD), has been shown to be not sufficient to weaken the magnetic braking at least at this stage of the disk formation. In this work, motivated by recent progress in the understanding of magnetic reconnection in turbulent environments, we appeal to the diffusion of magnetic field mediated by magnetic reconnection as an alternative mechanism for removing magnetic flux. We investigate numerically this mechanism during the later phases of the protostellar disk formation and show its high efficiency. By means of fully three-dimensional MHD simulations, we show that the diffusivity arising from turbulent magnetic reconnection is able to transport magnetic flux to the outskirts of the disk progenitor at timescales compatible with the collapse, allowing the formation of a rotationally supported disk around the protostar of dimensions similar to 100 AU, with a nearly Keplerian profile in the early accretion phase. Since MHD turbulence is expected to be present in protostellar disks, this is a natural mechanism for removing magnetic flux excess and allowing the formation of these disks. This mechanism dismisses the necessity of postulating a hypothetical increase of the ohmic resistivity as discussed in the literature. Together with our earlier work which showed that magnetic flux removal from molecular cloud cores is very efficient, this work calls for reconsidering the relative role of AD in the processes of star and planet formation.

Analysis of optimal control problems for the incompressible MHD equations and implementation in a finite element multiphysics code

This thesis deals with the study of optimal control problems for the incompressible Magnetohydrodynamics (MHD) equations. Particular attention to these problems arises from several applications in science and engineering, such as fission nuclear reactors with liquid metal coolant and aluminum casting in metallurgy. In such applications it is of great interest to achieve the control on the fluid state variables through the action of the magnetic Lorentz force. In this thesis we investigate a class of boundary optimal control problems, in which the flow is controlled through the boundary conditions of the magnetic field. Due to their complexity, these problems present various challenges in the definition of an adequate solution approach, both from a theoretical and from a computational point of view. In this thesis we propose a new boundary control approach, based on lifting functions of the boundary conditions, which yields both theoretical and numerical advantages. With the introduction of lifting functions, boundary control problems can be formulated as extended distributed problems. We consider a systematic mathematical formulation of these problems in terms of the minimization of a cost functional constrained by the MHD equations. The existence of a solution to the flow equations and to the optimal control problem are shown. The Lagrange multiplier technique is used to derive an optimality system from which candidate solutions for the control problem can be obtained. In order to achieve the numerical solution of this system, a finite element approximation is considered for the discretization together with an appropriate gradient-type algorithm. A finite element object-oriented library has been developed to obtain a parallel and multigrid computational implementation of the optimality system based on a multiphysics approach. Numerical results of two- and three-dimensional computations show that a possible minimum for the control problem can be computed in a robust and accurate manner.

Essential spectrum of systems of singular differential equations

In this paper we develop a new method to determine the essential spectrum of coupled systems of singular differential equations. Applications to problems from magnetohydrodynamics and astrophysics are given.

Kelvin-Helmholtz instabilities at the magnetic cavity boundary of comet 67P/Churyumov-Gerasimenko

We investigate the plasma environment of comet 67P/Churyumov-Gerasimenko, the target of the European Space Agency's Rosetta mission. Rosetta will rendezvous with the comet in 2014 at almost 3.5 AU and follow it all the way to and past perihelion at 1.3 AU. During its journey towards the inner solar system the comet's environment will significantly change. The interaction of the solar wind with a well developed neutral coma leads to the formation of an upstream bow shock and, closer to the comet, the inner shock separating the solar wind, with cometary pick-up ions mass-loaded, from the inner cometary ions which are dragged outward through abundant collisions and charge exchange with the expanding neutral gas. As a consequence the interplanetary magnetic field is prevented from penetrating the innermost region of the comet, the so-called magnetic cavity. We use our magnetohydrodynamics model BATSRUS (Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme) to simulate the solar wind - comet interaction. The model includes photoionization, ion-electron recombination, and charge exchange. Under certain conditions our model predicts an unstable plasma flow at the inner shock. We show that the plasma shear flow around the magnetic cavity can lead to Kelvin-Helmholtz instabilities. We investigate the onset of this phenomenon with change of heliocentric distance and furthermore show that a previously stable magnetic cavity boundary can become unstable when the neutral gas is predominately released from the dayside of the comet.

Plasma environment of a weak comet – Predictions for Comet 67P/Churyumov–Gerasimenko from multifluid-MHD and Hybrid models

The interaction of a comet with the solar wind undergoes various stages as the comet’s activity varies along its orbit. For a comet like 67P/Churyumov–Gerasimenko, the target comet of ESA’s Rosetta mission, the various features include the formation of a Mach cone, the bow shock, and close to perihelion even a diamagnetic cavity. There are different approaches to simulate this complex interplay between the solar wind and the comet’s extended neutral gas coma which include magnetohydrodynamics (MHD) and hybrid-type models. The first treats the plasma as fluids (one fluid in basic single fluid MHD) and the latter treats the ions as individual particles under the influence of the local electric and magnetic fields. The electrons are treated as a charge-neutralizing fluid in both cases. Given the different approaches both models yield different results, in particular for a low production rate comet. In this paper we will show that these differences can be reduced when using a multifluid instead of a single-fluid MHD model and increase the resolution of the Hybrid model. We will show that some major features obtained with a hybrid type approach like the gyration of the cometary heavy ions and the formation of the Mach cone can be partially reproduced with the multifluid-type model.

Force-free twisted magnetospheres of neutron stars

Context. The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. Aims. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. Methods. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. Results. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. Conclusions. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.

A new code for the Hall-driven magnetic evolution of neutron stars

Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.

Space charge limited current flow from a plane, plasma cathode,

"Prepared for the Air Force Ballistic Missile Division, Headquarters Air Research and Development Command, under contract AF 04(647)-309, Thermonuclear Propulsion Research."

On one-dimensional channel flow in the presence of a magnetic field.

"Prepared for the Air Force Ballistic Missile Division, Headquarters Air Research and Development Command, under contract AF 04(647)-309 advanced propulsion systems."

Radiation by plasma oscillations in a bounded plasma in a magnetic field.

At head of title: PRL-9-23.