BUG: airRecirculationRoom: reactingParcelFoam throws FPE for h/T
Summary
$FOAM_TUTORIALS/lagrangian/reactingParcelFoam/airRecirculationRoom
throws FPE for
Build : 34e226dfe3-20211220 OPENFOAM=2112 version=com
Arch : "LSB;label=32;scalar=64"
Exec : reactingParcelFoam -parallel
The same test case runs fine with v2106
.
Build : c15bfde3cb-20210624 OPENFOAM=2106
Arch : "LSB;label=32;scalar=64"
Exec : reactingParcelFoam -parallel
Steps to reproduce
cd $FOAM_TUTORIALS/lagrangian/reactingParcelFoam/airRecirculationRoom
./Allrun
A similar crashing case example can be found in EP1765, where the case throws negative temperature FATAL error.
Relevant logs and/or images
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 0.81706851, Final residual = 1.0147201e-08, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.76799964, Final residual = 1.5325339e-08, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.44172895, Final residual = 6.39482e-09, No Iterations 2
DILUPBiCGStab: Solving for CO2, Initial residual = 9.114497e-08, Final residual = 9.114497e-08, No Iterations 0
DILUPBiCGStab: Solving for O2, Initial residual = 9.1708675e-06, Final residual = 7.6550065e-11, No Iterations 1
DILUPBiCGStab: Solving for H2O, Initial residual = 4.273067e-08, Final residual = 4.273067e-08, No Iterations 0
[4] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[4] #1 Foam::sigFpe::sigHandler(int) at ??:?
[4] #2 ? in /lib64/libpthread.so.0
[4] #3 Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[4] #4 Foam::PBiCGStab::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[4] #5 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[4] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[4] #7 Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
[4] #8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
[4] #9 ? at ??:?
[4] #10 __libc_start_main in /lib64/libc.so.6
[4] #11 ? at /home/abuild/rpmbuild/BUILD/glibc-2.31/csu/../sysdeps/x86_64/start.S:122
[pikachu:38369:0:38369] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
Environment information
Problematic environment:
base0 = base
base1 = develop
api = 2112
patch = 0
HEAD = 34e226dfe3
version = com
compiler = Clang (system)
= clang version 9.0.1
mpi = SYSTEMOPENMPI
= mpirun (Open MPI) 1.10.7.0.5e373bf1fd
OS = Description: openSUSE Leap 15.3
opts = linux64ClangDPInt32Opt
Non-problematic environment:
base0 = release
base1 = v2106
api = 2106
patch = 0
HEAD = c15bfde3cb
version = v2106
compiler = Clang (system)
= clang version 9.0.1
mpi = SYSTEMOPENMPI
= mpirun (Open MPI) 1.10.7.0.5e373bf1fd
OS = Description: openSUSE Leap 15.3
opts = linux64ClangDPInt32Opt