CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

AMI simulation unstable within first few timesteps

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 24, 2014, 12:09
Default AMI simulation unstable within first few timesteps
New Member
A. Craig
Join Date: Feb 2014
Posts: 2
Rep Power: 0
craig0a is on a distinguished road
Hello everyone,

There are a couple of threads on this topic already, but after looking through them carefully, I am still not sure how to go about fixing my particular problem, since my simulation seems to be breaking down in a slightly different way .

Short: With my AMI set to not be rotating, my patch source/target weights are close to 1 when running, but my simulation is going unstable within a few timesteps. In Paraview, it looks like there is a build up of pressure/high velocity at the upstream surface of my AMI.


I am simulating a vertical axis wind turbine in 3D (again, like so many on this forum!). The geometry was made in Solidworks, modified in Blender to ensure the mesh was closed and to refine the triangulation, and exported as an ascii .obj file. In my case, I have three blades, a tower, and an "innerCylinderSmall" which encompasses just the blades, brushing the top surface of the tower.

I then followed the propellor tutorial as closely as possible. So my command flow is:

topoSet -dict system/removeRedundantZones.topoSetDict
topoSet -dict system/createInletOutletSets.topoSetDict
topoSet -dict system/createAMIFaces.topoSetDict
createPatch -overwrite
createBaffles -overwrite
mergeOrSplitBaffles -split -overwrite
transformPoints -scale "(0.00508 0.00508 0.00508)"

If I am understanding correctly, and this seems to be supported by my results, the "innerCylinderSmall" should be used in snappyHexMesh to create a more refined grid around my blades, but then it is removed and replaced with AMI surfaces. (So my AMI surface is a closed cylinder.)

Checking my results in ParaView, the mesh looks okay (I need to switch to an unstructured grid near the blades, but that can come later). The geometries seem to match up alright; ParaView recognizes both AMI1 and AMI2 as patches and they overlap exactly.

Running checkMesh, I get "Mesh OK". The cell sizes are about two orders of magnitude smaller than that used in the propeller tutorial. Again, AMI1 and AMI2 are recognized, each with 110444 faces, albeit one has 110682 points while the other has 110666 points.

When I try to run the simulation, setting the rotation to zero, running in water (like the propeller) at a Reynolds number of 600, for the first ~100 steps residuals, Courant number, patch/source target weights, etc. seem to be doing fine, but then in about 5 steps, everything explodes. (The number of steps to blow up seems to be the same regardless of initial time step size and I have been trying two orders of magnitude smaller than the propeller tutorial to match my smaller cell size.) Looking at the time sequence in Paraview, I am getting effectively unbounded growth in velocity and pressure at the upstream edge of my AMI surface (actually between blades of my geometry.)

I am guessing that this has something to do with how I am defining/meshing/setting boundary conditions on my AMI surface, but I don't know how to go about debugging this.

If anyone has the time to help me through this, I would greatly appreciate it. At present, I don't know what files/screenshots may be pertinent, so if you can help, please let me know what you would like to see.

Thank you!

A. Craig
craig0a is offline   Reply With Quote


ami, unstable

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
Problem with maximum number of Timesteps Phantalor CFX 5 July 13, 2013 07:15
Requiring some help regarding tornado simulation amarjogot Main CFD Forum 0 June 20, 2013 12:26
AMI simulation with VOF Ya_Squall2010 OpenFOAM Running, Solving & CFD 29 May 13, 2013 03:19
GUI crash and simulation engine still running RPJones FLOW-3D 2 November 9, 2010 08:18
3-D Contaminant Dispersal Simulation Apple L S Chan Main CFD Forum 1 December 23, 1998 10:06

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