CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   SIGFPE in reactingFoam(LES) (http://www.cfd-online.com/Forums/openfoam-bugs/87032-sigfpe-reactingfoam-les.html)

Hýr0 April 9, 2011 14:27

SIGFPE in reactingFoam(LES)
 
Hi,
I obtain a floating point exception when I run reactingFoam using a LesModel. From gdb:

if (turbulentReaction)
(gdb)
14 Cmix*sqrt(turbulence->muEff()/rho/turbulence->epsilon());
(gdb) n

Program received signal SIGFPE, Arithmetic exception.
0x01a9b232 in Foam::divide (res=..., f1=..., f2=...)
at fields/Fields/scalarField/scalarField.C:95
95 BINARY_OPERATOR(scalar, scalar, scalar, /, divide)

In 0 directory I don't assign a value for muEff, rho as well as for epsilon and I don't see the need to assign one in Les.

mevsfedora June 6, 2011 09:08

Quote:

Originally Posted by Hýr0 (Post 302904)
if (turbulentReaction)
(gdb)
14 Cmix*sqrt(turbulence->muEff()/rho/turbulence->epsilon());
(gdb) n

Program received signal SIGFPE, Arithmetic exception.
0x01a9b232 in Foam::divide (res=..., f1=..., f2=...)
at fields/Fields/scalarField/scalarField.C:95
95 BINARY_OPERATOR(scalar, scalar, scalar, /, divide)

In 0 directory I don't assign a value for muEff, rho as well as for epsilon and I don't see the need to assign one in Les.

Hi,

in reactingfoam you use the partially stirred reactor, which needs some fields like k or epsilon.

I'm new to Openfoam, but would like to use the PaSR for my case, also with LES turbulence model.
My Question is: Is there anyone, who already used reactingfoam with LES. How did you calculate epsilon?

Thanks for any answers.


All times are GMT -4. The time now is 16:13.