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

Beginner question OddEven Decoupling in SIMPLE Solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 9, 2005, 01:08
Default As part of my process of learn
  #1
brooksmoses
Guest
 
Posts: n/a
As part of my process of learning how to use OpenFOAM, I decided to try implementing a steady-flow solver with laminar flow and no turbulence. That is, essentially simpleFoam solver, but with the viscosity model and related bits from the icoFoam solver.

As a test case for the solver, I'm using the cavity problem from the first tutorial in the handbook, with the timestep variables adapted for a steady-state problem. Everything seems to be working, up until the point where I start looking at the results. I've got a nasty odd-even (checkerboard) instability going on.

Now, I know that theoretically this sort of thing happens because a central-differenced velocity discretization is decoupled from the pressure value in the center of the cell, and that this is generally best solved by providing some form of upwind biasing. However, this problem arises in my steady-state calculation, but it doesn't arise in the time-based calculations of the example....

Thus, my question: What's present in the icoFoam calculation that takes care of this that I may be missing? What is the best way to take care of the problem in my steady-state calculation? (Or perhaps: Is this something that ought not be happening in the first place, and I'd be better off looking for a bug?)

Thanks!
  Reply With Quote

Old   November 9, 2005, 07:30
Default What's present in the icoFoam
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Quote:
What's present in the icoFoam calculation that takes care of this that I may be missing? What is the best way to take care of the problem in my steady-state calculation?
You need some form of Rhie-Chow interpolation - have a look at the formulation of the pressure equation in icoFoam.

If you want to make a steady-state version of icoFoam, all you need to do is to get rid of the time derivative in the momentum equation and add implicit under-relaxation. Additionaly, you want SIMPLE instead of PISO (solve the pressure only once and under-relax it explicitly) and you're done.

Enjoy,

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   November 9, 2005, 20:38
Default Thanks for the advice! That s
  #3
brooksmoses
Guest
 
Posts: n/a
Thanks for the advice! That sounded like almost exactly what I'd done, so I went through and compared each line of my case to the icoFoam code, and made sure things matched up properly.

I'm embarassed to say what the problem was: I'd made a sign error on the viscosity term!

- Brooks
  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
question from ICEM beginner Jiuan CFX 1 March 6, 2009 19:41
How to write a solver beginner ivan_cozza OpenFOAM Running, Solving & CFD 23 January 31, 2008 16:17
LES beginner question Shuo Main CFD Forum 4 July 9, 2007 08:40
Question from a beginner Arun K Main CFD Forum 4 August 13, 2004 07:17
simple heat pipe flow-beginner mcm FLUENT 0 February 3, 2003 21:31


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