[ajnewman]

Hello,

I created a codedFixedBC with the following code

vector dir1 = vector(1,0,0);
vector dir2 = vector(0,1,0);
scalarField x = patch().Cf()&dir1;
scalarField y = patch().Cf()&dir2;
scalarField r = pow(pow(x,2) + pow(y,2),0.5);
operator == (vector(0,0,1)*r);

If it isn't evident the boundary condition is supposed to be linear growth in U_3 over a circular domain as a function of r

If I run the code for a single time step I can look at the nonuniform list<vector> for the patch at which this BC is applied. Here I can see a point having the value (0,0,0) which I would expect at the origin. If I examine the data in paraView using data along a line, or a point probe there is a non zero value of U_3 at (0,0,0). In fact the linear profile shows a distinct non zero plateau as the origin is approached.

Am I doing something dumb in paraView or is there actually a problem with my BC?

Any input is appreciated.

-Jensen