关于openfoam中如何实现公转加自转的运动？

[yanoo]

/--------------------------------- C++ -----------------------------------
========= |
\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\ / O peration | Website: https://openfoam.org
\ / A nd | Version: 9
\/ M anipulation |
*---------------------------------------------------------------------------*/
FoamFile
{
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
waiquan
{
type movingWallVelocity;
value uniform (0 0 0);
}
inlet
{
type movingWallVelocity;
value uniform (0 0 0);
}
neiquan
{
type movingWallVelocity;
value uniform (0 0 0);
}
outlet
{
type movingWallVelocity;
value uniform (0 0 0);
}
ball1
{
type rotatingWallVelocity;
relative yes;
origin (0 -0.01675 0);
axis (0 0 -1);
omega constant 837.76;
value uniform (0 0 0);
}
ball2
{
type rotatingWallVelocity;
relative yes;
origin (0 0.01675 0);
axis (0 0 -1);
omega constant 837.76;
value uniform (0 0 0);
}
ball3
{
type rotatingWallVelocity;
relative yes;
origin (-0.011844 -0.011844 0);
axis (0 0 -1);
omega constant 837.76;
value uniform (0 0 0);
}
ball4
{
type rotatingWallVelocity;
relative yes;
origin (-0.01675 0 0);
axis (0 0 -1);
omega constant 837.76;
value uniform (0 0 0);
}
ball5
{
type rotatingWallVelocity;
relative yes;
origin (-0.011844 0.011844 0);
axis (0 0 -1);
omega constant 837.76;
value uniform (0 0 0);
}
ball6
{
type rotatingWallVelocity;
relative yes;
origin (0.011844 0.011844 0);
axis (0 0 -1);
omega constant 837.76;
value uniform (0 0 0);
}
ball7
{
type rotatingWallVelocity;
relative yes;
origin (0.011844 -0.011844 0);
axis (0 0 -1);
omega constant 837.76;
value uniform (0 0 0);
}
ball8
{
type rotatingWallVelocity;
relative yes;
origin (0.01675 0 0);
axis (0 0 -1);
omega constant 837.76;
value uniform (0 0 0);
}
}
// ************************************************************************* //

[yanoo]

这个是我只加了动网格旋转速度，U边界条件用rotatingWallVelocity=0时的速度图，网格能正常旋转，但是流体公转速度似乎没有了

[李东岳]

是，rotatingWallVelocity跟movingWallVelocity冲突了。你有没有想过公转部分用SRF，改变参考系来实现

[yanoo]

@李东岳 感谢您的回复，SRF试过了可以正常转动，但是当我模型复杂一点就不适用了我还是想通过滑移网格来实现。我现在应该怎么修改实现呢？

[李东岳]

滑移网格还有自转的对应的rotatingWallVelocity跟movingWallVelocity冲突，这方面应该做算法的改动，但是我还没有思路怎么处理。

[yanoo]

@李东岳 这么难搞，有其他方法能实现这个公转加自转的模型吗。我看似乎可以参考Fluent中的设置方法，不知道openfoam中能不能实现？FLUENT动网格系列：公转&自转

[李东岳]

这个是两套滑移网格，不是rotatingVelocity。两套滑移网格不存在冲突，不过我没这么玩过，你可以试试。

[yanoo]

@李东岳 好的谢谢您，我再试试。

[tidedrinker]

可以参考这里

就结果来说，能实现，但我没看懂是怎么实现的。或者说，我不太理解这些参数是怎么起作用的，不理解内圈的旋转中心是怎么实现随外圈旋转而变化的。

[yanoo]

@tidedrinker 我针对这个案例Allrun后看不了结果，这两个box文件是什么意思好像要把他们合并？

[yanoo]

@tidedrinker 已经能正常查看结果，它应该是对内圈定义了两个旋转运动来实现公转加自转的，所以中心会随着外圈改变。这种跟上边Fluent两套滑移网格的方法类似，我再针对我的案例再试一下。感谢！

[jsy]

打扰一下，各位老师，我也有一个类似的问题，我希望用pisoFoam跑出tutorials/incompressible/SRFPimpleFoam/rotor2D案例的效果，模型没有怎么更改，可以运行但是不旋转，由于我需要使用CFDEM所以只能用pisoFoam。这片帖子上的rotatingWallVelocity我也试着做了一下，中间有一点，并没有随时间变化。有什么效果能在pisoFoam实现旋转的效果吗？

[yanoo]

@jsy SRF模型就是不旋转的，速度有旋转吗？你要想模型旋转就用滑移网格。

[jsy]

我的模型就是一个搅拌，中间有扇叶，滑移网格是设置dynamicMeshDict吧，pisoFoam求解器可以用吗，我是初学还没有接触过动网格。

[jsy]

我需要的效果就是这样，我就是不知道pisoFoam求解器可以跑吗，我试试吧，谢谢你。

[yanoo]

pisoFoam好像不支持动网格把？我也不太清楚，问下李老师把。@李东岳

[李东岳]

新版的pimpleFoam支持动网格， @jsy 你这个在标准算例里面就有，叫mixer2dAMI好像

[jsy]

谢谢老师，我现在的版本计算mixer2dAMI没有效果，我后期还需要使用CFDEM，使用pisoFoam，动网格好像不太适合。

[yanoo]

我用SRF模型定义Urel文件相对速度，一切正常，有自转公转速度合成方向。改用MRF模型，cellzoon为全部单元，定义U文件，发现我额外定义速度的边界上，附加不上MRF公转速度？改为Urel文件仍然不对。李老师能帮忙看下怎么回事吗，我真被这相对绝对速度搞混了。@李东岳

[李东岳]

我觉得你这个最好的方法就是SRF。

你MRF怎么定义的公转速度？