[ddemiroguzhan]

Dear All,
I am trying to make an analysis using overset functionality of OpenFOAM v17.06. 1 DoF is OK. 2 DoF is OK. But whenever I try using 2 DoF + angular rotation on the third axis (by the way the problem is 2D) the solution crashes. Log file of checkMesh and the error that OpenFOAM gives are attached. I know that the information is a bit insufficient but if anyone willing to go into the details, i can provide all the necessary data. By the way the solver is overPimpleDyMFoam

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1706 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : v1706
Arch : "LSB;label=32;scalar=64"
Exec : checkMesh
Date : Dec 07 2017
Time : 16:10:02
Host : "xxx"
PID : 179977
Case : /xxx/xxx/xxx/xxx/xxx/xxx/3dof_block_overset_mesh_3_rotation
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Time = 0

Mesh stats
points: 14578
internal points: 0
faces: 34459
internal faces: 15677
cells: 9057
faces per cell: 5.535607817
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 1

Overall number of cells of each type:
hexahedra: 4851
prisms: 4206
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: 2
The mesh has multiple regions which are not connected by any face.
<<Writing region information to "0/cellToRegion"
<<Writing region 0 with 4206 cells to cellSet region0
<<Writing region 1 with 4851 cells to cellSet region1

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
backandfront 18114 14578 ok (non-closed singly connected)
body 236 472 ok (non-closed singly connected)
farfield 198 400 ok (non-closed singly connected)
inlet 49 100 ok (non-closed singly connected)
outlet 49 100 ok (non-closed singly connected)
overset 136 272 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-20 -10 0) (20 10 1.000002775)
Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
Mesh has 2 solution (non-empty) directions (1 1 0)
***Number of edges not aligned with or perpendicular to non-empty directions: 1006
<<Writing 548 points on non-aligned edges to set nonAlignedEdges
Boundary openness (-1.698701896e-18 -1.941373595e-18 2.866559449e-18) OK.
Max cell openness = 2.190921877e-16 OK.
Max aspect ratio = 2.632205077 OK.
Minimum face area = 0.0005404535113. Maximum face area = 0.408163266. Face area magnitudes OK.
Min volume = 0.0005404535114. Max volume = 0.1649144533. Total volume = 834. Cell volumes OK.
Mesh non-orthogonality Max: 22.6836599 average: 3.265519369
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.4539633208 OK.
Coupled point location match (average 0) OK.

Failed 1 mesh checks.

End

Also the error that I have is:

Centre of rotation: (-744.8031209 -732.3031209 0)
Centre of mass: (-744.8031209 -732.3031209 0)
Orientation: (0.7668098028 0 0.6418743851 0 1 0 -0.6418743851 0 0.7668098028)
Linear velocity: (-1063020.895 -1063020.895 -0)
Angular velocity: (0 -12768045.14 0)
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/usr/lib64/libc.so.6"
#3 Foam::surfaceInterpolation::makeWeights() const at ??:?
#4 Foam::surfaceInterpolation::weights() const at ??:?
#5 Foam::surfaceInterpolation::makeNonOrthDeltaCoeffs () const at ??:?
#6 Foam::surfaceInterpolation::nonOrthDeltaCoeffs() const at ??:?
#7 Foam::fv::correctedSnGrad<double>::deltaCoeffs(Foa m::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#8 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#9 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double>(Foam::GeometricField< double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#10 Foam:atchDistMethods::Poisson::correct(Foam::Geo metricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) at ??:?
#11 void Foam::meshObject::movePoints<Foam::fvMesh>(Foam: bjectRegistry&) at ??:?
#12 Foam::fvMesh::movePoints(Foam::Field<Foam::Vector< double> > const&) at ??:?
#13 Foam::dynamicMotionSolverFvMesh::update() at ??:?
#14 Foam::dynamicOversetFvMesh::update() at ??:?
#15 ? at ??:?
#16 __libc_start_main in "/usr/lib64/libc.so.6"
#17 ? at ??:?
Floating point exception (core dumped)

Additionally, during the solution this problem occurs:

--> FOAM Warning :
From function static Foam::tetIndices Foam:olyMeshTetDecomposition::triangleTetIndices (const Foam:olyMesh&, Foam::label, Foam::label, Foam::label)
in file meshes/polyMesh/polyMeshTetDecomposition/polyMeshTetDecomposition.C at line 618
No base point for face 4, 4(749 748 2665 2666), vertices 4((7.779335446e-07 -1.39959937e-06 1.83904669e-05) (7.779329877e-07 -1.399598368e-06 1.874835471e-05) (-7.779874176e-07 1.399696298e-06 1.000018748) (-7.779868608e-07 1.399695296e-06 1.00001839)), produces a valid tet decomposition.

At the moment I can think of no solution. I have tried to play with fvSchemes file (interpolationScheme and laplacianScheme), fvSolution (the way cellDisplacement is solved), and oversetInterpolation( inverseDistance or cellVolumeWeight) but none of them worked for me. Thanks in advance.

Oguz

[vepe]

Hi Oguz,


it's been a while now that you posted this problem. Unfortunately i have the same problem... Is there any chance that you solved the problem meanwhile?


Thank you in advance!