obliqueshock

Calculates shock crossing values from Rankine-Hugoniot relations in the deHoffmann-Teller (dHT) frame.

obliqueshock.rankine(Tu, rhou, V, B, n, Vsh)

Call obliqueshock.rankine(Tu, rhou, V, B, n, Vsh) to compute the shock crossing values

Inputs:

Tu: upstream proton temperature [K] rhou: upstream proton number density [1/m3] V: 3-element upstream proton inflow velocity vector [m/s] B: 3-element upstream magnetic field vector [T] n: 3-element shock normal vector Vsh: 3-element shock velocity vector [m/s]

Returns:

The shock compression ratio X: Compression ratio from scipy.optimize X2: Compression ratio from the Newton method

Example:

obliqueshock.rankine(5e5, 1.0e6, [-750e3,0,0], [3.5355e-9,0,-3.5355e-9], [1,0,0], 0)

-> Computes the shock crossing for Tu = 500 kK, rhou = 1/cc, V = [-750,0,0]km/s (inflow plasma speed is 750 km/s in -X), B = [3.5355e,0,-3.5355e]nT (upstream magnetic field is 5 nT at 45 degree angle), n = [1,0,0], Vsh = 0 (shock front points in +X direction and is stationary in input frame)