CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

PrintStack error with Floating Point Exception rhoCentralFoam OF5.x

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 26, 2018, 15:25
Question PrintStack error with Floating Point Exception rhoCentralFoam OF5.x
  #1
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
Hi guys,

I am working on a compressible flow problem involving heat transfer since a couple of months and was initially using sonicFoam. However I have now started building the case in rhoCentralFoam to see if it gives me better and meaningful outcome (sonicFoam currently isn't).

So I got this error on my run during the third iteration:

Code:
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
fluxScheme: Kurganov

Starting time loop

Mean and max Courant Numbers = 0.756329 14929.4
Time = 3e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
ExecutionTime = 0.51 s  ClockTime = 1 s

Mean and max Courant Numbers = 0.774376 14933.4
Time = 6e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
ExecutionTime = 0.84 s  ClockTime = 1 s

#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::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5  ? at ??:?
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? at ??:?
The sqrt part got my attention and it made obvious that some variable under a square root is of value 0 or perhaps negative and that's why the error. But what exactly is it caused by? Can anyone help?

I want to try and use debugSwitch option in my controlDict but am confused about which one to use. I have opened a separate and related thread here: How to choose correct debugSwitch? (Floating point exception error)
deepbandivadekar is offline   Reply With Quote

Old   April 12, 2018, 07:33
Default
  #2
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
Bumping this thread as well in hope of grabbing attention of someone who knows how to deal with this.
deepbandivadekar is offline   Reply With Quote

Old   April 12, 2018, 08:41
Default
  #3
Member
 
Join Date: Sep 2014
Location: Germany
Posts: 88
Rep Power: 10
TobM is on a distinguished road
No Debug switch will help you in this case. Your maximum Courant number is way too big.
Decrease the time step size, or better, use adjustTimeStep with a reasonable maxCo.
TobM is offline   Reply With Quote

Old   April 12, 2018, 14:36
Default
  #4
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
Quote:
Originally Posted by TobM View Post
No Debug switch will help you in this case. Your maximum Courant number is way too big.
Decrease the time step size, or better, use adjustTimeStep with a reasonable maxCo.
Thanks for your inputs.

The time step right now is of the order of 10E-7. This already takes hours to finish (projected). I have tried reducing the time step all the way to 10E-12. This doesn't change the outcome in any way. I am convinced that the problem is not so simple.
deepbandivadekar is offline   Reply With Quote

Old   April 12, 2018, 16:07
Default
  #5
New Member
 
Yash Joshi
Join Date: Mar 2017
Posts: 5
Blog Entries: 1
Rep Power: 7
yashjo is on a distinguished road
Hi. What are the boundary conditions you are using and other case details?

It looks like somewhere the solver is encountering a division by zero.
yashjo is offline   Reply With Quote

Old   May 17, 2018, 14:08
Default Update:
  #6
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
Not sure how I missed the last reply to this thread.

Am still struggling with rhoCentralFoam and after trying few other ideas I'm back to this case. I have now different set of printStack errors.


Code:
[3] #0  Foam::error::printStack(Foam::Ostream&)[1] #0  Foam::error::printStack(Foam::Ostream&)[2] [7] #0  Foam::error::printStack(Foam::Ostream&)#0   at ??:?
Foam::error::printStack(Foam::Ostream&)[3] #1  Foam::sigFpe::sigHandler(int)[5] #0  Foam::error::printStack(Foam::Ostream&)[0] #0  Foam::error::printStack(Foam::Ostream&)[4] #0  Foam::error::printStack(Foam::Ostream&)[6] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[7] #1  Foam::sigFpe::sigHandler(int) at ??:?
[1] #1  Foam::sigFpe::sigHandler(int) at ??:?
[3] #2  ? at ??:?
[5] #1  Foam::sigFpe::sigHandler(int) at ??:?
[2] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
 in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&)[6] #1  Foam::sigFpe::sigHandler(int) at ??:?
[4] #1  Foam::sigFpe::sigHandler(int) at ??:?
[1] #2  ? at ??:?
[7] #2  ? at ??:?
[2] #2  ? at ??:?
[5] #2  ? at ??:?
[0] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[6] #2  ? at ??:?
[3] #4   in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[4] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[1] #4   in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[7] #4   at ??:?
[5] #4   at ??:?
[2] #4   at ??:?
[3] #5   at ??:?
[6] #4   at ??:?
[0] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)? at ??:?
[4] #4   at ??:?
[1] #5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[7] #5   at ??:?
[5] #5   at ??:?
[2] #5   at ??:?
[3] #6  __libc_start_mainFoam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[6] #5  ???? at ??:?
[0] #5   in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #7   at ??:?
[4] #5  ? at ??:?
[1] #6  __libc_start_main?? at ??:?
[7] #6  __libc_start_main? at ??:?
[5] #6  __libc_start_main at ??:?
[2] #6  __libc_start_main at ??:?
[6] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #7   at ??:?
[mae-hen-8163-db:09495] *** Process received signal ***
[mae-hen-8163-db:09495] Signal: Floating point exception (8)
[mae-hen-8163-db:09495] Signal code:  (-6)
[mae-hen-8163-db:09495] Failing at address: 0x3e900002517
[mae-hen-8163-db:09495] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fb3052eb4b0]
[mae-hen-8163-db:09495] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fb3052eb428]
[mae-hen-8163-db:09495] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fb3052eb4b0]
[mae-hen-8163-db:09495] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7fb3067a5728]
[mae-hen-8163-db:09495] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:09495] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:09495] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fb3052d6830]
[mae-hen-8163-db:09495] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:09495] *** End of error message ***
 in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7   in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #7   in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #7   at ??:?
[0] #6  __libc_start_main at ??:?
[4] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #7  ???? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #7   in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #7  ??? at ??:?
[mae-hen-8163-db:09499] *** Process received signal ***
[mae-hen-8163-db:09499] Signal: Floating point exception (8)
[mae-hen-8163-db:09499] Signal code:  (-6)
[mae-hen-8163-db:09499] Failing at address: 0x3e90000251b
[mae-hen-8163-db:09499] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7efcda8b94b0]
[mae-hen-8163-db:09499] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7efcda8b9428]
[mae-hen-8163-db:09499] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7efcda8b94b0]
[mae-hen-8163-db:09499] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7efcdbd73728]
[mae-hen-8163-db:09499] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:09499] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:09499] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7efcda8a4830]
[mae-hen-8163-db:09499] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:09499] *** End of error message ***
 at ??:?
[mae-hen-8163-db:09493] *** Process received signal ***
[mae-hen-8163-db:09493] Signal: Floating point exception (8)
[mae-hen-8163-db:09493] Signal code:  (-6)
[mae-hen-8163-db:09493] Failing at address: 0x3e900002515
[mae-hen-8163-db:09493] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f34318534b0]
[mae-hen-8163-db:09493] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f3431853428]
[mae-hen-8163-db:09493] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f34318534b0]
[mae-hen-8163-db:09493] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f3432d0d728]
[mae-hen-8163-db:09493] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:09493] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:09493] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f343183e830]
[mae-hen-8163-db:09493] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:09493] *** End of error message ***
 at ??:?
[mae-hen-8163-db:09497] *** Process received signal ***
[mae-hen-8163-db:09497] Signal: Floating point exception (8)
[mae-hen-8163-db:09497] Signal code:  (-6)
[mae-hen-8163-db:09497] Failing at address: 0x3e900002519
[mae-hen-8163-db:09497] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fd6cb47a4b0]
[mae-hen-8163-db:09497] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fd6cb47a428]
[mae-hen-8163-db:09497] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fd6cb47a4b0]
[mae-hen-8163-db:09497] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7fd6cc934728]
[mae-hen-8163-db:09497] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:09497] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:09497] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fd6cb465830]
[mae-hen-8163-db:09497] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:09497] *** End of error message ***
 at ??:?
[mae-hen-8163-db:09494] *** Process received signal ***
[mae-hen-8163-db:09494] Signal: Floating point exception (8)
[mae-hen-8163-db:09494] Signal code:  (-6)
[mae-hen-8163-db:09494] Failing at address: 0x3e900002516
[mae-hen-8163-db:09494] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f1a224a74b0]
[mae-hen-8163-db:09494] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f1a224a7428]
[mae-hen-8163-db:09494] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f1a224a74b0]
[mae-hen-8163-db:09494] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f1a23961728]
[mae-hen-8163-db:09494] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:09494] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:09494] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f1a22492830]
[mae-hen-8163-db:09494] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:09494] *** End of error message ***
 at ??:?
[mae-hen-8163-db:09498] *** Process received signal ***
[mae-hen-8163-db:09498] Signal: Floating point exception (8)
[mae-hen-8163-db:09498] Signal code:  (-6)
[mae-hen-8163-db:09498] Failing at address: 0x3e90000251a
[mae-hen-8163-db:09498] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f49e3cfd4b0]
[mae-hen-8163-db:09498] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f49e3cfd428]
[mae-hen-8163-db:09498] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f49e3cfd4b0]
[mae-hen-8163-db:09498] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f49e51b7728]
[mae-hen-8163-db:09498] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:09498] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:09498] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f49e3ce8830]
[mae-hen-8163-db:09498] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:09498] *** End of error message ***
 at ??:?
[mae-hen-8163-db:09496] *** Process received signal ***
[mae-hen-8163-db:09496] Signal: Floating point exception (8)
[mae-hen-8163-db:09496] Signal code:  (-6)
[mae-hen-8163-db:09496] Failing at address: 0x3e900002518
[mae-hen-8163-db:09496] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f978ea864b0]
[mae-hen-8163-db:09496] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f978ea86428]
[mae-hen-8163-db:09496] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f978ea864b0]
[mae-hen-8163-db:09496] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f978ff40728]
[mae-hen-8163-db:09496] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:09496] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:09496] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f978ea71830]
[mae-hen-8163-db:09496] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:09496] *** End of error message ***
 at ??:?
[mae-hen-8163-db:09492] *** Process received signal ***
[mae-hen-8163-db:09492] Signal: Floating point exception (8)
[mae-hen-8163-db:09492] Signal code:  (-6)
[mae-hen-8163-db:09492] Failing at address: 0x3e900002514
[mae-hen-8163-db:09492] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f6817a434b0]
[mae-hen-8163-db:09492] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f6817a43428]
[mae-hen-8163-db:09492] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f6817a434b0]
[mae-hen-8163-db:09492] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f6818efd728]
[mae-hen-8163-db:09492] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:09492] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:09492] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f6817a2e830]
[mae-hen-8163-db:09492] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:09492] *** End of error message ***

I don't even know where to start looking! Totally clueless.

fvSchemes and fvSolution are set exactly as the base case (which works). Boundary conditions include supersonic inlet and zeroGradient outlet.



Can anyone point me to something?
deepbandivadekar is offline   Reply With Quote

Old   May 18, 2018, 14:00
Default
  #7
ano
Member
 
ano
Join Date: Jan 2017
Location: Delft
Posts: 58
Rep Power: 8
ano is on a distinguished road
Quote:
Boundary conditions include supersonic inlet and zeroGradient outlet.
Perhaps it is better to use a wavetransmissive outlet to avoid numerical shocks. Perhaps it is a good idea to initialize the fields first with an approximate solution (from sonicFoam e.g.). It would be probably easier if you could give some more information on the physics and geometry.

Do you still get a residual of zero in your logfile?
ano is offline   Reply With Quote

Old   May 18, 2018, 17:22
Default
  #8
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
Quote:
Originally Posted by ano View Post
Perhaps it is better to use a wavetransmissive outlet to avoid numerical shocks.
Edit: My mistake. Yes, I do have pressure BC as waveTransmissive. zeroGradient is for velocity at outlet.

Quote:
Originally Posted by ano View Post
Perhaps it is a good idea to initialize the fields first with an approximate solution (from sonicFoam e.g.). It would be probably easier if you could give some more information on the physics and geometry.
That's an interesting proposition. I don't know how to do that though. Can you give me some directions?
Well it's a typical case of a blunt object in a supersonic flow. Thats's about it. I am using an axisymmetric mesh.


Quote:
Originally Posted by ano View Post
Do you still get a residual of zero in your logfile?
The residual zero is shown only for density and that is so because of the solver: diagonal, so those are expected to be zero. However, that was the previous case. As of now there's nothing in the log because the solver crashes right at the beginning.
deepbandivadekar is offline   Reply With Quote

Old   May 19, 2018, 05:53
Default
  #9
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
Does this part of error mean anything to you?

Code:
Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&)

If there's an invalid number inside the root (e.g. -ve) then it is likely that square root function isn't working.
deepbandivadekar is offline   Reply With Quote

Old   May 20, 2018, 15:30
Default
  #10
ano
Member
 
ano
Join Date: Jan 2017
Location: Delft
Posts: 58
Rep Power: 8
ano is on a distinguished road
From the few information you give, one could probably guess the following:

The square root error is most likely from the calculation of the speed of sound.

Code:
        volScalarField c("c", sqrt(thermo.Cp()/thermo.Cv()*rPsi));


As Greenshields et al. wrote in their paper about the solver
one problem are negative temperatures (see "Implementation of semi-discrete, non-staggered central schemes in a colocated, polyhedral, finite volume framework, for high-speed viscous flows").
They say it would be due to the fact that you calculate the temperature from the total energy minus kinetic energy and since the total energy and velocity used for the kinetic energy are calculated sequentially this difference could become negative. So this is probably the problem, since it affects the calculations of the square root in the speed of sound calculation.

Since you have constant specific heats (correct?), the problem in the sqrt() should be rPsi which includes psi


Code:
        volScalarField rPsi("rPsi", 1.0/psi);
which includes for the perfectGas Model the temperature (which is probably negative)

Code:
template<class Specie>
inline Foam::scalar Foam::perfectGas<Specie>::psi(scalar p, scalar T) const
{
    return 1.0/(this->R()*T);
}
But in the end this does not really help you besides knowing that your temperature is probably calculated negative somewhere. Since you got the problem directly at starting the simulation, you should try to either start with good initial conditions or give a less demanding problem which you then change stepwise to your real problem.

- For initializing the solution you could just use your results from sonicFoam (even if they don't make much sense) as starting point for rhoCentralFoam. Or perhaps rhoSimpleFoam or rhoPimpleFoam.
- start with an incompressible or subsonic totalPressure inlet boundary condition.
- start with a low pressure ratio between inlet and outlet and increase it over time. Set your pressure at the outlet initially to a value slightly below the one in the flow region. Calculate some steps, decrease the pressure a little bit, calculate again some steps. If this helps, you can later use a function doing this automatically for you.
- limit the pressure and/or temperature (e.g. to values above 1Pa, depending on your problem).
- try to use a coarser grid.
ano is offline   Reply With Quote

Old   May 21, 2018, 11:09
Default
  #11
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
Okay. So I tried your suggestions, except the following which I'm not sure how to implement.

Quote:
Originally Posted by ano View Post

- For initializing the solution you could just use your results from sonicFoam (even if they don't make much sense) as starting point for rhoCentralFoam. Or perhaps rhoSimpleFoam or rhoPimpleFoam.

These did not change the outcome. I was surprised, esp for subsonic case which I expected to be straight forward. So I think something entirely different is the cause. Similar errors throughout. Although, why coarser grid? Did you mean finer ?
Quote:

- start with an incompressible or subsonic totalPressure inlet boundary condition.
- start with a low pressure ratio between inlet and outlet and increase it over time. Set your pressure at the outlet initially to a value slightly below the one in the flow region. Calculate some steps, decrease the pressure a little bit, calculate again some steps. If this helps, you can later use a function doing this automatically for you.
- limit the pressure and/or temperature (e.g. to values above 1Pa, depending on your problem).
- try to use a coarser grid.
Nonetheless, I've decided to share the entire (modified) case folder. Please see attached files. This case is basically same as a tutorial.


Reproducing associated error right at the beginning for reference:

Code:
[3] #0  Foam::error::printStack(Foam::Ostream&)[7] #0  Foam::error::printStack(Foam::Ostream&)[6] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1  Foam::sigFpe::sigHandler(int) at ??:?
[7] #1  Foam::sigFpe::sigHandler(int)[2] #0  Foam::error::printStack(Foam::Ostream&)[5] [4] #0  Foam::error::printStack(Foam::Ostream&)[0] #0  Foam::error::printStack(Foam::Ostream&)#0  Foam::error::printStack(Foam::Ostream&)[1] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[7] #2  ? at ??:?
[6] #1  Foam::sigFpe::sigHandler(int) at ??:?
[4] #1  Foam::sigFpe::sigHandler(int) at ??:?
[3] #2  ? at ??:?
 at ??:?
 at ??:?
[1] #1  Foam::sigFpe::sigHandler(int)[5] #1  Foam::sigFpe::sigHandler(int) at ??:?
[2] #1  Foam::sigFpe::sigHandler(int) in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&)[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[6] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[4] #2  ? at ??:?
[0] #2  ? at ??:?
[7] #4   at ??:?
[1] #2  ? at ??:?
[5] #2  ? at ??:?
[2] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[3] #4   in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[4] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[5] #4   at ??:?
[2] #4   at ??:?
[7] #5   at ??:?
[6] #4   at ??:?
[1] #4   at ??:?
[0] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)?Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[3] #5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[5] #5   at ??:?
[7] #6  __libc_start_main at ??:?
[4] #5   at ??:?
[6] #5   at ??:?
[1] #5   at ??:?
[0] #5   at ??:?
[2] #5  ?????? at ??:?
[3] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[7] #7  ? at ??:?
[5] #6  __libc_start_main at ??:?
[6] #6  __libc_start_main at ??:?
[4] #6  __libc_start_main at ??:?
[0] #6  __libc_start_main at ??:?
[2] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #7  ? at ??:?
[1] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #7   in "/lib/x86_64-linux-gnu/libc.so.6"
[6] #7   in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #7   in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #7   in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #7  ? at ??:?
[mae-hen-8163-db:03457] *** Process received signal ***
[mae-hen-8163-db:03457] Signal: Floating point exception (8)
[mae-hen-8163-db:03457] Signal code:  (-6)
[mae-hen-8163-db:03457] Failing at address: 0x3e900000d81
[mae-hen-8163-db:03457] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f8082c154b0]
[mae-hen-8163-db:03457] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f8082c15428]
[mae-hen-8163-db:03457] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f8082c154b0]
[mae-hen-8163-db:03457] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f80840cf728]
[mae-hen-8163-db:03457] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:03457] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:03457] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f8082c00830]
[mae-hen-8163-db:03457] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:03457] *** End of error message ***
 in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #7  ???? at ??:?
[mae-hen-8163-db:03453] *** Process received signal ***
[mae-hen-8163-db:03453] Signal: Floating point exception (8)
[mae-hen-8163-db:03453] Signal code:  (-6)
[mae-hen-8163-db:03453] Failing at address: 0x3e900000d7d
[mae-hen-8163-db:03453] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fe350db84b0]
[mae-hen-8163-db:03453] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fe350db8428]
[mae-hen-8163-db:03453] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fe350db84b0]
[mae-hen-8163-db:03453] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7fe352272728]
[mae-hen-8163-db:03453] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:03453] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:03453] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fe350da3830]
[mae-hen-8163-db:03453] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:03453] *** End of error message ***
?? at ??:?
[mae-hen-8163-db:03455] *** Process received signal ***
[mae-hen-8163-db:03455] Signal: Floating point exception (8)
[mae-hen-8163-db:03455] Signal code:  (-6)
[mae-hen-8163-db:03455] Failing at address: 0x3e900000d7f
[mae-hen-8163-db:03455] [ 0]  at ??:?
[mae-hen-8163-db:03450] *** Process received signal ***
[mae-hen-8163-db:03450] Signal: Floating point exception (8)
[mae-hen-8163-db:03450] Signal code:  (-6)
[mae-hen-8163-db:03450] Failing at address: 0x3e900000d7a
[mae-hen-8163-db:03450] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f2dfb1cb4b0]
[mae-hen-8163-db:03450] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f2dfb1cb428]
[mae-hen-8163-db:03450] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f2dfb1cb4b0]
[mae-hen-8163-db:03450] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f2dfc685728]
[mae-hen-8163-db:03450] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:03450] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:03450] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f2dfb1b6830]
[mae-hen-8163-db:03450] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:03450] *** End of error message ***
/lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fdc534b14b0]
[mae-hen-8163-db:03455] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fdc534b1428]
[mae-hen-8163-db:03455] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fdc534b14b0]
[mae-hen-8163-db:03455] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7fdc5496b728]
[mae-hen-8163-db:03455] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:03455] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:03455] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fdc5349c830]
[mae-hen-8163-db:03455] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:03455] *** End of error message ***
 at ??:?
[mae-hen-8163-db:03454] *** Process received signal ***
[mae-hen-8163-db:03454] Signal: Floating point exception (8)
[mae-hen-8163-db:03454] Signal code:  (-6)
[mae-hen-8163-db:03454] Failing at address: 0x3e900000d7e
[mae-hen-8163-db:03454] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f717dab74b0]
[mae-hen-8163-db:03454] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f717dab7428]
[mae-hen-8163-db:03454] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f717dab74b0]
[mae-hen-8163-db:03454] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f717ef71728]
[mae-hen-8163-db:03454] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:03454] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:03454] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f717daa2830]
[mae-hen-8163-db:03454] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:03454] *** End of error message ***
 at ??:?
[mae-hen-8163-db:03451] *** Process received signal ***
[mae-hen-8163-db:03451] Signal: Floating point exception (8)
[mae-hen-8163-db:03451] Signal code:  (-6)
[mae-hen-8163-db:03451] Failing at address: 0x3e900000d7b
[mae-hen-8163-db:03451] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fabcf0314b0]
[mae-hen-8163-db:03451] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7fabcf031428]
[mae-hen-8163-db:03451] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7fabcf0314b0]
[mae-hen-8163-db:03451] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7fabd04eb728]
[mae-hen-8163-db:03451] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:03451] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:03451] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7fabcf01c830]
[mae-hen-8163-db:03451] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:03451] *** End of error message ***
 at ??:?
[mae-hen-8163-db:03452] *** Process received signal ***
[mae-hen-8163-db:03452] Signal: Floating point exception (8)
[mae-hen-8163-db:03452] Signal code:  (-6)
[mae-hen-8163-db:03452] Failing at address: 0x3e900000d7c
[mae-hen-8163-db:03452] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f94784cb4b0]
[mae-hen-8163-db:03452] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f94784cb428]
[mae-hen-8163-db:03452] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f94784cb4b0]
[mae-hen-8163-db:03452] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f9479985728]
[mae-hen-8163-db:03452] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:03452] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:03452] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f94784b6830]
[mae-hen-8163-db:03452] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:03452] *** End of error message ***
 at ??:?
[mae-hen-8163-db:03456] *** Process received signal ***
[mae-hen-8163-db:03456] Signal: Floating point exception (8)
[mae-hen-8163-db:03456] Signal code:  (-6)
[mae-hen-8163-db:03456] Failing at address: 0x3e900000d80
[mae-hen-8163-db:03456] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f83de8454b0]
[mae-hen-8163-db:03456] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f83de845428]
[mae-hen-8163-db:03456] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f83de8454b0]
[mae-hen-8163-db:03456] [ 3] /opt/openfoam5/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam4sqrtERNS_5FieldIdEERKNS_5UListIdEE+0x28)[0x7f83dfcff728]
[mae-hen-8163-db:03456] [ 4] rhoCentralFoam(_ZN4Foam4sqrtINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_+0x10a)[0x49663a]
[mae-hen-8163-db:03456] [ 5] rhoCentralFoam[0x42d71e]
[mae-hen-8163-db:03456] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f83de830830]
[mae-hen-8163-db:03456] [ 7] rhoCentralFoam[0x436899]
[mae-hen-8163-db:03456] *** End of error message ***
Attached Files
File Type: gz case.tar.gz (7.7 KB, 5 views)
deepbandivadekar is offline   Reply With Quote

Old   May 23, 2018, 13:10
Default
  #12
New Member
 
Khanh
Join Date: Dec 2015
Posts: 4
Rep Power: 9
khanhrovip is on a distinguished road
Hi Deep,
Did you solve your problem? I have the same.
khanhrovip is offline   Reply With Quote

Old   May 23, 2018, 15:29
Default
  #13
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
No. Still stuck with this.


I came across this thread and although I know I have more than enough swap as well as ram, tried running this in series instead and I get this error:
Code:
#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::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5  ? at ??:?
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? at ??:?
Floating point exception (core dumped)

I've a feeling that somewhere I have not set up the case correctly.
deepbandivadekar is offline   Reply With Quote

Old   May 23, 2018, 22:39
Default
  #14
New Member
 
Khanh
Join Date: Dec 2015
Posts: 4
Rep Power: 9
khanhrovip is on a distinguished road
I’m trying with my case, too. Thank you.
khanhrovip is offline   Reply With Quote

Old   May 24, 2018, 01:13
Default
  #15
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 727
Rep Power: 13
piu58 is on a distinguished road
I downloaded the case. Unfortunately, there is no mesh. Log files let me assume, that you got the mesh from "somewhere".

I recommend to start with a simpler case, which includes all of the physics at a very simple or at least a rather simple geometry. Use blockmesh for that. It would be best, if the case would be based on a reference solutions which includes at least most of the physics.

In one of your first posts you gave the Courant numbers, around 0.7 average (slightly the high from my experience) and more that 10 thousand at maximum, which is way to high. A decent mesh gives small differences between mean and max Co, so you assume at least one of the problems here.
__________________
Uwe Pilz
--
Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)
piu58 is offline   Reply With Quote

Old   May 25, 2018, 10:33
Default
  #16
Senior Member
 
Deep
Join Date: Oct 2017
Posts: 180
Rep Power: 7
deepbandivadekar is on a distinguished road
Quote:
Originally Posted by piu58 View Post
I downloaded the case. Unfortunately, there is no mesh. Log files let me assume, that you got the mesh from "somewhere".

I recommend to start with a simpler case, which includes all of the physics at a very simple or at least a rather simple geometry. Use blockmesh for that. It would be best, if the case would be based on a reference solutions which includes at least most of the physics.

In one of your first posts you gave the Courant numbers, around 0.7 average (slightly the high from my experience) and more that 10 thousand at maximum, which is way to high. A decent mesh gives small differences between mean and max Co, so you assume at least one of the problems here.

I've included the mesh file in the zip. It just needs to be exported to foam by gmshToFoam. zip only includes a clean case. This geometry itself I feel is pretty simple. Nothing complicated at all.
The first post is very old, things have changes thus far.
deepbandivadekar is offline   Reply With Quote

Reply

Tags
debugswitch, floating point exception, printstack, rhocentrafoam, sonicfoam

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
icoFoam floating point exception (8) leizhao512 OpenFOAM Running, Solving & CFD 7 November 1, 2018 12:43
rhoCentralFoam - kOmegaSST - Floating point exception DGW OpenFOAM Running, Solving & CFD 1 September 1, 2016 08:45
A floating point exception - SEM Model yansheng STAR-CCM+ 1 April 4, 2016 05:57
Floating point exception from twoPhaseEulerFoam openfoammaofnepo OpenFOAM Running, Solving & CFD 1 March 19, 2016 14:56
Floating point exception (core dumped) for GAMG solver yuhou1989 OpenFOAM Running, Solving & CFD 2 March 24, 2015 20:28


All times are GMT -4. The time now is 17:13.