BackLinks

Getting this plot using the BE_Field in AstroBEAR

As a measure of magnetic field strength in the pdf, I am going to use the magnetic pressure.

Magnetic pressure is an energy density, and in free space in cgs, is given by:

P_B = \frac{B^2}{8 \pi}

In the code, there is a field called 'BE_field' which is the magnetic energy (or pressure) field, given by:

BE\_field = \frac{1}{2} (B_x^2 + B_y^2 + B_z^2)

The B's in the above are in computational units (cu), implying -

BE\_field = \frac{1}{2} (||\vec{B}(cu)||^2) \equiv \frac{1}{2} B^2 (cu)

= \frac{1}{2} \frac{B^2(cgs)}{Bscale^2}

where

Bscale = \sqrt{4\pi *pscale}

Putting this into BE_field gives,

BE\_field = \frac{B^2(8\pi)^{-1}}{pscale} = \frac{P_B}{pscale} = P_B(cu)

aka - scaled magnetic pressure.

Therefore - using BE_field as is for my pdf - means I will be plotting scaled magnetic pressure vs. scaled number density.

What is pscale?

pscale = nscale(cm^{-3})*K_B(cgs)*tempscale(K)

For my problem set-up,

nscale=tempscale=1

which let's me do a cute trick -

P_B(cu) = \frac{P_B(cgs)}{nscale*K_B*tempscale}

which implies

P_B(cu)*nscale*tempscale = P_B(cu)[\frac{K}{cm^3}] = \frac{P_{B}(cgs)}{K_B}

That is, the LHS, which are the values I will be plotting, are in units of K/cm³, and therefore instead of calling it something like 'magnetic pressure in computational units', I can call it:

P_B K_B ^{-1} (\frac{K}{cm^3})

which is still a scaled magnetic pressure — but perhaps a more physically motivated one for readers..

This of course is plotted against number density, which again given nscale = 1, is already in units of cm^-3.

Code