[merlyn]

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

   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:

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/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

[gschaider]

Originally Posted by merlyn

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

   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:

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/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

The problem does not occur in groovyBC directly (although it is possible that it set a boundary to 0 and now your solver tries to divide by it) as it is not in the stack-trace. Have a look which values your solver divides by

Bernhard