CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   [groovyBC]error occurred when start to simulate convective BC (

merlyn March 28, 2011 15:11

[groovyBC]error occurred when start to simulate convective BC
Hi guys.

Previously, I set the boundary condition(surface temperature) as fixedValue, and everything was Ok.

Then, I tried to use groovyBC to modified 0/T, to make it a convective BC.
I installed groovyBC and ran demos with no problems.

But, after the modification of boundary as

        type                  groovyBC;
        variables              "htot=55.0;Tinf=423.0;k=0.53";
        valueExpression    "Tinf";
        fractionExpression  "1.0/(1.0 + k/(mag(Tinf-T)*htot))";

where, T is defined for the solid's internal temperature as:

internalField  uniform 298;
I ran my solver as curelaplacianFoam, which is only a little different from the laplacianFoam, error popped out as:

Time = 0.01

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam171/lib/linuxGccDPOpt/"
#1  Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linuxGccDPOpt/"
#2  Uninterpreted:
#3  Foam::divide(Foam::Field<double>&, double const&,  Foam::UList<double> const&) in  "/opt/openfoam171/lib/linuxGccDPOpt/"
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField,  Foam::volMesh> > Foam::operator/<Foam::fvPatchField,  Foam::volMesh>(Foam::dimensioned<double> const&,  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>  const&) in  "/home/merlin/OpenFOAM/merlin-1.7.1/applications/bin/linuxGccDPOpt/cureLaplacianFoam"
 in "/home/merlin/OpenFOAM/merlin-1.7.1/applications/bin/linuxGccDPOpt/cureLaplacianFoam"
#6  __libc_start_main in "/lib/"
 in "/home/merlin/OpenFOAM/merlin-1.7.1/applications/bin/linuxGccDPOpt/cureLaplacianFoam"
Floating point exception

Could anyone know how should I make it work? Basically, the case is I need a fixed value for Tinf, 423, and set the initial body as a uniform value, 298. Should not be a big deal.

Please help me, thank you very much in advance!!


All times are GMT -4. The time now is 10:37.