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

Convergence difficulties with 2d mixerGGI

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 3, 2009, 14:24
Default i am no expert of ggi but perh
  #1
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 200
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
i am no expert of ggi but perhaps you can replace p and pfinal brackets by

p ICCG 1e-06 0;
pFinal ICCG 1e-06 0;
louisgag is offline   Reply With Quote

Old   February 4, 2009, 03:55
Default Hi Louis, thanks for your s
  #2
New Member
 
Christina Smuda
Join Date: Mar 2009
Location: Germany
Posts: 12
Rep Power: 8
christinasmuda is on a distinguished road
Hi Louis,

thanks for your suggestion. I changed the solver to ICCG, but the convergence behaviour didn't change at all. Here is my case: rotor2Dggi.rar.

What else could I try in order to improve the convergence?

Regards,
Christina
christinasmuda is offline   Reply With Quote

Old   February 4, 2009, 14:09
Default Hello Christina, I noticed
  #3
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 200
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hello Christina,

I noticed your simulation endtime was 0.01 second. It would be surprising to see the velocity distribution you expect show up at 0.01s considering your rpm is 100.

I have relaxed your tolerances and reduced the nonOrtho correctors and am running your case at the moment. I'll update you on this later during the day.

Cheers,

-Louis
louisgag is offline   Reply With Quote

Old   February 4, 2009, 14:44
Default I was able to get some results
  #4
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 200
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
I was able to get some results that somewhat match the velocity distribution you were looking for by changing endtime to 1s, nCorrs to 0 and nonOrthocorrs to 0. Setting those to 0 might be extreme and you might want to play around with these values, but I set them to zero to speed up calculation and considering your mesh was very orthogonal.

Send me an email if you want the modified case along with the results (5mb).

Good luck,

-Louis


louisgag is offline   Reply With Quote

Old   February 6, 2009, 03:34
Default Hi Louis, thanks for sendin
  #5
New Member
 
Christina Smuda
Join Date: Mar 2009
Location: Germany
Posts: 12
Rep Power: 8
christinasmuda is on a distinguished road
Hi Louis,

thanks for sending the case.

I already tried running the case without any PISO correctors yesterday and I got the same flow field you sent to me. Unfortunately this flow field is wrong, in my opinion. I did the same calculations with fluent and with OpenFoam using sliding interfaces and both times I got the picture, I put in the link.

I'm still not really convinced about the solution you suggested. If I understand the PISO algorithm well, there's no sense in running it without any corrector steps, is there?

But as soon as I set nCorrectors to 1, I get a floating point error quickly. Setting the Courant Number to smaller values (0.25 or 0.1), the calculation is running for a longer time, but still ends up in a floating point error.

What is wrong with my set-up?

Any help is highly appreciated.

Thanks,
Christina
christinasmuda is offline   Reply With Quote

Old   February 6, 2009, 13:05
Default Dear Christina, After sendi
  #6
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 200
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Dear Christina,

After sending you the case yesterday, I also tried to increase the nCorrectors and ran into the same floating point error. It seems that the flow diverges suddenly at some point.

I question whether your mesh could be too small (being 3 mm wide), but checkMesh does not report any cell volume problem...

If I find the problem I'll let you know. Please let me know if you do.

Good luck in the meantime,

-Louis
louisgag 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
Convergence Difficulties sarath CFX 2 August 21, 2006 08:35
Convergence difficulties Freeman FLUENT 1 November 28, 2005 16:42
Post Processing difficulties Ashish CD-adapco 0 August 21, 2003 20:19
CAD difficulties Carlos Saraiva CD-adapco 1 January 9, 2003 12:09
Modeling Difficulties using Pro-Star satish CD-adapco 0 February 29, 2000 07:14


All times are GMT -4. The time now is 11:37.