CFD中文网

@Mania 怀疑fft变换时参数选错了，OF不大可能错，因为你说之前的都没问题。

@李东岳 确实如李老师所说，将固定壁面边界条件Omega的值设置为1e-6后，能算了！！！

@李东岳 谢谢李老师！我估计是边界条件的事，有些边界条件用的fixedValue，0，但是我改为10e-6后仍旧报错，我再试试。ps：帖子用的ubuntu系统，没中文输入法，无奈只能用英文发帖了

The error is as follows. However, if I input command "unset FOAM_SIGFPE" before "pimpleFoam", everything is OK. the OpenFOAM could run as usual. Could anybody help me for this stupid question?

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 1


PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
bounding k, min: 0 max: 0.00375 average: 0.00375
bounding omega, min: 0 max: 0.559 average: 0.559
RAS
{
    model           kOmegaSST;
    turbulence      on;
    printCoeffs     on;
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

No MRF models present

No finite volume options present
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleMomentumTransportModel<Foam::transportModel> > >, Foam::IncompressibleMomentumTransportModel<Foam::transportModel> >::F2() const at ??:?
#6  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleMomentumTransportModel<Foam::transportModel> > >, Foam::IncompressibleMomentumTransportModel<Foam::transportModel> >::F23() const at ??:?
#7  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleMomentumTransportModel<Foam::transportModel> > >, Foam::IncompressibleMomentumTransportModel<Foam::transportModel> >::correctNut() at ??:?
#8  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pimpleFoam"
Floating point exception (core dumped)

@Mania 因为之前踩过这个坑

就是renumberMesh后会重新生成一套网格数据，之前的还有个0文件夹网格数据。reconstructPar后，你将0文件夹删掉，再用paraFoam打开试试

是不是计算之前renumberMesh操作了？你只把renumberMesh后的网格和计算数据放一起，之前的不要放到一个目录里面，这样再用paraFoam打开试试。
我怀疑是乱码了

发个云图看看噻

因为前后网格没对应好。前后网格需要“完美”对应上，包括节点。其实你用cyclicAMI这个边界条件要好一点。

貌似有答案了。。。。

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ("libfvMotionSolvers.so");

motionSolver solidBody;

cellZone front_fan;

solidBodyMotionFunction rotatingMotion;

origin (2.5 2.5 0);
axis (0 0 1);
omega 3.2832; // rad/s
///////////////////////////////////////
dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ("libfvMotionSolvers.so");

motionSolver solidBody;

cellZone back_fan;

solidBodyMotionFunction rotatingMotion;

origin (7.5 2.5 0);
axis (0 0 1);
omega -3.2832; // rad/s

其中，dynamicMeshDict文件中的内容如下：
[0_1533984525943_dynamicMeshDict](正在上传 100%)

如2楼图所示，如果想要设置两个区域的动网格，而且这两个区域运动形式不一样（如：旋转的方向不同），请问如何在dynamicMeshDict文件中设置。一个区域的运动比较容易实现。谢谢！!

@东岳 据说速度没有原生的ubuntu快，但是我没有亲自测试，window界面不习惯，不过能算也行了

@东岳 哈哈 谢谢东岳老师！blueCFD编译foamMonitor 了，但是没有gnuplot所以没法运行（ps:虽然blueCFD安装完后有独立安装GNUPLOT）。现在官网提供blueCFD的安装exe文件了，还是很方便的，不妨一试。另外东岳老师在算例计算中，是如何看残差的？还是心中自有残差？

我想在blueCFD中使用foamMonitor命令，但是不成功，我没法了，有会的吗？谢谢了。^.^