Problem getting Parabolic profile in pipe flow
 

Dear all
I was trying to get a parabolic fully-developed profile for a flow through a simple tube, meshing an axisymmetric domain. Even though I was using the solver for reacting flow- the laminarSimpleSMOKE solver (an incompressible laminar flow solver with chemical reactions), as a first check I am not solving the species equations, just solving for the flow. I guess I am doing some mess in my pressure and velocity boundaries and thats why couldnt get the profile. The profile I get is attached. I boundary conditions I use are-

For P-
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 101325;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 101325;
}
axisym_f1
{
type wedge;
}
axisym_r1
{
type wedge;
}
outerwall
{
type zeroGradient;
}
}

For U-
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Umean 0.138475;

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform ($Umean 0 0);
}
outlet
{
type zeroGradient;
}
axisym_f1
{
type wedge;
}
axisym_r1
{
type wedge;
}
outerwall
{
type fixedValue;
value uniform (0 0 0);
}
}

For T-
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 298;


boundaryField
{
inlet
{
type fixedValue;
value uniform 298;
}
outlet
{
type zeroGradient;
}
axisym_f1
{
type wedge;
}
axisym_r1
{
type wedge;
}
outerwall
{
type fixedValue;
value uniform 300;
}
}

Could anyone please tell me if you experience similar problems? Please share your suggestions.
Thanks