CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   [groovyBC]error occurred when start to simulate convective BC (http://www.cfd-online.com/Forums/openfoam-solving/86626-groovybc-error-occurred-when-start-simulate-convective-bc.html)

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
Code:

  environment
    {
        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:
Code:

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

Time = 0.01

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#2  Uninterpreted:
#3  Foam::divide(Foam::Field<double>&, double const&,  Foam::UList<double> const&) in  "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
#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"
#5 
 in "/home/merlin/OpenFOAM/merlin-1.7.1/applications/bin/linuxGccDPOpt/cureLaplacianFoam"
#6  __libc_start_main in "/lib/libc.so.6"
#7 
 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!!

-Merlin


All times are GMT -4. The time now is 04:55.