有关SHM和多相流求解后处理时的异常情况

[吴建民]

我使用多相流求解器求解浓度场时，发现在障碍物附近对网格进行加密后，会出现这样的异常，我是用的snappyHexMesh绘制的网格，请问有没有老师之前遇到过类似的问题。是网格的问题还是求解器中的设置问题?

[吴建民]

感觉就是不同Level网格过渡处会有问题，我对背景网格加密过，或者取消边界层，问题都没解决

[吴建民]

上传一下我的snappyhexMeshDict文件，请各位老师帮忙看看```
code_text

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

#includeEtc "caseDicts/mesh/generation/snappyHexMeshDict.cfg"

castellatedMesh true;
snap            true;
addLayers       true;

geometry
{
    cylinder
    {
        type triSurfaceMesh;
        file "cylinder.stl";
    }

    refinementBox
    {
        type searchableBox;
        min  (0.0 0.0 0.08);
        max  (1.2 0.2 0.22);
    }
};

castellatedMeshControls
{


    maxLocalCells 200000;

    maxGlobalCells 5000000;

    minRefinementCells 10;

    maxLoadUnbalance 0.10;

    nCellsBetweenLevels 6;
    
    features
    (
      { file  "cylinder.eMesh"; level 3; }
    );

    refinementSurfaces
    {
        cylinder
        {
            level (3 3);
            patchInfo { type wall; }
        }
    }

    refinementRegions
    {
        refinementBox
        {
            mode    inside;
            level   1;
        }
    }

    insidePoint (0.4 0.1 0.15);
}

snapControls
{
    nSmoothPatch 3;

    tolerance 2.0;

    nSolveIter 30;

    nRelaxIter 5;
    
    nFeatureSnapIter 10;
    
    explicitFeatureSnap    true;
    
    implicitFeatureSnap    false;
    
    multiRegionFeatureSnap false;
}

addLayersControls
{
    layers
    {
        "(cylinder_stl_face<stlunit=MM>|cylinder_stl_top<stlunit=MM>)"
        {
            nSurfaceLayers 10;
        }
    }

    relativeSizes       true;
    expansionRatio      1.05;
    finalLayerThickness 0.8;
    minThickness        0.5;
}

meshQualityControls
{}

writeFlags
(
    // scalarLevels
    // layerSets
    // layerFields
);

mergeTolerance 1e-6;

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