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生成的网格却没有问题，说明求解器设置没有问题，所以错误出自哪里呢？
谢谢各位！ -
最后解决了吗。我现在还是会遇到这个问题，改变了边界条件还是这样的计算结果。
Create time Create mesh for time = 0 SIMPLE: convergence criteria field p tolerance 1e-05 field U tolerance 1e-05 field nuTilda tolerance 1e-05 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model SpalartAllmaras Selecting patchDistMethod meshWave RAS { RASModel SpalartAllmaras; turbulence on; printCoeffs on; sigmaNut 0.666667; kappa 0.41; Cb1 0.1355; Cb2 0.622; Cw2 0.3; Cw3 2; Cv1 7.1; Cs 0.3; } No MRF models present No finite volume options present Starting time loop momentumError momErr: turbulenceFields turbulenceFields1: storing fields: turbulenceProperties:nut turbulenceProperties:nuEff turbulenceProperties:k turbulenceProperties:R turbulenceProperties:I turbulenceProperties:L turbulenceProperties:devReff turbulenceProperties:epsilon Time = 0.01 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00113906, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.402392, No Iterations 1000 #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::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:? #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:? #5 Foam::GAMGSolver::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>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:? #8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:? #9 Foam::fvMatrix<double>::solve() in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/simpleFoam #10 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/simpleFoam #11 ? in /lib/x86_64-linux-gnu/libc.so.6 #12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6 #13 ? in ~/OpenFOAM/OpenFOAM-v2206/platforms/linux64GccDPInt32Opt/bin/simpleFoam 浮点异常 (核心已转储)
