[naveen k s]

Dear foamers,


I am trying to simulate the 2D plane poisuelle flow. I am using icoFoam to simulate the solution. These are my BC's in p and U field. The parallel plates are 2m long, 0.1m in height. Pressure at the inlet is 101.325 and outlet is zero. density is 1000 kg/m^3, kinematic viscosity is 0.01 m^2/s. For these conditions i am ending up with Umax of 12 m/s. Which is twice of the true value. and one more strange thing is with increase in simulation time Umax is also increasing. Anyone there to help me with this?



FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
top
{
type zeroGradient;
}

bottom
{
type zeroGradient;
}

inlet
{
type fixedValue;
value uniform 101.325;
}

outlet
{
type fixedValue;
value uniform 0;
}

frontAndBack
{
type empty;
}
}


U Field


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

internalField uniform (0 0 0);

boundaryField
{
top
{
type fixedValue;
value uniform (0 0 0);
}

bottom
{
type fixedValue;
value uniform (0 0 0);
}

inlet
{
type zeroGradient;
}

outlet
{
type zeroGradient;
}

frontAndBack
{
type empty;
}
}


blockMeshDict


FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}

ConvertToMeters 0.1;

vertices
(
(0 0 0)
(20 0 0)
(20 1 0)
(0 1 0)
(0 0 0.1)
(20 0 0.1)
(20 1 0.1)
(0 1 0.1)
);

blocks
(
hex ( 0 1 2 3 4 5 6 7 ) ( 20 20 1 ) simpleGrading ( 1 1 1 )
);

edges
(
);

boundary
(
top
{
type wall;
faces
(
(3 7 6 2)
);
}
bottom
{
type wall;
faces
(
(1 5 4 0)
);
}
inlet
{
type patch;
faces
(
(0 4 7 3)
);
}
outlet
{
type patch;
faces
(
(2 6 5 1)
);
}
frontAndBack
{
type empty;
faces
(
(0 3 2 1)
(4 5 6 7)
);
}
);

mergePatchPairs
(
);