-
Conserving Local Magnetic Helicity in Numerical Simulations
Authors:
Yossef Zenati,
Ethan T. Vishniac
Abstract:
Magnetic helicity is robustly conserved in systems with very large magnetic Reynolds numbers, including most systems of astrophysical interest, and unlike kinetic and magnetic energy is not dissipated at small scales. This plays a major role in suppressing the kinematic large-scale dynamo and may also be responsible for driving the large-scale dynamo through the magnetic helicity flux. Numerical s…
▽ More
Magnetic helicity is robustly conserved in systems with very large magnetic Reynolds numbers, including most systems of astrophysical interest, and unlike kinetic and magnetic energy is not dissipated at small scales. This plays a major role in suppressing the kinematic large-scale dynamo and may also be responsible for driving the large-scale dynamo through the magnetic helicity flux. Numerical simulations of astrophysical systems typically lack sufficient resolution to enforce global magnetic helicity over several dynamical times. In these simulations, magnetic helicity is lost either through numerical errors or through the action of an unrealistically large resistivity. Errors in the internal distribution of magnetic helicity are equally important and typically larger. Here we propose an algorithm for enforcing strict local conservation of magnetic helicity in the Coulomb gauge in numerical simulations so that their evolution more closely approximates that of real systems.
△ Less
Submitted 8 April, 2023; v1 submitted 10 June, 2021;
originally announced June 2021.
-
3D Turbulent Reconnection: Theory, Tests and Astrophysical Implications
Authors:
Alex Lazarian,
Gregory L. Eyink,
Amir Jafari,
Grzegorz Kowal,
Hui Li,
Siyao Xu,
Ethan T. Vishniac
Abstract:
Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines co…
▽ More
Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines constantly reconnect everywhere at all scales, making magnetic reconnection an intrinsic part of turbulent cascade. We note that this is inconsistent with the usual practice of regarding magnetic lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies, starting with Lazarian & Vishniac (1999), demonstrated that 3D turbulence makes magnetic reconnection fast and that these two processes are intrinsically connected. We discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence and demonstrate that in the presence of turbulence the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. This justifies an MHD-like treatment of magnetic reconnection at scales much larger than the relevant plasma scales. We discuss numerical and observational evidences supporting the turbulent reconnection model. In particular, we show that tearing reconnection is suppressed in 3D and, unlike the 2D case, the 3D reconnection induces turbulence that makes reconnection independent of resistivity. We show that turbulent reconnection dramatically affects the key astrophysical processes, e.g., star formation, turbulent dynamo, acceleration of cosmic rays. We provide criticism of the concept of "reconnection-mediated turbulence" and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity, and why the latter has fatal conceptual flaws.
△ Less
Submitted 3 January, 2020;
originally announced January 2020.
-
Direct Detection of Black Hole-Driven Turbulence in the Centers of Galaxy Clusters
Authors:
Yuan Li,
Marie-Lou Gendron-Marsolais,
Irina Zhuravleva,
Siyao Xu,
Aurora Simionescu,
Grant R. Tremblay,
Cassandra Lochhaas,
Greg L. Bryan,
Eliot Quataert,
Norman W. Murray,
Alessandro Boselli,
Julie Hlavacek-Larrondo,
Yong Zheng,
Matteo Fossati,
Miao Li,
Eric Emsellem,
Marc Sarzi,
Lev Arzamasskiy,
Ethan T. Vishniac
Abstract:
Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide ra…
▽ More
Supermassive black holes (SMBHs) are thought to provide energy that prevents catastrophic cooling in the centers of massive galaxies and galaxy clusters. However, it remains unclear how this "feedback" process operates. We use high-resolution optical data to study the kinematics of multi-phase filamentary structures by measuring the velocity structure function (VSF) of the filaments over a wide range of scales in the centers of three nearby galaxy clusters: Perseus, Abell 2597 and Virgo. We find that the motions of the filaments are turbulent in all three clusters studied. There is a clear correlation between features of the VSFs and the sizes of bubbles inflated by SMBH driven jets. Our study demonstrates that SMBHs are the main driver of turbulent gas motions in the centers of galaxy clusters and suggests that this turbulence is an important channel for coupling feedback to the environment. Our measured amplitude of turbulence is in good agreement with Hitomi Doppler line broadening measurement and X-ray surface brightness fluctuation analysis, suggesting that the motion of the cold filaments is well-coupled to that of the hot gas. The smallest scales we probe are comparable to the mean free path in the intracluster medium (ICM). Our direct detection of turbulence on these scales provides the clearest evidence to date that isotropic viscosity is suppressed in the weakly-collisional, magnetized intracluster plasma.
△ Less
Submitted 14 November, 2019;
originally announced November 2019.
-
Kelvin-Helmholtz versus Tearing Instability: What Drives Turbulence in Stochastic Reconnection?
Authors:
Grzegorz Kowal,
Diego A. Falceta-Gonçalves,
Alex Lazarian,
Ethan T. Vishniac
Abstract:
Over the last few years it became clear that turbulent magnetic reconnection and magnetized turbulence are inseparable. It was not only shown that reconnection is responsible for violating the frozen-in condition in turbulence, but also that stochastic reconnection in 3D generates turbulence by itself. The actual mechanism responsible for this driving is still unknown. Processes such tearing mode…
▽ More
Over the last few years it became clear that turbulent magnetic reconnection and magnetized turbulence are inseparable. It was not only shown that reconnection is responsible for violating the frozen-in condition in turbulence, but also that stochastic reconnection in 3D generates turbulence by itself. The actual mechanism responsible for this driving is still unknown. Processes such tearing mode or Kelvin-Helmholtz, among other plasma instabilities, could generate turbulence from irregular current sheets. We address the nature of driving mechanism for this process and consider a relative role of tearing and Kelvin-Helmholtz instabilities for the process of turbulence generation. In particular, we analyze the conditions for development of these two instabilities within three-dimensional reconnection regions. We show that both instabilities can excite turbulence fluctuations in reconnection regions. However, tearing mode has relatively slow growth rate, and at later times it becomes partially suppressed by transverse to the current sheet component of magnetic field, generated during the growth of turbulent fluctuations. On the contrary, the Kelvin-Helmholtz instability establishes quickly in the outflow region, and at later times it dominates the turbulence generation comparing to the contribution from tearing mode. Our results demonstrate that the tearing instability is subdominant compared to the the Kelvin-Helmholtz instability in terms of generation of turbulence in the 3D reconnection layers and therefore the self-driven reconnection is turbulent reconnection with tearing instability being important only at the initial stage of the reconnection.
△ Less
Submitted 20 January, 2020; v1 submitted 19 September, 2019;
originally announced September 2019.
-
Statistics of Reconnection-Driven Turbulence
Authors:
Grzegorz Kowal,
Diego A. Falceta-Gonçalves,
Alex Lazarian,
Ethan T. Vishniac
Abstract:
Magnetic reconnection is a process that changes magnetic field topology in highly conducting fluids. Within the standard Sweet-Parker model, this process would be too slow to explain observations (e.g. solar flares). In reality, the process must be ubiquitous as astrophysical fluids are magnetized and motions of fluid elements necessarily entail crossing of magnetic frozen-in field lines and magne…
▽ More
Magnetic reconnection is a process that changes magnetic field topology in highly conducting fluids. Within the standard Sweet-Parker model, this process would be too slow to explain observations (e.g. solar flares). In reality, the process must be ubiquitous as astrophysical fluids are magnetized and motions of fluid elements necessarily entail crossing of magnetic frozen-in field lines and magnetic reconnection. In the presence of turbulence, the reconnection is independent of microscopic plasma properties, and may be much faster than previously thought, as proposed in Lazarian & Vishniac (1999) and tested in Kowal et al. (2009, 2012). However, the considered turbulence in the Lazarian-Vishniac model was imposed externally. In this work we consider reconnection-driven magnetized turbulence in realistic three-dimensional geometry initiated by stochastic noise. We demonstrate through numerical simulations that the stochastic reconnection is able to self-generate turbulence through interactions between its outflows. We analyze the statistical properties of velocity fluctuations using power spectra and anisotropy scaling, which demonstrates that the reconnection produces Kolmogorov-like turbulence, compatible with Goldreich-Sridhar (1995) model. Anisotropy statistics are, however, strongly affected by the dynamics of reconnection outflows. Once the broad turbulent region is formed, the typical anisotropy scaling $l_\parallel \propto l_\perp^{2/3}$ is formed, especially for high resolution models, were the broader range of scales is available. The decay of reconnection outflows to turbulent-like fluctuations, characterized by different anisotropy scalings, strongly depends on $β$ plasma. Moreover, the estimated reconnection rates are weakly dependent on the resolution, suggesting that no external processes are required to make reconnection fast.
△ Less
Submitted 6 February, 2017; v1 submitted 11 November, 2016;
originally announced November 2016.
-
Inertial-Range Reconnection in Magnetohydrodynamic Turbulence and in the Solar Wind
Authors:
Cristian C. Lalescu,
Yi-Kang Shi,
Gregory L. Eyink,
Theodore D. Drivas,
Ethan T. Vishniac,
Alexander Lazarian
Abstract:
In situ spacecraft data on the solar wind show events identified as magnetic reconnection with outflows and apparent "`$X$-lines" $10^{3-4}$ times ion scales. To understand the role of turbulence at these scales, we make a case study of an inertial-range reconnection event in a magnetohydrodynamic (MHD) simulation. We observe stochastic wandering of field-lines in space, breakdown of standard magn…
▽ More
In situ spacecraft data on the solar wind show events identified as magnetic reconnection with outflows and apparent "`$X$-lines" $10^{3-4}$ times ion scales. To understand the role of turbulence at these scales, we make a case study of an inertial-range reconnection event in a magnetohydrodynamic (MHD) simulation. We observe stochastic wandering of field-lines in space, breakdown of standard magnetic flux-freezing due to Richardson dispersion, and a broadened reconnection zone containing many current sheets. The coarse-grain magnetic geometry is like large-scale reconnection in the solar wind, however, with a hyperbolic flux-tube or "$X$-line" extending over integral length-scales.
△ Less
Submitted 2 March, 2015;
originally announced March 2015.
-
Turbulent Reconnection and Its Implications
Authors:
Alex Lazarian,
Gregory L. Eyink,
Ethan T. Vishniac,
Grzegorz Kowal
Abstract:
Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes in magnetized plasmas. In most astrophysical environments the Reynolds numbers are large and therefore the transition to turbulence is inevitable. This turbulence must be taken into account for any theory of magnetic reconnection, since the initially laminar configurations can trans…
▽ More
Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes in magnetized plasmas. In most astrophysical environments the Reynolds numbers are large and therefore the transition to turbulence is inevitable. This turbulence must be taken into account for any theory of magnetic reconnection, since the initially laminar configurations can transit to the turbulence state, what is demonstrated by 3D high resolution numerical simulations. We discuss ideas of how turbulence can modify reconnection with the focus on the Lazarian & Vishniac (1999) reconnection model and present numerical evidence supporting the model and demonstrate that it is closely connected to the concept of Richardson diffusion and compatible with the Lagrangian dynamics of magnetized fluids. We point out that the Generalized Ohm's Law, that accounts for turbulent motion, predicts the subdominance of the microphysical plasma effects for a realistically turbulent media. We show that on of the most dramatic consequences of turbulence is the violation of the generally accepted notion of magnetic flux freezing. This notion is a corner stone of most theories dealing with magnetized plasmas and therefore its change induces fundamental shifts in accepted paradigms like turbulent reconnection entailing the diffusion process that is essential for understanding star formation. We argue, that at sufficiently high Reynolds numbers the process of tearing reconnection should transfer to turbulent reconnection. We discuss flares predicted by turbulent reconnection and relate them to solar flares and gamma ray bursts. We analyze solar observations, measurements in the solar wind or heliospheric current sheet, and show their correspondence with turbulent reconnection predictions. Finally, we discuss 1st Order Fermi acceleration as a natural consequence of the turbulent reconnection.
△ Less
Submitted 4 February, 2015;
originally announced February 2015.
-
Reconnection Studies Under Different Types of Turbulence Driving
Authors:
Grzegorz Kowal,
Alex Lazarian,
Ethan T. Vishniac,
Katarzyna Otmianowska-Mazur
Abstract:
We study a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian and Vishniac (1999) using three-dimensional direct numerical simulations. The model has been already successfully tested in Kowal et al. (2009) confirming the dependencies of the reconnection speed $V_{rec}$ on the turbulence injection power $P_{inj}$ and the injection scale $l_{inj}$ expressed b…
▽ More
We study a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian and Vishniac (1999) using three-dimensional direct numerical simulations. The model has been already successfully tested in Kowal et al. (2009) confirming the dependencies of the reconnection speed $V_{rec}$ on the turbulence injection power $P_{inj}$ and the injection scale $l_{inj}$ expressed by a constraint $V_{rec} \sim P_{inj}^{1/2} l_{inj}^{3/4}$ and no observed dependency on Ohmic resistivity. In Kowal et al. (2009), in order to drive turbulence, we injected velocity fluctuations in Fourier space with frequencies concentrated around $k_{inj}=1/l_{inj}$, as described in Alvelius (1999). In this paper we extend our previous studies by comparing fast magnetic reconnection under different mechanisms of turbulence injection by introducing a new way of turbulence driving. The new method injects velocity or magnetic eddies with a specified amplitude and scale in random locations directly in real space. We provide exact relations between the eddy parameters and turbulent power and injection scale. We performed simulations with new forcing in order to study turbulent power and injection scale dependencies. The results show no discrepancy between models with two different methods of turbulence driving exposing the same scalings in both cases. This is in agreement with the Lazarian and Vishniac (1999) predictions. In addition, we performed a series of models with varying viscosity $ν$. Although Lazarian and Vishniac (1999) do not provide any prediction for this dependence, we report a weak relation between the reconnection speed with viscosity, $V_{rec}\simν^{-1/4}$.
△ Less
Submitted 13 March, 2012;
originally announced March 2012.
-
Fast Magnetic Reconnection and Spontaneous Stochasticity
Authors:
Gregory L. Eyink,
Alex Lazarian,
Ethan T. Vishniac
Abstract:
Magnetic field-lines in astrophysical plasmas are expected to be frozen-in at scales larger than the ion gyroradius. The rapid reconnection of magnetic flux structures with dimensions vastly larger than the gyroradius requires a breakdown in the standard Alfvén flux-freezing law. We attribute this breakdown to ubiquitous MHD plasma turbulence with power-law scaling ranges of velocity and magnetic…
▽ More
Magnetic field-lines in astrophysical plasmas are expected to be frozen-in at scales larger than the ion gyroradius. The rapid reconnection of magnetic flux structures with dimensions vastly larger than the gyroradius requires a breakdown in the standard Alfvén flux-freezing law. We attribute this breakdown to ubiquitous MHD plasma turbulence with power-law scaling ranges of velocity and magnetic energy spectra. Lagrangian particle trajectories in such environments become "spontaneously stochastic", so that infinitely-many magnetic field-lines are advected to each point and must be averaged to obtain the resultant magnetic field. The relative distance between initial magnetic field lines which arrive to the same final point depends upon the properties of two-particle turbulent dispersion. We develop predictions based on the phenomenological Goldreich & Sridhar theory of strong MHD turbulence and on weak MHD turbulence theory. We recover the predictions of the Lazarian & Vishniac theory for the reconnection rate of large-scale magnetic structures. Lazarian & Vishniac also invoked "spontaneous stochasticity", but of the field-lines rather than of the Lagrangian trajectories. More recent theories of fast magnetic reconnection appeal to microscopic plasma processes that lead to additional terms in the generalized Ohm's law, such as the collisionless Hall term. We estimate quantitatively the effect of such processes on the inertial-range turbulence dynamics and find them to be negligible in most astrophysical environments. For example, the predictions of the Lazarian-Vishniac theory are unchanged in Hall MHD turbulence with an extended inertial range, whenever the ion skin depth $δ_i$ is much smaller than the turbulent integral length or injection-scale $L_i.$
△ Less
Submitted 17 August, 2011; v1 submitted 9 March, 2011;
originally announced March 2011.
-
Fast magnetic reconnection in three dimensional MHD simulations
Authors:
Bijia Pang,
Ue-Li Pen,
Ethan T. Vishniac
Abstract:
We present a constructive numerical example of fast magnetic reconnection in a three dimensional periodic box. Reconnection is initiated by a strong, localized perturbation to the field lines. The solution is intrinsically three dimensional, and its gross properties do not depend on the details of the simulations. $\sim 50%$ of the magnetic energy is released in an event which lasts about one Al…
▽ More
We present a constructive numerical example of fast magnetic reconnection in a three dimensional periodic box. Reconnection is initiated by a strong, localized perturbation to the field lines. The solution is intrinsically three dimensional, and its gross properties do not depend on the details of the simulations. $\sim 50%$ of the magnetic energy is released in an event which lasts about one Alfven time, but only after a delay during which the field lines evolve into a critical configuration. We present a physical picture of the process. The reconnection regions are dynamical and mutually interacting. In the comoving frame of these regions, reconnection occurs through an X-like point, analogous to Petschek reconnection. The dynamics appear to be driven by global flows, not local processes.
△ Less
Submitted 27 February, 2010;
originally announced March 2010.
-
Numerical Tests of Fast Reconnection in Weakly Stochastic Magnetic Fields
Authors:
G. Kowal,
A. Lazarian,
E. T. Vishniac,
K. Otmianowska-Mazur
Abstract:
We study the effects of turbulence on magnetic reconnection using 3D numerical simulations. This is the first attempt to test a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian & Vishniac (1999). This model predicts that weak turbulence, generically present in most of astrophysical systems, enhances the rate of reconnection by reducing the transverse sc…
▽ More
We study the effects of turbulence on magnetic reconnection using 3D numerical simulations. This is the first attempt to test a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian & Vishniac (1999). This model predicts that weak turbulence, generically present in most of astrophysical systems, enhances the rate of reconnection by reducing the transverse scale for reconnection events and by allowing many independent flux reconnection events to occur simultaneously. As a result the reconnection speed becomes independent of Ohmic resistivity and is determined by the magnetic field wandering induced by turbulence. To quantify the reconnection speed we use both an intuitive definition, i.e. the speed of the reconnected flux inflow, as well as a more sophisticated definition based on a formally derived analytical expression. Our results confirm the predictions of the Lazarian & Vishniac model. In particular, we find that Vrec Pinj^(1/2), as predicted by the model. The dependence on the injection scale for some of our models is a bit weaker than expected, i.e. l^(3/4), compared to the predicted linear dependence on the injection scale, which may require some refinement of the model or may be due to the effects like finite size of the excitation region. The reconnection speed was found to depend on the expected rate of magnetic field wandering and not on the magnitude of the guide field. In our models, we see no dependence on the guide field when its strength is comparable to the reconnected component. More importantly, while in the absence of turbulence we successfully reproduce the Sweet-Parker scaling of reconnection, in the presence of turbulence we do not observe any dependence on Ohmic resistivity, confirming that our reconnection is fast.
△ Less
Submitted 30 April, 2009; v1 submitted 11 March, 2009;
originally announced March 2009.
-
The Saturation Limit of the Magnetorotational Instability
Authors:
Ethan T. Vishniac
Abstract:
Simulations of the magnetorotational instability (MRI) in a homogeneous shearing box have shown that the asymptotic strength of the magnetic field declines steeply with increasing resolution. Here I model the MRI driven dynamo as a large scale dynamo driven by the vertical magnetic helicity flux. This growth is balanced by large scale mixing driven by a secondary instability. The saturated magne…
▽ More
Simulations of the magnetorotational instability (MRI) in a homogeneous shearing box have shown that the asymptotic strength of the magnetic field declines steeply with increasing resolution. Here I model the MRI driven dynamo as a large scale dynamo driven by the vertical magnetic helicity flux. This growth is balanced by large scale mixing driven by a secondary instability. The saturated magnetic energy density depends almost linearly on the vertical height of the typical eddies. The MRI can drive eddies with arbitrarily large vertical wavenumber, so the eddy thickness is either set by diffusive effects, by the magnetic tension of a large scale vertical field component, or by magnetic buoyancy effects. In homogeneous, zero magnetic flux, simulations only the first effect applies and the saturated limit of the dynamo is determined by explicit or numerical diffusion. The exact result depends on the numerical details, but is consistent with previous work, including the claim that the saturated field energy scales as the gas pressure to the one quarter power (which we interpret as an artifact of numerical dissipation). The magnetic energy density in a homogeneous shearing box will tend to zero as the resolution of the simulation increases, but this has no consequences for the dynamo or for angular momentum transport in real accretion disks. The claim that the saturated state depends on the magnetic Prandtl number may also be an artifact of simulations in which microphysical transport coefficients set the MRI eddy thickness. Finally, the efficiency of the MRI dynamo is a function of the ratio of the Alfvén velocity to the product of the pressure scale height and the local shear. As this approaches unity from below the dynamo reaches maximum efficiency.
△ Less
Submitted 5 February, 2009;
originally announced February 2009.
-
Numerical Studies of Weakly Stochastic Magnetic Reconnection
Authors:
G. Kowal,
A. Lazarian,
E. T. Vishniac,
K. Otmianowska-Mazur
Abstract:
We study the effects of turbulence on magnetic reconnection using three-dimensional numerical simulations. This is the first attempt to test a model of fast magnetic reconnection proposed by Lazarian & Vishniac (1999), which assumes the presence of weak, small-scale magnetic field structure near the current sheet. This affects the rate of reconnection by reducing the transverse scale for reconne…
▽ More
We study the effects of turbulence on magnetic reconnection using three-dimensional numerical simulations. This is the first attempt to test a model of fast magnetic reconnection proposed by Lazarian & Vishniac (1999), which assumes the presence of weak, small-scale magnetic field structure near the current sheet. This affects the rate of reconnection by reducing the transverse scale for reconnection flows and by allowing many independent flux reconnection events to occur simultaneously. We performed a number of simulations to test the dependencies of the reconnection speed, defined as the ratio of the inflow velocity to the Alfven speed, on the turbulence power, the injection scale and resistivity. Our results show that turbulence significantly affects the topology of magnetic field near the diffusion region and increases the thickness of the outflow region. We confirm the predictions of the Lazarian & Vishniac model. In particular, we report the growth of the reconnection speed proportional to ~ V^2, where V is the amplitude of velocity at the injection scale. It depends on the injection scale l as ~ (l/L)^(2/3), where L is the size of the system, which is somewhat faster but still roughly consistent with the theoretical expectations. We also show that for 3D reconnection the Ohmic resistivity is important in the local reconnection events only, and the global reconnection rate in the presence of turbulence does not depend on it.
△ Less
Submitted 11 December, 2008;
originally announced December 2008.
-
Glimm's Method for Relativistic Hydrodynamics
Authors:
J. K. Cannizzo,
N. Gehrels,
E. T. Vishniac
Abstract:
We present the results of standard one-dimensional test problems in relativistic hydrodynamics using Glimm's (random choice) method, and compare them to results obtained using finite differencing methods. For problems containing profiles with sharp edges, such as shocks, we find Glimm's method yields global errors ~1-3 orders of magnitude smaller than the traditional techniques. The strongest di…
▽ More
We present the results of standard one-dimensional test problems in relativistic hydrodynamics using Glimm's (random choice) method, and compare them to results obtained using finite differencing methods. For problems containing profiles with sharp edges, such as shocks, we find Glimm's method yields global errors ~1-3 orders of magnitude smaller than the traditional techniques. The strongest differences are seen for problems in which a shear field is superposed. For smooth flows, Glimm's method is inferior to standard methods. The location of specific features can be off by up to two grid points with respect to an exact solution in Glimm's method, and furthermore curved states are not modeled optimally since the method idealizes solutions as being composed of piecewise constant states. Thus although Glimm's method is superior at correctly resolving sharp features, especially in the presence of shear, for realistic applications in which one typically finds smooth flows plus strong gradients or discontinuities, standard FD methods yield smaller global errors. Glimm's method may prove useful in certain applications such as GRB afterglow shock propagation into a uniform medium.
△ Less
Submitted 8 February, 2008;
originally announced February 2008.
-
Magnetic Field Structure and Stochastic Reconnection in a Partially Ionized Gas
Authors:
A. Lazarian,
E. T. Vishniac,
J. Cho
Abstract:
We consider stochastic reconnection in a magnetized, partially ionized medium. Stochastic reconnection is a generic effect, due to field line wandering, in which the speed of reconnection is determined by the ability of ejected plasma to diffuse away from the current sheet along magnetic field lines, rather than by the details of current sheet structure. We consider the limit of weak stochastici…
▽ More
We consider stochastic reconnection in a magnetized, partially ionized medium. Stochastic reconnection is a generic effect, due to field line wandering, in which the speed of reconnection is determined by the ability of ejected plasma to diffuse away from the current sheet along magnetic field lines, rather than by the details of current sheet structure. We consider the limit of weak stochasticity, so that the mean magnetic field energy density is greater than either the turbulent kinetic energy density or the energy density associated with the fluctuating component of the field. We consider field line stochasticity generated through a turbulent cascade, which leads us to consider the effect of neutral drag on the turbulent cascade of energy. In a collisionless plasma, neutral particle viscosity and ion-neutral drag will damp mid-scale turbulent motions, but the power spectrum of the magnetic perturbations extends below the viscous cutoff scale. We give a simple physical picture of the magnetic field structure below this cutoff, consistent with numerical experiments. We provide arguments for the reemergence of the turbulent cascade well below the viscous cut-off scale and derive estimates for field line diffusion on all scales. We note that this explains the persistence of a single power law form for the turbulent power spectrum of the interstellar medium, from scales of tens of parsecs down to thousands of kilometers. We find that under typical conditions in the ISM stochastic reconnection speeds are reduced by the presence of neutrals, but by no more than an order of magnitude.
△ Less
Submitted 11 November, 2003;
originally announced November 2003.
-
A Numerical Gamma-Ray Burst Simulation Using Three-Dimensional Relativistic Hydrodynamics: The Transition from Spherical to Jet-like Expansion
Authors:
J. K. Cannizzo,
N. Gehrels,
E. T. Vishniac
Abstract:
We present the first unrestricted, three-dimensional relativistic hydrodynamical calculations of the blob of gas associated with the jet producing a gamma-ray burst. We investigate the deceleration phase of the blob corresponding to the time when afterglow radiation is produced, concentrating on the transition in which the relativistic beaming 1/gamma goes from being less than theta, where gamma…
▽ More
We present the first unrestricted, three-dimensional relativistic hydrodynamical calculations of the blob of gas associated with the jet producing a gamma-ray burst. We investigate the deceleration phase of the blob corresponding to the time when afterglow radiation is produced, concentrating on the transition in which the relativistic beaming 1/gamma goes from being less than theta, where gamma is the bulk Lorentz factor and theta is the angular width of the jet, to 1/gamma > theta. We study the time dependent evolution of the physical parameters associated with the jet, both parallel to the direction of motion and perpendicular to it. We calculate light curves for observers at varying angles with respect to the velocity vector of the blob, assuming optically thin emission that scales with the local pressure. Our main findings are that (i) gas ahead of the advancing blob does not accrete onto and merge with the blob material but rather flows around the blob, (ii) the decay light curve steepens at a time corresponding roughly to 1/gamma equals theta (in accord with earlier studies), and (iii) the rate of decrease of the forward component of momentum in the blob is well-fit by a simple model in which the gas in front of the blob exerts a drag force on the blob, and the cross sectional area of the blob increases quadratically with laboratory time (or distance).
△ Less
Submitted 4 October, 2003;
originally announced October 2003.
-
Ordinary and Viscosity-Damped MHD Turbulence
Authors:
Jungyeon Cho,
A. Lazarian,
Ethan T. Vishniac
Abstract:
We compare the properties of ordinary strong magnetohydrodynamic (MHD) turbulence in a strongly magnetized medium with the recently discovered viscosity-damped regime. We focus on energy spectra, anisotropy, and intermittency. Our most surprising conclusion is that in ordinary strong MHD turbulence the velocity and magnetic fields show different high-order structure function scalings. Moreover t…
▽ More
We compare the properties of ordinary strong magnetohydrodynamic (MHD) turbulence in a strongly magnetized medium with the recently discovered viscosity-damped regime. We focus on energy spectra, anisotropy, and intermittency. Our most surprising conclusion is that in ordinary strong MHD turbulence the velocity and magnetic fields show different high-order structure function scalings. Moreover this scaling depends on whether the intermittency is viewed in a global or local system of reference. This reconciles seemingly contradictory earlier results. On the other hand, the intermittency scaling for viscosity-damped turbulence is very different, and difficult to understand in terms of the usual phenomenological models for intermittency in turbulence. Our remaining results are in reasonable agreement with expectations. First, we find that our high resolution simulations for ordinary MHD turbulence show that the energy spectra are {\it compatible} with a Kolmogorov spectrum, while viscosity-damped turbulence shows a shallow $k^{-1}$ spectrum for the magnetic fluctuations. Second, a new numerical technique confirms that ordinary MHD turbulence exhibits Goldreich-Sridhar type anisotropy, while viscosity-damped MHD turbulence shows extremely anisotropic eddy structures. Finally, we show that many properties of incompressible turbulence for both the ordinary and viscosity-damped regimes carry over to the case of compressible turbulence.
△ Less
Submitted 3 July, 2003; v1 submitted 12 May, 2003;
originally announced May 2003.
-
Problems and Progress in Astrophysical Dynamos
Authors:
E. T. Vishniac,
A. Lazarian,
J. Cho
Abstract:
Astrophysical objects with negligible resistivity are often threaded by large scale magnetic fields. The generation of these fields is somewhat mysterious, since a magnetic field in a perfectly conducting fluid cannot change the flux threading a fluid element, or the field topology. Classical dynamo theory evades this limit by assuming that magnetic reconnection is fast, even for vanishing resis…
▽ More
Astrophysical objects with negligible resistivity are often threaded by large scale magnetic fields. The generation of these fields is somewhat mysterious, since a magnetic field in a perfectly conducting fluid cannot change the flux threading a fluid element, or the field topology. Classical dynamo theory evades this limit by assuming that magnetic reconnection is fast, even for vanishing resistivity, and that the large scale field can be generated by the action of kinetic helicity. Both these claims have been severely criticized, and the latter appears to conflict with strong theoretical arguments based on magnetic helicity conservation and a series of numerical simulations. Here we discuss recent efforts to explain fast magnetic reconnection through the topological effects of a weak stochastic magnetic field component. We also show how mean-field dynamo theory can be recast in a form which respects magnetic helicity conservation, and how this changes our understanding of astrophysical dynamos. Finally, we comment briefly on why an asymmetry between small scale magnetic and velocity fields is necessary for dynamo action, and how it can arise naturally.
△ Less
Submitted 31 May, 2002;
originally announced May 2002.
-
Magnetic Helicity Conservation and Astrophysical Dynamos
Authors:
Ethan T. Vishniac,
Jungyeon Cho
Abstract:
We construct a magnetic helicity conserving dynamo theory which incorporates a calculated magnetic helicity current. In this model the fluid helicity plays a small role in large scale magnetic field generation. Instead, the dynamo process is dominated by a new quantity, derived from asymmetries in the second derivative of the velocity correlation function, closely related to the `twist and fold'…
▽ More
We construct a magnetic helicity conserving dynamo theory which incorporates a calculated magnetic helicity current. In this model the fluid helicity plays a small role in large scale magnetic field generation. Instead, the dynamo process is dominated by a new quantity, derived from asymmetries in the second derivative of the velocity correlation function, closely related to the `twist and fold' dynamo model. The turbulent damping term is, as expected, almost unchanged. Numerical simulations with a spatially constant fluid helicity and vanishing resistivity are not expected to generate large scale fields in equipartition with the turbulent energy density. The prospects for driving a fast dynamo under these circumstances are uncertain, but if it is possible, then the field must be largely force-free. On the other hand, there is an efficient analog to the $α-Ω$ dynamo. Systems whose turbulence is driven by some anisotropic local instability in a shearing flow, like real stars and accretion disks, and some computer simulations, may successfully drive the generation of strong large scale magnetic fields, provided that $\partial_rΩ< \partial_θv_zω_θ>>0$. We show that this criterion is usually satisfied. Such dynamos will include a persistent, spatially coherent vertical magnetic helicity current with the same sign as $-\partial_rΩ$, that is, positive for an accretion disk and negative for the Sun. We comment on the role of random magnetic helicity currents in storing turbulent energy in a disordered magnetic field, which will generate an equipartition, disordered field in a turbulent medium, and also a declining long wavelength tail to the power spectrum. As a result, calculations of the galactic `seed' field are largely irrelevant.
△ Less
Submitted 18 December, 2000; v1 submitted 18 October, 2000;
originally announced October 2000.
-
The Generation of Magnetic Fields Through Driven Turbulence
Authors:
Jungyeon Cho,
Ethan T. Vishniac
Abstract:
We have tested the ability of driven turbulence to generate magnetic field structure from a weak uniform field using three dimensional numerical simulations of incompressible turbulence. We used a pseudo-spectral code with a numerical resolution of up to $144^3$ collocation points. We find that the magnetic fields are amplified through field line stretching at a rate proportional to the differen…
▽ More
We have tested the ability of driven turbulence to generate magnetic field structure from a weak uniform field using three dimensional numerical simulations of incompressible turbulence. We used a pseudo-spectral code with a numerical resolution of up to $144^3$ collocation points. We find that the magnetic fields are amplified through field line stretching at a rate proportional to the difference between the velocity and the magnetic field strength times a constant. Equipartition between the kinetic and magnetic energy densities occurs at a scale somewhat smaller than the kinetic energy peak. Above the equipartition scale the velocity structure is, as expected, nearly isotropic. The magnetic field structure at these scales is uncertain, but the field correlation function is very weak. At the equipartition scale the magnetic fields show only a moderate degree of anisotropy, so that the typical radius of curvature of field lines is comparable to the typical perpendicular scale for field reversal. In other words, there are few field reversals within eddies at the equipartition scale, and no fine-grained series of reversals at smaller scales. At scales below the equipartition scale, both velocity and magnetic structures are anisotropic; the eddies are stretched along the local magnetic field lines, and the magnetic energy dominates the kinetic energy on the same scale by a factor which increases at higher wavenumbers. We do not show a scale-free inertial range, but the power spectra are a function of resolution and/or the imposed viscosity and resistivity. Our results are consistent with the emergence of a scale-free inertial range at higher Reynolds numbers.
△ Less
Submitted 14 June, 2000; v1 submitted 27 March, 2000;
originally announced March 2000.
-
The Anisotropy of MHD Alfvénic Turbulence
Authors:
Jungyeon Cho,
Ethan T. Vishniac
Abstract:
We perform direct 3-dimensional numerical simulations for magnetohydrodynamic (MHD) turbulence in a periodic box of size $2π$ threaded by strong uniform magnetic fields. We use a pseudo-spectral code with hyperviscosity and hyperdiffusivity to solve the incompressible MHD equations. We analyze the structure of the eddies as a function of scale. A straightforward calculation of anisotropy in wave…
▽ More
We perform direct 3-dimensional numerical simulations for magnetohydrodynamic (MHD) turbulence in a periodic box of size $2π$ threaded by strong uniform magnetic fields. We use a pseudo-spectral code with hyperviscosity and hyperdiffusivity to solve the incompressible MHD equations. We analyze the structure of the eddies as a function of scale. A straightforward calculation of anisotropy in wavevector space shows that the anisotropy is scale-{\it independent}. We discuss why this is {\it not} the true scaling law and how the curvature of large-scale magnetic fields affects the power spectrum and leads to the wrong conclusion. When we correct for this effect, we find that the anisotropy of eddies depends on their size: smaller eddies are more elongated than larger ones along {\it local} magnetic field lines. The results are consistent with the scaling law $\tilde{k}_{\parallel} \sim \tilde{k}_{\perp}^{2/3}$ proposed by Goldreich and Sridhar (1995, 1997). Here $\tilde{k}_{\|}$ (and $\tilde{k}_{\perp}$) are wavenumbers measured relative to the local magnetic field direction. However, we see some systematic deviations which may be a sign of limitations to the model, or our inability to fully resolve the inertial range of turbulence in our simulations.
△ Less
Submitted 14 June, 2000; v1 submitted 27 March, 2000;
originally announced March 2000.
-
Thermal Equilibrium Curves and Turbulent Mixing in Keplerian Accretion Disks
Authors:
Pin-Gao Gu,
Ethan T. Vishniac,
John K. Cannizzo
Abstract:
We consider vertical heat transport in Keplerian accretion disks, including the effects of radiation, convection, and turbulent mixing driven by the Balbus-Hawley instability, in astronomical systems ranging from dwarf novae (DNe), and soft X-ray transients (SXTs), to active galactic nuclei (AGN). We propose a modified, anisotropic form of mixing-length theory, which includes radiative and turbu…
▽ More
We consider vertical heat transport in Keplerian accretion disks, including the effects of radiation, convection, and turbulent mixing driven by the Balbus-Hawley instability, in astronomical systems ranging from dwarf novae (DNe), and soft X-ray transients (SXTs), to active galactic nuclei (AGN). We propose a modified, anisotropic form of mixing-length theory, which includes radiative and turbulent damping. We also include turbulent heat transport, which acts everywhere within disks, regardless of whether or not they are stably stratified, and can move entropy in either direction. We have generated a series of vertical structure models and thermal equilibrium curves using the scaling law for the viscosity parameter $α$ suggested by the exponential decay of the X-ray luminosity in SXTs. We have also included equilibrium curves for DNe using an $α$ which is constant down to a small magnetic Reynolds number ($\sim 10^4$). Our models indicate that weak convection is usually eliminated by turbulent radial mixing. The substitution of turbulent heat transport for convection is more important on the unstable branches of thermal equilibrium S-curves when $α$ is larger. The low temperature turnover points $Σ_{max}$ on the equilibrium S-curves are significantly reduced by turbulent mixing in DNe and SXT disks. However, in AGN disks the standard mixing-length theory for convection is still a useful approximation when we use the scaling law for $α$, since these disks are very thin at the relevant radii. In accordance with previous work, we find that constant $α$ models give almost vertical S-curves in the $Σ-T$ plane and consequently imply very slow, possibly oscillating, cooling waves.
△ Less
Submitted 16 August, 1999;
originally announced August 1999.
-
Spectral Transition and Torque Reversal in X-ray Pulsar 4U 1626-67
Authors:
Insu Yi,
Ethan T. Vishniac
Abstract:
The accretion-powered, X-ray pulsar 4U 1626-67 has recently shown an abrupt torque reversal accompanied by a dramatic spectral transition and a relatively small luminosity change. The time-averaged X-ray spectrum during spin-down is considerably harder than during spin-up. The observed torque reversal can be explained by an accretion flow transition triggered by a gradual change in the mass accr…
▽ More
The accretion-powered, X-ray pulsar 4U 1626-67 has recently shown an abrupt torque reversal accompanied by a dramatic spectral transition and a relatively small luminosity change. The time-averaged X-ray spectrum during spin-down is considerably harder than during spin-up. The observed torque reversal can be explained by an accretion flow transition triggered by a gradual change in the mass accretion rate. The sudden transition to spin-down is caused by a change in the accretion flow rotation from Keplerian to sub-Keplerian. 4U 1626-67 is estimated to be near spin equilibrium with a mass accretion rate Mdot~2x10**16 g/s, Mdot decreasing at a rate ~6x10**14 g/s/yr, and a polar surface magnetic field of ~2b_p**{-1/2} 10^**12G where b_p is the magnetic pitch. During spin-up, the Keplerian flow remains geometrically thin and cool. During spin-down, the sub-Keplerian flow becomes geometrically thick and hot. Soft photons from near the stellar surface are Compton up-scattered by the hot accretion flow during spin-down while during spin-up such scattering is unlikely due to the small scale-height and low temperature of the flow. This mechanism accounts for the observed spectral hardening and small luminosity change. The scattering occurs in a hot radially falling column of material with a scattering depth ~0.3 and a temperature ~10^9K. The X-ray luminosity at energies >5keV could be a poor indicator of the mass accretion rate. We briefly discuss the possible application of this mechanism to GX 1+4, although there are indications that this system is significantly different from other torque-reversal systems.
△ Less
Submitted 17 May, 1999;
originally announced May 1999.
-
Galactic and Accretion Disk Dynamos
Authors:
Ethan T. Vishniac
Abstract:
Dynamos in astrophysical disks are usually explained in terms of the standard alpha-omega mean field dynamo model where the local helicity generates a radial field component from an azimuthal field. The subsequent shearing of the radial field gives rise to exponentially growing dynamo modes. There are several problems with this model. The exponentiation time for the galactic dynamo is hard to ca…
▽ More
Dynamos in astrophysical disks are usually explained in terms of the standard alpha-omega mean field dynamo model where the local helicity generates a radial field component from an azimuthal field. The subsequent shearing of the radial field gives rise to exponentially growing dynamo modes. There are several problems with this model. The exponentiation time for the galactic dynamo is hard to calculate, but is probably uncomfortably long. Moreover, numerical simulations of magnetic fields in shearing flows indicate that the presence of a dynamo does not depend on a non-zero average helicity. However, these difficulties can be overcome by including a fluctuating helicity driven by hydrodynamic or magnetic instabilities. Unlike traditional disk dynamo models, this `incoherent' dynamo does not depend on the presence of systematic fluid helicity or any kind of vertical symmetry breaking. It will depend on geometry, in the sense that the dynamo growth rate becomes smaller for very thin disks, in agreement with constraints taken from the study of X-ray novae. In this picture the galactic dynamo will operate efficiently, but the resulting field will have a radial coherence length which is a fraction of the galactic radius.
△ Less
Submitted 11 December, 1998;
originally announced December 1998.
-
Reconnection in a Weakly Stochastic Field
Authors:
A. Lazarian,
E. T. Vishniac
Abstract:
We examine the effect of weak, small scale magnetic field structure on the rate of reconnection in a strongly magnetized plasma. This affects the rate of reconnection by reducing the transverse scale for reconnection flows, and by allowing many independent flux reconnection events to occur simultaneously. Allowing only for the first effect and using Goldreich and Sridhar's model of strong turbul…
▽ More
We examine the effect of weak, small scale magnetic field structure on the rate of reconnection in a strongly magnetized plasma. This affects the rate of reconnection by reducing the transverse scale for reconnection flows, and by allowing many independent flux reconnection events to occur simultaneously. Allowing only for the first effect and using Goldreich and Sridhar's model of strong turbulence in a magnetized plasma with negligible intermittency, we find that the lower limit for the reconnection speed is the Alfven speed times the Lundquist number to the power (-3/16). The upper limit on the reconnection speed is typically a large fraction of Alfven speed. We argue that generic reconnection in turbulent plasmas will normally occur at close to this upper limit. The fraction of magnetic energy that goes directly into electron heating scales as Lundquist number to the power (-2/5) and the thickness of the current sheet scales as the Lundquist number to the power (-3/5). A significant fraction of the magnetic energy goes into high frequency Alfven waves. We claim that the qualitative sense of these conclusions, that reconnection is fast even though current sheets are narrow, is almost independent of the local physics of reconnection and the nature of the turbulent cascade. As the consequence of this the Galactic and Solar dynamos are generically fast, i.e. do not depend on the plasma resistivity.
△ Less
Submitted 19 January, 1999; v1 submitted 2 November, 1998;
originally announced November 1998.
-
Reconnection via the Tearing Instability
Authors:
A. Lazarian,
Ethan T. Vishniac
Abstract:
We discuss the role of tearing instabilities in magnetic reconnection. In three dimensions this instability leads to the formation of strong Alfvenic waves that remove plasma efficiently from the reconnection layer. As a result the instability proceeds at high rates while staying close to the linear regime. Our calculations show that for a resistive fluid the reconnection speed scales as the pro…
▽ More
We discuss the role of tearing instabilities in magnetic reconnection. In three dimensions this instability leads to the formation of strong Alfvenic waves that remove plasma efficiently from the reconnection layer. As a result the instability proceeds at high rates while staying close to the linear regime. Our calculations show that for a resistive fluid the reconnection speed scales as the product of the Alfven speed V_A over the magnetic Reynolds number to the power -0.3. In the limit of vanishing resistivity, tearing modes proceed at a non-zero rate, driven by the electron inertia term, giving rise to a reconnection speed V_A (c/ω_p L_x)^{3/5}, where ω_p is the plasma frequency and L_x is the transverse scale of the reconnection layer. Formally this solves the problem of fast reconnection, but in practice this reconnection speed is small.
△ Less
Submitted 16 April, 1998;
originally announced April 1998.
-
Disk Instabilities and Cooling Fronts
Authors:
Ethan T. Vishniac
Abstract:
Accretion disk outbursts, and their subsequent decline, offer a unique opportunity to constrain the physics of angular momentum transport in hot accretion disks. Recent work has centered on the claim by Cannizzo et al. that the exponential decay of luminosity following an outburst in black hole accretion disk systems is only consistent with a particular form for the dimensionless viscosity,…
▽ More
Accretion disk outbursts, and their subsequent decline, offer a unique opportunity to constrain the physics of angular momentum transport in hot accretion disks. Recent work has centered on the claim by Cannizzo et al. that the exponential decay of luminosity following an outburst in black hole accretion disk systems is only consistent with a particular form for the dimensionless viscosity, $α=35(c_s/rΩ)^{3/2}$. This result can be understood in terms of a simple model of the evolution of cooling fronts in accretion disks. In particular, the cooling front speed during decline is $\sim α_F c_{s,F}(c_{s,F}/rΩ)^{n}$, where $F$ denotes the position of the cooling front, and the exact value of $n$ depends on the hot state opacity, (although generally $n\approx 1/2$). Setting this speed proportional to $r$ constrains the functional form of $α$ in the hot phase of the disk, which sets it apart from previous arguments based on the relative durations of outburst and quiescence. However, it remains uncertain how well we know the exponent $n$. In addition, more work is needed to clarify the role of irradiation in these systems and its effect on the cooling front evolution.
△ Less
Submitted 17 February, 1998;
originally announced February 1998.
-
Reconnection in the ISM
Authors:
Ethan T. Vishniac,
A. Lazarian
Abstract:
We discuss the role of ambipolar diffusion for simple reconnection in a partially ionized gas, following the reconnection geometry of Parker and Sweet. When the recombination time is short the mobility and reconnection of the magnetic field is substantially enhanced as matter escapes from the reconnection region via ambipolar diffusion. Our analysis shows that in the interstellar medium it is th…
▽ More
We discuss the role of ambipolar diffusion for simple reconnection in a partially ionized gas, following the reconnection geometry of Parker and Sweet. When the recombination time is short the mobility and reconnection of the magnetic field is substantially enhanced as matter escapes from the reconnection region via ambipolar diffusion. Our analysis shows that in the interstellar medium it is the recombination rate that usually limits the rate of reconnection. We show that heating effects can reduce this rate by increasing the recombination time and raising the local ion pressure. In the colder parts of the ISM, when temperatures are or the order of 100K or less, we obtain a significant enhancement over the usual Sweet-Parker rate, but only in dense molecular clouds will the reconnection velocity exceed 0.001 times the Alfvén speed. The ratio of the ion orbital radius to the reconnection layer thickness is typically a few percent, except in dense molecular clouds where it can approach unity. We briefly discuss prospects for obtaining much faster reconnection speeds in astrophysical plasmas.
△ Less
Submitted 4 December, 1997;
originally announced December 1997.
-
Torque Reversal in Accretion-Powered X-ray Pulsars
Authors:
Insu Yi,
J. Craig Wheeler,
Ethan T. Vishniac
Abstract:
Accretion-powered X-ray pulsars 4U 1626-67, GX 1+4, and OAO 1657-415 have recently shown puzzling torque reversals. These reversals are characterized by short time scales, on the order of days, nearly identical spin-up and spin-down rates, and very small changes in X-ray luminosity. We propose that this phenomenon is the result of sudden dynamical changes in the accretion disks triggered by a gr…
▽ More
Accretion-powered X-ray pulsars 4U 1626-67, GX 1+4, and OAO 1657-415 have recently shown puzzling torque reversals. These reversals are characterized by short time scales, on the order of days, nearly identical spin-up and spin-down rates, and very small changes in X-ray luminosity. We propose that this phenomenon is the result of sudden dynamical changes in the accretion disks triggered by a gradual variation of mass accretion rates. These sudden torque reversals may occur at a critical accretion rate $\sim$ $10^{15}-10^{16}g~s^{-1}$ when the system makes a transition from (to) a primarily Keplerian flow to (from) a substantially sub-Keplerian, radial advective flow in the inner disk. For systems near spin equilibrium, the spin-up torques in the Keplerian state are slightly larger than the spin-down torques in the advective state, in agreement with observation. The abrupt reversals could be a signature of pulsar systems near spin equilibrium with the mass accretion rates modulated on a time scale of a year or longer near the critical accretion rate. It is interesting that cataclysmic variables and black hole soft X-ray transients change their X-ray emission properties at accretion rates similar to the pulsars' critical rate. We speculate that the dynamical change in pulsar systems shares a common physical origin with white dwarf and black hole accretion disk systems.
△ Less
Submitted 14 October, 1997; v1 submitted 28 April, 1997;
originally announced April 1997.
-
The Structure of Cooling Fronts in Accretion Disks
Authors:
Ethan T. Vishniac
Abstract:
Recent work has shown that the speed of the cooling front in soft X-ray transients may be an important clue in understanding the nature of accretion disk viscosity. In a previous paper (Vishniac and Wheeler 1996) we derived the scaling law for the cooling front speed. Here we derive a similarity solution for the hot inner part of disks undergoing cooling. This solution is exact in the limit of a…
▽ More
Recent work has shown that the speed of the cooling front in soft X-ray transients may be an important clue in understanding the nature of accretion disk viscosity. In a previous paper (Vishniac and Wheeler 1996) we derived the scaling law for the cooling front speed. Here we derive a similarity solution for the hot inner part of disks undergoing cooling. This solution is exact in the limit of a thin disk, power law opacities, and a minimum hot state column density which is an infinitesimal fraction of the maximum cold state density. For a disk of finite thickness the largest error is in the ratio of the mass flow across the cooling front to the mass flow at small radii. Comparison to the numerical simulations of Cannizzo et al. (1995) inidcates that the errors in the other parameters do not exceed $(c_{sF}/r_FΩ_F)^q$, that is, the ratio of the sound speed at the disk midplane to its orbital velocity, evaluated at the cooling front, to the qth power. Here $q\approx 1/2$. Its precise value is determined by the relevant hot state opacity law and the functional form of the dimensionless viscosity.
△ Less
Submitted 9 January, 1997;
originally announced January 1997.
-
Gravitational Instability in Collisionless Cosmological Pancakes
Authors:
Azita Valinia,
Paul R. Shapiro,
Hugo Martel,
Ethan T. Vishniac
Abstract:
The gravitational instability of cosmological pancakes composed of collisionless dark matter in an Einstein-de Sitter universe is investigated numerically to demonstrate that pancakes are unstable with respect to fragmentation and the formation of filaments. A ``pancake'' is defined here as the nonlinear outcome of the growth of a 1D, sinusoidal, plane-wave, adiabatic density perturbation. We ha…
▽ More
The gravitational instability of cosmological pancakes composed of collisionless dark matter in an Einstein-de Sitter universe is investigated numerically to demonstrate that pancakes are unstable with respect to fragmentation and the formation of filaments. A ``pancake'' is defined here as the nonlinear outcome of the growth of a 1D, sinusoidal, plane-wave, adiabatic density perturbation. We have used high resolution, 2D, N-body simulations by the Particle-Mesh (PM) method to study the response of pancakes to perturbation by either symmetric (density) or antisymmetric (bending or rippling) modes, with corresponding wavevectors k_s and k_a transverse to the wavevector k_p of the unperturbed pancake plane-wave. We consider dark matter which is initially ``cold'' (i.e. with no random thermal velocity in the initial conditions). We also investigate the effect of a finite, random, isotropic, initial velocity dispersion (i.e. initial thermal velocity) on the fate of pancake collapse and instability. Pancakes are shown to be gravitationally unstable with respect to all perturbations of wavelength l<l_p (where l_p= 2pi/k_p). These results are in contradiction with the expectations of an approximate, thin-sheet energy argument.
△ Less
Submitted 28 October, 1996;
originally announced October 1996.
-
The Speed of Cooling Fronts and the Functional Form of the Dimensionless Viscosity in Accretion Disks
Authors:
Ethan T. Vishniac,
J. Craig Wheeler
Abstract:
We examine the speed of inward traveling cooling fronts in accretion disks. We show that their speed is determined by the rarefaction wave that precedes them and is approximately $α_F c_{F} (H/r)^q$, where $α_F$ is the dimensionless viscosity, $c_{F}$ is the sound speed, $r$ is the radial coordinate, $H$ is the disk thickness, and all quantities are evaluated at the cooling front. The scaling ex…
▽ More
We examine the speed of inward traveling cooling fronts in accretion disks. We show that their speed is determined by the rarefaction wave that precedes them and is approximately $α_F c_{F} (H/r)^q$, where $α_F$ is the dimensionless viscosity, $c_{F}$ is the sound speed, $r$ is the radial coordinate, $H$ is the disk thickness, and all quantities are evaluated at the cooling front. The scaling exponent $q$ lies in the interval $[0,1]$, depending on the slope of the $(T,Σ)$ relation in the hot state. For a Kramer's law opacity and $α\propto (H/r)^n$, where $n$ is of order unity, we find that $q\sim 1/2$. This supports the numerical work of Cannizzo, Chen and Livio (1995) and their conclusion that $n\approx3/2$ is necessary to reproduce the exponential decay of luminosity in black hole X-ray binary systems. Our results are insensitive to the structure of the disk outside of the radius where rapid cooling sets in. In particular, the width of the rapid cooling zone is a consequence of the cooling front speed rather than its cause. We conclude that the exponential luminosity decay of cooling disks is probably compatible with the wave-driven dynamo model. It is not compatible with models with separate, constant values of $α$ for the hot and cold states.
△ Less
Submitted 29 March, 1996;
originally announced March 1996.
-
The Variability of Active Galactic Nuclei and the Radial Transport of Vertical Magnetic Flux
Authors:
Seok Jae Park,
Ethan T. Vishniac
Abstract:
We consider the radial buoyancy of vertical magnetic field lines in radiation and gas pressure dominated accretion disks. We find that in addition to radial drift driven by turbulent diffusion and biased by the global field geometry, there are buoyancy effects which tend to move magnetic flux outward. In gas pressure dominated disks the poloidal magnetic field will move outward at a rate compara…
▽ More
We consider the radial buoyancy of vertical magnetic field lines in radiation and gas pressure dominated accretion disks. We find that in addition to radial drift driven by turbulent diffusion and biased by the global field geometry, there are buoyancy effects which tend to move magnetic flux outward. In gas pressure dominated disks the poloidal magnetic field will move outward at a rate comparable to its inward advection. On the other hand, in a radiation pressure dominated disk the poloidal magnetic field will usually move outward faster than it is advected inward. This implies that the fields in disks in active galactic nuclei are generated at small radii by an internal disk dynamo. This conclusion can be avoided if the external field imposes a supersonic Alfvén speed within the disk without giving rise to interchange instabilities. In any case we note that variations in the mass transfer rate will lead directly to a modulation of the nonthermal emission from the disk system.
△ Less
Submitted 16 May, 1996; v1 submitted 26 February, 1996;
originally announced February 1996.
-
An Incoherent $α-Ω$ Dynamo in Accretion Disks
Authors:
Ethan T. Vishniac,
Axel Brandenburg
Abstract:
We use the mean-field dynamo equations to show that an incoherent alpha effect in mirror-symmetric turbulence in a shearing flow can generate a large scale, coherent magnetic field. We illustrate this effect with simulations of a few simple systems. In accretion disks, this process can lead to axisymmetric magnetic domains whose radial and vertical dimensions will be comparable to the disk heigh…
▽ More
We use the mean-field dynamo equations to show that an incoherent alpha effect in mirror-symmetric turbulence in a shearing flow can generate a large scale, coherent magnetic field. We illustrate this effect with simulations of a few simple systems. In accretion disks, this process can lead to axisymmetric magnetic domains whose radial and vertical dimensions will be comparable to the disk height. This process may be responsible for observations of dynamo activity seen in simulations of dynamo-generated turbulence involving, for example, the Balbus-Hawley instability. In this case the magnetic field strength will saturate at $\sim (h/r)^2$ times the ambient pressure in real accretion disks. The resultant dimensionless viscosity will be of the same order. In numerical simulations the azimuthal extent of the simulated annulus should be substituted for $r$. We compare the predictions of this model to numerical simulations previously reported by Brandenburg et al. (1995). In a radiation pressure dominated environment this estimate for viscosity should be reduced by a factor of $(P_{gas}/P_{radiation})^6$ due to magnetic buoyancy.
△ Less
Submitted 6 October, 1995;
originally announced October 1995.
-
Global Aspects of Elliptical Instability in Tidally Distorted Accretion Disks
Authors:
Dongsu Ryu,
Jeremy Goodman,
Ethan T. Vishniac
Abstract:
Tidally distorted accretion disks in binary star systems are subject to a local hydrodynamic instability which excites $m=1$ internal waves. This instability is three dimensional and approximately incompressible. We study the global aspects of this local instability using equations derived under the shearing sheet approximation, where the effects of the azimuthal variation along distorted orbita…
▽ More
Tidally distorted accretion disks in binary star systems are subject to a local hydrodynamic instability which excites $m=1$ internal waves. This instability is three dimensional and approximately incompressible. We study the global aspects of this local instability using equations derived under the shearing sheet approximation, where the effects of the azimuthal variation along distorted orbital trajectories are included in source terms which oscillate with local orbital phase. Linear analyses show that the excitation of the instability is essentially local, i.e. insensitive to radial boundary conditions. The region of rapid growth feeds waves into the region of slow or negligible growth, allowing the instability to become global. The global growth rate depends the maximum local growth rate, the size of the rapid growth region, and the local group velocity. We present an empirical expression for the global growth rate. We note that the local nature of the instability allows the excitation of waves with $m\ne 1$ when the local growth rate is large.
△ Less
Submitted 4 September, 1995;
originally announced September 1995.
-
Magnetohydrodynamic Turbulence in Accretion Disks
Authors:
Ethan T. Vishniac
Abstract:
Recent work on the structure of magnetic fields in a turbulent medium gives predictions for the properties of the magnetic flux tubes as a function of the Mach number and scale of the turbulence, and the resistivity and viscosity of the fluid. Here I discuss the implications of this work for accretion disks. I show that although accretion disk flux tubes are usually almost completely evacuated,…
▽ More
Recent work on the structure of magnetic fields in a turbulent medium gives predictions for the properties of the magnetic flux tubes as a function of the Mach number and scale of the turbulence, and the resistivity and viscosity of the fluid. Here I discuss the implications of this work for accretion disks. I show that although accretion disk flux tubes are usually almost completely evacuated, they are nevertheless less buoyant than previous estimates have suggested. I also note that vertical magnetic flux tends to be ejected from the outer edge of accretion disks, so evidence for continued magnetic activity in such systems should be interpreted as supporting the existence of dynamo activity.
△ Less
Submitted 18 August, 1995;
originally announced August 1995.
-
The Growth Rate of Tidally Excited Waves in Accretion Disks
Authors:
Ethan T. Vishniac,
Changsong Zhang
Abstract:
Accretion disks in close binary systems are subject to a tidally driven parametric instability which leads to the growth of internal waves near the outer edges of such disks (Goodman 1993). These waves are important in understanding the torque exerted on a disk by tidal forces and may play a role in the structure of the disk at small radii. Here we calculate the growth rate of this instability,…
▽ More
Accretion disks in close binary systems are subject to a tidally driven parametric instability which leads to the growth of internal waves near the outer edges of such disks (Goodman 1993). These waves are important in understanding the torque exerted on a disk by tidal forces and may play a role in the structure of the disk at small radii. Here we calculate the growth rate of this instability, including the effects of vertical structure and fluid compressibility. We find growth rates which are only slightly different from Goodman's original results, except that near the vertical resonance radius the growth rate can have an extremely broad and strong peak when the disk is stably stratified in the vertical direction. Higher order modes, in the sense of increasing number of vertical nodes, have similar growth rates. Our results differ from a previous calculation along these lines by Lubow et al. (1993). The difference is mostly due to their neglect of radial gradients in the tidally distorted streamlines.
△ Less
Submitted 8 August, 1995;
originally announced August 1995.
-
The Dynamics of Flux Tubes in a High Beta Plasma
Authors:
E. T. Vishniac
Abstract:
We suggest a new model for the structure of a magnetic field embedded high $β$ turbulent plasma, based on the popular notion that the magnetic field will tend to separate into individual flux tubes. We point out that interactions between the flux tubes will be dominated by coherent effects stemming from the turbulent wakes created as the fluid streams by the flux tubes. Balancing the attraction…
▽ More
We suggest a new model for the structure of a magnetic field embedded high $β$ turbulent plasma, based on the popular notion that the magnetic field will tend to separate into individual flux tubes. We point out that interactions between the flux tubes will be dominated by coherent effects stemming from the turbulent wakes created as the fluid streams by the flux tubes. Balancing the attraction caused by shielding effects with turbulent diffusion we find that flux tubes have typical radii comparable to the local Mach number squared times the large scale eddy length, are arranged in a one dimensional fractal pattern, have a radius of curvature comparable to the largest scale eddies in the turbulence, and have an internal magnetic pressure comparable to the ambient pressure. When the average magnetic energy density is much less than the turbulent energy density the radius, internal magnetic field and curvature scale of the flux tubes will be smaller than these estimates. Realistic resistivity does not alter the macroscopic properties of the fluid or the large scale magnetic field. In either case we show that the Sweet-Parker reconnection rate is much faster than an eddy turnover time. Realistic stellar plasmas are expected to either be in the ideal limit (e.g. the solar photosphere) or the resistive limit (most of the solar convection zone). All current numerical simulations of three dimensional MHD turbulence are in the viscous regime and are inapplicable to stars or accretion disks.
△ Less
Submitted 21 July, 1994;
originally announced July 1994.
-
Why Disks Shine: the Transport of Angular Momentum in Hot, Thin Disks
Authors:
E. T. Vishniac
Abstract:
I review recent work on the radial transport of angular momentum in ionized, Keplerian accretion disks. Proposed mechanisms include hydrodynamic and MHD local instabilities and long range effects mediated by wave transport. The most promising models incorporate the Velikhov-Chandrasekhar instability, caused by an instability of the magnetic field embedded in a differentially rotating disk. This…
▽ More
I review recent work on the radial transport of angular momentum in ionized, Keplerian accretion disks. Proposed mechanisms include hydrodynamic and MHD local instabilities and long range effects mediated by wave transport. The most promising models incorporate the Velikhov-Chandrasekhar instability, caused by an instability of the magnetic field embedded in a differentially rotating disk. This has the important feature that the induced turbulent motions necessarily transport angular momentum outward. By contrast, convective modes may transport angular momentum in either direction. Combining the magnetic field instability with an $α-Ω$ dynamo driven by internal waves leads to a model in which the dimensionless viscosity scales as $(H/r)^{4/3}$. However, this model has a phenomenology which is quite different from the $α$ disk model. For example, an active disk implies some source of excitation for the internal waves. In binary systems with a mass ratio of order unity the most likely exciting mechanism is a parametric instability due to tidal forces. This implies that in systems where the accretion stream is intermittent, like MV Lyrae or TT Ari, epochs when the mass flow is absent or very small will be epochs in which the disk shrinks and becomes relatively inactive and dark. This model also implies that forced vertical mixing is important, even in convectively stable disks. I discuss various observational tests of this model and the focus of current theoretical work.
△ Less
Submitted 28 September, 1993;
originally announced September 1993.
-
Does Dissipation in AGN Disks Couple to the Total Pressure?
Authors:
E. T. Vishniac
Abstract:
Recent work on the transport of angular momentum in accretion disks suggests that the Velikhov-Chandrasekhar instability, in which a large scale magnetic field generates small scale eddys in a shearing environment, may be ultimately responsible for this process. Although there is considerable controversy about the origin and maintenance of this field in accretion disks, it turns out that it is p…
▽ More
Recent work on the transport of angular momentum in accretion disks suggests that the Velikhov-Chandrasekhar instability, in which a large scale magnetic field generates small scale eddys in a shearing environment, may be ultimately responsible for this process. Although there is considerable controversy about the origin and maintenance of this field in accretion disks, it turns out that it is possible to argue, quite generally, using scaling arguments, that this process is sensitive to the total pressure in an AGN disk, rather than the pressure contributed by gas alone. We conclude that the resolution of the conceptual difficulties implied by the presence of strong thermal and viscous instabilities in radiation pressure and electron scattering dominated does not lie in models that couple the total dissipation rate to the gas pressure alone, or to some weighted mean of the gas and radiation pressures.
△ Less
Submitted 12 August, 1993;
originally announced August 1993.
-
On the Origin of the Galactic Magnetic Field
Authors:
Ethan T. Vishniac
Abstract:
The galactic magnetic field is commonly supposed to be due to a dynamo acting on some large scale seed field. A major difficulty with this idea is that estimates of reasonable seed field strengths tend to be quite low, on the order of $\sim10^{-20}$ gauss. Here we examine the contribution due to the flux entrained in winds from protostars formed in the first dynamo e-folding time of a galaxy's e…
▽ More
The galactic magnetic field is commonly supposed to be due to a dynamo acting on some large scale seed field. A major difficulty with this idea is that estimates of reasonable seed field strengths tend to be quite low, on the order of $\sim10^{-20}$ gauss. Here we examine the contribution due to the flux entrained in winds from protostars formed in the first dynamo e-folding time of a galaxy's existence. Using a minimal estimate of a protostellar magnetic field we find that if each protostar ejects a single current ring, sufficient to maintain flux freezing in the wind, than the large scale average dipole field from all such current rings will be at least 5 orders of magnitude larger than previous seed field estimates. Allowing for a reasonable amount of magnetic activity in protostars during an extended period of mass loss increases this to a dipole seed field of $\sim10^{-12}$ gauss. For the purposes of producing a seed field it is irrelevant whether or not this initial injection of flux takes place in a newly formed galactic disk, or in star forming proto-galactic clouds. The compression of this dipole field into a thin disk will lead to a large scale $B_r\sim 10^{-10.5}$ gauss. Initially, field strengths on smaller scales will be larger, but nowhere near current levels.
△ Less
Submitted 30 October, 1992;
originally announced October 1992.