Normalized to: Ruediger, G.
[1]
oai:arXiv.org:1610.05469 [pdf] - 1498509
Laboratory experiments and numerical simulations on magnetic
instabilities
Submitted: 2016-10-18
Magnetic fields of planets, stars and galaxies are generated by
self-excitation in moving electrically conducting fluids. Once produced,
magnetic fields can play an active role in cosmic structure formation by
destabilizing rotational flows that would be otherwise hydrodynamically stable.
For a long time, both hydromagnetic dynamo action as well as magnetically
triggered flow instabilities had been the subject of purely theoretical
research. Meanwhile, however, the dynamo effect has been observed in
large-scale liquid sodium experiments in Riga, Karlsruhe and Cadarache. In this
paper, we summarize the results of some smaller liquid metal experiments
devoted to various magnetic instabilities such as the helical and the azimuthal
magnetorotational instability, the Tayler instability, and the different
instabilities that appear in a magnetized spherical Couette flow. We conclude
with an outlook on a large scale Tayler-Couette experiment using liquid sodium,
and on the prospects to observe magnetically triggered instabilities of flows
with positive shear.
[2]
oai:arXiv.org:1601.03877 [pdf] - 1447663
Current helicity and electromotive force of magnetoconvection influenced
by helical background fields
Submitted: 2016-01-15
Motivated by the empirical finding that the known hemispheric rules for the
current helicity at the solar surface are not strict, the excitation of
small-scale current helicity by the influence of a large-scale helical magnetic
background fields on nonrotating magnetoconvection is demonstrated. It is shown
within a quasilinear analytic theory of driven turbulence and by nonlinear
simulations of magnetoconvection that the resulting small-scale current
helicity has the same sign as the large-scale current helicity while the ratio
of both pseudo-scalars is of the order of the magnetic Reynolds number of the
convection. The same models do not provide finite values of the small-scale
kinetic helicity. On the other hand, a turbulence-induced electromotive force
is produced including the diamagnetic pumping term as well as the eddy
diffusivity but no alpha effect. It is thus argued that the relations by
Pouquet & Patterson (1978) and Keinigs (1983) for the simultaneous existence of
small-scale current helicity and alpha effect do not hold for the considered
model of nonrotating magnetoconvection. Calculations for various values of the
magnetic Prandtl number demonstrate that for the considered diffusivities the
current helicity grows for growing magnetic Reynolds number which is not true
for the velocity of the diamagnetic pumping -- in agreement with the results of
the quasilinear analytical approximation.
[3]
oai:arXiv.org:1404.3562 [pdf] - 969762
The effective magnetic Prandtl number in magnetized and differentially
rotating stellar radiation zones
Submitted: 2014-04-14, last modified: 2015-03-24
With application to inner stellar radiative zones, a linear theory is used to
analyze the instability of a dipole-parity toroidal background field, in the
presence of density stratification, differential rotation, and realistically
small Prandtl numbers. The physical parameters are the normalized latitudinal
shear $a$ and the normalized field amplitude $b$. Only the solutions for the
wavelengths with the maximal growth rates are considered. If these scales are
combined to the radial values of velocity, one finds that the (very small)
radial velocity only depends slightly on $a$ and $b$, so that it can be used as
the free parameter of the eigenvalue system.
The resulting instability-generated tensors of magnetic diffusivity and eddy
viscosity are highly anisotropic. The eddy diffusivity in latitudinal direction
exceeds the eddy diffusivity in radial direction by orders of magnitude. Its
latitudinal profile shows a strong concentration toward the poles which is also
true for the effective viscosity which has been calculated via the angular
momentum transport of the instability pattern. The resulting effective magnetic
Prandtl number reaches values of $O(10^2)$, so that the differential rotation
decays much faster than the toroidal background field, which is {the} necessary
condition to explain the observed slow rotation of the early red-giant and
sub-giant cores by means of magnetic instabilities.
[4]
oai:arXiv.org:1404.4288 [pdf] - 1209006
The angular momentum transport by unstable toroidal magnetic fields
Submitted: 2014-04-16, last modified: 2014-10-30
We demonstrate with a nonlinear MHD code that angular momentum can be
transported due to the magnetic instability of toroidal fields under the
influence of differential rotation, and that the resulting effective viscosity
may be high enough to explain the almost rigid-body rotation observed in
radiative stellar cores. Only stationary current-free fields and only those
combinations of rotation rates and magnetic field amplitudes which provide
maximal numerical values of the viscosity are considered. We find that the
dimensionless ratio of the effective over molecular viscosity, $\nu_T/\nu$;,
linearly grows with the Reynolds number of the rotating fluid multiplied with
the square-root of the magnetic Prandtl number - which is of order unity for
the considered red sub-giant KIC 7341231.
For the considered interval of magnetic Reynolds numbers - which is
restricted by numerical constraints of the nonlinear MHD code - there is a
remarkable influence of the magnetic Prandtl number on the relative importance
of the contributions of the Reynolds stress and the Maxwell stress to the total
viscosity, which is magnetically dominated only for Pm $\gtrsim$ 0.5. We also
find that the magnetized plasma behaves as a non-Newtonian fluid, i.e. the
resulting effective viscosity depends on the shear in the rotation law. The
decay time of the differential rotation thus depends on its shear and becomes
longer and longer during the spin-down of a stellar core.
[5]
oai:arXiv.org:1408.3837 [pdf] - 1216347
The existence of the Lambda effect in the solar convection zone
indicated by SDO observations
Submitted: 2014-08-17
The empirical finding with data from the Solar Dynamics Observatory (SDO) of
positive (negative) horizontal Reynolds stress at the northern (southern)
hemisphere for solar giant cells (Hathaway et al. 2013) is discussed for its
consequences for the theory of the solar/stellar differential rotation. Solving
the nonlinear Reynolds equation for the angular velocity under neglect of the
meridional circulation we show that the horizontal Reynolds stress of the
northern hemisphere is always negative at the surface but it is positive in the
bulk of the solar convection zone by the action of the Lambda effect. The
Lambda effect, which describes the angular momentum transport of rigidly
rotating anisotropic turbulence and which avoids a rigid-body rotation of the
convection zones, is in horizontal direction of cubic power in the rotation
rate $\Omega$ and it is always equatorwards directed. Theories without Lambda
effect which may also provide the observed solar rotation law only by the
action of a meridional circulation lead to a horizontal Reynolds stress with
the opposite sign as observed.
[6]
oai:arXiv.org:1407.0240 [pdf] - 1284715
Diffusive MHD Instabilities: Beyond the Chandrasekhar Theorem
Submitted: 2014-07-01
The magnetohydrodynamic stability of axially unbounded cylindrical flows is
considered which contain a toroidal magnetic background field with the same
radial profile as the linear azimuthal velocity. Chandrasekhar (1956) has shown
for ideal fluids the stability of this configuration if the Alfven velocity of
the field equals the velocity of the background flow. It is demonstrated for
magnetized Taylor-Couette flows at the Rayleigh line, however, that for finite
diffusivity such flows become unstable against nonaxisymmetric perturbations
where the critical magnetic Reynolds number of the rotation rate does not
depend on the magnetic Prandtl number Pm if Pm much << 1. In order to study
this new diffusive azimuthal magnetorotational instability, flows and fields
with the same radial profile but with different amplitudes are considered. For
Pm << 1 the instability domain with the weakest fields and the slowest rotation
rates lies below the Chandrasekhar line of equal amplitudes for Alfven velocity
and rotation velocity. We find that then the lines of marginal instability
scale with the Reynolds number and the Hartmann number. The minimum values of
the field strength and the rotation rate which are needed for the instability
(slightly) grow for more and more flat rotation. Finally, the corresponding
electric current of the background field becomes so strong that the Tayler
instability (which even exists without rotation) also appears in the
bifurcation map at small Hartmann numbers displacing after all the azimuthal
magnetorotational instability.
[7]
oai:arXiv.org:1201.0652 [pdf] - 862289
The alpha-effect in a turbulent liquid-metal plane Couette flow
Submitted: 2012-01-03, last modified: 2014-02-10
We calculate the mean electromotive force in plane Couette flows of a
nonrotating conducting fluid under the influence of a large-scale magnetic
field for driven turbulence. A vertical stratification of the turbulence
intensity results in an alpha effect owing to the presence of horizontal shear.
Here we discuss the possibility of an experimental determination of the
components of the alpha tensor using both quasilinear theory and nonlinear
numerical simulations. For magnetic Prandtl numbers of the order of unity, we
find that in the high-conductivity limit the alpha effect in the direction of
the flow clearly exceeds the component in spanwise direction. In this limit,
alpha runs linearly with the magnetic Reynolds number Rm while in the
low-conductivity limit it runs with the product Rm*Re, where Re is the kinetic
Reynolds number so that for given Rm the alpha effect grows with decreasing
magnetic Prandtl number.
For the small magnetic Prandtl numbers of liquid metals, a common value for
the horizontal elements of the alpha tensor appears, which makes it unimportant
whether the alpha effect is measured in the spanwise or streamwise directions.
The resulting effect should lead to an observable voltage in both directions of
about 0.5 mV for magnetic fields of 1 kgauss and velocity fluctuations of about
1 m/s in a channel of 50 cm height (independent of its width).
[8]
oai:arXiv.org:1303.4621 [pdf] - 1165374
Astrophysical and experimental implications from the magnetorotational
instability of toroidal fields
Submitted: 2013-03-19, last modified: 2013-11-06
The interaction of differential rotation and toroidal fields that are
current-free in the gap between two corotating axially unbounded cylinders is
considered. It is shown that nonaxisymmetric perturbations are unstable if the
rotation rate and Alfv\'en frequency of the field are of the same order, almost
independent of the magnetic Prandtl number Pm. For the very steep rotation law
\Omega\propto R^{-2} (the Rayleigh limit) and for small Pm the threshold values
of rotation and field for this Azimuthal MagnetoRotational Instability (AMRI)
scale with the ordinary Reynolds number and the Hartmann number, resp. A
laboratory experiment with liquid metals like sodium or gallium in a
Taylor-Couette container has been designed on the basis of this finding. For
fluids with more flat rotation laws the Reynolds number and the Hartmann number
are no longer typical quantities for the instability. For the weakly nonlinear
system the numerical values of the kinetic energy and the magnetic energy are
derived for magnetic Prandtl numbers \leq 1. We find that the magnetic energy
grows monotonically with the magnetic Reynolds number Rm, while the kinetic
energy grows with Rm/\sqrt{Pm}. The resulting turbulent Schmidt number, as the
ratio of the `eddy' viscosity and the diffusion coefficient of a passive scalar
(such as lithium) is of order 20 for Pm=1, but for small Pm it drops to order
unity. Hence, in a stellar core with fossil fields and steep rotation law the
transport of angular momentum by AMRI is always accompanied by an intense
mixing of the plasma, until the rotation becomes rigid.
[9]
oai:arXiv.org:1201.2318 [pdf] - 1092873
Critical fields and growth rates of the Tayler instability as probed by
a columnar gallium experiment
Submitted: 2012-01-11, last modified: 2012-07-16
Many astrophysical phenomena (such as the slow rotation of neutron stars or
the rigid rotation of the solar core) can be explained by the action of the
Tayler instability of toroidal magnetic fields in the radiative zones of stars.
In order to place the theory of this instability on a safe fundament it has
been realized in a laboratory experiment measuring the critical field strength,
the growth rates as well as the shape of the supercritical modes. A strong
electrical current flows through a liquid-metal confined in a resting columnar
container with an insulating outer cylinder. As the very small magnetic Prandtl
number of the gallium-indium-tin alloy does not influence the critical Hartmann
number of the field amplitudes, the electric currents for marginal instability
can also be computed with direct numerical simulations. The results of this
theoretical concept are confirmed by the experiment. Also the predicted growth
rates of the order of minutes for the nonaxisymmetric perturbations are
certified by the measurements. That they do not directly depend on the size of
the experiment is shown as a consequence of the weakness of the applied fields
and the absence of rotation.
[10]
oai:arXiv.org:1202.1429 [pdf] - 1116444
The cross helicity at the solar surface by simulations and observations
Submitted: 2012-02-07, last modified: 2012-06-27
The quasilinear mean-field theory for driven MHD turbulence leads to the
result that the observed cross helicity <u \cdot b> may directly yield the
magnetic eddy diffusivity \eta_{T} of the quiet Sun. In order to model the
cross helicity at the solar surface, magnetoconvection under the presence of a
vertical large-scale magnetic field is simulated with the nonlinear MHD code
NIRVANA. The very robust result of the calculations is that < u_z b_z >\simeq 2
< u\cdot b> independent of the applied magnetic field amplitude. The
correlation coefficient for the cross helicity is about 10%. Of similar
robustness is the finding that the rms value of the magnetic perturbations
exceeds the mean-field amplitude (only) by a factor of five. The characteristic
helicity speed u_{\eta} as the ratio of the eddy diffusivity and the density
scale height for an isothermal sound velocity of 6.6 km/s proves to be 1 km/s
for weak fields. This value well coincides with empirical results obtained from
the data of the HINODE satellite and the Swedish 1-m Solar Telescope (SST)
providing the cross helicity component < u_z b_z >. Both simulations and
observations thus lead to a numerical value of \eta_{T} \simeq 10^12 cm^2 /s as
characteristic for the surface of the quiet Sun.
[11]
oai:arXiv.org:1203.4348 [pdf] - 1117417
Linear stability analysis of the Hall magnetorotational instability in a
spherical domain
Submitted: 2012-03-20, last modified: 2012-03-21
We investigate the stability of the Hall-MHD system and determine its
importance for neutron stars at their birth, when they still consist of
differentially rotating plasma permeated by extremely strong magnetic fields.
We solve the linearised Hall-MHD equations in a spherical shell threaded by a
homogeneous magnetic field. With the fluid/flow coupling and the Hall effect
included, the magnetorotational instability and the Hall effect are both acting
together. Results differ for magnetic fields aligned with the rotation axis and
anti-parallel magnetic fields. For a positive alignment of the magnetic field
the instability grows on a rotational time-scale for any sufficiently large
magnetic Reynolds number. Even the magnetic fields which are stable against the
MRI due to the magnetic diffusion are now susceptible to the shear-Hall
instability. In contrast, the negative alignment places strong restrictions on
the growth and the magnitude of the fields, hindering the effectiveness of the
Hall-MRI. While non-axisymmetric modes of the MRI can be suppressed by strong
enough rotation, there is no such restriction when the Hall effect is present.
The implications for the magnitude and the topology of the magnetic field of a
young neutron star may be significant.
[12]
oai:arXiv.org:1109.3345 [pdf] - 1084119
Suppression of the large-scale Lorentz force by turbulence
Submitted: 2011-09-15, last modified: 2012-01-03
The components of the total stress tensor (Reynolds stress plus Maxwell
stress) are computed within the quasilinear approximation for a driven
turbulence influenced by a large-scale magnetic background field. The
conducting fluid has an arbitrary magnetic Prandtl number and the turbulence
without the background field is assumed as homogeneous and isotropic with a
free Strouhal number St. The total large-scale magnetic tension is always
reduced by the turbulence with the possibility of a `catastrophic quenching'
for large magnetic Reynolds number Rm so that even its sign is reversed. The
total magnetic pressure is enhanced by turbulence with short correlation time
(`white noise') but it is reduced by turbulence with long correlation time.
Also in this case the sign of the total pressure may reverse but only for
special turbulences with sufficiently large St> 1.
The turbulence-induced terms of the stress tensor are suppressed by strong
magnetic fields. For the tension term this quenching grows with the square of
the Hartmann number of the magnetic field. For microscopic (i.e. small)
diffusivity values the magnetic tension term becomes thus highly quenched even
for field amplitudes much smaller than their equipartition value. In the
opposite case of large-eddy simulations the magnetic quenching is only mild but
then also the turbulence-induced Maxwell tensor components for weak fields
remain rather small.
[13]
oai:arXiv.org:1112.2103 [pdf] - 1092272
Experimental evidence for Tayler instability in a liquid metal column
Submitted: 2011-12-09
In the current-driven, kink-type Tayler instability (TI) a sufficiently
strong azimuthal magnetic field becomes unstable against non-axisymmetric
perturbations. The TI has been discussed as a possible ingredient of the solar
dynamo mechanism and a source of the helical structures in cosmic jets. It is
also considered as a size limiting factor for liquid metal batteries. We report
on a liquid metal TI experiment using a cylindrical column of the eutectic
alloy GaInSn to which electrical currents of up to 8 kA are applied. We present
results of external magnetic field measurements that indicate the occurrence of
the TI in good agreement with numerical predictions. The interference of TI
with the competing large scale convection, resulting from Joule heating, is
also discussed.
[14]
oai:arXiv.org:1110.4757 [pdf] - 1085047
Differential rotation and meridional flow on the lower zero age main
sequence: Reynolds stress versus baroclinic flow
Submitted: 2011-10-21
We study the variation of surface differential rotation and meridional flow
along the lower part of the zero age main sequence (ZAMS). We first compute a
sequence of stellar models with masses from 0.3 to 1.5 solar masses. We then
construct mean field models of their outer convection zones and compute
differential rotation and meridional flows by solving the Reynolds equation
with transport coefficients from the second order correlation approximation.
For a fixed rotation period of 2.5 d we find a strong dependence of the surface
differential rotation on the effective temperature with weak surface shear for
M dwarfs and very large values for F stars. The increase with effective
temperature is modest below 6000 K but very steep above 6000 K. The meridional
flow shows a similar variation with temperature but the increase with
temperature is not quite so steep. Both the surface rotation and the meridional
circulation are solar-type over the entire temperature range. We also study the
dependence of differential rotation and meridional flow on the rotation period
for masses. from 0.3 to 1.1 solar masses. The variation of the differential
rotation with period is weak except for very rapid rotation. The meridional
flow shows a systematic increase of the flow speed with the rotation rate.
Numerical experiments in which either the $\Lambda$ effect is dropped in the
Reynolds stress or the baroclinic term in the equation of motion is cancelled
show that for effective temperatures below 6000 K the Reynolds stress is the
dominant driver of differential rotation.
[15]
oai:arXiv.org:1107.2548 [pdf] - 1077940
Helicity and dynamo action in magnetized stellar radiation zones
Submitted: 2011-07-13
Helicity and \alpha effect driven by the nonaxisymmetric Tayler instability
of toroidal magnetic fields in stellar radiation zones are computed. In the
linear approximation a purely toroidal field always excites pairs of modes with
identical growth rates but with opposite helicity so that the net helicity
vanishes. If the magnetic background field has a helical structure by an extra
(weak) poloidal component then one of the modes dominates producing a net
kinetic helicity anticorrelated to the current helicity of the background
field.
The mean electromotive force is computed with the result that the \alpha
effect by the most rapidly growing mode has the same sign as the current
helicity of the background field. The \alpha effect is found as too small to
drive an \alpha^{2} dynamo but the excitation conditions for an \alpha\Omega
dynamo can be fulfilled for weak poloidal fields. Moreover, if the dynamo
produces its own \alpha effect by the magnetic instability then problems with
its sign do not arise. For all cases, however, the \alpha effect shows an
extremely strong concentration to the poles so that a possible \alpha\Omega
dynamo might only work at the polar regions. Hence, the results of our linear
theory lead to a new topological problem for the existence of large-scale
dynamos in stellar radiation zones on the basis of the current-driven
instability of toroidal fields.
[16]
oai:arXiv.org:1101.0535 [pdf] - 1042843
Helicity and alpha-effect by current-driven instabilities of helical
magnetic fields
Submitted: 2011-01-03, last modified: 2011-01-04
Helical magnetic background fields with adjustable pitch angle are imposed on
a conducting fluid in a differentially rotating cylindrical container. The
small-scale kinetic and current helicities are calculated for various field
geometries, and shown to have the opposite sign as the helicity of the
large-scale field. These helicities and also the corresponding $\alpha$-effect
scale with the current helicity of the background field. The $\alpha$-tensor is
highly anisotropic as the components $\alpha_{\phi\phi}$ and $\alpha_{zz}$ have
opposite signs. The amplitudes of the azimuthal $\alpha$-effect computed with
the cylindrical 3D MHD code are so small that the operation of an
$\alpha\Omega$ dynamo on the basis of the current-driven, kink-type
instabilities of toroidal fields is highly questionable. In any case the low
value of the $\alpha$-effect would lead to very long growth times of a dynamo
in the radiation zone of the Sun and early-type stars of the order of
mega-years.
[17]
oai:arXiv.org:1012.0720 [pdf] - 1042355
On radiation-zone dynamos
Submitted: 2010-12-03
It is shown that the magnetic current-driven (`kink-type') instability
produces flow and field patterns with helicity and even with \alpha-effect but
only if the magnetic background field possesses non-vanishing current helicity
\bar{\vec{B}}\cdot curl \bar{\vec{B}} by itself. Fields with positive
large-scale current helicity lead to negative small-scale kinetic helicity. The
resulting \alpha-effect is positive. These results are very strict for
cylindric setups without <i>z/I>-dependence of the background fields. The sign
rules also hold for the more complicated cases in spheres where the toroidal
fields are the result of the action of differential rotation (induced from
fossil poloidal fields) at least for the case that the global rotation is
switched off after the onset of the instability.
[18]
oai:arXiv.org:1011.5726 [pdf] - 1042197
The Tayler instability of toroidal magnetic fields in a columnar gallium
experiment
Submitted: 2010-11-26
The nonaxisymmetric Tayler instability of toroidal magnetic fields due to
axial electric currents is studied for conducting incompressible fluids between
two coaxial cylinders without endplates. The inner cylinder is considered as so
thin that even the limit of R_in \to 0 can be computed. The magnetic Prandtl
number is varied over many orders of magnitudes but the azimuthal mode number
of the perturbations is fixed to m=1. In the linear approximation the critical
magnetic field amplitudes and the growth rates of the instability are
determined for both resting and rotating cylinders. Without rotation the
critical Hartmann numbers do {\em not} depend on the magnetic Prandtl number
but this is not true for the growth rates. For given product of viscosity and
magnetic diffusivity the growth rates for small and large magnetic Prandtl
number are much smaller than those for Pm=1. For gallium under the influence of
a magnetic field at the outer cylinder of 1 kG the resulting growth time is 5
s. The minimum electric current through a container of 10 cm diameter to excite
the kink-type instability is 3.20 kA. For a rotating container both the
critical magnetic field and the related growth times are larger than for the
resting column.
[19]
oai:arXiv.org:1011.2684 [pdf] - 1041903
Magnetic fields of Ap stars as a result of an instability
Submitted: 2010-11-11
Ap star magnetism is often attributed to fossil magnetic fields which have
not changed much since the pre-main-sequence epoch of the stars. Stable
magnetic field configurations are known which could persist probably for the
entire main-sequence life of the star, but they may not show the complexity and
diversity exhibited by the Ap stars observed. We suggest that the Ap star
magnetism is not a result of stable configurations, but is the result of an
instability based on strong toroidal magnetic fields buried in the stars. The
highly nonaxisymmetric remainders of the instability are reminiscent of the
diversity of fields seen on Ap stars. The strengths of these remnant magnetic
fields is actually between a few per cent up to considerable fractions of the
internal toroidal field; this means field strengths of the order of kGauss
being compatible with what is observed. The magnetic fields emerge at the
surface rather quickly; rough estimates deliver time-scales of the order of a
few years. Since rotation stabilizes the instability, normal A stars may still
host considerable, invisible toroidal magnetic fields.
[20]
oai:arXiv.org:1010.3899 [pdf] - 1041395
Amplification and stability of magnetic fields and dynamo effect in
young A stars
Submitted: 2010-10-19
This study is concerned with the early evolution of magnetic fields and
differential rotation of intermediate-mass stars which may evolve into Ap
stars. We report on simulations of the interplay of differential rotation and
magnetic fields, the stability limits and non-linear evolution of such
configurations, and the prospects of dynamo action from the unstable cases. The
axisymmetric problem delivers a balance between field amplification and
back-reaction of the magnetic field on the differential rotation. The
non-axisymmetric case involves also the Tayler instability of the amplified
toroidal fields. We consider limits for field amplification and apply these to
young A stars. Apart from its application to Ap stars, the instability is
scrutinized for the fundamental possibility of a dynamo. We are not looking for
a dynamo as an explanation for the Ap star phenomenon. The kinetic helicity is
concentrated near the tangent cylinder of the inner sphere of the computational
domain and is negative in the northern hemisphere. This appears to be a
ubiquitous effect not special to the Tayler instability. The latter is actually
connected with a positive current helicity in the bulk of the spherical shell
giving rise to a small, but non-vanishing alpha-effect in non-linear evolution
of the instability.
[21]
oai:arXiv.org:1009.2315 [pdf] - 1034907
The pinch-type instability of helical magnetic fields
Submitted: 2010-09-13
To find out whether toroidal field can stably exist in galaxies the
current-driven instability of toroidal magnetic fields is considered under the
influence of an axial magnetic field component and under the influence of both
rigid and differential rotation. The MHD equations are solved in a simplified
model with cylindric geometry. We assume the axial field as uniform and the
fluid as incompressible. The stability of a toroidal magnetic field is strongly
influenced by uniform axial magnetic fields. If both field components are of
the same order of magnitude then the instability is slightly supported and
modes with m>1 dominate. If the axial field even dominates the most unstable
modes have again m>1 but the field is strongly stabilized. All modes are
suppressed by a fast rigid rotation where the m=1 mode maximally resists. Just
this mode becomes best re-animated for \Omega > \Omega^A (\Omega^A the Alfven
frequency) if the rotation has a negative shear. -- Strong indication has been
found for a stabilization of the nonaxisymmetric modes for fluids with small
magnetic Prandtl number if they are unstable for Pm=1. For rotating fluids the
higher modes with m>1 do not play an important role in the linear theory. In
the light of our results galactic fields should be marginally unstable against
perturbations with m<= 1. The corresponding growth rates are of the order of
the rotation period of the inner part of the galaxy.
[22]
oai:arXiv.org:1002.3037 [pdf] - 1025201
The exceptional Herbig Ae star HD101412: The first detection of resolved
magnetically split lines and the presence of chemical spots in a Herbig star
Submitted: 2010-02-16
We obtained high-resolution, high signal-to-noise UVES and a few lower
quality HARPS spectra revealing the presence of resolved magnetically split
lines. HD101412 is the first Herbig Ae star for which the rotational Doppler
effect was found to be small in comparison to the magnetic splitting. The
measured mean magnetic field modulus varies from 2.5 to 3.5kG, while the mean
quadratic field was found to vary in the range of 3.5 to 4.8kG. To determine
the period of variations, we used radial velocity, equivalent width, line
width, and line asymmetry measurements of variable spectral lines of several
elements, as well as magnetic field measurements. The most pronounced
variability was detected for spectral lines of He I and the iron peak elements,
whereas the spectral lines of CNO elements are only slightly variable. From
spectral variations and magnetic field measurements we derived a potential
rotation period P_rot=13.86d, which has to be proven in future studies with a
larger number of observations. It is the first time that the presence of
element spots is detected on the surface of a Herbig Ae/Be star. Our previous
study of Herbig Ae stars revealed a trend towards stronger magnetic fields for
younger Herbig Ae stars, confirmed by statistical tests. This is in contrast to
a few other (non-statistical) studies claiming that magnetic Herbig Ae stars
are progenitors of the magnetic Ap stars. New developments in MHD theory show
that the measured magnetic field strengths are compatible with a current-driven
instability of toroidal fields generated by differential rotation in the
stellar interior. This explanation for magnetic intermediate-mass stars could
be an alternative to a frozen-in fossil field.
[23]
oai:arXiv.org:1001.1895 [pdf] - 1019134
Dissipative Taylor-Couette flows under the influence of helical magnetic
fields
Submitted: 2010-01-12
The linear stability of MHD Taylor-Couette flows in axially unbounded
cylinders is considered, for magnetic Prandtl number unity. Magnetic fields
varying from purely axial to purely azimuthal are imposed, with a general
helical field parameterized by \beta=B_\phi/B_z. We map out the transition from
the standard MRI for \beta=0 to the nonaxisymmetric Azimuthal MagnetoRotational
Instability (AMRI) for \beta\to \infty. For finite \beta, positive and negative
wave numbers m, corresponding to right and left spirals, are no longer
identical. The transition from \beta=0 to \beta\to\infty includes all the
possible forms of MRI with axisymmetric and nonaxisymmetric modes. For the
nonaxisymmetric modes, the most unstable mode spirals in the opposite direction
to the background field. The standard (\beta=0) MRI is axisymmetric for weak
fields (including the instability with the lowest Reynolds number) but is
nonaxisymmetric for stronger fields. If the azimuthal field is due in part to
an axial current flowing through the fluid itself (and not just along the
central axis), then it is also unstable to the nonaxisymmetric Tayler
instability, which is most effective without rotation. For large \beta this
instability has wavenumber m=1, whereas for \beta\simeq 1 m=2 is most unstable.
The most unstable mode spirals in the same direction as the background field.
[24]
oai:arXiv.org:0908.1841 [pdf] - 1003435
The kink-type instability of toroidal stellar magnetic fields with
thermal diffusion
Submitted: 2009-08-13, last modified: 2009-12-04
The stability of toroidal magnetic fields in rotating radiative stellar zones
is studied for realistic values of both the Prandtl numbers. The two considered
models for the magnetic geometry represent fields with odd and even symmetry
with respect to the equator. In the linear theory in Boussinesq approximation
the resulting complex eigenfrequency (including growth rate and drift rate) are
calculated for a given radial wavenumber of a nonaxisymmetric perturbation with
m=1. The ratio of the Alfven frequency, \Omega_A, to the rate of the basic
rotation, \Omega, controls the eigenfrequency of the solution. For strong
fields with \Omega_A > \Omega the solutions do not feel the thermal diffusion.
The growth rate runs with \Omega_A and the drift rate is close to -\Omega so
that the magnetic pattern will rest in the laboratory system. For weaker fields
with \Omega_A < \Omega the growth rate strongly depends on the thermal
conductivity. For fields with dipolar parity and for typical values of the heat
conductivity the resulting very small growth rates are almost identical with
those for vanishing gravity. For fields with dipolar symmetry the differential
rotation of any stellar radiative zone (like the solar tachocline) is shown as
basically stabilizing the instability independent of the sign of the shear.
Finally, the current-driven kink-type instability of a toroidal background
field is proposed as a model for the magnetism of Ap stars. The recent
observation of a lower magnetic field treshold of about 300 Gauss for Ap stars
is understood as corresponding to the minimum magnetic field producing the
instability.
[25]
oai:arXiv.org:0911.5603 [pdf] - 1527983
Dissipative and nonaxisymmetric standard-MRI in Kepler disks
Submitted: 2009-11-30
Deviations from axial symmetry are necessary to maintain self-sustained
MRI-turbulence. We define the parameters region where nonaxisymmetric MRI is
excited and study dependence of the unstable modes structure and growth rates
on the relevant parameters. We solve numerically the linear eigenvalue problem
for global axisymmetric and nonaxisymmetric modes of standard-MRI in Keplerian
disks with finite diffusion. For small magnetic Prandtl number the microscopic
viscosity completely drops out from the analysis so that the stability maps and
the growth rates expressed in terms of the magnetic Reynolds number Rm and the
Lundquist number S do not depend on the magnetic Prandtl number Pm. The minimum
magnetic field for onset of nonaxisymmetric MRI grows with Rm. For given S all
nonaxisymmetric modes disappear for sufficiently high Rm. This behavior is a
consequence of the radial fine-structure of the nonaxisymmetric modes resulting
from the winding effect of differential rotation. It is this fine-structure
which presents severe resolution problems for the numerical simulation of MRI
at large Rm. For weak supercritical magnetic fields only axisymmetric modes are
unstable. Nonaxisymmetric modes need stronger fields and not too fast rotation.
If Pm is small its real value does not play any role in MRI.
[26]
oai:arXiv.org:0904.0132 [pdf] - 1001829
Eddy viscosity and turbulent Schmidt number by kink-type instability of
strong toroidal magnetic fields
Submitted: 2009-04-01
The potential of the nonaxisymmetric magnetic instability to transport
angular momentum and to mix chemicals is probed considering the stability of a
nearly uniform toroidal field between conducting cylinders with different
rotation rates. The fluid between the cylinders is assumed as incompressible
and to be of uniform density. With a linear theory the neutral-stability maps
for m=1 are computed. Rigid rotation must be subAlfvenic to allow instability
while for differential rotation with negative shear also an unstable domain
with superAlfvenic rotation exists. The rotational quenching of the magnetic
instability is strongest for magnetic Prandtl number Pm=1 and becomes much
weaker for Pm unequal 1. The effective angular momentum transport by the
instability is directed outwards(inwards) for subrotation(superrotation). The
resulting magnetic-induced eddy viscosities exceed the microscopic values by
factors of 10-100. This is only true for superAlfvenic flows; in the
strong-field limit the values remain much smaller. The same instability also
quenches concentration gradients of chemicals by its nonmagnetic fluctuations.
The corresponding diffusion coefficient remains always smaller than the
magnetic-generated eddy viscosity. A Schmidt number of order 30 is found as the
ratio of the effective viscosity and the diffusion coefficient. The magnetic
instability transports much more angular momentum than that it mixes chemicals.
[27]
oai:arXiv.org:0903.3831 [pdf] - 1001719
Tayler instability of toroidal magnetic fields in MHD Taylor-Couette
flows
Submitted: 2009-03-23
The nonaxisymmetric 'kink-type' Tayler instability (TI) of toroidal magnetic
fields is studied for conducting incompressible fluids of uniform density
between two infinitely long cylinders rotating around the same axis. It is
shown that for resting cylinders the critical Hartmann number for the unstable
modes does not depend on Pm. By rigid rotation the instability is suppressed
where the critical ratio of the rotation velocity and the Alfven velocity of
the field (only) slightly depends on the magnetic Prandtl number Pm. For Pm=1
the rotational quenching of TI takes its maximum. Rotation laws with negative
shear (i.e. d\Omega/dR<0) strongly destabilize the toroidal field if the
rotation is not too fast. For sufficiently high Reynolds numbers of rotation
the suppression of the nonaxisymmetric magnetic instability always dominates.
The angular momentum transport of the instability is anticorrelated with the
shear so that an eddy viscosity can be defined which proves to be positive. For
negative shear the Maxwell stress of the perturbations remarkably contributes
to the angular momentum transport. We have also shown the possibility of
laboratory TI experiments with a wide-gap container filled with fluid metals
like sodium or gallium. Even the effect of the rotational stabilization can be
reproduced in the laboratory with electric currents of only a few kAmp.
[28]
oai:arXiv.org:0901.0767 [pdf] - 20043
Magnetic pinch-type instability in stellar radiative zones
Submitted: 2009-01-07
The solar tachocline is shown as hydrodynamically stable against
nonaxisymmetric disturbances if it is true that no cos^{4}\theta term exists in
its rotation law. We also show that the toroidal field of 200 Gauss amplitude
which produces the tachocline in the magnetic theory of Ruediger & Kitchatinov
(1997) is stable against nonaxisymmetric MHD disturbances -- but it becomes
unstable for rotation periods slightly slower than 25 days. The instability of
such weak fields lives from the high thermal diffusivity of stellar radiation
zones compared with the magnetic diffusivity. The growth times, however, result
as very long (of order of 10\^5 rotation times). With estimations of the
chemical mixing we find the maximal possible field amplitude to be ~500 Gauss
in order to explain the observed lithium abundance of the Sun. Dynamos with
such low field amplitudes should not be relevant for the solar activity cycle.
With nonlinear simulations of MHD Taylor-Couette flows it is shown that for the
rotation-dominated magnetic instability the resulting eddy viscosity is only of
the order of the molecular viscosity. The Schmidt number as the ratio of
viscosity and chemical diffusion grows to values of ~20. For the majority of
the stellar physics applications, the magnetic-dominated Tayler instability
will be quenched by the stellar rotation.
[29]
oai:arXiv.org:0812.3790 [pdf] - 19698
Liquid metal experiments on the helical magnetorotational instability
Submitted: 2008-12-19
The magnetorotational instability (MRI) plays an essential role in the
formation of stars and black holes. By destabilizing hydrodynamically stable
Keplerian flows, the MRI triggers turbulence and enables outward transport of
angular momentum in accretion discs. We present the results of a liquid metal
Taylor-Couette experiment under the influence of helical magnetic fields that
show typical features of MRI at Reynolds numbers of the order 1000 and Hartmann
numbers of the order 10. Particular focus is laid on an improved experiment in
which split end caps are used to minimize the Ekman pumping.
[30]
oai:arXiv.org:0809.3345 [pdf] - 16515
Turbulent magnetic Prandtl numbers obtained with MHD Taylor-Couette flow
experiments
Submitted: 2008-09-19
The stability problem of MHD Taylor-Couette flows with toroidal magnetic
fields is considered in dependence on the magnetic Prandtl number. Only the
most uniform (but not current-free) field with B\_in = B\_out has been
considered. For high enough Hartmann numbers the toroidal field is always
unstable. Rigid rotation, however, stabilizes the magnetic (kink-)instability.
The axial current which drives the instability is reduced by the
electromotive force induced by the instability itself. Numerical simulations
are presented to probe this effect as a possibility to measure the turbulent
conductivity in a laboratory. It is shown numerically that in a sodium
experiment (without rotation) an eddy diffusivity 4 times the molecular
diffusivity appears resulting in a potential difference of ~34 mV/m. If the
cylinders are rotating then also the eddy viscosity can be measured. Nonlinear
simulations of the instability lead to a turbulent magnetic Prandtl number of
2.1 for a molecular magnetic Prandtl number of 0.01. The trend goes to higher
values for smaller Pm.
[31]
oai:arXiv.org:0803.1354 [pdf] - 1000664
The Tayler instability and the Hall effect in protoneutron stars
Submitted: 2008-03-10, last modified: 2008-09-10
Collapse calculations indicate that the hot newly born protoneutron stars
(PNS) rotate differentially so that strong toroidal magnetic field components
should exist in the outer crust where also the Hall effect appears to be
important when the Hall parameter \hat\beta=\omega_B \tau is of order unity.
The amplitudes of the induced toroidal magnetic fields are limited from above
by the Tayler instability. An important characteristic of the Hall effect is
its distinct dependence on the {\em sign} of the magnetic field. We find for
fast rotation that positive (negative) Hall parameters essentially reduce
(increase) the stability domain. It is thus concluded that the toroidal field
belts in PNS induced by their differential rotation should have different
amplitudes in both hemispheres which later are frozen in. Due to the effect of
magnetic suppression of the heat conductivity also the brightness of the two
hemispheres should be different. As a possible example for our scenario the
isolated neutron star RBS 1223 is considered which has been found to exhibit
different X-ray brightness at both hemispheres.
[32]
oai:arXiv.org:0808.0577 [pdf] - 15149
Stratorotational instability in MHD Taylor-Couette flows
Submitted: 2008-08-05
The stability of dissipative Taylor-Couette flows with an axial stable
density stratification and a prescribed azimuthal magnetic field is considered.
Global nonaxisymmetric solutions of the linearized MHD equations with toroidal
magnetic field, axial density stratification and differential rotation are
found for both insulating and conducting cylinder walls. Flat rotation laws
such as the quasi-Kepler law are unstable against the nonaxisymmetric
stratorotational instability (SRI). The influence of a current-free toroidal
magnetic field depends on the magnetic Prandtl number Pm: SRI is supported by
Pm > 1 and it is suppressed by Pm \lsim 1. For too flat rotation laws a smooth
transition exists to the instability which the toroidal magnetic field produces
in combination with the differential rotation. This nonaxisymmetric azimuthal
magnetorotational instability (AMRI) has been computed under the presence of an
axial density gradient. If the magnetic field between the cylinders is not
current-free then also the Tayler instability occurs and the transition from
the hydrodynamic SRI to the magnetic Tayler instability proves to be rather
complex. Most spectacular is the `ballooning' of the stability domain by the
density stratification: already a rather small rotation stabilizes magnetic
fields against the Tayler instability. An azimuthal component of the resulting
electromotive force only exists for density-stratified flows. The related
alpha-effect for magnetic SRI of Kepler rotation appears to be positive for
negative d\rho/dz <0.
[33]
oai:arXiv.org:0808.0403 [pdf] - 15116
Helical magnetorotational instability of Taylor-Couette flows in the
Rayleigh limit and for quasi-Kepler rotation
Submitted: 2008-08-04
The magnetorotational instability (MRI) of differential rotation under the
simultaneous presence of axial and azimuthal components of the (current-free)
magnetic field is considered. For rotation with uniform specific angular
momentum the MHD equations for axisymmetric perturbations are solved in a local
short-wave approximation. All the solutions are overstable for B_z \cdot B_\phi
\neq 0 with eigenfrequencies approaching the viscous frequency. For more flat
rotation laws the results of the local approximation do not comply with the
results of a global calculation of the MHD instability of Taylor-Couette flows
between rotating cylinders. -- With B_phi and B_z of the same order the
traveling-mode solutions are also prefered for flat rotation laws such as the
quasi-Kepler rotation. For magnetic Prandtl number Pm\to 0 they scale with the
Reynolds number of rotation rather than with the magnetic Reynolds number (as
for standard MRI) so that they can easily be realized in MHD laboratory
experiments. -- Regarding the nonaxisymmetric modes one finds a remarkable
influence of the ratio B_\phi /B_z only for the extrema. For B_\phi >> B_z and
for not too small Pm the nonaxisymmetric modes dominate the traveling
axisymmetric modes. For standard MRI with B_z >> B_\phi, however, the critical
Reynolds numbers of the nonaxisymmetric modes exceed the values for the
axisymmetric modes by many orders so that they are never prefered.
[34]
oai:arXiv.org:0711.0652 [pdf] - 6713
Energy distribution in nonaxisymmetric magnetic Taylor-Couette flow
Submitted: 2007-11-05
Azimuthal magnetorotational instability is a mechanism that generates
nonaxisymmetric field pattern. Nonlinear simulations in an infinite
Taylor-Couette system with current-free external field show, that not only the
linearly unstable mode m=1 appears, but also an inverse cascade transporting
energy into the axisymmetric field is possible. By varying the Reynolds number
of the flow and the Hartmann number for the magnetic field, we find that the
ratio between axisymmetric (m=0) and dominating nonaxisymmetric mode (m=1) can
be nearly free chosen. On the surface of the outer cylinder this mode
distribution appears similarly, but with weaker axisymmetric fields. We do not
find significant differences in the case that a constant current within the
flow is added.
[35]
oai:arXiv.org:0705.0883 [pdf] - 1042
Comment on "Helical MRI in magnetized Taylor-Couette flow"
Submitted: 2007-05-07, last modified: 2007-05-10
Liu et al. [Phys. Rev. E 74, 056302 (2006)] have presented a WKB analysis of
the helical magnetorotational instability (HMRI), and claim that it does not
exist for Keplerian rotation profiles. We show that if radial boundary
conditions are included, the HMRI can exist even for rotation profiles as flat
as Keplerian, provided only that at least one of the boundaries is sufficiently
conducting.
[36]
oai:arXiv.org:astro-ph/0702232 [pdf] - 89286
Sunspot models with bright rings
Submitted: 2007-02-08, last modified: 2007-02-12
A theoretical sunspot model is provided including magnetic suppression of the
diffusivities and also a strong stratification of density and temperature. Heat
diffusion alone with given magnetic field and zero mean flow only produces
(after a very long relaxation time) dark spots without any bright ring. Models
with full dynamics of both field and flow, however, provide rings and also the
observed correlation of ring temperature excess and the spot size. The rings
are formed as the result of heat transport by the resulting flow system and
increased thermal diffusivity due to reduced magnetic quenching around spots.
[37]
oai:arXiv.org:astro-ph/0701847 [pdf] - 88976
Stability of Toroidal Magnetic Fields in Rotating Stellar Radiation
Zones
Submitted: 2007-01-29
The questions of how strong magnetic fields can be stored in rotating stellar
radiative zones without being subjected to pinch-type instabilities and how
much radial mixing is produced if the fields are unstable are addressed. Linear
equations are derived for weak disturbances of magnetic and velocity fields
which are global in horizontal dimensions but short--scaled in radius. The
equations are solved to evaluate the stability of toroidal field patterns with
one or two latitudinal belts under the influence of a rigid basic rotation.
Hydrodynamic stability of latitudinal differential rotation is also considered.
The magnetic instability is essentially three--dimensional. It does not exist
in a 2D formulation with strictly horizontal disturbances on decoupled
spherical shells. Only stable (magnetically modified) r-modes are found in this
case. The instability recovers in 3D. The most rapidly growing modes for the
Sun have radial scales smaller than 1 Mm. The finite thermal conductivity makes
a strong destabilizing effect. The marginal field strength for the onset of the
instability in the upper part of the solar radiative zone is about 600 G. The
toroidal field can only slightly exceed this critical value for otherwise the
radial mixing produced by the instability would be too strong to be compatible
with the observed lithium abundance. Also the threshold for hydrodynamic
instability of differential rotation which exists in 2D is lowered in 3D. When
radial displacements are included, the value of 28% for critical shear is
reduced to 21%.
[38]
oai:arXiv.org:astro-ph/0701454 [pdf] - 88583
Destabilisation of hydrodynamically stable rotation laws by azimuthal
magnetic fields
Submitted: 2007-01-16
We consider the effect of toroidal magnetic fields on hydrodynamically stable
Taylor-Couette differential rotation flows. For current-free magnetic fields a
nonaxisymmetric m=1 magnetorotational instability arises when the magnetic
Reynolds number exceeds O (100). We then consider how this `azimuthal
magnetorotational instability' (AMRI) is modified if the magnetic field is not
current-free, but also has an associated electric current throughout the fluid.
This gives rise to current-driven Tayler instabilities (TI) that exist even
without any differential rotation at all. The interaction of the AMRI and the
TI is then considered when both electric currents and differential rotation are
present simultaneously. The magnetic Prandtl number Pm turns out to be crucial
in this case. Large Pm have a destabilizing influence, and lead to a smooth
transition between the AMRI and the TI. In contrast, small Pm have a
stabilizing influence, with a broad stable zone separating the AMRI and the TI.
In this region the differential rotation is acting to stabilize the Tayler
instabilities, with possible astrophysical applications (Ap stars). The growth
rates of both the AMRI and the TI are largely independent of Pm, with the TI
acting on the timescale of a single rotation period, and the AMRI slightly
slower, but still on the basic rotational timescale. The azimuthal drift
timescale is ~20 rotations, and may thus be a (flip-flop) timescale of stellar
activity between the rotation period and the diffusion time.
[39]
oai:arXiv.org:astro-ph/0701373 [pdf] - 88502
Structure and stability of the magnetic solar tachocline
Submitted: 2007-01-12
Rather weak fossil magnetic fields in the radiative core can produce the
solar tachocline if the field is poloidal and almost horizontal in the
tachocline region, i.e. if the field is confined within the core. This
particular field geometry is shown to result from a shallow (<~1Mm) penetration
of the meridional flow existing in the convection zone into the radiative core.
I.e., two conditions are crucial for a magnetic tachocline theory: (i) the
presence of meridional flow of a few meters per second at the base of the
convection zone, and (ii) a magnetic diffusivity inside the tachocline smaller
than 10^8 cm^2s^-1. Numerical solutions for both confined poloidal fields and
the resulting tachocline structures are presented. We find that the tachocline
thickness runs as <I>B</I>_p^-1/2 with the poloidal field amplitude falling
below 5% of the solar radius for <I>B</I>^p > 5 mG. The resulting toroidal
field amplitude inside the tachocline of about 100 G does not depend on the
<I>B</I>^p. The hydromagnetic stability of the tachocline is only briefly
discussed. For the hydrodynamic stability of latitudinal differential rotation
we found that the critical 29% of the 2D theory of Watson (1981) are reduced to
only 21% in 3D for marginal modes of about 6 Mm radial scale.
[40]
oai:arXiv.org:astro-ph/0611478 [pdf] - 86864
Theory of current instability experiments in magnetic Taylor-Couette
flows
Submitted: 2006-11-15, last modified: 2006-11-17
We consider the linear stability of dissipative MHD Taylor-Couette flow with
imposed toroidal magnetic fields. The inner and outer cylinders can be either
insulating or conducting; the inner one rotates, the outer one is stationary.
The magnetic Prandtl number can be as small as 10-5, approaching realistic
liquid-metal values. The magnetic field destabilizes the flow, except for
radial profiles of B$_\phi$(R) close to the current-free solution. The profile
with B$_{in}$=B$_{out}$ (the most uniform field) is considered in detail. For
weak fields the TC-flow is stabilized, until for moderately strong fields the
m=1 azimuthal mode dramatically destabilizes the flow again. There is thus a
maximum value for the critical Reynolds number. For sufficiently strong fields
(as measured by the Hartmann number) the toroidal field is always unstable,
even for Re=0.
The electric currents needed to generate the required toroidal fields in
laboratory experiments are a few kA if liquid sodium is used, somewhat more if
gallium is used. Weaker currents are needed for wider gaps, so a wide-gap
apparatus could succeed even with gallium. The critical Reynolds numbers are
only somewhat larger than the nonmagnetic values, so such an experiment would
require only modest rotation rates.
[41]
oai:arXiv.org:astro-ph/0607041 [pdf] - 316405
The traveling wave MRI in cylindrical Taylor-Couette flow: comparing
wavelengths and speeds in theory and experiment
Submitted: 2006-07-04, last modified: 2006-09-18
We study experimentally the flow of a liquid metal confined between
differentially rotating cylinders, in the presence of externally imposed axial
and azimuthal magnetic fields. For increasingly large azimuthal fields a
wave-like disturbance arises, traveling along the axis of the cylinders. The
wavelengths and speeds of these structures, as well as the field strengths and
rotation rates at which they arise, are broadly consistent with theoretical
predictions of such a traveling wave magnetorotational instability.
[42]
oai:arXiv.org:astro-ph/0604375 [pdf] - 81444
Do mean-field dynamos in nonrotating turbulent shear-flows exist?
Submitted: 2006-04-18
A plane-shear flow in a fluid with forced turbulence is considered. If the
fluid is electrically-conducting then a mean electromotive force (EMF) results
even without basic rotation and the magnetic diffusivity becomes a highly
anisotropic tensor. It is checked whether in this case self-excitation of a
large-scale magnetic field is possible (so-called WxJ-dynamo) and the answer is
NO. The calculations reveal the cross-stream components of the EMF
perpendicular to the mean current having the wrong signs, at least for small
magnetic Prandtl numbers. After our results numerical simulations with magnetic
Prandtl number of about unity have only a restricted meaning as the Prandtl
number dependence of the diffusivity tensor is rather strong. If, on the other
hand, the turbulence field is stratified in the vertical direction then a
dynamo active alpha-effect is produced. The critical magnetic Reynolds number
for such a self-excitation in a simple shear flow is slightly above 10 like for
the other - but much more complicated - flow patterns used in existing dynamo
experiments with liquid sodium or gallium.
[43]
oai:arXiv.org:astro-ph/0603417 [pdf] - 80626
Magnetic field confinement by meridional flow and the solar tachocline
Submitted: 2006-03-16
We show that the MHD theory that explains the solar tachocline by an effect
of the magnetic field can work with the decay modes of a fossil field in the
solar interior if the meridional flow of the convection zone penetrates
slightly the radiative zone beneath. An equatorward flow of about 10 m/s
penetrating to a maximum depth of 1000 km below the convection zone is able to
generate almost horizontal field lines in the tachocline region so that the
internal field is almost totally confined to the radiative zone. The theory of
differential solar rotation indeed provides meridional flows of about 10 m/s
and a penetration depth of < 1000 km for viscosity values that are
characteristic of a stable tachocline.
[44]
oai:arXiv.org:astro-ph/0512409 [pdf] - 78609
The angular momentum transport by the strato-rotational instability in
simulated Taylor-Couette flows
Submitted: 2005-12-15
Aims: To investigate the stability and angular momentum transport by the
strato-rotational instability in the nonlinear regime. Methods: The
hydrodynamic compressible equations are solved in a cartesian box in which the
outer cylinder is embedded. Gravity along the rotation axis leads to density
stratification. No-slip boundary conditions are used in the radial direction,
while free-slip conditions are used on the two ends of the cylinders. Results:
The strato-rotational instability is confirmed and the Reynolds stress is shown
to transport angular momentum away from the axis. However, the growth rate
decreases with increasing Reynolds number. This, as well as the presence of
boundaries renders this instability less relevant for astrophysical
applications.
[45]
oai:arXiv.org:astro-ph/0508041 [pdf] - 74930
New type of magneto-rotational instability in cylindrical Taylor-Couette
flow
Submitted: 2005-08-01
We study the stability of cylindrical Taylor-Couette flow in the presence of
combined axial and azimuthal magnetic fields, and show that adding an azimuthal
field profoundly alters the previous results for purely axial fields. For small
magnetic Prandtl numbers Pm, the critical Reynolds number Re_c for the onset of
the magneto-rotational instability becomes independent of Pm, whereas for
purely axial fields it scales as Pm^{-1}. For typical liquid metals, Re_c is
then reduced by several orders of magnitude, enough that this new design should
succeed in realizing this instability in the laboratory.
[46]
oai:arXiv.org:astro-ph/0507313 [pdf] - 74482
Taylor-Couette flow: MRI, SHI and SRI
Submitted: 2005-07-13
The linear stability theory of Taylor-Couette flows (unbounded in_z_) is
described including magnetic fields, Hall effect or a density stratification in
order to prepare laboratory experiments to probe the stability of differential
rotation in astrophysics. If an axial field is present then the
magnetorotational instability (MRI) is investigated also for small magnetic
Prandtl numbers. For rotating outer cylinder beyond the Rayleigh line
characteristic minima are found for magnetic Reynolds number of the order of 10
and for Lundquist numbers of order 1. The inclusion of extra toroidal
current-free fields leads to new oscillating solutions with rather low Reynolds
numbers and Hartmann numbers. The Hall effect establishes an unexpected
`shear-Hall instability' (SHI) where shear and magnetic field have opposite
signs. In this case even rotation laws increasing outwards may become unstable.
Recently global solutions have been found for the Taylor-Couette flows with
density stratified in axial directions (`stratorotational instability', SRI).
They exist beyond the Rayleigh line for Froude numbers of moderate order but
also not for too flat rotation laws with \hat\eta^2 < \hat\mu <\hat \eta.
[47]
oai:arXiv.org:astro-ph/0507123 [pdf] - 74292
The stability of MHD Taylor-Couette flow with current-free spiral
magnetic fields between conducting cylinders
Submitted: 2005-07-06
We study the magnetorotational instability in cylindrical Taylor-Couette
flow, with the (vertically unbounded) cylinders taken to be perfect conductors,
and with externally imposed spiral magnetic fields. The azimuthal component of
this field is generated by an axial current inside the inner cylinder, and may
be slightly stronger than the axial field. We obtain an instability beyond the
Rayleigh line, for Reynolds numbers of order 1000 and Hartmann numbers of order
10, and independent of the (small) magnetic Prandtl number. For experiments
with R_out = 2 R_in = 10 cm and Omega_out = 0.27 Omega_in, the instability
appears for liquid sodium for axial fields of ~20 Gauss and axial currents of
\~1200 A. For gallium the numbers are ~50 Gauss and ~3200 A. The vertical cell
size is about twice the cell size known for nonmagnetic experiments.
[48]
oai:arXiv.org:astro-ph/0506270 [pdf] - 73670
Three-dimensional stability of the solar tachocline
Submitted: 2005-06-13
The three-dimensional, hydrodynamic stability of the solar tachocline is
investigated based on a rotation profile as a function of both latitude and
radius. By varying the amplitude of the latitudinal differential rotation, we
find linear stability limits at various Reynolds numbers by numerical
computations. We repeated the computations with different latitudinal and
radial dependences of the angular velocity. The stability limits are all higher
than those previously found from two-dimensional approximations and higher than
the shear expected in the Sun. It is concluded that any part of the tachocline
which is radiative is hydrodynamically stable against small perturbations.
[49]
oai:arXiv.org:astro-ph/0504411 [pdf] - 72508
Differential rotation on the lower main sequence
Submitted: 2005-04-19
We compute the differential rotation of main sequence stars of the spectral
types F, G, K, and M by solving the equation of motion and the equation of
convective heat transport in a mean-field formulation. For each spectral type
the rotation rate is varied to study the dependence of the surface shear on
this parameter. The resulting rotation patterns are all solar-type. The
horizontal shear turns out to depend strongly on the effective temperature and
only weakly on the rotation rate. The meridional flow depends more strongly on
the rotation rate and has different directions in the cases of very slow and
very fast rotation, respectively.
[50]
oai:arXiv.org:astro-ph/0504173 [pdf] - 72270
Anti-solar differential rotation
Submitted: 2005-04-07
The differential rotation of anti-solar type detected by observations for
several stars may result from a fast meridional flow. The sufficiently
intensive meridional circulation may be caused by large-scale thermal
inhomogeneities or, perhaps, by tidal forcing from a companion star. First
results of simulations of the anti-solar rotation of a giant star with
magnetically induced thermal inhomogeneities are presented. Perspectives for
observational check of the theory are discussed.
[51]
oai:arXiv.org:astro-ph/0409721 [pdf] - 67802
A numerical MHD model for the solar tachocline
Submitted: 2004-09-29, last modified: 2005-03-09
There are successful approaches to explain the formation of the tachocline by
a poloidal magnetic field in the solar core. We present here the first MHD
simulations of the solar tachocline which self-consistently include the
meridional circulation. We show that the meridional flow significantly changes
the shape and the characteristics of the tachocline. We find that after the
inclusion of the meridional circulation, a tachocline can be formed even when
the poloidal field lines are crossing the boundary between the radiative zone
and the convection zone. We also discuss the effects of the magnetic Prandtl
number as well as of the magnetic Reynolds number on the properties of the
tachocline. The tachocline is much thinner for higher magnetic Reynolds numbers
and/or lower magnetic Prandtl numbers. We expect that a poloidal magnetic seed
field of around 1 G will be sufficient to produce the tachocline of the Sun.
However, the model requires the initial magnetic field to be in a narrow range
for satisfying tachocline solutions. The simulations including a stable
temperature gradient produce a shallower as well as slower meridional
circulation than the ones without it, as desired from the Lithium abundance at
the solar surface.
[52]
oai:arXiv.org:astro-ph/0407386 [pdf] - 66228
The influence of the Hall effect on the global stability of cool
protostellar disks
Submitted: 2004-07-19, last modified: 2005-02-10
The influence of the Hall effect on the global stability of cool Kepler disks
under the influence of an axial magnetic field is considered. For sufficiently
large magnetic Reynolds numbers Rm the magnetorotational instability (MRI)
exists in a finite interval of magnetic field amplitudes, $B$_{min} < $B$ <
$B$_{max}. For Kepler disks the pure MRI needs both rather high Rm
(representing the needed electrical conductivity) as well as $B$_{min} of order
0.1 G. The magnetic field pattern resulting from our global and linear
calculations is of quadrupolar parity. For magnetic fields $antiparallel$ to
the rotation axis the Hall effect reduces the minimum magnetic Reynolds number
by about one order of magnitude. The $B$_{min}, however, is even (sightly)
increased (see Fig. 6). For magnetic fields $parallel$ to the rotation axis the
Hall effect drives its own instability without the action of the Lorentz force.
The corresponding critical magnetic Reynolds number proves to be larger with
Hall effect (Rm ~ 10) than without Hall effect (Rm ~ 7) so that the Hall effect
for parallel fields even disturbs the formation of MHD-instability in cool
protoplanetary disks. If the disk is supercritical then the main result of the
Hall effect for positive fields is the strong reduction of the minimum magnetic
field amplitude which is necessary to start the instability. Observations must
show whether in star-forming regions the rotation axis and the magnetic field
orientation are correlated or are anticorrelated. If the magnetic fields are
high enough then our model predicts the dominance of fields antiparallel to the
rotation axis.
[53]
oai:arXiv.org:astro-ph/0501154 [pdf] - 70266
Stability of density-stratified viscous Taylor-Couette flows
Submitted: 2005-01-10, last modified: 2005-02-08
The stability of density-stratified viscous Taylor-Couette flows is
considered using the Boussinesq approximation but without any use of the
short-wave approximation. The flows which are unstable after the Rayleigh
criterion (\hat \mu<\hat \eta^2, with \hat \mu=\Omega_{out}/\Omega_{in} and
\hat \eta= R_{in}/R_{out}) now develop overstable axisymmetric Taylor vortices.
For the considered wide-gap container we find the nonaxisymmetric modes as the
most unstable ones. The nonaxisymmetric modes are unstable also beyond the
Rayleigh line. For such modes the instability condition seems simply to be
\hat\mu<1 as stressed by Yavneh, McWilliams & Molemaker (2001). However, we
never found unstable modes for too flat rotation laws fulfilling the condition
\hat \mu >\hat \eta. The Reynolds numbers rapidly grow to very high values if
this limit is approached (see Figs. 3 and 4).
Also striking is that the marginal stability lines for the higher $m$ do less
and less enter the region beyond the Rayleigh line so that we might have to
consider the stratorotational instability as a 'low-$m$ instability'. The
applicability of these results to the stability problem of accretion disks with
their strong stratification and fast rotation is shortly discussed.
[54]
oai:arXiv.org:astro-ph/0410729 [pdf] - 68578
Geodynamo alpha-effect derived from box simulations of rotating
magnetoconvection
Submitted: 2004-10-29
The equations for fully compressible rotating magnetoconvection are
numerically solved in a Cartesian box assuming conditions roughly suitable for
the geodynamo. The mean electromotive force describing the generation of mean
magnetic flux by convective turbulence in the rotating fluid is directly
calculated from the simulations, and the corresponding alpha-coefficients are
derived. Due to the very weak density stratification the alpha-effect changes
its sign in the middle of the box. It is positive at the top and negative at
the bottom of the convection zone. For strong magnetic fields we also find a
clear downward advection of the mean magnetic field. Both of the simulated
effects have been predicted by quasi-linear computations (Soward, 1979;
Kitchatinov and Ruediger, 1992). Finally, the possible connection of the
obtained profiles of the EMF with mean-field models of oscillating
alpha^2-dynamos is discussed.
[55]
oai:arXiv.org:astro-ph/0409246 [pdf] - 67327
Differential rotation of main sequence F stars
Submitted: 2004-09-10
The differential rotation of a 1.2 $M_\odot$ zero age MS star (spectral type
F8) is computed and the results are compared with those from a similar model of
the Sun. The rotation pattern is determined by solving the Reynolds equation
including the convective energy transport. The latter is anisotropic due to the
Coriolis force causing a horizontal temperature gradient of ~ 7 K between the
poles and the equator. Comparison of the transport mechanisms of angular
momentum (the eddy viscosity, the $\Lambda$-effect and the meridional flow)
shows that for the F star the $\Lambda$-effect is the most powerful transporter
for rotation periods of 7 d or less. In the limit of very fast rotation the
$\Lambda$-effect is balanced by the meridional flow alone and the rotation is
nearly rigid. The rotation pattern found for the F star is very similar to the
solar rotation law, but the horizontal shear is about twice the solar value. As
a function of the rotation period, the total equator-pole difference of the
angular velocity has a (slight) maximum at a period of 7 d and (slowly)
vanishes in both the limiting cases of very fast and very slow rotation. A
comparison of the solar models with those for the F-type star shows a much
stronger dependence of the differential surface rotation on the stellar
luminosity rather than on the rotation rate.
[56]
oai:arXiv.org:astro-ph/0407385 [pdf] - 66227
The eddy heat-flux in rotating turbulent convection
Submitted: 2004-07-19
The three components of the heat-flux vector $F =\rho C_p <u^\prime T^\prime>
are numerically computed for a stratified rotating turbulent convection using
the NIRVANA code in a flat box. The latitudinal component $F_\theta$ proves to
be negative (positive) in the northern (southern) hemisphere so that the heat
always flows towards the poles. As a surprise, the radial heat-flux $F_r$ peaks
at the equator rather than at the poles (Taylor numbers O(10^6)). The same
behavior is observed for the radial turbulence intensity $<u_r^{\prime 2}>$
which for \emph{free} turbulence is also believed to peak at the poles (see Eq.
(19) below). As we can show, however, the consequences of this unexpected
result (also obtained by Kaepylae, Korpi and Tuominen 2004) for the theory of
differential rotation are small as mainly the $F_\theta$ is responsible to
solve the `Taylor number puzzle'. In all our simulations the azimuthal
component $F_\phi$ proves to be negative so that the rotating turbulence
produces an westwards directed azimuthal heat-flux which should be observable.
Fluctuations with higher temperature are expected to be anticorrelated with
their own angular velocity fluctuations. We find this rotation-induced result
as understandable as the $F_\phi$ is closely related to the radial
\Lambda-effect which is known to be also negative in stratified and rapidly
rotating convection zones.
[57]
oai:arXiv.org:astro-ph/0207129 [pdf] - 50272
Differential rotation decay in the radiative envelopes of CP stars
Submitted: 2002-07-05
Stars of spectral classes A and late B are almost entirely radiative. CP
stars are a slowly rotating subgroup of these stars. It is possible that they
possessed long-lived accretion disks in their T Tauri phase. Magnetic coupling
of disk and star leads to rotational braking at the surface of the star.
Microscopic viscosities are extremely small and will not be able to reduce the
rotation rate of the core of the star. We investigate the question whether
magneto-rotational instability can provide turbulent angular momentum
transport. We illuminate the question whether or not differential rotation is
present in CP stars. Numerical MHD simulations of thick stellar shells are
performed. An initial differential rotation law is subject to the influence of
a magnetic field. The configuration gives indeed rise to magneto-rotational
instability. The emerging flows and magnetic fields transport efficiently
angular momentum outwards. Weak dependence on the magnetic Prandtl number
(~0.01 in stars) is found from the simulations. Since the estimated time-scale
of decay of differential rotation is 10^7-10^8 yr and comparable to the
life-time of A stars, we find the braking of the core to be an ongoing process
in many CP stars. The evolution of the surface rotation of CP stars with age
will be an observational challenge and of much value for verifying the
simulations.
[58]
oai:arXiv.org:astro-ph/0202299 [pdf] - 47793
Magnetoconvection and dynamo coefficients: II. Field-direction dependent
pumping of magnetic field
Submitted: 2002-02-15, last modified: 2002-06-28
We study the pumping of magnetic flux in three-dimensional compressible
magnetoconvection in the context of stellar dynamos. The simulation domain
represents a rectangular section from the lower part of a stellar convection
zone plus the underlying stably stratified layer, with a total depth of up to
five pressure scale heights. Once convection has attained a statistically
stationary state, a magnetic field is introduced. The magnetic field is
subsequently modified by the convective motions, and the resulting pumping
effects are isolated by calculating various coefficients of the expansion of
the electromotive force, uxb, in terms of components of the mean magnetic
field. The dependence of the pumping effects on rotation, latitude and other
parameters is studied. First numerical evidence is found for the existence of
pumping effects in the horizontal directions, unless the rotation axis
coincides with the vertical axis, as is the case on the poles. Evidence is
found that the pumping effects act differently on different components of the
mean magnetic field. Latitudinal pumping is mainly equatorward for toroidal
field, and can be poleward for poloidal field. Longitudinal pumping is mainly
retrograde for the radial field but prograde for the latitudinal field. The
pumping effect in the vertical direction is found to be dominated by the
diamagnetic effect, equivalent to a predominating downward advection with a
maximum speed in the turbulent case of about 10 percent of the rms convective
velocity. Where possible, an attempt is made to identify the physical origin of
the effect. Finally, some consequences of the results for stellar dynamos are
discussed.
[59]
oai:arXiv.org:astro-ph/0204308 [pdf] - 48853
Parity properties of an advection-dominated solar $\alpha^2\Om$-dynamo
Submitted: 2002-04-18
We have developed a high-precision code which solves the kinematic dynamo
problem both for given rotation law and meridional flow in the case of a low
eddy diffusivity of the order of $10^{11}$ cm$^2$/s known from the sunspot
decay. All our models work with an \alf-effect which is positive (negative) in
the northern (southern) hemisphere. It is concentrated in radial layers located
either at the top or at the bottom of the convection zone. We have also
considered an \alf-effect uniformly distributed in all the convection zone. In
the present paper the main attention is focused on i) the parity of the
solution, ii) the form of the butterfly diagram and iii) the phase relation of
the resulting field components. If the helioseismologically derived internal
solar rotation law is considered, a model without meridional flow of high
magnetic Reynolds number (corresponding to low eddy diffusivity) fails in all
the three issues in comparison with the observations. However, a meridional
flow with equatorial drift at the bottom of the convection zone of few meters
by second can indeed enforce the equatorward migration of the toroidal magnetic
field belts similar to the observed butterfly diagram but, the solution has
only a dipolar parity if the (positive) \alf-effect is located at the base of
the convection zone rather than at the top. We can, therefore, confirm the main
results of a similar study by Dikpati & Gilman (2001).
[60]
oai:arXiv.org:astro-ph/0110270 [pdf] - 45327
The stability of axisymmetric Taylor-Couette flow in hydromagnetics
Submitted: 2001-10-11
The linear marginal instability of an axisymmetric MHD Taylor-Couette flow of
infinite vertical extension is considered. The dependence of the flow stability
on magnetic Prandtl number, Pm, and gap-width between rotating cylinders is
investigated. There is an unstability for flows with resting outer cylinder
even without a magnetic field. Nevertheless, there are solutions with smaller
critical Reynolds numbers for certain (weak) magnetic fields. These solutions,
however, exist for not too small Pm or large enough gap. For hydrodynamically
stable flow, we found that a magnetic field amplitude always exsits where the
critical Reynolds number is minimal. In all calculated cases the minimal
critical Reynolds numbers are running with 1/Pm for small Pm and the critical
Reynolds numbers exceed value of 10^6 for Pm of sodium (10^{-5}) or gallium
(10^{-6}).