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

pimpleDyMFoam for VIV of a cylinder using LES

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

Like Tree1Likes
  • 1 Post By kkpal

Reply
 
LinkBack Thread Tools Display Modes
Old   October 9, 2015, 09:19
Default pimpleDyMFoam for VIV of a cylinder using LES
  #1
Senior Member
 
Join Date: Jan 2013
Posts: 118
Rep Power: 5
kkpal is on a distinguished road
dear foamers,
I am now running cases about flow induced vibration of a circular cylinder at Re=5000 with LES model. The solver is pimpleDyMFoam.
It runs well for the cases in which the amplitude of the cylinder is small. However, when it comes to the large amplitude cases, i.e., lock-in region, the following error message pumps out suddenly (by suddenly I mean the previous time steps are all good in terms of Cd, Cl, displacement, maxCo and so on ) and terminates the simulation.
Code:
   Courant Number mean: 0.0407862 max: 2.50865
Time = 109.01


Restraint verticalSpring:  attachmentPt - anchor (0 -0.257633 0) spring length 0.257633 force (-0 3.25133 -0)
Centre of mass: (0 -0.257633 2)
Linear velocity: (0 -0.849337 0)
Angular velocity: (2.0932e-06 -1.04789e-07 -2.56901e-08)
[4] processorPolyPatch::calcGeometry : Writing my 25 faces to OBJ file "/work/Re5000VIV/middle/processor4/procBoundary4to3throughfront_faces.obj"
[4] processorPolyPatch::calcGeometry : Dumping cell centre lines between corresponding face centres to OBJ file"/work/Re5000VIV/middle/processor4/procBoundary4to3throughfront_faceCentresConnections.obj"
[3] processorPolyPatch::calcGeometry : Writing my 25 faces to OBJ file "/work/Re5000VIV/middle/processor3/procBoundary3to4throughback_faces.obj"
[3] processorPolyPatch::calcGeometry : Dumping cell centre lines between corresponding face centres to OBJ file"/work/Re5000VIV/middle/processor3/procBoundary3to4throughback_faceCentresConnections.obj"
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 4 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
What does this error message mean?
kkpal is offline   Reply With Quote

Old   October 10, 2015, 09:52
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,748
Blog Entries: 39
Rep Power: 103
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quick answer: Looks like a physically impossible situation has occurred, probably where the mesh broke inside the cylinder.
In addition, this:
Code:
Courant Number mean: 0.0407862 max: 2.50865
is a dangerous situation, since this roughly means that the mass content of one cell can travel to 2-3 cells in a single time step, which would possibly explain why the mesh collapsed into itself.

When in doubt, save more time snapshots near the time of the crash, so that you can visually diagnose what the mesh looks like in each time step.
wyldckat is offline   Reply With Quote

Old   October 15, 2015, 21:20
Default
  #3
Member
 
Davi Barreira
Join Date: Apr 2014
Location: Fortaleza
Posts: 76
Rep Power: 4
davibarreira is on a distinguished road
Like Bruno said, it's probably that your mesh deformed too much. Look in Paraview in "surface with edges" mode, and see how the mesh cells are actually deforming. To solve this, you might increase the "outerDistance" in your dynamicMeshDict
davibarreira is offline   Reply With Quote

Old   October 17, 2015, 03:26
Default
  #4
Senior Member
 
Join Date: Jan 2013
Posts: 118
Rep Power: 5
kkpal is on a distinguished road
Dear Davi,
I think it might be associated with the mpi communications. I reconstructed the case at the last written time and run the case again with single processor, the same simulation run smoothly past the time at which mpirun stopped.
Besides, when I change the span-wise boundary condition of the computation domain from cyclic to slip, and run in parallel, the problem did not show up.

Quote:
Originally Posted by davibarreira View Post
Like Bruno said, it's probably that your mesh deformed too much. Look in Paraview in "surface with edges" mode, and see how the mesh cells are actually deforming. To solve this, you might increase the "outerDistance" in your dynamicMeshDict
davibarreira likes this.
kkpal is offline   Reply With Quote

Old   October 17, 2015, 10:34
Default
  #5
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,748
Blog Entries: 39
Rep Power: 103
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quick answer: From a blog post of mine: Notes about running OpenFOAM in parallel - see this post:
Quote:
stop when I run in parallel - post #28 provides a list of known reasons why a run might not go well in parallel
wyldckat is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Flow past a cylinder at Re 1e05 using LES, drag force coefficient is to low Scabbard Main CFD Forum 19 June 17, 2016 09:53
pimpleDyMFoam: moving cylinder in duct - crashes with high Courant Number khatereh OpenFOAM Running, Solving & CFD 1 July 18, 2015 12:54
adapt offset cylinder model for les turbulence shackman287 OpenFOAM 0 August 19, 2010 21:10
LES of a square cylinder gfilip OpenFOAM Running, Solving & CFD 1 June 24, 2010 12:33
3D LES simulation of a circular cylinder fpz Main CFD Forum 3 June 20, 2005 19:50


All times are GMT -4. The time now is 07:45.