CFD中文网

Gunther

@李东岳好的，李老师

Gunther

@李东岳李老师您好！这个链接失效了，麻烦补一下，谢谢

Gunther

cfd-online也有人提出了相同的问题，但Jasak表示没有问题，或许真的是安装版本的问题，请问用fe40的兄弟们有遇到相同问题吗？

如果可以正常运行，可以分享一下fe40安装包吗？

Gunther

@李东岳麻了，哭。。。我多么希望是我自己的问题

Gunther

@李东岳我要用到基于fe40的fsiFoam。这么说我的代码没有问题么？

Gunther

如图所示，第一个时间步网格会出现些许运动，之后就停止了。

Gunther

@Gunther 在 foam-extend-4.0 动网格未按照预期进行强迫运动中说：

@李东岳目前还没有在别的机器测试过，主要不确定是自己代码的原因还是向cfd-online所说是软件安装的问题。
麻烦您方便的时候看一下我贴出的代码有问题么，感谢感谢！我去找其他机器看看能不能动。

我换到其他机器也不能动，终端输出结果与贴出内容相同。

Gunther

@李东岳目前还没有在别的机器测试过，主要不确定是自己代码的原因还是向cfd-online所说是软件安装的问题。
麻烦您方便的时候看一下我贴出的代码有问题么，感谢感谢！我去找其他机器看看能不能动。

Gunther

@李东岳李老师您好！我说的是虚拟机fe40自带的tutorial，我运行后发现还是有问题。网格不动。

Gunther

在cfd-online找到相关的问题，说是程序安装的事，但是我改用dyfluid的虚拟机运行此算例，圆柱依旧无法运动。但是我确实看不出来哪里设置出现问题。

Dear Foamers,

I am trying to run the default wingMotion tutorial using sixDoFRigidBodyDisplacement in foam-extend 4.0. 
However, it seems to me that the body is only moving at the very first time step and then all motion ends. 
There is also no output regarding information about restraints, centre of mass etc. in the log file. 
I ran the tutorial in OpenFOAM 5.0 without a problem but can not get it to run in foam-extend 4.0.

Any ideas? Thanks in advance.

After some more tries I found that pimpleDyMFoam only moves the body at the first time step and then all movement is stopped. 
This is the same for sixDoFRigidBodyDisplacement or oscillatingDisplacement.

edit: I somehow messed the installation up. pimpleFoam worked excellent so I did not consider this first but after reinstallation the problem ist solved.

ps：我同时运行了教程自带的wingMotion sixDoFRigidBodyDisplacement算例，机翼也未发生运动。

Gunther

左侧为入口，右侧为出口；前一个圆柱在y方向做正弦运动，后一个圆柱固定。

计算域.png

dynamicMeshDict代码如下：

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

dynamicFvMesh      dynamicMotionSolverFvMesh;

motionSolverLibs ("libfvMotionSolver.so");

solver            displacementLaplacian;
//solver            velocityLaplacian;

diffusivity       inverseDistance (cylinder1);

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

0/pointDisplacement代码如下：

dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    frontAndBack
    {
        type            empty;
    }
    cylinder2
    {
         type            fixedValue;
         value           uniform (0 0 0);
    }
    cylinder1
    {
        type            oscillatingDisplacement;
        amplitude		(0 -1 0);
        omega			3.14; 
        value           uniform (0 0 0);
    }
    topAndBottom
    {
         type            symmetryPlane;
    }
    outlet
    {
         type            fixedValue;
         value           uniform (0 0 0);
    }
    inlet
    {
         type            fixedValue;
         value           uniform (0 0 0);      
    }
}

0/U代码如下：

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (1 0 0);

boundaryField
{
    frontAndBack
    {
        type            empty;
    }
    cylinder2
    {
        type            fixedValue;
        value           uniform (0 0 0); 
    }
    cylinder1
    {
        type            movingWallVelocity;
        value           uniform (0 0 0);
    }
    topAndBottom
    {
        type            symmetryPlane;
    }
    outlet
    {
        type            zeroGradient;
    }
    inlet
    {
        type            fixedValue;
        value           uniform (1 0 0);        
    }
}

现在出现的问题是，圆柱根本不动，以下是终端的输出结果：

Volume: new = 19.8431 old = 19.8431 change = 0 ratio = 0
Courant Number mean: 0.00403834 max: 0.0257354 velocity magnitude: 1.3391
Time = 3.052

GAMG:  Solving for cellDisplacementx, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG:  Solving for cellDisplacementy, Initial residual = 9.0384e-06, Final residual = 9.0384e-06, No Iterations 0
volume continuity errors : volume = 19.8431, max error = 2.21646e-12, sum local = 2.26977e-19, global = -4.50435e-20
BiCGStab:  Solving for Ux, Initial residual = 3.66548e-05, Final residual = 6.20396e-14, No Iterations 1
BiCGStab:  Solving for Uy, Initial residual = 0.000136469, Final residual = 3.21293e-13, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.000114889, Final residual = 5.46747e-06, No Iterations 10
time step continuity errors : sum local = 1.36508e-11, global = 4.94412e-13, cumulative = 1.47017e-06
DICPCG:  Solving for p, Initial residual = 9.3308e-06, Final residual = 9.74749e-07, No Iterations 18
time step continuity errors : sum local = 2.43382e-12, global = 4.80079e-13, cumulative = 1.47017e-06
ExecutionTime = 149.77 s  ClockTime = 149 s

麻烦各位大佬支支招！提前拜谢！

Gunther

@尚善若水谢谢大佬，笔芯

Gunther

@tidedrinker 非常感谢您分享您的经验！
我这个算例拿starccm跑过（二维网格，没有z轴厚度）结果并不理想，我也十分苦恼。

关于starccm网格导入openfoam，我个人的经验是先有starccm导出.ccm格式的三维网格，之后用ccm26ToFoam *.ccm操作就可以将三维网格导入openfoam了，值得注意的是，openfoam没有纯二维，都是通过设置empty边界条件实现二维的模拟，需要注意的是z方向网格数为1，这里会用到extrudeMesh的操作，具体内容见田东老师的视频
希望能够帮到你！