CFD中文网

是否是这样？应力张量我可以理解。我的疑问是速度梯度在笛卡尔坐标系下变换之后直接就变成在极坐标系中了速度梯度了？

具体问题是求粘弹流绕圆柱绕流，需设置应力张量在边界的值，在圆柱壁面应力张量的公式是极坐标系下，需要通过三步计算得到需要的结果：
1、笛卡尔坐标系下的速度梯度（张量）变换到极坐标系下的速度梯度。
2、通过（公式18）计算应力张量。
3、应力张量再变换到笛卡尔坐标系下。
请大家指点我应该用什么知识，应该有直接的变换公式吧？个人数学功底不扎实，不知去哪里找。。:upset:

圆柱半径0.01m，计算区域原点在圆心，x[-0.1 , 0.2]m，y[-0.08 , 0.08]m，入口速度0.01/s
运动粘度从0.02到2e-6,对应雷诺数为0.01到100；
使用二维网格， 转OpenFOAM网格厚度为0.003m。在计算阻力系数已考虑。

入口速度为固定值，压力为零梯度
圆柱上下边界使用的是固壁边界条件（这个设置有问题），速度固定为0，压力为0梯度。
出口速度为零梯度，压力为固定值0。

网格如下图，共12800个网格。

目前获得的数据有：
Re=100 ，Cd= 2.002967;
Re=10 ，Cd= 5.223205;
Re=0.1 ,Cd= 339.979;

这个与理论值不相符：

之后上下边界使用无滑移slip条件：
Re=0.1 ,Cd= 190.0602.

请大家指出问题，谢谢。