[vitokad]

Hello Foamers,

for my thesis I want to simulate a modified Savonius rotor. I created the mesh of the rotor and the surrounding environment with ICEM, I converted and merged (fluentMeshToFoam and mergeMeshes). I decomposed the domain into 4 subdomains (I have 4 cores), and I tried to start the solver in parallel (mpirun -np 4 pimpleDyMFoam -Parallel> log &).
The result was this:

[0] #0 Foam::error:rintStack(Foam::Ostream&)[1] #0 Foam::error:rintStack(Foam::Ostream&)[2] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[2] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #2 at ??:?
[2] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
[0] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
[1] #4 at ??:?
[2] #4 at ??:?
[0] #4 void Foam::divide<Foam::fvPatchField>(Foam::FieldField< Foam::fvPatchField, double>&, double const&, Foam::FieldField<Foam::fvPatchField, double> const&)void Foam::divide<Foam::fvPatchField>(Foam::FieldField< Foam::fvPatchField, double>&, double const&, Foam::FieldField<Foam::fvPatchField, double> const&)void Foam::divide<Foam::fvPatchField>(Foam::FieldField< Foam::fvPatchField, double>&, double const&, Foam::FieldField<Foam::fvPatchField, double> const&) at ??:?
[1] #5 at ??:?
[2] #5 at ??:?
[0] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[2] #6 at ??:?
[0] #6 at ??:?
[1] #6


[2] at ??:?
[2] #7 __libc_start_main[0] at ??:?
[0] #7 __libc_start_main[1] at ??:?
[1] #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #8 in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #8 in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #8


[0] at ??:?
[vito-Aspire-V5-571:03100] *** Process received signal ***
[vito-Aspire-V5-571:03100] Signal: Floating point exception (8)
[vito-Aspire-V5-571:03100] Signal code: (-6)
[vito-Aspire-V5-571:03100] Failing at address: 0x3e800000c1c
[vito-Aspire-V5-571:03100] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36ff0) [0x7fa1b456cff0]
[vito-Aspire-V5-571:03100] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7fa1b456cf77]
[vito-Aspire-V5-571:03100] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36ff0) [0x7fa1b456cff0]
[vito-Aspire-V5-571:03100] [ 3] /opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRK NS_5UListIdEE+0xc7) [0x7fa1b5866717]
[vito-Aspire-V5-571:03100] [ 4] pimpleDyMFoam(_ZN4Foam6divideINS_12fvPatchFieldEEE vRNS_10FieldFieldIT_dEERKdRKS4_+0x5b) [0x43c98b]
[vito-Aspire-V5-571:03100] [ 5] pimpleDyMFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7vol MeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_ 11dimensionedIdEERKS8_+0x126) [0x466c56]
[vito-Aspire-V5-571:03100] [ 6] pimpleDyMFoam() [0x4289bb]
[vito-Aspire-V5-571:03100] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7fa1b4557de5]
[vito-Aspire-V5-571:03100] [ 8] pimpleDyMFoam() [0x42bb01]
[vito-Aspire-V5-571:03100] *** End of error message ***
[1] at ??:?
[vito-Aspire-V5-571:03101] *** Process received signal ***
[vito-Aspire-V5-571:03101] Signal: Floating point exception (8)
[vito-Aspire-V5-571:03101] Signal code: (-6)
[vito-Aspire-V5-571:03101] Failing at address: 0x3e800000c1d
[vito-Aspire-V5-571:03101] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36ff0) [0x7f2082b9cff0]
[vito-Aspire-V5-571:03101] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7f2082b9cf77]
[vito-Aspire-V5-571:03101] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36ff0) [0x7f2082b9cff0]
[vito-Aspire-V5-571:03101] [ 3] /opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRK NS_5UListIdEE+0xc7) [0x7f2083e96717]
[vito-Aspire-V5-571:03101] [ 4] pimpleDyMFoam(_ZN4Foam6divideINS_12fvPatchFieldEEE vRNS_10FieldFieldIT_dEERKdRKS4_+0x5b) [0x43c98b]
[vito-Aspire-V5-571:03101] [ 5] pimpleDyMFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7vol MeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_ 11dimensionedIdEERKS8_+0x126) [0x466c56]
[vito-Aspire-V5-571:03101] [ 6] pimpleDyMFoam() [0x4289bb]
[vito-Aspire-V5-571:03101] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7f2082b87de5]
[vito-Aspire-V5-571:03101] [ 8] pimpleDyMFoam() [0x42bb01]
[vito-Aspire-V5-571:03101] *** End of error message ***
[2] at ??:?
[vito-Aspire-V5-571:03102] *** Process received signal ***
[vito-Aspire-V5-571:03102] Signal: Floating point exception (8)
[vito-Aspire-V5-571:03102] Signal code: (-6)
[vito-Aspire-V5-571:03102] Failing at address: 0x3e800000c1e
[vito-Aspire-V5-571:03102] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x36ff0) [0x7fc3f09f4ff0]
[vito-Aspire-V5-571:03102] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x37) [0x7fc3f09f4f77]
[vito-Aspire-V5-571:03102] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x36ff0) [0x7fc3f09f4ff0]
[vito-Aspire-V5-571:03102] [ 3] /opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKdRK NS_5UListIdEE+0xc7) [0x7fc3f1cee717]
[vito-Aspire-V5-571:03102] [ 4] pimpleDyMFoam(_ZN4Foam6divideINS_12fvPatchFieldEEE vRNS_10FieldFieldIT_dEERKdRKS4_+0x5b) [0x43c98b]
[vito-Aspire-V5-571:03102] [ 5] pimpleDyMFoam(_ZN4FoamdvINS_12fvPatchFieldENS_7vol MeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKNS_ 11dimensionedIdEERKS8_+0x126) [0x466c56]
[vito-Aspire-V5-571:03102] [ 6] pimpleDyMFoam() [0x4289bb]
[vito-Aspire-V5-571:03102] [ 7] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0x7fc3f09dfde5]
[vito-Aspire-V5-571:03102] [ 8] pimpleDyMFoam() [0x42bb01]
[vito-Aspire-V5-571:03102] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 3100 on node vito-Aspire-V5-571 exited on signal 8 (Floating point exception).


"vito-Aspire-V5-571" is my pc.

I don't know which is the error..any help is appreciated.

My case file: https://www.dropbox.com/s/i5s9rkyh41...to.tar.gz?dl=0


Best regards,

Vito

[ssss]

Does it run in serial?

There seems to be a division by zero but it's difficult to say

[vitokad]

I don't know, is the first time that happening something like that!

[maHein]

Hey,

I just had a look at your mesh. It looks like that your cyclicAMI boundary conditions are not overlapping. In particular, AMI2 does contain faces from your outlet boundary conditions.

You should have a look at these issues first.

Regards

[vitokad]

Hello maHein, i've updated the case.

When I try to start the solver (not in parallel), after a few number of iterations this is the result:

Courant Number mean: 6.13956e+14 max: 6.63231e+21
deltaT = 4.85945e-33
--> FOAM Warning :
From function Time:perator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 7 to 8 to distinguish between timeNames at time 0.000249669
Time = 0.00024966869

solidBodyMotionFunctions::rotatingMotion::transfor mation(): Time = 0.000249669 transformation: ((0 0 0) (1 (0 0 0.000784359)))
AMI: Creating addressing and weights between 200 source faces and 200 target faces
AMI: Patch source sum(weights) min/max/average = 0.999952, 1.00018, 1.00004
AMI: Patch target sum(weights) min/max/average = 0.999983, 1.00017, 1.00005
smoothSolver: Solving for Ux, Initial residual = 0.884486, Final residual = 9.41008e-07, No Iterations 58
smoothSolver: Solving for Uy, Initial residual = 0.863333, Final residual = 9.49771e-07, No Iterations 66
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8
at ??:?
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
at ??:?
Eccezione in virgola mobile (core dump creato)


The mesh now seems good. I think that there are some errors in fvschemes/fvsolution..

Someone can help me?

Thanks a lot

Vito


Case file: https://www.dropbox.com/s/0ewbspt795...lo.tar.gz?dl=0

[Aurelien Thinat]

Good morning Vito,

I had a look at your case. FvSolution and fvSchemes look good. I'll bet on a mesh problem : the velocity/turbulence are diverging at the end of the prism layer (see the attached picture).

You should try to smooth your mesh.

Anyway, do you have any article relative to this test case ? Velocity, pressure, Cd...

[vitokad]

Unfortunately I have not provided material relevant to this case, but only in such cases that I retrieved on the internet ..

Doing a checkMesh I get an error, I wrote a post about it in the sub-forum concerning the conversion of mesh from Ansys. However, by starting the simulation with only the command pimpleDyMFoam, it start..

[Aurelien Thinat]

Your mesh is globally not bad but the prism layers transition is really sharp : the cell volume jumps from 1 to 10 (or even more).

You should try to remesh locally at the trailing edge of the blade. Having a smooth transition near the blade could really help the solver.

[vitokad]

Well, I'd like to update you on the situation. Despite the error on the mesh, which gives me 4 interior points are not used, I was able to start the simulation in parallel. The problem was in the AMI interfaces: adding the string

preservePatches(AMI1 AMI2)

the simulation start in parallel without problem. If someone wants to know more about the case please write here.

Best regards,

Vito

[Aurelien Thinat]

Does it run until the convergence is achieve ?

I launched your case (without preserve patch option) in pimpleFoam (sliding mesh off), and it lead to a numerical divergence...

EDIT : It also diverges with simpleFoam. Both in single or parrallel run.

[vitokad]

I made some changes compared to the case that I had previously loaded, both in fvSchemes, and I create the rotor zone, which previously had not specified.