URGENT - Simulation setup using Reynolds Tensor Transport (R) - OpenFOAM

beluiz93 · November 29, 2016, 10:59

Hello guys!

I'm currently trying to simulate a LRR and SSG in a flow that is external to a sphere. Both of them, as you guys might now, are types of simulation that use Reynolds Stress Tensor Transport in your bases of calculation.

So, I added a "R" file in my 0/ folder already that comes with my Reynolds Tensor value.

I start the simulation and everything goes fine until around 8 seconds (when I use pimpleFoam) and 3 seconds (when I use pisoFoam). I don't know what it keeps getting wrong because the Courant looks fine, everything looks fine. Suddenly the simulation crashes and this happens:

#0 Foam::error:

rintStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#2 ? in "/lib64/libc.so.6"
#3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#6 Foam::fvMatrix<Foam::SymmTensor<double> >::solveSegregated(Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so"
#7 Foam::fvMatrix<Foam::SymmTensor<double> >::solve(Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so"
#8 Foam::fvMatrix<Foam::SymmTensor<double> >::solve() in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so"
#9 Foam::RASModels::LRR<Foam::IncompressibleTurbulenc eModel<Foam::transportModel> >::correct() in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so"
#10 ? in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
#11 __libc_start_main in "/lib64/libc.so.6"
#12 ? in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
Floating point exception

What is this fvMatrix error? I don't get that. I tried many fvSchemes like limitedLinear 1.0, backward, CrankNicolson and nothing works.

Please, I need help with this! Any answer will be great! Thanks in advance!

Bernardo Ribeiro

khedar · December 2, 2016, 07:51

floating point exception normally happens when you have 0 somewhere in a division. Probable reasons generally are bad mesh or bad boundary conditions.. either of these can lead to very high value of some flow quantity which leads to such error. your solver log (from start) would be more helpful. See for stuff like "bounding k/omega/U etc.."

November 29, 2016, 10:59	URGENT - Simulation setup using Reynolds Tensor Transport (R) - OpenFOAM	#1
beluiz93 New Member Bernardo Ribeiro Join Date: Aug 2016 Posts: 13 Rep Power: 10	Hello guys! I'm currently trying to simulate a LRR and SSG in a flow that is external to a sphere. Both of them, as you guys might now, are types of simulation that use Reynolds Stress Tensor Transport in your bases of calculation. So, I added a "R" file in my 0/ folder already that comes with my Reynolds Tensor value. I start the simulation and everything goes fine until around 8 seconds (when I use pimpleFoam) and 3 seconds (when I use pisoFoam). I don't know what it keeps getting wrong because the Courant looks fine, everything looks fine. Suddenly the simulation crashes and this happens: #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #2 ? in "/lib64/libc.so.6" #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so" #6 Foam::fvMatrix<Foam::SymmTensor<double> >::solveSegregated(Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so" #7 Foam::fvMatrix<Foam::SymmTensor<double> >::solve(Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so" #8 Foam::fvMatrix<Foam::SymmTensor<double> >::solve() in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so" #9 Foam::RASModels::LRR<Foam::IncompressibleTurbulenc eModel<Foam::transportModel> >::correct() in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so" #10 ? in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/pisoFoam" #11 __libc_start_main in "/lib64/libc.so.6" #12 ? in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/pisoFoam" Floating point exception What is this fvMatrix error? I don't get that. I tried many fvSchemes like limitedLinear 1.0, backward, CrankNicolson and nothing works. Please, I need help with this! Any answer will be great! Thanks in advance! Bernardo Ribeiro

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
OpenFOAM Training Jan-Apr 2017, Virtual, London, Houston, Berlin	cfd.direct	OpenFOAM Announcements from Other Sources	0	September 21, 2016 12:50
Reynolds number in 2D airfoil simulation	jck	OpenFOAM Running, Solving & CFD	2	June 3, 2016 04:56
how to calculate reynolds stresses from SA turbulence simulation?	hjj3118	FLUENT	0	May 3, 2016 00:39
Wave simulation using OpenFoam	Ellie	OpenFOAM	2	August 30, 2013 03:10
Convergence of Reynolds Stress Tensor / Fluctuations in OpenFOAM	Anne Lincke	OpenFOAM Running, Solving & CFD	0	July 17, 2012 11:50

December 2, 2016, 07:51		#2
khedar Senior Member khedar Join Date: Oct 2016 Posts: 111 Rep Power: 10	floating point exception normally happens when you have 0 somewhere in a division. Probable reasons generally are bad mesh or bad boundary conditions.. either of these can lead to very high value of some flow quantity which leads to such error. your solver log (from start) would be more helpful. See for stuff like "bounding k/omega/U etc.."