snappyHexMesh画网格

HIT_lly

使用snappyHexMesh画网格时，画出网格时没有什么问题但是在计算的时候会报错。有没有大佬知道这种情况该怎么解决

FOAM Warning : 
    From void Foam::cutFaceAdvect::cutPoints(Foam::label, Foam::scalar, Foam::DynamicList<Foam::Vector<double> >&)
    in file cellCuts/cutFace/cutFaceAdvect.C at line 927
    cutPoints = 4((0 0.0449362922089 0.385062422146) (0 0.0436853188202 0.385062441697) (0.01 0.0436853188202 0.385062441697) (0.01 0.0449362922089 0.385062422146)) for pts = 4((0 0.0449362922089 0.385062422146) (0 0.0436853188202 0.385062441697) (0.01 0.0436853188202 0.385062441697) (0.01 0.0449362922089 0.385062422146)), f - f0 = 4(23405 23403 38935 38937) and f0 = 2.22618323362e-13

李东岳

这是个warning。看起来好像是isoAdvector那面的算法？这个算法我不是太熟悉。回答也是乏力了。

HIT_lly

@李东岳谢谢东岳老师，这个我用interFoam也会报错，应该是snappyHexMesh切分网格的时候产生的问题。

李东岳

interFoam不需要cutFaceAdvect.C，你interFoam用的界面压缩的格式么？

HIT_lly

@李东岳东岳老师，上面那个报错是因为使用了isoAdvector算法。
使用interFoam进行计算也会进行报错，报错是这样的

FOAM Warning : 
    From Foam::triFace Foam::tetIndices::faceTriIs(const Foam::polyMesh&, bool) const
    in file meshes/polyMesh/polyMeshTetDecomposition/tetIndicesI.H at line 70
    No base point for face 7847, 4(5577 25968 25969 5578), produces a valid tet decomposition.

李东岳

你的interfoam的fvScheme怎么设置的

HIT_lly

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

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
default         Gauss upwind;
    div(rhoPhi,U)   Gauss upwind;
    div(U)          Gauss linear;
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss linear;
     div(phi,thermo:rho.water) Gauss upwind;
    div(phi,thermo:rho.air) Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
       div(rhoPhi,T)  Gauss upwind;
    div(rhoPhi,K)  Gauss upwind;

        div(phi,p)      Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(phi,omega) Gauss upwind;
    div(phi,U)      Gauss linearUpwind limited;

    energy          Gauss linearUpwind limited;
    div(phi,e)      $energy;
    div(phi,K)      $energy;
    div(phi,Ekp)    $energy;

    turbulence      Gauss linearUpwind limited;
    div(phi,k)      $turbulence;
    div(phi,omega)  $turbulence;
    div(rhoPhi,omega)  $turbulence;
div(rhoPhi,k)  $turbulence;


    div(phiv,p)     Gauss upwind;
    div((phi|interpolate(rho)),p) Gauss upwind;

}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method          meshWave;
}
oversetInterpolation
{
    method              inverseDistance;
    holeLayers          6;
    useLayer            3;
}

oversetInterpolationSuppressed
{
    grad(p_rgh);
    surfaceIntegrate(phiHbyA);
}

fluxRequired
{
    default         no;
    p_rgh;
    pcorr;
    alpha.water;
}


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

HIT_lly

@李东岳感谢东岳老师的回复，我的fvScheme是这样设置的，这用的是重叠网格技术。开始的isoadvector是我耦合的重叠网格算法，然后后面的compressibleinterFoam是openfoam 中自带的算法。

HIT_lly

之前snappyHexMesh画的3D网格进行重叠网格计算没有问题。后来为了节省计算成本用的2D网格，用的snappyHexMesh画的网格，然后使用extrudeMesh将3D网格变成2D网格后就出现问题了

CFD中文网

snappyHexMesh画网格