[Allankey]

Hi, I am trying to model a non uniform boundary condition at the inlet to represent fully developed laminar flow using CEL; but I am getting some strange results. I am following the methodology here:

http://www.arc.vt.edu/ansys_help/cfx...DefiCrea5.html


However, in my case the inlet is on the YZ plane, not in the XY plane in the example. I just noticed in the last sentence it says

"In this equation, x and y are defined as directions 1 and 2 (X and Y for Cartesian coordinate frames) respectively, in the selected reference coordinate frame."

Do you think I have a problem with referencing the coordinate where the radius is measured from?

I have attached an image of the resulting velocity across the inlet. As you can see it is not what I was aiming for i.e. fastest velocity in the center and 0 velocity at the wall.

Does anyone have any idea why this is happening or know of an alternative approach I can use to solve this problem?

Thanks in advance!

[ghorrocks]

For this approach to work you need to have the origin at the centre of the inlet boundary, and in the XY plane. If this is not the case you will need to account for the different origin location.

[Opaque]

Definitely, there is a problem with the use of the "r" CEL variable in your setup. As explained in the tutorial you pointed out, "r" refers to the radius measures about the Z axis. In your setup, the proper axis for you expression is X, and the radius must be computed accordingly.

Alternatives:

1 - Create your own radius variable. For example, Myr = sqrt(y^2 + z^2).. Assuming the inlet of interest is centered respect to the coordinate frame; otherwise, adjust for the location of the coordinate frame respect to the inlet, i.e. sqrt( (y - y_o)^2 + (z - z_o)^2 )

2 - Change the option from Cartesian Components to Cylindrical Components, set Tangential, and Radial components to 0, and use your expression using Radius instead of "r" for you Axial component. Do not forget to set "x" for you Axis Definition. NOTE: Radius depends on the Axis Definition, while "r" is always about "z" axis.

3 - Create a new coordinate frame centered on the inlet with "local z" aligned with "global x". Now "r" is about the "local z" and it should work as you expect.

Hope the above helps,

[Allankey]

1 - Create your own radius variable. For example, Myr = sqrt(y^2 + z^2).. Assuming the inlet of interest is centered respect to the coordinate frame; otherwise, adjust for the location of the coordinate frame respect to the inlet, i.e. sqrt( (y - y_o)^2 + (z - z_o)^2 )

Thank you both so much for your replies, that makes a lot more sense now! Will give it a try.

Thanks again