[Ahyar]

hello i'm trying to build a simulation with MRF on it. but when i tried to run the solution command, i've got this error



on the other hand, i've declared the cellZone on my snappyHexMeshDict

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

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

castellatedMesh true;
snap            true;
addLayers       false;

geometry
{
    inlet.stl
    {
        type triSurfaceMesh;
        name inlet;
    }
    tank.stl
    {
        type triSurfaceMesh;
        name tank;
    }
    impeller 
    {
        type triSurfaceMesh;
        file "impeller.stl";
    }
    MRF 
    {
        type searchableCylinder;
        point1 (0.0 0.0 -0.1);
        point2 (0.0 0.0 0.9);
        radius 0.4;
    }
}


castellatedMeshControls
{

    maxLocalCells 100000;

    maxGlobalCells 2000000;

    minRefinementCells 0;

    nCellsBetweenLevels 2;

    maxLoadUnbalance 0.10;

    features
    (
        {
            file "inlet.eMesh";
            level 0;
        }

        {
            file "tank.eMesh";
            level 0;
        }
        {
            file "impeller.eMesh";
            level 0;
        }

    );

    resolveFeatureAngle 30;

    refinementSurfaces
    {
        inlet
        {
            level (0 0);
        }
        tank
        {
            level (0 0);
        }
        impeller
        {
            level (0 0);
        }
        MRF 
        {
            level (1 1);
            cellZones cellMRFzone;
            faceZones faceMRFzone;
            cellZoneInside inside;
        }

    }

    refinementRegions
    {
        MRF
        {
            mode inside;
            levels ((1E15 1));
        }
        
    }

    locationInMesh (0.0 0.0 0.3); // Inside point

    allowFreeStandingZoneFaces true;
}

// Settings for the snapping.
snapControls
{

    nSmoothPatch 3;

    tolerance 1.0;

    nSolveIter 300;

    nRelaxIter 10;

    nFeatureSnapIter 10;

    implicitFeatureSnap true;

    explicitFeatureSnap false;

    multiRegionFeatureSnap true;
}

// Settings for the layer addition.
addLayersControls
{

    relativeSizes true;

    layers
    {
    }

    expansionRatio 1.0;

    finalLayerThickness 0.3;

    minThickness 0.25;

    nGrow 0;

    featureAngle 30;

    nRelaxIter 5;

    nSmoothSurfaceNormals 1;

    nSmoothNormals 3;

    nSmoothThickness 10;

    maxFaceThicknessRatio 0.5;

    maxThicknessToMedialRatio 0.3;

    minMedianAxisAngle 90;

    nBufferCellsNoExtrude 0;

    nLayerIter 50;

    nRelaxedIter 20;
}

meshQualityControls
{
    #include "meshQualityDict"

    relaxed
    {
        maxNonOrtho 75;
    }
}

writeFlags
(
);

mergeTolerance 1E-6;


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

and here's the MRFProperties if it matters

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

MRFImpeller
{
    cellZone    cellMRFzone;
    active      yes;

    // Fixed patches (by default they 'move' with the MRF zone)
    nonRotatingPatches ();

    origin    (0.0 0.0 0.0);
    axis      (0 1 0);
    omega     constant 12.56; // ~120 RPM
}

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

what should i do?
thanks in advance

[linnemann]

Hi

MRF 
{
     level (1 1);
     cellZones cellMRFzone;
     faceZones faceMRFzone;
     cellZoneInside inside;
}

Should be

MRF 
{
     level (1 1);
     cellZone cellMRFzone;
     faceZone faceMRFzone;
     cellZoneInside inside;
}

You can open the case in Paraview after meshing choose to select "Read zones" go to filter "Extract block" and choose the cellzone. That way you can inspect it has been created correctly. checkMesh command will also show if the cellZone has been created.

[Ahyar]

Originally Posted by linnemann

Hi

MRF 
{
     level (1 1);
     cellZones cellMRFzone;
     faceZones faceMRFzone;
     cellZoneInside inside;
}

Should be

MRF 
{
     level (1 1);
     cellZone cellMRFzone;
     faceZone faceMRFzone;
     cellZoneInside inside;
}

You can open the case in Paraview after meshing choose to select "Read zones" go to filter "Extract block" and choose the cellzone. That way you can inspect it has been created correctly. checkMesh command will also show if the cellZone has been created.

thank you so much, it really works