Showing 1–50 of 70 results for author: Blackman, E G

Thermal and magnetic evolution of an Earth-like planet with a basal magma ocean

Authors: Victor Lherm, Miki Nakajima, Eric G. Blackman

Abstract: Earth's geodynamo has operated for over 3.5 billion years. The magnetic field is currently powered by thermocompositional convection in the outer core, which involves the release of light elements and latent heat as the inner core solidifies. However, since the inner core nucleated no more than 1.5 billion years ago, the early dynamo could not rely on these buoyancy sources. Given recent estimates… ▽ More Earth's geodynamo has operated for over 3.5 billion years. The magnetic field is currently powered by thermocompositional convection in the outer core, which involves the release of light elements and latent heat as the inner core solidifies. However, since the inner core nucleated no more than 1.5 billion years ago, the early dynamo could not rely on these buoyancy sources. Given recent estimates of the thermal conductivity of the outer core, an alternative mechanism may be required to sustain the geodynamo prior to nucleation of the inner core. One possibility is a silicate dynamo operating in a long-lived basal magma ocean. Here, we investigate the structural, thermal, buoyancy, and magnetic evolution of an Earth-like terrestrial planet. Using modern equations of state and melting curves, we include a time-dependent parameterization of the compositional evolution of an iron-rich basal magma ocean. We combine an internal structure integration of the planet with energy budgets in a coupled core, basal magma ocean, and mantle system. We determine the thermocompositional convective stability of the core and the basal magma ocean, and assess their respective dynamo activity using entropy budgets and magnetic Reynolds numbers. Our conservative nominal model predicts a transient basal magma ocean dynamo followed by a core dynamo after 1 billion years. The model is sensitive to several parameters, including the initial temperature of the core-mantle boundary, the parameterization of mantle convection, the composition of the basal magma ocean, the radiogenic content of the planet, as well as convective velocity and magnetic scaling laws. We use the nominal model to constrain the range of basal magma ocean electrical conductivity and core thermal conductivity that sustain a dynamo. △ Less

Submitted 9 September, 2024; originally announced September 2024.

Structure and Dynamics of Magneto-Inertial, Differentially Rotating Laboratory Plasmas

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, D. R. Russell, S. Merlini, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise-time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set… ▽ More We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise-time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum, and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at $\approx 100$ km/s. The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$. As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value $\approx 20$ km/s. We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements. △ Less

Submitted 29 March, 2024; originally announced March 2024.

Comments: Submitted to the Journal of Plasma Physics. 22 pages, 9 figures

On The Structure of Plasma Jets in the Rotating Plasma Experiment

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, D. R. Russell, S. Merlini, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: Recent pulsed-power experiments have demonstrated the formation of astrophysically-relevant, differentially rotating plasmas [1]. Key features of the plasma flows are the discovery of a quasi-Keplerian rotation curve, the launching of highly-collimated angular-momentum-transporting axial jets, and a hollow density structure sustained by the centrifugal barrier effect. In this communication we disc… ▽ More Recent pulsed-power experiments have demonstrated the formation of astrophysically-relevant, differentially rotating plasmas [1]. Key features of the plasma flows are the discovery of a quasi-Keplerian rotation curve, the launching of highly-collimated angular-momentum-transporting axial jets, and a hollow density structure sustained by the centrifugal barrier effect. In this communication we discuss several features of the plasma structure in these experiments through order-of-magnitude models. First, we show that the observed rotation velocity would produce a centrifugal force strong enough to support the hollow density profile. Second, we show that the axial jet should diverge much faster than what was observed, were it not for a magnetized halo with 3T which surrounds the jet and exerts pressure on the interface.Finally, we discuss the temperature structure in the axial jet and plasma halo.We show that a 3T magnetic field would also suppress electron heat conduction,leading to the flat profile observed experimentally. We also find that the axial jet is efficiently radiatively cooled,whereas the halo is not, which would explain the thermal decoupling between the two regions. △ Less

Submitted 15 February, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: Second version with clarifications on text and figures, and explicit calculations regarding the omega-effect. Manuscript re-submitted to IEEE Transactions On Plasma Science. 4 Figures

Cooling and Instabilities in Colliding Radiative Flows with Toroidal Magnetic Fields

Authors: R. N. Markwick, A. Frank, E. G. Blackman, J. Carroll-Nellenback, S. V. Lebedev, D. R. Russell, J. W. D. Halliday, L. G. Suttle, P. M. Hartigan

Abstract: We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calcul… ▽ More We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calculation. We systematically vary plasma beta ($βべーた_m$) in the flows, the strength of the cooling ($Λらむだ_0$) and the exponent ($αあるふぁ$) of temperature-dependence of the cooling function. We find that for strong magnetic fields a counter-propagating jet called a "spine" is driven by pressure from shocked toroidal fields. The spines eventually become unstable and break apart. We demonstrate how formation and evolution of the spines depends on initial flow parameters and provide a simple analytic model that captures the basic features of the flow. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 14 pages, 16 figures. Submitted to MNRAS

Morphology of Shocked Lateral Outflows in Colliding Hydrodynamic Flows

Authors: R. N. Markwick, A. Frank, J. Carroll-Nellenback, E. G. Blackman, P. M. Hartigan, S. V. Lebedev, D. R. Russel, J. W. D. Halliday, L. G. Suttle

Abstract: Supersonic interacting flows occurring in phenomena such as protostellar jets give rise to strong shocks, and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross sectional radius of the colliding f… ▽ More Supersonic interacting flows occurring in phenomena such as protostellar jets give rise to strong shocks, and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. We introduce variations in the flow parameters of density, velocity, and cross sectional radius of the colliding flows %radius in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows (SLOs) being deflected from the horizontal when the flows are of differing cross-section. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Journal ref: Physics of Plasmas, Volume 29, Issue 10, id.102901, 14 pp (October 2022)

Ion and Electron Acoustic Bursts during Anti-Parallel Magnetic Reconnection Driven by Lasers

Authors: Shu Zhang, Abraham Chien, Lan Gao, Hantao Ji, Eric G. Blackman, Russ Follett, Dustin H. Froula, Joseph Katz, Chikang Li, Andrew Birkel, Richard Petrasso, John Moody, Hui Chen

Abstract: Magnetic reconnection converts magnetic energy into thermal and kinetic energy in plasma. Among numerous candidate mechanisms, ion acoustic instabilities driven by the relative drift between ions and electrons, or equivalently electric current, have been suggested to play a critical role in dissipating magnetic energy in collisionless plasmas. However, their existence and effectiveness during reco… ▽ More Magnetic reconnection converts magnetic energy into thermal and kinetic energy in plasma. Among numerous candidate mechanisms, ion acoustic instabilities driven by the relative drift between ions and electrons, or equivalently electric current, have been suggested to play a critical role in dissipating magnetic energy in collisionless plasmas. However, their existence and effectiveness during reconnection have not been well understood due to ion Landau damping and difficulties in resolving the Debye length scale in the laboratory. Here we report a sudden onset of ion acoustic bursts measured by collective Thomson scattering in the exhaust of anti-parallel magnetically driven reconnection using high-power lasers. The ion acoustic bursts are followed by electron acoustic bursts with electron heating and bulk acceleration. We reproduce these observations with 1D and 2D particle-in-cell simulations in which electron outflow jet drives ion-acoustic instabilities, forming double layers. These layers induce electron two-stream instabilities that generate electron acoustic bursts and energize electrons. Our results demonstrate the importance of ion and electron acoustic dynamics during reconnection when ion Landau damping is ineffective, a condition applicable to a range of astrophysical plasmas including near-Earth space, stellar flares, and black hole accretion engines. △ Less

Submitted 29 March, 2023; v1 submitted 26 September, 2022; originally announced September 2022.

Characterization of Quasi-Keplerian, Differentially Rotating, Free-Boundary Laboratory Plasmas

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, S. Merlini, D. R. Russell, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: We present results from pulsed-power driven differentially rotating plasma experiments designed to simulate physics relevant to astrophysical disks and jets. In these experiments, angular momentum is injected by the ram pressure of the ablation flows from a wire array Z pinch. In contrast to previous liquid metal and plasma experiments, rotation is not driven by boundary forces. Axial pressure gra… ▽ More We present results from pulsed-power driven differentially rotating plasma experiments designed to simulate physics relevant to astrophysical disks and jets. In these experiments, angular momentum is injected by the ram pressure of the ablation flows from a wire array Z pinch. In contrast to previous liquid metal and plasma experiments, rotation is not driven by boundary forces. Axial pressure gradients launch a rotating plasma jet upwards, which is confined by a combination of ram, thermal, and magnetic pressure of a surrounding plasma halo. The jet has subsonic rotation, with a maximum rotation velocity $23 \pm 3$ km/s. The rotational velocity profile is quasi-Keplerian with a positive Rayleigh discriminant $κかっぱ^2 \propto r^{-2.8\pm0.8}$ rad$^2$/s$^2$. The plasma completes $0.5 - 2$ full rotations in the experimental time frame ($\sim 150$ ns). △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 4 figures

Direct measurement of non-thermal electron acceleration from magnetically driven reconnection in a laboratory plasma

Authors: Abraham Chien, Lan Gao, Shu Zhang, Hantao Ji, Eric G. Blackman, William Daughton, Adam Stanier, Ari Le, Fan Guo, Russ Follett, Hui Chen, Gennady Fiksel, Gabriel Bleotu, Robert C. Cauble, Sophia N. Chen, Alice Fazzini, Kirk Flippo, Omar French, Dustin H. Froula, Julien Fuchs, Shinsuke Fujioka, Kenneth Hill, Sallee Klein, Carolyn Kuranz, Philip Nilson , et al. (2 additional authors not shown)

Abstract: Magnetic reconnection is a ubiquitous astrophysical process that rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy, and non-thermal energetic particles, including energetic electrons. Various reconnection acceleration mechanisms in different low-$βべーた$ (plasma-to-magnetic pressure ratio) and collisionless environments have been proposed theoretically and stu… ▽ More Magnetic reconnection is a ubiquitous astrophysical process that rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy, and non-thermal energetic particles, including energetic electrons. Various reconnection acceleration mechanisms in different low-$βべーた$ (plasma-to-magnetic pressure ratio) and collisionless environments have been proposed theoretically and studied numerically, including first- and second-order Fermi acceleration, betatron acceleration, parallel electric field acceleration along magnetic fields, and direct acceleration by the reconnection electric field. However, none of them have been heretofore confirmed experimentally, as the direct observation of non-thermal particle acceleration in laboratory experiments has been difficult due to short Debye lengths for \textit{in-situ} measurements and short mean free paths for \textit{ex-situ} measurements. Here we report the direct measurement of accelerated non-thermal electrons from low-$βべーた$ magnetically driven reconnection in experiments using a laser-powered capacitor coil platform. We use kiloJoule lasers to drive parallel currents to reconnect MegaGauss-level magnetic fields in a quasi-axisymmetric geometry. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that the mechanism of direct electric field acceleration by the out-of-plane reconnection electric field is at work. Scaled energies using this mechanism show direct relevance to astrophysical observations. Our results therefore validate one of the proposed acceleration mechanisms by reconnection, and establish a new approach to study reconnection particle acceleration with laboratory experiments in relevant regimes. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 13 pages, 6 figures

Cooling and Instabilities in Colliding Flows

Authors: R. N. Markwick, A. Frank, J. Carroll-Nellenback, B. Liu, E. G. Blackman, S. V. Lebedev, P. M. Hartigan

Abstract: Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature… ▽ More Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power law cooling function. The characteristic length and time scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently rapidly the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one and three dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 12 pages, 21 figures. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 2, pp.2266-2278 (Dec 2021)

Scaling of Turbulent Viscosity and Resistivity: Extracting a Scale-dependent Turbulent Magnetic Prandtl Number

Authors: Xin Bian, Jessica K. Shang, Eric G. Blackman, Gilbert W. Collins, Hussein Aluie

Abstract: Turbulent viscosity $νにゅー_t$ and resistivity $ηいーた_t$ are perhaps the simplest models for turbulent transport of angular momentum and magnetic fields, respectively. The associated turbulent magnetic Prandtl number $Pr_t\equiv νにゅー_t/ηいーた_t$ has been well recognized to determine the final magnetic configuration of accretion disks. Here, we present an approach to determining these ''effective transport'' coeffi… ▽ More Turbulent viscosity $νにゅー_t$ and resistivity $ηいーた_t$ are perhaps the simplest models for turbulent transport of angular momentum and magnetic fields, respectively. The associated turbulent magnetic Prandtl number $Pr_t\equiv νにゅー_t/ηいーた_t$ has been well recognized to determine the final magnetic configuration of accretion disks. Here, we present an approach to determining these ''effective transport'' coefficients acting at different length-scales using coarse-graining and recent results on decoupled kinetic and magnetic energy cascades [Bian & Aluie 2019]. By analyzing the kinetic and magnetic energy cascades from a suite of high-resolution simulations, we show that our definitions of $νにゅー_t$, $ηいーた_t$, and $Pr_t$ have power-law scalings in the ''decoupled range.'' We observe that $Pr_t\approx1 \text{~to~}2$ at the smallest inertial-inductive scales, increasing to $\approx 5$ at the largest scales. However, based on physical considerations, our analysis suggests that $Pr_t$ has to become scale-independent and of order unity in the decoupled range at sufficiently high Reynolds numbers (or grid-resolution), and that the power-law scaling exponents of velocity and magnetic spectra become equal. In addition to implications to astrophysical systems, the scale-dependent turbulent transport coefficients offer a guide for large eddy simulation modeling. △ Less

Submitted 2 July, 2021; originally announced July 2021.

On the shear-current effect: toward understanding why theories and simulations have mutually and separately conflicted

Authors: Hongzhe Zhou, Eric G. Blackman

Abstract: The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient $βべーた_{21}$ with a favorable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic (… ▽ More The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient $βべーた_{21}$ with a favorable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic ($βべーた^u_{21}$) and magnetic ($βべーた^b_{21}$) contributions. To resolve these discrepancies, we combine an analytical approach with simulations, and show that unlike $βべーた^b_{21}$, the kinetic SCE $βべーた^u_{21}$ has a strong dependence on the kinetic energy spectral index and can transit from positive to negative values at $\mathcal{O}(10)$ Reynolds numbers if the spectrum is not too steep. Conversely, $βべーた^b_{21}$ is always negative regardless of the spectral index and Reynolds numbers. For very steep energy spectra, the positive $βべーた^u_{21}$ can dominate even at energy equipartition $u_\text{rms}\simeq b_\text{rms}$, resulting in a positive total $βべーた_{21}$ even though $βべーた^b_{21}<0$. Our findings bridge the gap between the seemingly contradictory results from the second-order-correlation approximation (SOCA) versus the spectral-$τたう$ closure (STC), for which opposite signs for $βべーた^u_{21}$ have been reported, with the same sign for $βべーた^b_{21}<0$. The results also offer an explanation for the simulations that find $βべーた^u_{21}>0$ and an inconclusive overall sign of $βべーた_{21}$ for $\mathcal{O}(10)$ Reynolds numbers. The transient behavior of $βべーた^u_{21}$ is demonstrated using the kinematic test-field method. We compute dynamo growth rates for cases with or without rotation, and discuss opportunities for further work. △ Less

Submitted 13 September, 2021; v1 submitted 22 April, 2021; originally announced April 2021.

Comments: 17 pages, 12 figures; accepted by MNRAS

Influence of inhomogeneous stochasticity on the falsifiability of mean-field theories and examples from accretion disc modeling

Authors: Hongzhe Zhou, Eric G. Blackman

Abstract: Despite spatial and temporal fluctuations in turbulent astrophysical systems, mean-field theories can be used to describe their secular evolution. However, observations taken over time scales much shorter than dynamical time scales capture a system in a single state of its turbulence ensemble. Comparing with mean-field theory can falsify the latter only if the theory is additionally supplied with… ▽ More Despite spatial and temporal fluctuations in turbulent astrophysical systems, mean-field theories can be used to describe their secular evolution. However, observations taken over time scales much shorter than dynamical time scales capture a system in a single state of its turbulence ensemble. Comparing with mean-field theory can falsify the latter only if the theory is additionally supplied with a quantified precision. The central limit theorem provides appropriate estimates to the precision only when fluctuations contribute linearly to an observable and with constant coherent scales. Here we introduce an error propagation formula that relaxes both limitations, allowing for nonlinear functional forms of observables and inhomogeneous coherent scales and amplitudes of fluctuations. The method is exemplified in the context of accretion disc theories, where inhomogeneous fluctuations in the surface temperature are propagated to the disc emission spectrum--the latter being a nonlinear and non-local function of the former. The derived precision depends non-monotonically on emission frequency. Using the same method, we investigate how binned spectral fluctuations in telescope data change with the spectral resolving power. We discuss the broader implications for falsifiability of a mean-field theory. △ Less

Submitted 30 August, 2021; v1 submitted 12 December, 2020; originally announced December 2020.

Comments: 10 pages, 4 figures; accepted by MNRAS

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen… ▽ More Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

Persistent mysteries of jet engines, formation, propagation, and particle acceleration: have they been addressed experimentally?

Authors: Eric G. Blackman, Sergey V. Lebedev

Abstract: The physics of astrophysical jets can be divided into three regimes: (i) engine and launch (ii) propagation and collimation, (iii) dissipation and particle acceleration. Since astrophysical jets comprise a huge range of scales and phenomena, practicality dictates that most studies of jets intentionally or inadvertently focus on one of these regimes, and even therein, one body of work may be simply… ▽ More The physics of astrophysical jets can be divided into three regimes: (i) engine and launch (ii) propagation and collimation, (iii) dissipation and particle acceleration. Since astrophysical jets comprise a huge range of scales and phenomena, practicality dictates that most studies of jets intentionally or inadvertently focus on one of these regimes, and even therein, one body of work may be simply boundary condition for another. We first discuss long standing persistent mysteries that pertain the physics of each of these regimes, independent of the method used to study them. This discussion makes contact with frontiers of plasma astrophysics more generally. While observations theory, and simulations, and have long been the main tools of the trade, what about laboratory experiments? Jet related experiments have offered controlled studies of specific principles, physical processes, and benchmarks for numerical and theoretical calculations. We discuss what has been done to date on these fronts. Although experiments have indeed helped us to understand certain processes, proof of principle concepts, and benchmarked codes, they have yet to solved an astrophysical jet mystery on their own. A challenge is that experimental tools used for jet-related experiments so far, are typically not machines originally designed for that purpose, or designed with specific astrophysical mysteries in mind. This presents an opportunity for a different way of thinking about the development of future platforms: start with the astrophysical mystery and build an experiment to address it. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 37 pages, including 11 figures at the end; for New Astronomy Reviews special volume '100 Years of Jets' (eds. Rob Fender and Ralph Wijers)

Necessity of ventilation for mitigating virus transmission quantified simply

Authors: Eric G. Blackman, Gourab Ghoshal

Abstract: To mitigate the SARS-CoV-2 pandemic, officials have employed social distancing and stay-at-home measures, with increased attention to room ventilation emerging only more recently. Effective distancing practices for open spaces can be ineffective for poorly ventilated spaces, both of which are commonly filled with turbulent air. This is typical for indoor spaces that use mixing ventilation. While t… ▽ More To mitigate the SARS-CoV-2 pandemic, officials have employed social distancing and stay-at-home measures, with increased attention to room ventilation emerging only more recently. Effective distancing practices for open spaces can be ineffective for poorly ventilated spaces, both of which are commonly filled with turbulent air. This is typical for indoor spaces that use mixing ventilation. While turbulence initially reduces the risk of infection near a virion-source, it eventually increases the exposure risk for all occupants in a space without ventilation. To complement detailed models aimed at precision, minimalist frameworks are useful to facilitate order of magnitude estimates for how much ventilation provides safety, particularly when circumstances require practical decisions with limited options. Applying basic principles of transport and diffusion, we estimate the time-scale for virions injected into a room of turbulent air to infect an occupant, distinguishing cases of low vs. high initial virion mass loads and virion-destroying vs. virion-reflecting walls. We consider the effect of an open window as a proxy for ventilation. When the airflow is dominated by isotropic turbulence, the minimum area needed to ensure safety depends only on the ratio of total viral load to threshold load for infection. The minimalist estimates here convey simply that the equivalent of ventilation by modest sized open window in classrooms and workplaces significantly improves safety. △ Less

Submitted 18 August, 2020; v1 submitted 20 June, 2020; originally announced June 2020.

Comments: 8 pages, 1 table, 1 figure, submitted

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena throughout the Universe

Authors: H. Ji, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, L. -J. Chen, Y. Chen, A. Chien, D. Craig, J. Dahlin, W. Daughton, E. DeLuca, C. F. Dong, S. Dorfman, J. Drake, F. Ebrahimi , et al. (75 additional authors not shown)

Abstract: This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 6 pages, 1 figure, white paper submitted to both Plasma 2020 and Astro 2020 Decadal Surveys

Improving TBI Protection Measures and Standards for Combat Helmets

Authors: Eric G. Blackman

Abstract: How well do current combat helmets protect against Traumatic Brain Injury? To answer this question it is necessary to evaluate both the measures and standards of protection. We define measure as a physical test that the helmet must be subject to and the standard as the quantitative threshold of performance in this test that the helmet must satisfy to be acceptable. We find that both the measures a… ▽ More How well do current combat helmets protect against Traumatic Brain Injury? To answer this question it is necessary to evaluate both the measures and standards of protection. We define measure as a physical test that the helmet must be subject to and the standard as the quantitative threshold of performance in this test that the helmet must satisfy to be acceptable. We find that both the measures are inadequate and the standards are too low. Current combat helmets simply do not protect against TBI effectively as they have not been designed to do so. We explain this inadequacy of combat helmet protection against TBI, and how the situation can be dramatically improved both in the short and long term. △ Less

Submitted 14 March, 2020; originally announced March 2020.

Comments: Report of archival value from 2006-2007 (small update 2011) to the Defense Science Study Group (institute for Defense Analyses, Alexandria VA); 32 pages

Generalized quenching of large-scale magnetic dynamos in anisotropic flows

Authors: Hongzhe Zhou, Eric G. Blackman

Abstract: The buildup of small-scale magnetic helicity which accompanies the oppositely signed growth on large scales is central to conventional dynamical quenching theories of mean-field dynamos. However, the conventional formalism presumes isotropic turbulence and thereby excludes part of the magnetic Lorentz back-reaction. This renders it insufficient to predict the full quenching for general anisotropic… ▽ More The buildup of small-scale magnetic helicity which accompanies the oppositely signed growth on large scales is central to conventional dynamical quenching theories of mean-field dynamos. However, the conventional formalism presumes isotropic turbulence and thereby excludes part of the magnetic Lorentz back-reaction. This renders it insufficient to predict the full quenching for general anisotropic flows. To overcome this deficiency, we derive a new generalized quenching formalism that includes the full back-reacting Lorentz force, and a new "selective-damping-$τたう$" closure which conserves magnetic helicity. We apply the formalism to examples of $\bmαあるふぁ^2$ dynamos and show its predicted quenching for different cases of turbulence---isotropic helical, anisotropic helical, and anisotropic non-helical. It predicts stronger-than-conventional quenching in general, but reduces to the conventional case in the helical isotropic limit. △ Less

Submitted 9 June, 2020; v1 submitted 3 May, 2019; originally announced May 2019.

Comments: 16 pages, 2 figures; submitted to Journal of Plasma Physics

Solving the Riemann Problem for Realistic Astrophysical Fluids

Authors: Zhuo Chen, Matthew S. B. Coleman, Eric G. Blackman, Adam Frank

Abstract: We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the… ▽ More We present new methods to solve the Riemann problem both exactly and approximately for general equations of state (EoS) to facilitate realistic modeling and understanding of astrophysical flows. The existence and uniqueness of the new exact general EoS Riemann solution can be guaranteed if the EoS is monotone regardless of the physical validity of the EoS. We confirm that: (1) the solution of the new exact general EoS Riemann solver and the solution of the original exact Riemann solver match when calculating perfect gas Euler equations; (2) the solution of the new Harten-Lax-van Leer-Contact (HLLC) general EoS Riemann solver and the solution of the original HLLC Riemann solver match when working with perfect gas EoS; and (3) the solution of the new HLLC general EoS Riemann solver approaches the new exact solution. We solve the EoS with two methods, one is to interpolate 2D EoS tables by the bi-linear interpolation method, and the other is to analytically calculate thermodynamic variables at run-time. The interpolation method is more general as it can work with other monotone and realistic EoS while the analytic EoS solver introduced here works with a relatively idealized EoS. Numerical results confirm that the accuracy of the two EoS solvers is similar. We study the efficiency of these two methods with the HLLC general EoS Riemann solver and find that analytic EoS solver is faster in the test problems. However, we point out that a combination of the two EoS solvers may become favorable in some specific problems. Throughout this research, we assume local thermal equilibrium. △ Less

Submitted 2 November, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 42 pages, 19 figures, accepted by the Journal of Computational Physics. Corrected some grammatical errors and a typo in the latest version

Minimalist large scale dynamo from shear-driven inhomogeneity

Authors: F. Ebrahimi, E. G. Blackman

Abstract: We show that non-axisymmetric, non-helical perturbations in an unstratified shear flow produce a shear-plane averaged electromotive force (EMF) proportional to a spatially dependent kinetic helicity. This new "shear-driven $αあるふぁ$-effect" can amplify shear-plane averaged magnetic fields, but is undiscoverable when invoking the homogeneity assumption for dynamo coefficients. It is a candidate mechanism… ▽ More We show that non-axisymmetric, non-helical perturbations in an unstratified shear flow produce a shear-plane averaged electromotive force (EMF) proportional to a spatially dependent kinetic helicity. This new "shear-driven $αあるふぁ$-effect" can amplify shear-plane averaged magnetic fields, but is undiscoverable when invoking the homogeneity assumption for dynamo coefficients. It is a candidate mechanism underlying large scale field generation in astrophysical, laboratory, and simulated shear flows, and highlights the essentiality of correctly averaging dynamo coefficients and fields. △ Less

Submitted 5 May, 2019; v1 submitted 12 February, 2019; originally announced February 2019.

Comments: 4 pages

Calculating turbulent transport tensors by averaging single plume dynamics and application to dynamos

Authors: Hongzhe Zhou, Eric G. Blackman

Abstract: Transport coefficients in turbulence are comprised of correlation functions between turbulent fluctuations and efficient methods to calculate them are desirable. For example, in mean field dynamo theories used to model the growth of large scale magnetic fields of stars and galaxies, the turbulent electromotive force is commonly approximated by a series of tensor products of turbulent transport coe… ▽ More Transport coefficients in turbulence are comprised of correlation functions between turbulent fluctuations and efficient methods to calculate them are desirable. For example, in mean field dynamo theories used to model the growth of large scale magnetic fields of stars and galaxies, the turbulent electromotive force is commonly approximated by a series of tensor products of turbulent transport coefficients with successively higher order spatial derivatives of the mean magnetic field. One ingredient of standard models is the kinematic coefficient of the zeroth order term, namely the averaged kinetic pseudotensor $\bmαあるふぁ$, that converts toroidal to poloidal fields. Here we demonstrate an efficient way to calculate this quantity for rotating stratified turbulence, whereby the pre-averaged quantity is calculated for the motion of a single plume, and the average is then taken over an ensemble of plumes of different orientations. We calculate the plume dynamics in the most convenient frame, before transforming back to the lab frame and averaging. Our concise configuration space calculation gives essentially identical results to previous lengthier approaches. The present application exemplifies what is a broadly applicable method. △ Less

Submitted 28 November, 2018; v1 submitted 16 October, 2018; originally announced October 2018.

Comments: 6 pages; 1 figure; accepted by MNRAS Letters

Experimental Confirmation of the Standard Magnetorotational Instability Mechanism with a Spring-Mass Analogue

Authors: Derek M. H. Hung, Eric G. Blackman, Kyle J. Caspary, Erik P. Gilson, Hantao Ji

Abstract: The Magnetorotational Instability (MRI) has long been considered a plausibly ubiquitous mechanism to destabilize otherwise stable Keplerian flows to support radially outward transport of angular momentum. Such an efficient transport process would allow fast accretion in astrophysical objects such as stars and black holes to release copious kinetic energy that powers many of the most luminous sourc… ▽ More The Magnetorotational Instability (MRI) has long been considered a plausibly ubiquitous mechanism to destabilize otherwise stable Keplerian flows to support radially outward transport of angular momentum. Such an efficient transport process would allow fast accretion in astrophysical objects such as stars and black holes to release copious kinetic energy that powers many of the most luminous sources in the universe. But the standard MRI under a purely vertical magnetic field has heretofore never been directly measured despite numerous efforts over more than a decade. Here we report an unambiguous laboratory demonstration of the spring-mass analogue to the standard MRI by comparing motion of a spring-tethered ball within different rotating flows. The experiment corroborates the theory: efficient outward angular momentum transport manifests only for cases with a weak spring in quasi-Keperian flow. Our experimental method accomplishes this in a new way, thereby connecting solid and fluid mechanics to plasma astrophysics. △ Less

Submitted 11 December, 2018; v1 submitted 10 January, 2018; originally announced January 2018.

Comments: 7 pages; 4 figures; accepted by Communications Physics

Mass, energy, and momentum capture from stellar winds by magnetized and unmagnetized planets: implications for atmospheric erosion and habitability

Authors: Eric G. Blackman, John A. Tarduno

Abstract: The extent to which a magnetosphere protects its planetary atmosphere from stellar wind ablation depends upon how well it prevents energy and momentum exchange with the atmosphere and how well it traps otherwise escaping plasma. We focus on the former, and provide a formalism for estimating approximate upper limits on mass, energy, and momentum capture, and use them to constrain loss rates. Our ap… ▽ More The extent to which a magnetosphere protects its planetary atmosphere from stellar wind ablation depends upon how well it prevents energy and momentum exchange with the atmosphere and how well it traps otherwise escaping plasma. We focus on the former, and provide a formalism for estimating approximate upper limits on mass, energy, and momentum capture, and use them to constrain loss rates. Our approach quantifies a competition between the local deflection of incoming plasma by a planetary magnetic field and the increase in area for solar wind mass capture provided by this magnetosphere when the wind and planetary fields incur magnetic reconnection. Solar wind capture can be larger for a magnetized planet versus an unmagnetized planet with small ionosphere, even if the rate of energy transfer is less. We find this to be the case throughout Earth's history since the lunar forming impact-- the likely start of the geodynamo. Focusing of incoming charged particles near the magnetic poles can increase the local energy flux for a strong enough planetary field and weak enough wind, and we find that this may be the case for Earth's future. The future protective role of Earth's magnetic field would then depend on its ability to trap outgoing plasma. This competition between increased collection area and reduced inflow speed from a magnetosphere is likely essential in determining the net protective role of planetary fields and its importance in sustaining habitability, and can help explain the current ion loss rates from Mars, Venus and Earth. △ Less

Submitted 6 October, 2018; v1 submitted 2 January, 2018; originally announced January 2018.

Comments: 10 pages; 2 figures; accepted and in press, MNRAS

Derivation and precision of mean field electrodynamics with mesoscale fluctuations

Authors: Hongzhe Zhou, Eric G. Blackman, Luke Chamandy

Abstract: Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two ques… ▽ More Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two questions: (I) what are the appropriate generalized equations of MFE in the presence of mesoscale fluctuations? (II) how precise are theoretical predictions from MFE? We address both by first deriving the equations of MFE for different types of averaging, along with mesoscale correction terms that depend on the ratio of averaging scale to variation scale of the mean. We then show that even if these terms are small, predictions of MFE can still have a significant precision error. This error has an intrinsic contribution from the dynamo input parameters and a filtering contribution from differences in the way observations and theory are projected through the measurement kernel.Minimizing the sum of these contributions can produce an optimal scale of averaging that makes the theory maximally precise.The precision error is important to quantify when comparing to observations because it quantifies the resolution of predictive power. We exemplify these principles for galactic dynamos, comment on broader implications, and identify possibilities for further work. △ Less

Submitted 21 March, 2019; v1 submitted 6 October, 2017; originally announced October 2017.

Comments: 35 pages, 10 figures. Published in Journal of Plasma Physics; Mar. 21 2019: Typos of a sign in (3.6) and indices of $ξくしー$ in the paragraphs containing (B7)-(B9) corrected

Journal ref: Journal of Plasma Physics, 84(3) (2018)

Sustained turbulence and magnetic energy in non-rotating shear flows

Authors: Farrukh Nauman, Eric G. Blackman

Abstract: From numerical simulations, we show that non-rotating magnetohydrodynamic shear flows are unstable to finite amplitude velocity perturbations and become turbulent, leading to the growth and sustenance of magnetic energy, including large scale fields. This supports the concept that sustained magnetic energy from turbulence is independent of the driving mechanism for large enough magnetic Reynolds n… ▽ More From numerical simulations, we show that non-rotating magnetohydrodynamic shear flows are unstable to finite amplitude velocity perturbations and become turbulent, leading to the growth and sustenance of magnetic energy, including large scale fields. This supports the concept that sustained magnetic energy from turbulence is independent of the driving mechanism for large enough magnetic Reynolds numbers. △ Less

Submitted 8 March, 2017; v1 submitted 12 January, 2017; originally announced January 2017.

Comments: 6 pages, 9 figures, accepted by PRE

Journal ref: Phys. Rev. E 95, 033202 (2017)

Evolution of the magnetorotational instability on initially tangled magnetic fields

Authors: Pallavi Bhat, Fatima Ebrahimi, Eric G. Blackman, Kandaswamy Subramanian

Abstract: The initial magnetic field of previous magnetorotational instability (MRI) simulations has always included a significant system-scale component, even if stochastic. However, it is of conceptual and practical interest to assess whether the MRI can grow when the initial field is turbulent. The ubiquitous presence of turbulent or random flows in astrophysical plasmas generically leads to a small-scal… ▽ More The initial magnetic field of previous magnetorotational instability (MRI) simulations has always included a significant system-scale component, even if stochastic. However, it is of conceptual and practical interest to assess whether the MRI can grow when the initial field is turbulent. The ubiquitous presence of turbulent or random flows in astrophysical plasmas generically leads to a small-scale dynamo (SSD), which would provide initial seed turbulent velocity and magnetic fields in the plasma that becomes an accretion disc. Can the MRI grow from these more realistic initial conditions? To address this we supply a standard shearing box with isotropically forced SSD generated magnetic and velocity fields as initial conditions, and remove the forcing. We find that if the initially supplied fields are too weak or too incoherent, they decay from the initial turbulent cascade faster than they can grow via the MRI. When the initially supplied fields are sufficient to allow MRI growth and sustenance, the saturated stresses, large-scale fields, and power spectra match those of the standard zero net flux MRI simulation with an initial large scale vertical field. △ Less

Submitted 11 October, 2017; v1 submitted 21 September, 2016; originally announced September 2016.

Comments: Published in MNRAS, 6 pages, 5 figures

Vector cylindrical harmonics for low-dimensional convection models

Authors: Douglas H. Kelley, Eric G. Blackman

Abstract: Approximate empirical models of thermal convection can allow us to identify the essential properties of the flow in simplified form, and to produce empirical estimates using only a few parameters. Such "low-dimensional" empirical models can be constructed systematically by writing numerical or experimental measurements as superpositions of a set of appropriate basis modes, a process known as Galer… ▽ More Approximate empirical models of thermal convection can allow us to identify the essential properties of the flow in simplified form, and to produce empirical estimates using only a few parameters. Such "low-dimensional" empirical models can be constructed systematically by writing numerical or experimental measurements as superpositions of a set of appropriate basis modes, a process known as Galerkin projection. For three-dimensional convection in a cylinder, those basis modes should be vector-valued, mutually orthogonal, and defined in cylindrical coordinates. Here we construct such a basis set and demonstrate that it has these desired properties and boundary conditions when the exact constraint of incompressibility is relaxed. We show its use for representing sample simulation data and point out its potential for low-dimensional convection models. △ Less

Submitted 6 October, 2016; v1 submitted 16 May, 2016; originally announced May 2016.

Comments: 9 pages; 8 figures, (revised to include fundamental improvements that result from eliminating the constraint of incompressibility), submitted to Physical Review Fluids

Large scale dynamo action precedes turbulence in shearing box simulations of the magnetorotational instability

Authors: Pallavi Bhat, Fatima Ebrahimi, Eric G. Blackman

Abstract: We study the dynamo generation (exponential growth) of large scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal ($x$-$y$) averaging, we demonstrate the presence of large scale fields when either horizontal or vertical ($y$-$z$) averaging is employed. By computing planar average… ▽ More We study the dynamo generation (exponential growth) of large scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal ($x$-$y$) averaging, we demonstrate the presence of large scale fields when either horizontal or vertical ($y$-$z$) averaging is employed. By computing planar averaged fields and power spectra, we find large scale dynamo action in the early MRI growth phase---a previously unidentified feature. Fast growing horizontal low modes and fiducial vertical modes over a narrow range of wave numbers amplify these planar averaged fields in the MRI growth phase, before turbulence sets in. The large scale field growth requires linear fluctuations but not nonlinear turbulence (as defined by mode-mode coupling) and grows as a direct global mode of the MRI. Only by vertical averaging, can it be shown that the growth of horizontal low wavenumber MRI modes directly feed-back to the initial vertical field providing a clue as to why the large scale vertical field sustains against turbulent diffusion in the saturation regime. We compute the terms in the planar averaged mean field equations to identify the individual contributions to large scale field growth for both vertical and horizontal averaging. The large scale fields obtained from such vertical averaging are found to compare well with global cylindrical simulations and quasilinear analytical analysis from a previous study by Ebrahimi \& Blackman. We discuss the potential implications of these new results for understanding large scale MRI dynamo saturation and turbulence. △ Less

Submitted 11 May, 2016; v1 submitted 9 May, 2016; originally announced May 2016.

Comments: 12 pages, 9 figues, submitted to MNRAS, comments are welcome

Shearing box simulations in the Rayleigh unstable regime

Authors: Farrukh Nauman, Eric G. Blackman

Abstract: We study the stability properties of Rayleigh unstable flows both in the purely hydrodynamic and magnetohydrodynamic (MHD) regimes for two different values of the shear $q=2.1, 4.2$ ($q = - d\lnΩおめが/ d\ln r$) and compare it with the Keplerian case $q=1.5$. We find that the $q>2$ regime is unstable both in the hydrodynamic and in the MHD limit (with an initially weak magnetic field). In this regime, t… ▽ More We study the stability properties of Rayleigh unstable flows both in the purely hydrodynamic and magnetohydrodynamic (MHD) regimes for two different values of the shear $q=2.1, 4.2$ ($q = - d\lnΩおめが/ d\ln r$) and compare it with the Keplerian case $q=1.5$. We find that the $q>2$ regime is unstable both in the hydrodynamic and in the MHD limit (with an initially weak magnetic field). In this regime, the velocity fluctuations dominate the magnetic fluctuations. In contrast, in the $q<2$ (magnetorotational instability (MRI)) regime the magnetic fluctuations dominate. This highlights two different paths to MHD turbulence implied by the two regimes, suggesting that in the $q>2$ regime the instability produces primarily velocity fluctuations that cause magnetic fluctuations, with the causality reversed for the $q<2$ MRI unstable regime. We also find that the magnetic field correlation is increasingly localized as the shear is increased in the Rayleigh unstable regime. In calculating the time evolution of spatial averages of different terms in the MHD equations, we find that the $q>2$ regime is dominated by terms which are nonlinear in the fluctuations, whereas for $q<2$, the linear terms play a more significant role. △ Less

Submitted 23 January, 2017; v1 submitted 16 July, 2015; originally announced July 2015.

Comments: accepted by MNRAS: 10 pages, 14 figures, 1 table; revised and expanded

Motivation and challenge to capture both large scale and local transport in next generation accretion theory

Authors: Eric G. Blackman, Farrukh Nauman

Abstract: Accretion disc theory is less developed than stellar evolution theory although a similarly mature phenomenological picture is ultimately desired. While the interplay of theory and numerical simulations has amplified community awareness of the role of magnetic fields in angular momentum transport, there remains a long term challenge to incorporate insight gained from simulations back into improving… ▽ More Accretion disc theory is less developed than stellar evolution theory although a similarly mature phenomenological picture is ultimately desired. While the interplay of theory and numerical simulations has amplified community awareness of the role of magnetic fields in angular momentum transport, there remains a long term challenge to incorporate insight gained from simulations back into improving practical models for comparison with observations. Here we emphasize the need to incorporate the role of non-local transport more precisely. To show where large scale transport would fit into the theoretical framework and how it is currently missing, we review why the wonderfully practical approach of Shakura-Sunyaev (1973,SS73) is necessarily a mean field theory, and one which does not include large scale transport. Observations of coronae and jets combined with the interpretation of results even from shearing box simulations of the magnetorotational instability (MRI) suggest that a significant fraction of disc transport is indeed non-local. We show that the Maxwell stresses in saturation are dominated by large scale contributions and the physics of MRI transport is not fully captured by a viscosity. We also clarify the standard physical interpretation of the MRi as it applies to shearing boxes. Computational limitations have so far focused most attention toward local simulations but the next generation of global simulations should help to inform improved mean field theories. Mean field accretion theory and mean field dynamo theory should in fact be unified into a single theory that predicts the time evolution of spectra and luminosity from separate disc, corona, and outflow contributions. Finally, we note that any mean field theory has a finite predictive precision that needs to be quantified when comparing the predictions to observations. △ Less

Submitted 4 July, 2017; v1 submitted 1 January, 2015; originally announced January 2015.

Comments: 31 Pages, 4 figures, published in Journal of Plasma Physics (for Special Issue commemorating Ya. B. Zeldovich); this version updated to correct mathematical typos in the first paragraph of section 4.4.1, identified in a post-publication corrigendum

Journal ref: J. Plasma Phys. 81 (2015) 395810505

Sensitivity of the Magnetorotational Instability to the shear parameter in stratified simulations

Authors: Farrukh Nauman, Eric G. Blackman

Abstract: The magnetorotational instability (MRI) is a shear instability and thus its sensitivity to the shear parameter $q = - d\lnΩおめが/d\ln r $ is of interest to investigate. Motivated by astrophysical disks, most (but not all) previous MRI studies have focused on the Keplerian value of $ q=1.5$. Using simulation with 8 vertical density scale heights, we contribute to the subset of studies addressing the the… ▽ More The magnetorotational instability (MRI) is a shear instability and thus its sensitivity to the shear parameter $q = - d\lnΩおめが/d\ln r $ is of interest to investigate. Motivated by astrophysical disks, most (but not all) previous MRI studies have focused on the Keplerian value of $ q=1.5$. Using simulation with 8 vertical density scale heights, we contribute to the subset of studies addressing the the effect of varying $q$ in stratified numerical simulations. We discuss why shearing boxes cannot easily be used to study $q>2$ and thus focus on $q<2$. As per previous simulations, which were either unstratified or stratified with a smaller vertical domain, we find that the $q$ dependence of stress for the stratified case is not linear, contrary to the Shakura-Sunyaev model. We find that the scaling agrees with \cite{1996MNRAS.281L..21A} who found it to be proportional to the shear to vorticity ratio $q/(2-q)$. We also find however, that the shape of the magnetic and kinetic energy spectra are relatively insensitive to $q$ and that the ratio of Maxwell stress to magnetic energy ratio also remains nearly independent of $q$. This is consistent with a theoretical argument in which the rate of amplification of the azimuthal field depends linearly on $q$ and the turbulent correlation time $τたう$ depends inversely on $q$. As such, we measure the correlation time of the turbulence and find that indeed it is inversely proportional to $q$. △ Less

Submitted 30 October, 2014; v1 submitted 8 September, 2014; originally announced September 2014.

Comments: 8 pages, 13 figures, 1 table, version accepted by MNRAS

Experimental demonstration of an inertial collimation mechanism in nested outflows

Authors: R. Yurchak, A. Ravasio, A. Pelka, S. Pikuz Jr., E. Falize, T. Vinci, M. Koenig, B. Loupias, A. Benuzzi-Mounaix, M. Fatenejad, P. Tzeferacos, D. Q. Lamb, E. G. Blackman

Abstract: Interaction between a central outflow and a surrounding wind is common in astrophysical sources powered by accretion. Understanding how the interaction might help to collimate the inner central outflow is of interest for assessing astrophysical jet formation paradigms. In this context, we studied the interaction between two nested supersonic plasma flows generated by focusing a long pulse high-ene… ▽ More Interaction between a central outflow and a surrounding wind is common in astrophysical sources powered by accretion. Understanding how the interaction might help to collimate the inner central outflow is of interest for assessing astrophysical jet formation paradigms. In this context, we studied the interaction between two nested supersonic plasma flows generated by focusing a long pulse high-energy laser beam onto a solid target. A nested geometry was created by shaping the energy distribution at the focal spot with a dedicated phase plate. Optical and X-ray diagnostics were used to study the interacting flows. Experimental results and numerical hydrodynamic simulations indeed show the formation of strongly collimated jets. Our work experimentally confirms the "shock-focused inertial confinement" mechanism proposed in previous theoretical astrophysics investigations. △ Less

Submitted 19 March, 2014; originally announced March 2014.

Comments: 6 pages, 5 figures, 1 table; accepted by Physical Review Letters

On characterizing nonlocality and anisotropy for the magnetorotational instability

Authors: Farrukh Nauman, Eric G. Blackman

Abstract: The extent to which angular momentum transport in accretion discs is primarily local or non-local and what determines this is an important avenue of study for understanding accretion engines. Taking a step along this path, we analyze simulations of the magnetorotational instability (MRI) by calculating energy and stress power spectra in stratified isothermal shearing box simulations in several new… ▽ More The extent to which angular momentum transport in accretion discs is primarily local or non-local and what determines this is an important avenue of study for understanding accretion engines. Taking a step along this path, we analyze simulations of the magnetorotational instability (MRI) by calculating energy and stress power spectra in stratified isothermal shearing box simulations in several new ways. We divide our boxes in two regions, disc and corona where the disc is the MRI unstable region and corona is the magnetically dominated region. We calculate the fractional power in different quantities, including magnetic energy and Maxwell stresses and find that they are dominated by contributions from the lowest wave numbers. This is even more dramatic for the corona than the disc, suggesting that transport in the corona region is dominated by larger structures than the disc. By calculating averaged power spectra in one direction of $k$ space at a time, we also show that the MRI turbulence is strongly anisotropic on large scales when analyzed by this method, but isotropic on small scales. Although the shearing box itself is meant to represent a local section of an accretion disc, the fact that the stress and energy are dominated by the largest scales highlights that the locality is not captured within the box. This helps to quantify the intuitive importance of global simulations for addressing the question of locality of transport, for which similar analyses can be performed. △ Less

Submitted 9 April, 2014; v1 submitted 17 March, 2014; originally announced March 2014.

Comments: 7 pages, 4 figures, 1 table, version accepted by MNRAS

Journal ref: MNRAS, Volume 441, Issue 3, p.1855-1860 07/2014

Ribbons characterize magnetohydrodynamic magnetic fields better than lines: a lesson from dynamo theory

Authors: Eric G. Blackman, Alexander Hubbard

Abstract: Blackman & Brandenburg argued that magnetic helicity conservation in dynamo theory can in principle be captured by diagrams of mean field dynamos when the magnetic fields are represented by ribbons or tubes, but not by lines. Here we present such a schematic ribbon diagram for the $αあるふぁ^2$ dynamo that tracks magnetic helicity and provides distinct scales of large scale magnetic helicity, small scale… ▽ More Blackman & Brandenburg argued that magnetic helicity conservation in dynamo theory can in principle be captured by diagrams of mean field dynamos when the magnetic fields are represented by ribbons or tubes, but not by lines. Here we present such a schematic ribbon diagram for the $αあるふぁ^2$ dynamo that tracks magnetic helicity and provides distinct scales of large scale magnetic helicity, small scale magnetic helicity, and kinetic helicity involved in the process. This also motivates our construction of a new "2.5 scale" minimalist generalization of the helicity-evolving equations for the αあるふぁ^2 dynamo that separately allows for these three distinct length scales while keeping only two dynamical equations. We solve these equations and, as in previous studies, find that the large scale field first grows at a rate independent of the magnetic Reynolds number R_M before quenching to an R_M dependent regime. But we also show that the larger the ratio of the wavenumber where the small scale current helicity resides to that of the forcing scale, the earlier the non-linear dynamo quenching occurs, and the weaker the large scale field is at the turnoff from linear growth. The harmony between the theory and the schematic diagram exemplifies a general lesson that magnetic fields in MHD are better visualized as two-dimensional ribbons (or pairs of lines) rather than single lines. △ Less

Submitted 9 May, 2014; v1 submitted 13 March, 2014; originally announced March 2014.

Comments: 11 pages, 10 figure panels, version accepted by MNRAS

Magnetic Helicity and Large Scale Magnetic Fields: A Primer

Authors: Eric G. Blackman

Abstract: Magnetic fields of laboratory, planetary, stellar, and galactic plasmas commonly exhibit significant order on large temporal or spatial scales compared to the otherwise random motions within the hosting system. Such ordered fields can be measured in the case of planets, stars, and galaxies, or inferred indirectly by the action of their dynamical influence, such as jets. Whether large scale fields… ▽ More Magnetic fields of laboratory, planetary, stellar, and galactic plasmas commonly exhibit significant order on large temporal or spatial scales compared to the otherwise random motions within the hosting system. Such ordered fields can be measured in the case of planets, stars, and galaxies, or inferred indirectly by the action of their dynamical influence, such as jets. Whether large scale fields are amplified in situ or a remnant from previous stages of an object's history is often debated for objects without a definitive magnetic activity cycle. Magnetic helicity, a measure of twist and linkage of magnetic field lines, is a unifying tool for understanding large scale field evolution for both mechanisms of origin. Its importance stems from its two basic properties: (1) magnetic helicity is typically better conserved than magnetic energy; and (2) the magnetic energy associated with a fixed amount of magnetic helicity is minimized when the system relaxes this helical structure to the largest scale available. Here I discuss how magnetic helicity has come to help us understand the saturation of and sustenance of large scale dynamos, the need for either local or global helicity fluxes to avoid dynamo quenching, and the associated observational consequences. I also discuss how magnetic helicity acts as a hindrance to turbulent diffusion of large scale fields, and thus a helper for fossil remnant large scale field origin models in some contexts. I briefly discuss the connection between large scale fields and accretion disk theory as well. The goal here is to provide a conceptual primer to help the reader efficiently penetrate the literature. △ Less

Submitted 17 March, 2014; v1 submitted 4 February, 2014; originally announced February 2014.

Comments: 32 pages, 18 figures (many with multiple panels); this version updated with a number of small corrections made in the proofs preparation stage; in press by Space Science Reviews and to eventually appear in associated ISSI book volume on "Multiscale Structure Formation and Dynamics in Cosmic Plasmas"; A. Balogh, Editor

Resilience of helical fields to turbulent diffusion II: direct numerical simulations

Authors: Pallavi Bhat, Eric G. Blackman, Kandaswamy Subramanian

Abstract: The recent study of Blackman and Subramanian (Paper I) indicates that large scale helical magnetic fields are resilient to turbulent diffusion in the sense that helical fields stronger than a critical value, decay on slow (~resistively mediated), rather than fast ($\sim$ turbulent) time scales. This gives more credence to potential fossil field origin models of the magnetic fields in stars, galaxi… ▽ More The recent study of Blackman and Subramanian (Paper I) indicates that large scale helical magnetic fields are resilient to turbulent diffusion in the sense that helical fields stronger than a critical value, decay on slow (~resistively mediated), rather than fast ($\sim$ turbulent) time scales. This gives more credence to potential fossil field origin models of the magnetic fields in stars, galaxies and compact objects. We analyze a suite of direct numerical simulations (DNS) of decaying large scale helical magnetic fields in the presence of non-helical turbulence to further study the physics of helical field decay. We study two separate cases: (1) the initial field is large enough to decay resistively, is tracked until it transitions to decay fast, and the critical large scale helical field at that transition is sought; (2) the case of Paper I, wherein there is a critical initial helical field strength below which the field undergoes fast decay right from the beginning. In case (1), both DNS and solutions of the two scale model (from Paper 1), reveal that the transition energy, $E_{c1}$, is independent of the scale of the turbulent forcing, within a small range of $\Rm$. We also find that the kinetic alpha, $αあるふぁ_K$, is subdominant to magnetic alpha, $αあるふぁ_M$, in the DNS, justifying an assumption in the two scale model. For case (2), we show exact solutions of two scale model in the limit of $ηいーた\rightarrow 0$ in fully helical case, leading to the transition energy, $E_{c2} = (k_1/\kf)^2 M_{eq}$, where $k_1$ and $\kf$ are the large scale and small turbulent forcing scale respectively and $M_{eq}$ is the equipartition magnetic energy. The DNS in this case agree qualitatively with the two scale model but the $R_M$ currently achievable, is too small to satisfy a condition $3/R_M << (k_1/k_f)^2$, necessary to robustly reveal the transition, $E_{c2}$ (Abridged). △ Less

Submitted 2 October, 2013; originally announced October 2013.

Comments: 15 pages, 16 figures, submitted to MNRAS

A study of large scale dynamo growth rates from numerical simulations and implications for mean field theories

Authors: Kiwan Park, Eric G. Blackman, Kandaswamy Subramanian

Abstract: Understanding large scale magnetic field growth in turbulent plasmas in the magnetohydrodynamic limit is a goal of magnetic dynamo theory. In particular, assessing how well large scale helical field growth and saturation in simulations matches that predicted by existing theories is important for progress. Using numerical simulations of isotropically forced turbulence without large scale shear with… ▽ More Understanding large scale magnetic field growth in turbulent plasmas in the magnetohydrodynamic limit is a goal of magnetic dynamo theory. In particular, assessing how well large scale helical field growth and saturation in simulations matches that predicted by existing theories is important for progress. Using numerical simulations of isotropically forced turbulence without large scale shear with the implications, we focus on several aspects of this comparison that have not been previously tested: (1) Leading mean field dynamo theories which break the field into large and small scales predict that large scale helical field growth rates are determined by the difference between kinetic helicity and current helicity with no dependence on the non-helical energy in small scale magnetic fields. Our simulations show that the growth rate of the large scale field from fully helical forcing is indeed unaffected by the presence or absence of small scale magnetic fields amplified in a precursor non-helical dynamo. However, because the precursor non helical dynamo in our simulations produced fields that were strongly sub-equipartition with respect to the kinetic energy, we cannot yet rule out the potential influence of stronger non- helical small scale fields. (2) We have identified two features in our simulations which cannot be explained by the most minimalist versions of two-scale mean field theory: (i) fully helical small scale forcing produces significant non-helical large scale magnetic energy and (ii) the saturation of the large scale field growth is time-delayed with respect to what minimalist theory predicts. We comment on desirable generalizations to the theory in this context and future desired work. △ Less

Submitted 13 May, 2013; v1 submitted 9 May, 2013; originally announced May 2013.

Comments: 8 pages, 9 figures. This will be published soon(Phys. Rev. E)

On deriving flux freezing in magnetohydrodynamics by direct differentiation

Authors: Eric G. Blackman

Abstract: The magnetic flux freezing theorem or Alfvén's theorem, is a basic principle of ideal magnetohydrodynamics (MHD), a commonly used approximation to describe the aspects of astrophysical and laboratory plasmas. The theorem states that the magnetic flux---the integral of magnetic field penetrating a surface---is conserved in time as that surface is distorted in in time by fluid motions. Pedagogues of… ▽ More The magnetic flux freezing theorem or Alfvén's theorem, is a basic principle of ideal magnetohydrodynamics (MHD), a commonly used approximation to describe the aspects of astrophysical and laboratory plasmas. The theorem states that the magnetic flux---the integral of magnetic field penetrating a surface---is conserved in time as that surface is distorted in in time by fluid motions. Pedagogues of MHD commonly derive flux freezing without showing how to take the material derivative of a general flux integral and/or assuming a vanishing field divergence from the outset. Here I avoid these shortcomings and derive flux freezing by direct differentiation, explicitly using a Jacobian to transform between the evolving field-penetrating surface at different times. The approach is instructive for its generality and helps elucidate the role of magnetic monopoles in breaking flux freezing. The paucity of appearances of this derivation in standard MHD texts suggests that its pedagogic value is underappreciated. △ Less

Submitted 1 February, 2013; v1 submitted 15 January, 2013; originally announced January 2013.

Comments: 7 pages, version accepted by European Journal of Physics (minor edits to previous version)

Particle-in-cell simulations of particle energization via shock drift acceleration from low Mach number quasi-perpendicular shocks in solar flares

Authors: Jaehong Park, Chuang Ren, Jared C. Workman, Eric G. Blackman

Abstract: Low Mach number, high beta fast mode shocks can occur in the magnetic reconnection outflows of solar flares. These shocks, which occur above flare loop tops, may provide the electron energization responsible for some of the observed hard X-rays and contemporaneous radio emission. Here we present new 2D particle-in-cell simulations of low Mach number/high beta quasi-perpendicular shocks. The simula… ▽ More Low Mach number, high beta fast mode shocks can occur in the magnetic reconnection outflows of solar flares. These shocks, which occur above flare loop tops, may provide the electron energization responsible for some of the observed hard X-rays and contemporaneous radio emission. Here we present new 2D particle-in-cell simulations of low Mach number/high beta quasi-perpendicular shocks. The simulations show that electrons above a certain energy threshold experience shock-drift-acceleration. The transition energy between the thermal and non-thermal spectrum and the spectral index from the simulations are consistent with some of the X-ray spectra from RHESSI in the energy regime of $E\lesssim 40\sim 100$ keV. Plasma instabilities associated with the shock structure such as the modified-two-stream and the electron whistler instabilities are identified using numerical solutions of the kinetic dispersion relations. We also show that the results from PIC simulations with reduced ion/electron mass ratio can be scaled to those with the realistic mass ratio. △ Less

Submitted 13 January, 2013; v1 submitted 20 October, 2012; originally announced October 2012.

Journal ref: ApJ, 765, 147, 2013

On the resilience of helical magnetic fields to turbulent diffusion and the astrophysical implications

Authors: Eric G. Blackman, Kandaswamy Subramanian

Abstract: The extent to which large scale magnetic fields are susceptible to turbulent diffusion is important for interpreting the need for in situ large scale dynamos in astrophysics and for observationally inferring field strengths compared to kinetic energy. By solving coupled equations for magnetic energy and magnetic helicity in a system initiated with isotropic turbulence and an arbitrarily helical la… ▽ More The extent to which large scale magnetic fields are susceptible to turbulent diffusion is important for interpreting the need for in situ large scale dynamos in astrophysics and for observationally inferring field strengths compared to kinetic energy. By solving coupled equations for magnetic energy and magnetic helicity in a system initiated with isotropic turbulence and an arbitrarily helical large scale field, we quantify the decay rate of the latter for a bounded or periodic system. The energy associated with the non-helical magnetic field rapidly decays by turbulent diffusion, but the decay rate of the helical component depends on whether the ratio of its magnetic energy to the turbulent kinetic energy exceeds a critical value given by M_{1,c} =(k_1/k_2)^2, where k_1 and k_2 are the wave numbers of the large and forcing scales. Turbulently diffusing helical fields to small scales while conserving magnetic helicity requires a rapid increase in total magnetic energy. As such, only when the helical fields are sub-critical can they so diffuse. When super-critical, the large scale helical field decays slowly, at a rate determined by microphysical dissipation even when macroscopic turbulence is present. Amplification of small scale magnetic helicity abates the turbulent diffusion. Two implications are that: (1) Standard arguments supporting the need for in situ large scale dynamos based on the otherwise rapid turbulent diffusion of large scale fields require re-thinking since only the non-helical field is so diffused in a closed system. Boundary terms could however provide potential pathways for rapid change of the large scale helical field. (2) Since M_{1,c} <<1 for k_1 << k_2, the presence of long-lived ordered large scale helical fields, as in extragalactic jets, does not guarantee that the magnetic field dominates the kinetic energy. △ Less

Submitted 16 November, 2015; v1 submitted 11 September, 2012; originally announced September 2012.

Comments: published in MNRAS (in this replacement, the missing .bbl file has been added)

Journal ref: Blackman, E.G., & Subramanian, K., 2013, MNRAS, 429, 1398

Toward reducing impact induced brain injury: Lessons from a computational study of army and football helmet pads

Authors: W. C. Moss, M. J. King, E. G. Blackman

Abstract: We use computational simulations to compare the impact response of different football and U.S. Army helmet pad materials. We conduct experiments to characterize the material response of different helmet pads. We simulate experimental helmet impact tests performed by the U.S. Army to validate our methods. We then simulate a cylindrical impactor striking different pads. The acceleration history of t… ▽ More We use computational simulations to compare the impact response of different football and U.S. Army helmet pad materials. We conduct experiments to characterize the material response of different helmet pads. We simulate experimental helmet impact tests performed by the U.S. Army to validate our methods. We then simulate a cylindrical impactor striking different pads. The acceleration history of the impactor is used to calculate the Head Injury Criterion for each pad. We conduct sensitivity studies exploring the effects of pad composition, geometry, and material stiffness. We find that: (1) The football pad materials do not outperform the currently used military pad material in militarily-relevant impact scenarios; (2) Optimal material properties for a pad depend on impact energy; and (3) Thicker pads perform better at all velocities. Our analysis suggests that by using larger helmet shells with correspondingly thicker pads, impact-induced traumatic brain injury may be significantly reduced. Keywords: helmet, pad, head injury, traumatic brain injury, head injury criterion, impact △ Less

Submitted 12 April, 2012; originally announced April 2012.

Comments: 24 pages, 12 figures, 1 table

Report number: LLNL-JRNL-490182-DRAFT

Journal ref: Computer Methods in Biomechanics and Biomedical Engineering (2012)

Particle-in-cell simulations of particle energization from low Mach number fast mode shocks

Authors: Jaehong Park, Jared C. Workman, Eric G. Blackman, Chuang Ren, Robert Siller

Abstract: Astrophysical shocks are often studied in the high Mach number limit but weakly compressive fast shocks can occur in magnetic reconnection outflows and are considered to be a site of particle energization in solar flares. Here we study the microphysics of such perpendicular, low Mach number collisionless shocks using two-dimensional particle-in-cell (PIC) simulations with a reduced ion/electron ma… ▽ More Astrophysical shocks are often studied in the high Mach number limit but weakly compressive fast shocks can occur in magnetic reconnection outflows and are considered to be a site of particle energization in solar flares. Here we study the microphysics of such perpendicular, low Mach number collisionless shocks using two-dimensional particle-in-cell (PIC) simulations with a reduced ion/electron mass ratio and employ a moving wall boundary method for initial generation of the shock. This moving wall method allows for more control of the shock speed, smaller simulation box sizes, and longer simulation times than the commonly used fixed wall, reflection method of shock formation. Our results, which are independent of the shock formation method, reveal the prevalence shock drift acceleration (SDA) of both electron and ions in a purely perpendicular shock with Alfvén Mach number $M_A=6.8$ and ratio of thermal to magnetic pressure $βべーた=8$. We determine the respective minimum energies required for electrons and ions to incur SDA. We derive a theoretical electron distribution via SDA that compares to the simulation results. We also show that a modified two-stream instability due to the incoming and reflecting ions in the shock transition region acts as the mechanism to generate collisionless plasma turbulence that sustains the shock. △ Less

Submitted 17 May, 2012; v1 submitted 29 February, 2012; originally announced March 2012.

Journal ref: Phys. Plasmas 19, 062904 (2012)

Simulations of a Magnetic Fluctuation Driven Large Scale Dynamo and Comparison with a Two-scale Model

Authors: Kiwan Park, E. G. Blackman

Abstract: Models of large scale (magnetohydrodynamic) dynamos (LSD) which couple large scale field growth to total magnetic helicity evolution best predict the saturation of LSDs seen in simulations. For the simplest so called "αあるふぁ2" LSDs in periodic boxes, the electromotive force driving LSD growth depends on the difference between the time-integrated kinetic and current helicity associated with fluctuations… ▽ More Models of large scale (magnetohydrodynamic) dynamos (LSD) which couple large scale field growth to total magnetic helicity evolution best predict the saturation of LSDs seen in simulations. For the simplest so called "αあるふぁ2" LSDs in periodic boxes, the electromotive force driving LSD growth depends on the difference between the time-integrated kinetic and current helicity associated with fluctuations. When the system is helically kinetically forced (KF), the growth of the large scale helical field is accompanied by growth of small scale magnetic (and current) helicity which ultimately quench the LSD. Here, using both simulations and theory, we study the complementary magnetically forced(MF) case in which the system is forced with an electric field that supplies magnetic helicity. For this MF case, the kinetic helicity becomes the back-reactor that saturates the LSD. Simulations of both MF and KF cases can be approximately modeled with the same equations of magnetic helicity evolution, but with complementary initial conditions. A key difference between KF and MF cases is that the helical large scale field in the MF case grows with the same sign of injected magnetic helicity, whereas the large and small scale magnetic helicities grow with opposite sign for the KF case. The MF case can arise even when the thermal pressure is approximately smaller than the magnetic pressure, and requires only that helical small scale magnetic fluctuations dominate helical velocity fluctuations in LSD driving. We suggest that LSDs in accretion discs and Babcock models of the solar dynamo are actually MF LSDs. △ Less

Submitted 10 January, 2012; v1 submitted 3 January, 2012; originally announced January 2012.

Comments: 12 pages, 34 figures

Modeling Poynting flux vs. kinetic-energy dominated jets

Authors: Martín Huarte-Espinosa, Adam Frank, Eric G. Blackman

Abstract: We present 3D-MHD AMR simulations of Poynting flux dominated (PFD) jets formed by injection of magnetic energy. We compare their evolution with a hydrodynamic jet which is formed by injecting kinetic energy with the same energy flux than the PFD jets. We predict characteristic emission distributions for each of these jets. Current-driven perturbations in PFD jets are amplified by both cooling and… ▽ More We present 3D-MHD AMR simulations of Poynting flux dominated (PFD) jets formed by injection of magnetic energy. We compare their evolution with a hydrodynamic jet which is formed by injecting kinetic energy with the same energy flux than the PFD jets. We predict characteristic emission distributions for each of these jets. Current-driven perturbations in PFD jets are amplified by both cooling and rotation for the regimes studied: Shocks and thermal pressure support are weakened by cooling, making the jets more susceptible to kinking. Rotation amplifies the toroidal magnetic field which also exacerbates the kink instability. △ Less

Submitted 3 November, 2011; originally announced November 2011.

Comments: 2 pages, 1 figure to appear in the proceedings of the "Magnetic Fields in the Universe III" meeting, Zakopane, PL

Comparison Between Turbulent Helical Dynamo Simulations and a Nonlinear Three-Scale Theory

Authors: Kiwan Park, E. G. Blackman

Abstract: Progress toward understanding principles of nonlinear growth and saturation of large scale magnetic fields has emerged from comparison of theoretical models that incorporate the evolution of magnetic helicity with numerical simulations for problems that are more idealized than expected in astrophysical circumstances, but still fully non-linear. We carry out a new comparison of this sort for the ma… ▽ More Progress toward understanding principles of nonlinear growth and saturation of large scale magnetic fields has emerged from comparison of theoretical models that incorporate the evolution of magnetic helicity with numerical simulations for problems that are more idealized than expected in astrophysical circumstances, but still fully non-linear. We carry out a new comparison of this sort for the magnetic field growth from forced isotropic helical turbulence in a periodic box. Previous comparisons be- tween analytic theory and simulations of this problem have shown that a two-scale model compares well with the simulations in agreeing that the driver of large scale field growth is the difference between kinetic and current helicities associated with the small scale field, and that the backreaction that slows the growth of the large scale field as the small scale current helicity grows. However, the two-scale model artificially re- stricts the small scale current helicity to reside at the same scale as the driving kinetic helicity. In addition, previous comparisons have focused on the late time saturation regime, and less on the early-time growth regime. Here we address these issues by comparing a three scale model to new simulations for both early and late time growth regimes. We find that the minimalist extension to three scales provides a better model for the field evolution than the two-scale model because the simulations show that the small scale current helicity does not reside at the same scale as that of the driving kinetic helicity. The simulations also show that the peak of the small scale current helicity migrates toward lower wave numbers as the growth evolves from the fast to saturated growth regimes. △ Less

Submitted 21 June, 2012; v1 submitted 30 August, 2011; originally announced August 2011.

Comments: This paper has been withdrawn by the author due to some typos

Not Much Helicity is Needed to Drive Large Scale Dynamos

Authors: Jonathan Pietarila Graham, Eric G. Blackman, Pablo D. Mininni, Annick Pouquet

Abstract: Understanding the in situ amplification of large scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity f_h,C has not been previously well quantified.… ▽ More Understanding the in situ amplification of large scale magnetic fields in turbulent astrophysical rotators has been a core subject of dynamo theory. When turbulent velocities are helical, large scale dynamos that substantially amplify fields on scales that exceed the turbulent forcing scale arise, but the minimum sufficient fractional kinetic helicity f_h,C has not been previously well quantified. Using direct numerical simulations for a simple helical dynamo, we show that f_h,C decreases as the ratio of forcing to large scale wave numbers k_F/k_min increases. From the condition that a large scale helical dynamo must overcome the backreaction from any non-helical field on the large scales, we develop a theory that can explain the simulations. For k_F/k_min>8 we find f_h,C< 3%, implying that very small helicity fractions strongly influence magnetic spectra for even moderate scale separation. △ Less

Submitted 21 June, 2012; v1 submitted 15 August, 2011; originally announced August 2011.

Comments: 5 pages, 4 figures

Journal ref: Physical Review E 85, 066406 (2012)

Fast Shocks From Magnetic Reconnection Outflows

Authors: Jared C. Workman, Eric G. Blackman, Chuang Ren

Abstract: Magnetic reconnection is commonly perceived to drive flow and particle acceleration in flares of solar, stellar, and astrophysical disk coronae but the relative roles of different acceleration mecha- nisms in a given reconnection environment are not well understood. We show via direct numerical simulations that reconnection outflows produce weak fast shocks, when conditions for fast recon- nection… ▽ More Magnetic reconnection is commonly perceived to drive flow and particle acceleration in flares of solar, stellar, and astrophysical disk coronae but the relative roles of different acceleration mecha- nisms in a given reconnection environment are not well understood. We show via direct numerical simulations that reconnection outflows produce weak fast shocks, when conditions for fast recon- nection are met and the outflows encounter an obstacle. The associated compression ratios lead to a Fermi acceleration particle spectrum that is significantly steeper than the strong fast shocks commonly studied, but consistent with the demands of solar flares. While this is not the only acceleration mechanism operating in a reconnection environment, it is plausibly a ubiquitous one. △ Less

Submitted 4 September, 2011; v1 submitted 30 April, 2011; originally announced May 2011.

Anderson and Stoner Published White Dwarf Mass Limits Before Chandrasekhar

Authors: Eric G. Blackman

Abstract: In their engaging recountals of Chandrasekhar's extraordinary career (Physics Today, vol 63, Issue 12, Dec 2010), neither Dyson nor Wali mention that Chandrasekhar was the third person not the first, to publish a white dwarf mass limit incorporating a relativistic treatment of degenerate electrons. As it has become a common misconception that Chandrasekhar was the first, a clarifying reminder on t… ▽ More In their engaging recountals of Chandrasekhar's extraordinary career (Physics Today, vol 63, Issue 12, Dec 2010), neither Dyson nor Wali mention that Chandrasekhar was the third person not the first, to publish a white dwarf mass limit incorporating a relativistic treatment of degenerate electrons. As it has become a common misconception that Chandrasekhar was the first, a clarifying reminder on this historical point is warranted. In short, the white dwarf mass limit widely attributed to Chandrasekhar (1931) should be the specific white dwarf mass limit calculated for a polytrope. The insight that a relativistic treatment of degeneracy leads to the existence of a white dwarf mass limit first appeared in papers of W. Anderson (1929) and E.C. Stoner (1930) for a uniform density star. Accordingly, Chandrasekhar (1931) cites Stoner (1930) and points out that the polytrope white dwarf mass limit is less than Stoner's uniform density white dwarf mass limit by about 20%. △ Less

Submitted 9 March, 2011; v1 submitted 7 March, 2011; originally announced March 2011.

Comments: 1 page letter, submitted to Physics Today (Dec 25, 2010), (in v2: added a reference)

Incorporating Human Body Mass in Standards of Helmet Impact Protection against Traumatic Brain Injury

Authors: Eric G. Blackman

Abstract: Impact induced traumatic brain injury (ITBI) describes brain injury from head impact not necessarily accompanied by skull fracture. For sufficiently abrupt head impact decelerations, ITBI results from brain tissue stress incurred as the brain crashes into the inside of the skull wall, displacing the surrounding cerebral spinal fluid (CSF). Proper helmet cushioning can damp the impact force and r… ▽ More Impact induced traumatic brain injury (ITBI) describes brain injury from head impact not necessarily accompanied by skull fracture. For sufficiently abrupt head impact decelerations, ITBI results from brain tissue stress incurred as the brain crashes into the inside of the skull wall, displacing the surrounding cerebral spinal fluid (CSF). Proper helmet cushioning can damp the impact force and reduce ITBI. But force is mass times acceleration and commonly used helmet blunt impact standards are based only on acceleration thresholds. Here I show how this implies that present standards overestimate the minimum acceleration onset for ITBI by implicitly assuming that the brain is mechanically decoupled from the body. I quantify how an arbitrary orientation of the body with respect to impact direction increases the effective mass that should be used in calculating the required damping force and injury threshold accelerations. I suggest a practical method to incorporate the body mass and impact angle into ITBI helmet standards and provide direction for further work. △ Less

Submitted 18 April, 2009; originally announced April 2009.

Comments: 12 Pages 1 figure, submitted to Journal of Biomechanics

Formation of Episodic Magnetically Driven Radiatively Cooled Plasma Jets in the Laboratory

Authors: F. Suzuki-Vidal, S. V. Lebedev, A. Ciardi, S. N. Bland, J. P. Chittenden, G. N. Hall, A. Harvey-Thompson, A. Marocchino, C. Ning, C. Stehle, A. Frank, E. G. Blackman, S. C. Bott, T. Ray

Abstract: We report on experiments in which magnetically driven radiatively cooled plasma jets were produced by a 1 MA, 250 ns current pulse on the MAGPIE pulsed power facility. The jets were driven by the pressure of a toroidal magnetic field in a ''magnetic tower'' jet configuration. This scenario is characterized by the formation of a magnetically collimated plasma jet on the axis of a magnetic ''bubbl… ▽ More We report on experiments in which magnetically driven radiatively cooled plasma jets were produced by a 1 MA, 250 ns current pulse on the MAGPIE pulsed power facility. The jets were driven by the pressure of a toroidal magnetic field in a ''magnetic tower'' jet configuration. This scenario is characterized by the formation of a magnetically collimated plasma jet on the axis of a magnetic ''bubble'', confined by the ambient medium. The use of a radial metallic foil instead of the radial wire arrays employed in our previous work allows for the generation of episodic magnetic tower outflows which emerge periodically on timescales of ~30 ns. The subsequent magnetic bubbles propagate with velocities reaching ~300 km/s and interact with previous eruptions leading to the formation of shocks. △ Less

Submitted 1 April, 2009; originally announced April 2009.

Comments: 6 pages, 5 figures. Accepted for publication in Astrophysics & Space Science

Journal ref: Astrophys.Space Sci.322:19-23,2009