CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Community Contributions

[OLAFLOW] OLAFOAM baseWaveFlume ,adding slope resulted in Simulation Termination

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By Phicau

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   September 4, 2017, 07:11
Default OLAFOAM baseWaveFlume ,adding slope resulted in Simulation Termination
  #1
New Member
 
Ahmed
Join Date: Aug 2017
Location: Germany
Posts: 2
Rep Power: 0
Siddiki is on a distinguished road
Dear All,
I am new with OpenFOAM(v4.1) and OLAFOAM even with forums, so pardon my irrelevancy if that is so.
I want to simulate waves hitting on a dike + its wave run up so i picked up base wave flume tutorial, it worked fine , then i changed the blockMeshDict by bringing two of the right points closer towards origin to model a slope of 1/5 (which is close to what i need).

The simulation ran for two iterations and then the following statements occurred and the simulation got completed:

blockMesh meshing...
Preparing 0 folder...
Setting the fields...
Running...
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#6 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#7 ? at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? at ??:?
./runCase: line 16: 2570 Floating point exception(core dumped) olaFoam > olaFoam.log
Simulation complete.

I tried checkMesh to re assure and it seemed :

Build : 4.1
Exec : checkMesh
Date : Sep 04 2017
Time : 12:19:21
Host : "siddiqui"
PID : 3137
Case : /home/siddiqui/OpenFOAM/siddiqui-4.1/OLAFOAM/tutorials/OF_OF+_latest/baseWaveFlume
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 17892
internal points: 0
faces: 35195
internal faces: 17305
cells: 8750
faces per cell: 6
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 8750
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
inlet 70 142 ok (non-closed singly connected)
outlet 70 142 ok (non-closed singly connected)
bottom 125 252 ok (non-closed singly connected)
atmosphere 125 252 ok (non-closed singly connected)
frontAndBack 17500 17892 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (0 -0.2 0) (10 0 1)
Mesh has 2 geometric (non-empty/wedge) directions (1 0 1)
Mesh has 2 solution (non-empty) directions (1 0 1)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (4.21952e-18 -3.44198e-15 4.18793e-17) OK.
Max cell openness = 2.15085e-16 OK.
Max aspect ratio = 10.56 OK.
Minimum face area = 0.00057551. Maximum face area = 0.016. Face area magnitudes OK.
Min volume = 0.000115102. Max volume = 0.000227755. Total volume = 1.5. Cell volumes OK.
Mesh non-orthogonality Max: 78.6458 average: 62.4877
*Number of severely non-orthogonal (> 70 degrees) faces: 7784.
Non-orthogonality check OK.
<<Writing 7784 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 1.78278 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

The blockMesh.log is :
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Creating block mesh from
"/home/siddiqui/OpenFOAM/siddiqui-4.1/OLAFOAM/tutorials/OF_OF+_latest/baseWaveFlume/system/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Reading patches section

Creating block mesh topology

Reading physicalType from existing boundary file

Default patch type set to empty

Check topology

Basic statistics
Number of internal faces : 0
Number of boundary faces : 6
Number of defined boundary faces : 6
Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 1
Block 0 cell size :
i : 0.04 .. 0.04
j : 0.2 .. 0.2
k : 0.0142857 .. 0.0142857


Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (0 -0.2 0) (10 0 1)
nPoints: 17892
nCells: 8750
nFaces: 35195
nInternalFaces: 17305
----------------
Patches
----------------
patch 0 (start: 17305 size: 70) name: inlet
patch 1 (start: 17375 size: 70) name: outlet
patch 2 (start: 17445 size: 125) name: bottom
patch 3 (start: 17570 size: 125) name: atmosphere
patch 4 (start: 17695 size: 17500) name: frontAndBack

End

As far as what i get from different sources is that i am having some divide by zero issue somewhere and i don't know how to trace it .
I have not changed anything in alpha.water, p_rgh , U and alpha.water.org
Please help me out with this.Let me know where am i loosing.
Thankyou so much. and sorry for such long post.
The blockTopology is attached.
BlockTopology_Slope.jpg
Siddiki is offline   Reply With Quote

 


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
Help in FSI : termination due to maximum simulation time kqlee CFX 5 August 16, 2012 03:14


All times are GMT -4. The time now is 23:56.