[fantasteak]

Hi foamers,

I'm trying to set up a mesh for a combustion model and have a problem in placing a cricular inlet on the front wall of the rectangular burner. After searching for a while I think topoSet and createPatch are the utilities I need, but I don't get it work.
So how can I define a circular geometry in topoSet? After cylinderToCell I'm very baffled.
Maybe someone knows how to slove this problem?

Thank you very much in advance!

Here's my blockMeshDict:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 1;

vertices
(
    (0 0 0)		//0, Ursprung, FRONT
    (0.2 0 0)		//1, x_1-Wert, FRONT
    (0.2 0.2 0)		//2, y_1-Wert, FRONT
    (0 0.2 0)		//3, z_1-Wert, FRONT
    (0 0 2.5)		//4, Ursprung, BACK
    (0.2 0 2.5)		//5, x_2-Wert, BACK
    (0.2 0.2 2.5)	//6, y_2-Wert, BACK
    (0 0.2 2.5)		//7, z_2-Wert, BACK
    
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (35 35 251) simpleGrading (((1 0.23 1) (1 0.54 1) (1 0.23 1)) ((1 0.23 1) (1 0.54 1) (1 0.23 1)) 1)
   
);

edges
(
);

boundary
(

    front
    {
        type patch;
        faces
        (
            (0 3 2 1)
	);
	
    }

    back
    {
	type patch;
	faces
	(
	    (4 5 6 7)
	);
    }

   Seiten
    {
	type wall;
	faces
	(
	   (0 1 5 4)
	   (3 2 6 7)
	   (1 5 6 2)
	   (0 4 7 3)
	);
    }
);


mergePatchPairs
(
);

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

topoSetDict:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      topoSetDict;
}

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

actions
(
    // front
{
        name    roundCells;
        type    cellSet;
        action  new;
         source cylinderToCell;
    sourceInfo
    {
       p1       (0.088 0.088 0);   // start point on cylinder axis
       p2       (0.112 0.112 0);   // end point on cylinder axis
       radius   0.012;
    }
}

{
        name    frontSet;
        type    faceSet;
        action  new;
        source  cellToFace;
        sourceInfo
        {
            set     roundCells;
            option  all;
        }
    }

{
        name    frontSet;
        type    faceSet;
        action  subset;
        source  boxToFace;
        sourceInfo
        {
            box (0.088 0.088 0)(0.112 0.112 0);
        }
}
	
{
        name    roundCells1;
        type    cellSet;
        action  new;
        source  cylinderToCell;
        sourceInfo
        {
           p1       (0.088 0.088 0);  // start point on cylinder axis
           p2       (0.112 0.112 0);   // end point on cylinder axis
           radius   0.012;
        }
}

{
        name    frontSet1;
        type    faceSet;
        action  new;
        source  cellToFace;
        sourceInfo
        {
            set     roundCells1;
            option  all;
        }
    } 
	
{
        name    frontSet1;
        type    faceSet;
        action  subset;
        source  boxToFace;
        sourceInfo
        {
            box (0.088 0.088 0)(0.112 0.112 0);
        }
}
   
{
        name    front;
        type    faceZoneSet;
        action  new;
        source  setToFaceZone;
        sourceInfo
        {
            faceSet frontSet;
        }
}

{
        name    front1;
        type    faceZoneSet;
        action  new;
        source  setToFaceZone;
        sourceInfo
        {
            faceSet frontSet1;
        }
    }
    
);


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

createPatchDict:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      createPatchDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Do a synchronisation of coupled points after creation of any patches.
// Note: this does not work with points that are on multiple coupled patches
//       with transformations (i.e. cyclics).
pointSync true;

// Patches to create.
patches
(
{
	name inletA;
	patchInfo
	{
	type cyclicAMI;
	neighbourPatch inletB;
	coupleGroup inlet;
	}
	constructFrom set;
	set front;
}

{
	name inletB;
	patchInfo
	{
	type cyclicAMI;
	neighbourPatch inletA;
	}
	constructFrom set;
	set front1;
}

{
	name outlet;
	patchInfo
	{
	type empty;
	}
	constructFrom patches;
	patches (back);
}
);

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

[avila.vc]

Hi, it's most likely not the best solution, but if it takes too long for someone to give you the right answer and for you to manage to do your circular inlet, I would recommend to create your mesh with other software and then export to OpenFOAM. I usually do that with ICEM (Ansys) or Gmsh(free software), I think they are easier since they have a graphical interface to set up geometry and meshes.

Best luck

[Antimony]

Hi,

You are creating patches from faceZoneSets instead of faceSets. Have you tried it with the faceSets?

Cheers,
Antimony

[fantasteak]

Hi,
thanks for your answers.
Yes I tried with the faceSets, but the result remains the same.
And when deactivating the subsets I get this error:

--> FOAM FATAL ERROR:
Face 1193 specified in set frontSet is not an external face of the mesh.
This application can only repatch existing boundary faces.

    From function createPtch.exe
    in file createPatch.C at line 749.

FOAM exiting

Using other software is unfortunately not really an option...