[jherb]

I would like to add an equation like the following to my solver:





is a 3x3 matrix, and are 3-dim vectors. All are fields (volVectorField, volTensorField), so are functions of the location (cells).


Does a build-in method exists in OpenFOAM to solve such an equation (actually solving this equation for all cells)?


Thank you for any help.

[agustinvo]

Hi


just write

y=B/A

When doing this with a tensor A and a vector B, it does

y = inv(A) * B

Please check here.

[dlahaye]

See https://www.linkedin.com/pulse/openf...menico-lahaye/ and references cited.

Kind wishes, Domenico Lahaye.

[jherb]

Originally Posted by agustinvo

Hi


just write

y=B/A

When doing this with a tensor A and a vector B, it does

y = inv(A) * B

Please check here.

Thank you. This works.

[syavash]

Originally Posted by jherb

Thank you. This works.

Dear Joachim,

Would you give some elaboration on what the type of vector and matrix is?
I have a similar problem in which I have defined the matrix of coefficient as

scalarSquareMatrix A(3, 3, 0);

and the coefficient vector and the unknown vector as

scalarDiagonalMatrix Coeff(3, 0);
scalarDiagonalMatrix xx(3, 0);

But when I input the following line in my code

xx = Coeff/A;

I get a bunch of errors. I suppose I should change the types of the matrix and vectors but I am not sure how.

Syavash

[jherb]

This is the principle implementation I use:

    volVectorField rightSide
    (
        something
    );

    volTensorField leftSide
    (
        I
    );

    volVectorField result = inv(leftSide) & rightSide;