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打开也没有问题
但使用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生成的网格却没有问题,说明求解器设置没有问题,所以错误出自哪里呢?
谢谢各位!
2019年9月23日 02:37
1/2
2019年10月24日 10:17