obliqueshock_nif

Script calculates shock crossing values from Rankine-Hugoniot relations in the normal incidence frame (NIF).

obliqueshock_nif.rankine(Tu, rhou, V, B)

Call obliqueshock_nif.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]

Note: Give the plasma velocity relative to the shock speed, which is set to zero in the script.

Returns:

The shock compression ratio X: Plasma compression ratio XB: Magnetic field compression ratio

Example:

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

-> 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),