[teofil]

I would like to perform a very simple thermal analysis using laplacianFOAM, but I don't understand well the results, mostly if in comparison with other codes like CalculiX and Elmer (have a look to the attached pictures).

For the 0/T file boundaries conditions I used:

internalField uniform 20;

boundaryField
{
heatSource
{
type fixedGradient;
gradient uniform 520.0208;
value 20;
}

"auto1|auto2|auto3|auto4" (paths for natural convectio...)
{
type codedMixed;

refValue uniform 0;
refGradient uniform 0;
valueFraction uniform 1;

name freeConvection;
code
#{
float tAmb = 20;
float h = 5.246;
this->refGrad() = -h*(this->db().time().value() - tAmb);
this->refValue() = 20;
#};
codeOptions
#{
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude
#};
codeInclude
#{
#include "fvCFD.H"
#include <cmath>
#include <iostream>
#};
}
}

Could anyone help me to understand how I could set-up the correct boundaries?

Alex

[Sophia]

Hi Alex,

Have you solved your problem? I am working on something similar and I am very much a beginner but I do see something that doesn't convince me. In the free convection code you write this->db().time().value().. but that gives you a value of time not the value of the temperature that you need.

Let me know if you get somewhere!

Sophia