[Fabf]

Hi,

little newbie-issue regarding use of class Foam::vectorTools. I want to calculate the angle in deg between two vectors (g and U). Therefore I write

const volScalarField Phib_=vectorTools::degAngleBetween(g_, U_);

To call this function the beginning of the file reads

#include "vectorTools.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
using namespace Foam::vectorTools;
    
namespace Foam
{
namespace compressible
{
namespace RASModels
{

There Error says the function can not be called. So I am quit sure I do not use the right namespace. But how should I call the function?

Note: Application is a turbulence model for buoyancy.

[alexeym]

Hi,

There is no definition of the function for volVectorFields, it is defined only for Vectors.

So either you implement the function for volVectorFields, or iterate over U_ and call the function with Vector arguments.

[Fabf]

Oh Okay. So I tried

const volScalarField DegAngle = acos (g.component(0) * U.component(0) + g.component(1) * U.component(1) + g.component(2) * U.component(2) )/ (normalize(g_)*normalize(U_)))*180.0 / PI;

I can not call the normalize function. Also with quaternion::normalize it did not work. Any idea?

[alexeym]

Hi,

And why do you thinks normalize is defined for vector fields (there is documentation, you know)?

What you are trying to do?

[Fabf]

... I messed it up. Time to go home for today..

Wanted magnitude. So now I calculate the angle between two vectors with

const volScalarField DegAngle = (acos((g_.component(0) * U_.component(0) + g_.component(1) * U_.component(1) + g_.component(2) * U_.component(2)) / (mag(g_)*mag(U_))))*180.0 / constant::mathematical::pi;

and this works.

Sorry.