Scientific interests


How to make astrophysical jets?

Astrophysical jets are defined as highly collimated beams of matter moving with high velocity. It is now generally accepted that magnetic fields are responsible for both acceleration and collimation of the flow.

As a general property observed in all astrophysical jet sources, is the signature for the existence of an accretion disk. This holds for jets sources of different scale of energy and spatial extent:
Active Galactic Nuclei,
Micro Quasars
(jets from compact stars), and
Protostellar Objects.

We apply numerical MHD simulations to calculate the time-dependent structure of the disk-jet magnetosphere, and the dynamics of the plasma in the formation region of magnetized jets.

With jet formation we denote the accelaration and collimation of a disk wind into a jet flow, with jet launching we denote the process which lifts the accreting material into the direction vertical to the disk, producing a disk wind.

Below we present (in chronological order) example results of our time-dependent MHD jet formation simulations using the ZEUS-3D code and the PLUTO code. These simulations follow the general framework of Blandford & Payne (1982) and Ouyed & Pudritz (1997). For details we refer to our published papers.

Model scenario of a magnetised young star-disk system (Fendt 1994, Fendt & Camenzind 1995).


MHD simulations: flaring disk-jet magnetospheres causing time-variable mass fluxes (2009)

This paper investigates time-dependent evolution of a superposed disk magnetic field with a stellar dipolar magnetosphere. (see Fendt, 2009, 692, 346 pdf-file).
Field distributions with aligned and anti-aligned magnetic axis of disk and stellar field are investigated. We use our version of ZEUS-3D extended for physical magnetic diffusivity, thus being able to follow reconnection events. The X-point of where anit-aligned disk field and stellar field meet along the disk surface gives rise to strong reconnection events with strong flares emerging across the jet magnetosphere. The flares evolve similar to coronal mass ejections changing the mass flux and velocity profile across the jet and also changing the mass flux number value by a factor of two to four with in short time. Our hypothesis is is that these flares initiate the observed jet knots. For our magnetic diffusivity model the time scales somehow fit.

The following gif-animations show the long-term evolution of the flow (use xanim).

Large frames, initial time steps (~3MB) Medium size frames, long run (~15MB)

Fewer frames, long run (~34MB) Large frames, long run (~35MB)


MHD simulations: disk magnetization profile and jet collimation (2006)

This paper investigates the interrelation between the profile of the accretion disk magnetic field profile (and the disk wind density profile) and the degree of jet collimation. (see Fendt, 2006, ApJ 651, 272, pdf-file).
Disk density profile and magnetic flux profile are prescribed as power law profiles. Both variables can be unifoed in the disk magnetisation profile. We find a unique relation between the magnetization profile power law index and the jet degeree of collimation (measured by the ratio of axial to lateral mass fluxes). Flat magnetic field / magnetization profiles collimate the outflow to a higher degree. Steep profiles such as a dipolar-like field distribution or a so-called X-wind do not produce highly collimated jets.


MHD simulations: magnetic diffusivity and jet collimation (2002)

In another paper we investigated the relation between the magnitude of jet magnetic diffusivity and the degree of jet collimation. (see Fendt & Cemeljic, 2002, A&A 395, 1045, pdf-file).
Diffusive jets are generally less collimated and there seems to be a critical value for the diffusivity (eta) above which the "jet" remains uncollimated. The jet velocities become faster with increasing diffusivity but the bow shock propagates slower (see density contours below with normalized eta = 0, 0.01, 0.1)


MHD simulations: dipolar magnetosphere (1999)

Using the ZEUS 3D code in the axisymmetry option the evolution of a stellar dipolar-type magnetic field interacting with an accretion disk is calculated. The boundary condition is an inflow from a Keplerian disk (as in Ouyed and Pudritz 1997). In the movies the accretion disk is at the lower boundary. The initial magnetic field is locked in a rigidly rotating stellar surface and the disk. The size of the domain is 20x20 inner disk radii. Shown are density (colors) and poloidal field lines (black).

In the first simulation, the star is at rest. 100 Keplerian periods of the inner disk are calculated. A bubble forms disrupting the dipole. A disk wind accelerates and slowly collimates, indicating a possible final stationary state (Fendt & Elstner, 1999, A&A 349, L61, pdf-file ).

Dipolar field, star at rest (~2MB)

In the second example the star rotates with a corotation radius at the inner disk radius. The initial field structure is a force-free dipole quenched along the equatorial plane. A two-component outflow (disk wind and stellar wind) is formed which is uncollimated. More than 2500 Keplerian periods of the inner disk are calculated in order to obtain a quasi-stationary final state (run S2). A larger stellar wind mass flow rate stabilizes the flow along the axis (run L5) (see Fendt & Elstner, 2000, A&A 363, 208 pdf-file). The following gif-animations show the long-term evolution of the flow (use xanim).

Small frames, long run (~11MB) Small frames, final time steps (~3MB)
Large frames, final time steps (~9MB) Medium size frames, long run (~25MB)
Large frames, long run (~35MB)


MHD jet formation simulations (1999)

Using the ZEUS-3D code in the axisymmetry option we calculated the jet formation from an accretion disk into a hydrostatic disk corona. The model of Ouyed and Pudritz (1997) is applied and used as reference for further simulations, e.g. of dipolar-type magnetic fields interacting with an accretion disk (below).

The jet runs from left (accretion disk) to the right (hydrostatic corona). Shown are density (colors) and poloidal field lines (black). The jet ACCELERATES and COLLIMATES.

Movie of the jet simulation (~1MB)

©2010 Christian Fendt - last modified July 2010