[bepps]

Hi folks, I am trying to develop a codedFixedValue boundary condition in the attached file, and I need help with some API things:

-- How do I interpolate the velocity field at a given (x,y,z) point? (line 157)

-- How do I get the (x,y,z) coordinates of the locations on the patch where the boundary condition is to be applied? (lines 178--182)

-- How do I prescribe the boundary condition? (line 194)

-- Do I need to do anything special to enable this boundary condition to be used in simulations on multiple processors?

My desired b.c. applies to swirling flows, such as a propeller, where the pressure field should vary radially according to the normal momentum equation, dp/dr = rho * Vt^2 / r, where Vt is the tangential (swirl) component of velocity. I need to interpolate the velocity at several points on the outlet, compute the desired pressure distribution, p(r), and then apply that pressure at the actual boundary points in the mesh. The attached code has all this theory implemented, but I am struggling with the OpenFOAM 11 API. Can anyone please help me fill in the blanks? Thanks in advance!

[bjnieuwboer]

Hi,

You could take a look at http://www.wolfdynamics.com/wiki/prog_codestreamBC.pdf. There they implement a parabolic velocity profile using codeStream or codedFixedVale.

I hope this helps you.