CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Vector Fields and Operators on the Boundary (https://www.cfd-online.com/Forums/openfoam-programming-development/84880-vector-fields-operators-boundary.html)

 Noggin February 10, 2011 14:04

Vector Fields and Operators on the Boundary

Hi Foamers,

just wondering if anyone knows if N Vector Fields are readily available in OpenFoam 1.7.0.

I want to have a geometric field with componants of, say, vectors with 5 elements. I hoping to have such a vector field which includes related algebra. I'm thinking that a geometricField<List<5>,...> would be ok to store the data but would not have the algebra provided?

I tried creating an N vector copying and pasting the vector definition, but it seems the VectorSpace class requires a constant int.

I'm a little muddled about what the best approach would be. If there is a better solution, and any hints about the above muchly appreciated.

Similarly, is an NxN matrix field, to make the Jacobian matrix geometric field possible?

Finally, is it possible, and if so, how would one go about applying the standard field operators only on the boundary conditions?

Many thanks.

 Noggin February 14, 2011 05:37

Just in case anyone else wants to do it, to make an N Vector geometric field, there is a bunch of copying and pasting to be done. Unfortunately it is a constant size, 5 in my case.
Start off by making a new 5 vector: copy vector.C/H, Vector.C/H dimensionedVector, and making new versions of them - I changed the .C files to .H files as well otherwise the wmake system seemed not to find dependencies.
Then update the vectorFieldFwd.H, vectorField.H,vectorIOField.C/H.
Next update everything in the volFields directory. (All this is done in a local directory so as not to recompile the whole source - OpenFOAM very good code to enable this!)
So far, so good, but the boundaries need to be individually updated as well. I have only included the calculated boundary so far which included a copy and paste of the fvPatchFields/Fwd.H/C but not fvPatchField, and the calculated boundary directory.

Not sure if this is a useful solution generally.

Next to try and get the operators to work on it...

Still not sure about operators operating exclusively on boundaries ...

 kathrin_kissling February 15, 2011 03:43

Hi Noggin,

I'm not really sure what you are planning to to, but I think everything you want is there in OpenFoam-1.6-ext. (There src/VectorN)