[ttjs3]

Dear Foamers,
I am trying to set up a new injection model. I need to decide the injection position based on the time varying flow rate through an orifice in the domain.

In controldict we can set up a function object to calculate the flow rate through a cross section of the orifice and write it out to a file. But is there a way to pass the flow rate as a parameter into my new injection model so that I can use it to calculate the parcel injection position?

In Fortran we can just store it in a variable and pass it on across subroutines. But it's much less straightforward in c++. I am really stumped. Could someone help me out? Thanks in advance.

Callum

[anishtain4]

Hi,

If you are developing your own solver I don't think there's a way to pass a function object variable to the solver as function objects will be compiled during the run time and solver is not aware of them. I think the best way would be to read the averaging area from an input dictionary and then implement that part in your solver. Averaging a part in your solver should not be hard.

[ttjs3]

Hiya, think I probably give the solution I found myself for this problem in case someone else has the need to do something similar.

Basically I just put the corresponding parts of code used to calculate the flow in function objects into my injection model. The functions that I put into my codes are:

sampledSurfaceFaces
getFieldValues
combineFields

The codes for these functions can be found in the original OpenFOAM codes. They just need to be adapted to the new code where they are used.

In addition, the following header files are also needed:
#include "interpolationCellPoint.H"
#include "ListListOps.H"
#include "Pstream.H"

Hope this can be helpful to other users. Cheers.