CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

[groovyBC]error occurred when start to simulate convective BC

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 27, 2011, 20:15
Default [groovyBC]error occurred when start to simulate convective BC
  #1
New Member
 
Wei Ma
Join Date: Feb 2011
Posts: 14
Rep Power: 6
merlyn is on a distinguished road
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
merlyn is offline   Reply With Quote

Old   March 29, 2011, 05:24
Default
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by merlyn View Post
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
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
gschaider is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
TimeVaryingMappedFixedValue irishdave OpenFOAM Running, Solving & CFD 28 May 28, 2015 13:37
Transient DPM-About Start Stop times Julie FLUENT 3 January 24, 2014 08:39
Start star ccm+ remotely with mpiexec Dagnir CD-adapco 0 November 24, 2010 13:33
Convective / Conductive Heat Transfer in Hypersonic flows enigma Main CFD Forum 2 November 1, 2009 23:53
Lets start the public domain CFD-Project! Heinz Wilkening Main CFD Forum 3 March 11, 1999 23:55


All times are GMT -4. The time now is 20:30.