CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Six DoF - Overset in OpenFOAM

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 7, 2017, 08:21
Question Six DoF - Overset in OpenFOAM
  #1
New Member
 
Oğuzhan Demir
Join Date: Dec 2016
Posts: 19
Rep Power: 7
ddemiroguzhan is on a distinguished road
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
Quote:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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:

Quote:
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:
Quote:
--> 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
ddemiroguzhan is offline   Reply With Quote

Old   January 24, 2019, 05:36
Default
  #2
New Member
 
Join Date: Nov 2018
Posts: 1
Rep Power: 0
vepe is on a distinguished road
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!
vepe is offline   Reply With Quote

Reply

Tags
overpimpledymfoam, overset, sixdof

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Getting Started with OpenFOAM wyldckat OpenFOAM 25 August 14, 2022 13:55
Map of the OpenFOAM Forum - Understanding where to post your questions! wyldckat OpenFOAM 10 September 2, 2021 05:29
Overset Mesh error : no startcell found ! Cris ZC STAR-CCM+ 4 September 9, 2019 10:43
How to contribute to the community of OpenFOAM users and to the OpenFOAM technology wyldckat OpenFOAM 17 November 10, 2017 15:54
OpenFOAM Training Jan-Apr 2017, Virtual, London, Houston, Berlin cfd.direct OpenFOAM Announcements from Other Sources 0 September 21, 2016 11:50


All times are GMT -4. The time now is 10:57.