athenakit.physics

Physics modules for computing derived quantities in Newtonian and general-relativistic regimes.

GRMHD

Derived variables and helper functions for general-relativistic MHD in stationary (Kerr) spacetimes.

athenakit.physics.grmhd.cks_geometry(x, y, z, a)

athenakit.physics.grmhd.cks_to_sks_tens_con(axx, axy, axz, ayx, ayy, ayz, azx, azy, azz, x, y, z, a)

athenakit.physics.grmhd.cks_to_sks_vec_con(ax, ay, az, x, y, z, a)

athenakit.physics.grmhd.cks_to_sks_vec_cov(a_x, a_y, a_z, x, y, z, a)

athenakit.physics.grmhd.functions(a)

athenakit.physics.grmhd.lower_vector(at, ax, ay, az, g_tt, g_tx, g_ty, g_tz, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz)

athenakit.physics.grmhd.norm_to_coord(uut, uux, uuy, uuz, alpha, g_tx, g_ty, g_tz)

athenakit.physics.grmhd.normal_lorentz(uux, uuy, uuz, g_xx, g_xy, g_xz, g_yy, g_yz, g_zz)

athenakit.physics.grmhd.r_horizon(a)

athenakit.physics.grmhd.three_to_four_field(bbx, bby, bbz, ut, ux, uy, uz, u_x, u_y, u_z)

athenakit.physics.grmhd.variables(f, a)

Calculate all the variables used in the GRMHD simulation

Parameters:

• f (function) – A function that takes a string as input and returns the value of the variable

• a (float) – The spin parameter of the black hole

Returns:

v – A dictionary containing all the variables

Return type:

dict

Metric

Kerr metric utilities (Boyer-Lindquist and Kerr-Schild coordinates).

athenakit.physics.metric.kerr_schild_metric_and_inverse(x, y, z, a)

athenakit.physics.metric.kerr_schild_radius(x, y, z, a)

SNR

Analytic solutions for supernova remnant evolution.

class athenakit.physics.snr.SNR_ABC(A, B=1.0, C=1.0, gamma=1.6666666666666667, E0=1e+51, n0=1.0, mu=0.618, psi_x=0.2)

Bases: object

__init__(A, B=1.0, C=1.0, gamma=1.6666666666666667, E0=1e+51, n0=1.0, mu=0.618, psi_x=0.2)

func(x, Y)

class athenakit.physics.snr.SNR_SedovTaylor(gamma=1.6666666666666667, E0=1e+51, n0=1.0, mu=0.618)

Bases: object

__call__(xs, ts)

Call self as a function.

__init__(gamma=1.6666666666666667, E0=1e+51, n0=1.0, mu=0.618)

func(Y, x)

post_solve(ts)

set_K()

solve(xs)

class athenakit.physics.snr.SNR_WhiteLong(tau1, mucl, gamma=1.6666666666666667, E0=1e+51, n0=1.0, mu=0.618)

Bases: object

__call__(xs, ts)

Call self as a function.

__init__(tau1, mucl, gamma=1.6666666666666667, E0=1e+51, n0=1.0, mu=0.618)

func(x, Y)

post_solve(ts)

set_K()

solve(xs)

class athenakit.physics.snr.SNR_evo(n=1, M=3, E=1, mu=0.618, gamma=1.6666666666666667, config=True)

Bases: object

__init__(n=1, M=3, E=1, mu=0.618, gamma=1.6666666666666667, config=True)

config(t=array([1.00000000e-04, 1.03128316e-04, 1.06354496e-04, 1.09681601e-04, 1.13112788e-04, 1.16651313e-04, 1.20300535e-04, 1.24063916e-04, 1.27945027e-04, 1.31947552e-04, 1.36075289e-04, 1.40332154e-04, 1.44722187e-04, 1.49249555e-04, 1.53918552e-04, 1.58733611e-04, 1.63699300e-04, 1.68820332e-04, 1.74101565e-04, 1.79548012e-04, 1.85164842e-04, 1.90957383e-04, 1.96931134e-04, 2.03091762e-04, 2.09445114e-04, 2.15997219e-04, 2.22754295e-04, 2.29722754e-04, 2.36909207e-04, 2.44320476e-04, 2.51963593e-04, 2.59845810e-04, 2.67974609e-04, 2.76357701e-04, 2.85003044e-04, 2.93918840e-04, 3.03113550e-04, 3.12595900e-04, 3.22374888e-04, 3.32459793e-04, 3.42860186e-04, 3.53585937e-04, 3.64647222e-04, 3.76054540e-04, 3.87818715e-04, 3.99950910e-04, 4.12462638e-04, 4.25365773e-04, 4.38672559e-04, 4.52395623e-04, 4.66547988e-04, 4.81143084e-04, 4.96194760e-04, 5.11717301e-04, 5.27725435e-04, 5.44234355e-04, 5.61259726e-04, 5.78817704e-04, 5.96924951e-04, 6.15598650e-04, 6.34856522e-04, 6.54716840e-04, 6.75198452e-04, 6.96320794e-04, 7.18103909e-04, 7.40568469e-04, 7.63735792e-04, 7.87627861e-04, 8.12267350e-04, 8.37677640e-04, 8.63882844e-04, 8.90907830e-04, 9.18778243e-04, 9.47520530e-04, 9.77161967e-04, 1.00773068e-03, 1.03925568e-03, 1.07176689e-03, 1.10529514e-03, 1.13987227e-03, 1.17553107e-03, 1.21230540e-03, 1.25023015e-03, 1.28934130e-03, 1.32967597e-03, 1.37127243e-03, 1.41417017e-03, 1.45840988e-03, 1.50403355e-03, 1.55108448e-03, 1.59960730e-03, 1.64964807e-03, 1.70125428e-03, 1.75447489e-03, 1.80936041e-03, 1.86596292e-03, 1.92433614e-03, 1.98453546e-03, 2.04661800e-03, 2.11064268e-03, 2.17667025e-03, 2.24476338e-03, 2.31498667e-03, 2.38740677e-03, 2.46209240e-03, 2.53911443e-03, 2.61854596e-03, 2.70046235e-03, 2.78494135e-03, 2.87206312e-03, 2.96191033e-03, 3.05456825e-03, 3.15012480e-03, 3.24867066e-03, 3.35029934e-03, 3.45510729e-03, 3.56319397e-03, 3.67466194e-03, 3.78961698e-03, 3.90816818e-03, 4.03042803e-03, 4.15651256e-03, 4.28654141e-03, 4.42063797e-03, 4.55892950e-03, 4.70154722e-03, 4.84862648e-03, 5.00030684e-03, 5.15673224e-03, 5.31805113e-03, 5.48441658e-03, 5.65598646e-03, 5.83292359e-03, 6.01539588e-03, 6.20357648e-03, 6.39764396e-03, 6.59778248e-03, 6.80418197e-03, 7.01703829e-03, 7.23655342e-03, 7.46293569e-03, 7.69639990e-03, 7.93716762e-03, 8.18546731e-03, 8.44153460e-03, 8.70561248e-03, 8.97795155e-03, 9.25881025e-03, 9.54845510e-03, 9.84716096e-03, 1.01552113e-02, 1.04728984e-02, 1.08005237e-02, 1.11383983e-02, 1.14868426e-02, 1.18461873e-02, 1.22167735e-02, 1.25989528e-02, 1.29930878e-02, 1.33995527e-02, 1.38187331e-02, 1.42510267e-02, 1.46968439e-02, 1.51566076e-02, 1.56307542e-02, 1.61197336e-02, 1.66240098e-02, 1.71440614e-02, 1.76803818e-02, 1.82334800e-02, 1.88038809e-02, 1.93921257e-02, 1.99987727e-02, 2.06243975e-02, 2.12695939e-02, 2.19349740e-02, 2.26211693e-02, 2.33288310e-02, 2.40586305e-02, 2.48112606e-02, 2.55874352e-02, 2.63878910e-02, 2.72133877e-02, 2.80647085e-02, 2.89426612e-02, 2.98480792e-02, 3.07818214e-02, 3.17447741e-02, 3.27378510e-02, 3.37619944e-02, 3.48181763e-02, 3.59073989e-02, 3.70306958e-02, 3.81891330e-02, 3.93838098e-02, 4.06158599e-02, 4.18864524e-02, 4.31967930e-02, 4.45481252e-02, 4.59417313e-02, 4.73789339e-02, 4.88610967e-02, 5.03896263e-02, 5.19659730e-02, 5.35916329e-02, 5.52681486e-02, 5.69971109e-02, 5.87801607e-02, 6.06189899e-02, 6.25153435e-02, 6.44710211e-02, 6.64878784e-02, 6.85678294e-02, 7.07128478e-02, 7.29249692e-02, 7.52062927e-02, 7.75589832e-02, 7.99852734e-02, 8.24874655e-02, 8.50679342e-02, 8.77291280e-02, 9.04735724e-02, 9.33038717e-02, 9.62227117e-02, 9.92328623e-02, 1.02337180e-01, 1.05538610e-01, 1.08840192e-01, 1.12245057e-01, 1.15756437e-01, 1.19377664e-01, 1.23112175e-01, 1.26963513e-01, 1.30935333e-01, 1.35031404e-01, 1.39255613e-01, 1.43611969e-01, 1.48104605e-01, 1.52737785e-01, 1.57515906e-01, 1.62443501e-01, 1.67525247e-01, 1.72765966e-01, 1.78170632e-01, 1.83744372e-01, 1.89492477e-01, 1.95420401e-01, 2.01533769e-01, 2.07838382e-01, 2.14340223e-01, 2.21045463e-01, 2.27960464e-01, 2.35091788e-01, 2.42446202e-01, 2.50030685e-01, 2.57852435e-01, 2.65918874e-01, 2.74237657e-01, 2.82816678e-01, 2.91664078e-01, 3.00788252e-01, 3.10197859e-01, 3.19901828e-01, 3.29909369e-01, 3.40229977e-01, 3.50873446e-01, 3.61849876e-01, 3.73169684e-01, 3.84843611e-01, 3.96882735e-01, 4.09298482e-01, 4.22102632e-01, 4.35307336e-01, 4.48925126e-01, 4.62968923e-01, 4.77452054e-01, 4.92388263e-01, 5.07791724e-01, 5.23677054e-01, 5.40059328e-01, 5.56954091e-01, 5.74377375e-01, 5.92345715e-01, 6.10876161e-01, 6.29986298e-01, 6.49694260e-01, 6.70018750e-01, 6.90979055e-01, 7.12595063e-01, 7.34887289e-01, 7.57876886e-01, 7.81585671e-01, 8.06036141e-01, 8.31251499e-01, 8.57255673e-01, 8.84073340e-01, 9.11729948e-01, 9.40251743e-01, 9.69665789e-01, 1.00000000e+00]))

momr(t)

pres(t)

r(t)

v(t)

class athenakit.physics.snr.SNR_mpe(tau1, mucl, B=1.0, C=1.0, gamma=1.6666666666666667, E0=1e+51, n0=1.0, mu=0.618, psi_x=0.2)

Bases: object

__call__(xs, ts)

Call self as a function.

__init__(tau1, mucl, B=1.0, C=1.0, gamma=1.6666666666666667, E0=1e+51, n0=1.0, mu=0.618, psi_x=0.2)

call(xs, ts)

func(x, Y)

post_solve(ts)

set_K()

solve(xs)

class athenakit.physics.snr.SedovTaylor(E, rho, gamma=1.6666666666666667)

Bases: object

__init__(E, rho, gamma=1.6666666666666667)

r_s(t)

v_s(t)