(mu/2pi)log(b/a) for a=coreradius, b=shellradius