|
[Sponsors] |
createPatch produces error "faces with negative volume decomposition" |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 25, 2013, 05:39 |
createPatch produces error "faces with negative volume decomposition"
|
#1 |
New Member
Danesh S
Join Date: Jul 2013
Location: Bochum, Germany
Posts: 27
Rep Power: 12 |
Hello All,
I just rebuild a closed cylindrical mesh into a cyclic cylindrical mesh. For the purpose of reducing the number etc. I used createPatch, respective a createPatchdict. Everything worked fine so far, but now I get this error when trying to run pisoFoam: Code:
Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type LESModel Selecting LES turbulence model Smagorinsky #0 Foam::error::printStack(Foam::Ostream&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib64/libc.so.6" #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::operator/(double const&, Foam::tmp<Foam::Field<double> > const&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #5 Foam::surfaceInterpolation::makeDeltaCoeffs() const in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #6 Foam::surfaceInterpolation::deltaCoeffs() const in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" #7 Foam::incompressible::LESModel::LESModel(Foam::word const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so" #8 Foam::incompressible::LESModels::Smagorinsky::Smagorinsky(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so" #9 Foam::incompressible::LESModel::adddictionaryConstructorToTable<Foam::incompressible::LESModels::Smagorinsky>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so" #10 Foam::incompressible::LESModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so" #11 Foam::incompressible::turbulenceModel::addturbulenceModelConstructorToTable<Foam::incompressible::LESModel>::NewturbulenceModel(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so" #12 Foam::incompressible::turbulenceModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so" #13 in "/opt/net/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/pisoFoam" #14 __libc_start_main in "/lib64/libc.so.6" #15 at /home/abuild/rpmbuild/BUILD/glibc-2.15/csu/../sysdeps/x86_64/elf/start.S:116 Gleitkomma-Ausnahme Code:
Create polyMesh for time = 0 Time = 0 Mesh stats points: 456000 faces: 1319760 internal faces: 1272240 cells: 432000 boundary patches: 6 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 432000 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. Topological cell zip-up check OK. Face-face connectivity OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface topology Bounding box aeussereWand 18000 18240 ok (non-closed singly connected) (-1 0 -1) (1 5 1) innereWand 18000 18240 ok (non-closed singly connected) (-0.5 0 -0.5) (0.5 5 0.5) cyc_half0 2880 3025 ok (non-closed singly connected) (-1 0 0) (1 0 1) cyc_half1 2880 3025 ok (non-closed singly connected) (-1 0 -1) (1 0 0) cyc_half10 2880 3025 ok (non-closed singly connected) (-1 5 0) (1 5 1) cyc_half11 2880 3025 ok (non-closed singly connected) (-1 5 -1) (1 5 0) Checking geometry... Overall domain bounding box (-1 0 -1) (1 5 1) Mesh (non-empty, non-wedge) directions (1 1 1) Mesh (non-empty) directions (1 1 1) Boundary openness (1.00823e-17 -8.79755e-18 -1.68404e-16) OK. Max cell openness = 2.53899e-16 OK. Max aspect ratio = 5.92719 OK. Minimum face area = 0.000278357. Maximum face area = 0.00174528. Face area magnitudes OK. Min volume = 1.85571e-05. Max volume = 3.59782e-05. Total volume = 11.7796. Cell volumes OK. Mesh non-orthogonality Max: 0.000138361 average: 2.95756e-06 Non-orthogonality check OK. Face pyramids OK. Max skewness = 0.0127689 OK. Coupled point location match (average 0) OK. ***Error in face tets: 11520 faces with low quality or negative volume decomposition tets. <<Writing 11520 faces with low quality or negative volume decomposition tets to set lowQualityTetFaces Min/max edge length = 0.0130896 0.0666667 OK. All angles in faces OK. Face flatness (1 = flat, 0 = butterfly) : average = 1 min = 1 All face flatness OK. Cell determinant (wellposedness) : minimum: 0.260484 average: 3.33264 Cell determinant check OK. Concave cell check OK. Failed 1 mesh checks. Do you have any suggestions concerning this? Is it a bug? Why do I have these negative cells now, since I did not change the order of vertices for the faces? Help is indeed welcome Cheers, Danesh Edit: "Gleitkomma-Ausnahme" is german for floating point exception :-D |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] mesh airfoil NACA0012 | anand_30 | OpenFOAM Meshing & Mesh Conversion | 13 | March 7, 2022 17:22 |
how to set periodic boundary conditions | Ganesh | FLUENT | 15 | November 18, 2020 06:09 |
[snappyHexMesh] How to define to right point for locationInMesh | Mirage12 | OpenFOAM Meshing & Mesh Conversion | 7 | March 13, 2016 14:07 |
[Commercial meshers] negative volume in gambit | sofie1 | OpenFOAM Meshing & Mesh Conversion | 6 | January 18, 2012 02:34 |
remeshing due to negative volume error | Doginal | CFX | 1 | August 21, 2011 21:50 |