[anu.raghun]

I am trying to simulate laminar flow over a delta wing in OpenFoam. I keep getting this error with blockMesh though, could someone help me out?

--> FOAM FATAL ERROR:
face 0 in patch 1 does not have neighbour cell face: 4(4 3 5 4)

From function Foam::labelList Foam:olyMesh::facePatchFaceCells(const faceList&, const labelListList&, const faceListList&, Foam::label) const
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 118.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam:olyMesh::facePatchFaceCells(Foam::List<Foam ::face> const&, Foam::List<Foam::List<int> > const&, Foam::List<Foam::List<Foam::face> > const&, int) const at ??:?
#3 Foam:olyMesh::setTopology(Foam::List<Foam::cellS hape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::List<int>&, Foam::List<int>&, int&, int&, Foam::List<Foam::cell>&) at ??:?
#4 Foam:olyMesh:olyMesh(Foam::IOobject const&, Foam::Xfer<Foam::Field<Foam::Vector<double> > > const&, Foam::List<Foam::cellShape> const&, Foam::List<Foam::List<Foam::face> > const&, Foam::List<Foam::word> const&, Foam::PtrList<Foam::dictionary> const&, Foam::word const&, Foam::word const&, bool) at ??:?
#5 Foam::blockMesh::createTopology(Foam::IOdictionary const&, Foam::word const&) at ??:?
#6 Foam::blockMesh::blockMesh(Foam::IOdictionary const&, Foam::word const&) at ??:?
#7 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/blockMesh"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/blockMesh"
Aborted (core dumped)

This is what my blockMeshDict looks like:

convertToMeters 0.1;

vertices
(
(0 0 0) //0
(-0.5 0 -2) //1
(0.5 0 -2) //2
(0 -2 0) //3
(-0.5 -2 -2) //4
(0.5 -2 -2) //5
);

blocks
(
hex (0 1 2 0 3 4 5 3) (50 1 200) simpleGrading (0.2 1 1)
);

edges
(
);

boundary
(
axis
{
type empty;
faces
(
(0 3 3 0)
);
}
plate
{
type plate;
faces
(
(4 3 5 4)
);
}
inlet1
{
type patch;
faces
(
(4 1 0 3)
);
}
inlet2
{
type patch;
faces
(
(3 0 2 5)
);
}
top
{
type symmetry;
faces
(
(1 2 0 1)
);
}
outlet
{
type patch;
faces
(
(1 4 5 2)
);
}
);

mergePatchPairs
(
);

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

I appreciate any help!

[Antimony]

Hi,

It would seem that the face definitions for the "top" and "plate" patches have to follow a different order.

For plate, change (4 3 5 4) to (3 4 5 3) and for top, change (1 2 0 1) to (0 1 2 0).

As to why this is, it might be the case that when you have a wedge, OF is fine when the repeated vertex in the hex definition appears twice in a boundary patch definition.

Hope this helps.

Cheers,
Antimony

[anu.raghun]

Thank you so much! That worked perfectly.