[Lorg]

Dear Foamers,

I have got a problem using buoyantSimpleFoam for a heat transfer case. Whenever I use a velocity boundary condition for the two inlets higher than a certain value (the maximum value that works is 13.2) I get the following error code at fourth time step:

Starting time loop

Time = 1

DILUPBiCGStab:  Solving for h, Initial residual = 1, Final residual = 0.00170966, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00557072, No Iterations 7
time step continuity errors : sum local = 0.10991, global = 0.00504564, cumulative = 0.00504564
rho max/min : 1.17397 1.00626
DILUPBiCGStab:  Solving for omega, Initial residual = 0.00769473, Final residual = 3.32288e-06, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 1, Final residual = 0.00646544, No Iterations 1
ExecutionTime = 12.99 s  ClockTime = 14 s

Time = 2

DILUPBiCGStab:  Solving for h, Initial residual = 0.707512, Final residual = 0.00134701, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.0188145, Final residual = 0.000139613, No Iterations 4
time step continuity errors : sum local = 10.5002, global = -0.508423, cumulative = -0.503377
rho max/min : 5.68128 1.03367
DILUPBiCGStab:  Solving for omega, Initial residual = 0.029101, Final residual = 7.18413e-06, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 0.960528, Final residual = 0.00728786, No Iterations 1
ExecutionTime = 18.35 s  ClockTime = 19 s

Time = 3

DILUPBiCGStab:  Solving for h, Initial residual = 0.638894, Final residual = 0.0021587, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.011095, Final residual = 9.57394e-05, No Iterations 4
time step continuity errors : sum local = 6.01821, global = 1.52077, cumulative = 1.01739
rho max/min : 3.06482 0.945442
DILUPBiCGStab:  Solving for omega, Initial residual = 0.0956394, Final residual = 0.00013444, No Iterations 1
bounding omega, min: -59.9638 max: 30099.9 average: 71.9794
DILUPBiCGStab:  Solving for k, Initial residual = 0.205061, Final residual = 0.000447495, No Iterations 1
ExecutionTime = 23.77 s  ClockTime = 25 s


Time = 4


DILUPBiCGStab:  Solving for h, Initial residual = 0.197362, Final residual = 0.000685917, No Iterations 1
[0] [1] #0  Foam::error::printStack(Foam::Ostream&)[2] #0  Foam::error::printStack(Foam::Ostream&)#0  Foam::error::printStack(Foam::Ostream&) at ??:?
[2] #1  Foam::sigFpe::sigHandler(int) at ??:?
[1] #1  Foam::sigFpe::sigHandler(int) at ??:?
[0] #1  Foam::sigFpe::sigHandler(int) at ??:?
[2] #2  ? at ??:?
[0] #2  ? at ??:?
[1] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
[2] #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[0] #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[1] #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[0] #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[2] #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[1] #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[2] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[0] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[1] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[0] #7   at ??:?
[2] #7   at ??:?
[1] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[0] #8   in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[2] #8   in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[1] #8  Foam::fvMatrix<double>::solve()Foam::fvMatrix<double>::solve()Foam::fvMatrix<double>::solve() in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[0] #9   in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[2] #9   in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[1] #9  ??? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[0] #10  __libc_start_main in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[2] #10  __libc_start_main in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[1] #10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #11   in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #11   in "/lib/x86_64-linux-gnu/libc.so.6"
[1] #11  ???[lars-B450M-DS3H:20606] *** Process received signal ***
 in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[lars-B450M-DS3H:20606] Signal: Floating point exception (8)
[lars-B450M-DS3H:20606] Signal code:  (-6)
[lars-B450M-DS3H:20606] Failing at address: 0x3e80000507e
[lars-B450M-DS3H:20606] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f43a4531f20]
[lars-B450M-DS3H:20606] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f43a4531e97]
[lars-B450M-DS3H:20606] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f43a4531f20]
[lars-B450M-DS3H:20606] [ 3] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xae)[0x7f43a593c5be]
[lars-B450M-DS3H:20606] [ 4] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x90a)[0x7f43a593feda]
[lars-B450M-DS3H:20606] [ 5] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x4af)[0x7f43a594232f]
[lars-B450M-DS3H:20606] [ 6] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x17d)[0x7f43a9804ced]
[lars-B450M-DS3H:20606] [ 7] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x15b)[0x56486602f70b]
[lars-B450M-DS3H:20606] [ 8] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x117)[0x56486602fa97]
[lars-B450M-DS3H:20606] [ 9] buoyantSimpleFoam(+0x33f8a)[0x564865fdcf8a]
[lars-B450M-DS3H:20606] [10] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f43a4514b97]
[lars-B450M-DS3H:20606] [11] buoyantSimpleFoam(+0x36e7a)[0x564865fdfe7a]
[lars-B450M-DS3H:20606] *** End of error message ***
 in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[lars-B450M-DS3H:20608] *** Process received signal ***
[lars-B450M-DS3H:20608] Signal: Floating point exception (8)
[lars-B450M-DS3H:20608] Signal code:  (-6)
[lars-B450M-DS3H:20608] Failing at address: 0x3e800005080
[lars-B450M-DS3H:20608] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f602c8a3f20]
[lars-B450M-DS3H:20608] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f602c8a3e97]
[lars-B450M-DS3H:20608] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f602c8a3f20]
[lars-B450M-DS3H:20608] [ 3] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xae)[0x7f602dcae5be]
[lars-B450M-DS3H:20608] [ 4] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x90a)[0x7f602dcb1eda]
[lars-B450M-DS3H:20608] [ 5] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x4af)[0x7f602dcb432f]
[lars-B450M-DS3H:20608] [ 6] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x17d)[0x7f6031b76ced]
[lars-B450M-DS3H:20608] [ 7] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x15b)[0x5556ea2e970b]
[lars-B450M-DS3H:20608] [ 8] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x117)[0x5556ea2e9a97]
[lars-B450M-DS3H:20608] [ 9] buoyantSimpleFoam(+0x33f8a)[0x5556ea296f8a]
[lars-B450M-DS3H:20608] [10] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f602c886b97]
[lars-B450M-DS3H:20608] [11] buoyantSimpleFoam(+0x36e7a)[0x5556ea299e7a]
[lars-B450M-DS3H:20608] *** End of error message ***
 in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/buoyantSimpleFoam"
[lars-B450M-DS3H:20607] *** Process received signal ***
[lars-B450M-DS3H:20607] Signal: Floating point exception (8)
[lars-B450M-DS3H:20607] Signal code:  (-6)
[lars-B450M-DS3H:20607] Failing at address: 0x3e80000507f
[lars-B450M-DS3H:20607] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f9525c34f20]
[lars-B450M-DS3H:20607] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f9525c34e97]
[lars-B450M-DS3H:20607] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f9525c34f20]
[lars-B450M-DS3H:20607] [ 3] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xae)[0x7f952703f5be]
[lars-B450M-DS3H:20607] [ 4] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x90a)[0x7f9527042eda]
[lars-B450M-DS3H:20607] [ 5] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x4af)[0x7f952704532f]
[lars-B450M-DS3H:20607] [ 6] /opt/openfoam7/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x17d)[0x7f952af07ced]
[lars-B450M-DS3H:20607] [ 7] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x15b)[0x55ec2854a70b]
[lars-B450M-DS3H:20607] [ 8] buoyantSimpleFoam(_ZN4Foam8fvMatrixIdE5solveEv+0x117)[0x55ec2854aa97]
[lars-B450M-DS3H:20607] [ 9] buoyantSimpleFoam(+0x33f8a)[0x55ec284f7f8a]
[lars-B450M-DS3H:20607] [10] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f9525c17b97]
[lars-B450M-DS3H:20607] [11] buoyantSimpleFoam(+0x36e7a)[0x55ec284fae7a]
[lars-B450M-DS3H:20607] *** End of error message ***

I tried converting the case to simpleFoam by changing the solver in the controlDict file, adapting the p_rgh file to p and changing the entries in the fvSolution file for p_rgh to p and it worked just fine using higher velocities.
Does anyone have any idea why the simulation crashes with buoyantSImpleFoam at higher velocities? Is there something wrong with the boundary conditions? Please find the case attached in case you want to try it on your own pc.
I would be very happy if anyone could help me since I have been struggling with this problem for a long time now!

Best regards,

Lorg

[Lorg]

If it helps there is a table of all boundary conditions on one glance for the better overview attached.
Cheers!