[charitonas]

Hello Foamers,

I am using rhoCentralFoam 2.1.x and 2.2.2. I am trying to add a new stress tensor instead of tauMC.
The usual tensor is:

volTensorField tauMC("tauMC", muEff*dev2(Foam::T(fvc::grad(U))));

so now i want to check some features of the solver and now i am trying to change it to:

volTensorField tauMC("tauMC", muEff*dev2(Foam::T(fvc::grad(5*U))));

If i put a value higher than 5 then the simulation crash with the below message:

[1] #0 Foam::error::printStack(Foam::Ostream&)[0] #0 Foam::error::printStack(Foam::Ostream&) in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int) in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigHandler(int) in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2 in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2 in "/usr/lib/libc.so.6"
[1] #3 Foam::ePsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/usr/lib/libc.so.6"
[0] #3 Foam::ePsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[1] #4 Foam::ePsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[0] #4 Foam::ePsiThermo<Foam::pureMixture<Foam::sutherlandTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[1] #5 in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
[0] #5 

[1] in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
[1] #6 __libc_start_main[0] in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
[0] #6 __libc_start_main in "/usr/lib/libc.so.6"
[1] #7 in "/usr/lib/libc.so.6"
[0] #7 

[1] in "/home/jlight/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/rhoCentralFoam"

I assume that when there is a huge difference between the U and the new vector
(5*U) then it crushes. Does anyone know if there is any way to increase the "tolerance" so i can put a very high value for the stress tensor?

Thanks