[Muerio]

Hello,

I'm currently trying to run refineMesh in parallel after completing the castellated mesh step in snappyHexMesh.
My goal is to refine the mesh in the x direction only.

My refineMeshDict_xRefine looks like this:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      refineMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

set xRefine;

globalCoeffs
{
    tan1 (1 0 0);
    tan2 (0 1 0);
}

(
    tan1
);

useHexTopology  true;

geometricCut    false;

writeMesh       false;

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

When running

runParallel refineMesh -overwrite -noFunctionObjects -dict system/refineMeshDict_xRefine

it results in errors like this one for every processor:

[1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] Cut (8143888 2468835) on face (0.055919118 0.18161516 0.0017083333) of coupled patch procBoundary1to0 is not consistent with coupled cut (1860044 8143888)
[1] 
[1]     From function void Foam::cellCuts::syncProc()
[1]     in file meshCut/cellCuts/cellCuts.C at line 256.
[1] 
FOAM parallel run exiting

Changing the switches for useHexTopology and geometricCut results in the same errors.

Does anyone know how to prevent this error?

[huazhongdk@gmail.com]

Originally Posted by Muerio

Hello,

I'm currently trying to run refineMesh in parallel after completing the castellated mesh step in snappyHexMesh.
My goal is to refine the mesh in the x direction only.

My refineMeshDict_xRefine looks like this:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      refineMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

set xRefine;

globalCoeffs
{
    tan1 (1 0 0);
    tan2 (0 1 0);
}

(
    tan1
);

useHexTopology  true;

geometricCut    false;

writeMesh       false;

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

When running

runParallel refineMesh -overwrite -noFunctionObjects -dict system/refineMeshDict_xRefine

it results in errors like this one for every processor:

[1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] Cut (8143888 2468835) on face (0.055919118 0.18161516 0.0017083333) of coupled patch procBoundary1to0 is not consistent with coupled cut (1860044 8143888)
[1] 
[1]     From function void Foam::cellCuts::syncProc()
[1]     in file meshCut/cellCuts/cellCuts.C at line 256.
[1] 
FOAM parallel run exiting

Changing the switches for useHexTopology and geometricCut results in the same errors.

Does anyone know how to prevent this error?

hi Muerio,
you can try to run snappyHexMesh with only setting castellatedMesh as true, then use refineMesh to refine the mesh. After that, remove the cellLevel and pointLevel files in constant/polyMesh, and run snappyHexMesh again with snap as true. I solve the similar problem by doing these steps, hope it will help you.

1. blockMesh -dict ./system/blockMeshDict
2. snappyHexMesh -dict ./system/snappyHexMeshDict -overwrite
3. topoSet -dict ./system/topoSetDict
4. refineMesh -dict ./system/refineMeshDict -overwrite
5. cd ./constant/polyMesh && rm cellLevel && rm pointLevel
6. snappyHexMesh -dict ./system/snappyHexMeshDict.snappy -overwrite