Dear All,

In a stenosed pipe flow problem, I specify parabolic velocity profile u=2*(1−(r^2)), where r=sqrt(y^2+z^2) as seen below.

<BOUNDARYREGIONS>
            <B ID="0"> C[7] </B>  <!-- Inlet -->
            <B ID="1"> C[9] </B>  <!-- Outlet -->
            <B ID="2"> C[8] </B>  <!-- Wall -->
        </BOUNDARYREGIONS>

        <BOUNDARYCONDITIONS>
            <REGION REF="0">
                <D VAR="u" VALUE="2*(1-(rad(y,z)^2))" />
                <D VAR="v" VALUE="0" />
                <D VAR="w" VALUE="0" />
                <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
            </REGION>
            <REGION REF="1">
                <N VAR="u" VALUE="0" />
                <N VAR="v" VALUE="0" />
                <N VAR="w" VALUE="0" />
                <D VAR="p" VALUE="0" />
            </REGION>
            <REGION REF="2">
                <D VAR="u" VALUE="0" />
                <D VAR="v" VALUE="0" />
                <D VAR="w" VALUE="0" />
                <N VAR="p" USERDEFINEDTYPE="H" VALUE="0" />
            </REGION>
        </BOUNDARYCONDITIONS>

However, since u=0 at wall of the pipe due to no slip boundary condition definition, the calculated value of parabolic velocity profile at the wall (u = 1.5 ) can not be achieved. Is there any way to overcome this problem?

Regards,
Kamil