PlasmaPy / hack-week

2022 Hack Week material...slides,notebooks, and more
BSD 2-Clause "Simplified" License
9 stars 41 forks source link

Implement the relaxation rates (collision frequencies) in the NRL formulary #23

Closed pheuer closed 2 years ago

pheuer commented 2 years ago

The NRL formulary includes relaxation rates/collision frequencies (pg. 31 in this version) for several collision processes. These are expressed as functions of an integral Phi, which itself is a function of the ratio between the bulk flow energy (1/2 mv^2) and the thermal energy (kT) of the plasma: x.

The next page (pg. 32) expresses these four frequencies for different combinations of ions and electrons in both the fast (x>>1) and slow (x<<1) limits (in which the integral of x is approximately constant). However, numerically it is easy enough to just evaluate the integral Phi and get one expression for both regimes.

It would be great to have the integral Phi and the four expressions for the relaxation rates $\nu$ in the PlasmaPy collisions package! In the fast beam limit these are particularly useful for estimating collisionality in shock experiments.

image

JaydenR2305 commented 2 years ago

Hi Peter. I'll give this a try