对于在柱坐标中进行模拟的代码,我需要流出产生的速度 v(R) 和 v(z)。
我模拟了这样的流出,但使用球坐标(圆柱坐标现在不起作用),所以我在 r 方向vr(r,theta)
和 theta 方向都有速度vt(r,theta)
。在我的数据文件中,我计算了R
, z
, vR
,vz
但这些数据是二维数组,取决于r
和theta
- 所以我不知道我是否可以使用它们。
我以前从未做过插值,但我知道我可以使用我认为可以执行类似操作的工具更改数据的坐标:
R, z = np.linspace()
[Rm, zm ] = meshgrid(R,z)
def Sph2cyl()
r = np.sqrt(Rm**2 + zm**2)
theta = np.arctan2(Rm,zm)
[rc, theta_c] = Sph2cyl(Rm,zm)
vrc = Interp2(r,theta,vr,rc,theta_c) # I know name of the function is different in python.
我想用类似的东西我可以用圆柱坐标来表达我的径向球面速度,但它没有给我圆柱半径R
或的速度z
。
会像这样的转变:
vR = vr * R/rc
vz = vr * z/rc
vrc
插值和足够了rc
吗?