我想绘制圆柱对称磁场中正电荷的运动。

我假设一个围绕 z 轴的圆柱体，磁场顺时针方向。B 场的大小为 6 T，与 z 轴的距离 R 为 3 m。带电粒子沿 z 轴正方向发射，能量为 2 MeV。

我不确定如何正确模拟这个 B 场。我正在考虑在圆柱坐标中创建 B 场，圆柱从 0 到 2pi：

theta = numpy.linspace(0, 2*numpy.pi, 360)
x = R*numpy.cos(theta)
y = R*numpy.sin(theta)

Bx = B0*(numpy.cos(numpy.arctan2(y,x)
By = B0*(-numpy.sin(numpy.arctan2(y,x)))
Bz = 0

然后创建一个向量B=[Bx, By, Bz]，我将根据该向量使用洛伦兹力计算时间跨度 t 的加速度。

但我想我正在绕圈子。是否有另一种方法来创建圆柱对称磁场？