[carserg92]

Hello everyone,

I made a plane in OpenFoam and I want to apply rhoCentralFoam as solver. But when I do this it happens the next error:

#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::heThermo<Foam:: psiThermo, Foam:: pureMixture<Foam::constTransport<Foam::species::th ermo<Foam::hConstThermo<Foam:: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#4 Foam:: psiThermo::addfvMeshConstructorToTable<Foam::hePsi Thermo<Foam:: psiThermo, Foam:: pureMixture<Foam::constTransport<Foam::species::th ermo<Foam::hConstThermo<Foam:: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#5 Foam::autoPtr<Foam:: psiThermo> Foam::basicThermo::New<Foam:: psiThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
#6 Foam:: psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7
at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
at ??:?


Could anyone help me please? It's for a proyect and I need it. I will be very thanfuly for it.

Thank you.

I add the data of the proyect. Thank you very much

Attachment 39585

[T.D.]

Hi try to check your mesh first. It failed 7 checks:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 682480
faces: 1961600
internal faces: 1878400
cells: 640000
faces per cell: 6
boundary patches: 11
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 639988
prisms: 0
wedges: 12
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
*Number of regions: 6
The mesh has multiple regions which are not connected by any face.
<<Writing region information to "0/cellToRegion"
<<Writing region 0 with 192000 cells to cellSet region0
<<Writing region 1 with 64000 cells to cellSet region1
<<Writing region 2 with 64000 cells to cellSet region2
<<Writing region 3 with 64000 cells to cellSet region3
<<Writing region 4 with 128000 cells to cellSet region4
<<Writing region 5 with 128000 cells to cellSet region5

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
cilindro 6400 6560 ok (non-closed singly connected)
estabilizadorv 6400 6560 ok (non-closed singly connected)
cola 8000 8081 ok (non-closed singly connected)
alaizq 8000 8081 ok (non-closed singly connected)
alader 8000 8081 ok (non-closed singly connected)
estabhorizq 9600 9602 ok (closed singly connected)
estabhorder 9600 9602 ok (closed singly connected)
sharkletsizq 8000 8081 ok (non-closed singly connected)
sharkletsder 8000 8075 ok (non-closed singly connected)
punta 8000 8081 ok (non-closed singly connected)
defaultFaces 3200 3362 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-14.36 0 -1.37) (14.36 38.65 7.81)
Mesh (non-empty, non-wedge) directions (0 0 0)
Mesh (non-empty) directions (0 0 0)
***Number of edges not aligned with or perpendicular to non-empty directions: 1592360
<<Writing 682480 points on non-aligned edges to set nonAlignedEdges
***Boundary openness (4.27990244461409e-16 -0.0216521575582463 -6.15926392209936e-16) possible hole in boundary description.
***Open cells found, max cell openness: 1, number of open cells 1600
<<Writing 1600 non closed cells to set nonClosedCells
Minimum face area = 3.0577969933228e-09. Maximum face area = 0.0514085168221845. Face area magnitudes OK.
***Zero or negative cell volume detected. Minimum negative volume: -0.00272489720430368, Number of negative volume cells: 325808
<<Writing 325808 zero volume cells to set zeroVolumeCells
Mesh non-orthogonality Max: 180 average: 91.4114929588283
*Number of severely non-orthogonal (> 70 degrees) faces: 25074.
***Number of non-orthogonality errors: 956857.
<<Writing 981931 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 1956722 faces are incorrectly oriented.
<<Writing 1000364 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 381.743726683058, 212 highly skew faces detected which may impair the quality of the results
<<Writing 212 skew faces to set skewFaces
Coupled point location match (average 0) OK.

Failed 7 mesh checks.

Originally Posted by carserg92

Hello everyone,

I made a plane in OpenFoam and I want to apply rhoCentralFoam as solver. But when I do this it happens the next error:

#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::heThermo<Foam:: psiThermo, Foam:: pureMixture<Foam::constTransport<Foam::species::th ermo<Foam::hConstThermo<Foam:: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#4 Foam:: psiThermo::addfvMeshConstructorToTable<Foam::hePsi Thermo<Foam:: psiThermo, Foam:: pureMixture<Foam::constTransport<Foam::species::th ermo<Foam::hConstThermo<Foam:: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#5 Foam::autoPtr<Foam:: psiThermo> Foam::basicThermo::New<Foam:: psiThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
#6 Foam:: psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7
at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
at ??:?


Could anyone help me please? It's for a proyect and I need it. I will be very thanfuly for it.

Thank you.

I add the data of the proyect. Thank you very much

Attachment 39585

[carserg92]

hi! thanks for the reply. I minimize the problem to a cylinder and a cone together. I do the checkMesh and everything is ok but the next error continuea appearing when I apply the rhoCentralFoam:

#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::heThermo<Foam:: psiThermo, Foam:: pureMixture<Foam::constTransport<Foam::species::th ermo<Foam::hConstThermo<Foam:: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#4 Foam:: psiThermo::addfvMeshConstructorToTable<Foam::hePsi Thermo<Foam:: psiThermo, Foam:: pureMixture<Foam::constTransport<Foam::species::th ermo<Foam::hConstThermo<Foam:: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#5 Foam::autoPtr<Foam:: psiThermo> Foam::basicThermo::New<Foam:: psiThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
#6 Foam:: psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7
at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
at ??:?

Could you help me please? Thank you so much!! I add the new mesh Proyecto2.zip

[T.D.]

Hi ,

I checked your attached final file. Your problem is due to the bad definition of pressure boundary condition (BC). If you change only the pressure BC at the outlet to:
outlet
{
type zeroGradient;
}

or to

outlet
{
type fixedValue;
value uniform 1;
}
The simulations run fine.

So, Try to pose well the Physics for your Boundary Conditions.

Good Luck !
Regards,

T.D.

Originally Posted by carserg92

hi! thanks for the reply. I minimize the problem to a cylinder and a cone together. I do the checkMesh and everything is ok but the next error continuea appearing when I apply the rhoCentralFoam:

#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::heThermo<Foam:: psiThermo, Foam:: pureMixture<Foam::constTransport<Foam::species::th ermo<Foam::hConstThermo<Foam:: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#4 Foam:: psiThermo::addfvMeshConstructorToTable<Foam::hePsi Thermo<Foam:: psiThermo, Foam:: pureMixture<Foam::constTransport<Foam::species::th ermo<Foam::hConstThermo<Foam:: perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#5 Foam::autoPtr<Foam:: psiThermo> Foam::basicThermo::New<Foam:: psiThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
#6 Foam:: psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7
at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
at ??:?

Could you help me please? Thank you so much!! I add the new mesh Attachment 39643