The Cosmic Engine: How Matter Conquers a Black Hole's Magnetic Shield

The Core Discovery

By modeling a black hole of mass $M$ with a high spin parameter ($a = 0.99$), 2D HARM numerical simulations demonstrate that accreting matter effectively crushes the Meissner expulsion. It drags magnetic lines of force across the event horizon, fundamentally reshaping the environment of the abyss.

Why It Matters: The Universe's Plumbing

This discovery explains the mechanics of the universe’s most violent engines. It shows black holes don't just swallow matter; they act as cosmic looms. By weaving magnetic fields, they create "magnetically ejected disks"—intermittent equatorial outflows that spit material back into space.

The Arena

The simulated domain stretched from $0.65 R_g$ out to a vast $10^3 R_g$.

• Initialization: A non-magnetized Bondi inflow of plasma.
• The Challenge: An axially symmetric Wald uniform magnetic field acting as a wall against the infalling matter.

The Tug-of-War

In cases of strong magnetization (pressure ratio $\beta \approx 1$), a critical shift occurred:

• The mass accretion rate ($\dot{M}$) didn't just stabilize—it was suppressed.
• This suppression was driven by the sudden birth of equatorial outflows.

The Four-Phase Evolution

The team observed a distinct evolution, peaking at saturation around $t = 180$ (in geometrical units).

• Key Outcome: The once-uniform magnetic field was forced into a "split-monopole" topology.
• Result: This creates a high-pressure current sheet at the equator where magnetic reconnection—the snapping and realigning of field lines—fires "plasmoids" (high-energy blobs of plasma) radially outward.

Technical Foundations

To maintain numerical stability in vacuum-dominated regions, the simulation relied on:

• A grid resolution of $600 \times 512$ points.
• A maintained density floor for scenarios where $\beta < 10^{-4}$.

Model Limitations

1. Dimensional Constraint: The use of a 2D axially symmetric scheme means certain three-dimensional instabilities remain unexplored.
2. Physics Omission: The model does not yet account for radiative cooling.
3. Driver of Events: The observed reconnection events are currently driven by numerical resistivity rather than a modeled physical one.