||March 28, 2011 15:11
[groovyBC]error occurred when start to simulate convective BC
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
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"
#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"
#6 __libc_start_main in "/lib/libc.so.6"
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!!