# Beginner question OddEven Decoupling in SIMPLE Solver

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

 November 9, 2005, 01:08 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!

November 9, 2005, 07:30
What's present in the icoFoam
#2
Senior Member

Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,762
Rep Power: 21
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

 November 9, 2005, 20:38 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

 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 Jiuan CFX 1 March 6, 2009 19:41 ivan_cozza OpenFOAM Running, Solving & CFD 23 January 31, 2008 16:17 Shuo Main CFD Forum 4 July 9, 2007 08:40 Arun K Main CFD Forum 4 August 13, 2004 07:17 mcm FLUENT 0 February 3, 2003 21:31

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