# Solution convergence problem - poor mesh?

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

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.78e-3 s. Problem with this is that using the adjustTimeStep in controlDict to keep maxCo low gives me timesteps on the order of 1e-7. 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 5-sided cells has anything to do with the problem, or possibly my boundary layer, which has 50 inflated layers. Turbulence model is k-w 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
Attached Images
 stator.jpg (89.8 KB, 27 views) statorggi.png (33.0 KB, 21 views) rotor.jpg (96.9 KB, 39 views) rotorblade.jpg (94.6 KB, 25 views) rotorbladele.jpg (79.6 KB, 23 views)

 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.6-ext>)?
- Anyone have any experience with the new rotatingMotion in OF2.0? Perhaps I should try that approach.

Thanks,
Shawn
Attached Images
 B2_LE.jpg (19.8 KB, 20 views)

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post RameshK Main CFD Forum 7 July 17, 2010 14:13 Bedotto NUMECA 1 March 18, 2010 00:40 geryes FLUENT 0 February 25, 2010 17:32 hammam CFX 3 August 7, 2007 18:33 limseokmin FLUENT 3 November 14, 2004 13:43