[sahm]

Hi FOAMers

I have made a solver that part of it does solid heat diffusion. I want to add anisotropic thermal conductivity to this code, however part of the code has a convection boundary condition that requires to use the thermal conductivity normal to the boundary condition.

my problem is when I am defining a dot product of thermal conductivity (lambda) and the normal vector of the boundary surface:

vectorField NrmlVctr = this->patch().nf(); // vector normal to boundary

And Lambda is defined as a volVectorField in the createFields file.

This line gives an error:

this->valueFraction() = 1.0 / ( 1.0 + (lambda.boundaryField() & NrmlVctr)/(ConvCoeff_*mag(delta)) );

The code does not compile and gives an error that says operator& is not defined for these types of operands. Can someone help me resolve this problem?

Thanks.

[sahm]

Ok, it seems I got that problem kinda sorted out. I had to have it defined as:

(lambda.boundaryField()[patchindex].patchInternalField() & SurfaceNormalVector)

However, it seems the problem is more complicated than that. I need to define the thermal conductivity as a Tensor and modify heat diffusion equation in the solver based on that.
Do I need to change the form of heat diffusion equation? Or the discretization in OpenFOAM works based on that?

Thanks.

[sahm]

One more question. I am trying to redefine the convective boundary condition for Anisotropic material. How am I supposed to implement it? I have asked this question before for isotropic material, and someone helped me out. But that method used mixed patch boundary condition type. How can I use a similar thing?

Thanks.