CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   boundary cond correct in U, not in paraview (https://www.cfd-online.com/Forums/openfoam-solving/134530-boundary-cond-correct-u-not-paraview.html)

ajnewman April 30, 2014 19:32

boundary cond correct in U, not in paraview
 
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


All times are GMT -4. The time now is 14:03.