[Hisuku]

Hi there,

I inherited a project from a classmate I am trying to get running and currently only get crashes for nearly everything I try. The case is a bubble column in a pipe which flows around a simple obstacle (another smaller pipe going transversal trough the other). My classmate used twoPhaseEulerFoam and extended the example given in the official examples tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn.

The error I get is most of the time the same but it comes in two flavors (I will paste both outputs at the end of this text).

What I already tried:
- variations of deltaT and writeInterval in controlDict
- variations of U.* -> 0.5, 1 and 2 m/2
- other meshes
- putting my meshes in the official sample with only the necessary adjustments

Out of all only the smallest flowspeed with very small deltaT and a fine mesh managed to hold out a while but failed a few hours ago as well and is also way to slow (0.26sec out of 5 I need in somewhat around 9 days simulation time).
All simulations run parallelised on 16 cores.

If you need further details, let me know. Ill paste in anythign we need to discuss. I am still somewhat new to openFoam (but have enough fun with it to make it a study project) and even more new to this forum so please be kind ;D

error variation 1:

[2] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[2] #1 Foam::sigFpe::sigHandler(int) at ??:?
[2] #2 ? in /lib64/libc.so.6
[2] #3 ? in /lib64/libm.so.6
[2] #4 Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
[2] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
[2] #6 Foam::dragModels::SchillerNaumann::CdRe() const at ??:?
[2] #7 Foam::dragModel::Ki() const at ??:?
[2] #8 Foam::dragModel::K() const at ??:?
[2] #9 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:?
[2] #10 Foam::twoPhaseSystem::Kd() const at ??:?
[2] #11 ? at ??:?
[2] #12 __libc_start_main in /lib64/libc.so.6
[2] #13 ? at ??:?
[node787:30705] *** Process received signal ***
[node787:30705] Signal: Floating point exception (8)
[node787:30705] Signal code: (-6)
[node787:30705] Failing at address: 0x1059a000077f1
[node787:30705] [ 0] /lib64/libc.so.6(+0x36400)[0x2b55c7789400]
[node787:30705] [ 1] /lib64/libc.so.6(gsignal+0x37)[0x2b55c7789387]
[node787:30705] [ 2] /lib64/libc.so.6(+0x36400)[0x2b55c7789400]
[node787:30705] [ 3] /lib64/libm.so.6(+0x7750)[0x2b55c7242750]
[node787:30705] [ 4] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UL istIdEERKd+0x31)[0x2b55c6811901]
[node787:30705] [ 5] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam3powINS _12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14Geometr icFieldIdT_T0_EEEERKS7_RKNS_11dimensionedIdEE+0x1e 2)[0x2b55c1e9fe22]
[node787:30705] [ 6] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam10dragModels15SchillerNaumann4CdReEv+0x1c1)[0x2b55c2bbcfb1]
[node787:30705] [ 7] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel2KiEv+0xec)[0x2b55c2b9a3fc]
[node787:30705] [ 8] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel1KEv+0x21)[0x2b55c2b9a701]
[node787:30705] [ 9] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam23Blende dInterfacialModelINS_9dragModelEE1KEv+0x43f)[0x2b55c28f2f4f]
[node787:30705] [10] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam14twoPha seSystem2KdEv+0x1e)[0x2b55c28d8eae]
[node787:30705] [11] twoPhaseEulerFoam[0x433333]
[node787:30705] [12] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b55c7775555]
[node787:30705] [13] twoPhaseEulerFoam[0x43fae0]
[node787:30705] *** End of error message ***


error variation 2:

[15] --------------------------------------------------------------------------
A process has executed an operation involving a call to the
"fork()" system call to create a child process. Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your job may hang, crash, or produce silent
data corruption. The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.

The process that invoked fork was:

Local host: [[21830,1],15] (PID 17905)

If you are *absolutely sure* that your application will successfully
and correctly survive a call to fork(), you may disable this warning
by setting the mpi_warn_on_fork MCA parameter to 0.
--------------------------------------------------------------------------
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
[15] #1 Foam::sigFpe::sigHandler(int) at ??:?
[15] #2 ? in /lib64/libpthread.so.0
[15] #3 ? in /lib64/libm.so.6
[15] #4 Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
[15] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
[15] #6 Foam::dragModels::SchillerNaumann::CdRe() const at ??:?
[15] #7 Foam::dragModel::Ki() const at ??:?
[15] #8 Foam::dragModel::K() const at ??:?
[15] #9 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:?
[15] #10 Foam::twoPhaseSystem::Kd() const at ??:?
[15] #11 ? at ??:?
[15] #12 __libc_start_main in /lib64/libc.so.6
[15] #13 ? at ??:?
[node774:17905] *** Process received signal ***
[node774:17905] Signal: Floating point exception (8)
[node774:17905] Signal code: (-6)
[node774:17905] Failing at address: 0x1059a000045f1
[node774:17905] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b6e7ba85630]
[node774:17905] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x2b6e7ba854fb]
[node774:17905] [ 2] /lib64/libpthread.so.0(+0xf630)[0x2b6e7ba85630]
[node774:17905] [ 3] /lib64/libm.so.6(+0x7750)[0x2b6e7b565750]
[node774:17905] [ 4] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UL istIdEERKd+0x31)[0x2b6e7aae0811]
[node774:17905] [ 5] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam3powINS _12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14Geometr icFieldIdT_T0_EEEERKS7_RKNS_11dimensionedIdEE+0x1f 1)[0x2b6e7947ad91]
[node774:17905] [ 6] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam10dragMo dels15SchillerNaumann4CdReEv+0x150)[0x2b6e79b998b0]
[node774:17905] [ 7] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam9dragMod el2KiEv+0x102)[0x2b6e79b8e632]
[node774:17905] [ 8] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam9dragMod el1KEv+0x21)[0x2b6e79b8ea11]
[node774:17905] [ 9] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam23Blende dInterfacialModelINS_9dragModelEE1KEv+0x456)[0x2b6e79b5e646]
[node774:17905] [10] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam14twoPha seSystem2KdEv+0x1e)[0x2b6e79b3750e]
[node774:17905] [11] twoPhaseEulerFoam[0x438773]
[node774:17905] [12] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b6e7bcb4555]
[node774:17905] [13] twoPhaseEulerFoam[0x44ac5b]
[node774:17905] *** End of error message ***
--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 15 with PID 17905 on node node774 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

[SHY22]

Originally Posted by Hisuku

Hi there,

I inherited a project from a classmate I am trying to get running and currently only get crashes for nearly everything I try. The case is a bubble column in a pipe which flows around a simple obstacle (another smaller pipe going transversal trough the other). My classmate used twoPhaseEulerFoam and extended the example given in the official examples tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn.

The error I get is most of the time the same but it comes in two flavors (I will paste both outputs at the end of this text).

What I already tried:
- variations of deltaT and writeInterval in controlDict
- variations of U.* -> 0.5, 1 and 2 m/2
- other meshes
- putting my meshes in the official sample with only the necessary adjustments

Out of all only the smallest flowspeed with very small deltaT and a fine mesh managed to hold out a while but failed a few hours ago as well and is also way to slow (0.26sec out of 5 I need in somewhat around 9 days simulation time).
All simulations run parallelised on 16 cores.

If you need further details, let me know. Ill paste in anythign we need to discuss. I am still somewhat new to openFoam (but have enough fun with it to make it a study project) and even more new to this forum so please be kind ;D

error variation 1:

[2] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[2] #1 Foam::sigFpe::sigHandler(int) at ??:?
[2] #2 ? in /lib64/libc.so.6
[2] #3 ? in /lib64/libm.so.6
[2] #4 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
[2] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
[2] #6 Foam::dragModels::SchillerNaumann::CdRe() const at ??:?
[2] #7 Foam::dragModel::Ki() const at ??:?
[2] #8 Foam::dragModel::K() const at ??:?
[2] #9 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:?
[2] #10 Foam::twoPhaseSystem::Kd() const at ??:?
[2] #11 ? at ??:?
[2] #12 __libc_start_main in /lib64/libc.so.6
[2] #13 ? at ??:?
[node787:30705] *** Process received signal ***
[node787:30705] Signal: Floating point exception (8)
[node787:30705] Signal code: (-6)
[node787:30705] Failing at address: 0x1059a000077f1
[node787:30705] [ 0] /lib64/libc.so.6(+0x36400)[0x2b55c7789400]
[node787:30705] [ 1] /lib64/libc.so.6(gsignal+0x37)[0x2b55c7789387]
[node787:30705] [ 2] /lib64/libc.so.6(+0x36400)[0x2b55c7789400]
[node787:30705] [ 3] /lib64/libm.so.6(+0x7750)[0x2b55c7242750]
[node787:30705] [ 4] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UL istIdEERKd+0x31)[0x2b55c6811901]
[node787:30705] [ 5] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam3powINS _12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14Geometr icFieldIdT_T0_EEEERKS7_RKNS_11dimensionedIdEE+0x1e 2)[0x2b55c1e9fe22]
[node787:30705] [ 6] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam10dragModels15SchillerNaumann4CdReEv+0x1c1)[0x2b55c2bbcfb1]
[node787:30705] [ 7] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel2KiEv+0xec)[0x2b55c2b9a3fc]
[node787:30705] [ 8] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleEulerianInterfacialModels.so(_ZNK4F oam9dragModel1KEv+0x21)[0x2b55c2b9a701]
[node787:30705] [ 9] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam23Blende dInterfacialModelINS_9dragModelEE1KEv+0x43f)[0x2b55c28f2f4f]
[node787:30705] [10] /software/openfoam/v1806/OpenFOAM-v1806/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam14twoPha seSystem2KdEv+0x1e)[0x2b55c28d8eae]
[node787:30705] [11] twoPhaseEulerFoam[0x433333]
[node787:30705] [12] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b55c7775555]
[node787:30705] [13] twoPhaseEulerFoam[0x43fae0]
[node787:30705] *** End of error message ***


error variation 2:

[15] --------------------------------------------------------------------------
A process has executed an operation involving a call to the
"fork()" system call to create a child process. Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your job may hang, crash, or produce silent
data corruption. The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.

The process that invoked fork was:

Local host: [[21830,1],15] (PID 17905)

If you are *absolutely sure* that your application will successfully
and correctly survive a call to fork(), you may disable this warning
by setting the mpi_warn_on_fork MCA parameter to 0.
--------------------------------------------------------------------------
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[15] #1 Foam::sigFpe::sigHandler(int) at ??:?
[15] #2 ? in /lib64/libpthread.so.0
[15] #3 ? in /lib64/libm.so.6
[15] #4 Foam:ow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
[15] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:ow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
[15] #6 Foam::dragModels::SchillerNaumann::CdRe() const at ??:?
[15] #7 Foam::dragModel::Ki() const at ??:?
[15] #8 Foam::dragModel::K() const at ??:?
[15] #9 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:?
[15] #10 Foam::twoPhaseSystem::Kd() const at ??:?
[15] #11 ? at ??:?
[15] #12 __libc_start_main in /lib64/libc.so.6
[15] #13 ? at ??:?
[node774:17905] *** Process received signal ***
[node774:17905] Signal: Floating point exception (8)
[node774:17905] Signal code: (-6)
[node774:17905] Failing at address: 0x1059a000045f1
[node774:17905] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b6e7ba85630]
[node774:17905] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x2b6e7ba854fb]
[node774:17905] [ 2] /lib64/libpthread.so.0(+0xf630)[0x2b6e7ba85630]
[node774:17905] [ 3] /lib64/libm.so.6(+0x7750)[0x2b6e7b565750]
[node774:17905] [ 4] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam3powERNS_5FieldIdEERKNS_5UL istIdEERKd+0x31)[0x2b6e7aae0811]
[node774:17905] [ 5] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTurbulenceModels.so(_ZN4Foam3powINS _12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14Geometr icFieldIdT_T0_EEEERKS7_RKNS_11dimensionedIdEE+0x1f 1)[0x2b6e7947ad91]
[node774:17905] [ 6] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam10dragMo dels15SchillerNaumann4CdReEv+0x150)[0x2b6e79b998b0]
[node774:17905] [ 7] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam9dragMod el2KiEv+0x102)[0x2b6e79b8e632]
[node774:17905] [ 8] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam9dragMod el1KEv+0x21)[0x2b6e79b8ea11]
[node774:17905] [ 9] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam23Blende dInterfacialModelINS_9dragModelEE1KEv+0x456)[0x2b6e79b5e646]
[node774:17905] [10] /software/openfoam/v2112/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libcompressibleTwoPhaseSystem.so(_ZNK4Foam14twoPha seSystem2KdEv+0x1e)[0x2b6e79b3750e]
[node774:17905] [11] twoPhaseEulerFoam[0x438773]
[node774:17905] [12] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b6e7bcb4555]
[node774:17905] [13] twoPhaseEulerFoam[0x44ac5b]
[node774:17905] *** End of error message ***
--------------------------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 15 with PID 17905 on node node774 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

Hi! Have you solved this problem? I'm using the reactingTwoPhaseEulerFoam and I meet the same problem like you. Could you give me some suggestion about how to fixed this trouble?

[Hisuku]

Sadly not really. In my university there is only one person capable of openfoam at another chair. But he only does single phase simulation and was ill most of the time.

As far as I understood there seems to be a weakness with the dragModel in use (specified in phaseProperties). Changing the 'air in water' entry to the IshiiZuber like

(air in water)

{

type IshiiZuber;

residualRe 1e-3;

swarmCorrection

{

type none;

residualAlpha Tomiyama;

l 1.75;

}

}

as described in https://www.cemf.ir/wp-content/uploa...ble-Column.pdf seems to somewhat work better depending on the alpha.air and how it is distributed over the inlet but the error seems to focus on the 'water in air' entry bellow it. Entry #6 of every error output I had always calls the used drag model for 'water in air'. The article I linked mentions other models which are available for twoPhaseEulerFoam. I didn't manage to find any clue about what is supposed to work better. If you get any insight in that I would appreciate it if you could tell me.

[SHY22]

Originally Posted by Hisuku

Sadly not really. In my university there is only one person capable of openfoam at another chair. But he only does single phase simulation and was ill most of the time.

As far as I understood there seems to be a weakness with the dragModel in use (specified in phaseProperties). Changing the 'air in water' entry to the IshiiZuber like

(air in water)

{

type IshiiZuber;

residualRe 1e-3;

swarmCorrection

{

type none;

residualAlpha Tomiyama;

l 1.75;

}

}

as described in https://www.cemf.ir/wp-content/uploa...ble-Column.pdf seems to somewhat work better depending on the alpha.air and how it is distributed over the inlet but the error seems to focus on the 'water in air' entry bellow it. Entry #6 of every error output I had always calls the used drag model for 'water in air'. The article I linked mentions other models which are available for twoPhaseEulerFoam. I didn't manage to find any clue about what is supposed to work better. If you get any insight in that I would appreciate it if you could tell me.

Thanks for your reply! I'll try it on my case and I'll get back to you again if I find anything! Have a great time at work!