Contact
CoCalc Logo Icon
StoreFeaturesDocsShareSupport News Sign UpSign In
| Download
Project: BHLectures
Views: 612
Kernel: SageMath 8.0.beta0

Meridional plots in Kerr spacetime for a/m=0.99a/m=0.99

%display latex

O'Neill exponential coordinates:

x,y = var('x y') r = 1/2*ln(x^2 + y^2) costh2 = y^2/(x^2+y^2) sinth2 = x^2/(x^2+y^2)

Function ff defining the ergoregion by f<0f<0:

a = 0.99 f = r^2 -2*r + a^2*costh2 f
14log(x2+y2)2+0.980100000000000y2x2+y2log(x2+y2)\renewcommand{\Bold}[1]{\mathbf{#1}}\frac{1}{4} \, \log\left(x^{2} + y^{2}\right)^{2} + \frac{0.980100000000000 \, y^{2}}{x^{2} + y^{2}} - \log\left(x^{2} + y^{2}\right)
ergo = region_plot(f < 0, (x,-8, 8), (y, -5, 5), incol='lightgray', bordercol='grey', axes_labels=[r'$\mathrm{e}^{r/m}\sin\theta$', r'$\mathrm{e}^{r/m}\cos\theta$']) ergo += text(r'$\mathscr{E}^+$', (1.05*e^2, 0.5*e), color='grey', fontsize=20) ergo += text(r'$\mathscr{E}^-$', (1.5, 0.8), color='grey', fontsize=20)

Various remarkable surfaces:

Rp = exp(1 + sqrt(1-a^2)) Rm = exp(1 - sqrt(1-a^2)) Hp = circle((0,0), Rp, color='black', thickness=2) + \ text(r'$\mathscr{H}$', (0.84*Rp, 0.77*Rp), color='black', fontsize=20) Hm = circle((0,0), Rm, color='green', thickness=2) + \ text(r'$\mathscr{H}_{\rm in}$', (0.75*Rm, 0.9*Rm), color='green', fontsize=20) R0 = circle((0,0), 1, color='darkorange', linestyle='dotted', thickness=3) + \ text(r'$r\!=\!0$', (1.1,-1.), color='darkorange', fontsize=16) sing = circle((1,0), 0.1, color='red', fill=True) + \ circle((-1,0), 0.1, color='red', fill=True) rminf = circle((0,0), 0.1, edgecolor='black', facecolor='white', fill=True) region_label = text(r'${\rm I}$', (-2.3*Rp, Rp), fontsize=20) + \ text(r'${\rm II}$', (-0.72*Rp, 0.5*Rp), fontsize=20) + \ text(r'${\rm III}$', (-0.4*Rm, 0.5*Rm), fontsize=20) graph = ergo + Hp + Hm + R0 + sing + rminf + region_label

Carter time machine:

ft = (r^2+a^2)*(r^2+a^2*costh2) + 2*a^2*r*sinth2
tmachine = region_plot(ft < 0, (x,-2, 2), (y, -2, 2), incol='yellow', bordercol='gold') graph += tmachine show(graph, aspect_ratio=1)
Image in a Jupyter notebook
graph.save("ker_ergo_a99.pdf", aspect_ratio=1)