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算例,机翼也未发生运动。