CFD中文网

有没有大佬可以指点一下

两位前辈我现在在多轴组合运动上存在一些问题，想请教一下你们是怎么设置这种组合运动的吗？如果方便能给我一个类似案例设置文件吗？我的问题是一个球体在一个圆柱流域中公转加自转的组合运动。我目前尝试了设置，但是计算出现了发散，我转速很小的情况下库朗数居然会增加到几万。@卡洛 @冷空气LQW

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2406                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh dynamicMultiMotionSolverFvMesh;


// Specification of rotation around centre
rotation1
{
    solidBodyMotionFunction  rotatingMotion;
    rotatingMotionCoeffs
    {
        origin      (0 0 0);
        axis        (0 1 0);
        omega       5; 
    }
}


dynamicMultiMotionSolverFvMeshCoeffs
{
    domain
    {
        solver solidBody;
        cellZone domain;
        solidBodyCoeffs
        {
            $rotation1;
        }
    }

    ball
    {
        cellZone ball;

  
        solver solidBody;

        solidBodyCoeffs
        {
            solidBodyMotionFunction  multiMotion;
            multiMotionCoeffs
            {
                rotation_1
                {
                    $rotation1;
                }
                rotation_2
                {

                    solidBodyMotionFunction rotatingMotion;
                    rotatingMotionCoeffs
                    {
                        origin      (0 0 0.018);
                        axis        (0 1 0);
                        omega      -50; 
                    }
                }
            }
        }
    }
}


// ************************************************************************* //

模型介绍：一个轴承模型，腔体内部有钢球。轴承内圈旋转，外圈静止，球绕轴线公转加绕自心自转。目前一种方法是通过给流体域一个旋转速度，再给球壁面加自转壁面速度。

SRF模型能正确模拟球壁面附近两个速度合成矢量，但实际模型更复杂还包含其它静止区域，需要使用MRF或滑移网格方法。但是在MRF和滑移网格simpleFoam、pimpleFoam只能对绝对速度文件U进行设置，在旋转壁面上无法实现两个速度的合成。针对MRF和动网格能对旋转壁面进行相对速度的设置吗？像SRF那样设置Urel文件？

另外我发现在Fluent中似乎可以直接设置旋转壁面相对某个运动区域设置一个相对速度值。如下图，但我并没有尝试过不确定，我希望在OpenFOAM中进行类似设置。希望有大佬可以给我一些建议，非常感谢！！

最后附上我SRF模型和MRF模型的案例设置文件
MRF 副本.zip 本.zip)
SRF副本.zip

@李东岳 我的模型两侧还有静止区域，只是这个简单模型还没加上，所以SRF是不适用的。

@李东岳 是的。 全部网格应该也行把 ，我看有个案例就是把全部网格MRF了。

我的模型两侧还有静止区域，只是这个简单模型还没加上，所以SRF是不适用的。

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

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

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

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

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

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

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

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

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

/--------------------------------- 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);
}
}
// ************************************************************************* //

dynamicMeshDict

一个圆柱形腔体内部有8个球体障碍物。球体在内部公转加自转来搅动腔体内流体。目前思路是采用滑移网格让腔体公转运动，为球体壁面加自转速度来模拟自转运动。但是我发现在我U文件中如果设置球壁面为rotatingWallVelocity时，整个腔体内流体公转速度就没了，但是网格却是在运动的。如果使用movingWallVelocity就能正常有速度，请问是哪里出现了问题？