
[Sponsors] 
November 10, 2011, 21:08 
Solution convergence problem  poor mesh?

#1 
Member
Shawn
Join Date: Oct 2011
Posts: 56
Rep Power: 7 
I'm trying to figure out why I can't converge my solution  yes, I know, lots of possibilities here. I'm wondering if it has to do with my mesh. My simulation is a vertical axis turbine, using the GGI interface and pimpleDyMFoam. WI previously used CFX to do the simulations, and found appropriate grid size, dimensions, boundary layers, timestep so that incremental changes in the simulations parameters changed the solution by about ~1%. To do 3D cases, I need to move to openFoam so that I can run on more processors. In order to compare the inputs (timestep, grid, etc.), I need to compare a 2D openFoam solution to a 2D CFX solution. I got a small (smaller, simpler geometry with quite coarse mesh) test case running in openFoam and have moved to my full sized mesh and am running into convergence/timestep problems. I would like to be able to run my cases with timesteps equal to about 1 degree of rotation, this is based on CFX results giving about a 1.5%2.0% change in both power coefficient and maximum thrust coefficient between steps equivalent to 2.5 degrees, 1.0 degrees and 0.5 degrees. For my first test, this corresponds to a timestep of 2.78e3 s. Problem with this is that using the adjustTimeStep in controlDict to keep maxCo low gives me timesteps on the order of 1e7. On top of that, I have a lot of numerical instability so that my velocity (inlet conditon is about 10 m/s) climbs past 10e6, but seems to come back down, slowly, eventually. I set up my grid so that I would have 360 equal divisions on my ggi interface, one cell incrementation per "ideal" timestep. I was in contact with someone else who did a similar simulation, linnemann, and he used a much different grid  structured and hex. I'm wondering if my unstructured, grid and 5sided cells has anything to do with the problem, or possibly my boundary layer, which has 50 inflated layers. Turbulence model is kw SST. I ran his case which seemed to converge well, even with a fairly high (50) limite on maxCo. I've tried increasing the number of corrector steps (nCorrectors, nOuterCorrectors, nNonOrthoganalCorrectors), relaxation, solver method (GAMG, PCG, ...). Any tips are welcome.
Thanks 

November 10, 2011, 22:47 

#2 
Member
Shawn
Join Date: Oct 2011
Posts: 56
Rep Power: 7 
I should add that the linnemann mesh was much more coarse around the blaedes (300 nodes vs 1000 nodes) with fewer boundary layers (5 vs 50).
I will post my controlDict/fvSolution/fvSchemes files a little later, had to leave for other work. I was also wondering if maybe I should start with a simpleSolver and get a steady state solution at my initial blade position to use as initial conditions, since the very first iteration of my transient solution seems to generate ridiculously higher velocity pockets/pressure near the blade surfaces. I have read that sometimes there can be some serious sensitivity to the initial conditons. I suppose it's worth a try. Shawn 

November 12, 2011, 20:31 

#3 
Member
Shawn
Join Date: Oct 2011
Posts: 56
Rep Power: 7 
Turns out that I just get high velocity (and gradient) at the leading edge of the blade. Given the high velocity, the segregated PISO solver and the small cell size, it's very easy for my Courant number to rise unless the timesteps are VERY small.
Increasing the number of nOuterCorrectors seemed to help, as did changing my scheme to upwind. I'm going to experiment with the number of various correctors and the schemes a little more.  In lieu of that, are there any coupled/fully implicit solvers available that allow for a rotor/stator interface (similar to my using pimpleDyMFoam with the ggi <OF1.6ext>)?  Anyone have any experience with the new rotatingMotion in OF2.0? Perhaps I should try that approach. Thanks, Shawn 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
First order in time and Central Difference Convergence problem  RameshK  Main CFD Forum  7  July 17, 2010 14:13 
History Convergence: Graphical problem  Bedotto  NUMECA  1  March 18, 2010 00:40 
How to use a solution as an inlet boundary condition in another problem in fluent  geryes  FLUENT  0  February 25, 2010 17:32 
transition model convergence problem  hammam  CFX  3  August 7, 2007 18:33 
convergence problem  limseokmin  FLUENT  3  November 14, 2004 13:43 