CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Savonius pimpleDyMFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 8, 2014, 14:50
Default Savonius pimpleDyMFoam
  #1
New Member
 
Vito Fasano
Join Date: Apr 2014
Posts: 14
Rep Power: 9
vitokad is on a distinguished road
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
vitokad is offline   Reply With Quote

Old   September 8, 2014, 16:57
Default
  #2
Senior Member
 
anonymous
Join Date: Aug 2014
Posts: 205
Rep Power: 10
ssss is on a distinguished road
Does it run in serial?

There seems to be a division by zero but it's difficult to say
ssss is offline   Reply With Quote

Old   September 8, 2014, 17:02
Default
  #3
New Member
 
Vito Fasano
Join Date: Apr 2014
Posts: 14
Rep Power: 9
vitokad is on a distinguished road
I don't know, is the first time that happening something like that!
vitokad is offline   Reply With Quote

Old   September 9, 2014, 10:57
Default
  #4
Member
 
Join Date: Jun 2012
Posts: 76
Rep Power: 11
maHein is on a distinguished road
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
maHein is offline   Reply With Quote

Old   September 9, 2014, 18:35
Default
  #5
New Member
 
Vito Fasano
Join Date: Apr 2014
Posts: 14
Rep Power: 9
vitokad is on a distinguished road
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

Last edited by vitokad; September 10, 2014 at 04:18.
vitokad is offline   Reply With Quote

Old   September 11, 2014, 05:48
Default
  #6
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 13
Aurelien Thinat is on a distinguished road
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...
Attached Images
File Type: jpg divergence.jpg (15.9 KB, 28 views)
Aurelien Thinat is offline   Reply With Quote

Old   September 11, 2014, 06:02
Default
  #7
New Member
 
Vito Fasano
Join Date: Apr 2014
Posts: 14
Rep Power: 9
vitokad is on a distinguished road
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..
vitokad is offline   Reply With Quote

Old   September 11, 2014, 07:07
Default
  #8
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 13
Aurelien Thinat is on a distinguished road
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.
Aurelien Thinat is offline   Reply With Quote

Old   September 16, 2014, 03:45
Default
  #9
New Member
 
Vito Fasano
Join Date: Apr 2014
Posts: 14
Rep Power: 9
vitokad is on a distinguished road
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
vitokad is offline   Reply With Quote

Old   September 16, 2014, 08:05
Default
  #10
Senior Member
 
Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 13
Aurelien Thinat is on a distinguished road
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.
Aurelien Thinat is offline   Reply With Quote

Old   September 16, 2014, 09:30
Default
  #11
New Member
 
Vito Fasano
Join Date: Apr 2014
Posts: 14
Rep Power: 9
vitokad is on a distinguished road
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.
vitokad is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
pimpleDyMFoam error message laurentb OpenFOAM Running, Solving & CFD 7 May 13, 2015 05:48
pimpleDymFoam for tidal turbine Jackie Chen OpenFOAM 6 August 18, 2014 11:09
pimpleDyMFoam issue giovanidiniz OpenFOAM Running, Solving & CFD 1 July 5, 2013 07:25
pimpleDyMFoam samiam1000 OpenFOAM 2 September 19, 2012 10:11
Error with pimpleDyMFoam samiam1000 OpenFOAM 2 June 11, 2012 06:21


All times are GMT -4. The time now is 20:33.