[deepbandivadekar]

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:

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]

Bumping this thread as well in hope of grabbing attention of someone who knows how to deal with this.

[TobM]

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.

[deepbandivadekar]

Originally Posted by TobM

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.

[yashjo]

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.

[deepbandivadekar]

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.

[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?

[ano]

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?

[deepbandivadekar]

Originally Posted by ano

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.

Originally Posted by ano

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.

Originally Posted by ano

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]

Does this part of error mean anything to you?

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.

[ano]

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.

        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

        volScalarField rPsi("rPsi", 1.0/psi);

which includes for the perfectGas Model the temperature (which is probably negative)

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.

[deepbandivadekar]

Okay. So I tried your suggestions, except the following which I'm not sure how to implement.

Originally Posted by ano

- 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 ?

- 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:

[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 ***

[khanhrovip]

Hi Deep,
Did you solve your problem? I have the same.

[deepbandivadekar]

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:

#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.

[khanhrovip]

I’m trying with my case, too. Thank you.

[piu58]

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.

[deepbandivadekar]

Originally Posted by piu58

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.