[sippycup]

I have been using the tutorials to become familiar with OpenFOAM and have been trying to recreate a cfdesign "turbulent flow over heated cylinder" case to compare different cfd packages. I know how to use FoamX, paraview etc, but my experience with flow is pretty minimal... I'm mostly the computer guy that's trying to get it to work.

I've created a satisfactory 2d (3d in foam but 2d for my purposes) mesh: it's a pipe cut in half lengthwise, with half a circle cut out where the "heated cylinder" would is. I therefore have six patches: left (needs to be a velocity inlet at 50 m/s), right (pressure outlet at 0 pa) down and cylinder (symmetryplanes), and up (wallFunctions).

I could not find an explicit "velocity inlet" in FoamX for the left side. From the example "cavity" case in the turbFoam directory I decided to make it a "wallFunction" and defined 50 for the x coordinate.

I calculated k and epsilon values for my particular case and input them in FoamX. They are 9.375 and 188.7 respectively. I think the delta_t required for the courant number is something like 0.003 s.

What I'm having problems with is the p field file. I get the following error: Create mesh for

time = 0
Reading field p
--> FOAM FATAL IO ERROR : attempt to read beyond EOF

I'm guessing that it would try to read beyond the end of file if it wasn't finding everything it needed, but that's a pretty feeble guess.

I tried to compare the cavity p and my file for differences, and the only thing I see is that I defined a p=0 for the right patch, all of the cavity patches are zeroGradient, and my internalField is nonuniform, whereas the cavity one is set to 0. I'm going to append my p file to the end of this message. Any help at all would be greatly appreciated.

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

// Field Dictionary

FoamFile
{
version 2.0;
format ascii;

root "/home/flowloop/OpenFOAM/flowloop-1.4/run/tutorials/turbFoam";
case "2.20";
instance "0";
local "";

class volScalarField;
object p;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


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

internalField nonuniform;

boundaryField
{
down
{
type symmetryPlane;
}

right
{
type fixedValue;
value uniform 0;
}

up
{
type zeroGradient;
}

left
{
type zeroGradient;
}

cylinder
{
type symmetryPlane;
}

defaultFaces
{
type empty;
}
}


// ************************************************** *********************** //