Bicoherence in electrostatic turbulence driven by high magnetohydrodynamic activity in Tokamak Chauffage Alfven Bresilien

During some discharges in Tokamak Chauffage Alfven Bresilien [R. M. O. Galvao et al., Plasma Phys. Controlled Fusion 43, 1181 (2001)] high magnetohydrodynamic activity may appear with a peaked frequency spectrum. Whenever this peak occurs, the ambient broadband electrostatic turbulence is remarkably modified, synchronizing into the dominant magnetic fluctuation frequency and presenting high bicoherence in the whole plasma edge with a maximum bicoherence inside the plasma. A phenomenological model is introduced to investigate this driven turbulence bicoherence, consisting of nonlinearly coupled phase-randomized drift modes with time-periodic external driving at the dominant magnetic fluctuation frequency. The bicoherence spectrum of this model can mimic features of the experimental results. (C) 2009 American Institute of Physics. [DOI: 10.1063/1.3099701]

Electrostatic turbulence driven by high magnetohydrodynamic activity in Tokamak Chauffage Alfven Bresilien

In Tokamak Chauffage Alfven Bresilien [R. M. O. Galvao , Plasma Phys. Controlled Fusion 43, 1181 (2001)], high magnetohydrodynamic (MHD) activity may appear spontaneously or during discharges with a voltage biased electrode inserted at the plasma edge. The turbulent electrostatic fluctuations, measured by Langmuir probes, are modulated by Mirnov oscillations presenting a dominant peak with a common frequency around 10 kHz. We report the occurrence of phase locking of the turbulent potential fluctuations driven by MHD activity at this frequency. Using wavelet cross-spectral analysis, we characterized the phase and frequency synchronization in the plasma edge region. We introduced an order parameter to characterize the radial dependence of the phase-locking intensity. (c) 2008 American Institute of Physics.

MAGNETOHYDRODYNAMIC SIMULATIONS OF RECONNECTION AND PARTICLE ACCELERATION: THREE-DIMENSIONAL EFFECTS

Magnetic fields can change their topology through a process known as magnetic reconnection. This process in not only important for understanding the origin and evolution of the large-scale magnetic field, but is seen as a possibly efficient particle accelerator producing cosmic rays mainly through the first-order Fermi process. In this work we study the properties of particle acceleration inserted in reconnection zones and show that the velocity component parallel to the magnetic field of test particles inserted in magnetohydrodynamic (MHD) domains of reconnection without including kinetic effects, such as pressure anisotropy, the Hall term, or anomalous effects, increases exponentially. Also, the acceleration of the perpendicular component is always possible in such models. We find that within contracting magnetic islands or current sheets the particles accelerate predominantly through the first-order Fermi process, as previously described, while outside the current sheets and islands the particles experience mostly drift acceleration due to magnetic field gradients. Considering two-dimensional MHD models without a guide field, we find that the parallel acceleration stops at some level. This saturation effect is, however, removed in the presence of an out-of-plane guide field or in three-dimensional models. Therefore, we stress the importance of the guide field and fully three-dimensional studies for a complete understanding of the process of particle acceleration in astrophysical reconnection environments.

VELOCITY FIELD OF COMPRESSIBLE MAGNETOHYDRODYNAMIC TURBULENCE: WAVELET DECOMPOSITION AND MODE SCALINGS

We study compressible magnetohydrodynamic turbulence, which holds the key to many astrophysical processes, including star formation and cosmic-ray propagation. To account for the variations of the magnetic field in the strongly turbulent fluid, we use wavelet decomposition of the turbulent velocity field into Alfven, slow, and fast modes, which presents an extension of the Cho & Lazarian decomposition approach based on Fourier transforms. The wavelets allow us to follow the variations of the local direction of the magnetic field and therefore improve the quality of the decomposition compared to the Fourier transforms, which are done in the mean field reference frame. For each resulting component, we calculate the spectra and two-point statistics such as longitudinal and transverse structure functions as well as higher order intermittency statistics. In addition, we perform a Helmholtz-Hodge decomposition of the velocity field into incompressible and compressible parts and analyze these components. We find that the turbulence intermittency is different for different components, and we show that the intermittency statistics depend on whether the phenomenon was studied in the global reference frame related to the mean magnetic field or in the frame defined by the local magnetic field. The dependencies of the measures we obtained are different for different components of the velocity; for instance, we show that while the Alfven mode intermittency changes marginally with the Mach number, the intermittency of the fast mode is substantially affected by the change.

Turbulence in collisionless plasmas: statistical analysis from numerical simulations with pressure anisotropy

In recent years, we have experienced increasing interest in the understanding of the physical properties of collisionless plasmas, mostly because of the large number of astrophysical environments (e. g. the intracluster medium (ICM)) containing magnetic fields that are strong enough to be coupled with the ionized gas and characterized by densities sufficiently low to prevent the pressure isotropization with respect to the magnetic line direction. Under these conditions, a new class of kinetic instabilities arises, such as firehose and mirror instabilities, which have been studied extensively in the literature. Their role in the turbulence evolution and cascade process in the presence of pressure anisotropy, however, is still unclear. In this work, we present the first statistical analysis of turbulence in collisionless plasmas using three-dimensional numerical simulations and solving double-isothermal magnetohydrodynamic equations with the Chew-Goldberger-Low laws closure (CGL-MHD). We study models with different initial conditions to account for the firehose and mirror instabilities and to obtain different turbulent regimes. We found that the CGL-MHD subsonic and supersonic turbulences show small differences compared to the MHD models in most cases. However, in the regimes of strong kinetic instabilities, the statistics, i.e. the probability distribution functions (PDFs) of density and velocity, are very different. In subsonic models, the instabilities cause an increase in the dispersion of density, while the dispersion of velocity is increased by a large factor in some cases. Moreover, the spectra of density and velocity show increased power at small scales explained by the high growth rate of the instabilities. Finally, we calculated the structure functions of velocity and density fluctuations in the local reference frame defined by the direction of magnetic lines. The results indicate that in some cases the instabilities significantly increase the anisotropy of fluctuations. These results, even though preliminary and restricted to very specific conditions, show that the physical properties of turbulence in collisionless plasmas, as those found in the ICM, may be very different from what has been largely believed.

Analytical solutions for Tokamak equilibria with reversed toroidal current

In tokamaks, an advanced plasma confinement regime has been investigated with a central hollow electric current with negative density which gives rise to non-nested magnetic surfaces. We present analytical solutions for the magnetohydrodynamic equilibria of this regime in terms of non-orthogonal toroidal polar coordinates. These solutions are obtained for large aspect ratio tokamaks and they are valid for any kind of reversed hollow current density profiles. The zero order solution of the poloidal magnetic flux function describes nested toroidal magnetic surfaces with a magnetic axis displaced due to the toroidal geometry. The first order correction introduces a poloidal field asymmetry and, consequently, magnetic islands arise around the zero order surface with null poloidal magnetic flux gradient. An analytic expression for the magnetic island width is deduced in terms of the equilibrium parameters. We give examples of the equilibrium plasma profiles and islands obtained for a class of current density profile. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624551]

Nonlocal magnetorotational instability

An analytical theory of the nonlocal magnetorotational instability (MRI) is developed for the simplest astrophysical plasma model. It is assumed that the rotation frequency profile has a steplike character, so that there are two regions in which it has constant different values, separated by a narrow transition layer. The surface wave approach is employed to investigate the MRI in this configuration. It is shown that the main regularities of the nonlocal MRI are similar to those of the local instability and that driving the nonaxisymmetric MRI is less effective than the axisymmetric one, also for the case of the nonlocal instability. The existence of nonlocal instabilities in nonmagnetized plasma is predicted. (c) 2008 American Institute of Physics.

EVOLUTION AND LIFETIME OF TRANSIENT CLUMPS IN THE TURBULENT INTERSTELLAR MEDIUM

We study the evolution of dense clumps and provide an argument that the existence of the clumps is not limited by their crossing times. We claim that the lifetimes of the clumps are determined by turbulent motions on a larger scale, and we predict the correlation of clump lifetime with column density. We use numerical simulations to successfully test this relation. In addition, we study the morphological asymmetry and the magnetization of the clumps as functions of their masses.

Thixocasting of an A356 alloy: Fluidity, porosity distribution and thermomechanical fatigue behavior

An extensive set of experiments was performed on a semi-solid A356 alloy in order to assess its flow behavior, mechanical properties, microstructural evolution and porosity level. Three different microstructural conditioning techniques (raw material preparation) were employed: deformation recrystallization, magnetohydrodynamic stirring and low temperature pouring. Measurement of microstructural parameters such as Al-alpha particle size, shape factor, contiguity and entrapped liquid showed a relative equivalency among the various conditioning techniques. It was found that the strongest influence on semi-solid slurry fluidity is exerted by the mould temperature. Tensile properties and porosity levels were measured on a demonstration part produced with different slurry ingate velocities. Results showed similar strength levels among all thixocast samples, a strong correlation between elongation and pore volume fraction and porosity levels much lower than the typical figure for permanent mould or die cast Al-Si alloys. Finally, thermomechanical fatigue tests results were much more favorable to the semi-solid material when compared with the conventionally cast alloy, a result attributed to lower porosity, spheroidal shape of the Al-alpha phase, and refined Si eutectic particles. alpha 2007 Elsevier B.V. All rights reserved.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS. II. THE EFFECTS OF A TILTED MAGNETOSPHERE AND PLANETARY INTERACTIONS

Based on our previous work, we investigate here the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such a configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully three-dimensional (3D) approach. We perform 3D numerical magnetohydrodynamic simulations of stellar winds and study the effects caused by different model parameters, namely the misalignment angle theta(t), the stellar period of rotation, the plasma-beta, and the heating index.. Our simulations take into account the interplay between the wind and the stellar magnetic field during the time evolution. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern. Furthermore, we obtain that by increasing theta(t), the wind velocity increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our 3D, time-dependent wind models allow us to study the interaction of a magnetized wind with a magnetized extrasolar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet`s magnetic field lines and produce electron cyclotron radiation at radio wavelengths. The power released in the interaction depends on the planet`s magnetic field intensity, its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like planet orbiting at 0.05 AU presents a radio power that is similar to 5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. For three selected simulations, we find a variation of the radio power of a factor 1.3-3.7, depending on theta(t). Moreover, we extend the investigation done in Vidotto et al. and analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the timescale tau(w) for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau(w) similar or equal to 100 Myr, for a stellar magnetosphere tilted by theta(t) = 30 degrees, tau(w) ranges from similar to 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction on tau(w) might occur for even larger misalignment angles and/or different wind parameters.

DIFFUSION OF MAGNETIC FIELD AND REMOVAL OF MAGNETIC FLUX FROM CLOUDS VIA TURBULENT RECONNECTION

The diffusion of astrophysical magnetic fields in conducting fluids in the presence of turbulence depends on whether magnetic fields can change their topology via reconnection in highly conducting media. Recent progress in understanding fast magnetic reconnection in the presence of turbulence reassures that the magnetic field behavior in computer simulations and turbulent astrophysical environments is similar, as far as magnetic reconnection is concerned. This makes it meaningful to perform MHD simulations of turbulent flows in order to understand the diffusion of magnetic field in astrophysical environments. Our studies of magnetic field diffusion in turbulent medium reveal interesting new phenomena. First of all, our three-dimensional MHD simulations initiated with anti-correlating magnetic field and gaseous density exhibit at later times a de-correlation of the magnetic field and density, which corresponds well to the observations of the interstellar media. While earlier studies stressed the role of either ambipolar diffusion or time-dependent turbulent fluctuations for de-correlating magnetic field and density, we get the effect of permanent de-correlation with one fluid code, i.e., without invoking ambipolar diffusion. In addition, in the presence of gravity and turbulence, our three-dimensional simulations show the decrease of the magnetic flux-to-mass ratio as the gaseous density at the center of the gravitational potential increases. We observe this effect both in the situations when we start with equilibrium distributions of gas and magnetic field and when we follow the evolution of collapsing dynamically unstable configurations. Thus, the process of turbulent magnetic field removal should be applicable both to quasi-static subcritical molecular clouds and cores and violently collapsing supercritical entities. The increase of the gravitational potential as well as the magnetization of the gas increases the segregation of the mass and magnetic flux in the saturated final state of the simulations, supporting the notion that the reconnection-enabled diffusivity relaxes the magnetic field + gas system in the gravitational field to its minimal energy state. This effect is expected to play an important role in star formation, from its initial stages of concentrating interstellar gas to the final stages of the accretion to the forming protostar. In addition, we benchmark our codes by studying the heat transfer in magnetized compressible fluids and confirm the high rates of turbulent advection of heat obtained in an earlier study.

THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF MAGNETIZED WINDS OF SOLAR-LIKE STARS

By means of self-consistent three-dimensional magnetohydrodynamics (MHD) numerical simulations, we analyze magnetized solar-like stellar winds and their dependence on the plasma-beta parameter (the ratio between thermal and magnetic energy densities). This is the first study to perform such analysis solving the fully ideal three-dimensional MHD equations. We adopt in our simulations a heating parameter described by gamma, which is responsible for the thermal acceleration of the wind. We analyze winds with polar magnetic field intensities ranging from 1 to 20 G. We show that the wind structure presents characteristics that are similar to the solar coronal wind. The steady-state magnetic field topology for all cases is similar, presenting a configuration of helmet streamer-type, with zones of closed field lines and open field lines coexisting. Higher magnetic field intensities lead to faster and hotter winds. For the maximum magnetic intensity simulated of 20 G and solar coronal base density, the wind velocity reaches values of similar to 1000 km s(-1) at r similar to 20r(0) and a maximum temperature of similar to 6 x 10(6) K at r similar to 6r(0). The increase of the field intensity generates a larger ""dead zone"" in the wind, i.e., the closed loops that inhibit matter to escape from latitudes lower than similar to 45 degrees extend farther away from the star. The Lorentz force leads naturally to a latitude-dependent wind. We show that by increasing the density and maintaining B(0) = 20 G the system recover back to slower and cooler winds. For a fixed gamma, we show that the key parameter in determining the wind velocity profile is the beta-parameter at the coronal base. Therefore, there is a group of magnetized flows that would present the same terminal velocity despite its thermal and magnetic energy densities, as long as the plasma-beta parameter is the same. This degeneracy, however, can be removed if we compare other physical parameters of the wind, such as the mass-loss rate. We analyze the influence of gamma in our results and we show that it is also important in determining the wind structure.

Rotation effect on geodesic and zonal flow modes in tokamak plasmas with isothermal magnetic surfaces

The electrostatic geodesic mode oscillations are investigated in rotating large aspect ratio tokamak plasmas with circular isothermal magnetic surfaces. The analysis is carried out within the magnetohydrodynamic model including heat flux to compensate for the non-adiabatic pressure distribution along the magnetic surfaces in plasmas with poloidal rotation. Instead of two standard geodesic modes, three geodesic continua are found. The two higher branches of the geodesic modes have a small frequency up-shift from ordinary geodesic acoustic and sonic modes due to rotation. The lower geodesic continuum is a newzonal flowmode (geodesic Doppler mode) in plasmas with mainly poloidal rotation. Limits to standard geodesic modes are found. Bifurcation of Alfven continuum by geodesic modes at the rational surfaces is also discussed. Due to that, the frequency of combined geodesic continuum extends from the poloidal rotation frequency to the ion-sound band that can have an important role in suppressing plasma turbulence.

Integrable maps with non-trivial topology: application to divertor configurations

We explore a method for constructing two-dimensional area-preserving, integrable maps associated with Hamiltonian systems, with a given set of fixed points and given invariant curves. The method is used to find an integrable Poincare map for the field lines in a large aspect ratio tokamak with a poloidal single-null divertor. The divertor field is a superposition of a magnetohydrodynamic equilibrium with an arbitrarily chosen safety factor profile, with a wire carrying an electric current to create an X-point. This integrable map is perturbed by an impulsive perturbation that describes non-axisymmetric magnetic resonances at the plasma edge. The non-integrable perturbed map is applied to study the structure of the open field lines in the scrape-off layer, reproducing the main transport features obtained by integrating numerically the magnetic field line equations, such as the connection lengths and magnetic footprints on the divertor plate.

High-frequency extensions of magnetorotational instability in astrophysical plasmas

High-frequency extensions of magnetorotational instability driven by the Velikhov effect beyond the standard magnetohydrodynamic (MHD) regime are studied. The existence of the well-known Hall regime and a new electron inertia regime is demonstrated. The electron inertia regime is realized for a lesser plasma magnetization of rotating plasma than that in the Hall regime. It includes the subregime of nonmagnetized electrons. It is shown that, in contrast to the standard MHD regime and the Hall regime, magnetorotational instability in this subregime can be driven only at positive values of dln Omega/dlnr, where Omega is the plasma rotation frequency and r is the radial coordinate. The permittivity of rotating plasma beyond the standard MHD regime, including both the Hall regime and the electron inertia regime, is calculated.