foam-extend-4.0 动网格未按照预期进行强迫运动
-
左侧为入口，右侧为出口；前一个圆柱在y方向做正弦运动，后一个圆柱固定。
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麻烦各位大佬支支招！提前拜谢！
-
在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 在 foam-extend-4.0 动网格未按照预期进行强迫运动 中说：
@李东岳 目前还没有在别的机器测试过，主要不确定是自己代码的原因还是向cfd-online所说是软件安装的问题。
麻烦您方便的时候看一下我贴出的代码有问题么，感谢感谢！我去找其他机器看看能不能动。我换到其他机器也不能动，终端输出结果与贴出内容相同。
-
如图所示，第一个时间步网格会出现些许运动，之后就停止了。
2024年11月11日 16:19
3/12
2024年11月12日 14:01