pointwise输出二维网格到openfoam无法计算



  • 大家好,我用pointwise输出了一个翼型的网格,使用checkMesh检查没有错误

    Mesh stats
        points:           28336
        internal points:  0
        faces:            56210
        internal faces:   27874
        cells:            14014
        faces per cell:   6
        boundary patches: 3
        point zones:      0
        face zones:       0
        cell zones:       0
    
    Overall number of cells of each type:
        hexahedra:     14014
        prisms:        0
        wedges:        0
        pyramids:      0
        tet wedges:    0
        tetrahedra:    0
        polyhedra:     0
    
    Checking topology...
        Boundary definition OK.
        Cell to face addressing OK.
        Point usage OK.
        Upper triangular ordering OK.
        Face vertices OK.
        Number of regions: 1 (OK).
    
    Checking patch topology for multiply connected surfaces...
        Patch               Faces    Points   Surface topology                  
        airFoil2D           154      308      ok (non-closed singly connected)  
        farfield            154      308      ok (non-closed singly connected)  
        frontAndBack        28028    28336    ok (non-closed singly connected)  
    
    Checking geometry...
        Overall domain bounding box (-4.99867 -5.67032 0) (6.78888 5.66637 3)
        Mesh (non-empty, non-wedge) directions (1 1 0)
        Mesh (non-empty) directions (1 1 0)
        All edges aligned with or perpendicular to non-empty directions.
        Boundary openness (-1.85094e-17 2.46792e-17 -2.22841e-16) OK.
        Max cell openness = 4.12267e-15 OK.
        Max aspect ratio = 170.167 OK.
        Minimum face area = 4.43442e-08. Maximum face area = 1.67744.  Face area magnitudes OK.
        Min volume = 1.33032e-07. Max volume = 0.378026.  Total volume = 313.461.  Cell volumes OK.
        Mesh non-orthogonality Max: 17.8429 average: 1.07748
        Non-orthogonality check OK.
        Face pyramids OK.
        Max skewness = 1.9585 OK.
        Coupled point location match (average 0) OK.
    
    Mesh OK.
    
    End
    
    

    我用paraFoam打开也没有问题
    naca.png
    但使用simpleFoam运算就会出现问题:

    smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 156.2, No Iterations 1000
    smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 4.72215, No Iterations 1000
    GAMG:  Solving for p, Initial residual = 1, Final residual = 2.66928e+147, No Iterations 1000
    time step continuity errors : sum local = 1.37947e+146, global = 9.17194e+129, cumulative = 9.17194e+129
    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::sigFpe::sigHandler(int) at ??:?
    #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
    #4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
    #5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
    #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
    #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
    #8  Foam::fvMatrix<double>::solve() at ??:?
    #9  Foam::incompressible::RASModels::SpalartAllmaras::correct() at ??:?
    #10  ? at ??:?
    #11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #12  ? at ??:?
    Floating point exception (core dumped)
    

    我用Construct2D或者snappyHexmesh生成的网格却没有问题,说明求解器设置没有问题,所以错误出自哪里呢?
    谢谢各位!


Log in to reply